Quantitative measures of adiabatic lapse rate, equilibrium line altitude (ELA), and accumulation-area ratio (AAR) are important to understand the hydrological processes and conduct hydrological modeling in a highly glaciated watershed. We present a detailed analysis of temperature data from 21 climatic stations, hypsometric analyses of glacier distributions, and a method to analyze ablation gradients and runoff curves concurrently to quantify these parameters for the watersheds of the Upper Indus Basin (UIB), with 15,062 km2 of glacierized area and an elevation range of 361–8,611 m. We show that the ELA varies considerably from one watershed to another, implying a highly variable upper elevation limit up to which melting of snows and glaciers takes place throughout the basin. This is in sharp contrast to the assertions made by previous researchers. We show that the ELA is as low as 4,840 m in Astore watershed and it is as high as 6,200 m in Shyok watershed. In accordance with the variation of ELA, the AAR also varies considerably from one watershed to another. It is as low as 0.10 in Gilgit and as high as 0.65 in Zanskar watersheds. We ascribe 15–20% uncertainty to these estimates of ELA and AAR in UIB.
INTRODUCTION
In mountainous river basins, where melt water from seasonal and perennial snow and ice cover forms an integral component of the river discharge, an understanding of the hydrological processes and hydrologic modeling requires quantitative models of altitudinal variations of temperature. When such basins have significant glacierized areas, the equilibrium line altitude (ELA) is a useful parameter to depict the link between glacier mass balance and climate change. The ELA is widely used to infer present and past climatic conditions (e.g., Andrews 1975; Porter 1975, 1977) and is an important climate descriptor for a glacierized watershed's summer air temperature and winter precipitation. Another parameter of great importance in the descriptions of hydrological processes on glaciers is the accumulation-area ratio (AAR), which is the ratio of the area of accumulation zone to the total glaciated area. This ratio is necessary to perform glacier mass balance calculations (e.g., see Appendix B in Mukhopadhyay & Khan 2015). In most cases, particularly in rugged and remote terrains, these parameters are largely unknown due to difficulties in obtaining sufficient field-based observations that lead to the estimations of these parameters. In this paper, we provide the estimates of these parameters for the Upper Indus Basin (UIB).
UIB is among the most melt water-dependent river basins worldwide. The Upper Indus River originates at an elevation of about 5,166 m above sea level (asl) in the remote region of western Tibet and flows in a general northwest direction between the Zanskar, Greater Himalayas and the Karakoram Mountains all the way to the Hindu Kush Mountains. Thus, UIB straddles two very high altitude mountain ranges that contain numerous peaks with elevations greater than 6,000 m asl and several notable peaks in the range of 7,000–8,000 m asl such as Rakaposhi, Nanga Parbat, and K2 (Mount Godwin Austen). According to the elevation data derived from the Shuttle Radar Topography Mission (SRTM), launched by the National Aeronautics and Space Administration (NASA) in 2005, the elevation in the basin ranges from 361 m at the basin outlet (the Tarbela reservoir intake) to 8,569 m asl (Khan et al. 2014) with a mean elevation of 4,654 m asl. The highest point in the basin is the peak of Mount Godwin Austen (K2) at 8,611 m, the second highest peak in the world (the difference is due to errors in the SRTM data; see Mukhopadhyay & Singh 2011). The mean elevation refers to the elevation from where half of the basin area is above and the other half is below. In addition to the high altitude nature of the terrain, UIB is the abode of some remarkable glaciers such as Siachen, Masherbrum, Panmah, Baltoro, Biafo, Chogo Lungma, Batura, Hispar, and Rimo glaciers, among others, and surrounding snow covered mountain peaks and slopes. According to Bajracharya & Shrestha (2011), there are 11,413 glaciers in UIB with an area of 15,062 km2. Consequently melt water from glaciers, and perennial and seasonal snow cover, forms the dominant constituent of the river flows in UIB.
Mukhopadhyay & Khan (2015) have shown that glacial melt far outweighs snowmelt in the rivers draining the Karakoram and Zanskar ranges. According to their estimates, in the Karakoram, the annual glacial melt proportion varies from 43 to 50% whereas snowmelt varies from 27 to 31%. On the other hand, snowmelt dominates over glacial melt in the rivers draining the western Greater Himalayas and the Hindu Kush. Here, the snowmelt percentage in river discharge varies from 31 to 53%, whereas that of glacial melt ranges from 16 to 30%. In the main stem of the Upper Indus River, the snowmelt fraction in most cases is slightly greater than the glacial melt fraction. In the main stem, the snowmelt percentage ranges from 35 to 44% whereas the glacial melt percentage ranges from 25 to 36%. The Upper Indus River just upstream of the Tarbela Reservoir carries annual flows constituted of 70% melt water, of which 26% is contributed by glacial melts and 44% by snowmelts.
In this paper, we present models of altitudinal variation of temperature in UIB where the variations of elevation are remarkably great. Archer (2004) has examined this problem. However, the data used in his study are mostly from the climatic stations that are located at valley floor elevations. Since then, further data have become available not only from these stations but also from several stations placed at higher elevations more recently. These datasets provide an opportunity to reexamine this problem, to better define the variation of temperature as a function of space, season, and altitude within this river basin, which in turn may help to better refine the snowmelt models of UIB developed previously (e.g., Mukhopadhyay & Dutta 2010; Tahir et al. 2011).
The ELA is best determined by careful measurements of specific balance quantities at many points on the glacier surface, so that reasonably reliable isolines of zero balance may be drawn (Braithwaite & Müller 1980). Such measurements have been carried out on relatively few glaciers in the Karakoram Range and, even then, for only a few years (e.g., Wake 1989; Hewitt 2005). In this paper, we develop an innovative approach using hydrological methods to estimate the ELA within various watersheds of UIB. We also provide the estimations of the AAR in these watersheds from estimates of glacierized areas and the ELA.
DATA
We have collected temperature data from the existing climatic stations maintained and operated by two official authorities of Pakistan, one in India and another in China (Table 1). Therefore, these data are from reliable sources and are of good quality as described below.
List of climatic stations from where data are collected and used in this investigation
. | . | . | Coordinates . | . | . | . | |
---|---|---|---|---|---|---|---|
Station . | Watershed . | Altitude (m) . | Latitude . | Longitude . | Period of record . | No. of years . | Source* . |
Astore | Astore | 2,168 | 35.3665 | 74.8647 | 1954–2010 | 57 | 1 |
Rama | Astore | 3,179 | 35.4548 | 74.7765 | 1999–2010 | 12 | 2 |
Rattu | Astore | 2,718 | 35.1615 | 74.7852 | 1999–2010 | 12 | 2 |
Burzil | Outside | 4,239 | 34.8994 | 75.0791 | 1999–2010 | 12 | 2 |
Gilgit | Gilgit | 1,460 | 35.9206 | 74.3271 | 1951–2010 | 60 | 1 |
Ushkore | Gilgit | 3,051 | 36.0271 | 73.4148 | 1999–2010 | 12 | 2 |
Yasin | Gilgit | 3,280 | 36.4505 | 73.2942 | 1999–2010 | 12 | 2 |
Gupis | Gilgit | 2,156 | 36.1789 | 73.4386 | 1955–2010 | 56 | 1 |
Shandur | Outside | 3,712 | 36.0878 | 72.5474 | 1994–2010 | 17 | 2 |
Zani Pass | Outside | 3,839 | 36.3518 | 72.1691 | 1994–2010 | 17 | 2 |
Hunza | Hunza | 2,374 | 36.3220 | 74.6461 | 2007–2010 | 4 | 1 |
Naltar | Hunza | 2,898 | 36.1679 | 74.1753 | 1999–2012 | 12 | 2 |
Khunjerab | Hunza | 4,630 | 36.8119 | 75.3321 | 1999–2013 | 12 | 2 |
Ziarat | Hunza | 3,020 | 36.8290 | 74.4178 | 1999–2014 | 12 | 2 |
Shigar | Shigar | 2,367 | 35.6982 | 75.4814 | 1996–2010 | 15 | 2 |
Hushey | Shyok | 3,075 | 35.3424 | 76.1395 | 1994–2010 | 17 | 2 |
Deosai | Upper Indus | 4,149 | 35.003756 | 75.592112 | 1994–2010 | 17 | 2 |
Bunji | Upper Indus | 1,372 | 35.6462 | 74.6292 | 1953–2010 | 58 | 1 |
Skardu | Upper Indus | 2,317 | 35.2864 | 75.5630 | 1952–2010 | 59 | 1 |
Leh | Upper Indus | 3,505 | 34.164 | 77.587 | 1876–1990 | 106 | 3 |
Demchok | Upper Indus | 4,278 | 32.500 | 80.080 | 1979–2010 | 31 | 4 |
. | . | . | Coordinates . | . | . | . | |
---|---|---|---|---|---|---|---|
Station . | Watershed . | Altitude (m) . | Latitude . | Longitude . | Period of record . | No. of years . | Source* . |
Astore | Astore | 2,168 | 35.3665 | 74.8647 | 1954–2010 | 57 | 1 |
Rama | Astore | 3,179 | 35.4548 | 74.7765 | 1999–2010 | 12 | 2 |
Rattu | Astore | 2,718 | 35.1615 | 74.7852 | 1999–2010 | 12 | 2 |
Burzil | Outside | 4,239 | 34.8994 | 75.0791 | 1999–2010 | 12 | 2 |
Gilgit | Gilgit | 1,460 | 35.9206 | 74.3271 | 1951–2010 | 60 | 1 |
Ushkore | Gilgit | 3,051 | 36.0271 | 73.4148 | 1999–2010 | 12 | 2 |
Yasin | Gilgit | 3,280 | 36.4505 | 73.2942 | 1999–2010 | 12 | 2 |
Gupis | Gilgit | 2,156 | 36.1789 | 73.4386 | 1955–2010 | 56 | 1 |
Shandur | Outside | 3,712 | 36.0878 | 72.5474 | 1994–2010 | 17 | 2 |
Zani Pass | Outside | 3,839 | 36.3518 | 72.1691 | 1994–2010 | 17 | 2 |
Hunza | Hunza | 2,374 | 36.3220 | 74.6461 | 2007–2010 | 4 | 1 |
Naltar | Hunza | 2,898 | 36.1679 | 74.1753 | 1999–2012 | 12 | 2 |
Khunjerab | Hunza | 4,630 | 36.8119 | 75.3321 | 1999–2013 | 12 | 2 |
Ziarat | Hunza | 3,020 | 36.8290 | 74.4178 | 1999–2014 | 12 | 2 |
Shigar | Shigar | 2,367 | 35.6982 | 75.4814 | 1996–2010 | 15 | 2 |
Hushey | Shyok | 3,075 | 35.3424 | 76.1395 | 1994–2010 | 17 | 2 |
Deosai | Upper Indus | 4,149 | 35.003756 | 75.592112 | 1994–2010 | 17 | 2 |
Bunji | Upper Indus | 1,372 | 35.6462 | 74.6292 | 1953–2010 | 58 | 1 |
Skardu | Upper Indus | 2,317 | 35.2864 | 75.5630 | 1952–2010 | 59 | 1 |
Leh | Upper Indus | 3,505 | 34.164 | 77.587 | 1876–1990 | 106 | 3 |
Demchok | Upper Indus | 4,278 | 32.500 | 80.080 | 1979–2010 | 31 | 4 |
*Sources: 1 = Pakistan Meteorological Department; 2 = Water and Power Development Authority, Pakistan; 3 = Indian Meteorological Department (IMD) (at Leh, data for 1970–1977 are missing); 4 = China Meteorological Administration. All stations' daily data (except Leh and Demchok) have been obtained, and a monthly average has been estimated from these datasets.
We have collected temperature data from 21 climatic stations (Table 1). The valley-based stations within the Pakistani part of the basin are maintained by the Pakistan Meteorological Department (PMD) and have a long period of record (>55 years), and the high altitude stations are maintained by the Water and Power Development Authority of Pakistan (WAPDA) for its snow and ice project, with a shorter period of record (Table 1). We have obtained daily maximum and minimum temperatures from PMD and WAPDA. PMD collects climatic data in accordance with the guidelines prescribed by the World Meteorological Organization (WMO). WAPDA has auto-recording climatic stations, and those data are collected in accordance with WMO guidelines too. In addition, we have collected temperature data from the Indian Meteorological Department for one station (Leh) within the Indian part of the basin and also from the China Meteorological Administration for one station (Demchok) within the Tibetan part of UIB.
The other datasets we have used include data related to topography, snow covered areas, glacier cover, and river flows. These data are described in Khan et al. (2014), Mukhopadhyay (2012, 2013), and Mukhopadhyay & Khan (2014a, 2014b, 2015).
THE MAJOR WATERSHEDS OF UIB
Major rivers and watersheds of UIB with locations of the climatic stations from where data have been collected and used in the present study.
Major rivers and watersheds of UIB with locations of the climatic stations from where data have been collected and used in the present study.
Major watersheds of UIB
Watershed . | Watershed area (km2) . | Glacial cover RGI-3.0 (km2)1 . | Glacier extents (range in m) . |
---|---|---|---|
Shyok | 33,470 | 7,696 ± 770 | 3,221–7,803 |
Hunza | 13,734 | 3,843 ± 384 | 2,299–7,809 |
Shigar | 7,040 | 2,121 ± 212 | 2,716–8,569 |
Gilgit | 12,713 | 836.3 ± 84 | 2,750–6,992 |
Zanskar | 14,861 | 1,189 ± 119 | 3,927–6,509 |
Shingo | 10,502 | 733 ± 73 | 3,579–7,103 |
Astore | 3,988 | 541 ± 54 | 2,897–8,020 |
Kharmong (UIB)* | 44,667 | 601 ± 60 | 4,438–6,507 |
Tarbela (UIB)** | 26,053 | 776± 78 | 2,602–8,069 |
Watershed . | Watershed area (km2) . | Glacial cover RGI-3.0 (km2)1 . | Glacier extents (range in m) . |
---|---|---|---|
Shyok | 33,470 | 7,696 ± 770 | 3,221–7,803 |
Hunza | 13,734 | 3,843 ± 384 | 2,299–7,809 |
Shigar | 7,040 | 2,121 ± 212 | 2,716–8,569 |
Gilgit | 12,713 | 836.3 ± 84 | 2,750–6,992 |
Zanskar | 14,861 | 1,189 ± 119 | 3,927–6,509 |
Shingo | 10,502 | 733 ± 73 | 3,579–7,103 |
Astore | 3,988 | 541 ± 54 | 2,897–8,020 |
Kharmong (UIB)* | 44,667 | 601 ± 60 | 4,438–6,507 |
Tarbela (UIB)** | 26,053 | 776± 78 | 2,602–8,069 |
*Kharmong (UIB): Watershed draining the headwater catchments of Upper Indus River and those along its course up to Kharmong (before its confluence with Shyok, Zanskar and Dras Nala watersheds are excluded).
**Tarbela (UIB): Watershed draining the catchments of the Upper Indus River along its course from Kharmong up to the point of inflow into the Tarbela Reservoir (all major tributary watersheds are excluded).
Data from the Randolph Glacier Inventory version 3.0, released in 2013 (http://www.glims.org/RGI/). These areas are estimated for glaciers only. The accuracy is assumed to be 90% (error bars at 10%).
The extents of glaciers in the elevation range have been derived from the SRTM DEM by overlying the RGI 3.0 glacier maps on it.
Percentages of surface areas in five elevation zones (hypsometric bands) for the nine major watersheds of UIB
. | Hypsometric bands (areas are in percentages of total areas) . | |||||
---|---|---|---|---|---|---|
Watershed . | < 1,000 m . | 1,000–2,500 m . | 2,500–3,500 m . | 3,500–5,300 m . | 5,300–6,500 m . | > 6,500 m . |
Tarbela | 0.15 | 2.89 | 10.36 | 64.57 | 21.93 | 0.1 |
Shyok | 0 | 0.84 | 3.73 | 44.73 | 50.2 | 0.5 |
Hunza | 0 | 3.64 | 16.75 | 69.9 | 5.82 | 3.89 |
Shigar | 0 | 4.26 | 8.52 | 65.34 | 19.89 | 1.99 |
Gilgit | 0 | 4.64 | 14.24 | 78.58 | 2.44 | 0.1 |
Zanskar | 0 | 0 | 1.48 | 69.98 | 28.26 | 0.28 |
Shingo | 0 | 0 | 8.57 | 87.13 | 2.66 | 1.64 |
Astore | 0 | 2.51 | 18.8 | 76.48 | 1.25 | 0.96 |
Kharmong (UIB)* | 0 | 0.04 | 2.19 | 72.77 | 24.85 | 0.15 |
. | Hypsometric bands (areas are in percentages of total areas) . | |||||
---|---|---|---|---|---|---|
Watershed . | < 1,000 m . | 1,000–2,500 m . | 2,500–3,500 m . | 3,500–5,300 m . | 5,300–6,500 m . | > 6,500 m . |
Tarbela | 0.15 | 2.89 | 10.36 | 64.57 | 21.93 | 0.1 |
Shyok | 0 | 0.84 | 3.73 | 44.73 | 50.2 | 0.5 |
Hunza | 0 | 3.64 | 16.75 | 69.9 | 5.82 | 3.89 |
Shigar | 0 | 4.26 | 8.52 | 65.34 | 19.89 | 1.99 |
Gilgit | 0 | 4.64 | 14.24 | 78.58 | 2.44 | 0.1 |
Zanskar | 0 | 0 | 1.48 | 69.98 | 28.26 | 0.28 |
Shingo | 0 | 0 | 8.57 | 87.13 | 2.66 | 1.64 |
Astore | 0 | 2.51 | 18.8 | 76.48 | 1.25 | 0.96 |
Kharmong (UIB)* | 0 | 0.04 | 2.19 | 72.77 | 24.85 | 0.15 |
*Kharmong (UIB): Watershed draining the headwater catchments of the Upper Indus River and those along its course up to Kharmong (before its confluence with Shyok, Zanskar and Shingo watersheds are excluded).
As noted above, UIB spans three mountain ranges, namely the Greater Himalayas, the Karakoram, and the Hindu Kush. According to Winiger et al. (2005), these three mountain ranges have to be divided into different climatic regions. Accordingly, three watersheds, namely Hunza, Shigar, and Shyok, cover the western, central, and eastern parts of the Karakoram; Gilgit watershed covers the Hindu Kush and four watersheds, namely Astore, Shingo, Zanskar, and Kharmong astride the western, central, west-central, and eastern sections of the Western Greater Himalayas. The remaining watershed, called Tarbela, lies along the valley between the ridges of the Karakoram, the Western Greater Himalayas, and the Hindu Kush.
ALTITUDINAL VARIATION OF TEMPERATURE

Monthly average temperatures at the climatic stations: (a) Gilgit, Gupis, Yasin, Ushkore, Shandur and Zani Pass; (b) Naltar, Ziarat, Khunjerab and Hunza; (c) Skardu, Shigar, Hushey, Leh and Demchok; and (d) Astore, Rama, Rattu, Burzil and Bunji. Elevations and period of record for the stations are given in Table 1.
Monthly average temperatures at the climatic stations: (a) Gilgit, Gupis, Yasin, Ushkore, Shandur and Zani Pass; (b) Naltar, Ziarat, Khunjerab and Hunza; (c) Skardu, Shigar, Hushey, Leh and Demchok; and (d) Astore, Rama, Rattu, Burzil and Bunji. Elevations and period of record for the stations are given in Table 1.
The pair that produces the most reasonable results is selected as the representative value for a particular watershed. Consistency in the calculated values, optimum fluctuations from day to day or from month to month, and values oscillating around the standard environmental lapse rate of 0.00649 °C m−1 are the criteria used to judge the values that are most reasonable and consistent.
Lapse rates are calculated for the individual months and for the individual watersheds within the basin. Table 4 provides an internally consistent set of values of λa for different months within the major watersheds of UIB.
Monthly average lapse rates in the individual watersheds of UIB
. | Gilgit . | Hunza . | Shigar . | Shyok . | Astore . | Shingo . | Kharmong . |
---|---|---|---|---|---|---|---|
Jan | 0.00591 | 0.00602 | 0.00641 | 0.00403 | 0.00577 | 0.00709 | 0.00840 |
Feb | 0.00663 | 0.00631 | 0.00668 | 0.00587 | 0.00619 | 0.00824 | 0.00819 |
Mar | 0.00690 | 0.00655 | 0.00716 | 0.00730 | 0.00657 | 0.00955 | 0.00899 |
Apr | 0.00663 | 0.00649 | 0.00747 | 0.00743 | 0.00658 | 0.00922 | 0.01026 |
May | 0.00618 | 0.00660 | 0.00743 | 0.00738 | 0.00677 | 0.00912 | 0.00962 |
Jun | 0.00556 | 0.00661 | 0.00771 | 0.00680 | 0.00658 | 0.00861 | 0.00908 |
Jul | 0.00605 | 0.00625 | 0.00734 | 0.00598 | 0.00675 | 0.00825 | 0.00995 |
Aug | 0.00596 | 0.00596 | 0.00746 | 0.00594 | 0.00647 | 0.00756 | 0.01017 |
Sep | 0.00538 | 0.00618 | 0.00736 | 0.00494 | 0.00619 | 0.00751 | 0.00905 |
Oct | 0.00575 | 0.00666 | 0.00787 | 0.00446 | 0.00626 | 0.00729 | 0.00986 |
Nov | 0.00561 | 0.00656 | 0.00749 | 0.00346 | 0.00617 | 0.00727 | 0.00956 |
Dec | 0.00547 | 0.00641 | 0.00689 | 0.00357 | 0.00588 | 0.00707 | 0.00933 |
. | Gilgit . | Hunza . | Shigar . | Shyok . | Astore . | Shingo . | Kharmong . |
---|---|---|---|---|---|---|---|
Jan | 0.00591 | 0.00602 | 0.00641 | 0.00403 | 0.00577 | 0.00709 | 0.00840 |
Feb | 0.00663 | 0.00631 | 0.00668 | 0.00587 | 0.00619 | 0.00824 | 0.00819 |
Mar | 0.00690 | 0.00655 | 0.00716 | 0.00730 | 0.00657 | 0.00955 | 0.00899 |
Apr | 0.00663 | 0.00649 | 0.00747 | 0.00743 | 0.00658 | 0.00922 | 0.01026 |
May | 0.00618 | 0.00660 | 0.00743 | 0.00738 | 0.00677 | 0.00912 | 0.00962 |
Jun | 0.00556 | 0.00661 | 0.00771 | 0.00680 | 0.00658 | 0.00861 | 0.00908 |
Jul | 0.00605 | 0.00625 | 0.00734 | 0.00598 | 0.00675 | 0.00825 | 0.00995 |
Aug | 0.00596 | 0.00596 | 0.00746 | 0.00594 | 0.00647 | 0.00756 | 0.01017 |
Sep | 0.00538 | 0.00618 | 0.00736 | 0.00494 | 0.00619 | 0.00751 | 0.00905 |
Oct | 0.00575 | 0.00666 | 0.00787 | 0.00446 | 0.00626 | 0.00729 | 0.00986 |
Nov | 0.00561 | 0.00656 | 0.00749 | 0.00346 | 0.00617 | 0.00727 | 0.00956 |
Dec | 0.00547 | 0.00641 | 0.00689 | 0.00357 | 0.00588 | 0.00707 | 0.00933 |
Station pairs and period of records.
Gilgit: Gilgit–Schandur (1995–2005); Hunza: Naltar–Khunjerab (1999–2010); Shigar: Shigar–Khunjerab (1999–2010).
Shyok: Skardu–Hushey (1999–2010); Astore: Astore–Burzil (1999–2010); Shingo: Skardu–Deosai (1999–2010).
Kharmong: Leh–Demchok (1978–1990).
Typically, values of λa decrease from the dry adiabatic lapse rate of 0.0098 °C m−1 depending on the humidity, and are close to the global mean lapse rate of 0.00649 °C m−1, also known as the environmental lapse rate. However, in mountainous terrains, λa can show great diurnal and seasonal variations and can depart significantly from the standard environmental lapse rate. For example, Archer (2003) has noted that within the northwestern part of UIB, the value of λa ranges from 0.0065–0.0075 °C m−1. Jain et al. (2008) observed λa to vary from 0.0060–0.0074 °C m−1 in the Sutlej River basin of the Western Himalayas. Singh (1991) observed that in the Western Himalayas, calculated values of λa also depend on the elevations of the stations from which measured temperatures are used in the estimation of λa. He observed a decrease in values of λa with increasing station elevations. Kattel et al. (2015) show that the annual cycle of the λa on the northern and southern slopes of the Himalayas exhibits distinct seasonal patterns controlled by the elevation-dependent moisture content of the atmosphere and solar radiation.
It should be noted at this point that lapse rates may not remain constant with time, given the evidence for elevation-dependent warming (EDW) with more rapid rise, historical and predicted, at higher elevations (Qin et al. 2009; Mountain Research Initiative (MRI) EDW Working Group 2015). MRI (2015) notes, for example, that over the past 20 years temperatures above 4,000 m asl have warmed nearly 75% faster than temperatures in areas below 2,000 m asl. However, Qin et al. (2009) note that ‘it is found… that the warming rate increases with respect to elevation below roughly 5,000 m ASL and such an increase tendency disappears for higher elevations (>5,000 m)’. This variation in tendency with respect to EDW may well also be reflected in lapse rates. More rapid increases in temperature at higher elevations, both historic and projected, have led or will lead to a decrease in lapse rates.
Our estimations indicate that lapse rates are lower in the Shyok watershed compared to those in the Shigar and Hunza watersheds. This is most likely, as noted above, due to a general decrease in elevation from eastern to western Karakoram. For this reason, these lapse rates, for the purpose of any hydrologic modeling, should be used relative to the elevations of the base stations used in their derivations. The highest lapse rates are in the eastern watersheds of the basin over the Western Greater Himalayas, likely due to greater aridity in this part of the basin compared to the northwestern parts (see Immerzeel et al. 2015). The lower autumn/winter lapse rate in the Shyok watershed compared to the Kharmong watershed is possibly due to the fact that the Shyok watershed straddles eastern Karakoram, where the general slope is south facing, whereas the Kharmong watershed spans the western Higher Himalayas, where the general slope is north facing. Therefore, during autumn and winter, altitudinal variation of temperature in the western Himalayas is greater than that in the eastern Karakoram. Kattel et al. (2015) have shown the control of the northern and southern slopes of the Himalayas on temperature lapse rates.
In general, within a watershed, the lapse rate increases from January to April/May and then remains fairly constant during the summer months, and then starts to decline from October with the commencement of autumn and continues to decline in the winter. In other words, the greatest altitudinal variation in temperature occurs during the snow and glacier melting season. Immerzeel et al. (2015) show the spatial variation of precipitation over UIB. According to their analysis, the highest precipitation is found along the monsoon-influenced southern Himalayan arc, while in the north of the Himalayan range the precipitation decreases quickly towards the vast dry area in the northeastern part of UIB (Shyok watershed). On the other hand, the northwestern part of the basin (the Hunza and Shigar watersheds) receives considerable amounts of precipitation at high altitude. Thayyen & Dimri (2014) demonstrate that manifestations of the presence or absence of moisture is the single most important factor determining the temperature distribution along the higher Himalayan slopes, driven by the orographic forcing. These authors also suggest that ‘Diverse slope environmental lapse rate of temperature for different altitudinal sections of monsoon and cold-arid systems of the Himalaya amply demonstrate that the use of environmental lapse rate or any other temperature lapse rate values arbitrarily for extrapolating the temperature to the Himalayan cryospheric systems is not appropriate’.

Monthly FLA or 0 °C isotherm within individual watersheds of UIB for: (a) Gilgit, Hunza, Shigar and Shyok; and (b) Astore, Shongo and Kharmong.
Monthly FLA or 0 °C isotherm within individual watersheds of UIB for: (a) Gilgit, Hunza, Shigar and Shyok; and (b) Astore, Shongo and Kharmong.
ESTIMATION OF ELAs
Schematic illustration of the concept of ablation gradient in relation to ELA and altitudinal variation of accumulation.
Schematic illustration of the concept of ablation gradient in relation to ELA and altitudinal variation of accumulation.
The ablation gradient can be measured from detailed field studies. Unfortunately, such studies are lacking in UIB. The only glaciers for which the ablation gradient has been established from field investigation are the Biafo glacier (Hewitt et al. 1989) and the Baltoro glacier (Mihalcea et al. 2006). According to the results of the study conducted by Hewitt et al. (1989), an estimate of , as noted by Yu et al. (2013). In the present study, we test three values of αg: 1 m : 0.0025 m; 1 m : 0.005 m, and 1 m : 0.01 m since it is largely an unknown parameter.
Highly generalized diagram showing zones of origin of two distinct types of melt water (M1 and M2) in UIB (developed from information given in Mukhopadhyay & Khan 2014a, 2014b). M1 melt water originates from seasonal snows below the zone of ablation of glacierized surface (approximately below 3,500/3,000 m). M2 melt water originates from both seasonal snows and glaciers from the elevation range that straddles the zone of ablation of the glacierized surface.
Highly generalized diagram showing zones of origin of two distinct types of melt water (M1 and M2) in UIB (developed from information given in Mukhopadhyay & Khan 2014a, 2014b). M1 melt water originates from seasonal snows below the zone of ablation of glacierized surface (approximately below 3,500/3,000 m). M2 melt water originates from both seasonal snows and glaciers from the elevation range that straddles the zone of ablation of the glacierized surface.
Separated hydrographs showing M1, M2 and base flow components and total annual hydrographs developed from the long-term averages of the flow records at the gauging stations near the outlet of the watersheds (for locations of the gauging stations and period of records, see Mukhopadhyay & Khan 2014a). The M2 peaks of July and August are used to calculate the runoff gradients.
Separated hydrographs showing M1, M2 and base flow components and total annual hydrographs developed from the long-term averages of the flow records at the gauging stations near the outlet of the watersheds (for locations of the gauging stations and period of records, see Mukhopadhyay & Khan 2014a). The M2 peaks of July and August are used to calculate the runoff gradients.
Ablation gradients and runoff gradients for the Astore watershed, illustrating the derivation of the ELA of the watershed.
Ablation gradients and runoff gradients for the Astore watershed, illustrating the derivation of the ELA of the watershed.
Hypsometric data for distribution of glacierized areas within the individual watersheds of UIB. Note that the cumulative glacierized area increases with elevation.
Hypsometric data for distribution of glacierized areas within the individual watersheds of UIB. Note that the cumulative glacierized area increases with elevation.
Figure 8 shows the ablation curve calculated by assuming a value of as discussed above, and the runoff curve calculated for the M2 zone of the Astore watershed, as an example. Let at any elevation the amount of runoff generated be denoted by
and the amount of melt water produced be denoted by
. Mostly due to loss in the ground and due to evaporation,
. Let an ablation curve and the runoff curve intersect at elevation
(see Figure 8). The following relationships are noted:
for, , which is not possible, but for
, and at
.


We have calculated the ablation and runoff curves for each of the major watersheds of UIB (Figure 1) to derive the ELA within the respective watersheds for both of the months of July and August. Table 5 provides the results of the calculations. The lower and upper limits of the ELA for a month correspond to the lower and upper limits of the values of αg selected (which are 1:0.0025 m and 1:0.01 m, respectively). From the calculations of the ELA for the months of July and August, we take the averages of the lower and upper limits of the ELA as given in Table 5.
Estimated ranges of the ELA for the months of July and August within the individual watersheds of UIB
. | ELA July (m) . | ELA August (m) . | Average ELA (m) . | |||
---|---|---|---|---|---|---|
Watersheds . | Lower limit . | Upper limit . | Lower limit . | Upper limit . | Lower limit . | Upper limit . |
Gilgit | 5,264 | 5,675 | 5,380 | 5,740 | 5,322 | 5,707 |
Hunza | 5,280 | 5,480 | 5,280 | 5,480 | 5,280 | 5,480 |
Shigar | 5,080 | 5,330 | 5,080 | 5,330 | 5,080 | 5,330 |
Shyok | 6,080 | 6,200 | 6,080 | 6,200 | 6,080 | 6,200 |
Kharmong | 5,280 | 5,510 | 5,310 | 5,520 | 5,295 | 5,515 |
Zanskar | 5,250 | 5,500 | 5,180 | 5,480 | 5,215 | 5,490 |
Shingo | 5,020 | 5,170 | 5,020 | 5,170 | 5,020 | 5,170 |
Astore | 4,780 | 5,100 | 4,900 | 5,140 | 4,840 | 5,120 |
. | ELA July (m) . | ELA August (m) . | Average ELA (m) . | |||
---|---|---|---|---|---|---|
Watersheds . | Lower limit . | Upper limit . | Lower limit . | Upper limit . | Lower limit . | Upper limit . |
Gilgit | 5,264 | 5,675 | 5,380 | 5,740 | 5,322 | 5,707 |
Hunza | 5,280 | 5,480 | 5,280 | 5,480 | 5,280 | 5,480 |
Shigar | 5,080 | 5,330 | 5,080 | 5,330 | 5,080 | 5,330 |
Shyok | 6,080 | 6,200 | 6,080 | 6,200 | 6,080 | 6,200 |
Kharmong | 5,280 | 5,510 | 5,310 | 5,520 | 5,295 | 5,515 |
Zanskar | 5,250 | 5,500 | 5,180 | 5,480 | 5,215 | 5,490 |
Shingo | 5,020 | 5,170 | 5,020 | 5,170 | 5,020 | 5,170 |
Astore | 4,780 | 5,100 | 4,900 | 5,140 | 4,840 | 5,120 |
Different researchers have also estimated various values of ELA in this terrain (Table 6) using other methods. For example, Gardelle et al. (2013) used approximately 2000 Landsat images to estimate the ELA in the Karakoram and Hindu Kush mountains. Our calculations are in general agreement with the values listed in Table 6. However, our estimates of ELA for each of the individual watersheds within UIB show that with the exception of the Shyok watershed, the ELA is in the general range of 5,000–5,500 m. In the case of the Shyok watershed, the ELA is in the range of 6,000–6,200 m. This higher range of ELA in this watershed, compared to the other watersheds of the basin, can be ascribed to the greater aridity of this watershed compared to the other watersheds of the basin (Immerzeel et al. 2015, Figure 5).
ELA in different parts of UIB as estimated variously
Region . | ELA (m) . | Source . |
---|---|---|
Karakoram (central and western) | 4,200–4,800 | Young & Hewitt (1993) |
Karakoram–southern slopes | 4,600–5,300 | Owen & Benn (2005) |
Karakoram | 4,750–5,310 | Gardelle et al. (2013) |
Karakoram (central) | 5,200 | Mayer et al. (2010) |
Hindu Kush | 4,890–5,210 | Gardelle et al. (2013) |
Nanga Parbat, Western Himalayas | 4,500–4,700 | Owen & Benn (2005) |
Zanskar Range, Western Himalayas | 5,500–5,800 | Osmaston (1994) |
Region . | ELA (m) . | Source . |
---|---|---|
Karakoram (central and western) | 4,200–4,800 | Young & Hewitt (1993) |
Karakoram–southern slopes | 4,600–5,300 | Owen & Benn (2005) |
Karakoram | 4,750–5,310 | Gardelle et al. (2013) |
Karakoram (central) | 5,200 | Mayer et al. (2010) |
Hindu Kush | 4,890–5,210 | Gardelle et al. (2013) |
Nanga Parbat, Western Himalayas | 4,500–4,700 | Owen & Benn (2005) |
Zanskar Range, Western Himalayas | 5,500–5,800 | Osmaston (1994) |
ESTIMATION OF AARs
We have extracted AARs for various watersheds, using lower and upper limits of ELAs, RGI v 3.0 glacier inventory, and the SRTM DEM. First, we have selected the lower limit of the ELA in a watershed and estimated the glacier area above it, using RGI and DEM data, and then divided it by the inventory-based total glacierized area of the same watershed. This computation has resulted in defining the upper limit of AAR. Similarly, we have used the upper limit of the ELA to estimate the lower limit of the AAR. Table 7 gives the list of AAR values within the individual watersheds.
Estimated AAR in various watersheds of UIB
. | AAR . | |
---|---|---|
Watersheds . | Lower limit . | Upper limit . |
Gilgit | 0.04 | 0.10 |
Hunza | 0.27 | 0.40 |
Shigar | 0.30 | 0.43 |
Shyok | 0.14 | 0.16 |
Kharmong | 0.38 | 0.60 |
Zanskar | 0.30 | 0.65 |
Shingo | 0.40 | 0.50 |
Astore | 0.33 | 0.40 |
. | AAR . | |
---|---|---|
Watersheds . | Lower limit . | Upper limit . |
Gilgit | 0.04 | 0.10 |
Hunza | 0.27 | 0.40 |
Shigar | 0.30 | 0.43 |
Shyok | 0.14 | 0.16 |
Kharmong | 0.38 | 0.60 |
Zanskar | 0.30 | 0.65 |
Shingo | 0.40 | 0.50 |
Astore | 0.33 | 0.40 |
Note: The lower limit of the AAR in the Gilgit watershed most likely has a much larger uncertainty than that present in the other estimates.
DISCUSSION
In this paper we have developed an innovative approach using hydrological methods to estimate the ELA for the principal watersheds of UIB. We have developed this method after establishing the FLA in the watersheds from calculations of adiabatic lapse rates and using the concepts of the ablation gradient and runoff curve. Due to the presence of intrinsic uncertainties in the input data, we ascribe 15–20% errors to the estimates provided in this study. No analytical expression can be developed linking the independent and dependent variables. Therefore, straightforward application of an error propagation formula cannot be applied to estimate these errors directly. For this reason, these error estimates are approximations based on assumed errors in the input variables that include topographic data from the SRTM, the glacierized areas from RGI database, temperature estimates, and so on.
The estimations of the ELAs and AARs, presented in this study, can be used in improved hydrological modelling and glacier mass balance studies in UIB. Due to the absence of regional ELA estimates in the UIB, some of the previous studies have adopted ELA estimates from other western Himalayan regions for glacial mass balance analysis and future forecasts of climate change impact on the Karakoram glaciers (e.g., Chaturvedi et al. 2014). However, our estimations show that Karakoram, the Hindu Kush and the Himalayas have different ELA and AAR values. Therefore, use of the estimates given in the present study can be useful for estimating uncertainties in future climate impact analysis due to variability in ELAs and AARs. Using the ELA and AAR values provided in the current study will provide more improved and accurate results in future glacier mass balance studies and hydrologic modeling to make predictions on water availability within UIB.
The ablation gradient and runoff curve have inverse relationships in an X-Z space. The ablation gradient is a theoretical curve that is mostly a rate function of temperature that increases with decreasing elevation. The runoff curve, on the other hand, is the curve that is based on the observed runoff amount, which increases with increasing elevation. Greater runoff is generated at higher elevation because of two principal reasons. First, the glacierized surface within a watershed increases with elevation (Figure 7). Secondly, even though a principal form of energy that causes snow and ice melting is atmospheric temperature, it is a poor proxy to the total energy budget on snow and glacier surfaces due to its decreasing strength with elevation. Kattel et al. (2015) have shown that a considerable amount of solar radiation at higher elevations causes the atmospheric warming in the Himalayas. Therefore, solar radiation is the prime source of energy to drive the melting process at elevations < the FLA. In the following discussion, we show that in general solar radiation increases with elevation in UIB.
Point estimates of solar radiation for the months of July and August as a function of elevation in (a) Shigar and (b) Zanskar watersheds of UIB as examples to illustrate the general increase in solar radiation with increase in elevation.
Point estimates of solar radiation for the months of July and August as a function of elevation in (a) Shigar and (b) Zanskar watersheds of UIB as examples to illustrate the general increase in solar radiation with increase in elevation.
Our calculations show that the ELA varies considerably from one watershed to another within UIB. This in turn implies that the upper elevation limit up to which melting of snows and glaciers takes place is not constant throughout the basin. This is in sharp contrast to the assertions made by previous researchers. For example, Young & Hewitt (1993) have considered 4,800 m as the upper elevation limit of the zone of melt production. On the other hand, Hewitt (2011) has considered this elevation to be 4,650 m. Yu et al. (2013) have set 5,000 m as the highest elevation from which melt water originates to supply to river runoff within UIB. We show that the upper elevation limit of production of M2 melt water can be as low as 4,840 m in the Astore watershed to as high as 6,200 m in the Shyok watershed. In the western and central Karakoram, this elevation ranges from 5,080 to 5,480 m. Thus, in general, the ELA is in the range of ∼5,000–5,500 m.
In accordance with the variation of ELA, the AAR also varies considerably from one watershed to another. This finding has a significant bearing on the determination of loss or gain of glacial storage from the slopes of the trend lines of August flows (see Mukhopadhyay & Khan 2014b; Mukhopadhyay et al. 2015). For example, Mukhopadhyay & Khan (2014b) assumed that 30–40% of the glacierized area of the Shigar watershed belongs to the zone of ablation. Our present estimates show that 57–70% of the glacierized area of this watershed is in the zone of ablation. Therefore, the loss of glacier storage in Shigar watershed is actually much higher than the estimates given by Mukhopadhyay & Khan (2014b). Our calculations indicate that in general, the AAR in UIB can range from as low as 0.10 ± 0.04 in Gilgit watershed to as high as 0.65 ± 0.07 in Zanskar watershed.
The variations in the ELA and AAR are essentially due to the spatial variations of the FLA during the summer ablation season. These variations are principally caused by variations of altitude (terrain characteristics) and differential atmospheric moisture availability in different parts of the river basin.
CONCLUSIONS
Adiabatic lapse rates, ELA, and accumulation-area ratio are three important parameters for an understanding of the hydrological processes and hydrological modeling of mountainous watersheds with a significant presence of glaciers. UIB is such a river basin, where elevation and glacierized areas are highly significant. Thus far, quantitative measures of these parameters of this river basin have been sketchy or lacking. Here we provide a detailed analysis of temperature data from 21 climatic stations, hypsometric analyses of glacier distributions, and a new hydrological method to analyze the ablation gradient and the runoff curve concurrently to quantify these parameters for the constituent watersheds of the basin.
The 0 °C isotherm or FLA in all of the watersheds is at its lowest elevation in January (1,896–2,735 m). It starts to recede to higher elevations from February and attains the maximum elevation in July in all watersheds except Astore, where the maximum elevation is attained in August. The maximum FLA ranges from 5,239 to 6,263 m. After attainment of maximum elevation in either July or August, the FLA starts to descend to lower elevations with the onset of autumn. Accordingly, there is a great spatial variation of FLA in the basin, possibly controlled by the humidity and terrain characteristics in the different parts of the basin.
Peak discharge in the rivers occurs in July and August when the ablation zone attains its maximum extent. The elevation of the upper limit of the ablation zone or ELA during July and August ranges from 4,840 to 6,200 m. Thus, there is a great altitudinal range within the basin up to which melting of snows and glaciers occurs. This is in contrast to the assumption used in several previous investigations (e.g., Young & Hewitt 1993; Yu et al. 2013) where the upper limit of the zone of ablation has been considered to be in the range of 4,800–5,000 m. The fraction of glacierized area that is subject to ablation also varies greatly from one watershed to another as, shown by the variation in the AAR. In general, the AAR ranges from 0.27 to 0.65. However, in two watersheds, namely Gilgit and Shyok, the AAR is relatively low.
ACKNOWLEDGEMENTS
The authors are immensely grateful to the WAPDA of Pakistan for making the river discharge and climatic data available for research. The authors are also thankful to the PMD for providing long-term climatic data. The authors express gratitude to Dr Ghulam Rasul, Chief Meteorologist, PMD, Islamabad, Pakistan, Mr Muhammad Riaz, Chief Meteorologist, PMD, Lahore, Mr Danial Hashmi, Senior Officer, WAPDA, Lahore, Pakistan, and Mrs Shan Aman Rana, Department of Economics, London School of Economics and Political Sciences, UK, for their help and support in data collection. The authors are also thankful to Dr Dong Jun Seo of the Department of Civil and Environmental Engineering of the University of Texas at Arlington for supporting A.K. to hold a visiting research scholar position for this collaborative research that has resulted during A.K.'s stay in Texas, USA.