Increased Nitrogen Loading Boosts Summer Phytoplankton Growth by Alterations in Resource and Zooplankton Control: A Mesocosm Study

Abstract

The effectiveness of controlling nitrogen (N) to manage eutrophication of aquatic ecosystems remains debated. To understand the mechanisms behind phytoplankton growth in shallow lakes (resource and grazing effects) under contrasting N loading scenarios, we conducted a 70-days mesocosm experiment in summer. The mesocosms contain natural plankton communities deriving from a 10-cm layer of lake sediment and 450 L of lake water. We also added two juvenile crucian carp (Carassius carassius) in each mesocosm to simulate presence of the prevailing omni-benthivorous fish in subtropical lakes. Our results showed that N addition increased not only water N levels but also total phosphorus (TP) concentrations, which together elevated the phytoplankton biomass and caused strong dominance of cyanobacteria. Addition of N significantly lowered the herbivorous zooplankton to phytoplankton biomass ratio and promoted the phytoplankton yield per nutrient (Chl-a: TP or Chl-a: TN ratio), indicating that summer N addition likely also favored phytoplankton growth through reduced grazing by zooplankton. Accordingly, our study indicates that summer N loading may boost eutrophication via both changes in resource and grazing control in shallow lakes. Thus, alleviation of eutrophication in shallow eutrophic lakes requires a strategic approach to control both nutrients (N and P) appropriately.

Introduction

Nitrogen (N) and phosphorus (P) are key limiting factors determining autotroph growth in both terrestrial and aquatic ecosystems (Elser et al., 2007; Conley et al., 2009). The N:P ratio of reactive N and P anthropogenic inputs to the biosphere has dramatically increased during the past half century, inducing significant alterations in the organisms and ecosystems of Earth (Peñulas et al., 2012; Peñulas et al., 2013; Tong et al., 2020).

In shallow eutrophic lakes, in-lake TN:TP ratios are often low in summer due to the high denitrification and strong P release at high temperatures (Xu et al., 2010; Xu et al., 2021; Søndergaard et al., 2013; Søndergaard et al., 2017; Shatwell and Köhler, 2019). From a resource availability perspective, excess N loading could directly increase the summer phytoplankton biomass via alleviation of temporal N deficiency, and this has been confirmed by numerous bioassays and mesocosm experiments (Xu et al., 2010; He et al., 2015; Kolzau et al., 2018). Nitrogen loading may also influence the internal P loading via feedback mechanisms between phytoplankton blooms and sediment P release (Xie et al., 2003; Gao et al., 2014; Chen et al., 2018). That is, a higher phytoplankton biomass due to summer N enrichment augments the P release from P-rich sediments via various mechanisms including increasing diffusion, bacteria activity, pH, and likely also degradation of settled algae (Xie et al., 2003; Gao et al., 2014; Chen et al., 2018; Ma et al., 2018; Ma et al., 2021). This released P, in turn, hinders a change of the system towards P limitation at persistent external N loading and, therefore, it creates a positive feedback loop benefitting phytoplankton growth (Figure 1). Nonetheless, high doses of nitrate, may also prevent P release by improving oxidation of the sediment in shallow lakes (Andersen, 1982; Boström and Pettersson, 1982). For instance, a long-term investigation of 22 Danish shallow lakes revealed absence of P release if in-lake nitrate concentrations exceeded about 0.5 mg L⁻¹ (Andersen, 1982). However, this pattern was not supported in studies of four Danish shallow eutrophic lakes where summer nitrate additions increased the P release from the sediment (Jensen and Andersen, 1992). Thus, the dual effect of oxidized nitrogen on sediment P release makes prediction of the effects of summer N loading in shallow lakes uncertain.

FIGURE 1

Excess N loading may also favor phytoplankton growth via reducing the fitness and biomass of herbivorous zooplankton (Trommer et al., 2017; Lorenz et al., 2019). For instance, mesocosm experiments found that N addition caused stoichiometric shifts in the biomass of primary producers towards higher N:P and C:P ratios, which lowered cladoceran and total zooplankton biomasses but favored the dominance of rotifers (Lorenz et al., 2019). In shallow lakes, although P has traditionally been considered the primary nutrient causing harmful cyanobacterial blooms, recent findings from laboratory studies and lakes throughout the world demonstrate that N can be important in controlling the timing, density, and toxicity of some non-diazotrophic cyanobacterial blooms (reviewed by Gobler et al., 2016). Beaulieu et al. (2013) assessed cyanobacteria blooms in 1,147 lakes and reservoirs with differing trophic status across the U.S. and found that the best predictors of these events were TN and water temperature. Multi-year observations from western Lake Erie in United States demonstrate that microcystin concentrations peak in parallel with inorganic N concentrations and are significantly lower during years of reduced inorganic nitrogen loading and concentrations (Gobler et al., 2016). Thus, in cyanobacteria-dominated shallow lakes, it can be inferred that summer N loading might create high cyanobacterial biomass and cyanotoxin concentrations. As numerous studies have recorded that bloom-forming cyanobacteria generally are a poor food resource for both cladocerans and copepods due to their low nutritional value, e.g., insufficient amounts of essential lipids like PUFAs and sterols (Caramujo et al., 2008; Martin-Creuzburg et al., 2008), poor manageability (long filaments or aggregate colonies), and toxicity (DeMott and Moxter, 1991; Kirk and Gilbert, 1992), the feedback between high cyanobacteria blooms and low zooplankton grazing driven by summer N loadings may produce another “vicious cycle” that expectedly boosts eutrophication (Figure 1).

