9 December 2014 Rapid glacier retreat in the Naimona’Nyi region, western Himalayas, between 2003 and 2013
Author Affiliations +
J. of Applied Remote Sensing, 8(1), 083508 (2014). doi:10.1117/1.JRS.8.083508
The Northern Hemisphere experienced its warmest temperatures in the early 21st century, and this clearly increased glacier melting on the Tibetan Plateau. This study analyzed the glacier change in the area, their terminus positions and surface elevations over the last decade (2003 to 2013) in the Naimona’Nyi region of the western Himalayas, by comparing remote sensing data and differential global position system (dGPS) data from in situ surveys. The results show an accelerating glacier retreat over the past decade in this region. The area covered by glaciers was reduced by a total of 13.2±0.0022  km2 over this period. The terminus of the Naimona’Nyi main glacier studied in the paper retreated by over 191±35  m. We also compared the ice cloud and elevation satellite elevation data with in situ measured decimeter accuracy dGPS elevation data, thus providing changes in glacier surface elevations in different periods. The maximum measured glacier thinning rate reached to 0.58±0.06  m/a between 2009 and 2013. In addition to the dependence on elevation and glacier size, we found a larger retreat on the north facing slopes than on the south facing slopes in the region. Meteorological data show that glacier changes within the study area can probably be attributed to the observed rapid temperature rise.
Zhu, Tian, Wang, Wang, and Cui: Rapid glacier retreat in the Naimona’Nyi region, western Himalayas, between 2003 and 2013



Glacier change is the best indicator of climate change, and it rapidly responds to changes in precipitation and temperature.1,2 Glacier shrinkage in the high Himalayas and the Tibetan Plateau has drawn worldwide attention to the changing climate. In addition, glacial melt water provides an important freshwater resource and significantly contributes to global sea-level rise.34.5.6

The Himalayas hold more than one-fourth of the glacier area in all of “High Asia” and most of the glaciers have, in recent decades, been rapidly retreating and losing ice.78.9 Great efforts have been made to assess the shrinkage of glaciers in the eastern and central Himalayas of China.1011.12.13.14 However, less research has been done in the western Himalaya due to the problem of difficult access. Studies on glacier change will provide valuable information related to climate change both regionally and globally.15 In addition, glaciers within our study area are important water resources for production and local habitants.16

Earlier studies on the variability of the Naimona’Nyi glacier have been carried out by a few researchers. Yao et al. estimated the glacier mass balance from stake measurements,16 while Wang et al. and Ye et al. calculated glacier areas and volume losses before 2003 based on remote sensing data.17,18 Kehrwald et al. analyzed the ice core data from the Naimona’Nyi main glacier and found that the surface of this glacier had melted due to the warming climate.3 Since the results of those researches are earlier works, glacier changes in Naimona’Nyi area should be further monitored because the past decade was the warmest decade of the Northern Hemisphere19 and more new data is available.

A variety of methods has been applied to monitor glacial changes. The most commonly used methods compare data from a variety of images (e.g., remote sensing, digital elevation models, and topographic maps) in different periods.9,12 Compared with glacier area change, glacier surface elevation for the most recent changes is more difficult to investigate. The traditional way to measure glacier thinning is direct annual measurements. However, annually, in situ measurements are difficult because of the harsh working conditions. Laser altimetry data acquired by the geoscience laser altimeter system (GLAS) carried on the ice cloud and elevation satellite (ICESat) have now been proven to provide accurate data for glacier surface elevations.20,21 However, the satellite tracks in different years do not overlap and do not allow for easy annual comparisons for these data. Combining the ICESat data with in situ decimeter-accuracy differential global position system (dGPS) measurements of the glacier may provide a better approach for monitoring changes in glacier thickness with a high accuracy. Fortunately, the Naimona’Nyi main glacier in the research region is covered by ICESat data and is now available for this kind of comparison.

This study discusses how these glaciers have changed over the last decade in the Naimona’Nyi region of the western Himalaya, based on both ICESat data and site geodetic measurements, including the glacier area, the glacier’s terminal position, and specific glacier surface elevation changes, in order to determine their responses to climate change.


Study Area

Mt. Naimona’Nyi lies in the western Himalayas, southwest of the Tibetan Plateau (Fig. 1). The peak reaches an elevation of 7694 m, the highest mountain in the western Himalayas. A large number of glaciers were formed here because of the high elevation. Earlier work shows that the equilibrium line altitude (ELA) ranges between 5100 and 6000 m, with the ELA on the northern slopes 400 m higher than that on the southern slopes.22 The Naimona’Nyi main glacier, the largest glacier in the research region, extends over 10 km and has a width of approximately 3 km. The maximum thickness of this main glacier exceeds 250 m; it was measured in 2008 by ground penetrating radar and the result has not been published anywhere. Its total area exceeds 10km2.

