One of the most important physical characteristics driving lifecycle events in lakes is stratification. Already subtle variations in the timing of stratification onset and break-up (phenology) are known to have major ecological effects, mainly by determining the availability of light, nutrients, carbon and oxygen to organisms. Despite its ecological importance, historic and future global changes in stratification phenology are unknown. Here, we used a lake-climate model ensemble and long-term observational data, to investigate changes in lake stratification phenology across the Northern Hemisphere from 1901 to 2099. Under the high-greenhouse-gas-emission scenario, stratification will begin 22.0 ± 7.0 days earlier and end 11.3 ± 4.7 days later by the end of this century. It is very likely that this 33.3 ± 11.7 day prolongation in stratification will accelerate lake deoxygenation with subsequent effects on nutrient mineralization and phosphorus release from lake sediments. Further misalignment of lifecycle events, with possible irreversible changes for lake ecosystems, is also likely.

This paper describes the rationale and the protocol of the first component of the third simulation round of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a, http://www.isimip.org, last access: 2 November 2023) and the associated set of climate-related and direct human forcing data (CRF and DHF, respectively). The observation-based climate-related forcings for the first time include high-resolution observational climate forcings derived by orographic downscaling, monthly to hourly coastal water levels, and wind fields associated with historical tropical cyclones. The DHFs include land use patterns, population densities, information about water and agricultural management, and fishing intensities. The ISIMIP3a impact model simulations driven by these observation-based climate-related and direct human forcings are designed to test to what degree the impact models can explain observed changes in natural and human systems. In a second set of ISIMIP3a experiments the participating impact models are forced by the same DHFs but a counterfactual set of atmospheric forcings and coastal water levels where observed trends have been removed. These experiments are designed to allow for the attribution of observed changes in natural, human, and managed systems to climate change, rising CH₄ and CO₂ concentrations, and sea level rise according to the definition of the Working Group II contribution to the IPCC AR6.

Lake thermal dynamics have been considerably impacted by climate change, with potential adverse effects on aquatic ecosystems. To better understand the potential impacts of future climate change on lake thermal dynamics and related processes, the use of mathematical models is essential. In this study, we provide a comprehensive review of lake water temperature modeling. We begin by discussing the physical concepts that regulate thermal dynamics in lakes, which serve as a primer for the description of process‐based models. We then provide an overview of different sources of observational water temperature data, including in situ monitoring and satellite Earth observations, used in the field of lake water temperature modeling. We classify and review the various lake water temperature models available, and then discuss model performance, including commonly used performance metrics and optimization methods. Finally, we analyze emerging modeling approaches, including forecasting, digital twins, combining process‐based modeling with deep learning, evaluating structural model differences through ensemble modeling, adapted water management, and coupling of climate and lake models. This review is aimed at a diverse group of professionals working in the fields of limnology and hydrology, including ecologists, biologists, physicists, engineers, and remote sensing researchers from the private and public sectors who are interested in understanding lake water temperature modeling and its potential applications.

Empirical evidence demonstrates that lakes and reservoirs are warming across the globe. Consequently, there is an increased need to project future changes in lake thermal structure and resulting changes in lake biogeochemistry in order to plan for the likely impacts. Previous studies of the impacts of climate change on lakes have often relied on a single model forced with limited scenario-driven projections of future climate for a relatively small number of lakes. As a result, our understanding of the effects of climate change on lakes is fragmentary, based on scattered studies using different data sources and modelling protocols, and mainly focused on individual lakes or lake regions. This has precluded identification of the main impacts of climate change on lakes at global and regional scales and has likely contributed to the lack of lake water quality considerations in policy-relevant documents, such as the Assessment Reports of the Intergovernmental Panel on Climate Change (IPCC). Here, we describe a simulation protocol developed by the Lake Sector of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) for simulating climate change impacts on lakes using an ensemble of lake models and climate change scenarios for ISIMIP phases 2 and 3. The protocol prescribes lake simulations driven by climate forcing from gridded observations and different Earth system models under various representative greenhouse gas concentration pathways (RCPs), all consistently bias-corrected on a 0.5^∘ × 0.5^∘ global grid. In ISIMIP phase 2, 11 lake models were forced with these data to project the thermal structure of 62 well-studied lakes where data were available for calibration under historical conditions, and using uncalibrated models for 17 500 lakes defined for all global grid cells containing lakes. In ISIMIP phase 3, this approach was expanded to consider more lakes, more models, and more processes. The ISIMIP Lake Sector is the largest international effort to project future water temperature, thermal structure, and ice phenology of lakes at local and global scales and paves the way for future simulations of the impacts of climate change on water quality and biogeochemistry in lakes.

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.

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 and bottom waters persist up to a month after ice‐off, and that theory predicts a salinity gradient of 1.3–1.4 g kg⁻¹ 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.