Environmental drivers of phytoplankton community dynamics in an agriculturally-influenced tributary in the lower Great Lakes

Phytoplankton community composition in tributaries differs from that in their receiving waters, due to light limitation from suspended particles and other factors such as nutrient availability and temperature. This study was designed to manipulate light levels in early, mid, and late summer to determine the combined effects of light attenuation and naturally varying nutrient availability on phytoplankton community composition in an agriculturally-influenced tributary of the lower Great Lakes. In all trials, in situ microcosm experiments show that phytoplankton abundance increased under three light attenuation treatments (60 %, 75 %, and 85 % attenuation) relative to time-zero, but higher light attenuation reduced total phytoplankton abundance relative to controls. Highest phytoplankton diversity in terms of richness and evenness occurred in September (late summer), and across all three trials was lowest under the highest light attenuation treatments (85 %). Phytoplankton community composition followed a normal seasonal shift from diatoms dominating in June (early summer), followed by cyanobacteria dominating in mid to late summer. In general, lower light levels (especially 85 % attenuation) corresponded with an increased dominance of cyanobacteria. These findings support the hypothesis that phytoplankton abundance and diversity vary with light and nutrient availability and that light attenuation promotes the shift from buoyant cyanobacteria to other taxa more tolerant of low light levels.

Aerobic anoxygenic phototrophs play important roles in nutrient cycling within cyanobacterial Microcystis bloom microbiomes

BACKGROUND

During the bloom season, the colonial cyanobacterium Microcystis forms complex aggregates which include a diverse microbiome within an exopolymer matrix. Early research postulated a simple mutualism existing with bacteria benefitting from the rich source of fixed carbon and Microcystis receiving recycled nutrients. Researchers have since hypothesized that Microcystis aggregates represent a community of synergistic and interacting species, an interactome, each with unique metabolic capabilities that are critical to the growth, maintenance, and demise of Microcystis blooms. Research has also shown that aggregate-associated bacteria are taxonomically different from free-living bacteria in the surrounding water. Moreover, research has identified little overlap in functional potential between Microcystis and members of its microbiome, further supporting the interactome concept. However, we still lack verification of general interaction and know little about the taxa and metabolic pathways supporting nutrient and metabolite cycling within Microcystis aggregates.

RESULTS

During a 7-month study of bacterial communities comparing free-living and aggregate-associated bacteria in Lake Taihu, China, we found that aerobic anoxygenic phototrophic (AAP) bacteria were significantly more abundant within Microcystis aggregates than in free-living samples, suggesting a possible functional role for AAP bacteria in overall aggregate community function. We then analyzed gene composition in 102 high-quality metagenome-assembled genomes (MAGs) of bloom-microbiome bacteria from 10 lakes spanning four continents, compared with 12 complete Microcystis genomes which revealed that microbiome bacteria and Microcystis possessed complementary biochemical pathways that could serve in C, N, S, and P cycling. Mapping published transcripts from Microcystis blooms onto a comprehensive AAP and non-AAP bacteria MAG database (226 MAGs) indicated that observed high levels of expression of genes involved in nutrient cycling pathways were in AAP bacteria.

CONCLUSIONS

Our results provide strong corroboration of the hypothesized Microcystis interactome and the first evidence that AAP bacteria may play an important role in nutrient cycling within Microcystis aggregate microbiomes. Video Abstract.

Microbiome processing of organic nitrogen input supports growth and cyanotoxin production of Microcystis aeruginosa cultures

Nutrient-induced blooms of the globally abundant freshwater toxic cyanobacterium Microcystis cause worldwide public and ecosystem health concerns. The response of Microcystis growth and toxin production to new and recycled nitrogen (N) inputs and the impact of heterotrophic bacteria in the Microcystis phycosphere on these processes are not well understood. Here, using microbiome transplant experiments, cyanotoxin analysis, and nanometer-scale stable isotope probing to measure N incorporation and exchange at single cell resolution, we monitored the growth, cyanotoxin production, and microbiome community structure of several Microcystis strains grown on amino acids or proteins as the sole N source. We demonstrate that the type of organic N available shaped the microbial community associated with Microcystis, and external organic N input led to decreased bacterial colonization of Microcystis colonies. Our data also suggest that certain Microcystis strains could directly uptake amino acids, but with lower rates than heterotrophic bacteria. Toxin analysis showed that biomass-specific microcystin production was not impacted by N source (i.e. nitrate, amino acids, or protein) but rather by total N availability. Single-cell isotope incorporation revealed that some bacterial communities competed with Microcystis for organic N, but other communities promoted increased N uptake by Microcystis, likely through ammonification or organic N modification. Our laboratory culture data suggest that organic N input could support Microcystis blooms and toxin production in nature, and Microcystis-associated microbial communities likely play critical roles in this process by influencing cyanobacterial succession through either decreasing (via competition) or increasing (via biotransformation) N availability, especially under inorganic N scarcity.

Effective removal of nitrate in water by continuous-flow electro-dialysis ion exchange membrane bioreactor (CF-EDIMB): Performance optimization and microbial analysis

The use of nitrogen fertilizer has been causing nitrate pollution in groundwater, and there is an urgent need for efficient approach to remove nitrate from groundwater. In our job, a novel continuous-flow electrodialysis ion exchange membrane bioreactor system (CF-EDIMB) was set up to remove nitrate (NO3-) from water for the first time. Nitrate removal was positively dependent on water chamber HRT and voltage; voltage had significant effect on the water chamber effluent pH; acetate utilization efficiency was closely correlated with acetate dosage. The optimal conditions forecasted through response surface method (RSM) were given as follows: water chamber HRT was 20 h, biological chamber HRT was 24 h, voltage was 6.65 V and acetate dosage was 454.99 mg/L, dedicating to nitrate removal of 81.90% (83.70% in prediction), water chamber effluent pH of 7.10 (7.00 in prediction) and acetate utilization efficiency of 92.87% (96.51% in prediction). Meanwhile, microorganisms are crucial for nitrate removal, and the microbial community was not sensitive to the variation of acetate dosage. The microbial analysis results indicated that when CF-EDIMB system was operated for 20 d, the sulfate-reducing bacteria Sediminibacterium appeared in the biological chamber, and the effluent sulfate concentration of biological chamber was decreased. During the whole operation, Thauera was the dominant genus. Denitrifying functional genes nirS presented a better expression than the gene narG, and there was no accumulation of nitrite.

Transcriptionally active nitrogen fixation and biosynthesis of diverse secondary metabolites by Dolichospermum and Aphanizomenon-like Cyanobacteria in western Lake Erie Microcystis blooms

Cyanobacterial harmful algal blooms (cyanoHABs) in the western basin of Lake Erie are dominated by microcystin producing Microcystis spp., but other cyanobacterial taxa that coexist in these communities may play important roles in production of toxins and shaping bloom dynamics and community function. In this study, we used metagenomic and metatranscriptomic data from the 2014 western Lake Erie cyanoHAB to explore the genetic diversity and biosynthetic potential of cyanobacteria belonging to the Anabaena, Dolichospermum, Aphanizomenon (ADA) clade. We reconstructed two near-complete metagenome-assembled genomes from two distinct ADA clade species, each containing biosynthetic gene clusters that encode novel and known secondary metabolites, including those with toxic and/or known taste and odor properties, that were transcriptionally active. However, neither ADA metagenome-assembled genome contained genes encoding guanitoxins, anatoxins, or saxitoxins, which are known to be produced by ADA. The ADA cyanobacteria accounted for most of the metagenomic and metatranscriptomic reads from nitrogen fixation genes, suggesting they were the dominant N-fixers at the times and stations sampled. Despite their relatively low abundance, our results highlight the possibility that ADA taxa could influence the water quality and ecology of Microcystis blooms, although the extent of these impacts remains to be quantified.

Short-term Lake Erie algal bloom prediction by classification and regression models

The recent outbreaks of harmful algal blooms in the western Lake Erie Basin (WLEB) have drawn tremendous attention to bloom prediction for better control and management. Many weekly to annual bloom prediction models have been reported, but they only employ small datasets, have limited types of input features, build linear regression or probabilistic models, or require complex process-based computations. To address these limitations, we conducted a comprehensive literature review, complied a large dataset containing chlorophyll-a index (from 2002 to 2019) as the output and a novel combination of riverine (the Maumee & Detroit Rivers) and meteorological (WLEB) features as the input, and built machine learning-based classification and regression models for 10-d scale bloom predictions. By analyzing the feature importance, we identified 8 most important features for the HAB control, including nitrogen loads, time, water levels, soluble reactive phosphorus load, and solar irradiance. Here, both long- and short-term nitrogen loads were for the first time considered in HAB models for Lake Erie. Based on these features, the 2-, 3-, and 4-level random forest classification models achieved an accuracy of 89.6%, 77.0%, and 66.7%, respectively, and the regression model achieved an R² value of 0.69. In addition, long-short term memory (LSTM) was implemented to predict temporal trends of four short-term features (N, solar irradiance, and two water levels) and achieved the Nash-Sutcliffe efficiency of 0.12-0.97. Feeding the LSTM model predictions for these features into the 2-level classification model reached an accuracy of 86.0% for predicting the HABs in 2017-2018, suggesting that we can provide short-term HAB forecasts even when the feature values are not available.