Both the westerly circulation and the southwest monsoons affect the Naimona’Nyi region. The nearest meteorological station (Purang: 30°17′N, 81°15′E; the elevation of 3900 m), which is only 20 km away from the Naimona’Nyi main glacier, shows an annual average precipitation of 169 mm and average air temperature of 3°C. However, precipitation could be much greater in the high mountain regions, perhaps between 760 and 1000 mm.16


Data and Methods


Optical Remote Sensing Data


Data preprocessing and glacier delineation

This study used four Landsat enhanced Thematic Mapper Plus (ETM+) images and one from the Landsat operational land image (OLI) provided by the United States Geological Survey23 to map the glaciers. Except for imaging time, these images have the same path/row (144/39) and spatial resolution (30 m). In order to delineate the glaciers and reduce the influence of seasonal snow and cloud cover, we selected images in late summer and early autumn. Because no image from this season was available in 2013, we used images from December 2013 instead. Because there was no heavy snow prior to December according to meteorological data in this region and light snow could not last, the seasonal snow influence is weak and was ignored here. The spatial resolutions for both Landsat ETM+ and Landsat OLI optical remote sensing images have been determined to be 30 and 15 m for the panchromatic band. Because of a sensor problem, a gas stripe appears in the Landsat ETM+ images after May 31, 2003. In order to ensure the data accuracy, we repaired the data, registered the new data with the high quality Landsat OLI image, and constrained the error within one pixel. Global digital elevation model data from National Aeronautics and Space Administration have a resolution of 30 m24 (The Land Processes Distributed Active Archive Center). All image projections use the WGS84-UTM coordinates.

Image fusion helps to improve image resolution and the identification of terrestrial objects,25,26 and is commonly used for monitoring glaciers by optical remote sensing.27 The selection of fusion methods depends on the features discernible in the fused images and the purpose of the fusion. Studies have shown that the principal component analysis transform algorithm for Landsat ETM+ panchromatic and multispecial images provides the optimal fusion,28 which has been used in this work. Our study fused the false color images synthesized from the ETM+ 543 band and the panchromatic images to produce multispectral images with a resolution of 15 m. The fused images were used to extract glacier boundaries.

A variety of methods have been used to determine glacier boundaries from optical remote sensing images, including: supervised classification, unsupervised classification, Normalized Difference Snow Index, Band Ratio, visual interpretation, and so on. Among them, visual interpretation has the highest accuracy, especially for discrimination between snow accumulation and debris areas with a low probability of error. Band ratio is also a suitable method in practice.29 Considering the existence of the seasonal snow, shadow areas and small clouds, we did not use these automatic identification methods, which could affect the accuracy. Instead, we primarily relied on visual interpretation and used the TM4/TM5 band ratio as complementary to delineate the glaciers in the Naimona’Nyi region. A threshold was set to 2.3 to separate glaciers from nonglacier covered areas.


Uncertainty in delineating glacier outlines

Previous experience suggests that spatial resolution and registration errors provide two key causes of uncertainty. The uncertainty can be estimated by Eq. (1):30,31


where UT represents the uncertainty of the glacier terminus, λ is the original spatial resolution of a single image, and ε is a single image registration uncertainty.

The various sources of uncertainty are complex. However, Eq. (1) provides a simple and an effective way to estimate glacier terminus uncertainty in practice and has been used in this work. Results indicated that the glacier terminus measurement uncertainties are 26, 25, and 24 m, respectively, for 2003 to 2006, 2006 to 2009, and 2009 to 2013. The total UT among the four images was 35 m during 2003 to 2013.

The errors in glacial terminus result in uncertainties in the delineated glacier areas. Ye et al. obtained the glacier area uncertainty (UA) by using Eq. (2):31



Therefore, we estimated the uncertainty of the glacier area in the Naimona’Nyi region to be 0.0022km2.


Surface Elevation Data


Ice cloud and elevation satellite and high precision differential GPS measurements

The GLAS on ICESat can monitor the elevation change of a glacier surface, especially on the polar ice sheet.32 GLAS offers 15 products for the entire Earth, which can be downloaded from the International Snow and Ice Data Center (NSIDC).33 This paper used the latest “release 33-level 2” of GLAS14 products, which have the same parameters except for the dates. Compared with earlier versions, “release 33” incorporates a variety of optimizations and improvements on the fast line-of-sight atmospheric analysis of spectral hypercubes, parameter corrections, accuracy improvements, elevation treatments, and so on.

