Distribution of glaciofluvial sediment within and on the surface of a high arctic valley glacier: Marthabreen,Svalbard

Marthabreen is a 7·8 km long valley glacier in SW Spitsbergen. The glacier is partially covered by a layer of angular debris derived from rockfall in its accumulation area, pierced in places by pinnacles and ridges of glaciofluvial sediment. These concentrations of glaciofluvial sediment fall into three categories: (1) debris pinnacles; (2) longitudinal sediment dykes; (3) longitudinal ridge accumulations. Debris pinnacles are slabs of sediment (predominantly sands, gravels and cobbles) elevated to the glacier surface along thrusts. Longitudinal sediment dykes are low (<0·5 m high) ridges of debris melting out of vertical sediment dykes within the body of the glacier. They are composed of a range of facies including sands, granule gravels, pebble gravels and diamiction. These dykes are sub-parallel to the longitudinal foliation on the glacier and form during folding of the stratification. Longitudinal ridge accumulations are higher (>1 m high) ridges of sorted sand and gravels which are not associated with penetrative ice structures. Their occurrence downglacier of sediment dykes and debris pinnacles suggests that they originate as supraglacial or englacial channels or tunnels filled by sediment derived from the dykes or thrusts. The presence of large quantities of glaciofluvial sediment on the surface of Marthabreen does not imply englacial water flow at high levels within the glacier, but is related to ice deformational processes such as thrusting and folding of debris-rich stratification. Copyright © 1999 John Wiley & Sons, Ltd.     

Englacial and subglacial water flow at Skálafellsjökull,Iceland derived from ground penetrating radar,in situ Glacsweb probe and borehole water level measurements

We reconstruct englacial and subglacial drainage at Skálafellsjökull, Iceland, using ground penetrating radar (GPR) common offset surveys, borehole studies and Glacsweb probe data. We find that englacial water is not stored within the glacier (water content ~0–0.3%). Instead, the glacier is mostly impermeable and meltwater is able to pass quickly through the main body of the glacier via crevasses and moulins. Once at the glacier bed, water is stored within a thin (1 m) layer of debris‐rich basal ice (2% water content) and the till. The hydraulic potential mapped across the survey area indicates that when water pressures are high (most of the year), water flows parallel to the margin, and emerges 3 km down glacier at an outlet tongue. GPR data indicates that these flow pathways may have formed a series of braided channels. We show that this glacier has a very low water‐storage capacity, but an efficient englacial drainage network for transferring water to the glacier bed and, therefore, it has the potential to respond rapidly to changes in melt‐water inputs. © 2015 The Authors. Earth Surface Processes and Landforms published by 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.     

Heterogeneous water storage and thermal regime of supraglacial ponds on debris‐covered glaciers

The water storage and energy transfer roles of supraglacial ponds are poorly constrained, yet they are thought to be important components of debris‐covered glacier ablation budgets. We used an unmanned surface vessel (USV) to collect sonar depth measurements for 24 ponds to derive the first empirical relationship between their area and volume applicable to the size distribution of ponds commonly encountered on debris‐covered glaciers. Additionally, we instrumented nine ponds with thermistors and three with pressure transducers, characterizing their thermal regime and capturing three pond drainage events. The deepest and most irregularly‐shaped ponds were those associated with ice cliffs, which were connected to the surface or englacial hydrology network (maximum depth = 45.6 m), whereas hydrologically‐isolated ponds without ice cliffs were both more circular and shallower (maximum depth = 9.9 m). The englacial drainage of three ponds had the potential to melt ~100 ± 20 × 10³ kg to ~470 ± 90 × 10³ kg of glacier ice owing to the large volumes of stored water. Our observations of seasonal pond growth and drainage with their associated calculations of stored thermal energy have implications for glacier ice flow, the progressive enlargement and sudden collapse of englacial conduits, and the location of glacier ablation hot‐spots where ponds and ice cliffs interact. Additionally, the evolutionary trajectory of these ponds controls large proglacial lake formation in deglaciating environments. Copyright © 2017 John Wiley & Sons, Ltd.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 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.     

