A lake data set for the Tibetan Plateau from the 1960s, 2005, and 2014

A Comprehensive Dataset of Surface Water Quality Spanning 1940-2023 for Empirical and ML Adopted Research

Assessment and monitoring of surface water quality are essential for food security, public health, and ecosystem protection. Although water quality monitoring is a known phenomenon, little effort has been made to offer a comprehensive and harmonized dataset for surface water at the global scale. This study presents a comprehensive surface water quality dataset that preserves spatio-temporal variability, integrity, consistency, and depth of the data to facilitate empirical and data-driven evaluation, prediction, and forecasting. The dataset is assembled from a range of sources, including regional and global water quality databases, water management organizations, and individual research projects from five prominent countries in the world, e.g., the USA, Canada, Ireland, England, and China. The resulting dataset consists of 2.82 million measurements of eight water quality parameters that span 1940 - 2023. This dataset can support meta-analysis of water quality models and can facilitate Machine Learning (ML) based data and model-driven investigation of the spatial and temporal drivers and patterns of surface water quality at a cross-regional to global scale.

Dynamic evolution trend and driving mechanisms of water conservation in the Yellow River Basin, China

Water conservation (WC) is a critical ecological service function in the Yellow River Basin (YRB). There is currently a lack of detailed exploration of WC development processes and the impact mechanisms of driving factors at spatiotemporal scales in the YRB. By collecting data on DEM, land use, soil, meteorology, reservoirs, and observed discharge, this study established a large-scale WC model using the soil and water assessment tool (SWAT). The abrupt change test, empirical orthogonal function (EOF), wavelet analysis, hierarchical partitioning analysis (HPA), geodetectors, and aridity index were employed to analyze the multi-spatiotemporal characteristics and driving forces of WC calculated using the water balance method. The results are as follows: (1) The average WC among the YRB was 9.11 mm (74.68 × 108 m3) from 1960 to 2020. Pasture and forests contributed to 48.65% and 22.05% of the average annual WC, respectively. (2) WC exhibited four forms: less/more in the YRB, more in the southeast (northwest), and less in the northwest (southeast). (3) Forests and pastures in land use had higher average WC capacity, while Gansu, Shaanxi, and Qinghai ranked in the top three for average WC among the nine provinces. (4) Precipitation was the major driving force affecting WC variations, with the interaction between precipitation and actual evapotranspiration being the most significant. (5) Drought was a significant cause of negative WC. Protecting and managing crucial WC areas was essential for improving the ecological environment. This research elucidates the driving forces of WC in the YRB, providing scientific support for improving regional WC and promoting sustainable development.

Establishing a quantitative assessment methodology framework of water conservation based on the water balance method under spatiotemporal and different discontinuous ecosystem scales

Water conservation (WC) is an essential terrestrial ecosystem service that mitigates surface runoff and replenishes groundwater, which has received considerable attention under the dual pressures of climate change and human activity. However, there is insufficient understanding of the trends in WC changes on temporal (annual, monthly, daily), spatial, and ecosystem scales. This study proposed a quantitative assessment methodology framework (QAMF) for analyzing the spatiotemporal variation of WC under different discontinuous ecosystems. The QAMF mainly used models and methods such as the hydrological model (SWAT), calibration and uncertainty program (SWAT-CUP), WC calculation formula (water balance method), and spatial analysis method (empirical orthogonal function and wavelet analysis). It was applied to the source region of the Yellow River (SRYR), where the ecological landscape pattern underwent varying degrees of degradation, and WC capacity decreased. The results show that: Firstly, the constructed SWAT in the SRYR had high accuracy, and the proposed formula for calculating WC was suitable for multi-temporal scale analysis of WC in spatially distributed discontinuous basins. Secondly, the annual and monthly WC were respectively 81.00-184.13 mm and -28.58-107.64 mm, and daily WC was positive during extreme precipitation periods and negative during dry periods. The regulating effect of WC was fully reflected on the daily scale, partially reflected on the monthly scale, and absent on the annual scale. Third, the crucial WC area was mainly distributed in the southeast, and there was a significant primary yearly cycle of WC in the SRYR. Finally, different ecosystems exhibited different WC capabilities, and protecting the diversity of ecosystems played an essential role in maintaining and improving the WC function in the SRYR. This project has great scientific significance and technological support for scientifically evaluating the WC capacity in the SRYR.

Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets

In this study, multiple remote sensing data were used to quantitatively evaluate the contributions of surface water, soil moisture and groundwater to terrestrial water storage (TWS) changes in five groundwater resources zones of Inner Mongolia (GW_I, GW_II, GW_III, GW_IV and GW_V), China. The results showed that TWS increased at the rate of 2.14 mm/a for GW_I, while it decreased at the rate of 4.62 mm/a, 5.89 mm/a, 2.79 mm/a and 2.62 mm/a for GW_II, GW_III, GW_IV and GW_V during 2003-2021. Inner Mongolia experienced a widespread soil moisture increase with the rate of 4.17 mm/a, 2.13 mm/a, 1.20 mm/a, 0.25 mm/a and 1.36 mm/a for the five regions, respectively. Significant decreases were detected for regional groundwater storage (GWS) with the rate of 2.21 mm/a, 6.76 mm/a, 6.87 mm/a, 3.01 mm/a, and 4.14 mm/a, respectively. Soil moisture was the major contributor to TWS changes in GW_I, which accounted 58% of the total TWS changes. Groundwater was the greatest contributor to TWS changes in other four regions, especially GWS changes, which accounted for 76% TWS changes in GW_IV. In addition, this study found that the role of surface water was notable for calculating regional GWS changes.

