• Adel, M. M., 2002: Man-made climatic changes in the Ganges basin. Int. J. Climatol., 22 , 9931016.

  • Aizen, V. B., E. M. Aizen, J. M. Melack, and J. Dozier, 1997: Climatic and hydrologic changes in the Tien Shan, central Asia. J. Climate, 10 , 13931404.

    • Search Google Scholar
    • Export Citation
  • Archer, D. R., 2003: Contrasting hydrological regimes in the Indus Basin. J. Hydrol., 274 , 198210.

  • Archer, D. R., 2004: Hydrological implications of spatial and altitudinal variation in temperature in the Upper Indus Basin. Nord. Hydrol, 35 , 3. 209222.

    • Search Google Scholar
    • Export Citation
  • Archer, D. R., and H. J. Fowler, 2004: Spatial and temporal variations in precipitation in the Upper Indus Basin, global teleconnections and hydrological implications. Hydrol. Earth Syst. Sci., 8 , 4761.

    • Search Google Scholar
    • Export Citation
  • Briffa, K. R., T. J. Osborn, F. H. Schweingruber, I. C. Harris, P. D. Jones, S. G. Shiyatov, and E. A. Vaganov, 2001: Low frequency temperature variations from northern tree ring density network. J. Geophys. Res., 106 , 29292941.

    • Search Google Scholar
    • Export Citation
  • Cook, E. R., P. J. Krusic, and P. D. Jones, 2003: Dendroclimatic signals in long tree-ring chronologies from the Himalayas of Nepal. Int. J. Climatol., 23 , 707732.

    • Search Google Scholar
    • Export Citation
  • Dai, A., K. E. Trenberth, and T. R. Karl, 1999: Effects of clouds, soil moisture, precipitation, and water vapor on diurnal temperature range. J. Climate, 12 , 24512473.

    • Search Google Scholar
    • Export Citation
  • Easterling, D. R., and Coauthors, 1997: Maximum and minimum temperature trends for the globe. Science, 277 , 364367.

  • Folland, C. K., and Coauthors, 2001a: Observed climate variability and change. Climate Change 2001: The Scientific Basis, J. T. Houghton et al., Eds., Cambridge University Press, 99–182.

    • Search Google Scholar
    • Export Citation
  • Folland, C. K., and Coauthors, 2001b: Global temperature change and its uncertainties since 1861. Geophys. Res. Lett., 28 , 26212624.

  • Fowler, H. J., and D. R. Archer, 2005: Hydro-climatological variability in the Upper Indus Basin and implications for water resources. Regional Hydrological Impacts of Climatic Change—Impact Assessment and Decision Making, T. Wagener et al., Eds., IAHS Publication 295, 131–138.

    • Search Google Scholar
    • Export Citation
  • Gallo, K. P., D. R. Easterling, and T. C. Peterson, 1996: The influence of land use/land cover on climatological values of the diurnal temperature range. J. Climate, 9 , 29412944.

    • Search Google Scholar
    • Export Citation
  • Giorgi, F., J. W. Hurrell, M. R. Marinucci, and M. Beniston, 1997: Elevation dependency of the surface climate change signal: A model study. J. Climate, 10 , 288296.

    • Search Google Scholar
    • Export Citation
  • Hasnain, S. I., 1999: Report of Working Group on Himalayan Glaciology. International Commission on Snow and Ice, 22 pp.

  • Hasnain, S. I., 2002: Himalayan glaciers meltdown: Impact on South Asian rivers. FRIEND 2002—Regional Hydrology: Bridging the Gap between Research and Practice, H. A. J. van Lanen and S. Demuth, Eds., IAHS Publication 274, 417–423.

    • Search Google Scholar
    • Export Citation
  • Hewitt, K., 1998: Glaciers receive a surge of attention in the Karakoram Himalaya. Eos, Trans. Amer. Geophys. Union, 79 , 104105.

  • Hewitt, K., 2005: The Karakoram anomaly? Glacier expansion and the “elevation effect,” Karakoram, Himalaya. Mt. Res. Dev., 25 , 332340.

    • Search Google Scholar
    • Export Citation
  • Hewitt, K., and G. J. Young, 1993: The Snow and Ice Hydrology Project: A Pakistan–Canada research and training programme. Snow and Glacier Hydrology: Proceedings of International Symposium, Kathmandu, Nepal, 16–21 November 1992, G. J. Young, Ed., IAHS/AISH Publication 218, 49–58.

    • Search Google Scholar
    • Export Citation
  • Hingane, L. S., K. R. Kumar, and B. V. R. Murty, 1985: Long term trends of surface air temperature in India. J. Climatol., 5 , 521528.

    • Search Google Scholar
    • Export Citation
  • Hu, Z. Z., S. Yang, and R. G. Wu, 2003: Long-term climate variations in China and global warming signals. J. Geophys. Res., 108 .4614, doi:10.1029/2003JD003651.

    • Search Google Scholar
    • Export Citation
  • Jones, P. D., and A. Moberg, 2003: Hemispheric and large-scale surface air temperature variations: An extensive revision and an update to 2001. J. Climate, 16 , 206223.

    • Search Google Scholar
    • Export Citation
  • Jones, P. D., M. New, D. E. Parker, S. Martin, and I. G. Rigor, 1999: Surface air temperature and its changes over the past 150 years. Rev. Geophys., 37 , 173199.

    • Search Google Scholar
    • Export Citation
  • Karl, T. R., and Coauthors, 1993: A new perspective on recent global warming— Asymmetric trends of daily maximum and minimum temperature. Bull. Amer. Meteor. Soc., 74 , 10071023.

    • Search Google Scholar
    • Export Citation
  • Kothyari, U. C., and V. P. Singh, 1996: Rainfall and temperature trends in India. Hydrol. Processes, 10 , 357372.

  • Kumar, K. R., and L. S. Hingane, 1988: Long term variation of surface air temperature at major industrial cities of India. Climatic Change, 13 , 287307.

    • Search Google Scholar
    • Export Citation
  • Kumar, K. R., K. K. Kumar, and G. B. Pant, 1994: Diurnal asymmetry of surface temperature trends over India. Geophys. Res. Lett., 21 , 677680.

    • Search Google Scholar
    • Export Citation
  • Lal, M., and H. Harasawa, 2001: Future climate change scenarios for Asia as inferred from selected coupled atmosphere–ocean global climate models. J. Meteor. Soc. Japan, 79 , 219227.

    • Search Google Scholar
    • Export Citation
  • Lal, M., H. Harasawa, and D. Murdiyarso, 2001: Asia. Climate Change 2001: Impacts, Adaptation and Vulnerability, J. J. McCarthy et al., Eds., Cambridge University Press, 533–590.

    • Search Google Scholar
    • Export Citation
  • Liu, X., and B. Chen, 2000: Climatic warming in the Tibetan Plateau during recent decades. Int. J. Climatol., 20 , 17291742.

  • Mann, M. E., and P. D. Jones, 2003: Global surface temperatures over the past two millennia. Geophys. Res. Lett., 30 .1820, doi:10.1029/2003GL017814.

    • Search Google Scholar
    • Export Citation
  • Mastny, L., cited. 2000: Melting of earth’s ice cover reaches new high. World Watch Institute, Washington, DC. [Available online at http://www.worldwatch.org/alerts/000306.html .].

  • Miehe, G., 1992: Recent climatic changes in the Karakoram as indicated by vegetation features and geomorphological processes. Abstracts, Int. Symp. on the Karakoram and Kunlun Mountains, Kashi, China, UNESCO, 110.

  • Pan, Z., R. W. Arritt, E. S. Takle, W. J. Gutowski, C. J. Anderson, and M. Segal, 2004: Altered hydrologic feedback in a warming climate introduces a “warming hole.”. Geophys. Res. Lett., 31 .L17109, doi:10.1029/2004GL020528.

    • Search Google Scholar
    • Export Citation
  • Ramanathan, V., and Coauthors, 2005: Atmospheric brown clouds: Impacts on South Asian climate and hydrological cycle. Proc. Natl. Acad. Sci. USA, 102 , 53265333.

    • Search Google Scholar
    • Export Citation
  • Sen, K. K., R. L. Senwal, U. Rana, S. Nautiyal, R. K. Maikuri, K. S. Rao, and K. G. Saxena, 2002: Patterns and implications of land use/cover change: A case study in Pranmati watershed (Garhawal Himalaya, India). Mt. Res. Dev., 22 , 5662.

    • Search Google Scholar
    • Export Citation
  • Sharma, K. P., B. Moore, and C. J. Vorosmarty, 2000: Anthropogenic, climatic and hydrologic trends in the Kosi basin, Himalaya. Climatic Change, 47 , 141165.

    • Search Google Scholar
    • Export Citation
  • Shrestha, A. B., C. Wake, P. A. Mayewski, and J. Dibb, 1999: Maximum temperature trends in the Himalaya and its vicinity: An analysis based on temperature records from Nepal for the period 1971–94. J. Climate, 12 , 27752786.

    • Search Google Scholar
    • Export Citation
  • Shrestha, M. L., and A. B. Shrestha, 2004: Recent trends and potential climate change impacts on glacier retreat/glacier lakes in Nepal and potential adaptation measures. OECD Global Forum on Sustainable Development: Development and Climate Change, Paris, France, OECD, 5–14.

    • Search Google Scholar
    • Export Citation
  • Singh, N., and N. A. Sontakke, 2002: On climatic fluctuations and environmental changes of the Indo-Gangetic plains, India. Climatic Change, 52 , 287313.

    • Search Google Scholar
    • Export Citation
  • Singh, P., and N. Kumar, 1997: Impact assessment of climate change on the hydrological response of a snow and glacier melt runoff dominated Himalayan river. J. Hydrol., 193 , 316330.

    • Search Google Scholar
    • Export Citation
  • Singh, P., and L. Bengtsson, 2005: Impact of warmer climate on melt and evaporation for the rainfed, snowfed and glacierfed basins in the Himalayan region. J. Hydrol., 300 , 140154.

    • Search Google Scholar
    • Export Citation
  • Stone, D. A., and A. J. Weaver, 2003: Factors contributing to diurnal temperature range trends in twentieth and twenty-first century simulations of the CCCma coupled model. Climate Dyn., 20 , 435445.

    • Search Google Scholar
    • Export Citation
  • Thorne, P. W., and Coauthors, 2003: Probable causes of late twentieth century tropospheric temperature trends. Climate Dyn., 21 , 573591.

    • Search Google Scholar
    • Export Citation
  • Watson, R. T., and the Core Writing Team, Eds., 2001: Climate Change 2001: Synthesis Report. Cambridge University Press, 398 pp.

  • World Glacier Monitoring Service, cited. 2002: Glacier mass balance data 2000/01. [Available online at www.geo.unizh.ch/wgms/ .].

  • World Meteorological Organization, 1970: Guide to Hydrometeorological Practices. 2d ed. WMO No. 168, Tech. Paper 82.

  • Yadav, R. R., W-K. Park, J. Singh, and B. Dubey, 2004: Do the western Himalayas defy global warming? Geophys. Res. Lett., 31 .L17201, doi:10.1029/2004GL020201.

    • Search Google Scholar
    • Export Citation
  • Yamada, T., T. Shiraiwa, H. Iida, T. Kadota, T. Watanabe, B. Rana, Y. Ageta, and H. Fushimi, 1992: Fluctuations of the glaciers from the 1970s to 1989 in the Khumbu, Shorong and Langtang regions, Nepal Himalayas. Bull. Glacier Res., 10 , 1119.

    • Search Google Scholar
    • Export Citation
  • View in gallery
    Fig. 1.

    The Upper Indus Basin showing relief, rivers, and stations used in the analysis. Inset shows continental location of study area.

  • View in gallery
    Fig. 2.

    Trend in diurnal temperature range at Gilgit, Skardu, and Astore: (a) annual and (b) June to November.

  • View in gallery
    Fig. 3.

    Trend in winter (DJF) maximum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

  • View in gallery
    Fig. 4.

    Trend in winter (DJF) minimum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

  • View in gallery
    Fig. 5.

    Trend in summer (JJA) maximum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

  • View in gallery
    Fig. 6.

    Trend in summer (JJA) minimum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

  • View in gallery
    Fig. 7.

    Trend in winter (DJF) mean temperature at Srinagar from 1894 to 2000 with fitted linear trend shown for periods 1894–1960 and 1961–2000 and overall trend shown for the whole period (dashed line). Separate mean temperature time series are given for each of the months from December to April, with fitted linear trend lines.

  • View in gallery
    Fig. 8.

    Seasonal variations in DTR at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring.

  • View in gallery
    Fig. 9.

    Seasonal variations in minimum temperature at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring. Note that different scales are used for the y axis to highlight the differences between minimum temperatures in different seasons.

  • View in gallery
    Fig. 10.

    Seasonal variations in maximum temperature at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring. Note that different scales are used for the y axis to highlight the differences between maximum temperatures in different seasons.

  • View in gallery
    Fig. 11.

    Relationship between temperature and runoff at (left) Gilgit and the river Hunza at Dainyor Bridge and (right) Skardu and the river Shyok at Yugo. (a) Temperature trend from June to August; (b) relationship between seasonal temperature and runoff (mm) from June to August; and (c) runoff trend from June to August.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1703 664 39
PDF Downloads 1110 382 31

Conflicting Signals of Climatic Change in the Upper Indus Basin

H. J. FowlerWater Resource Systems Research Laboratory, School of Civil Engineering and Geosciences, University of Newcastle, Newcastle upon Tyne, United Kingdom

Search for other papers by H. J. Fowler in
Current site
Google Scholar
PubMed
Close
and
D. R. ArcherJBA Consulting, Skipton, North Yorkshire, United Kingdom

Search for other papers by D. R. Archer in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Temperature data for seven instrumental records in the Karakoram and Hindu Kush Mountains of the Upper Indus Basin (UIB) have been analyzed for seasonal and annual trends over the period 1961–2000 and compared with neighboring mountain regions and the Indian subcontinent. Strong contrasts are found between the behavior of winter and summer temperatures and between maximum and minimum temperatures. Winter mean and maximum temperature show significant increases while mean and minimum summer temperatures show consistent decline. Increase in diurnal temperature range (DTR) is consistently observed in all seasons and the annual dataset, a pattern shared by much of the Indian subcontinent but in direct contrast to both GCM projections and the narrowing of DTR seen worldwide. This divergence commenced around the middle of the twentieth century and is thought to result from changes in large-scale circulation patterns and feedback processes associated with the Indian monsoon.

The impact of observed seasonal temperature trend on runoff is explored using derived regression relationships. Decreases of ∼20% in summer runoff in the rivers Hunza and Shyok are estimated to have resulted from the observed 1°C fall in mean summer temperature since 1961, with even greater reductions in spring months. The observed downward trend in summer temperature and runoff is consistent with the observed thickening and expansion of Karakoram glaciers, in contrast to widespread decay and retreat in the eastern Himalayas. This suggests that the western Himalayas are showing a different response to global warming than other parts of the globe.

Corresponding author address: Dr. H. J. Fowler, WRSRL School of Civil Engineering and Geosciences, Cassie Building, University of Newcastle, Newcastle upon Tyne NE1 7RU, United Kingdom. Email: h.j.fowler@ncl.ac.uk

Abstract

Temperature data for seven instrumental records in the Karakoram and Hindu Kush Mountains of the Upper Indus Basin (UIB) have been analyzed for seasonal and annual trends over the period 1961–2000 and compared with neighboring mountain regions and the Indian subcontinent. Strong contrasts are found between the behavior of winter and summer temperatures and between maximum and minimum temperatures. Winter mean and maximum temperature show significant increases while mean and minimum summer temperatures show consistent decline. Increase in diurnal temperature range (DTR) is consistently observed in all seasons and the annual dataset, a pattern shared by much of the Indian subcontinent but in direct contrast to both GCM projections and the narrowing of DTR seen worldwide. This divergence commenced around the middle of the twentieth century and is thought to result from changes in large-scale circulation patterns and feedback processes associated with the Indian monsoon.

The impact of observed seasonal temperature trend on runoff is explored using derived regression relationships. Decreases of ∼20% in summer runoff in the rivers Hunza and Shyok are estimated to have resulted from the observed 1°C fall in mean summer temperature since 1961, with even greater reductions in spring months. The observed downward trend in summer temperature and runoff is consistent with the observed thickening and expansion of Karakoram glaciers, in contrast to widespread decay and retreat in the eastern Himalayas. This suggests that the western Himalayas are showing a different response to global warming than other parts of the globe.

Corresponding author address: Dr. H. J. Fowler, WRSRL School of Civil Engineering and Geosciences, Cassie Building, University of Newcastle, Newcastle upon Tyne NE1 7RU, United Kingdom. Email: h.j.fowler@ncl.ac.uk

1. Introduction

Instrumental records show a systematic increase in global mean temperature (Folland et al. 2001a, b), with global mean temperature increasing at a rate of 0.07°C decade−1 over the last century (Jones and Moberg 2003). In addition, the 1990s were the warmest decade, and 1998 the warmest year since the start of the global mean temperature record in 1856 (Jones and Moberg 2003). However, the warming has not been globally uniform. High northern latitudes have been particularly affected, with reconstructions of mean surface temperature over the past two millennia suggesting that the late twentieth century warmth is unprecedented (Mann and Jones 2003) and attributed to the anthropogenic forcing of climate (Thorne et al. 2003).

In most parts of the world there have been differential changes in daily maximum and minimum temperatures, resulting in both a narrowing of diurnal temperature range (DTR) and an increase in mean temperature (Karl et al. 1993; Easterling et al. 1997; Jones et al. 1999). However, there are suggestions that the western Himalaya region is showing a different response to global warming (Kumar et al. 1994; Yadav et al. 2004), with an increase in DTR and a cooling of mean temperature in some seasons, possibly as a result of local forcing factors.

Meteorological data from the western Himalaya and trans-Himalaya mountains have been little studied. However, here we report on instrumental records of temperature data obtained for the Upper Indus Basin (UIB) for 1900–2000. The UIB (Fig. 1) is a key area for climatological studies as it is located at the boundary between tropical and continental climatic influences and has quite different climatic controls from the eastern or Greater Himalaya. Indian monsoon air masses, which bring significant rainfall on the southern margin, penetrate infrequently across the Himalayan mountain divide to the trans-Himalayan regions of the Karakoram and Hindu Kush Mountains. Precipitation in these northern ranges is concentrated in winter and spring months (Archer and Fowler 2004) and carried on westerly disturbances originating predominantly from Mediterranean and Caspian Sea regions. Winter precipitation provides the principal source for accumulation on UIB glaciers in the greatest area of perennial ice outside the polar regions (22 000 km2) and upon melting, for runoff in the river Indus.

In the UIB, summer runoff (July to September) is highly correlated with winter precipitation in middle-elevation basins dependent on an ephemeral snowpack (Fowler and Archer 2005). However, summer temperature and runoff are negatively correlated on these snow-fed catchments where increased temperature results in increased evaporative loss and, since snow cover volume is limiting, in reduced runoff (Singh and Bengtsson 2005).

In contrast, summer runoff volume in major high-elevation Indus tributaries such as the rivers Hunza and Shyok (Fig. 1), fed by glaciers and permanent snowpack, is uncorrelated with winter precipitation but highly correlated with summer mean temperature (Archer 2003). Linear regression analysis by Archer (2003) suggests that a 1°C rise in mean summer temperature arising from climate change would result in an increase of 17% in summer runoff for the river Shyok (basin area 65 025 km2) and a 16% increase for the river Hunza (basin area 13 925 km2), respectively. Together these two catchments contribute more than 25% of the inflow to Tarbela Dam, which is the main controlling structure for the Indus Basin Irrigation System, one of the world’s largest integrated irrigation networks.

These significant, but opposing, responses to variations in summer temperatures indicate potentially important practical impacts of temperature trend on river runoff and consequently upon local rural livelihoods in the mountains and the productivity of agriculture on the Indus plains. In addition, temperature trend (along with precipitation) will influence glacier mass balance. While most of the world’s mountain glaciers have been shrinking for at least the last 30 years (World Glacier Monitoring Service 2002), including the neighboring Greater Himalaya (Hasnain 1999; Mastny 2000; Shrestha and Shrestha 2004). Hewitt (1998) reports the widespread expansion of larger glaciers in the central Karakoram, accompanied by an exceptional number of glacier surges. This contrast in the behavior of glaciers suggests a pattern of climatic change in Karakoram and Hindu Kush different from that in the Greater Himalaya.

The main focus of this study is the investigation of trend in seasonal temperature series over the period from 1900 to 2000 for seven instrumental records in the UIB. However, before this is addressed, the extent to which individual temperature records or groups of records are representative of the region in which they are located has been investigated.

2. Data

The analysis investigates temperature data for six trans-Himalayan stations. Srinagar, in the neighboring Jhelum Basin south of the Himalayan divide, has also been included for comparison. The stations cover a distance east to west of about 300 km and from north to south of over 200 km. Several data sources have been used, with measurement locations shown in Fig. 1 and listed in Table 1.

  1. For the prepartition period up to 1947, climate statistics for Gilgit and Skardu were derived from daily maximum and minimum temperature records in 6-monthly daily weather reports published by the India Meteorological Department (IMD). Records commence in 1900 and 1903 for Skardu and Gilgit, respectively, and are essentially complete from 1905 to 1935 and thereafter intermittent until 1947.

  2. For the postpartition period, the Pakistan Meteorological Department (PMD) provided daily records of maximum and minimum temperature for Gilgit and Skardu to 1999, and monthly records for Astore, Dir, Drosh, and Bunji to 1997. With the exception of Dir (starting in 1967) these records commence in the 1950s.

  3. Monthly mean temperature records for Srinagar (1894–2000) were obtained from the Climate Research Unit at the University of East Anglia, United Kingdom.

Validation was based primarily on the inspection of records as station history files were not available. Interviews with PMD staff indicated that standard meteorological practice, established by IMD in 1891, had been adopted by PMD. However, site changes had occurred during partition at Gilgit and Skardu. The principal validation check of records was through double mass curve analysis, which was carried out by season and separately for maximum and minimum temperatures. No evidence of heterogeneity was found in any of the maximum temperature series. However, two years of inconsistent winter minima were identified at Drosh and these were substituted by regression with other stations. In addition, evidence of heterogeneity was identified in the winter minimum record at Bunji with a break of slope in the double mass curve in 1977; the pre-1977 record was adjusted to be consistent with the later record. While Bunji has been included in time series analysis, results must be considered of lower reliability than other stations.

The network has recently been supplemented by more than 20 automatic observations stations at elevations up to 4700 m (Hewitt and Young 1993) and although record lengths are still insufficient for time series analysis, they provide a basis for assessing altitudinal variations in temperature (Archer 2004).

3. Spatial correlation of temperature

Mountainous regions are said to require a much greater density of climate stations than neighboring flatlands to achieve the same reliability of areal estimates. In this most rugged of world landscapes, the existing network of weather stations falls short by more than an order of magnitude of the minimum requirement recommended by World Meteorological Organization (1970). The climatological network in the UIB is also biased by the predominant location of stations on valley floors. To assess regionally representative trends, it is necessary to establish the extent to which seasonal variations are correlated across the region. Archer (2004) investigated such spatial relationships by drawing up correlation matrices of annual, seasonal, and monthly temperature between weather stations in the region. A summary of the results for the key period of ablation, divided between the 3-month periods April–June and July–September is shown in Table 2.

Conclusions drawn from this table and from the related analysis are as follows:

  1. In spite of the separation of stations in distance and by major topographic barriers, seasonal temperatures and their anomalies are remarkably conservative. Positive correlation coefficients occur between all stations, with better correlation over shorter distances. The high correlation between temperature at Srinagar, south of the main Himalayas, and the northern valley stations of Gilgit and Skardu is particularly notable.

  2. Correlation of mean temperature is much higher during spring and summer than during winter months. It is suspected that in winter, under the prevailing influence of the Tibetan anticyclone, more local conditions such as valley temperature inversions prevail.

  3. On average, the station with the best seasonal correlation with other stations is Astore. For the spring and summer season the correlation coefficient for mean temperature is greater than 0.75 with all stations except Drosh (July–September).

  4. Correlation of monthly mean temperature at high-altitude stations with short records at Kunjerab (4733 m) and Shandur (3719 m) with Gilgit and Skardu yielded on average eight months with correlation coefficients greater than 0.8.

The high correlation in mean temperature between valley stations separated by considerable horizontal distances and intervening mountain barriers strongly implies that these correlations will be reflected over much shorter vertical distances. This view is supported by the high correlation coefficients demonstrated between valley stations and short records from high-altitude stations (Archer 2004) and also by the strong observed correlation between valley temperature stations and catchment runoff that is generated at higher elevations (Archer 2003). It is concluded that the temperature time series at valley stations can be used to provide a regional picture of seasonal and year-to-year variation in temperature.

4. Analysis of temporal temperature change

a. Recent trends in the UIB, 1961–2000

Trends in annual, seasonal, and monthly mean, maximum, and minimum temperature were investigated for stations at Dir, Drosh, Bunji, Gilgit, Astore, and Skardu for the period 1961–99. For Srinagar, only mean temperature for the period 1961–2000 was available for investigation. The records were analyzed by fitting a linear least squares trend line to the annual deviation from the mean and assessing the significance of trend using Student’s t test.

A summary of the results of the analysis of annual and seasonal trend in mean, maximum, and minimum temperature is provided in Table 3 as a basis for illustrating common patterns and regional differences. Increases (+) and decreases (–) are shown as change in °C decade−1.

1) Mean temperature