Properties of natural supraglacial debris in relation to modelling sub‐debris ice ablation

As debris‐covered glaciers become a more prominent feature of a shrinking mountain cryosphere, there is increasing need to successfully model the surface energy and mass balance of debris‐covered glaciers, yet measurements of the processes operating in natural supraglacial debris covers are sparse. We report measurements of vertical temperature profiles in debris on the Ngozumpa glacier in Nepal, that show: (i) conductive processes dominate during the ablation season in matrix‐supported diamict; (ii) ventilation may be possible in coarse surface layers; (iii) phase changes associated with seasonal change have a marked effect on the effective thermal diffusivity of the debris. Effective thermal conductivity determined from vertical temperature profiles in the debris is generally ~30% higher in summer than in winter, but values depend on the volume and phase of water in the debris. Surface albedo can vary widely over small spatial scales, as does the debris thickness. Measurements indicate that debris thickness is best represented as a probability density function with the peak debris thickness increasing down‐glacier. The findings from Ngozumpa glacier indicate that the probability distribution of debris thickness changes from positively skewed in the upper glacier towards a more normal distribution nearer the terminus. Although many of these effects remain to be quantified, our observations highlight aspects of spatial and temporal variability in supraglacial debris that may require consideration in annual or multi‐annual distributed modelling of debris‐covered glacier surface energy and mass balance. Copyright © 2012 John Wiley & Sons, Ltd.     

The influence of ground ice distribution on geomorphic dynamics since the Little Ice Age in proglacial areas of two cirque glacier systems

Holocene glaciers have contributed to an abundance of unstable sediments in mountainous environments. In permafrost environments, these sediments can contain ground ice and are subject to rapid geomorphic activity and evolution under condition of a warming climate. To understand the influence of ground ice distribution on this activity since the Little Ice Age (LIA), we have investigated the Pierre Ronde and Rognes proglacial areas, two cirque glacier systems located in the periglacial belt of the Mont Blanc massif. For the first time, electrical resistivity tomography, temperature data loggers and differential global positioning systems (dGPS) are combined with historical documents and glaciological data analysis to produce a complete study of evolution in time and space of these small landsystems since the LIA. This approach allows to explain spatial heterogeneity of current internal structure and dynamics. The studied sites are a complex assemblage of debris‐covered glacier, ice‐rich frozen debris and unfrozen debris. Ground ice distribution is related to former glacier thermal regime, isolating effect of debris cover, water supply to specific zones, and topography. In relation with this internal structure, present dynamics are dominated by rapid ice melt in the debris‐covered upper slopes, slow creep processes in marginal glacigenic rock glaciers, and weak, superficial reworking in deglaciated moraines. Since the LIA, geomorphic activity is mainly spatially restricted within the proglacial areas. Sediment exportation has occurred in a limited part of the former Rognes Glacier and through water pocket outburst flood and debris flows in Pierre Ronde. Both sites contributed little sediment supply to the downslope geomorphic system, rather by episodic events than by constant supply. In that way, during Holocene and even in a paraglacial context as the recent deglaciation, proglacial areas of cirque glaciers act mostly as sediment sinks, when active geomorphic processes are unable to evacuate sediment downslope, especially because of the slope angle weakness. Copyright © 2014 John Wiley & Sons, Ltd.     

Morphology and evolution of supraglacial hummocks on debris-covered Himalayan glaciers

