Introduction

Most glaciers in the Tien Shan have retreated at accelerating rates over the 20th century (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006; Reference Li, Zhu, Zhang, Pei, Qin and ZhouLi and others, 2006; Reference Kutuzov and ShahgedanovaKutuzov and Shahgedanova, 2009). This retreat is particularly evident in the northern Tien Shan ranges, which are located at the transition between the mild temperate conditions of southern Siberia and the more severe continental climate of High Asia. Retreating glaciers have been reported from the central Tien Shan regions too, but with lower recession rates. Among the glacierized ranges, the Akshiirak massif has attracted the most attention because of its extensive glaciation, comprising 178 glaciers with a total area of 317.6 km² (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006). Although most of the glaciers retreated over the 20th century, a prominent advance including surge-type movements has been reported for subperiods in the massif (Reference Dolgushin and OsipovaDolgushin and Osipova, 1982; Reference Solomina, Barry and BodnyaSolomina and others, 2004). The overall balance of the glacierized area was negative in the last century: 82% of glaciers retreated, 14% were in a steady state and only 4% advanced (Reference Dyurgerov, Liu and XieDyurgerov and others, 1995). The rate of glacierized area reduction was relatively low from 1943 to 1977 (~0.5 km² a^-1) but increased to ~1.4 km² a^-1 in 1977-2003 (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006).

Whereas the distribution and extent of glaciers in the central Tien Shan is relatively well known, information about the recent glacier fluctuations is comparatively poor, with the exception of the Akshiirak massif. The first record of glacier recession in this region was reported for the largest valley glacier in the Akshiirak massif, named Petrov by Reference KaulbarsKaulbars (1875) after the second lieutenant of his expedition. In 1869, Kaulbars described the area around the glacier terminus with a small lake in front of the terminal moraine. Davydov visited the area in 1911 and drew a sketch map of the glacier terminus with a lake within the terminal moraines. The outline is a rough approximation of the site and unfortunately does not help us to reconstruct the glacier terminus position. Reference DavydovDavydov’s (1927) report from his second visit to the area suggests that Petrov glacier was in a steady state between 1911 and 1925. Reference Kalesnik and EpsteinKalesnik and Epstein (1935) noted that the glacier terminus was located by the third inner ridge of the terminal moraine complex in 1932. Relics of this moraine probably form a part of the recent shoreline at the lake dam. An aerial survey carried out on behalf of the central government of the Soviet Union in Central Asia on 16 September 1943 revealed a retreat of Petrov glacier which terminated in a proglacial lake. According to Reference AvsiukAvsiuk (1953), the length of the glacier had decreased by 1.4 m a^-1 since 1911. However, this interpretation is very tentative as short glacier advances could have occurred during this period (Reference BondarevBondarev, 1963). Aerial photographs taken in 1957 show further recession of the glacier terminus, which retreated by about 400 m (Reference BondarevBondarev, 1963). Between 1957 and 1977 the length of the glacier decreased considerably (Reference Sevast’yanov and FuntikovSevast’yanov and Funtikov, 1981; Reference KuzmichenokKuzmichenok, 1990). Fluctuations of the glacier at the end of the 20th century were described briefly by Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others (2006) and Reference Janský, Engel, Šobr, Beneš, Špacek and YerokhinJanský and others (2009). However, no field observation of the glacier has been available since the beginning of the 21st century.

This paper presents results of monitoring changes in ice- cover extent of Petrov glacier in 2003-09. A synopsis of the historical evolution of this site was published by Reference Janský, Engel, Šobr, Beneš, Špacek and YerokhinJanský and others (2009). In this paper, we present results of geomorphological and bathymetric investigations carried out on Petrov glacier and Petrov lake from 2006 to 2009. Our first aim is to reconstruct the recent fluctuations of Petrov glacier. The second objective is a detailed description of recent changes in the proglacial area including volume changes of Petrov lake.

Study Area

