Applicability of an ultra-long-range terrestrial laser scanner to monitor the mass balance of Muz Taw Glacier,Sawir Mountains, China

FeiTeng Wang , ChunHai Xu ,2, ZhongQin Li *, Muhammad Naveed Anjum ,2, Lin Wang

1. State Key Laboratory of Cryospheric Science/Tien Shan Glaciological Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China

2. University of Chinese Academy of Sciences, Beijing 100049, China

1　Introduction

Glacier mass-balance provides significant information in understanding climate-glacier processes and its consequences such as water resource estimation and global sea level rise (IPCC, 2013). Numerous studies have been performed on glacier area, length,thickness (volume) and mass balance by using remotely sensed observations (Li KMet al., 2011; Gardelleet al., 2013; Shangguanet al., 2015) and in-situ measurements (Zempet al., 2015). Mass-balance can be derived by using in-situ direct glaciological measurements and repeated geodetic surveys. The direct glaciological method monitors mass balance using ablation stakes and snow pits and net accumulation is measured by digging pits in the accumulated areas and stakes are drilled into ablation of the glacier surface (Cogleyet al., 2011). Although the glaciological method provides quantitative results at high temporal resolution (Zempet al., 2015), it is usually limited to individual glaciers. Direct measurements are time consuming and labor intensive. Thus, at present, only 41 glaciers have been studied for ongoing mass-balance series with more than 30 years of observations(WGMS, 2017).

Compared with traditional glaciological measurements, the geodetic method provides glacier volume changes by repeated mapping from ground, air or spaceborne surveys and subsequent differences of multi-temporal digital elevation models (DEMs) of glacier surface topography (Zempet al., 2013). The geodetic mass-balance extends the glaciological sample both in space and time, which has become a widely used method to calculate glacier mass balance and evaluate the contribution of glaciers to runoff and periodically validate and calibrate long-term glaciological mass-balance series (Cogley, 2009; Cogleyet al., 2011; Zempet al., 2013). Whereas, the limited spatiotemporal resolution of existing DEMs usually influences the calculation of accurate geodetic mass balance. A new generation of terrestrial lidar techniques has offered dense point clouds and high resolution and precision DEMs of glacier surface (Fischeret al., 2016).

In China, glacier mass-balance monitoring mainly focus on Tien Shan, Qilian Shan and the Tibetan Plateau (Li ZQet al., 2011; Wanget al., 2017; Zhuet al.,2017). Even though glaciers in Sawir Mountains have undergone an acceleration of mass loss and terminal retreat with the consequence of floods, these glaciers have received little attention at present (Huaiet al.,2015). The present study aims to provide a detailed and precise mass balance and front variation monitoring of Muz Taw Glacier using an ultra-long-range terrestrial laser scanner (TLS). We then use Muz Taw Glacier as a case study to discuss the applicability of TLS systems to monitor the mass balance of mountain glaciers. Thus, the present study provides valuable data for poorly understood mass balance and front variation of mountain glaciers.

2　Study area

The Sawir Mountains are located south of the Ertix River (Wang, 1988). According to the first Glacier Inventory of China (GIC), there were a total of 21 glaciers with a total area of 16.84 km2in Chinese Sawir Mountains (Shi, 2005). Due to accelerated shrinkage, glacier area reduced by 42.74% at a rate of 0.14 km2/a from 1959 to 2013 (Huaiet al., 2015).With the growth of numerous water resource management problems (e.g., flood protection, water supply,and operation of hydroelectric facilities) in this region, timely and correct assessment of glaciers in the Sawir Mountains is crucial to fill the gap, and it has an extremely important significance for Jeminay's economic development and people's livelihood.