Spatiotemporal patterns of organic carbon burial over the last century in Lake Qinghai, the largest lake on the Tibetan Plateau

Lakes are important carbon sinks in terrestrial environments. However, the estimation of the global lake carbon sink has large uncertainty. Data from plateau and remote lakes are rare, and most studies of carbon sequestration in large lakes have been based on single or a few sediment cores. Here, twenty-five sediment cores were collected by grid sampling covering Lake Qinghai, the largest lake on the Tibetan Plateau. Age models were established by combining radionuclide ²¹⁰Pb¹³⁷Cs dating with magnetic susceptibility chronostratigraphy of sediment cores. Furthermore, the spatiotemporal variations of the organic carbon burial rate (OCBR) over the past century were investigated. The OCBR of Lake Qinghai has increased significantly since the 1990s in association with warm-humid climates, increased nutrient supply and, enhanced land-use changes. The spatial distributions of OCBR were insignificant during the 1900s-1960s and 1960s-1990s and then shifted to a pattern of high values occurring in the southwestern lake areas during the post-1990s period. The spatial distribution of OCBR was mainly determined by sediment mass accumulation rate, primary production, and potential mineralization. The average OCBR in the all sediment cores showed no correlation with water depth (12-30 m) and was within one standard deviation of the whole-lake average value for most cores. These results suggest that the average OCBR of a sediment core in a relatively flat lakebed can generally represent the whole-basin level of Lake Qinghai. The average OCBR was 22.5 ± 5.5 g m^-2 yr^-1, which is close to those values reported previously for lakes of boreal forest and taiga but significantly higher than those reported for tundra lakes. Our findings highlight that the remote lakes on the Tibetan Plateau have great carbon sequestration potential in sediments and may act as a significant natural carbon sink.

Application and recent progress of inland water monitoring using remote sensing techniques

Hyperspectral remote sensing, which retrieves the water quality parameters by direct high-resolution analysis of the electromagnetic spectrum reflected from the water surface, has been widely applied for inland water quality detection. Such a new approach provides an opportunity to generate real-time data from water with the noncontact method, largely improving working efficiency. By summarizing the development and current applications of hyperspectral remote sensing, we compare the relative merits of varying remote sensing platforms, popular inversion models, and the application of hyperspectral monitoring of chlorophyll-a (Chl-a), transparency, total suspended solids (TSS), colored dissolved organic matter (CDOM), phycocyanin (PC), total phosphorus (TP), and total nitrogen (TN) water quality parameters. Most studies have focused on spaceborne remote sensing, which is usually used to monitor large waterbodies for Chl-a and other water quality parameters with optical properties; semiempirical, bio-optical, and semianalytical models are frequently used. With the rapid development of aerospace technology and near-surface remote sensing, the spectral resolution of remote sensing imaging technology has been dramatically improved and has begun to be applied to small waterbodies. In the future, the multiplatform linkage monitoring approach may become a new research direction. Advanced computer technology has also enabled machine learning models to be applied to water quality parameter inversion, and machine learning models have higher robustness than the three commonly used models mentioned above. Although nitrogen and phosphorus, with nonoptical properties, have also received attention and research from some scholars in recent years, the uncertainty of their mechanisms makes it necessary to maintain a cautious attitude when treating such research.

Lake volume variation in the endorheic basin of the Tibetan Plateau from 1989 to 2019

Lake storage change serves as a unique indicator of natural climate change on the Tibetan Plateau (TP). However, comprehensive lake storage data, especially for lakes smaller than 10 km2, are still lacking in the region. In this dataset, we completed a census of annual relative lake volume (RLV) for 976 lakes, which are larger than 1 km2, on the endorheic basin of the Tibetan Plateau (EBTP) during 1989-2019 using Landsat imagery and digital terrain models. Our method first identifies individual lakes, determines their analysis extents and calculates annual lake area from Landsat imagery. It then derives lake area-elevation relationship, estimates lake surface elevation, and calculates RLV. Validation and comparison with several existing datasets indicate our data are more reliable and comprehensive. Our study complements existing lake datasets by providing a complete and long-term lake water volume change data for the region.

Extraction of Water Body Information from Remote Sensing Imagery While Considering Greenness and Wetness Based on Tasseled Cap Transformation

