Hawliau Cyffredinol / General rightsCopyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.• Users may download and print one copy of any publication from the public portal for the purpose of private study or research.• You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal ? Take down policyIf you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

Heat uptake is a key variable for understanding the Earth system response to greenhouse gas forcing. Despite the importance of this heat budget, heat uptake by inland waters has so far not been quantified. Here we use a unique combination of global-scale lake models, global hydrological models and Earth system models to quantify global heat uptake by natural lakes, reservoirs, and rivers. The total net heat uptake by inland waters amounts to 2.6 ± 3.2 ×10 20 J over the period 1900-2020, corresponding to 3.6% of the energy stored on land. The overall uptake is dominated by natural lakes (111.7%), followed by reservoir warming (2.3%). Rivers contribute negatively (-14%) due to a decreasing water volume. The thermal energy of water stored in artificial reservoirs exceeds inland water heat uptake by a factor ∼10.4. This first quantification underlines that the heat uptake by inland waters is relatively small, but non-negligible. Plain Language Summary Human-induced emissions of CO 2 and other greenhouse gases cause energy accumulation in the Earth system. Oceans trap most of this excess energy, thereby largely buffering the warming of the atmosphere. However, the fraction of excess energy stored in lakes, reservoirs, and rivers is currently unknown, despite the high heat capacity of water. Here we quantify this human-induced heat storage, and show that it amounts up to 3.6% of the energy stored on land, while covering 2.58% of the land surface. The increase in heat storage from 1900 to 2020 is dominated by warming of lakes. The thermal heat contained in the water stored in man-made reservoirs is about ten times larger. Our study overall highlights the importance of inland waters-next to oceans, ice and land-for buffering atmospheric warming, especially on regional scale.

Abstract. Questions related to historical and future water resources and scarcity have been addressed by several macroscale hydrological models. One of these models is the Variable Infiltration Capacity (VIC) model. However, further model developments were needed to holistically assess anthropogenic impacts on global water resources using VIC. Our study developed VIC-WUR, which extends the VIC model using (1) integrated routing, (2) surface and groundwater use for various sectors (irrigation, domestic, industrial, energy, and livestock), (3) environmental flow requirements for both surface and groundwater systems, and (4) dam operation. Global gridded datasets on sectoral demands were developed separately and used as an input for the VIC-WUR model. Simulated national water withdrawals were in line with reported Food and Agriculture Organization (FAO) national annual withdrawals (adjusted R2 > 0.8), both per sector and per source. However, trends in time for domestic and industrial water withdrawal were mixed compared with previous studies. Gravity Recovery and Climate Experiment (GRACE) monthly terrestrial water storage anomalies were well represented (global mean root-mean-squared error, RMSE, values of 1.9 and 3.5 mm for annual and interannual anomalies respectively), whereas groundwater depletion trends were overestimated. The implemented anthropogenic impact modules increased simulated streamflow performance for 370 of the 462 anthropogenically impacted Global Runoff Data Centre (GRDC) monitoring stations, mostly due to the effects of reservoir operation. An assessment of environmental flow requirements indicates that global water withdrawals have to be severely limited (by 39 %) to protect aquatic ecosystems, especially with respect to groundwater withdrawals. VIC-WUR has potential for studying the impacts of climate change and anthropogenic developments on current and future water resources and sector-specific water scarcity. The additions presented here make the VIC model more suited for fully integrated worldwide water resource assessments.

Abstract. 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.

Abstract. 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. The protocol prescribes lake simulations driven by climate forcing from gridded observations and different Earth system models under various Representative Greenhouse Gas Concentration Pathways, 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 for nearly 17,500 lakes using uncalibrated models and forcing data from the global grid where lakes are present. 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.