Functional responses of zooplankton communities to depth, trophic status, and ion content in mine pit lakes

Functional traits can be used to identify various ecosystem processes that are influenced by natural and anthropogenic factors. We tested the hypothesis that the species composition and functional feeding traits of zooplankton are influenced by environmental gradients, such as depth, trophic status, and ion content. This hypothesis was tested in three lignite strip mine lakes in Central Poland. Two shallow, polymictic lakes were reclaimed (1994 and 2008), and the third deep, stratified lake has been reclaimed during the flooding process, since 2011. Samples were collected monthly in 2016 and 2017. The functional responses of zooplankton communities to environmental variables were tested through redundancy analyses. Both eutrophy and oligotrophy contributed to the taxonomic heterogeneity and functional homogeneity of zooplankton. In eutrophic condition small microphagous rotifers and small-sized crustaceans dominated, while oligotrophy and high ion concentrations favored dominance of large microphagous and stationary/suspended feeders. In turn, the depth gradient of a newly emerging, stratified lake resulted in taxonomic and functional heterogeneity of zooplankton. Our results complement existing knowledge about zooplankton in post-industrial areas, and also provide new information on the mechanisms of formation and functioning of mine pit lake ecosystems.


Introduction
Functional traits are used to explain variation in the structure and function of zooplankton communities in aquatic ecosystems. The application and interpretation of functional traits requires a comprehensive ecosystem approach, where the diversity and distribution of traits in communities of co-existing species are associated with ecosystem structure and environmental processes (Barnett & Beisner, 2007;Vogt et al., 2013;Pomerleau et al., 2015;Obertegger & Flaim, 2018;Sodré & Bozelli, 2019;Simões et al., 2020). The functional diversity of zooplankton is usually analyzed by grouping species with similar morphological, physiological and ecological traits. Unlike conventional analyses of taxonomic diversity, this approach does not focus on individual species (Naeem & Wright, 2003;Boersma et al., 2016). Therefore, communities with similar functionalities are not necessarily composed of the same species, but they are characterized by the same pattern of functional traits, such as body size, feeding strategy and habitat utilisation (Barnett & Beisner, 2007;Oh et al., 2017;Sodré & Bozelli, 2019;Simões et al., 2020).
According to community assembly theory, the composition, structure and functionality of zooplankton communities are determined by environmental filters and competition (Booth & Swanton, 2002). Environmental filters are abiotic factors (e.g. temperature, depth, water trophy, salinity) that eliminate nonadapted species and, by forming groups with common ecological tolerance, lead to the homogenization of traits. In contrast, competition acts as a factor of trait divergency that leads to the ecological differentiation or homogenization of co-existing zooplankton species. In the competition for common resources, the more efficient competitor wins, and the number of species decreases. However, there are also species in the community whose functional traits allow them to use a different food resource. Consequently, the presence of different functional groups in an ecosystem compensates and the community becomes functionally and often taxonomically heterogeneous. According to Cornwell et al. (2006), both environmental filters and competition occur simultaneously, but they have different consequences for the ecological functionality of communities. Chesson (2000) highlights that both environmental filters and competition can lead to trait similarity and dissimilarity, and the final result depends on the robustness of ecosystem mechanisms.
The analyses of functional traits support the identification of ecosystem phenomena that are not revealed by taxonomic analyses (Barnett & Beisner, 2007;Parravicini et al., 2014).
Most studies on zooplankton functional traits have been conducted in natural ecosystems using long-term datasets. However, there is a general scarcity of observations concerning the functional diversity of artificial reservoirs with a much shorter history, including mine pit lakes whose number continues to increase in reclaimed mining sites. Pit lakes occupy basins created by strip mining (Nixdorf et al., 2005;Gammons et al., 2009;Sloss, 2013). The flooding of excavated pits is often used to reclaim former mines to remedy environmental disturbance resulting from strip mining (Larondelle & Haase, 2012;Blanchette & Lund, 2016;Stachowski et al., 2018), and it delivers social and environmental benefits (McCullough & Lund, 2006;Lei et al., 2016;Lechner et al., 2017). Newly constructed aquatic ecosystems quickly integrate into the local environment and landscape, and biological succession processes occur in these ecosystems (Kalin et al., 2001;Gammons et al., 2009;Goździejewska et al. 2019b). The rate and direction of transformations are influenced by the morphometric and geochemical properties of excavated pits and the hydrological regime (Kumar et al., 2009;Larondelle & Haase, 2012;Sloss, 2013). The geologic structure of the area and the groundwater quality are the main factors determining the primary composition and chemical properties of water. Similarly to natural ecosystems, the land-use pattern in the direct catchment, the quality of inflows from the catchment, and the type and intensity of lake use play an important role in the development of mine pit lakes (Axler et al., 1998;Kleeberg et al., 2008). The sensitivity of these relatively young environments can be analyzed based on selected functional traits of the colonizing species to expand our understanding of numerous ecological relationships and to identify potential threats (El-Bassat & Taylor, 2007;Marszelewski et al., 2017). This is an important consideration for preserving water resources created in former industrial areas that are devoid of natural lakes.
This study analyzed zooplankton structure in three mine pit lakes (Janiszew, Przykona and Władysławów) which had been formed as a result of a land reclamation project initiated by the Adamów Coal Mine in former lignite strip mines in Central Poland. The first limnological studies of lakes Janiszew and Przykona were carried out by Marszelewski et al. (2017), whereas Lake Władysławów was previously analyzed by Goździejewska et al. (2019b) and Skrzypczak et al. (2019). These studies revealed differences in the rate and direction of ecological succession in the analyzed lakes, which can be attributed to the lakes' age, morphometric traits, and catchment characteristics. Two environmental characteristics, depth and trophic state, were identified as potential determinants of lake functionality. The depth of the lake determined the water mixing dynamics, thermal and oxygen conditions, while the trophic level determined nutrients concentration. The aim of this study was to investigate the impact of depth and trophic gradients on planktonic communities. We hypothesized that the species composition and functional feeding traits of zooplankton would be different in lakes vary of depth and trophic status. We also noted the importance of geochemical properties, such as ion concentrations in the water. We expected that an increase in water trophy will result in taxonomic and functional homogenization of the zooplankton. In turn, an increase of depth gradient will favor the heterogeneity of these zooplankton characteristics. This study complements existing knowledge and provides new information about the functioning of food webs in mine pit lakes.

