A multi‐dimensional analysis of pro‐glacial landscape change at Sólheimajökull,southern Iceland

Pro‐glacial landscapes are some of the most active on Earth. Previous studies of pro‐glacial landscape change have often been restricted to considering either sedimentological, geomorphological or topographic parameters in isolation and are often mono‐dimensional. This study utilized field surveys and digital elevation model (DEM) analyses to quantify planform, elevation and volumetric pro‐glacial landscape change at Sólheimajökull in southern Iceland for multiple time periods spanning from 1960 to 2010. As expected, the most intense geomorphological changes persistently occurred in the ice‐proximal area. During 1960 to 1996 the pro‐glacial river was relatively stable. However, after 2001 braiding intensity was higher, channel slope shallower and there was a shift from overall incision to aggradation. Attributing these pro‐glacial river channel changes to the 1999 jökulhlaup is ambiguous because it coincided with a switch from a period of glacier advance to that of glacier retreat. Furthermore, glacier retreat (of ~40 m yr^?1) coincided with ice‐marginal lake development and these two factors have both altered the pro‐glacial river channel head elevation. From 2001 to 2010 progressive increase in channel braiding and progressive downstream incision occurred; these together probably reflecting stream power due to increased glacier ablation and reduced sediment supply due to trapping of sediment by the developing ice‐marginal lake. Overall, this study highlights rapid spatiotemporal pro‐glacial landscape reactions to changes in glacial meltwater runoff regimes, glacier terminus position, sediment supply and episodic events such as jökuhlaups. Recognizing the interplay of these controlling factors on pro‐glacial landscapes will be important for understanding the geological record and for landscape stability assessments. Copyright © 2014 John Wiley & Sons, Ltd.     

Controls on englacial sediment deposition during the November 1996 jökulhlaup,Skeiđarárjökull,Iceland

This paper presents sedimentary evidence for rapid englacial debris entrainment during jökulhlaups. Previous studies of jökulhlaup sedimentology have focused predominantly on proglacial impact, rather than depositional processes within glaciers. However, observations of supraglacial floodwater outbursts suggest that englacial sediment emplacement is possible during jökulhlaups. The November 1996 jökulhlaup from Skei?arárjökull, Iceland presented one of the first opportunities to examine englacial flood deposits in relation to former supraglacial outlets. Using observations from Skei?arárjökull, this paper identifies and explains controls on the deposition of englacial flood sediments and presents a qualitative model for englacial jökulhlaup deposition. Englacial jökulhlaup deposits were contained within complex networks of upglacier‐dipping fractures. Simultaneous englacial deposition of fines and boulder‐sized sediment demonstrates that englacial fracture discharge had a high transport capacity. Fracture geometry was an important control on the architecture of englacial jökulhlaup deposits. The occurrence of pervasively frozen flood deposits within Skei?arárjökull is attributed to freeze‐on by glaciohydraulic supercooling. Floodwater, flowing subglacially or through upglacier‐dipping fractures, would have supercooled as it was raised to the surface faster than its pressure‐melting point could increase as glaciostatic pressure decreased. Evidence for floodwater contact with the glacier bed is supported by the ubiquitous occurrence of sheared diamict rip‐ups and intra‐clasts of basal ice within jökulhlaup fractures, deposited englacially some 200–350 m above the bed of Skei?arárjökull. Evidence for fluidal supercooled sediment accretion is apparent within stratified sands, deposited englacially at exceptionally high angles of rest in the absence of post‐depositional disturbance. Such primary sediment structures cannot be explained unless sediment is progressively accreted to opposing fracture walls. Ice retreat from areas of former supraglacial outbursts revealed distinct ridges characterized by localized upwellings of sediment‐rich floodwater. These deposits are an important addition to current models of englacial sedimentation and demonstrate the potential for post‐jökulhlaup landform development. Copyright © 2001 John Wiley & Sons, Ltd.     

The impact of pro‐glacial lakes on downstream sediment delivery in Norway