NSIDC gave the central geographic coordinates of each footprint (footprint: diameter 70 m), but did not include the center point’s elevation. In this paper, we regarded the footprint’s elevation as the center point’s elevation, and calibrated the elevation through verification points, so that we could match the field dGPS survey and the footprint of ICESat. This approach basically compares the elevation changes from ICESat data and dGPS data for specific footprints on the glacier. First, we extract the center geographic coordinates and elevation for each ICEsat footprint. Subsequently, we determine the center point position by dGPS, remeasure the current elevation, and subtract the ICEsat elevation from the remeasured dGPS elevation. Thus, the glacier surface elevation change during a specific period can be determined.

ICESat provides the elevation data for the period from 2003 to 2009. Our study obtained the surface elevations of the Naimona’Nyi glacier in 2013 by dGPS remeasurements for each specific footprint. We made the dGPS measurement on September 12, 2013, using the NavCom Starfire3050 differential GPS receiver from Unistrong Company. We tested the measurement precision by comparing the dGPS data at a National first-class geographic reference point in Tibet and the vertical error was within ±0.06m. The GLAS14 elevation data are based on TOPEX/Poseidon as a reference datum, while the SF3050 satellite differential measurement was based on WGS-84. Therefore, we combined the two kinds of reference data into WGS-84.


Data qualification for ice cloud and elevation satellite

Although the footprint of ICESat proved to be an accurate data source for glacier surface elevations, the elevation is the average value of the footprint (diameter: 70 m) not the center point of the footprint. As a result, the elevation provided in ICESat may be doubtful on rough surface. The part of the glacier surface where we measured on the Naimona’Nyi main glacier is large and quite flat, with a slope lower than 6 deg. This will constrain the uncertainty in elevation, especially within the scope of a footprint. However, the calibration is still necessary. Previous studies have used decimeter accuracy dGPS differential GPS (dGPS) to verify the difference of ICESat data,34,35 with results showing that the deviation is closely related to spatial location. We examined these differences in the Naimona’Nyi region using an in situ dGPS survey for specific footprints. We selected 18 footprints to estimate the difference between the center point’s elevation and footprint data in the flat nonglacier terrain near the Naimona’Nyi glacier. The result shows that the maximum difference, the minimum difference, and average difference in this area are 3.6, 1.3, and 2.2 m, respectively (Fig. 2).

Fig. 1

Location of the Naimona’Nyi glaciers in the western Himalayas. The inner rectangle shows the glaciers studied in this work.


Fig. 2

The differences between footprint elevation and center point elevation of 18 verification points.


We then used the average difference value to calibrate ICESat elevations so as to remove the uncertainty of the elevation values for the glacier surface.




Changes in Glacier Area

We estimated the changes in glacier area in the Naimona’Nyi region from satellite images. Our results show that the total glacier coverage in the region reached 79.6km2 in 2003, 76.3km2 in 2006, 72.3km2 in 2009, and 66.3km2 in 2013 during the three intervening intervals (Table 1). This is equivalent of annual area change rates of 1.3%, 1.7%, and 2.1%, respectively. An earlier work using advanced spaceborne thermal emission and reflection radiometer images showed that glaciers in this region covered 79.3km2 in 2003,18 matching well with our results. Slight differences might come from data sources, interpretation methods and timing. Results show that glacier area in 2013 only covered 83% of the glacier area in 2003, or a shrinkage of 13.2km2 over the last decade. Glacier retreating rates varied in different periods, with the most rapid retreat during 2009 to 2013 (Table 1).

Table 1

Glacier area changes in the Naimona’Nyi region between 2003 and 2013.

TimeArea(km2)Area change (km2)Annual change rate (km2/a)Annual change rate (%)Glacier terminus change (m)Annual change (m/a)
Total glacier area2003 to 200679.5 to 76.3−3.2±0.0015−1.1−1.3
2006 to 200976.3 to 72.3−4.0±0.0015−1.3−1.7
2009 to 201372.3 to 66.3−6.0±0.0015−1.5−2.1
2003 to 201379.5 to 66.3−13.2±0.0022−1.3−1.6
Main glacier2003 to 200615.4 to 14.6−0.8±0.0015−0.3−1.8−82−27.3
2006 to 200914.6 to 14.90.3±0.00150.10.6−58−19.3
2009 to 201314.9 to 14.6−0.3±0.0015−0.1−0.5−51−12.8
2003 to 201315.4 to 14.6−0.8±0.0022−0.1−0.5−191−19.1

Glacier area changes between the northern and the southern slopes in the Naimona’Nyi area are also compared [Fig. 3(a)]. The glaciers on the southern slopes have lost 3.4km2 in area over 10 years, which is only 35% of that on the northern slopes. Also, the temporal patterns of retreat are different between the two slopes. The glacier coverage of the southern slopes changed by 0.2km2, 3.9km2, and 0.3km2 during the periods of 2003 to 2006, 2006 to 2009, and 2009 to 2013, respectively. During the same periods, the glacier area on the northern slopes changed by 3.4km2, 0.1km2, and 6.2km2. Compared with the consistent retreats on the northern slopes, the southern slopes showed only a moderate retreat or even slight advances.

