Water temperature forecasting in lakes and reservoirs is a valuable tool to manage crucial freshwater resources in a changing and more variable climate, but previous efforts have yet to identify an optimal modelling approach. Here, we demonstrate the first multi-model ensemble (MME) reservoir water temperature forecast, a forecasting method that combines individual model strengths in a single forecasting framework. We developed two MMEs: a three-model process-based MME and a five-model MME that includes process-based and empirical models to forecast water temperature profiles at a temperate drinking water reservoir. Our results showed that the five-model MME improved forecast performance by 8-30% relative to individual models and the process-based MME, as quantified using an aggregated probabilistic skill score. This increase in performance was due to large improvements in forecast bias in the five-model MME, despite increases in forecast uncertainty. High correlation among the process-based models resulted in little improvement in forecast performance in the process-based MME relative to the individual process-based models. The utility of MMEs is highlighted by two results: 1) no individual model performed best at every depth and horizon (days in the future), and 2) MMEs avoided poor performances by rarely producing the worst forecast for any single forecasted period (<6% of the worst ranked forecasts over time). This work presents an example of how existing models can be combined to improve water temperature forecasting in lakes and reservoirs and discusses the value of utilising MMEs, rather than individual models, in operational forecasts.

Small freshwater reservoirs are ubiquitous and likely play an important role in global greenhouse gas (GHG) budgets relative to their limited water surface area. However, constraining annual GHG fluxes in small freshwater reservoirs is challenging given their footprint area and spatially and temporally variable emissions. To quantify the GHG budget of a small (0.1 km2) reservoir, we deployed an eddy covariance system in a small reservoir located in southwestern Virginia, USA over two years to measure carbon dioxide (CO2) and methane (CH4) fluxes near-continuously. Fluxes were coupled with in situ sensors measuring multiple environmental parameters. Over both years, we found the reservoir to be a large source of CO2 (633-731 g CO2-C m-2 yr-1) and CH4 (1.02-1.29 g CH4-C m-2 yr-1) to the atmosphere, with substantial sub-daily, daily, weekly, and seasonal timescales of variability. For example, fluxes were substantially greater during the summer thermally-stratified season as compared to the winter. In addition, we observed significantly greater GHG fluxes during winter intermittent ice-on conditions as compared to continuous ice-on conditions, suggesting GHG emissions from lakes and reservoirs may increase with predicted decreases in winter ice-cover. Finally, we identified several key environmental variables that may be driving reservoir GHG fluxes at multiple timescales, including, surface water temperature and thermocline depth followed by fluorescent dissolved organic matter. Overall, our novel year-round eddy covariance data from a small reservoir indicate that these freshwater ecosystems likely contribute a substantial amount of CO2 and CH4 to global GHG budgets, relative to their surface area.

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). To assess the current state of freshwater forecasting and identify opportunities for future progress, we synthesized freshwater forecasting papers published in the past five 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 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 five 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 necessitates 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.

TITLE : Embedding communication concepts in forecasting training increases students’ understanding of ecological uncertainty Submitted as an Article to Ecosphere , Eco-Education TrackAUTHOR LIST: Whitney M. Woelmera*, Tadhg N. Moorea,b11Present address: School of Biological, Earth and Environmental Sciences, University College Cork, Cork, Ireland, Mary E. Loftona, R. Quinn Thomasa,b, and Cayelan C. CareyaaDepartment of Biological Sciences, Virginia Tech, Blacksburg, VA, USAbDepartment of Forest Resources and Environmental Conservation, Virginia Tech, Blacksburg, VA, USA*Corresponding author: [email protected] RESEARCH STATEMENT : This study collected and analyzed human subject data and was approved by the Virginia Tech Institutional Review Board (19-669) and the Carleton College Institutional Review Board (19-20 065). Data for this study have been anonymized and aggregated and can be found at Woelmer (2023) along with all code to reproduce the analysis and figures within this study.Woelmer, W. 2023. Wwoelmer/module8_public_ecosphere: Ecosphere submission March 2023 (v1.0). Zenodo. https://doi.org/10.5281/zenodo.7733965KEYWORDS : active learning, ecology education, ecological forecast, Macrosystems EDDIE, R Shiny, teaching modules, translational ecology, undergraduate curricula, visualization literacyABSTRACT : Communicating and interpreting uncertainty in ecological model predictions is notoriously challenging, motivating the need for new educational tools which introduce ecology students to core concepts in uncertainty communication. Ecological forecasting, an emerging approach to estimate future states of ecological systems with uncertainty, provides a relevant and engaging framework for introducing uncertainty communication to undergraduate students, as forecasts can be used as decision support tools for addressing real-world ecological problems and are inherently uncertain. To provide critical training on uncertainty communication and introduce undergraduate students to the use of ecological forecasts for guiding decision-making, we developed a hands-on teaching module within the Macrosystems EDDIE (Environmental Data-Driven Inquiry and Exploration; MacrosystemsEDDIE.org) educational program. Our module used an active learning approach by embedding forecasting activities in an R Shiny application to engage introductory students in data science, ecological modeling, and forecasting without needing advanced computational or programming skills. Pre- and post-module assessment data from >250 undergraduate ecology students indicate that the module significantly increased students’ ability to interpret forecast visualizations with uncertainty, identify different ways to communicate forecast uncertainty for diverse users, and correctly define ecological forecasting terms. Specifically, students were more likely to describe visual, numeric, and probabilistic methods of uncertainty communication following module completion. Students were also able to identify more benefits of ecological forecasting following module completion, with the key benefits of using forecasts for prediction and decision-making most commonly described. These results show promise for introducing ecological model uncertainty, data visualizations, and forecasting into undergraduate ecology curricula via software-based learning, which can increase students’ ability to engage and understand complex ecological concepts.

