Skip to main content


  • Access
  • Cited by 30
  • Cited by
    This article has been cited by the following publications. This list is generated based on data provided by CrossRef.

    Kulkarni, Anil V. and Pratibha, S. 2018. Science and Geopolitics of The White World. p. 183.

    Khan, Abul Amir Pant, N. C. and Ravindra, Rasik 2018. Science and Geopolitics of The White World. p. 161.

    Dobreva, Iliyana Bishop, Michael and Bush, Andrew 2017. Climate–Glacier Dynamics and Topographic Forcing in the Karakoram Himalaya: Concepts, Issues and Research Directions. Water, Vol. 9, Issue. 6, p. 405.

    Brahmbhatt, Rupal M. Bahuguna, I. M. Rathore, B. P. Kulkarni, Anil V. Shah, Rajesh D. Rajawat, A. S. and Kargel, Jeffrey S. 2017. Significance of glacio-morphological factors in glacier retreat: a case study of part of Chenab basin, Himalaya. Journal of Mountain Science, Vol. 14, Issue. 1, p. 128.

    Kulkarni, Anil V. Nayak, Shailesh and Pratibha, S. 2017. Observed Climate Variability and Change over the Indian Region. p. 193.

    King, Owen Quincey, Duncan J. Carrivick, Jonathan L. and Rowan, Ann V. 2017. Spatial variability in mass loss of glaciers in the Everest region, central Himalayas, between 2000 and 2015. The Cryosphere, Vol. 11, Issue. 1, p. 407.

    Garg, Purushottam Kumar Shukla, Aparna Tiwari, Reet Kamal and Jasrotia, Avtar Singh 2017. Assessing the status of glaciers in part of the Chandra basin, Himachal Himalaya: A multiparametric approach. Geomorphology, Vol. 284, p. 99.

    Mir, Riyaz Ahmad and Majeed, Zahid 2016. Frontal recession of Parkachik Glacier between 1971-2015, Zanskar Himalaya using remote sensing and field data. Geocarto International, p. 1.

    Nainwal, H.C. Banerjee, Argha Shankar, R. Semwal, Prabhat and Sharma, Tushar 2016. Shrinkage of Satopanth and Bhagirath Kharak Glaciers, India, from 1936 to 2013. Annals of Glaciology, Vol. 57, Issue. 71, p. 131.

    Miles, Evan S. Pellicciotti, Francesca Willis, Ian C. Steiner, Jakob F. Buri, Pascal and Arnold, Neil S. 2016. Refined energy-balance modelling of a supraglacial pond, Langtang Khola, Nepal. Annals of Glaciology, Vol. 57, Issue. 71, p. 29.

    Banerjee, Argha and Azam, Mohd Farooq 2016. Temperature reconstruction from glacier length fluctuations in the Himalaya. Annals of Glaciology, Vol. 57, Issue. 71, p. 189.

    Gaddam, Vinay Kumar Kulkarni, Anil V. and Kumar Gupta, Anil 2016. Estimation of glacial retreat and mass loss in Baspa basin, Western Himalaya. Spatial Information Research, Vol. 24, Issue. 3, p. 257.

    Vijay, Saurabh and Braun, Matthias 2016. Elevation Change Rates of Glaciers in the Lahaul-Spiti (Western Himalaya, India) during 2000–2012 and 2012–2013. Remote Sensing, Vol. 8, Issue. 12, p. 1038.

    Chand, Pritam Sharma, Milap Chand and Prasad, Ram Nagesh 2016. Climate Change, Glacier Response, and Vegetation Dynamics in the Himalaya. p. 155.

    Chand, Pritam and Sharma, Milap Chand 2016. Monitoring Frontal Changes of Shah Glacier in the Ravi Basin, Himachal Himalaya (India) from 1965 to 2013. National Academy Science Letters, Vol. 39, Issue. 2, p. 109.

    Raj, K. Babu Govindha and Kumar, K. Vinod 2016. Inventory of Glacial Lakes and its Evolution in Uttarakhand Himalaya Using Time Series Satellite Data. Journal of the Indian Society of Remote Sensing, Vol. 44, Issue. 6, p. 959.

    Sharma, Parmanand Patel, Lavkush K Ravindra, Rasik Singh, Ajit K, Mahalinganathan and Thamban, Meloth 2016. Role of debris cover to control specific ablation of adjoining Batal and Sutri Dhaka glaciers in Chandra Basin (Himachal Pradesh) during peak ablation season. Journal of Earth System Science, Vol. 125, Issue. 3, p. 459.

    2016. Recent Climate Change Impacts on Mountain Glaciers. p. 152.

    Thakuri, Sudeep Salerno, Franco Bolch, Tobias Guyennon, Nicolas and Tartari, Gianni 2016. Factors controlling the accelerated expansion of Imja Lake, Mount Everest region, Nepal. Annals of Glaciology, Vol. 57, Issue. 71, p. 245.

    Shukla, Aparna and Qadir, Junaid 2016. Differential response of glaciers with varying debris cover extent: evidence from changing glacier parameters. International Journal of Remote Sensing, Vol. 37, Issue. 11, p. 2453.




      • Send article to Kindle

        To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India
        Available formats
        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your Dropbox account. Find out more about sending content to Dropbox.

        The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India
        Available formats
        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your Google Drive account. Find out more about sending content to Google Drive.

        The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India
        Available formats
Export citation



We investigated area changes in glaciers covering an area of ∼200 km2 in the Tista basin, Sikkim, Eastern Indian Himalaya, between ∼1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 ± 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 ± 0.82, 2.56 ± 0.61 and 2.28 ± 2.01 km2 for the periods 1989–97, 1997–2004/05 and 2004–2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987–2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0°C rise in the summer minimum temperature and a 2.0°C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.