Muz Taw Glacier (47.06°N, 85.56°E; glacier ID 5A259C0001) is a northeast-orientated valley glacier with an area of 4.27 km2and a length of 3.7 km according to the first GIC (Shi, 2005), located at the northern side of the Sawir Mountains (Figure 1a). The elevation from the glacier terminus to the highest point ranged from 3,090 m to 3,800 m a.s.l.. The snow line and ice volume of this glacier were about 3,310 m a.s.l. and 0.28 km3, respectively. The average ice thickness was 66 m. The mass-balance monitoring program was implemented beginning in 2014(Figure 1b). From field observations, there was no debris on the glacier surface. Muz Taw Glacier has been dramatically shrinking due to extensive melting over the past several decades.

3　Data and methods

3.1　Riegl VZ®-6000 ultra-long-range TLS

Traditional TLS systems use class 1 laser beams with wavelengths around 1,500 nm, with low reflectance over snow and ice cover, and the possible scanning distance is limited to a few hundred meters(Deemset al., 2013), which restricts the application of TLS surveys in the cryosphere (Rabatelet al., 2008).The Riegl VZ®-6000 ultra-long-range TLS operates with class 3B laser beams (wavelengths around 1,064 nm), and the system is exceptionally well-suited for measuring snowy and icy terrain. Faster surveys(up to 222,000 measurements/s) are possible even at long range (＞6,000 m) with unprecedented high accuracy and precision (RIEGL Laser Measurement System, 2014a).

The Riegl VZ®-6000 uses near-infrared laser signals to capture repeated dense point clouds of the glacier surface terrain, based on the principle of time-offlight measurement with echo digitization and online waveform processing (RIEGL Laser Measurement Systems, 2013). The laser pulse transmitter emits a beam of laser pulse, which reaches the target object and then reflected back to the laser receiver, and the pulse is used to calculate the time difference of the pulse (RIEGL Laser Measurement Systems, 2014a).

3.2　Field measurements

On June 2, 2015 and July 25, 2016, Riegl VZ®-6000 TLS measurements of Muz Taw Glacier were performed with two scanning stations (Figure 1c), coincident with same day of direct glaciological measurements. Each scanning station was fixed with a GPS-leveling point in order to prohibit ground motion. The instrument was mounted on a tripod placed on stable bedrock for each scan (Figure 2). Three-dimensional (3 D)coordinates of the two scanning stations were obtained using the real-time kinematic (RTK) global positioning system. The coordinated system of RTK surveys was the World Geodetic System 1984 (WGS84).

Figure 1　An over view of the study area. (a) Location of Muz Taw Glacier; (b) Measuring point network in 2014;and (c) S1 and S2 stations at the terminus of the glacier (image source: Google Earth)

Figure 2　TLS survey of Muz Taw Glacier in July 2016 with the Riegl VZ®-6000 terrestrial laser scanner

For TLS work, both data quality and acquisition time should be considered, thus surveying parameters are of significant importance. To avoid range ambiguity and to capture point clouds of the glacier surface to the maximum extent, coarse scanning was first performed with vertical and horizontal angles range of 60°-120° from zenith and 0°-360°, respectively, and corresponding laser pulse repetition rate was set as 50 kHz. Based on coarse scanning, laser pulse repetition rate was selected as 30 kHz for the implementation of fine scanning of the glacier surface terrain region for two scanning stations. We should note that the overlap percentage of two scans was more than 30% to satisfy the basic requirements for multi-station adjustment (MSA) (Mukupaet al., 2017). Surveying parameters are listed in Table 1.

3.3　Point cloud data processing

Here, we use RiSCAN PRO®v1.81 software to perform the point cloud data processing (RIEGL Laser Measurement Systems, 2014b), which includes:

1 Direct georeferencing. An additional sensor(e.g., navigation sensors, total station) is used to transform the laser scanner's system (Mukupaet al., 2017).In this study, GPS-RTK surveying points of each scanning station were used to transform the Scanner's Own Coordinate System (SOCS) into a Global Coordinate System (GLCS). The theory can be illustrated as (Lichtiet al., 2005)

2 Point cloud registration. The position of each scan was fixed after direct georeferencing. Whereas,due to the influence of orientation, the point clouds of the overlapped areas cannot coincide completely. The iterative closest point (ICP) algorithm was used to complete multi-station adjustment (MSA) (Besl and McKay, 1992; Zhang, 1994). The orientation of each scan was iteratively modified to calculate the best fit among the points patches based on least-squares minimization of residuals until point cloud layers of each scan coincided (Besl and McKay, 1992).