To obtain a better mechanistic understanding of the effects of summer N loading on the resource and grazing control of phytoplankton growth in shallow eutrophic lakes, we conducted a 70-days mesocosm experiment involving contrasting N loadings at the shore of Lake Taihu (China), one of China’s best studied lakes and known for its serious non-diazotrophic cyanobacterial (mainly Microcystis spp.) blooms (Guo et al., 2007; Qin et al., 2019). Previous bioassay experiments and field monitoring studies have revealed that the phytoplankton in Lake Taihu is generally N limited or N, P co-limited during summer and fall (Xu et al., 2010; Xu et al., 2021; Paerl et al., 2011). Thus, our hypothesis was that with increasing N loading in summer, the phytoplankton biomass would increase due to the alleviation of N deficiency and higher sediment P release. We also expected that N addition would increase the biomass and dominance of grazing-resistant cyanobacteria and thereby reduce zooplankton herbivory, which, in turn, would favor phytoplankton growth via release of grazing control.

2 Methods and Materials

2.1 Experimental Design

The 70-days outdoor mesocosm experiment was conducted from June 8 to August 17, 2020, at Taihu Lake Laboratory Ecosystem Research Station (TLLER), located in Meiliang Bay on the northern edge of Lake Taihu, China. On May 30, six cylindrical, high-density polyethylene tanks (95 cm height × 87 cm upper diameter × 77 cm bottom diameter) were filled with a 10 cm layer of lake sediment and 450 L of water collected from Meiliang Bay. The sediments had previously been sieved (mesh size: 1.7 mm) to remove large invertebrates and mixed to ensure uniformity. The prepared sediment contained natural quantities of organism resting stages, which allowed zooplankton to hatch and supported community development during the experimental period. The water was also screened using a 64-μm mesh filter to remove macroinvertebrates and inorganic particles before use. To simulate the prevailing omni-benthivorous fish in subtropical lakes, two juvenile crucian carp (Carassius carassius) were added to each mesocosm. Each carp had an average wet weight of 1.95 ± 0.07 g and an average body length of 5.17 ± 0.13 cm. In each mesocosm, the fish approximated a natural density of 66 kg ha⁻¹, which is often recorded in Chinese subtropical or tropical lakes, for example Lake Taihu (Mao et al., 2020).

When the experiment began, nitrate (KNO₃) was supplied weekly to three mesocosms throughout the experiment. The daily N loading was equal to 130 μg N L⁻¹ day⁻¹, which was within the range of external loading in Lake Taihu (Paerl et al., 2011). The total amount of N added in each of the N treatments was equal to 4.095 g during the experiment. The remaining three mesocosms functioned as controls.

2.2 Sampling and Analytical Methods

2.2.1 Water Physico-Chemical Properties and Chlorophyll a

Samples were collected biweekly (every 14th day) during the experiment. Prior to sampling, pH and dissolved oxygen (DO) concentrations were measured in situ by a portable instrument (YSI 9500, YSI Inco, Yellow Springs, United States). Thereafter, a tube sampler (8 cm diameter, 64 cm length) was used to collect 2-L depth-integrated water samples from each mesocosm. These water samples were subsequently taken to the laboratory for analyses of nutrient and chlorophyll a (Chl-a) concentrations.

Total suspended solids (TSS) was determined from 100 to 200 ml water samples filtered through pre-combusted (450°C for 2 h) and pre-weighed GF/C filters, which were then dried to a constant weight at 60°C for 48 h. After determining TSS, the filters were combusted in a muffle furnace at 550°C for 2 h, cooled in a desiccator, and finally weighed to determine the level of inorganic suspended solid (ISS). The organic fraction (OSS) was calculated by subtracting ISS from TSS. Total nitrogen (TN), total phosphorus (TP), total dissolved nitrogen (TDN), and total dissolved phosphorus (TDP) were determined through colorimetry after digestion with K₂S₂O₈ and NaOH solution (APHA, 1998). Particulate nitrogen (PN) or particulate phosphorus (PP) was calculated by the differences of TN and DTN and TP and DTP, separately. We calculated TN: TP and PN: PP ratios (by weight) to assess the nutrient limitation of phytoplankton (Redfield, 1958; Guildford and Hecky, 2000). We calculated DTP: TP ratios (by weight) to explain the potential mechanism of P release (Ma et al., 2021). Chl-a concentrations were measured spectrophotometrically from matter retained on a GF/C filter and extracted in a 90% (v/v) acetone/water solution for 24 h. No correction was carried out for pheophytin interference (SEPA, 2002).

2.2.2 Sediment P Release

We measured the sediment P profile on July 6 (day 28 of the experiment) and August 17 (day 70 of the experiment). After water sampling, a plexiglass cylinder with an internal diameter of 3 cm and a length of 60 cm was used to extract a sediment core in the middle of each mesocosm. We sliced off the top 4 cm of each sediment core and measured its dry weight after drying it in the shade. Total-P of the sediment was measured according to the methods introduced by Ruban et al. (1999). P forms were also extracted according to Hupfer et al. (1995), Rydin (2000), and Meis et al. (2012). Dried sediment (0.5 g) was added to 50-ml centrifuge tubes and sequentially extracted with 1) 1 M NH₄Cl (NH₄Cl-P); 2) 0.11 M NaHCO₃ and 0.11 M Na₂S₂O₄ (BD-P); 3) 0.1 M NaOH (NaOH-P); and 4) 0.5 M HCl (HCl-P); 5) the residual P (Residual-P) was calculated by subtracting all the extractable P forms from Total-P. For each mesocosm, we calculated the P concentration (mg m⁻²) using the formula:where P_t is the P content (mg kg⁻¹) on the sampling date t (Day 28 or Day 70), DW_t is the dry weight of sediments on date t, and Area represents the sampling area of the sediments.

2.2.3 Phytoplankton, Zooplankton, and Top-Down Control

Besides measurements of Chl a, phytoplankton communities were recorded at the end of the experiment (Day 70). For this purpose, a 500 ml depth-integrated water sample from each mesocosm was treated with 5 ml Lugol’s iodine solution and sedimented for 48 h. The supernatant was removed and the residue collected and examined at 100–400x magnification for enumeration. Phytoplankton was identified to genus level following Hu and Wei (2006) using, as far as possible, recent taxonomic revisions (Guiry and Guiry, 2020). The biomass of common phytoplankton taxa was calculated based on cell size measurements of at least 30 cells of each taxon using formulae for geometric shapes approximating cell forms (Zhang and Huang, 1991). For less common taxa, biomass calculations were based on fewer measurements; otherwise, at least 1,000 cells of each taxon were counted per sample.

