The Tibetan Plateau Glacier Data –TPG2013 is a glacial coverage data on the Tibetan Plateau around 2013. 128 Landsat 8 Operational Land Imager (OLI) images were selected with 30-m spatial resolution, for comparability with previous and current glacier inventories. Besides, about 20 images acquired in 2014 were used to complete the full coverage of the TP. The most frequent year in this period was defined as the reference year for the mosaic image: i.e. 2013. Glacier outlines were digitized on-screen manually from the 2013 image mosaic, relying on false-colour image composites (RGB by bands 654), which allowed us to distinguish ice/snow from cloud. Debris-free ice was distinguished from the debris and debris-covered ice by its higher reflectance. Debris-covered ice was not delineated in this data. [To minimize the effects of snow or cloud cover on glacierized areas, high-resolution (30 m spatial resolution and 4-day repetition cycle) images were also used for reference in glacier delineation from the Chinese satellites HJ-1A and HJ-1B, which were launched on Sep.6th 2008. Both carried as payload two 4-band CCD cameras with swath width 700 km (360 km per camera). All HJ-1A/1B data in 2012, 2013 and 2014 (65 scenes, Fig.S1, Table S1) were from China Centre for Resources Satellite Data and Application (CRESDA; http://www.cresda.com/n16/n92006/n92066/n98627/index.html). Each scene was orthorectified with respect to the 30m-resolution digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) and Landsat images.] The delineated glacier outlines were compared with band-ratio (e.g. TM3/TM5) results, and validated by overlapping them onto Google Earth imagery, SRTM DEM, topographic maps and corresponding satellite images. Topographic maps from the 1970s and all available satellite images (including Google EarthTM imagery and HJ-1A/1B satellite data) were used as base reference data. For areas with mountain shadows and snow cover, they were verified by different methods using data from different seasons. For glaciers in deep shadow, Google EarthTM imagery from different dates was used as the reference for manual delineation. Steep slopes or headwalls were also excluded in the TPG2013. Areas that appeared in any of these sources to have the characteristics of exposed ground/basement/bed rock were manually delineated as non-glacier, and were also cross-checked with CGI-1 and CGI-2. Steep hanging glaciers were included in TPG2013 if they were identifiable on images in all three epochs (i.e. TPG1976, TPG2001, and TPG2013). The accuracy of manual digitization was controlled within one half-pixel. All glacier areas were calculated on the WGS84 spheroid in an Albers equal-area map projection centred at (95°E, 30°N) with standard parallels at 15°N and 65°N. Our results showed that the relative deviation of manual interpretation was less than 3.9%.

Geladandong region is an important and typical source region of great rivers and lakes in the Qinghai Tibet Plateau. This data set provides DEM covering glaciers in the source region of the Yangtze River and Selin Co with different time scales and resolutions to calculate the seasonal and decadal changes of glacier surface elevation in the source region. This data set includes seven 5-meter resolution TanDEM-X data from July 2016 to 2017, which can be used to calculate the seasonal change of glacier surface elevation; it includes one KH-9 DEM with a resolution of 30m in 1976, five TanDEM-X with a resolution of 30m in 2011, one TanDEM-X in 2014 and three TanDEM-X in 2017 with a resolution of 30m. The data can be used to calculate the change of glacier surface elevation during 1976-2000, 2000-20112011-2017. At the same time, Landsat ETM data are used to extract the glacier outline in 1976and we divide it according to the RGI6.0; The right figure shows the spatial and temporal coverage information of the data set, and the base figure is the orthophoto corrected kh-9 image.

Snow over sea ice controls the energy budgets, affects the sea ice growth/melting, and thus has essential climatic effects. Snow depth, one of the fundamental properties of snow cover, is essential for understanding of the rapid change in Antarctic climate and for sea ice thickness estimation. Passive microwave radiometer can be used for basin-scale snow depth estimation in daily scale, however, previous published methods applied for Antarctic snow depth shows clear underestimation, which limits their further application. Here, we construct a new and robust linear regression equation for snow depth retrieval using microwave radiometers by including lower frequencies, and we produce the snow depth product over Antarctic sea ice from 2002 to 2020 from AMSR-E, AMSR-2, SSMIS based on this method. A regression analysis using 7 years of Operation IceBridge (OIB) airborne snow depth measurements shows that the gradient ratio (GR) calculated using brightness temperatures in vertical polarized 37 and 19 GHz, i.e., GR(37/7), is the optimal one for deriving Antarctic snow depth with an root mean square deviation (RMSD) of 8.92 cm and a correlation coefficient of -0.64, the related equation coefficients are then derived. GR(37/19) is used to retrieve snow depth from SSMIS data to fill the observation gaps between AMSR-E and AMSR-2, and the estimated snow depth is corrected for the consistence with these from AMSR-E/2. An averaged uncertainty of 3.81 cm is found based on a Gaussian error propagation, which accounts for 12% of the estimated mean snow depth. The evaluation of proposed method with in-situ measurements from Australian Antarctic Data Centre shows that the proposed method outperforms the previous available method, with a mean difference of 5.64 cm and an RMSD of 13.79 cm, comparing to -14.47 cm and 19.49 cm. Comparison to shipborne observations from Antarctic Sea Ice Processes and Climate indicates that the proposed method shows slight better performance than previous method (RMSDs of 16.85 cm and 17.61 cm, respectively); and comparable performances in growth and melting seasons suggests that the proposed method can still be used in the melting season. We generate a complete snow depth product over Antarctic sea ice from 2002 to 2020 in daily scale, and negative trends can be found in all sea sectors and seasons. This dataset can be further used in the reanalysis data evaluation, sea ice thickness estimation, climate model and other aspects.

The “Long-term series of daily global snow depth” was produced using the passive microwave remote sensing data. The temporal range is 1980~2018, and the coverage is the global land. The spatial resolutions is 25,067.53 m and the temporal resolution is daily. A dynamic brightness temperature gradient algorithm was used to derive snow depth. In this algorithm, the spatial and temporal variations of snow characteristics were considered and the spatial and seasonal dynamic relationships between the temperature difference between 18 GHz and 36 GHz and the measured snow depth were established. The long-term sequence of satellite-borne passive microwave brightness temperature data used to derive snow depth came from three sensors (SMMR, SSM/I and SSMI/S), and there is a certain system inconsistency among them. So, the inter-sensor calibration was performed to improve the temporal consistency of these brightness temperature data before snow depth derivation. The accuracy analysis shows that the relative deviation of Eurasia snow depth data is within 30%. The data are stored as a txt file every day, each file is a 1383*586 snow depth matrix, and each snow depth represents a 25,067.53m* 25,067.53m grid. The projection of this data is EASE-Grid, and following is the file header which describes the projection detail. File header: ncols 1383 nrows 586 xllcorner -17334193.54 yllcorner -7344787.75 cellsize 25,067.53 NODATA_value -1

The multi-decadal lake number and area changes in China during 1960s–2020 are derived from historical topographic maps and >42151 Landsat satellite images, including lakes as fine as ≥1 km^2 in size for the past 60 years (1960s, 1970s, 1990, 1995, 2000, 2005, 2010, 2015, 2020). From the 1960s to 2020, the total number of lakes (≥ 1 km ^ 2) in China increased from 2127 to 2621, and the area expanded from 68537 km ^ 2 to 82302 km ^ 2.

Ec-earth-heihe USES the output of the global model of ec-earth as the driving field to simulate the 6-hour data of the Heihe river basin in 2006-2080 under the scenarios of 1980-2005 and RCP4.5.Spatial scope: the grid center of the simulation area is located at (40.30n, 99.50e), the horizontal resolution is 3 km, and the number of simulated grid points in the model is 161 (meridional) X 201 (zonal). Projection: LAMBERT conformal projection, two standard latitudes of 30N and 60N. Time range: from January 1, 1980 to December 31, 2010, with an interval of 6 hours. Description of file contents: monthly storage by grads without format.Except the maximum and minimum temperature as the daily scale, the other variables are all 6-hour data. MATLAB can be used to read, visible tmax_erain_xiong_heihe.m file description. Data description of heihe river basin: 1) Anemometer west wind (m/s) abbreviation usurf 2) Anemometer south wind(m/s), abbreviation vsurf 3) Anemometer temperature (deg K) abbreviation tsurf 4) maximal temperature (deg K) abbreviation tmax 5) minimal temperature (deg K) abbreviated tmin 6) Anemom specific humidity (g/kg) abbreviation qsurf 7) Accumulated precipitation (mm/hr) abbreviation precip 8) Accumulated evaporation (mm/hr) abbreviation evap 9) Accumulated sensible heat (watts/m**2/hr) abbreviation sensible 10) Accumulated net infrared radiation (watts/m * * 2 / hr) abbreviation netrad File name definition: Abbreviation-ec-earth-6hour，YTD For example, precip-ec-earth-6hour.198001，Is the data of 6-hour precipitation in January, 1980 (1) historical 6-hour data driven by the ec-earth global climate model from 1980 to 2005 (2) produce 6-hour data of heihe river basin under the scenario of RCP 4.5 for the global climate model ec-earth from 2006 to 2080

