skip to main content


Title: Lake thermal structure drives interannual variability in summer anoxia dynamics in a eutrophic lake over 37 years
Abstract. 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.  more » « less
Award ID(s):
1517823 1753639 1759865 1753657 2025982
NSF-PAR ID:
10220713
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Hydrology and Earth System Sciences
Volume:
25
Issue:
2
ISSN:
1607-7938
Page Range / eLocation ID:
1009 to 1032
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Dataset to run a 37-year simulation (1979-2015) of the Lake Mendota lake ecosystem using the vertical 1D GLM-AED2 model. The focus of this modeling study is on determining the drivers of year-to-year variability in the spatial and temporal extent of hypolimnetic anoxia. 
    more » « less
  2. Abstract. Hypolimnetic oxygen depletion during summer stratification in lakes can lead to hypoxic and anoxic conditions. Hypolimnetic anoxia is a water quality issue with many consequences, including reduced habitat for cold-water fish species, reduced quality of drinking water, and increased nutrient and organic carbon (OC) release from sediments. Both allochthonous and autochthonous OC loads contribute to oxygen depletion by providing substrate for microbial respiration; however, their relative contributions to oxygen depletion across diverse lake systems remain uncertain. Lake characteristics, such as trophic state, hydrology, and morphometry, are also influential in carbon-cycling processes and may impact oxygen depletion dynamics. To investigate the effects of carbon cycling on hypolimnetic oxygen depletion, we used a two-layer process-based lake model to simulate daily metabolism dynamics for six Wisconsin lakes over 20 years (1995–2014). Physical processes and internal metabolic processes were included in the model and were used to predict dissolved oxygen (DO), particulate OC (POC), and dissolved OC (DOC). In our study of oligotrophic, mesotrophic, and eutrophic lakes, we found autochthony to be far more important than allochthony to hypolimnetic oxygen depletion. Autochthonous POC respiration in the water column contributed the most towards hypolimnetic oxygen depletion in the eutrophic study lakes. POC water column respiration and sediment respiration had similar contributions in the mesotrophic and oligotrophic study lakes. Differences in terms of source of respiration are discussed with consideration of lake productivity and the processing and fates of organic carbon loads.

     
    more » « less
  3. Abstract

    Lake water clarity, phytoplankton biomass, and hypolimnetic oxygen concentration are metrics of water quality that are highly degraded in eutrophic systems. Eutrophication is linked to legacy nutrients stored in catchment soils and in lake sediments. Long lags in water quality improvement under scenarios of nutrient load reduction to lakes indicate an apparent ecosystem memory tied to the interactions between water biogeochemistry and lake sediment nutrients. To investigate how nutrient legacies and ecosystem memory control lake water quality dynamics, we coupled nutrient cycling and lake metabolism in a model to recreate long‐term water quality of a eutrophic lake (Lake Mendota, Wisconsin, USA). We modeled long‐term recovery of water quality under scenarios of nutrient load reduction and found that the rates and patterns of water quality improvement depended on changes in phosphorus (P) and organic carbon storage in the water column and sediments. Through scenarios of water quality improvement, we showed that water quality variables have distinct phases of change determined by the turnover rates of storage pools—an initial and rapid water quality improvement due to water column flushing, followed by a much longer and slower improvement as sediment P pools were slowly reduced. Water clarity, phytoplankton biomass, and hypolimnetic dissolved oxygen differed in their time responses. Water clarity and algal biomass improved within years of nutrient reductions, but hypolimnetic oxygen took decades to improve. Even with reduced catchment loading, recovery of Lake Mendota to a mesotrophic state may require decades due to nutrient legacies and long ecosystem memory.

     
    more » « less
  4. Green Lake is the deepest natural inland lake in Wisconsin, with a maximum depth of about 72 meters. In the early 1900s, the lake was believed to have very good water quality (low nutrient concentrations and good water clarity) with low dissolved oxygen (DO) concentrations occurring in only the deepest part of the lake. Because of increased phosphorus (P) inputs from anthropogenic activities in its watershed, total phosphorus (TP) concentrations in the lake have increased; these changes have led to increased algal production and low DO concentrations not only in the deepest areas but also in the middle of the water column (metalimnion). The U.S. Geological Survey has routinely monitored the lake since 2004 and its tributaries since 1988. Results from this monitoring led the Wisconsin Department of Natural Resources (WDNR) to list the lake as impaired because of low DO concentrations in the metalimnion, and they identified elevated TP concentrations as the cause of impairment. As part of this study by the U.S. Geological Survey, in cooperation with the Green Lake Sanitary District, the lake and its tributaries were comprehensively sampled in 2017–18 to augment ongoing monitoring that would further describe the low DO concentrations in the lake (especially in the metalimnion). Empirical and process-driven water-quality models were then used to determine the causes of the low DO concentrations and the magnitudes of P-load reductions needed to improve the water quality of the lake enough to meet multiple water-quality goals, including the WDNR’s criteria for TP and DO. Data from previous studies showed that DO concentrations in the metalimnion decreased slightly as summer progressed in the early 1900s but, since the late 1970s, have typically dropped below 5 milligrams per liter (mg/L), which is the WDNR criterion for impairment. During 2014–18 (the baseline period for this study), the near-surface geometric mean TP concentration during June–September in the east side of the lake was 0.020 mg/L and in the west side was 0.016 mg/L (both were above the 0.015-mg/L WDNR criterion for the lake), and the metalimnetic DO minimum concentrations (MOMs) measured in August ranged from 1.0 to 4.7 mg/L. The degradation in water quality was assumed to have been caused by excessive P inputs to the lake; therefore, the TP inputs to the lake were estimated. The mean annual external P load during 2014–18 was estimated to be 8,980 kilograms per year (kg/yr), of which monitored and unmonitored tributary inputs contributed 84 percent, atmospheric inputs contributed 8 percent, waterfowl contributed 7 percent, and septic systems contributed 1 percent. During fall turnover, internal sediment recycling contributed an additional 7,040 kilograms that increased TP concentrations in shallow areas of the lake by about 0.020 mg/L. The elevated TP concentrations then persisted until the following spring. On an annual basis, however, there was a net deposition of P to the bottom sediments. Empirical models were used to describe how the near-surface water quality of Green Lake would be expected to respond to changes in external P loading. Predictions from the models showed a relatively linear response between P loading and TP and chlorophyll-a (Chl-a) concentrations in the lake, with the changes in TP and Chl-a concentrations being less on a percentage basis (50–60 percent for TP and 30–70 percent for Chl-a) than the changes in P loading. Mean summer water clarity, quantified by Secchi disk depths, had a greater response to decreases in P loading than to increases in P loading. Based on these relations, external P loading to the lake would need to be decreased from 8,980 kg/yr to about 5,460 kg/yr for the geometric mean June–September TP concentration in the east side of the lake, with higher TP concentrations than in the west side, to reach the WDNR criterion of 0.015 mg/L. This reduction of 3,520 kg/yr is equivalent to a 46-percent reduction in the potentially controllable external P sources (all external sources except for precipitation, atmospheric deposition, and waterfowl) from those measured during water years 2014–18. The total external P loading would need to decrease to 7,680 kg/yr (a 17-percent reduction in potentially controllable external P sources) for near-surface June–September TP concentrations in the west side of the lake to reach 0.015 mg/L. Total external P loading would need to decrease to 3,870–5,320 kg/yr for the lake to be classified as oligotrophic, with a near-surface June–September TP concentration of 0.012 mg/L. Results from the hydrodynamic water-quality model GLM–AED (General Lake Model coupled to the Aquatic Ecodynamics modeling library) indicated that MOMs are driven by external P loading and internal sediment recycling that lead to high TP concentrations during spring and early summer, which in turn lead to high phytoplankton production, high metabolism and respiration, and ultimately DO consumption in the upper, warmer areas of the metalimnion. GLM–AED results indicated that settling of organic material during summer might be slowed by the colder, denser, and more viscous water in the metalimnion and thus increase DO consumption. Based on empirical evidence from a comparison of MOMs with various meteorological, hydrologic, water quality, and in-lake physical factors, MOMs were lower during summers, when metalimnetic water temperatures were warmer, near-surface Chl-a and TP concentrations were higher, and Secchi depths were lower. GLM–AED results indicated that the external P load would need to be reduced to about 4,060 kg/yr, a 57-percent reduction from that measured in 2014–18, to eliminate the occurrence of MOMs less than 5 mg/L during more than 75 percent of the years (the target provided by the WDNR). Large reductions in external P loading are expected to have an immediate effect on the near-surface TP concentrations and metalimnetic DO concentrations in Green Lake; however, it may take several years for the full effects of the external-load reduction to be observed because internal sediment recycling is an important source of P for the following spring. 
    more » « less
  5. Abstract

    Absence of dissolved oxygen (anoxia) in the hypolimnion of lakes can eliminate habitat for sensitive species and may induce the release of sediment‐bound phosphorus. Lake anoxia generally results from decomposition of organic matter, which is exacerbated by high nutrient loads. Total phosphorus (TP) in lakes is regulated by static aspects of the lake’s watershed, but lake TP can be readily increased by human activities. In some low‐nutrient lakes, basin morphometry may induce naturally occurring anoxia. The occurrence of natural anoxia is especially important to consider in lake water quality assessments that compare observed conditions to expected reference conditions. To investigate the occurrence of natural vs. anthropogenically influenced anoxia, we constructed a logistic regression model to calculate the probability of low‐nutrient lakes (TP < 15 µg/L) developing aerial anoxic extent ≥10% by testing the predictive potential of variables related to basin morphometry, depths of lake thermal strata, epilimnetic TP, and dissolved organic carbon (DOC). Maximum lake depth and the proportion of lake area under the top of the metalimnion were the most important variables to predict the likelihood of hypolimnetic anoxia, which correctly predicted anoxic condition in 84% of lakes (Model 1). Adding TP as a third variable to Model 1 produced a significantly improved model (Model 2) but the prediction success rate was comparable (86%). We also present a model for lakes with limited bathymetric data, which predicts anoxia with 81% accuracy based on maximum lake depth and mean thermocline depth at peak stratification. DOC was relatively low (4.3 ± 1.5 mg/L [mean ± SD]) in the study lakes and its inclusion did not improve model performance. In Model 1, lakes with an anoxic extent ≥10% of lake area had significantly higher epilimnetic TP than lakes with oxic hypolimnia, regardless of prediction category or success. Our results indicate that including TP as a variable helps refine models based on morphometry alone, but lake morphometry and stratification dynamics are the most important factors in the development of anoxic extent in low‐nutrient temperate lakes. Our approach informs studies concerned with identifying key factors that influence regime shifts in a variety of ecosystems.

     
    more » « less