Thick supraglacial debris layers often have an undulating, hummocky topography that influences the lateral transport of debris and meltwater and provides basins for supraglacial ponds. The role of ablation and other processes associated with supraglacial debris in giving rise to this hummocky topography is poorly understood. Characterizing hummocky topography is a first step towards understanding the feedbacks driving the evolution of debris-covered glacier surfaces and their potential impacts on mass balance, hydrology and glacier dynamics. Here we undertake a geomorphological assessment of the hummocky topography on five debris-covered glaciers in the Everest region of the central Himalaya. We characterize supraglacial hummocks through statistical analyses of their vertical relief and horizontal geometry. Our results establish supraglacial hummocks as a distinct landform. We find that a typical hummock has an elongation ratio of 1.1:1 in the direction of ice flow, length of 214 ± 109 m and width of 192 ± 88 m. Hummocky topography has a greater amplitude across-glacier (15.4 ± 10.9 m) compared to along the glacier flow line (12.6 ± 8.3 m). Consequently, hummock slopes are steeper in the across-glacier direction (8.7 ± 4.3°) than in the direction of ice flow (5.6 ± 4.0°). Longer, wider and higher-amplitude hummocks are found on larger glaciers. We postulate that directional anisotropy in the hummock topography arises because, while the pattern of differential ablation driving topography evolution is moderated by processes including the gravitational redistribution of debris across the glacier surface, it also inherits an orientation preference from the distribution of englacial debris in the underlying ice. Our morphometric data inform future efforts to model these interactions, which should account for additional factors such as the genesis of supraglacial ponds and ice cliffs and their impact on differential ablation.     

Interrelations among pyroclastic surge,pyroclastic flow,and lahars in Smith Creek valley during first minutes of 18 May 1980 eruption of Mount St. Helens,USA

A devastating pyroclastic surge and resultant lahars at Mount St. Helens on 18 May 1980 produced several catastrophic flowages into tributaries on the northeast volcano flank. The tributaries channeled the flows to Smith Creek valley, which lies within the area devastated by the surge but was unaffected by the great debris avalanche on the north flank. Stratigraphy shows that the pyroclastic surge preceded the lahars; there is no notable “wet” character to the surge deposits. Therefore the lahars must have originated as snowmelt, not as ejected water-saturated debris that segregated from the pyroclastic surge as has been inferred for other flanks of the volcano. In stratigraphic order the Smith Creek valley-floor materials comprise (1) a complex valley-bottom facies of the pyroclastic surge and a related pyroclastic flow, (2) an unusual hummocky diamict caused by complex mixing of lahars with the dry pyroclastic debris, and (3) deposits of secondary pyroclastic flows. These units are capped by silt containing accretionary lapilli, which began falling from a rapidly expanding mushroom-shaped cloud 20 minutes after the eruption's onset. The Smith Creek valley-bottom pyroclastic facies consists of (a) a weakly graded basal bed of fines-poor granular sand, the deposit of a low-concentration lithic pyroclastic surge, and (b) a bed of very poorly sorted pebble to cobble gravel inversely graded near its base, the deposit of a high-concentration lithic pyroclastic flow. The surge apparently segregated while crossing the steep headwater tributaries of Smith Creek; large fragments that settled from the turbulent surge formed a dense pyroclastic flow along the valley floor that lagged behind the front of the overland surge. The unusual hummocky diamict as thick as 15 m contains large lithic clasts supported by a tough, brown muddy sand matrix like that of lahar deposits upvalley. This unit contains irregular friable lenses and pods meters in diameter, blocks incorporated from the underlying dry and hot pyroclastic material that had been deposited only moments earlier. The hummocky unit is the deposit of a high-viscosity debris flow which formed when lahars mingled with the pyroclastic materials on Smith Creek valley floor. Overlying the debris flow are voluminous pyroclastic deposits of pebbly sand cut by fines-poor gas-escape pipes and containing charred wood. The deposits are thickest in topographic lows along margins of the hummocky diamict. Emplaced several minutes after the hot surge had passed, this is the deposit of numerous secondary pyroclastic flows derived from surge material deposited unstably on steep valley sides.     

Subglacial drumlins and englacial fractures at the surge-type glacier,Múlajökull,Iceland

