1. Introduction
Glaciers are sensitive objects that can indicate global climatic changes. Correct evaluation of glacier area and volume change has wide practical application in water-resource, water-supply and hydropower assessments. Despite the low amount of precipitation and extremely dry climate in central Asia, the Tien Shan hold one of the greatest concentrations of glacial ice in the mid-latitudes, and they constitute a vital source of water for more than 100 million people living in this region. Over the last 150 years, since the end of the ‘Little Ice Age’, and notably since the 1970s, the glaciers of central Asia have tended to retreat rapidly due to an increase in air temperatures and changes in precipitation partitioning (more rain than snow at high elevations) that have caused the glaciers to have a mainly negative mass balance, which has also changed river runoff regimes (Reference GlazirinGlazirin, 1996; Reference Aizen, Aizen, Melack and DozierAizen and others, 1997, 2006; Reference GlantzGlantz, 1999; Reference Denisov, Agaltseva and PakDenisov and others, 2000; Reference Agrawala, Barlow, Cullen and LyonAgrawala and others, 2001; Reference Aizen and CarpenterAizen, 2003). An increase in air temperature and melting of glaciers has a significant impact on desertification of the central Asian lowlands, on glacier outburst floods and on mudflows in the upper valleys (Reference Aizen, Aizen, Melack and DozierAizen and others, 1997, 2006; Reference SarsembekovSarsembekov, 2000). Furthermore, the stability of the central Asian hydrological systems relies on the stream-flow buffering capacity of the glacial ice.
Changes in the glacial coverage area in high mountains were poorly documented until aerial photography and particularly satellite imagery provided new techniques for remotely sensed glacier monitoring. Since the mid-1980s, the role of satellite glacier monitoring has been widely recognized. World Glacier Monitoring Service (WGMS; Reference Haeberli, Cihlar and BarryHaeberli and others, 2000) and Global Land Ice Measurements from Space (GLIMS; Reference BishopBishop and others, 2004) are two international projects aimed at studying glacier fluctuations. Presently, as access to remote-sensing data has become less complex, many researchers have used it to study land-surface objects, but not all researches are as accurate as could be. Recent research based on remote-sensing data revealed significant glacier recession (Reference Khromova, Dyurgerov and BarryKhromova and others, 2003) which seemed speculative and thus required further validation.
The main objective of our research is precise evaluation of the Tien Shan glacial recession using modern remote-sensing techniques in combination with earlier geodetic and meteorological data, and aerial photographs. We therefore compiled >50 years of meteorological data, >100 years of data from land-surface topographical surveys, topographic maps and modern data from on-site global positioning system (GPS) observations, in addition to remotely sensed data from the Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Shuttle Radar Topography Mission (SRTM). Together, these have provided unique information for interpreting climatic events and glaciological consequences. Furthermore, the authors report on numerous field studies they have completed in the Tien Shan in the last 25 years, particularly in the Akhsiirak and Ala Archa areas, which allow validation of past and present glacier positions on the remotely sensed images with great accuracy.
This research did not set itself the task of comparing glacier changes in different mountain systems of the world. Obviously, there are differences and similarities in alpine glacier fluctuation. For example, glaciers in the Alps are more vulnerable and retreat faster than Tien Shan, Pamir or Tibetan glaciers due to lower absolute elevations (Reference Paul, Kääb, Maisch, Kellenberger and HaeberliPaul and others, 2004 ). Glaciers of Tibet, the Himalaya and central Asia can be more stable and may even advance (Reference LiuLiu and others, 2006; Reference Narama, Shimamura, Nakayama and AbdrakhmatovNarama and others, 2006). Therefore, an accurate estimation of the real glacier-covered area in the past and at present is very important in understanding the impact of global climate change on water resources in arid and semi-arid regions.
2. Area of Research
For our analysis we selected two Tien Shan glacierized basins located approximately 800 km from each other that have different climatic regimes (Fig. 1). Both basins have rich glaciological data for the same period of time and have been the focus of many studies since the 19th century (Reference VenukovVenukov, 1861; Reference KaulbarsKaulbars, 1875).
Akshiirak glacierized massif
The Akshiirak glacierized massif (41˚40’–42˚05’ N, 78˚00’–78˚32’ E), the second largest glacierized massif in the Tien Shan, is composed of 178 glaciers with an area of 371.6 km2 (in 2003) and an altitudinal extent of 3600–5000 m a.s.l. (Figs 1 and 2). Large valley glaciers form 87% of the Akshiirak glacierized area. The Akshiirak subcontinental glacier system receives 350 mm of annual precipitation, with 45% occurring in three summer months. Akshiirak glaciers feed two large central Asian rivers: the Narin river, the headwaters of the Syr Dar’ya river, and the Saridjaz river, a tributary of the Tarim river in northwestern China.
Ala Archa glacierized basin
The Ala Archa glacierized basin is located at the northern edge of the Tien Shan mountain system (42˚24’–42˚36’ N, 74˚24’–74˚34’ E) and contains 48 glaciers covering an area of 36.31 km2 (2003). Ala Archa glaciers are spread through an altitude range of 3300–4800 m a.s.l. (Figs 1 and 2). Approximately 83% of the Ala Archa glacier area consists of large valley glaciers, with about 76% of the total glacier-covered areas located between 3700 and 4100m a.s.l.
The Ala Archa glaciers receive 700 mm of annual precipitation, mainly during spring to summer months (48% from April to June) (Reference AizenAizen, 1988a, b; Reference Aizen, Aizen, Glazirin and LoaicigaAizen and others, 2000), and they feed the closed drainage basin of the Chu river, which is the major irrigation and water-supply source for northern Kyrgyzstan and southern Kazakhstan.
3. Data and Methods
Tacheometric surveying, aerial photographs and topographic maps
Akshiirak
The first tacheometric surveying of Akshiirak glaciers was accomplished by expeditions of the Russian Imperial Geographical Society in 1869 (Reference KaulbarsKaulbars, 1875). At the beginning of the 20th century, tacheometric and terrestrial photogrammetric measurements were repeated by expeditions of the Russian Academy of Sciences (Reference Vorob’evVorob’ev, 1935; Reference AvsyukAvsyuk, 1953; Reference Zabirov and KnijnikovZabirov and Knijnikov, 1962). These measurements were applied only to the ablation area of the large valley glaciers in the northwestern part of Akshiirak massif. One of the present authors (V.A. Kuzmichenok), using aerial photographs, carried out the first comprehensive study of Akshiirak glacier area and volume changes between 1943 and 1977. The 1943 and 1977 glacier boundaries were transferred to intermediate 1 : 10 000 scale topographic maps by traditional stereophotogrammetry using the same set of ground-control points (GCPs) that were used to produce the 1 :25000 scale military topographic maps. Two digital elevation models (DEMs) with 100 m resolution were also generated by point-by-point stereophotogrammetric measurements for estimation of surface changes. The DEM vertical accuracies of 2.8 m for 1943 and 2.7 m for 1977 were estimated by repeat measurements and comparison along map sheet boundaries (Reference KuzmichenokKuzmichenok, 1985). The results were published in a 1 :50000 map (Reference KuzmichenokKuzmichenok, 1990a, Reference Kuzmichenokb, Reference Kuzmichenok1991), which contained generalized 1943 and 1977 glacier boundaries and 10 m surface change contour lines. The present positions of seven Akshiirak glacier termini were corrected using large-scale aerial photographs from 1995, acquired from AeroMap US Inc.
Ala Archa
The first tacheometric surveying of the Ala Archa glaciers termini positions were conducted in the mid-19th century (Reference VenukovVenukov, 1861). Several tacheometric surveys also took place during the 20th century (Agricultural Report, 1915; Reference BeznosovBeznosov, 1916; Reference KorzjenewskiKorzjenewski, 1933; Reference VisnevskiVisnevski, 1937; Reference MarichekMarichek, 1950; Reference FreifeldFreifeld, 1952; Reference Il’inIl’in, 1954; Reference Aizen and KalinichenkoAizen and Kalinichenko, 1979; Reference Aizen, Maksimov and SolodovAizen and others, 1983). These measurements did not cover all the glaciers in the studied areas, and therefore were used only for intermediate updating of the position of some glacier boundaries. Thus, complete information about all glacier area changes in the Ala Archa basin was obtained using aerial photographs (1963, 1981) and satellite remotely sensed data for the last 40 years only. An issue of the Glacier Inventory of the USSR (Katalog Lednikov SSSR, 1973) contains information on glacier-covered areas in the Ala Archa basin, which was generated from aerial photographs of 1963. Change in glacier-covered area between 1963 and 1981 was evaluated by Reference AizenAizen (1984).
Graphical renditions of the Akshiirak and Ala Archa tacheometric surveying during the period 1861–1950 were registered to 1 : 25 000 scale modern topographic maps by affine transformation using common relief features (moraines, rocks, etc.) and, in some cases, survey benchmarks determined on site. The estimated transformation accuracy is 30–40m and up to 100 m for the earliest surveys.
Remote-sensing data and DEMs
The present glacier-covered area and the glacier boundaries in Akshiirak and Ala Archa were determined from ASTER images acquired on 17 August (Ala Archa) and 18 August (Akshiirak) 2003 (Fig. 2). Both images were clear of clouds and revealed glacier surfaces near the end of the ablation period when glacier boundaries are most visible. The raw L1A ASTER images were orthorectified with orthobase photogrammetric software. GCPs were collected from the topographic maps of 1977 (Akshiirak) and 1963 (Ala Archa) and identified on the ASTER images. The nominal vertical accuracy of the maps is 3.3 m and horizontal accuracy is 5 m. For the orthorectification, 15 m resolution DEMs were generated for Aksiirak and Ala Archa from the 30 1 : 25000 scale topographic maps that cover all the glaciers and surrounding areas. The 10 m contour lines and spot elevations were manually digitized and processed with the ANUDEM algorithm available in the ArcGIS 9 software package. The orthorectification root-mean-square error (rmse) of the 28 GCPs of the Akshiirak ASTER image is 9 m, and for the 20 GCPs of the Ala Archa ASTER image it is 10 m. To achieve maximum accuracy for the glacier boundaries, manual digitizing on false-color composites of visible/near-infrared (VNIR) bands was used. True hardware-enabled stereo viewing with nadir 3N and backward-looking 3B bands was applied to delineate glaciers in problem areas (debris-covered termini and shadows). The accuracy of digitized 2003 glacier boundaries was confirmed by 2002 on-site GPS measurements on seven glacier termini in the Akshiirak massif.
To estimate the glacier volume changes in Akshiirak from 1977 to 2000 (Surazakov and Aizen, in press), we used an unedited version of 3 arcsec SRTM data acquired in the C-band and processed by NASA–Jet Propulsion Laboratory (JPL), USA (Reference Rabus, Eineder, Roth and BamlerRabus and others, 2003), and the 1977 DEM. For accurate comparison with the 1977 DEM, the SRTM data were transformed from the World Geodetic System 1984 (WGS84) to the ‘Pulkovo’ Russian 1942 datum with 2 m accuracy. Accuracy of the surface change measurements was estimated on flat depressions surrounding Akshiirak massif and glacier-free outcrops. More than 110000 SRTM points were compared to the 1977 DEM to validate the accuracy. The standard deviation of differences is 6.3m on slopes <258. The glacier surface-elevation changes on steep slopes (>258) were linearly interpolated using a triangular irregular network (TIN) with the SRTM points of flat areas and the glacier boundaries as zero change. After correction for penetration of the radar signal into snow and ice (Reference Rignot, Echelmeyer and KrabillRignot and others, 2001), the estimated overall error of glacier surface change from 1977 to 2000 is 8.2 m. The modern rates of surface change were confirmed by comparing SRTM data with 2003 Ice, Cloud and land Elevation Satellite (ICESat) laser altimetry data (Surazakov and Aizen, in press).
4. Results
Central Tien Shan (Akshiirak glacierized massif)
In 1943, Akshiirak glacier covered 424.7 km2. From 1943 to 1977 the glacial area shrank to 406.8 km2 (−4.2%), and it decreased to 371.6 km2 by 2003 (−8.7%). In total, Akshiirak glaciers have lost 12.5% of their surface area from 1943 to 2003 (Table 1a). The rate of area reduction increased from 0.53 ±0.0004 km2 a- 1 during the first period (1943−77) to 1.35 ± 0.001 km2 a- 1 during the second period (1977−2003). The rate of glacier thinning also changed from 0.24± 0.12 m a- 1 to 0.69± 0.37 m a- 1. The volume of Akshiirak glaciers has been reduced 9.7±0.01 km3 since 1943 (Table 1b), which is approximately equal to 29% of initial glacier volume estimated as difference between the glacier surface elevation in 1943 and the radio-echo sounding glacier thickness measurements in 1986−87 (Kuzmichenok, 1990, Reference Kuzmichenok1996).
However, the glacier recession during the first period was not uniform. For example, 7 of the 178 Akshiirak glaciers advanced between 1943 and 1977 and an increase in surface elevation was observed on 32 glaciers. In the second period (1977–2003), all Akshiirak glaciers experienced recession. The four glaciers with areas <1 km2 completely disappeared by 2003, and one glacier has been removed by Kumtor Mining Co. to access a gold-bearing lode (Figs 2–5).
We cannot represent all Akshiirak massif glacier behavior in this paper, so only a few of the 178 glaciers were selected for discussion (Table 2; Figs 2 and 3):
Bordoo glacier (1.25 km2) advanced 270 m from 1943 to 1955 and then retreated 700m by 2003.
Sari Tor glacier (3.05 km2) retreated 70 m from 1932 (Reference Kuzmichenok and KasenovKuzmichenok and Kasenov, 2002) to 1943. From 1943 to the mid-1970s this glacier was stable but it retreated 310m from 1977 to 2003.
Davidova glacier (10.75 km2) terminus was advancing before 1932 (Reference Kuzmichenok and KasenovKuzmichenok and Kasenov, 2002), but by the mid-1940s it was retreating (Figs 3 and 5). The second advance of the glacier terminus was recorded in 1978 when the glacier terminus developed a 40m high steep front edge and advanced 30 m relative to 1977, while its surface elevation increased 50 m (Fig. 3b). However, the glacier terminus did not reach the 1956 margin. From the end of the 1980s the retreat of the glacier terminus accelerated, particularly on the right side of the tongue (down-ice view), and by 2000 the surface elevation of the glacier tongue had decreased 70 m. By the mid-1990s, the right side of the glacier terminus had been buried under a bank of waste from the gold mining identified by SRTM data as an increase in surface elevation (Fig. 5).
Lysyi glacier (4.17 km2) has retreated 690m during the past 70 years, and the average surface elevation lowered 32 m from 1943 to 2000. The surface of this glacier has also been covered by a bank of waste from the gold mining and by blocks of ice that were removed from a hanging glacier to accommodate gold mining.
Petrova glacier (65.33 km2) is one of the largest glaciers in the Akshiirak massif, and its glacier terminus descends to its moraine lake. From 1869 to 2003, this glacier gradually retreated >2.5 km, and the surface of the glacier terminus has lowered >100 m. The average surface elevation lowered 19.3 m from 1943 to 2000.
Dvoinoi-Levyi glacier (2.1 km2) receded 380 m from 1943 to 2003, separating from Petrova glacier, and the average surface elevation lowered 24.1 m from 1943 to 2000.
Bezimyanniy glacier (4.83 km2) has exhibited some surging behavior. The glacier terminus stagnated before 1940, and by the mid-1950s it had advanced 340 m. Since then it has receded. By 2003 the elevation of the middle part of the glacier terminus had risen 35 m, although the elevation of the area above this point decreased 40–45 m (Fig. 4).
Northern Tien Shan (Ala Archa glacierized basin)
From the mid-19th century to the beginning of the 20th century, the termini of the Ala Archa glaciers retreated 1.0 km on average from <2800 m a.s.l. to 3100 m a.s.l. (Reference VenukovVenukov, 1861; Reference LangwagenLangwagen, 1908). They then stagnated approximately at the same level until 1910–15 (Agricultural Report, 1915; Reference KorzjenewskiKorzjenewski, 1933). However, these results were based on information for large glaciers. According to the tacheometric surveying (Reference VisnevskiVisnevski, 1937; Reference MarichekMarichek, 1950; Reference FreifeldFreifeld, 1952), from 1915 to the end of the 1940s the glaciers retreated another 0.5 km, reaching 3200 m a.s.l. Comparing these data with the aerial photographs of 1963, we found that Ala Archa glacier margins were stationary during this period. In 1963, Ala Archa glaciers covered 42.83 km2. By 1981, the glacier area had shrunk to 40.62 km2 (5.2% from 1963) and glacier termini had retreated to 3500 ma.s.l. By 2003, the glaciers had lost another 10.6% of their area, shrinking to 36.31 km2 (Table 3; Fig. 6). Table 4 presents the recession of the terminus of Golubina glacier, the largest glacier in the Ala Archa basin, from 1861 to 2003.
Golubina glacier (5.8 km2) occupied an altitudinal range of 3300–4400ma.s.l. (2003) and has a north-northwest orientation (Figs 6 and 7). In 1861, the Golubina glacier terminus lay in a deep canyon at 2900 ma.s.l. (Reference VenukovVenukov, 1861). From 1861 until 1913, the glacier terminus retreated to 3050 m a.s.l. In 1949, the glacier terminus was at 3150m a.s.l and by 1963 it reached 3250m a.s.l. Between 1958 and 1982, when glaciological observations and tachymetric surveys on Golubina glacier were conducted on a regular basis, small annual oscillations of the terminus position up to 80 m occurred, and growth of the glacier surface elevations of 25–40m occurred along the central axial line (1971, 1974, 1976) (Reference Aizen, Maksimov and SolodovAizen and others, 1983). From 1958 to 1972, the Golubina glacier mass balance was positive (Reference AizenAizen, 1988b), but since 1973 its mass balance has been predominantly negative (Reference Aizen, Maksimov and SolodovAizen and others, 1983, 1988a, b). Consequently, the glacier began to retreat more rapidly, and from 1981 to 2003 it retreated 260m and thinned 25–30m on average along the longitudinal section, while glacier area shrank 0.2 km2.
Thus, from 1861 to 2003, glaciers in the Ala Archa basin retreated 1.5 km on average. In the last 40 years from 1963 to 2003, the total Ala Archa glacier-covered area shrank 6.52 km2, or 15.2%. From 1963 to 1981, nine small glaciers <1 km2 in area totally disappeared (Reference AizenAizen, 1984).
5. Discussion
Our research presents a long-term (about 140 years) estimation of the northern and central Tien Shan glacier changes based on a complex of geodetic ground surveys and remote-sensing data, which quantified the spatial and temporal regime of the glaciers (Tables 1–4). There has been a definite trend of glacier recession over the last 140 years and especially since the mid-1970s, which indicates an abrupt climate-change effect. The Tien Shan, located in the center of the large Eurasian continent, have spring–summer precipitation maxima, which preserve glaciers from intensive melt during the ablation period by increasing the glaciers’ surface albedo. At the same time, the main factor controlling the glacier regime is the impact of air temperature which affects the type of precipitation, the duration and the intensity of snow and ice melt throughout altitudinal belts. The modern increase of air temperature, which is also observed in the Tien Shan’s alpine areas, extends the period and intensity of melt and the glacier recession (Reference Aizen, Aizen, Melack and DozierAizen and others, 1997).
Analysis of data from the station located closest to the Ala Archa glaciers shows no significant trends in observed annual precipitation or warm-season (May–September) temperature for the period 1913–2000 or 1963–2000 at Baitik station (1580 m a.s.l.), Ala Archa basin (Fig. 8a). However, a significant negative trend in annual precipitation is observed during the period 1963–81 (29% of average), while air temperature significantly increased during the 20 year period 1981–2000 (0.93˚C) (Fig. 8a; Table 3), which extends the period of glacier ablation and consequently accelerates glacier recession (Fig. 8c; Table 4).
Analysis of data from the station located closest to the Akshiirak glaciers shows there is no significant trend in annual precipitation in the central Tien Shan for the last 60 years (1943–2003) of instrumental observations (Tien Shan station (3614 m a.s.l.)) (Fig. 8b). However, warm-season air temperature increased by 1.0˚C for the same period (Table 1). The increase in air temperature, especially observable since the 1970s, accelerated the recession of the Tien Shan glaciers (Fig. 8b; Table 1).
Our results, showing a decrease of 8.6% in the Akshiirak glacier area from 1977 to 2003, are in accordance with Reference LiuLiu and others (2006) and Reference Narama, Shimamura, Nakayama and AbdrakhmatovNarama and others (2006) who estimated the glacier recession in the northern and eastern Tien Shan between 1970 and 2004. However, our results do not correspond with the –23% estimate by Reference Khromova, Dyurgerov and BarryKhromova and others (2003), particularly for the Akshiirak glacier area for the same period. We suppose that several methodological issues may contribute to the discrepancy in the results: (1) Reference Khromova, Dyurgerov and BarryKhromova and others (2003) used a map developed by Reference KuzmichenokKuzmichenok (1990a, b) to estimate glacier area changes from 1977 to 2001. However, the 1990 map (Reference KuzmichenokKuzmichenok, 1990a,Reference Kuzmichenokb) cannot be used for such analysis due to restrictions on large-scale thematic maps published for civilian use in the former USSR: there is no coordinate information on the map, and accurate georegistration is not possible. Also, during compilation of the 1 : 50 000 map from original 1 : 10 000 maps, significant non-systematic horizontal errors were introduced across the map. (2) Manual delineation of glaciers is subject to personal judgment as mentioned by Reference Khromova, Dyurgerov and BarryKhromova and others (2003) concerning ‘thin ice and snow patches’. To minimize this type of error, in our study the same person, one of the authors, interpreted the 1943, 1977 and 2003 data. (3) In the ASTER image used in Reference Khromova, Dyurgerov and BarryKhromova and others’ (2003) estimation, glaciers are almost completely covered by fresh snow, which complicates glacier boundary identification. (4) There was no assessment of the glacier area covered by debris. The six large Akshiirak glaciers have significant debris cover (around 1 km2 in total) that we identified on large-scale aerial photographs of 1943 and 1977, and by stereo viewing of 3N and 3B bands of the ASTER image of 2003.
6. Conclusion
By using ground measurements and remote-sensing data, we have determined that continuous glacial recession has occurred over the last 140 years in the Ala Archa and Akshiirak glacial regions of the Tien Shan. The glaciers retreated up to 3 km between the 1860s and 2003, the surface of their ablation areas lowered >100m, and the glacier area shrank by 15.2% in the Ala Archa over the last 40 years. The Akshiirak glacierized massif lost about 10 km3 of glacier ice between 1943 and 2000. However, the overall trend of glacier recession is variable. On a broader view, the fact that most of the Tien Shan glaciers have thinned over the study period does support the observation that, in general, glacial retreat accelerated from the mid-1970s to the present. The rate of glacier recession is about 3% higher in the northern peripheral Tien Shan ranges than in the central Tien Shan, which may be explained by severe continental climate and exclusively summer maximum precipitation in the central Tien Shan, while maximum precipitation in the northern Tien Shan is in spring and the first month of summer. Hence, continuous increase of warm-season air temperatures in the northern and central Tien Shan since the mid-1970s without increase of precipitation may further accelerate glacier recession and intensify desertification processes in central Asia and northwestern China.
Acknowledgements
This work is a part of research supported by the US National Science Foundation grant (ATM-0233583). The authors express their sincere gratitude to D. Kasenov, K. Aitimbetov and R. Piligrim (Kumtor Mining Co.) for assistance in organizing the 2002 field research, and H. Brecher for careful editorial work and valuable comments and suggestions.