Water, as an important part of ecosystems, is also an important topic in the field of remote sensing. Shadows and dense vegetation negatively affect most traditional methods used to extract water body information from remotely sensed images. As a result, extracting water body information with high precision from a wide range of remote sensing images which contain complex ground-based objects has proved difficult. In the present study, a method used for extracting water body information from remote sensing imagery considers the greenness and wetness of ground-based objects. Ground objects with varied water content and vegetation coverage have different characteristics in their greenness and wetness components obtained by the Tasseled Cap transformation (TCT). Multispectral information can be output as brightness, greenness, and wetness by Tasseled Cap transformation, which is widely used in satellite remote sensing images. Hence, a model used to extract water body information was constructed to weaken the influence of shadows and dense vegetation. Jiangsu and Anhui provinces are located along the Yangtze River, China, and were selected as the research area. The experiment used the wide-field-of-view (WFV) sensor onboard the Gaofen-1 satellite to acquire remotely sensed photos. The results showed that the contours and spatial extent of the water bodies extracted by the proposed method are highly consistent, and the influence of shadow and buildings is minimized; the method has a high Kappa coefficient (0.89), overall accuracy (92.72%), and user accuracy (88.04%). Thus, the method is useful in updating a geographical database of water bodies and in water resource management.

Geochemical and Statistical Analyses of Trace Elements in Lake Sediments from Qaidam Basin, Qinghai-Tibet Plateau: Distribution Characteristics and Source Apportionment

The safety of lake ecosystems on the Qinghai-Tibet Plateau (QTP) has attracted increasing attention, owing to its unique location and ecological vulnerability. Previous studies have shown that the aquatic systems on the QTP have been polluted to varying degrees by trace elements. However, little is known of the distribution and sources of trace elements in lakes in the northeast QTP. Here, 57 sediment samples were collected from six lakes (Dasugan Lake, Xiaoqaidam Lake, Kreuk Lake, Toson Lake, Gahai Lake and Xiligou Lake) in the Qaidam Basin, northeast QTP, and the trace elements (V, Cr, Ni, Cu, Zn, As, Ba, Tl, Cd, Pb, and U) were analyzed. The results indicated that Ba, Zn, V, and Cr had a higher content and a wider distribution relative to the other tested elements. Correlation coefficient matrix results showed that the trace elements in the study area were strongly correlated, revealing their source of similarity. Self-organizing maps (SOM, an artificial neural network algorithm) results indicated that the degree of pollution in Xiaoqaidam Lake was the highest, and that of Dasugan Lake was the lowest. Furthermore, all sampling points were clustered into four categories through K-means clustering. The positive matrix factorization (PMF) results indicated that atmospheric deposition and anthropogenic inputs were the main trace elements sources in these lakes, followed by traffic emission and geological sources. Collectively, trace elements of six lakes in Qinghai-Tibet Plateau presented high-content and pollution characteristics. This research provides a scientific basis for better water environment management and ecological protection on the QTP.

Monitoring Water Quality of the Haihe River Based on Ground-Based Hyperspectral Remote Sensing

The Haihe River is a typical sluice-controlled river in the north of China. The construction and operation of sluice dams change the flow and other hydrological factors of rivers, which have adverse effects on water, making it difficult to study the characteristics of water quality change and water environment control in northern rivers. In recent years, remote sensing has been widely used in water quality monitoring. However, due to the low signal-to-noise ratio (SNR) and the limitation of instrument resolution, satellite remote sensing is still a challenge to inland water quality monitoring. Ground-based hyperspectral remote sensing has a high temporal-spatial resolution and can be simply fixed in the water edge to achieve real-time continuous detection. A combination of hyperspectral remote sensing devices and BP neural networks is used in the current research to invert water quality parameters. The measured values and remote sensing reflectance of eight water quality parameters (chlorophyll-a (Chl-a), phycocyanin (PC), total suspended sediments (TSS), total nitrogen (TN), total phosphorus (TP), ammonia nitrogen (NH4-N), nitrate-nitrogen (NO3-N), and pH) were modeled and verified. The results show that the performance R2 of the training model is above 80%, and the performance R2 of the verification model is above 70%. In the training model, the highest fitting degree is TN (R2 = 1, RMSE = 0.0012 mg/L), and the lowest fitting degree is PC (R2 = 0.87, RMSE = 0.0011 mg/L). Therefore, the application of hyperspectral remote sensing technology to water quality detection in the Haihe River is a feasible method. The model built in the hyperspectral remote sensing equipment can help decision-makers to easily understand the real-time changes of water quality parameters.

Long-Term Lake Area Change and Its Relationship with Climate in the Endorheic Basins of the Tibetan Plateau

Lakes are sensitive indicators of climate change in the Tibetan Plateau (TP), which have shown high temporal and spatial variability in recent decades. The driving forces for the change are still not entirely clear. This study examined the area change of the lakes greater than 1 km2 in the endorheic basins of the Tibetan Plateau (EBTP) using Landsat images from 1990 to 2019, and analysed the relationships between lake area and local and large-scale climate variables at different geographic scales. The results show that lake area in the EBTP has increased significantly from 1990 to 2019 at a rate of 432.52 km2·year−1. In the past 30 years, lake area changes in the EBTP have mainly been affected by local climate variables such as precipitation and temperature. At a large scale, Indian Summer Monsoon (ISM) has correlations with lake area in western sub-regions in the Inner Basin (IB). While Atlantic Multidecadal Oscillation (AMO) has a significant connection with lake area, the North Atlantic Oscillation (NAO) does not. We also found that abnormal drought (rainfall) brought by the El Niño/La Niña events are significantly correlated with the lake area change in most sub-regions in the IB.

