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.

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.

The past frozen soil map of the Tibetan Plateau was based on a small number of temperature station observations and used a classification system based on continuity. This data set used the geographically weighted regression model (GWR) to synthesize MODIS surface temperature, leaf area index, snow cover ratio and multimodel soil moisture forecast products of the National Meteorological Information Center through spatiotemporal reconstruction. In addition, precipitation observations of more than 40 meteorological stations, the precipitation products of FY2 satellite observations and the multiyear average temperature observation data of 152 meteorological stations from 2000 to 2010 were integrated to simulate the average temperature data of the Tibetan Plateau, and the permafrost thermal condition classification system was used to classify permafrost into several types: Very cold, Cold, Cool, Warm, Very warm, and Likely thawing. The map shows that, after deducting lakes and glaciers, the total area of permafrost on the Tibetan Plateau is approximately 1,071,900 square kilometers. Verification shows that this map has higher accuracy. It can provide support for future planning and design of frozen soil projects and environmental management.

Qinghai Tibet Plateau is the largest permafrost area in the world. At present, some permafrost distribution maps have been compiled. However, due to the limited data sources, unclear standards, insufficient verification and lack of high-quality spatial data sets, there is great uncertainty in drawing Permafrost Distribution Maps on TP. Based on the improved medium resolution imaging spectrometer (MODIS) surface temperature (LSTS) model of 1 km clear sky mod11a2 (Terra MODIS) and myd11a2 (Aqua MODIS) product (reprocessing version 5) in 2003-2012, the data set simulates the distribution of permafrost and generates the permafrost map of Qinghai Tibet Plateau. The map was verified by field observation, soil moisture content and bulk density. Permafrost attributes mainly include: seasonally frozen ground, permafrost and unfrozen ground. The data set provides more detailed data of Permafrost Distribution and basic data for the study of permafrost in the Qinghai Tibet Plateau.

This dataset contains five types of boundaries. 1. TPBoundary_ 2500m: Based on ETOPO5 Global Surface Relief, ENVI+IDL was used to extract data at an elevation of 2500m within the longitude (65~105E) and latitude (20~45N) range in the Tibetan Plateau. 2. TPBoundary_ 3000m: Based on ETOPO5 Global Surface Relief, ENVI+IDL was used to extract data at an elevation of 3000m within the longitude (65~105E) and latitude (20~45N) range in the Tibetan Plateau. 3. TPBoundary_ HF (high_frequency): This boundary is defined according to 2 previous studies. Bingyuan Li (1987) had a systematic discussion on the principles for determining the extent of the Tibetan Plateau and the specific boundaries. From the perspective of the formation and basic characteristics of the Tibetan Plateau, he proposed the basic principles for determining the extent of the Tibetan Plateau based on the geomorphological features, the plateau surface and its altitude, while considering the integrity of the mountain. Yili Zhang (2002) determined the extent and boundaries of the Tibetan Plateau based on the new results of research in related fields and years of field practice. He combined information technology methods to precisely locate and quantitatively analyze the extent and boundary location of the Tibetan Plateau, and concluded that the Tibetan Plateau in China extends from the Pamir Plateau in the west to the Hengduan Mountains in the east, from the southern edge of the Himalayas in the south to the northern side of the Kunlun-Qilian Mountains in the north. On April 14, 2017, the Ministry of Civil Affairs of the People's Republic of China issued the Announcement on Adding Geographical Names for Public Use in the Southern Tibetan Region (First Batch), adding six geographical names in the southern Tibetan region, including Wo’gyainling, Mila Ri, Qoidêngarbo Ri, Mainquka, Bümo La, and Namkapub Ri. ４. TPBoundary_ New (2021): Along with the in-depth research on the Tibetan Plateau, the improvement of multidisciplinary research and understanding inside and outside the plateau, and the progress of geographic big data and Earth observation science and technology, the development of the 2021 version of the Tibetan Plateau boundary data by Yili Zhang and et al. was completed based on the comprehensive analysis of ASTER GDEM and Google Earth remote sensing images. The range boundary starts from the northern foot of the West Kunlun Mountain-Qilian Mountain Range in the north and reaches the southern foot of the Himalayas and other mountain ranges in the south, with a maximum width of 1,560 km from north to south; from the western edge of the Hindu Kush Mountains and the Pamir Plateau in the west to the eastern edge of the Hengduan Mountains and other mountain ranges in the east, with a maximum length of about 3,360 km from east to west; the latitude and longitude range is 25°59′30″N~40°1′0″N, 67°40′37″E~104°40′57″E, with a total area of 3,083,400km2 and an average altitude of about 4,320m. Administratively, the Tibetan Plateau is distributed in nine countries, including China, India, Pakistan, Tajikistan, Afghanistan, Nepal, Bhutan, Myanmar, and Kyrgyzstan. ５. TPBoundary_ Rectangle: The rectangle was drawn according to the range of Lon (63~105E) and Lat (20~45N). The data are in latitude and longitude projection WGS84. As the basic data, the boundary of the Tibetan Plateau can be used as a reference basis for various geological data and scientific research on the Tibetan Plateau.