Table 1　Riegl VZ®-6000 TLS surveying parameters of Muz Taw

3 Point cloud vacuation and filtering. We used an octree algorithm to achieve point cloud data vacuation. In particular, the algorithm was applied to the registered scans to generate points with equal spacing(Schnabel and Klein, 2006; Perroyet al., 2010). Then a terrain filter was used to remove noise and nonground points (RIEGL Laser Measurement Systems,2014b).

In the support of ArcGIS10.1 software, TLS-derived DEMs were generated using natural neighbor interpolation with a resampled resolution of 1 m.

3.4　Glaciological mass balance measurements

The glaciological mass balance measurements of Muz Taw Glacier were performed in 2014. About 20 ablation stakes were drilled into the glacier surface evenly at different elevation bands using a stream drill(Figure 1c). Generally, a spatial distribution of singlepoint ablation or accumulation, snow density (if there is snow cover) were determined by stakes and snow pits at intervals of a month from the beginning of May to early September each year. Therefore, the measured programs include the stake vertical height over the glacier surface, thickness of superimposed ice,thickness and density of each snow/firn layer at individual stakes. By using in-situ measured densities of snow and ice, stake and snow pit changing values were converted into mass (water equivalent) to calculate single point mass balance (the summation of snow/firn, superimposed ice and glacier ice mass balance). Then point observations can be extrapolated to glacier-wide mass balance contour-line and isoline methods (Xie and Liu, 2010).

3.5　Geodetic mass balance

3.5.1Calculation of geodetic mass balance

The geodetic mass balance was calculated by the subtraction of two high-resolution DEMs. The volume change ΔVwas determined by summing the elevation change Δhiat the individual pixelrof different time periods (Zempet al., 2013).

whereNis the number of pixels covering the glacier at the maximum extent, andris the pixel size (1m×1m).

To estimate glacier mass balance (Bgeo), observed geodetic volume change must be converted to specific mass balance (m w.e.) by using the following Equations(4)and(5):

whereρpresents the average density of volume changes of 850±60 kg/m3as recommended by Huss(2013),Sis the mean glacier area over the PoR, the two survey datest1andt2, assuming a linear change(Thibertet al., 2008; Zempet al., 2013).

3.5.2Uncertainty analysis

Uncertainties in the geodetic mass balance include systematic or random errors (Zempet al.,2013), those include: (1) point cloud data acquisition errors, including terrain and atmospheric environment (wind and moisture); (2) data processing errors and DEM creation,i.e., registration (multi-station adjustment), point cloud filtering and interpolation(smoothing terrain information) (Gabbudet al., 2015;Hartzellet al., 2015; Fischeret al., 2016); and (3)spatial bias of the two DEMs and sampling density of snow and ice (Nuth and kääb, 2011; Zempet al.,2013; Thomsonet al., 2017).

Point clouds of the glacier surface should be captured in sunny, windless weather to remove the influence of wind and moisture. The octree algorithm built the topological relationship of scattered points to reduce the point cloud, in order to preserve terrain information as much as possible. To eliminate the spatial bias of the two DEMs, scanning point clouds of 2015 were set as the reference; a relative registration selects the ground control points of the non-glaciated stable terrain to match the point clouds of the two periods.

Uncertainties of the average TLS-derived glacier surface elevation changes for individual glaciers were quantified based on geostatistical analysis methods(Rolstadet al., 2009), which can be approximated by

whereSis the glacierized area of the subtracted DEMs, which is same as theSdefined by Equation(5)at the annual scale,Scoris the area over which there are spatial correction errors in the subtraction of the two DEMs, and we conservatively assumeScor=S(Fischeret al., 2016). TheσΔhTLSis the standard deviation for TLS-derived bedrock (non-glaciated ground) elevation changes.