Simulation of the Dynamic Water Storage and Its Gravitational Effect in the Head Region of Three Gorges Reservoir Using Imageries of Gaofen-1

The construction of a high-resolution dynamic water storage model, driven by the mass load of the huge water storage of the Three Gorges Reservoir (TGR), is the necessary basic data for accurately simulating changes in the geophysical field, e.g., gravity, crustal deformation, and stress. However, previously established models cannot meet the needs of accurately simulating the impoundment effects of TGR, because these models were simplified and approximated and did not consider the variation of river boundaries caused by water level changes. In this study, we combined high-resolution Gaofen-1 (GF-1) satellite imageries and real-time water level in front of the dam and extracted 31 river boundaries of the head region of TGR between the lowest (145 m) and the highest (175 m) impoundment stages based on the Normalized Differential Water Index (NDWI) and threshold segmentation from Otsu method. Developed dynamic water storage model based on higher-resolution GF-1 data can show the true river boundary changes more exactly, especially in local areas. Compared to the previous approximate models, the model that we constructed accurately depicts the boundary distribution information of the different impoundment stages. Moreover, we simulated TGR-induced gravitational effects based on the high-precision forward modeling of the dynamic water storage model (i.e., considering changes of dynamic water area and water level). The theoretical modelled results are consistent with in situ gravity measurements with the difference mainly within 10 μGal. Our results indicate that water storage variations of TGR mainly affect the gravity field response within 1000 m of the reservoir bank with its maximum amplitude up to several hundred μGal. The dynamic water storage and its simulation results of gravitational effects can effectively eliminate the impact of surface water load driven by the TGR under human control and greatly improve the signal-to-noise ratio of regional gravity observational data. Thus, this work will be beneficial in the application of geophysical and geodetic monitoring aimed to comprehensively track the local and regional geological structural stability, e.g., artificial reservoir induced earthquake and landslide.

Recent Abnormal Hydrologic Behavior of Tibetan Lakes Observed by Multi-Mission Altimeters

Inland lakes in the Tibetan Plateau (TP) with closed catchments and minimal human disturbance are an important indicator of climate change. However, the examination of changes in the spatiotemporal patterns of Tibetan lakes, especially water level variations, is limited due to inadequate access to measurements. This obstacle has been improved by the development of satellite altimetry observations. The more recent studies revealed that the trend of central TP to grow decreased or reversed between 2010 and 2016. However, thus far, this trend has not been investigated to determine whether this pattern would last for the following years. This study aims to combine the traditional (launched before 2010, e.g., TOPEX/POSEIDON, ERS-1, ERS-2, Jason-1/-2, and Envisat) and recently advanced (launched after 2010, e.g., SARAL and Sentinel-3) altimetry observations to understand the Tibetan lake changes further in recent years. Therefore, we acquired information on the continuous lake level changes in Tibetan lakes using the lake level sequence integration method based on multisource altimetry satellites. The results revealed that water level changes in 22 examined lakes showed abrupt rises in 2016–2018, but the onsets and magnitudes of the rises varied among the lakes. During the study period, the water levels of the lakes (except Nam Co) revealed a drastic rising tendency with a mean rate of 0.74 m/a, which was remarkably higher than the average rate of water level rise over the period 2010–2015 (approximately 0.28 m/a). Specifically, the water level of the nine lakes in the Northern TP (NTP) displayed a significant rising trend, with an average rate of 0.82 m/a. In the Central TP (CTP), the lake level changes were generally divided into two categories. The water levels for the lakes in the Western CTP rose rapidly, while, in the Eastern CTP, the lake water levels rose slowly, with an average rising rate less than 0.40 m/a. The water levels for the lakes in the Northeastern TP (NETP) and Northwestern TP (NWTP) kept a stable rising tendency. According to the results of the climate analysis, the spatial differences of the lake level rise rates were primarily caused by the spatial and temporal changes of precipitation over the TP.

Achieving Higher Resolution Lake Area from Remote Sensing Images Through an Unsupervised Deep Learning Super-Resolution Method

Lakes have been identified as an important indicator of climate change and a finer lake area can better reflect the changes. In this paper, we propose an effective unsupervised deep gradient network (UDGN) to generate a higher resolution lake area from remote sensing images. By exploiting the power of deep learning, UDGN models the internal recurrence of information inside the single image and its corresponding gradient map to generate images with higher spatial resolution. The gradient map is derived from the input image to provide important geographical information. Since the training samples are only extracted from the input image, UDGN can adapt to different settings per image. Based on the superior adaptability of the UDGN model, two strategies are proposed for super-resolution (SR) mapping of lakes from multispectral remote sensing images. Finally, Landsat 8 and MODIS (moderate-resolution imaging spectroradiometer) images from two study areas on the Tibetan Plateau in China were used to evaluate the performance of UDGN. Compared with four unsupervised SR methods, UDGN obtained the best SR results as well as lake extraction results in terms of both quantitative and visual aspects. The experiments prove that our approach provides a promising way to break through the limitations of median-low resolution remote sensing images in lake change monitoring, and ultimately support finer lake applications.