The interaction between drumlins and overriding glacier ice is not well studied, largely due to the difficulty of identifying and accessing suitable active subglacial environments. The surge-type glacier Múlajökull, in central Iceland, overlies a known field of actively forming drumlins and therefore provides a rare opportunity to investigate the englacial structures that have developed in association with ice flow over the subglacial drumlins. In this study detailed ground-penetrating radar surveys are combined with field observations to identify clear sets of up-glacier and down-glacier dipping fractures at Múlajökull's margin. These are interpreted as conjugate shear planes or P- and R-type Reidel shears that developed and filled with saturated sediment derived from the glacier bed, during a previous surge. The fracture sets exhibit focused spatial distributions that are influenced by the subglacial topography. In particular, down-glacier dipping fractures are strongly focused over drumlin stoss slopes. These fractures, although well developed at depth, were mostly unable to transmit basal water and sediment up to the glacier surface during the surge cycle. In contrast, up-glacier dipping fractures formed over drumlin lee sides and in more gently sloping swales, and more frequently connected to the glacier surface, providing a pathway for the evacuation of basal water and water-saturated sediment. The study suggests that the subglacial drumlins under Múlajökull's margin have influenced the nature and distribution of englacial fractures, which could potentially contribute to spatial variations in basal water pressure during a surge. BGS © UKRI 2018     

Periodic drainage of ice-dammed lakes as a result of variations in glacier velocity

Previous discussions of the catastrophic drainage of ice-dammed lakes have centred on mechanisms where characteristics of the lake are crucial to drainage initiation, for example dam flotation or tunnel formation at a critical lake depth. This paper describes a mechanism for lake drainage where drainage initiation depends on the characteristics of the glacier and is independent of the characteristics of the lake. Prediction of this mechanism must be based on glacier dynamics, whereas the mechanisms most commonly discussed previously are best predicted primarily on the basis of lake evolution. An ice-dammed lake at the margin of the glacier Solheimajokull, in southern Iceland, was observed to drain rapidly into the sub- or englacial drainage system, supplying water and debris to the bed or interior of the glacier. Geomorphological evidence suggests that the lake drains and refills periodically, discharging up to 13300 m³ of water into the glacier-hydrological system. The depth of the maximum lake is insufficient to cause either flotation of the ice margin or tunnel opening by plastic deformation of the ice, and we suggest that sudden drainage is related to ice-bed separations associated with specific glacier flow states rather than to a critical lake depth threshold. This mechanism of lake drainage has implications for conditions at the glacier bed, for the development of basal ice and for the entrainment of debris into the glacier, as well as for the prediction of potentially hazardous catastrophic drainage events and jokulhlaups from ice-dammed lakes.     

Subglacial observations from øksfjordjøkelen,North Norway

Stepped bedrock topography at the snout of a small outlet glacier from Øksfjordjøkelen, North Norway, produces an extensive subglacial cavity system which stretches some 70m across and 100m up-glacier, giving access beneath ice ≤50 m thick. Inside the cavity, regelation ice, clean glacier ice and deforming basal ice have been observed. Samples were taken and basal debris concentrations at the glacier sole were found to vary between 0.005 and 15.38 per cent by volume. The basal ice velocity has been determined using a linear variable differential transformer attached to an analogue recorder, and also by means of measured displacements of ice crews and clasts embedded in the basal ice. Velocities were found to differ both spatially and temporally from a maximum of 2.55 mm h¹ to a minimum of 0.3 mm h^?1. The measurements and observations, which have been related to present theory, show how spatially averaged values for a number of variables could lead to inaccuracies in predicted erosion values, certainly at a local scale. On the exposed foreland, jointcontrolled lee-side faces provide evidence for extensive subglacial plucking (here taken to mean the removal of preloosened bed material and/or material resulting from bed failure). Indeed, in the cavity the early stages of removal of joint-controlled blocks by ice deformation along joints have been observed. The importance of debris-rich basal ice is shown in the formation of large striations (up to 500cm × 16cm × 2cm) present on the foreland.     

