Progress and opportunities in advancing near-term forecasting of freshwater quality.

Near-term freshwater forecasts, defined as sub-daily to decadal future predictions of a freshwater variable with quantified uncertainty, are urgently needed to improve water quality management as freshwater ecosystems exhibit greater variability due to global change. Shifting baselines in freshwater ecosystems due to land use and climate change prevent managers from relying on historical averages for predicting future conditions, necessitating near-term forecasts to mitigate freshwater risks to human health and safety (e.g., flash floods, harmful algal blooms) and ecosystem services (e.g., water-related recreation and tourism). To assess the current state of freshwater forecasting and identify opportunities for future progress, we synthesized freshwater forecasting papers published in the past 5 years. We found that freshwater forecasting is currently dominated by near-term forecasts of water quantity and that near-term water quality forecasts are fewer in number and in the early stages of development (i.e., non-operational) despite their potential as important preemptive decision support tools. We contend that more freshwater quality forecasts are critically needed and that near-term water quality forecasting is poised to make substantial advances based on examples of recent progress in forecasting methodology, workflows, and end-user engagement. For example, current water quality forecasting systems can predict water temperature, dissolved oxygen, and algal bloom/toxin events 5 days ahead with reasonable accuracy. Continued progress in freshwater quality forecasting will be greatly accelerated by adapting tools and approaches from freshwater quantity forecasting (e.g., machine learning modeling methods). In addition, future development of effective operational freshwater quality forecasts will require substantive engagement of end users throughout the forecast process, funding, and training opportunities. Looking ahead, near-term forecasting provides a hopeful future for freshwater management in the face of increased variability and risk due to global change, and we encourage the freshwater scientific community to incorporate forecasting approaches in water quality research and management.

Anoxia decreases the magnitude of the carbon, nitrogen, and phosphorus sink in freshwaters.

Oxygen availability is decreasing in many lakes and reservoirs worldwide, raising the urgency for understanding how anoxia (low oxygen) affects coupled biogeochemical cycling, which has major implications for water quality, food webs, and ecosystem functioning. Although the increasing magnitude and prevalence of anoxia has been documented in freshwaters globally, the challenges of disentangling oxygen and temperature responses have hindered assessment of the effects of anoxia on carbon, nitrogen, and phosphorus concentrations, stoichiometry (chemical ratios), and retention in freshwaters. The consequences of anoxia are likely severe and may be irreversible, necessitating ecosystem-scale experimental investigation of decreasing freshwater oxygen availability. To address this gap, we devised and conducted REDOX (the Reservoir Ecosystem Dynamic Oxygenation eXperiment), an unprecedented, 7-year experiment in which we manipulated and modeled bottom-water (hypolimnetic) oxygen availability at the whole-ecosystem scale in a eutrophic reservoir. Seven years of data reveal that anoxia significantly increased hypolimnetic carbon, nitrogen, and phosphorus concentrations and altered elemental stoichiometry by factors of 2-5× relative to oxic periods. Importantly, prolonged summer anoxia increased nitrogen export from the reservoir by six-fold and changed the reservoir from a net sink to a net source of phosphorus and organic carbon downstream. While low oxygen in freshwaters is thought of as a response to land use and climate change, results from REDOX demonstrate that low oxygen can also be a driver of major changes to freshwater biogeochemical cycling, which may serve as an intensifying feedback that increases anoxia in downstream waterbodies. Consequently, as climate and land use change continue to increase the prevalence of anoxia in lakes and reservoirs globally, it is likely that anoxia will have major effects on freshwater carbon, nitrogen, and phosphorus budgets as well as water quality and ecosystem functioning.

Near-term phytoplankton forecasts reveal the effects of model time step and forecast horizon on predictability.