Petrov glacier and its proglacial lake are situated in the Akshiirak massif (41°40′-42^°05′N; 78^°00′-78°35′E; 938 km²), central Tien Shan. The massif consists of three mountain ranges with a southwest orientation. The ranges are elevated 1200-1400 m above the surrounding relief. Well-developed glacial landforms and extensive glaciation covering 371.6 km² (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006) are characteristic features of the massif. The study area is located on the northwestern slope of the Akshiirak massif, which rises above a gently sloping valley of the Kumtor river. Petrov glacier terminates in a wide valley with relatively low and slightly inclined slopes. The proglacial Petrov lake (3741 m a.s.l.) stretches 2.9 km from the glacier terminus towards the west-northwest, where it is dammed by a system of terminal moraines. A surface outflow is situated in the northern part of the moraine (Fig. 1). The moraine dam comprises till of three to five stages of glaciation, with the most prominent stage correlated with the Little Ice Age period (Reference BondarevBondarev, 1963; Reference SavoskulSavoskul, 1995). A number of thermokarst lakes in the dam document the presence of massive ice in the recently exposed moraines. Degradation of the buried ice is accompanied by relief sinking, mass movements and potholing (Reference Janský, Engel, Šobr, Beneš, Špacek and YerokhinJanský and others, 2009). On the distal side of the moraine complex the glaciofluvial sediments cover the bottom of the Kumtor valley, which is up to 1.5 km wide. From the north, a tributary hanging valley runs into the right slope of the main valley 100 m above Petrov lake. The slopes of the valley around the lake consist of Carboniferous limestone and sandstone, while the eastern part of the valley consists mainly of Lower Carboniferous gneiss and schist (Reference IvanovIvanov, 1964). Proterozoic granitoids surround the accumulation area of Petrov glacier in the central part of the Akshiirak massif.

Fig. 1. Location map of glaciers in the study area, showing glacier numbers and flowlines. The dark-grey textured glaciers drain into the Kumtor river, whereas glaciers plotted in white are situated beyond the watershed. Distribution of glaciers after Reference KhromovaKhromova (2005), numerical data after Reference SuslovSuslov (1973).

The glacierized area is dominated by Petrov glacier, which is a large dendritic valley glacier (Fig. 1). The glacier is characterized by an accumulation-area ratio value of 0.6, long constituent flow units, low average slope angle (<10°) and an even longitudinal profile. The glacier is not considered a surge-type glacier, unlike Northern Karasai or Bordu in the southern part of the Akshiirak (Reference BondarevBondarev, 1963; Reference Dolgushin and OsipovaDolgushin and Osipova, 1982). Two small valley glaciers fill the bottom of deeply incised tributary valleys located north of Petrov lake (Fig. 1;Table 1). These glaciers have irregular longitudinal profiles, higher slope angles (10-15°) and shorter response times than Petrov glacier. The equilibrium line-altitude (ELA) in the region ranged from 4150 to 4340 m a.s.l. between 1985 and 1989 (Reference Dyurgerov, Liu and XieDyurgerov and others, 1995).

Table 1. Morphologic characteristics of glaciers in study area

The study area falls within the zone of mountain permafrost. The lower limit of continuous permafrost in the Akshiirak region is located at 3500-3600 m a.s.l. (Reference Dikich, Bakov, Koshoev, Melnikova and DikichDikich and others, 1991; Reference Gorbunov, Seversky and TitkovGorbunov and others, 1996). Based on ground temperature measurements in nearly 50 boreholes, the permafrost is up to 370 m thick (Reference Ermolin, Nemov and PopovErmolin and others, 1989). Since the lower limit of continuous mountain permafrost is well below the ELA, we expect that glaciers in the study area are mostly cold-based. The sporadic data on thermal regimes of glaciers in nearby ranges (Reference Dyurgerov, Liu and XieDyurgerov and others, 1995) support this view.

The region is characterized by a continental climate with marked seasonal variation. Climate conditions are affected by the position of the Akshiirak massif, which represents an obstacle for western transport of moist Atlantic air masses. In the zone of predominantly westerly winds, precipitation mainly falls in the spring and summer months, while the winters are cold and dry. Mean annual precipitation (MAP) reaches 291 mm, with 45% occurring during the more humid season from June to August (Tien Shan meteorological station data, 1930-97). Mean January, July and annual air temperatures (MAAT) are -21.7°C, 4.5°C and -7.6°C, respectively (1930-97). According to continuous meteorological data collected before the replacement of the Tien Shan meteorological station in 1997, air temperature has increased since 1930, but no significant trend in precipitation has occurred (Reference Solomina, Barry and BodnyaSolomina and others, 2004; Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2007).

Materials and Methods

Remote-sensing data and maps