Modelling the hydrological response of debris‐free and debris‐covered glaciers to present climatic conditions in the semiarid Andes of central Chile

We apply the process‐based, distributed TOPKAPI‐ETH glacio‐hydrological model to a glacierized catchment (19% glacierized) in the semiarid Andes of central Chile. The semiarid Andes provides vital freshwater resources to valleys in Chile and Argentina, but only few glacio‐hydrological modelling studies have been conducted, and its dominant hydrological processes remain poorly understood. The catchment contains two debris‐free glaciers reaching down to 3900 m asl (Bello and Yeso glaciers) and one debris‐covered avalanche‐fed glacier reaching to 3200 m asl (Piramide Glacier). Our main objective is to compare the mass balance and runoff contributions of both glacier types under current climatic conditions. We use a unique dataset of field measurements collected over two ablation seasons combined with the distributed TOPKAPI‐ETH model that includes physically oriented parameterizations of snow and ice ablation, gravitational distribution of snow, snow albedo evolution and the ablation of debris‐covered ice. Model outputs indicate that while the mass balance of Bello and Yeso glaciers is mostly explained by temperature gradients, the Piramide Glacier mass balance is governed by debris thickness and avalanches and has a clear non‐linear profile with elevation as a result. Despite the thermal insulation effect of the debris cover, the mass balance and contribution to runoff from debris‐free and debris‐covered glaciers are similar in magnitude, mainly because of elevation differences. However, runoff contributions are distinct in time and seasonality with ice melt starting approximately four weeks earlier from the debris‐covered glacier, what is of relevance for water resources management. At the catchment scale, snowmelt is the dominant contributor to runoff during both years. However, during the driest year of our simulations, ice melt contributes 42 ± 8% and 67 ± 6% of the annual and summer runoff, respectively. Sensitivity analyses show that runoff is most sensitive to temperature and precipitation gradients, melt factors and debris cover thickness. Copyright © 2016 John Wiley & Sons, Ltd.     

Interannual variability in glacier contribution to runoff from a high‐elevation Andean catchment: understanding the role of debris cover in glacier hydrology

We present a field‐data rich modelling analysis to reconstruct the climatic forcing, glacier response, and runoff generation from a high‐elevation catchment in central Chile over the period 2000–2015 to provide insights into the differing contributions of debris‐covered and debris‐free glaciers under current and future changing climatic conditions. Model simulations with the physically based glacio‐hydrological model TOPKAPI‐ETH reveal a period of neutral or slightly positive mass balance between 2000 and 2010, followed by a transition to increasingly large annual mass losses, associated with a recent mega drought. Mass losses commence earlier, and are more severe, for a heavily debris‐covered glacier, most likely due to its strong dependence on snow avalanche accumulation, which has declined in recent years. Catchment runoff shows a marked decreasing trend over the study period, but with high interannual variability directly linked to winter snow accumulation, and high contribution from ice melt in dry periods and drought conditions. The study demonstrates the importance of incorporating local‐scale processes such as snow avalanche accumulation and spatially variable debris thickness, in understanding the responses of different glacier types to climate change. We highlight the increased dependency of runoff from high Andean catchments on the diminishing resource of glacier ice during dry years.     

Englacial vs lacustrine origin of volcanic table mountains: evidence from iceland