The glaciers of the Himalaya contribute to a varying degree to the overall river runoff in South Asia. The highest contribution is from the Indus River which originates in the northwest Himalaya (Immerzeel and others, 2010). The contribution of glacier meltwater to the overall runoff, which is higher in the upper reaches of the rivers (Kaser and others, 2010), constitutes the main water source for micro and mini hydroelectric projects, irrigation, and domestic and adventure tourism for different regions (Krishna, 2005). However, detailed information about the importance of glacier melt is still rare, though precise data about glacier area and volume, their changes over time and the causes of these changes exist. While glacier area changes have been comparatively well documented in the western Indian Himalaya (Kulkarni and others, 2011; Bolch and others, 2012), little information is available for the eastern Himalaya including the Sikkim Himalayan glaciers. Many glaciers in the Himalaya are covered with supraglacial debris that alters their response to climate change (Scherler and others, 2011; Bolch and others, 2012). While a thick debris cover is expected to reduce glacier retreat (Scherler and others, 2011), several recent studies have shown that debris-covered glaciers have lost significant mass in recent decades or are clearly retreating despite thick debris cover (Bolch and others, 2011; Kääb and others, 2012; Nuimura and others, 2012; Bhambri and others, 2013). This has been attributed to enhanced ice melt on glacial lakes and ice cliffs (Sakai and others, 2000; Bolch and others, 2012).

An extensive glacier inventory for Sikkim Himalaya using Indian Remote-sensing Satellite IRS-1C/1D LISS III (Linear Imaging Self-Scanning sensor III) data of 1 January 1997 mapped 84 glaciers occupying an area of 440.30 km2 and estimated a volume of 140 km3 (Bahuguna and others, 2001). For the volume estimate, the area–volume scaling method suggested by Meier (1961) and UNESCO/International Association of Scientific Hydrology (1970) was adopted. In the inventory, the majority of glaciers were smaller than 5 km2 and only two were larger than 20 km2. The average altitude of the snowline was ∼6000 m a.s.l. in January 1990. The Sikkim glaciers comprise valley glaciers with elevations ranging from 3720 to ∼8400 m a.s.l. In addition, permanent snowfields covering an area of 251.22 km2 were mapped (Bahuguna and others, 2001). Many glaciers in Sikkim are covered by debris. The largest and most well-known glacier of Sikkim is Zemu Glacier (Raina and others, 1973) which is situated below Mount Kanchenjunga (8586 m a.s.l.) and was first mapped during an expedition of German researchers (Finsterwalder, 1935). Changme Khangpu, one of the Sikkim glaciers, has been studied using Si-32 and Pb-210 chronology, and the results showed (1) a retreat rate of 6 m a−1 between 1976 and 1981, (2) an accumulation rate of 0.7 m a−1, and (3) a basal flow of 0.01 m a−1 compared to a surface flow of 40 m a−1, indicating that the deeper ice near the bedrock travels much more slowly than the surface ice (Shukla and others, 1983; Nijampurkar, 1985; Nijampurkar and others, 1985). Another study estimated the average rate of glacier retreat as 13 m a−1, for 26 glaciers in Sikkim, using satellite images between 1976 and 2005 (Raina, 2009). Snow studies using the normalized difference snow index (NDSI) method showed a maximum snow-cover extent of 50% in February, with accumulation peaks in October and February. The snow extent was comparatively high (35–40%) even in the summer months, indicating different snow accumulation and ablation patterns in Sikkim as compared to the western Himalaya (Basnett and others, 2011). However, there is still very little information on glacier changes in Sikkim. The remote and rugged terrain of Sikkim Himalaya and the politically restricted boundaries inhibit the study of mass-balance and glacier fluctuations. In addition, supraglacial and moraine-dammed lakes, which are common in Sikkim (Raj and others 2012), influence glacier changes.

The aims of this paper are therefore (1) to present a detailed and long-term analysis of glacier area changes in Sikkim Himalaya, and (2) to understand the influence of debris cover and moraine-dam lakes on glacier retreat.

Study Area

