We comprehensively estimated water volume changes for 1132 lakes larger than 1 km2. Overall, the water mass stored in the lakes increased by 169.7±15.1 Gt (3.9±0.4 Gt yr-1) between 1976 and 2019, mainly in the Inner-TP (157.6±11.6 or 3.7±0.3 Gt yr-1). A substantial increase in mass occurred between 1995 and 2019 (214.9±12.7 Gt or 9.0±0.5 Gt yr-1), following a period of decrease (-45.2±8.2 Gt or -2.4±0.4 Gt yr-1) prior to 1995. A slowdown in the rate of water mass increase occurred between 2010 and 2015 (23.1±6.5 Gt or 4.6±1.3 Gt yr-1), followed again by a high value between 2015 and 2019 (65.7±6.7 Gt or 16.4±1.7 Gt yr-1). The increased lake-water mass occurred predominately in glacier-fed lakes (127.1±14.3 Gt) in contrast to non-glacier-fed lakes (42.6±4.9 Gt), and in endorheic lakes (161.9±14.0 Gt) against exorheic lakes (7.8±5.8 Gt) over 1976−2019.

Vegetation functional type (PFT) is a combination of large plant species according to the ecosystem function and resource utilization mode of plant species. Each planting functional type shares similar plant attributes, which simplifies the diversity of plant species into the diversity of plant function and structure.The concept of vegetation-functional has been advocated by ecologists especially ecosystem modelers.The basic assumption is that globally important ecosystem dynamics can be expressed and simulated through limited vegetative functional types.At present, vegetation-functional model has been widely used in biogeographic model, biogeochemical model, land surface process model and global dynamic vegetation model. For example, the land surface process model of the national center for atmospheric research (NCAR) in the United States has changed the original land cover information into the applied vegetation-functional map (Bonan et al., 2002).Functional vegetation has been used in the dynamic global vegetation model (DGVM) to predict the changes of ecosystem structure and function under the global change scenario. 1. Functional classification system of vegetation 1 Needleleaf evergreen tree, temperate 2 Needleleaf evergreen tree, boreal 3 Needleleaf deciduous tree 4 Broadleaf evergreen tree, tropical 5 Broadleaf evergreen tree, temperate 6 Broadleaf deciduous tree, tropical 7 Broadleaf deciduous tree, temperate 8 Broadleaf deciduous tree, boreal 9 Broadleaf evergreen shrub, temperate 10 Broadleaf deciduous shrub, temperate 11 Broadleaf deciduous shrub, boreal 12 C3 grass, arctic 13 C3 grass 14 C4 grass 15 Crop 16 Permanent wetlands 17 Urban and built-up lands 18 Snow and ice 19 Barren or sparsely vegetated lands 20 Bodies of water 2. Drawing method China's 1km vegetation function map is based on the climate rules of land cover and vegetation function conversion proposed by Bonan et al. (Bonan et al., 2002).Ran et al., 2012).MICLCover land cover map is a blend of 1:100000 data of land use in China in 2000, the Chinese atlas (1:10 00000) the type of vegetation, China 1:100000 glacier map, China 1:10 00000 marshes and MODIS land cover 2001 products (MOD12Q1) released the latest land cover data, using IGBP land cover classification system.The evaluation shows that it may be the most accurate land cover map on the scale of 1km in China.Climate data is China's atmospheric driven data with spatial resolution of 0.1 and temporal resolution of 3 hours from 1981 to 2008 developed by he jie et al. (2010).The data incorporates Princeton land-surface model driven data (Sheffield et al., 2006), gewex-srb radiation data (Pinker et al., 2003), TRMM 3B42 and APHRODITE precipitation data, and observations from 740 meteorological stations and stations under the China meteorological administration.According to the evaluation results of RanYouhua et al. (2010), GLC2000 has a relatively high accuracy in the current global land cover data set, and there is no mixed forest in its classification system. Therefore, the mixed forest in the MICLCover land cover diagram USES GLC2000 (Bartholome and Belward, 2005).The information in xu wenting et al., 2005) was replaced.The data can be used in land surface process model and other related researches.