Aerial photographs and satellite images were used to determine glacier terminus changes since the mid-20th century (Table 2). The paper prints of the aerial photographs and a sheet of 1 : 25 000 topographic map based on aerial survey from 1977 were scanned at 600 dpi and registered using coordinates of control points. A Trimble GeoXH GPS receiver (Reference TrimbleTrimble, 2009) was used to determine coordinates for the points. Averaging and accuracy-based logging were applied to ensure horizontal root-mean-square (rms) position accuracy <3.5 m. The Landsat and QuickBird data were georeferenced using 7-11 control points. The pixel resolution of the data and registration rms errors are summarized in Table 2. The uncertainty associated with glacier terminus changes between two dates was calculated as the sum of the resolution-related uncertainty and uncertainty due to errors in registration of dataset pairs (Reference Li, Zhu, Zhang, Pei, Qin and ZhouLi and others, 2006). The partial uncertainties are estimated using

Table 2. Overview of maps, aerial photographs and satellite imageries used in the study

where R _a and R _b are the resolution or registration errors of images a and b, respectively (Reference Li, Zhu, Zhang, Pei, Qin and ZhouLi and others, 2006). The uncertainties in the area extent are determined using

where A = x ², x is the image resolution and d is the uncertainty related to the spatial resolution of the dataset (Reference Hall, Bayr, Schoner, Bindschadler and ChienHall and others, 2003). The calculated uncertainties are listed in Table 3.

Field mapping and geophysical soundings

Combined GPS/geodetic mapping of glacier limits and landforms in the proglacial area was carried out on 411 July 2006, 12-13 August 2008 and 5-9 September 2009. A Trimble GeoXH GPS was used to delimit accessible boundaries of the glaciers and to determine coordinates for the base points of geodetic measurements with rms accuracy <3.5 m. The calving terminus of Petrov glacier was located using a pulsed laser rangefinder of a Leica TCR 705 total station in non-prism mode. Though the declared accuracy of non-prism mode measurement is 3 mm (+ 2 ppm × measuring distance) the overall accuracy of measurement was determined by the GPS positioning of the base points. The collected positional data were analysed in conjunction with remote-sensing imagery using a GIS (MapInfo Professional 10.0) in order to quantify length and area fluctuations of Petrov glacier and its lake.

Ground-penetrating radar (GPR) was used at the terminus of Petrov glacier to acquire information on glacier thickness and subglacial topography. GPR data were collected along a profile that was perpendicular to the centre flowline of Petrov glacier and lay on average 130 m from the glacier terminus. The location of the profiles was chosen so that measurements described the whole width of the glacier. An unshielded 50MHz RTA (Rough Terrain Antenna) and MAL? CU-II control unit (MALå GeoScience, 2005) were used for data collection. The signal acquisition time was set to 860 ns, and scan spacing was about 0.07 m. GPR soundings were done on 8 September 2009 in favourable dry conditions. The collected data were processed and interpreted using the Reflex-Win software version 4.5. The depth axis of profiles was calculated from the time axis and the wave velocity for glacier ice (0.167 m ns^-1). The permittivity for glacier ice corresponds with the value of 3.17 (Reference Davis and AnnanDavis and Annan, 1989).

Table 3. Uncertainty associated with the Petrov glacier terminus change

Results

Fluctuations of Petrov glacier

Two major periods of Petrov glacier fluctuations may be distinguished since the first investigations in 1869. From 1869 to 1977, the terminus of Petrov glacier seems to have undergone relatively small annual oscillations showing similar changes in both northern and southern sections (Fig. 2). During this period, the glacier retreated at a mean rate of 13.6m a^-1. At the end of the 1970s, the recession accelerated substantially and the land-based and calving sections of the Petrov glacier terminus retreated at different rates. Since 1977, linear retreat rates have not dropped below 33.6 ± 1.9 m a^-1, with a mean value of 40.6 ± 0.5 m a^-1. The surface area of the glacier terminus decreased at a rate of 0.07 ± 0.0002 km² a^-1 over the period 1977-2009.

Fig. 2. Positions of the Petrov glacier terminus over the 20th century. The map documents the expansion of the proglacial lake from 1943 to 2003. QuickBird image (2003).

