skip to main content


Title: Uncertainty in projections of future lake thermal dynamics is differentially driven by lake and global climate models
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.  more » « less
Award ID(s):
1933016 1933102
NSF-PAR ID:
10440847
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
PeerJ
Volume:
11
ISSN:
2167-8359
Page Range / eLocation ID:
e15445
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. Abstract. Empirical evidence demonstrates that lakes and reservoirs are warming acrossthe globe. Consequently, there is an increased need to project futurechanges in lake thermal structure and resulting changes in lakebiogeochemistry in order to plan for the likely impacts. Previous studies ofthe impacts of climate change on lakes have often relied on a single modelforced with limited scenario-driven projections of future climate for arelatively small number of lakes. As a result, our understanding of theeffects of climate change on lakes is fragmentary, based on scatteredstudies using different data sources and modelling protocols, and mainlyfocused on individual lakes or lake regions. This has precludedidentification of the main impacts of climate change on lakes at global andregional scales and has likely contributed to the lack of lake water qualityconsiderations in policy-relevant documents, such as the Assessment Reportsof the Intergovernmental Panel on Climate Change (IPCC). Here, we describe asimulation protocol developed by the Lake Sector of the Inter-SectoralImpact Model Intercomparison Project (ISIMIP) for simulating climate changeimpacts on lakes using an ensemble of lake models and climate changescenarios for ISIMIP phases 2 and 3. The protocol prescribes lakesimulations driven by climate forcing from gridded observations anddifferent Earth system models under various representative greenhouse gasconcentration pathways (RCPs), all consistently bias-corrected on a0.5∘ × 0.5∘ global grid. In ISIMIP phase 2, 11 lakemodels were forced with these data to project the thermal structure of 62well-studied lakes where data were available for calibration underhistorical conditions, and using uncalibrated models for 17 500 lakesdefined for all global grid cells containing lakes. In ISIMIP phase 3, thisapproach was expanded to consider more lakes, more models, and moreprocesses. The ISIMIP Lake Sector is the largest international effort toproject future water temperature, thermal structure, and ice phenology oflakes at local and global scales and paves the way for future simulations ofthe impacts of climate change on water quality and biogeochemistry in lakes. 
    more » « less
  3. Abstract Water temperature, ice cover, and lake stratification are important physical properties of lakes and reservoirs that control mixing as well as bio-geo-chemical processes and thus influence the water quality. We used an ensemble of vertical one-dimensional hydrodynamic lake models driven with regional climate projections to calculate water temperature, stratification, and ice cover under the A1B emission scenario for the German drinking water reservoir Lichtenberg. We used an analysis of variance method to estimate the contributions of the considered sources of uncertainty on the ensemble output. For all simulated variables, epistemic uncertainty, which is related to the model structure, is the dominant source throughout the simulation period. Nonetheless, the calculated trends are coherent among the five models and in line with historical observations. The ensemble predicts an increase in surface water temperature of 0.34 K per decade, a lengthening of the summer stratification of 3.2 days per decade, as well as decreased probabilities of the occurrence of ice cover and winter inverse stratification by 2100. These expected changes are likely to influence the water quality of the reservoir. Similar trends are to be expected in other reservoirs and lakes in comparable regions. 
    more » « less
  4. Abstract The physical processes of heat exchange between lakes and the surrounding atmosphere are important in simulating and predicting terrestrial surface energy balance. Latent and sensible heat fluxes are the dominant physical process controlling ice growth and decay on the lake surface, as well as having influence on regional climate. While one-dimensional lake models have been used in simulating environmental changes in ice dynamics and water temperature, understanding the seasonal to daily cycles of lake surface energy balance and its relationship to lake thermal properties, atmospheric conditions, and how those are represented in models is still an open area of research. We evaluated a pair of one-dimensional lake models, Freshwater Lake (FLake) and the General Lake Model (GLM), to compare modeled latent and sensible heat fluxes against observed data collected by an eddy covariance tower during a 1-yr period in 2017, using Lake Mendota in Madison, Wisconsin, as our study site. We hypothesized transitional periods of ice cover as a leading source of model uncertainty, and we instead found that the models failed to simulate accurate values for large positive heat fluxes that occurred from late August into late December. Our results ultimately showed that one-dimensional models are effective in simulating sensible heat fluxes but are considerably less sensitive to latent heat fluxes than the observed relationships of latent heat flux to environmental drivers. These results can be used to focus future improvement of these lake models especially if they are to be used for surface boundary conditions in regional numerical weather models. Significance Statement While lakes consist of a small amount of Earth’s surface, they have a large impact on local climate and weather. A large amount of energy is stored in lakes during the spring and summer, and then removed from lakes before winter. The effect is particularly noticeable in high latitudes, when the seasonal temperature difference is larger. Modeling this lake energy exchange is important for weather models and measuring this energy exchange is challenging. Here we compare modeled and observed energy exchange, and we show there are large amounts of energy exchange happening in the fall, which models struggle to capture well. During periods of partial ice coverage in early winter, lake behavior can change rapidly. 
    more » « less
  5. Abstract

    Ice‐off dates on lakes are some of the longest phenological records in the field of ecology, and some of the best evidence of long‐term climatic change. However, there has been little investigation as to whether the date of ice‐off on a lake impacts spring and summer ecosystem dynamics. Here, I analyzed 274 years of long‐term data from eight north temperate lakes in two climate zones to address whether lakes have ecological memory of ice‐off in the subsequent summer. Five metrics were investigated: epilimnion temperatures, hypolimnion temperatures, hypolimnetic oxygen drawdown, water clarity, and spring primary productivity. The response of the metrics to ice‐off date were variable across latitude and lake type. The northern set of lakes stratified quickly following ice‐off, and early ice‐off years resulted in significantly warmer hypolimnetic temperatures. Oxygen depletion in the hypolimnion was not impacted by ice‐off date, likely because in late ice‐off years the lakes did not fully mix. In the southern lakes, ice‐off date was not correlated to the onset of stratification, with the latter being a more dominant control on hypolimnetic temperature and oxygen. The implications of these findings is that as ice‐off date trends earlier in many parts of the world, the lakes that will likely experience the largest changes in spring and summer ecosystem properties are the lakes that currently have the longest duration of lake ice. In considering a future with warmer winters, these results provide a starting point for predicting how lake ecosystem properties will change with earlier ice‐off.

     
    more » « less