Heterogeneity in supraglacial debris thickness and its role in glacier mass changes of the Mount Gongga

In the Tibetan Plateau, many glaciers have extensive covers of supraglacial debris in their ablation zones, which affects glacier response to climate change by altering ice melting and spatial patterns of mass loss. Insufficient debris thickness data make it difficult to analyze regional debris-cover effects. Maritime glaciers of the Mount Gongga have been characterized by a substantial reduction in glacier area and ice mass in recent decades. The thermal property of the debris layer estimated from remotely sensed data reveals that debris-covered glaciers are dominant in this region, on which the proportion of debris cover to total glacier area varies from 1.74% to 53.0%. Using a physically-based debris-cover effect assessment model, we found that although the presence of supraglacial debris has a significant insulating effect on heavily debris-covered glaciers, it accelerates ice melting on ~10.2% of total ablation zone and produces rapid wastage of ~25% of the debris-covered glaciers, leading to the similar mass losses between the debris-covered and debris-free glaciers. Widespread debris cover also facilitates the development of active terminus regions. Regional differences in debris-cover effects are apparent, highlighting the importance of debris cover for understanding glacier mass changes in the Tibetan Plateau and other mountain ranges around the world.     

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.     

Processes at the margins of supraglacial debris cover: Quantifying dirty ice ablation and debris redistribution

Current glacier ablation models have difficulty simulating the high-melt transition zone between clean and debris-covered ice. In this zone, thin debris cover is thought to increase ablation compared to clean ice, but often this cover is patchy rather than continuous. There is a need to understand ablation and debris dynamics in this transition zone to improve the accuracy of ablation models and the predictions of future debris cover extent. To quantify the ablation of partially debris-covered ice (or ‘dirty ice’), a high-resolution, spatially continuous ablation map was created from repeat unmanned aerial systems surveys, corrected for glacier flow in a novel way using on-glacier ablation stakes. Surprisingly, ablation is similar (range ~ 5 mm w.e. per day) across a wide range of percentage debris covers (~ 30–80%) due to the opposing effects of a positive correlation between percentage debris cover and clast size, countered by a negative correlation with albedo. Once debris cover becomes continuous, ablation is significantly reduced (by 61.6% compared to a partial debris cover), and there is some evidence that the cleanest ice (<~ 15% debris cover) has a lower ablation than dirty ice (by 3.7%). High-resolution feature tracking of clast movement revealed a strong modal clast velocity where debris was continuous, indicating that debris moves by creep down moraine slopes, in turn promoting debris cover growth at the slope toe. However, not all slope margins gain debris due to the removal of clasts by supraglacial streams. Clast velocities in the dirty ice area were twice as fast as clasts within the continuously debris-covered area, as clasts moved by sliding off their boulder tables. These new quantitative insights into the interplay between debris cover characteristics and ablation can be used to improve the treatment of dirty ice in ablation models, in turn improving estimates of glacial meltwater production. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

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.     

Evolution of debris cover on glaciers of the Eastern Alps,Austria, between 1996 and 2015

Debris cover on glaciers is an important component of glacial systems as it influences climate–glacier dynamics and thus the lifespan of glaciers. Increasing air temperatures, permafrost thaw and rock faces freshly exposed by glacier downwasting in accumulation zones result in increased rockfall activity and debris input. In the ablation zone, negative mass balances result in an enhanced melt-out of englacial debris. Glacier debris cover thus represents a clear signal of climate warming in mountain areas. To assess the temporal development of debris on glaciers of the Eastern Alps, Austria, we mapped debris cover on 255 glaciers using Landsat data at three time steps. We applied a ratio-based threshold classification technique and analysed glacier catchment characteristics to understand debris sources better. Across the Austrian Alps, debris cover increased by more than 10% between 1996 and 2015 while glaciers retreated in response to climate warming. Debris cover distribution shows significant regional variability, with some mountain ranges being characterised by mean debris cover on glaciers of up to 75%. We also observed a general rise of the mean elevation of debris cover on glaciers in Austria. The debris cover distribution and dynamics are highly variable due to topographic, lithological and structural settings that determine the amount of debris delivered to and stored in the glacier system. Despite strong variation in debris cover, all glaciers investigated melted at increasing rates. We conclude that the retarding effects of debris cover on the mass balance and melt rate of Austrian glaciers is strongly subdued compared with other mountain areas. The study indicates that, if this trend continues, many glaciers in Austria may become fully debris covered. However, since debris cover seems to have little impact on melt rates, this would not lead to prolonged existence of debris-covered ice compared with clean ice glaciers.     