The content is the daily runoff observation record of the outlet weir of the Pailugou basin. The spatial range of Pailugou: 38.529-38.558N, 100.286-100.536E. Data dates include May 1, 2013 to September 5, 2013. The unit is m3/day.

Soil data is important both on a global scale and on a local scale, and due to the lack of reliable soil data, land degradation assessments, environmental impact studies, and sustainable land management interventions have received significant bottlenecks . Affected by the urgent need for soil information data around the world, especially in the context of the Climate Change Convention, the International Institute for Applied Systems Analysis (IIASA) and the Food and Agriculture Organization of the United Nations (FAO) and the Kyoto Protocol for Soil Carbon Measurement and FAO/International The Global Agroecological Assessment Study (GAEZ v3.0) jointly established the Harmonized World Soil Database version 1.2 (HWSD V1.2). Among them, the data source in China is the second national land in 1995. Investigate 1:1,000,000 soil data provided by Nanjing Soil. The resolution is 30 seconds (about 0.083 degrees, 1km). The soil classification system used is mainly FAO-90. The core soil system unit unique verification identifier: MU_GLOBAL-HWSD database soil mapping unit identifier, connected to the GIS layer. MU_SOURCE1 and MU_SOURCE2 source database drawing unit identifiers SEQ-soil unit sequence in the composition of the soil mapping unit; The soil classification system utilizes the FAO-7 classification system or the FAO-90 classification system (SU_SYM74 resp. SU_SYM90) or FAO-85 (SU_SYM85). The main fields of the soil property sheet include: ID (database ID) MU_GLOBAL (Soil Unit Identifier) ​​(Global) SU_SYMBOL soil drawing unit SU_SYM74 (FAO74 classification); SU_SYM85 (FAO85 classification); SU_SYM90 (name of soil in the FAO90 soil classification system); SU_CODE soil charting unit code SU_CODE74 soil unit name SU_CODE85 soil unit name SU_CODE90 soil unit name DRAINAGE (19.5); REF_DEPTH (soil reference depth); AWC_CLASS(19.5); AWC_CLASS (effective soil water content); PHASE1: Real (soil phase); PHASE2: String (soil phase); ROOTS: String (depth classification to the bottom of the soil); SWR: String (soil moisture content); ADD_PROP: Real (specific soil type in the soil unit related to agricultural use); T_TEXTURE (top soil texture); T_GRAVEL: Real (top gravel volume percentage); (unit: %vol.) T_SAND: Real (top sand content); (unit: % wt.) T_SILT: Real (surface layer sand content); (unit: % wt.) T_CLAY: Real (top clay content); (unit: % wt.) T_USDA_TEX: Real (top layer USDA soil texture classification); (unit: name) T_REF_BULK: Real (top soil bulk density); (unit: kg/dm3.) T_OC: Real (top organic carbon content); (unit: % weight) T_PH_H2O: Real (top pH) (unit: -log(H+)) T_CEC_CLAY: Real (cation exchange capacity of the top adhesive layer soil); (unit: cmol/kg) T_CEC_SOIL: Real (cation exchange capacity of top soil) (unit: cmol/kg) T_BS: Real (top level basic saturation); (unit: %) T_TEB: Real (top exchangeable base); (unit: cmol/kg) T_CACO3: Real (top carbonate or lime content) (unit: % weight) T_CASO4: Real (top sulfate content); (unit: % weight) T_ESP: Real (top exchangeable sodium salt); (unit: %) T_ECE: Real (top conductivity). (Unit: dS/m) S_GRAVEL: Real (bottom crushed stone volume percentage); (unit: %vol.) S_SAND: Real (bottom sand content); (unit: % wt.) S_SILT: Real (bottom sludge content); (unit: % wt.) S_CLAY: Real (bottom clay content); (unit: % wt.) S_USDA_TEX: Real (bottom USDA soil texture classification); (unit: name) S_REF_BULK: Real (bottom soil bulk density); (unit: kg/dm3.) S_OC: Real (underlying organic carbon content); (unit: % weight) S_PH_H2O: Real (bottom pH) (unit: -log(H+)) S_CEC_CLAY: Real (cation exchange capacity of the underlying adhesive layer soil); (unit: cmol/kg) S_CEC_SOIL: Real (cation exchange capacity of the bottom soil) (unit: cmol/kg) S_BS: Real (underlying basic saturation); (unit: %) S_TEB: Real (underlying exchangeable base); (unit: cmol/kg) S_CACO3: Real (bottom carbonate or lime content) (unit: % weight) S_CASO4: Real (bottom sulfate content); (unit: % weight) S_ESP: Real (underlying exchangeable sodium salt); (unit: %) S_ECE: Real (underlying conductivity). (Unit: dS/m) The database is divided into two layers, with the top layer (T) soil thickness (0-30 cm) and the bottom layer (S) soil thickness (30-100 cm). For other attribute values, please refer to the HWSD1.2_documentation documentation.pdf, The Harmonized World Soil Database (HWSD V1.2) Viewer-Chinese description and HWSD.mdb.

