Exploring nutrient and light limitation of algal production in a shallow turbid reservoir.

Harmful algal blooms are increasingly recognized as a threat to the integrity of freshwater reservoirs, which serve as water supplies, wildlife habitats, and recreational attractions. While algal growth and accumulation is controlled by many environmental factors, the relative importance of these factors is unclear, particularly for turbid eutrophic systems. Here we develop and compare two models that test the relative importance of vertical mixing, light, and nutrients for explaining chlorophyll-a variability in shallow (2-3 m) embayments of a eutrophic reservoir, Jordan Lake, North Carolina. One is a multiple linear regression (statistical) model and the other is a process-based (mechanistic) model. Both models are calibrated using a 15-year data record of chlorophyll-a concentration (2003-2018) for the seasonal period of cyanobacteria dominance (June-October). The mechanistic model includes a novel representation of vertical mixing and is calibrated in a Bayesian framework, which allows for data-driven inference of important process rates. Both models show that chlorophyll-a concentration is much more responsive to nutrient variability than mixing, light, or temperature. While both models explain approximately 60% of the variability in chlorophyll-a, the mechanistic model is more robust in cross-validation and provides a more comprehensive assessment of algal drivers. Overall, these models indicate that nutrient reductions, rather than changes in mixing or background turbidity, are critical to controlling cyanobacteria in a shallow eutrophic freshwater system.

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.

Bayesian mechanistic modeling characterizes Gulf of Mexico hypoxia: 1968-2016 and future scenarios.

The hypoxic zone in the northern Gulf of Mexico is among the most dramatic examples of impairments to aquatic ecosystems. Despite having attracted substantial attention, management of this environmental crisis remains challenging, partially due to limited monitoring to support model development and long-term assessments. Here, we leverage new geostatistical estimates of hypoxia derived from nearly 150 monitoring cruises and a process-based model to improve characterization of controlling mechanisms, historic trends, and future responses of hypoxia while rigorously quantifying uncertainty in a Bayesian framework. We find that November-March nitrogen loads are important controls of sediment oxygen demand, which appears to be the major oxygen sink. In comparison, only ~23% of oxygen in the near-bottom region appears to be consumed by net water column respiration, which is driven by spring and summer loads. Hypoxia typically exceeds 15,600 km2 in June, peaks in July, and declines below 10,000 km2 in September. In contrast to some previous Gulf hindcasting studies, our simulations demonstrate that hypoxia was both severe and worsening prior to 1985, and has remained relatively stable since that time. Scenario analysis shows that halving nutrient loadings will reduce hypoxia by 37% with respect to 13,900 km2 (1985-2016 median), while a +2°C change in water temperature will cause a 26% hypoxic area increase due to enhanced sediment respiration and reduced oxygen solubility. These new results highlight the challenges of achieving hypoxia reduction targets, particularly under warming conditions, and should be considered in ecosystem management.

A space-time geostatistical model for probabilistic estimation of harmful algal bloom biomass and areal extent.

Harmful algal blooms (HABs) have been increasing in intensity worldwide, including the western basin of Lake Erie. Substantial efforts have been made to track these blooms using in situ sampling and remote sensing. However, such measurements do not fully capture HAB spatial and temporal dynamics due to the limitations of discrete shipboard sampling over large areas and the effects of clouds and winds on remote sensing estimates. To address these limitations, we develop a space-time geostatistical modeling framework for estimating HAB intensity and extent using chlorophyll a data sampled during the HAB season (June-October) from 2008 to 2017 by five independent monitoring programs. Based on the Bayesian information criterion for model selection, trend variables explain bloom northerly and easterly expansion from Maumee Bay, wind effects over depth, and variability among sampling methods. Cross validation results demonstrate that space-time kriging explains over half of the variability in daily, location-specific chlorophyll observations, on average. Conditional simulations provide, for the first time, comprehensive estimates of overall bloom biomass (based on depth-integrated concentrations) and surface areal extent with quantified uncertainties. These new estimates are contrasted with previous Lake Erie HAB monitoring studies, and deviations among estimates are explored and discussed. Overall, results highlight the importance of maintaining sufficient monitoring coverage to capture bloom dynamics, as well as the benefits of the proposed approach for synthesizing data from multiple monitoring programs to improve estimation accuracy while reducing uncertainty.

Long-Term Phosphorus Loading and Springtime Temperatures Explain Interannual Variability of Hypoxia in a Large Temperate Lake.

Anthropogenic eutrophication has led to the increased occurrence of hypoxia in inland and coastal waters around the globe. While low dissolved oxygen conditions are known to be driven primarily by nutrient loading and water column stratification, the relative importance of these factors and their associated time scales are not well understood. Here, we explore these questions for Lake Erie, a large temperate lake that experiences widespread annual summertime hypoxia. We leverage a three-decade data set of summertime hypoxic extent (1985-2015) and examine the role of seasonal and long-term nutrient loading, as well as hydrometeorological conditions. We find that a linear combination of decadal total phosphorus loading from tributaries and springtime air temperatures explains a high proportion of the interannual variability in average summertime hypoxic extent (R2 = 0.71). This result suggests that the lake responds primarily to long-term variations in phosphorus inputs, rather than springtime or annual loading as previously assumed, which is consistent with internal phosphorus loading from lake sediments likely being an important contributing mechanism. This result also demonstrates that springtime temperatures have a substantial impact on summertime hypoxia, likely by impacting the timing of onset of thermal stratification. These findings imply that management strategies based on reducing tributary phosphorus loading would take several years to reap full benefits, and that projected future increases in temperatures are likely to exacerbate hypoxia in Lake Erie and other temperate lakes.