In contrast to many studies indicating unprecedented warming in the twentieth century (e.g., Jones and Moberg 2003), about half the records in the UIB show a reduction in mean annual temperature since the 1960s and more strongly since the mid-1970s. The trend is statistically significant (p < 0.05) at Bunji and Dir, with cooling rates of –0.52 and –0.26°C decade−1, respectively, since 1961 (Table 3a), and is related to significant cooling of mean temperature during spring [March to May (MAM)] and summer [June to August (JJA)].

All records show summer cooling. This is large and statistically significant at some stations, with rates of –0.30, –0.99, –0.35 and –0.38°C decade−1 at Astore, Bunji, Dir, and Gilgit, respectively (Table 3a). Similarly, a cooling of spring mean temperature is seen at four out of seven stations, although this is only significant at two stations, and the Skardu record exceptionally shows significant warming (Table 3a). Change in fall [September to November (SON)] mean temperature over the same period is comparable to that of spring.

In contrast, in winter [December to February (DJF)] a warming of mean temperature is seen in most records from 1961 to 2000. There is evidence of a significant winter warming trend at Skardu and Srinagar (p < 0.05), Gilgit (p < 0.10), and Drosh, with an insignificant trend seen at other stations.

2) Maximum and minimum temperature

Since 1961 there has been an increase in annual maximum temperature at all stations except Bunji and Astore (Table 3b). Significant increases in maxima occur mainly in the months from November to April. However, there has also been a contrasting and large regionwide reduction in annual minimum temperature (Table 3c). This has resulted in the most striking change in the UIB since 1961: the large increase in DTR at all stations. Figure 2a shows the large increase in annual DTR at Skardu, Gilgit, and Dir during 1961–2000. This increase is strongest at the northern stations, with Skardu, Gilgit, Bunji, and Dir showing positive trends in DTR in all months as well as all seasons. The largest change in DTR is found during the period from June to November, summer and fall seasons, when all stations in the UIB show a positive trend (see Fig. 2b).