The effect of single versus dual nutrient decreases on phytoplankton growth rates, community composition, and Microcystin concentration in the western basin of Lake Erie

The primary management strategy for minimizing harmful algal blooms (HABs) in Lake Erie has been to reduce springtime loading of phosphorus (P) to the lake. However, some studies have shown that the growth rate and toxin content for the HABs-causing cyanobacterium Microcystis also respond to the availability of dissolved inorganic nitrogen (N). This evidence is based on both observational studies that correlate bloom development with changes in N forms and concentrations in the lake, and experiments in which P and/or N are added at concentrations in excess of those present in the lake. The goal of this study was to determine whether a combined decrease in N and P concentrations from ambient levels in Lake Erie could limit the development of HABs more than a reduction in P concentration only. To directly test the impact of P-only versus dual N and P concentration decreases on phytoplankton in the western basin of Lake Erie, we evaluated changes in growth rate, community composition, and microcystin (MC) concentration through eight bioassay experiments performed from June through October 2018, which encompassed the normal Lake Erie Microcystis-dominated HAB season. Our results showed that during the first five experiments covering June 25 to August 13, the P-only and the dual N and P decrease treatments had similar effects. However, when ambient N became scarce later in the season, the N and P decrease treatments resulted in negative growth rates for cyanobacteria, whereas -P only decreases did not. During low ambient N conditions, dual nutrient decreases lowered the prevalence of cyanobacteria among the total phytoplankton community and decreased microcystin concentrations. The results presented here complement previous experimental work on Lake Erie and suggest that dual nutrient control could be an effective management strategy to decrease microcystin production during the bloom and even possibly diminish or shorten the duration of the bloom based on creating nutrient limiting conditions sooner in the HAB growing season.

Nitrogen Limitation of Intense and Toxic Cyanobacteria Blooms in Lakes within Two of the Most Visited Parks in the USA: The Lake in Central Park and Prospect Park Lake

The Lake in Central Park (LCP) and Prospect Park Lake (PPL) in New York City (NYC), USA, are lakes within two of the most visited parks in the USA. Five years of nearshore sampling of these systems revealed extremely elevated levels of cyanobacteria and the toxin, microcystin, with microcystin levels averaging 920 µg L−1 and chlorophyll a from cyanobacterial (cyano-chla) populations averaging 1.0 × 105 µg cyano-chla L−1. Both lakes displayed elevated levels of orthophosphate (DIP) relative to dissolved inorganic nitrogen (DIN) during summer months when DIN:DIP ratios were < 1. Nutrient addition and dilution experiments revealed that N consistently limited cyanobacterial populations but that green algae were rarely nutrient limited. Experimental additions of public drinking water that is rich in P and, to a lesser extent N, to lake water significantly enhanced cyanobacterial growth rates in experiments during which N additions also yielded growth enhancement. Collectively, this study demonstrates that the extreme microcystin levels during blooms in these highly trafficked lakes represent a potential human and animal health threat and that supplementation of these artificial lakes with public drinking water to maintain water levels during summer may promote the intensity and N limitation of blooms.

The crucial influence of trophic status on the relative requirement of nitrogen to phosphorus for phytoplankton growth

Clarifying the pattern of relative nitrogen (N)-to-phosphorus (P) requirements for phytoplankton growth is of great significance for eutrophication mitigation and aquatic system management. The relative N-to-P requirement for phytoplankton growth is considered an essential trait determining species dominance within ecosystems and explaining phytoplankton response to nutrient availability. These requirements vary with environmental trophic status, though this variation remains unclear. Here, we evaluated the relative N-to-P requirements under different absolute nutrient levels using previous and current experimental data on eight phytoplankton species (three studied by us and five extrapolated from previous studies). Results showed that relative N-to-P requirements for phytoplankton growth decreased as absolute nutrient levels increased. Thus, N may be crucial for enhancing phytoplankton growth under low nutrient conditions, whereas P may be the primary limiting factor of phytoplankton growth under sufficient nutrient conditions. This result applies to single species as well as species assemblages, which are independent of species shifts occurring along water N:P gradients. The response observed in our large trophic status gradient may help elucidate the relative importance of N and P reductions in mitigating the impact of eutrophication on ecosystems.

Quantification of microcystin production and biodegradation rates in the western basin of Lake Erie

Cyanobacterial biomass forecasts currently cannot predict the concentrations of microcystin, one of the most ubiquitous cyanotoxins that threaten human and wildlife health globally. Mechanistic insights into how microcystin production and biodegradation by heterotrophic bacteria change spatially and throughout the bloom season can aid in toxin concentration forecasts. We quantified microcystin production and biodegradation during two growth seasons in two western Lake Erie sites with different physicochemical properties commonly plagued by summer Microcystis blooms. Microcystin production rates were greater with elevated nutrients than under ambient conditions and were highest nearshore during the initial phases of the bloom, and production rates were lower in later bloom phases. We examined biodegradation rates of the most common and toxic microcystin by adding extracellular stable isotope-labeled microcystin-LR (1 μg L-1), which remained stable in the abiotic treatment (without bacteria) with minimal adsorption onto sediment, but strongly decreased in all unaltered biotic treatments, suggesting biodegradation. Greatest biodegradation rates (highest of -8.76 d-1, equivalent to the removal of 99.98% in 18 h) were observed during peak bloom conditions, while lower rates were observed with lower cyanobacteria biomass. Cell-specific nitrogen incorporation from microcystin-LR by nanoscale imaging mass spectrometry showed that a small percentage of the heterotrophic bacterial community actively degraded microcystin-LR. Microcystin production and biodegradation rates, combined with the microcystin incorporation by single cells, suggest that microcystin predictive models could be improved by incorporating toxin production and biodegradation rates, which are influenced by cyanobacterial bloom stage (early vs. late bloom), nutrient availability, and bacterial community composition.

Quantifying the contribution of direct runoff and baseflow to nitrogen loading in the Western Lake Erie Basins

Soluble nitrogen is highly mobile in soil and susceptible to leaching. It is important to identify nitrogen transport pathways so that the sources can be efficiently targeted in environment management. This study quantified the contribution of direct runoff and baseflow to nitrate + nitrite loading by separating flow and nitrate + nitrite concentration measurements into two periods depending on whether only baseflow was present or not using baseflow separation methods. When both direct runoff and baseflow were present in streamflow, their nitrate + nitrite concentrations were assumed based on the hydrological reasoning that baseflow does not change rapidly, and streamflow mostly consists of direct runoff within a rainfall event. For this study, we obtained and investigated daily flow and nitrate + nitrite concentration observations made at the outlets of 22 watersheds located in the Western Lake Erie area. Results showed that baseflow was responsible for 26 to 77% of the nitrate + nitrite loads. The relative nitrate + nitrite load contributions of direct runoff and baseflow substantially varied with the sizes of drainage areas and agricultural land uses. Increases in drainage areas tend to prolong the travel time of surface runoff and thus help its reinfiltration into soil, which then could increase the baseflow contribution. In addition, the artificial drainage networks common in the agricultural fields of the study areas would promote the drainage of nutrient-laden excess water from soils. Such findings suggest the need for environmental management customized considering nitrogen transport pathways.

Dynamic Phycobilin Pigment Variations in Diazotrophic and Non-diazotrophic Cyanobacteria Batch Cultures Under Different Initial Nitrogen Concentrations

Increased anthropogenic nutrient loading has led to eutrophication of aquatic ecosystems, which is the major cause of harmful cyanobacteria blooms. Element stoichiometry of cyanobacteria bloom is subject to nutrient availabilities and may significantly contribute to primary production and biogeochemical cycling. Phycobilisome is the antenna of the photosynthetic pigment apparatus in cyanobacteria, which contains phycobilin pigments (PBPs) and linker proteins. This nitrogen (N)-rich protein complex has the potential to support growth as a N-storage site and may play a major role in the variability of cyanobacteria N stoichiometry. However, the regulation of PBPs during bloom formation remains unclear. We investigated the temporal variation of N allocation into PBPs and element stoichiometry for two ubiquitous cyanobacteria species, Microcystis aeruginosa and Dolichospermum flos-aquae, in a batch culture experiment with different initial N availabilities. Our results indicated that the N allocation into PBPs is species-dependent and tightly regulated by the availability of nutrients fueling population expansion. During the batch culture experiment, different nutrient uptake rates led to distinct stoichiometric imbalances of N and phosphorus (P), which substantially altered cyanobacteria C: N and C: P stoichiometry. Microcystis invested cellular N into PBPs and exhibited greater flexibility in C: N and C: P stoichiometry than D. flos-aquae. The dynamics of such N-rich macromolecules may help explain the N stoichiometry variation during a bloom and the interspecific difference between M. aeruginosa and D. flos-aquae. Our study provides a quantitative understanding of the elemental stoichiometry and the regulation of PBPs for non-diazotrophic and diazotrophic cyanobacteria blooms.

