Iceberg calving, one of the key process of Antarctic mass balance, has been regarded as an important variable in fine monitoring the changes of ice shelves. The authors used multi-source remote sensing data near early August of each year from 2005 to 2020, including ENVISAT ASAR (WSM) images from 2005 to 2011, Terra/Aqua MODIS 7-2-1 band composite images from 2012 to 2014, Landsat-8 OLI 4-3-2 band composite images from 2013 to 2020, and Sentinel-1 SAR (EW) images from 2015 to 2020, to generate annual circum-Antarctic image mosaics after pre-processing. Next, combining MEaSUREs ice velocity dataset, grounding line, ice thickness dataset (Bedmap 2 and Bedmachine), spatial calculation and map digitization techniques were applied to extract all annual calving events larger than 1 km² that occurred on the Antarctic ice shelves from August 2005 to August 2020. Also, their area, thickness, mass and calving recurrence cycle were calculated to derive the annual iceberg calving dataset of the Antarctic ice shelves (2005-2020). This dataset contains the distribution of 15-year annual calving events, along with the attributes of each individual calving event including calving year, length, area, average thickness, mass, and recurrence interval. This dataset can directly reflect the magnitude characteristics and distribution of Antarctic iceberg calving in different years, which fills the gap of fine monitoring dataset of iceberg calving and provides fundamental data for subsequent research on calving mechanism and mass balance of Antarctic ice shelf-ice sheet system.
QI Mengzhen, LIU Yan, CHENG Xiao, HUI Fengming, CHEN Zhuoqi
Qiangyong glacier: 90.23 °E, 28.88° N, 4898 m asl. The surface is bedrock. The record contains data of 1.5 m temperature, 1.5 m humidity, 2 m wind speed, 2 m wind orientation, surface temperature, etc. Data from the automated weather station was collected using USB equipment at 19:10 on August 6, 2019, with a recording interval of 10 minutes, and data was downloaded on December 20, 2020. There is no missing data but a problem with the wind speed data after 9:30 on July 14, 2020 (most likely due to damage to the wind vane). Jiagang glacier: 88.69°E, 30.82°N, 5362 m asl. The surface is rubble and weeds. The records include 1.5 meters of temperature, 1.5 meters of humidity, 2 meters of wind speed, 2 meters of wind direction, surface temperature, etc. The initial recording time is 15:00 on August 9, 2019, and the recording interval is 1 minute. The power supply is mainly maintained by batteries and solar panels. The automatic weather station has no internal storage. The data is uploaded to the Hobo website via GPRS every hour and downloaded regularly. At 23:34 on January 5, 2020, the 1.5 meter temperature and humidity sensor was abnormal, and the temperature and humidity data were lost. The data acquisition instrument will be retrieved on December 19, 2020 and downloaded to 19:43 on June 23, 2020 and 3:36 on September 25, 2020. Then the temperature and humidity sensors were replaced, and the observations resumed at 12:27 on December 21. The current data consists of three segments (2019.8.9-2020.6.30; 2020.6.23-2020.9.25; 2020.12.19-2020.12.29), Some data are missing after inspection. Some data are duplicated in time due to recording battery voltage, which needs to be checked. The meteorological observation data at the front end of Jiagang mountain glacier are collected by the automatic weather station Hobo rx3004-00-01 of onset company. The model of temperature and humidity probe is s-thb-m002, the model of wind speed and direction sensor is s-wset-b, and the model of ground temperature sensor is s-tmb-m006. The meteorological observation data at the front end of Jianyong glacier are collected by the US onset Hobo u21-usb automatic weather station. The temperature and humidity probe model is s-thb-m002, the wind speed and direction sensor model is s-wset-b, and the ground temperature sensor model is s-tmb-m006.
The data in the form of .xlsx store the meteorological varialbes observed on the East Rongbuk glacier from May to July. Two sheets, named "Surface_energy_budget" and "Cycle", respectivley, are included. In the sheet of "surface_energy_budget", the meteorological variables are as follows: Four-component radiations (incident solar radiation, reflected shortwave radiation, incoming longwave radiation, outgoing longwave radiation)、wind speed and direction, air temperature and relative humidity, cloud index, south Asian summer monsoon and albedo. In addition, net shortwave radiation, net longwave radiation, net radiation, sensible heat, latent heat and subsurface heat are also included. Energy fluxes are in unit of W m-2. The sheet of "Cycle" stores the diurnal cycle of the meteorological variables mentioned above. In the first line, the prefixes of "1"、"2" and “3” indicate three observational periods, i.e., "1" represents days from 1 - 28 May, "2" represents the period between 29 May 16 June and "3" indicates time episode from 17 June to 22 July.
This data set is the physical property data of Hengduan Mountain Glacier, reflecting the temperature condition of Hengduan Mountain Glacier. It was observed on Baishui No.1 glacier on the east slope of Yulong Mountain and dagongba glacier on the west slope of Gongga Mountain by the comprehensive scientific investigation team of Qinghai Tibet Plateau of Chinese Academy of Sciences from 1982 to 1984. The temperature field location, altitude, drilling information, ice surface condition, sampling time, sampling depth and measured temperature of Baishui No. 1 glacier on the east slope of Yulong and dagongba glacier on the west slope of Gongga are recorded in detail in the data, which are obtained from field investigation and calculation. At the same time, the velocity data of dagongba glacier and the surface strain rate, normal strain rate and its error and principal strain rate at 4700m of Baishui No.1 glacier in Yulongshan are available. This data is of great significance to the study of temperature and movement of glacial active layer in Hengduan Mountain area.
The melting observation of Hengduan Moutain glacier is mainly carried out on Hailuogou Glacier on the east slope of Gongga and the large and small Gongba glacier on the west slope of Gongga. In addition, some ablation observations have been made on Baishui 1 glacier on the east slope of Yulong. According to the melting observation of the four glaciers in the above two mountains, there are some regional representativeness, which makes them reflect the basic situation of melting glaciers in Hengduan Mountain. This data set records the glacier ablation data of different time and different places: from June to August 1982, the Glacier No. 1 in Baishui on the east slope of Yulong mountain was observed at the altitude of 4200m, 4600m and 4800m. From August 27, 1982 to the end of August 1983, the annual measured data of different heights of Hailuogou Glacier tongue on the east slope of Gongga Mountain were collected. From July 12, 1982 to August 6, 1983, the observation data of Gongba glacier melting on the west slope of Gongga Mountain were recorded.
This data is the statistics of the glaciers and their types in Hengduan Mountain area, the information of each glacier, and the data of some glacier snow lines and related parameters in China. The data includes eight tables, which are glacier Statistics (measured data) of Hengduan Mountains, glacier Statistics (measured data) of Hengduan Mountains, glacier types (measured data), basic characteristics (measured data) of some glacier recharge areas in Gongga Mountain, AAR value and avalanche area (measured data) of some glaciers in Gongga mountain, and ice field in Gongga mountain Data statistics of Sichuan (measured data), thickness measurement statistics of 4 glaciers in Gongga Mountain (measured data), snow line data of some glaciers in China and related parameters (data statistics).
The data set is a record of glacier distribution in Hoh Xil region, including three tables: the distribution of modern glaciers in various mountain areas in Hoh Xil region, the distribution of modern glaciers in various river basins in Hoh Xil region, and the distribution of modern glaciers in different mountain height segments in Hoh Xil region. Hoh Xil, located in the hinterland of the Qinghai Tibet Plateau, has an average altitude of more than 5000m and a very cold climate. According to the catalogue of China's glaciers and the author's re statistics on the 1 / 100000 topographic map, 437 modern glaciers are developed in the whole region, covering an area of 1552.39 square kilometers, with ice reserves of 162.8349 cubic kilometers, becoming an important source of water supply for many rivers and lakes in the region. Through this data set, we can know more about the distribution of glaciers in this area.
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) These data main included the GPR-surveyed ice thickness of six typical various-sized glaciers in 2016-2018; the GlabTop2-modeled ice thickness of the entire UIB sub-basins, discharge data of the hydrological stations, and related raw & derived data. 2) Data sources and processing methods: We compared the plots and profiles of GPR-surveyed ice bed elevation with the GlabTop2-simulated results and selected the optimal parametric scheme, then simulated the ice thickness of the whole UIB basin and assessed its hydrological effect. These processed results were stored as tables and tif format， 3) Data quality description: The simulated ice thickness has a spatial resolution of 30 m, and has been verified by the GPR-surveyed ice thickness for the MD values were less than 10 m. The maximum error of the GPR-measured data was 230.2 ± 5.4 m, within the quoted glacier error at ± 5%. 4) Synthesizing knowledge of the ice thickness and ice reserves provides critical information for water resources management and regional glacial scientific research, it is also essential for several other fields of glaciology, including hydrological effect, regional climate modeling, and assessment of glacier hazards.
This dataset includes annual mosaics of Antarctic ice velocity derived from Landsat 8 images between December, 2013 and April, 2019, which was updated in 2020 in order to produce multi-year annual ice velocity mosaics and improve the quality of products including non-local means (NLM) filter, and absolute calibration using rock outcrops data. The resulting Version 2 of the mosaics offer reduced local errors, improved spatial resolution as described in the README file.
The data are collected from the automatic weather station (AWS, Campbell company) in the moraine area of the 24K glacier in the Southeast Tibet Plateau, Chinese Academy of Sciences. The geographic coordinates are 29.765 ° n, 95.712 ° E and 3950 m above sea level. The data include daily arithmetic mean data of air temperature (℃), relative humidity (%), wind speed (M / s), net radiation (w / m2), water vapor pressure (kPa) and air pressure (mbar). In the original data, an average value was recorded every 30 minutes before October 2018, and then an average value was recorded every 10 minutes. The temperature and humidity are measured by hmp155a temperature and humidity probe. The net radiation probe is nr01, the atmospheric pressure sensor probe is ptb210, and the wind speed sensor is 05103. These probes are 2 m above the ground. Data quality: the data has undergone strict quality control. The original abnormal data of 10 minutes and 30 minutes are removed first, and then the arithmetic mean of each hour is calculated. Finally, the daily value is calculated. If the number of hourly data is less than 24, the data is removed, and the corresponding date data in the data table is empty. In addition to the lack of some parameter data due to the thick snow and low temperature in winter and spring, the data can be used by scientific researchers who study climate, glacier and hydrology through strict quality control.
The data set contains the stable oxygen isotope data of ice core from 1864 to 2006. The ice core was obtained from Noijinkansang glacier in the south of Southern Tibetan Plateau, with a length of 55.1 meters. Oxygen isotopes were measured using a MAT-253 mass spectrometer (with an analytical precision of 0.05 ‰) at the Key Laboratory of CAS for Tibetan Environment and Land Surface Processes, China. Data collection location: Noijinkansang glacier (90.2 ° e, 29.04 ° n, altitude: 5950 m)
The ages of glacial traces of the last glacial maximum, Holocene and little ice age in the Westerlies and monsoon areas were determined by Cosmogenic Nuclide (10Be and 26Al) exposure dating method to determine the absolute age sequence of glacial advance and retreat. The distribution of glacial remains is investigated in the field, the location of moraine ridge is determined, and the geomorphic characteristics of moraine ridge are measured. According to the geomorphic location and weathering degree of glacial remains, the relationship between the new and the old is determined, and the moraine ridge of the last glacial maximum is preliminarily determined. The exposed age samples of glacial boulders on each row of moraine ridges were collected from the ridge upstream. This data includes the range of glacier advance and retreat in Karakoram area during climate transition period based on 10Be exposure age method.
Among many indicators reflecting climate and environmental change, the stable isotope index of ice core is an indispensable parameter in the study of ice core record, and is one of the most reliable and effective ways to recover the past climate change. Ice core accumulation is a direct record of precipitation on glaciers, and high resolution ice core records ensure the continuity of precipitation records. Therefore, ice core records provide an effective means to recover precipitation changes. The isotope and accumulation of ice cores drilled from the Qinghai Tibet Plateau can be used to reconstruct the changes of temperature and precipitation, which is a good record of climate and environment. This data set provides stable isotope records of hushe ice core in Karakoram area and provides data support for the study of climate change in Qinghai Tibet Plateau.
XU Baiqing, WANG Mo
The coverage time of glacier runoff data set in the five major river source areas of the Qinghai Tibet Plateau is from 1971 to 2015, and the time resolution is year by year, covering the source areas of five major rivers (Yellow River source, Yangtze River source, Lancang River source, Nu River source, Yarlung Zangbo River source). The data is based on multi-source remote sensing and measured data. The glacier runoff data is simulated by using the daily scale meteorological data of five major river source areas and their surrounding meteorological stations, the global vegetation products of umd-1km, the igbp-dis soil database, the first and second glacier catalogue data, and the distributed hydrological model vic-cas coupled with the glacier module is used to simulate the glacier runoff data. The simulation results are verified by the site measured data to enhance the quality control. Data indicators include: Glacier runoff (rate of glacier runoff:%), total runoff (mm / a), snow runoff (rate of snow runoff:%), and rainfall runoff rate (rainfall runoff rate:%).
The data involved three periods of geodetic glacier mass storage change of three Rongbuk glaciers and its debris-covered ice in the Rongbuk Catchment from 1974-2016 (unit: m w.e. a-1). It is stored in the ESRI vector polygon format. The data sets are composed of three periods of glacier surface elevation difference between 1974-2000，2000-2016 and 1974-2006, i.e. DHPRISM2006-DEM1974（DH2006-1974）、DHSRTM2000-DEM1974（DH2000-1974）、DHASTER2016-SRTM2000（DH2016-2000）. DH2006-1974 was surface elevation change between ALOS/PRISMDEM(PRISM2006) and DEM1974, i.e. the DEM1974 was subtracted from PRISM2006, DH2006-1974 =PRISM2006 – DEM1974. The PRISM2006 was generated from stereo pairs of ALOS/PRISM on 4 Dec. 2006. The earlier historical DEM (DEM1974, spatial resolution 25m) was derived from 1:50,000 topographic maps in October 1974(DEM1974,spatial resolution 25m). The uncertainty in the ice free areas of DHPRISM2006-DEM1974 was ±0.24 m a-1. DHSRTM2000-DEM1974（DH2000-1974）was surface elevation change between SRTM DEM(SRTM2000) and DEM1974. The uncertainty in the ice free areas of DHSRTM2000-DEM1974 was ±0.13 m a-1. DHASTER2016-SRTM2000（DH2016-2000）was the surface elevation change between ASTER DEM2016 and SRTM DEM(SRTM2000). The uncertainty in the ice free areas of DHASTER2016-SRTM2000 was ±0.08 m a-1. Glacier-averaged annual mass balance change (m w.e.a-1) was averaged annually for each glacier, which was calculated by DH2006-1974/DH2000-1974/DH2016-2000, glacier coverage area and ice density of 850 ± 60 kg m−3. The attribute data includes Glacier area by Shape_Area (m2), EC2000-1974/EC2016-2000/EC2006-1974, i.e. Glacier-averaged surface elevation change in each period(m a-1), MB2000-1974/ MB2016-2000/MB2006-1974, i.e. Glacier-averaged annual mass balance in each period (m w.e.a-1), and MC2000-1974/ MC2016-2000/MC2006-1974,Glacier-averaged annual mass change in each period(m3 w.e.a-1), Uncerty_EC is the maximum uncertainty of glacier surface elevation change（m a-1）、Uncerty_MB, is the maximum uncertainty of glacier mass balance（m w.e. a-1），Uncerty_MC, is the maximum uncertainty of glacier mass change（m3w.e. a-1）。 MinUnty_EC，is the minimum uncertainty of glacier surface elevation change，MinUnty_MB，is the minimum uncertainty of glacier mass balance（m w.e. a-1），MinUnty_MC is the minimum uncertainty of glacier mass change（m3 w.e. a-1.The data sets could be used for glacier change, hydrological and climate change studies in the Himalayas and High Mountain Asia.
The data involved geodetic glacier mass change of 71pieces of glaciers during 2000-2014 in the east of the Yigongzangbu, Southeast Tibetan Plateau. It is stored in the ESRI vector polygon format.Glacier-averaged mass balance (m w.e.a-1) was calculated by the surface elevation difference between 2000-2014 ( Dh2000-2014)、glacier coveraged vector data (CGI2/TPG1976/RGI6.0) and ice density of 850 ± 60 kg m−3. Dh2000-2014 is obtained from surface elevation change by D-InSAR technique from a pair of TSX / TDx SAR images on February 7, 2014 and SRTM DEM. CGI2/TPG1976/RGI6.0 were used to extract glacier boundary and GLIMS-ID. SRTM DEM is the reference DEM and datum DEM with spatial resolution 30m. The attribute data includes GLIMS-ID, Area,EC_m_a-1,,MB_m w.e.a-1, MC_m3 w.e.a-1, MC_Gt.a-1, Uncerty_EC, Uncerty_MB, UT_MCm3w.e. a-1. Respectively, EC_m_a-1,,is the glacier-averaged annual elevation change during 2000-2014(m a-1),MB_m w.e.a-1, is glacier-averaged annual mass balance during 2000-2014(m w.e.a-1), MC_m3 w.e.a-1, is glacier-averaged annual mass change during 2000-2014 (m3 w.e.a-1), MC_Gt.a-1,is glacier-averaged annual mass change during 2000-2014 (Gt.a-1)Uncerty_EC is the uncertainty of glacier surface elevation change（±m a-1）、Uncerty_MB ,is the uncertainty of glacier mass balance（±m w.e. a-1），UT_MCm3w.e. a-1, is the uncertainty of glacier mass change（±m3w.e. a-1）。The data sets could be used for glacier change, hydrological and climate change studies in the southeast of Tibetan Plateau.
The data set involved geodetic annual glacier-averaged mass balance and mass change data atMt.Xixiabangma areasin the Himalayas from 1974 to 2017. It is stored in the ESRI vector polygon format and is composed of two periods, which includes surface elevation difference between 1974-2000 (DH1974-2000, from KH-9 DEM1974 and SRTM DEM2000), surface elevation difference between 2000-2017(DH2000-2017, by DinSAR techniquesfrom SRTM DEM2000 and TSX/TDX data in 2017). KH-9 DEM is a DEM of the study area in 1974, which was generated from three scenes of optical stereo pairs from KH-9. Geodetic glacier mass change was calculated by DH above, glacier cover vector data from TPG1976/CGI2/RGI6.0 with ice density of 850 ± 60 kg m−3. The attribute data included: GLIMSId means the glacier code from GLIMS data base, Area（km2）is the glacier area by km2, area_m2 is glacier area by (m2）, the glacier name, EC74_2000, the surface elevation change rate from 1974 to 2000(m a-1), EC00_2017, the surface elevation change rate from 2000 to 2017 (m a-1), MB74_2000, the geodetic glacier mass balance between 1974 and 2000（m w.e. a-1），MB00_2017, the geodetic glacier mass balance between 2000 and 2017（m w.e. a-1）.MC74_2000, the geodetic glacier mass change from 1974 to 2000 (m3w.e. a-1), MC00_2017, the geodetic glacier mass change from 2000 to 2017(m3 w.e. a-1). Ut_EC74_00 is the uncertainty of glacier surface elevation change（m a-1） in 1974-2000、Ut_MB74_00, is the uncertainty of glacier mass balance for each glacier（m w.e. a-1）in 1974-2000，Ut_MC74_00, is the uncertainty of glacier mass change for each glacier（m3w.e. a-1）in 1974-2000. Ut_EC00_17，is the uncertainty of glacier surface elevation change in 2000-2017（m a-1），Ut_MB00_17，is the uncertainty of glacier mass balance for each glacier in 2000-2017（m w.e. a-1），Ut_MC00_17 is the uncertainty of glacier mass change for each glacier in 2000-2017（m3 w.e. a-1）.This data set is used for the study glaciers melting and its hydrological effects in the Central Himalayas.It also could be used in studies of climatic change and disasters research in the Himalayas.
The data set involved geodetic annual glacier-averagedmass balance and mass change data at the Ponkar area in Nepal on the Southern slope of the Himalayas from 1974 to 2014. It is stored in the ESRI vector polygon format and is composed of two periods, which includes surface elevation difference between 1974-2000 (DH1974-2000, from KH-9 DEM1974 and SRTM DEM2000), surface elevation difference between 2000-2014 (DH2000-2014,by DinSAR techniques from SRTM DEM2000 and TSX/TDX data in 2014). KH-9 DEM is a DEM of the study area in 1974, which was generated from three scenes of optical stereo pairs from KH-9. Geodetic glacier mass change was calculated by DH above, glacier cover vector data from TPG1976/CGI2/RGI6.0 with ice density of 850 ± 60 kg m−3. The attribute data included: GLIMSId means the glacier code from GLIMS data base, the glacier_area（m2）、Area（km2）, the glacier name, EC74_2000, the surface elevation change rate from 1974 to 2000(m a-1), EC00_2014, the surface elevation change rate from 2000 to 2014 (m a-1), MB74_2000, the geodetic glacier mass balance between 1974 and 2000（m w.e. a-1），MB00_2014, the geodetic glacier mass balance between 2000 and 2014（m w.e. a-1）.MC74_2000, the geodetic glacier mass change from 1974 to 2000 (m3w.e. a-1), MC00_2014, the geodetic glacier mass change from 2000 to 2014(m3w.e. a-1). Ut_EC74_00 is the uncertainty of glacier surface elevation change（m a-1） in 1974-2000、Ut_MB74_00, is the uncertainty of glacier mass balance for each glacier（m w.e. a-1）in 1974-2000，Ut_MC74_00, is the uncertainty of glacier mass change for each glacier（m3w.e. a-1）in 1974-2000. Ut_EC00_14，is the uncertainty of glacier surface elevation change in 2000-2014（m a-1），Ut_MB00_14，is the uncertainty of glacier mass balance for each glacier in 2000-2014（m w.e. a-1），Ut_MC00_14 is the uncertainty of glacier mass change for each glacier in 2000-2014（m3 w.e. a-1）. This data set is used for the study glaciers melting and its hydrological effects in Ponkar area in Nepal in the Southern slope of the Himalayas. It also could be used in studies of climatic change and disasters research in the Himalayas.
The data involved two periods of geodetic glacier mass storage change of Naimona’Nyi glaciers in the western of Himalaya from 1974-2013 (unit: m w.e. a-1). It is stored in the ESRI vector polygon format. The data sets are composed of two periods of glacier surface elevation difference between 1974-2000 and 2000-2013, i.e. DHSRTM2000-DEM1974（DH2000-1974）、DHTanDEM2013-SRTM2000（DH2013-2000）. DH2000-1974 was surface elevation change between SRTM2000 and DEM1974, i.e. the earlier historical DEM (DEM1974, spatial resolution 25m) was derived from 1:50,000 topographic maps in October 1974(DEM1974,spatial resolution 25m). The uncertainty in the ice free areas of DH2000-1974 was ±0.13 m a-1. The surface elevation difference between 2000-2013 (DH2000-2013, by DinSAR techniques from SRTM DEM2000 and TSX/TDX data on Oct.17th in 2013) The uncertainty in the ice free areas of DH2013-2000 was ±0.04 m a-1. Glacier-averaged annual mass balance change (m w.e.a-1) was averaged annually for each glacier, which was calculated by DH2000-1974/DH2013-2000, glacier coverage area and ice density of 850 ± 60 kg m−3. The attribute data includes Glacier area by Shape_Area (m2), EC74_00, EC00_13, i.e. Glacier-averaged surface elevation change in 1974-2000 and 2000-2013(m a-1), MB74_00, MB00_13 i.e. Glacier-averaged annual mass balance in 1974-2000 and 2000-2013 (m w.e.a-1), and MC74_00, MC00_13, Glacier-averaged annual mass change in 1974-2000 and 2000-2013 (m3 w.e.a-1), Uncerty_MB, is the uncertainty of glacier-averaged annual mass balance（m w.e. a-1）， Uncerty_MC, is the Maximum uncertainty of glacier-averaged annual mass change（m3 w.e. a-1）. The data sets could be used for glacier change, hydrological and climate change studies in the Himalayas and High Mountain Asia.