Top-down and bottom-up effects on zooplankton size distribution in a deep stratified lake

Trophic interactions in the pelagic area of lakes and the opposing effects of fish feeding (top-down) and phytoplankton biomass (bottom-up) on zooplankton communities are central topics in limnology. We hypothesized that zooplankton size distributions should be a more sensitive approach to disentangle top-down and bottom-up effects than the commonly measured zooplankton biomass. We examined zooplankton size distributions from 148 samples collected during summer months in the upper and lower pelagic layers of a deep mesotrophic lake among 13 years of sampling. Top-down effects, namely fish size and biomass, and bottom-up effects, including water temperature and total phosphorus and chlorophyll a concentrations, were considered. To add robustness to our analyses, we expressed the zooplankton size distributions as size spectra based on log-binning, as continuous size spectra and by the size diversity, a measure that has been developed to mimic taxonomic diversity indices. Among numerous regressions tested, significant top-down or bottom-up effects could rarely been detected. Our results indicate that the overall zooplankton size distribution was not significantly affected by fish predation and lake productivity measured as total phosphorus or chlorophyll a concentration. However, we found negative correlations between fish biomass and the preferred zooplankton prey, including Bosmina longirostris, Daphnia cucullata and nauplii in the upper depth layer. However, due to their small body size, low biomass and therefore relative small contribution to the zooplankton size distribution, predation on preferred zooplankton species did not translate into a statistically significant modification of the entire size spectrum. Consequently, the size spectrum seems to be relatively robust against predation effects, but might reflect the lake-wide energy availability and transfer efficiency in the food web.