Simplified Method for the Determination of Total Kjeldahl Nitrogen in Wastewater

The increase in total Kjeldahl nitrogen (TKN) concentrations is correlated with increases in cyanobacterial bloom biomass. Standard methods for the measurement of TKN are tedious, costly, time-consuming and involve the use of hazardous catalysts, such as mercury, high temperatures and significant amounts of toxic acids and bases. Since TKN plays a pivotal role in influencing algal blooms, there is an urgent need to develop simpler, safer and more accurate methods for the determination of TKN. The simplified TKN method (s-TKN™) developed by Hach® offers several advantages over the traditional TKN method, including eliminating the use of mercury, requiring low sample and reagent volumes and being cost-efficient and user-friendly. This communication presents preliminary results comparing the efficacy of s-TKN™ and the standard method, using commonly used primary standards and waste, estuarine and lake water matrices. For all primary standards analyzed, the s-TKN™ method exhibited good accuracy across a wide range of concentrations. The repeatability for the glycine–para-toluene sulfonic acid (Gly-PTSA) standard using the s-TKN™ method was 4.1% at the highest concentrations analyzed, with overall repeatability across concentrations comparable to the standard EPA method. For wastewater, estuarine and lake matrices, a good correlation (r2 = 0.9917) between the two methods and no statistical difference in the values (p > 0.05) were obtained between two methods. Preliminary studies indicate that the s-TKN™ method has the potential to reduce the expenditure associated with the cost of analysis and has the potential to be a safer and cheaper alternative, while providing comparable analytical results to the standard method.

Modeling the sensitivity of cyanobacteria blooms to plausible changes in precipitation and air temperature variability

Many recent studies have attributed the observed variability of cyanobacteria blooms to meteorological drivers and have projected blooms with worsening societal and ecological impacts under future climate scenarios. Nonetheless, few studies have jointly examined their sensitivity to projected changes in both precipitation and temperature variability. Using an Integrated Assessment Model (IAM) of Lake Champlain's eutrophic Missisquoi Bay, we demonstrate a factorial design approach for evaluating the sensitivity of concentrations of chlorophyll a (chl-a), a cyanobacteria surrogate, to global climate model-informed changes in the central tendency and variability of daily precipitation and air temperature. An Analysis of Variance (ANOVA) and multivariate contour plots highlight synergistic effects of these climatic changes on exceedances of the World Health Organization's moderate 50 μg/L concentration threshold for recreational contact. Although increased precipitation produces greater riverine total phosphorus loads, warmer and drier scenarios produce the most severe blooms due to the greater mobilization and cyanobacteria uptake of legacy phosphorus under these conditions. Increases in daily precipitation variability aggravate blooms most under warmer and wetter scenarios. Greater temperature variability raises exceedances under current air temperatures but reduces them under more severe warming when water temperatures exceed optimal values for cyanobacteria growth more often. Our experiments, controlled for wind-induced changes to lake water quality, signal the importance of larger summer runoff events for curtailing bloom growth through reductions of water temperature, sunlight penetration and stratification. Finally, the importance of sequences of wet and dry periods in generating cyanobacteria blooms motivates future research on bloom responses to changes in interannual climate persistence.

Evaluation of the Potential Release Risk of Internal N and P from Sediments—A Preliminary Study in Two Freshwater Reservoirs in South China

Growing evidence has demonstrated the influence of internal nitrogen (N) and phosphorus (P) on harmful algae blooms in eutrophic freshwater ecosystems. However, the main controlling factors for internal N and P release risks, and whether these factors vary as environmental conditions change, remains poorly understood. We evaluated potential release risks of N and P from sediments in two freshwater reservoirs in Beihai City, southern China, by evaluating apparent nutrient fluxes during simulated static incubation experiments at two temperatures (15 °C and 25 °C). Sediments were analyzed to determine their basic properties as well as N and P fractions. Results showed that the main controlling factors of the apparent fluxes in dissolved total P, soluble reactive P, total N, and ammonium were related to sediment adsorption properties, redox properties, and microbial-mediated properties (e.g., water-extractable P, total inorganic N, redox-sensitive P, total organic carbon, organic P). The primary controlling factors for apparent N and P fluxes were dependent on the form of N and P and changed with temperature. The results suggest that care should be taken when simply using total N and P contents in sediments to evaluate their internal nutrient release risks.

Dynamic Responses of Endosymbiotic Microbial Communities Within Microcystis Colonies in North American Lakes to Altered Nitrogen, Phosphorus, and Temperature Levels

The toxic cyanobacterium, Microcystis, is a pervasive cyanobacterial harmful algal bloom (CHAB) - forming genus that naturally occurs in colonies that harbor diverse microbiomes of heterotrophic bacteria. While the effects of nutrient loading and climatic warming on CHABs are well-known, little is known regarding how these environmental drivers alter the structural and functional potential of the microbial assemblages associated with blooms that, in turn, may impact cyanobacterial growth. Here, we used next-generation sequencing of 16S ribosomal rRNA genes to characterize the dynamics of the bacterial assemblages within Microcystis colonies in two temperate North American lakes: Lake Erie and Lake Agawam (NY, United States) and quantified their responses to experimentally increased levels of nitrogen (N), phosphorus (P) and temperature. Across experiments, Microcystis populations were consistently and significantly promoted by N and, to a lesser extent, elevated temperature (p < 0.05). In contrast, bacterial assemblages within Microcystis colonies were more resilient to environmental perturbations, with the relative abundance of 7–16% of amplicon sequence variants changing and several individual taxa displaying significant (p < 0.05) increases and decreases in relative abundance, primarily in response to elevated temperature and to a lesser extent, N. In contrast to individual taxa, community diversity was not significantly altered by individual treatments during experiments but rather was inversely correlated with the intensity of Microcystis blooms (p < 0.001). While predicted metabolic function was even less impacted by environmental drivers than microbial diversity, the predicted abundance of nitrogenase (nifH), alkaline phosphatase (phoX), and urease (ure) genes significantly increased in response to N but decreased in response to increased temperature (p < 0.05). Collectively, the resilience of microbial community structure and function within colonies suggests they may support the ability of Microcystis to persist through short-term fluctuations in environmental conditions by supplying essential nutrients.

Novel high throughput sequencing - fluorometric approach demonstrates Microcystis blooms across western Lake Erie are promoted by grazing resistance and nutrient enhanced growth

Cyanobacterial harmful algal blooms (CHABs) are a global public health threat. While CHABs are often promoted by nutrients, an important and often overlooked influence on bloom dynamics is zooplankton grazing. In the present study, zooplankton grazing and nutrient enrichment experiments were combined with next generation sequencing and fluorometric analyses to quantify differential grazing and nutrient effects on specific cyanobacterial genera across the western basin of Lake Erie. Grazing by two different sized daphnids, Daphnia magna and Daphnia pulex, was compared to protozooplankton grazing effects assessed via a dilution approach at sites within the Maumee and Sandusky Bays where Planktothrix, Microcystis, Synechococcus, and Dolichospermum were the dominant genera. Daphnid grazing significantly reduced Synechococcus net growth rates at most sites as well as Planktothrix net growth in Sandusky Bay and Dolichospermum in Maumee Bay. Dilution resulted in significant growth increase of Synechococcus at half of the sites and Planktothrix at most sites evidencing substantial grazing pressure by the protozooplankton community on these genera. In contrast, Microcystis populations were largely unaffected by daphnids and protozooplankton grazing but benefitted from nutrient enrichment more than other CHAB genera. When diatoms were present in moderate abundance, grazing rates by daphnids on diatoms were significantly greater than grazing rates on cyanobacteria. The novel approach used in this study established differences in grazing pressure and nutrient effects on differing taxa and revealed that, while many taxa were grazed by multiple classes of zooplankton (e.g. Planktothrix, Synechococcus, Dolichospermum, diatoms), the lack of grazing pressure on Microcystis coupled with nutrient-enhanced growth in western Lake Erie promotes the occurrence of CHABs of this genus.

The Lake Erie HABs Grab: A binational collaboration to characterize the western basin cyanobacterial harmful algal blooms at an unprecedented high-resolution spatial scale

Monitoring of cyanobacterial bloom biomass in large lakes at high resolution is made possible by remote sensing. However, monitoring cyanobacterial toxins is only feasible with grab samples, which, with only sporadic sampling, results in uncertainties in the spatial distribution of toxins. To address this issue, we conducted two intensive "HABs Grabs" of microcystin (MC)-producing Microcystis blooms in the western basin of Lake Erie. These were one-day sampling events during August of 2018 and 2019 in which 100 and 172 grab samples were collected, respectively, within a six-hour window covering up to 2,270 km² and analyzed using consistent methods to estimate the total mass of MC. The samples were analyzed for 57 parameters, including toxins, nutrients, chlorophyll, and genomics. There were an estimated 11,513 kg and 30,691 kg of MCs in the western basin during the 2018 and 2019 HABs Grabs, respectively. The bloom boundary poses substantial issues for spatial assessments because MC concentration varied by nearly two orders of magnitude over very short distances. The MC to chlorophyll ratio (MC:chl) varied by a factor up to 5.3 throughout the basin, which creates challenges for using MC:chl to predict MC concentrations. Many of the biomass metrics strongly correlated (r > 0.70) with each other except chlorophyll fluorescence and phycocyanin concentration. While MC and chlorophyll correlated well with total phosphorus and nitrogen concentrations, MC:chl correlated with dissolved inorganic nitrogen. More frequent MC data collection can overcome these issues, and models need to account for the MC:chl spatial heterogeneity when forecasting MCs.