Stability of supraglacial debris

Rock debris on the surface of ablating glaciers is not static, and is often transported across the ice surface as relief evolves during melt. This supraglacial debris transport has a strong influence on the spatial distribution of melt, and is implicated in the formation of hummocky glacial topography in deglaciated terrain. Furthermore, as ice‐dammed lakes and ice‐cored slopes become increasingly common in deglaciating watersheds, there is rising concern about hazards to humans and infrastructure posed by mass‐wasting of ice‐cored debris. The existing quantitative framework for describing these debris transport processes is limited, making it difficult to account for transport in mass balance, hazard assessment, and landscape development models. This paper develops a theoretical framework for assessing slope stability and gravitational mass transport in a debris‐covered ice setting. Excess water pressure at the interface between ablating ice and lowering debris is computed by combining Darcy's law with a meltwater balance. A limit‐equilibrium slope stability analysis is then applied to hypothetical debris layers with end‐member moisture conditions derived from a downslope meltwater balance that includes production and seepage. The resulting model system constrains maximum stable slope angles and lengths that vary with debris texture, thickness, and the rate of meltwater production. Model predictions are compared with field observations and with digital elevation model (DEM)‐derived terrain metrics from two modern debris‐covered glaciers on Mount Rainier, USA. Copyright © 2017 John Wiley & Sons, Ltd.     

The evolution of brittle and ductile structures at the surface of a partly debris-covered,rapidly thinning and slowly moving glacier in 1998–2012 (Pasterze Glacier,Austria)

Many glaciers in alpine regions are currently rapidly receding and thinning at historically unobserved rates causing changes in the velocity field and in normal and shear stresses affecting the surface expression of structures within the ice. We studied the distribution of brittle and ductile structures at the surface of Pasterze Glacier during a 14-year period by analysing orthophotos and digital elevation models of five stages (1998, 2003, 2006, 2009 and 2012). A structural glaciological mapping key was applied. Normal faults, strike-slip faults, en échelon structures (systematic stepping of fractures), thrust faults, and band ogives were distinguished. Results indicate substantial deceleration and glacier thinning in 1998–2012. Glacier thinning was not homogenous over time related to the uneven distribution of supraglacial debris causing differential ablation or the selective ablation effects of subglacial water channels. Peculiar supraglacial features observed are circular collapse structures with concentric crevasses which form when the ice between the surface and the roof of water channels decreases. The total length of brittle structures increased from 38.4 km to 56.9 km whereas the extent of the glacier tongue decreased by 25%. The fracture density doubled from 0.009 to 0.018 m/m². Areas of the glacier tongue which were up to 100 m away from the nearest brittle structure increased by 16%. The visual appearance of thrust faults shifted upglacier due to decreasing glacier velocity causing horizontal shortening or due to exhumation of faults that did not previously extend to the surface. A large number of brittle structures are progressively independent from glacier motion. Our study suggests that glacier tongues which are in a state of rapid decay and thinning are prone to fracturing due to normal fault formation and glacier disintegration. Water further increases ablation rates substantially if rather large amounts drain through supra-, en- or subglacial water channels. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Geomorphological evolution of a debris-covered glacier surface

