Effect of changes in the fractal structure of a littoral zone in the course of lake succession on the abundance, body size sequence and biomass of beetles

Study area and field studies

The study covered 40 dystrophic lakes located in northern Poland: the South Baltic Coastland, the West Pomeranian, Olsztyn, Mrągowo Lakelands and the Suwałki Lakeland (Fig. 1, Appendix S1). The lakes chosen for the study were of various surface areas and having a floating peat mat growth belonging to Sphagnum sp. that varied in size. The lakes represented various succession stages in a disharmonic series—from oligo- to polyhumic lakes. The oligohumic lakes were inhabited by Juncus bulbosus, Eleocharis palustris, Phragmites australis, Typha angustifolia, Typha latiffolia, Sphagnum sp., Lobelia dortmanna, Isoëtes lacustris, Drosera rotundifolia, D. anglica, D. intermedia, Sparganium angustifolium, Lycopodiella inundata. In polyhumic lakes, the dominant species included Sphagnum sp., Oxycoccus quadripelatus, Andromeda polifonia, Eriophorum vaginatum, Calluna vulgaris, Erica tetralix, Empetrum nigrum. The lakes were arranged a priori according to the seral stages of succession (Bloechl et al., 2010), and the point of reference was assumed to be the coverage of a lake’s surface by a Sphagnum mat (0–75%) (Appendix S1). The percentage cover by a Sphagnum mat in each lake was calculated in the GIS system, supported by ArcMap 9.3.1 software. For charting the research sites, data made available through Geoportal 2 in the WMS format were employed.

The study was carried out from 2002 to 2014, in spring, summer and autumn. The fauna samples were collected with a dip net from a surface area of about 1 m². In the compacted Sphagnum mat environment, a sample comprised 10 subsamples (aggregated afterwards), which were collected using a 0.1 m² strainer. Sampling sites were chosen so as to represent the biggest array of the littoral habitat diversity and areas of individual lakes. Thus, five different littoral components (habitats) were identified: (1) the Sphagnum mat, (2) sparse macrophyte zone, (3) dense macrophyte zone, (4) sandy bottom zone and (5) pockets and ponds within a Sphagnum mat. The vegetation cover was assessed using Braun-Blanquet (1964) phytosociological records. All lakes were identified with respect to the surface area, extent of the lake’s surface covered with a Sphagnum mat, which corresponds to succession stages (1–3), percent share of individual habitats in the littoral zone (1–4) and a seral stage (1–40) (Appendix S1). It was assumed that these habitats, due to different degrees of complexity, represent different fractal dimensions, all within an interval of 2.0 >D >3.0 (Tokeshi & Arakaki, 2012). A schematic representation of fractals different in size can be seen in Fig. 2. The smallest size (1) corresponds to the sandy bottom zone, while the largest one (5) is assigned to the Sphagnum mat. The intermediate sizes belong to: (2) pockets and ponds within the Sphagnum mat, (3) sparse macrophyte zone, and (4) dense macrophyte zone. Field studies were conducted in accordance with field-study-approvals: Certificate of RDOS in Olsztyn (WOPN.070.21.2018.AKI), WNP Certificate (PNE 08-070/ 18) and ZPK Certificate (060.1.2018).