Cyanobacterial bloom phenology in Saginaw Bay from MODIS and a comparative look with western Lake Erie

Saginaw Bay and western Lake Erie basin (WLEB) are eutrophic catchments in the Laurentian Great Lakes that experience annual, summer-time cyanobacterial blooms. Both basins share many features including similar size, shallow depths, and equivalent-sized watersheds. They are geographically close and both basins derive a preponderance of their nutrient supply from a single river. Despite these similarities, the bloom phenology in each basin is quite different. The blooms in Saginaw Bay occur at the same time and place and at the same moderate severity level each year. The WLEB, in contrast, exhibits far greater interannual variability in the timing, location, and severity of the bloom than Saginaw Bay, consistent with greater and more variable phosphorus inputs. Saginaw Bay has bloom biomass that corresponds to relatively mild blooms in WLEB, and also has equivalent phosphorus loads. This result suggests that if inputs of P into the WLEB were reduced to similarly sized loads as Saginaw Bay the most severe blooms would be abated. Above 500 t P input, which occur in WLEB, blooms increase non-linearly indicating any reduction in P-input at the highest inputs levels currently occurring in the WLEB, would yield disproportionately large reductions in cyanobacterial bloom intensity. As the maximum phosphorus loads in Saginaw Bay lie just below this inflection point, shifts in the Saginaw Bay watershed toward greater agriculture uses and less wetlands may substantially increase the risk of more intense cyanobacterial blooms than presently occur.

Elucidating controls on cyanobacteria bloom timing and intensity via Bayesian mechanistic modeling

The adverse impacts of harmful algal blooms (HABs) are increasing worldwide. Lake Erie is a North American Great Lake highly affected by cultural eutrophication and summer cyanobacterial HABs. While phosphorus loading is a known driver of bloom size, more nuanced yet crucial questions remain. For example, it is unclear what mechanisms are primarily responsible for initiating cyanobacterial dominance and subsequent biomass accumulation. To address these questions, we develop a mechanistic model describing June-October dynamics of chlorophyll a, nitrogen, and phosphorus near the Maumee River outlet, where blooms typically initiate and are most severe. We calibrate the model to a new, geostatistically-derived dataset of daily water quality spanning 2008-2017. A Bayesian framework enables us to embed prior knowledge on system characteristics and test alternative model formulations. Overall, the best model formulation explains 42% of the variability in chlorophyll a and 83% of nitrogen, and better captures bloom timing than previous models. Our results, supported by cross validation, show that onset of the major midsummer bloom is associated with about a month of water temperatures above 20 °C (occurring 19 July to 6 August), consistent with when cyanobacteria dominance is usually reported. Decreased phytoplankton loss rate is the main factor enabling biomass accumulation, consistent with reduced zooplankton grazing on cyanobacteria. The model also shows that phosphorus limitation is most severe in August, and nitrogen limitation tends to occur in early autumn. Our results highlight the role of temperature in regulating bloom initiation and subsequent loss rates, and suggest that a 2 °C increase could lead to blooms that start about 10 days earlier and grow 23% more intense.

Roles of Nutrient Limitation on Western Lake Erie CyanoHAB Toxin Production

Cyanobacterial harmful algal bloom (CyanoHAB) proliferation is a global problem impacting ecosystem and human health. Western Lake Erie (WLE) typically endures two highly toxic CyanoHABs during summer: a Microcystis spp. bloom in Maumee Bay that extends throughout the western basin, and a Planktothrix spp. bloom in Sandusky Bay. Recently, the USA and Canada agreed to a 40% phosphorus (P) load reduction to lessen the severity of the WLE blooms. To investigate phosphorus and nitrogen (N) limitation of biomass and toxin production in WLE CyanoHABs, we conducted in situ nutrient addition and 40% dilution microcosm bioassays in June and August 2019. During the June Sandusky Bay bloom, biomass production as well as hepatotoxic microcystin and neurotoxic anatoxin production were N and P co-limited with microcystin production becoming nutrient deplete under 40% dilution. During August, the Maumee Bay bloom produced microcystin under nutrient repletion with slight induced P limitation under 40% dilution, and the Sandusky Bay bloom produced anatoxin under N limitation in both dilution treatments. The results demonstrate the importance of nutrient limitation effects on microcystin and anatoxin production. To properly combat cyanotoxin and cyanobacterial biomass production in WLE, both N and P reduction efforts should be implemented in its watershed.

Long-term nitrate removal in three riparian buffers: 21 years of data from the Bear Creek watershed in central Iowa, USA

Riparian buffers are a conservation practice that increases vegetation diversity on the agricultural landscape while providing environmental benefits. This study specifically focused on the ability of riparian buffers to remove nitrate from shallow groundwater. There are many studies that assessed nitrate removal within buffers, but not many have a long-term, continuous data set that can analyze for variation in nitrate removal rates over time. Here we report on 21 years of nitrate well data, from 1996 through 2017, for three buffers in the Bear Creek watershed in central Iowa. These buffers are named using abbreviations to help keep landowners anonymous (e.g. RN, RS, and ST). Studied buffers RS and ST showed greater nitrate reduction (or removal) after 10 and 6 years of its establishment, respectively. Buffer RN did not experience a significant nitrate removal increase with time, but instead had higher nitrate removal rates when compared to buffers RS and ST of 10.3 g NO₃^--N m^-1 day^-1 from the start of this study. From this data, we suggest that past land management played a major role in the responses observed. RN had previously been established in cool-season grasses for grazing before being converted to a buffer, while RS and ST had been managed in a corn and soybean rotation. RN was thought to have higher denitrification immediately with increased labile soil carbon input and enhanced soil aggregation due to the grassland perennials, while buffer vegetation establishment increased soil carbon inputs and soil aggregation over time for RS and ST. These nitrate removal trends would not have been observed without access to long-term, continuous data. This study highlighted the importance of long-term data sets and the need to assess conservation practices over time to determine their longevity and efficiency with time.

The role of phosphorus and nitrogen on chlorophyll a: Evidence from hundreds of lakes

The effect of nutrients on phytoplankton biomass in lakes continues to be a subject of debate by aquatic scientists. However, determining whether or not chlorophyll a (CHL) is limited by phosphorus (P) and/or nitrogen (N) is rarely considered using a probabilistic method in studies of hundreds of lakes across broad spatial extents. Several studies have applied a unified CHL-nutrient relationship to determine nutrient limitation, but pose a risk of ecological fallacy because they neglect spatial heterogeneity in ecological contexts. To examine whether or not CHL is limited by P, N, or both nutrients in hundreds of lakes and across diverse ecological settings, a probabilistic machine learning method, Bayesian Network, was applied. Spatial heterogeneity in ecological context was accommodated by the probabilistic nature of the results. We analyzed data from 1382 lakes in 17 US states to evaluate the cause-effect relationships between CHL and nutrients. Observations of CHL, total phosphorus (TP), and total nitrogen (TN) were discretized into three trophic states (oligo-mesotrophic, eutrophic, and hypereutrophic) to train the model. We found that although both nutrients were related to CHL trophic state, TP was more related to CHL than TN, especially under oligo-mesotrophic and eutrophic CHL conditions. However, when the CHL trophic state was hypereutrophic, both TP and TN were important. These results provide additional evidence that P-limitation is more likely under oligo-mesotrophic or eutrophic CHL conditions and that co-limitation of P and N occurs under hypereutrophic CHL conditions. We also found a decreasing pattern of the TN/TP ratio with increasing CHL concentrations, which might be a key driver for the role change of nutrients. Previous work performed at smaller scales support our findings, indicating potential for extension of our findings to other regions. Our findings enhance the understanding of nutrient limitation at macroscales and revealed that the current debate on the limiting nutrient might be caused by failure to consider CHL trophic state. Our findings also provide prior information for the site-specific eutrophication management of unsampled or data-limited lakes.

Spatiotemporal Variability in Phytoplankton Bloom Phenology in Eastern Canadian Lakes Related to Physiographic, Morphologic, and Climatic Drivers