The data set mainly includes the ice observation frequency (ICO) of north temperate lakes in four periods from 1985 to 2020, as well as the location, area and elevation of the lakes. Among them, the four time periods are 1985-1998 (P1), 1999-2006 (P2), 2007-2014 (P3) and 2015-2020 (P4) respectively, in order to improve the "valid observation" times in the calculation period and improve the accuracy. The ICO of the four periods is calculated by the ratio of "icing" times and "valid observation" times counted by all Landsat images in each period. Other lake information corresponds to the HydroLAKEs data set through the "hylak_id" column in the table. In addition, the data only retains about 30000 lakes with an area of more than 1 square kilometer, which are valid for P1-P4 observation. The data set can reflect the response of Lake icing to climate change in recent decades.

The dataset is a nearly 36-year (1983.7-2018.12) high-resolution (3 h, 10 km) global SSR (surface solar radiation) dataset, which can be used for hydrological modeling, land surface modeling and engineering application. The dataset was produced based on ISCCP-HXG cloud products, ERA5 reanalysis data, and MODIS aerosol and albedo products with an improved physical parameterization scheme. Validation and comparisons with other global satellite radiation products indicate that our SSR estimates were generally better than those of the ISCCP flux dataset (ISCCP-FD), the global energy and water cycle experiment surface radiation budget (GEWEX-SRB), and the Earth's Radiant Energy System (CERES). This SSR dataset will contribute to the land-surface process simulations and the photovoltaic applications in the future. The unit is W/㎡, instantaneous value.

The long-term sequence data set of lake areas on the Tibetan Plateau contains area data of 364 lakes with areas greater than 10 square kilometers from 1970s to 2013. Based on Landsat images, Landsat data in October are mainly used, and one data is taken every three years to reduce seasonal variation and make the available data reach the maximum. The data set is extracted by the NDWI Water Index, and each lake undergoes manual visual inspection and edition. The data set can be used to study lake change, lake water balance and climate change on the Tibetan Plateau. Data type: Vector data. Projection: WGS84.

Gridded climatic datasets with fine spatial resolution can potentially be used to depict the climatic characteristics across the complex topography of China. In this study we collected records of monthly temperature at 1153 stations and precipitation at 1202 stations in China and neighboring countries to construct a monthly climate dataset in China with a 0.025° resolution (~2.5 km). The dataset, named LZU0025, was designed by Lanzhou University and used a partial thin plate smoothing method embedded in the ANUSPLIN software. The accuracy of LZU0025 was evaluated based on three aspects: (1) Diagnostic statistics from the surface fitting model during 1951–2011. The results indicate a low mean square root of generalized cross validation (RTGCV) for the monthly air temperature surface (1.06 °C) and monthly precipitation surface (1.97 mm1/2). (2) Error statistics of comparisons between interpolated monthly LZU0025 with the withholding of climatic data from 265 stations during 1951–2011. The results show that the predicted values closely tracked the real true values with values of mean absolute error (MAE) of 0.59 °C and 70.5 mm, and standard deviation of the mean error (STD) of 1.27 °C and 122.6 mm. In addition, the monthly STDs exhibited a consistent pattern of variation with RTGCV. (3) Comparison with other datasets. This was done in two ways. The first was via comparison of standard deviation, mean and time trend derived from all datasets to a reference dataset released by the China Meteorological Administration (CMA), using Taylor diagrams. The second was to compare LZU0025 with the station dataset in the Tibetan Plateau. Taylor diagrams show that the standard deviation, mean and time trend derived from LZU had a higher correlation with that produced by the CMA, and the centered normalized root-mean-square difference for this index derived from LZU and CMA was lower. LZU0025 had high correlation with the Coordinated Energy and Water Cycle Observation Project (CEOP) - Asian Monsoon Project, (CAMP) Tibet surface meteorology station dataset for air temperature, despite a non-significant correlation for precipitation at a few stations. Based on this comprehensive analysis, we conclude that LZU0025 is a reliable dataset. LZU0025, which has a fine resolution, can be used to identify a greater number of climate types, such as tundra and subpolar continental, along the Himalayan Mountain. We anticipate that LZU0025 can be used for the monitoring of regional climate change and precision agriculture modulation under global climate change.

