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 dataset uses daily temperature data from SMMR (1978-1987), SSM/I (1987-2009) and SSMIS (2009-2015). It is generated by the dual-index (TB, 37v, SG) freeze-thaw discrimination algorithm. The classification results include the frozen surface, the thawed surface, the deserts and water bodies. The data coverage is the main part of China’s mainland, with a spatial resolution of 25.067525 km via the EASE-Grid projection method, and it is stored in ASCIIGRID format. All the ASCII files in this data set can be opened directly with a text program such as Notepad. Except for the head file, the body content is numerically characterized by the freeze/thaw status of the surface soil: 1 for frozen, 2 for thawed, 3 for desert, and 4 for precipitation. If you want to use the icon for display, we recommend using the ArcView + 3D or Spatial Analyst extension module for reading; in the process of reading, a grid format file will be generated, and the displayed grid file is the graphical expression of the ASCII file. The read method comprises the following. [1] Add the 3D or Spatial Analyst extension module to the ArcView software and then create a new View. [2] Activate View, click File menu, and select the Import Data Source option. When the Import Data Source selection box pops up, select ASCII Raster in the Select import file type box. When the dialog box for selecting the source ASCII file automatically pops up, click to find any ASCII file in the data set, and then press OK. [3] Type the name of the Grid file in the Output Grid dialog box (it is recommended that a meaningful file name is used for later viewing) and click the path to store the Grid file, press OK again, and then press Yes (to select integer data) and Yes (to put the generated grid file into the current view). The generated files can be edited according to the Grid file standard. This completes the process of displaying an ASCII file into a Grid file. [4] In the batch processing, the ASCIGRID command of ARCINFO can be used to write AML files, and then use the Run command to complete the process in the Grid module: Usage: ASCIIGRID <in_ascii_file> <out_grid> {INT | FLOAT}. The production of this data is supported by the following Natural Science Foundation Projects: Environmental and Ecological Science Data Center of West China (90502010), Land Data Assimilation System of West China (90202014) and Active and Passive Microwave Radiation Transmission Simulation and Radiation Scattering Characteristics of the Frozen Soil (41071226).

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.

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.