The uncertaintyσBTLSin the homogenized mass balance results from errors in the geodetic mass change (Husset al., 2009):

whereσρis the uncertainty of density.is the mean of TLS-derived geodetic elevation changes and the related uncertainty relies on the accuracy of the DEMs. Derived values are given in Table 2.

Table 2　Density (ρ) and its uncertainty (σρ) according to Huss (2013), TLS-derived glacier surface elevation changes (ΔhTLS),standard deviation of surface elevation changes (σΔhTLS) in stable bedrock area, and uncertainty of mass balance (σBTLS)

4　Results and discussion

4.1　Glacier front variation (length changes)

Glacier front variation data are a key element for evaluating the regional representativeness of the few glaciological measurements at spatiotemporal scales(Zempet al., 2015). Using Landsat TM/ETM+ images, previous studies show that the glacier terminus had an observed retreat of 269.19 m from 1977 to 2013, at a retreat rate of 7.28 m/a (Huaiet al., 2015).However, the low spatial resolution images usually limited the accuracy of glacier front variations. As in-situ measurements of glacier terminus were poor,in this work, based on the visual interpretation of the TLS-derived high-resolution point clouds, we extracted six stable landmarks and calculated their corresponding distance to Muz Taw Glacier terminus, and those landmarks distributed evenly at the front of the glacier (Figure 3). Observed glacier front variations ranged from 8.4 to 10.6 m, with an average value of 9.3 m from June 2, 2015 to July 25, 2016 (Table 3),and the annual retreat rate was 8.1 m, which demonstrated an accelerated terminus retreat rate in recent years.

Figure 3　Landmarks and their distance to glacier terminus were extracted based on the visual interpretation of the TLS-derived high-resolution point clouds

4.2　Glacier surface elevation changes and geodetic mass balance

The map of TLS-derived geodetic elevation changes revealed a serious ice loss at the ablation area(Figure 4) and the glacier tongue showed more negative ice loss than the upper part. The mean elevation change was -4.55 m for the ablation area between June 2, 2015 and July 25, 2016. Elevation changes measured along the main flow line allow more detailed insights into the characteristics of glacier behavior. The elevation differences showed well changing rule, that was to say, the ice thinning value slackened when the length increased (the distance from glacier terminus) (Figure 5), which agrees well with field observations. Figure 4 illustrates the distribution of individual ablation stakes (black crosses) and TLS-derived geodetic elevation changes(black digitals) were calculated at the corresponding stakes, which closely match direct glaciological elevation changes (stake height over glacier surface)from in-situ measurements. Considering the uncertainties of the two measuring methods, glacier ice movement and the influence of terrain factors, differences in surface elevation changes of the point measurements were neglected for this work, indicating that TLS-derived glacier surface elevation changes are reliable.

Table 3　Statistical front variations of Muz Taw Glacier

Figure 4　Spatial distribution of TLS-derived annual surface elevation changes at the ablation area of the Muz Taw Glacier.Black crosses and letters present the location of ablation stakes,red and black digitals in the boxes show the elevation changes of corresponding stakes (points) derived by the glaciological methods and TLS, respectively

For the conversion of observed glacier elevation changes into the mass budget a density of ice/firn/snow needs to be assumed. Here, a recommended density assumption of 850±60 kg/m3is used and the density was provided by Huss (2013) using an empirical firn densification model with idealized surface mass balance forcing. Using Equations(3)-(5), the mass balance of the ablation area was determined, and we selected 375 stable points in the non-glaciated ground to evaluate the uncertainty of TLS-derived elevation changes and geodetic mass balance (Equations(6),(7)). The mass balance uncertainty of the observed period was ±0.378 m w.e.. and the calculated balance was -3.864±0.378 m w.e.. Compared to the mass balance, the uncertainty was very small,which indicates a high accuracy of TLS-derived DEMs.

Figure 5　Longitudinal profile of glacier surface elevation changes for the observed period

4.3　Applicability and recommendation of the TLS system to monitor the mass balance