Glacier recessions caused by climate change may uncover pro‐glacial lakes that form important sedimentation basins regulating the downstream sediment delivery. The impact of modern pro‐glacial lakes on fluvial sediment transport from three different Norwegian glaciers: Nigardsbreen, Engabreen and Tunsbergdalsbreen, and their long‐term development has been studied. All of these lakes developed in modern times in overdeepened bedrock basins. The recession of Nigardsbreen uncovered a 1.8 km long and on average 15 m deep pro‐glacial lake basin during 1937 to 1968. Since then the glacier front has been situated entirely on land, and the sediment input and output of the lake has been measured. The suspended sediment transport into and out of the lake averaged 11 730 t yr^?1 and 2340 t yr^?1 respectively. Thus, 20% remained in suspension at the outlet. The measured mean annual bedload supplied to the lake was 11 800 t yr^?1, giving a total transport of 23 530 t yr^?1 which corresponds to a specific sediment yield of 561 t km^?2 yr^?1. A 1.9 km long and up to 90 m deep pro‐glacial lake basin downstream from Engabreen glacier was uncovered during 1890 to 1944. The average suspended sediment load delivered from the glacier during the years 1970–1981 amounted to 12 375 t yr^?1and the transport out of the lake was 2021 t yr^?1, giving an average of 16% remaining in suspension. The mean annual bedload was 8000 t yr^?1, thus the total transport was 20 375 t yr^?1, giving a specific sediment yield of 566 t km^?2 yr^?1. For Tunsbergdalsbreen glacier, measurements in the early 1970s indicated that the suspended sediment transport was on average 44 000 t yr^?1. From 1987 to 1993 the recession of the glacier uncovered a small pro‐glacial lake, 0.3 km long and around 9 m deep. Downstream from this, the suspended sediment load measured in 2009 was 28 000 t yr^?1, indicating that as much as 64% remained in suspension. Flow velocity, grain size of sediment, and morphology of the lake are important factors controlling the sedimentation rate in the pro‐glacial lakes. A survey of the sub‐glacial morphology of Tunsbergdalsbreen revealed that there are several overdeepened basins beneath the glacier. The largest is 4 km long and 100 m deep. When the glacier melts back they will become lakes and act as sedimentation basins. Despite an expected increase in sediment yield from the glacier, little sediment will pass these lakes and downstream sediment delivery will be reduced markedly. Beneath Nigardsbreen there was only a small depression that may form a lake and the sediment delivery will not be significantly affected. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Ice fracturing during jökulhlaups: implications for englacial floodwater routing and outlet development

Theoretical studies of glacial outburst floods (jökulhlaups) assume that: (i) intraglacial floodwater is transported efficiently in isolated conduits; (ii) intraglacial conduit enlargement operates proportionally to increasing discharge; (iii) floodwater exits glaciers through pre‐existing ice‐marginal outlets; and (iv) the morphology and positioning of outlets remains fixed during flooding. Direct field observations, together with historical jökulhlaup accounts, confirm that these theoretical assumptions are not always correct. This paper presents new evidence for spatial and temporal changes in intraglacial floodwater routing during jökulhlaups; secondly, it identifies and explains the mechanisms controlling the position and morphology of supraglacial jökulhlaup outlets; and finally, it presents a conceptual model of the controls on supraglacial outbursts. Field observations are presented from two Icelandic glaciers, Skeiðarárjökull and Sólheimajökull. Video footage and aerial photographs, taken before, during and after the Skeiðarárjökull jökulhlaup and immediately after the Sólheimajökull jökulhlaup, reveal changes in floodwater routing and the positioning and morphology of outlets. Field observations confirm that glaciers cannot transmit floodwater as efficiently as previously assumed. Rapid increases in jökulhlaup discharge generate basal hydraulic pressures in excess of ice overburden. Under these circumstances, floodwater can be forced through the surface of glaciers, leading to the development of a range of supraglacial outlets. The rate of increase in hydraulic pressure strongly influences the type of supraglacial outlet that can develop. Steady increases in basal hydraulic pressure can retro‐feed pre‐existing englacial drainage, whereas transient increases in pressure can generate hydraulic fracturing. The position and morphology of supraglacial outlets provide important controls on the spatial and temporal impact of flooding. The development of supraglacial jökulhlaup outlets provides a new mechanism for rapid englacial debris entrainment. Copyright © 2000 John Wiley & Sons, Ltd.     

Glacier retreat and landform production on an overdeepened glacier foreland: the debris‐charged glacial landsystem at Kvíárjökull,Iceland

