Explore the vast range of titles available.

The concentration of oxygen is fundamental to lake water quality and ecosystem functioning through its control over habitat availability for organisms, redox reactions, and recycling of organic material. In many eutrophic lakes, oxygen depletion in the bottom layer (hypolimnion) occurs annually during summer stratification. The temporal and spatial extent of summer hypolimnetic anoxia is determined by interactions between the lake and its external drivers (e.g., catchment characteristics, nutrient loads, meteorology) as well as internal feedback mechanisms (e.g., organic matter recycling, phytoplankton blooms). How these drivers interact to control the evolution of lake anoxia over decadal timescales will determine, in part, the future lake water quality. In this study, we used a vertical one-dimensional hydrodynamic–ecological model (GLM-AED2) coupled with a calibrated hydrological catchment model (PIHM-Lake) to simulate the thermal and water quality dynamics of the eutrophic Lake Mendota (USA) over a 37 year period. The calibration and validation of the lake model consisted of a global sensitivity evaluation as well as the application of an optimization algorithm to improve the fit between observed and simulated data. We calculated stability indices (Schmidt stability, Birgean work, stored internal heat), identified spring mixing and summer stratification periods, and quantified the energy required for stratification and mixing. To qualify which external and internal factors were most important in driving the interannual variation in summer anoxia, we applied a random-forest classifier and multiple linear regressions to modeled ecosystem variables (e.g., stratification onset and offset, ice duration, gross primary production). Lake Mendota exhibited prolonged hypolimnetic anoxia each summer, lasting between 50–60 d. The summer heat budget, the timing of thermal stratification, and the gross primary production in the epilimnion prior to summer stratification were the most important predictors of the spatial and temporal extent of summer anoxia periods in Lake Mendota. Interannual variability in anoxia was largely driven by physical factors: earlier onset of thermal stratification in combination with a higher vertical stability strongly affected the duration and spatial extent of summer anoxia. A measured step change upward in summer anoxia in 2010 was unexplained by the GLM-AED2 model. Although the cause remains unknown, possible factors include invasion by the predacious zooplankton Bythotrephes longimanus. As the heat budget depended primarily on external meteorological conditions, the spatial and temporal extent of summer anoxia in Lake Mendota is likely to increase in the near future as a result of projected climate change in the region.

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 modeling 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. We found 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: (a) no individual model performed best at every depth and horizon (days in the future), and (b) 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 utilizing MMEs, rather than individual models, in operational forecasts. Key Points Aggregated lake temperature forecast skill was higher for multi‐model ensemble (MME) forecasts than individual model forecasts Including baseline empirical models (day‐of‐year, persistence) with process models improved MME forecast performance MME forecasts improved forecast skill by “hedging,” as no individual model performed best at all horizons or depths

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.

The quality of lake ice is of uppermost importance for ice safety and under-ice ecology, but its temporal and spatial variability is largely unknown. Here we conducted a coordinated lake ice quality sampling campaign across the Northern Hemisphere during one of the warmest winters since 1880 and show that lake ice during 2020/2021 commonly consisted of unstable white ice, at times contributing up to 100% to the total ice thickness. We observed that white ice increased over the winter season, becoming thickest and constituting the largest proportion of the ice layer towards the end of the ice cover season when fatal winter drownings occur most often and light limits the growth and reproduction of primary producers. We attribute the dominance of white ice before ice-off to air temperatures varying around the freezing point, a condition which occurs more frequently during warmer winters. Thus, under continued global warming, the prevalence of white ice is likely to substantially increase during the critical period before ice-off, for which we adjusted commonly used equations for human ice safety and light transmittance through ice. Under continued global warming, lakes will increasingly be covered by white ice, in particular towards the end of the ice cover season when fatal winter drownings occur most often and light limits the growth and reproduction of primary producers.