The last three decades may be divided into five periods: three periods (1977-80, 1990-95 and 2008-09) marked by a rapid retreat of the glacier terminus, and two intermediate periods during which the rate of recession remained almost constant (Fig. 3). Between 1977 and 1980 the mean retreat rate for Petrov glacier reached 51.6 ±5.4 m a^-1, which is the second highest value over the period 1869-2009 (Table 4). The more intensive fluctuations occurred in the calving part of the glacier, with a rate of recession exceeding 66.1 ± 5.4 m a^-1. Between 1980 and 1990 the retreat rate of Petrov glacier decreased to 34.9 ±4.1 ma^-1, and the southern flow unit retreated at a mean rate 1.5 times higher than the northern flow unit. The recession of Petrov glacier increased considerably from 1990 to 1995, reaching the highest recorded retreat rate (54.1 ± 8.2 m a^-1). A subsequent period is marked by a relatively steady retreat rate, which increased again after 2008 for the northern flow unit. Over the period 1995-2009 the northern terminus of Petrov glacier retreated at a rate up to three times higher than the southern terminus.

Fig. 3. Mean annual variations of the Petrov glacier length (solid line) and area (grey shading) since 1943. Dotted and dashed lines represent length changes of the southern and northern flow units, respectively.

The increased retreat rate of Petrov glacier in the last two decades coincides with observed changes in adjacent valley glaciers N. 371 and N. 372. These glaciers retreated from 1943 to 1977 by as much as 250 m and slightly advanced in the 1980s (Reference MeinersMeiners, 1997). Between 1990 and 2009 the lengths of the glaciers decreased by up to 100 ± 40 m, representing a mean retreat rate of 5 ± 2ma^-1. The retreat rate accelerated during the period 2006-09 to 11.0 ± 1.6 m a ¹, which is the highest observed rate since 1957.

Table 4. Length and area changes of Petrov glacier and its constituent flow units

Table 5. Area and volume changes of Petrov lake

Petrov lake enlargement

A small proglacial lake that first appeared in the 1930s extended slowly until the late 1970s (Table 5). After 1978 the surface enlargement increased and remained relatively high for the next three decades. Slightly higher surface changes occurred in the periods 1978-80 and 1995-2003, when the lake area increased yearly by ~0.1 km². The progress of the lake volume changes had a similar pattern to the largest increase after 1978. An accelerated increase in the volume was recorded in 2003-06 when the glacier terminus retreated into the deepest part of the lake basin. After a minor decrease of the volume change during 2006-08, the lake grew by (2.2 ± 0.009) × 10⁶ m³ from 2008 to 2009.

The eastern shoreline of the lake changed considerably between 2003 and 2006. The southern section of the Petrov glacier terminus retreated from the lake, and a till sheet and glaciofluvial fan became the main landforms of the lake shore. In the same period, the northern part of Petrov lake expanded into the deepest part (70 m) of the lake basin (Reference Janský, Engel, Šobr, Beneš, Špacek and YerokhinJanský and others, 2009). The bottom of the basin is wide and relatively flat in this area (Fig. 4). The retreat of Petrov glacier after 2006 led to an expansion of the lake into a shallower part of the basin in front of the calving terminus, which is characterized by convergence of contours towards the central line.

Fig. 4. Evolution of Petrov lake between 1911 and 2009. Dark and light grey areas represent lakes and glaciers, respectively. Dotted textures delimit moraines. Modified after Reference DavydovDavydov (1927), Reference AvsiukAvsiuk (1953) and Reference Sevast’yanov and FuntikovSevast’yanov and Funtikov (1981).

Morphology of the Petrov glacier terminus

The different behaviour of the northern and southern sections of the terminus of Petrov glacier has affected the shape and surface morphology of the terminus. The steep calving terminus has the form of a cliff with numerous crevasses that extend hundreds of metres back into the glacier (Fig. 5). The rugged surface indicates relatively high activity of the flow unit which is manifested in frequent calving. The higher mass turnover led to the rapid retreat of the calving terminus and to the ice front lowering by about 5 m between 2006 and 2009. In contrast, the landterminating glacier front is marked by predominantly moderate slope angle and lower crevasse density. The smoother ice surface may be directly related to a lower tensile stress and mass exchange on the southern flow unit. The marginal parts of the terminus were frozen to the ground, as observed in subglacial cavities in central and southern parts of the terminus. The GPR data show a distinctive reflection boundary corresponding to an ice-bed boundary of Petrov glacier. The amplitude of reflections is smaller in the area below the northern flow unit, implying transition from ice to subglacial till rather than to solid bedrock. Applying a wave velocity of 0.167 mns^-1 for the depth calculation, the thickness of the northern flow unit exceeds 65m and the maximum depth of the southern flow unit is 38 m (Fig. 6).