The China Meteorological Forcing Dataset (CMFD) is a high spatial-temporal resolution gridded near-surface meteorological dataset that was developed specifically for studies of land surface processes in China. The dataset was made through fusion of remote sensing products, reanalysis dataset and in-situ observation data at weather stations. Its record starts from January 1979 and keeps extending (currently up to December 2018) with a temporal resolution of three hours and a spatial resolution of 0.1°. Seven near-surface meteorological elements are provided in CMFD, including 2-meter air temperature, surface pressure, specific humidity, 10-meter wind speed, downward shortwave radiation, downward longwave radiation and precipitation rate.

This data is originated from the 1:100,000 national basic geographic database, which was open freely for public by the National Basic Geographic Information Center in November 2017. The boundary of the Qinghai-Tibet Plateau was spliced and clipped as a whole, so as to facilitate the study on the Qinghai-Tibet plateau. This data set is the 1:100,000 administrative boundaries of the qinghai-tibet plateau, including National_Tibet_line、 Province_Tibet、City_Tibet、County_Tibet_poly and County_Tibet_line. Administrative boundary layer (County_Tibet_poly) property name and definition: Item Properties Describe Example PAC Administrative division code 513230 NAME The name of the County line name Administrative boundary layer (BOUL) attribute name and definition: Item Properties Describe Example GB classification code 630200 Administrative boundary layer (County_Tibet_line) attribute item meaning: Item Properties Describe Example GB 630200 Provincial boundary GB 640200 Prefectural, municipal and state administrative boundaries GB 650201 county administrative boundaries (determined)

This data set comprises the plateau soil moisture and soil temperature observational data based on the Tibetan Plateau, and it is used to quantify the uncertainty of model products of coarse-resolution satellites, soil moisture and soil temperature. The observation data of soil temperature and moisture on the Tibetan Plateau (Tibet-Obs) are from in situ reference networks at four regional scales, which are the Nagqu network of cold and semiarid climate, the Maqu network of cold and humid climate, and the Ali network of cold and arid climate，and Pali network. These networks provided representative coverage of different climates and surface hydrometeorological conditions on the Tibetan Plateau. - Temporal resolution: 1hour - Spatial resolution: point measurement - Measurement accuracy: soil moisture, 0.00001; soil temperature, 0.1 °C; data set size: soil moisture and temperature measurements at nominal depths of 5, 10, 20, 40 - Unit: soil moisture, cm ^ 3 cm ^ -3; soil temperature, °C

A multi-layer soil particle-size distribution dataset (sand, silt and clay content), based on USDA (United States Department of Agriculture) standard for regional land and climate modelling in China. was developed The 1:1,000,000 scale soil map of China and 8595 soil profiles from the Second National Soil Survey served as the starting point for this work. We reclassified the inconsistent soil profiles into the proper soil type of the map as much as possible because the soil classification names of the map units and profiles were not quite the same. The sand, silt and clay maps were derived using the polygon linkage method, which linked soil profiles and map polygons considering the distance between them, the sample sizes of the profiles, and soil classification information. For comparison, a soil type linkage was also generated by linking the map units and soil profiles with the same soil type. The quality of the derived soil fractions was reliable. Overall, the map polygon linkage offered better results than the soil type linkage or the Harmonized World Soil Database. The dataset, with a 1-km resolution, can be applied to land and climate modelling at a regional scale. Data characteristics： projection：projection Coverage: China Resolution: 0.00833 (about 1 km) Data format: FLT, TIFF Value range: 0%-100% Document describing： Floating point raster files include: Sand1. FLT, clay1. FLT -- surface (0-30cm) sand, clay content. Sand2. FLT, clay2. FLT -- content of sand and clay in the bottom layer (30-100cm). PSD. HDR -- header file: Ncols - the number of columns Nrows- rows Xllcorner - latitude in the lower left corner Yllcorner - longitude of the lower left corner Cellsize - cellsize NODATA_value - a null value byteorder - LSBFIRST, Least Significant Bit First TIFF raster files include: Sand1. Tif, clay1. Tif - surface (0-30cm) sand, clay content. Sand2. Tif, clay2. Tif - bottom layer (30-100cm) sand, clay content.

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.