The GAME/Tibet project conducted a short-term pre-intensive observing period (PIOP) at the Amdo station in the summer of 1997. From May to September 1998, five consecutive IOPs were scheduled, with approximately one month per IOP. More than 80 scientific workers from China, Japan and South Korea went to the Tibetan Plateau in batches and carried out arduous and fruitful work. The observation tests and plans were successfully completed. After the completion of the IOP in September, 1998, five automatic weather stations (AWS), one Portable Atmospheric Mosonet (PAM), one boundary layer tower and integrated radiation observatory (Amdo) and nine soil temperature and moisture observation stations have been continuously observed to date and have obtained extremely valuable information for 8 years and 6 months consecutively (starting from June 1997). The experimental area is located in Nagqu, in northern Tibet, and has an area of 150 km × 200 km (Fig. 1), and observation points are also established in D66, Tuotuohe and the Tanggula Mountain Pass (D105) along the Qinghai-Tibet Highway. The following observation stations (sites) are set up on different underlying surfaces including plateau meadows, plateau lakes, and desert steppe. (1) Two multidisciplinary (atmosphere and soil) observation stations, Amdo and NaquFx, have multicomponent radiation observation systems, gradient observation towers, turbulent flux direct measurement systems, soil temperature and moisture gradient observations, radiosonde, ground soil moisture observation networks and multiangle spectrometer observations used as ground truth values for satellite data, etc. (2) There are six automatic weather stations (D66, Tuotuohe, D105, D110, Nagqu and MS3608), each of which has observations of wind, temperature, humidity, pressure, radiation, surface temperature, soil temperature and moisture, precipitation, etc. (3) PAM stations (Portable Automated Meso - net) located approximately 80 km north and south of Nagqu (MS3478 and MS3637) have major projects similar to the two integrated observation stations (Amdo and NaquFx) above and to the wind, temperature and humidity turbulence observations. (4) There are nine soil temperature and moisture observation sites (D66, Tuotuohe, D110, WADD, NODA, Amdo, MS3478, MS3478 and MS3637), each of which has soil temperature measurements of 6 layers and soil moisture measurement of 9 layers. (5) A 3D Doppler Radar Station is located in the south of Nagqu, and there are seven encrypted precipitation gauges in the adjacent (within approximately 100 km) area. The radiation observation system mainly studies the plateau cloud and precipitation system and serves as a ground true value station for the TRMM satellite. The GAME-Tibet project seeks to gain insight into the land-atmosphere interaction on the Tibetan Plateau and its impact on the Asian monsoon system through enhanced observational experiments and long-term monitoring at different spatial scales. After the end of 2000, the GAME/Tibet project joined the “Coordinated Enhanced Observing Period (CEOP)” jointly organized by two international plans, GEWEX (Global Energy and Water Cycle Experiment) and CL IVAR (Climate Change and Forecast). The Asia-Australia Monsoon Project (CAMP) on the Tibetan Plateau of the Global Coordinated Enhanced Observation Program (CEOP) has been started. The data set contains POP data for 1997 and IOP data for 1998. Ⅰ. The POP data of 1997 contain the following. 1. Precipitation Gauge Network (PGN) 2. Radiosonde Observation at Naqu 3. Analysis of Stable Isotope for Water Cycle Studies 4. Doppler radar observation 5. Large-Scale Hydrological Cycle in Tibet (Link to Numaguchi's home page) 6. Portable Automated Mesonet (PAM) [Japanese] 7. Ground Truth Data Collection (GTDC) for Satellite Remote Sensing 8. Tanggula AWS (D105 station in Tibet) 9. Syamboche AWS (GEN/GAME AWS in Nepal) Ⅱ. The IOP data of 1998 contain the following. 1. Anduo （1） PBL Tower, 2） Radiation, 3） Turbulence SMTMS 2. D66 （1） AWS （2） SMTMS （3） GTDC （4) Precipitation 3. Toutouhe （1） AWS （2） SMTMS （3 ）GTDC 4. D110 （1） AWS （2） SMTMS (3) GTDC (4) SMTMS 5. MS3608 （1） AWS （2） SMTMS （3） Precipitation 6. D105 （1） Precipitation (2) GTDC 7. MS3478(NPAM) （1） PAM （2） Precipitation 8. MS3637 （1） PAM （2） SMTMS （3） Precipitation 9. NODAA （1） SMTMS (2) Precipitation 10. WADD （1） SMTMS （2） Precipitation （3） Barometricmd 11. AQB （1） Precipitation 12. Dienpa (RS2) （1） Precipitation 13. Zuri （1） Precipitation （2） Barometricmd 14. Juze （1） Precipitation 15. Naqu hydrological station （1） Precipitation 16. MSofNaqu （1） Barometricmd 16. Naquradarsite （1）Radar system （2） Precipitation 17. Syangboche [Nepal] （1） AWS 18. Shiqu-anhe （1） AWS （2） GTDC 19. Seqin-Xiang （1） Barometricmd 20. NODA （1）Barometricmd （2） Precipitation （3) SMTMS 21. NaquHY （1） Barometricmd （2） Precipitation 22. NaquFx(BJ) （1） GTDC（2) PBLmd (3) Precipitation 23. MS3543 （1） Precipitation 24. MNofAmdo （1） Barometricmd 25. Mardi （1） Runoff 26. Gaize （1） AWS （2） GTDC （3） Sonde A CD of the data GAME-Tibet POP/IOP dataset cd （vol. 1) GAME-Tibet POP/IOP dataset cd （vol. 2)

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.

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.

This data set is collected from the supplementary information part of the paper: Yao, T. , Thompson, L. , & Yang, W. . (2012). Different glacier status with atmospheric circulations in tibetan plateau and surroundings. Nature Climate Change, 1580, 1-5. This paper report on the glacier status over the past 30 years by investigating the glacial retreat of 82 glaciers, area reductionof 7,090 glaciers and mass-balance change of 15 glaciers. This data set contains 8 tables, the names and content are as follows: Data list: The data name list of the rest tables; t1: Distribution of Glaciers in the TP and surroundings; t2: Data and method for analyzing glacial area reduction in each basin; t3: Glacial area reduction during the past three decades from remote sensing images in the TP and surroundings; t4: Glacial length fluctuationin the TP and surroundings in the past three decades; t5: Detailed information on the glaciers for recent mass balance measurement in the TP and surroundings; t6: Recent annual mass balances in different regions in the TP; t7: Mass balance of Long-time series for the Qiyi, Xiaodongkemadi and Kangwure Glaciers in the TP. See attachments for data details: Supplementary information.pdf, Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings.pdf.

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