This dataset is land surface phenology estimated from 16 days composite MODIS NDVI product (MOD13Q1 collection6) in the Three-River-Source National Park from 2001 to 2020. The spatial resolution is 250m. The variables include Start of Season (SOS) and End of Season (EOS). Two phenology estimating methods were used to MOD13Q1, polynomial fitting based threshold method and double logistic function based inflection method. There are 4 folders in the dataset. CJYYQ_phen is data folder for source region of the Yangtze River in the national park. HHYYQ_phen is data folder for source region of Yellow River in the national park. LCJYYQ_phen is data folder for source region of Lancang River in the national park. SJY_phen is data folder for the whole Three-River-Source region. Data format is geotif. Arcmap or Python+GDAL are recommended to open and process the data.

This data is derived from the National Basic Geographic Information Resources Catalogue Service System, which provides 11 million national basic geographic databases free of charge by the National Basic Geographic Information Center in November 2017. We have spliced and cut the source of the three rivers as a whole, so as to facilitate the use of the study of the source area of the three rivers. This data set is composed of 1:1 million administrative boundary layers (BOUA) and administrative boundary line layers (BOUL) in Sanjiangyuan area. Names and definitions of BOUA attribute items: Attribute Item Description Fill in Example PAC Administrative Division Code 513230 NAME Name Rangtang County Names and definitions of BOUL attribute items: Attribute Item Description Fill in Example GB National Standard Classification Code 630200 The meaning of BOUL attribute items: Attribute Item Code Description GB 630200 Provincial Boundary GB 640200 District, Municipal and State Administrative Region GB 650201 County administrative boundaries (determined)

This dataset includes one scene acquired on (yy-mm-dd) 2012-07-25, covering the natural oasis eco-hydrology experimental area in the lower reaches of the Heihe River Basin. This datum contains panchromatic and multi-spectral bands, with spatial resolution of 0.6 m and 2.4 m, respectively. The data product level of this image is Level 2A. QuickBird dataset was acquired through purchase.

On 19 August 2012, a Leica ALS70 airborne laser scanner boarded by the Y-12 aircraft was used to obtain the point cloud data. Leica ALS70 airborne laser scanner has unlimited numbers of returns intensities measurements including the first, second, third return intensities. The wavelength of laser light is 1064 nm. The absolute flight altitude is 2900 m with the point cloud density 1 point per square meter. Aerial LiDAR-DSM was obtained through parameter calibration, automatic classification of point cloud density and manual editing.