Fig. 5. Terminus of Petrov glacier in July 2005. The northern part of the glacier borders the eastern shore of Petrov lake whereas the southern flow unit terminates on the former glacier bed 30–70m from the shoreline. Photo: Z. Engel.

Fig. 6. A 50MHz GPR section for the terminal part of Petrov glacier. The most prominent reflection corresponds to the lower limit of the glacier.

The effects of uneven ablation were observed on the glacier terminus during 2005-09. The repeated measurements of the terminus position confirmed a remarkable retreat and surface lowering of the glacier sections exposed to direct insolation, which were 130-240 m further up the valley, and up to 5 m lower, in 2009 than in 2006. In contrast, debris-covered parts of the glacier terminus retreated by up to 135 m, and the surface lowering was smaller than 2 m according to repeat elevation measurements. The thickness of the debris cover ranged from tenths of a metre to about 1.5 m. The total amount of iceberg calving from the northern glacier terminus accounted for the total mass loss of (17.7 ± 0.4) × 10⁶ m³ from 2003 to 2009. The mean annual calving flux increased from (2.5 ± 0.06) × 10⁶m³ in 2003-06 to (4.6 ± 0.09) × 10⁶m³ in 2008-09. Subaerial calving was observed during our fieldwork in 2006 and 2009, when small sections of the cliff were calved off once every day or two, forming icebergs up to 8 m long. The presence of a number of larger icebergs along the cliff indicates that subaqueous calving may contribute to the disintegration of the glacier terminus.

Climate and discharge variations

Analyses of meteorological data from the Tien Shan station show a statistically significant increase in WST (r =0.323, p = 0.06) and decrease in MAP (r =-0.343, p = 0.04) in the period 1930-97. In contrast, no significant trend was observed in the period 1998-2009. The air-temperature record indicates that WST increased by 0.8°C between the periods 1930-97 and 1998-2009 (Fig. 7). At the same time, MAAT increased by 1.8°C, and the period of positive air temperatures extended from 3.6 to 4.2 months. However, the inferred increase in temperature is complicated by the relocation of the meteorological station. MAP increased by 60mm between the two time series. The increase of precipitation was limited to the warm season only. During the remaining months, precipitation slightly decreased.

Fig. 7. Mean annual precipitation (bars), mean annual (bold line) and warm-season temperature (grey line) at the Tien Shan meteorological station. Dotted line shows the annual retreat rate of Petrov glacier. MAP, WST and MAAT values are reported for periods 1930–97 and 1998–2009.

Discharge recordings of the Kumtor river have not revealed any trend for the period 2001-09 (Fig. 8). A positive correlation (r = 0.74, p = 0.02) between discharge and WST exists, but no statistically significant correlation was found between discharge and precipitation. A maximum discharge is seen in the late warm season (August), 1 month after the warmest month and 2 months after the maximum precipitation (Fig. 9). The delay of the glacier discharge results from the retaining of precipitation in the watershed (net glacier storage) during the accumulation season.

Fig. 8. Annual discharge of the Kumtor river (Q), precipitation (P) and warm-season temperatures (WST) at the Tien Shan meteorological station 2001–09.

Fig. 9. Precipitation (bars) and monthly discharge for the Kumtor drainage area.

Discussion

