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.

Water level drawdowns are increasingly common in lakes and reservoirs worldwide as a result of both climate change and water management. Drawdowns can have direct effects on physical properties of a waterbody (e.g., by altering stratification and light dynamics), which can interact to modify the waterbody’s biology and chemistry. However, the ecosystem-level effects of drawdown remain poorly characterized in small, thermally-stratified reservoirs, which are common in many regions of the world. Here, we intensively monitored a small eutrophic reservoir for two years, including before, during, and after a month-long drawdown that reduced total reservoir volume by 36%. During drawdown, stratification strength (maximum buoyancy frequency) and surface phosphate concentrations both increased, contributing to a substantial surface phytoplankton bloom. The peak in phytoplankton biomass was followed by cascading changes in surface water chemistry associated with bloom degradation, with sequential peaks in dissolved organic carbon, dissolved carbon dioxide, and ammonium concentrations that were up to an order of magnitude higher than the previous year. Dissolved oxygen concentrations substantially decreased in the surface waters during drawdown (to 41% saturation), which was associated with increased total iron and manganese concentrations. Combined, our results illustrate how changes in water level can have cascading effects on coupled physical, chemical, and biological processes. As climate change and water management continue to increase the frequency of drawdowns in lakes worldwide, our results highlight the importance of characterizing how water level variability can alter complex in-lake ecosystem processes, thereby affecting water quality.

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.

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.

The biogeochemical cycles of iron (Fe) and manganese (Mn) in lakes and reservoirs have predictable seasonal trends, largely governed by stratification dynamics and redox conditions in the hypolimnion. However, short-term (i.e., sub-weekly) trends in Fe and Mn cycling are less well-understood, as most monitoring efforts focus on longer-term (i.e., monthly to yearly) time scales. The potential for elevated Fe and Mn to degrade water quality and impact ecosystem functioning, coupled with increasing evidence for high spatiotemporal variability in other biogeochemical cycles, necessitates a closer evaluation of the short-term Fe and Mn cycling dynamics in lakes and reservoirs. We adapted a UV-visible spectrophotometer coupled with a multiplexor pumping system and PLSR modeling to generate high spatiotemporal resolution predictions of Fe and Mn concentrations in a drinking water reservoir (Falling Creek Reservoir, Vinton, VA, USA) equipped with a hypolimnetic oxygenation (HOx) system. We quantified hourly Fe and Mn concentrations during two distinct transitional periods: reservoir turnover (Fall 2020) and initiation of the HOx system (Summer 2021). Our sensor system was able to successfully predict mean Fe and Mn concentrations as well as capture sub-weekly variability, ground-truthed by traditional grab sampling and laboratory analysis. During fall turnover, hypolimnetic Fe and Mn concentrations began to decrease more than two weeks before complete mixing of the reservoir occurred, with rapid equalization of epilimnetic and hypolimnetic Fe and Mn concentrations in less than 48 hours after full water column mixing. During the initiation of hypolimnetic oxygenation in Summer 2021, we observed that Fe and Mn were similarly affected by physical mixing in the hypolimnion, but displayed distinctly different responses to oxygenation, as indicated by the rapid oxidation of soluble Fe but not soluble Mn. This study demonstrates that Fe and Mn concentrations are highly sensitive to shifting DO and stratification and that their dynamics can substantially change on hourly to daily time scales in response to these transitions.