Phytoplankton bloom monitoring in freshwaters is a challenging task, particularly when biomass is dominated by buoyant cyanobacterial communities that present complex spatiotemporal patterns. Increases in bloom frequency or intensity and their earlier onset in spring were shown to be linked to multiple anthropogenic disturbances, including climate change. The aim of the present study was to describe the phenology of phytoplankton blooms and its potential link with morphological, physiographic, anthropogenic, and climatic characteristics of the lakes and their watershed. The spatiotemporal dynamics of near-surface blooms were studied on 580 lakes in southern Quebec (Eastern Canada) over a 17-year period by analyzing chlorophyll-a concentrations gathered from MODIS (Moderate Resolution Imaging Spectroradiometer) satellite images. Results show a significant increase by 23% in bloom frequency across all studied lakes between 2000 and 2016. The first blooms of the year appeared increasingly early over this period but only by 3 days (median date changing from 6 June to 3 June). Results also indicate that high biomass values are often reached, but the problem is seldom extended to the entire lake surface. The canonical correlation analysis between phenological variables and environmental variables shows that higher frequency and intensity of phytoplankton blooms and earlier onset date occurred for smaller watersheds and higher degree-days, lake surface area, and proportion of urban zones. This study provides a regional picture of lake trophic state over a wide variety of lacustrine environments in Quebec, a detailed phenology allowing to go beyond local biomass assessments, and the first steps on the development of an approach exploiting regional trends for local pattern assessments.

The Composition and Function of Microbiomes Within Microcystis Colonies Are Significantly Different Than Native Bacterial Assemblages in Two North American Lakes

The toxic cyanobacterium Microcystis is one of the most pervasive harmful algal bloom (HAB) genera and naturally occurs in large colonies known to harbor diverse heterotrophic bacterial assemblages. While colony-associated microbiomes may influence Microcystis blooms, there remains a limited understanding of the structure and functional potential of these communities and how they may be shaped by changing environmental conditions. To address this gap, we compared the dynamics of Microcystis-attached (MCA), free-living (FL), and whole water (W) microbiomes during Microcystis blooms using next-generation amplicon sequencing (16S rRNA), a predictive metagenome software, and other bioinformatic approaches. Microbiomes were monitored through high resolution spatial-temporal surveys across two North American lakes, Lake Erie (LE) and Lake Agawam (LA; Long Island, NY, United States) in 2017, providing the largest dataset of these fractions to date. Sequencing of 126 samples generated 7,922,628 sequences that clustered into 7,447 amplicon sequence variants (ASVs) with 100% sequence identity. Across lakes, the MCA microbiomes were significantly different than the FL and W fractions being significantly enriched in Gemmatimonadetes, Burkholderiaceae, Rhizobiales, and Cytophagales and depleted of Actinobacteria. Further, although MCA communities harbored > 900 unique ASVs, they were significantly less diverse than the other fractions with diversity inversely related to bloom intensity, suggesting increased selection pressure on microbial communities as blooms intensified. Despite taxonomic differences between lakes, predicted metagenomes revealed conserved functional potential among MCA microbiomes. MCA communities were significantly enriched in pathways involved in N and P cycling and microcystin-degradation. Taxa potentially capable of N2-fixation were significantly enriched (p < 0.05) and up to four-fold more abundant within the MCA faction relative to other fractions, potentially aiding in the proliferation of Microcystis blooms during low N conditions. The MCA predicted metagenomes were conserved over 8 months of seasonal changes in temperature and N availability despite strong temporal succession in microbiome composition. Collectively, these findings indicate that Microcystis colonies harbor a statistically distinct microbiome with a conserved functional potential that may help facilitate bloom persistence under environmentally unfavorable conditions.

Effects of Detritivores on Nutrient Dynamics and Corn Biomass in Mesocosms

(1) Background: Strategies aimed at managing freshwater eutrophication should be based on practices that consider cropland invertebrates, climatic change, and soil nutrient cycling. Specifically, detritivores play a crucial role in the biogeochemical processes of soil through their consumptive and burrowing activities. Here, we evaluated the effectiveness of increasing detritivore abundance as a strategy for nutrient management under varied rainfall. (2) Methods: We manipulated soil macroinvertebrate abundance and rainfall amount in an agricultural mesocosms. We then measured the phosphorus, nitrogen, and carbon levels within the soil, corn, invertebrates, and soil solution. (3) Results: Increasing detritivore abundance in our soil significantly increased corn biomass by 2.49 g (p < 0.001), reduced weed growth by 18.2% (p < 0.001), and decreased soil solution nitrogen and total organic carbon (p < 0.05) and volume by 31.03 mL (p < 0.001). Detritivore abundance also displayed a significant interaction effect with rainfall treatment to influence soil total P (p = 0.0019), total N (p < 0.001), and total C (p = 0.0146). (4) Conclusions: Soil detritivores play an important role in soil nutrient cycling and soil health. Incorporating soil macroinvertebrate abundance into management strategies for agricultural soil may increase soil health of agroecosystems, preserve freshwater ecosystems, and protect the valuable services they both provide for humans.

Assessment of hydrology and nutrient losses in a changing climate in a subsurface-drained watershed

Studies assessing the impact of subsurface drains on hydrology and nutrient yield in a changing climate are limited, specifically for Western Lake Erie Basin. This study aimed to evaluate the impact of changing climate on hydro-climatology and nutrient loadings in agricultural subsurface-drained areas on a watershed in northeastern Indiana. The study was conducted using a hydrologic model - the Soil and Water Assessment Tool (SWAT) - under two different greenhouse gas emission scenarios (RCP 4.5 and RCP 8.5). Based on analysis, annual subsurface drain flow totals could increase by 70% with respect to the baseline by the end of the 21st century. Surface runoff could increase by 10 to 140% and changes are expected to be greater under RCP 8.5. Soluble phosphorus yield over the basin in a year via subsurface drains could decrease by 30 to 60% under either emission scenarios. Annual total soluble phosphorus yield (soluble phosphorus loading to stream) from subsurface drains and surface runoff could vary from 0.041 to 0.058 kg/ha under RCP 4.5 and 0.035 to 0.064 kg/ha under RCP 8.5 by the end of the 21st century while the values from the baseline model were 0.051 kg/ha. This was attributable to the fact that future climate could have a greater increase in surface runoff than subsurface drain flow based on analysis of the different climate scenarios. Outputs from individual climate model data rather than ensembles provided a band of influence of watershed responses, while outputs from different timelines provided details for evaluating management practice suitability with respect to anticipated differences in climate. Results provide valuable information for stakeholders and policy makers for planning management practices to protect water quality.

A synoptic assessment of water quality in two Great Lakes connecting channels

We conducted a probabilistic water quality assessment of two Great Lakes connecting channels, the St. Marys River, and the Lake Huron-Lake Erie Corridor (HEC) in 2014-2015. We compared the condition of the channels to each other and to the up- and down-river Great Lakes with data from an assessment of the Great Lakes nearshore conducted in 2015. We assessed the condition of each channel as good, fair, or poor by applying the most protective water quality thresholds for the down-channel lake. Condition in the St. Marys River rated mostly fair for total phosphorus (TP) and mostly good for chlorophyll a, and area-weighted mean concentrations were intermediate to nearshore Lake Superior and Lake Huron. A large proportion of the area of the St. Marys River was in poor condition for water clarity based on Secchi depth; while nearshore Lakes Superior and Huron were mostly in good condition for water clarity. Area-weighted mean concentrations of TP and chlorophyll a in the HEC were more like nearshore Lake Huron than Lake Erie. For those indicators, most of the area of the HEC was rated good. The HEC appears more degraded when Lake Huron thresholds are applied rather than Lake Erie thresholds. Appropriate thresholds for the connecting channels should align with assessment objectives and be at least as protective as thresholds for the down-channel lake. Future iterations of this assessment will allow evaluation of water quality trends in the connecting channels.

Effect of micronutrients on algae in different regions of Taihu, a large, spatially diverse, hypereutrophic lake

Eutrophication or excessive nutrient richness is an impairment of many freshwater ecosystems and a prominent cause of harmful algal blooms. It is generally accepted that nitrogen and phosphorus nutrients are the primary causative factor, however, for systems subject to large anthropogenic perturbation, this may no longer be true, and the role of micronutrients is often overlooked. Here we report a study on Lake Tai (Taihu), a large, spatially diverse and hypereutrophic lake in China. We performed small-scale mesocosm nutrient limitation bioassays using boron, iron, cobalt, copper, molybdenum, nitrogen and phosphorus on phytoplankton communities sampled from different locations in Taihu to test the relative effects of micronutrients on in situ algal assemblages. In addition to commonly-used methods of chemical and biological analysis (including algal phytoplankton counting), we used flow cytometry coupled with data-driven analysis to monitor changes to algal assemblages. We found statistically significant effects of limitation or co-limitation for boron, cobalt, copper and iron. For copper at one location chlorophyll-a was over four times higher for amendment with copper, nitrogen and phosphorous than for the latter two alone. Since copper is often proposed as amendment for the environmental management of harmful algal blooms, this result is significant. We have three primary conclusions: first, the strong effects for Cu that we report here are mutually consistent across chlorophyll-a results, count data, and results determined from a data-driven approach to flow cytometry. Given that we cannot rule out a role for a Fe-Cu homeostatic link in causing these effects, future research into MNs and how they interact with N, P, and other MNs should be pursued to explore new interventions for effective management of HABs. Second, in view of the stimulatory effect that Cu exhibited, management of HABs with Cu as an algal biocide may not always be advisable. Third, our approach to flow cytometry offers data confirming our results from chemical and biological analysis, however also holds promise for future development as a high-throughput tool for use in understanding changes in algal assemblages. The results from this study concur with a small and emerging body of literature suggesting that the potential role of micronutrients in eutrophication requires further consideration in environmental management.