The rapid retreat of Petrov glacier in the period 1943-57 described above is consistent with intense recession of large dendritic glaciers in the Tien Shan. In the Akshiirak massif, most of the valley glaciers retreated by up to 50 m a^-1 during this period (Reference BondarevBondarev, 1963). Increased rate of glacier retreat was also reported from the nearby Terskey-Alatoo (Reference Kutuzov and ShahgedanovaKutuzov and Shahgedanova, 2009), the Khan Tengri massif (Reference KoshoevKoshoev, 1986) and the Kyrgyz Range (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006). From 1957 to 1977, the retreat rate of Petrov glacier decreased, coincident with reduced recession of glaciers throughout the northern and central Tien Shan (Reference Makarevich, Liu, Dyurgerov, Liu and XieMakarevich and Liu, 1995). According to Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others (2006), mean retreat rates of glaciers in the Akshiirak ranged from 5 to 19 m a^-1 over this period. Despite the overall glacier recession in the 1960s and the first half of the 1970s, many glaciers in the northern Tien Shan were in equilibrium, or even advancing slightly (Reference MaksimovMaksimov, 1980; Reference Yermolov, Solodkov and FirsovYermolov and others, 1986; Reference Vilesov, Makarevich and PolyakovVilesov and others, 1993). In contrast, the recession of ?r?mqi glacier No. 1, eastern Tien Shan, was most rapid from the 1960s to the mid-1970s, when the retreat rate decreased substantially (Reference Yao, Wang, Liu, Pu, Shen and LuYao and others, 2004). These differences result from more variable climate conditions in peripheral mountain areas (Reference Aizen and AizenAizen and Aizen, 1994).

In the late 1970s, a period of minor glacier recession in the study area changed to a rapid retreat. The observed rapid retreat event of Petrov glacier supports the fluctuation scheme described by Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others (2006) in the Kyrgyz Range where Golubina glacier retreated by 18ma^-1 in 1976-78. Although the glacier retreat in the study area had decreased in 1980, it remained significantly higher than in the period 1957-77. The rapid glacier retreat observed in the study area in the 1980s and 1990s coincides with increased glacier recession that was reported from numerous Tien Shan ranges (Reference Aizen and AizenAizen and Aizen, 1994; Vilesov and Uvarov 2001; Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006; Reference Li, Zhu, Zhang, Pei, Qin and ZhouLi and others, 2006; Reference BolchBolch, 2007; Reference Niederer, Bilenko, Ershova, Hurni, Yerokhin and MaselliNiederer and others, 2008) including the Akshiirak massif (Reference Khromova, Dyurgerov and BarryKhromova and others, 2003;Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006) and nearby Terskey-Alatoo (Reference Narama, Shimamura, Nakayama and AbdrakhmatovNarama and others, 2006; Reference Kutuzov and ShahgedanovaKutuzov and Shahgedanova, 2009). The highest values of the Petrov glacier retreat rates are in the 1990s, which is consistent with the observations made in the Akshiirak massif by Reference Khromova, Dyurgerov and BarryKhromova and others (2003). Accelerated recession was also reported from the northern Tien Shan where the annual reduction of glacier area doubled after 1986 (Reference Niederer, Bilenko, Ershova, Hurni, Yerokhin and MaselliNiederer and others, 2008). In the central and eastern Tien Shan, the recession rate also increased but not so much (Reference Li, Zhu, Zhang, Pei, Qin and ZhouLi and others, 2006; Reference Narama, Shimamura, Nakayama and AbdrakhmatovNarama and others, 2006). Since the beginning of the 21st century, the glaciers in the study area have increasingly retreated. The described trend is in agreement with observations of Narama and others (2010b), who reported slightly higher recession rates in 2000-07 than in previous decades within the northern and central Tien Shan.

A number of studies indicate that the spatial and temporal variations of glaciers in the Tien Shan over the 20th century were controlled by climate changes. However, there is disagreement about the main climate factor that controls glacier fluctuations. With the exception of Reference BondarevBondarev (1963), who related glacial fluctuations primarily to changes in annual precipitation, authors mostly considered air temperature as the main factor controlling glacier mass balance (Reference BakovBakov, 1983; Reference Glazirin and FaizrachmanovGlazirin and Faizrachmanov, 1988; Reference Aizen and AizenAizen and Aizen, 1994; Reference Solomina, Barry and BodnyaSolomina and others, 2004; Reference Han, Ding, Ye, Liu and JiaoHan and others, 2006). According to Reference Aizen, Aizen, Melack and DozierAizen and others (1997), an increase in air temperature extends the period and intensity of melt and increases glacier recession. Recent studies of the climate data from the Tien Shan station showed that air temperatures observed between May and September were the main factors controlling the fluctuations of Akshiirak glaciers in the second half of the 20th century (Reference Aizen, Kuzmichenok, Surazakov and AizenAizen and others, 2006). Results presented in this study confirm this view, as the WST from 1943 to 1995 correlates with the Petrov glacier retreat rate (r=-0.922, p =0.03). A strong negative correlation implies the importance of the ablation-season warming for glacier retreat. However, the glacier fluctuations were also closely interconnected with precipitation (r =0.943, p = 0.02). The statistically significant correlation between the climate data and glacier retreat was not found for the period 1999-2009. The lack of a significant relation results from the short time series and limited number of described recession periods.

