Abstract
Lake surface temperatures are projected to increase under climate change, which could trigger shifts in the future distribution of thermally sensitive aquatic species. Of particular concern for lake ecosystems are when temperatures increase outside the range of natural variability, without analogue either today or in the past. However, our knowledge of when such no-analogue conditions will appear remains uncertain. Here, using daily outputs from a large ensemble of SSP3-7.0 Earth system model projections, we show that these conditions will emerge at the surface of many northern lakes under a global warming of 4.0â°C above pre-industrial conditions. No-analogue conditions will occur sooner, under 2.4â°C of warming, at lower latitudes, primarily due to a weaker range of natural variability, which increases the likelihood of the upper natural limit of lake temperature being exceeded. Similar patterns are also projected in subsurface water, with no-analogue conditions occurring first at low latitudes and occurring last, if at all, at higher latitudes. Our study suggests that global warming will induce changes across the water column, particularly at low latitudes, leading to the emergence of unparalleled climates with no modern counterparts, probably affecting their habitability and leading to rearrangements of freshwater habitats this century.
Similar content being viewed by others
Main
Lake surface water temperature can increase considerably in response to greenhouse warming1,2,3, leading to unprecedented shifts in lake thermal habitats and consequently in lake ecosystem functioning4,5,6. The implications of lake warming are particularly severe when temperatures increase above the thermal tolerance limit of key aquatic species7,8,9. One potential severe consequence of long-term future warming in lakes is the emergence of thermal conditions that are wholly outside the range of historical variability and not analogous to conditions observed either today or in the past10,11. In other words, future climate change could lead to the emergence of no-analogue conditions (that is, when lake temperatures become increasingly new) and likewise exceed temperatures that organisms will be able to adapt to. Future lake temperatures with no current analogue (that is, not found today or in the past) and the disappearance of some extant climates could lead to the extinction of some species and the redistribution of others at the global scale12,13. Projecting the emergence of no-analogue lake temperatures is thus fundamental for adaptation, planning and future risk assessments.
A key challenge for projecting the emergence of no-analogue conditions is to partition the uncertainties introduced by natural variability of local climates and, specifically, to assess when temperatures will increase outside the upper natural variability range. With climate projections using an individual realization of an Earth system model, one is faced with the challenge of deconvolving the forced signal from natural variability, with this becoming particularly pronounced at local scales. With large-ensemble simulations, where the only difference between the individual simulations lies in the initial conditions, the forced signal can be identified through the ensemble mean14. Thereby, this method facilitates a clear identification of natural variability as deviations from the ensemble mean. Previous studies have investigated the emergence of no-analogue conditions in both terrestrial and marine ecosystems11,15,16. In lakes, however, the use of ensemble modelling is less common17,18,19âparticularly the use of large ensembles20âleading to a dearth of knowledge regarding the emergence of no-analogue conditions in one of Earthâs most critical ecosystems. In this Article, we aim to fill this knowledge gap by calculating the emergence of no-analogue conditions in lakes using projections from a 100-member ensemble of historical (1850â2014) and future (2015â2100) simulations conducted with a state-of-the-art Earth system model: the Community Earth System Model version 2 Large Ensemble (CESM2-LE)21, which includes a fully coupled lake model (see Methods).
Using the CESM2-LE, we investigated the forced (anthropogenic) and natural variability components of depth-resolved lake temperature changes during the open-water, ice-free (that is, warm) season (see Methods), enabling identification of: (1) the forced climate change signal itself; and (2) a timescale over which such changes will emerge over natural variability. Conceptually, no-analogue conditions represent unprecedented thermal environments that have not been observed within the historical record. Understanding these conditions is important as they signal a shift beyond the adaptive capacities of current aquatic ecosystems, potentially leading to ecological disruptions. The emergence of no-analogue conditions is subsequently communicated relative to global warming levels (calculated from CESM2-LE). In Fig. 1, we show for illustrative purposes the method for calculating the emergence of no-analogue conditions at the surface of one modelled lake, but the same method is then applied to all simulated lakes, encompassing both surface and interior conditions. In this study, we excluded all model grid points at which lakes experienced at least 8âmonths of annual ice cover during the historical (1850â1900) period (Supplementary Fig. 1). The CESM2-LE lake temperature projections analysed in this study have been validated against in situ and satellite data, as well as against other process-based lake models (Supplementary Figs. 2â10 and Supplementary Table 1).
a,b, Shown for one example lake (situated at 42.9°âN and 111.3°âW) is the method used to estimate the emergence of no-analogue conditions. In a, each grey line represents the simulated warm season surface water temperature from the individual ensemble members (nâ=â100). The black line indicates the ensemble mean. The red lines represent two standard deviations (\(\sigma\)) around the mean (\(\mu\)). The blue horizontal lines represent the natural variability range in lake surface temperature, defined as 2\(\sigma\) around the \(\mu\) of the ensemble (1850â1900). The emergence of no-analogue conditions, shown by the vertical dashed line, is estimated as the time when lake temperatures will permanently exceed the natural variability range. b, Global average surface air temperature anomalies (relative to 1850â1900) from the ensemble mean, which the emergence of no-analogue conditions (a) is compared against.
No-analogue conditions in lake surface waters
The emergence of no-analogue conditions at the surface of lakes is projected to occur under different levels of future global warming, with strong regional dependence (Fig. 2a). These conditions indicate a considerable shift in lake thermal conditions that surpasses historical variability, potentially leading to notable ecological and biological impacts. Low-latitude lakes are projected to experience no-analogue conditions first (that is, under lower levels of global warming; Fig. 2 and Supplementary Fig. 11). This suggests that these ecosystems, which are already adapted to narrower thermal ranges, may be more vulnerable to early changes, highlighting the importance of considering regional differences in climate change impacts. To investigate further the influence of lake location on the emergence of no-analogue conditions, we grouped lakes by the thermal regions in which they are found22. Within each thermal region, we calculated the average and standard deviation of the global mean surface air temperatures under which no-analogue conditions occur at the lake surface (Fig. 2c). Our results suggest that no-analogue conditions at the lake surface will emerge under lower levels of global warming within the warmest thermal lake regions (for example, northern warm, tropical hot and southern hot lakes). For the colder lakes situated in, for example, northern frigid or southern temperate regions, the emergence of no-analogue conditions is projected to occur only after global warming reaches ~4.0â°C (Fig. 2c). Our analysis also suggests that a proportion of the studied lakesâprimarily those situated at higher latitudesâare unlikely to experience no-analogue warm season temperature conditions during this century (Fig. 2a), at least under the warming scenario considered.
a, Projected emergence of no-analogue conditions at the lake surface under global warming. This illustrates the magnitude of projected global warming that our simulations suggest will result in no-analogue conditions. b, Distribution of lake thermal regions worldwide. c, Emergence of no-analogue conditions across the thermal regions during the warm season. The data represent meansâ±âs.d. NC, northern cool; NF, northern frigid; NH, northern hot; NT, northern temperate; NW, northern warm; SH, southern hot; ST, southern temperate; SW, southern warm; TH, tropical hot.
Role of anthropogenic forcing and natural variability
The emergence of no-analogue conditions in lakes is determined by both the anthropogenic trend and the range of natural variability in lake water temperature, both of which vary with latitude (Supplementary Fig. 12) and thus across lake thermal regions (Fig. 3). It is the ratio of these quantities that determines when no-analogue conditions emerge. When the ratio of the trend to the natural variability is high, no-analogue conditions will occur sooner as lake temperatures can reach the natural upper bounds earlier. In contrast, if this ratio is low, no-analogue conditions will occur later. In Fig. 3, we compare the trend (2000â2100) and natural variability (1850â1900) in lake surface water temperature. Our projections suggest that during this century lakes situated in colder climates (northern frigid and northern cool) will experience strong surface temperature trends ~1.5 times higher than lakes situated in warmer regions (for example, northern hot and tropical hot). However, these cold lakes also experience the largest natural variability, which is approximately twice that projected in the climatologically warmer lakes. Both the trend and natural variability in lake surface water temperature closely follow those projected for surface air temperature (Supplementary Fig. 13). This is expected, given the degree to which surface air temperature variations influence lake surface temperature. However, the amplified response of lake surface temperature to air temperature changes is also influenced by lake depth, with deeper lakes typically shown to experience greater long-term change23,24. Other atmospheric drivers can also influence lake surface temperatures, including wind speed25 and solar radiation26, among others27, all of which are included in our model projections (see Methods). There is also a lag between air and water temperature, particularly for large lakes28, which can result in a partial decoupling of air and water temperature at seasonal timescales. These factors, as well as the influence of lake morphometry29 and the presence of winter ice cover1, account for some of the variability in airâwater temperature relationships and partly explain why natural variability and warming rates can differ somewhat between air and surface water temperature (Supplementary Fig. 13). Overall, our results suggest that although warmer tropical lakes experience lower temperature trends, their level of background natural variability is proportionally even lower, so they can experience relatively early emergence of no-analogue conditions.
a,b, Rate of change in lake surface water temperature across the studied lakes from 2000â2100 (a) and differences in these trends across lake thermal regions (b). c,d, Natural variability (1850â1900) in surface water temperatures worldwide (c) and differences in this variability across thermal regions (d). The trends in surface water temperature were calculated from an average across the 100 ensemble members using a linear regression model. The data in b and d represent meansâ±âs.d.
Novel conditions in subsurface waters
Understanding the emergence of no-analogue conditions at the lake surface will be critical for assessing the impacts of climate change on lake ecosystems, particularly for organisms and aquatic species that cannot shift their vertical distributions to track a suitable thermal habitat. However, some aquatic species have adapted to light and thermal conditions at depth and hence will not be directly influenced by surface water temperature changes. It is therefore important to also consider the emergence of no-analogue conditions in subsurface layers, and to identify when the natural range of variations in temperature will be exceeded throughout the water column. We investigated differences in the emergence of no-analogue conditions by depth in our global simulations. For consistency, we initially focused on the subsurface layers ranging from 0â8.5âm depth. Lakes that are 8.5âm deep (that is, in terms of average depth) or shallower represent 96% of the global lake distribution30. In some lake thermal regions, our simulations suggest clear differences across depths in the emergence of no-analogue conditions (Fig. 4). The most apparent variations with depth were found in northern lakes (northern frigid to northern hot) where the emergence of no-analogue conditions is projected to occur much later for deeper waters. For example, in northern temperate lakes, no-analogue conditions at the lake surface will emerge under a projected global warming of 2.7â±â0.6â°C. However, at a depth of ~8.5âm, no-analogue conditions become emergent only for a global temperature increase of 4.5â°C (corresponding to ~2100âce) (Fig. 4). Several lakes in this thermal regimen are unlikely to experience no-analogue conditions in deeper waters before 2100âce under the projected levels of global warming considered in this study. However, in tropical hot and many southern lake thermal regions (southern hot, southern warm and southern temperate), there are only marginal differences in the projected emergence of no-analogue conditions with varying depth (Fig. 4fâi). In tropical hot lakes, no-analogue conditions will emerge under a global warming of 2.5â±â0.5â°C at the surface or 2.6â±â0.7â°C at a depth of 8.5âm.
aâi, Vertical profiles of the calculated emergence of no-analogue conditions across thermal lake regions. Subsurface patterns are shown for lakes in northern frigid (a), northern cool, (b) northern temperate (c), northern warm (d), northern hot (e), tropical hot (f), southern hot (g), southern warm (h) and southern temperate regions (i). For each thermal region, we selected lakes that were at least 8.5âm deep and calculated the average (solid line) and s.d. (dashed lines) in the emergence of no-analogue conditions during the warm season.
Our simulations also suggest clear differences in the emergence of no-analogue conditions in bottom water temperature (Fig. 5a). Notably, warmer low-latitude lakes experience no-analogue conditions first and cooler lakes situated at higher latitudes (that is, cooler lake thermal regions) experience no-analogue conditions at depth last, if at all. An ordinary least-squares analysis (sample sizeâ=â1,084) suggested that, in lakes with a depth of 8.5âm (and in which no-analogue conditions emerge), 87.7% of the variance in the emergence of no-analogue conditions can be explained by the ratio of the trend to the natural variability. Our simulations suggested that the calculated trends in bottom water temperature during the twenty-first century will be highest at low latitudes and the natural variability will be lowest at low latitudes (Fig. 5 and Supplementary Fig. 14). The trend in deep-water temperature is primarily determined by the trend in surface water temperature and, critically, the warming signal that is transferred vertically within a water body. The latter is influenced by the frequency of overturning days during the ice-free season. Specifically, the trend and natural variability in bottom water temperatures follow closely those projected at the surface in tropical hot lakes and shallow polymictic lakes worldwide. In our simulations, 65% of lakes (with a depth of 8.5âm; 590 lakes) in the tropics overturn for more than 182âd annually. However, the trend of bottom water in lakes that stratify has a weak relationship with the trends in surface water temperature. In northern temperate and northern cool lakes with a mean depth of 8.5âm and where no-analogue conditions of the bottom water are projected to emerge this century, the frequency of overturning days explains 73.1% of the variance (sample sizeâ=â294) in the trend of bottom water temperature, which outweighs the impacts of the trend in surface water. These lakes are dimictic and experience ice coverage annually, and the amount of downward heat penetration increases along with the frequency of overturning across lakes. In lakes that experience intermittent mixing or stratify permanently, bottom waters will, to a large extent, be shielded from much of the influence of surface temperature, although deep mixing31,32 and the downward diffusion of heat33,34 are able to introduce heat into deeper water. The slow increase in bottom water temperature in these lakes, however, would lead to a delayed emergence of no-analogue conditions at depth.
a, Projected emergence of no-analogue conditions by bottom water temperature increase under global warming. Lakes where no-analogue conditions are unlikely to occur this century, under the global warming levels considered, are shown in black. b, Rate of change (2000â2100) in bottom water temperature during the warm season, calculated from the average of the 100-member ensemble. c, Natural variability in bottom water temperature during the warm season, estimated as twice the standard deviation (\(\sigma\)) around the mean (\(\mu\)) of the large ensemble between 1850 and 1900. d, Percentage of overturning (that is, mixed) days during the warm season (1850â1900), calculated as the average of the 100-member ensemble.
Impact of no-analogue conditions on lake ecosystems
The ecological responses expected from the emergence of no-analogue conditions at the surface or at depth are likely to be idiosyncratic and will depend on attributes such as species adaptive capacity, the ability of aquatic species to migrate and the disruption of ecological interactions. Nevertheless, the emergence of no-analogue lake temperatures over the twenty-first century implies that there will be a substantial disruption of the normal lake conditions that have shaped lake ecosystems in the past. Lake ecosystems are already being threatened by several anthropogenic stressors35, including over-fishing36, eutrophication37,38 and water abstraction39,40. Climate change exacerbates the degradation of lake ecosystems by pushing organisms to adapt to warmer waters. In addition, extreme events such as lake heatwaves are expected to increase under persistent anthropogenic perturbations41, with potentially devastating effects on lake ecosystems. Indeed, the occurrence of extreme events in recent years has highlighted the serious consequences of reaching historically unprecedented conditions7,9,42,43,44.
Although our study focuses on the emergence of no-analogue conditions relative to natural variability, it is crucial to consider the absolute tolerances of organisms inhabiting lakes. Even temporary deviations in temperature beyond these absolute tolerances can have ecological impacts, including stress, mortality and disruption of ecosystem services7,8. For instance, species adapted to specific thermal niches may face critical thresholds beyond which survival and reproductive success are compromised, regardless of whether such conditions are unprecedented within the historical context of variability. Cold water fish species such as lake trout (Salvelinus namaycush), which have optimal environmental temperatures of less than ~12â°C, could be particularly susceptible to no-analogue conditions, or even to changes that are not considered to be outside of the bounds of natural variability. In lakes experiencing no-analogue warming, temperatures could exceed these limits for a sustained duration, potentially leading to decreased growth rates, increased susceptibility to disease and higher mortality rates. Similarly, for aquatic plants such as the water lily (Nymphaea odorata), which thrive in temperatures up to ~30â°C, mean warming of lakes might push temperatures beyond this threshold, impacting their photosynthesis and growth. Several uncertainties must also be considered when evaluating these impacts, including that speciesâ thermal tolerances vary45,46, genetic and phenotypic variability can influence adaptability47,48 and ecological interactions may amplify or mitigate effects49,50. For example, changes in predatorâprey dynamics, competition for resources and symbiotic relationships can influence how species respond to altered thermal conditions. Shifts in community composition and trophic interactions may have cascading effects throughout the ecosystem, altering nutrient cycling and productivity. Additionally, the effects of no-analogue conditions may be influenced by habitat complexity, further complicating predictions of ecological responses. Thus, although our study provides a framework for identifying no-analogue conditions, the actual ecological impacts will depend on a multitude of factors, requiring comprehensive assessment for accurate predictions and mitigation strategies.
The simulations employed in this study demonstrate the value of scrutinizing climate change perturbations on ecosystem thermal environments across the water column with respect to surface conditions, which has been the focus of previous studies exploring future warming impacts on lakes. Our results thus shed light on a complex picture of the impacts of climate change on lake ecosystems by adding the vertical dimension and recognizing the importance of natural variability. We find that ongoing climate change will induce changes across the water column, particularly in low-latitude lakes, leading to the emergence of unparalleled climates with no modern counterparts. If anthropogenic emissions continue to rise, we project that the natural variability in water temperature will be exceeded and no-analogue conditions will become ubiquitous in lakes this century, probably affecting their habitability. Assuming that aquatic organisms are adapted to the current thermal environment, no-analogue conditions may lead to substantial rearrangements of freshwater habitats this century.
Methods
Lake temperature projections
We evaluated the influence of climate change on the thermal environment of lakes by analysing simulations from CESM2-LE21. Most notably, we investigated lake temperature projections generated by the Lake, Ice, Snow, and Sediment Simulator (LISSS)51, which is a one-dimensional (1D) thermodynamic lake model embedded within and coupled to the land surface module of CESM (that is, the Community Land Model; version 5)52. The CESM2-LE lake simulations have been used, and described in detail, in previous studies20,53. LISSS is a thermal diffusion modelâit is an adapted version of the Hostetler 1D modelâand was developed specifically for simulating lakes in Earth system models. Compared with other lake models, LISSS has lower computational expense and requires minimal lake-specific calibration. More complex models, such as the General Lake Model54 or SimStrat55, improve the characterization of deep lakes, but also often require lake-specific data or calibration and have not yet been integrated into Earth system models. LISSS was developed to accurately simulate lakes in Earth system models, but also to be computationally efficient. In brief, LISSS simulates the thermal environment of lakes by balancing the surface energy budget and solving the 1D thermal diffusion equation for lake water temperature, then it simulates vertical mixing, including diffusion of wind-driven eddies, convective mixing, molecular diffusion and enhanced mixing due to unresolved 3D processes. LISSS has been tested extensively in previous studies investigating lake responses to climate change20,51,56,57. LISSS has been shown to provide an accurate representation of water temperature, surface fluxes and the thickness of ice and snow in previous studies for which observational data were available20,51,58,59,60.
LISSS was used to simulate depth-resolved lake water temperatures at a longitudeâlatitude resolution of ~1° by ~1° and a 30âmin time step. Lakes simulated within each grid in CESM2-LE are based on the depth and surface area of all known lakes in that region61,62. The projections thus represent a typical lake for each grid cell57,63,64. Each lake grid has ten vertical layers in the simulation. A piecewise cubic Hermite interpolation method65 was applied to interpolate lake water temperatures onto desired depths. The CESM2-LE simulations include historical projections from 1850â2014 and future (2015â2100) projections simulated under the Shared Socioeconomic Pathways forcing scenario SSP3-7.0 (ref. 66). An ensemble approach was used in CESM2-LE including 100 members, with each member initialized from a set of different combinations of initial states of the atmosphere and ocean21.
Although previous studies have validated lake temperature simulations from LISSS20,51,56,57, here we performed an additional validation of the outputs from CESM2-LE (Supplementary Figs. 2â10 and Supplementary Table 1). Most notably, we compared the lake temperature simulations with satellite-derived surface water temperatures from ARC-Lake version 3.0 (ref. 67) and depth-resolved in situ temperatures from ref. 68. All comparisons were made between the observations and simulations by extracting the water temperature time series from the grid in which each lake was located. In the case of satellite-derived lake surface temperatures, we highlight that the obtained value is sensitive to the skin temperature of the water, which is the temperature of a layer <0.1âmm thick from which thermal radiation is emitted by the lake. Thus, the satellite datum is an estimate of this skin temperature, which may differ from the temperature as measured by a thermometer, or as simulated by a lake model, a few centimetres below the airâwater interface. Typically, the difference between skin and sub-skin lake surface temperature is of the order of a few tenths of a degree69,70,71. This temperature difference is commonly known as the cool skin effect72,73, the magnitude of which is influenced by airâwater surface heat fluxes70,71. Satellite lake surface temperatures have been used to quantify worldwide aspects of lake thermal dynamics22,74,75. Regarding our depth-resolved lake temperature observations, these are mostly only available for lakes in the Northern Hemisphere, but also include a few tropical systems68. In turn, our depth-resolved lake temperature simulations are not as well validated at low latitude but do provide a reasonable validation of our modelled data when available. To verify that our simulations were able to simulate multidecadal variations in lake temperature (and would therefore also be informative with respect to future climate change impacts), we compared the CESM2-LE simulations with long-term (1910â2018) in situ water temperature observations from Wörthersee, Austria (Supplementary Fig. 8), available from ref. 76.
We also compared the CESM2-LE water temperature simulations with global-scale projections from the Intersectoral Impact Model Intercomparison Project phase 2b (ISIMIP2b) Lake Sector18. Specifically, we repeated the model validation exercise described above for CESM2-LE, but for the ISIMIP2b projections, then compared the validation summary statistics. Here we compared the CESM2-LE simulations with those simulated by more complex models, including: (1) SimStrat-UoG, a modified version of the buoyancy-extended kâε model SimStrat, which is a kâε turbulence model with buoyancy and internal seiche parameterization77,78; (2) LAKE, an extended one-dimensional model of thermodynamic, hydrodynamic and biogeochemical processes in lakes79,80 that has been tested extensively with respect to thermal and ice regimens under contrasting climate conditions81,82,83,84; and (3) VIC-LAKE, a 1D model derived from the Variable Infiltration Capacity Macroscale Hydrologic Model85. These models were included in ISIMIP2b to simulate lakes worldwide at a 0.5° by 0.5° grid resolution using climate projections from GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR and MIROC5 as input. The ISIMIP2b Lake Sector simulations have been used extensively to investigate the impact of climate change on lakes57,63,64. The model validation exercise demonstrated that despite the spatial resolution of CESM2-LE being coarser (1° by 1°) than that of ISIMIP2b (0.5° by 0.5°), the simulations provide comparable performance for the lakes tested. We also note that in the lakes tested the standard deviations of the observations and simulations from CESM2-LE are comparable at the lake surface, but the difference typically (although not always) increases in deeper water.
The CESM2-LE simulations were also used in this study to estimate the number of overturning (that is, mixed) days in lakes worldwide. Although there is no universal definition of stratification/mixing in lakes86,87, here we define an overturning day as the day that has at least one full water column overturning during a 24-h period.
Emergence of no-analogue conditions
The emergence of no-analogue conditions was estimated from the 100-member CESM2-LE as the time when an increase in lake temperature is permanently above its natural variability. In this study, we considered the natural variability in lake temperature to be twice the standard deviation (\(\sigma\)) around the mean (\(\mu\)) of the ensemble between 1850 and 1900. This period was chosen to represent the pre-industrial period when the influence of anthropogenic emissions was considered minimal. Thus, we took lake temperatures of all of the ensembles (nâ=â100) and all of the years between 1850 and 1900 as a sample set that was formed under a pre-industrial climate. The natural variability range was defined as twice the standard deviation around the mean of the pre-industrial sample set. Along with continuous warming, no-analogue conditions emerge when \(\mu\)âââ2\(\sigma\) (the lower red solid line in Fig. 1a) of the 100-member ensemble permanently exceeds the natural variability range (that is, when the upper threshold is exceeded). The calculated emergence of no-analogue conditions was subsequently communicated relative to global warming levels, which provided a way to reference change in climatic conditions independent of the emissions scenario used in the simulations.
In this study, we only considered the open-water, ice-free (that is, warm) season when estimating the emergence of no-analogue conditions (Supplementary Fig. 1). We define the warm season as January to March in the Southern Hemisphere (90°âSâ23.5°âS) and July to September in the Northern Hemisphere (23.5°âNâ90°âN). In tropical lakes (23.5°âSâ23.5°âN), the warm season refers to the entire year. However, we note that our definition of a warm season does not strictly include the warmest months in all of the studied lakes, with some lakes experiencing, for example, higher surface water temperatures in June compared with September. In this study we also excluded all model grid points where lakes experienced at least 8âmonths of annual ice cover during the 1850â1900 period (Supplementary Fig. 1). Finally, the lake temperature for the warm season was calculated for each lake and used to estimate the natural variability and emergence of no-analogue conditions.
To investigate across-lake differences in the emergence of no-analogue conditions, we grouped the studied lakes according to the lake thermal region in which they are found22. This classification allowed us to better understand the varying capacities of different lake types to cope with thermal changes, and the potential onset of no-analogue conditions serves as an indicator of when and where these ecosystems may face unprecedented thermal stress.
Limitations of the model projections
Although we consider our results to be robust and believe that they bridge an important knowledge gap in climate change assessments, there are some limitations to consider when interpreting our key findings. First, our simulations are for representative lakes for each CESM2-LE grid. Individual lakes within a grid will behave in a different way from the typical lake considered as, for example, lake surface area and depth are known to influence both lake surface and bottom water temperatures due to differences in their thermal inertia23,29 and thus could also influence the computed emergence of no-analogue conditions. We also note that as our projections were generated with a 1D process-based lake model, which largely represents lake mean conditions88,89, horizontal features in lakes (for example, lateral advection) and the intra-lake responses to climate change will not have been captured23,24.
In this study, we also calculated the emergence of no-analogue conditions based on the results of SimStrat-UoG model simulation from the local lake sector of ISIMIP. Instead of simulating virtual lakes based on average depth in individual grids, this project configured the lake model with actual bathymetry for 59 lakes. Of these, 56 had corresponding lake grids in CESM2-LE. The calculation of emergence of no-analogue conditions for the ISIMIP local lake simulation was based on the single realization as the meteorological forcing for the lake simulation was achieved from the output of the GFDL-ESM2M model under the RCP8.5 scenario in CMIP5. In the calculation, the natural variability range was calculated as twice the standard deviation around the mean of water temperature during 1850â1900. The forcing signal was identified as the fourth-order polynomial fit of the raw data and the natural variability (that is, the envelope denoted by the two red solid lines in Fig. 1a) was calculated as the two standard deviations of temperature anomalies (by retracting the fourth-order polynomial fit values from the raw data) in a sliding window of 30âyears. After calculating the emergence of no-analogue conditions, we compared the year of no-analogue conditions emerging against the magnitude of projected global warming. To maintain consistency in the calculation method, we applied the same method to individual members of the CESM2-LE. The results indicated that the emergence of no-analogue conditions in the surface water estimated by the ISIMIP local project is close to the CESM2-LE-based estimations, except for four lakes (that is, Annie, Biel, Kinneret and Sau) in the northern mid- and high latitudes that have a difference of about 1â°C in warming magnitude. Notably, despite the difference in the model configuration of depth, the emergence of no-analogue conditions projected by the CESM2-LE shares the same vertical pattern with the emergence estimated by the ISIMIP local project, which uses the actual depths of lakes. For example, in Lake Kivu, which is in the tropics and has a maximum depth of more than 490âm, the CESM2-LE projects the same pattern of emergence as the ISIMIP local project for the upper 50âm layer. Indeed, there exists a distinct difference in one of the 56 lakes (that is, Lake Rotorua, New Zealand), in which the vertical patterns of emergence are opposite to the two simulations. Although the model configuration of virtual lake depth induces bias in the projection of emergence of no-analogue conditions, the CESM2-LE-based projections are in the same range of uncertainty as the estimations based on the simulation of actual lakes.
Another limitation of the lake temperature simulations is that the light attenuation coefficient (Kd) does not vary temporally but is assumed to be constant for a given lake (but does vary from one lake to another). Although this is common in 1D global lake simulations18,57,59, it means that changes in water transparency during the twentieth and twenty-first century are not considered in the historic and future projections of lake temperature. Transparency can either increase or decrease in the future, as it has during the historic period, due to changes inâamong other thingsâland management practices and precipitation patterns90. These changes in transparency can either amplify or suppress lake surface temperatures under climate change91. Although it is unclear how transparency will change during the twenty-first century, we now discuss how some of the potential changes could influence our long-term projections of water temperature. Regarding surface water temperature, stronger light attenuation can lead to higher temperatures in spring and summer (due to higher absorption of solar radiation) but result in faster cooling in autumn due to the presence of a shallower thermocline92,93.
As suggested by our results, the influence of climatic change on bottom water temperature will differ across lakes depending on their mixing regimen, with, for example, bottom temperatures experiencing a stronger response to climatic warming in lakes that frequently mix. An increase in transparency in these lakes will probably result in an even greater increase in bottom temperature, due to an increase in the depth at which solar radiation can penetrate through the water column. Conversely, a decrease in transparency could result in a weaker relationship between air and bottom temperature, or even a cooling of the deeper layers, due to well-mixed lakes potentially shifting to a more stable mixing regimen94,95. However, such transparency-induced alterations in mixing regimens are most likely to occur in marginal lakes (that is, those that historically transition between two mixing regimen types)75,95. In seasonally stratifying lakes, a change in transparency could influence how solar radiation is absorbed in the water column, essentially altering the vertical partitioning of heat. A scenario of increasing water transparency could result in further amplification of bottom warming compared with those projected in this study91, leading to an earlier emergence of no-analogue conditions. Conversely, a decrease in transparency could result in cooler bottom water temperatures in summer due to thermal shielding42,91,96,97, thus leading to a later emergence of no-analogue conditions. The modelling protocol followed in this study assumed no temporal change in water transparency nor changes in topography (for example, wind sheltering), both of which can influence lake temperature98, but instead focused solely on the more robust future projections of changes in climate.
It is important to consider that our study was based on simulations from a single lake model coupled to an Earth system model. Recent studies have highlighted the benefit of investigating climate change impacts in lakes with the use of lake model ensembles17,18. Indeed, multi-lake model ensemble simulations are increasingly used to obtain robust assessments of freshwater ecosystem responses to climate change. A lake model ensemble approach was not feasible within CESM2-LE, given its large computational requirements. Future studies could benefit from considering multiple lake model simulations when investigating lake responses to climate change within a large ensemble. Despite these limitations, our results provide an important step forward in understanding changes in lake thermal conditions within a warming world.
Data availability
The model outputs of CESM2-LE are available at https://www.cesm.ucar.edu/projects/community-projects/LENS2/data-sets.html. The main results of this work are openly accessible at https://zenodo.org/records/11427870 (ref. 99). ISIMIP2b local lake simulations are available at https://data.isimip.org/10.48364/ISIMIP.563533. ISIMIP2b global lake simulations are available at https://data.isimip.org/10.48364/ISIMIP.931371. The in situ vertical temperature profile data used in this study are available from https://doi.org/10.6073/pasta/73eb330b4be65faadf38eb2288665848. Satellite-derived lake surface water temperature data from ARC-Lake are available from http://www.laketemp.net/home_ARCLake/targets_phase3.php.
Code availability
The analysis and visualization codes are available from https://github.com/geohuanglei/TOE_LAKE_Temperature and https://zenodo.org/records/11427870 (ref. 99).
Change history
24 July 2024
In the version of the article initially published, in ref. 89, the surname of the first author was incorrect and has now been amended to Ramón in the HTML and PDF versions of the article.
References
Austin, J. A. & Colman, S. M. Lake Superior summer water temperatures are increasing more rapidly than regional air temperatures: a positive iceâalbedo feedback. Geophys. Res. Lett. 34, L06604 (2007).
Schneider, P. & Hook, S. J. Space observations of inland water bodies show rapid surface warming since 1985. Geophys. Res. Lett. 37, L22405 (2010).
OâReilly, C. M. et al. Rapid and highly variable warming of lake surface waters around the globe. Geophys. Res. Lett. 42, 10773â10781 (2015).
Cline, T. J., Bennington, V. & Kitchell, J. F. Climate change expands the spatial extent and duration of preferred thermal habitat for Lake Superior fishes. PLoS ONE 8, e62279 (2013).
Jonsson, T. & Setzer, M. A freshwater predator hit twice by the effects of warming across trophic levels. Nat. Commun. 6, 5992 (2015).
Kraemer, B. et al. Climate change drives widespread shifts in lake thermal habitat. Nat. Clim. Change 11, 521â529 (2021).
Till, A. et al. Fish die-offs are concurrent with thermal extremes in north temperate lakes. Nat. Clim. Change 9, 637â641 (2019).
Kelly, S. et al. Warming winters threaten peripheral Arctic charr populations of Europe. Climatic Change 163, 599â618 (2020).
Tye, S. P. et al. Climate warming amplifies the frequency of fish mass mortality events across north temperate lakes. Limnol. Oceanogr. Lett. 7, 510â519 (2022).
Williams, J. W. & Jackson, S. T. Novel climates, no-analog communities, and ecological surprises. Front. Ecol. Environ. 5, 475â482 (2007).
Mora, C. et al. The projected timing of climate departure from recent variability. Nature 502, 183â187 (2013).
Comte, L. & Olden, J. D. Climatic vulnerability of the worldâs freshwater and marine fishes. Nat. Clim. Change 7, 718â722 (2017).
Woolway, R. I. & Maberly, S. C. Climate velocity in inland standing waters. Nat. Clim. Change 10, 1124â1129 (2020).
Deser, C. et al. Insights from Earth system model initial-condition large ensembles and future prospects. Nat. Clim. Change 10, 277â286 (2020).
Schlunegger, S. et al. Emergence of anthropogenic signals in the ocean carbon cycle. Nat. Clim. Change 9, 719â725 (2019).
Silvy, Y. et al. Human-induced changes to the global ocean water masses and their time of emergence. Nat. Clim. Change 10, 1030â1036 (2020).
Moore, T. N. et al. LakeEnsemblR: an R package that facilitates ensemble modelling of lakes. Environ. Model. Softw. 143, 105101 (2021).
Golub, M. et al. A framework for ensemble modelling of climate change impacts on lakes worldwide: the ISIMIP Lake Sector. Geosci. Model. Dev. 15, 4597â4623 (2022).
La Fuente, S. et al. Ensemble modeling of global lake evaporation under climate change. J. Hydrol. 631, 130647 (2024).
Huang, L. et al. Emerging unprecedented lake ice loss in climate change projections. Nat. Commun. 13, 5798 (2022).
Rodgers, K. B. et al. Ubiquity of human-induced changes in climate variability. Earth Syst. Dyn. 12, 1393â1411 (2021).
Maberly, S. C. et al. Global lake thermal regions shift under climate change. Nat. Commun. 11, 1232 (2020).
Woolway, R. I. & Merchant, C. J. Amplified surface temperature response of cold, deep lakes to interannual air temperature variability. Sci. Rep. 7, 4130 (2017).
Calamita, E. et al. On the role of local depth and latitude on surface warming heterogeneity in the Laurentian Great Lakes. Inland Waters 11, 208â222 (2021).
Woolway, R. I. et al. Northern Hemisphere atmospheric stilling accelerates lake thermal responses to a warming world. Geophys. Res. Lett. 46, 11983â11992 (2019).
Schmid, M. et al. Excess warming of a Central European lake driven by solar brightening. Water Resour. Res. 52, 8103â8116 (2016).
Edinger, J. E., Duttweiler, D. W. & Geyer, J. C. The response of water temperatures to meteorological conditions. Water Resour. Res. 4, 1137â1143 (1968).
Toffolon, M. et al. On the use of averaged indicators to assess lakesâ thermal response to changes in climatic conditions. Environ. Res. Lett. 15, 034060 (2020).
Kraemer, B. et al. Morphometry and average temperature affect lake stratification responses to climate change. Geophys. Res. Lett. 42, 4981â4988 (2015).
Messager, M. et al. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nat. Commun. 7, 13603 (2016).
Livingstone, D. M. An example of the simultaneous occurrence of climate-driven âsawtoothâ deep-water warming/cooling episodes in several Swiss lakes. Verh. Int. Ver. Limnol. 26, 822â826 (1997).
Peters, F. et al. Modeling 50 years of historical temperature profiles in a large central European lake. Limnol. Oceanogr. 47, 186â197 (2002).
Ambrosetti, W. & Barbanti, L. Deep water warming in lakes: an indicator of climatic change. J. Limnol. 58, 1â9 (1999).
Verburg, P. & Hecky, R. E. The physics of the warming of Lake Tanganyika by climate change. Limnol. Oceanogr. 54, 2418â2430 (2009).
Jenny, J.-P. et al. Scientistâs warning to humanity: rapid degradation of the worldâs large lakes. J. Great Lakes Res. 46, 686â702 (2020).
Allan, J. D. et al. Overfishing of inland waters. Bioscience 55, 1041â1051 (2005).
Moss, B. et al. Allied attack: climate change and eutrophication. Inland Waters 1, 101â105 (2011).
Foley, B., Jones, I. D., Maberly, S. C. & Rippey, B. Long-term changes in oxygen depletion in a small temperate lake: effects of climate change and eutrophication. Freshw. Biol. 57, 278â289 (2012).
Wurtsbaugh, W. A. et al. Decline of worldâs saline lakes. Nat. Geosci. 10, 816â821 (2017).
Schulz, S. Climate change or irrigated agricultureâwhat drives the water level decline of Lake Urmia. Sci. Rep. 10, 236 (2020).
Woolway, R. I. et al. Lake heatwaves under climate change. Nature 589, 402â407 (2021).
Jankowski, T., Livingstone, D. M., Bührer, H., Forster, R. & Niederhauser, P. Consequences of the 2003 European heat wave for lake temperature profiles, thermal stability, and hypolimnetic oxygen depletion: implications for a warmer world. Limnol. Oceanogr. 51, 815â819 (2006).
De Eyto, E. et al. Response of a humic lake ecosystem to an extreme precipitation event: physical, chemical, and biological implications. Inland Waters 6, 483â498 (2016).
Bartosiewicz, M., Laurion, I., Clayer, F. & Maranger, R. Heat-wave effects on oxygen, nutrients, and phytoplankton can alter global warming potential of gases emitted from a small shallow lake. Environ. Sci. Technol. 50, 6267â6275 (2016).
Angilletta, M. J. Jr., Steury, T. D. & Sears, M. W. Temperature, growth rate, and body size in ectotherms: fitting pieces of a life-history puzzle. Integr. Comp. Biol. 44, 498â509 (2004).
Pörtner, H. O. & Farrell, A. P. Physiology and climate change. Science 322, 690â692 (2008).
Hoffmann, A. A. & Sgrò, C. M. Climate change and evolutionary adaptation. Nature 470, 479â485 (2011).
Reusch, T. B. H. et al. Ecosystem recovery after climatic extremes enhanced by genotypic diversity. Proc. Natl Acad. Sci. USA 102, 2826â2831 (2005).
Walther, G. R. et al. Ecological responses to recent climate change. Nature 416, 389â395 (2002).
Parmesan, C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 37, 637â669 (2006).
Subin, Z. M., Riley, W. J. & Mironov, D. An improved lake model for climate simulations: model structure, evaluation, and sensitivity analyses in CESM1. J. Adv. Model. Earth Syst. 4, M02001 (2012).
Lawrence, D. M. et al. The Community Land Model Version 5: description of new features, benchmarking, and impact of forcing uncertainty. J. Adv. Model. Earth Syst. 11, 4245â4287 (2019).
Woolway, R. I. et al. Lake ice will be less for recreation and transportation under future warming. Earths Future 10, e2022EF002907 (2022).
Bruce, L. C. et al. A multi-lake comparative analysis of the General Lake Model (GLM): stress-testing across a global observatory network. Environ. Model. Softw. 102, 274â291 (2018).
Perroud, M. et al. Simulation of multiannual thermal profiles in deep Lake Geneva: a comparison of one-dimensional lake models. Limnol. Oceanogr. 54, 1574â1594 (2009).
Butcher, J. B. et al. Sensitivity of lake thermal and mixing dynamics to climate change. Climatic Change 129, 295â305 (2015).
Grant, L. et al. Attribution of global lake systems change to anthropogenic forcing. Nat. Geo. 14, 849â854 (2021).
Hu, C. et al. Trends in evaporation of a large subtropical lake. Theor. Appl. Climatol. 129, 159â170 (2017).
Wang, W. et al. Global lake evaporation accelerated by changes in surface energy allocation in a warmer climate. Nat. Geosci. 11, 410â414 (2018).
Xiao, K. et al. Evaporation from a temperate closed-basin lake and its impact on present, past, and future water level. J. Hydrol. 561, 59â75 (2018).
Lehner, B. & Döll, P. Development and validation of a global database of lakes, reservoirs and wetlands. J. Hydrol. 296, 1â22 (2004).
Kourzeneva, E., Asensio, H., Martin, E. & Faroux, S. Global gridded dataset of lake coverage and lake depth for use in numerical weather prediction and climate modelling. Tellus A Dyn. Meteorol. Oceanogr. 64, 15640 (2012).
Vanderkelen, I. et al. Global heat uptake by inland waters. Geophys. Res. Lett. 47, e2020GL087867 (2020).
Woolway, R. I. et al. Phenological shifts in lake stratification under climate change. Nat. Comm. 12, 2318 (2021).
Fritsch, F. N. & Butland, J. A method for constructing local monotone piecewise cubic interpolants. SIAM J. Sci. Comput. 5, 300â304 (1984).
OâNeill, B. C. et al. The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev. 9, 3461â3482 (2016).
Merchant, C. & MacCallum, S. Lake Surface Water Temperature ARC-Lake v3 (1995â2012). University of Reading dataset https://doi.org/10.17864/1947.186 (2018).
Pilla, R. M. et al. Deeper waters are changing less consistently than surface waters in a global analysis of 102 lakes. Sci. Rep. 10, 20514 (2020).
Hook, S. J. et al. Retrieval of lake bulk and skin temperatures using Along-Track Scanning Radiometer (ATSR-2) data: a case study using Lake Tahoe, California. J. Atmos. Ocean. Technol. 20, 534â548 (2003).
Wilson, R. C., Hook, S. J., Schneider, P. & Schladow, G. S. Skin and bulk temperature difference at lake Tahoe: a case study on lake skin effect. J. Geophys. Res. Atmos. 118, 346 (2013).
Hondzo, M. et al. Measurement and scaling of lake surface skin temperatures. Geophys. Res. Lett. 49, e2021GL093226 (2022).
Fairall, C. W. et al. Cool-skin and warm-layer effects on sea surface temperature. J. Geophys. Res. 101, 1295â1308 (1996).
Minnett, P. J., Smith, M. & Ward, B. Measurements of the oceanic thermal skin effect. Deep Sea Res. Pt II Top. Stud. Oceanogr. 58, 861â868 (2011).
Fichot, C. et al. Assessing change in the overturning behaviour of the Laurentian Great Lakes using remotely sensed lake surface water temperatures. Remote Sens. Environ. 235, 111427 (2019).
Woolway, R. I. & Merchant, C. J. Worldwide alteration of lake mixing regimes in response to climate change. Nat. Geosci. 12, 271â276 (2019).
Woolway, R. I. et al. Substantial increase in minimum lake surface temperatures under climate change. Climatic Change 155, 81â94 (2019).
Goudsmit, G. H., Burchard, H., Peeters, F. & Wüest, A. Application of k-â turbulence models to enclosed basins: the role of internal seiches. J. Geophys. Res. 107, 3230 (2002).
Gaudard, A., RÃ¥man VinnÃ¥, L., Bärenbold, F., Schmid, M. & Bouffard, D. Toward an open access to high-frequency lake modelling and statistics data for scientists and practitionersâthe case of Swiss lakes using Simstrat v2.1. Geosci. Model Dev. 12, 3955â3974 (2019).
Stepanenko, V. et al. LAKE 2.0: a model for temperature, methane, carbon dioxide and oxygen dynamics in lakes. Geosci. Model Dev. 9, 1977â2006 (2016).
Stepanenko, V., Repina, I. A., Ganbat, G. & Davaa, G. Numerical simulation of ice cover of saline lakes. Izv. Atmos. Ocean Phys. 55, 129â138 (2019).
Stepanenko, V. et al. A one-dimensional model intercomparison study of thermal regime of a shallow, turbid midlatitude lake. Geosci. Model Dev. 6, 1337â1352 (2013).
Stepanenko, V. et al. Simulation of surface energy fluxes and stratification of a small boreal lake by a set of one-dimensional models. Tellus A Dyn. Meteorol.Oceanogr. 66, 21389 (2014).
Thiery, W. et al. LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models. Tellus A Dyn. Meteorol. Oceanogr. 66, 21390 (2014).
Guseva, S. et al. Multimodel simulation of vertical gas transfer in a temperate lake. Hydrol. Earth Syst. Sci. 24, 697â715 (2020).
Bowling, L. C. & Lettenmaier, D. P. Modeling the effects of lakes and wetlands on the water balance of arctic environments. J. Hydrometeorol. 11, 276â295 (2010).
Gray, E. et al. Wide-spread inconsistency in estimation of lake mixed depth impacts interpretation of limnological processes. Water Res. 168, 115136 (2020).
Wilson, H. L. et al. Variability in epilimnion depth estimations in lakes. Hydrol. Earth Syst. Sci. 24, 5559â5577 (2020).
Ulloa, H. N., Winters, K. B., Wuest, A. & Bouffard, D. Differential heating drives downslope flows that accelerate mixed-layer warming in ice-covered waters. Geophys. Res. Lett. 46, 13872â13882 (2019).
Ramón, C.L., Ulloa, H. N., Doda, T., Winters, K. B. & Bouffard, D. Bathymetry and latitude modify lake warming under ice. Hydrol. Earth Syst. Sci. 25, 1813â1825 (2021).
Rose, K. C. et al. Annual precipitation regulates spatial and temporal drivers of lake water clarity. Ecol. Appl. 27, 632â643 (2017).
Rose, K. C., Winslow, L. A., Read, J. S. & Hansen, G. J. A. Climate-induced warming of lakes can be either amplified or suppressed by trends in water clarity. Limnol. Oceanogr. Lett. 1, 44â53 (2016).
Hocking, G. C. & StraÅ¡kraba, M. The effect of light extinction on thermal stratification in reservoirs and lakes. Int. Rev. Hydrobiol. 84, 535â556 (1999).
Heiskanen, J. J. et al. Effects of water clarity on lake stratification and lakeâatmosphere heat exchange. J. Geophys. Res. Atmos. 120, 7412â7428 (2015).
Shatwell, T., Adrian, R. & Kirillin, G. Planktonic events may cause polymicticâdimictic regime shifts in temperate lakes. Sci. Rep. 6, 24361 (2016).
Shatwell, T. et al. Future projections of temperature and mixing regime of European temperate lakes. Hydrol. Earth Syst. Sci. 23, 1533â1551 (2019).
Pilla, R. M. et al. Browning-related decreases in water transparency lead to long-term increases in surface water temperature and thermal stratification in two small lakes. J. Geophys. Res. Biogeosci. 123, 1651â1665 (2018).
Bartosiewicz, M. et al. Hot tops, cold bottoms: synergistic climate warming and shielding effects increase carbon burial in lakes. Limnol. Oceanogr. Lett. 4, 132â144 (2019).
Tanentzap, A. J. et al. Cooling lakes while the world warms: effects of forest regrowth and increased dissolved organic matter on the thermal regime of a temperate, urban lake. Limnol. Oceanogr. 53, 404â410 (2007).
Lei, H. Emergence of no-analogue conditions in lakes this century. Zenodo https://doi.org/10.5281/zenodo.11427869 (2024).
Acknowledgements
This work was supported by a UK Research and Innovation Natural Environment Research Council Independent Research Fellowship (NE/T011246/1) and a Natural Environment Research Council grant (reference number NE/X019071/1; UK Earth Observation Climate Information Service) (both awarded to R.I.W.), as well as the National Natural Science Foundation of China (number 42201049) and World Premier International Research Center Initiative of the Ministry of Education, Culture, Sports, Science and Technology, Japan (to K.B.R.). A.T. and S.-S.L. were supported by the Institute for Basic Science (project code IBS-R028-D1). L.H. was supported by the Second Tibetan Plateau Scientific Expedition and Research (grant number 2019QZKK0202). The CESM2-LE simulations presented here were conducted through a partnership between the Institute for Basic Science Center for Climate Physics in South Korea and the CESM group at the National Center for Atmospheric Research in the United States, representing a broad collaborative effort between scientists from both centres. The simulations were conducted on the Institute for Basic Science Center for Climate Physics supercomputer, Aleph. We also acknowledge support from KREONET. For their roles in producing, coordinating and making available the ISIMIP climate scenarios, we acknowledge support from the ISIMIP cross-sectoral science team. For producing the lake model simulations, we thank the ISIMIP Lake Sector modellers and ISIMIP Lake Sector coordinators.
Author information
Authors and Affiliations
Contributions
R.I.W. and L.H. initiated the project and led the data analysis and design of visualizations. L.H. provided the model projections. S.-S.L., K.B.R. and A.T. were instrumental in producing the large-scale simulations. L.H. and R.Y. developed the protocol for determining the emergence of no-analogue conditions. R.I.W. led the writing of the manuscript. All authors participated in discussions, revisions and the final production of this manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Geoscience thanks Dirk Schlabing and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Tamara Goldin, in collaboration with the Nature Geoscience team.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Figs. 1â14 and Supplementary Table 1.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Huang, L., Woolway, R.I., Timmermann, A. et al. Emergence of lake conditions that exceed natural temperature variability. Nat. Geosci. 17, 763â769 (2024). https://doi.org/10.1038/s41561-024-01491-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41561-024-01491-5