Due to high temporal-spatial resolution and ultralong-range of Riegl VZ®-6000, previous studies have shown that this technique can provide accurate mass balance data of very small glaciers in the Swiss Alps(Fischeret al., 2016) and surface melting values of Alpine glaciers at the seasonal and hourly scales(Gabbudet al., 2015). What's more, TLS can work without contact with objects and in a rugged environ-ment, thus allowing measurements of dangerous or inaccessible areas of the glacier surface. Although we had inserted more than 20 ablation stakes on the Muz Taw Glacier surface, they were distributed evenly at different elevation bands. Some areas, such as the terrain of firn basin is precipitous, and the surface topography of the tributary glacier (right branch) is also dangerous (Figure 1c). Therefore, it is not easy for us to insert ablation stakes and dig snow pits in those inaccessible areas. However, glacial terrain and insufficient spatial distribution of stakes have significant influence on the calculation of glacier-wide mass balance. TLS technology makes it possible for us to monitor those inaccessible regions, so we recommend the application of Riegl VZ®-6000 TLS for future mass balance monitoring of Muz Taw Glacier. To better monitor the mass balance, several recommendations need to be followed to improve the results of glacier measurements: (1) good weather (dry and windless atmosphere) is of great importance for field surveys; (2) stable and visual angles of the scanning stations are necessary between surveys, which could minimize registration errors among the scans of the same time; (3) adequate and visible areas of non-glaciated stable terrain are required to guarantee a good quality of the relative registration of different temporal point clouds; (4) TLS measurements should be implemented at the end of ablation season to decrease the influence of fresh snow cover; (5) the overlap percentage of each scanning station of the same period should be more than 30% to achieve an ideal mosaic of multi-station adjustment; and (6) monitoring of glacier density and glaciological mass balance at Muz Taw Glacier should continue. Long-term mass-balance series measurements are needed, especially for important regions where there is a shortage of in-situ monitoring data. The TLS-derived geodetic results can validate the reliability of glaciological balance and complement the glaciological method.

5　Conclusion and outlook

The Riegl VZ®-6000 ultra-long-range TLS system is well suited for surveying snow- and ice-covered terrain and has a ranging accuracy better than 10 mm.The system can capture surface elevation changes and the mass of mountain glaciers at high spatiotemporal resolution, so the quality of TLS-derived geodetic net mass balance is very good. In the present paper, we have shown the applicability of an ultra-long-range terrestrial laser scanner to monitor mass balance of Muz Taw Glacier.

The TLS-derived high spatial resolution point clouds determined the front variation (glacier retreat)of Muz Taw Glacier of 9.3 m from June 2, 2015 to July 25, 2016. The geodetic elevation change was 4.55 m at the ablation area and the calculated balance was-3.864±0.378 m w.e.. Compared with in-situ measurements, the glaciological elevation change of individual stakes and the TLS-derived geodetic elevation change of corresponding points matched closely.Mean uncertainty in the TLS-derived mass balance was ±0.378 m w.e.. The uncertainty is attributed to data acquisition errors, data-processing and DEM creation errors. Thus, in further studies, we should take into account recommendations mentioned in Section 4.3 to generate high quality data and minimize the uncertainty of TLS-derived results.

This research was supported by the National Natural Science Foundation of China (41601076, 41471058 and 91425303), the "Light of West China" program for Talent Introduction of Chinese Academy. We would like to thank the Tien Shan Glaciological Station for field surveys.

Besl PJ, McKay ND, 1992. A method for registration of 3-D shapes.IEEE Transactions on Pattern Analysis and Machine Intelligence,14(2): 239-256. DOI: 10.1109/34.121791.

Cogley JG, 2009. Geodetic and direct mass-balance measurements:comparison and joint analysis. Annals of Glaciology, 50(50):96-100. DOI: 10.3189/172756409787769744.

Cogley JG, Hock R, Rasmussen LA,et al., 2011. Glossary of Glacier Mass Balance and Related Terms. IHP-VII Technical Documents in Hydrology No. 86, IACS Contribution No. 2. Paris:UNESCO-IHP.