Figure 3 clearly shows the regionwide increase in winter maximum temperature over the period 1961–99 at stations in the UIB. There have been smaller changes in winter minimum temperature, with small increases in minima at some sites and decreases in minima at others (Fig. 4). Since temperature maxima have been increasing at a greater rate than temperature minima, even where temperature minima are also increasing, there has been a resulting increase in winter DTR at all stations, and an overall winter warming.

In contrast, in summer months there has been a regionwide reduction in temperature maxima (Fig. 5; Table 3c). Perhaps more surprising, however, is the concomitant reduction in summer temperature minima (Fig. 6). This is statistically significant (p < 0.05) at all stations except Drosh, and there are high rates of cooling, with reductions of −1.11°C decade−1 at Bunji (Table 3c) and over −0.4°C decade−1, at all stations. These reductions are at their highest in June. Summer cooling of mean temperature therefore arises as a result of reduction in both minimum and maximum temperature, although the former has by far the larger effect. Since summer minimum temperatures show greater rates of cooling than maximum temperatures, there has also been a significant increase in summer DTR.

In spring and fall, a similar cooling trend is also observed in temperature minima (Table 3c). Combination with a significant upward trend in temperature maxima at many stations in fall (Table 3b) results in an increasing DTR in spring and fall.

b. Long-term change at Gilgit, Skardu, and Srinagar