Fig. 3

Changes in glacier areas in Naimona’ Nyi region with (a) different slopes and (b) elevation.


Approximately 89% of the glacier areas’ loss during the period between 2003 and 2013 is within the elevation zone between 5500 and 6500 m. However, on the lower zones of the glaciers, up to 96% of the glacier area below 5500 m disappeared. Between 5500 and 6000 m, the glacier area lost 23%, and between 6000 and 6500 m the glacier area lost 11%. Above 6500 m, no obvious glacier area change is found [Fig. 3(b)].


Change in the Area and Terminal Position of the Naimona’Nyi Main Glacier

The Naimona’Nyi main glacier is the largest single glacier in this region and accounts for 20% of the total glacier area. This study addresses changes of this main glacier, which lost 5.3% of its area between 2003 and 2013. Despite the overall retreating trend, our study found this glacier had expanded 0.3km2 during the period between 2006 and 2009. The most rapid glacier retreat took place between 2003 and 2006, with a retreat rate of 0.3km2/a, three times larger than the average between 2003 and 2013.

The glacier’s terminal position may respond differently to climate change than the glacier area change. The glacier’s terminal position retreated by 82, 58, and 51 m during the intervals of 2003 to 2006, 2006 to 2009, and 2009 to 2013. During the past 10 years, the main glacier’s terminus retreated 191 m, with an annual average rate of 19.1m/a. The largest retreat of the glacier terminus occurred between 2003 and 2006, in accordance with the glacier area change. Unlike the zero changes in glacier area between 2006 and 2013, the glacier terminus still showed a considerable rate of retreat (Fig. 4).

Fig. 4

Changes of the main glacier’s outline and glacier terminal position (a) between 2003 and 2013. (Color dots refer to ICESat footprints in different times).



Surface Elevation Changes of the Naimona’Nyi Main Glacier

The coverage of the ICESat orbit only allows the study of glacier surface elevation changes on specific parts of the Naimona’Nyi main glacier. A total of 74 ICESat footprints from different years (Table 2) were remeasured by dGPS on the glacier surface in 2013. Then we compared the remeasured dGPS data with ICESat footprint elevations as described in the early section to obtain changes of the glacier surface elevations. However, because the data came from different periods, we present the annual change rates during specific periods (Table 2). For the ICEsat points in 2003, the average measured glacier elevation change was 2.7m on average between 2003 and 2013, with an annual rate of 0.27m/a. The largest rate of elevation change was 0.58m/a during the period of 2009 to 2013.

Table 2

Average changes of the surface elevation on the main glacier.

Time periodAverage change (m)Annual change rate (m/a)Number of footprints
2003 to 2013−2.7−0.27±0.0610
2005 to 2013−3.6−0.45±0.0611
2006 to 2013−2.6−0.37±0.0623
2007 to 2013−2.3−0.39±0.0611
2009 to 2013−2.3−0.58±0.0619
Total: 74

We found a slight changing rate of glacier surface elevation with different altitude zones. We divided the glacier surface into seven elevation zones based on an elevation intervals of 50 m, and then compared the glacier surface elevation changes in each group (Fig. 5). The maximum glacier elevation change rate of 0.53m/a was within the elevation of 5750 to 5800 during the period of 2009 to 2013. A slight decrease in the elevation change rate was noted between the lower and the upper portions of the glacier, as shown in the black dashed line in Fig. 5.

Fig. 5

Surface elevation changes in different elevations of the main glacier.





Glacier Retreat

Since the early 20th century, glaciers on the Tibetan Plateau have generally been retreating, and this trend has accelerated since the 1990s.36 Our work on the Naimona’Nyi glaciers in the western Himalayan shows that glacier retreat increased during the last decade. The glacier retreat rate between 2003 and 2013 is five times of that in the previous 30 years (1973 to 2003).18 The observed annual retreat rate of 1.6%/a for the Naimona’Nyi glaciers significantly exceeds the 0.42%/a for the entire western Himalayas, 0.41%/a for the middle Himalayas, and 0.57%/a for the eastern Himalayas (1980–2010).15 Glacier retreat in the Himalayas appears to be accelerating over the past decade, which matches conclusions from other regions.37,38 This work also found that the main glacier in the Naimona’Nyi region remains relatively inert compared with the majority of other glaciers in this region. The retreat rate for the Naimona’Nyi main glacier remains lower than most glaciers on the Tibetan Plateau, but seems to be roughly similar to those in the West Kunlun region.39 This may be related to the fact that small glaciers are more sensitive to climate change, and small glaciers account for much of the glacier area loss.