The impact of cation concentration on Microcystis (cyanobacteria) scum formation

Cyanobacterial scums at the surface of the lakes are potentially harmful phenomena with increasing occurrence in the last decades, and the causes that lead to their formation are still an unresolved issue. In order to better understand what triggers the scums, we investigated the effect of several Mg²⁺ and Ca²⁺ ion concentrations in promoting them in eight Microcystis aeruginosa strains. The possibility to prevent scum formation by using the ion chelator EDTA was also explored. We found that in some strains the cell aggregation takes place under lower ion source concentrations (20 mM MgSO₄ or CaCl₂), while in others this phenomenon does not occur even at 60 mM concentration. The scum formation correlated to the amount of extracellular polymeric substances (between 234 and 351 µg/cell). EDTA failed to prevent the scum formation in most strains, and in turn it caused cell lysis followed by the release of cellular content into the culture medium. We emphasize the relevance of these results for cyanobacterial scum formation in the environment and we also suggest that controlling the salinity of the medium (by manipulating the ion concentration) is a potentially efficient method for biomass harvesting in large ponds/tanks.

Reduced forms of nitrogen are a driver of non-nitrogen-fixing harmful cyanobacterial blooms and toxicity in Lake Erie

Western Lake Erie (WLE) experiences anthropogenic eutrophication and annual, toxic cyanobacterial blooms of non-nitrogen (N) fixing Microcystis. Numerous studies have shown that bloom biomass is correlated with an increased proportion of soluble reactive phosphorus loading from the Maumee River. Long term monitoring shows that the proportion of the annual Maumee River N load of non-nitrate N, or total Kjeldahl nitrogen (TKN), has also increased significantly (Spearman's ρ = 0.68, p = 0.001) over the last few decades and is also significantly correlated to cyanobacterial bloom biomass (Spearman's ρ = 0.64, p = 0.003). The ratio of chemically reduced N to oxidized N (TKN:NO3) concentrations was also compared to extracted chlorophyll and phycocyanin concentrations from all weekly sampling stations within WLE from 2009 to 2015. Both chlorophyll (Spearman's ρ = 0.657, p < 0.0001) and phycocyanin (Spearman's ρ = 0.714, p < 0.0001) were significantly correlated with TKN:NO3. This correlation between the increasing fraction of chemically reduced N from the Maumee River and increasing bloom biomass demonstrates the urgent need to control N loading, in addition to current P load reductions, to WLE and similar systems impacted by non-N-fixing, toxin-producing cyanobacteria.

A Hybrid Lagrangian–Eulerian Particle Model for Ecosystem Simulation

Current numerical methods for simulating biophysical processes in aquatic environments are typically constructed in a grid-based Eulerian framework or as an individual-based model in a particle-based Lagrangian framework. Often, the biogeochemical processes and physical (hydrodynamic) processes occur at different time and space scales, and changes in biological processes do not affect the hydrodynamic conditions. Therefore, it is possible to develop an alternative strategy to grid-based approaches for linking hydrodynamic and biogeochemical models that can significantly improve computational efficiency for this type of linked biophysical model. In this work, we utilize a new technique that links hydrodynamic effects and biological processes through a property-carrying particle model (PCPM) in a Lagrangian/Eulerian framework. The model is tested in idealized cases and its utility is demonstrated in a practical application to Sandusky Bay. Results show the integration of Lagrangian and Eulerian approaches allows for a natural coupling of mass transport (represented by particle movements and random walk) and biological processes in water columns which is described by a nutrient-phytoplankton-zooplankton-detritus (NPZD) biological model. This method is far more efficient than traditional tracer-based Eulerian biophysical models for 3-D simulation, particularly for a large domain and/or ensemble simulations.

Accuracy of data buoys for measurement of cyanobacteria, chlorophyll, and turbidity in a large lake (Lake Erie, North America): implications for estimation of cyanobacterial bloom parameters from water quality sonde measurements

Microcystin (MCY)-producing harmful cyanobacterial blooms (cHABs) are an annual occurrence in Lake Erie, and buoys equipped with water quality sondes have been deployed to help researchers and resource managers track cHABs. The objective of this study was to determine how well water quality sondes attached to buoys measure total algae and cyanobacterial biomass and water turbidity. Water samples were collected next to two data buoys in western Lake Erie (near Gibraltar Island and in the Sandusky subbasin) throughout summers 2015, 2016, and 2017 to determine correlations between buoy sonde data and water sample data. MCY and nutrient concentrations were also measured. Significant (P < 0.001) linear relationships (R² > 0.75) occurred between cyanobacteria buoy and water sample data at the Gibraltar buoy, but not at the Sandusky buoy; however, the coefficients at the Gibraltar buoy differed significantly across years. There was a significant correlation between buoy and water sample total chlorophyll data at both buoys, but the coefficient varied considerably between buoys and among years. Total MCY concentrations at the Gibraltar buoy followed similar temporal patterns as buoy and water sample cyanobacterial biomass data, and the ratio of MCY to cyanobacteria-chlorophyll decreased with decreased ambient nitrate concentrations. These results suggest that buoy data are difficult to compare across time and space. Additionally, the inclusion of nitrate concentration data can lead to more robust predictions on the relative toxicity of blooms. Overall, deployed buoys with sondes that are routinely cleaned and calibrated can track relative cyanobacteria abundance and be used as an early warning system for potentially toxic blooms.

Long-term Cropping Effects on Partitioning of Water Flow and Nitrate Loss between Surface Runoff and Tile Drainage

Surface runoff and tile drainage are the main pathways for water movement and entry of agricultural nitrate into water resources. The objective of this 5-yr study was to characterize the partitioning of water flow and nitrate loss between these pathways for a humid-temperate Brookston clay loam soil under 54 to 59 yr of consistent cropping and fertilization. Cropping treatments included monoculture corn ( L., MC), continuous bluegrass ( L.) sod (CS), and a corn-oat-alfalfa ( L.)-alfalfa rotation (RC-RO-RA1-RA2). Fertilization treatments included annual fertilizer addition (F) and no fertilizer addition (NF). Tile drainage and surface runoff occurred primarily during the nongrowing season (November-April), and they were highly correlated with the mean saturated hydraulic conductivity of the near-surface soil profile. Tile drainage accounted for 69 to 90% of cumulative water flow and 79 to 96% of cumulative nitrate loss from fertilized rotation and CS, whereas surface runoff accounted for the majority of the nitrate losses in MC (i.e., 75-93% of water flow and 65-96% of nitrate loss). Cumulative nitrate losses were highest in the RC-F (152 kg N ha), RC-NF (101 kg N ha), RA2-F (121 kg N ha), and RA2-NF (75 kg N ha) plots, and these high losses are attributed to N mineralization from the plowed alfalfa and fertilization (if applicable). Fertilization increased cumulative nitrate loss in tile drainage from all treatments, whereas no fertilization increased cumulative nitrate loss in surface runoff from the rotation. Cropping system and fertilization on clay loam soil changed how water flow and nitrate loss were partitioned between tile drainage and surface runoff.

Internal phosphorus loading from sediments causes seasonal nitrogen limitation for harmful algal blooms

It is proposed that the internal loading of phosphorus (P) from sediments plays an important role in seasonal nitrogen (N) limitation for harmful algal blooms (HABs), although there is a lack of experimental evidence. In this study, an eutrophic bay from the large and shallow Lake Taihu was studied for investigating the contribution of internal P to N limitation over one-year field sampling (February 2016 to January 2017). A prebloom-bloom period was identified from February to August according to the increase in Chla concentration in the water column, during which the ratio of total N to total P (TN/TP) exponentially decreased with month from 43.4 to 7.4. High-resolution dialysis (HR-Peeper) and diffusive gradients in thin films (DGT) analysis showed large variations in the vertical distribution of mobile P (SRP and DGT-labile P) in sediments, resulting in the SRP diffusion flux at the sediment-water interface ranging from -0.01 to 6.76mg/m²/d (minus sign denotes downward flux). Significant and linear correlations existed between SRP and soluble Fe(II) concentrations in pore water, reflecting that the spatial-temporal variation in mobile P was controlled by microbe-mediated Fe redox cycling. Mass estimation showed that the cumulative flux of SRP from sediments accounted for 54% of the increase in TP observed in the water column during the prebloom-bloom period. These findings are supported by the significantly negative correlation (p<0.01) observed between sediment SRP flux and water column TN/TP during the same period. Overall, these results provide solid evidence for the major role of internal P loading in causing N limitation during the prebloom-bloom period.

Nitrogen limitation, toxin synthesis potential, and toxicity of cyanobacterial populations in Lake Okeechobee and the St. Lucie River Estuary, Florida, during the 2016 state of emergency event