Glacier recession and landform development in a debris‐charged glacial landsystem characterized by an overdeepening is quantified using digital photogrammetry, digital elevation model (DEM) construction and mapping of the Icelandic glacier Kvíárjökull for the period 1945–2003. Melting of ice‐cores is recorded by surface lowering rates of 0·8 m yr^–1 (1945–1964), 0·3 m yr^–1 (1964–1980), 0·015 m yr^–1 (1980–1998) and 0·044 m yr^–1 (1998–2003). The distribution/preservation of pushed and stacked ice‐cored moraine complexes are determined by the location of the long‐term glacial drainage network in combination with retreat from the overdeepening, into which glacifluvial sediment is being directed and where debris‐rich ice masses are being reworked and replaced by esker networks produced in englacial meltwater pathways that bypassed the overdeepening and connected to outwash fans prograding over the snout. Recent accelerated retreat of Kvíárjökull, potentially due to increased mass balance sensitivity, has made the snout highly unstable, especially now that the overdeepening is being uncovered and the snout flooded by an expanding pro‐glacial, and partially supraglacial, lake. This case study indicates that thick sequences of debris‐charged basal ice/controlled moraine have a very low preservation potential but ice‐cored moraine complexes can develop into hummocky moraine belts in de‐glaciated terrains because they are related to the process of incremental stagnation, which at Kvíárjökull has involved periodic switches from transport‐dominant to ablation‐dominant conditions. Glacier recession is therefore recorded temporally and spatially by two suites of landforms relating to two phases of landform production which are likely typical for glaciers occupying overdeepenings: an early phase of active, temperate recession recorded by push moraines and lateral moraines and unconfined pro‐glacial meltwater drainage; and a later phase of incremental stagnation and pitted outwash head development initiated by the increasing topographic constraints of the latero‐frontal moraine arc and the increasing importance of the overdeepening as a depo‐centre. Copyright © 2012 John Wiley & Sons, Ltd.     

Propagation of a subglacial flood wave during the initiation of a jökulhlaup

Abstract

Observations from the jökulhlaup from Grímsvötn in Vatnajökull, south-eastern Iceland, in 1996 indicate that the jökulhlaup was initiated by the movement of a localised pressure wave that travelled 50 km in 10 h from Grimsvötn to the terminus, forming a subglacial pathway along the glacier bed. Shortly after this wave reached the terminus, the jökulhlaup was flowing at a high discharge through a tunnel that would have needed much longer time to form by ice melting as assumed in existing theories of jökulhlaups. Frozen sediments formed in crevasses and frazil ice on the surface of the flood waters indicate the flow of supercooled water in the terminus region, demonstrating that the rate of heat transfer from subglacial flood water to the overlying ice is greatly underestimated in current theories.     

Sedimentation in arctic proglacial lakes: Mittivakkat Glacier,south‐east Greenland

Several sediment cores were collected from two proglacial lakes in the vicinity of Mittivakkat Glacier, south‐east Greenland, in order to determine sedimentation rates, estimate sediment yields and identify the dominant sources of the lacustrine sediment. The presence of varves in the ice‐dammed Icefall Lake enabled sedimentation rates to be estimated using a combination of X‐ray photography and down‐core variations in ¹³⁷Cs activity. Sedimentation rates for individual cores ranged between 0·52 and 1·06 g cm⁻² year⁻¹, and the average sedimentation rate was estimated to be 0·79 g cm⁻² year⁻¹. Despite considerable down‐core variability in annual sedimentation rates, there is no significant trend over the period 1970 to 1994. After correcting for autochthonous organic matter content and trap efficiency, the mean fine‐grained minerogenic sediment yield from the 3·8 km² basin contributing to the lake was estimated to be 327 t km⁻² year⁻¹. Cores were also collected from the topset beds of two small deltas in Icefall Lake. The deposition of coarse‐grained sediment on the delta surface was estimated to total in excess of 15 cm over the last c. 40 years. In the larger Lake Kuutuaq, which is located about 5 km from the glacier front and for which the glacier represents a smaller proportion of the contributing catchment, sedimentation rates determined for six cores collected from the centre of the lake, based on their ¹³⁷Cs depth profiles, were estimated to range between 0·05 and 0·11 g cm⁻² year⁻¹, and the average was 0·08 g cm⁻² year⁻¹. The longer‐term (c. 100–150 years) average sedimentation rate for one of the cores, estimated from its unsupported ²¹⁰Pb profile, was 0·10–0·13 g cm⁻² year⁻¹, suggesting that sedimentation rates in this lake have been essentially constant over the last c. 100–150 years. The average fine‐grained sediment yield from the 32·4 km² catchment contributing to the lake was estimated to be 13 t km⁻² year⁻¹. The ¹³⁷Cs depth profiles for cores collected from the topset beds of the delta of Lake Kuutuaq indicate that in excess of 27 cm of coarse‐grained sediment had accumulated on the delta surface over the last approximately 40 years. Caesium‐137 concentrations associated with the most recently deposited (uppermost) fine‐grained sediment in both Icefall Lake and Lake Kuutuaq were similar to those measured in fine‐grained sediment collected from steep slopes in the immediate proglacial zone, suggesting that this material, rather than contemporary glacial debris, is the most likely source of the sediment deposited in the lakes. This finding is confirmed by the ¹³⁷Cs concentrations associated with suspended sediment collected from the Mittivakkat stream, which are very similar to those for proglacial material. Copyright © 2000 John Wiley & Sons, Ltd.     