Despite the growing use of Aquatic Ecosystem Models for lake modeling, there is currently no widely applicable framework for their configuration, calibration, and evaluation. Calibration is generally based on direct data comparison of observed versus modeled state variables using standard statistical techniques, however, this approach may not give a complete picture of the model's ability to capture system‐scale behavior that is not easily perceivable in observations, but which may be important for resource management. The aim of this study is to compare the performance of “naïve” calibration and a “system‐inspired” calibration, an approach that augments the standard state‐based calibration with a range of system‐inspired metrics (e.g., thermocline depth, metalimnetic oxygen minima), to increase the coherence between the simulated and natural ecosystems. A coupled physical‐biogeochemical model was applied to a focal site to simulate two key state‐variables: water temperature and dissolved oxygen. The model was calibrated according to the new system‐inspired modeling convention, using formal calibration techniques. There was an improvement in the simulation using parameters optimized on the additional metrics, which helped to reduce uncertainty predicting aspects of the system relevant to reservoir management, such as the occurrence of the metalimnetic oxygen minima. Extending the use of system‐inspired metrics when calibrating models has the potential to improve model fidelity for capturing more complex ecosystem dynamics. Key Points We assessed the use of system‐inspired metrics in a novel approach to calibrating Aquatic Ecosystem Models (AEMs) The use of system‐inspired metrics in calibration improved model performance relative to traditional calibration methods Implementation of system‐inspired metrics has the potential to greatly improve model prediction of complex ecosystem dynamics

Here we draw attention to the potential for pelagic bloom-forming cyanobacteria to have substantial effects on nutrient cycling and ecosystem resilience across a wide range of lakes. Specifically, we hypothesize that cyanobacterial blooms can influence lake nutrient cycling, resilience, and regime shifts by tapping into pools of nitrogen (N) and phosphorus (P) not usually accessible to phytoplankton. The ability of many cyanobacterial taxa to fix dissolved N 2 gas is a well-known potential source of N, but some taxa can also access pools of P in sediments and bottom waters. Both of these nutrients can be released to the water column via leakage or mortality, thereby increasing nutrient availability for other phytoplankton and microbes. Moreover, cyanobacterial blooms are not restricted to high nutrient (eutrophic) lakes: blooms also occur in lakes with low nutrient concentrations, suggesting that changes in nutrient cycling and ecosystem resilience mediated by cyanobacteria could affect lakes across a gradient of nutrient concentrations. We used a simple model of coupled N and P cycles to explore the effects of cyanobacteria on nutrient dynamics and resilience. Consistent with our hypothesis, parameters reflecting cyanobacterial modification of N and P cycling alter the number, location, and/or stability of model equilibria. In particular, the model demonstrates that blooms of cyanobacteria in low-nutrient conditions can facilitate a shift to the high-nutrient state by reducing the resilience of the low-nutrient state. This suggests that cyanobacterial blooms warrant attention as potential drivers of the transition from a low-nutrient, clear-water regime to a high-nutrient, turbid-water regime, a prediction of particular concern given that such blooms are reported to be increasing in many regions of the world due in part to global climate change.

Environmental research requires understanding nonlinear ecological dynamics that interact across multiple spatial and temporal scales. The analysis of long‐term and high‐frequency sensor data combined with simulation modeling enables interpretation of complex ecological phenomena, and the computational skills needed to conduct these analyses are increasingly being integrated into graduate student training programs in ecology. Despite its importance, however, computational literacy—that is, the ability to harness the power of computer technologies to accomplish tasks—is rarely taught in undergraduate ecology classrooms, representing a major gap in training students to tackle complex environmental challenges. Through our experience developing undergraduate curricula in long‐term and high‐frequency data analysis and simulation modeling for two environmental science pedagogical initiatives, Project EDDIE (Environmental Data‐Driven Inquiry and Exploration) and Macrosystems EDDIE, we have found that students often feel intimidated by computational tasks, which is compounded by the lack of familiarity with software (e.g., R) and the steep learning curves associated with script‐based analytical tools. The use of prepackaged, flexible modules that introduce programming as a mechanism to explore environmental datasets and teach inquiry‐based ecology, such as those developed for Project EDDIE and Macrosystems EDDIE, can significantly increase students’ experience and comfort levels with advanced computational tools. These types of modules in turn provide great potential for empowering students with the computational literacy needed to ask ecological questions and test hypotheses on their own. As continental‐scale sensor observatory networks rapidly expand the availability of long‐term and high‐frequency data, students with the skills to manipulate, visualize, and interpret such data will be well‐prepared for diverse careers in data science, and will help advance the future of open, reproducible science in ecology. Ecology increasingly requires the use of advanced computational techniques including simulation modeling and the analysis of high‐frequency data, yet teaching these techniques in undergraduate ecology courses is not common. Here, we synthesize our experiences developing and assessing teaching modules that bring “big data” into undergraduate ecology courses. These experiences have emphasized the importance, challenges, and opportunities of explicitly embedding computational literacy skills (e.g., data manipulation and visualization using the R programming language) into undergraduate ecology curricula, and we use data from student assessments to support our observations.

