These datasets include mean annual ground temperature (MAGT) at the depth of zero annual amplitude (approximately 3 m to 25 m), active layer thickness (ALT), the probability of the permafrost occurrence, and the new permafrost zonation based on hydrothermal condition for the period of 2000-2016 in the Northern Hemisphere with an 1-km resolution by integrate unprecedentedly large amounts of field data (1,002 boreholes for MAGT and 452 sites for ALT) and multisource geospatial data, especially remote sensing data, using statistical learning modelling with an ensemble strategy, and thus more accurate than previous circumpolar maps.

Effective evaluation of future climate change, especially prediction of future precipitation, is an important basis for formulating adaptation strategies. This data is based on the RegCM4.6 model, which is compatible with multi-model and different carbon emission scenarios: CanEMS2 (RCP 45 and RCP85), GFDL-ESM2M (RCP2.6, RCP4.5, RCP6.0 and RCP8.5), HadGEM2-ES (RCP2.6, RCP4.5 And RCP8.5), IPSL-CM5A-LR (RCP2.6, RCP4.5, RCP6.0 and RCP8.5), MIROC5 (RCP2.6, RCP4.5, RCP6.0 and RCP8.5). The future climate data (2007-2099) has 21 sets, with a spatial resolution at 0.25 degrees and the temporal resolution at 3 hours (or 6 hours), daily and yearly scales.

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.

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.

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

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

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 field observation platform of the Tibetan Plateau is the forefront of scientific observation and research on the Tibetan Plateau. The land surface processes and environmental changes based comprehensive observation of the land-boundary layer in the Tibetan Plateau provides valuable data for the study of the mechanism of the land-atmosphere interaction on the Tibetan Plateau and its effects. This dataset integrates the 2005-2016 hourly atmospheric, soil hydrothermal and turbulent fluxes observations of Qomolangma Atmospheric and Environmental Observation and Research Station, Chinese Academy of Sciences (QOMS/CAS), Southeast Tibet Observation and Research Station for the Alpine Environment, CAS (SETORS), the BJ site of Nagqu Station of Plateau Climate and Environment, CAS (NPCE-BJ), Nam Co Monitoring and Research Station for Multisphere Interactions, CAS (NAMORS), Ngari Desert Observation and Research Station, CAS (NADORS), Muztagh Ata Westerly Observation and Research Station, CAS (MAWORS). It contains gradient observation data composed of multi-layer wind speed and direction, temperature, humidity, air pressure and precipitation data, four-component radiation data, multi-layer soil temperature and humidity and soil heat flux data, and turbulence data composed of sensible heat flux, latent heat flux and carbon dioxide flux. These data can be widely used in the analysis of the characteristics of meteorological elements on the Tibetan Plaetau, the evaluation of remote sensing products and development of the remote sensing retrieval algorithms, and the evaluation and development of numerical models.

China's second glacier inventory uses the high-resolution Landsat TM/ETM+ remote sensing satellite data as the main glacier boundary data source and extracts the data source with the latest global digital elevation model, SRTM V4, as the glacier attribute, using the current international ratio threshold segmentation method to extract the glacier boundary in bare ice areas. The ice ridge extraction algorithm is developed to extract the glacier ice ridge, and it is used for the segmentation of a single glacier. At the same time, the international general algorithm is used to calculate the glacier attributes, so that the vector data and attribute data that contain the glacier information of the main glacier regions in west China are obtained. Compared with some field GPS field measurement data and higher resolution remote sensing images (such as from QuickBird and WorldView), the glacial vector data in the second glacier inventory data set of China have higher positioning accuracy and can meet the requirements for glacial data in national land, water conservancy, transportation, environment and other fields. Glacier inventory attributes: Glc_Name, Drng_Code, FCGI_ID, GLIMS_ID, Mtn_Name, Pref_Name, Glc_Long, Glc_Lati, Glc_Area, Abs_Accu, Rel_Accu, Deb_Area, Deb_A_Accu, Deb_R_Accu, Glc_Vol_A, Glc_Vol_B, Max_Elev, Min_Elev, Mean_Elev, MA_Elev, Mean_Slp, Mean_Asp, Prm_Image, Aux_Image, Rep_Date, Elev_Src, Elev_Date, Compiler, Verifier. For a detailed data description, please refer to the second glacier inventory data description.