Constraints on the dynamics of subglacial basalt eruptions from geological and geochemical observations at Kverkfjöll, NE-Iceland

The Kverkfjöll area, NE Iceland is characterised by subglacial basalt pillow lavas erupted under thick ice during the last major glaciation in Iceland. The water contents of slightly vesiculated glassy rims of pillows in six localities range from 0.85±0.03 to 1.04±0.03 wt %. The water content measurements allow the ice thickness to be estimated at between 1.2 and 1.6 km, with the range reflecting the uncertainty in the CO₂ and water contents of the melt. The upper estimates agree with other observations and models that the ice thickness in the centre of Iceland was 1.5–2.0 km at the time of the last glacial maximum. Many of the pillows in the Kverkfjöll area are characterised by vesiculated cores (40–60% vesicles) surrounded by a thick outer zone of moderately vesicular basalt (15–20% vesicles). The core contains ～1 mm diameter spherical vesicles distributed uniformly. This observation suggests a sudden decompression and vesiculation of the still molten core followed by rapid cooling. The cores are attributed to a jökulhlaup in which melt water created by the eruption is suddenly released reducing the environmental pressure. Mass balance and solubility relationships for water allow a pressure decrease to be calculated from the observed change of vesicularity of between 4.4 and 4.7 MPa depressurization equivalent to a drop in the water level in the range 440–470 m. Consideration of the thickness of solid crust around the molten cores at the time of the jökulhlaup indicates an interval of 1–3 days between pillow emplacement and the jökulhlaup. Upper limits for ice melting rates of order 10^?3 m/s are indicated. This interpretation suggests that jökulhlaups can reactivate eruptions.     

Outburst susceptibility assessment of moraine‐dammed lakes in Western Himalaya using an analytic hierarchy process

Glacier retreat results in the formation and expansion, and sometimes outburst, of moraine‐dammed lakes worldwide. Sudden outburst floods from such lakes have caused enormous damage to settlements and infrastructure located downstream. Such lakes located in the Himalayan region are highly prone to outburst floods due to climatic conditions and geotectonic settings. In this study, multi‐temporal Landsat images from 2002–2014, digital elevation models (DEMs), geomorphic analysis and modelling were used to assess the changes in glacial lakes and the outburst susceptibility of moraine‐dammed lakes in the Chandra–Bhaga basin of the north‐western Indian Himalaya. An inventory of lakes was developed using satellite data, thematic maps and ground‐based investigations for the Chandra–Bhaga basin. The total area of all glacial lakes (size >5000 m²) increased by 47% from 2002 to 2014, with a pronounced increase of 57% for moraine‐dammed lakes. Sixteen moraine‐dammed lakes were identified and assessed for outburst susceptibility using the analytic hierarchy process (AHP). Forty‐one reported glacial lake outburst flood (GLOF) events from moraine‐dammed lakes in Himalayan regions were analysed, culminating in the identification of 11 critical factors for assessing outburst susceptibility using the AHP, including those related to the lake area and change, surrounding terrain characteristics, dam geometry, regional seismicity and rainfall history. The past three GLOF events in the Himalayan region were used to validate the method and to classify moraine‐dammed lakes as having very high, high, medium or low outburst susceptibility. Eight lakes classified as very high and high outburst susceptibility should be further investigated in detail. The proposed AHP‐based approach is suitable for first‐order identification of critical lakes for prioritising future detailed investigation and monitoring of moraine‐dammed glacial lakes in the Himalayan region. Copyright © 2017 John Wiley & Sons, Ltd.     