Globally significant quantities of carbon (C), nitrogen (N), and phosphorus (P) enter freshwater reservoirs each year. These inputs can be buried in sediments, respired, taken up by organisms, emitted to the atmosphere, or exported downstream. While much is known about reservoir-scale biogeochemical processing, less is known about spatial and temporal variability of biogeochemistry within a reservoir along the continuum from inflowing streams to the dam. To address this gap, we examined longitudinal variability in surface water biogeochemistry (C, N, and P) in two small reservoirs throughout a thermally stratified season. We sampled total and dissolved fractions of C, N, and P, as well as chlorophyll-a from each reservoir’s major inflows to the dam. We found that heterogeneity in biogeochemical concentrations was greater over time than space. However, dissolved nutrient and organic carbon concentrations had high site-to-site variability within both reservoirs, potentially as a result of shifting biological activity or environmental conditions. When considering spatially explicit processing, we found that certain locations within the reservoir, most often the stream–reservoir interface, acted as “hotspots” of change in biogeochemical concentrations. Our study suggests that spatially explicit metrics of biogeochemical processing could help constrain the role of reservoirs in C, N, and P cycles in the landscape. Ultimately, our results highlight that biogeochemical heterogeneity in small reservoirs may be more variable over time than space, and that some sites within reservoirs play critically important roles in whole-ecosystem biogeochemical processing.

Arboreal herbivory is rare among mammals. The few species with this lifestyle possess unique adaptions to overcome size-related constraints on nutritional energetics. Sloths are folivores that spend most of their time resting or eating in the forest canopy. A three-toed sloth will, however, descend its tree weekly to defecate, which is risky, energetically costly and, until now, inexplicable. We hypothesized that this behaviour sustains an ecosystem in the fur of sloths, which confers cryptic nutritional benefits to sloths. We found that the more specialized three-toed sloths harboured more phoretic moths, greater concentrations of inorganic nitrogen and higher algal biomass than the generalist two-toed sloths. Moth density was positively related to inorganic nitrogen concentration and algal biomass in the fur. We discovered that sloths consumed algae from their fur, which was highly digestible and lipid-rich. By descending a tree to defecate, sloths transport moths to their oviposition sites in sloth dung, which facilitates moth colonization of sloth fur. Moths are portals for nutrients, increasing nitrogen levels in sloth fur, which fuels algal growth. Sloths consume these algae-gardens, presumably to augment their limited diet. These linked mutualisms between moths, sloths and algae appear to aid the sloth in overcoming a highly constrained lifestyle.

Ecosystems around the globe are experiencing changes in both the magnitude and fluctuations of environmental conditions 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 Forecasting Lake and Reservoir Ecosystems (FLARE), 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.81°C, mean 7‐day RMSE of 1.15°C, and mean 35‐day RMSE of 1.94°C. Aggregated across the year, daily data assimilation yielded the most skillful forecasts at 1‐ to 7‐day‐ahead horizons, but weekly data assimilation resulted in the most skillful forecasts at 8‐ to 35‐day‐ahead horizons. Within a year, forecasts with weekly data assimilation consistently outperformed forecasts with daily data assimilation after the 8‐day forecast horizon during mixed spring/autumn periods and 5‐ to 14‐day‐ahead horizons during the summer‐stratified period, depending on depth. 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.