Zooplankton were sampled on days 0, 28, 56, and 70 of the experiment. A 10 L depth-integrated water sample was filtered through a 64 μm net and the retained material was preserved in Lugol’s iodine solution. The zooplankton included cladocerans, copepods, and rotifers that were counted at magnifications between×40 and ×100. Species were identified according to Chiang and Du (1979) and Shen and Du (1979). Copepod nauplii were counted without further taxonomic distinction. Where possible, up to 20 individuals of each taxon were measured. Biomass (dry weight) was estimated by length–weight relationships (Bottrell 1976; McCauley et al., 1984).

To evaluate the grazing pressure of zooplankton on the phytoplankton, we calculated three top-down metrics: the ratios of zooplankton biomass to Chl-a concentration, Chl-a to TP, and Chl-a to TN ratios. We also calculated the density-weighted body mass of the zooplankton community and the biomass ratio of cladocerans to copepods to assess the predation strength of fish on zooplankton as large-sized individuals and taxa with low mobility (e.g., cladocerans) were more sensitive to fish predation (Jackson et al., 2007; Jeppesen et al., 2011). These metrics (when combined) are considered to be good indicators of top-down control and are widely used in long-term monitoring, cross-system lake comparison, and mesocosm studies (Gyllström et al., 2005; Jackson et al., 2007; Jeppesen et al., 2011; Liu et al., 2018; He et al., 2018, 2021).

2.3 Data Analysis

All data processes and statistical analyses were carried out in R (R Core Team, 2020). We studied the influence of N addition (fixed effects) on nutrient concentrations and ratios, pH, DO, Chl-a, zooplankton biomass, and top-down metrics using generalized linear mixed models with time and mesocosm identity as random effects (GLMM; Bolker et al., 2009). For each response variable, we assessed the relative fit of the two different candidate models (fixed effect: N addition and null) and selected the best model based on the lowest AIC value from the default anova function. We assumed Gaussian error distributions for all response variables and fitted modes using the glmmTMB function from the glmmTMB package (Brooks et al., 2017). We also used generalized linear mixed models to compare the differences of sediment P contents between both N scenarios at both sampling dates. To explore whether the difference of P contents between both N scenarios changed temporally, we assessed the relative fit of the five different candidate models (fixed effect: N addition * Time, N addition + Time, N addition, Time and null) for each P fraction and selected the best model based on the lowest AIC value from the default anova function. Phytoplankton biomass at the end of experiment was compared by applying Student’s t test (function t. test) after normality (function shapiro. test) and homoscedasticity (function bartlett. test) tests.

3 Results

3.1 Water Physico-Chemical Properties

Model selection by GLMMs confirmed significant positive effects of N addition on the TN、TDN, and TP concentrations, while its effect on the TDP concentration was not significant (Table 1; Figure 2). Model selection by GLMMs confirmed that N addition significantly promoted water pH, DO, and DTP: TP (Table 1; Figure 2). N addition did not significantly elevate TN: TP and PN: PP ratios as only the intercept term was included in the best GLMM models (Table 1; Figure 2). During the experiement, the TN: TP ratios were within the range 9–23 and PN: PP ratios were close to the therotical Redfield ratio 7 (Redfield, 1958) (Figure 2).

FIGURE 2

High turbidity was observed in both treatments, with TSS concentrations generally exceeding 20 mg L⁻¹ during the experiment. N addition significantly increased TSS, OSS, and ISS concentrations (Table 1; Figure 2).

3.2 Sediment P Release

The concentrations of Total-P and all P fractions except Residual-P at the end of the experiment (Day 70) were significantly lower than on Day 28 (Table 2; Figure 3). Moreover, the significant impact of the interactions of N addition and time on Total-P, NH₄Cl-P, NaOH-P, and HCl-P means that the magnitude of temporal declines was significantly greater in the N addition than in the control treatments (Table 2; Figure 3). N addition also significantly reduced BD-P, but its interaction with time was not significant (Table 2; Figure 3). No significance of N addition, time, and their interactions for Residual-P was observed in our study (Table 2; Figure 3).

FIGURE 3

3.3 Phytoplankton

N addition significantly boosted phytoplankton growth as both the Chl-a concentration (Table 1) and phytoplankton biomass (biovolume) (t-test, p = 0.008) were much higher in the N mesocosms than in the controls (Figures 4A,B). Meanwhile, we found a significant positive linear relation between phytoplankton Chl-a and biomass at the end of the experiment, which indicated that chlorophyll a could be used as an estimate of phytoplankton biomass in our study (Figure 4C). At the end of the experiment, large differences in phytoplankton communities between the two treatments were observed (Figure 4D). In the control mesocosms, the phytoplankton was dominated by diatoms (mainly Gyrosigma and Surirella), while filamentous cyanobacteria, e.g., Oscillatoria and Anabaena, dominated in the N treatment (Figure 4D). The differences in cyanobacteria and diatom dominance (%) between the two treatments were significant (t-test, p = 0.027 for cyanobacteria; p = 0.029 for diatoms; Figure 4D).

FIGURE 4

3.4 Zooplankton and Grazing Control

A total of 21 zooplankton taxa were recorded in the two treatments across the entire experimental period, including six cladoceran species (Ceriodaphnia quadrangular, Chydorus sphaericus, Moina sp., Bosmina sp. Alona sp., and Diaphanosoma sp.), two copepods (Mesocyclops leuckarti and Sinocalanus dorrii), and 13 rotifers (mainly Conochilus unicornis and Brachionus spp.). In our study, Mesocyclops leuckarti was the dominant species in both treatments. As the adult and most copepodid stages of this species can be carnivorous, we separately calculated total (all species included) and herbivorous zooplankton (with elimination of copepodites and adults of M. leuckarti) biomass in our study. Model selection by GLMMs revealed that no significant effect of N addition on cladoceran, copepod, and total and herbivorous zooplankton biomass occurred, but rotifer biomass was significantly higher in the N treatments than in the controls (Table 1; Figure 5).

