The evaporative loss from global lakes (natural and artificial) is a critical component of the terrestrial water and energy balance. However, the evaporation volume of these water bodies—from the spatial distribution to the long-term trend—is as of yet unknown. Here, using satellite observations and modeling tools, we quantified the evaporation volume from 1.42 million global lakes from 1985 to 2018. We find that the long-term average lake evaporation is 1500 ± 150 km³ year⁻¹ and it has increased at a rate of 3.12 km³ year⁻¹. The trend attributions include an increasing evaporation rate (58%), decreasing lake ice coverage (23%), and increasing lake surface area (19%). While only accounting for 5% of the global lake storage capacity, artificial lakes (i.e., reservoirs) contribute 16% to the evaporation volume. Our results underline the importance of using evaporation volume, rather than evaporation rate, as the primary index for assessing climatic impacts on lake systems.

Surface water reservoirs are increasingly being relied upon to meet rising demands in the context of growing population and changing climate. However, the amount of water available in reservoirs (and the corresponding trends) have not been well quantified at the global scale. Here we use satellite observations to estimate the storage variations of 7245 global reservoirs from 1999 to 2018. Total global reservoir storage has increased at a rate of 27.82 ± 0.08 km³/yr, which is mainly attributed to the construction of new dams. However, the normalized reservoir storage (NS)—the ratio of the actual storage to the storage capacity—has declined by 0.82 ± 0.01%. The decline of NS values is especially pronounced in the global south, while the global north mainly exhibits an NS increase. With predicted decreasing runoff and increasing water demand, these observed diminishing storage returns of reservoir construction will likely persist into the future.

Knowledge of reservoir bathymetry is essential for many studies on terrestrial hydrological and biogeochemical processes. However, there are currently no cost-effective approaches to derive reservoir bathymetry at the global scale. This study explores the potential of generating high-resolution global bathymetry using elevation data collected by the 532-nm Advanced Topographic Laser Altimeter System (ATLAS) onboard the Ice, Cloud, and Land Elevation Satellite (ICESat-2). The novel algorithm was developed and tested using the ICESat-2 airborne prototype, the Multiple Altimeter Beam Experimental Lidar (MABEL), with Landsat-based water classifications (from 1982 to 2017). MABEL photon elevations were paired with Landsat water occurrence percentiles to establish the elevation-area (E-A) relationship, which in turn was applied to the percentile image to obtain partial bathymetry over the historic dynamic range of reservoir area. The bathymetry for the central area was projected to achieve the full bathymetry. The bathymetry image was then embedded onto the digital elevation model (DEM). Results were validated over Lake Mead against survey data. Results over four transects show coefficient of determination (R²) values from 0.82 to 0.99 and root-mean-square error (RMSE) values from 1.18 to 2.36 m. In addition, the E-A and elevation-storage (E-S) curves have RMSEs of 1.56 m and 0.08 km³, respectively. Over the entire dynamic reservoir area, the derived bathymetry agrees very well with independent survey data, except for within the highest and lowest percentile bands. With abundant overpassing tracks and high spatial resolution, the newly launched ICESat-2 should enable the derivation of bathymetry over an unprecedented number of reservoirs.

Area-volume-elevation (AVE) curves are critical for reservoir operation rules. However, such curves are not publicly available for most global reservoirs. Here, we present a framework to derive reservoir AVE curves from TanDEM-X data, using Lake Mead (~600 km ² ) as an example. First, the maximum water extent from 1984 to 2018-provided by the global surface water (GSW) dataset-was used as a mask to obtain the TanDEM-X data. Then, the TanDEM-X water indication mask (WAM) was applied to the extracted TanDEM-X data to obtain the visible bathymetry, which represents the topography between the maximum extent (according to GSW) and the water extent from WAM. Last, the AVE curve was generated by integrating the volume values from the top to bottom layers. TanDEM-X also captures the elevation values of the transitional waters, which are defined as the difference between the highest and lowest water levels. The transitional waters were obtained by thresholding amplitude and coherence images, and their elevations were then added to the visible bathymetry to extend the AVE curves with an elevation range extending from 344-369 m to 341-369 m. Validation results against in situ lidar survey values suggest a high-accuracy of elevation-area (E-A) relationships with R ² values of >0.99 and NRMSE values from 2.11% to 2.45%, and elevation-volume (E-V) relationships with R ² values of 1 and NRMSE values from 1.11% to 1.29%. Results also show that TanDEM-X data can capture the interannual variations due to multiple acquisitions, and that the elevation measurements for the lake shore areas are reliable.

The changing climate and intensifying human activities have made an impact on the hydrological processes in the upper Yangtze River (UYR), but quantifying their effects remains uncertain. This study used the Budyko framework to investigate the response of runoff (Q) to climate change and human activities during 1956–2017 and evaluate the impacts of human activities, including land use/cover change, water use, dam construction, and vegetation change, on watershed characteristic. Results show that climate change is the dominant driver of Q variations in the Wujiang River (WJR), Jialing River (JLR), and Jinsha River (JSR) watersheds, with contributions of 58.6%, 66.9%, and 67.6%, respectively. However, in Mingjiang River (MJR) and UYR watersheds, human activities contribute more to Q variations with 55.2% and 51.2%, respectively. Human activities play important roles in variation of watershed characteristics, and they can explain 22%, 26%, 36%, 25%, and 53% of the watershed character change in UYR, WJR, JLR, MJR, and JSR, respectively. This study conducts a comprehensive analysis of the causes of Q change in UYR, and provides a new perspective to explore the effects of specific human activities on watershed characteristics.