Deems JS, Painter TH, Finnegan DC, 2013. Lidar measurement of snow depth: a review. Journal of Glaciology, 59(215): 467-479.DOI: 10.3189/2013JoG12J154.

Fischer M, Huss M, Kummert M,et al., 2016. Application and validation of long-range terrestrial laser scanning to monitor the mass balance of very small glaciers in the Swiss Alps. The Cryosphere,10(3): 1279-1295. DOI: 10.5194/tc-10-1279-2016.

Gabbud C, Micheletti N, Lane SN, 2015. Lidar measurement of surface melt for a temperate Alpine glacier at the seasonal and hourly scales. Journal of Glaciology, 61(299): 963-974. DOI:10.3189/2015JoG14J226.

Gardelle J, Berthier E, Arnaud Y,et al., 2013. Region-wide glacier mass balances over the Pamir-Karakoram-Himalaya during 1999-2011. The Cryosphere, 7(4): 1263-1286. DOI: 10.5194/tc-7-1263-2013.

Hartzell PJ, Gadomski PJ, Glennie CL,et al., 2015. Rigorous error propagation for terrestrial laser scanning with application to snow volume uncertainty. Journal of Glaciology, 61(230): 1147-1158.DOI: 10.3189/2015JoG15J031.

Huai BJ, Li ZQ, Wang FT,et al., 2015. Variation of glaciers in the Sawuer Mountain within Chinese territory during 1959～2013.Journal of Glaciology and Geocryology, 37(5): 1141-1149. DOI:10.7522/j.isnn.1000-0240.2015.0128. (in Chinese)

Huss M, Bauder A, Funk M, 2009. Homogenization of long-term mass-balance time series. Annals of Glaciology, 50(50): 198-206.DOI: 10.3189/172756409787769627.

Huss M, 2013. Density assumptions for converting geodetic glacier volume change to mass change. The Cryosphere, 7(3): 877-887.DOI: 10.5194/tc-7-877-2013.

Intergovernmental Panel on Climate Change (IPCC), 2013. Climate Change 2013: the Physical Science Basis. Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press.

Li KM, Li ZQ, Gao WY,et al., 2011. Recent glacial retreat and its effect on water resources in eastern Xinjiang. Chinese Science Bulletin, 56(33): 3596-3604. DOI: 10.1007/s11434-011-4720-8.

Li ZQ, Li HL, Chen YN, 2011. Mechanisms and simulation of accelerated shrinkage of continental glaciers: a case study of Urumqi Glacier No. 1 in Eastern Tianshan, Central Asia. Journal of Earth Science, 22(4): 423-430. DOI: 10.1007/s12583-011-0194-5.

Lichti DD, Gordon SJ, Tipdecho T, 2005. Error models and propagation in directly georeferenced terrestrial laser scanner networks.Journal of Surveying Engineering, 131(4): 135-142. DOI:10.1061/(ASCE)0733-9453(2005)131:4(135).

Mukupa W, Roberts GW, Hancock CM,et al., 2017. A review of the use of terrestrial laser scanning application for change detection and deformation monitoring of structures. Survey Review, 49(353):99-116. DOI: 10.1080/00396265.2015.1133039.

Nuth C, Kääb A, 2011. Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change. The Cryosphere, 5(1): 271-290. DOI: 10.5194/tc-5-271-2011.

Perroy RL, Bookhagen B, Asner GP,et al., 2010. Comparison of gully erosion estimates using airborne and ground-based LiDAR on Santa Cruz Island, California. Geomorphology, 118(3-4):288-300. DOI: 10.1016/j.geomorph.2010.01.009.

Rabatel A, Deline P, Jaillet S,et al., 2008. Rock falls in high-alpine rock walls quantified by terrestrial LiDAR measurements: a case study in the Mont Blanc area. Geophysical Research Letters,35(10): L10502. DOI: 10.1029/2008GL033424.

