The freezing / thawing state of near surface soil represents the dormancy and activity of land surface processes. This alternation of freezing and thawing phases can cause a series of complex surface process trajectory mode mutations, and affect the water cycle processes such as soil hydrothermal characteristics, surface runoff and groundwater recharge, and also affect climate change through water and energy cycle mechanism. This data set is based on AMSR-E and amsr2 passive microwave data, using discriminant algorithm to prepare global near earth surface freeze-thaw state (spatial resolution: 0.25 °; time span: 2002-2019), data storage type: 8-bit unsigned integer (file type:. HDF5) 5) Among them: 0: water body and missing data; 1: frozen soil; 2: thawed soil; 3: precipitation; 15: perennial snow and ice sheet. It can be used to analyze the spatial distribution and trend of the global freeze-thaw cycle, such as the start / end date, freezing / thawing duration, freezing range and other indicators. It can provide data support for understanding the interaction mechanism between land surface freeze-thaw cycle and water and energy exchange process under the background of global change. For detailed naming and missing of data, please refer to the data description.
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.
The “Long-term series of daily global snow depth” was produced using the passive microwave remote sensing data. The temporal range is 1979~2017, 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
Snow is a significant component of the ecosystem and water resources in high-mountain Asia (HMA). Therefore, accurate, continuous, and long-term snow monitoring is indispensable for the water resources management and economic development. The present study improves the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua satellites 8 d (“d” denotes “day”) composite snow cover Collection 6 (C6) products, named MOD10A2.006 (Terra) and MYD10A2.006 (Aqua), for HMA with a multistep approach. The primary purpose of this study was to reduce uncertainty in the Terra–Aqua MODIS snow cover products and generate a combined snow cover product. For reducing underestimation mainly caused by cloud cover, we used seasonal, temporal, and spatial filters. For reducing overestimation caused by MODIS sensors, we combined Terra and Aqua MODIS snow cover products, considering snow only if a pixel represents snow in both the products; otherwise it is classified as no snow, unlike some previous studies which consider snow if any of the Terra or Aqua product identifies snow. Our methodology generates a new product which removes a significant amount of uncertainty in Terra and Aqua MODIS 8 d composite C6 products comprising 46 % overestimation and 3.66 % underestimation, mainly caused by sensor limitations and cloud cover, respectively. The results were validated using Landsat 8 data, both for winter and summer at 20 well-distributed sites in the study area. Our validated adopted methodology improved accuracy by 10 % on average, compared to Landsat data. The final product covers the period from 2002 to 2018, comprising a combination of snow and glaciers created by merging Randolph Glacier Inventory version 6.0 (RGI 6.0) separated as debris-covered and debris-free with the final snow product MOYDGL06*. We have processed approximately 746 images of both Terra and Aqua MODIS snow containing approximately 100 000 satellite individual images. Furthermore, this product can serve as a valuable input dataset for hydrological and glaciological modelling to assess the melt contribution of snow-covered areas. The data, which can be used in various climatological and water-related studies, are available for end users at https://doi.org/10.1594/PANGAEA.901821 (Muhammad and Thapa, 2019).
The data sets include four sets of data obtained from the Scanning Multi-channel Microwave Radiometer (SMMR), Special Sensor Microwave Imager (SSM/I) and the Special Sensor Microwave Imager Sounder (SSMIS) sensors using passive microwave remote sensing inversion. SMMR was aboard the Nimbus-7 satellite, and its working period was from October 26, 1978 to July 8, 1987. Since July 1987, the data provided by the SSM/I and the SSMIS aboard the US Defense Meteorological Satellite Program (DMSP) satellite group have been used. The first three data sets contain sea ice concentration data, covering the Antarctic region with a spatial resolution of 25 km: (1) The data were obtained from Nimbus-7 SMMR and DMSP SSM/I-SSMIS Version 1 by applying the NASA Team algorithm inversion. The temporal coverage is from November 1978 to February 2017, with a temporal resolution of one month. A bin file is stored every month. (2) The data source is the same as the first set. The temporal coverage is from 1978-10-26 to 2017-2-28. The temporal resolution is two days, and the spatial resolution is 25 km. A folder was stored every year, and a bin file was stored every other day. (3) The data were obtained from near-real-time DMSP SSMIS by applying the NASA Team algorithm inversion. The temporal coverage is from 2015-1-1 to 2018-2-3, and the temporal resolution is one day. A bin file is stored every day. Each file consists of a 300-byte file title (data time information, projection pattern, file name) and a 316*332 matrix. The fourth set of data is the sea ice coverage and sea ice area time series. The temporal coverage is from November 1978 to December 2017. This data set is a time series sequence of sea ice coverage and sea ice area in the Antarctic. The temporal resolution is one month, and an ASCII file is stored every month. Each file consists of a file title (time, data type), a 39*1 sea ice cover matrix and a 39*1 sea ice area matrix. For further details on the data, please visit the US Ice and Snow Data Center NSIDC website - Data Description http://nsidc.org/data/NSIDC-0051; http://nsidc.org/data/NSIDC-0081; http://nsidc.org/data/G02135
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.
This dataset is the spatial distribution map of the marshes in the source area of the Yellow River near the Zaling Lake-Eling Lake, covering an area of about 21,000 square kilometers. The data set is classified by the Landsat 8 image through an expert decision tree and corrected by manual visual interpretation. The spatial resolution of the image is 30m, using the WGS 1984 UTM projected coordinate system, and the data format is grid format. The image is divided into five types of land, the land type 1 is “water body”, the land type 2 is “high-cover vegetation”, the land type 3 is “naked land”, and the land type 4 is “low-cover vegetation”, and the land type 5 is For "marsh", low-coverage vegetation and high-coverage vegetation are distinguished by vegetation coverage. The threshold is 0.1 to 0.4 for low-cover vegetation and 0.4 to 1 for high-cover vegetation.
This dataset contains the glacier outlines in Qilian Mountain Area in 2019. The dataset was produced based on classical band ratio criterion and manual editing. Chinese GF series images collected in 2019 were used as basic data for glacier extraction. Google images and Map World images were employed as reference data for manual adjusting. The dataset was stored in SHP format and attached with the attributions of coordinates, glacier ID and glacier area. Consisting of 1 season, the dataset has a spatial resolution of 2 meters. The accuracy is about 1 pixel (±2 meter). The dataset directly reflects the glacier distribution within the Qilian Mountain in 2018, and can be used for quantitative estimation of glacier mass balance and the quantitative assessment of glacier change’s impact on basin runoff.
The global Cryosat-2 GDR dataset is generated by the European Space Agency (ESA); it has a temporal coverage from 2010 to 2016 and covers the globe. On April 8, 2010, the ESA launched the Cryosat-2 high-tilt polar orbit satellite. The satellite is equipped with an SAR Interferometer Radar Altimeter (SIRAL), which is mainly used to monitor polar ice thickness and sea ice thickness changes, and, furthermore, to study the effects of melting polar ice on global sea level rise and that of global climate change on Antarctic ice thickness. The altimeter operates in the Ku-band and at a frequency of 13.575 GHz, it includes three measurement modes. One is a low-resolution altimeter measurement mode (LRM) that points to the subsatellite point to obtain all surface observations for land, sea, and ice sheets; its processing is similar to ENVISAT/RA-2, with an orbital resolution of 5 to 7 km. The second is the Synthetic Aperture Radar (SAR) measurement mode, which is mainly used to improve the accuracy and resolution of sea ice observations; it can make the resolution along the orbit reach approximately 250 m. The third is the Interferometric Synthetic Aperture Radar (InSAR), which is mainly used to improve the accuracy of areas with complex terrain such as the edges of ice sheets or ice shelves. The CryoSat -2/SIRAL data products mainly include 0-level data, 1b-level data, 2-level data and high-level data. The Cryosat-2/SIRAL products consist of two files: an XML head file (.HDR) and a data product file (.DBL). The HDR file is an auxiliary ASCII file for fast identification and retrieval of the data files. 1b-level products are stored separately according to the measurement modes, and the data recording formats of different modes are also different. Each waveform in LRM mode and SAR mode has 128 sampling points, while that in SARIn mode has 512 sampling points. 2-level GDR products are available for most scientific applications, including measurement time, geographic location, altitude, and more. In addition, the altitude information in GDR products has been obtained through instrumental calibration, transmission delay corrections, geometric corrections, and geophysical corrections (such as atmospheric corrections and tidal corrections). The GDR products are single global full-track data, that is, the measurement results of the three modes. After different processing, they are combined in chronological order; thereby, the data recording formats are unified. The data in the three modes use different waveform retracking algorithms to obtain altitude values. In the latest updated Baseline C data, the LRM mode data use three algorithms: Refined CFI, UCL and Refined OCOG.