There exists a need to advance our understanding of debris-covered glacier surfaces over relatively short timescales due to rapid, climatically induced areal expansion of debris cover at the global scale, and the impact debris has on mass balance. We applied unpiloted aerial vehicle structure-from-motion (UAV-SfM) and digital elevation model (DEM) differencing with debris thickness and debris stability modelling to unravel the evolution of a 0.15 km² region of the debris-covered Miage Glacier, Italy, between June 2015 and July 2018. DEM differencing revealed widespread surface lowering (mean 4.1 ± 1.0 m a^-1; maximum 13.3 m a^-1). We combined elevation change data with local meteorological data and a sub-debris melt model, and used these relationships to produce high resolution, spatially distributed maps of debris thickness. These maps were differenced to explore patterns and mechanisms of debris redistribution. Median debris thicknesses ranged from 0.12 to 0.17 m and were spatially variable. We observed localized debris thinning across ice cliff faces, except those which were decaying, where debris thickened. We observed pervasive debris thinning across larger, backwasting slopes, including those bordered by supraglacial streams, as well as ingestion of debris by a newly exposed englacial conduit. Debris stability mapping showed that 18.2–26.4% of the survey area was theoretically subject to debris remobilization. By linking changes in stability to changes in debris thickness, we observed that slopes that remain stable, stabilize, or remain unstable between periods almost exclusively show net debris thickening (mean 0.07 m a^-1) whilst those which become newly unstable exhibit both debris thinning and thickening. We observe a systematic downslope increase in the rate at which debris cover thickens which can be described as a function of the topographic position index and slope gradient. Our data provide quantifiable insights into mechanisms of debris remobilization on glacier surfaces over sub-decadal timescales, and open avenues for future research to explore glacier-scale spatiotemporal patterns of debris remobilization. © 2020 The Authors. Earth Surface Processes and Landforms published by 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.     

Organic carbon fluxes of a glacier surface: A case study of Foxfonna,a small Arctic glacier

Arctic glaciers are rapidly responding to global warming by releasing organic carbon (OC) to downstream ecosystems. The glacier surface is arguably the most biologically active and biodiverse glacial habitat and therefore the site of important OC transformation and storage, although rates and magnitudes are poorly constrained. In this paper, we present measurements of OC fluxes associated with atmospheric deposition, ice melt, biological growth, fluvial transport and storage (in superimposed ice and cryoconite debris) for a supraglacial catchment on Foxfonna glacier, Svalbard (Norway), across two consecutive years. We found that in general atmospheric OC input (averaging 0.63 ± 0.25 Mg a^-1 total organic carbon, i.e. TOC, and 0.40 ± 0.22 Mg a^-1 dissolved organic carbon, i.e. DOC) exceeded fluvial OC export (0.46 ± 0.04 Mg a^-1 TOC and 0.36 ± 0.03 Mg a^-1 DOC). Early in the summer, OC was mobilised in snowmelt but its release was delayed by temporary storage in superimposed ice on the glacier surface. This delayed the export of 28.5% of the TOC in runoff. Biological production in cryoconite deposits was a negligible potential source of OC to runoff, while englacial ice melt was far more important on account of the glacier's negative ice mass balance (–0.89 and –0.42 m a^-1 in 2011 and 2012, respectively). However, construction of a detailed OC budget using these fluxes shows an excess of inputs over outputs, resulting in a net retention of OC on the glacier surface at a rate that would require c. 3 years to account for the OC stored as cryoconite debris. © 2018 John Wiley & Sons, Ltd.     

An enhanced temperature index model for debris-covered glaciers accounting for thickness effect

Debris-covered glaciers are increasingly studied because it is assumed that debris cover extent and thickness could increase in a warming climate, with more regular rockfalls from the surrounding slopes and more englacial melt-out material. Debris energy-balance models have been developed to account for the melt rate enhancement/reduction due to a thin/thick debris layer, respectively. However, such models require a large amount of input data that are not often available, especially in remote mountain areas such as the Himalaya, and can be difficult to extrapolate. Due to their lower data requirements, empirical models have been used extensively in clean glacier melt modelling. For debris-covered glaciers, however, they generally simplify the debris effect by using a single melt-reduction factor which does not account for the influence of varying debris thickness on melt and prescribe a constant reduction for the entire melt across a glacier.In this paper, we present a new temperature-index model that accounts for debris thickness in the computation of melt rates at the debris-ice interface. The model empirical parameters are optimized at the point scale for varying debris thicknesses against melt rates simulated by a physically-based debris energy balance model. The latter is validated against ablation stake readings and surface temperature measurements. Each parameter is then related to a plausible set of debris thickness values to provide a general and transferable parameterization. We develop the model on Miage Glacier, Italy, and then test its transferability on Haut Glacier d’Arolla, Switzerland.The performance of the new debris temperature-index (DETI) model in simulating the glacier melt rate at the point scale is comparable to the one of the physically based approach, and the definition of model parameters as a function of debris thickness allows the simulation of the nonlinear relationship of melt rate to debris thickness, summarised by the Østrem curve. Its large number of parameters might be a limitation, but we show that the model is transferable in time and space to a second glacier with little loss of performance. We thus suggest that the new DETI model can be included in continuous mass balance models of debris-covered glaciers, because of its limited data requirements. As such, we expect its application to lead to an improvement in simulations of the debris-covered glacier response to climate in comparison with models that simply recalibrate empirical parameters to prescribe a constant across glacier reduction in melt.     

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.     

