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.

HiWATER: Dataset of hydrometeorological observation network (eddy covariance system of A’rou Superstation, 2013)

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.

HiWATER: Dataset of hydrometeorological observation network (eddy covariance system of A’rou Superstation, 2014)

This data set contains the eddy correlativity observation data from January 19, 2014 to December 31, 2014 of the super station at the upper reaches of heihe hydrometeorological observation network.The station is located in caoban village, aru township, qilian county, qinghai province.The longitude and latitude of the observation point are 100.4643e, 38.0473n and 3033m above sea level.The rack height of the vortex correlativity meter is 3.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked with red font, in which data is missing from January 1 to January 18 due to collector damage. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), stability Z/L (dimensionless), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al.(2018), and for observation data processing, please refer to Liu et al.(2011).

HiWATER: Dataset of hydrometeorological observation network (eddy covariance system of A’rou Superstation, 2015)

This data set contains the eddy correlograph observation data from January 1, 2015 to December 31, 2015 at the alou superstation, upstream of the heihe hydrometeorological observation network.The station is located in caoban village, aru township, qilian county, qinghai province.The longitude and latitude of the observation point are 100.4643e, 38.0473n and 3033m above sea level.The rack height of the vortex correlativity meter is 3.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked with red font, in which the calibration data of Li7500A of the eddy current system on April 16-17 is missing;When 10Hz data is missing due to a problem with the memory card storage data (9.20-10.21,11.3-11.18), the data is replaced by the 30-min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al. (2018), and for observation data processing, please refer to Liu et al. (2011).

HiWATER: Dataset of hydrometeorological observation network (eddy covariance system of A’rou Superstation, 2016)

This data set contains the observation data of vorticity correlativity from January 1, 2016 to December 31, 2016 of the super station at the upper reaches of heihe hydrometeorological observation network.The station is located in caoban village, aru township, qilian county, qinghai province.The longitude and latitude of the observation point are 100.4643e, 38.0473n and 3033m above sea level.The rack height of the vortex correlativity meter is 3.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked with red font, in which Li7500A calibration data of vorticity system from April 30 to May 1 is missing;When 10Hz data is missing due to a problem with the memory card storage data (11.9-11.24), the data will be replaced by 30min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al. (2018), and for observation data processing, please refer to Liu et al. (2011).

HiWATER: Dataset of hydrometeorological observation network (eddy covariance system of A’rou Superstation, 2017)

This data set contains the eddy correlativity observation data from January 1, 2017 to December 31, 2017 at the super station at the upper reaches of heihe hydrometeorological observation network.The station is located in caoban village, aru township, qilian county, qinghai province.The longitude and latitude of the observation point are 100.4643e, 38.0473n and 3033m above sea level.The rack height of the vortex correlativity meter is 3.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked with red font, in which the calibration data of the vortex system Li7500A from April 13 to April 14 is missing;When 10Hz data is missing due to a problem with the storage card (2.17-2.23, 3.3-4.12), the data will be replaced by the 30-min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al. (2018), and for observation data processing, please refer to Liu et al. (2011).