Introduction
Trophic interactions in the pelagic area of lakes are one of the cornerstones of classical and modern limnology. Many scientific studies centred around the opposing effects of fish feeding (top-down) and phytoplankton biomass (bottom-up) on zooplankton communities (Brooks and Dodson 1965;Carpenter et al. 1985;McQueen et al. 1986), while recent studies have also recognised the importance of the conditions at which top-down and bottom-up effects dominate (Carpenter et al. 2016). Zooplankton are an integral part of aquatic ecosystems, playing a crucial role in linking primary producers and higher trophic levels such as fish. In turn, zooplankton communities are sensitive to changes of their resources and their predators, and therefore mirror a balance of food web processes by their body size distribution and taxonomic composition (Mills and Schiavone 1982;Hansson et al. 2007).
However, detecting bottom-up and top-down effects on zooplankton becomes complicated because the zooplankton community is not a homogenous group, but is composed of several taxonomic and functional groups (e.g., cladocerans vs. copepods; filter-feeders vs. particulate feeders; herbivorous vs. omnivorous or carnivorous species). Accordingly, both the selective feeding of zooplankton on phytoplankton and the selective feeding of zooplanktivorous fish on zooplankton may modify the relative densities of zooplankton taxa and subsequently the strength of competitive interactions within the zooplankton community. It is therefore not surprising that studies on pelagic food chains have not identified a single 'master' mechanism, which drives the zooplankton community composition and biomass in the pelagic area of all lakes.
In contrast, it might be possible that the zooplankton size distribution is a more sensitive reflection of the strength of bottom-up or top-down effects in pelagic layers. The body size of individuals is a simple metric to collect, while it is a key characteristic correlated with physiological and ecological properties. The size of an organism scales with traits such as metabolic rate, growth (Marquet et al. 2005;Woodward et al. 2005), life span (Peters 1983;Calder 1984) and trophic position. In aquatic ecosystems, smaller organisms are typically more abundant than larger ones (Kerr 1974;Sprules et al. 1991;Thiebaux and Dickie 1993). This negative relationship between size and abundance emerges from organism traits and ecological interactions (Sprules and Barth 2016) and is often referred to as a community size spectrum (SS). The SS describes the relationship between the logarithm of size and the logarithm of abundance, which is usually linear and can therefore be characterized by slopes and intercepts of linear regressions (Sheldon et al. 1972;Kerr and Dickie 2001). The intercept reflects the overall abundance of organism in the community and the slope mirrors the relative abundances of small and large organism (Ahrens and Peters 1991). Therefore, SS are a highly effective approach to summarize and compare the size structure of aquatic communities (Cottingham 1999;Cózar et al. 2003) and might be useful to evaluate resource availability for and selective predation on zooplankton communities.
A combination of abiotic and biotic factors influences zooplankton community size structure. Growth of zooplankton populations reflects factors regulating the metabolism, including temperature (Tordesillas et al. 2016) and food availability (Liu et al. 2015). Rising temperatures favour small organisms (Atkinson et al. 2003;Havens et al. 2015;Andersen et al. 2016) and therefore result in a shift towards smaller body size (Gardner et al. 2011), reflected in more negative SS slopes and lower size diversity, whereas greater resource availability leads to an increase in relative abundance of large organisms and therefore flatter SS slopes (Gaedke et al. 2004;Sprules and Barth 2016). Specifically, increases in total phosphorus (TP) lead to increases of overall zooplankton community abundances via higher food availability (Hanson and Peters 1984;Pinto-Coelho et al. 2005). However, top-down predation effects from fish (Jeppesen et al. 1997;Zhang et al. 2013) strongly affect biomass and body size distribution of zooplankton communities (Brooks and Dodson 1965;Badosa et al. 2007), leading to a reduced size diversity (Brucet et al. 2010). It is challenging to elucidate, which of these factors exhibit the greatest effect on trophic levels and therefore population dynamics and size structure, as the impact of each occurs simultaneously (Shurin et al. 2012;Liu et al. 2020) and interactive effects of these factors can be very complex.
In the present study, we used zooplankton SS to evaluate the strength of bottom-up and top-down directed effects on zooplankton size distributions in a deep mesotrophic lake. We examined the zooplankton size distribution from 148 samples collected exclusively during summer months among 13 years in the upper and lower pelagic layers of a deep mesotrophic lake. To add robustness to our analyses, we expressed the size distributions as SS based on log-binning, as continuous SS and by the size diversity, a measure that has been developed to mimic taxonomic diversity indices (Brucet et al. 2006(Brucet et al. , 2010(Brucet et al. , 2017Quintana et al. 2008Quintana et al. , 2016. We expected that both annually varying predation by the dominant pelagic planktivorous fish and resource availability via phytoplankton biomass would affect the year-to-year variation of the zooplankton size distributions. Therefore, we explored whether TP, chlorophyll a (chl a) concentrations and water temperature modified the size distribution of zooplankton. In turn, we tested whether the size distribution of fish and their biomasses as obtained by two different fishing methods had an effect on zooplankton size distributions. The lake is dominated by coregonid fishes, which show strong variations in year-class strength and annual growth rates resulting in strongly differing abundances between subsequent years (Helminen and Sarvala 1994;Auvinen et al. 2000;Viljanen et al. 2004;Mehner et al. 2011). Accordingly, strong inter-annual effects on their zooplankton prey can be expected. We hypothesized that the size distributions of zooplankton would show strong differences between years, which reflect the relative contributions of fish predation and resource availability. Therefore, we aimed to demonstrate that the size distribution of pelagic zooplankton communities is a superior and easy-to-generate metric that informs about the major structuring forces in the pelagic area of deep lakes.

Study site and sampling
This study was conducted in Lake Stechlin, located approximately 120 km north of Berlin (Germany). This mesotrophic lake covers about 4.25 km 2 with mean and maximum depths of 22.3 m and 69 m, respectively (see Fig. 1a, Koschel and Adams 2003). A total of 13 fish species inhabit the lake, but the pelagic fish community is dominated by two sympatric species of ciscoes, Coregonus albula and C. fontanae (more than 95% of fish abundance; Mehner and Schulz 2002;Schulz and Freyhof 2003;Helland et al. 2007). Furthermore, the average stock exploitation rate by the single commercial fishery is 0.08 (Wanke et al. 2017), suggesting no strong effect of fish removal on fish abundance differences between years. We estimated counts and individual body wet mass (in g) of fish using a pelagic midwater trawl annually at the end of June between 2005 and 2018 (n = 13 years, data for 2016 missing due to net malfunction). Trawling is commonly employed for SS studies on fish, but its inherent size selectivity (Klein et al. 2019) may result in an overrepresentation of medium-sized fish. Therefore, simultaneously to trawling, we also performed hydroacoustic surveys to obtain pelagic fish mean mass (in g) and biomass (kilogram per hectare, see Fig. 1b).

Crustacean assessment
Crustacean plankton was collected by vertical tows in the upper depth layer from 0 to 22 m and the lower depth layer from 22.1 to 65 m, using a cone-shaped closing net (mesh size 90 lm, opening 27 cm 2 ; length 1.2 m, Hydrobios, Kiel, Germany) at the deepest location of the lake (NE basin, 69.5 m). Since fish data were exclusively collected in summer, we similarly considered only zooplankton samples collected during summer months between 2005 and 2018 (13 years, excl. 2016). Therefore, from June to August, the lake was sampled twice a month except 2014 (one sampling in August), 2015 (each month sampled only once) and 2018 (one sampling in August) resulting in a total of 148 zooplankton samples (74 samplings 9 2 depth layers). Samples were preserved in a 4% sugarformaldehyde solution. Taxa were identified according to Flößner (1972Flößner ( , 2000, Kiefer and Fryer (1978), Einsle (1993) andLieder (1996). A sub-sample containing at least 100 individuals of the dominating species or groups was counted (abundance in Ind m -3 ) using an inverted microscope at 60-fold magnification. Mean length of zooplankton was estimated by measuring approximately 15 individuals per each frequent group (species level in most cases; see Table 1 in Appendix for list) and then transformed to carbon body mass (lg C Ind -1 ) using species-specific allometric equations from published studies (Bottrell et al. 1976;Kasprzak 1984). To examine the taxon-specific proportions from the entire zooplankton community in the upper and lower depth layers, we calculated the biomass of each taxon group for each depth layer. Abundance (Ind m -3 ) was multiplied with the carbon body mass for each taxon and the group-specific annual average biomass was determined from the individual samples.

Fish assessment: trawling
A pelagic trawl with 28/20/10 mm mesh size, 10 mm mesh size in the cod-end and a total length of 14.8 m (stretched on land) was utilized. The net with an opening area of approximately 10 m 2 (opening width 3.5 m) was towed by a boat (length 7 m, width 2 m), which was driven by a 60 hp engine over four longitudinal transects in the deepest lake basins (see Fig. 1b). Trawling speed (mean ± SD) was 6.5 ± 0.6 km h -1 (1.8 ± 0.2 m s -1 ), while towed distance, as the product of trawling speed and trawling time, ranged between 500 and 1600 m (mean 840 m) with an average towing time (± SD) of 7.8 min (± 1.9). During each survey, a total of four hauls were conducted at approximately 12 m ± 0.65 m, 15 m ± 1.09 m, 25 m ± 2.08 m and 32 m ± 0.90 m (mean ± SD), whereby the two deeper hauls had shorter trawling times and slightly slower hauling speeds. Due to limited spatial extension of water layers with more than 35 m depth, deeper hauls could not be performed. The actual sampling depth was recorded with a diving computer attached on the trawl's head rope. The depth variability during each tow never exceeded 3 m. At the end of each haul, the trawl was quickly lifted by hydraulic winches at speeds comparable to the trawling speeds, thus preventing escapement of fish from the net. Fish were counted, and individually measured to determine total length (TL, nearest centimetre) and body mass (wet mass, wm, nearest gram). Only subsamples were measured and weighed when catches were too large. The abundance of fish (Ind 1000 m -3 ) was estimated from the towed distance and the opening area of the trawl. Due to the shorter hauled distance in the deeper layers, the overall size distribution per year calculated as the total sum of all fish caught by the four hauls would have been biased towards the size of fish in the upper water layers. Therefore, we employed a weighting factor to calculate a weighted mean abundance and size distribution across the four depth layers. The weighting factor (WF) was calculated as the ratio between the maximum volume fished by one among the four hauls in this year and the volume fished in the respective depth. Accordingly, the WF equals the abundance and size distribution per depth layer to identical fished volumes. Some size distribution analyses require fish numbers per size class as integers; then the weighting factor was rounded to the nearest half integer (see Table 2 in Appendix for an example).

Fish assessment: hydroacoustics
Hydroacoustic surveys were completed as a series of transects using a SIMRAD (Kongsberg) EY60 splitbeam hydroacoustic unit (120 kHz, circular transducers, beam width 7°9 7°, pulse duration 128 ms, ping rate 3 pings s -1 ). After calibration with standard spheres provided by the manufacturer, the surveys were conducted during complete darkness starting 1 h after sunset every June between the 6th and 29th. Surveys were performed with vertical beaming along 20 transects (total distance about 12 km) across this tri-basin lake, with transect length ranging from 635 to 1332 m (see Fig. 1b). Data were stored in a computer, processed and analysed by the postprocessing Sonar 5 Pro software (version 6.0.4-R340, Balk and Lindem 2017). The lower target strength (TS) threshold (as a measure of fish body length in dB) was set to -55 dB for the echogram corresponding to fish of an approximate total length of 4.2 cm and 0.42 g wm. Parameters for amplitude echograms were set 6 dB lower (-61 dB) to accept targets out of the half power edge of the sound beam. The TS (dB) of the targets was converted into fish total length (cm) and body wm (g) using the equation obtained for Coregonus spp. from Lake Stechlin by Mehner (2006a, b): The records from all 20 sampled transects were combined and averaged for each year. Fish biomass (kg ha -1 ) was obtained by echo-integration, by using the size-frequency distribution of the detected individual targets per transect.

Environmental descriptor assessment
Physical and chemical data were recorded at the deepest site of the lake, simultaneously with the crustacean data. Vertical profiles of water temperature (°C) were determined every 5 m between the surface and 65 m, while TP (mg L -1 ) concentrations were measured from samples at 0, 5, 10, 40, 60 and 65 m depth photometrically. TP was measured after wet digestion of unfiltered subsamples in an autoclave (potassium peroxodisulfate, 30 min, 134°C). TP content was determined by flow analysis and spectrometric detection (Foss FIAstar 5000 Analyzer). Chl a (mg L -1 ) concentration was measured at 0, 5 and 10 m depth using photometric determination (Hitachi spectrophotometer U-2900). The arithmetic means of temperature and TP per year were determined for the upper (0-20 m) and lower (21-65 m) layers separately, while the mean chl a concentration was only based on samples from the upper layer. However, since chl a concentration in shallower depths may have an effect on processes of the entire water column, the determined chl a mean was also utilized as a predictor for dependent variables in the lower layer.

Data analyses
Zooplankton abundances (Ind m -3 ) were classified into a total of 17 body mass size classes, ranging from 0.1 to 25 lg C Ind -1 (log 2 range: -3.36 to 4.63) and logarithmically transformed (base 2) for each sampling in the upper and lower depth layer (total of 148). The arithmetic mean zooplankton abundances per size class based on all samplings in the upper or lower depth layers for each year were determined and log transformed. Therefore, zooplankton SS were calculated for each depth layer separately per year (total of 26 SS, 2 for each of the 13 years). We utilized the LT (log-transform) method where the log 2 mean counts were aggregated in log 2 size bins of constant width and plotted against the midpoint of the size interval. Ordinary least square regression (OLS) method was used to estimate the SS slope and intercept (Rice and Gislason 1996). The intercept of size spectra is a general measure of total organism abundance. Normality of data was visually examined and tested with a Shapiro-Wilk normality test.
For scrutiny, we also used the maximum likelihood estimation when fitting zooplankton size spectra (Edwards et al. 2017). Here individual zooplankton body masses were not binned, but utilized in continuous form. The exponent of the individual size distribution (b), known to be positively related to the slope of the SS, was determined for each summer sampling in the upper and lower layers. Then the arithmetic mean of exponent b and the confidence interval were determined for each depth layer per year from these samplings. Furthermore, we determined zooplankton size diversity (l) as based on carbon body mass (lg C), calculated from each sample and then averaged for each depth layer per year following Quintana et al. (2008), employing the Shannon diversity expression but adapted for body size as a continuous variable. Finally, we calculated the total biomass and median body mass of only those zooplankton taxa that are commonly included in Coregonus diet, which encompasses mainly cladocerans, including Daphnia spp., Bosmina spp., Chydorus sphaericus, Diaphanosoma brachyrum, Ceriodaphnia quandrangula and nauplii  in each depth layer as a response variable. However, we did not calculate SS for the limited prey taxa as these only represent a fraction of the sizes of the entire community. Therefore, one cannot expect a negative relationship between size and abundance as predicted by SS theory for the entire zooplankton community.
Fish size distribution from trawling was comparably described as continuous SS using maximum likelihood estimation of b, which is the recommended method for fitting fish SS (Edwards et al. 2017). This was based on depth-weighted fish abundances. In addition, the arithmetic mean biomass (g 1000 m -3 , depth-weighted average from four trawl hauls) was calculated for each year. Finally, fish biomass as kilogram per hectare (kg ha -1 ) and mean mass (in g) of fish were calculated from the hydroacoustic surveys.
To identify relationships between zooplankton community parameters (size distribution) and size or biomass of fish and mean TP, chl a and temperature, we utilized multiple regression analyses (see Table 1). Response variables tested in the multiple regressions were zooplankton slope and intercept (estimated by OLS) and zooplankton mean size diversity index, both per depth layer and year, annual summer averages of exponent b (estimated by MLE) of the SS, median body mass (lg C) and mean abundance (Ind m -3 ) of potential prey taxa in the upper and lower layer. As independent variables we utilized the exponent b of the fish SS (estimated by MLE) and the mean fish biomass (g 1000 m -3 , both from trawling), as well as mean mass (g) and biomass (kg ha -1 ) as based on hydroacoustic surveys. Water temperature and mean chl a and TP concentrations were also included as environmental and productivity predictors, respectively, in all regressions. Only predictor variables that were not strongly cross-correlated (Pearson correlation coefficient below 0.5, and Variance Inflation Factor (VIF) below 2) were included in the analyses (Belsley et al. 1980), as multicollinearity can result in inflated variance among predictors in the model. Automatic stepwise backwards selection and the Akaike information criterion (AIC) were used to select the variables in the final model by accounting for increasing complexity (as seen in Akaike 1974). The most parsimonious model obtained was the combination of variables that produced the lowest AIC. Residual partial plots for each significant variable were plotted to show the relationship between the predictor and the response variables. All statistical Specifically the MLE's of fish and zooplankton size spectra were performed by the provided R code and package (sizeSpectra) from Edwards et al. (2017) and Edwards (2020).

Results
Zooplankton community structure, abundance and biomass (2005-2018, excl. 2016) During the summer months, the zooplankton community of Lake Stechlin was numerically dominated by copepod nauplii and Bosmina longirostris.  Table 3 in Appendix). However, total zooplankton community biomass in the upper layer was clearly dominated by the calanoid copepod Eudiaptomus gracilis (10-30%), with Daphnia cucullata (8-24%) and nauplii (12-25%) on the second and third rank. In the lower layer, Eurytemora lacustris (27-76%) and Bosmina longirostris (4-47%) contributed most to the overall zooplankton community (see Table 2). The body masses of zooplankton ranged from Chydorus sphaericus (0.09 lg C Ind -1 ) in the upper layer to Cyclops spp. (24.80 lg C Ind -1 ) in the lower depth layer. The mean (± SD) body length and mass of the entire zooplankton communities were higher in the lower (x 0.54 ± 0.01 mm; x = 1.42 ± 0.08 lg C bm) than in the upper depth layer (x 0.50 ± 0.007 mm; x = 0.95 ± 0.03 lg C bm).
Zooplankton spectral slope and shape and MLE b The zooplankton size distributions substantially deviated from a linear negative log-log relationship, with only six out of the 26 (2 layers 9 13 years) linear regressions being significant (p \ 0.05, see Table 4 in Appendix). In most of the years, log size and log mean abundance exhibited a bimodal pattern (see Figures 1-4 in Appendix). The highest mean abundance per body mass size-class was often observed in the log 2 body mass classes -1.75 and -1.25 (0.20-0.41 lg C) dominated by Bosmina longirostris and nauplii and log 2 body mass class 1.25 and 1.75 (1.68-3.35 lg C) dominated by Eudiaptomus gracilis, Daphnia cucullata and Eurytemora lacustris (see Table 3 in Appendix). Zooplankton spectral intercepts differed significantly between both depth layers. The mean intercept in the upper layer was significantly higher compared to the lower layer (t test: t (24) = -15.09, p \ 0.0001). In contrast, mean spectral slopes were identical between upper and lower layers (t-test: t (26) = -0.0005, p = 0.99). The mean exponent b of the SS estimated by MLE was slightly higher in the upper than in the lower depth layer (t test: t (155) = 0.68, p = 0.5). Zooplankton mean size diversity was significantly higher in the lower than in the upper layer (t test: t (120) = 3.52, p \ 0.0006, see Table 3).

Fish abundance and range between 2005 and 2018
Fish size distributions varied strongly between years, while the individual wet mass of fish sampled during trawling ranged from 0.018 to 103 g. Some years were characterised by numerous small young fish (less than 4 g wet mass, e.g.

Effect of fish abundance and environmental descriptors on zooplankton community
Among the 48 combinations for 6 response variables and the 4 different fish predictors plus the 3 environmental predictors used in all analyses in two depth layers, the majority (45 out of 48) of the regressions were not significant (see Table 5 in Appendix). Only the mean zooplankton spectral intercept (estimated by OLS) was significantly positively related to the exponent b of the fish SS from trawling, and was positively related to TP and negatively related to chl a in the upper depth layer (R 2 = 0.64: F = 5.5, df = 9, p = 0.01; see Fig. 4). No similar correlations were found for the lower layer. Fish biomass (kg ha -1 ) based on hydroacoustics had a significantly negative effect on median body mass of potential fish prey taxa in the upper (R 2 = 0.3, F (11) = 4.8, p = 0.049) and lower depth layers (R 2 = 0.53, F (11) = 12.69, p = 0.004, see Fig. 5).

Discussion
In our study, we aimed to determine top-down and bottom-up drivers for zooplankton size distribution in the pelagic area of the mesotrophic Lake Stechlin during summer months across 13 years of sampling. We tested seven commonly used parameters to describe the zooplankton community size structure and density from several successive years, by expecting that potentially strong effects from fish size or biomass and the environmental predictors would become visible on all or many of the six zooplankton parameters. Surprisingly, however, we found almost no effect of annual differences in fish size, biomass, TP, chl a and temperature on the zooplankton size structure. This result seems to be robust because we expressed the zooplankton size structure as linearly binned SS, as continuous SS and as size diversity. Furthermore, the fish size and biomass data were obtained by two independent approaches, midwater trawling and hydroacoustics. We conclude that feeding by the zooplanktivorous coregonids in Lake Stechlin does not affect the size distribution of the entire pelagic zooplankton community. However, the few significant regressions indicate that there is a limited top-down effect on the density and size of those zooplankton taxa that constitute the major prey of the coregonids.
Our analyses have shown that fish size distributions as based on the exponent b of the continuous SS obtained from trawling were significantly positively correlated with the spectral intercept of the zooplankton community in the upper layer of the lake. High b b Fig. 2 Mean abundance (Ind m -3 ) of zooplankton taxa grouped as: (1) Bosmina spp., (2) Chydorus spp. and Diaphanosoma spp., (3) calanoids, (4) copepodits and nauplii, and (5) Daphnia spp. for 2005-2018 (excl. 2016) in the a upper and b lower depth layer of Lake Stechlin exponents indicate a dominance of larger fish in the size distribution. Therefore, years with a relatively high number of medium and large-sized fish (e.g., in 2014 and 2015) were associated with a high density of zooplankton organisms (reflected by the intercept of the zooplankton SS). In contrast, years characterised OLS means size spectrum estimated by ordinary least squares, while MLE indicates fit of continuous size spectrum by maximum likelihood Fig. 3 Depth-weighted counts of fish distributed into body wet mass (g) classes as obtained from four trawling hauls in Lake Stechlin for each sampling year by numerous small fish (low exponent b, e.g. in 2012 and 2017) were associated with lower zooplankton densities. This positive relationship between zooplankton density and the contribution of larger coregonids to the size distribution may indicate that zooplankton density might be top-down controlled by fish feeding only in years with a dominance of small coregonids. Furthermore, a direct feeding effect of fish was reflected by the negative correlations between hydroacoustically estimated fish biomass and the median body mass of zooplankton prey taxa. These negative correlations were found in both the upper and lower layers. These results suggest that fish feeding might affect size structure and density of the zooplankton community, but with a stronger effect in the upper layers, and with a stronger effect on those cladocerans that constitute the prey of the coregonids .
Cladocerans were more frequently found in the upper than in the lower water layer. Coregonids feed mainly during the day in the deeper hypolimnetic layers where light intensity is sufficient to allow for selective feeding (Gjelland et al. 2009;Mehner et al. 2010) and while migrating between the lower and upper hypolimnion close to the thermocline at low light conditions during dusk and dawn (Scheuerell and Schindler 2003;Gjelland et al. 2009). Therefore, coregonids in Lake Stechlin perform diel vertical migration for efficient zooplankton feeding, as the great majority of the cladoceran zooplankton in Lake Stechlin concentrate between the surface and 25 m depth (Kasprzak and Schwabe 1987). Our data were collected during summer months when the higher water temperatures that fish experience close to the thermocline, may result in higher daily basal energy requirements. It was already suggested by Mehner (2006a) that hungry juvenile coregonids prolong their stay in the shallower water layers, despite higher predation risk, to cope with the increasing energy demands during summer. Therefore, they may control the zooplankton density and size in the upper layers during years with a dominance of juvenile coregonids. Previous studies in Lake Stechlin suggested that the top-down control of the zooplankton biomass by fish predation was unlikely, as trophic transfer efficiency between zooplankton and fish was determined to be relatively low (6%, Schulz et al. 2004). Although prey taxa dominated the upper layer in numbers, we found that prey and non-prey taxa contributed equal amounts of biomass to the zooplankton community in the upper layer. Furthermore, it was shown in a previous study that specifically in summer months pelagic coregonids feed preferentially upon Bosmina longirostris and B. coregoni , whilst these taxa contribute only 9% to the total zooplankton community biomass in the upper layer. Therefore, even if they are eaten by fish, their removal has only a minimal impact on the size distribution of the entire zooplankton community. Even more drastically it can be seen in the lower layer, where prey taxa of coregonids likewise dominate in numbers but they contribute only by 28% to the biomass of the zooplankton community. In contrast, the much bigger Eurytemora lacustris, which is not fed upon by the fish, contributes 44% to the community biomass alone. The data obtained here support that the size distribution of the entire zooplankton community is not top-down controlled by fish feeding. There are signs of negative correlations between zooplankton prey size and biomass of smaller coregonids, but these effects are obviously not strong enough to induce traceable changes in the size structure of the entire zooplankton community. These results support that pelagic trophic cascades are weak in lakes with lower productivity (Elser and Goldman 1991;Persson and Diehl 1992;Sarnelle 1992;Mehner et al. 2008).
As the productivity of Lake Stechlin has increased during the last decade (see Figures 9, 10 in Appendix), we expected that chl a and TP concentrations would likewise affect zooplankton size distributions across the years. However, we found only a weak positive effect by TP on zooplankton density (reflected by the intercept of the size spectra), whereas chl a was negatively correlated with the spectral intercept in the Fig. 5 Partial residuals of the dependent variable median body mass of zooplankton taxa constituting coregonid prey and the independent variables fish biomass (kg ha -1 ) from hydroacoustic surveys in the upper (a) and lower (b) depth layers of Lake Stechlin upper layer. This result is in part counterintuitive because we would have expected a positive correlation between zooplankton resources (phytoplankton biomass expressed as chl a) and zooplankton density. However, a recent long-term study from Lake Stechlin found evidence for a decreasing efficiency in trophic coupling between zooplankton and phytoplankton caused by an increasing dominance of filamentous cyanobacteria in the phytoplankton community (Selmeczy et al. 2019). Cyanobacteria are a known inadequate food source for zooplankton (Gulati and DeMott 1997) and therefore their dominance reduced the resource use efficiency and biomass of zooplankton. This dynamics may explain the significant negative relationship of chl a and zooplankton spectral intercept found in the upper layer in our study. However, these correlations became only visible in the multiple regression when the exponent of the fish SS was likewise included. Therefore, another plausible explanation would be that in years with a dominance of small coregonids, the zooplankton density (primarily cladocerans) is diminished, and the lower cladoceran density cascades down to a higher phytoplankton biomass in the upper layer of Lake Stechlin. Overall, however, there is a long-term signal that the increasing TP concentrations may increase the zooplankton densities via bottom-up control. In a recent comparative study across more than 200 lakes it has been clearly shown that chl a concentrations are more strongly bottom-up correlated with fish production rates, than fish feeding rates affect chl a concentrations negatively in the sense of top-down control (Bartrons et al. 2020).
There are several insights from the results that might be important for food web studies beyond Lake Stechlin. Although coregonids are strictly planktivorous and may achieve high biomasses in deep lakes, they seem to have only a limited capacity to induce pelagic trophic cascades down to phytoplankton. There are two reasons for their weak predation effects. First, coregonids have a relatively small gape size, and hence prefer small zooplankton prey such as Bosmina spp. . In turn, they do not feed upon the large calanoid and cyclopoid copepods, which often dominate in the colder hypolimnetic layers of deep lakes. This strong size and species selectivity prevents that the total zooplankton density is massively affected by coregonid predation. Second, there is a spatial disconnect between coregonids and their preferred prey. While the coregonids prefer the colder hypolimnetic zones with temperatures around 10-12°C and ascend into warmer water only during the night (Mehner et al. 2010), cladocerans grow best in the warm epilimnetic layers (Helland et al. 2007). Therefore, fish have only a limited time window during which they can feed upon their preferred prey. This short feeding period may explain why there is only a weak top-down effect of coregonids even on their preferred prey taxa.
Ultimately, we have to revise our expectation that zooplankton size distributions are sensitive to the feeding effect of fish and to environmental predictors. Although species-and size-selective feeding of coregonids is likely and has effects on the densities of the preferred taxa, these effects do not translate into a statistically significant modification of the entire SS. Instead, the SS seems to be relatively robust against predation effects, and seems to reflect the lake-wide energy availability and transfer efficiency in the food web (Mehner et al. 2018). This suggests that the classical pelagic food chain from fish to phytoplankton is only part of the lake food web. Creating SS across the majority of the organismal groups in both pelagic and benthic habitats might therefore allow inferring the major drivers of food web configurations. Furthermore, additional information on phytoplankton size distribution and taxonomy would allow a better mechanistic understanding of the interaction between phytoplankton and zooplankton communities. In our study, we used exclusively the chl a concentration as an indicator for phytoplankton biomass. However, many zooplankton species feed in a size-and speciesselective way. For example, juvenile daphniids usually avoid green filamentous algae, such as Ulothrix, whereas larger individuals include it in their diet (Chen et al. 2011). Shifts in phytoplankton composition and size structure, such as higher biomass of filamentous cyanobacteria induced by global warming (Selmeczy et al. 2019), may modify the strength of top-down effects of zooplankton on phytoplankton (Ersoy et al. 2017).
Our results are based on correlational evidence, which not necessarily allows implications of causal relationships. In particular, data on stomach content analyses or stable isotopes of fish would allow a better understanding of the strength of trophic links that ultimately affect the size distributions (Brose 2010;Boukal 2014). Furthermore, despite having included data from 13 successive years, in the framework of multiple linear regressions this is a low number of observations, and hence is associated with low statistical power (Hair et al. 2014). Accordingly, the overall surprisingly low number of significant multiple linear regressions here has to be considered cautiously. However, our study on pelagic trophic interactions and drivers of zooplankton SS may inspire others with significantly longer datasets to repeat a critical reconsideration of top-down and bottom-up effects on zooplankton size structure in stratified lakes.