Black carbon and mercury in the surface sediments of Selin Co, central Tibetan Plateau: Covariation with total carbon

Tibetan Plateau (TP) is an important geographical region for investigating the long-range transport of pollutants as limited emission sources exist in this region. In this study, based on analysis of 61 surface samples, we report the spatial distribution and concentrations of BC, Hg, total organic carbon (TOC) and inorganic carbon (IC) in surface sediments of Selin Co, the largest lake in central Tibet. The mean BC and Hg concentrations were 0.62 ± 0.34 mg/g and 32.03 ± 9.88 ng/g (range: 0.03-1.47 mg/g and 13.83-51.81 ng/g respectively), which were lower than the values from other lakes in the Himalayan-Tibetan Plateau (HTP). BC and Hg exhibited similar spatial distribution in the surface sediments. Similarly, the mean TOC and IC were 2.19 ± 1.46% and 3.13 ± 1.07% (range: 0.0007-7.78% and 0.30-5.30% respectively). BC/TOC ratio, as well as char/soot ratio, suggests biomass burning as a major source of BC in the sediments via the influence of long-range transport. The positive correlation between the concentrations of BC and Hg suggests similar emission sources or transport pathway. Concentrations of BC and Hg were higher in fine grain particles (size <~50 μm) which were capable of transport and deposit in the deeper part of the lake, as suggested by a significant relationship between water depth and particle size. This study elucidates the extent of pollution in very recent ages and also could serve as the basis for paleo-environmental studies in future.

Subpixel Mapping of Surface Water in the Tibetan Plateau with MODIS Data

This article presents a comprehensive subpixel water mapping algorithm to automatically produce routinely open water fraction maps in the Tibetan Plateau (TP) with the Moderate Resolution Imaging Spectroradiometer (MODIS). A multi-index threshold endmember extraction method was applied to select the endmembers from MODIS images. To incorporate endmember variability, an endmember selection strategy, called the combined use of typical and neighboring endmembers, was adopted in multiple endmember spectral mixture analysis (MESMA), which can assure a robust subpixel water fractions estimation. The accuracy of the algorithm was assessed at both the local scale and regional scale. At the local scale, a comparison using the eight pairs of MODIS/Landsat 8 Operational Land Imager (OLI) water maps demonstrated that subpixels water fractions were well retrieved with a root mean square error (RMSE) of 7.86% and determination coefficient (R2) of 0.98. At the regional scale, the MODIS water fraction map in October 2014 matches well with the TP lake data set and the Global Lake and Wetland Database (GLWD) in both latitudinal and longitudinal distribution. The lake area estimation is more consistent with the reference TP lake data set (difference of −3.15%) than the MODIS Land Water Mask (MOD44W) (difference of −6.39%).

Estimation of satellite-derived lake water surface temperatures in the western Mediterranean: Integrating multi-source, multi-resolution imagery and a long-term field dataset using a time series approach

Lake surface water temperature (LSWT) is a key parameter to help study the environmental and ecological impacts of climate change. In this work, we measured the LSWT of 1 natural and 23 artificial lakes located on the island of Sardinia in the western Mediterranean, which is a region where changes in climate are projected to have significant impacts. By integrating multi-source and multi-resolution datasets of the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat and long-term in situ temperature observations, we detected, measured, and analysed the LSWT trends during the period of 2000-2018 across all the investigated lakes. Methodologically, we demonstrated that a simplified approached based on Planck's equation for Landsat thermal infrared (TIR) data could be a valid alternative to radiative transfer equation retrieval methods for the retrieval of LSWT without loss of accuracy. Moreover, we demonstrated that rescaled and independently validated MOD112A-derived LSWT showed good accuracy, efficiently filled the spatial and temporal gaps in long-term in situ LSWT, and could be used for long-term LSWT trend detection and measurement. All 24 lakes showed an annual warming trend of +0.010 °C/y, warming winter trend of +0.013 °C/y, and cooling summer trend of -0.038 °C/y during the period of 2000-2018. This study demonstrated that the measured trend rates could be explained by and were strongly correlated with the climatology of Italy for the 2000-2018 period. Finally, we demonstrated the key role and the importance of the availability of long-term in situ temperature datasets. The approach used in this study is up-scalable to other medium to low-resolution TIR sensors as well as to other long-term monitoring sites, such as LTER-Italy, LTER-Europe, or ILTER sites.

Reconstructing Six Decades of Surface Temperatures at a Shallow Lake