Study area
The study was conducted in three mine pit lakes located in the vicinity of the Turek-Adamów lignite strip mine: Janiszew (JAN), Przykona (PRZ) and Władysławów (WLA) (Fig. 1). The lakes vary in age, depth and basin shape based on the mining operations that were present. They also differ in the type of direct catchment and form of use, depending on ownership rights and land-use types prescribed in the local zoning plan (Ró _ zkowski et al., 2010;Marszelewski et al., 2017;Goździejewska et al., 2019b; Table 1). PRZ is the oldest (1994), largest and shallowest of the lakes, and it occupies the former opencast mine in Adamów (OM-A; Fig. 1). PRZ was supplied with water from the Teleszyna River and groundwater removed from the mine. This polymictic lake has a regularly shaped basin, a flat bottom, gently sloping shores and high water mixing dynamics at average depth (4.3 m) (Marszelewski et al., 2017;Goździejewska et al., 2019b). The lake became municipal property in 2004, and it has experienced considerable recreational pressure since that time. Residential buildings, recreational facilities, bathing areas, an angling base, and a slip for launching motor boats have been developed in the lake's direct catchment (Ró _ zkowski et al., 2010). Goździejewska et al. (2019b) reported an increase in water trophy of PRZ.
In 2008, the basin of JAN was filled with water removed from the opencast mine in Koźmin (OM-K; Fig. 1). The morphometric traits and the water mixing regime in JAN are similar to those noted in PRZ, and the lake's water turns over from top to bottom during the year ( Supplementary Fig. S1). JAN has an average depth of 6.5 m. It is surrounded by fallow land and fields directly along its shoreline. The lake and its immediate surroundings are owned by the mine, and fishing and other recreational activity are prohibited (Marszelewski et al., 2017;Goździejewska et al., 2019b).
WLA occupies the former Władysławów strip mine which has been reclaimed by long-term flooding since 2011 (Ró _ zkowski et al., 2010;Stachowski et al., 2018). The lake does not have an outflow, and it is supplied mainly with groundwater, rain water and catchment water. Unlike PRZ and JAN, WLA is a deep lake whose maximum depth was estimated at 25 m during the study (the maximum target depth is 37 m, with an estimated average of 21.5 m). Steep slopes with an inclination angle of up to 30% account for around 80% of the shore line. The remaining part of the shoreline features gentle slopes. The lake is characterized by thermal and oxygen stratification in summer ( Supplementary Fig. S1). The flooding process has not yet been completed, the lake constitutes the property of the mine, and it is not open to the public.
Sampling and analytical procedure Zooplankton were sampled every four weeks between July and October 2016, and between May and September 2017. A single sampling sites was located in the deepest part of each lake. In shallow, polymictic lakes PRZ and JAN, samples were collected between the surface and the bottom at a depth of 0, 2, 4 and 6 m (PRZ) and 0, 3, 5 and 8 m (JAN), and combined into single samples for each lake. This layer was defined as the epilimnion due to the polymictic and homothermic character of the entire water column. In the deep, stratified lake WLA, zooplankton samples were collected at three depths: 0-5 m, 6-10 m and 15-20 m. These layers were formed within the reach of the epilimnion, metalimnion and hypolimnion, respectively, during the summer stagnation period. All samples from the three lakes were collected with a 5 l Patalas trap. Samples of 20 l from each layer were passed through a plankton net with 30 lm mesh size, preserved in Lugol's solution, and fixed in 2% formalin solution. During the experiment, a total of 45 zooplankton samples were collected, including 9 PRZ, 9 JAN, and 27 WLA.
The taxonomic variability of zooplankton was evaluated with the dominance index (D) (Kasprzak and Niedbała, 1981), diversity index (Shannon's index H') and evenness index (Pielou's index J'). The taxonomic diversity and similarity of zooplankton was calculated based on mean species abundance (ind l -1 ) colonizing the epilimnion layer of the lakes and was compared with the use of MVSP 3.22 software (Kovach, 2015).
To assess the functional variability of zooplankton, we used traits based on body size and feeding strategies. On this basis, rotifer and crustacean species were classified into four groups: small microphagous (SMF), large microphagous (LMF), raptorials (RAP), and stationary/suspended (STA/SUS) feeders, with the use of the methods described by ), Bertani et al. (2012, Vogt et al. (2013), and Moreira et al. (2016).
The ratio of trophic groups (GR') in the zooplankton was determined to characterize ecosystem functionality of the lakes (Smith et al., 2009;. GR' was calculated for the mean biomass (mg l -1 ) of each functional groups in the entire water column, from the surface to the bottom of each lake, with the following formula: The values of GR' ranged from -1 to 1. Values below 0 indicate the dominance of microphagous feeders, and values higher than 0 denote the dominance of raptorial feeders.
To determine the response of zooplankton trophic groups to environmental variables, we used a grazing rate coefficient. The grazing rates for each feeder type were calculated as a function of the biomass of zooplankton feeding on phytoplankton (Lampert, 1988): where: c is the grazing rate (% day -1 ) and W is the biomass of zooplankton trophic groups (lg l -1 ).

Statistical analyses
The overall differences in the physical and chemical parameters of water and zooplankton parameters across the studied lakes were determined by nonparametric analyses (Statistica 13.0 for Windows, Statsoft, Tulsa). Specifically, we used non-parametric ANOVA and Kruskal-Wallis test to determine statistically significant differences in the physical and chemical parameters of water and the taxonomic and functional structure of zooplankton between lakes (H, P B 0.05). Correlations between zooplankton species and environmental variables were calculated with the Spearman's rank-order correlation (P B 0.05).
The influence of the trophic, depth, and ion gradients on the spatiotemporal variation in the distribution of zooplankton taxa and functional groups was determined with the use of multivariate ordination techniques in CANOCO 4.56 software (ter Braak & Š milauer, 2002). Two analyses that differed in the range of data were performed. The first analysis evaluated the influence of the trophic gradient, productivity and ions on the zooplankton in the polymictic zone of lakes (epilimnion). The abundance of zooplankton species were log (x ? 1) transformed, and taxa with minimum 2% relative abundance in a given lake were used in further analysis. A preliminary detrended correspondence analysis (DCA with detrending by segments) was performed to estimate gradient length in standard deviation (SD) units. Gradient length was shorter than 3 (2.475 SD), and subsequent numerical analyses involved techniques that were based on the underlying linear speciesresponse model (Jongman et al., 1995). Variance inflation factors (VIF) were applied to identify correlated variables, and the results were used to eliminate selected variables from further analyses. A total of eighteen environmental variables were included in preliminary analyses, and six of those were eliminated. SDT, PO 4 3-, NO 3 -, EC, Ca 2? and color were not included in the analyzed dataset because these parameters were strongly correlated with the selected variables.
The second analysis examined the effects of depth (temperature, oxygen concentration and density of water) on the zooplankton in the stratified Lake WLA, and species with a minimum relative abundance of 1% were used in the multivariate ordination analysis. A total of seven environmental variables were included in analyses.
The significance of canonical axes and environmental variables to explain the variance of the community was tested using Monte Carlo simulations with 499 permutations. The forward selection option in CANOCO was then used to identify environmental variables that significantly influenced the distribution of zooplankton.
A comparison of zooplankton structure based on Sorensen's coefficient revealed qualitative differences in each lake across seasons (Fig. 2). In a dendrogram based on the values of the Bray-Curtis dissimilarity index, the zooplankton groups in the studied lakes  were divided into two clusters, revealing the differences between the seasons of 2016 and 2017 in PRZ. The structure of zooplankton in PRZ was characterized by 57% similarity to JAN in 2017, and by 32% similarity to WLA in 2017. The zooplankton community in JAN was most stable and similar (63%) during the entire study (Fig. 2). The zooplankton trophic groups were most abundant in PRZ lake, in order: SMF (0.03-0.32 mg l -1 ), LMF (1.81-31.28 mg l -1 ) and raptorial (0.29-9.43 mg l -1 ). The biomass of STA/SUS was highest in JAN (0.12-3.48 mg l -1 ; Fig. 3). The differences in the mean biomass of SMF, LMF and RAP functional groups in the lakes were statistically significant (P B 0.05; Fig. 3). Seasonal variations in GR' values ranged from -0.392 to -0.929 in PRZ and from -0.222 to -0.732 in JAN (Fig. 4). Large microphagous feeders in lakes PRZ and LMF and STA/SUS groups in JAN had the highest proportions of total zooplankton biomass during the entire study ( Supplementary Fig. S2). The broadest range of GR' values was noted in WLA between the epilimnion in August (-0.923) and the hypolimnion in October (0.985) 2016 (Fig. 4), which resulted from temporal  Fig. S2).
The mean zooplankton grazing rates ranged from 18.9% day -1 (JAN) to 78.8% day -1 (PRZ), but only remained stable during both years of the study in JAN. The grazing rate in Lake WLA decreased from 57% in 2016 to 7% day -1 in 2017, and it increased from 18% to 128% day -1 in PRZ in the corresponding years. SMF group was characterized by the lowest grazing rate which ranged from 1% (WLA) to 6% day -1 (PRZ) on average (Fig. 5). The abundance and total grazing efficiency of the LMF group peaked in PRZ in 2017 (1453 ind. 1 -1 and 419% day -1 on average; Fig. 5) when the lake was massively colonized by Keratella quadrata and B. longirostris. In the remaining lakes, the grazing rate of the LMF group ranged from 37% to 75% day -1 on average, and a minimum value of around 6% day -1 was noted in WLA in 2017. The representative raptorials, A. ovalis and Polyarthra spp., dominated in terms of zooplankton biomass and abundance, mostly in WLA and PRZ (Fig. 5, Table 2).
In shallow polymictic lakes JAN and PRZ water temperature changed in a uniform manner between the surface and the bottom under the influence of seasonal changes in air temperature. In the stratified WLA lake, thermal stratification and oxygen distribution also changed with depth ( Supplementary Fig. S1). During the summer stagnation period, a thermocline was observed between the depth of 7 m and 9 m, and temperature and oxygen saturation gradients between the epilimnion and the hypolimnion were determined at 13.9°C and 2.3 mg O 2 l -1 , respectively. Metalimnetic oxygen maxima with a positive heterograde oxygen curve (11.3 mg O 2 l -1 , 107% saturation) were observed in July 2017 ( Supplementary Fig. S1).

Primary predictors of zooplankton variability
In a redundancy analysis of the polymictic layer, five of the 12 environmental variables (TN, TH, SUS, Chla and BOD) explained 67.4% (P \ 0.05) of total variance in zooplankton data. Forward selection of explanatory variables demonstrated that the best predictors zooplankton assemblage were TN (16%, P = 0.002), followed by TH (8%, P = 0.004), SUS (8%, P = 0.010), Chl-a (8%, P = 0.026) and BOD (5%, P = 0.046; Supplementary Table S1). The Monte Carlo permutation test revealed significant correlations between species and environmental factors along the first two RDA axes (F = 2.415, P \ 0.002). The correlation coefficients and the explained variation for axis 1 and axis 2 were determined at 0.97 and 0.912, and 21.5% and 15.5%, respectively (Supplementary Table S2). On the first axis, the highest correlation between environmental variables and sample location was determined for TN (r = 0.79). Axis 2 was correlated with Chl-a (r = 0.36).
All predictors of zooplankton variability, indicating water trophy (TN, TP), productivity (Chl-a, suspension) and organic matter content (BOD, TOC) were correlated in PRZ (Fig. 6). The mean values of TN and TP were significant highest in PRZ (P B 0.05; Table 3), and they were correlated with BOD (r = 0.64 and r = 0.51, respectively; P \ 0.05). The concentration of Chl-a, turbidity and water color were also significantly higher in PRZ than in JAN and WLA (Table 3). The RDA biplot for species and environmental variables revealed that Bosmina longirostris (LMF) and Asplanchna priodonta Gosse, 1850 (RAP) were positively correlated with TN; Daphnia cucullata (LMF) and Keratella quadrata (LMF) were positively correlated with Chl-a; and Brachionus angularis Gosse, 1851 (SMF), Filinia longiseta (RAP) were correlated with BOD (Fig. 7). Total hardness (8%, P = 0.004) was the second predictor of explanatory variables in zooplankton communities in RDA. Total hardness is an indicator of ion content (in addition to EC, chlorides and calcium), and it was correlated with sampling sites in WLA and JAN (Fig. 6). The differences in the mean values of the above parameters in the lakes were statistically significant (P B 0.05; Table 3). The RDA biplot for species and environmental variables revealed that most taxa, in particular LMF (Daphnia cucullata, Diaphanosoma brachyurum, Chydorus sphaericus (Müller, 1785), Eubosmina longicornis kessleri (Uljanin, 1874), and Keratella quadrata) and SMF (e.g. Keratella cochlearis, Pompholyx sulcata Hudson, 1885 and nauplii), were negatively correlated with TH. In turn, raptorial species Gastropus stylifer Imhof, 1891, Polyarthra longiremis and Ascomorpha ovalis were positively correlated with TH and Cl (Fig. 7).
An RDA relating to the depth gradient of WLA revealed that 3 of the 7 environmental variables (DO, %DO and EC) explained 59.5% (P \ 0.05) of the variation in the distribution of zooplankton taxa and   Table S3). The speciesenvironment correlations were significant along first two RDA axes (Monte Carlo test, F = 3.152, P \ 0.002). Axis 1 accounted for 37.2% and axis 2 for 9.5% of the explained variation (Supplementary Table S4). The highest correlation between environmental variables and sampling location was observed for %DO (r = -0.46) on the first axis. Axis 2 was correlated with T (r = -0.50). The RDA biplot for species and environmental variables indicated that taxa such as Polyarthra spp. (RAP), Diaphanosoma brachyurum (LMF) and Thermocyclops crassus (Fischer, 1853) (RAP) were positively correlated with epilimnion conditions: temperature, DO and Chl-a; whereas Ascomorpha ovalis (RAP), Gastropus stylifer (RAP) and Keratella tecta (Gosse, 1851) (SMF) were positively correlated with %DO and EC, and were characteristic of the metalimnion. In turn, copepods Cyclops vicinus Uljanin, 1875 (RAP) and nauplii (SMF) were negatively correlated with temperature, oxygen concentration and productivity, and they were identified in the hypolimnion (Fig. 8).

Discussion
Our study revealed that the water trophy, depth and ion concentrations influenced the diversity of species and functional structure of zooplankton. The chemical and zooplankton indices confirmed the differences in trophy between the polymictic lakes JAN and PRZ. The high trophic status of PRZ can be attributed to eutrophication that was accelerated by anthropogenic factors. In comparison with the results reported by Marszelewski et al. (2017), the average values of TSI SD, TSI Chl-a and TSI ROT increased by 36%, 25% and 5%, respectively, over a period of only 5 years. As a result, zooplankton abundance increased more than 20-fold on average, with a high proportion of trophiclevel indicator species, including rotifers Keratella quadrata, Brachionus angularis and Trichocerca similis (Ejsmont-Karabin, 2012) and crustaceans Daphnia cucullata, Bosmina longirostris, Chydorus sphaericus, Diaphanosoma brachyurum and Thermocyclops sp. (Ejsmont-Karabin & Karabin, 2013). It should be noted that the some functional traits of above-mentioned bioindicator species were previously correlated with water trophy. For example Manca (2011) andOh et al. (2017) positively correlated the morphological traits (trophi type) of the genera Keratella and Brachionus and Polyarthra and Trichocerca with high trophic-level parameters and productivity (Chl-a). In turn, the smaller body size of Cladocera and Copepoda indicator species (Bosmina, Eubosmina, Chydorus, Thermocyclops) were previously correlated with higher lake trophy (Barnett & Beisner, 2007;Sodré & Bozelli, 2019). Therefore, high water trophy in PRZ led to the prevalence of species characterized by small size and different trophic type, and it also supported colonization by many eutrophic species, confirmed by RDA analysis. The above contributed to high taxonomic diversity but low lake functionality, indicated by low gradient related to the ratio of zooplankton trophic groups (GR'). A dominance of LMF was observed throughout the entire study, and it is indicative of ecosystem disturbance and functional depletion (homogenization) of PRZ (Hampton, 2005;Hogsten et al., 2009). According to Petchey (2009), prolonged and unidirectional use of food resources can compromise an ecosystem's resistance and contribute to its degradation.
PRZ is subjected to high levels of anthropogenic pressure (water recreation, fishing, possible sewage inflows) and is characterized by small volume; therefore, eutrophication proceeded rapidly. The analyzed lake's productivity increased 11-fold during the twoyear study (2016 and 2017), which prompted a rapid response from both zooplankton taxa (Bray-Curtis dissimilarity index) and feeding guilds. Obertegger and Manca (2011) demonstrated that fluctuations in water trophy significantly influence competition among rotifers. Microphagous rotifers dominate in lakes characterized by higher trophy, whereas raptorials are most prevalent when trophic characteristic are low. These variations can be attributed to the availability of food sources. According to Galir Balkic et al. (2018), small filtering zooplankton such as the smallest rotifers generally play a minor and relatively stable role in algal grazing, even during rapid fluctuations in environmental conditions, probably due to a high share of detritus in their diet. However, during dense phytoplankton blooms (Chl-a = 59.9 lg l -1 ) in PRZ in 2017, the grazing efficiency of SMF increased by 50% on average, whereas their abundance remained similar relative to 2016. The above observation was attributed to the fact that the small Keratella cochlearis had been replaced by Copepoda nauplii which grazed more efficiently on algal blooms. According to Kim et al. (2000), copepod nauplii are more effective filter feeders than many Rotifera species, including K. cochlearis; and this difference is particularly pronounced when both groups graze on algae. In the analyzed case, this taxonomic replacement occurred during dense phytoplankton blooms (Chl-a [ 20 lg l -1 ; Kawecka & Eloranta, 1994) with a dominance ([ 50% quantitative proportion) of Microcystis aeruginosa (Kützing) Kützing 1846, Aphanizomenon sp. and Pseudanabaena mucicola (Naumann & Hub.-Pest.) Schwabe 1964 cyanobacteria (Table 1). Zhao et al. (2008) observed that high copepod functionality is associated with selective feeding on alternative food sources, such as diatoms, detritus and green algae, during cyanobacterial blooms, whereas Rotifera and Cladocera (non-selective filter-feeders with identical feeding preferences) can be limited or eliminated.
Cyanobacteria forming large mucilaginous colonies or large filaments with aerotopes (such as the above species) are generally not a suitable food source for rotifers (Colina et al., 2015), but large populations of LMF (Keratella quadrata, Brachionus angularis) and RAP (Polyarthra spp., Trichocerca similis) were identified in PRZ. According to Soares et al. (2010) and Galir Balkić et al. (2018), this could be the result of adaptation to cyanobacteria in their diet (Brachionus), grazing on alternative food sources such as bacteria and detritus (K. quadrata), or switching to animal food resource, such as protozoa (T. similis). Tõnno et al. (2016) showed resistance and/or feeding selectivity of small Cladocera species, including Ch. sphaericus towards cyanobacterial colonies, which explains/allows their presence during bloom. These effects were observed in PRZ in 2017 when algivorous species dominated other trophic groups in terms of grazing efficiency, which points to strong interactions with the food web.
The oligotrophic conditions in JAN occupied the opposite end of the gradient relative to the eutrophic PRZ, which could be attributed to the absence of anthropogenic influences in JAN which is not open for public use. This comparison revealed that the content of ions in the water (TH, Ca 2? , EC, Cl -), resulting from the geochemical properties of excavated pits, exert a significant direct and indirect influence on ecosystem, including the trophic state of lakes.
The concentration of calcium cations plays a major role in food webs by promoting optimal water hardness (1.5-4 mval l -1 and [ 2 mg Ca 2? l -1 ) (Wurts & Durborow, 1992;Hessen et al., 2000;Korosi et al., 2012;Boyd et al., 2016). Water hardness in JAN and the remaining lakes was within the upper range of the cited values, whereas the concentration of calcium cations exceeded the critical minimum threshold 20-to 40-fold. Therefore, these factors could not directly disrupt life processes (osmoregulation, growth, reproduction) of planktonic crustaceans, and similar observations were made by Hessen et al. (1995), Waervågen et al. (2002), Ashforth &Yan (2008), andJeziorski et al. (2012). However, calcium concentration in the oligomesotrophic JAN was nearly two-fold higher, and chlorophyll a concentration was approximately six-fold lower than in the eutrophic PRZ lake. The above indicates that Ca 2? could be the main regulator of phytoplankton productivity and species composition which influenced the availability of phosphorus, the principal nutrient (Korosi et al., 2012;Boyd et al., 2016). In JAN, the phosphorus balance could also have been affected by stoneworts (Characeae) which densely colonized the lake bottom in both shallow and deep zones. According to Kufel & Ozimek (1994), these macroalgae are characteristic of hard mesotrophic waters with low phosphorus content. Phosphorus is also more rapidly accumulated by Chara spp. than by other aquatic plants (macrophytes and phytoplankton), thus eliminating or inhibiting the development of these trophic groups through competition for food or allelopathy (Kufel & Kufel, 2002;Lürling et al., 2006). Low phosphorus levels in JAN contributed to the abundance of planktonic algae of the genus Dinobryon (Chrysophyceae). The diet of representatives of this mixotrophic genus may also include bacteria, therefore they often dominate in oligotrophic water bodies, forming large and numerous colonies (Sommer et al., 2003;Kamjunke et al., 2007;Padisák et al., 2009). Regardless of its phenology, the genus Dinobryon (which occurs in cool seasons of the year, including spring; Heinze et al., 2013) was observed during the entire study in JAN, and its abundance peaked in August 2017. The above also indicates that low phosphorus concentration was a permanent, rather than a seasonal characteristic of JAN.
Lower water trophy in JAN than in PRZ led to a four-fold decrease in total zooplankton abundance. Ecosystems with low trophic characteristic (oligotrophic lakes, flowing waters) are generally characterized by abundance of Rotifera species which have a short life cycle and low nutritional requirements in comparison with crustaceans (Baranyi et al., 2002;Goździejewska et al., 2010;Goździejewska et al., 2016). However, rotifers were least abundant in JAN. The above could have directly resulted from the absence of flagellates such as Dinobryon in rotifer diets. In turn, a high concentration of dissolved calcium in JAN and the availability of large Dinobryon colonies probably contributed to the development and survival of large planktonic crustaceans which were abundantly represented by Ceriodaphnia quadrangula and Eudiaptomus graciloides. These species are highly effective filtering feeders and grazers, and they effectively competed for food resources with SMF, and, possibly, directly eliminated SMF from the environment by mechanical interference them in filtering chambers (Sommer et al., 2003;Balvert et al., 2009).
The abundance of many Rotifera species in JAN could also be limited by high electrolytic conductivity (EC). The EC of natural fresh water bodies is usually positively correlated with their trophic status (concentrations of nutrients and chlorophyll a) (Bērziņ š & Pejler, 1989;Diaz et al., 2007;Brito et al., 2011). In man-made reservoirs, including mine pit lakes, EC is often correlated with the presence of dissociated carbonate, chloride and sulfate ions or heavy metals (Molenda, 2014;Pociecha & Bielańska-Grajner, 2015). The RDA indicated strong correlations between EC and all parameters responsible for EC (Ca 2? , TH, Cl -). Electrolytic conductivity was bound by a negative correlation with TP and Chl-a, which indicates that EC strongly influenced the functioning of plankton communities in JAN. The quantitative composition of species in JAN was dominated by K. cochlearis (SMF) and Polyarthra spp. (RAP)-rotifers that are capable of colonizing multiple environments, including mine pit lakes with high EC (EC [ 600 ls l -1 ; Bielańska-Grajner & Gładysz, 2010; Pociecha & Bielańska-Grajner, 2015). Therefore, ion concentrations in JAN appear to play an important role, by favouring specific functional traits of zooplankton, such as large body size in Cladocera (Ceriodaphnia) and Copepoda (Eudiaptomus). Ion concentrations is also responsible for competition in the food web (elimination of rotifers) and the homogenization of phytoplankton resources in a reservoir with permanently low TP levels. As a result, the trophic structure of zooplankton was based on the codominance of LMF and STA/SUS throughout the entire study. The above was responsible for the low trophic gradient of zooplankton feeding guild taxa, and it testifies to the low functionality of the JAN ecosystem (Moreira et al., 2016).
The observed changes in environmental conditions, which were reflected by planktonic structures in WLA, resulted from early dynamic ecological succession processes as well as depth gradient. According to Gammons et al. (2009) and Audet et al. (2013), the abundance and diversity of zooplankton and phytoplankton are a function of the time elapsed since a lake's formation. WLA was characterized by a predominance of Pyrrophyta (genera Peridinium and Ceratium), which are often regarded as pioneer species in mine pit lakes younger than 10 years (Gammons et al., 2009). These freshwater dinoflagellates have a vertical distribution pattern, and they accumulate at depths that are optimal for photosynthesis (Regel et al., 2004;Darki & Krakhmalnyi, 2019). As a result, oxygen levels increase in the metalimnion (at a depth of approx. 8-9 m in this study) as an effect of intensive photosynthesis and favorable light conditions. The variations in phytoplankton structure and biomass in stratified lakes induced by a vertical gradient of light, temperature and density (EC) significantly differentiate habitats and the structure of zooplankton communities (Obertegger & Flaim, 2018). Oxygen saturation is yet another important factor that changes with depth (Karpowicz et al., 2020). In WLA, most of the dominant zooplankton species colonized the epi-metalimnion. In 2016, Diaphanosoma brachyurum (Liévin, 1848), a large filter feeder, was the dominant species mostly in the polymictic layer (epilimnion); whereas in 2017, the availability of Peridinium and Ceratium as food resources led to colonization by efficient raptorial consumers Ascomorpha ovalis, Polyarthra spp. and Gastropus stylifer, mostly in metalimnetic conditions (Stelzer, 1998). The above changes in zooplankton community can be also attributed to top-down biotic factors associated with fish pressure on planktonic crustaceans. In the future, there are plans to stock WLA with fish for recreational purposes (Stachowski et al., 2018). However, the reservoir has been already colonized by fish as a result of natural processes and uncontrolled stocking (Skrzypczak et al., 2019). Thermophilic fish most readily graze in the epimetalimnion where they prey on large and slow organisms (Cladocera). Therefore, migrating crustaceans often find shelter in deeper parts of the pelagic zone (hypolimnion) during the day (Gauthier et al., 2014, Karpowicz et al., 2020. However, these layers may be less accessible due to decreased availability of nutrients (lower primary production) and progressing oxygen depletion. Oxygen stress can affect both taxonomic and functional structure of zooplankton as small species replace large ones. These processes are intensified with hypoxic conditions (DO \ 2 mg l -1 ) (Karpowicz et al., 2020). In the hypolimnion of WLA, oxygen concentrations at a depth of 25 m and temperature of around 6.5°C ranged from 5.4 to 7.7 mg l -1 . These conditions did not compromise the activity of grazing raptorials represented by numerous populations of large-sized copepods Cyclops vicinus and C. scutifer G.O. Sars, 1863.
This study confirmed that the depth gradient significantly contributed to variations in habitats and food niches in the young ecosystem of WLA. However, these observations were not confirmed by the values of biodiversity indices, and the species diversity of zooplankton was lowest in WLA relative to JAN and PRZ. In turn, the highest value of GR' in WLA could be attributed to differences in trophic characteristic in two successive years of the study and between depth zones. The predominance of raptorials in WLA is indicative to low water trophy levels and the scarcity of phytoplankton, which means low density and/or dominance of large algae species, unavailable to filter feeders (Obertegger & Manca, 2011).

Conclusions
The taxonomic diversity and functionality of zooplankton in the mine pit lakes were determined by the selective influence of environmental conditions. High primary production contributed to taxonomic heterogeneity, but it led to the functional homogeneity of plankton communities with a dominance of cyanobacteria and functionally specialized rotifers and crustaceans (microphagous rotifers and small-sized crustaceans). Low primary production, oligotrophy and high ion concentrations contributed to high species diversity of zooplankton and the functional simplification of the plankton community with a predominance of Dinobryon colonies and large-sized Cladocera (LMF) and Calanoida (STA/SUS). In this case, geochemical properties were a strong and stable influencing factor which led to the development of permanent competitive and allelopathic (Characeae) relationships in plankton communities and consequently their homogenization. The ecosystem of the deep stratified lake was characterized by greater spatial heterogeneity of habitats and food niches due to a vertical gradient of temperature, light, oxygen saturation and water density. Large thermophilic filter feeders and raptorial rotifers grazed in the epimetalimnion. Large-sized Copepoda that are resistant to oxygen stress found refuge against predatory fish in the hypolimnion.
Post mine lakes are a unique phenomenon of the formation of ''new'' ecosystems. It is possible to observe ecosystem changes related to the early stages of biological succession in these lakes. Our analyses of zooplankton functional traits showed the stable effects of oligotrophy and the dynamics of formation of the new ecosystem. They are also an excellent tool for monitoring threats in aquatic ecosystems resulting from anthropogenic pressure.