FIGURE 5

No significant effects of N addition on zooplankton body mass and cladoceran to copepod biomass ratios were observed (Table 1; Figures 6A,B). However, N addition significantly reduced both total zooplankton and herbivorous zooplankton biomass to Chl-a ratios and increased the chlorophyll a to TP and TN ratios (Table 1; Figures 6C–F).

FIGURE 6

4. Discussion

Understanding the ecological mechanisms that underlie eutrophication in response to nutrient loading has long been an urgent objective in limnological research (Smith and Schindler, 2009). We found that summer N addition increased not only the N levels but also P availability for phytoplankton growth. We also found that N addition lowered the zooplankton-phytoplankton biomass ratio and promoted the phytoplankton yield, expressed as the chlorophyll a to TP and TN ratios, indicating that N addition favored phytoplankton growth through reduction of zooplankton grazing. Accordingly, our study suggests that in eutrophic shallow lakes, N loading may increase the phytoplankton biomass in summer via both resource and grazer effects. Our study thus showed that decreasing N loading may lead to reduced eutrophication in summer in shallow lakes, which is consistent with long-term field investigations conducted in temperate shallow waterbodies, e.g., Lake Müggelsee in Germany and shallow meres in England (Moss et al., 2005; Shatwell and Köhler, 2019).

As there was no external P loading in our study, the enhanced water TP concentrations induced by N loading can only be attributed to sediment P release, which was partially supported by the greater temporal reductions of sediment Total-P in the N addition treatments than in the control (Figure 3). Promotion by N addition of TP concentrations has been widely observed in eutrophic shallow waterbodies where sediments serve as an important P source (Jensen and Andersen, 1992; Xie et al., 2003; Sagrario et al., 2005; He et al., 2015; Ma et al., 2018). The mechanisms behind include increases in diffusion, bacteria activity, and pH and in the degradation of settled algae via the feedback between an N-induced high algal biomass and sediment P (Gao et al., 2014; Ma et al., 2018; 2021). Some of these mechanisms were also observed in our study. For instance, we found a clear decline in the contribution of TDP to TP after Day 42 and promotion of pH by N addition (Figure 2), indicating a high potential of diffusion (Xie et al., 2003; Chen et al., 2018; Ma et al., 2021) and desorption of PO₄^3- from suspended surface sediments (Koski-Vähälä and Hartikainen 2001; Søndergaard et al., 2003; Jin et al., 2006; Gao et al., 2014). This was also evidenced by the high reduction of NH₄Cl-P and NaOH-P caused by N addition in our study (Figure 2; Figure 3). Moreover, in our study, sediment P release may also have been enhanced by increased strength of sediment resuspension (Søndergaard et al., 2003; He et al., 2015) as the ISS concentrations in the N addition treatments were higher than in the control mesocosms (Figure 2). A possible explanation of the stronger resuspension by N loading, through speculative, is that fish had grown bigger in the N addition mesocosms due to the higher food production, as seen in the study by Wang et al. (2017). However, it is important to note that the reductions of Total-P in top 4 cm sediment were much larger than the increases of TP in the water. This means that only a small fraction of P in the sediment surface was released into the overlying water and most of the P was transported to the deeper sediment layers, likely due to the increased DO (Figure 2) and potential oxidized thickness by high NO₃⁻ addition (Jensen and Andersen, 1992). Overall, due to the simultaneous increase of water TP concentrations, persistent N addition did not significantly increase either the TN: TP or the PN: PP ratios, with phytoplankton growth likely being co-limited by N and P during the experimental period (Figure 2). The external N loading caused high internal P loading, thereby snowballing the eutrophication process via increased resources, as shown in Figure 1.

In our study, N addition caused a dominance replacement of diatoms by cyanobacteria (Figure 4D), as also seen in similar previous mesocosm experiments and lake investigations (Beaulieu et al., 2013; He et al., 2015; Gobler et al., 2016). The increased cyanobacteria dominance caused by N loading may be attributed to the increased competitive advantage by the shading-resistant cyanobacteria (e.g., Oscillatoria) in turbid and nutrient-enriched environments (Havens et al., 2003). Interestingly, N addition increased, not reduced, the dominance of Anabaena, an N₂-fixing taxon, in our study. This contrasted with the common view that declining N inputs would favor N₂-fixing cyanobacteria in lakes (Schindler et al., 2008), but concurred with observation from studies of a large number of Danish shallow lakes (Jensen et al., 1994; Jeppesen et al., 2005), suggesting that the mechanisms underlying the occurrence of N₂-fixing cyanobacteria in lakes are complex (Kolzau et al., 2018; Shatwell and Köhler, 2019). In our study, both TN: TP and PN: PP ratios were not obviously changed by N addition (Figure 2), which indicates that the growth of Anabaena in N addition mesocosms was unlikely caused by N deficiency.

Nutrient availability, however, was not the only reason why N addition increased the summer phytoplankton growth as it could not explain the observed elevation of the Chl-a: TP and Chl-a: TN ratios following N addition. The increased amount of chlorophyll a per unit of N and P suggests that factors other than resources affected the accumulation of phytoplankton biomass. Combined with the low total and herbivorous zooplankton to phytoplankton biomass ratios in the treatments with N addition, the results indicate that phytoplankton was also favored by the reduced herbivory. Similarly, Bergström et al. (2021) found much lower zooplankton to phytoplankton biomass ratios after whole-lake fertilization with N along a gradient of lake DOC concentrations. It could be argued that the increased Chl-a to TP and TN ratios by N addition were not due to changes in zooplankton herbivory, but instead were a result of changes of cellular pigment contents between phytoplankton groups (Nicholls and Dillon, 1978). However, this is not the case as we found a close relationship between phytoplankton biomass and Chl a concentration (Figure 4C).