Glacier elevation and glacier orientation affect their rates of glacier area change. This has been found in the middle of Tienshan and in western Nyainqentanglha.40,41 In the Naimona’Nyi region, the lower portions of the glaciers show the fastest retreats, and glaciers below 5500 m have almost entirely disappeared. The maximum melt zone in the Naimona’Nyi region lies in the 5500 to 6000 m elevation range, which is similar to that in the Mt.Qomolangma National Nature Reserve.42 More glaciers in the Naimona’Nyi region lie on the northern slopes, including the Naimona’Nyi main glacier. Roughly, the glaciers on the northern slope should retreat less than those on the southern slope which receives more direct sunlight. However, in the Naimona’Nyi region, glacial retreat shows the opposite relationship with glacier orientation. Glaciers on the southern slopes retreated less than those on the northern slopes, and some even advanced during some periods. This is possibly because most of the south slope glaciers are located at higher elevations, which could also explain the different retreat patterns between the two slopes of the Naimona’Nyi glaciers.

The Naimona’Nyi main glacier terminal position has been constantly retreating over the period of this study. However, no obvious trend in retreat rate has been observed. This pattern differs from the glacier area change. For the Naimona’Nyi main glacier, the retreat rate was 19.1m/a between 2003 and 2013 based on remote sensing image, and significantly exceeds the retreat rate of 7.8m/a measured by Yao et al. during 2004 to 2006.16 The difference likely relates to the study methods,42 in that Yao et al. only measured part of the glacier terminus.16 Large annual variations of the glacier terminal position and somewhat different time periods of study may also prevent accurate direct comparisons. However, this retreat rate falls below that determined for the Panu glacier (1997 to 2007: 24.7m/a) and the Jicong Pu glacier (1977 to 2003: 48m/a).43,44

Our study also presented glacier surface elevation changes. Glacier thinning rates clearly vary with different periods, but they generally show rapid thinning since 2009. Compared with the results from other regions, the thinning rate remains moderate for the Naimona’Nyi main glacier. For instance, Zhang et al. report the East Rongbuk Glacier to be thinning at an estimated rate of 0.86m/a;45 Kääb et al. report thinning of 0.21±0.05m/a in the Kush-Karakoram-Himalaya region;20 and Wang et al. report 0.15m/a of thinning on the Urumqi No. 1 Glacier.46 Furthermore, our study found that glacier thinning rates varied with location on the glacier, especially with elevation. Below 5950 m, the thinning rate varied over a large range. Above that elevation, the thinning rate shows an apparent decreasing trend with elevation, and this trend held for all of our measurement periods.


Climate Controls Changes on the Naimona’Nyi Glacier

Climate changes, especially changes in air temperature and precipitation, directly relate to glacier changes.1,2 The following presents the climate change trend available from the Purang meteorological station (1973 to 2013), the station nearest to our study area, and compared it with glacier change.

The annual temperatures and summer temperatures increased at the rate of 0.43°C/10a and 0.27°C/10a, respectively, at the Purang station over the past 40 years, which are much greater than the entire Tibetan Plateau and the global average.47 However, annual precipitation and winter precipitation show slight decreasing trends, with decreases of 0.44 and 0.21mm/a, respectively (Fig. 6). Earlier studies showed that summer temperatures and winter precipitation provide the primary controls affecting these glacier changes.48,49 The annual variations of air temperature and precipitation did not match with the annual glacier changes, implying that the glacial change results more from long-term climate variations. Our study concludes that the accelerating retreat of the Naimona’Nyi glaciers mainly results from the significant climate warming and less so from the slight precipitation changes.

Fig. 6

Annual precipitation, winter precipitation (January to March) and annual averaged temperature, summer averaged temperature (June to August) changes at Purang between 1973 and 2013. (Dashed vertical lines indicate the period of 2003 to 2013).




The investigation examined glacier changes in the Naimona’Nyi region over the last decade (2003 to 2013) based on satellite images, ICESat data, and dGPS data. Our efforts concentrated on changes in glaciers’ areas, terminal positions, and surface elevations.

The area covered by glacier ice in the Naimona’Nyi region rapidly decreased, with a total area decrease of 16.5% over the past decade. The retreat rate shows an increasing trend, with the most rapid retreat occurring in the last 4 years. The terminal position of the main glacier retreated 191±35m in total in the past decade.

This study used in situ surveys on the main glacier surface and comparisons with the ICESat data to present precise changes in glacier surface elevations. The results show the glacier thinning rate is accelerating. The largest measured glacier thinning rate is 0.58±0.06m/a between 2009 and 2013, and the glacier surface elevation changes spatially vary on the glacier.