RIEGL Laser Measurement Systems, 2013. Preliminary Data Sheet,07.05.2013; RIEGL VZ®-6000 - 3D Ultra long range terrestrial laser scanner with online waveform processing. Horn, Austria:RIEGL Laser Measurement Systems.

RIEGL Laser Measurement Systems, 2014a. 3D terrestrial laser scanner Riegl VZ®-4000 / Riegl VZ®-6000 General Description and Data Interfaces. Horn, Austria: RIEGL Laser Measurement Systems.

RIEGL Laser Measurement Systems, 2014b. RiSCAN PRO®- Version 1.8.1. Horn, Austria: Riegl Laser Measurement Systems.

Rolstad C, Haug T, Denby B, 2009. Spatially integrated geodetic glacier mass balance and its uncertainty based on geostatistical analysis:application to the western Svartisen ice cap, Norway. Journal of Glaciology, 55(192): 666-680. DOI: 10.3189/002214309789470950.

Schnabel R, Klein R, 2006. Octree-based point-cloud compression. In:Proceedings of the 3rd Eurographics/IEEE VGTC Conference On Point-based Graphics. Boston, MA, USA: ACM, 111-120. DOI:10.2312/SPBG/SPBG06/111-120.

Shangguan DH, Bolch T, Ding YJ,et al., 2015. Mass changes of Southern and Northern Inylchek Glacier, Central Tian Shan,Kyrgyzstan, during ～1975 and 2007 derived from remote sensing data. The Cryosphere, 9(2): 703-717. DOI: 10.5194/tc-9-703-2015.

Shi YF, 2005. Concise China Glacier Inventory. Shanghai: Shanghai Popular Science Press, pp. 101-105.

Thibert E, Blanc R, Vincent C,et al., 2008. Glaciological and volumetric mass-balance measurements: error analysis over 51 years for Glacier de Sarennes, French Alps. Journal of Glaciology, 54(186):522-532. DOI: 10.3189/002214308785837093.

Thomson LI, Zemp M, Copland L,et al., 2017. Comparison of geodetic and glaciological mass budgets for White Glacier, Axel Heiberg Island, Canada. Journal of Glaciology, 63(237): 55-66. DOI:10.1017/jog.2016.112.

Wang S, Yao TD, Tian LD,et al., 2017. Glacier mass variation and its effect on surface runoff in the Beida River catchment during 1957-2013. Journal of Glaciology, 63(239): 523-534. DOI:10.1017/jog.2017.13.

Wang ZT, 1988. New statistical figures and distribution feature of glaciers on the various Mountains in China. Arid Land Geography,11(3): 8-14. DOI: 10.13826/j.cnki.cn65-1103/x.1988.03.002.

World Glacier Monitoring Service (WGMS), 2017. Global Glacier Change Bulletin No. 2 (2014-2015). In: Zemp M, Nussbaumer SU,Gärtner-Roer I,et al. (eds.), ICSU(WDS)/IUGG(IACS)/UNEP/UNESCO/WMO, World Glacier Monitoring Service,Zurich, pp.1-7. DOI: 10.5904/wgms-fog-2017-10.

Xie ZC, Liu CH, 2010. Introduction to Glaciology. Shanghai Popular Science Press, Shanghai, pp. 1-490.

Zemp M, Thibert E, Huss M,et al., 2013. Reanalysing glacier mass balance measurement series. The Cryosphere, 7(4): 1227-1245.DOI: 10.5194/tc-7-1227-2013.

Zemp M, Frey H, Gärtner-Roer I,et al., 2015. Historically unprecedented global glacier decline in the early 21st century. Journal of Glaciology, 61(228): 745-762. DOI: 10.3189/2015JoG15J017.

Zhang ZY, 1994. Iterative point matching for registration of free-form curves and surfaces. International Journal of Computer Vision,13(2): 119-152. DOI: 10.1007/BF01427149.

Zhu ML, Yao TD, Yang W,et al., 2017. Differences in mass balance behavior for three glaciers from different climatic regions on the Tibetan Plateau. Climate Dynamics. DOI: 10.1007/s00382-017-3817-4. (in Press)