Proglacial sediment trapping in recently formed Silt Lake,upper Lillooet Valley,Coast Mountains,British Columbia

The sedimentology of proglacial Silt Lake was assessed by lake sediment coring and monitoring of lacustrine processes during a late‐summer period of high glacier melt to characterize sediment delivery from the heavily glacierized catchment and investigate the sediment trapping dynamics of this upland lake. A complete varve chronology was established for a distal basin of the lake which was exposed by Lillooet Glacier retreat between 1947 and 1962. The varve record showed decreasing sedimentation rates in the basin while the glacier retreated, and as the lake became free of ice contact in the early 1970s. Although recession has continued over recent decades, and glacier proximity to the lake has, therefore, continued decreasing, lacustrine sedimentation rates are now accelerating due to changing basin morphometry caused by delta progradation. Over shorter time scales, lake sedimentation patterns respond to changing runoff conditions, including late‐summer glacier melt intensity, intra‐annual flooding events, diumal runoff fluctuations, and within‐lake turbidity currents. Turbidity currents included quasi‐regular flows during high diurnal discharges and an episodic flushing of temporarily stored sediment from the sandur or delta at a time of low stage. Suspended sediment yield to Silt Lake is estimated to exceed 10³ Mg km^?2 a^?1, a magnitude that surpasses previous local and regional yield estimates for the glacierized headwaters of the Lillooet River valley. Since Silt Lake currently traps a significant prooportion of that upland sediment supply, and the trapping efficiency of the basin has been variable at decadal time scales, the formation and continued development of Lilt Lake has likely had a significant influence on downstream sediment delivery. Lacustrine sediment‐based proxies of long‐term hydroclimatic variability being developed in glacially distal settings should include provisions for dynamic sediment trapping effects in upstream water bodies that often form in the active proglacial environment. Copyright © 2008 John Wiley & Sons, Ltd.     

Reconstructing historic Glacial Lake Outburst Floods through numerical modelling and geomorphological assessment: Extreme events in the Himalaya

Recession of high‐mountain glaciers in response to climatic change frequently results in the development of moraine‐dammed glacial lakes. Moraine dam failure is often accompanied by the release of large volumes of water and sediment, termed a Glacial Lake Outburst Flood (GLOF). Chukhung Glacier is a small (~3 km²) receding valley glacier in Mt. Everest (Sagarmatha) National Park, Nepal. Unlike many Himalayan glaciers, which possess a thick mantle of supraglacial debris, its surface is relatively clean. The glacier terminus has receded 1.3 km from its maximum Holocene position, and in doing so provided the space for an ice‐contact moraine‐dammed lake to develop. The lake had a maximum volume of 5.5 × 10⁵ m³ and drained as a result of breaching of the terminal moraine. An estimated 1.3 × 10⁵ m³ of material was removed from the terminal moraine during breach development. Numerical dam‐breach modelling, implemented within a Generalised Likelihood Uncertainty Estimation (GLUE) framework, was used to investigate a range of moraine‐dam failure scenarios. Reconstructed outflow peak discharges, including failure via overtopping and piping mechanisms, are in the range 146–2200 m³ s^‐1. Results from two‐dimensional hydrodynamic GLOF modelling indicate that maximum local flow depths may have exceeded 9 m, with maximum flow velocities exceeding 20 m s^‐1 within 700 m of the breach. The floodwaters mobilised a significant amount of material, sourced mostly from the expanding breach, forming a 300 m long and 100 m wide debris fan originating at the breach exit. moraine‐dam. These results also suggest that inundation of the entire floodplain may have been achieved within ten minutes of initial breach development, suggesting that debris fan development was rapid. We discuss the key glaciological and geomorphological factors that have determined the evolution of a hazardous moraine‐dammed lake complex and the subsequent generation of a GLOF and its geomorphological impact. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Formation,mechanism and significance of alluvial-dammed lakes in Golmud River catchment,north-eastern Qinghai-Tibetan Plateau