Glacier retreat correlates somewhat with slope, orientation, and area. Northward facing glaciers retreated faster than south facing ones. Low elevation zones of glaciers lost more area than the glacier’s high elevation zones. As previously recognized, small glaciers did not survive climate warming. By comparisons with meteorological data, we conclude that the Naimona’Nyi glacier retreat mostly results from long-term air temperature increases in the region and is less influenced from slight precipitation changes.


We thank two anonymous reviewers for constructive comments and suggestions so that we can revise and significantly improve this paper. This work is founded by the “Strategic Priority Research Program (B)” of the Chinese Academy of Sciences (Grant No. XDB 03030100), and the National Natural Science Foundation of China (Grant Nos. 41025002, 41371086, and 41101072). Major Special Project-the China High-Resolution Earth Observation System (30-Y30B13-9003-14/16-01).



T. P. Barnettet al., “Potential impacts of a warming climate on water availability in snow-dominated regions,” Nature 438, 303–309 (2005).NATUAS0028-0836http://dx.doi.org/10.1038/nature04141Google Scholar


J. Oerlemans, “Extracting a climate signal from 169 glacier records,” Science 308, 675–677 (2005).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1107046Google Scholar


N. M. Kehrwaldet al., “Mass loss on Himalayan glacier endangers water resources,” Geophys. Res. Lett. 35(22), 1–6 (2008).GPRLAJ0094-8276http://dx.doi.org/10.1029/2008GL035556Google Scholar


W. W. ImmerzeelL. P. H. BeekM. F. P. Bierkens, “Climate change will affect the Asian water towers,” Science 328(5984), 1382–1385 (2010).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1183188Google Scholar


M. F. Meieret al., “Glaciers dominate eustatic sea-level rise in the 21st century,” Science 317(5841), 1064–1067 (2007).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1143906Google Scholar


J. Thomaset al., “Recent contributions of glaciers and ice caps to sea level rise,” Nature 482, 514–518 (2012).NATUAS0028-0836http://dx.doi.org/10.1038/nature10858Google Scholar


Z. T. WangH. C. Su, “Glaciers in the world and china: distribution and their significance as water resources,” J. Glaciol. Geocryol. 25(5), 498–503 (2003).Google Scholar


A. V. Kulkarniet al., “Glacial retreat in Himalaya using Indian remote sensing satellite data,” Curr. Sci. 92(1), 69–74 (2007).CUSCAM0011-3891http://dx.doi.org/10.1117/12.694004Google Scholar


T. Cheet al., “Changes in glaciers and glacial lakes and the identification of dangerous glacial lakes in the Pumqu River Basin, Xizang (Tibet),” Adv. Meteorol. 2014, 1–8 (2014).1687-9309http://dx.doi.org/10.1155/2014/903709Google Scholar


S. D. Kaspariet al., “Snow accumulation rate on Qomolangma (Mount Everest), Himalaya: synchronicity with sites across the Tibetan Plateau on 50–100 year timescales,” J. Glaciol. 54(185), 343–352 (2008).JOGLAO0022-1430http://dx.doi.org/10.3189/002214308784886126Google Scholar


Z. G. Liet al., “Monitoring glacial variations based on remote sensing in the Luozha region, eastern Himalayas, 1980–2007,” Geogr. Res. 30(5), 939–952 (2011).http://dx.doi.org/10.11821/yj2011050018Google Scholar


L. L. Maet al., “Recent area and ice volume change of Kangwure glacier in the middle of Himalayas,” Chin. Sci. Bull. 55(20), 2088–2096 (2010).CSBUEF1001-6538http://dx.doi.org/10.1007/s11434-010-3211-7Google Scholar


Q. H. Yeet al., “Monitoring glacier and supra-glacier lakes from space in Mt. Qomolangma region of the Himalayas on the Tibetan Plateau in China,” J. Mt. Sci. 6(3), 211–220 (2009).1672-6316http://dx.doi.org/10.1007/s11629-009-1016-4Google Scholar


K. FujitaT. Nuimura, “Spatially heterogeneous wastage of Himalayan glaciers,” Proc. Natl. Acad. Sci. 108(34), 14011–14014 (2011).PMASAX0096-9206http://dx.doi.org/10.1073/pnas.1106242108Google Scholar


T. D. Yaoet al., “Different glacier status with atmospheric circulation in Tibetan Plateau and surroundings,” Nat. Clim. Change. 2(9), 663–667 (2012).NCCACZ1758-678Xhttp://dx.doi.org/10.1038/nclimate1580Google Scholar


T. D. Yaoet al., “Recent rapid retreat of the Naimona’Nyi glacier in southwestern Tibet Plateau,” J. Glaciol. Geocryol. 29(4), 503–508 (2007).Google Scholar