Though the changes of the Petrov glacier length are controlled in particular by climate, differences between the retreat rates of the northern and southern sections of the glacier terminus imply modification of the fluctuations by local environmental factors. According to our results, the differences in retreat rates started to increase during the period 1999-2003 and increased considerably after 2003. The bathymetric map (Fig. 4) shows that the glacier terminus has retreated under different conditions since the 1990s. The northern flow unit started to retreat from a relatively shallow and less inclined part of the lake basin into its entrenched deepest part. In contrast, the southern glacier section retreated back into the shallow eastern end of the lake basin and became land-terminating by 2003. As the southern terminus retreated from the lake and started to move back across the till sheet, its recession decreased substantially. Since the end of the 1990s, the average retreat rate of the southern land-terminating terminus has been two to three times lower than the retreat rate of the northern terminus (in a proglacial lake). The contrast in retreat rates is further enhanced by the fact that the southern flow unit is about half as thick as the northern flow unit and is partly frozen to its bed. The difference between the retreat rates of the land-based and calving glacier sections is in agreement with the notion that calving into proglacial lakes accelerates glacier advances and retreat rates (Reference Warren and AniyaWarren and Aniya, 1999; Reference Naruse and SkvarcaNaruse and Skvarca, 2000; Reference Tsutaki, Nishimura, Yoshizawa and SugiyamaTsutaki and others, 2011). Similar differences of retreat rates were described in the nearby Terskey-Alatoo, where a number of proglacial lakes were formed in the second half of the 20th century. The central part of Popov glacier retreated in a lake environment by 190m between 1956 and 1965, while sections outside the lake retreated by only 70 m (Reference Kutuzov and ShahgedanovaKutuzov and Shahgedanova, 2009). The same authors further noted an 80 m difference in retreat between lake-calving and land-terminating sections of south Ashu-Tor glacier in 1977-2006.

The observed rate of calving flux is difficult to link with the mass loss of calving glaciers in the Tien Shan owing to the limited availability of data. Reference Mayer, Lambrecht, Hagg, Helm and ScharrerMayer and others (2008) found that the calving flux from southern Inylchek glacier, central Tien Shan, to glacier-dammed Merzbacher lake is 12 × 10⁶m³, which is substantially more than the estimated calving flux of Petrov glacier. This is probably related to the higher horizontal surface velocity of the larger Inylchek glacier (up to 0.8 m d^-1) compared with the much slower (up to 0.1 md^-1) ice front of Petrov glacier (Reference BakovBakov, 1983; Reference Dyurgerov, Liu and XieDyurgerov and others, 1995; Reference Mayer, Lambrecht, Hagg, Helm and ScharrerMayer and others, 2008). Although existing estimates of the calving flux are inherently variable (Reference Warren and AniyaWarren and Aniya, 1999; Reference Van der VeenVan der Veen, 2002), our estimates of total mass loss due to calving from Petrov glacier are consistent with computed wastage of lake-calving glaciers in the Swiss Alps and Alaska. Reference Funk and RothlisbergerFunk and Rothlisberger (1989) reported calving fluxes of six glaciers in the Alps in the range 5.7-11.6 × 10⁶m³. According to Reference Boyce, Motyka and TrufferBoyce and others (2007), the mean annual mass loss of Mendenhall Glacier, Coast Mountains, Alaska, was 0.75 × 10⁷m³w.e. over the period 2002-05. In contrast, the mean annual calving flux reported for glaciers in Patagonia was two orders of magnitude higher at the end of the 20th century (Reference Warren and AniyaWarren and Aniya, 1999; Reference Koppes, Conway, Rasmussen and ChernosKoppes and others, 2011). The higher calving flux of these glaciers, similar to tidewater glaciers, may be attributed to their size and terminus position in deep lakes (Reference VenterisVenteris, 1999).