In our study, the decline of zooplankton grazing by N addition could not be attributed to increased fish predation as both indicators, average zooplankton body mass and cladoceran to copepod biomass ratio, showed no obvious changes (Figure 6). Thus, changes in phytoplankton properties driven by N addition are the most likely explanation of the decreased zooplankton grazing. Supporting this view, an obvious N-driven substitution of diatoms by filamentous cyanobacteria was observed (Figure 4D). The latter are generally considered as low-quality food for both copepods and cladocerans (Caramujo et al., 2008; Martin-Creuzburg et al., 2008). For instance, laboratory feeding experiments revealed that when the freshwater harpacticoid copepod Attheyella trispinosa grew on cyanobacteria (Leptolyngbya foveolarum and Cylindrospermum stagnale), its survival was reduced, and the development from the first copepodite stage to adult lasted 26% longer than when they grew on diatoms (Achnanthes lanceolate and Nitzschia perminuta) (Caramujo et al., 2008). Moreover, cyanobacteria, e.g., Oscillatoria and Anabaena, may also resist zooplankton grazing because of poor manageability due to their long filaments and potential toxicity (DeMott and Moxter, 1991; Koski et al., 1999; Engström-Öst et al., 2015). Although feeding experiments have confirmed that Mesocyclops and other copepods may feed size-selectively on filaments and discriminate against toxic strains (DeMott and Moxter, 1991: Kumar and Ramakrishna Rao, 1999; Sarma et al., 2013), it is energetically costly and compensated for by a reduction in population size (DeMott and Moxter, 1991; Koski et al., 1999; Engström-öst et al., 2015).

In recent decades, increasing TN: TP ratios resulting from an imbalance between N and P supplies have been observed in many freshwater shallow lakes (Tong et al., 2020). However, the necessity of controlling N inputs remains a subject of debate (Schindler et al., 2008; Conley et al., 2009; Paerl et al., 2011; Shatwell and Köhler, 2019). Our results from a shallow eutrophic mesocosm study clearly showed that N inputs disproportionally increased the summer phytoplankton biomass (both Chl a and Chl a to TP or TN ratios) due to the combined effects of resource and grazing control. Our study therefore implies that N loading may not only physiologically support phytoplankton growth, it also modifies the abiotic and biotic environment and thereby creates conditions that are beneficial for phytoplankton growth. Thus, alleviation of eutrophication in shallow eutrophic lakes requires a balanced and strategic approach to appropriately control both nutrients (N and P).

References

• 1

  AndersenJ. M. (1982). Effect of Nitrate Concentration in Lake Water on Phosphate Release from the Sediment. Water Res.16, 1119–1126. 10.1016/0043-1354(82)90128-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  APHA (1998). Standard Methods for the Examination of Water and Waste Water. 20th ed.Washington, DC: American Public Health Association.

  • Google Scholar
  • View reference in article

• 3

  BeaulieuM.PickF.Gregory-EavesI. (2013). Nutrients and Water Temperature Are Significant Predictors of Cyanobacterial Biomass in a 1147 Lakes Data Set. Lim58, 1736–1746. 10.4319/lo.2013.58.5.1736

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BergströmA.-K.DeiningerA.JonssonA.KarlssonJ.VredeT. (2021). Effects of Nitrogen Enrichment on Zooplankton Biomass and N:P Recycling Ratios across a DOC Gradient in Northern-Latitude Lakes. Hydrobiologia848, 4991–5010. 10.1007/s10750-021-04689-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BolkerB. M.BrooksM. E.ClarkC. J.GeangeS. W.PoulsenJ. R.StevensM. H. H.et al (2009). Generalized Linear Mixed Models: A Practical Guide for Ecology and Evolution. Trends Ecol. Evol.24, 127–135. 10.1016/j.tree.2008.10.008

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BoströmB.PetterssonK. (1982). Different Patterns of Phosphorus Release from lake Sediments in Laboratory Experiments. Hydrobiologia91-92, 415–429. 10.1007/PL00020032

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BottrellH. C. (1976). A Review of Some Problems in Zooplankton Production Studies. Norw. J. Zool.24, 419–456.

  • Google Scholar
  • View reference in article

• 8

  BrooksM. E.KristensenK.BenthemK. J. v.MagnussonA.BergC. W.NielsenA.et al (2017). glmmTMB Balances Speed and Flexibility Among Packages for Zero-Inflated Generalized Linear Mixed Modeling. R. Journal9, 378–400. 10.32614/rj-2017-066

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  CaramujoM. J.BoschkerH. T.AdmiraalW. I. M. (2008). Fatty Acid Profiles of Algae Mark the Development and Composition of Harpacticoid Copepods. Freshw. Biol53, 77–90. 10.1111/j.1365-2427.2007.01868.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ChenM.DingS.ChenX.SunQ.FanX.LinJ.et al (2018). Mechanisms Driving Phosphorus Release during Algal Blooms Based on Hourly Changes in Iron and Phosphorus Concentrations in Sediments. Water Res.133, 153–164. 10.1016/j.watres.2018.01.040

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  ChiangS. C.DuN. S. (1979). Fauna Sinica; Crustacea: Freshwater Cladocera. Beijing: Science press.

  • Google Scholar
  • View reference in article