As climate and land use increase the variability of many ecosystems, forecasts of ecological variables are needed to inform management and use of ecosystem services. In particular, forecasts of phytoplankton would be especially useful for drinking water management, as phytoplankton populations are exhibiting greater fluctuations due to human activities. While phytoplankton forecasts are increasing in number, many questions remain regarding the optimal model time step (the temporal frequency of the forecast model output), time horizon (the length of time into the future a prediction is made) for maximizing forecast performance, as well as what factors contribute to uncertainty in forecasts and their scalability among sites. To answer these questions, we developed near-term, iterative forecasts of phytoplankton 1-14 days into the future using forecast models with three different time steps (daily, weekly, fortnightly), that included a full uncertainty partitioning analysis at two drinking water reservoirs. We found that forecast accuracy varies with model time step and forecast horizon, and that forecast models can outperform null estimates under most conditions. Weekly and fortnightly forecasts consistently outperformed daily forecasts at 7-day and 14-day horizons, a trend that increased up to the 14-day forecast horizon. Importantly, our work suggests that forecast accuracy can be increased by matching the forecast model time step to the forecast horizon for which predictions are needed. We found that model process uncertainty was the primary source of uncertainty in our phytoplankton forecasts over the forecast period, but parameter uncertainty increased during phytoplankton blooms and when scaling the forecast model to a new site. Overall, our scalability analysis shows promising results that simple models can be transferred to produce forecasts at additional sites. Altogether, our study advances our understanding of how forecast model time step and forecast horizon influence the forecastability of phytoplankton dynamics in aquatic systems and adds to the growing body of work regarding the predictability of ecological systems broadly.

Using near-term forecasts and uncertainty partitioning to inform prediction of oligotrophic lake cyanobacterial density.

Near-term ecological forecasts provide resource managers advance notice of changes in ecosystem services, such as fisheries stocks, timber yields, or water quality. Importantly, ecological forecasts can identify where there is uncertainty in the forecasting system, which is necessary to improve forecast skill and guide interpretation of forecast results. Uncertainty partitioning identifies the relative contributions to total forecast variance introduced by different sources, including specification of the model structure, errors in driver data, and estimation of current states (initial conditions). Uncertainty partitioning could be particularly useful in improving forecasts of highly variable cyanobacterial densities, which are difficult to predict and present a persistent challenge for lake managers. As cyanobacteria can produce toxic and unsightly surface scums, advance warning when cyanobacterial densities are increasing could help managers mitigate water quality issues. Here, we fit 13 Bayesian state-space models to evaluate different hypotheses about cyanobacterial densities in a low nutrient lake that experiences sporadic surface scums of the toxin-producing cyanobacterium, Gloeotrichia echinulata. We used data from several summers of weekly cyanobacteria samples to identify dominant sources of uncertainty for near-term (1- to 4-week) forecasts of G. echinulata densities. Water temperature was an important predictor of cyanobacterial densities during model fitting and at the 4-week forecast horizon. However, no physical covariates improved model performance over a simple model including the previous week's densities in 1-week-ahead forecasts. Even the best fit models exhibited large variance in forecasted cyanobacterial densities and did not capture rare peak occurrences, indicating that significant explanatory variables when fitting models to historical data are not always effective for forecasting. Uncertainty partitioning revealed that model process specification and initial conditions dominated forecast uncertainty. These findings indicate that long-term studies of different cyanobacterial life stages and movement in the water column as well as measurements of drivers relevant to different life stages could improve model process representation of cyanobacteria abundance. In addition, improved observation protocols could better define initial conditions and reduce spatial misalignment of environmental data and cyanobacteria observations. Our results emphasize the importance of ecological forecasting principles and uncertainty partitioning to refine and understand predictive capacity across ecosystems.

Increased adoption of best practices in ecological forecasting enables comparisons of forecastability.

Near-term iterative forecasting is a powerful tool for ecological decision support and has the potential to transform our understanding of ecological predictability. However, to this point, there has been no cross-ecosystem analysis of near-term ecological forecasts, making it difficult to synthesize diverse research efforts and prioritize future developments for this emerging field. In this study, we analyzed 178 near-term (≤10-yr forecast horizon) ecological forecasting papers to understand the development and current state of near-term ecological forecasting literature and to compare forecast accuracy across scales and variables. Our results indicated that near-term ecological forecasting is widespread and growing: forecasts have been produced for sites on all seven continents and the rate of forecast publication is increasing over time. As forecast production has accelerated, some best practices have been proposed and application of these best practices is increasing. In particular, data publication, forecast archiving, and workflow automation have all increased significantly over time. However, adoption of proposed best practices remains low overall: for example, despite the fact that uncertainty is often cited as an essential component of an ecological forecast, only 45% of papers included uncertainty in their forecast outputs. As the use of these proposed best practices increases, near-term ecological forecasting has the potential to make significant contributions to our understanding of forecastability across scales and variables. In this study, we found that forecastability (defined here as realized forecast accuracy) decreased in predictable patterns over 1-7 d forecast horizons. Variables that were closely related (i.e., chlorophyll and phytoplankton) displayed very similar trends in forecastability, while more distantly related variables (i.e., pollen and evapotranspiration) exhibited significantly different patterns. Increasing use of proposed best practices in ecological forecasting will allow us to examine the forecastability of additional variables and timescales in the future, providing a robust analysis of the fundamental predictability of ecological variables.