The duration of temperature records at Gilgit and Skardu is unique for currently operating trans-Himalayan stations. Despite discontinuity in the temperature record arising from site changes, comparison has been made of trend in maximum, minimum, and mean monthly temperatures and diurnal temperature range for the periods 1905–35 and 1958–2000 and mean temperature only for Srinagar.

There are significant features in common between the UIB stations and Srinagar during the later period. All records show a warming trend in winter months related to increases in maximum temperature and cooling in summer months related to decreases in minimum temperature.

Both the significant correlation in seasonal temperature between Srinagar and UIB stations (Table 2) and the similarity of the seasonal pattern of change with Gilgit and Skardu suggest that the homogeneous record at Srinagar may provide an analog of long-term climatic change for the UIB.

Over the longer period from 1894 to 2000, there has been a small but significant increase in annual mean temperature at Srinagar (p < 0.01) of 0.07°C decade−1. During the winter season (DJF), prior to 1960 the temperature increase was significant (p < 0.05) at 0.10°C decade−1 (Fig. 7a). However, since 1960 the rate of increase has accelerated to 0.51°C decade−1. If individual months are further examined then large and statistically significant increases in mean temperature are seen in the months from December to April (p < 0.05) (Figs. 7b–f), with smaller increases during other months. This suggests that the rate of increase of winter temperatures in the UIB region may also have accelerated in the period since 1960, a possible signal of global warming.

More remarkable, however, is the sudden divergence in temperature trend between the high latitudes and the western Himalayan region in the latter half of the twentieth century. Climate proxies and instrumental records in the northern high latitudes indicate unprecedented warming (Jones and Moberg 2003), but instrumental records from some UIB stations suggest a cooling of mean temperature, annually and in all seasons except winter. Figure 8 shows the striking increase in DTR at Skardu and Gilgit in all seasons since the middle of the twentieth century. This has been caused by a pronounced cooling of minimum temperatures since 1960, particularly in summer and autumn (Fig. 9), coupled to a slight increase in maximum temperatures in all seasons since 1960 (Fig. 10). This is in contrast to the narrowing of DTR seen in the rest of the globe since 1950 (Karl et al. 1993; Easterling et al. 1997). However, during 1900–35 DTR in the western Himalayas remains constant or shows a decreasing trend. This is akin to trends seen in the remainder of the northern high latitudes since 1950 (e.g., Karl et al. 1993).

5. Discussion

While there are clearly differences between stations in the amount and even the direction of change, several common patterns emerge, especially when the series are decomposed into seasonal groupings or minimum and maximum series. Trends are less evident in annual and mean temperature series because of differences in behavior between seasons or between maxima and minima. The observed trends in the UIB may be compared with those in neighboring parts of the Himalayas, central Asia, and the Indian subcontinent. The most significant common patterns in the UIB with respect to the 1961 to 2000 period are the following:

  1. Decrease in summer mean temperature (seven of seven sites, four significant)

  2. Decrease in summer minimum temperature (six of six sites, five significant)

  3. Decrease in fall minimum temperature (six of six sites, four significant)

  4. Increase in fall maximum temperature (four of six sites, three significant)

  5. Increase in winter mean temperature (five of seven sites, three significant)

  6. Increase in winter maximum temperature (five of six sites, three significant)

  7. Decrease in annual minimum temperature (six of six sites, four significant)

  8. Increase in annual maximum temperature (four of six sites, two significant)