This study investigated a series of dammed lakes and downstream-adjacent alluvial fans in the upstream to middle reaches of the Golmud River in the eastern Kunlun Mountain, on the north-eastern Qinghai-Tibetan Plateau (QTP). An optically stimulated luminescence (OSL) chronology shows the sediments of five dammed lakes developed from c. 45–40, 30–25, 18–14, and 12–8 ka, corresponding to MIS 3b, late MIS 3a, Last Deglaciation, and early Holocene, respectively. The remote sensing data show these dammed lakes have a total area of 109.4 km², with the lake volume of more than 4.0 km³. Symmetric alluvial fans from north–south tributary valleys produced OSL ages of c. 61–52, 42–31, 26–20, and 16–10 ka, corresponding to glaciation periods: the MIS 3c and MIS 3a, MIS 2, and the Last Deglaciation. This suggests that glacial activity is responsible for the alluvial fan development, where dammed rivers occurred first, but lake formation did not take place synchronously until later periods of strong hydrologic activity, resulting from northward intrusions of the Indian summer monsoon (ISM) or glacier melt. Thus, the blocking pattern is that river valleys were dammed during periods of glacial activity and lakes formed during wet periods. The lake formation and subsequent drainage may have resulted in: (i) impeded headwater incision and strengthening of downstream dissection; (ii) enriched the halite and potash in the distal Qarhan Salt Lake through hydrologic and hydrochemical processes of abundant water input, the salt lake expansion, salt redissolution from playa and final resedimentation during later dry periods. The alluvial-dammed lake pattern in the mountain-basin systems of eastern Kunlun Mountain offers a model for assessing the linkages between monsoon dynamics, geomorphic processes and distal salt lake evolutions in other arid regions.     

OSL dating of past lake levels for a large dammed lake in southern Tibet and determination of possible controls on lake evolution

Fluvio‐lacustrine terraces along Phung Chu (river) on the central southern Tibetan Plateau indicate that a large palaeo‐dammed‐lake formerly existed in this area. Based on landscape survey, optically stimulated luminescence (OSL) dating and sedimentary analyses, this research shows that the Phung Chu was blocked and a dammed‐lake over 2500 km² in size formed before 30 ka ago. OSL dating analysis suggests the fluvio‐lacustrine sediments were well bleached and yield accurate age estimates for two lake drainage events. The first drainage event took place after 30 ka, resulted in river incision and formed a high terrace at 50 m height from the present river level. The second drainage happened after 3.7 ka, resulted in further river incision and formed the second terrace at 25 m height from the present river level. According to the distribution of the fluvio‐lacustrine sediments, active normal faults (particularly the Kharta Fault) in this region and the high gradient slopes after Phung Chu enters the Yö Ri gorge, seismically‐induced landsliding is regarded as highly likely to have been the cause of river blockage and associated formation of a dammed‐lake, although glacial damming is also a possible cause. The volume of drainages from this dammed‐lake may have led to catastrophic flooding and analogous modern lakes represent significant geo‐hazard risks to down‐river human settlements. As dammed‐lakes are special phases in fluvial evolution, often involving river blockage, breakthrough and drastic catchment change, these processes can reveal how tectonic or climatic events modify landforms. However, such tectonic‐derived landform changes can also impact palaeo‐climate of the region. Thus this study has added new evidence regarding the evolutionary history of a dammed lake including its formation, duration, extent and final drainage, which is crucial for understanding its general landscape process mechanisms and for better assessing geo‐hazard risks. Copyright © 2016 John Wiley & Sons, Ltd.     

Rapid expansion of glacial lakes caused by climate and glacier retreat in the Central Himalayas