A document by R. Quinn Thomas. Click on the document to view its contents.

Mary Lofton1,2a, Tadhg Moore1b, Whitney Woelmer1c, Quinn Thomas1,2,3d, Cayelan Carey1,2e1Department of Biological Sciences, Virginia Tech. 926 West Campus Drive, Blacksburg, Virginia, 24061, United States2Center for Ecosystem Forecasting, Virginia Tech. 1015 Life Science Circle, Blacksburg, VA 24061, United States3Department of Forest Resources and Environmental Conservation, Virginia Tech. 310 West Campus Drive, Blacksburg, Virginia, 24061, United Statesa Postdoctoral Research Associate; [email protected] Postdoctoral Research Associate; [email protected] Ph.D. Candidate; [email protected] Associate Professor; [email protected] Professor; [email protected]

Ecosystems around the globe are experiencing increased variability due to land use and climate change. In response, ecologists are increasingly using near-term, iterative ecological forecasts to predict how ecosystems will change in the future. To date, many near-term, iterative forecasting systems have been developed using high temporal frequency (minute to hourly resolution) data streams for assimilation. However, this approach may be cost-prohibitive or impossible for forecasting ecological variables that lack high-frequency sensors or have high data latency (i.e., a delay before data are available for modeling after collection). To explore the effects of data assimilation frequency on forecast skill, we developed water temperature forecasts for a eutrophic drinking water reservoir and conducted data assimilation experiments by selectively withholding observations to examine the effect of data availability on forecast accuracy. We used  in-situ sensors, manually collected data, and a calibrated water quality ecosystem model driven by forecasted weather data to generate future water temperature forecasts using FLARE (Forecasting Lake And Reservoir Ecosystems), an open-source water quality forecasting system. We tested the effect of daily, weekly, fortnightly, and monthly data assimilation on the skill of 1 to 35-day-ahead water temperature forecasts. We found that forecast skill varied depending on the season, forecast horizon, depth, and data assimilation frequency, but overall forecast performance was high, with a mean 1-day-ahead forecast root mean square error (RMSE) of 0.94°C, mean 7-day RMSE of 1.33°C, and mean 35-day RMSE of 2.15°C. Aggregated across the year, daily data assimilation yielded the most skillful forecasts at 1-7-day-ahead horizons, weekly data assimilation resulted in the most skillful forecasts at 8-35-day-ahead horizons. Within a year, daily to fortnightly data assimilation substantially outperformed monthly data assimilation in the stratified summer period, whereas all data assimilation frequencies resulted in skillful forecasts across depths in the mixed spring/autumn periods for shorter forecast horizons. Our results suggest that lower-frequency data (i.e., weekly) may be adequate for developing accurate forecasts in some applications, further enabling the development of forecasts broadly across ecosystems and ecological variables without high-frequency sensor data.

Freshwater ecosystems provide vital services, yet are facing increasing risks from global change. In particular, lake thermal dynamics have been altered around the world as a result of climate change, necessitating a predictive understanding of how climate will continue to alter lakes in the future as well as the associated uncertainty in these predictions. Numerous sources of uncertainty affect projections of future lake conditions but few are quantified, limiting the use of lake modeling projections as management tools. To quantify and evaluate the effects of two potentially important sources of uncertainty, lake model selection uncertainty and climate model selection uncertainty, we developed ensemble projections of lake thermal dynamics for a dimictic lake in New Hampshire, USA (Lake Sunapee). Our ensemble projections used four different climate models as inputs to five vertical one-dimensional (1-D) hydrodynamic lake models under three different climate change scenarios to simulate thermal metrics from 2006 to 2099. We found that almost all the lake thermal metrics modeled (surface water temperature, bottom water temperature, Schmidt stability, stratification duration, and ice cover, but not thermocline depth) are projected to change over the next century. Importantly, we found that the dominant source of uncertainty varied among the thermal metrics, as thermal metrics associated with the surface waters (surface water temperature, total ice duration) were driven primarily by climate model selection uncertainty, while metrics associated with deeper depths (bottom water temperature, stratification duration) were dominated by lake model selection uncertainty. Consequently, our results indicate that researchers generating projections of lake bottom water metrics should prioritize including multiple lake models for best capturing projection uncertainty, while those focusing on lake surface metrics should prioritize including multiple climate models. Overall, our ensemble modeling study reveals important information on how climate change will affect lake thermal properties, and also provides some of the first analyses on how climate model selection uncertainty and lake model selection uncertainty interact to affect projections of future lake dynamics.

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, seven-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.

The National Ecological Observatory Network (NEON)’s standardized monitoring program provides an unprecedented opportunity for comparing the predictability of ecosystems. To harness the power of NEON data for examining environmental predictability, we scaled a near-term, iterative water temperature forecasting system to all six conterminous NEON lakes. We generated 1 to 35-day ahead forecasts using a process-based hydrodynamic model that was updated with observations as they became available. Forecasts were more accurate than a null model up to 35-days ahead among lakes, with an aggregated 1-day ahead RMSE (root-mean square error) of 0.60℃ and 35-days ahead RMSE of 2.17℃. Water temperature forecast accuracy was positively associated with lake depth and water clarity, and negatively associated with catchment size and fetch. Our results suggest that lake characteristics interact with weather to control the predictability of thermal structure. Our work provides some of the first probabilistic forecasts of NEON sites and a framework for examining continental-scale 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 to 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 which 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.