a. Summer temperatures

The consistent decline in summer mean and minimum temperatures is perhaps surprising in view of the frequently observed and reported retreat of neighboring Himalayan glaciers, which has generally been ascribed to an increase in summer temperatures (Yamada et al. 1992; Hasnain 2002). Thus, Shrestha et al. (1999) derived a consistent upward trend in summer temperatures (June to September) for the period 1977–94, ranging from 0.14°C for low-elevation stations to significant increases of 0.62°C for the Nepal Himalayas and 1.09°C for trans-Himalayan stations. However their analysis was based on maximum temperatures, while minimum temperatures, graphed but not analyzed, show steady or falling trends. The decline in minimum temperatures was confirmed by Sharma et al. (2000) for several stations in Nepal. Only two of six stations in the UIB show an upward trend in maximum summer temperatures, but the rates of decrease in the remaining four stations is much less than for minimum temperatures.

Cook et al. (2003) reexamined a longer Kathmandu mean temperature record and compared it with a 0.5° latitude–longitude gridded dataset based on records from neighboring northern India; both showed a cooling trend in the monsoon (June to September) as well as in premonsoon (February to May) for the period 1901–95. In addition, Hingane et al. (1985) found that northwest India shows a significant decreasing trend in mean temperature in summer (−0.05°C decade−1), confirmed by Kumar et al. (1994) for minimum temperature at the same decreasing rate while maximum temperatures have been steady. Summer cooling trends during the last few decades of the twentieth century have also been documented for the Tibetan Plateau (Liu and Chen 2000), central China (Hu et al. 2003), and central Asia (Briffa et al. 2001).

Cook et al. (2003) attributed this summer cooling to a coincident increase in monsoonal rainfall in Nepal after 1950, which would tend to suppress temperatures through increased daytime cloud cover. In the UIB, increasing summer rainfall since 1961 (Archer and Fowler 2004) and a coincident increase in daytime cloud cover may also provide an explanation for the summer cooling of maximum temperature. It does not, however, provide an explanation for the large decreases in minimum temperature seen in the summer. The decrease in summer temperature minima could be explained by a coincident decrease in cloud cover during the night, although currently no data are available to assess this.

Other research has argued for an aerosol induced cooling over much of South Asia and the Indian subcontinent. Ramanathan et al. (2005) suggest that during the “Indian dry season,” from October to May, the haze from the Asian Brown Cloud reduces the surface receipt of solar radiation and precipitation efficiency, causing the observed cooling of sea surface temperatures over the North Indian Ocean. This, in turn, has reduced monsoonal summer precipitation by ∼5% in India. However, it is unlikely that aerosol induced cooling offers an explanation for the observed summer cooling in the UIB as during the monsoon the effect of the Asian Brown Cloud is diminished in India. Additionally, the UIB region has seen a recent increase in summer precipitation (Archer and Fowler 2004).

Finally, there is the possibility that upward trends in summer temperature are more marked at higher elevations as shown by the analysis of Shrestha et al. (1999) for Nepal; this would imply an increased lapse rate. A tendency for warming trends to increase with elevation has been found in both the Tien Shan of central Asia (Aizen et al. 1997) and the Tibetan Plateau (Liu and Chen 2000). There is no evidence of elevation dependency of warming (cooling) in our dataset although this would not necessarily be expected as the range in station elevation is only 1000 m.

b. Winter temperatures

The UIB shows a consistent increase in winter mean temperatures, with increasing trends since 1961 of 0.17° and 0.38°C decade−1 at Gilgit and Skardu, respectively. The increase in winter maximum temperature is large and statistically significant (p < 0.05) at Gilgit, Skardu, and Dir, with increases of 0.27°, 0.55°, and 0.51°C decade−1, respectively. However, winter minima only show positive trend at three of six stations and none is significant. The greatest warming in neighboring regions has also been in winter temperature. For India, Kumar et al. (1994) derived winter warming trends of 0.09°C decade−1 for maximum temperature in an all-India temperature series and 0.08°C decade−1 for the northwest region for 1901– 87. Shrestha et al. (1999) found increases in winter maxima of 0.61°, 0.90°, and 1.24°C decade−1 for all-Nepal, Himalayan, and trans-Himalayan stations, respectively. Cook et al. (2003) found indications of warming since 1960 at Kathmandu only in the period from October to January. Additionally, decadal winter increases of 0.32°C in mean temperature have been suggested for the Tibetan Plateau by Liu and Chen (2000).

Change in winter temperature is less hydrologically significant than change during the summer as the most hydrologically active areas above an elevation of 4000 m remain above the freezing level even for maximum temperatures from November to March (Archer 2004). However, average winter slope lapse rates of 0.70°C (100 m)−1 (Archer 2004) suggest a shift in the freezing level of up to 300 m and a decrease in the ratio of solid to total precipitation at lower elevations. The perception of decreased winter precipitation and increased duration of opening of passes in the Karakoram and Hindu Kush through the winter may therefore be explained by an increase in winter temperature rather than by the decrease in winter precipitation suggested by Miehe (1992).

c. Asymmetry of temperature trends

The most striking temperature change in the UIB is the recent large increase in DTR observed in all seasons and in the annual dataset. Contrasting changes in maximum and minimum temperature for Gilgit and Skardu shown in Figs. 9 and 10 suggest that this divergence commenced around the middle of the twentieth century. This change is not unique to the high-elevation UIB but has also been observed over much of the Indian subcontinent (Kumar et al. 1994; Yadav et al. 2004). In contrast, in many other parts of the world and especially at higher latitudes, minimum temperatures have increased at a higher rate than maximum temperatures, resulting in a narrowing of DTR since 1950 (Karl et al. 1993; Easterling et al. 1997; Jones et al. 1999). Several explanations have been put forward for the increase in DTR on the Indian subcontinent.

First, it has been proposed by Cook et al. (2003) that change in DTR may result from change in precipitation regime. Cloud cover associated with precipitation reduces incoming solar radiation and traps more outgoing longwave radiation, hence altering surface temperature, and up to 80% of DTR variance can be explained by cloud and precipitation records (Dai et al. 1999). Given the observed increase in DTR, one would expect a decrease in precipitation, with warmer days and cooler nights. However, Archer and Fowler (2004) found an increase in annual precipitation in the UIB from 1961 to 2000 with a statistically significant (p < 0.05) upward trend in precipitation during both winter and summer seasons.

Simple monthly correlation of various indices of temperature with precipitation is shown in Table 4. There is a significant (p < 0.05) negative relationship between maximum temperature and precipitation in summer and autumn months, suggesting that increased precipitation may induce a cooling effect on maximum temperature. A negative relationship is found between DTR and precipitation in all months, suggesting that an increase in precipitation should cause a decrease in DTR. Therefore, although change in summer precipitation regime provides a possible explanation for reduction in summer maximum temperature in the UIB, it does not explain the increase in DTR.

Second, change in snow-covered area and thus regional albedo may cause change in DTR. Snow-covered area is mainly dependent on winter precipitation. Low winter precipitation can therefore lower the regional albedo and cause more solar radiation to be absorbed at the surface, causing increased surface melting. This positive feedback may result in increased daytime maxima (Giorgi et al. 1997). D. R. Archer (2001, unpublished manuscript) illustrated this feedback for the UIB by examining correlation between winter half-year precipitation and spring temperature. Significant negative correlation was found only for the months of March and April and did not persist into the summer. Additionally, as noted above, winter precipitation has shown an increasing trend in the UIB, significant (p < 0.05) at three stations (Archer and Fowler 2004). Snow albedo feedback is therefore unlikely to have caused the observed change in DTR. Besides, Kumar et al. (1994) report that the temperature asymmetry is shared with northwest India where there is no question of snow albedo effects.