Using structure‐from‐motion to create glacier DEMs and orthoimagery from historical terrestrial and oblique aerial imagery

Increased resolution and availability of remote sensing products, and advancements in small‐scale aerial drone systems, allows observations of glacial changes at unprecedented levels of detail. Software developments, such as structure‐from‐motion (SfM), now allow users an easy and efficient method to generate three‐dimensional (3D) models and orthoimages from aerial or terrestrial datasets. While these advancements show promise for current and future glacier monitoring, many regions still suffer a lack of observations from earlier time periods. We report on the use of SfM to extract spatial information from various historic imagery sources. We focus on three geographic regions, the European Alps, high Arctic Norway and the Nepal Himalayas. We used terrestrial field photographs from 1896, high oblique aerial photographs from 1936 and aerial handheld photographs from 1978 to generate digital elevation models (DEMs) and orthophotos of the Rhone glacier, Brøggerhalvøya and the lower Khumbu glacier, respectively. Our analysis shows that applying SfM to historic imagery can generate high quality models using only ground control points. Limited camera/orientation information was largely reproduced using self‐calibrated model data. Using these data, we calculated mean ground sampling distances across each site which demonstrates the high potential resolution of resulting models. Vertical errors for our models are ±5.4 m, ±5.2 m and ±3.3 m. Differencing shows similar patterns of thinning at lower Rhone (European Alps) and Brøggerhalvøya (Norway) glaciers, which have mean thinning rates of 0.31 m a^?1 (1896–2010) to 0.86 m a^?1 (1936–2010) respectively. On these clean ice glaciers thinning is highest in the terminus region and decreasing up‐glacier. In contrast to these glaciers, uneven topography, exposed ice‐cliffs and debris cover on the Khumbu glacier create a highly variable spatial distribution of thinning. The mean thinning rate for the Khumbu study area was found to be 0.54 ± 0.9 m a^?1 (1978–2015). Copyright © 2017 John Wiley & Sons, Ltd.     

Occurrence,evolution and ice content of ice-debris complexes in the Ak-Shiirak,Central Tien Shan revealed by geophysical and remotely-sensed investigations

Rock glaciers and large ice-debris complexes are common in many mountain ranges and are especially prominent in semi-arid mountains such as the Andes or the Tien Shan. These features contain a significant amount of ice but their occurrence and evolution are not well known. Here, we present an inventory of the ice-debris complexes for the Ak-Shiirak, Tien Shan's second largest glacierised massif, and a holistic methodology to investigate two characteristic and large ice-debris complexes in detail based on field investigations and remote sensing analysis using Sentinel-1 SAR data, 1964 Corona and recent high resolution stereo images. Overall, we found 74 rock glaciers and ice-debris complexes covering an area of 11.2 km² (3.2% of the glacier coverage) with a mean elevation of about 3950 m asl. Most of the complexes are located south-east of the main ridge of Ak-Shiirak. Ground penetrating radar (GPR) measurements reveal high ice content with the occurrence of massif debris-covered dead-ice bodies in the parts within the Little Ice Age glacier extent. These parts showed significant surface lowering, in some places exceeding 20 m between 1964 and 2015. The periglacial parts are characterised by complex rock glaciers of different ages. These rock glaciers could be remnants of debris-covered ice located in permafrost conditions. They show stable surface elevations with no or only very low surface movement. However, the characteristics of the fronts of most rock glacier parts indicate slight activity and elevation gains at the fronts slight advances. GPR data indicated less ice content and slanting layers which coincide with the ridges and furrows and could mainly be formed by glacier advances under permafrost conditions. Overall, the ice content is decreasing from the upper to the lower part of the ice-debris complexes. Hence, these complexes, and especially the glacier-affected parts, should be considered when assessing the hydrological impacts of climate change. © 2018 John Wiley & Sons, Ltd.     