Y. T. Wanget al., “Monitoring glacier volume change based on Multi-source DEM and multi-temporal remote sensing images,” J. Glaciol. Geocryol. 32(10), 126–132 (2010).Google Scholar


Q. H. Yeet al., “Glacier variations in the Naimona’Nyi region, Western Himalaya, in the last three decades,” Ann. Glaciol. 43, 385–389 (2006).ANGLDN0260-3055http://dx.doi.org/10.3189/172756406781812032Google Scholar


IPCC, “Climate change 2013: the physical science basis,” Fifth assessment report Stockholm, Sweden, Cambridge University Press, Cambridge, United Kingdom, and New York (2013).Google Scholar


A. Kääbet al., “Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas,” Nature 488, 485–498 (2012).NATUAS0028-0836http://dx.doi.org/10.1038/nature11324Google Scholar


N. Neckelet al., “Glacier mass changes on the Tibetan Plateau 2003–2009 derived from ICESat laser altimetry measurements,” Environ. Res. Lett. 9, 1–7 (2014).ERLNAL1748-9326http://dx.doi.org/10.1088/1748-9326/9/1/014009Google Scholar


Y. F. Shi, “Glacier distribution in the flow area,” Chapter 4 in Concise Glacier Inventory of China, Y. F. Shi, Ed., pp. 84–89, Shanghai Popular Science Press, Shanghai, China (2005).Google Scholar


United States Geological Survey, “Landsat Enhanced Thematic Mapper,” http://glovis.usgs.gov/ (14 July 2003).Google Scholar


National Aeronautics and Space Administration, “Advanced Spaceborne Theemal Emission and Reflection Radiometer Global Digital Elevation Model,” http://reverb.echo.nasa.gov (17 October 2011).Google Scholar


R. S. GautamD. SinghA. Mittal, “Application of principal component analysis and information fusion technique to detect hotspots in NOAA/AVHRR images of Jharia coalfield, India,” J. Appl. Remote Sens. 1(1), 013523 (2007).1931-3195http://dx.doi.org/10.1117/1.2771256Google Scholar


A. Y. Shiet al., “Multispectral and panchromatic image fusion based on improved bilateral filter,” J. Appl. Remote. Sens. 5(1), 053542 (2011).1931-3195http://dx.doi.org/10.1117/1.3616010Google Scholar


T. BolchT. D. YaoS. C. Kang, “A glacier inventory for the western Nyainqentanglha Range and the Nam Co Basin, Tibetan, and glacier changes 1976–2009,” Cryosphere 4, 419–433 (2010).1994-0416http://dx.doi.org/10.5194/tc-4-419-2010Google Scholar


L. P. YangD. S. XiaF. H. Chen, “On fusion algorithms of Landsat 7 ETM+ pan-multispectral images,” J. Lanzhou. U. Nat. Sci. 43, 7–11 (2007).LCTHAF0455-2059Google Scholar


T. Bolch, “Climate change and glacier retreat in northern TienShan (Kazakhstan/Kyrgyzstan) using remote sensing data,” Global Planet Change 56, 1–12 (2007).GPCHE40921-8181http://dx.doi.org/10.1016/j.gloplacha.2006.07.009Google Scholar


W. SilverioJ. M. Jaquet, “Glacial cover mapping (1987–1996) of the Cordillera Blanca (Peru) using satellite imagery,” Remote Sens. Environ. 95(3), 342–350 (2005).RSEEA70034-4257http://dx.doi.org/10.1016/j.rse.2004.12.012Google Scholar


Q. H. Yeet al., “Monitoring glacier variations on Geladandong Mountain, central Tibetan Plateau, from 1969 to 2002 using remote-sensing and GIS technologies,” J. Glaciol. 52(179), 537–545 (2006).JOGLAO0022-1430http://dx.doi.org/10.3189/172756506781828359Google Scholar


T. Bolchet al., “Mass loss of Greenland’s glaciers and ice caps 2003–2008 revealed from ICESat laser altimetry data,” Geophys. Res. Lett. 40, 875–881 (2013).GPRLAJ0094-8276http://dx.doi.org/10.1002/grl.v40.5Google Scholar


National Snow and Ice Data Center, “Geoscience Laser Altimeter System,” http://nsidc.org/ (13 January 2003).Google Scholar


C. B. Fan, Space-Borne Laser Altimeter Theory and Application, Wuhan University Press, Wuhan (2007).Google Scholar


H. J. Wenet al., “Accuracy assessment of ICESat laser altimeter data using GPS measurements,” Geo. Inf. Sci. Wuhan Univ. 36(3), 262–266 (2011).Google Scholar


Y. F. ShiS. Y. Liu, “Estimation on the response of glaciers in China to the global warming in the 21st century,” Chin. Sci. Bull. 45(7), 668–672 (2000).CSBUEF1001-6538http://dx.doi.org/10.1007/BF02886048Google Scholar


