The dataset includes soil physical and chemical attributes: pH value, organic matter fraction, cation exchange capacity, root abundance, total nitrogen (N), total phosphorus (P), total potassium (K), alkali-hydrolysable N, available P, available K, exchangeable H+, Al3+, Ca2+, Mg2+, K+ , Na+, horizon thickness, soil profile depth, sand, silt and clay fractions, rock fragment, bulk density, porosity, structure, consistency and soil color. Quality control information (QC) was provided. The resolution is 30 arc-seconds (about 1 km at the equator). The vertical variation of soil property was captured by eight layers to the depth of 2.3 m (i.e. 0- 0.045, 0.045- 0.091, 0.091- 0.166, 0.166- 0.289, 0.289- 0.493, 0.493- 0.829, 0.829- 1.383 and 1.383- 2.296 m) for convenience of use in the Common Land Model and the Community Land Model (CLM). 1.THSCH.nc: Saturated water content of FCH 2.PSI_S.nc: Saturated capillary potential of FCH 3.LAMBDA.nc: Pore size distribution index of FCH 4.K_SCH.nc: Saturate hydraulic conductivity of FCH 5.THR.nc: Residual moisture content of FGM 6.THSGM.nc: Saturated water content of FGM 7.ALPHA.nc: The inverse of the air-entry value of FGM 8.N.nc: The shape parameter of FGM 9.L.nc: The pore-connectivity parameter of FGM 10.K_SVG.nc: Saturated hydraulic conductivity of FGM 11.TH33.nc: Water content at -33 kPa of suction pressure, or field capacity 12.TH1500.nc: Water content at -1500 kPa of suction pressure, or permanent wilting point

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 dataset is the monthly precipitation data of China, with a spatial resolution of 0.0083333 ° (about 1km) and a time range of 1901.1-2021.12. The data format is NETCDF, i.e.. Nc format. This dataset is generated in China through the Delta spatial downscaling scheme based on the global 0.5 ° climate dataset released by CRU and the global high-resolution climate dataset released by WorldClim. In addition, 496 independent meteorological observation point data are used for verification, and the verification results are reliable. This data set covers the main land areas in China (including Hong Kong, Macao and Taiwan), excluding islands and reefs in the South China Sea. In order to facilitate storage, the data are all int16 type and stored in nc files, with precipitation units of 0.1mm. NC data can be mapped using ArcMAP software; Matlab software can also be used for extraction processing. Matlab has released the function to read and store nc files. The read function is ncread, and switch to the nc file storage folder. The statement is expressed as: ncread ('XXX.nc ',' var ', [i j t], [leni lenj lent]), where XXX.nc is the file name, and is the string required' '; Var is from XXX The variable name read in NC. If it is a string, '' is required; i. J and t are the starting row, column and time of the read data respectively, and leni, lenj and lent i are the length of the read data in the row, column and time dimensions respectively. In this way, this function can be used to read in any region and any time period in the study area. There are many commands about NC data in the help of Matlab, which can be viewed. WGS84 is recommended for data coordinate system.

As the main parameter in the land surface energy balance, surface temperature indicates the degree of land-atmosphere energy and water transfer and is widely used in research on climatology, hydrology and ecology. In the study of frozen soil, climate is one of the decisive factors for the existence and development of frozen soil. The surface temperature is the main climatic factor affecting the distribution of frozen soil and affects the occurrence, development and distribution of frozen soil. It is the upper boundary condition for modelling frozen soil and is significant to the study of hydrological processes in cold regions. The data set was based on the DEM and observation station data of the Tibetan Plateau Engineering Corridor and analysed the changing trend of surface temperature on the Tibetan Plateau from 2000 to 2014. Using the surface temperature data products MOD11A1/A2 and MYD11A1/A2 of MODIS aboard Terra and Aqua, the surface temperature information under cloud cover was reconstructed based on the spatio-temporal information of the images. The reconstruction information and surface temperature representativeness problems were analysed using information obtained from 8 sites, including the Kunlun Mountains (wetland, grassland), Beiluhe (grassland, meadow), Kaixinling (meadow, grassland), and Tanggula Mountain (meadow, wetland). According to the correlation coefficient (R2), root-mean-square error (RMSE), mean absolute error (MAE) and mean deviation (MBE), the following results were obtained: (1) the reconstruction accuracy of MODIS surface temperature under cloud cover is higher when it is based on spatio-temporal information; (2) the weighted average representation is the best when generalizing four observations of Terra and Aqua. By analysing the reconstruction of MODIS surface temperature information and representativeness problems, the average annual MODIS surface temperature data of the Tibetan Plateau and the engineering corridor from 2000 to 2010 were obtained. According to the data set, the surface temperature from 2000 to 2010 also experienced volatile rising trends from 2000 to 2010, which is basically consistent with the changing trend of the climate change in the permafrost regions of the Tibetan Plateau and the Qinghai-Tibet Engineering Corridor.

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.

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.

This data set is based on the evaluation of existing land cover data and the evidence theory，including a 1:100,000 land use map for the year 20 2000、a 1:1,000,000 vegetation map、a 1:1,000,000 swamp-wetland map, a glacier map and a Moderate-Resolution Imaging Spectroradiometer land cover map for China in 2001 (MODIS2001) were merged，Finally, the decision is made based on the principle of maximum trust, and a new 1KM land cover data of China in 2000 with IGBP classification system is produced. The new land cover data not only maintain the overall accuracy of China's land use data, but also supplement the information of vegetation types and vegetation seasons in China's vegetation map, update China's wetland map, add the latest information of China's glacier map, and make the classification system more general.

Based on the analysis of brgdgts and hydrogen isotopes of leaf wax in lake sediments from Tengchong Qinghai (tcqh) in Yunnan Province, this study shows for the first time the high-resolution annual average temperature change history of low latitude land since the last glacial period (since the last 88000 years). According to the annual average temperature of South Asia established by tcqh core, there are two warm periods of 88000-71000 years and 45000-22000 years in this region, and the temperature range is about 2-3 ° C. Since the Holocene, the temperature has been increasing for about 1-2 years ° C。

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