Third, in the lowlands of the Indian subcontinent, land-use changes including urbanization (Kumar and Hingane 1988), deforestation (Kothyari and Singh 1996; Yadav et al. 2004), intensification of agriculture (Singh and Sontakke 2002), and widespread irrigation (Adel 2002) have been associated with change in temperature. Urbanized areas often show a narrower DTR than nearby rural areas (Gallo et al. 1996) as urbanization increases nighttime temperature minima through the urban heat island effect. Land-use change such as deforestation may reduce temperature minima. Yadav et al. (2004) suggest that land-use change, particularly deforestation, has contributed to a widening of the DTR in northern India. The removal of tree canopies causes direct heating of the soil during the daytime to increase. This is then followed by the intense cooling of these soils by terrestrial radiation loss at night, enhanced by the fact that they are moisture deficient due to surface runoff increases associated with deforestation.

Development activities leading to widespread deforestation do appear to offer a plausible explanation for the increase in DTR in northern India with, for example, 50% of the areal expansion of agricultural land in the Garhwal region during 1963–93 coming from previously forested area (Sen et al. 2002). However, the same trend of widening DTR is also seen in the neighboring region of the UIB where none of the same developmental pressures are relevant. Here, natural vegetation is sparse and limited either by low rainfall at lower elevations or by low temperature at higher elevations. The UIB is also sparsely populated and neither the limited urbanization nor land-use change on valley floors is considered sufficient to have initiated significant changes in thermal regime. The upward trend in DTR in the UIB is thus difficult to explain but the regionwide increase in DTR in the western Himalayas suggests that changes in atmospheric circulation patterns are likely to be responsible.

GCM future climate simulations for the region from the Intergovernmental Panel on Climate Change (IPCC) suggest more pronounced increases in minimum rather than maximum temperatures over Asia related to increased monsoon cloudiness, hence an increase in mean temperature annually and during all seasons and a decrease in annual and winter DTR (Lal et al. 2001; Watson et al. 2001). During the summer an increase in DTR is simulated, but this results from a more pronounced increase in maximum temperature relative to the minimum temperature (Lal and Harasawa 2001) rather than the reduction in both summer maximum and minimum temperatures seen since 1961 in the UIB. GCM integrations also suggest that variations in DTR are much more sensitive to changes in feedbacks than in direct forcings, with variations arising mostly from changes in clouds and in soil moisture (Stone and Weaver 2003). It is possible that global warming has induced process coupling in this region that occurs on scales too small for current GCMs to resolve well. Using a regional climate model, Pan et al. (2004) were able to simulate a local minimum of warming in the central United States (a “warming hole”) associated with changes in low-level circulations caused by regional-scale feedback processes under global warming. This was absent from GCM simulations of the same region, but this region is known to have experienced a cooling trend in the last 25 yr of the twentieth century. A relationship between circulation changes during the Northern Hemisphere winter and winter DTR has also been examined by Easterling et al. (1997) for the European region. They found a bipolar pattern resulting from stronger westerly flows, with Iberia experiencing increased DTR and decreased DTR over northern Europe.

While explanation remains problematic, the shared response with the Indian subcontinent suggests a regional or latitudinal effect associated with large-scale climatic processes and their feedbacks, and related to changes in the monsoon circulation. Recent modeling work by Ramanathan et al. (2005) suggests that aerosol cooling brought about by the Asian Brown Cloud in the nonmonsoonal months from October to May has caused the observed cooling of sea surface temperatures over the North Indian Ocean. This, in turn, weakens the monsoon circulation and causes a southward shift, reducing rainfall in India by ∼5%. In the UIB, summer precipitation has increased but this may be due to the incursion of more westerly weather rather than monsoonal incursions. The relationship with circulation patterns will be investigated in further work.

d. Impacts on hydrological regime and glaciological response

The overall impact of climate change on glacier mass balance and hydrological regime require the joint consideration of change in both precipitation and the energy budget. The emphasis here is on the energy budget as represented by temperature. Trend in temperature will have both short- and long-term effects on glaciers and runoff.

Short-term effects of seasonal temperature are revealed by the analysis of current variability. A seasonal mean temperature range for both spring (April to June) and summer (July to September) of 4° to 5°C since 1961 provides a basis for the assessment of impact on runoff. Observed summer runoff on high-altitude glacier-fed catchments shows significant positive correlation with summer temperatures. For the river Hunza a 1°C fall below the current mean temperature—the approximate fall in summer mean temperature since 1961—gives a predicted 16% decrease in summer runoff (r = 0.69); the equivalent decrease for the river Shyok is 17% (r = 0.73) using derived regression equations (Fig. 11). An equivalent figure of 16%–18% increase for a 1°C temperature rise was identified from modeling studies on western Himalayan catchments by Singh and Kumar (1997). Predicted runoff decreases for a 1°C fall in the spring months are more striking, 48% and 62% for the rivers Hunza (r = 0.73) and Shyok (r = 0.53), respectively, albeit with much smaller runoff totals.

On average, climate stations in the UIB show a reduction in mean summer temperature of 1.2°C over the period 1961 to 2000 (Fig. 11a), representing a fall in average summer runoff of ∼20% on the Hunza and Shyok (Fig. 11b). Similarly, during the more sensitive spring months a 40-yr mean temperature reduction of 0.3°C yields a fall in spring runoff of 15% to 20%. These results are confirmed for the river Hunza at Dainyor where there is a significant downward trend in observed runoff for spring and summer of 46% and 35%, respectively, over the record period (Fig. 11c). The absence of significant trend for the river Shyok (Fig. 11c) may be explained by regional variation in climatic trends. The nearest climate station to the Shyok catchment at Skardu has an upward trend in spring and a less steep downward trend in summer than other UIB stations (Table 3a).

In contrast, temperature and runoff are negatively correlated on middle-altitude (snow-fed) catchments although the regression relationships are nonsignificant and confounded by the much stronger influence of winter precipitation (Archer 2003). Singh and Bengtsson (2005) suggest through modeling studies that increased temperature results in increased evaporative loss and—since snow cover volume is limiting—in reduced runoff, with an estimated decrease of 9% for a 1°C rise in temperature for the Satluj basin in the western Himalayas.

These varied responses to changes in the driving variables mean that caution must be exercised in the prediction of runoff response to climate change, especially in complex catchments such as the main stem of the river Indus where a mix of climatic controls impact upon hydrological regimes.

In the longer term, a reduction in summer temperature implies reduced ablation and a positive mass balance of glaciers, given constant precipitation. In addition, climate stations in the Karakoram show consistent positive trends in winter precipitation, averaging a 7% increase per decade for the period since 1961 (Archer and Fowler 2004). In contrast, the observed higher temperatures in winter and associated decrease in the ratio of solid to total precipitation is expected to have an effect only on the lower reaches of glaciers. While the lag time of large Karakoram glaciers may be many decades, the reported expansion and thickening of many glaciers (Hewitt 1998, 2005) is consistent with the observed climatic changes. Karakoram glaciers offer a significant contrast with those in the eastern Himalayas (Hasnain 1999).

6. Conclusions