D. A. Carloet al., “Glacier shrinkage driven by climate during half a century (1954–2007) in the Ortles-Cevedale group (Stelvio National Park, Lombardy, Italian Alps),” Theor. Appl. Clim. 116, 169–190 (2013).TACLEK1434-4483http://dx.doi.org/10.1007/s00704-013-0938-5Google Scholar


S. DuttaA. L. RamanathanA. Linda, “Glacier fluctuation using satellite data in Beas basin, 1972–2006, Himachal Pradesh, India,” J. Earth Syst. 121(5), 1105–1112 (2012).0253-4126http://dx.doi.org/10.1007/s12040-012-0219-1Google Scholar


C. X. LiT. B. YangH. Z. Tian, “Variation of West Kunlun mountains glacier during 1990–2011,” Prog. Geogr. 32(4), 548–559 (2013).0309-1325http://dx.doi.org/10.11820/dlkxjz.2013.04.007Google Scholar


B. L. Liet al., “Glacier change over the past four decades in the Middle Chinese Tien Shan,” J. Glaciol. 52(178), 1–9 (2006).JOGLAO0022-1430http://dx.doi.org/10.3189/172756506781828557Google Scholar


X. Wanget al., “Glacier temporal spatial change characteristics in Western Nyainqentanglha range, Tibetan Plateau 1977–2010,” J. Chin. Univ. Geosci. 37(5), 1082–1092 (2012).JCUGEXGoogle Scholar


Y. Nieet al., “Monitoring glacier change based on remote sensing in the Mt. Qomolangma National Nature Preserve, 1976–2006,” Acta Geogr. Sin. 65(1), 13–28 (2010).Google Scholar


T. Cheet al., “Monitoring glacier and associated glacial lakes on the east slopes of mount Xixabangma from remote sensing images,” J. Glaciol. Geocryol. 27(6), 801–805 (2005).Google Scholar


S. C. Kanget al., “Glacier retreating dramatically on the Mt. Nyainqentanglha during the last 40 years,” J. Glaciol. Geocryol. 39(6), 869–873 (2007).Google Scholar


T. Zhanget al., “Ice thickness observation and landform study of East Rongbuk glacier, Mt. Qomolangma,” J. Glaciol. Geocryol. 34(5), 1059–1066 (2012).Google Scholar


P. Y. Wanget al., “Changes of ice-thickness and volume for representative glaciers in Tienshan Mountains in the past 50 years,” Acta Geogr. Sin. 67(7), 929–940 (2012).http://dx.doi.org/10.11821/xb201207006Google Scholar


B. Wanget al., “Tibetan Plateau warming and precipitation changes in East Asia,” Geophys. Res. Lett. 35, 1–5 (2008).GPRLAJ0094-8276http://dx.doi.org/10.1029/2008GL033375Google Scholar


R. Bonannoet al., “Glacier response to current climate change and future scenario in the northwestern Italian Alps,” Region. Environ. Change 14(2), 633–643 (2014).1436-3798http://dx.doi.org/10.1007/s10113-013-0523-6Google Scholar


N. L. Wanget al., “Variations in equilibrium line elevation of the Qiyi Glacier, Qilian Mountains, over the past 50 years,” Chin. Sci. Bull. 55(32), 3107–3115 (2010).CSBUEF1001-6538http://dx.doi.org/10.1007/s11434-010-4167-3Google Scholar


Dayun Zhu is a PhD candidate in the School of Geographical Science, Southwest University, China. He is mainly engaged in research on snow remote sensing, glaciers, and climate change.

Lide Tian is a professor in the Institute of Tibetan Plateau Research, Chinese Academy of Sciences, China. He is mainly engaged in the research of third pole environment, Tibetan Plateau environment changes, and land surface processes.

Jianli Wang is a professor in the School of Geographical Science, Southwest University, China. He is mainly engaged in the research of sedimentology and sedimentary environment evolution, landform, and geographic information systems.

Yong Wang is a professor in the School of Geographical Science, Southwest University, China. He is mainly engaged in the research of paleoclimate and environmental pollution.

Jiangpeng Cui is a PhD candidate in the Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences. He is mainly engaged in the research of isotope and climate change.

Dayun Zhu, Lide Tian, Jianli Wang, Yong Wang, Jiangpeng Cui, "Rapid glacier retreat in the Naimona’Nyi region, western Himalayas, between 2003 and 2013," Journal of Applied Remote Sensing 8(1), 083508 (9 December 2014). http://dx.doi.org/10.1117/1.JRS.8.083508
Submission: Received ; Accepted

Climate change

Temperature metrology


Remote sensing

Earth observing sensors

Image fusion


Back to Top