Lake Okeechobee, FL, USA, has been subjected to intensifying cyanobacterial blooms that can spread to the adjacent St. Lucie River and Estuary via natural and anthropogenically-induced flooding events. In July 2016, a large, toxic cyanobacterial bloom occurred in Lake Okeechobee and throughout the St. Lucie River and Estuary, leading Florida to declare a state of emergency. This study reports on measurements and nutrient amendment experiments performed in this freshwater-estuarine ecosystem (salinity 0–25 PSU) during and after the bloom. In July, all sites along the bloom exhibited dissolved inorganic nitrogen-to-phosphorus ratios < 6, while Microcystis dominated (> 95%) phytoplankton inventories from the lake to the central part of the estuary. Chlorophyll a and microcystin concentrations peaked (100 and 34 μg L^-1, respectively) within Lake Okeechobee and decreased eastwards. Metagenomic analyses indicated that genes associated with the production of microcystin (mcyE) and the algal neurotoxin saxitoxin (sxtA) originated from Microcystis and multiple diazotrophic genera, respectively. There were highly significant correlations between levels of total nitrogen, microcystin, and microcystin synthesis gene abundance across all surveyed sites (p < 0.001), suggesting high levels of nitrogen supported the production of microcystin during this event. Consistent with this, experiments performed with low salinity water from the St. Lucie River during the event indicated that algal biomass was nitrogen-limited. In the fall, densities of Microcystis and concentrations of microcystin were significantly lower, green algae co-dominated with cyanobacteria, and multiple algal groups displayed nitrogen-limitation. These results indicate that monitoring and regulatory strategies in Lake Okeechobee and the St. Lucie River and Estuary should consider managing loads of nitrogen to control future algal and microcystin-producing cyanobacterial blooms.

Interactions between nitrogen form, loading rate, and light intensity on Microcystis and Planktothrix growth and microcystin production

The toxin-producing, bloom-forming cyanobacterial genera Microcystis and Planktothrix require fixed nitrogen (N), such as nitrate, ammonium, or organic N (e.g., urea) for growth and production of microcystins (MC). Bioavailable N can enter lakes in pulses via tributary discharge and through in-lake recycling, which can maintain low N concentrations. Additionally, light intensity has been suggested to play a role in MC production. This study examined how three forms of N (nitrate, ammonium, and urea) interacted with N loading rate (one large pulse vs. many small pulses) and light intensity to stimulate Microcystis and Planktothrix growth and MC production using nutrient enrichment experiments. Enrichments of nitrate, ammonium, and urea resulted in greater cyanobacterial biovolumes and MC concentrations than phosphorus-only enrichments, and there was no difference between pulse (100 μmol/L) and press treatments (8.3 μmol/L every 4 h). Analysis of mcyD transcripts showed significant up-regulation within 4 h of ammonium and urea enrichment. High light intensities (300 μmol photons/m²/s) with N enrichment resulted in greater cyanobacterial biovolumes and MC concentrations than lower light intensities (30 and 3 μmol photons/m²/s). Overall, the results suggest Microcystis and Planktothrix can use many forms of N and that high light intensities enhance MC production during elevated N concentrations. Moreover, the results here further demonstrate the importance of considering N, as well as P, in management strategies aimed at mitigating cyanobacterial blooms.

A Nitrogen Physical Input-Output Model for Illinois

Nitrogen (N) presents an important challenge for sustainability. Human intervention in the global nitrogen cycle has been pivotal in in providing goods and services to society. However, release of N beyond its intended societal use has many negative health and environmental consequences. Several systems modeling approaches have been developed to understand the trade-offs between the beneficial and harmful effects of N. These efforts include life cycle modeling, integrated management practices and sustainability metrics for individuals and communities. However, these approaches do not connect economic and ecological N flows in physical units throughout the system, which could better represent these trade-offs for decision-makers. Physical Input-Output Table (PIOT) based models present a viable complementary solution to overcome this limitation. We developed a N-PIOT for Illinois representing the interdependence of sectors in 2002, using N mass units. This allows studying the total N flow required to produce a certain amount of N in the final product. An Environmentally Extended Input Output (EEIO) based approach was used to connect the physical economic production to environmental losses; allowing quantification of total environmental impact to support agricultural production in Illinois. A bottom up approach was used to develop the N-PIOT using Material Flow Analysis (MFA) tracking N flows associated with top 3 commodities (Corn, Soybean and Wheat). These three commodities cover 99% of N fertilizer use in Illinois. The PIOT shows that of all the N inputs to corn production the state exported 68% of N embedded in useful products, 9% went to animal feed manufacturing and only 0.03% was consumed directly within the state. Approximately 35% of N input to soybean farming ended up in animal feed. Release of N to the environment was highest from corn farming, at about 21.8% of total N fertilizer inputs, followed by soybean (9.2%) and wheat farming (4.2%). The model also allowed the calculation of life cycle N use efficiency for N based on physical flows in the economy. Hence, PIOTs prove to be a viable tool for developing a holistic approach to manage disrupted biogeochemical cycles, since these provide a detailed insight into physical flows in economic systems and allow physical coupling with ecological N flows.

Ecophysiological Examination of the Lake Erie Microcystis Bloom in 2014: Linkages between Biology and the Water Supply Shutdown of Toledo, OH

Annual cyanobacterial blooms dominated by Microcystis have occurred in western Lake Erie (U.S./Canada) during summer months since 1995. The production of toxins by bloom-forming cyanobacteria can lead to drinking water crises, such as the one experienced by the city of Toledo in August of 2014, when the city was rendered without drinking water for >2 days. It is important to understand the conditions and environmental cues that were driving this specific bloom to provide a scientific framework for management of future bloom events. To this end, samples were collected and metatranscriptomes generated coincident with the collection of environmental metrics for eight sites located in the western basin of Lake Erie, including a station proximal to the water intake for the city of Toledo. These data were used to generate a basin-wide ecophysiological fingerprint of Lake Erie Microcystis populations in August 2014 for comparison to previous bloom communities. Our observations and analyses indicate that, at the time of sample collection, Microcystis populations were under dual nitrogen (N) and phosphorus (P) stress, as genes involved in scavenging of these nutrients were being actively transcribed. Targeted analysis of urea transport and hydrolysis suggests a potentially important role for exogenous urea as a nitrogen source during the 2014 event. Finally, simulation data suggest a wind event caused microcystin-rich water from Maumee Bay to be transported east along the southern shoreline past the Toledo water intake. Coupled with a significant cyanophage infection, these results reveal that a combination of biological and environmental factors led to the disruption of the Toledo water supply. This scenario was not atypical of reoccurring Lake Erie blooms and thus may reoccur in the future.

Spatial and temporal variability in the nitrogen cyclers of hypereutrophic Lake Taihu

Harmful cyanobacterial blooms (cyanoHABs) are a major threat to freshwater ecosystems worldwide. Evidence suggests that both nitrogen and phosphorus are important nutrients in the development and proliferation of blooms, yet much less is known about nitrogen cycling dynamics in these systems. To assess the potential nitrogen cycling function of the cyanoHAB community, surface water samples were collected in Lake Tai (Taihu), China over a 5-month bloom event in 2014. The expression of six nitrogen cycling genes (nifH, hzsA, nxrB, nrfA, amoA, nosZ) was surveyed using a targeted microarray with probes designed to provide phylogenetic information. N-Cycling gene expression varied spatially across Taihu, most notably near the mouth of the Dapu River. Expression of nifH was observed across the lake and attributable to both Proteobacteria and Cyanobacteria: Proteobacteria were major contributors to nifH signal near shore. Other N transformations such as anaerobic ammonia oxidation and denitrification were evident in the surface waters as well. Observations in this study highlight the potential importance of heterotrophic bacteria in N-cycling associated with cyanoHABs.

Characterizing nutrient distributions and fluxes in a eutrophic reservoir, Midwestern United States

Harmful algal blooms are increasingly common in aquatic ecosystems and have been linked to runoff from agricultural land. This study investigated the internal nutrient (i.e., phosphorus (P) and nitrogen (N)) dynamics of a eutrophic reservoir in the Midwestern United States to constrain the potential for sedimentary nutrients to stimulate harmful algal blooms. The spatial distribution of nutrients in the water column (soluble reactive P (SRP), nitrate/nitrite-N (NO_x-N), and ammonium-N (NH₄⁺-N)) and sediments (total P, total carbon (C), total N, and organic matter (OM)) were quantified and mapped. Water column nutrients varied spatially and temporally, with generally higher concentrations near the dam wall during normal lake levels. The upper portion of the lake, near the inlet, was sampled during a flood event and had overall higher nutrient concentrations and lower chlorophyll levels compared to normal lake level samples. Mean sedimentary total P (936mg/kg) was ~30% higher in the reservoir than the surrounding upland soils, with the highest concentrations near the dam wall (1661mg/kg) and a significant positive correlation found between sedimentary total P, total C, and OM. Additionally, 15 intact sediment cores were manipulated ex situ to examine mechanisms of nutrient flux across the sediment-water interface (SWI) that may trigger algal blooms. Core treatment conditions included advection (i.e., simulating potential nutrient fluxes during wind events through sediment resuspension) and diffusion. Core experiments indicated both advective and diffusive conditions at the SWI may trigger the flux of nutrients important for algal growth from lake sediments, with diffusion contributing both N and P to the water column, while intense advection increased water column N, but decreased P. Release of P to the water column may be more diffusion-driven than advection-driven, whereas N release to the water column appears to be both diffusion- and advection-driven.