Although there are variations between the stations analyzed, in terms of the amount and direction of change in temperature, common patterns emerge that can be compared with neighboring mountain regions and the Indian subcontinent:

  1. Mean and minimum summer temperatures show a consistent cooling in the period since 1961, a characteristic shared with northwest India and, at least, for lower-level stations in Nepal. Summer maxima show little clear trend.

  2. Winter mean and maximum temperatures show statistically significant increases. Although winter minima show no significant trend, the large positive changes in the maxima are sufficient to tip the balance of change in the mean annual temperature to significantly positive at the most easterly stations of Skardu and Srinagar.

  3. The most striking recent change is the large increase in diurnal temperature range, observed at all seasons and in the annual dataset. Long records at Gilgit and Skardu suggest that the trend of divergence commenced in the middle of the twentieth century. This change is in common with much of the Indian subcontinent but in direct contrast to most other parts of the globe where a narrowing of DTR has been observed. Explanations are problematic but appear to be associated with the coupling of regional feedbacks such as aerosol cooling to large-scale climatic processes.

  4. Using linear relationships between spring and summer temperatures and runoff in the rivers Hunza and Shyok and the cooling regional temperature trend produced a predicted reduction of 20% in runoff over the record period. This predicted fall was exceeded by actual runoff decreases on the river Hunza. The absence of equivalent decline in the river Shyok can be explained by regional variations in temperature trend, with easterly stations in the UIB showing more positive temperature trend.

  5. Summer temperature reductions and positive trend in winter precipitation imply reduced ablation and increased accumulation of Karakoram glaciers. These climatic changes are consistent with the observed thickening and expansion of glaciers in the UIB region, in contrast to widespread retreat and decay in the eastern Himalayas.

  6. The hydrological and glaciological impact of change in temperature is highly dependent upon the season in which the change occurs. Summer cooling has a widespread impact at all elevations whereas the effects of winter warming are restricted to lower elevations only.

  7. The viability of temperature measurements at valley stations as indices of change at higher elevations is confirmed both by the strong links between temperature and runoff and by the strong spatial correlation in seasonal temperature across the region.

Acknowledgments

The work was initiated while the second author was a volunteer with Voluntary Service Overseas (VSO) and employed by the German Agency for Technical Development (GTZ). The authors wish to thank colleagues of both agencies for the opportunity and their support. Particular thanks are due to GTZ counterpart Dr. Juan Jose Victoria. Support was also provided by officials at WAPDA and at the Pakistan Meteorological Department, and the help of the director, Shaukat Ali Awan, is acknowledged for supply of data and helpful discussion. The authors also thank Chris Kilsby (Newcastle University), Phil Jones (Climatic Research Unit, University of East Anglia), and three anonymous reviewers for comments that helped to improve the paper.