Lake surface water temperature (LSWT) plays a fundamental role in the lake energy budget. However, direct observations of LSWT require considerable effort for acquisition and hence are rare relative to a large number of lakes. In lakes where LSWT has not been covered sufficiently by in situ measurements, remote sensing and lake modeling can be used to produce a fine spatio-temporal record of LSWTs. In our study, the Moderate-Resolution Imaging Spectroradiometer (MODIS) LSWT was used to compare with in situ data at the overpass times over the six sites in Lake Chaohu, a large shallow lake in China. MODIS-derived LSWT reflected the variation of lake surface temperature well, with a correlation coefficient of 0.96 and a cool bias of 1.25 °C. The bias was modified by an “Upper Envelop” smoothing method and then employed to evaluate the general lake model (GLM) performance, a one-dimensional hydrodynamic model. The GLM simulations showed good performance compared with MODIS LSWT data at an interannual time scale. A 57-year record of simulated LSWT was hindcast by the well-calibrated GLM for Lake Chaohu. The results showed that LSWT decreased by 0.08 °C/year from 1960 to 1981 and then increased by 0.05 °C/year. These trends were most likely caused by a cooling effect of decreased surface incident solar radiation and a warming effect of reduced wind speed. Our study promoted the use of MODIS-derived LSWT as an alternative data source, and then combined with a numerical model for inland water surface temperature, and also further provided an understanding of climate warming effect on such a shallow eutrophic lake. Key points: (1) Moderate-Resolution Imaging Spectroradiometer (MODIS) lake water surface temperature (LSWT) was validated with real-time in situ data collected at Lake Chaohu with high accuracy; (2) MODIS LSWT was modified by the bias correction and employed to evaluate a one-dimensional lake model at interannual and intraannual scale; The LSWT hindcast by a well-calibrated model at Lake Chaohu decreased by 0.08 °C/year from 1960 to 1981 and increased by 0.05 °C/year from 1982 to 2016.

Climate change leads to a doubling of turbidity in a rapidly expanding Tibetan lake

Recent climate change is causing most lakes on the Tibetan Plateau to grow at an unprecedented rate. Changes in the physical properties and water storage of the lakes are now relatively well documented. Yet the impacts on their water quality remain poorly understood. Turbidity is a well-established optical water-quality indicator related to suspended particulate matter concentration which can affect vertical light attenuation and ecosystem functioning. Here, we use remotely sensed data to assess the seasonal and long-term variations in turbidity in Siling Lake, one of the fastest growing lakes on the Tibetan Plateau, and to identify potential driving mechanisms of this change. The lake experiences two distinct peaks of turbidity during the year: one in August (warm season) caused by the seasonal influx of sediments from the Zagya Zangbo River, and one in December (cold season) caused by the wind-driven resuspension of sediments along the lakes' shorelines. The analysis further revealed a persistent increasing trend that doubled the average lake turbidity between 2000 and 2017. Evidence suggests this rise in turbidity results from a climate-driven increase in sediment supply from the Zagya Zangbo River, and from sediment resuspension associated with the erosion of shorelines recently submerged during the rapid expansion of the lake (paleoshorelines). Our results highlight the vulnerability of the Tibetan Lakes' water quality to climate change.

A Two-Stage Fusion Framework to Generate a Spatio–Temporally Continuous MODIS NDSI Product over the Tibetan Plateau

The Tibetan Plateau (TP) is an important component of the global environmental system, on which the snow cover greatly affects the regional climate and ecology. Moderate resolution imaging spectroradiometer (MODIS) snow cover products have been demonstrated to be appropriate for investigating the snow cover over the TP. However, they are subject to cloud obscuration, and the TP’s extremely complex terrain makes the snow monitoring difficult. Therefore, in this paper, we propose a two-stage spatio–temporal fusion framework for the cloud removal of MODIS C6 snow products, including an adjusted Terra and Aqua combination (TAC) and a spatio–temporal fusion based on Gaussian kernel function and error correction (STF-GKF-EC). To the best of our knowledge, this is the first time that a spatio–temporally continuous daily 500-m MODIS normalized difference snow index (NDSI) product has been generated for the TP, which greatly improves the spatial and temporal resolutions of the current snow cover products. The main stage, STF-GKF-EC, adaptively weights the spatial and temporal correlations by the Gaussian kernel function, and further takes the rapid changes of snow cover into consideration through the error correction. The experiments indicated that STF-GKF-EC removes clouds completely, achieving an overall accuracy (OA) and mean absolute error (MAE) of 91.48% and 3.88, respectively. Based on the cloud-removed results, during 2001–2017, as far as the intra-annual variation is concerned, a large proportion of the snow cover appears between October and May, with a peak in February/March, and the variation is mainly controlled by temperature. For the inter-annual variation, an obvious increasing trend of 0.68/year for NDSI is observed before 2005, followed by a slight decreasing trend of 0.16/year, in which precipitation is a better explanation factor than temperature.

Assessment of Water Storage Change in China’s Lakes and Reservoirs over the Last Three Decades