In total, 499 samples were collected. Subsequently, the collected samples were described according to the chosen environmental parameters: habitat (1—Sphagnum mat, 2—sparse macrophytes, 3—dense macrophytes, 4—sandy bottom, 5—pockets and ponds within Sphagnum mat), place (location) (1–ecotone, 2—pockets and ponds within the Sphagnum mat, 3—compacted Sphagnum mat), bottom (soil substrate) (1—sand, 2—sand and mud, 3—Sphagnum), depth (1—0–10 cm, 2—0–20 cm, 3—20–40 cm, 4—40–60 cm), nympheids (1—none, 2—present), elodeids (1—none, 2—present), (detritus (1—none, 2—scarce and fine, 3—abundant and fine, 4—abundant and coarse), debris (1 –none, 2—present), fractal dimension (1–5), stage—seral stage of succession (1–40).

Ecological and statistical analyses

Species diversity was calculated using: S—number of species, N—number of individuals and D—percent share. All collected beetles were divided into five classes of different body size (where a body size meant the total length of an organisms), i.e., 1—very small beetles (<3.0 mm), 2—small beetles (3.0–5.0 mm), 3—medium beetles (5.1–10.0 mm), 4—large beetles (10.1–20.0 mm) and 5—very large beetles (>20.0 mm). In addition, the body weight was measured within each of the five classes. The ranges of body weight values of beetles in the distinguished body size classes are shown in Table 1. Three functional groups were distinguished in the trophic structure of beetles: predators (families: Gyrinidae, Dytiscidae and Noteridae), polyphages (Haliplidae) and saprophages (Helophoridae, Hydrophilidae and Hydraenidae) (Appendix S2). Collected adult water beetles were identified using standard identification keys (Galewski, 1990; Galewski & Tranda, 1978; Friday, 1988; Hansen, 1987; Holmen, 1987).

Because the samples of water beetle fauna were collected several times (taking into account the phenological aspect) from the same lakes and from the habitats distinguished within these lakes, we used a GLM (Generalized Linear Model) for repeated measures ANOVA (Hocking, 1996) to determine the significance of differences in the number of beetles within each functional group in the distinguished habitats, and to identify dependencies between the type of habitats, body size, abundance and biodiversity of Coleoptera. First, we checked the assumptions of normality (the Shapiro–Wilk test) and homogeneity of variances (the Levene’s test), respectively. The GLM repeated measure models were calculated on the basis of Type III sums of squares so as to take the unbalanced design into account. Significant results were tested for pair-wise comparisons by the Tukey’s HSD post-hoc tests. Dependent variables (abundance and number of species) were transformed where necessary to fulfil the requirements of parametric tests (Saint-Germain et al., 2007; Cremona, Planas & Lucotte, 2008).

In four separate analyses, we calculated the Pearson’s correlation coefficients (r_p) in order to determine correlations between: (1) body size and counts of beetles in the distinguished habitats, (2) cover of the Sphagnum mat and counts of beetles, (3) cover of the Sphagnum mat and functional groups, and (4) species and environmental variables.

Linear regression analysis was performed to determine the influence of: (1) the body size on the abundance of beetles, (2) fractal dimension on the mean body size on, (3) fractal dimension on the mean weight, and (4) the seral stage of succession on the mean body size. Multidimensional correspondence analysis (MCA) (Clausen, 1998) served to determine dependencies between the abundance of water beetles within the identified body size classes, the share of the distinguished habitats in the littoral zones and in the analysed lakes. MCA has been used in similar studies, e.g., Usseglio-Polatera, Bournaud & Tachet, 2000; Obolewski, Glińska-Lewczuk & Kobus, 2009. The average density of body size classes in a sample was adjusted by the value corresponding to the total contribution of a given littoral zone type in individual lakes, thus giving the weighted average. The analysis included two dimensions, of which the first explained the biggest part of the general chi-squared statistics (% of inertia), whereas the inclusion of the other dimension increased the percentage of explained inertia. Similarities in the fauna between particular lakes were also analysed by the single linkage method for the hierarchical clustering of objects. The distance measure is an Euclidean distance. The results are presented in the form of a dendrogram drawn for 40 lakes.

Relationships between the presence of beetles and the analysed environmental parameters in individual lakes were determined with the help of a principal component analysis, PCA. All calculations were performed in Statistica, ver. 13.0 (StatSoft, Tulsa, USA).

General description of the collected material

The collected material comprised 10,139 specimens representing 124 species classified to seven families (Gyrinidae, Haliplidae, Noteridae, Dytiscidae, Hydraenidae, Hydrochidae and Hydrophilidae) (Appendix S1). Most species (55.6%) belonged to the two smallest body size classes. In total, they composed 84.75% of all collected specimens. The remaining beetles represented the consecutively higher body size classes, of which the largest one had the fewest representatives (Table 1). Most beetles (28.8% of the entire material) were collected in pockets and ponds of the Sphagnum mat. It was also there that the highest species richness was determined (91 spp.). The species diversity was lower in the dense macrophytes zone (26.6%; 74 spp.), sparse macrophytes zone (19.64%; 66 spp.), compacted Sphagnum mat (18.9%; 64 spp.), and finally in the sandy bottom (5.9%; 39 spp.). Very small beetles, which represented the first body size class, demonstrated the highest species richness in the compacted Sphagnum mat (41.5% of the whole collected material) (Fig. 3), while being the most numerous in the dense macrophyte zone. Small beetles (class 2 of the body size) were most numerously caught in the pockets and ponds within the Sphagnum mat (73.4%). The presence of very large beetles was notable in the sparse macrophytes zone, pockets and ponds in the Sphagnum mat and in the sandy bottom zone, although they made up a small share of the collected material (around 1%).

In the trophic structure, the highest percentage (71.9%) was achieved by predatory beetles, represented mainly by Dytiscidae (Fig. 4). Although their species diversity was similar in all the habitats, they were most numerous in the dense macropytes zone (91.3%) and in the sandy bottom zone (84.5%). Most saprophages (47.2%) were noted in the compacted Sphagnum mat, while the number of species representing polyphages was low in all of the habitats, with the highest score (six spp.) determined in the compacted Sphagnum mat. The GLM Repeated Measure Anova results (F_(2,499) = 6.74, p = 0.0089) indicate significant differences in counts of the distinguished trophic groups in the particular habitats. The counts of predatory beetles in the habitats in question did not differ significantly (the Tukey’s HSD test, p = 0.18). Conversely, differences in the counts of predators and polyphages were significant (p = 0.0003).

Dependencies between types of habitats, body size and abundance as well as biodiversity of Coleoptera

The results of the GLM repeated measure Anova analysis point to a significant influence of the synergistic effect between the type of habitat and body size class on the abundance of beetles (F_(16,1996) = 7.25, p = 0.0000002). Significant differences (the Tukey’s HSD post-hoc test) between the analysed subclasses are illustrated in Fig. 5A (cf. Appendix S3). Moreover, a significant influence of the synergistic effect was determined between the indicated factors (habitat and body size class) and the number of species (F_(16,1996) = 16.22, p = 0.00002) (Appendix S3). Results of the GLM repeated measure also point to a significant influence of the type of habitat on the number of beetles F_(4,499) = 6.75, p = 0.00003 (Fig. 5B). Significant differences (the Tukey’s HSD post-hoc test) were determined between the compacted Sphagnum mat and sparse macrophytes zone (p = 0.031), dense macrophytes zone and sparse macrophytes zone (p = 0.00014), dense macrophytes zone and pockets (p = 0.0479). Moreover, significant differences were demonstrated in the numbers of beetles representing the distinguished body size classes (F_(4,1996) = 77.61, p = 0.00002). The Tukey’s post-hoc test reveals significant differences in the counts of beetles representing all the classes of body size (p = 0.00002), except 3 and 4 class (p = 0.77). Both factors, the type of habitat (F_(4,499) = 11.04, p = 0.000003) and body size (F_(4,1996) = 375.61, p = 0.000001), also had a significant impact on the number of species determined. With respect to the species diversity, significant differences were found between the Sphagnum mat and sparse macrophytes (p = 0.0048), and between the pockets and ponds in the Sphagnum mat and sandy bottom zone (p = 0.0046), and between the Sphagnum mat and sparse macrophytes (p = 0.00018) and dense macrophytes (p = 0.0014). Significant differences in the number of species were shown between all body size classes (p = 0.00002), except 1 and 3 (p = 0.24).

Negative correlations were observed in all the habitats between the body size and the number of beetles. It was only in the sandy bottom zone that these correlations were not significant (p = 0.47) (Table 2). Similarly, negative correlations appeared between the body size and species richness in the compacted Sphagnum mat (p = 0.01) and in the sparse macrophytes zone (p = 0.0008) (Table 2).

Dependencies between the type of habitats, body size and abundance of Coleoptera in lakes

Our analysis of the body size structure of beetles in the fauna populating every lake demonstrated significant differences in the counts of beetles representing different body size classes (x² = 171.18, df = 156; p = 0.0001). The plotted diagram points to lakes with a more strongly developed Sphagnum mat concentrating small organisms (two smallest body size classes) (Fig. 6). In the dendrogram (Fig. 7), these lakes create centrally located clusters (2–7). In the lakes less densely overgrown with Sphagnum sp., with the littoral zone either bare or weakly overgrown by sparse macrophytes, a larger contribution of organisms that belonged to the other body size classes, especially middle-size ones, was notable (clusters 8–9) (Figs. 6 and 7). Lakes included in cluster 1 are characterized by a minimal share of the smallest beetles (prevalent are the beetles that belong to the second body size class). The lakes Żabionek and Wielkie Gacno were distinguished by a demonstrably large contribution of big beetles.

Presence of beetles in lakes versus the fractal structure

Our analysis of the correlations confirmed a strong positive correlation between the development of a Sphagnum mat in the littoral zone of each lake and the abundance of beetles which belong to the first body size class (r_p = 0.7025, p = 0.000013), while the correlation with the number of the largest beetles was negative (r_p =  − 0.21, p = 0.04). Both releationships were linear: LOG (Body size 1) = 0.9045 + 0.4736 * LOG (Sphagnum mat%); LOG (Body size 5) = 0.1531–0.101 * LOG (Sphagnum mat%). The remaining correlations were not significant statistically (Table 3).

A negative relationship between the size of body and the number of organisms in lakes emerged from our analysis conducted on the level of samples (r_p = 0.7046, p = 0.00001). This correlation was linear as well (Fig. 8A). The relationship between the weight and abundance of beetles was similar in character (Fig. 8B).

The results of the linear regression analysis showed that as the fractal dimension increased, the mean body size (r_P =  − 0.86, p = 0.037) and the mean weight of Coleoptera (r_P =  − 0.91, p = 0.032) decreased (Fig. 9). In turn, the relationship between the total biomass of beetles and the fractal dimension was not significant (r_P =  − 0.31, p = 0.623).

Our analysis performed for all the lakes (arranged according to the seral stages of succession), which took into account the extent of the Sphagnum mat cover and the shares of the other habitats with different fractal dimensions, also indicated negative correlations between the seral stage of succession and the mean body size (r_P =  − 0.465, p = 0.0034) (Fig. 10).

The principal component analysis of the variables representing the parameters of habitats and trophic groups distinguished in our study (Fig. 11) suggests that the first axis corresponding to the highest own values most strongly corresponds with the variables ‘predators’ and ‘saprophages’, while the second axis shows the strongest correlations with the variables ‘body size’, ‘weight’, ‘N—abundance’, ‘place’, ‘depth’, ‘stage’ and ‘fractal dimension’. Positive correlations were determined between ‘abundance’ and fractal dimension (r_p = 0.26, p = 0.003), ‘abundance’ and ‘stage’ (r_p = 0.06, p = 0.001), ‘saprophages’ and ‘detritus’ (r_p = 0.75, p = 0.004), ‘polyphages’ and ‘fractal dimension’ (r_p = 0.81, p = 0.003), while negative correlations appeared between ‘abundance’ and ‘depth’ (r_p =  − 0.41, p = 0.003), ‘saprophages’ and ‘depth’ (r_p =  − 0.85, p = 0.001), ‘abundance’ and ‘weight’ and ‘body size’ (r_p =  − 0.15, p = 0.008) and between ‘predators” and place (r_p =  − 0.85, p = 0.004) and between ‘abundance’ and ‘nympheids’ (r_p =  − 0.78, p = 0.006).

Influence of the habitat and architecture on the richness, abundance and body size of Coleoptera

This article is an outcome of 12 years of field studies, during which time we acquired a very rich amount of comparative material (124 species) from 40 lakes. Among the numerous references dedicated to water beetles, there are few which document a larger number of lakes (e.g., Nilsson & Söderberg, 1996; Buczyński, 2012).

In our study, same as Linzmeier & Ribeiro-Costa (2011), we observed the highest species diversity and density of beetles representing two smallest body size classes, mostly of the genus Hydroporus. According to Ulrich (2007), empirical studies show the highest species richness among the medium-size species, while simulation models point to the smallest species. A contrary opinon is held by Scheffer & Van Nes (2006) and Scheffer et al. (2015), who maintain that intermediate body sizes in animal assemblages are found extremely rarely because coevolution in competitive systems favours the co-occurence of species either very similar to or very different from one another.

Thus, it seems to be the matter of assigning the criteria distinguishing body size categories. At the same time, we demonstrated very scarce presence of the largest species, a finding also reported by Nilsson & Söderberg (1996), Ulrich (2007) and Linzmeier & Ribeiro-Costa (2011).

Similarly to Nilsson & Söderberg (1996), Saint-Germain et al. (2007) and Linzmeier & Ribeiro-Costa (2011), we revealed negative correlations between the body size (and the individual body mass) and the abundance of beetles, as well as between the body size and the species richness. According to Nilsson & Söderberg (1996), a negative relationship between the abundance and body size is typical for most aquatic organisms. On the other hand, Saint-Germain et al. (2007) emphasise the negative correlation between the abundance of species and body size.

For clarification, we made an in-depth analysis of the dependencies between the body size and the number of beetles relative to the aspects defining the complex nature of each habitat. Such an approach enjoys a long tradition, as suggested by Ulrich (2007), although it is rarely implemented in hydrobiological research (Cremona, Planas & Lucotte, 2008; Tokeshi & Arakaki, 2012; Désamoré et al., 2018). We demonstrated significant differences between the abundance and species richness of beetles representing different body size classes in the particular habitats. Figure 5 shows that the highest abundance and species richness occurred in the Sphagnum mat and mostly with respect to the smallest beetles (body size classes 1 and 2), mainly of the genera Anacaena and Hydroporus. In the same habitat, the least numerous were large beetles and the largest ones, such as Dytiscus or Cybister, were completely absent. A reverse situation was observed in the sandy bottom zone and in the pockets and ponds within the Sphagnum mat, where the smallest beetles (class 1) were the least numerous, while the largest ones were more abundant than elsewhere. Medium-size beetles (classes 2 and 3), e.g., Noterus, Agabus, and Ilybius, were numerous in all habitats. These relationships are explained perfectly well by the fractal structure. Should we take into consideration an increase in the fractal dimension from the smallest one (the less complex form of a habitat) in the sandy bottom zone to the highest (the more complex form of a habitat with increasingly small structural elements) in the Sphagnum mat, then our observations are in accord with the ones reported by other hydrobiologists, who conclude that the species diversity and density increase as the fractal dimensions increases (e.g., Verberk et al., 2001; Tokeshi & Arakaki, 2012). Meanwhile, small spaces between leaves in a Sphagnum mat inhibit the presence of very large beetles, a conclusion supported by the results provided by Tokeshi & Arakaki (2012). Similar conclusions were drawn by Scheffer et al. (2015), who suggest that the maximum body size is limited by the available space in which beetles could move.

In our study, negative correlations between the body size, weight or abundance of water beetles were determined in particular habitats. Considering the total biomass of beetles in individual habitats, we were unable to identify any significant differences, which again agrees with the conclusions drawn by Tokeshi & Arakaki (2012), namely that biomass does not change in the fractal dimensions. An increase in the density of beetles is offset by a decrease in the individual size of the body, which is concordant with the results of Tokeshi & Arakaki’s research (2012).

A very compact structure of vegetation creates niches for small beetles, offering them egg-laying safety, food for larvae (much humus) and a shelter from predators (Verberk et al., 2001). How is the co-occurrence of these smallest beetles possible if, according to Scheffer et al. (2015), specimens of the same size compete with one another most strongly? Some hydrobiologists, e.g., Scheffer & Van Nes (2006) and Scheffer et al. (2015), draw attention to the fact that these organisms create self-organising assemblages, the presence of which is the result of co-evolution between potential competitors.

In the trophic structure we analysed we discovered that 70% of beetles were predators, being the most numerous groups in all habitats. According to Bloechl et al. (2010), the number of predatory beetles depends on the amount of their prey. In our study, the smallest quantitative impact of predatory beetles (same as the number of saprophages) was identified in the Sphagnum mat. There they represented classes of organisms with small body size, e.g., Hydroporus, which do not limit the presence of other water beetles, as they use other food resources, e.g., zooplankton, smaller insects, like mayflies, or insect eggs (Frelik & Pakulnicka, 2015; Frelik, Koszałka & Pakulnicka, 2016; Perissinotto, Bird & Bilton, 2016). According to Scheffer & Van Nes (2006), Pakulnicka et al. (2013) and Scheffer et al. (2015), what happens here is the evident division of functions, which relies on the principle of minimising similarity (being similar albeit different), as this minimizes competitiveness. As Scheffer et al. (2015) maintain, these mechanisms ensure a certain measure of redundancy of similar species in the environment, which is essential for the functioning of ecosystems during unfavourable changes.

Implications of the disharmonic succession of lakes on coleopteran fauna

The species diversity and density of organisms are the properties most often applied to measuring both biodiversity and quality of the environment, also in the context of succession changes (Kowalik, 1968; Ranta, 1985; Kordylas, 1990; Linzmeier & Ribeiro-Costa, 2011). According to Vamosi, Naydani & Vamosi (2007), species richness and body size change predictably along the spatial gradient, whereas Linzmeier & Ribeiro-Costa (2011) suggest that changes in the body size can be an important indicator, while Saint-Germain et al. (2007) claim that body mass is a better measure than abundance.

Lake succession is a phenomenon that occurs in a given place and over a certain period of time, and with time (in subsequent stages) habitat conditions, mostly shaped by macrophytes, change as well. For particular organisms, optimal living conditions appear here and now, and then they disappear. Thus, changes in the primary structure of lakes result in changes in the secondary structure, formed by assemblages of various organisms. Many ecologists, including Siemann, Haarstad & Tilman (1999), Brown (2003), White et al. (2007) and Linzmeier & Ribeiro-Costa (2011), underline that body size is correlated with numerous morphological, physiological and ecological features, such as an ability to disperse, metabolic capacity, digestion capacity, reproduction rate and duration of a generation, as well as biodiversity. Scheffer & Van Nes (2006) claim that these characteristics act as specific mechanisms which prevent competition between species, especially ones with similar body size.

We have demonstrated that body size and body mass are also useful measurements in investigations into the mechanisms of succession of beetles in dystrophic lakes. Communities of lake-dwelling beetles representing different seral stages are characterised by diverse shares of species representing different body size classes. This is a consequence of changes in the fractal structure, that is the representativeness of habitats (fractals) with different fractal dimensions. Hence, in young lakes (see the Dendrogram), where the littoral zone is very modestly overgrown with plants, we noted very few smallest beetles associated with peats (tyrphophiles), e.g., Anacaena, Helophorus or Hydroporus, which agrees with the results reported by Bloechl et al. (2010). However, these beetles were very numerous in mature lakes, where a Sphagnum mat dominates. Thus, we determined negative correlations between the body size and the percentage of an area overgrown with Sphagnum moss in the littoral zone of lakes, as well as between the body size (and mean weight) and the abundance of beetles in lakes. Also, we observed that as the fractal dimension of lakes increased, the mean body size (and mean weight) decreased. As spaces between the habitat’s components, i.e., plants, decrease in size and become more and more complex, smaller organisms are evidently preferred, even though larger ones may occur in such habitats as well, because they can push away or move plant stems, small shoots or leaves, an observation also made by McAbendroth et al. (2005). This could be the reason why large beetles were found in the pockets and ponds within the Sphagnum mat, and in the ecotone between the Sphagnum mat and the open water surface.

We support the view of other researchers, e.g., McAbendroth et al. (2005), Bloechl et al. (2010), Linzmeier & Ribeiro-Costa (2011), and believe that large species gain benefits during the early stages of succession, while smaller and specialised ones are at an advantage in the final stages. We have shown (see Fig. 11) that the abundance of beetles is correlated with the age of lakes (the seral stage) (contrary to Bloechl et al. (2010), and with the fractal dimension (McAbendroth et al., 2005). These relationships are connected with the development of plants, which offer a shelter from predators, places for laying eggs and food supply. Like Bloechl et al. (2010) and Cremona, Planas & Lucotte (2008), we demonstrated a relationship between polyphages (Haliplidae) and submerged, more compact plants, while a negative correlation was proven with nympheids, which restrain the growth of food (algae) eaten by polyphages. In turn, the presence of saprophages is correlated with detritus and drifters (which additionally contribute to a growth in the fractal dimension by creating microhabitats that can be used by smaller organisms), which has been confirmed by Verberk et al. (2001) and Cremona, Planas & Lucotte (2008). Moreover, we found out that the abundant presence of predatory beetles in all lakes, and especially among loose vegetation, is favoured by the availability of potential prey, a conclusion also supported by the findings made by Bloechl et al. (2010), whereas a more compact structure is clearly a barrier to the occurrence of the largest organisms. Another factor that limits the number of beetles is the depth. This suggestion is confirmed by the research completed by Brittain (1978), Pakulnicka, Górski & Bielecki (2015b), Pakulnicka et al. (2013) and Pakulnicka et al. (2015a). Among the lakes which we analysed, the ones that contained dark brown water with acidic reaction (Pakulnicka & Zawal, 2018) there were no fish or amphibians, which are usually considered as one of the major factors limiting the presence of water beetles (Eriksson et al., 1980; Solé & Rödder, 2010).