• 12

  ConleyD. J.PaerlH. W.HowarthR. W.BoeschD. F.SeitzingerS. P.HavensK. E.et al (2009). Controlling Eutrophication: Nitrogen and Phosphorus. Science323, 1014–1015. 10.1126/science.1167755

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DeMottW. R.MoxterF. (1991). Foraging Cyanobacteria by Copepods: Responses to Chemical Defense and Resource Abundance. Ecology72, 1820–1834. 10.2307/1940981

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  ElserJ. J.BrackenM. E. S.ClelandE. E.GrunerD. S.HarpoleW. S.HillebrandH.et al (2007). Global Analysis of Nitrogen and Phosphorus Limitation of Primary Producers in Freshwater, Marine and Terrestrial Ecosystems. Ecol. Lett.10 (12), 1135–1142. 10.1111/j.1461-0248.2007.01113.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  Engström-ÖstJ.BrutemarkA.VehmaaA.MotwaniN. H.KatajistoT. (2015). Consequences of a Cyanobacteria Bloom for Copepod Reproduction, Mortality and Sex Ratio. J. Plankton Res.37, 388–398. 10.1093/plankt/fbv004

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  GaoY.CornwellJ. C.StoeckerD. K.OwensM. S. (2014). Influence of Cyanobacteria Blooms on Sediment Biogeochemistry and Nutrient Fluxes. Limnol. Oceanogr.59, 959–971. 10.4319/lo.2014.59.3.0959

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  GoblerC. J.BurkholderJ. M.DavisT. W.HarkeM. J.JohengenT.StowC. A.et al (2016). The Dual Role of Nitrogen Supply in Controlling the Growth and Toxicity of Cyanobacterial Blooms. Harmful Algae54, 87–97. 10.1016/j.hal.2016.01.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  Gonzalez SagrarioM. A.JeppesenE.GomàJ.SondergaardM.JensenJ. P.LauridsenT.et al (2005). Does High Nitrogen Loading Prevent Clear-Water Conditions in Shallow Lakes at Moderately High Phosphorus Concentrations?Freshw. Biol50 (1), 27–41. 10.1111/j.1365-2427.2004.01290.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  GuildfordS. J.HeckyR. E. (2000). Total Nitrogen, Total Phosphorus, and Nutrient Limitation in Lakes and Oceans: Is There a Common Relationship?Limnol. Oceanogr.45, 1213–1223. 10.4319/lo.2000.45.6.1213

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  GuiryM. D.GuiryG. M. (2020). Algaebase. Galway: World- wide Electronic PublicationNational University of Ireland. Available at: http://www.algaebase.org (Accessed November 07, 2020).

  • Google Scholar
  • View reference in article

• 21

  GuoL. (2007). Doing battle with the green monster of Taihu Lake. Science.317, 1166. 10.1126/science.317.5842.1166

  • CrossRef
  • Google Scholar

• 22

  GyllströmM.HanssonL.-A.JeppesenE.CriadoF. G.GrossE.IrvineK.et al (2005). The Role of Climate in Shaping Zooplankton Communities of Shallow Lakes. Limnol. Oceanogr.50, 2008–2021. 10.4319/lo.2005.50.6.2008

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  HavensK. E.JamesR. T.EastT. L.SmithV. H. (2003). N:P Ratios, Light Limitation, and Cyanobacterial Dominance in a Subtropical Lake Impacted by Non-Point Source Nutrient Pollution. Environ. Pollut.122, 379–390. 10.1016/s0269-7491(02)00304-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  HeH.JinH.JeppesenE.LiK.LiuZ.ZhangY. (2018). Fish-Mediated Plankton Responses to Increased Temperature in Subtropical Aquatic Mesocosm Ecosystems: Implications for lake Management. Water Res.144, 304–311. 10.1016/j.watres.2018.07.055

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  HeH.KangY.LiuZ. (2015). Nitrogen Inputs Enhance Phytoplankton Growth during Sediment Resuspension Events: A Mesocosm Study. Hydrobiologia744, 297–305. 10.1007/s10750-014-2085-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  HeH.NingX.ChenK.LiQ.LiK.LiuZ.et al (2021). Intraguild Predation Dampens Trophic Cascades in Shallow Aquatic Mesocosms in the Subtropics: Implications for Lake Restoration by Biomanipulation. Freshw. Biol.66, 1571–1580. 10.1111/fwb.13739

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HuH. J.WeiY. X. (2006). The Freshwater Algae of China: Systematics, Taxonomy and Ecology. Beijing: Science Press. (in Chinese).

  • Google Scholar

• 28

  HupferM.GächterR.GiovanoliR. (1995). Transformation of Phosphorus Species in Settling Seston and during Early Sediment Diagenesis. Aquat. Sci.57, 305–324. 10.1007/BF00878395

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  JacksonL. J.LauridsenT. L.SøndergaardM.JeppesenE. (2007). A Comparison of Shallow Danish and Canadian Lakes and Implications of Climate Change. Freshw. Biol52, 1782–1792. 10.1111/j.1365-2427.2007.01809.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  JensenH. S.AndersenF. O. (1992). Importance of Temperature, Nitrate, and pH for Phosphate Release from Aerobic Sediments of Four Shallow, Eutrophic Lakes. Limnol. Oceanogr.37, 577–589. 10.4319/lo.1992.37.3.0577

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  JensenJ. P.JeppesenE.OlrikK.KristensenP. (1994). Impact of Nutrients and Physical Factors on the Shift from Cyanobacterial to Chlorophyte Dominance in Shallow Danish Lakes. Can. J. Fish. Aquat. Sci.51, 1692–1699. 10.1139/f94-170

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  JeppesenE.JensenJ. P.SondergaardM.LauridsenT. L. (2005). Response of Fish and Plankton to Nutrient Loading Reduction in Eight Shallow Danish Lakes with Special Emphasis on Seasonal Dynamics. Freshw. Biol50, 1616–1627. 10.1111/j.1365-2427.2005.01413.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  JeppesenE.NõgesP.DavidsonT. A.HabermanJ.NõgesT.BlankK.et al (2011). Zooplankton as Indicators in Lakes: a Scientific-Based Plea for Including Zooplankton in the Ecological Quality Assessment of Lakes According to the European Water Framework Directive (WFD). Hydrobiologia676, 279–297. 10.1007/s10750-011-0831-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  JinX.WangS.PangY.Chang WuF. (2006). Phosphorus Fractions and the Effect of pH on the Phosphorus Release of the Sediments from Different Trophic Areas in Taihu Lake, China. Environ. Pollut.139, 288–295. 10.1016/j.envpol.2005.05.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  KirkK. L.GilbertJ. J. (1992). Variation in Herbivore Response to Chemical Defenses: Zooplankton Foraging on Toxic Cyanobacteria. Ecology73, 2208–2217. 10.2307/1941468

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  KolzauS.DolmanA. M.VossM.WiednerC. (2018). The Response of Nitrogen Fixing Cyanobacteria to a Reduction in Nitrogen Loading. Internat Rev. Hydrobiol103, 5–14. 10.1002/iroh.201601882

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  KoskiM.EngströmJ.ViitasaloM. (1999). Reproduction and Survival of the Calanoid Copepod Eurytemora Affinis Fed with Toxic and Non-Toxic Cyanobacteria. Mar. Ecol. Prog. Ser.186, 187–197. 10.3354/meps186187

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  Koski-VähäläJ.HartikainenH. (2001). Assessment of the Risk of Phosphorus Loading due to Resuspended Sediment. J. Environ. Qual.30, 960–966. 10.2134/jeq2001.303960x

  • CrossRef
  • Google Scholar