Lakes and reservoirs are essential elements of the hydrological and biochemical cycles, considered sentinels of global climate change. However, comprehensive quantifications of their water storage changes (∆V) at a large spatiotemporal scale are still rare. Here, we integrated a global surface water dataset and SRTM digital elevation models, both available from Google Earth Engine platform at a spatial resolution of 30 m, to evaluate ∆V for a total of 760 lakes and reservoirs across China at an annual timescale since 1984. The results indicated that (1) the aggregated water storage went through a slight increase of 41.5 Gt (1.7 Gt/yr) during 1985–2005, a significant decrease of 100 Gt (-20.6 Gt/yr) during 2005–2009, and then increased by 136.3 Gt (21.3 Gt/yr) during 2009–2015. (2) The increasing trend was largely attributed to lakes and reservoirs in the Tibetan Plateau Lake Zone, and the decreasing trend was mainly due to the North and Northwest Lake Zone, with little variations observed for the Northeast and Southwest Lake Zones. (3) Qinghai lake was associated with the largest increase (18.3 Gt) and Poyang lake presented the largest decline (-9.2 Gt). The results can help advance our understanding of the impact of climate change and improve future projection.

A New Digital Lake Bathymetry Model Using the Step-Wise Water Recession Method to Generate 3D Lake Bathymetric Maps Based on DEMs

The availability of lake bathymetry maps is imperative for estimating lake water volumes and their variability, which is a sensitive indicator of climate. It is difficult, if not impossible, to obtain bathymetric measurements from all of the thousands of lakes across the globe due to costly labor and/or harsh topographic regions. In this study, we develop a new digital lake bathymetry model (DLBM) using the step-wise water recession method (WRM) to generate 3-dimensional lake bathymetric maps based on the digital elevation model (DEM) alone, with two assumptions: (1) typically, the lake’s bathymetry is formed and shaped by geological processes similar to those that shaped the surrounding landmasses, and (2) the agent rate of water (the thickness of the sedimentary deposit proportional to the lake water depth) is uniform. Lake Ontario and Lake Namco are used as examples to demonstrate the development, calibration, and refinement of the model. Compared to some other methods, the estimated 3D bathymetric maps using the proposed DLBM could overcome the discontinuity problem to adopt the complex topography of lake boundaries. This study provides a mathematically robust yet cost-effective approach for estimating lake volumes and their changes in regions lacking field measurements of bathymetry, for example, the remote Tibetan Plateau, which contains thousands of lakes.

A data set of inland lake catchment boundaries for the Qiangtang Plateau

A catchment is the basic unit for studying hydrologic cycle processes and associated climate change impacts. Accurate catchment delineation is essential in the field of hydrology, environment, and meteorology. Traditionally, catchment delineation is most easily carried out where the outflow area can be easily determined because of a well-defined outlet. The obstacle of the current study is to determine accurately the catchment boundary of lakes that are internally draining and, therefore, lack a well-defined outflow (i.e. inland lakes). This study describes a catchment delineation method which demarcated all the catchments of the lakes in the Qiangtang Plateau, especially for the inland lakes and their closed catchments. Lake catchment boundaries determined for the Qiangtang Plateau provide a significant advancement for water resource and climate change evaluation and agriculture production in the area.

Water Storage Variations in Tibet from GRACE, ICESat, and Hydrological Data

The monitoring of water storage variations is essential not only for the management of water resources, but also for a better understanding of the impact of climate change on hydrological cycle, particularly in Tibet. In this study, we estimated and analyzed changes of the total water budget on the Tibetan Plateau from the Gravity Recovery And Climate Experiment (GRACE) satellite mission over 15 years prior to 2017. To suppress overall leakage effect of GRACE monthly solutions in Tibet, we applied a forward modeling technique to reconstruct hydrological signals from GRACE data. The results reveal a considerable decrease in the total water budget at an average annual rate of −6.22 ± 1.74 Gt during the period from August 2002 to December 2016. In addition to the secular trend, seasonal variations controlled mainly by annual changes in precipitation were detected, with maxima in September and minima in December. A rising temperature on the plateau is likely a principal factor causing a continuous decline of the total water budget attributed to increase melting of mountain glaciers, permafrost, and snow cover. We also demonstrate that a substantial decrease in the total water budget due to melting of mountain glaciers was partially moderated by the increasing water storage of lakes. This is evident from results of ICESat data for selected major lakes and glaciers. The ICESat results confirm a substantial retreat of mountain glaciers and an increasing trend of major lakes. An increasing volume of lakes is mainly due to an inflow of the meltwater from glaciers and precipitation. Our estimates of the total water budget on the Tibetan Plateau are affected by a hydrological signal from neighboring regions. Probably the most significant are aliasing signals due to ground water depletion in Northwest India and decreasing precipitation in the Eastern Himalayas. Nevertheless, an integral downtrend in the total water budget on the Tibetan Plateau caused by melting of glaciers prevails over the investigated period.

A long-term dataset of lake surface water temperature over the Tibetan Plateau derived from AVHRR 1981-2015

