This biophysical permafrost zonation map was produced using a rule-based GIS model that integrated a new permafrost extent, climate conditions, vegetation structure, soil and topographic conditions, as well as a yedoma map. Different from the previous maps, permafrost in this map is classified into five types: climate-driven, climate-driven/ecosystem-modified, climate-driven/ecosystem protected, ecosystem-driven, and ecosystem-protected. Excluding glaciers and lakes, the areas of these five types in the Northern Hemisphere are 3.66×106 km2, 8.06×106 km2, 0.62×106 km2, 5.79×106 km2, and 1.63×106 km2, respectively. 81% of the permafrost regions in the Northern Hemisphere are modified, driven, or protected by ecosystems, indicating the dominant role of ecosystems in permafrost stability in the Northern Hemisphere. Permafrost driven solely by climate occupies 19% of permafrost regions, mainly in High Arctic and high mountains areas, such as the Qinghai-Tibet Plateau.

DEM is the English abbreviation of Digital Elevation Model, which is the important original data of watershed topography and feature recognition.DEM is based on the principle that the watershed is divided into cells of m rows and n columns, the average elevation of each quadrilateral is calculated, and then the elevation is stored in a two-dimensional matrix.Since DEM data can reflect local topographic features with a certain resolution, a large amount of surface morphology information can be extracted through DEM, which includes slope, slope direction and relationship between cells of watershed grid cells, etc..At the same time, the surface flow path, river network and watershed boundary can be determined according to certain algorithm.Therefore, to extract watershed features from DEM, a good watershed structure pattern is the premise and key of the design algorithm. Elevation data map 1km data formed according to 1:250,000 contour lines and elevation points in China, including DEM, hillshade, Slope and Aspect maps. Data set projection: Two projection methods: Equal Area projection Albers Conical Equal Area (105, 25, 47) Geodetic coordinates WGS84 coordinate system

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 distribution of lakes in space and its change over time are closely related to agricultural, environmental and ecological issues, and are critical factors for human socio-economic development. In the past decades, satellite based remote sensing has been developed rapidly to provide essential data sources for monitoring temporal lakes dynamics with its advantage of rapidness, wide coverage, and lower cost. This dataset was produced from Landsat images using the automated water detection method (Feng et al, 2015). We collected 96,278 Landsat images (about 25 terabytes) that acquired since 2000 with less than 80% cloud contamination in the arid region of central Asia and Tibetan Plateau. Water is detected in each of the image and then aggregated to monthly temporal resolution by taking advantage of the high-performance processing capability and large data storage provided by Global Land Cover Facility (GLCF) at University of Maryland. The results are validated systematically and quantitatively using manually interpreted dataset, which consists of a set of locations collected by a stratified random sampling strategy to effectively represent different spatial-temporal distributions in the region. The validation suggests high accuracy of the results (overall accuracy: 99.45(±0.59); user accuracy: 85.37%±(3.74); produce accuracy: 98.17(±1.05)).

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

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.

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 dataset includes the monthly precipitation data with 0.0083333 arc degree (~1km) for China from Jan 1901 to Dec 2020. The data form belongs to NETCDF, namely .nc file. The unit of the data is 0.1 mm. The dataset was spatially downscaled from CRU TS v4.02 with WorldClim datasets based on Delta downscaling method. The dataset was evaluated by 496 national weather stations across China, and the evaluation indicated that the downscaled dataset is reliable for the investigations related to climate change across China. The dataset covers the main land area of China, including Hong Kong, Macao and Taiwan regions, and excluding islands and reefs in South China Sea.

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.