The data set includes ASTER GDEM data and its Mosaic. ASTER Global DEM (ASTER GDEM) is a Global digital elevation data product jointly released by NASA and Japan's ministry of economy, trade and industry (METI) on June 29, 2009. The DEM data is based on the observation results of NASA's new earth observation satellite TERRA.It is produced by the ASTER(Advanced Space borne Thermal Emission and Reflection Radio meter) sensor, which collects 1.3 million stereo image data, covering more than 99% of the earth's land surface.The data has a horizontal accuracy of 30 m (95% confidence) and an elevation accuracy of 7-14 m (95% confidence).This data is the third global elevation data, which is significantly higher than previous SRTM3 DEM and GTOPO30 data. We from NASA's web site (http://wist.echo.nasa.gov/api) to download the data of heihe river basin, and through the data center to distribute.The data distributed by the center completely retains the original appearance of the data without any modification to the data.If users need details about ASTER GDEM preparation process, please refer to the data documents of metadata connections, or visit http://www.ersdac.or.jp/GDEM/E/3.html or directly from https://lpdaac.usgs.gov/ reading and ASTER Global DEM related documents. ASTER GDEM is divided into several data blocks of 1×1 degree in distribution, and the distribution format is zip compression format. Each compressed file includes three files. The file naming format is as follows: ASTGTM_NxxEyyy_dem.tif ASTGTM_NxxEyyy_num.tif reademe.pdf Where xx is the starting latitude and yyy is the starting longitude._dem. Tif is the dem data file, _num. Tif is the data quality file, and reademe is the data description file. In order to facilitate users to use the data, on the basis of the fractional ASTER GDEM data, we splice fractional SRTM data to prepare the ASTER GDEM Mosaic map of the black river basin, which retains all the original features of ASTER GDEM without any resamulation. This data includes two files: heihe_aster_gdem_mosaic_dem.img Heihe_Aster_GDEM_Mosaic_num. Img The data is stored in the format of Erdas image, where the file _dem.img is the dem data file and the file _num. Img is the data quality file.

The SRTM sensor has two bands, namely C-band and X-band. The SRTM we are using now comes from the C-band. The publicly released SRTM digital elevation products include DEM data at three different resolutions:     * SRTM1 covers only the continental United States, with a spatial resolution of 1s;     * SRTM3 data covers the world with a spatial resolution of 3s. This is the most widely used dataset. The elevation reference of SRTM3 is the geoid of EGM96 and the horizontal reference is WGS84. The nominal absolute elevation accuracy is ± 16m, and the absolute plane accuracy is ± 20m.     * SRTM30 data also covers the world, with a resolution of 30s. There are multiple versions of SRTM data. The early SRTM data was completed by NASA's "JPL" (Jet Propulsion Laboratory) ground data processing system (GDPS). The data is called SRTM3- 1. The National Geospatial Intelligence Agency has further processed the data, and the lack of data has been significantly improved. The data is called SRTM3-2. This dataset is mainly the fourth version of SRTM terrain data obtained by CIAT (International Center for Tropical Agriculture) using a new interpolation algorithm. This method better fills the SRTM 90 data hole. The interpolation algorithm comes from Reuter et al. (2007). The data of SRTM is organized as follows: every 5 latitude and longitude grids is divided into a file, which are divided into 24 rows (-60 to 60 degrees) and 72 columns (-180 to 180 degrees). The file naming rule is srtm_XX_YY.zip, where XX indicates the number of columns (01-72), and YY indicates the number of rows (01-24). The resolution of the data is 90 m. Data use: SRTM data uses a 16-bit value to represent the elevation value (-/ + / 32767 meters), the maximum positive elevation is 9000 meters, and the negative elevation (12,000 meters below sea level). -32767 standard for empty data.

This data set includes the daily averages of the temperature, pressure, relative humidity, wind speed, precipitation, global radiation, P2.5 concentration and other meteorological elements observed by the Qomolangma Station for Atmospheric and Environmental Observation and Research from 2005 to 2016. The data are aimed to provide service for students and researchers engaged in meteorological research on the Tibetan Plateau. The precipitation data are observed by artificial rainfall barrel, the evaporation data are observed by Φ20 mm evaporating pan, and all the others are daily averages and ten-day means obtained after half hour observational data are processed. All the data are observed and collected in strict accordance with the Equipment Operating Specifications, and some obvious error data are eliminated when processing the generated data.

This data set is based on the evaluation of existing land cover data and the evidence theory，including a 1:100,000 land use map for the year 20 2000、a 1:1,000,000 vegetation map、a 1:1,000,000 swamp-wetland map, a glacier map and a Moderate-Resolution Imaging Spectroradiometer land cover map for China in 2001 (MODIS2001) were merged，Finally, the decision is made based on the principle of maximum trust, and a new 1KM land cover data of China in 2000 with IGBP classification system is produced. The new land cover data not only maintain the overall accuracy of China's land use data, but also supplement the information of vegetation types and vegetation seasons in China's vegetation map, update China's wetland map, add the latest information of China's glacier map, and make the classification system more general.

There are many lakes in the Qinghai Tibet Plateau. The glacial phenology and duration of lakes in this region are very sensitive to regional and global climate change, so they are used as the key indicators of climate change research, especially the comparative study of the three polar environmental changes of the earth. However, due to its poor natural environment and sparse population, there is a lack of conventional field measurement of lake ice phenology. The lake ice was monitored with a resolution of 500 meters by using the normalized difference snow index (NDSI) data of MODIS. The traditional snow map algorithm is used to detect the lake daily ice amount and coverage under the condition of sunny days, and the lake daily ice amount and coverage under the condition of cloud cover are re determined through a series of steps based on the spatiotemporal continuity of the lake surface conditions. Through time series analysis, 308 lakes larger than 3km2 are identified as effective records of lake ice range and coverage, forming a daily lake ice range and coverage data set, including 216 lakes.

Near surface atmospheric forcing data were produced by using Wether Research and Forecasting (WRF) model over the Heihe River Basin at hourly 0.05 * 0.05 DEG resolution, including the following variables: 2m temperature, surface pressure, water vapor mixing ratio, downward shortwave & upward longwave radiation, 10m wind field and the accumulated precipitation. The forcing data were validated by observational data collected by 15 daily Chinese Meteorological Bureau conventional automatic weather station (CMA), a few of Heihe River eco-hydrological process comprehensive remote sensing observation (WATER and HiWATER) site hourly observations were verified in different time scales, draws the following conclusion: 2m surface temperature, surface pressure and 2m relative humidity are more reliable, especially 2m surface temperature and surface pressure, the average errors are very small and the correlation coefficients are above 0.96; correlation between downward shortwave radiation and WATER site observation data is more than 0.9; The precipitation agreed well with observational data by being verified based on rain and snow precipitation two phases at yearly, monthly, daily time scales . the correlation coefficient between rainfall and the observation data at monthly and yearly time scales were up to 0.94 and 0.84; the correlation between snowfall and observation data at monthly scale reached 0.78, the spatial distribution of snowfall agreed well with the snow fractional coverage rate of MODIS remote sensing product. Verification of liquid and solid precipitation shows that WRF model can be used for downscaling analysis in complex and arid terrain of Heihe River Basin, and the simulated data can meet the requirements of watershed scale hydrological modeling and water resources balance. The data for 2000-2012 was provided in 2013. The data for 2013-2015 was updated in 2016. The data for 2016-2018 was updated in 2019. The data for 2019-2021 was updated in 2021.

The Antarctic ice sheet is one of the largest potential sources of global sea level rise. Accurately determining the mass budget of the ice sheet is the key to understand the dynamic changes of the Antarctic ice sheet. It is very important to understand the evolution process of the ice sheet and accurately predict the future global sea level rise. Based on the MEaSUREs Antarctic groundingline and the basin boundaries, we discretize the groundingline, combine the MEaSUREs and RAMP annual ice velocity data from 1985 to 2015 with the BedMachine ice thickness data, and vectorially calculate the ice discharge at each flux gate of the groundingline. We use the surface mass balance data of RACMO2.3p2 model to spatially calculate the surface mass balance of each basin, and combined it with the ice discharge results to obtain the Antarctic ice sheet mass balance data set (1985-2015). The data set includes the mass balance results of each basin of the Antarctic ice sheet in the year 1985, 2000 and 2015, and the annual ice velocity data, ice thickness and annual ice discharge corresponding to the location of each flux gate. The data set realizes the fine evaluation of ice flux at the groundingline, and reflect the changes and spatial distribution characteristics of the mass balance of each basin of the Antarctic ice sheet in recent 30 years. It provides basic data for the subsequent fine change evaluation and prediction of the mass balance of the Antarctic ice sheet and the exploration of the mechanism of ice sheet loss.

The mass loss of the Greenland ice sheet has been the main contributor to global sea level rise in recent decades. Under the trend of global warming, the Greenland ice sheet is melting faster. It is of great scientific significance to explore the causes of mass loss and its response to climate change. Based on the MEaSUREs Greenland groundingline and the basin boundaries, we discretize the groundingline, combine the MEaSUREs annual ice velocity data from 1985 to 2015 with the BedMachine v3 ice thickness data, and vectorially calculate the ice discharge at each flux gate of the groundingline. We use the surface mass balance data of RACMO2.3p2 model to spatially calculate the surface mass balance of each basin, and combined it with the ice discharge results to obtain the Greenland ice sheet mass balance data set (1985-2015). The data set includes the mass balance results of each basin of the Greenland ice sheet in the year 1985, 2000 and 2015, and the annual ice velocity data, ice thickness and annual ice discharge corresponding to the location of each flux gate. The data set realizes the fine evaluation of ice flux at the groundingline, and reflect the changes and spatial distribution characteristics of the mass balance of each basin of the Greenland ice sheet in recent 30 years. It provides basic data for the subsequent fine change evaluation and prediction of the mass balance of the Greenland ice sheet and the exploration of the mechanism of ice sheet loss.

Based on AVHRR-CDR SR products, a daily cloud-free snow cover extent dataset with a spatial resolution of 5 km from 1981 to 2019 was prepared by using decision tree classification method. Each HDF4 file contains 18 data elements, including data value, data start date, longitude and latitude, etc. At the same time, to quickly preview the snow distribution, the daily file contains the snow area thumbnail, which is stored in JPG format. This data set will be continuously supplemented and improved according to the real-time satellite remote sensing data and algorithm update (up to may 2019), and will be fully open and shared.

The data of this study is mainly based on Google Earth Engine big data cloud processing platform. Sentinel-2 of The Three River Headwater region, Pul and Yukon River Basins in 2017 is selected as the basic data, STRM-DEM and Global Surface Water are used as auxiliary data. AWEIn，AWEIs，WI2015，MNDWI，NDWI and other index threshold extraction are selected to obtain seasonal water body and permanent water body according to annual water frequency(spatial resolution 10m). This water data product provides effective basic data for high spatial-temporal resolution water body change and permafrost hydrological analysis.

By applying Supply-demand Balance Analysis, the water resource supply and demand of the whole river basin and each county or district were calculated and used to evaluate the vulnerability of the water resources system of the basin. The IPAT equation was used to set a future water resource demand scenario to establish the scenario by setting variables such as future population growth rate, economic growth rate, and unit GDP water consumption. By taking 2005 as the base year and using assorted forecasting data of population size and economic scale, the future water demand scenarios of various counties and cities from 2010 to 2050 were predicted. By applying the basic structure of the HBV conceptual hydrological model of the Swedish Hydrometeorological Institute, a model of the variation tendency of the basin under climate change was designed. The glacial melting scenario was used as the model input to construct the runoff scenario under climate change. According to the national regulations of the water resources allocation of the basin, a water distribution plan was set up to calculate the water supply comprehensively. Considering of the supply and demand situation, the water resource system vulnerability was evaluated by the water shortage rate. By calculating the (grain production) land pressure index of the major counties and cities in the basin, the balance of supply and demand of land resources under the climate change, glacial melt and population growth scenarios was analyzed, and the vulnerability of the agricultural system was evaluated. The Miami formula and HANPP model were used to calculate the human appropriation of net primary biomass and primary biomass in the major counties and cities for the future, and the vulnerability of ecosystems from the perspective of supply and demand balance was assessed.

The Antarctic ice sheet elevation data were generated from radar altimeter data (Envisat RA-2) and lidar data (ICESat/GLAS). To improve the accuracy of the ICESat/GLAS data, five different quality control indicators were used to process the GLAS data, filtering out 8.36% unqualified data. These five quality control indicators were used to eliminate satellite location error, atmospheric forward scattering, saturation and cloud effects. At the same time, dry and wet tropospheric, correction, solid tide and extreme tide corrections were performed on the Envisat RA-2 data. For the two different elevation data, an elevation relative correction method based on the geometric intersection of Envisat RA-2 and GLAS data spot footprints was proposed, which was used to analyze the point pairs of GLAS footprints and Envisat RA-2 data center points, establish the correlation between the height difference of these intersection points (GLAS-RA-2) and the roughness of the terrain relief, and perform the relative correction of the Envisat RA-2 data to the point pairs with stable correlation. By analyzing the altimetry density in different areas of the Antarctic ice sheet, the final DEM resolution was determined to be 1000 meters. Considering the differences between the Prydz Bay and the inland regions of the Antarctic, the Antarctic ice sheet was divided into 16 sections. The best interpolation model and parameters were determined by semivariogram analysis, and the Antarctic ice sheet elevation data with a resolution of 1000 meters were generated by the Kriging interpolation method. The new Antarctic DEM was verified by two kinds of airborne lidar data and GPS data measured by multiple Antarctic expeditions of China. The results showed that the differences between the new DEM and the measured data ranged from 3.21 to 27.84 meters, and the error distribution was closely related to the slope.