The formation and development of glacial lakes in mountainous
regions is one of the consequences of glacier recession. Such lakes may drain
partially or completely when the stability of their dams is disturbed or as
a consequence of impacts. We present a case study from the Central Asian
mountain range of Tien Shan – a north-oriented tributary of the Adygine
Valley, where the retreat of a polythermal glacier surrounded by permafrost
has resulted in the formation of several generations of lakes. The aim of
this study was to analyse the past development of different types of glacial
lakes influenced by the same glacier, to project the site's future
development, and to evaluate the outburst susceptibility of individual lakes
with an outlook for expected future change. We addressed the problem using a
combination of methods, namely bathymetric, geodetic and geophysical on-site
surveys, satellite images and digital elevation model analysis, and modelling
of glacier development. Based on this case of the glacial lakes being of
varied age and type, we demonstrated the significance of glacier ice in lake
development. Lake 3, which is in contact with the glacier terminus, has
changed rapidly over the last decade, expanding both in area and depth and
increasing its volume by more than 13 times (7800 to 106 000 m
Glacial lakes often appear as gems gleaming amid harsh mountain
environments. At the same time, they can pose a serious threat to downstream
settlements and infrastructure. If we focus on Asian mountain ranges, the
now common term “glacial lake outburst flood” (GLOF), first used for the
Himalayan region, is ubiquitous. Mainly due to the development of
satellite-based sensor technologies, the regions under the scrutiny of
investigators have increased from the Himalayas (Benn et al., 2012; Fujita
et al., 2009; Shrestha et al., 2010) and Karakoram (Chen et al., 2010;
Haemmig et al., 2014) to the Tibetan Plateau (Liu et al., 2014; Wang et al.,
2013; Zhang et al., 2017), Pamir (Mergili and Schneider, 2011), and Tien Shan
(Bolch et al., 2011; Engel et al., 2012; Narama et al., 2017; Sorg et al.,
2012). In the territory of Kyrgyzstan, there are about 2000 glacial lakes
(
Some glacial lakes form and drain within a relatively short time (Erokhin et al., 2017). Others exist for years and decades without major change. Lastly, there are lakes that are in a phase of enlargement. The expansion of glacial lakes broadly correlates with climatic warming trends and negative glacier mass balance patterns (Lei et al., 2012; Mergili et al., 2013; Zhao et al., 2015). In addition, a recent study has shown that the expansion of glacial lakes under a dry continental climate regime can be closely related to thermokarst processes in permafrost as well (Li et al., 2014). This means that lakes enlarge when in contact with a retreating glacier tongue by the melting of ice in the lake basin bed or around the sides or by filling due to increased inflow from a glacier. Other factors that play a role in lake evolution include the erosion of surface drainage channels, the formation, expansion, or blockage of subsurface drainage channels, dam morphology changes, or slope movements adjacent to the lake.
Outbursts from mountain lakes occur almost annually in Kyrgyzstan (in recent years, these included Lake Merzbacher, 2017, Lake Chelektor, 2017, and Lake Aksai, 2015), and given the generally increasing use of mountain valleys, greater losses can be expected as a consequence (Dussaillant et al., 2010). Once an outburst is triggered, a flood (often evolving into a hyperconcentrated flow or debris flow) will threaten areas in the lower parts of a valley. To reduce the risk, an assessment framework (see, e.g. the GAPHAZ Technical Guidance Document; GAPHAZ, 2017) can be applied, which includes the following steps: (i) identification of potentially dangerous lakes, (ii) detailed evaluation of the hazard, and (iii) application of mitigating measures. Our study includes the second step – a detailed evaluation of the outburst susceptibility of selected lakes that had already been identified as potentially dangerous (Erokhin and Zaginaev, 2016). As the need for field data for such an assessment has been generally acknowledged, here, we present an assessment based on a combination of on-the-spot and remotely obtained data.
Although the study site is not particularly suited for the purpose of only monitoring glacier retreat, we did find it most appropriate in terms of observing the site in a more holistic way. The glacier itself is considered to be representative of smaller glaciers, which are subject to relatively rapid retreat rates (Narama et al., 2010) and generally have shorter response times to changes in climatic conditions compared to larger ones (Wang et al., 2014). This allowed us to follow a larger development period at a smaller timescale. Also, the existence of glacial lakes documenting glacier shrinkage was convenient – being at various stages of development, they were seen as able to help us to better understand the evolution of different types of lakes over time. Due to the occurrence of lake outbursts in the region, the high elevation of the site, and its proximity to populated areas, an outburst susceptibility assessment strategy, with an estimation of future change, was considered to be desirable and of practical relevance.
Based on a detailed study of one particular site, we looked at the issue of glacier development and the associated, actively evolving glacier forefield and its consequences on downstream areas. We have used the term “ice-debris complex” (as coined by Bolch et al., 2019) to describe the glacier and the debris landforms with varying shares of buried ice in front of it. In order to address the problems associated with the rapid processes that operate in mountainous regions, we incorporated perspectives of future scenarios for the development of the site, and the related outburst susceptibility, into our study. Because our intention was to promote a more complex approach to the issue of susceptibility assessment, our study included three particular objectives: (i) to analyse the past evolution of different types of glacial lakes influenced by the same glacier, (ii) to provide insight into the probable future evolution of the site, and (iii) to evaluate the outburst susceptibility of individual lakes, with an outlook for possible future change. As a consequence of our objectives, the structure of the paper is as follows: first, we present an analysis of the past development of the local setting, following both the glacier and the proglacial lakes; then, we report on the potential future conditions for the site; and, finally, we include evaluations of the lakes' outburst susceptibilities and estimations of how they will likely evolve under changing climatic conditions. In our opinion, to analyse glacier retreat and lake formation without touching on the topic of consequences (present and near future) would be as incomplete as evaluating outburst susceptibility without knowing the past evolution of the site.
The research site is situated in northern Tien Shan, specifically Kyrgyz
Ala-Too. The highest peak of this east–west-oriented range is Pik Semenova
Tienshanskogo, at 4895 m a.s.l. Although the Kyrgyz Ala-Too is not as
severely glaciated as the ranges in central Tien Shan, the currently
glacier-free parts of the valleys do show signs of the former glacier
extent. The main Ala Archa valley currently has
Representation of different lake types in Kyrgyz Ala-Too, northern
Tien Shan, in 2017. Only lakes with a minimum area of 1500 m
The development of glaciation in the Ala Archa Basin is monitored, in the
long term, mainly due to its position in proximity to the Kyrgyz capital,
Bishkek, and the popularity of the valley with tourists visiting the
National Park. According to Aizen et al. (2006), there was a change in the
glaciated area of the Ala Archa watershed by
The Adygine watershed has an area of 39.6 km
Overview of the Adygine ice-debris complex
The upper part of the complex involves a polythermal glacier (2.8 km
Basic parameters of the three largest lakes of the site.
A Holocene moraine complex, at an elevation of 3450 m a.s.l., forms the lower part of the Adygine complex. It consists of several parts of varying age (Shatravin, 2000), some of them being identifiable based on visual inspection only. The western part is composed of creeping ice-rich debris, with a typical relief of arcuate ridges and furrows in its terminal part; ice lenses are present and occasionally exposed (Fig. 1b). From the east, a smaller ice-debris matrix moves under the influence of gravity and adjoins the central landform. The main body is a so-called debris-derived rock glacier that formed below the glacier terminus. Lying probably within the discontinuous permafrost zone, this landform has a prominently oversteepened front that is almost 700 m wide. In its western part, thermokarst processes are manifested in subsidence craters, cracks, and sagging, and numerous thermokarst lakes and exposed ice lenses (Fig. 1c) can be found here. The eastern part is considered to be of an older generation, with a rather flattened surface and stable lakes with lower turbidity.
Climatic conditions in the area are continental, characterised by relatively
low precipitation and high annual and daily air temperature fluctuations.
This is clearly evident at the Ala Archa station (2200 m a.s.l.), which
reports a 22
A research station, built at the site in 2008 next to the Lake 2 surface outflow, provided an adequate location for field measurements to be taken. Climatic data and ground temperatures were obtained from an automated weather station situated near the Adygine glacier terminus at 3560 m a.s.l. The hydrological regime of the lakes was monitored with pressure sensors installed in Lake 1 (2012–2015), Lake 2 (2007–2017), and Lake 3 (2012–2017).
The glacier retreat has been monitored since 2006, during the ablation
season, by means of a geodetic survey carried out with a Leica TCR 705 total
station (infrared beam on a reflective prism, accuracy: 0.005 m). The survey
was focused solely on the terminal part of the tongue, where the retreat is
most significant and can result in the exposure of overdeepenings (potential
future lakes). This time series was supplemented with glacier limits
acquired from aerial images taken in 1988 and 1962 (scale of the survey:
The selected key locations at the site (the lakes' surroundings and their dams) were surveyed by electrical resistivity tomography and spontaneous polarisation, in order to detect buried glacier remnants, frozen ground, and possible seepage routes. In order to acquire information on glacier thickness and bed topography for adjusting the model of future glacier extent and identifying potential areas of lake formation, ground-penetrating radar (GPR) was used. The location of the profiles is shown in Fig. 1a.
The method of electrical resistivity tomography (ERT) consists of measuring the resistivity of the subsurface by means of a number of electrodes located along a profile. Interpretation of such measured data was performed using the software Res2Dinv (Loke and Barker, 1995), which provided 2-D inversion of the measured resistance data and calculated the 2-D cross section under the profile. The resistive section of the observed profile approximates the actual distribution of the individual resistive layers in the depth range given by the maximum spacing of the electrodes. We used the ARS-200E model (GF Instruments, CR), with 48 electrodes at 4.4 m increments. The maximum range of the electrodes was 206.8 m, with a depth range of about 50 m.
The spontaneous polarisation (SP) method measures the natural electrical potential of the rock environment; in the case of lake dams, the filtering potential that arises from water filtration through the porous environment can be tested. Measurements were made with a GEOTOR I with offset compensation. Non-polarisable electrodes were used to minimise the effect of transient resistance between the measuring electrode and the geological environment. The measurements were carried out using potential variation, in which the electrical potential of a given point on the profile is measured relative to a fixed reference point. The measurement steps were around 10 m, and the location of every fifth point was recorded using GPS.
GPR was used to acquire information on the glacier thickness and subglacial topography. GPR data were collected along the central longitudinal profile, for five transverse and three short, connecting profiles. An unshielded 50 MHz rough terrain antenna and MALÅ CU-II control unit were used for data collection. The signal acquisition time was set to 860 ns and the scan spacing was about 0.07 m. The collected data were processed and interpreted using Reflex-Win software version 4.5.
The glacier mass balance model – glacier evolution runoff model, or GERM – applied in this study has been well established for simulating glacier mass balance from climate data, based on a sophisticated degree-day approach (Hock, 1999; Huss et al., 2008). In our case, the future glacier mass balances of the Adygine glacier until 2050 were calculated with the GERM, forced by empirically downscaled scenarios' data for daily air temperature and precipitation, as described below. As the glacier mass balance was not measured at the Adygine glacier, we used the mass balance data from the nearby Golubina glacier (WGMS, 2017) for validating the simulated mass balances of the GERM against measured mass balances for the period 1981–1994. Mass balances were validated for both winter and summer. In addition, we were able to use discharge data from the Adygine River (from gauging near the confluence with the Ala Archa River) for 1960–1987 to validate the modelled discharge of the GERM with observed data.
Despite the generally poor availability of data in Central Asia, daily air temperature means and precipitation totals were gathered from two stations close to the study area: Alplager (2130 m a.s.l., 1978–2010) and Baityk (1580 m a.s.l., 1914–1979). The meteorological station data were used to calibrate and validate the downscaling model (multivariate regression and analogue method; Benestad, 2004; Benestad et al., 2007). Reanalysis data from NCEP/NCAR (Kalnay et al., 1996) were used as predictor fields in order to downscale long-term datasets of precipitation and air temperature for the study region. For scenario simulations until 2050, three different SRES scenarios – A1B, B1, and A2 (IPCC, 2007) – were performed by the ECHAM5/MPI-OM coupled atmosphere–ocean model (Röckner et al., 2003). Although the SRES scenario data were replaced by the Representative Concentration Pathways approach of the Intergovernmental Panel on Climate Change (IPCC, 2013), the uncertainty ranges of future large-scale temperature and precipitation changes provided by the two approaches are rather similar and smaller compared to the uncertainty range of climate data coming from further downscaling. We argue that, for the purpose of assessing future glacier extent, in particular for the future development of glacial lakes, these scenarios are sufficiently realistic.
In order to estimate the potential locations for glacial lake formation, the data on glacier bed topography acquired by GPR were used (see Supplement). We also consulted the digital elevation model (DEM) data for the glacier surface from the Shuttle Radar Topography Mission, which has a resolution of 30 m (Farr et al., 2007). The combined information on future glacier extent and the topography of the exposed area was analysed in GIS by determining flow directions and identifying sinks (i.e. overdeepenings).
The presented lake outburst susceptibility assessment is qualitative and was adapted to accommodate the distinctive regional features – dry climate, the common occurrence of buried ice and permafrost. We based the assessment procedure on knowledge we gained through fieldwork in the region, focused especially on the hydrological systems of proglacial lakes. The specificity of the region is the origin of the lakes – the most common type is lakes formed in intramorainic depressions (Table 1). Many lake hazard evaluations are based on cases of moraine dams with a distinct morphometry; however, this type of lake does not occur that frequently in the region, representing only 2 % of the potentially dangerous ones (Erokhin and Zaginaev, 2016). Nevertheless, inspiration was drawn from assessments by ICIMOD (ICIMOD, 2011), Allen et al. (2016), Frey et al. (2010), and Huggel et al. (2004). Hazard assessments based solely on remotely sensed data are appropriate when applied over a large area and as a first step; however, at least for lakes labelled as potentially dangerous, additional evaluations, which would include on-site investigation, need to be carried out.
In this paper, the overall susceptibility is introduced as being a
combination of conditioning factors (i.e. a lake's inherent characteristics
determining the lake's susceptibility to bursting and causing flooding) and
the presence of possible triggers that would have the capacity to cause an
outburst. To assess a lake's inherent susceptibility, the following
parameters were selected and a simple qualitative rating suggested (Table 3):
The
Inherent susceptibility parameters and ratings of lakes. The scheme on the right shows possible combinations of parameter ratings (H – high, M – medium, L – low) and the resulting susceptibility. Note that this simplifying table serves as a guideline for the susceptibility estimate and that even medium resulting susceptibility means a lake is dangerous and may burst and cause damage.
The determination of the degree to which a lake is susceptible to outburst is followed by an assessment of potential triggers. The means of determining whether a trigger has the potential to cause outburst in the specific lake or not may be distant (satellite imagery, DEM), as is the case for the first group of triggers in Table 4, or can be on-site, as is necessary for inner dam stability or subsurface channel functioning, for example.
Triggers with potential to cause lake outburst and their description.
The fall of mass into a lake is one of the most common triggers of lake outburst (Ding and Liu, 1992; Emmer and Cochachin, 2013; Falatkova, 2016), resulting in an impact wave that can overtop or destabilise the dam. Following a widely used procedure, the distance between a lake and a slope and slope steepness are the parameters that determine the trigger potential (Alean, 1985; Fischer et al., 2012; Noetzli et al., 2003). In the case of calving, a lake's contact with a cliff-like glacier terminus is a crucial condition.
The development and stability of a glacial lake can be affected by the melting of buried ice within a dam. During stagnation of a debris-covered glacier, an ice-cored moraine dam can be formed. These cores present a weak point in a dam, as their ablation changes the dam's inner structure (Richardson and Reynolds, 2000). Besides that, dam stability is also a function of the overall frozen or unfrozen condition of the dam material.
Not only landslide dams are prone to failure caused by internal erosion (piping). Cases of GLOF due to piping in moraine dams have also been observed (Xu, 1988). Moraines consist of unconsolidated heterogeneous material and therefore may be prone to seepage and piping as a conduit grows and a dam is weakened (Awal et al., 2011). Another dam-stability-decreasing phenomenon – earthquakes – has been confirmed as being a primary trigger of several lake outbursts (Clague and Evans, 2000; Lliboutry et al., 1977).
In the case of a lake not having surface drainage, and meltwater from a glacier being routed into the lake, there is the potential for a significant hydrostatic pressure increase at times of higher melt rates of snow or ice. A sudden air temperature rise could therefore cause a disturbance in the lake's hydrological balance, with possible consequences for dam stability. An increase in hydrostatic pressure can also be caused by the blockage of a subsurface channel, which is a very unpredictable, and still rather poorly understood, phenomenon. Monitoring changes in lake morphometry due to frequent sliding of lake walls or disruptions in the hydrological regime of a lake can serve as a tool for determining channel blockage potential.
Changes in lake outburst susceptibility result from a complex process associated with the interaction of geomorphological processes and the behaviour of the glacio-hydrological system. Development of the englacial and subsurface routing of meltwater is especially crucial; while being difficult to determine, this was beyond the scope of this paper. The presented estimate of susceptibility change is linked to climatic changes (rising MAAT) only and involves impacts accompanying deglaciation that have been observed in other mountain ranges, such as glacier retreat, altered glacier melt rate, slope instability (debuttressing effect), permafrost degradation, and buried ice melting (Allen et al., 2016; Frey et al., 2010; Haeberli et al., 2017). The general assumptions for future climatic conditions, glacier retreat, and altered melt rates are supported by the glacier mass balance model (Sect. 3.3).
The terminus position of the investigated glacier has changed significantly
since the 1960s. Over the past five decades, the furthermost part of the
glacier terminus has retreated
Retreat of the Adygine glacier tongue between 1962 and 2017. 1–7 are the proglacial lakes. The arrow marks the rock outcrop (riegel). Base image: Worldview-2 (2011).
There are numerous lakes of varying size, type, and stage of development in the study area. They form a three-level cascade, and at least some of the lakes are hydrologically connected.
The lowest part of the cascade is formed from about a dozen thermokarst
lakes, situated on a large morainic, perennially frozen platform at an
elevation of
Lake 2, the largest lake at the study site, began forming in about 1960,
after the western glacier tongue retreated behind a rock riegel at 3540 m a.s.l. Since then, the lake area has extended, together with further
terminus recession, reaching 32 700 m
Electrical resistivity results along Profile 1 and Profile 2 for Lake 2 and Profile 3 for Lake 3. Resistivity values of materials depend on grain size, pore size, water content, salinity, temperature, and phase resistivity increases due to the change from electrically conductive water to non-conductive ice (Kneisel et al., 2008). For location of the profiles (P1–P3), see Fig. 1a.
Since 2004, several lakes have formed in the proximity of the retreating glacier terminus, filling rather shallow depressions. The lakes' development is closely linked to the amount and routing of glacier meltwater as well as to basin stability. Our data suggest the importance of direct contact with the terminus, as only those lakes showed an increase in area. Lake 4 is an example of one that was enlarging slightly until it lost contact with the glacier in 2011–2012. Most lakes formed, enlarged, or remained stable over a period of several years and finally drained by subsurface routes (for detailed information on the lakes' morphometric changes, see the Supplement). An exception is Lake 6, which seemed stable over the monitoring period, until it suddenly drained in 2015 but was filled again the following year. The area exposed following the glacier retreat is underlain by buried ice lenses (visual inspection; Fig. 1a) within frozen ground (documented by the measured resistivity; Fig. 3). A possible explanation for the lake's sudden drainage might be permafrost degradation caused by warming by the lake water (Shur and Jorgensen, 2007).
Although these proglacial lakes are highly unstable, the hazard they pose is
negligible, as only a small amount of water (up to 4000 m
Bathymetric map and area changes of Lake 3 between 2007 and 2017. Lake depth data obtained during a field survey at the site carried out in 2017.
At the level of Lake 1, there are several small glacial lakes of varying age (Fig. 5), which were formed by thermokarst processes. As these lakes only retain small volumes of water, and no significant areal development has been observed, a closer investigation from the outburst susceptibility point of view was not carried out. However, a few lakes situated in the youngest generation of a moraine recently uncovered the basin sides formed by buried glacier ice (Fig. 5, left). These lakes are expected to develop further in the future, either by enlarging or draining through newly opened subsurface channels.
Thermokarst lakes of varying age in the moraine complex below
the Adygine glacier. These small lakes are situated at
According to all three modelled scenarios, an increase in air temperature in
the area is highly likely in future decades, with an estimated higher
temperature increase in spring and summer and a moderate temperature
increase in autumn and winter. Scenarios of annual precipitation showed a
reduction under all three scenarios. In particular, summer and autumn
depicted a decrease in precipitation. Winter and spring displayed a minor
reduction, or even a slight increase, in precipitation totals (for winter
under A2, for spring under B1). Consequently, future mass balances simulated
for the Adygine glacier by the GERM are negative, and runoff will be
altered. Low-pass-filtered values of the mass balance of the Adygine glacier
were negative throughout the simulation period (Fig. 6a). Scenario A1B
showed significantly less negative mass balances between 2015 and 2035
compared to the other scenarios. In this period, Scenario A2 was
characterised by the highest decadal variation, with filtered values of
between
Glacier mass balance
Due to the negative mass balances (Fig. 6a), all scenarios resulted in a
reduction in glacier area. Under Scenario A1B, the glacier area reduced
only slightly in the first two decades (2015–2035), when compared to the
other two scenarios (Fig. 6b). After 2035, Scenario A1B showed accelerated
areal shrinkage (resulting from the accelerated mass loss). All three
scenarios simulated a significant decrease in glacier area by 2050,
specifically leaving only 73.2 % (A1B), 62.3 % (A2), and 55.6 % (B1) of
the 2012 glacier area. The model simulations indicate (Fig. 7) that, under
Scenario B1, the glacier will disintegrate into several parts, whereas,
under Scenario A2, only one smaller part will remain. The modelled glacier
retreat implies the potential for new lakes to form in the exposed area. The
glacier bed topography (Fig. A2), based on GPR profiling, revealed several
overdeepenings. Under Scenario A1B, three new lakes, with a total area of
Modelled future glacier extent and potential areas for lake formation for 2050 under Scenario A1B, Scenario A2, and Scenario B1.
As a result of greater negative glacier mass balances, total runoff from the Adygine glacier catchment is expected to rise. According to Scenario A1B and Scenario A2, the total runoff will increase by 19.7 % and 25 %, respectively, between 2010 and 2050 (Fig. 8a). However, Scenario B1 showed a rather sensitive reaction of the glacier dynamics to climate change, with a clear peak in discharge values around 2020. This was followed by a decline and overall stagnation in the following decade, resulting in no significant trend in the modelled period.
Obviously, runoff for individual months will change in the future. The most pronounced changes were seen at the beginning and end of an ablation period – specifically, in April, May, and October. Slight increases in runoff are also expected in June and September. However, for the core period of the ablation season – July and August – no significant changes were modelled. Separating total runoff into its two main components – melt from snow and melt from ice – for both today and the future (Fig. 8b) highlighted a significant temporal shift for the snow component. The time of peak runoff from snowmelt, which was in mid-June for the first decade of 2010–2020, will change to mid-May in the decade 2040–2050. This will coincide with a slightly earlier start of ice melt at the glacier, as the snow cover will have depleted earlier as well, thus exposing bare ice.
The resulting lakes' inherent susceptibility to outburst.
The inherent outburst susceptibility of the studied lakes, based on the lakes' characteristics, determined here by five parameters, is summarised in Table 5. Lake 1 and Lake 3 were marked as having medium susceptibility, whereas the largest lake of the site, Lake 2, was considered to be less inclined to burst and cause a flood. Several areas were determined as being a potential source of mass movement that could reach the lakes (Fig. 9). These are mainly the steep upper part of the glacier with visible deformations and cracks and also the eastern lateral moraine with buried ice. Both surface and subsurface hydrological connections among the lakes indicate routing of excess water in the case of rapid melting or lake drainage. As shown by the model projections, and also indicated by the position of the ELA, the glacier is expected to shrink significantly in the coming decades. Sorg et al. (2014) confirmed that future susceptibility is likely to change due to glacier retreat and the accompanying consequences. For instance, the slopes may become more prone to destabilisation due to permafrost degradation, in particular those that are south-oriented (i.e. exposed to direct solar radiation).
Geomorphological conditions and potential triggers of lake outburst at the Adygine ice-debris complex. View from the north (Google Earth, 2018).
Lake 1 has a medium susceptibility to cause flooding, mainly due to its origin (intramorainic depression), subsurface drainage, and the potential for enlargement. Outburst could be triggered by a fall of mass from the north-oriented slope of unconsolidated debris, which contains buried ice. The weak point is the subsurface drainage channel that leads through the moraine complex. Changes in the capacity and, above all, the blockage of such a channel would lead to the lake filling up and subsequent outburst. That could also be initiated by the outburst or overflow of Lake 3. Its outburst susceptibility is expected to increase in the future, due to the intensified melting of buried and exposed ice near the lake basin and destabilised slopes in its vicinity.
Lake 2 has a low susceptibility to cause flooding because of its stable surface drainage, distance from the glacier terminus and the associated ending of lake basin development. As there are no proximal unstable slopes (only rock outcrops), outburst is not likely to be triggered by a displacement wave from the fall of mass into the lake. A possible weak point is represented by buried ice and water seepage in the western part of the dam that could lead to changes in the subsurface drainage network. A change in outburst susceptibility is expected to be negligible in the future, as surface drainage would be able to deal with possible greater inflow, and the lake stability will not be threatened by unstable slopes.
Lake 3 has a medium susceptibility to cause flooding, due to it being in contact with the glacier, having relatively unstable subsurface drainage, and having the potential to grow both in area and depth. Despite the contact with the glacier tongue, a dam overflow by impact wave from calving is not probable because the terminus is not crevassed and a calving front has not developed. However, the lake is within reach of an ice avalanche from the steep glacier slope and also a landslide from the adjacent, unstable lateral moraine. The eastern part of the dam is prone to failure, due to the presence of buried ice. Overflow of the dam caused by increased inflow from the glacier may lead to progressive dam erosion. Outburst susceptibility is expected to increase in combination with increasing lake volume and hydrostatic pressure on the dam. Also, destabilised slopes in the lake's proximity will threaten its stability.
Individually, the lakes present a certain amount of threat. Nevertheless, as these lakes are interconnected, the overall hazard represented by the site should be considered. In this case, a chain reaction is a very probable scenario. The lowest lake (Lake 1) is subject to changes in its basin morphometry and thus in the functioning of its subsurface drainage channel. In the case of a rapid increase in drainage from Lake 2, or outburst from Lake 3, this lake would be hit, although it is rather unclear how this still-evolving, subsurface-drained basin would react. In regard to similar cases in the region, it can be expected that the lake will fill up, with subsequent rapid drainage.
The GPR data show a two-layered structure of the glacier, with prevailing cold ice and a temperate layer along the glacier bed. The terminus is frozen to the ground, and cold ice dominates the glacier base 450 m up-valley. At mid-length, on a narrow glacier snout around 3700 m a.s.l., the ice thickness increases to more than 50 m, and temperate ice appears at the glacier base. Temperate ice spreads along the base over the next 1 km and terminates below cold ice in the accumulation area. The cold–temperate transition along most of this section is located 30 to 50 m below the glacier's surface. The upper part of the glacier (above 3900 m a.s.l.) consists mainly of ice below the pressure-melting point and is thus connected to the permafrost environment in bedrock.
The presence of both types of ice confirms the polythermal regime of the Adygine glacier that was originally assumed only for large valley glaciers in the Tien Shan region (Dyurgerov et al., 1995). The distribution of temperate ice is similar to that of the characteristic polythermal glaciers described by Copland and Sharp (2001) and Etzelmüller and Hagen (2005). Following the classification scheme by thermal structure proposed by Irvine-Fynn et al. (2011), the Adygine glacier corresponds to a thermal configuration with the presence of temperate ice near the bed in the ablation area. The smaller extent of warm ice in the accumulation zone implies a different temperature structure than that of the Tuyuksu glacier (the Zailiyskiy Alatau), which has a maximum thickness of warm ice in the accumulation zone, where the cold–temperate transition reaches the glacier surface (Nosenko et al., 2016).
Lake 2 is a typical example of a lake formed when a glacier retreats over a rock step. The lake's growth was coupled with glacier recession until the terminus reached a higher elevation than the lake outflow. After the loss of contact, the lake stabilised and its formation ended. Lake 3 had a similar genesis and is currently in the stage of growth. Thanks to intensive thermal exchange between lake water and ice, a glacier terminus recedes by melting, and sometimes even by calving, in favour of lake area and depth. As a lake loses contact with a glacier tongue, some development can continue due to the melting of residual ice in the lake bed or dam. The extent of changes is, however, considerably reduced.
The case of Lake 1, formed in an intramorainic depression, demonstrates the
difference in development. The lake's development was not terminated or
stabilised after loss of contact with the glacier. Its further evolution was
linked to buried ice, its exposure and melting, and the degradation of
permafrost. The depression formed in the 1960s and filled with water two
decades later; yet even today, the lake is subject to changes in basin
morphology and drainage. Such lakes do not grow to large dimensions like
typical moraine-dammed lakes (e.g. Tsho Rolpa and Dig Tsho); however, they
can hold sufficient volumes of water to cause an outburst flood. In many
cases, such lakes are even non-stationary, filling up only when a subsurface
channel gets blocked, often in times of increased inflow due to snow or ice
melting (Erokhin et al., 2017). The 2012 outburst from the Teztor Lake,
situated in a neighbouring valley at a similar altitude to Adygine, serves
as an example. When filled, the lake has a volume of 70 000 m
Huss et al. (2014) described, in detail, the uncertainties of glacier and runoff modelling in general. In our case, the GERM was calibrated using the summer and winter mass balance data of Golubina glacier from 1981 to 1994. Whereas simulated winter mass balances for Golubina showed no consistent trend of under- or overestimation, the modelled summer mass balance slightly overestimated the real conditions from 1981 to 1984 and generally underestimated it from 1985 to 1994. The performance of the model for simulating discharge is described by the Nash–Sutcliffe efficiency (Nash and Sutcliffe, 1970), which is 0.64 based on monthly mean values of the Adygine catchment for 1960–1987. The biases and uncertainties in the downscaled data for precipitation over the complex mountain topography remain a weakness and need further improvements in the future. Also, the lack of a DEM from 1960, as well as the glacier area from this time, attenuated the accuracy of the calibration of the GERM simulations.
Both the area and volume of glaciers all over Tien Shan are expected to
decrease throughout the 21st century (Sorg et al., 2012). Aizen et al. (2007b) also predicted significant glacier degradation linked to the ELA
shifting to higher altitudes because rising air temperatures will not be
balanced out by sufficient increases in precipitation totals. Sorg et al. (2014) used the GERM for glacier and runoff modelling in Chon Kemin Valley
(Zailiyskiy and Kungey Alatau, Tien Shan), where the glaciers are expected
to have vanished by 2080 under the more pessimistic scenarios (dry–warm,
wet–warm). The future runoff results of Chon Kemin (Sorg et al., 2014) are
very similar to those of our study area, with “warm” scenarios expecting
peak runoff in the 2020s and the main change in spring runoff (increase)
being caused by higher winter precipitation and enhanced snowmelt in
spring. A general prediction for runoff change in Tien Shan is a possible
increase in the near future, followed by a steady decline until the end of
the 21st century (Sorg et al., 2012). However, different regions in
Tien Shan may exhibit varying responses and uncertainties connected to
future runoff modelling. An advanced state of glacier degradation was
presented by Huss and Fischer (2016), focusing on small-sized glaciers (
Our approach to the identification of potential locations for future lake formation is based on glacier bed topography derived from in situ data. For obvious reasons, it is not possible to apply such a method to the whole region. Numerical models are a good substitute and are widely used in large areas, from the Peruvian Andes (Colonia et al., 2017) to the Swiss Alps (Linsbauer et al., 2012) and the Himalayas–Karakoram (Linsbauer et al., 2016). The commonly used model is GlabTop, introduced by Linsbauer et al. (2009) and Paul and Linsbauer (2012), with an automated version (GlabTop2) introduced by Frey et al. (2014). The model calculates glacier ice thickness, based on a DEM, and glacier outlines and branch lines; the automated version (GlabTop2) avoids manual delineation of the lines. A similar concept and results (Frey et al., 2014) were obtained from a model presented in Huss and Farinotti (2012). Potential sites with overdeepenings can also be identified, or confirmed, based on the three morphological criteria (glacier surface features) introduced by Frey et al. (2010). This is a more laborious approach, which includes manual analysis of the DEMs and high-resolution satellite/aerial imagery. If we apply these morphological criteria to the Adygine glacier, one location with a distinct narrowing, and also a change in slope, can be identified. This coincides with our results pointing to the presence of overdeepenings at the location of the current glacier tongue. Nevertheless, as Haeberli et al. (2016b) correctly pointed out, there are still significant limitations to understanding the principles of the depth erosion by glaciers and thus to modelling glacier bed overdeepenings. Moreover, the estimation of debris cover volume left after glacier recession is a challenge, as excessive debris can smooth topographic irregularities and lead to the formation of outwash plains, and not lakes (Linsbauer et al., 2016). We are also aware of the limitations to our approach. Due to the low resolution of the GPR-derived ice thickness raster, the resulting extent of potential sites for lake formation must be interpreted with caution.
The presented evaluation strategy combined knowledge from the field with the
benefits of remotely obtained data, providing a good insight into the
interaction of the site's elements (a moraine, buried ice, perennially frozen
ground, meltwater channels, a glacier, lakes) and a broader spatial and
temporal perspective, thanks to the airborne imagery. Field investigations
are necessary to produce a realistic assessment of lake outburst
susceptibility. McKillop and Clague (2007) mentioned that knowledge of
seepage and lake bathymetry and its changes is important information for
understanding dam hydraulic conditions. Bolch et al. (2008) highlighted
subsurface glacier meltwater routing, which we also believe to be extremely
important, together with the buried ice distribution. At the Adygine
ice-debris complex, the ERT technique helped us to discover an extensive ice
core near Lake 2, buried under a considerably thick (8–10 m) layer of
mostly frozen debris. A similar structure was observed, for example, at the
moraine dam of the Thulagi glacier lake, Himalayas (Pant and Reynolds, 2000).
The seepage routes and hydrological processes within moraine dams were also
examined with SP and ERT in Thompson et al. (2012) at Miage glacier, Italian
Alps. The resistivity values attributed to materials detected at our study
site are consistent with the common resistivity value ranges of the different
materials summarised in Kneisel (2006). Frozen
ground can exhibit a wide range of resistivity values, varying from 10
Remote sensing has proved to be a very useful tool in the field of environmental studies and is the basis for many hazard assessments (e.g. Allen et al., 2016; Bolch et al., 2008; Huggel et al., 2004; McKillop and Clague, 2007). Comprehensive hazard assessments elaborated as a basis for wide usage include, for example, the GAPHAZ Technical Guidance Document (GAPHAZ, 2017), the strategy prepared by ICIMOD, focused on the high-mountain region of the Himalayas–Karakoram (ICIMOD, 2011), and a useful summary of various lake outburst hazard assessments is also provided in Emmer and Vilimek (2013).
For a first, large-scale outburst hazard assessment, the trigger is often
represented solely by a fall of mass into a lake (e.g. Allen et al., 2016).
The presented procedure aims to encompass all trigger factors that could
have the capacity to lower a dam's or lake basin's stability and cause lake
drainage. However, there are still some limitations to a full understanding
of the varying mechanisms of lake outbursts. We did not include permafrost
degradation as a separate trigger, but we did consider it to be one of the
main factors that can increase the outburst susceptibility in future, as
presented for the Alps by Haeberli et al. (2017). The presence of permafrost
at the site is supported by other studies (e.g. Marchenko et al., 2007), the
global permafrost model of Gruber (2012), and the site's negative MAAT
of recent years. Although the measurement period is rather short
(2008–2013), the long-term MAAT of the site may, in fact, be even lower,
given the low value (MAAT of
Glacier retreat in our site's case may not lead to new ice- or snow-avalanche-starting zones, but the rising MAAT could influence the glacier thermal regime and cause a shift in certain parts from cold- to warm-based, thereby lowering the friction between the glacier bed and steep ice masses (Huggel et al., 2008). Areas exposed after glacier recession often consist of unconsolidated material, which, if steep, may become a potential starting zone for landslides or debris flows (Haeberli et al., 2017). Already exposed steep slopes may become destabilised and fail due to the melting of interstitial ice or, where present, the degradation of mountain permafrost (Huggel et al., 2010). The formation of new lakes at higher elevations could increase the hazard of the site due to a possible cascade effect, where even a small volume release could trigger the outburst of a lower-lying lake. Such situations are expected to occur in the future, and the potentially arising challenges should be addressed in advance (Haeberli et al., 2016a).
This paper aimed to describe the impact of glacier retreat on its forefield, to highlight the differences in evolution and outburst susceptibility of proglacial lakes, and to promote a more field-based, regionally specified approach to hazard and susceptibility assessment.
The main results concerning the studied site are as follows:
The glacier is subject to relatively rapid recession, which is comparable
to that of other glaciers in Tien Shan. Since 1962, it has retreated, on
average, by 11.5 m a As a consequence, glacier runoff will change in the coming decades.
Specifically, the beginning of the ablation season is likely to occur earlier,
compared to the current situation, and the peak runoff from melting snow
will very likely shift from mid-June to mid-May, resulting in a slightly
earlier start to glacier ice melting. A three-level cascade of glacial lakes has formed in the glacier
forefield, with the largest (Lake 2) retaining a volume of about 200 000 m Lake 1 and Lake 3 were categorised as having medium outburst susceptibility,
mainly due to their subsurface drainage and the presence of buried ice in
their vicinities. Despite its large size, the outburst susceptibility of
Lake 2 was assessed as low. The destabilisation of steep slopes, exposure
and melting of buried ice, and changes in glacier runoff may further
deteriorate the stability of Lake 1 and Lake 3, leading to an increase in
outburst susceptibility in the future.
A very similar development in other glacier complexes in this region is expected, with runoff peaks occurring at present or in the near future, accompanied by permafrost degradation, the gradual melting of buried ice, and the formation of new lakes. This will result in new outburst hazards that could arise in a relatively short time. One of the main outputs of the study is an analysis of lake development dynamics, which indicate a lake's potential to grow and its related outburst susceptibility. As we have documented, lake size is not as important a factor, in terms of susceptibility, and we would rather emphasise the influence of buried ice in relation to the presence of perennially frozen ground and the unpredictability of subsurface drainage channels. In terms of future susceptibility change, we would like to draw attention to the influence of glacier runoff. The expected extension of the ablation season, and the related activity in glacial lakes, may lead to increased outburst susceptibility. The prolonged period of the melting season in relation to buried ice will result in a greater potential for the alteration of the hydrological network of drainage channels. Further research should be dedicated to the hydrological connections in proglacial morainic environments and the interactions of subsurface water flow with permafrost and buried ice, as there is still an insufficient understanding of these processes.
Data used in this study are available at
The supplement related to this article is available online at:
KF, MŠ, and BJ carried out the fieldwork and processed the data obtained. KF designed the outburst susceptibility assessment, produced the figures, and wrote the manuscript, with contributions from AN and WS for the modelling part. ZE provided results from the GPR survey; VB provided results from the ERT and SP surveys. AN and WS ran the glacier evolution model and HH supervised the EURAS-CLIMPACT project and initiated the concept of this joint paper. All authors contributed to improving drafts of the paper.
The authors declare that they have no conflict of interest.
Great thanks are due to Sergey Erokhin for providing support essential for successful field work, valuable material, and consultation of the hazard situation in Kyrgyz Ala-Too. We would like to thank Vitalii Zaginaev for reliable help with field work activities, Stefan Reisenhofer for performing the downscaling of climate data, and Matthias Huss for providing the GERM glacier model. The work connected to future development of the Adygine glacier was a part of the EURAS-CLIMPACT project. We highly appreciate the work of the reviewers and editors; special thanks go to Wilfried Haeberli, who provided us with many constructive comments and suggestions which helped to improve the manuscript.
This paper was edited by Richard Gloaguen and reviewed by Wilfried Haeberli, Praveen Thakur, and one anonymous referee.