Data content: this data set is the historical archived satellite data of the domestic high score series (GF1 / 2 / 3 / 4) in the key river and lake research areas of the Qinghai Tibet Plateau from 2015 to 2020, which can cover the typical river and lake areas for effective monitoring. The time range of the data is from 2015 to 2020. Data source and processing method: the data are level 1 products. After equalizing radiation correction, the changes affecting the sensors are corrected by the equalizing functions of different detectors. Some data are based on the Landsat 8 images in the same period as the base map, and control points are selected for geometric correction of the images. Then, orthophoto correction is carried out based on DEM data, and band fusion processing is carried out for the corresponding data. Data quality description: the Gaofen series satellites are processed by the China Resources Satellite Application Center. There are raw data received by the satellite ground receiving station of the Chinese Academy of Sciences and processed products at all levels. Among them, level 1a (pre-processing level radiometric correction image product): image data processed by data analysis, uniform radiometric correction, noise removal, MTFC, CCD splicing, band registration, etc; And provide RPC files for satellite direct attitude orbit data production. Refer to the data website of China Resources Satellite Application Center for details. Data application achievements and prospects: the data are domestic high-resolution data with high resolution, which can be used to monitor the changes of the Qinghai Tibet Plateau as a water tower in Asia and the generated images, and test the accuracy of other data in the region

The dataset is the remote sensing image data ofGF-1 satellite in the Qinghai-Tibet engineering corridor obtained by China High Resolution Earth Observation Center. After the fusion processing of multispectral and panchromatic bands, the image data with a spatial resolution of 2 m is obtained. In the process of obtaining ground vegetation information, the classification technology of combining object-oriented computer automatic interpretation and manual interpretation is adopted, The object-oriented classification technology is to collect adjacent pixels as objects to identify the spectral elements of interest, make full use of high-resolution panchromatic and multispectral data space, texture and spectral information to segment and classify, and output high-precision classification results or vectors. In actual operation, the image is automatically extracted by eCognition software. The main processes are image segmentation, information extraction and accuracy evaluation. After verification with the field survey, the overall extraction accuracy is more than 90%.

The data set includes annual mass balance of Naimona’nyi glacier (northern branch) from 2008 to 2018, daily meteorological data at two automatic meteorological stations (AWSs) near the glacier from 2011 to 2018 and monthly air temperature and relative humidity on the glacier from 2018 to 2019. In the end of September or early October for each year , the stake heights and snow-pit features (snow layer density and stratigraphy) are manually measured to derive the annual point mass balance. Then the glacier-wide mass balance was then calculated （Please to see the reference). Two automatic weather stations (AWSs, Campbell company) were installed near the Naimona’nyi Glacier. AWS1, at 5543 m a. s.l., recorded meteorological variables from October 2011 at half hourly resolution, including air temperature (℃), relative humidity (%), and downward shortwave radiation (W m-2) . AWS2 was installed at 5950 m a.s.l. in October 2010 at hourly resolution and recorded wind speed (m/s), air pressure (hPa), precipitation (mm). Data quality: the quality of the original data is better, less missing. Firstly, the abnormal data in the original records are removed, and then the daily values of these parameters are calculated. Two probes (Hobo MX2301) which record air temperature and relative humidity was installed on the glacier at half hour resolution since October 2018. The observed meteorological data was calculated as monthly values. The data is stored in Excel file. It can be used by researchers for studying the changes in climate, hydrology, glaciers, etc.

(1) Data content: data set of precipitation field of the three poles (Arctic, Antarctic and Qinghai Tibet Plateau) in the past millennium; (2) Data source and processing method: the data is independently produced by the author and is produced by assimilating the precipitation proxy data in the three polar regions through the paleoclimate data assimilation method; (3) Data quality description: there is a high degree of spatial-temporal consistency between the data set and the precipitation data sets measured by multiple instruments (correlation coefficient is above 0.35, P < 0.001; Nash efficiency coefficient is above 0.3). In addition, the correlation coefficient with multiple precipitation data series reconstructed based on proxy data is between 0.2 and 0.6 (P < 0.001); (4) It can be used to study the temporal and spatial changes of precipitation in the past millennium in the three polar regions.