Glacial lake outburst floods are among the most serious natural hazards in the Himalayas. Such floods are of high scientific and political importance because they exert trans‐boundary impacts on bordering countries. The preparation of an updated inventory of glacial lakes and the analysis of their evolution are an important first step in assessment of hazards from glacial lake outbursts. Here, we report the spatiotemporal developments of the glacial lakes in the Poiqu River basin, a trans‐boundary basin in the Central Himalayas, from 1976 to 2010 based on multi‐temporal Landsat images. Studied glacial lakes are classified as glacier‐fed lakes and non‐glacier‐fed lakes according to their hydrologic connection to glacial watersheds. A total of 119 glacial lakes larger than 0.01 km² with an overall surface area of 20.22 km² (±10.8%) were mapped in 2010, with glacier‐fed lakes being predominant in both number (69, 58.0%) and area (16.22 km², 80.2%). We found that lakes connected to glacial watersheds (glacier‐fed lakes) significantly expanded (122.1%) from 1976 to 2010, whereas lakes not connected to glacial watersheds (non‐glacier‐fed lakes) remained stable (+2.8%) during the same period. This contrast can be attributed to the impact of glaciers. Retreating glaciers not only supply meltwater to lakes but also leave space for them to expand. Compared with other regions of the Hindu Kush Himalayas (HKH), the lake area per glacier area in the Poiqu River basin was the highest. This observation might be attributed to the different climate regimes and glacier status along the HKH. The results presented in this study confirm the significant role of glacier retreat on the evolution of glacial lakes. Copyright © 2014 John Wiley & Sons, Ltd.     

Late‐summer thermal regime of a small proglacial lake

This study was motivated by an interest in understanding the potential effects of climate change and glacier retreat on late summer water temperatures in alpine areas. Fieldwork was carried out between July and September 2007 at Place Lake, located below Place Glacier in the southern Coast Mountains of British Columbia. Place Lake has an area of 72 000 m², a single inlet and outlet channel, and an approximate residence time of 4 days. Warming between the inlet and outlet of the lake ranged up to 3 °C and averaged 1.8 °C, which exceeds the amount of warming that occurred over the 1 km reach of Place Creek between the lake outlet and tree line. Over a 23‐day period, net radiation totalled about 210 MJ·m^–2, with sensible heat flux adding another 56 MJ m^‐2. The latent heat flux consumed about 8% of the surface heat input. The dominant heat sink was the net horizontal advection associated with lake inflow and outflow. Early in the study period, temperatures between the surface and 6‐m depth were dominantly at or above 4 °C and were generally neutral to thermally stable, whereas temperatures decreased with depth below 6 m and exhibited irregular sub‐diurnal variations. The maximum outflow temperature of almost 7 °C occurred in this period. We hypothesize that turbidity currents associated with cold, sediment‐laden glacier discharge formed an underflow and influenced temperatures in the deeper portion of the lake but did not mix with the upper layers. Later in the study period, the lake was dominantly well mixed with some near‐surface stability associated with nocturnal cooling. Further research is required to examine the combined effects of sediment concentrations and thermal processes on mixing in small proglacial lakes to make projections of the consequences of glacier retreat on alpine lake and stream temperatures. Copyright © 2011 John Wiley & Sons, Ltd.     

Hazardous processes and events from glacier and permafrost areas: lessons from the Chilean and Argentinean Andes

Glacier and permafrost hazards such as glacial‐lake outburst floods and rock–ice avalanches cause significant socio‐economic damages worldwide, and these processes may increase in frequency and magnitude if the atmospheric temperature rises. In the extratropical Andes nearly 200 human deaths were linked to these processes during the twentieth century. We analysed bibliographical sources and satellite images to document the glacier and permafrost dynamics that have caused socio‐economic damages in this region in historic time (including glacial lake outburst floods, ice and rock–ice avalanches and lahars) to unravel their causes and geomorphological impacts. In the extratropical Andes, at least 15 ice‐dammed lakes and 16 moraine‐dammed lakes have failed since the eighteenth century, causing dozens of floods. Some floods rank amongst the largest events ever recorded (5000 × 10⁶ m³ and 229 × 10⁶ m³, respectively). Outburst flood frequency has increased in the last three decades, partially as a consequence of long‐term (decades to centuries) climatic changes, glaciers shrinkage, and lake growth. Short‐term (days to weeks) meteorological conditions (i.e. intense and/or prolonged rainfall and high temperature that increased meltwater production) have also triggered outburst floods and mass movements. Enormous mass failures of glaciers and permafrost (> 10 × 10⁶ m³) have impacted lakes, glaciers, and snow‐covered valleys, initiating chain reactions that have ultimately resulted in lake tsunamis and far‐reaching (> 50 km) flows. The eruption of ice‐covered volcanoes has also caused dozens of damaging lahars with volumes up to 45 × 10⁶ m³. Despite the importance of these events, basic information about their occurrence (e.g. date, causes, and geomorphological impact), which is well established in other mountain ranges, is absent in the extratropical Andes. A better knowledge of the processes involved can help to forecast and mitigate these events. Copyright © 2014 John Wiley & Sons, Ltd.     