Lake surface water temperature (LSWT) is of vital importance for hydrological and meteorological studies. The LSWT ground measurements in the Tibetan Plateau (TP) were quite scarce because of its harsh environment. Thermal infrared remote sensing is a reliable way to calculate historical LSWT. In this study, we present the first and longest 35-year (1981-2015) daytime lake-averaged LSWT data of 97 large lakes (>80 km2 each) in the TP using the 4-km Advanced Very High Resolution Radiometer (AVHRR) Global Area Coverage (GAC) data. The LSWT dataset, taking advantage of observations from NOAA's afternoon satellites, includes three time scales, i.e., daily, 8-day-averaged, and monthly-averaged. The AVHRR-derived LSWT has a similar accuracy (RMSE = 1.7 °C) to that from other data products such as MODIS (RMSE = 1.7 °C) and ARC-Lake (RMSE = 2.0 °C). An inter-comparison of different sensors indicates that for studies such as those considering long-term climate change, the relative bias of different AVHRR sensors cannot be ignored. The proposed dataset should be, to some extent, a valuable asset for better understanding the hydrologic/climatic property and its changes over the TP.

Monitoring 40-Year Lake Area Changes of the Qaidam Basin, Tibetan Plateau, Using Landsat Time Series

Areal changes of high-altitude inland lakes on the Qaidam Basin (QB) of the Tibetan Plateau are reliable indicators of climate change and anthropogenic disturbance. Due to the physical difficulties to access, our knowledge of the spatial patterns and processes of climatic and human impacts on the Basin has been limited. Focusing on lake area changes, this study used long-term Landsat images to map the patterns of lakes and glaciers in 1977, 1990, 2000, and 2015, and to monitor the spatially explicit changes of lakes between 1977 and 2015. Results revealed that the total number of lakes (area > 0.5 km2) increased by 18, while their total area expanded by 29.8%, from 1761.5 ± 88.1 km2 to 2285.9 ± 91.4 km2. Meanwhile, glaciers have decreased in area by 259.16 km2 in the past four decades. The structural equation model (SEM) was applied to examine the integrative effects of natural and anthropogenic factors on lake area. Precipitation change exhibited the most significant influence on lake area in the QB from 1977 to 2000, while human activities also played an important role in the expansion of lakes in the QB in the period 2000–2015. In particular, extensive exploitation of salt lakes as mining resources resulted in severe changes in lake area and landscape. The continuously expanding salt lakes inundated the road infrastructure nearby, posing great threats to road safety. This study shed new light on the impacts of recent environmental changes and human interventions on lakes in the Qaidam Basin, which could assist policy-making for protecting the lakes and for strengthening the ecological improvement of this vast, arid basin.

A Modified Empirical Retracker for Lake Level Estimation Using Cryosat-2 SARin Data

Satellite radar altimetry is an important technology for monitoring water levels, but issues related to waveform contamination restrict its use for rivers, narrow reservoirs, and small lakes. In this study, a novel and improved empirical retracker (ImpMWaPP) is presented that can derive stable inland lake levels from Cryosat-2 synthetic aperture radar interferometer (SARin) waveforms. The retracker can extract a robust reference level for each track to handle multi-peak waveforms. To validate the lake levels derived by ImpMWaPP, the in situ gauge data of seven lakes in the Tibetan Plateau are used. Additionally, five existing retrackers are compared to evaluate the performance of the proposed ImpMWaPP retracker. The results reveal that ImpMWaPP can efficiently process the multi-peak waveforms of the Cryosat-2 SARin mode. The root-mean-squared errors (RMSEs) obtained by ImpMWaPP for Qinghai Lake, Nam Co, Zhari Namco, Ngoring Lake, Longyangxia Reservoir, Bamco, and Dawa Co are 0.085 m, 0.093 m, 0.109 m, 0.159 m, 0.573 m, 0.087 m, and 0.122 m, respectively. ImpMWaPP obtains the lowest mean RMSE (0.175 m) over the seven lakes, indicating that it extracts lake levels well during icing and no-ice periods, and is more suitable for lakes frozen in winter.

A comprehensive data set of lake surface water temperature over the Tibetan Plateau derived from MODIS LST products 2001-2015

Lake surface water temperature (LSWT) is sensitive to long-term changes in thermal structure of lakes and regional air temperature. In the context of global climate change, recent studies showed a significant warming trend of LSWT based on investigating 291 lakes (71% are large lakes, ≥50 km² each) globally. However, further efforts are needed to examine variation in LSWT at finer regional spatial and temporal scales. The Tibetan Plateau (TP), known as ‘the Roof of the World’ and ‘Asia’s water towers’, exerts large influences on and is sensitive to regional and even global climates. Aiming to examine detailed changing patterns and potential driven mechanisms for temperature variations of lakes across the TP region, this paper presents the first comprehensive data set of 15-year (2001–2015) nighttime and daytime LSWT for 374 lakes (≥10 km² each), using MODIS (Moderate Resolution Imaging Spectroradiometer) Land Surface Temperature (LST) products as well as four lake boundary shapefiles (i.e., 2002, 2005, 2009, and 2014) derived from Landsat/CBERS/GaoFen-1 satellite images. The data set itself reveals significant information on LSWT and its changes over the TP and is an indispensable variable for numerous applications related to climate change, water budget analysis (particularly lake evaporation), water storage changes, glacier melting and permafrost degradation, etc.