Hydrology and dynamics of a polythermal (mostly cold) High Arctic glacier

To improve our understanding of the interactions between hydrology and dynamics in mostly cold glaciers (in which water flow is limited by thermal regime), we analyse short‐term (every two days) variations in glacier flow in the ablation zone of polythermal John Evans Glacier, High Arctic Canada. We monitor the spatial and temporal propagation of high‐velocity events, and examine their impacts upon supraglacial drainage processes and evolving subglacial drainage system structure. Each year, in response to the rapid establishment of supraglacial–subglacial drainage connections in the mid‐ablation zone, a ‘spring event’ of high horizontal surface velocities and high residual vertical motion propagates downglacier over two to four days from the mid‐ablation zone to the terminus. Subsequently, horizontal velocities fall relative to the spring event but remain higher than over winter, reflecting channelization of subglacial drainage but continued supraglacial meltwater forcing. Further transient high‐velocity events occur later in each melt season in response to melt‐induced rising supraglacial meltwater inputs to the glacier bed, but the dynamic response of the glacier contrasts with that recorded during the spring event, with the degree of spatial propagation a function of the degree to which the subglacial drainage system has become channelized. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Hummock corridors in the south‐central sector of the Fennoscandian ice sheet,morphometry and pattern

Subglacial conditions strongly influence the flow of ice‐sheets, in part due to the availability of melt water. Contemporary ice sheets are retreating and are affected by increased melting as climate warms. The south Swedish uplands (SSU) were deglaciated during the relatively warm Bølling‐Allerød interval, and by studying the glacial landforms there it is possible to increase the understanding of the subglacial environment during this period of warming. Across the study area, vast tracts of hummocks have long been recognized. However, recent mapping shows a pattern of elongated zones of hummocks radially oriented, hereafter referred to as ‘hummock corridors’. Morphometric parameters were measured on the hummock corridors using a 2 m horizontal resolution digital elevation model. Corridor width varies between 0.2 and 4.9 km and their length between 1.5 and 11.8 km. A majority of hummock corridors are incised in drumlinised till surfaces. The pattern of hummock corridors shows a clear relation to the overall ice‐flow. Further, hummock corridors do not follow topographic gradients, and in at least one place an esker overlies hummocks on the corridor floor. The lateral spacing of hummock corridors and corridor morphology are similar to tunnel valleys, eskers and glaciofluvial corridors reported elsewhere. Such relationships support a subglacial genesis of the corridors in the SSU by water driven by the subglacial hydraulic gradient and that hummock corridors are forms that can be identified as tunnel valleys and glaciofluvial corridors (GFC). Ages were assigned to hummock‐corridor cross‐sections from a deglacial reconstruction of the Fennoscandian Ice Sheet. By comparing the frequency of corridors per age interval with climate variations from a Greenland ice core, we hypothesize that an increase in the number of corridors is related to the Bølling‐Allerød warming, indicating a higher rate of delivery of surface melt water to the bed at this time. Copyright © 2017 John Wiley & Sons, Ltd.     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

Flow prediction in high altitude ungauged catchments: A case study in the Italian Alps (Pantano Basin,Adamello Group)

We develop a simple model to evaluate the daily flow discharges in the ablation season for the 11 km² Pantano basin in the Retiche Italian Alps, based upon the data gathered during a three years field campaign. The Pantano basin embeds the Venerocolo debris covered and the Avio debris free glaciers, covering 2.14 km² in the Adamello Group, where the widest Italian glacier Adamello is located. First, degree-day models based upon air temperature are tuned to calculate snow and ice melt at daily scale. Glaciers’ meteorological data are collected from an automatic weather station (AWS), operating on the glacier during summer 2007. The melt factors in the debris covered areas of the glacier are estimated against debris thickness, using a data driven parameterization. The flow discharge from the catchment is estimated using semi distributed flow routing for the ablation seasons of four years, from 2006 to 2009. The predicted discharges are compared to those derived from inverse reservoir's routing at the Benedetto lake, catching the basin outflow. The proposed approach is valuable as a tool to investigate the hydrology of poorly gauged glacierized areas, including those with debris covered ice, widely diffused and yet poorly understood. Pending accurate parameterization the approach is usable for water resources evaluation and for long term assessment of the climate change impact on the glacierized areas within the Alps.