The data on Petrov lake development document substantial increase of the surface area at the end of the 1970s, and continued enlargement during subsequent decades. From 1978 to 2006, the surface area increased at a mean annual rate of 2%, experiencing comparable enlargement to that of the proglacial lake Dzhuuku-5 in the Terskey-Alatoo. The surface area of this moraine-dammed lake at 3830ma.s.l. grew from 0.133 km² in 1976 (Reference Sevast’yanovSevast’yanov, 1980) to 0.254km² in 2005, i.e. at a mean annual rate of 1.6%. The regular annual increase (1.0-1.5%) of Petrov lake since 1999 is in contrast to the unstable development of glacial lakes in the Terskey-Alatoo. Area changes of glacier lakes presented by Reference Narama, Duishonakunov, Kääb, Daiyrov and AbdrakhmatovNarama and others (2010a) indicate strong interannual surface variations ranging from about -25% to 65%. The difference between surface changes of Petrov lake and that of lakes in the Terskey-Alatoo may be attributed to different volumes and phases of lake development. The surface area of the lakes noted by Reference Narama, Duishonakunov, Kääb, Daiyrov and AbdrakhmatovNarama and others (2010a) is one or two orders of magnitude smaller than the area of Petrov lake, and their sensitivity to annual changes in climate conditions is limited.

Summary and Conclusions

Petrov glacier had undergone relatively small annual oscillations until 1977, retreating at a mean rate of 13.6m a^-1. At the end of the 1970s, the recession increased substantially and an accelerated retreat followed during the next three decades. Three periods of rapid glacier retreat (1977-80, 1990-95 and 2008-09) separated by two intervals with a reduced retreat occurred during this time. The rate of recession was highest between 1990 and 1995, when a mean annual retreat of 54.1 ± 8.2 m occurred along the centre flowline of the glacier. Between 1977 and 2009, the glacier retreated at a mean rate of 40.6 ± 0.5 ma^-1 and the surface area of the glacier terminus decreased at a rate of 0.07 ± 0.0002 km² a^-1. The increased retreat rate of Petrov glacier during this period coincides with the accelerated glacier wastage reported from various regions in the Tien Shan. The observed length changes of the glacier correlate with warm-season air temperatures recorded in the second half of the 20th century. A negative correlation between warm-season temperatures and retreat rate indicates the importance of the ablation-season warming for glacier retreat. Between the periods 1930-97 and 1998-2009, the mean temperature and the length of the ablation season increased by 0.8°C and by 0.6 months, respectively. However, the correlation between climate conditions and glacier retreat was not found for the period 1999-2009.

The lake-terminating northern section of the Petrov glacier terminus has retreated at rates up to three times higher than the land-terminating southern part since 1995. The difference in retreat rate may be attributed to intense calving of the northern section in a proglacial environment and to subfreezing conditions in the marginal parts of the southern part of the glacier base. The calving accounted for a total mass loss of (17.7 ±0.4) × 10⁶m³ from the calving glacier terminus in 2003-09. The mean annual calving flux reached typical values reported for alpine lake-calving glaciers. The calving terminus of Petrov glacier is more than 65 m thick according to GPR measurements, whereas the maximum thickness of the southern flow unit is about 38 m. The surface lowering of the glacier terminus ranged from 2 m in debris-covered areas to 5 m in sections exposed to direct insolation from 2006 to 2009.

The fluctuations of Petrov glacier controlled the growth of the area and volume of Petrov lake. The mean annual changes of the lake surface area were low until the late 1970s, ranging up to 0.03 km² a^-1. After 1978, the surface enlargement increased and remained high for the next three decades. The mean annual enlargement by 0.04-0.1 km² is in agreement with the lake surface changes in the central Tien Shan over the past three decades. The mean annual volume increase ranged from (1.3 ± 0.012) × 10⁶m³ to (2.2 ± 0.009) × 10⁶m³ in 1978-2009. The volume enlargement increased considerably in 2003-06, when the calving glacier terminus retreated from the deepest part of the lake basin. The decrease in the depth of the lake basin and rising bed topography below the calving glacier terminus detected by the GPR survey imply that the increase in volume will not be as fast in the near future. The retreat rate of Petrov glacier is likely to be affected by the decrease in the depth of the lake, hence a slower retreat may be expected in the near future. This piece of knowledge was unknown to authors of hazard assessments from Petrov lake outbursts (e.g. Reference Janský, Šobr and EngelJanský and others, 2010) and it should stimulate further studies.