Late Pleistocene outburst floods from Issyk Kul,Kyrgyzstan?

Elevated shorelines and lake sediments surrounding Issyk Kul, the world's second largest mountain lake, record fluctuating lake levels during Quaternary times. Together with bathymetric and geochemical data, these markers document alternating phases of lake closure and external drainage. The uppermost level of lake sediments requires a former damming of the lake's western outlet through the Boam gorge. We test previous hypothesised ice or landslide dam failures by exploring possible links between late Quaternary lake levels and outbursts. We review and recompile the chronology of reported changes in lake site, and offer new ages of abandoned shorelines using ¹⁴C in bivalve and gastropod shells, and plant detritus, as well as sand lenses in delta and river sediments using Infrared Stimulated Luminescence. Our dates are consistent with elevated lake levels between ~45 ka and 22 ka. Cosmogenic ¹⁰Be and ²⁶Al exposure ages of fan terraces containing erratic boulders (>3 m) downstream of the gorge constrain the timing of floods to 20.5–18.5 ka, postdating a highstand of Issyk Kul. A flow‐competence analysis gives a peak discharge of >10⁴ m³ s^–1 for entraining and transporting these boulders. Palaeoflood modelling, however, shows that naturally dammed lakes unconnected to Issyk Kul could have produced such high discharges upon sudden emptying. Hence, although our data are consistent with hypotheses of catastrophic outburst floods, average lake‐level changes of up to 90 mm yr^–1 in the past 150 years were highly variable without any outbursts, so that linking lake‐level drops to catastrophic dam breaks remains ambiguous using sedimentary archives alone. This constraint may readily apply to other Quaternary lakes of that size elsewhere. Nonetheless, our reconstructed Pleistocene floods are among the largest reported worldwide, and motivate further research into the palaeoflood hydrology of Central Asia. Copyright © 2017 John Wiley & Sons, Ltd.     

A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

Retreating ice: research in pro‐glacial areas matters

Glacier forefields are landscapes in transition from glacial to non‐glacial conditions; this implies intense geomorphic, hydrological and ecological dynamics with important on‐ and off‐site effects. This special issue collects 13 papers covering recent research in both (sub‐)polar and alpine pro‐glacial environments that focus on (i) pro‐glacial sediment sources, (ii) pro‐glacial rivers, (iii) pro‐glacial lakes, (iv) ground water and ice, and (v) the development of soil and vegetation in its interplay with morphodynamics. Advances in mapping, surveying and geophysical techniques form the basis for research perspectives related to the historical evolution of pro‐glacial areas, the understanding of complex interactions of multiple processes, and the effects of continued glacier recession. Copyright © 2015 John Wiley & Sons, Ltd.     

Flocculation: a key process in the sediment flux of a large,glacier‐fed lake

This study investigates the consequences of flocculation for sediment flux in glacier‐fed Lillooet Lake, British Columbia based on density, fractal dimension, in situ profiles of sediment concentration and size distribution, and settling velocity equations presented in the literature. Sediment flux attributed to macroflocs during the late spring and summer accounts for a significant portion of sediment flux in the lake, equivalent to at least one‐quarter of the average annual sediment flux. Fine sediment is reaching the lake floor faster in flocs than occurs if settling as individual grains. This flux varies both spatially and temporally over the observation period, suggesting a link between deposition via flocculation and the properties of bottom sediments. Macrofloc flux increased through June, reached a peak during July, and then declined into August. Macrofloc flux was greatest in the distal end of the first basin, approximately 10 km from the point of inflow. Relatively high excess densities (～0·1 g cm^–3 at 500 µm) for flocs in situ are consistent with a composition dominated by inorganic primary particles. Microlaminations within Lillooet Lake varves have been linked by earlier workers to discharge events, and the action of turbidity currents, emanating from the Lillooet River. While turbidity currents undoubtedly occur in Lillooet Lake, these results demonstrate flocculation as an adjunct process linking discharge, lake level, macrofloc flux, bulk density and microlaminations. In situ measurements of sediment settling velocity in glacier‐fed lakes are required to better constrain flux rates, and permit comparison between flocculation in lacustrine environments with existing studies of estuarine, marine and fluvial flocculation. Copyright © 2009 John Wiley & Sons, Ltd.