We focus on Sikkim Himalaya where many debris-covered glaciers are present. Sikkim is a part of the Eastern Indian Himalayan ranges and lies between Nepal and Bhutan (Fig. 1). It is a mountainous state of India, extending ∼114 km from north to south and 64 km from east to west, with a total geographical area of 7096 km2. The state is situated between 27°00′46″ and 28°07′48″ N and 88°00′58″ and 88°55′25″ E. The mountainous topography ranges from 300 to 8598 m a.s.l., covering subtropical to alpine ecozones, within a short horizontal distance. The state receives an annual precipitation of 2000–4000 mm (, while the Asian monsoon contributes >80% of precipitation during the summer months (Bookhagen and Burbank, 2010). Most of the geographical area of Sikkim is drained by the Tista River, a tributary of the Brahmaputra River.

Fig. 1. Location map of Sikkim. The region is located between Nepal and Bhutan. The map shows locations of glaciers in the Tista basin, Sikkim. The meteorological data of Gangtok were analysed.

Data and Methods

Data sources

Landsat Thematic Mapper (TM) (Level 1T) and Indian Remote-sensing Satellite (IRS) LISS III satellite images were analysed to map and monitor the glacier changes (Tables 1 and 2). The LISS III images were co-registered to the Landsat images and used for delineating glacier boundaries. In addition, the Landsat TM (Level 1T) and topographic maps of the Survey of India were used as reference data. The Landsat images were downloaded from the US Geological Survey (USGS) website (, and the topographic map was acquired from the Sikkim State Council of Science and Technology and the Department of Mines, Minerals and Geology, Government of Sikkim. The glacier extents from the topographic maps, surveyed in 1962/63 and 1966/67, indicate a general recession. However, we did not use these outlines for a quantitative assessment, as glacier boundaries on the topographic maps are uncertain due to errors in the identification of debris-covered glacier tongues and snow-fields (Bhambri and Bolch, 2009).

Table 1. Satellite and sensor characteristics

Table 2. List of satellite data used in the investigation

In the western Himalaya where the summer monsoon is of less influence than in the east, the glaciers are usually mapped using satellite images of the July–September season, as the seasonal snow cover is at a minimum and the glaciers are fully exposed (Kulkarni, 1991). However, in Sikkim, in the eastern Himalaya where little precipitation occurs during winter, images from November–December and January have been selected based on minimal seasonal snow and cloud cover (Basnett and Kulkarni, 2012). The images were supplemented with the data of the March and October season. In addition, multi-date scenes for the same year were used to overcome issues related to shadows.

Glacier mapping and analysis

We manually delineated the glaciers and moraine-dam lakes using false-colour composites with band combinations of red, near-infrared (NIR) and shortwave infrared (SWIR), and a true-colour composite of red, green and NIR, of Landsat TM/ETM+ and LISS III images. We also used contrast techniques such as brightness and linear stretching, and a shaded relief map derived from the Shuttle Radar Topography Mission digital elevation model (SRTM DEM) (90 m; to aid the manual interpretation of glacial features. In addition, we referred the Survey of India maps for the identification, location and interpretation of glaciers on satellite images. We also used the NDSI (NDSI = (green – SWIR)/(green + SWIR)) applied on an Advanced Wide Field Sensor (AWiFS) image of the years 2004 and 2005 to eliminate non-glacier areas under mountain shadow conditions (cf. Kulkarni and others, 2006). Glacier boundaries in the accumulation area were demarcated by ‘ice divides’, which were visually detected and manually digitized based on the shaded relief map.

The glacier area was determined for four different periods: 1989/90, 1997, 2004 and 2009/10. The 1990 glacier boundary was used as a ‘master inventory’ because of good image quality. As the oldest data layer, it was used to estimate glacier retreat for the following years. The changes were estimated only in the lower region, as very little or no change was observed in the upper regions. Moreover, since the upper regions including the accumulation area were generally snow-covered on the available images, it was difficult to identify small changes in those areas.

Estimation of uncertainties

Estimating the uncertainty is crucial to assess the significance of the results and to avoid misinterpretation. The accuracy of the outlines depends, typically, on the resolution of the images used, the conditions during image acquisition (seasonal snow, shadow) and the contrast between the glacier and its surroundings (DeBeer and Sharp, 2007). The latter is especially low for debris-covered glaciers. The most accurate way to assess glacier outlines would be high-resolution imagery (Paul and others, 2013). However, such data were, as in most studies in remote mountainous areas, not available for our study region. A widely used and recommended approach is to estimate the uncertainty using a buffer of up to one pixel size to the glacier margin (Granshaw and Fountain, 2006; Bolch and others, 2010; Paul and others, 2013). In this study, uncertainties in the estimation of ‘glacier extent’ and ‘retreat’ were considered due to (1) the spatial resolution of satellite data and (2) the lack of a distinct difference between glacial and non-glacial areas caused by a debris cover. These are called ‘Mapping’ and ‘Interpretation’ uncertainties. We estimate the mapping uncertainty as


where N is the number of pixels along the glacier boundary and A is the area of the pixel.

The interpretation uncertainty was estimated manually by mapping possible alternative glacier boundaries (Fig. 2). This technique was applied for the first time while interpreting glacier boundaries around the snout, especially for debris-covered glaciers. Alternative snout positions were estimated using a visual interpretation technique, based on shape, size, pattern, tone, texture and association (Kulkarni, 1991). Two extreme positions were identified at (a) ‘G max’ and (b) ‘G min’ snout positions as shown in Figure 2. These were used to estimate uncertainty in the interpretation:


where G max is the maximum possible glacier extent and G min is the minimum possible glacier extent.

However, the uncertainty of glacier shrinkage is likely to be overestimated when calculating uncertainty based on the method mentioned above, because, for each time period, the vast majority of the area changes are usually restricted to the snout. Hence, to estimate uncertainties in area changes between two time periods, the changes in the areal extent near the snout are considered. These uncertainties in glacier area changes can be written as

where N S is the number of pixels along the boundary of the snout where changes occurred.

Fig. 2. The identification of a glacier terminus can sometimes be difficult using satellite imagery, if glaciers are covered by debris. In this study, alternative snout positions have been estimated using a ‘visual interpretation technique’, based on tone, texture and colour. This image shows a true-colour composite with a band combination of green, red and NIR. Two extreme positions were identified at (a) ‘Gmax’ and (b) ‘Gmin’ snout positions. This was used to estimate uncertainty in the interpretation. A satellite image of LISS III, IRS-P6, 25 March 2004, path 107 row 52, is shown.


Glacier area in Sikkim

In this study, 39 glaciers were mapped and monitored, to estimate the loss in glacier area (Fig. 1). Though the Sikkim glacier inventory by Bahuguna and others (2001) mapped 84 glaciers covering 440.31 km2, in this study not all the glaciers could be mapped due to a cloud and snow cover. Nevertheless, the 39 glaciers mapped occupy an area of 333.0 ± 9.1 km2 in 1989/90, which represents ∼75% of the glacierized area in the Tista basin. The overall loss in glacier area in 20 years, from 1989/90 to 2009/10, was 6.9 ± 1.5 km2 (∼2%). Zemu Glacier (Fig. 1, glacier No. 3) covered 37% (∼131 km2) of the investigated glacier area and showed no significant change of the terminus in the 1990, 1997, 2004 and 2010 images. This is possibly due to its large size, steep slope, glacier dynamics and thick debris cover at the tongue (Venkatesh and others, 2011). Zemu Glacier was not included in the change analysis, in order to understand the fluctuations of small and medium-sized glaciers. The remaining 38 glaciers covered 201.9 ± 7.0 km2 in 1989/90 and showed an overall loss of 6.9 ± 1.5 km2 (∼3%) until 2010 (∼0.33 km2 a−1).The rate of change was −0.16 ± 0.10% a−1 from 1990 to 2010, with an increasing rate after 1997 (Tables 3 and 4).

Table 3. Total glacier area for 38 glaciers in the Tista basin

Table 4. Area loss for 38 glaciers in the Tista basin

A good correlation (r 2 = 0.25) was observed between glacier size and the relative area loss for debris-free glaciers, indicating a lower loss for large glaciers. However, the inclusion of debris-covered glaciers weakened the correlation (r 2 = 0.10) between the size and area loss (Fig. 3), suggesting a strong influence of debris on the glacier retreat. Therefore, we carried out a more detailed analysis for debris-free glaciers (‘clean’), debris-covered glaciers without lakes (‘debris’), and debris-covered glaciers with moraine-dammed lakes (‘debris with lakes’) to estimate the influences of these parameters on glacier retreat.

Fig. 3. The influence of glacial size on retreat for debris-covered and clean glaciers (debris-free) between 1989 and 2010. The clean glaciers show a high loss in glacier area for small-sized glaciers, but the debris-covered glaciers show variable retreat and are independent of glacier size.

Glacier comparison: ‘clean’, ‘debris’ and ‘debris with lakes’

In 1989/90, 12 out of 38 glaciers had no debris cover (‘clean’), 12 were covered by supraglacial debris (‘debris’), while 14 were debris-covered and associated with moraine-dam lakes (‘debris with lakes’). Between 1989/90 and 2010, we observed (1) an increase in the area of debris cover by 6.5 ± 1.4 km2 and (2) the formation and expansion of supraglacial lakes on two glaciers (Fig. 4c–f). These lakes joined and developed moraine-dam lakes on two debris-covered glaciers (Fig. 5), indicating a transformation in the characteristics of the glaciers, as debris cover has promoted the development of lakes. These moraine-dam lakes developed at the glacier termini at certain altitudes (∼6300 m a.s.l.) on the medium-sized ‘debris-covered glaciers’ (Fig. 4f). A field photograph showing moraine-dam lakes at the terminus of Tista Khangse Glacier (No. 30 in Fig. 1) is shown in Figure 6.

Fig. 4. The distribution and number of ‘clean’, ‘debris’ and ‘debris with lakes’ glaciers: (a) ‘clean’ glaciers (1990) (n = 12); (b) ‘clean’ glaciers (2010) (n = 12); (c) ‘debris-covered glaciers without lakes’ (1990) (n = 10); (d) ‘debris-covered glaciers without lakes’ (2010) (n = 10); (e) ‘debris-covered glaciers with lakes’ (1990) (n = 14); (f) ‘debris-covered glaciers without lakes’ (2010) (n = 16). (f) shows an increase in the number of lakes on debris-covered glaciers over 20 years, between 1990 and 2010.

Fig. 5. Evolution and coalescence of a supraglacial lake and the formation of a moraine dam: (a, b) show no frontal change between 1990 and 1997; (b) shows the evolution of a supraglacial lake and (c) the coalescence of a supraglacial lake, which occupies glacier area between two lateral moraines; and (d) shows the formation of a moraine dam lake leading to glacier area loss. The four images used are a false-colour composite with a band combination of red, NIR and SWIR.

Fig. 6. Moraine-dam lake (Khangchung Chho) at the terminus of Tista Khangse Glacier (No. 30 in Fig. 1).

Glacier area changes

We observed a large difference in the area changes of (1) clean (mean size 3.2 km2), (2) debris (mean size 6 km2) and (3) debris-covered glaciers with lakes (mean size 6.4 km2) between 1989/90 and 2009/10. The 12 ‘clean’ glaciers lost 5.3 ± 1.0% of their area, the 10 ‘debris-covered glaciers’ lost 3.9 ± 1.3%, and the 16 ‘debris with lakes’ lost 9.0 ± 1.4% (Fig. 7). Hence, glaciers with ‘debris’ appear to shrink at a lower rate than ‘debris-free’ glaciers, but the ‘debris with lakes’ glaciers lose area at a much higher rate. This is about 1.7 times higher than ‘clean’ and 2.3 times higher than ‘debris-covered glaciers’, suggesting that the debris cover retards melting while the presence of lakes on debris-covered glaciers accelerates glacier retreat (Fig. 7).

Fig. 7. The influence of debris-covered and moraine-dam lakes on glacier retreat. The numbers of glaciers used to estimate the percentage loss in area are 12, 10 and 16 for ‘clean’, ‘debris’ and ‘debris with lakes’, respectively. The figure suggests that debris cover can retard retreat. However, debris cover can also lead to the formation of lakes on moraines and this can accelerate retreat.

Altitude zone analysis

The 38 glaciers in this study ranged from ∼4200 to ∼7600 m, with an average altitude of ∼5780 m. The mean altitude was (1) 5783 m for ‘clean’, (2) 5610 m for ‘debris-covered’ and (3) 5825 m for ‘debris with lakes’ glaciers. The altitude distribution pattern (Fig. 8) showed that despite similar glacier altitudes of ‘clean’ and ‘debris with lakes’, the ‘debris with lakes’ glaciers lost area at a higher rate (Fig. 7), indicating an accelerated retreat of debris-covered glaciers when associated with lakes. We also observed a large difference in the altitude location of glaciers with ‘debris’ versus ‘debris with lakes’. The glaciers with ‘debris’ are located ∼200 m lower than ‘debris with lakes’ glaciers (Fig. 9), but showed a three times lower retreat (Fig. 7). This indicates that the presence of a water body influences and accelerates melting, while altitudinal differences have a clearly lower influence. Moreover, though the mean size of ‘debris’ glaciers and ‘debris with lakes’ was comparatively similar, a higher area loss was observed on ‘debris with lakes’ glaciers, suggesting the influence of lakes on glacier retreat.

Fig. 8. The glaciers associated with lakes showed a higher loss in glacier area for all three observational periods, suggesting the influence of lakes on glacier retreat.

Fig. 9. Area–altitude distribution of glaciers in Sikkim. The ‘debris-covered glaciers with lakes’ are distributed at higher altitudes than ‘debris-covered glaciers without lakes’. The ‘clean’ and ‘debris-covered glaciers with lakes’ are distributed almost equally in the higher regions.

Lake area changes

In 2010, 16 lakes were associated with debris-covered glaciers, and 7 were associated with debris-free glaciers. The 16 debris-covered glaciers and the 7 debris-free glaciers associated with lakes both showed a two times higher loss in glacier area than the glaciers without lakes. Thus, of the 38 glaciers studied in this investigation, moraine-dam lakes were associated with 23 and these glaciers showed a five times greater loss in glacier area, suggesting the influence of lakes on glacier retreat (Fig. 8). The areal extent of the individual lakes ranged between 0.04 and 2.10 km2, and the overall lake area increased from 6.6 ± 0.8 km2 to 9.6 ± 1.1 km2, between 1990 and 2010. This increase in area was due to the expansion of moraine-dammed lakes, but also to the formation and coalescence of supraglacial lakes on two debris-covered glaciers. An example of lake expansion, near the terminus of South Lhonak Glacier (glacier No. 11 in Fig. 1), is shown in Figure 10.

Fig. 10. The moraine-dammed lake at the terminus of South Lhonak Glacier (glacier No. 11 in Fig. 1) has expanded from 0.50 ± 0.05 km2 to 1.08 ± 0.07 km2 between 1990 and 2010. The expansion of the lake has reduced the glacier area by 4.39 ± 0.71%.


Influence of glacier size on glacier retreat

In this study, the relationship between the glacier size and the area loss was strong for ‘clean’ glaciers. These glaciers showed a higher loss in glacial area for smaller glaciers in comparison with larger ones, which is consistent with the results in many other mountain regions in the Himalaya (e.g. the Baspa River and the Bhagirati and Alknanda River basins (Kulkarni and Alex, 2003; Bhambri and others, 2011)) and elsewhere in the world (e.g. Paul and others, 2004; Bolch and others, 2010). However, the correlation between size and loss in area was weak for debris-covered glaciers and even weaker for ‘debris with lakes’ glaciers, indicating a different reaction for the latter. Moreover, the loss in area for the ‘debris with lakes’ glaciers was twice as high as for the ‘debris’ glaciers (Fig. 7), though they were located almost at the same altitude (Fig. 9) and were of similar size (∼6 km2). This suggests that glacier retreat is affected much more by the presence of a water body on debris-covered glaciers than by altitude. Therefore, in the following subsections, detailed discussions are carried out to understand the influence of debris cover and lakes on glacier retreat.

Influence of debris cover and glacial lakes on glacier variations

Thick debris cover on glaciers leads, in general, to less frontal retreat or even stable tongues, as it lowers the rate of surface melting, as in Zemu Glacier and many other glaciers in the Himalaya (Tangborn and Rana, 2000; Bolch and others, 2008; Nainwal and others, 2008; Racoviteanu and others, 2008a,b; Bhambri and others, 2011; Scherler and others, 2011). However, this study suggests that debris-covered glaciers can have a large area loss despite a thick debris cover, if associated with lakes. This fact is important for Sikkim Himalaya as the majority of debris-covered glaciers in this region have moraine-dam lakes (Fig. 4f). Initially, as observed in this study, some debris-covered glaciers had stagnant fronts (Fig. 5), but over a period of time, lakes developed (1) near the terminus and (2) on the middle regions of the glacier tongues. These lakes, when formed at the terminus, merged and developed into moraine-dam lakes. In addition, the supraglacial lakes formed on the interiors of the glacier tongues, united and extended their area onto the lateral margins of the glaciers. When the boundaries of the lakes touch the lateral boundaries of the glacier, they cease to be a part of the glacier. The ice beyond the lake then becomes inactive and is separated from the glacier, leading to a loss in glacier area. This phenomenon was well observed in Changsang Glacier (17.9% debris-covered; Fig. 5), where between 1990 and 1997 a supraglacial lake was part of the glacier. After 2004, however, it was no longer part of the glacier area, as the lake periphery collided with the lateral boundaries of the glacier, resulting in a loss of glacial area, caused by the formation of lakes and dead ice mounds. Once a stable lake is formed, it leads to the transfer of latent heat flux at the interface of water and ice (Sakai and others, 2000). In addition, a lake can absorb more solar radiation due to the lower albedo of water. These energy exchange processes cause the melting of ice below the lake water and finally its disappearance, if the lake is large and stable over a long period. This situation was also observed at several other debris-covered glaciers in the Himalaya and described for Imja Lake in the Mount Everest region (Fujita and others, 2009; Watanabe and others, 2009).

An increase in area of the debris cover on glaciers has also been detected elsewhere in the Himalaya (Bolch and others, 2008; Bhambri and others, 2011; Kamp and others 2011), indicating increased debris production, which reduces glacier transport capacity and leads to negative mass balances (Benn and others, 2012; Bolch and others, 2012).

Debris cover can cause differential rates of melting. Near the terminus, a thick debris cover retards melting, but at the upper reaches of the glacier tongue the thinner debris cover causes higher rates of melting. This causes a lessening, or even reversal, of the surface gradient (cf. Bolch and others 2011; Benn and others 2012) and leads to the formation of supraglacial lakes (Bennett and Glasser, 2009). These lakes then become prone to glacier melt (Sakai and Fujita, 2010; Benn and others, 2012). They may be temporary and drain if they reach the englacial drainage system (‘perched lakes’; Gulley and Benn, 2007; Benn and others, 2012). However, if the lakes reach the base level, individual lakes join together and form a larger lake which can be dammed behind the terminal moraine (‘base-level lakes’), as observed, for example, at Jonsung Glacier, Sikkim (glacier 10 in Fig. 1).

Impact of lake formation

Moraine-dam lakes have a strong impact on glacier dynamics. The lakes at the terminus of glaciers modify the stress regime of the glacier ice in contact. The presence of the lakes increases the terminus disintegration by the process of calving, and also transmits thermal energy to the ice front, accelerating melting which leads to a loss in glacier area (Komori, 2008; Röhl, 2008; Sakai and others, 2009; Gardelle and others, 2011). Detailed investigations in the Mount Everest area show a significant larger retreat and volume loss for debris-covered Imja Glacier, where a lake developed after the 1960s. This was not true for adjacent glaciers with similar characteristics but without any, or with very small, glacier lakes (Bajracharja and Mool, 2009; Bolch and others, 2011). This is observed in this study where moraine-dam lakes have substantially influenced glacier retreat, with an overall deglaciation of 0.20 km2 per glacier, between 1989 and 2010. Likewise, in the Himalaya, studies have reported a much larger lake growth in the eastern Himalaya than in the central and drier northwest regions (Xin and others, 2008; Gardelle and others, 2011).

Climate considerations

The mean annual air temperature (MAAT) record at Gangtok station (located at 1812 m a.s.l. (Fig. 1); has 24 years of continuous and complete long-term records) showed a rise of 1.03°C (0.04°C a−1) between 1987 and 2011. The summer (April–October) minimum temperature also recorded a rise of ∼1°C, but a significant increase of 2.04°C was observed in the winter (November–March) minimum temperature. The winter maximum temperature rose by 1.13°C as opposed to a 0.17°C rise in the summer maximum temperature, suggesting warmer winters (Fig. 11). The precipitation had no clear trend due to the high variability in rainfall patterns, but, on average, there was a slight decrease in winter precipitation (Fig. 12).

Fig. 11. Trends in mean annual, summer and winter temperatures from 1987 to 2011 at Gangtok, Sikkim. The summer months are April–October and the winter months are November–March. The temperature trends are extracted from the mean monthly temperature records at Gangtok station.

Fig. 12. Variable trends in summer and winter rainfall between 1987 and 2011 at Gangtok. The rainfall data are extracted from the monthly total rainfall records. The summer rainfall is highly variable.

The climate observations showed a rise in MAAT and a decline in the winter precipitation. These observations are similar to weather station data observed in Nepal Himalaya (Shrestha and others, 2000). The Sikkim glaciers are predominantly monsoon-affected and are more sensitive to temperature changes than winter-accumulation-type glaciers, as the temperature increase directly reduces solid precipitation (Fujita, 2008). One reason for the loss in glacier area and the development of lakes on the Sikkim glaciers could be the rise in MAAT over recent decades and, especially, higher winter and summer minimum temperatures, leading to a prolonged melting season, accompanied by a slight decrease in winter precipitation (Figs 11 and 12). Precipitation may fall as rain on the higher regions of glaciers and reduce the accumulation area, leading to an upward shift in the 0°C isotherm line. This causes a high rate of retreat at lower elevations (Racoviteanu and others, 2008b; Bolch and others, 2012).

However, the climate record is too short-term for confirmation, as the glacier response time, which can be some decades for larger glaciers, has to be taken into account. In addition, Gangtok station and most of the regional weather stations are located below 2000 m in Sikkim and may not provide applicable comparisons to the glacier environment. This is especially true for precipitation, which is highly variable, especially in mountainous terrain, so trends can differ significantly (Bhambri and others, 2011). Thus, setting up a glacier observatory would fill the data gap.


Remote-sensing data were used to estimate the glacier area in the Tista basin for the years 1989/90, 1997, 2004/05 and 2009/10. The ‘mapping’ and ‘interpretation’ uncertainties were estimated using perimeters, spatial resolution of sensors and considering uncertainties in identifying glacier snouts. The investigation has shown an overall reduction of 6.9 ± 1.5 km2 (∼0.16% a−1) in the glacier area between 1989 and 2010. The loss in glacier area in Sikkim accelerated after 1997. Moreover, the development and expansion of lakes and supraglacial lakes on debris-covered glaciers was observed during this period. The number of ‘debris with lakes’ glaciers, the area of debris cover and the size of moraine-dammed lakes has also increased and this could further influence mass loss. Analysis of the climate data indicates that the loss in glacier area is possibly caused by increasing temperatures. However, the reaction of glaciers to climate change is not uniform, as ‘clean’, ‘debris’ and ‘debris with lakes’ glaciers showed variations in area loss, though located in almost the same altitude ranges. Most ‘debris-covered glaciers’ in Sikkim have moraine-dam lakes at their termini, and showed a greater loss in area than ‘clean’ glaciers. Though some ‘debris-covered glaciers’ showed stable fronts, the gradual development and coalescence of supraglacial lakes led to the formation of moraine-dam lakes at the terminus. This study suggests that ‘debris cover’ on glaciers can enhance the development of glacial lakes. As a consequence, the retreat of debris-covered glaciers associated with lakes is clearly higher than that of debris-free glaciers.


We gratefully acknowledge funding from the Divecha Centre for Climate Change, Indian Institute of Science, Bangalore. We thank two anonymous reviewers for helpful comments. T. Bolch acknowledges funding through Deutsche Forschungsgemeinschaft (DFG). The mean monthly temperature and precipitation data from 1987 to 2011, of Gangtok station were procured from the India Meteorological Department.


IM Bahuguna , AV Kulkarni , ML Arrawatia and DG Shresta (2001) Glacier atlas of Tista basin (Sikkim Himalaya). (SAC/RESA/MWRG-GLI/SN/16/2001) Space Application Centre, Ahmedabad
SR Bajracharya and P Mool (2009) Glaciers, glacial lakes and glacial lake outburst floods in the Mount Everest region, Nepal. Ann. Glaciol., 50(53), 8186 (doi: 10.3189/172756410790595895)
S Basnett and AV eds. Kulkarni (2012) Monitoring the seasonal snow cover in Sikkim Himalaya using remote sensing techniques. In ML Arrawatia and S Tambe eds. Climate change in Sikkim: patterns, impacts and initiatives. Information and Public Relations Department, Government of Sikkim, Gangtok, 6980
S Basnett , AV Kulkarni , BP Rathore , ML Arrawatia and DG Shrestha (2011) Monitoring Sikkim Himalayan cryosphere. In B Shrestha and B Bajracharya eds. International Symposium on Benefiting from Earth Observation: Bridging the Data Gap for Adaptation to Climate Change in the Hindu Kush–Himalayan Region, 4– 6 October 2010, Kathmandu, Nepal. International Centre for Integrated Mountain Development, Kathmandu
DI Benn and 9 others (2012) Response of debris-covered glaciers in the Mount Everest region to recent warming, and implications for outburst flood hazards. Earth-Sci. Rev., 114(1–2), 156174 (doi: 10.1016/j.earscirev.2012.03.008)
MR Bennett and NF Glasser (2009) Glacial geology: ice sheets and landforms, 2nd edn. Wiley-Blackwell, Chichester
R Bhambri and T Bolch (2009) Glacier mapping: a review with special reference to the Indian Himalayas. Progr. Phys. Geogr., 33(5), 672704 (doi: 10.1177/0309133309348112)
R Bhambri , T Bolch , RK Chaujar and SC Kulshreshtha (2011) Glacier changes in the Garwal Himalaya, India, from 1968 to 2006 based on remote sensing. J. Glaciol., 57(203), 543556 (doi: 10.3189/002214311796905604)
R Bhambri , T Bolch , P Kawishwar , DP Dobhal , D Srivastrava and B Pratap (2013) Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram. Cryosphere, 7, 13851398 (doi: 10.5194/tc-7-1385-2013)
T Bolch , M Buchroithner , T Pieczonka and A Kunert (2008) Planimetric and volumetric glacier changes in the Khumbu Himal, Nepal, since 1962 using Corona, Landsat TM and ASTER data. J. Glaciol., 54(187), 592600 (doi: 10.3189/002214308786570782)
T Bolch , B Menounos and R Wheate (2010) Landsat-based inventory of glaciers in western Canada, 1985–2005. Remote Sens. Environ., 114(1), 127137 (doi: 10.1016/j.rse.2009.08.015)
T Bolch , T Pieczonka and DI Benn (2011) Multi-decadal mass loss of glaciers in the Everest area (Nepal Himalaya) derived from stereo imagery. Cryosphere, 5(2), 349358 (doi: 10.5194/tc-5-349-2011)
T Bolch and 10 others (2012) The state and fate of Himalayan glaciers. Science, 336(6079), 310314 (doi: 10.1126/science.1215828)
B Bookhagen and DW Burbank (2010) Toward a complete Himalayan hydrological budget: spatiotemporal distribution of snowmelt and rainfall and their impact on river discharge. J. Geophys. Res., 115(F3), F03019 (doi: 10.1029/2009JF001426)
CM DeBeer and MJ Sharp (2007) Recent changes in glacier area and volume within the southern Canadian Cordillera. Ann. Glaciol., 46, 215221 (doi: 10.3189/172756407782871710)
R Finsterwalder (1935) On the map of the Zemu Glacier. Himalayan J., 7, 125138
K Fujita (2008) Effect of precipitation seasonality on climatic sensitivity of glacier mass balance. Earth Planet. Sci. Lett., 276(1–2), 1419 (doi: 10.1016/j.epsl.2008.08.028)
K Fujita , A Sakai , T Nuimura , S Yamaguchi and RR Sharma (2009) Recent changes in Imja Glacial Lake and its damming moraine in the Nepal Himalaya revealed by in situ surveys and multitemporal ASTER imagery. Environ. Res. Lett., 4, 045205 (doi: 10.1088/1748-9326/4/4/045205)
J Gardelle , Y Arnaud and E Berthier (2011) Contrasted evolution of glacial lakes along the Hindu Kush Himalaya mountain range between 1990 and 2009. Global Planet. Change, 75(1–2), 4755 (doi: 10.1016/j.gloplacha.2010.10.003)
FD Granshaw and AG Fountain (2006) Glacier change (1958–1998) in the North Cascades National Park Complex, Washington, USA. J. Glaciol., 52(177), 251256 (doi: 10.3189/172756506781828782)
J Gulley and DI Benn (2007) Structural control of englacial drainage systems in Himalayan debris-covered glaciers. J. Glaciol., 53(182), 399412 (doi: 10.3189/002214307783258378)
WW Immerzeel , LPH Van Beeke and MFP Bierkens (2010) Climate change will affect the Asian water towers. Science, 328(5984), 13821385 (doi: 10.1126/science.1183188)
A Kääb , E Berthier , C Nuth , J Gardelle and Y Arnaud (2012) Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature, 488(7412), 495498 (doi: 10.1038/nature11324)
U Kamp , M Byrne and T Bolch (2011) Glacier fluctuations between 1975 and 2008 in the Greater Himalaya Range of Zanskar, southern Ladakh. J. Mt. Sci., 8(3), 374389 (doi: 10.1007/s11629-011-2007-9)
G Kaser , M Grosshauser and B Marzeion (2010) Contribution potential of glaciers to water availability in different climate regimes. Proc. Natl Acad. Sci. USA (PNAS), 107(47), 20 22320 227 (doi: 10.1073/pnas.1008162107)
J Komori (2008) Recent expansions of glacial lakes in the Bhutan Himalayas. Quat. Int., 184(1), 177186 (doi: 10.1016/j.quaint.2007.09.012)
AP Krishna (2005) Snow and glacier cover assessment in the high mountains of Sikkim Himalaya. Hydrol. Process., 19(12), 23752383 (doi: 10.1002/hyp.5890)
AV Kulkarni (1991) Glacier inventory in Himachal Pradesh using satellite data. J. Indian Soc. Remote Sens., 19(3), 195203
AV Kulkarni and S Alex (2003) Estimation of recent glacial variations in Baspa basin using remote sensing technique. J. Indian Soc. Remote Sens., 31(2), 8190 (doi: 10.1007/BF03030775)
AV Kulkarni , SK Singh , P Mathur and VD Mishra (2006) Algorithm to monitor snow cover using AWiFS data of RESOURCESAT-1 for the Himalayan region. Int.J. Remote Sens., 27(12), 24492457
AV Kulkarni , BP Rathore , SK Singh and IM Bahuguna (2011) Understanding changes in the Himalayan cryosphere using remote sensing techniques. Int. J. Remote Sens., 32(3), 601615 (doi: 10.1080/01431161.2010.517802)
MF Meier (1961) Distribution and variations of glaciers in the United States exclusive of Alaska. IASH Publ. 54 (General Assembly of Helsinki 1960 – Snow and Ice), 420429
HC Nainwal , BDS Negi , M Chaudhary , KS Sajwan and A Gaurav (2008) Temporal changes in rate of recession: evidence from Satopanth and Bhagirath Kharak glaciers, Uttarakhand, using Total Station Survey. Current Sci., 94(5), 653660
VN Nijampurkar (1985) Basal flow rate of Changme–Khangpu glacier, Sikkim Himalaya based on 32Si and 210Pb chronology. Proc. Indian Acad. Sci., Sect. A, 94(2), 8389 (doi: 10.1007/BF02871940)
VN Nijampurkar , N Bhandari , DV Borole and U Bhattacharya (1985) Radiometric chronology of Changme-Khangpu Glacier, Sikkim. J. Glaciol., 31(107), 2833
T Nuimura , K Fujita , S Yamaguchi and RR Sharma (2012) Elevation changes of glaciers revealed by multitemporal digital elevation models calibrated by GPS survey in the Khumbu region, Nepal Himalayas, 1992–2008. J. Glaciol., 58(210), 648656 (doi: 10.3189/2012JoG11J061)
F Paul , A Kääb , M Maisch , T Kellenberger and W Haeberli (2004) Rapid disintegration of Alpine glaciers observed with satellite data. Geophys. Res. Lett., 31(21), L21402 (doi: 10.1029/2004GL020816)
F Paul and 19 others (2013) On the accuracy of glacier outlines derived from remote-sensing data. Ann. Glaciol., 54(63 Pt 1), 171182 (doi: 10.3189/2013AoG63A296)
AE Racoviteanu , Y Arnaud , MW Williams and J Ordoñez (2008a) Decadal changes in glacier parameters in the Cordillera Blanca, Peru, derived from remote sensing. J. Glaciol., 54(186), 499510
AE Racoviteanu , MW Williams and RG Barry (2008b) Optical remote sensing of glacier characteristics: a review with focus on the Himalaya. Sensors, 8(5), 33553383 (doi: 10.3390/s8053355)
VK Raina (2009) Himalayan glaciers: a state-of-art review of glacial studies, glacial retreat and climate change. (SPPI Reprint Series, MoEF Discussion Paper) Ministry of Environment & Forests, Government of India/G.B. Pant Institute of Himalayan Environment and Development, New Delhi/Kosi-Katarmal, Almora
VK Raina , V Bhattachorya and S Pattnaik (1973) ‘Zemu glacier’. Rec. Geol. Surv. India, 105(2), 95106
BKG Raj , VK Kumar and SN Remya (2012) Remote sensing-based inventory of glacial lakes in Sikkim Himalaya: semi-automated approach using satellite data. Geomat. Natur. Hazards Risk, 4(3), 241253 (doi: 10.1080/19475705.2012.707153)
K Röhl (2008) Characteristics and evolution of supraglacial ponds on debris-covered Tasman Glacier, New Zealand. J. Glaciol., 54(188), 867880 (doi: 10.3189/002214308787779861)
A Sakai and K Fujita (2010) Correspondence. Formation conditions of supraglacial lakes on debris-covered glaciers in the Himalaya. J. Glaciol., 56(195), 177181 (doi: 10.3189/002214310791190785)
A Sakai , N Takeuchi , K Fujita and M Nakawo (2000) Role of supraglacial ponds in the ablation process of a debris-covered glacier in the Nepal Himalayas. IAHS Publ. 264 (Symposium at Seattle 2000 – Debris-Covered Glaciers), 119130
A Sakai , K Nishimura , T Kadota and N Takeuchi (2009) Onset of calving at supraglacial lakes on debris-covered glaciers of the Nepal Himalaya. J. Glaciol., 55(193), 909917 (doi: 10.3189/002214309790152555)
D Scherler , B Bookhagen and MR Strecker (2011) Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nature Geosci., 4(3), 156159 (doi: 10.1038/ngeo1068)
AB Shrestha , CP Wake , JE Dibb and PA Mayewski (2000) Precipitation fluctuations in the Nepal Himalaya and its vicinity and relationship with some large scale climatological parameters. Int. J. Climatol., 20(3), 317327
PN Shukla , N Bhandari , VN Nijampurkar , DK Rao , VMK Puri and S Sharma (1983) Ice accumulation rate in Changme Khangpu glacier, Sikkim. Proc. Indian Acad. Sci., Sect. A, 92(3), 255260 (doi: 10.1007/BF02854594)
W Tangborn and B Rana (2000) Mass balance and runoff of the partially debris-covered Langtang Glacier, Nepal. IAHS Publ. 264 (Symposium at Seattle 2000 – Debris-Covered Glaciers), 99108
UNESCO/International Association of Scientific Hydrology (IASH) (1970) Perennial ice and snow masses: a guide for compilation and assemblage of data for a world inventory, ed. F Müller . (Technical Papers in Hydrology 1) UNESCO/IASH, Paris
TN Venkatesh , AV Kulkarni and J Srinivasan (2011) Relative effect of slope and equilibrium line altitude on the retreat of Himalayan glaciers. Cryos. Discuss., 5(5), 25712604 (doi: 10.5194/tcd-5-2571-2011)
T Watanabe , D Lamsal and JD Ives (2009) Evaluating the growth characteristics of a glacial lake and its degree of danger of outburst flooding: Imja Glacier, Khumbu Himal, Nepal. Nor. Geogr. Tidsskr., 63(4), 255267 (doi: 10.1080/00291950903368367)
W Xin , S Liu , W Guo and J Xu (2008) Assessment and simulaton of glacier lake outburst floods for Longbasaba and Pida Lakes, China. Mt. Res. Dev., 28(3), 310317 (doi: 10.1659/mrd.0894)