Detailed facies analysis of hyaloclastites and associated lavas from eight table mountains and similar "hyaloclastite volcanoes" in the Icelandic rift zone contradict a rapid and continuous, "monogenetic", entirely subglacial evolution of most volcanoes studied. The majority of the exposed hyaloclastite deposits formed in large, stable lakes as indicated by widespread, up to 300-m-thick, continuous sections of deep water, shallow water and emergent facies. Salient features include extensively layered or bedded successions comprising mainly debris flow deposits, turbidites, base surge and fallout deposits consisting of texturally and compositionally variable, slightly altered hyaloclastites, as well as sheet and pillow lavas. In contrast, chaotic assemblages of coarser-grained, more poorly sorted and more strongly palagonitized hyaloclastite tuffs and breccias, as well as scoria and lava are interpreted to have formed under sub- or englacial conditions in small, chimney-like ice cavities or ice-bound lakes. Irregularly shaped and erratically arranged hyaloclastite bodies produced at variable water levels appear to have resulted mainly from rapid changes of the eruptive environment due to repeated build-up and drainage of ice-bound lakes as well as the restricted space between the ice walls. We distinguish a "deep water" facies formed during high water levels of the lake, a hydroclastic shallow water and emergent facies (leakage of the lake or growth of the volcano above the water surface). Our model implies the temporary existence of large, stable lakes in Iceland probably formed by climatically induced ice melting. The highly complex edifices of many table mountains and similar volcanoes were constructed during several eruptive periods in changing environments characterized by contrasting volcanic and sedimentary processes. Received: 10 June 1997 / Accepted: 28 July 1998     

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.     

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.     

Ice flow‐unit influence on glacier structure,debris entrainment and transport

The links between structural glaciology, glacial debris entrainment and transport have been established in a number of different glacier settings. Here we document the structural evolution of a temperate Alpine valley glacier (Vadrec del Forno, Switzerland) and demonstrate that individual flow units within the glacier have very different structural and debris characteristics. The glacier consists of a broad accumulation area with multiple basins feeding a relatively narrow tongue and is formed from six distinct flow units. Each flow unit has its own characteristic structural assemblage. Flow units that narrow rapidly down‐glacier are dominated by primary stratification that has evolved into longitudinal foliation. In contrast, wider flow units preferentially develop an axial planar foliation. Glacier structure plays a limited role in the entrainment of debris, which is more strongly influenced by ice‐marginal rockfall and avalanche inputs onto the glacier surface. However, once entrained, glacier structure controls the reorientation and redistribution of debris within the ice mass. By taking a whole‐glacier approach to describing glacier structure and debris transport, we conclude that individual flow units are unique with regard to structure and debris transfer. Copyright © 2013 John Wiley & Sons, Ltd.     

Implications of decadal to century scale glacio‐hydrological change for water resources of the Hood River basin,OR, USA

In glacier‐fed rivers, melting of glacier ice sustains streamflow during the driest times of the year, especially during drought years. Anthropogenic and ecologic systems that rely on this glacial buffering of low flows are vulnerable to glacier recession as temperatures rise. We demonstrate the evolution of glacier melt contribution in watershed hydrology over the course of a 184‐year period from 1916 to 2099 through the application of a coupled hydrological and glacier dynamics model to the Hood River basin in Northwest Oregon, USA. We performed continuous simulations of glaciological processes (mass accumulation and ablation, lateral flow of ice and heat conduction through supra‐glacial debris), which are directly linked with seasonal snow dynamics as well as other key hydrologic processes (e.g. evapotranspiration and subsurface flow). Our simulations show that historically, the contribution of glacier melt to basin water supply was up to 79% at upland water management locations. We also show that supraglacial debris cover on the Hood River glaciers modulates the rate of glacier recession and progression of dry season flow at upland stream locations with debris‐covered glaciers. Our model results indicate that dry season (July to September) discharge sourced from glacier melt started to decline early in the 21st century following glacier recession that started early in the 20th century. Changes in climate over the course of the current century will lead to 14–63% (18–78%) reductions in dry season discharge across the basin for IPCC emission pathway RCP4.5 (RCP8.5). The largest losses will be at upland drainage locations of water diversions that were dominated historically by glacier melt and seasonal snowmelt. The contribution of glacier melt varies greatly not only in space but also in time. It displays a strong decadal scale fluctuations that are super‐imposed on the effects of a long‐term climatic warming trend. This decadal variability results in reversals in trends in glacier melt, which underscore the importance of long‐time series of glacio‐hydrologic analyses for evaluating the hydrological response to glacier recession. Copyright © 2016 John Wiley & Sons, Ltd.