According to the global soil map. Net standard, the 0-1m soil depth is divided into 5 layers: 0-5cm, 5-15cm, 15-30cm, 30-60cm and 60-100cm. According to the principle of soil landscape model, the spatial distribution data products of soil silt content in different layers are made by using the digital soil mapping method. The American system classification is used as the standard of soil particle classification. The source data of this data set comes from the soil profile data integrated by the major research plan integration project of Heihe River Basin (soil data integration and soil information product generation of Heihe River Basin, 91325301). Scope: Heihe River Basin; Projection: WGS · 1984 · Albers; Spatial resolution: 100M; Data format: TIFF; Dataset content: hh_silt_layer1.tif: 0-5cm soil silt content; hh_silt_layer2.tif: 5-15cm soil silt content; hh_silt_layer3.tif: 15-30cm soil silt content; hh_silt_layer4.tif: 30-60cm soil silt content; hh_silt_layer5.tif：60-100cm soil silt content;

According to the global soil map. Net standard, the 0-1m soil depth is divided into 5 layers: 0-5cm, 5-15cm, 15-30cm, 30-60cm and 60-100cm. According to the principle of soil landscape model, the spatial distribution data products of soil clay content in different layers are made by using the digital soil mapping method. The American system classification is used as the standard of soil particle classification. The source data of this data set comes from the soil profile data integrated by the major research plan integration project of Heihe River Basin (soil data integration and soil information product generation of Heihe River Basin, 91325301). Scope: Heihe River Basin; Projection: WGS · 1984 · Albers; Spatial resolution: 100M; Data format: TIFF; Dataset content: hh_clay_layer1.tif: 0-5cm soil clay content; hh_clay_layer2.tif: 5-15cm soil clay content; hh_clay_layer3.tif: 15-30cm soil clay content; hh_clay_layer4.tif: 30-60cm soil clay content; hh_clay_layer5.tif: 60-100cm soil clay content;

This dataset includes five scenes, covering the artificial oasis eco-hydrology experimental area of the Heihe River Basin, which were acquired on (yy-mm-dd) 2012-04-05, 2012-04-21, 2012-05-07, 2012-06-24, 2012-07-10. The data were all acquired around 11:50 (BJT) with data product of Level 2. Landsat ETM+ dataset was downloaded from http://glovis.usgs.gov/.

This dataset includes one scene acquired on (yy-mm-dd hh:mm, BJT) 2012-07-06 06:30, covering the artificial oasis eco-hydrology experimental area of the Heihe River Basin. This datum was acquired at Stripmap-Quad mode with product level of SLC, and this image includes VV, VH, HH and HV polarization with a spatial resolution of 8 m. Radarsat-2 dataset was acquired from the Institute of Remote Sensing and Digital Earth of Chinese Academy of Sciences (Courtesy: Dr. Chen Quan).

The aerosol optical thickness data of the Arctic Alaska station is based on the observation data products of the atmospheric radiation observation plan of the U.S. Department of energy at the Arctic Alaska station. The data coverage time is updated from 2017 to 2019, with the time resolution of hour by hour. The coverage site is the northern Alaska station, with the longitude and latitude coordinates of (71 ° 19 ′ 22.8 ″ n, 156 ° 36 ′ 32.4 ″ w). The source of the observed data is retrieved from the radiation data observed by mfrsr instrument. The characteristic variable is aerosol optical thickness, and the error range of the observed inversion is about 15%. The data format is NC format. The aerosol optical thickness data of Qomolangma station and Namuco station in the Qinghai Tibet Plateau is based on the observation data products of Qomolangma station and Namuco station from the atmospheric radiation view of the Institute of Qinghai Tibet Plateau of the Chinese Academy of Sciences. The data coverage time is from 2017 to 2019, the time resolution is hour by hour, the coverage sites are Qomolangma station and Namuco station, the longitude and latitude coordinates are (Qomolangma station: 28.365n, 86.948e, Namuco station Mucuo station: 30.7725n, 90.9626e). The source of the observed data is retrieved from the radiation data observed by mfrsr instrument. The characteristic variable is aerosol optical thickness, and the error range of the observed inversion is about 15%. The data format is TXT.

The data set includes soil bulk density data of representative soil samples collected from July 2012 to August 2013 in the Heihe River Basin. The first soil survey was conducted in 2012. After the representativeness evaluation of collected samples, we conducted an additional sampling in 2013. These samples are representative enough to represent the soil variation in the Heihe River Basin, of which the soil variation in each landscape could be accounted for. The sampling depths in field refer to the sampling specification of Chinese Soil Taxonomy, in which soil samples were taken from genetic soil horizons.

The data set includes soil organic carbon concentrations data of representative soil samples collected from July 2012 to August 2013 in the Heihe River Basin. The first soil survey was conducted in 2012. After the representativeness evaluation of collected samples, we conducted an additional sampling in 2013. These samples are representative enough to represent the soil variation in the Heihe River Basin, of which the soil variation in each landscape could be accounted for. The sampling depths in field refer to the sampling specification of Chinese Soil Taxonomy, in which soil samples were taken from genetic soil horizons.