REFERENCES

  • Adel, M. M., 2002: Man-made climatic changes in the Ganges basin. Int. J. Climatol., 22 , 9931016.

  • Aizen, V. B., E. M. Aizen, J. M. Melack, and J. Dozier, 1997: Climatic and hydrologic changes in the Tien Shan, central Asia. J. Climate, 10 , 13931404.

    • Search Google Scholar
    • Export Citation
  • Archer, D. R., 2003: Contrasting hydrological regimes in the Indus Basin. J. Hydrol., 274 , 198210.

  • Archer, D. R., 2004: Hydrological implications of spatial and altitudinal variation in temperature in the Upper Indus Basin. Nord. Hydrol, 35 , 3. 209222.

    • Search Google Scholar
    • Export Citation
  • Archer, D. R., and H. J. Fowler, 2004: Spatial and temporal variations in precipitation in the Upper Indus Basin, global teleconnections and hydrological implications. Hydrol. Earth Syst. Sci., 8 , 4761.

    • Search Google Scholar
    • Export Citation
  • Briffa, K. R., T. J. Osborn, F. H. Schweingruber, I. C. Harris, P. D. Jones, S. G. Shiyatov, and E. A. Vaganov, 2001: Low frequency temperature variations from northern tree ring density network. J. Geophys. Res., 106 , 29292941.

    • Search Google Scholar
    • Export Citation
  • Cook, E. R., P. J. Krusic, and P. D. Jones, 2003: Dendroclimatic signals in long tree-ring chronologies from the Himalayas of Nepal. Int. J. Climatol., 23 , 707732.

    • Search Google Scholar
    • Export Citation
  • Dai, A., K. E. Trenberth, and T. R. Karl, 1999: Effects of clouds, soil moisture, precipitation, and water vapor on diurnal temperature range. J. Climate, 12 , 24512473.

    • Search Google Scholar
    • Export Citation
  • Easterling, D. R., and Coauthors, 1997: Maximum and minimum temperature trends for the globe. Science, 277 , 364367.

  • Folland, C. K., and Coauthors, 2001a: Observed climate variability and change. Climate Change 2001: The Scientific Basis, J. T. Houghton et al., Eds., Cambridge University Press, 99–182.

    • Search Google Scholar
    • Export Citation
  • Folland, C. K., and Coauthors, 2001b: Global temperature change and its uncertainties since 1861. Geophys. Res. Lett., 28 , 26212624.

  • Fowler, H. J., and D. R. Archer, 2005: Hydro-climatological variability in the Upper Indus Basin and implications for water resources. Regional Hydrological Impacts of Climatic Change—Impact Assessment and Decision Making, T. Wagener et al., Eds., IAHS Publication 295, 131–138.

    • Search Google Scholar
    • Export Citation
  • Gallo, K. P., D. R. Easterling, and T. C. Peterson, 1996: The influence of land use/land cover on climatological values of the diurnal temperature range. J. Climate, 9 , 29412944.

    • Search Google Scholar
    • Export Citation
  • Giorgi, F., J. W. Hurrell, M. R. Marinucci, and M. Beniston, 1997: Elevation dependency of the surface climate change signal: A model study. J. Climate, 10 , 288296.

    • Search Google Scholar
    • Export Citation
  • Hasnain, S. I., 1999: Report of Working Group on Himalayan Glaciology. International Commission on Snow and Ice, 22 pp.

  • Hasnain, S. I., 2002: Himalayan glaciers meltdown: Impact on South Asian rivers. FRIEND 2002—Regional Hydrology: Bridging the Gap between Research and Practice, H. A. J. van Lanen and S. Demuth, Eds., IAHS Publication 274, 417–423.

    • Search Google Scholar
    • Export Citation
  • Hewitt, K., 1998: Glaciers receive a surge of attention in the Karakoram Himalaya. Eos, Trans. Amer. Geophys. Union, 79 , 104105.

  • Hewitt, K., 2005: The Karakoram anomaly? Glacier expansion and the “elevation effect,” Karakoram, Himalaya. Mt. Res. Dev., 25 , 332340.

    • Search Google Scholar
    • Export Citation
  • Hewitt, K., and G. J. Young, 1993: The Snow and Ice Hydrology Project: A Pakistan–Canada research and training programme. Snow and Glacier Hydrology: Proceedings of International Symposium, Kathmandu, Nepal, 16–21 November 1992, G. J. Young, Ed., IAHS/AISH Publication 218, 49–58.

    • Search Google Scholar
    • Export Citation
  • Hingane, L. S., K. R. Kumar, and B. V. R. Murty, 1985: Long term trends of surface air temperature in India. J. Climatol., 5 , 521528.

    • Search Google Scholar
    • Export Citation
  • Hu, Z. Z., S. Yang, and R. G. Wu, 2003: Long-term climate variations in China and global warming signals. J. Geophys. Res., 108 .4614, doi:10.1029/2003JD003651.

    • Search Google Scholar
    • Export Citation
  • Jones, P. D., and A. Moberg, 2003: Hemispheric and large-scale surface air temperature variations: An extensive revision and an update to 2001. J. Climate, 16 , 206223.

    • Search Google Scholar
    • Export Citation
  • Jones, P. D., M. New, D. E. Parker, S. Martin, and I. G. Rigor, 1999: Surface air temperature and its changes over the past 150 years. Rev. Geophys., 37 , 173199.

    • Search Google Scholar
    • Export Citation
  • Karl, T. R., and Coauthors, 1993: A new perspective on recent global warming— Asymmetric trends of daily maximum and minimum temperature. Bull. Amer. Meteor. Soc., 74 , 10071023.

    • Search Google Scholar
    • Export Citation
  • Kothyari, U. C., and V. P. Singh, 1996: Rainfall and temperature trends in India. Hydrol. Processes, 10 , 357372.

  • Kumar, K. R., and L. S. Hingane, 1988: Long term variation of surface air temperature at major industrial cities of India. Climatic Change, 13 , 287307.

    • Search Google Scholar
    • Export Citation
  • Kumar, K. R., K. K. Kumar, and G. B. Pant, 1994: Diurnal asymmetry of surface temperature trends over India. Geophys. Res. Lett., 21 , 677680.

    • Search Google Scholar
    • Export Citation
  • Lal, M., and H. Harasawa, 2001: Future climate change scenarios for Asia as inferred from selected coupled atmosphere–ocean global climate models. J. Meteor. Soc. Japan, 79 , 219227.

    • Search Google Scholar
    • Export Citation
  • Lal, M., H. Harasawa, and D. Murdiyarso, 2001: Asia. Climate Change 2001: Impacts, Adaptation and Vulnerability, J. J. McCarthy et al., Eds., Cambridge University Press, 533–590.

    • Search Google Scholar
    • Export Citation
  • Liu, X., and B. Chen, 2000: Climatic warming in the Tibetan Plateau during recent decades. Int. J. Climatol., 20 , 17291742.

  • Mann, M. E., and P. D. Jones, 2003: Global surface temperatures over the past two millennia. Geophys. Res. Lett., 30 .1820, doi:10.1029/2003GL017814.

    • Search Google Scholar
    • Export Citation
  • Mastny, L., cited. 2000: Melting of earth’s ice cover reaches new high. World Watch Institute, Washington, DC. [Available online at http://www.worldwatch.org/alerts/000306.html .].

  • Miehe, G., 1992: Recent climatic changes in the Karakoram as indicated by vegetation features and geomorphological processes. Abstracts, Int. Symp. on the Karakoram and Kunlun Mountains, Kashi, China, UNESCO, 110.

  • Pan, Z., R. W. Arritt, E. S. Takle, W. J. Gutowski, C. J. Anderson, and M. Segal, 2004: Altered hydrologic feedback in a warming climate introduces a “warming hole.”. Geophys. Res. Lett., 31 .L17109, doi:10.1029/2004GL020528.

    • Search Google Scholar
    • Export Citation
  • Ramanathan, V., and Coauthors, 2005: Atmospheric brown clouds: Impacts on South Asian climate and hydrological cycle. Proc. Natl. Acad. Sci. USA, 102 , 53265333.

    • Search Google Scholar
    • Export Citation
  • Sen, K. K., R. L. Senwal, U. Rana, S. Nautiyal, R. K. Maikuri, K. S. Rao, and K. G. Saxena, 2002: Patterns and implications of land use/cover change: A case study in Pranmati watershed (Garhawal Himalaya, India). Mt. Res. Dev., 22 , 5662.

    • Search Google Scholar
    • Export Citation
  • Sharma, K. P., B. Moore, and C. J. Vorosmarty, 2000: Anthropogenic, climatic and hydrologic trends in the Kosi basin, Himalaya. Climatic Change, 47 , 141165.

    • Search Google Scholar
    • Export Citation
  • Shrestha, A. B., C. Wake, P. A. Mayewski, and J. Dibb, 1999: Maximum temperature trends in the Himalaya and its vicinity: An analysis based on temperature records from Nepal for the period 1971–94. J. Climate, 12 , 27752786.

    • Search Google Scholar
    • Export Citation
  • Shrestha, M. L., and A. B. Shrestha, 2004: Recent trends and potential climate change impacts on glacier retreat/glacier lakes in Nepal and potential adaptation measures. OECD Global Forum on Sustainable Development: Development and Climate Change, Paris, France, OECD, 5–14.

    • Search Google Scholar
    • Export Citation
  • Singh, N., and N. A. Sontakke, 2002: On climatic fluctuations and environmental changes of the Indo-Gangetic plains, India. Climatic Change, 52 , 287313.

    • Search Google Scholar
    • Export Citation
  • Singh, P., and N. Kumar, 1997: Impact assessment of climate change on the hydrological response of a snow and glacier melt runoff dominated Himalayan river. J. Hydrol., 193 , 316330.

    • Search Google Scholar
    • Export Citation
  • Singh, P., and L. Bengtsson, 2005: Impact of warmer climate on melt and evaporation for the rainfed, snowfed and glacierfed basins in the Himalayan region. J. Hydrol., 300 , 140154.

    • Search Google Scholar
    • Export Citation
  • Stone, D. A., and A. J. Weaver, 2003: Factors contributing to diurnal temperature range trends in twentieth and twenty-first century simulations of the CCCma coupled model. Climate Dyn., 20 , 435445.

    • Search Google Scholar
    • Export Citation
  • Thorne, P. W., and Coauthors, 2003: Probable causes of late twentieth century tropospheric temperature trends. Climate Dyn., 21 , 573591.

    • Search Google Scholar
    • Export Citation
  • Watson, R. T., and the Core Writing Team, Eds., 2001: Climate Change 2001: Synthesis Report. Cambridge University Press, 398 pp.

  • World Glacier Monitoring Service, cited. 2002: Glacier mass balance data 2000/01. [Available online at www.geo.unizh.ch/wgms/ .].

  • World Meteorological Organization, 1970: Guide to Hydrometeorological Practices. 2d ed. WMO No. 168, Tech. Paper 82.

  • Yadav, R. R., W-K. Park, J. Singh, and B. Dubey, 2004: Do the western Himalayas defy global warming? Geophys. Res. Lett., 31 .L17201, doi:10.1029/2004GL020201.

    • Search Google Scholar
    • Export Citation
  • Yamada, T., T. Shiraiwa, H. Iida, T. Kadota, T. Watanabe, B. Rana, Y. Ageta, and H. Fushimi, 1992: Fluctuations of the glaciers from the 1970s to 1989 in the Khumbu, Shorong and Langtang regions, Nepal Himalayas. Bull. Glacier Res., 10 , 1119.

    • Search Google Scholar
    • Export Citation

Fig. 1.
Fig. 1.

The Upper Indus Basin showing relief, rivers, and stations used in the analysis. Inset shows continental location of study area.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 2.
Fig. 2.

Trend in diurnal temperature range at Gilgit, Skardu, and Astore: (a) annual and (b) June to November.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 3.
Fig. 3.

Trend in winter (DJF) maximum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 4.
Fig. 4.

Trend in winter (DJF) minimum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 5.
Fig. 5.

Trend in summer (JJA) maximum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 6.
Fig. 6.

Trend in summer (JJA) minimum temperature at UIB stations from 1961 to 1999. Panels show temperature series with fitted linear trend line (dashed) given for each series.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 7.
Fig. 7.

Trend in winter (DJF) mean temperature at Srinagar from 1894 to 2000 with fitted linear trend shown for periods 1894–1960 and 1961–2000 and overall trend shown for the whole period (dashed line). Separate mean temperature time series are given for each of the months from December to April, with fitted linear trend lines.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 8.
Fig. 8.

Seasonal variations in DTR at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 9.
Fig. 9.

Seasonal variations in minimum temperature at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring. Note that different scales are used for the y axis to highlight the differences between minimum temperatures in different seasons.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 10.
Fig. 10.

Seasonal variations in maximum temperature at Skardu and Gilgit. The dashed lines represent a fitted linear trend for the periods 1905–35 and 1955–99: (a) summer and autumn and (b) winter and spring. Note that different scales are used for the y axis to highlight the differences between maximum temperatures in different seasons.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Fig. 11.
Fig. 11.

Relationship between temperature and runoff at (left) Gilgit and the river Hunza at Dainyor Bridge and (right) Skardu and the river Shyok at Yugo. (a) Temperature trend from June to August; (b) relationship between seasonal temperature and runoff (mm) from June to August; and (c) runoff trend from June to August.

Citation: Journal of Climate 19, 17; 10.1175/JCLI3860.1

Table 1.

Climate stations used in analysis.

Table 1.
Table 2.

Correlation coefficient (r) of spring mean temperature (April to June)—numbers above table (top left to bottom right) diagonal—and summer mean temperature (July to September)—numbers below table diagonal. Bold figures: significance 0.01. Italic: significance 0.05. Summarized from Archer (2004).

Table 2.
Table 3.

Trend in annual and seasonal temperature from 1961 to 1999 for (a) mean, (b) mean maximum, and (c) mean minimum, showing change in °C decade−1. Bold p < 0.05. Bold italic p < 0.10.

Table 3.
Table 4.

Correlation between monthly precipitation and various indices of monthly temperature at Astore, 1961–97: (A) maximum temperature, (B) minimum temperature, (C) mean temperature, and (D) diurnal temperature range. Bold p < 0.05.

Table 4.
Save