Impacts of the 2014 severe drought on the Microcystis bloom in San Francisco Estuary

The increased frequency and intensity of drought with climate change may cause an increase in the magnitude and toxicity of freshwater cyanobacteria harmful algal blooms (CHABs), including Microcystis blooms, in San Francisco Estuary, California. As the fourth driest year on record in San Francisco Estuary, the 2014 drought provided an opportunity to directly test the impact of severe drought on cyanobacteria blooms in SFE. A field sampling program was conducted between July and December 2014 to sample a suite of physical, chemical, and biological variables at 10 stations in the freshwater and brackish reaches of the estuary. The 2014 Microcystis bloom had the highest biomass and toxin concentration, earliest initiation, and the longest duration, since the blooms began in 1999. Median chlorophyll a concentration increased by 9 and 12 times over previous dry and wet years, respectively. Total microcystin concentration also exceeded that in previous dry and wet years by a factor of 11 and 65, respectively. Cell abundance determined by quantitative PCR indicated the bloom contained multiple potentially toxic cyanobacteria species, toxic Microcystis and relatively high total cyanobacteria abundance. The bloom was associated with extreme nutrient concentrations, including a 20-year high in soluble reactive phosphorus concentration and low to below detection levels of ammonium. Stable isotope analysis suggested the bloom varied with both inorganic and organic nutrient concentration, and used ammonium as the primary nitrogen source. Water temperature was a primary controlling factor for the bloom and was positively correlated with the increase in both total and toxic Microcystis abundance. In addition, the early initiation and persistence of warm water temperature coincided with the increased intensity and duration of the Microcystis bloom from the usual 3 to 4 months to 8 months. Long residence time was also a primary factor controlling the magnitude and persistence of the bloom, and was created by a 66% to 85% reduction in both the water inflow and diversion of water for agriculture during the summer. We concluded that severe drought conditions can lead to a significant increase in the abundance of Microcystis and other cyanobacteria, as well as their associated toxins.

It Takes Two to Tango: When and Where Dual Nutrient (N & P) Reductions Are Needed to Protect Lakes and Downstream Ecosystems

Preventing harmful algal blooms (HABs) is needed to protect lakes and downstream ecosystems. Traditionally, reducing phosphorus (P) inputs was the prescribed solution for lakes, based on the assumption that P universally limits HAB formation. Reduction of P inputs has decreased HABs in many lakes, but was not successful in others. Thus, the "P-only" paradigm is overgeneralized. Whole-lake experiments indicate that HABs are often stimulated more by combined P and nitrogen (N) enrichment rather than N or P alone, indicating that the dynamics of both nutrients are important for HAB control. The changing paradigm from P-only to consideration of dual nutrient control is supported by studies indicating that (1) biological N fixation cannot always meet lake ecosystem N needs, and (2) that anthropogenic N and P loading has increased dramatically in recent decades. Sediment P accumulation supports long-term internal loading, while N may escape via denitrification, leading to perpetual N deficits. Hence, controlling both N and P inputs will help control HABs in some lakes and also reduce N export to downstream N-sensitive ecosystems. Managers should consider whether balanced control of N and P will most effectively reduce HABs along the freshwater-marine continuum.

Comparative bioavailability of ammonium, nitrate, nitrite and urea to typically harmful cyanobacterium Microcystis aeruginosa

Phosphorus is generally considered as the prime limiting nutrient responsible for cyanobacterial blooms. However, recent research is drawing attention to the importance of bioavailable nitrogen (N) in freshwater eutrophication. This study investigated the bioavailability of NO3(-)-N, NO2(-)-N, NH4(+)-N and Urea-N under different concentrations of 1.2, 3.6 and 6.0mgL(-1) to Microcystis aeruginosa. Overall, Urea-N ranked the first in promoting M. aeruginosa growth, followed by NO3(-)-N and NO2(-)-N. However, the algal growth cultured in NH4(+)-N was depressed under test N levels. The bioavailability of N to M. aeruginosa was seriously influenced by both N forms and N concentrations (p<0.01). Total N concentrations in Urea-N treatment decreased the fastest, which were corresponding with the μ values of M. aeruginosa. The high enzymic activities of nitrate reductase, nitrite reductase and glutamine synthetase indicated that the decomposition process for urea is effective, which contributed in N assimilation and utilization in M. aeruginosa cells.

The dual role of nitrogen supply in controlling the growth and toxicity of cyanobacterial blooms

Historically, phosphorus (P) has been considered the primary limiting nutrient for phytoplankton assemblages in freshwater ecosystems. This review, supported by new findings from Lake Erie, highlights recent molecular, laboratory, and field evidence that the growth and toxicity of some non-diazotrophic blooms of cyanobacteria can be controlled by nitrogen (N). Cyanobacteria such as Microcystis possess physiological adaptations that allow them to dominate low-P surface waters, and in temperate lakes, cyanobacterial densities can be controlled by N availability. Beyond total cyanobacterial biomass, N loading has been shown to selectively promote the abundance of Microcystis and Planktothrix strains capable of synthesizing microcystins over strains that do not possess this ability. Among strains of cyanobacteria capable of synthesizing the N-rich microcystins, cellular toxin quotas have been found to depend upon exogenous N supplies. Herein, multi-year observations from western Lake Erie are presented demonstrating that microcystin concentrations peak in parallel with inorganic N, but not orthophosphate, concentrations and are significantly lower (p<0.01) during years of reduced inorganic nitrogen loading and concentrations. Collectively, this information underscores the importance of N as well as P in controlling toxic cyanobacteria blooms. Furthermore, it supports the premise that management actions to reduce P in the absence of concurrent restrictions on N loading may not effectively control the growth and/or toxicity of non-diazotrophic toxic cyanobacteria such as the cosmopolitan, toxin-producing genus, Microcystis.

Nutrient-Controlled Niche Differentiation of Western Lake Erie Cyanobacterial Populations Revealed via Metatranscriptomic Surveys

Although toxic cyanobacterial blooms in western Lake Erie threaten drinking water supplies and are promoted by nutrient loading, the precise nutrient regime that selects specific cyanobacteria populations is poorly understood. Here, we assess shifts in cyanobacterial abundances and global gene-expression patterns in response to natural and manipulated gradients in nitrogen and phosphorus to identify gene pathways that facilitate dominance by different cyanobacteria. Gradients in soluble reactive phosphorus shaped cyanobacterial communities and elicited the largest transcriptomic responses. Under high-P conditions (closest to the mouth of the Maumee River), Anabaena and Planktothrix were the dominant cyanobacterial populations, and experimental P and ammonium enrichment promoted nitrogen fixation gene (nifH) expression in Anabaena. For Microcystis, experimental additions of P up-regulated genes involved in phage defense, genomic rearrangement, and nitrogen acquisition but led to lower abundances. Within offshore, low-P regions of the western basin of Lake Erie, Microcystis up-regulated genes associated with P scavenging (pstSCAB, phoX) and dominated cyanobacterial communities. Experimental additions of ammonium and urea did not alter Microcystis abundances but did up-regulate protease inhibitors (aer and mcn gene sets) and microcystin synthetase genes (mcy), with urea enrichment yielding significant increases in microcystin concentrations. Our findings suggest that management plans that reduce P loads alone may not significantly reduce the risk of cyanobacterial blooms in western Lake Erie but rather may promote a shift among cyanobacterial populations (Microcystis, Anabaena, and Planktothrix) toward a greater dominance by toxic strains of Microcystis.

Effects of increasing nitrogen and phosphorus concentrations on phytoplankton community growth and toxicity during Planktothrix blooms in Sandusky Bay, Lake Erie

Sandusky Bay experiences annual toxic cyanobacterial blooms dominated by Planktothrix agardhii/suspensa. To further understand the environmental drivers of these events, we evaluated changes in the growth response and toxicity of the Planktothrix-dominated blooms to nutrient amendments with orthophosphate (PO4) and inorganic and organic forms of dissolved nitrogen (N; ammonium (NH4), nitrate (NO3) and urea) over the bloom season (June - October). We complemented these with a metagenomic analysis of the planktonic microbial community. Our results showed that bloom growth and microcystin (MC) concentrations responded more frequently to additions of dissolved N than PO4, and that the dual addition of NH4 + PO4 and Urea + PO4 yielded the highest MC concentrations in 54% of experiments. Metagenomic analysis confirmed that P. agardhii/suspensa was the primary MC producer. The phylogenetic distribution of nifH revealed that both heterocystous cyanobacteria and heterotrophic proteobacteria had the genetic potential for N2 fixation in Sandusky Bay. These results suggest that as best management practices are developed for P reductions in Sandusky Bay, managers must be aware of the negative implications of not managing N loading into this system as N may significantly impact cyanobacterial bloom size and toxicity.