(1) Data content: Millennial temperature (near-surface air temperature anomaly based on the millinnial mean）datasets over the three poles, e.g., Arctic, Antarctic, and Qinghai-Tibet Plateau; (2) Data sources and processing methods: These datasets were produced by the authors themselves using the paleoclimate data assimilationand approach based on climatic proxies over the three poles; (3) Description of data quality: There are high spatio-temporal consistency between these datasets and several instrumental gridded temperature datasets (correlation coefficient above 0.6, p <0.001; Nash efficiency coefficient above 0.5). In addition, the correlations between these datasets and several proxy-based temperature series are between 0.4 and 0.8 (p <0.001). (4) Data application achievements and prospects: These datasets can be used to investigate the temporal and spatial variations in temperature over the three poles during the past millennium.

The basic data of hydrometeorology, land use and DEM were collected through the National Meteorological Information Center, the hydrological Yearbook, the China Statistical Yearbook and the Institute of geographical science and resources of the Chinese Academy of Sciences. The distributed time-varying gain hydrological model (DTVGM) with independent intellectual property rights is adopted for modeling, and the Qinghai Tibet Plateau is divided into 10937 sub basins with a threshold of 100 square kilometers. The daily flow data of 14 flow stations in Heihe River, Yarlung Zangbo River, Yangtze River source, Yellow River source, Yalong River, Minjiang River and Lancang River Basin were selected to draft and verify the model. The daily scale Naxi efficiency coefficient is above 0.7 and the correlation coefficient is above 0.8. The actual evaporation simulation is basically consistent with the station observation published by the Meteorological Bureau. The model simulates the water cycle process from 1998 to 2017. After verification, the spatial and temporal distribution of the actual evaporation (including soil evaporation and plant transpiration) on the 0.01 degree daily scale in the whole Tibetan Plateau is given.

The basic data of hydrometeorology, land use and DEM were collected through the National Meteorological Information Center, the hydrological Yearbook, the China Statistical Yearbook and the Institute of geographical science and resources of the Chinese Academy of Sciences. The distributed time-varying gain hydrological model (DTVGM) with independent intellectual property rights is adopted for modeling, and the Qinghai Tibet Plateau is divided into 10937 sub basins with a threshold of 100 square kilometers. The daily flow data of 14 flow stations in Heihe River, Yarlung Zangbo River, Yangtze River source, Yellow River source, Yalong River, Minjiang River and Lancang River Basin were selected to draft and verify the model. The daily scale Naxi efficiency coefficient is above 0.7 and the correlation coefficient is above 0.8. The model simulates the water cycle process from 1998 to 2017, and gives the spatial and temporal distribution of 0.01 degree daily scale runoff in the whole Qinghai Tibet Plateau.

The global high-resolution simulated near sea surface temperature precipitation SST data set from 1990 to 2020 is from the latest cmip6 project. Cmip6 is the sixth climate model comparison program organized by the world climate research project (WCRP). Original data source: https://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6 。 The data set includes the global near ocean surface temperature (TMP), precipitation (PR) and sea surface temperature (TOS). The air temperature and precipitation data include the rectangular combination of shared social economic path (SSP) and representative concentration path (RCP) of four different experimental scenarios of scenario MIP in cmip6. (1) Ssp126: upgrade rcp2.6 scenario based on ssp1 (low forcing scenario) (radiation forcing will reach 2.6w/m2 in 2100). (2) Ssp245: upgrade rcp4.5 scenario based on SSP2 (moderate forcing scenario) (radiation forcing will reach 4.5 w / m2 in 2100). (3) Ssp370: a new rcp7.0 emission path based on ssp3 (medium forcing scenario) (radiation forcing will reach 7.0 w / m2 in 2100). (4) Ssp585: upgrade rcp8.5 scenario based on ssp5 (high forcing scenario) (ssp585 is the only SSP scenario that can make radiation forcing reach 8.5 w / m2 in 2100). SST data provides ssp126 scenario data.

The three pole soil microbial post-treatment products in typical years collected the distribution and analysis results of soil samples from the north and south polar regions from 2005 to 2006 and the distribution and analysis results of soil samples from the Qinghai Tibet Plateau in 2015. Through sorting and summarizing, the post-processing data products of soil microorganisms in the three pole region are obtained. The data format is excel, which is convenient for users to view. Among them, the collection time of samples from the north and south polar regions was from December 13, 2005 to December 8, 2006, including 52 samples from three regions in the Arctic (Spitsbergen slijeringa, Spitsbergen vestpynten, and Alexandra fjord Highlands), and 171 samples from five regions in the Antarctic (Mitchell Peninsula, Casey station main power house, Robinson ridge, herring Island, browning Peninsula); The Qinghai Tibet Plateau was collected from July 1 to July 15, 2015, including meadow, grassland and desert ecosystems. There were 18 sampling points in total, and the number of samples at each sampling point was 3-5. The precipitation, air temperature and drought degree of the sampling point are estimated from the meteorological information for reference. The soil surface samples were collected and stored in liquid nitrogen, then transported back to the Sydney Laboratory for extraction by fastprep DNA kit. The extracted DNA samples were amplified with the 16S rRNA gene fragment using 27F (5'-gagttttgatcntggctca-3') and 519r (5'-gtnttacngcgckctg-3'). The amplified fragments were sequenced by 454 method, and the original data were analyzed by mothur software. The sequences with poor sequencing quality were first removed, and then the chimeric sequences were sequenced and removed. After that, the similarity between sequences is calculated. Sequences with a similarity of more than 97% are clustered into one OTU, and OTU representative sequences are defined. The OTU representative sequences were aligned with the Silva database, and were identified to belong to the first level when the reliability was greater than 80%.

We propose an algorithm for ice crack identification and detection using u-net network, which can realize the automatic detection of Antarctic ice cracks. Based on the data of sentinel-1 EW from January to February every year, in order to suppress the speckle noise of SAR image, the probabilistic patch based weights (ppb) algorithm is selected for filtering, and then representative samples are selected and input into the u-net network for model training, and the ice cracks are predicted according to the trained model. Taking five typical ice shelves（Amery、Fimbul、Nickerson、Shackleton、Thwaiters) in Antarctica as an example, the average accuracy of classification results can reach 94.5%, of which the local accuracy of fissure area can reach 78.6%, and the recall rate is 89.4%.

We propose an algorithm for ice fissure identification and detection using u-net network, which can realize the automatic detection of ice fissures of Typical Glaciers in Greenland ice sheet. Based on the data of sentinel-1 IW from July and August every year, in order to suppress the speckle noise of SAR image, the probabilistic patch based weights (ppb) algorithm is selected for filtering, and then the representative samples are selected and input into the u-net network for model training, and the ice cracks are predicted according to the trained model. Taking two typical glaciers in Greenland (Jakobshavn and Kangerdlussuaq) as examples, the average accuracy of classification results can reach 94.5%, of which the local accuracy of fissure area can reach 78.6%, and the recall rate is 89.4%.

This dataset provided the long-term elevation change rates data of the GrIS in three different periods using the ICESat data (February 2003 to October 2009), the Cryosat-2 data (August 2010 to October 2018) and the ICESat-2 data (October 2018 to December 2020) respectively. The dataset is named by the data.The data revealed that the elevation change rates of the GrIS were -12.19±3.81 cm/yr, -19.70±3.61 cm/yr and -23.39±3.06 cm/yr in the three different periods. In general, the obtained results agree with the trends discovered by other studies that were also derived from satellite altimetry data. This dataset provides the basic data for research into the impact of climate change over the GrIS.