• 39

  KumarR. A. M.Ramakrishna RaoT. (1999). Demographic Responses of Adult Mesocyclops Thermocyclopoides (Copepoda, Cyclopoida) to Different Plant and Animal Diets. Freshw. Biol42 (3), 487–501. 10.1046/j.1365-2427.1999.00485.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LiuZ.HuJ.ZhongP.ZhangX.NingJ.LarsenS. E.et al (2018). Successful Restoration of a Tropical Shallow Eutrophic Lake: strong Bottom-Up but Weak Top-Down Effects Recorded. Water Res.146, 88–97. 10.1016/j.watres.2018.09.007

  • CrossRef
  • Google Scholar

• 41

  LorenzP.TrommerG.StiborH. (2019). Impacts of Increasing Nitrogen:Phosphorus Ratios on Zooplankton Community Composition and whitefish ( Coregonus Macrophthalmus ) Growth in a Pre‐Alpine Lake. Freshw. Biol.64, 1210–1225. 10.1111/fwb.13296

  • CrossRef
  • Google Scholar

• 42

  MaS.-N.WangH.-J.WangH.-Z.ZhangM.LiY.BianS.-J.et al (2021). Effects of Nitrate on Phosphorus Release from Lake Sediments. Water Res.194, 116894. 10.1016/j.watres.2021.116894

  • CrossRef
  • Google Scholar

• 43

  MaS. N.WangH. J.WangH. Z.LiY.LiuM.LiangX. M.et al (2018). High Ammonium Loading Can Increase Alkaline Phosphatase Activity and Promote Sediment Phosphorus Release: A Two-Month Mesocosm experiment. Water Res.145, 388–397. 10.1016/j.watres.2018.08.043

  • CrossRef
  • Google Scholar

• 44

  MaoZ.GuX.CaoY.ZhangM.ZengQ.ChenH.et al (2020). The Role of Top-Down and Bottom-Up Control for Phytoplankton in a Subtropical Shallow Eutrophic Lake: Evidence Based on Long-Term Monitoring and Modeling. Ecosystems23, 1449–1463. 10.1007/s10021-020-00480-0

  • CrossRef
  • Google Scholar

• 45

  Martin-CreuzburgD.von ElertE.HoffmannK. H. (2008). Nutritional Constraints at the Cyanobacteria-Daphnia Magna Interface: The Role of Sterols. Limnol. Oceanogr.53, 456–468. 10.4319/lo.2008.53.2.0456

  • CrossRef
  • Google Scholar

• 46

  McCauleyE. (1984) “The Estimation of the Abundance and Biomass of Zooplankton in Samples,” in A Manual on Methods for the Assessment of Secondary Productivity in Fresh Waters. Editors DowningJ.RiglerF. H. (Oxford: Blackwell Scientific), 228–265.

  • Google Scholar

• 47

  MeisS.SpearsB. M.MaberlyS. C.O’MalleyM. B.PerkinsR. G. (2012). Sediment Amendment with Phoslock in Clatto Reservoir (Dundee, UK): Investigating Changes in Sediment Elemental Composition and Phosphorus Fractionation. J. Environ. Manage.93, 185–193. 10.1016/j.jenvman.2011.09.015

  • CrossRef
  • Google Scholar

• 48

  MossB.BarkerT.StephenD.WilliamsA. E.BalaylaD. J.BekliogluM.et al (2005). Consequences of Reduced Nutrient Loading on a Lake System in a Lowland Catchment: Deviations from the Norm?Freshw. Biol50, 1687–1705. 10.1111/j.1365-2427.2005.01416.x

  • CrossRef
  • Google Scholar

• 49

  NichollsK. H.DillonP. J. (1978). An Evaluation of Phosphorus-Chlorophyll-Phytoplankton Relationships for Lakes. Int. Revue Ges. Hydrobiol. Hydrogr.63 (2), 141–154. 10.1002/iroh.19780630203

  • CrossRef
  • Google Scholar

• 50

  PaerlH. W.XuH.McCarthyM. J.ZhuG.QinB.LiY.et al (2011). Controlling Harmful Cyanobacterial Blooms in a Hyper-Eutrophic lake (Lake Taihu, China): The Need for a Dual Nutrient (N & P) Management Strategy. Water Res.45, 1973–1983. 10.1016/j.watres.2010.09.018

  • CrossRef
  • Google Scholar

• 51

  PeñuelasJ.PoulterB.SardansJ.CiaisP.Van Der VeldeM.BoppL.et al (2013). Human-Induced Nitrogen-Phosphorus Imbalances Alter Natural and Managed Ecosystems across the globe. Nat. Commun.4, 1–10. 10.1038/ncomms3934

  • CrossRef
  • Google Scholar

• 52

  PeñuelasJ.SardansJ.Rivas-ubachA.JanssensI. A. (2012). The Human-Induced Imbalance between C, N and P in Earth's Life System. Glob. Change Biol.18, 3–6. 10.1111/j.1365-2486.2011.02568.x

  • CrossRef
  • Google Scholar

• 53

  QinB.PaerlH. W.BrookesJ. D.LiuJ.JeppesenE.ZhuG.et al (2019). Why Lake Taihu Continues to be Plagued with Cyanobacterial Blooms Through 10 years (2007–2017) Efforts. Sci. Bull.64, 354–356. 10.1016/j.scib.2019.02.008

  • CrossRef
  • Google Scholar