This data is obtained by spatial interpolation and permafrost simulation through the surface temperature at 0 cm of nine stations in and outside the source area of the upper reaches of Heihe River. In the figure, 1 represents seasonal frozen soil and 2 represents permafrost. The data is in TIFF format, WGS-84 is used for projection, and the spatial range is 37.7263n-39.0976n, 98.5769e-101.1608e.

The dataset contains all individual glacial storage (unit: km3) over the Qinghai-Tibetan Plateau in 1970s and 2000s. It is sourced from the resultant data of the paper entitled "Consolidating the Randolph Glacier Inventory and the Glacier Inventory of China over the Qinghai-Tibetan Plateau and Investigating Glacier Changes Since the mid-20th Century". The first draft of this paper has been completed and is planned to be submitted to Earth System Science Data journal. The baseline glacier inventories in 1970s and 2000s are the Randolph Glacier Inventory 4.0 dataset, and the Glacier Inventory of China, respectively. Based on the individual glacial boundaries extracted from the above-mentioned two datasets, the grid-based bedrock elevation dataset (https://www.ngdc.noaa.gov/mgg/global/global.html, DOI: 10.7289/v5c8276m), and the glacier surface elevation obtained by a slope-dependent method, the individual glacier volumes in 1970s and 2000s are then calculated. In addition, the calculated results of individual glacier volumes in this study have been compared and verified with the existent results of several glacier volumes, relevant remote sensing datasets, and the global glacier thickness dataset based on the average of multiple glacier model outputs (https://www.research-collection.ethz.ch/handle/20.500.11850/315707, doi:10.3929/ethz-b-000315707), and the errors in the calculations have also been quantified. The established dataset in this study is expected to provide the data basis for the future regional water resources estimation and glacier ablation-involved researches. Moreover, the acquisition of the data also provides a new idea for the future glacier storage estimation.

The data set contains all single glacial reserves (in KM3) in the Tibetan Plateau of 1970s and 2000s. This data set comes from the result data of the paper entitled "consolidating the Randolph glacier inventory and the glacier inventory of China over the Qinghai titanium plate and investigating glacier changes since the mid-20th century". The first draft of this paper has been completed and is planned to be submitted to earth system science data. The 1970s basic glacier catalog data in the dataset is extracted from Randolph glacier Inventory data set, 2000s basic glacial catalogue is from China's second glacial catalogue data set. Based on the glacial boundary extracted from the two data sets and combined with the grid based bedrock elevation data set (https://www.ngdc.noaa.gov/mgg/global/global.html, DOI: 10.7289/v5c8276m) and the glacial table obtained by a slope dependent method Based on the surface elevation data set, the single glacier reserves in the two catalogues are calculated. In addition, the calculation results of single glacier reserves obtained in this study have been compared and verified with the calculation results of partial glacier reserves, relevant remote sensing data sets, and the global glacier thickness data set based on the average of multiple glacier model sets in multiple directions, and the errors in the calculation results have also been quantified. The establishment of the data set is expected to provide the data basis for the future regional water resources estimation and glacier ablation research, and the acquisition of the data also provides a new idea for the future glacier reserves research.

"Coupling and Evolution of Hydrological-Ecological-Economic Processes in Heihe River Basin Governance under the Framework of Water Rights" (91125018) Project Data Convergence-The documents of the west of Taolai River water conservancy team project plan 1. Data summary: The documents of the west of Taolai River water conservancy team project plan 2. Data content: Taolai River water conservancy team project plan, including the project plan of reservoir irrigation and drainage in the west of the river region

This dataset contains the flux measurements from the A’rou Superstation eddy covariance system (EC) in the upper reaches of the Heihe hydrometeorological observation network from 29 December, 2012, to 19 November, 2013. The site (100.464° E, 38.047° N) was located in the Alpine grassland, Caodaban village, A’rou Country, Qilian County in Qilian Province. The elevation is 3033 m. The EC was installed at a height of 3.5 m, and the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (CSAT3 & Li7500A) was 0.15 m. The raw data acquired at 10 Hz were processed using the Edire post-processing software (University of Edinburgh, http://www.geos.ed.ac.uk/abs/research/micromet/EdiRe/), including the spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. The observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC), as proposed by Foken and Wichura [1996]: class 1 (level 0: Δst<30 and ITC<30), class 2 (level 1: Δst<100 and ITC<100), and class 3 (level 2: Δst>100 and ITC>100), which represent high-, medium-, and low-quality data, respectively. In addition to the above processing steps, the half-hourly flux data were screened using a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data collected before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 3% of the 30 min raw record; and (4) data were rejected at night when the friction velocity (u*) was less than 0.1 m/s. There were 48 records per day, and the missing data were replaced with -6999. Suspicious data were marked in red. The data from 14 February to 28 February were missing due to a storage problem of CF card in the datalogger. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m^3), CO2 mass density (CO2, mg/m^3), friction velocity (ustar, m/s), stability (z/L), sensible heat flux (Hs, W/m2), latent heat flux (LE, W/m^2), carbon dioxide flux (Fc, mg/ (m^2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xls format. For more information, please refer to Liu et al. (2018) (for hydrometeorological observation network or sites information), Liu et al. (2011) (for data processing) in the Citation section.