The Antarctic McMurdo Dry Valleys ice velocity product is based on the Antarctic Ice Sheet Velocity and Mapping Project (AIV) data product, which is post-processed with advanced algorithms and numerical tools. The product is mapped using Sentinel-1/2/Landsat data and provides uniform, high-resolution (60m) ice velocity results for McMurdo Dry Valleys, covering the period from 2015 to 2020.
JIANG Liming JIANG Liming JIANG Liming
Based on the data of GF-1 and GF-2 in China, the freeze-thaw disaster distribution data of Qinghai Tibet project corridor is produced by using the deep learning classification method and manual visual interpretation and correction. The geographical range of the data is 40km along the Xidatan Anduo section of Qinghai Tibet highway. The data include the distribution data of thermokast lakes and the distribution data of thermal melting landslides. The dataset can provide data basis for the research of freeze-thaw disaster and engineering disaster prevention and reduction in Qinghai Tibet engineering corridor. The spatial distribution of freezing and thawing disasters within 40km along the Xidatan-Anduo section of Qinghai Tibet highway is self-made based on the domestic GF-2 image data. Firstly, the deep learning method is used to extract the mud flow terrace block from GF-2 data; Then, ArcGIS is used for manual editing.
NIU Fujun, LUO Jing LUO Jing
Pine Island Glacier, Swett Glacier, etc. are distributed in the basins of the Antarctic Ice Sheet 21 and 22, which is one of the areas with the most severe melting in the Southwest Antarctica. This dataset first uses Cryosat-2 data (August 2010 to October 2018) to establish a plane equation in each regular grid, taking into account terrain items, seasonal fluctuations, backscattering coefficients, wave front width, lifting rails and other factors, and calculates the elevation change of ice cover surface in the grid through least square regression. In addition, we used ICESat-2 data (October 2018 to December 2020) to calculate the surface elevation change during the two periods by obtaining the elevation difference at the intersection of satellite lifting orbits in each regular grid. The spatial resolution of surface elevation change data in two periods is 5km × 5km, the file format is GeoTIFF, the projection coordinate is polar stereo projection (EPSG 3031), and it is named by the name of the satellite altimetry data used. The data can be opened using ArcMap, QGIS and other software. The results show that the average elevation change rate of the region from 2010 to 2018 is -0.34 ± 0.08m/yr, which belongs to the area with severe melting. The annual average elevation change rate from October 2018 to November 2020 is -0.38 ± 0.06m/yr, which is in an intensified state compared with CryoSat-2 calculation results.
YANG Bojin , HUANG Huabing , LIANG Shuang , LI Xinwu
This data set includes five periods of lake transparency data, including 1995, 2002, 2005, 2010 and 2015. The data sources are: Landsat 5, Landsat 7 and Landsat 8. Method of use: It is convenient to measure the spectral reflectance. On the basis of analyzing the relationship between the spectral reflectance and the transparency measured synchronously, the semi empirical method is used to select the best band combination, establish the transparency algorithm of Qinghai Tibet Plateau lakes, and obtain the water transparency. The verification of measured points shows that the relative error of water transparency estimation is 35%.
SONG Kaishan
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%.
LI Xinwu , LIANG Shuang , YANG Bojin , ZHAO Jingjing
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%.
LI Xinwu , LIANG Shuang , YANG Bojin , ZHAO Jingjing
Global solar radiation and diffuse horizontal solar radiation at Dome C (Antarctica) are measured by radiation sensors (pyranometers CM22, Kipp & Zonen Inc., The Netherlands), and water vapor pressure (hPa) at the ground are obtained from the IPEV/PNRA Project “Routine Meteorological Observation at Station Concordia”, http://www.climantartide.it. This dataset includes hourly solar radiation and its absorbing and scattering losses caused by the absorbing and scattering atmospheric substances (MJ m-2, 200-3600 nm), and the albedos at the top of the atmosphere and the surface. The above solar radiations are calculated by using an empirical model of global solar radiation (Bai, J.; Zong, X.; Lanconelli, C.; Lupi, A.; Driemel, A.; Vitale, V.; Li, K.; Song, T. 2022. Long-Term Variations of Global Solar Radiation and Its Potential Effects at Dome C (Antarctica). Int. J. Environ. Res. Public Health, 19, 3084. https://doi.org/10.3390/ijerph19053084). The observed global solar radiation and meteorological parameters are available at https://doi.org/10.1594/PANGAEA.935421. The data set can be used to study solar radiation and its attenuation at Dome C, Antarctica.
BAI Jianhui
Global solar radiation at Qomolangma station (The Tibetan Plateau) is measured by radiation sensor (pyranometers CM22, Kipp & Zonen Inc., The Netherlands), and water vapor pressure (hPa) at the ground is measured by HMP45C-GM (Vaisala Inc., Vantaa, Finland). This dataset includes hourly solar radiation and its absorbing and scattering losses caused by the absorbing and scattering atmospheric substances (MJ m-2, 200-3600 nm), and the albedos at the top of the atmosphere and the surface. The above solar radiations are calculated by using an empirical model of global solar radiation (Bai, J.; Zong, X.; Ma, Y.; Wang, B.; Zhao, C.; Yang, Y.; Guang, J.; Cong, Z.; Li, K.; Song, T. 2022. Long-Term Variations in Global Solar Radiation and Its Interaction with Atmospheric Substances at Qomolangma. Int. J. Environ. Res. Public Health, 19, 8906. https://doi.org/10.3390/ijerph19158906). The observed global solar radiation and meteorological variables are available at https://data.tpdc.ac.cn/zh-hans/data/b9ab35b2-81fb-4330-925f-4d9860ac47c3/. The data set can be used to study solar radiation and its attenuation at Qomolangma region.
BAI Jianhui
The fractional snow cover (FSC) is the ratio of snow cover area (SCA) to unit pixel area. The data set is made by bv-blrm snow area proportional linear regression empirical model; The source data used are mod09ga 500m global daily surface reflectance products and mod09a1 500m 8-day synthetic global surface reflectance products; The production platform uses Google Earth engine; The data range is global, the data preparation time is from 2000 to 2021, the spatial resolution is 500 meters, and the temporal resolution is year by year. This set of data can provide quantitative information of snow cover distribution for regional climate simulation and hydrological models.
MA Yuan
Fractional Vegetation Cover (FVC) refers to the percentage of the vertical projected area of vegetation to the total area of the study area. It is an important indicator to measure the effectiveness of ecological protection and ecological restoration. It is widely used in the fields of climate, ecology, soil erosion and so on. FVC is not only an ideal parameter to reflect the productivity of vegetation, but also can play a good role in evaluating topographic differences, climate change and regional ecological environment quality. This research work is mainly to post process two sets of glass FVC data, and give a more reliable vegetation coverage of the circumpolar Arctic Circle (north of 66 ° n) and the Qinghai Tibet Plateau (north of 26 ° n to 39.85 °, east longitude 73.45 ° to 104.65 °) in 2013 and 2018 through data fusion, elimination of outliers and clipping.
YE Aizhong
NDVI reflects the background effects of plant canopy, such as soil, wet ground, snow, dead leaves, roughness, etc., and is related to vegetation cover. It is one of the important parameters to reflect the crop growth and nutrient information. According to this parameter, the N demand of crops in different seasons can be known, which is an important guide to the reasonable application of N fertilizer. Correct NDVI (C-NDVI) is the value of NDVI after excluding the influence of climate elements (temperature, precipitation, etc.) on NDVI. Taking precipitation as an example, studies on the lag effect of precipitation on vegetation growth show that the lag time of precipitation effects varies in different regions due to differences in vegetation composition and soil types. In this study, we post-processed the MODIS NDVI data and firstly correlated the NDVI value of the current month with the precipitation of the current month, the average value of the precipitation of the current month with that of the previous month, and the average value of the precipitation of the current month with that of the previous two months to determine the optimal lag time. The NDVI was regressed on precipitation and air temperature to obtain the correlation coefficients, and then the corrected NDVI values were calculated by the difference between the MODIS NDVI and the NDVI regressed on climate factors. We corrected NDVI using climate data to give reliable vegetation correction indices for the circum-Arctic Circle (range north of 66°N) and the Tibetan Plateau (range 26°N to 39.85°N and 73.45°E to 104.65°E) for 2013 and 2018. The spatial resolution of the data is 0.5 degrees and the temporal resolution is monthly values.
YE Aizhong
This dataset consists of four files including (1) Lake ice thickness of 16 large lakes measured by satellite altimeters for 1992-2019 (Altimetric LIT for 16 large lakes.xlsx); (2) Daily lake ice thickness and lake surface snow depth of 1,313 lakes with an area > 50 km2 in the Northern Hemisphere modeled by a one-dimensional remote sensing lake ice model for 2003-2018 (in NetCDF format); (3) Future lake ice thickness and surface snow depth for 2071-2099 modeled by the lake ice model with a modified ice growth module (table S1.xlsx); (4) A lookup table containing lake IDs, names, locations, and areas. This daily lake ice and snow thickness dataset could provide a benchmark for the estimation of global lake ice and snow mass, thereby improving our understanding of the ecological and economical significance of freshwater ice as well as its response to climate change.
LI Xingdong, LONG Di, HUANG Qi, ZHAO Fanyu
Precipitation over the Tibetan Plateau (TP) known as Asia's water tower plays a critical role in regional water and energy cycles, largely affecting water availability for downstream countries. Rain gauges are indispensable in precipitation measurement, but are quite limited in the TP that features complex terrain and the harsh environment. Satellite and reanalysis precipitation products can provide complementary information for ground-based measurements, particularly over large poorly gauged areas. Here we optimally merged gauge, satellite, and reanalysis data by determining weights of various data sources using artificial neural networks (ANNs) and environmental variables including elevation, surface pressure, and wind speed. A Multi-Source Precipitation (MSP) data set was generated at a daily timescale and a spatial resolution of 0.1° across the TP for the 1998‒2017 period. The correlation coefficient (CC) of daily precipitation between the MSP and gauge observations was highest (0.74) and the root mean squared error was the second lowest compared with four other satellite products, indicating the quality of the MSP and the effectiveness of the data merging approach. We further evaluated the hydrological utility of different precipitation products using a distributed hydrological model for the poorly gauged headwaters of the Yangtze and Yellow rivers in the TP. The MSP achieved the best Nash-Sutcliffe efficiency coefficient (over 0.8) and CC (over 0.9) for daily streamflow simulations during 2004‒2014. In addition, the MSP performed best over the ungauged western TP based on multiple collocation evaluation. The merging method could be applicable to other data-scarce regions globally to provide high quality precipitation data for hydrological research. The latitude and longitude of the left bottom corner across the TP, the number of rows and columns, and grid cells information are all included in each ASCII file.
HONG Zhongkun , LONG Di
Aiming at the 179000 km2 area of the pan three rivers parallel flow area of the Qinghai Tibet Plateau, InSAR deformation observation is carried out through three kinds of SAR data: sentinel-1 lifting orbit and palsar-1 lifting orbit. According to the obtained InSAR deformation image, it is comprehensively interpreted in combination with geomorphic and optical image features. A total of 949 active landslides below 4000m above sea level were identified. It should be noted that due to the difference of observation angle, sensitivity and observation phase of different SAR data, there are some differences in the interpretation of the same landslide with different data. The scope and boundary of the landslide need to be corrected with the help of ground and optical images. The concept of landslide InSAR recognition scale is different from the traditional spatial resolution and mainly depends on the deformation intensity. Therefore, some landslides with small scale but prominent deformation characteristics and strong integrity compared with the background can also be interpreted (with SAR intensity map, topographic shadow map and optical remote sensing image as ground object reference). The minimum interpretation area can reach several pixels. For example, a highway slope landslide with only 4 pixels is interpreted with reference to the highway along the Nujiang River.
YAO Xin
This data is a high-resolution soil freeze/thaw (F/T) dataset in the Qinghai Tibet Engineering Corridor (QTEC) produced by fusing sentinel-1 SAR data, AMSR-2 microwave radiometer data, and MODIS LST products. Based on the newly proposed algorithm, this product provides the detection results of soil F/T state with a spatial resolution of 100 m on a monthly scale. Both meteorological stations and soil temperature stations were used for results evaluation. Based on the ground surface temperature data of four meteorological stations provided by the national meteorological network, the overall accuracy of soil F/T detection products achieved 84.63% and 77.09% for ascending and descending orbits, respectively. Based on the in-situ measured 5 cm soil temperature data near Naqu, the average overall accuracy of ascending and descending orbits are 78.58% and 76.66%. This high spatial resolution F/T product makes up traditional coarse resolution soil F/T products and provides the possibility of high-resolution soil F/T monitoring in the QTEC.
ZHOU Xin , LIU Xiuguo , ZHOU Junxiong , ZHANG Zhengjia , CHEN Qihao , XIE Qinghua
Soil moisture is an important boundary condition of earth-atmosphere exchanges, and it has been defined as an essential climate variable by GCOS. Vegetation optical depth is a physical variable to measure the attenuation of vegetation in microwave radiative transfer model, and it has been proved to be a good indicator of vegetation water content and biomass. This dataset uses the multi-channel collaborative algorithm (MCCA) to retrieve both soil moisture and polarized vegetation optical depth with SMAP brightness temperature. The algorithm uses a self-constraint relationship between land parameters and an analytical relationship between brightness temperature at different channels to perform the retrieval process. The MCCA does not depend on other auxiliary data on vegetation properties and can be applied to a variety of satellites. The soil moisture product from this dataset includes the soil moisture content in the unfrozen period and the liquid water content in the frozen period. Both horizontal- and vertical-polarization vegetation optical depth are retrieved. So far as we know, it is the first polarization-dependent vegetation optical depth product at L-band. This dataset was validated by 19 dense soil moisture observation networks (9 core validation sites used by SMAP team and 13 sites not used by them), and the widely used soil climate analysis network (SCAN). It was found that ubRMSE (unbiased root mean square error) of MCCA retrieved soil moisture is generally smaller than that of other SMAP products.
ZHAO Tianjie, PENG Zhiqing , YAO Panpan, SHI Jiancheng
Based on long-term series Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover products, daily snow cover products without data gaps at 500 m spatial resolution over the Tibetan Plateau from 2002 to 2021 were generated by employing a Hidden Markov Random Field (HMRF) modeling technique. This HMRF framework optimally integrates spectral, spatiotemporal, and environmental information together, which not only fills data gaps caused by frequent clouds, but also improves the accuracy of the original MODIS snow cover products. In particular, this technology incorporates solar radiation as an environmental contextual information to improve the accuracy of snow identification in mountainous areas. Validation with in situ observations and snow cover derived from Landsat-8 OLI images revealed that these new snow cover products achieved an accuracy of 98.31% and 92.44%, respectively. Specifically, the accuracy of the new snow products is higher during the snow transition period and in complex terrains with higher elevations as well as sunny slopes. These gap-free snow cover products effectively improve the spatiotemporal continuity and the low accuracy in complex terrains of the original MODIS snow products, and is thus the basis for the study of climate change and hydrological cycling in the TP.
HUANG Yan , XU Jianghui
This data set is a code file set of TCA (triple collision analysis) algorithm, which is used to generate the global daily-scale soil moisture fusion dataset from 2011 to 2018.
XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, XIE Qiuxia, JIA Li , HU Guangcheng
The vegetation type map was created by the random forest (RF) classification approach, based on 319 ground-truth samples, combined with a set of input variables derived from the visible, infrared, and thermal Landsat-8 images. According to vegetation characteristics, four types include alpine swamp meadow (ASM), alpine meadow (AM), alpine steppe (AS), and alpine desert (AD) were classified in this map. Based on a spatial resolution of 30 m, the map can provide more detailed vegetation information.
ZHOU Defu, ZOU Defu, ZOU Defu, Zhao Lin, ZHAO Lin, Liu Guangyue, LIU Guangyue, Du Erji, DU Erji, LI Zhibin , LI Zhibin, Wu Tonghua, WU Xiaodong, CHEN Jie CHEN Jie
This data is the land cover data at 30m resolution of Southeast Asia in 2015. The data format of the data is NetCDF, and the variable name is "land cover type". The data was obtained by mosaicing and extracting the From-GLC data. Several land cover types, such as snow and ice that do not exist in Southeast Asia were eliminated.The legend were reintegrated to match the new data. The data provide information of 8 land cover types: cropland, forest, grassland, shrub, wetland, water, city and bare land. The overall accuracy of the data is 71% (Gong et al., 2019). The data can provide the land cover information of Southeast Asia for hydrological models and regional climate models.
LIU Junguo
Photosynthetically active radiation (PAR) is fundamental physiological variable driving the process of material and energy exchange, and is indispensable for researches in ecological and agricultural fields. In this study, we produced a 35-year (1984-2018) high-resolution (3 h, 10 km) global grided PAR dataset with an effective physical-based PAR model. The main inputs were cloud optical depth from the latest International Satellite Cloud Climatology Project (ISCCP) H-series cloud products, the routine variables (water vapor, surface pressure and ozone) from the ERA5 reanalysis data, aerosol from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) products and albedo from Moderate Resolution Imaging Spectroradiometer (MODIS) product after 2000 and CLARRA-2 product before 2000. The grided PAR products were evaluated against surface observations measured at seven experimental stations of the SURFace RADiation budget network (SURFRAD), 42 experimental stations of the National Ecological Observatory Network (NEON), and 38 experimental stations of the Chinese Ecosystem Research Network (CERN). The instantaneous PAR was validated at the SURFRAD and NEON, and the mean bias errors (MBEs) and root mean square errors (RMSEs) are 5.6 W m-2 and 44.3 W m-2, and 5.9 W m-2 and 45.5 W m-2, respectively, and correlation coefficients (R) are both 0.94 at 10 km scale. When averaged to 30 km, the errors were obviously reduced with RMSEs decreasing to 36.3 W m-2 and 36.3 W m-2 and R both increasing to 0.96. The daily PAR was validated at the SURFRAD, NEON and CERN, and the RMSEs were 13.2 W m-2, 13.1 W m-2 and 19.6 W m-2, respectively at 10 km scale. The RMSEs were slightly reduced to 11.2 W m-2, 11.6 W m-2, and 18.6 W m-2 when upscaled to 30 km. Comparison with the other well-known global satellite-based PAR product of the Earth's Radiant Energy System (CERES) reveals that our PAR product was a more accurate dataset with higher resolution than the CRERS. Our grided PAR dataset would contribute to the ecological simulation and food yield assessment in the future.
TANG Wenjun
Based on the data of Gaogao No. 1 and No. 2 in China from 2019 to 2020, the freeze-thaw disaster distribution data of Qinghai Tibet project corridor is produced by using the deep learning classification method and manual visual interpretation and correction. The geographical range of the data is 40km along the Xidatan Anduo section of Qinghai Tibet highway. The data include the distribution data of thermal melting lakes and ponds and the distribution data of thermal melting landslides. The data set can provide data basis for the research of freeze-thaw disaster and engineering disaster prevention and reduction in Qinghai Tibet engineering corridor. The spatial distribution of freezing and thawing disasters within 40km along the Xidatan Anduo section of Qinghai Tibet highway is self-made based on the domestic gaogao-2 image data. Firstly, the deep learning method is used to extract the mud flow terrace block from Gaogao No. 2 data; Then, ArcGIS is used for manual editing. During the production process, the operators are required to strictly abide by the operation specifications, and a special person is responsible for the quality review.
NIU Fujun, LUO Jing LUO Jing
Surface melting is the primary reason that affects the mass balance of Greenland ice sheet. At the same time, ice and snow have high albedo, and ice sheet surface melting will cause the difference of radiation energy budget, and then affects the energy exchange between sea-land-air. The high-resolution ice sheet surface melting product provides important information support for the study of Greenland ice sheet surface melting and its response to global climate change. This dataset combined microwave radiometer product and optical albedo product, the daily, winter (June-August) averages and July averages of the former are used for layer-stacking, then Gram-Schmidt Spectral Sharpening was adapted to fuse the layer-stacking results with MODIS GLASS albedo product. The spatial resolution of fusion-results has been downscaled from 25 km to 0.05˚. By employing a threshold-based melt detection approach for each fusion-results pixel, Greenland ice sheet surface melt daily product for 1985, 2000, 2015 (DSSMIS) was generated. The spatial resolution of DSSMIS is higher than that of published data sets at home and abroad. Combined with the advantages of radiometer and albedo data, the spatial details characteristics are enhanced and consistent with the extraction range of the original radiometer products, effectively reducing the noise of the radiometer. DSSMIS’s data type is integer, where 1 is melted, 0 is not melted, 255 is masked area besides Greenland ice sheet, and the data set is stored as *.nc.
WEI Siyi,
This data includes the land cover data of Central Asia, South Asia and Indochina Peninsula in the from 1992 to 2020 with a spatial resolution of 300mLand cover data includes 10 primary categories, which are combined from the secondary categories of the original data. The data source is the surface coverage product CCI-LC of ESA, where the spatial distribution of cropland, built-up land, and water for the land cover data from 1992 to 2020. Combined with the Tsinghua university global land cover data (FROM GLC, 30 m grid), NASA MODIS global land cover data (MCD12Q1, 500 m grid), the United States Geological Survey (USGS global land data (GFSAD30, 30 m), Japanese global forest data (PALSAR/PALSAR - 2, 25 m), the training sample dataset of land cover interpretation were built from the consistent areas of multiple products. The Google Earth Engine and random forest algorithm were used to correct the cropland, built-up land, and water of temporal CCI-LC data. Using the high resolution images in Google Earth at 2019 and 2020, the accuracy of change areas of cropland, built-up land, and water was validated by the stratified random sampling. A total of 3,600 land parcels were selected from 1,200 land parcels of the three land cover types, indicating that the accuracy of our corrected product increased in the range of 11% to 26% for the change areas compared to the CCI-LC product.
XU Erqi
Snow, ice, and glaciers have the highest albedo of any part of Earth's surface. The increase in melting of the polar ice sheet results in a rapid and sequential decrease in albedo and subsequently influences the global energy balance. The hydrological system derived from surface melt and basal meltwater will affect the dynamic stability of ice sheet and therefore mass balance. The dataset combined microwave radiometer product and optical albedo product, the daily, winter (June-August) averages and July averages of the former are used for layer-stacking, then Gram-Schmidt Spectral Sharpening was adapted to fuse the layer-stacking results with MODIS GLASS albedo product. The spatial resolution of fusion-results has been downscaled from 25 km to 0.05˚. By employing a threshold-based melt detection approach for each fusion-results pixel, Antarctic ice sheet surface melt daily product for 1985-1986, 2000-2001, 2015-2016 (DSSMIS) was generated. The spatial resolution of DSSMIS is higher than that of published data sets at home and abroad. Combined with the advantages of radiometer and albedo data, the spatial details characteristics are enhanced and consistent with the extraction range of the original radiometer products, effectively reducing the noise of the radiometer. It better reflects the melting gradient of mountainous area, groundline area and ice shelf over time, DSSMIS has a higher accuracy. DSSMIS’s data type is integer, where 1 is melted, 0 is not melted, 255 is masked area besides Antarctic ice sheet, and the data set is stored as *.nc.
WEI Siyi,
Under the funding of the first project (Development of Multi-scale Observation and Data Products of Key Cryosphere Parameters) of the National Key Research and Development Program of China-"The Observation and Inversion of Key Parameters of Cryosphere and Polar Environmental Changes", the research group of Zhang, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, developed the snow depth downscaling product in the Qinghai-Tibet Plateau. The snow depth downscaling data set for the Tibetan Plateau is derived from the fusion of snow cover probability dataset and Long-term snow depth dataset in China. The sub-pixel spatio-temporal downscaling algorithm is developed to downscale the original 0.25° snow depth dataset, and the 0.05° daily snow depth product is obtained. By comparing the accuracy evaluation of the snow depth product before and after downscaling, it is found that the root mean square error of the snow depth downscaling product is 0.61 cm less than the original product. The details of the product information of the Downscaling of Snow Depth Dataset for the Tibetan Plateau (2000-2018) are as follows. The projection is longitude and latitude, the spatial resolution is 0.05° (about 5km), and the time is from September 1, 2000 to September 1, 2018. It is a TIF format file. The naming rule is SD_yyyyddd.tif, where yyyy represents year and DDD represents Julian day (001-365). Snow depth (SD), unit: centimeter (cm). The spatial resolution is 0.05°. The time resolution is day by day.
YAN Dajiang, MA Ning, MA Ning, ZHANG Yinsheng
This data set is hyperspectral observation data of typical vegetation along Sichuan Tibet Railway in September 2019, using the airborne spectrometer of Dajiang M600 resonon imaging system. Including the hyperspectral data observed in the grassland area of Lhasa in 2019, with its own latitude and longitude. The hyperspectral survey was mainly sunny. Before flight, whiteboard calibration was carried out; when data were collected, there was a target (that is, the standard reflective cloth suitable for the grass), which was used for spectral calibration; there were ground mark points (that is, letters with foam plates), and the longitude and latitude coordinates of each mark were recorded for geometric precise calibration. The DN value recorded by Hyperspectral camera of UAV can be converted into reflectivity by using Spectron Pro software. Hyperspectral data is used to extract spectral characteristics of different vegetation types, vegetation classification, inversion of vegetation coverage and so on.
ZHOU Guangsheng, JI Yuhe, LV Xiaomin, SONG Xingyang
The spectral characteristics of different land use types are mainly determined by spectrograph in the surface spectral data set of Qinghai Tibet Plateau. The measured ground features are mainly divided into woodland, (Alpine) shrub, (Alpine) grassland, wetland, cultivated land and bare land. It includes the field observation points in Lhasa, Linzhi, Shigatse, Ali and Naqu. The spectral characteristics of forests were measured based on the different growth stages of vegetation; The spectral characteristics of grassland were measured based on different coverage; The spectral characteristics of cultivated land were measured based on the main crop types, rape flowers and highland barley; The measurements of wetlands were conducted on the rivers, low-lying valleys and lakes; The measurements of bare lands were conducted on the desert, Gobi and roads, which have no vegetation cover. The measurement conducted from July to August in 2019, and the data is daily observation data. The data set can provide a reference for the field verification of remote sensing interpretation.
FENG Xiaoming
This dataset includes the Antarctica ice sheet mass balance estimated from satellite gravimetry data, April 2002 to December 2019. The satellite measured gravity data mainly come from the joint NASA/DLR mission, Gravity Recovery And Climate Exepriment (GRACE, April 2002 to June 2017), and its successor, GRACE-FO (June 2018 till present). Considering the ~1-year data gap between GRACE and GRACE-FO, we extra include gravity data estimated from GPS tracking data of ESA's Swarm 3-satellite constellation. The GRACE data used in this study are weighted mean of CSR, GFZ, JPL and OSU produced solutions. The post-processing includes: replacing GRACE degree-1, C20 and C30 spherical harmonic coefficients with SLR estimates, destriping filtering, 300-km Gaussian smoothing, GIA correction using ICE6-G_D (VM5a) model, leakage reduction using forward modeling method and ellipsoidal correction.
C.K. Shum
Land surface temperature (LST) is a key variable for high temperature and drought monitoring and climate and ecological environment research. Due to the sparse distribution of ground observation stations, thermal infrared remote sensing technology has become an important means of quickly obtaining ground temperature over large areas. However, there are many missing and low-quality values in satellite-based LST data because clouds cover more than 60% of the global surface every day. This article presents a unique LST dataset with a monthly temporal resolution for China from 2003 to 2017 that makes full use of the advantages of MODIS data and meteorological station data to overcome the defects of cloud influence via a reconstruction model. We specifically describe the reconstruction model, which uses a combination of MODIS daily data, monthly data and meteorological station data to reconstruct the LST in areas with cloud coverage and for grid cells with elevated LST error, and the data performance is then further improved by establishing a regression analysis model. The validation indicates that the new LST dataset is highly consistent with in situ observations. For the six natural subregions with different climatic conditions in China, verification using ground observation data shows that the root mean square error (RMSE) ranges from 1.24 to 1.58 K, the mean absolute error (MAE) varies from 1.23 to 1.37 K and the Pearson coefficient (R2) ranges from 0.93 to 0.99. The new dataset adequately captures the spatiotemporal variations in LST at annual, seasonal and monthly scales. From 2003 to 2017, the overall annual mean LST in China showed a weak increase. Moreover, the positive trend was remarkably unevenly distributed across China. The most significant warming occurred in the central and western areas of the Inner Mongolia Plateau in the Northwest Region, and the average annual temperature change is greater than 0.1K (R>0:71, P<0:05), and a strong negative trend was observed in some parts of the Northeast Region and South China Region. Seasonally, there was significant warming in western China in winter, which was most pronounced in December. The reconstructed dataset exhibits significant improvements and can be used for the spatiotemporal evaluation of LST in high-temperature and drought-monitoring studies. More detail please refer to Zhao et al (2020). doi.org/10.5281/zenodo.3528024
MAO Kebiao
The Land Surface Temperature in China dataset contains land surface temperature data for China (about 9.6 million square kilometers of land) during the period of 2003-2017, in Celsius, in monthly temporal and 5600 m spatial resolution. It is produced by combing MODIS daily data(MOD11C1 and MYD11C1), monthly data(MOD11C3 and MYD11C3) and meteorological station data to reconstruct real LST under cloud coverage in monthly LST images, and then a regression analysis model is constructed to further improve accuracy in six natural subregions with different climatic conditions.
MAO Kebiao
This dataset contains measurements of L-band brightness temperature by an ELBARA-III microwave radiometer in horizontal and vertical polarization, profile soil moisture and soil temperature, turbulent heat fluxes, and meteorological data from the beginning of 2016 till August 2019, while the experiment is still continuing. Auxiliary vegetation and soil texture information collected in dedicated campaigns are also reported. This dataset can be used to validate the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite based observations and retrievals, verify radiative transfer model assumptions and validate land surface model and reanalysis outputs, retrieve soil properties, as well as to quantify land-atmosphere exchanges of energy, water and carbon and help to reduce discrepancies and uncertainties in current Earth System Models (ESM) parameterizations. ELBARA-III horizontal and vertical brightness temperature are computed from measured radiometer voltages and calibrated internal noise temperatures. The data is reliable, and its quality is evaluated by 1) Perform ‘histogram test’ on the voltage samples (raw-data) of the detector output at sampling frequency of 800 Hz. Statistics of the histogram test showed no non-Gaussian Radio Frequency Interference (RFI) were found when ELBAR-III was operated. 2) Check the voltages at the antenna ports measured during sky measurements. Results showed close values. 3) Check the instrument internal temperature, active cold source temperature and ambient temperature. 3) Analysis the angular behaviour of the processed brightness temperatures. -Temporal resolution: 30 minutes -Spatial resolution: incident angle of observation ranges from 40° to 70° in step of 5°. The area of footprint ranges between 3.31 m^2 and 43.64 m^2 -Accuracy of Measurement: Brightness temperature, 1 K; Soil moisture, 0.001 m^3 m^-3; Soil temperature, 0.1 °C -Unit: Brightness temperature, K; Soil moisture, m^3 m^-3; Soil temperature, °C/K
BOB Su, WEN Jun
The dataset include ground-based passive microwave brightness temperature, multi-angle brightness temperature, ten-minute 4-component radiation and snow temperature, daily snow pit data and hourly meteorological data observed at Altay base station（lon：88.07、lat: 44.73）from November 27, 2015 to March 26, 2016. Daily snow pit parameters include: snow stratification, stratification thickness, density, particle size, temperature. These data are stored in five NetCDF files: TBdata. nc, TBdata-multiangle. nc, ten-minute 4 component radiation and snow temperature. nc, hourly meteorological and soil data. nc and daily snow pit data.nc. TBdata. nc is brightness temperature at 3 channels for both polarizations automatically collected by a six-channel dual polarized microwave radiometer RPG-6CH-DP. The contents include Year, month, day, hour, minute, second, Tb1h, Tb1v, Tb18h, Tb18v, Tb36h, Tb36v, incidence angle, azimuth angle. TBdata-multiangle.nc is 7 groups of multi-angle brightness temperatures at 3 channels for both polarizations. The contents include Year, month, day, hour, minute, second, Tb1h, Tb1v, Tb18h, Tb18v, Tb36h, Tb36v, incidence angle, azimuth angle. The ten-minute 4 component radiation and snow temperature.nc contains 4 component radiation and layered snow temperatures. The contents include Year, month, day, hour, minute, SR_DOWN, SR_UP, LR_DOWN, LR_UP, T_Sensor, ST_0cm, ST_5cm, ST_15cm, ST_25cm, ST_35cm, ST_45cm, ST_55cm. The hourly meteorological and soil data.nc contains hourly weather data and layered soil data. The contents include Year, month, day, hour, Tair, Wair, Pair, Win, SM_10cm, SM_20cm, Tsoil_5cm, Tsoil_10cm, Tsoil_15 cm, Tsoil_20cm. The daily snow pit data.nc. is manual snow pit data. The observation time was 8:00-10:100 am local time. The contents include Year, month, day, snow depth, thickness_layer1, thickness_layer2, thickness_layer3, thickness_layer4, thickness_layer5, thickness_layer6, Long_layer1, Short_layer1, Long_layer2, Short_layer2, Long_layer3, Short_layer3, Long_layer 4, Short_layer4, Long_layer5, Short_layer5, Long_layer6, Short_layer 6, Stube, Snow shovel_0-10, Snow shovel _10-20, Snow shovel _20-30, Snow shovel _30-40, Snow shovel _40-50, Snow fork_5, Snow fork _10, Snow fork _15, Snow fork_20, Snow fork_25, Snow fork_30, Snow fork_35, Snow fork_40, Snow fork_45, Snow fork_50, shape1, shape2, shape3, shape4, shape5,
DAI Liyun
The Tibetan Plateau Glacier Data –TPG2017 is a glacial coverage data on the Tibetan Plateau from selected 210 scenes of Landsat 8 Operational Land Imager (OLI) images with 30-m spatial resolution from 2013 to 2018, among of which 90% was in 2017 and 85% in winter. Therefore, 2017 was defined as the reference year for the mosaic image. Glacier outlines were digitized on-screen manually from the 2017 image mosaic, relying on false-colour image composites (RGB by bands 654), which allowed us to distinguish ice/snow from cloud. Debris-free ice was distinguished from the debris and debris-covered ice by its higher reflectance. Debris-covered ice was not delineated in this data. The delineated glacier outlines were compared with band-ratio (e.g. TM3/TM5) results, and validated by overlapping them onto Google Earth imagery, SRTM DEM, topographic maps and corresponding satellite images. For areas with mountain shadows and snow cover, they were verified by different methods using data from different seasons. For glaciers in deep shadow, Google EarthTM imagery from different dates was used as the reference for manual delineation. Steep slopes or headwalls were also excluded in the TPG2017. Areas that appeared in any of these sources to have the characteristics of exposed ground/basement/bed rock were manually delineated as non-glacier, and were also cross-checked with CGI-1 and CGI-2. Steep hanging glaciers were included in TPG2017 if they were identifiable on images in all other three epochs (i.e. TPG1976, TPG2001, and TPG2013). The accuracy of manual digitization was controlled within one half-pixel. All glacier areas were calculated on the WGS84 spheroid in an Albers equal-area map projection centred at (95°E, 30°N) with standard parallels at 15°N and 65°N. Our results showed that the relative deviation of manual interpretation was less than 3.9%.
YE Qinghua
The permafrost stability map was created based on the classification system proposed by Guodong Cheng (1984), which mainly depended on the inter-annual variation of deep soil temperature. By using the geographical weighted regression method, many auxiliary data was fusion in the map, such as average soil temperature, snow cover days, GLASS LAI, soil texture and organic from SoilGrids250, soil moisture products from CLDAS of CMA, and FY2/EMSIP precipitation products. The permafrost stability data spatial resolution is 1km and represents the status around 2010. The following table is the permafrost stability classification system. The data format is Arcgis Raster.
RAN Youhua
This dataset includes component temperatures measured by the thermal imager at the Mixed Forest and Sidaoqiao stations between 23 July and 18 August, 2014. The Mixed Forest (101.1335 °E, 41.9903 °N, 874 m.a.s.l.) and Sidaoqiao (101.1374 °E, 42.0012 °N, 873 m.a.s.l.) stations were located in the downstream of the Heihe River basin, Dalaihubu Town, Ejin Banner, Inner Mongolia. At the Mixed Forest station, a Testo 890-2 thermal imager (Testo Inc., Germany) with a resolution of 640 × 480 pixels was employed to acquire brightness temperature images. The imager was manually operated from a 10-m height platform of the tower between 10:00-16:00 (China Standard Time, CST) with an observation interval of 1-h on cloudless days. On August 4th observations were acquired between 11:00 and 17:00 at an interval of 10-min to match observations acquired with an airborne TIR imager. The ground based imager was pointed to five viewing directions (southeast-SE, east-E, northeast-NE, northwest-NW, and southwest-SW) and was inclined 25°–45° below the horizon depending on viewing direction. At Sidaoqiao station, a Testo 875-2i imager (Testo Inc., Germany) with a resolution of 160 × 120 pixels was manually operated from a 10-m high platform to acquire brightness temperature images in directions SW, SE, NE, and NW. Depending on the targets in each viewing direction, the imager was inclined to 30°–45° below the horizon. Observations at Sidaoqiao and Mixed Forest stations were almost synchronous. Furthermore, visible images were taken simultaneously with the aforementioned two TIR imagers (2048 × 1536 pixels for Testo 890-2 and 640 × 480 pixels for Testo 875-2i).
LI Mingsong , MA Jin
This dataset includes component temperatures measured by the thermal infrared (TIR) radiometers at the Mixed Forest and Sidaoqiao stations between 22 July, 2014 and 19 July, 2016. The Mixed Forest (101.1335 °E, 41.9903 °N, 874 m.a.s.l.) and Sidaoqiao (101.1374 °E, 42.0012 °N, 873 m.a.s.l.) stations were located in the downstream of the Heihe River basin, Dalaihubu Town, Ejin Banner, Inner Mongolia. At the Mixed Forest station, two TIR radiometers (SI-111, Apogee Instruments Inc., USA) connected to a data logger (CR800, Campbell Scientific Inc., USA) measured component temperatures of the sunlit canopy and shaded canopy. TIR radiometers were mounted horizontally at 5 m height on iron rods just south and north of a tree and pointed to its canopy. The distance from the sensor to the canopy was ~1 m. At the Sidaoqiao station, two SI-111 TIR radiometers connected to a CR800 data logger measured component temperatures of the soil and shrub. The first sensor pointed from 2 m height under a viewing zenith angle of 45° to bare soil; the second sensor was mounted at 1-m height and pointed horizontally into the shrub canopy.
ZHOU Ji, LI Mingsong , MA Jin
This dataset is blended by two other sets of data, snow cover dataset based on optical instrument remote sensing with 1km spatial resolution on the Qinghai-Tibet Plateau (1989-2018) produced by National Satellite Meteorological Center, and near-real-time SSM/I-SSMIS 25km EASE-grid daily global ice concentration and snow extent (NISE, 1995-2018) provided by National Snow and Ice Data Center (NSIDC, U.S.A). It covers the time from 1995 to 2018 (two periods, from January to April and from October to December) and the region of Qinghai-Tibet Plateau (17°N-41°N, 65°E-106°E) with daily product, which takes equal latitude and longitude projection with 0.01°×0.01° spatial resolution, and characterizes whether the ground is covered by snow. The input data sources include daily snow cover products generated by NOAA/AVHRR, MetOp/AVHRR, and alternative to AVHRR taken from TERRA/MODIS corresponding observation, and snow extent information of NISE derived from observation by SSM/I or SSMIS of DMSP satellites. The processing method of data collection is as following: first, taking 1km snow cover product from optical instruments as initial value, and fully trusting its snow and clear sky without snow information; then, under the aid of sea-land template with relatively high resolution, replacing the pixels or grids where is cloud coverage, no decision, or lack of satellite observation, by NISE's effective terrestrial identification results. For some water and land boundaries, there still may be a small amount of cloud coverage or no observation data area that can’t be replaced due to the low spatial resolution of NISE product. Blended daily snow cover product achieves about 91% average coincidence rate of snow and non-snow identification compared to ground-based snow depth observation in years. The dataset is stored in the standard HDF4 files each having two SDSs of snow cover and quality code with the dimensions of 4100-column and 2400-line. Complete attribute descriptions is written in them.
ZHENG Zhaojun, CAO Guangzhen
Snow cover dataset is produced by snow and cloud identification method based on optical instrument observation data, covering the time from 1989 to 2018 (two periods, from January to April and from October to December) and the region of Qinghai-Tibet Plateau (17°N-41°N, 65°E-106°E) with daily product, which takes equal latitude and longitude projection with 0.01°×0.01° spatial resolution, and characterizes whether the ground under clear sky or transparent thin cloud is covered by snow. The input data sources include AVHRR L1 data of NOAA and MetOp serials of satellites, and L1 data corresponding to AVHRR channels taken from TERRA/MODIS. Decision Tree algorithm (DT) with dynamic thresholds is employed independent of cloud mask and its cloud detection emphasizes on reserving snow, particularly under transparency cirrus. It considers a variety of methods for different situations, such as ice-cloud over the water-cloud, snow in forest and sand, thin snow or melting snow, etc. Besides those, setting dynamic threshold based on land-surface type, DEM and season variation, deleting false snow in low latitude forest covered by heavy aerosol or soot, referring to maximum monthly snowlines and minimum snow surface brightness temperature, and optimizing discrimination program, these techniques all contribute to DT. DT discriminates most snow and cloud under normal circumstances, but underestimates snow on the Qinghai-Tibet Plateau in October. Daily product achieves about 95% average coincidence rate of snow and non-snow identification compared to ground-based snow depth observation in years. The dataset is stored in the standard HDF4 files each having two SDSs of snow cover and quality code with the dimensions of 4100-column and 2400-line. Complete attribute descriptions is written in them.
ZHENG Zhaojun, CHU Duo
On 1 August 2012, Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Harbin Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area. WIDAS includes a CCD camera with a spatial of resolution 0.08 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 0.4 m), and a thermal image camera with a spatial resolution of 2 m. The CCD camera data are recorded in DN values processed by mosaic and orthorectification.
XIAO Qing, Wen Jianguang
On 3 August 2012, Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Harbin Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area (5×5 km). WIDAS includes a CCD camera with a spatial resolution of 0.08 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 0.4 m), and a thermal image camera with a spatial resolution of 2 m. The CCD camera data are recorded in DN values processed by mosaic and orthorectification.
XIAO Qing, Wen Jianguang
On 26 July 2012, Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Harbin Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area (5×5 km). WIDAS includes a CCD camera with a spatial resolution of 0.2 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 1 m), and a thermal image camera with a spatial resolution of 4.8 m. The CCD camera data are recorded in DN values processed by mosaic and orthorectification.
XIAO Qing, Wen Jianguang
On 25 August and 28 August, 2012, a RCD30 camera of Leica Company boarded on the Y-12 aircraft was used to obtain CCD image. RCD30 camera has a focal length of 80 mm and four bands including red, green, blue and near-infrared bands. The absolute flight altitude is 4800 and 5200 m, and ground sample distance is 6-19 cm. The product includes TIF images and exterior orientation elements.
XIAO Qing, Wen Jianguang
On 19 August 2012, a RCD30 camera of Leica Company boarded on the Y-12 aircraft was used to obtain the CCD image. RCD30 camera has a focal length of 80 mm and four bands including red, green, blue and near-infrared bands. The absolute flight altitude is 2900 m and ground sample distance is 10 cm. The data includes TIF images and exterior orientation elements.
XIAO Qing, Wen Jianguang
On 2 August 2012, Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Harbin Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area (30×30 km). WIDAS includes a CCD camera with a spatial resolution of 0.26 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 1.3 m), and a thermal image camera with a spatial resolution of 6.3 m. The CCD camera data are recorded in DN values processed by mosaic and orthorectification.
XIAO Qing, Wen Jianguang
On 25 August 2012, a RCD30 camera of Leica Company boarded on the Y-12 aircraft was utilized to obtain the optical remote sensing data. RCD30 camera has a focal length of 80 mm and four bands including red, green, blue and near-infrared bands. The absolute flight altitude is 5200 m and ground sample distance is 6-19 cm. The product includes TIF images and exterior orientation elements.
XIAO Qing, Wen Jianguang
On 25 August 2012, a RCD30 camera of Leica Company boarded on the Y-12 aircraft was used to obtain the CCD image. RCD30 camera has a focal length of 80 mm and four bands including red, green, blue and near-infrared bands. The absolute flight altitude is 4800 m and 5200 m, and ground sample distance is 8-19 cm. The product includes TIF images and exterior orientation elements.
XIAO Qing, Wen Jianguang
This is the vegetation index (NDVI) for Maduo County in July, August and September of 2016. It is obtained through calculation based on the multispectral data of GF-1. The spatial resolution is 16 m. The GF-1 data are processed by mosaicking, projection coordinating, data subsetting and other methods. The maximum synthesis is then conducted every month in July, August, and September.
LI Fei, Fei Li, Zhijun Zhang
This is the 1976, 1991, 2000, and 2010 vector data set of glaciers and glacial lakes in the Boqu Basin in Central Himalaya based on Landsat satellite images. The data source is from Landsat remote images. 1976: LM21510411975306AAA05, LM21510401976355AAA04 1991: LT41410401991334XXX02, LT41410411991334XXX02 2000: LE71410402000279SGS00, LE71400412000304SGS00, LE71410402000327EDC00, LE71410412000327EDC00 2010: LT51400412009288KHC00, LT51410402009295KHC00, LT51410412009311KHC00, LT51410402011237KHC00. The boundaries of glaciers and glacial lakes are extracted manually from the various remote sensing images. The extraction error of the boundaries of glaciers and glacial lakes is estimated to be 0.5 pixels. Data file: Glacial_1976: Glacier vector data in 1976 Glacial_1991: Glacier vector data in 1991 Glacial_2000: Glacier vector data in 2000 Glacial_2010: Glacier vector data in 2010 Glacial_Lake_1976: Glacial lake vector data in 1976年 Glacial_Lake_1991: Glacial lake vector data in 1991 Glacial_Lake_2000: Glacial lake vector data in 2000 Glacial_Lake_2010: Glacial lake vector data in 2010 The glacial lake vector data fields include Number, name, latitude and longitude, altitude, area, orientation, type of glacial lake, length, width, and distance from the glacier.
WANG Weicai
The Tibetan Plateau Glacier Data –TPG2013 is a glacial coverage data on the Tibetan Plateau around 2013. 128 Landsat 8 Operational Land Imager (OLI) images were selected with 30-m spatial resolution, for comparability with previous and current glacier inventories. Besides, about 20 images acquired in 2014 were used to complete the full coverage of the TP. The most frequent year in this period was defined as the reference year for the mosaic image: i.e. 2013. Glacier outlines were digitized on-screen manually from the 2013 image mosaic, relying on false-colour image composites (RGB by bands 654), which allowed us to distinguish ice/snow from cloud. Debris-free ice was distinguished from the debris and debris-covered ice by its higher reflectance. Debris-covered ice was not delineated in this data. [To minimize the effects of snow or cloud cover on glacierized areas, high-resolution (30 m spatial resolution and 4-day repetition cycle) images were also used for reference in glacier delineation from the Chinese satellites HJ-1A and HJ-1B, which were launched on Sep.6th 2008. Both carried as payload two 4-band CCD cameras with swath width 700 km (360 km per camera). All HJ-1A/1B data in 2012, 2013 and 2014 (65 scenes, Fig.S1, Table S1) were from China Centre for Resources Satellite Data and Application (CRESDA; http://www.cresda.com/n16/n92006/n92066/n98627/index.html). Each scene was orthorectified with respect to the 30m-resolution digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) and Landsat images.] The delineated glacier outlines were compared with band-ratio (e.g. TM3/TM5) results, and validated by overlapping them onto Google Earth imagery, SRTM DEM, topographic maps and corresponding satellite images. Topographic maps from the 1970s and all available satellite images (including Google EarthTM imagery and HJ-1A/1B satellite data) were used as base reference data. For areas with mountain shadows and snow cover, they were verified by different methods using data from different seasons. For glaciers in deep shadow, Google EarthTM imagery from different dates was used as the reference for manual delineation. Steep slopes or headwalls were also excluded in the TPG2013. Areas that appeared in any of these sources to have the characteristics of exposed ground/basement/bed rock were manually delineated as non-glacier, and were also cross-checked with CGI-1 and CGI-2. Steep hanging glaciers were included in TPG2013 if they were identifiable on images in all three epochs (i.e. TPG1976, TPG2001, and TPG2013). The accuracy of manual digitization was controlled within one half-pixel. All glacier areas were calculated on the WGS84 spheroid in an Albers equal-area map projection centred at (95°E, 30°N) with standard parallels at 15°N and 65°N. Our results showed that the relative deviation of manual interpretation was less than 3.9%.
YE Qinghua
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