Iron and manganese fluxes across the sediment-water interface in a drinking water reservoir.

The development of low dissolved oxygen (DO) concentrations in the hypolimnion of drinking water reservoirs during thermal stratification can lead to the reduction of oxidized, insoluble iron (Fe) and manganese (Mn) in sediments to soluble forms, which are then released into the water column. As metals degrade drinking water quality, robust measurements of metal fluxes under changing oxygen conditions are critical for optimizing water treatment. In this study, we conducted benthic flux chamber experiments in summer 2018 to directly quantify Fe and Mn fluxes at the sediment-water interface under different DO and redox conditions of a eutrophic drinking water reservoir with an oxygenation system (Falling Creek Reservoir, Vinton, VA, USA). Throughout the experiments, we monitored DO, oxidation-reduction potential (ORP), water temperature, and pH in the chambers and compared the metal fluxes in the chambers with time-series of fluxes calculated using a hypolimnetic mass balance method. Our results showed that metal fluxes were highly variable during the monitoring period and were sensitive to redox conditions in the water column at the sediment-water interface. The time-series changes in fluxes and relationship to redox conditions are suggestive of "hot moments", short time periods of intense biogeochemical cycling. Although the metal concentrations and fluxes are specific to this site, the approaches for examining relationships between metals, oxygen concentrations and overall redox conditions can be applied by water utilities to improve water quality management of Fe and Mn.

Metalimnetic oxygen minima alter the vertical profiles of carbon dioxide and methane in a managed freshwater reservoir.

Metalimnetic oxygen minimum zones (MOMs) commonly develop during the summer stratified period in freshwater reservoirs because of both natural processes and water quality management. While several previous studies have examined the causes of MOMs, much less is known about their effects, especially on reservoir biogeochemistry. MOMs create distinct redox gradients in the water column which may alter the magnitude and vertical distribution of dissolved methane (CH4) and carbon dioxide (CO2). The vertical distribution and diffusive efflux of CH4 and CO2 was monitored for two consecutive open-water seasons in a eutrophic reservoir that develops MOMs as a result of the operation of water quality engineering systems. During both summers, elevated concentrations of CH4 accumulated within the anoxic MOM, reaching a maximum of 120 µM, and elevated concentrations of CO2 accumulated in the oxic hypolimnion, reaching a maximum of 780 µM. Interestingly, the largest observed diffusive CH4 effluxes occurred before fall turnover in both years, while peak diffusive CO2 effluxes occurred both before and during turnover. Our data indicate that MOMs can substantially change the vertical distribution of CH4 and CO2 in the water column in reservoirs, resulting in the accumulation of CH4 in the metalimnion (vs. at the sediments) and CO2 in the hypolimnion.

Chaoborus spp. Transport CH4 from the Sediments to the Surface Waters of a Eutrophic Reservoir, But Their Contribution to Water Column CH4 Concentrations and Diffusive Efflux Is Minor.

Chaoborus spp. (midge larvae) live in the anoxic sediments and hypolimnia of freshwater lakes and reservoirs during the day and migrate to the surface waters at night to feed on plankton. It has recently been proposed that Chaoborus take up methane (CH4) from the sediments in their tracheal gas sacs, use this acquired buoyancy to ascend into the surface waters, and then release the CH4, thereby serving as a CH4 "pump" to the atmosphere. We tested this hypothesis using diel surveys and seasonal monitoring, as well as incubations of Chaoborus to measure CH4 transport in their gas sacs at different depths and times in a eutrophic reservoir. We found that Chaoborus transported CH4 from the hypolimnion to the lower epilimnion at dusk, but the overall rate of CH4 transport was minor, and incubations revealed substantial variability in CH4 transport over space and time. We calculated that Chaoborus transport ∼0.1 mmol CH4 m-2 yr-1 to the epilimnion in our study reservoir, a very low proportion (<1%) of total CH4 diffusive flux during the summer stratified period. Our data further indicate that CH4 transport by Chaoborus is sensitive to water column mixing, Chaoborus density, and Chaoborus species identity.