Based on the existing natural hole data of 15 active layer depth monitoring sites in the Qinghai-Tibet Engineering Corridor, the active layer depth distribution map of the Qinghai-Tibet Engineering Corridor was simulated using the GIPL2.0 frozen soil model. The model required synthesis of a temperature data set of time series. The temperature data were divided into two phases according to the time spans, which were 1980-2009 and 2010-2015. The data of the first phase were from the Chinese meteorological driving data set (http://dam. Itpcas.ac.cn/rs/?q=data#CMFD_0.1), and the data of the second phase was the application of MODIS surface temperature products (MOD11A1/A2 and MYD11A1/A2) with a spatial resolution of 1 km. In addition, the soil type data required by the model came from the China Soil Database (V1.1) and have a resolution of 1 km. At the same time, the topography was also considered. The research area was classified into 88 types based on the measured soil thermophysical parameters and land cover types, and then the simulation was performed. The simulation results were compared with the field measured data. The results showed that they were highly consistent, and the correlation coefficient reached 0.75. In alpine areas, the average depth of the active layer is below 2.0 m. However, in the river valleys, the average depth of the active layer is above 4.0 m. In the high plain area, the depth of the active layer is usually between 3.0 m and 4.0 m.

The data set was produced based on the SRTM DEM data collected by Space Shuttle Radar terrain mission in 2016, the reference data such as river, lake and other water system auxiliary data , using the arcgis hydrological model to analyze and extract the river network. There are 12 sub-basins over the Tibet Plateau, including AmuDayra、Brahmaputra、Ganges、Hexi、Indus、Inner、Mekong、Qaidam、Salween、Tarim、Yangtze、Yellow. The outer boundary is based on the 2500-metre contour line and national boundaries.

Soil bulk density, porosity, water content, water characteristic curve, saturated hydraulic conductivity, particle analysis, infiltration rate, and sampling point location information in the upper reaches of the Heihe River Basin. 1. The data is for 2014 supplementary sampling for 2012, using the ring knife to take the original soil; 2. The soil bulk density is the dry bulk density of the soil and is measured by the drying method. The original ring-shaped soil sample collected in the field was thermostated at 105 ° C for 24 hours in an oven, and the soil dry weight was divided by the soil volume (100 cubic centimeters) , unit: g/cm 3 . 3. Soil porosity is obtained according to the relationship between soil bulk density and soil porosity; 4. Soil infiltration analysis data set, the data is the field experimental measurement data from 2013 to 2014. 5. The infiltration data is measured by “MINI DISK PORTABLE TENSION INFILTROMETER”, and the approximate saturated hydraulic conductivity under a certain negative pressure is obtained. 6. Soil particle size data was measured at the Grain Granulation Laboratory of the Key Laboratory of the Ministry of Education of Lanzhou University. The measuring instrument is a Malvern laser particle size analyzer MS2000. 7. The saturated hydraulic conductivity is measured according to the enamel hair self-made instrument of Yi Yanli (2009). The Marioot bottle was used to maintain the head during the experiment; at the same time, the Ks measured at the time was converted to the Ks value at 10 °C for analysis and calculation. 8. Soil water content data is measured using ECH2O, including 5 layers of soil water content and soil temperature. 9. The water characteristic curve is measured by the centrifuge method: the undisturbed soil of the ring cutter collected in the field is placed in a centrifuge, and each of the speeds is measured at 0, 310, 980, 1700, 2190, 2770, 3100, 5370, 6930, 8200, 11600. The secondary rotor weight is obtained.

Mean annual ground temperature (MAGT) at a depth of zero annual amplitude and permafrost thermal stability type are fundamental importance for engineering planning and design, ecosystem management in permafrost region. This dataset is produced by integrating remotely sensed freezing degree-days and thawing degree-days, snow cover days, leaf area index, soil bulk density, high-accuracy soil moisture data, and in situ MAGT measurements from 237 boreholes for the 2010s (2005-2015) on the Tibetan Plateau (TP) by using an ensemble learning method that employs a support vector regression (SVR) model based on distance-blocked resampling training data with 200 repetitions. Validation of the new permafrost map indicates that it is probably the most accurate of all available maps at present. The RMSE of MAGT is approximately 0.75 °C and the bias is approximately 0.01 °C. This map shows that the total area of permafrost on the TP is approximately 115.02 (105.47-129.59) *104 km2. The areas corresponding to the very stable, stable, semi-stable, transitional, and unstable types are 0.86*104 km2, 9.62*104 km2, 38.45*104 km2, 42.29*104 km2, and 23.80*104 km2, respectively. This new dataset is available for evaluate the permafrost change in the future on the TP as a baseline. More details can be found in Ran et al., (2020) that published at Science China Earth Sciences.