• 54

  R Core Team(2020). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available at: https://www.R-project.org/ (Accessed September 25, 2020).

  • Google Scholar

• 55

  RedfieldA. C. (1958). The Biological Control of Chemical Factors in the Environment. Am. Sci.46, 230A–221A.

  • Google Scholar

• 56

  RubanV.López-SánchezJ. F.PardoP.RauretG.MuntauH.QuevauvillerP. (1999). Selection and Evaluation of Sequential Extraction Procedures for the Determination of Phosphorus Forms in Lake Sediment. J. Environ. Monit.1, 51–56. 10.1039/A807778I

  • CrossRef
  • Google Scholar

• 57

  RydinE. (2000). Potentially Mobile Phosphorus in Lake Erken Sediment. Water Res.34, 2037–2042. 10.1016/S0043-1354(99)00375-9

  • CrossRef
  • Google Scholar

• 58

  SarmaS. S. S.Jiménez-ContrerasJ.FernándezR.NandiniS.García-GarcíaG. (2013). Functional Responses and Feeding Rates ofMesocyclops pehpeiensisHu (Copepoda) Fed Different Diets (Rotifers, Cladocerans, Alga and Cyanobacteria). J. Nat. Hist.47, 841–852. 10.1080/00222933.2012.747636

  • CrossRef
  • Google Scholar

• 59

  SchindlerD. W.HeckyR. E.FindlayD. L.StaintonM. P.ParkerB. R.PatersonM. J.et al (2008). Eutrophication of Lakes Cannot Be Controlled by Reducing Nitrogen Input: Results of a 37-year Whole-Ecosystem experiment. Proc. Natl. Acad. Sci.105, 11254–11258. 10.1073/pnas.0805108105

  • CrossRef
  • Google Scholar

• 60

  SEPA (2002). Analytical Methods for Water and Wastewater Monitor. 4th edn. Beijing: Chinese Environmental Science Press. (in Chinese).

  • Google Scholar

• 61

  ShatwellT.KöhlerJ. (2019). Decreased Nitrogen Loading Controls Summer Cyanobacterial Blooms without Promoting Nitrogen‐Fixing Taxa: Long‐Term Response of a Shallow lake. Limnol. Oceanogr64, S166–S178. 10.1002/lno.11002

  • CrossRef
  • Google Scholar

• 62

  ShenJ. R.DuN. S. (1979). Fauna Sinica, Crustacea, Freshwater Copepoda. Beijing: Science Press. (in Chinese).

  • Google Scholar

• 63

  SmithV. H.SchindlerD. W. (2009). Eutrophication Science: Where Do We Go from Here?Trends Ecol. Evol.24, 201–207. 10.1016/j.tree.2008.11.009

  • CrossRef
  • Google Scholar

• 64

  SøndergaardM.BjerringR.JeppesenE. (2013). Persistent Internal Phosphorus Loading during Summer in Shallow Eutrophic Lakes. Hydrobiologia710, 95–107. 10.1007/s10750-012-1091-3

  • CrossRef
  • Google Scholar

• 65

  SøndergaardM.JensenJ. P.JeppesenE. (2003). Role of Sediment and Internal Loading of Phosphorus in Shallow Lakes. Hydrobiologia506-509, 135–145. 10.1023/B:HYDR.0000008611.12704.dd

  • CrossRef
  • Google Scholar

• 66

  SøndergaardM.LauridsenT. L.JohanssonL. S.JeppesenE. (2017). Nitrogen or Phosphorus Limitation in Lakes and its Impact on Phytoplankton Biomass and Submerged Macrophyte Cover. Hydrobiologia795, 35–48. 10.1007/s10750-017-3110-x

  • CrossRef
  • Google Scholar

• 67

  TongY.WangM.PeñuelasJ.LiuX.PaerlH. W.ElserJ. J.et al (2020). Improvement in Municipal Wastewater Treatment Alters lake Nitrogen to Phosphorus Ratios in Populated Regions. Proc. Natl. Acad. Sci. USA117, 11566–11572. 10.1073/pnas.1920759117

  • CrossRef
  • Google Scholar

• 68

  TrommerG.PoxleitnerM.LorenzP.BitzilekisE.GogaladzeA.SchultesS.et al (2017). Altered Food-Web Dynamics under Increased Nitrogen Load in Phosphorus Deficient Lakes. Aquat. Sci.79, 1009–1021. 10.1007/s00027-017-0551-2

  • CrossRef
  • Google Scholar

• 69

  WangH.-J.XiaoX.-C.WangH.-Z.LiY.YuQ.LiangX.-M.et al (2017). Effects of High Ammonia Concentrations on Three Cyprinid Fish: Acute and Whole-Ecosystem Chronic Tests. Sci. Total Environ.598, 900–909. 10.1016/j.scitotenv.2017.04.070

  • CrossRef
  • Google Scholar

• 70

  XieL.XieP.LiS.TangH.LiuH. (2003). The Low TN:TP Ratio, a Cause or a Result of Microcystis Blooms?Water Res.37, 2073–2080. 10.1016/S0043-1354(02)00532-8

  • CrossRef
  • Google Scholar

• 71

  XuH.McCarthyM. J.PaerlH. W.BrookesJ. D.ZhuG.HallN. S.et al (2021). Contributions of External Nutrient Loading and Internal Cycling to Cyanobacterial Bloom Dynamics in Lake Taihu, China: Implications for Nutrient Management. Limnol. Oceanogr.66, 1492–1509. 10.1002/lno.11700

  • CrossRef
  • Google Scholar

• 72

  XuH.PaerlH. W.QinB.ZhuG.GaoaG. (2010). Nitrogen and Phosphorus Inputs Control Phytoplankton Growth in Eutrophic Lake Taihu, China. Limnol. Oceanogr.55, 420–432. 10.4319/lo.2010.55.1.0420

  • CrossRef
  • Google Scholar

• 73

  ZhangZ. S.HuangX. F. (1991). Methods in Freshwater Plankton Study. Beijing: Science Press. (in Chinese).

  • Google Scholar