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. more » 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. « less
; ; ; ; ; ; ;
Award ID(s):
1517823 1753639 1759865 1753657 2025982
Publication Date:
Journal Name:
Hydrology and Earth System Sciences
Page Range or eLocation-ID:
1009 to 1032
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 themore »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.« less
  2. Changes in mixing regimes and CO2 availability may promote harmful cyanobacterial blooms in polymictic lakes and ponds globally, but the underlying mechanisms still remain unclear. We integrated results from a natural experiment comprising an average-wet year (2011) and one with heat waves (2012), a long-term meteorological dataset (1960–2010), historical phosphorus concentrations and corresponding sedimentary pigment records, to determine, on different temporal scales, the mechanistic controls of cyanobacterial blooms in a eutrophic polymictic lake. Intense warming in 2012 was associated with: 1) increased stability of the water column with buoyancy frequencies exceeding 40 cph at the surface, 2) high phytoplankton biomassmore »in spring (up to 125 mg WW L-1), 3) reduced downward transport of heat and 4) persistently depleted epilimnetic CO2 concentrations. CO2 depletion was effectively maintained by intense uptake by phytoplankton (influx up to 30 mmol m-2 d-1) in combination with reduced carbon inputs from the watershed during dry periods. Under eutrophic conditions these effects triggered massive bloom of buoyant cyanobacteria (up to 300 mg WW L-1). Complementary evidence from polynomial regression modelling using long-term datasets revealed that warming is the most important predictor of cyanobacterial abundance during the second half of the last century explaining 78% of the observed positive trend, whereas phosphorus concentration explained only 10% thereof. Together the results from the interannual comparison and the multi-decadal record indicate that hotter and drier climates increase water column stratification and decrease CO2 availability in eutrophic polymictic lakes. This combination catalyzes blooms of buoyant cyanobacteria.« less
  3. Anthropogenic freshwater salinization affects thousands of lakes worldwide, and yet little is known about how salt loading may shift timing of lake stratification and spring mixing in dimictic lakes. Here, we investigate the impact of salinization on mixing in Lakes Mendota and Monona, Wisconsin, by deploying under-ice buoys to record salinity gradients, using an analytical approach to quantify salinity thresholds that prevent spring mixing, and running an ensemble of vertical one-dimensional hydrodynamic lake models (GLM, GOTM, and Simstrat) to investigate the long-term impact of winter salt loading on mixing and stratification. We found that spring salinity gradients between surface andmore »bottom waters persist up to a month after ice-off, and that theory predicts a salinity gradient of 1.3–1.4 g kg-1 would prevent spring mixing. Numerical models project that salt loading delays spring mixing and increases water column stability, with ramifications for oxygenation of bottom waters, biogeochemistry, and lake habitability.« less
  4. Climate change is expected to decrease ice coverage and thickness globally while increasing the variability of ice coverage and thickness on midlatitude lakes. Ice thickness affects physical, biological, and chemical processes as well as safety conditions for scientists and the general public. Measurements of ice thickness that are both temporally frequent and spatially extensive remain a technical challenge. Here new observational methods using repurposed soil moisture sensors that facilitate high spatial–temporal sampling of ice thickness are field tested on Lake Mendota in Wisconsin during the winter 2015/16 season. Spatial variability in ice thickness was high, with differences of 10 cmmore »of ice column thickness over 1.05 km of horizontal distance. When observational data were compared with manual measurements and model output from both the Freshwater Lake (FLake) model and General Lake Model (GLM), ice thickness from sensors matches manual measurements, whereas GLM and FLake results showed a thinner and thicker ice layer, respectively. The FLake-modeled ice column temperature effectively remained at 0°C, not matching observations. We also show that daily ice dynamics follows the expected linear function of ice thickness growth/melt, improving confidence in sensor accuracy under field conditions. We have demonstrated a new method that allows low-cost and high-frequency measurements of ice thickness, which will be needed both to advance winter limnology and to improve on-ice safety.

    « less
  5. China has made a concerted effort to successfully improve water quality of rivers, but lake water quality has not improved. Lakes require controls on both catchment external nutrient loads and in-lake internal loads, where nature-based solutions are coupled with engineered systems to achieve the United Nations Sustainable Development Goals (SDGs).