Initial results from a distributed,physically based model of glacier hydrology

This paper describes the development and testing of a distributed, physically based model of glacier hydrology. The model is used to investigate the behaviour of the hydrological system of Haut Glacier d'Arolla, Valais, Switzerland. The model has an hourly time-step and three main components: a surface energy balance submodel, a surface flow routing submodel and a subglacial hydrology submodel. The energy balance submodel is used to calculate meltwater production over the entire glacier surface. The surface routing submodel routes meltwater over the glacier surface from where it is produced to where it either enters the subglacial hydrological system via moulins or runs off the glacier surface. The subglacial hydrology submodel calculates water flow in a network of conduits, which can evolve over the course of a melt season simulation in response to changing meltwater inputs. The main model inputs are a digital elevation model of the glacier surface and its surrounding topography, start-of-season snow depth distribution data and meteorological data. Model performance is evaluated by comparing predictions with field measurements of proglacial stream discharge, subglacial water pressure (measured in a borehole drilled to the glacier bed) and water velocities inferred from dye tracer tests. The model performs best in comparison with the measured proglacial stream discharges, but some of the substantial features of the other two records are also reproduced. In particular, the model results show the high amplitude water pressure cycles observed in the borehole in the mid-melt season and the complex velocity/discharge hysteresis cycles observed in dye tracer tests. The results show that to model outflow hydrographs from glacierized catchments effectively, it is necessary to simulate spatial and temporal variations in surface melt rates, the delaying effect of the surface snowpack and the configuration of the subglacial drainage system itself. The model's ability to predict detailed spatial and temporal patterns of subglacial water pressures and velocities should make it a valuable tool for aiding the understanding of glacier dynamics and hydrochemistry. © 1998 John Wiley & Sons, Ltd.     

Using synchronous fluorescence spectroscopy and principal components analysis to monitor dissolved organic matter dynamics in a glacier system

The molecular characteristics of dissolved organic matter (DOM) reflect both its source material and its biogeochemical history. In glacial systems, DOM characteristics might be expected to change over the course of a melt season as changes in the glacier drainage system cause the mobilization of DOM from different OM pools. To test this hypothesis we used Principal Components Analysis (PCA) of synchronous fluorescence spectra to detect and describe changes in the DOM in meltwater from a glacier system in the Coast Mountains of northern British Columbia, Canada. For most of the melt season, the dominant component of subglacially routed meltwater DOM is characterized by a tyrosine‐like fluorophore. This DOM component is most likely derived from supraglacial snowmelt. During periods of high discharge, a second component of DOM is present which is humic in character and similar to DOM sampled from a nearby non‐glacial stream. This DOM component is inferred to be derived from a moss‐covered soil environment that has been glacially overrun. It is probably entrained into glacial melt waters when the supraglacial meltwater flux exceeds the capacity of the principal subglacial drainage channels and water floods areas of the glacier bed that are normally isolated from the subglacial drainage system. Another source of DOM also appears to be mobilized during periods of high air temperatures. It is characterized by both humic and proteinaceous fluorophores and may be derived from the drainage of supraglacial cryoconite holes. Copyright © 2009 John Wiley & Sons, Ltd.     

Basal sediment evacuation by subglacial meltwater: suspended sediment transport from Haut Glacier d'Arolla,Switzerland

Proglacial suspended sediment transport was monitored at Haut Glacier d'Arolla, Switzerland, during the 1998 melt season to investigate the mechanisms of basal sediment evacuation by subglacial meltwater. Sub‐seasonal changes in relationships between suspended sediment transport and discharge demonstrate that the structure and hydraulics of the subglacial drainage system critically influenced how basal sediment was accessed and entrained. Under hydraulically inefficient subglacial drainage at the start of the melt season, sediment availability was generally high but sediment transport increased relatively slowly with discharge. Later in the melt season, sediment transport increased more rapidly with discharge as subglacial meltwater became confined to a spatially limited network of channels following removal of the seasonal snowpack from the ablation area. Flow capacity is inferred to have increased more rapidly with discharge within subglacial channels because rapid changes in discharge during highly peaked diurnal runoff cycles are likely to have been accommodated largely by changes in flow velocity. Basal sediment availability declined during channelization but increased throughout the remainder of the monitored period, resulting in very efficient basal sediment evacuation over the peak of the melt season. Increased basal sediment availability during the summer appears to have been linked to high diurnal water pressure variation within subglacial channels inferred from the strong increase in flow velocity with discharge. Basal sediment availability therefore appears likely to have been increased by (1) enhanced local ice‐bed separation leading to extra‐channel flow excursions and[sol ]or (2) the deformation of basal sediment towards low‐pressure channels due to a strong diurnally reversing hydraulic gradient between channels and areas of hydraulically less‐efficient drainage. Copyright © 2005 John Wiley & Sons, Ltd.     

Glacial conditions that contribute to the regeneration of Fountain Glacier proglacial icing,Bylot Island,Canada

Proglacial icings are one of the most common forms of extrusive ice found in the Canadian Arctic. However, the icing adjacent to Fountain Glacier, Bylot Island, is unique due to its annual cycle of growth and decay, and perennial existence without involving freezing point depression of water due to chemical characteristics. Its regeneration depends on the availability of subglacial water and on the balance between ice accretion and hydro‐thermal erosion. The storage and conduction of the glacial meltwater involved in the accretion of the icing were analyzed by conducting topographic and ground penetrating radar surveys in addition to the modelling of the subglacial drainage network and the thermal characteristics of the glacier base. The reflection power analysis of the geophysical data shows that some areas of the lower ablation zone have a high accumulation of liquid water, particularly beneath the centre part of the glacier along the main supraglacial stream. A dielectric permittivity model of the glacier – sediment interface suggests that a considerable portion of the glacier is warm based; allowing water to flow through unfrozen subglacial sediments towards the proglacial outwash plain. All these glacier‐related characteristics contribute to the annual regeneration of the proglacial icing and allow for portions of the icing to be perennial. Copyright © 2013 John Wiley & Sons, Ltd.     

Evidence for seasonal subglacial outburst events at a polythermal glacier,Finsterwalderbreen, Svalbard

Bulk runoff and meteorological data suggest the occurrence of two meltwater outburst events at Finsterwalderbreen, Svalbard, during the 1995 and 1999 melt seasons. Increased bulk meltwater concentrations of Cl^? during the outbursts indicate the release of snowmelt from storage. Bulk meltwater hydrochemical data and suspended sediment concentrations suggest that this snowmelt accessed a chemical weathering environment characterized by high rock:water ratios and long rock–water contact times. This is consistent with a subglacial origin. The trigger for both the 1995 and 1999 outbursts is believed to be high rates of surface meltwater production and the oversupply of meltwater to areas of the glacier bed that were at the pressure melting point, but which were unconnected to the main subglacial drainage network. An increase in subglacial water pressure to above the overburden pressure lead to the forcing of a hydrological connection between the expanding subglacial reservoir and the ice‐marginal channelized system. The purging of ice blocks from the glacier during the outbursts may indicate the breach of an ice dam during connection. Although subglacial meltwater issued continually from the glacier terminus via a subglacial upwelling during both melt seasons, field observations showed outburst meltwaters were released solely via an ice‐marginal channel. It is possible that outburst events are a seasonal phenomenon at this glacier and reflect the periodic drainage of meltwaters from the same subglacial reservoir from year to year. However, the location of this reservoir is uncertain. A 100 m high bedrock ridge traverses the glacier 6·5 km from its terminus. The overdeepened area up‐glacier from this is the most probable site for subglacial meltwater accumulation. Copyright © 2001 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.     

Hydrochemistry of meltwaters draining a polythermal‐based,high‐Arctic glacier,south Svalbard: II. Winter and early Spring

The hydrochemistry of naled and upwelling water sampled from the forefields of Finsterwalderbreen, Svalbard, during spring are used for the first time to infer the hydrology of overwinter meltwaters at a polythermal‐based glacier. Hydrochemical variations in naled are explained in terms of different water sources and their chemical alteration during freezing. Two water sources to naled are identified: surficially routed snowmelt and subglacial water. Naled that results from the freezing of the former is enriched in atmospherically derived ions such as Na⁺ and Cl⁻, and is believed to be formed during winter warm periods. Naled of subglacial origin contains relatively high proportions of crustally derived solute. It reflects the freezing of subglacial meltwaters that continue to issue from a subterranean upwellling during winter. An increasing dominance of SO²⁻₄ Mg²⁺, Na⁺ and Cl⁻ in subglacial naled with increasing distance from the upwelling reflects the progressive freezing of this water body and the associated removal of Ca²⁺ and HCO by calcite precipitation. These spatial trends are accentuated by the leaching of soluble ions from the naled close to its source by subsequent upwelling waters. The chemistry of spring upwelling waters, also of subglacial origin, strongly reflects this process. Meltwater produced by geothermal heating of glacier basal ice is believed to be the principal source of water to the subglacial drainage system during winter. Solute acquisition by this meltwater is limited by a scarcity of proton suppliers. Evolution of this dilute meltwater carries an imprint of ion exchange processes. Some stored subglacial water from the end of the previous ablation season may supplement the basal meltwater component in early winter. Copyright © 2000 John Wiley & Sons, Ltd.     

Seasonal changes in the morphology of the subglacial drainage system,Haut Glacier d'Arolla,Switzerland

Dye tracing techniques were used to investigate the glacier-wide pattern of change in the englacial/subglacial drainage system of Haut Glacier d'Arolla during the ablation seasons of 1990 and 1991. Analysis of breakthrough curve characteristics indicate that over the course of a melt season, a system of major channels developed by headward growth at the expense of a hydraulically inefficient distributed system. By the end of the melt season, this channel system extended at least 3·3 km from the snout of the 4 km long glacier and drained the bulk of supraglacially derived meltwater passing through the glacier. The upper limit of the channel system closely followed the retreating snowline up-glacier. Rates of headward channel growth reached c. 65 m d⁻¹, although these rates decreased in the upper 1 km of the glacier where snowline retreat exposed a patchy firn aquifer. It appears that the removal of snow (with its high albedo and significant water storage capacity) from the glacier surface resulted in a dramatic increase in the volume of runoff into moulins, and in the peakedness of daily runoff cycles. This induced transient high water pressures within the distributed drainage system, which caused it to evolve rapidly into a channelised system. It is therefore likely that, at a local scale, channel growth occurred down-glacier from moulins, and that the overall up-glacier-directed pattern of channel formation was caused by the retreating snowline exposing new moulins and crevasses to inputs of ice-derived meltwater. Damping of diurnal melt inputs by storage in the firm aquifer accounts for the slowing of channel growth in the upper glacier. © 1998 John Wiley & Sons, Ltd.     

Spatio‐temporal tracer variability in the glacier melt end‐member — How does it affect hydrograph separation results?

Geochemical and isotopic tracers were often used in mixing models to estimate glacier melt contributions to streamflow, whereas the spatio‐temporal variability in the glacier melt tracer signature and its influence on tracer‐based hydrograph separation results received less attention. We present novel tracer data from a high‐elevation catchment (17 km², glacierized area: 34%) in the Oetztal Alps (Austria) and investigated the spatial, as well as the subdaily to monthly tracer variability of supraglacial meltwater and the temporal tracer variability of winter baseflow to infer groundwater dynamics. The streamflow tracer variability during winter baseflow conditions was small, and the glacier melt tracer variation was higher, especially at the end of the ablation period. We applied a three‐component mixing model with electrical conductivity and oxygen‐18. Hydrograph separation (groundwater, glacier melt, and rain) was performed for 6 single glacier melt‐induced days (i.e., 6 events) during the ablation period 2016 (July to September). Median fractions (±uncertainty) of groundwater, glacier melt, and rain for the events were estimated at 49±2%, 35±11%, and 16±11%, respectively. Minimum and maximum glacier melt fractions at the subdaily scale ranged between 2±5% and 76±11%, respectively. A sensitivity analysis showed that the intraseasonal glacier melt tracer variability had a marked effect on the estimated glacier melt contribution during events with large glacier melt fractions of streamflow. Intra‐daily and spatial variation of the glacier melt tracer signature played a negligible role in applying the mixing model. The results of this study (a) show the necessity to apply a multiple sampling approach in order to characterize the glacier melt end‐member and (b) reveal the importance of groundwater and rainfall–runoff dynamics in catchments with a glacial flow regime.     

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.     

Seasonal changes in ground‐penetrating radar signature observed at a polythermal glacier,Bylot Island,Canada

In recent years, ground‐penetrating radar (GPR) has been increasingly used for characterization of subglacial and englacial environments at polythermal glaciers. The geophysical method is able to exploit the dielectric difference between water, air, sediment and ice, allowing delineation of subsurface hydrological, thermal and structural conditions. More recent GPR research has endeavoured to examine temporal change in glaciers, in particular the distribution of the cold ice zone at polythermal glaciers. However, the exact nature of temporal change that can be identified using GPR has not been fully examined. This research presents the results of three GPR surveys conducted over the course of a summer ablation season at a polythermal glacier in the Canadian Arctic. A total of approximately 30 km of GPR profiles were collected in 2002 repeatedly covering the lower 2 km of Stagnation Glacier, Bylot Island (72°58′ N 78°22′ W). Comparison between profiles indicated changes in the radar signature, including increased noise, appearance and disappearance of englacial reflections, and signal attenuation in the latter survey. Further, an area of chaotic returns in up‐glacier locations, which was interpreted to be a wet temperate ice zone, showed marked recession over the course of the ablation season. Combining all the temporal changes that were detected by GPR, results indicate that a polythermal glacier may exhibit strongly seasonal changes in hydrological and thermal characteristics throughout the ice body, including the drainage of 17 000 m³ of temporarily stored intra‐glacial meltwater. It is also proposed that the liquid water content in the temperate ice zone of polythermal glaciers can be described as a fraction of a specific retention capacity. Copyright © 2006 John Wiley & Sons, Ltd.     

Trickle or treat: The dynamics of nutrient export from polar glaciers

Cold‐based polar glacier watersheds contain well‐defined supraglacial, ice‐marginal, and proglacial elements that differ in their degree of hydrologic connectivity, sources of water (e.g., snow, ice, and/or sediment pore water), meltwater residence times, allochthonous and autochthonous nutrient, and sediment loads. We investigated 11 distinct hydrological units along the supraglacial, ice marginal, and proglacial flow paths that drain Joyce Glacier in the McMurdo Dry Valleys of Antarctica. We found that these units play unique and important roles as sources and/or sinks for dissolved inorganic nitrogen and dissolved inorganic phosphorus and for specific fractions of dissolved organic matter (DOM) as waters are routed from the glacier into nutrient‐poor downstream ecosystems. Changes in nutrient export from the glacial system as a whole were observed as the routing and residence times of meltwater changed throughout the melt season. The concentrations of major ions in the proglacial stream were inversely proportional to discharge, such that there was a relatively constant “trickle” of these solutes into downstream ecosystems. In contrast, NO₃^? concentrations generally increased with discharge, resulting in delivery of episodic pulses of dissolved inorganic nitrogen‐rich water (“treats”) into those same ecosystems during high discharge events. DOM concentrations or fluorescence did not correlate with discharge rate, but high variability in DOM concentrations or fluorescence suggests that DOM may be exported downstream as episodic treats, but with spatial and/or temporal patterns that remain poorly understood. The strong, nutrient‐specific responses to changes in hydrology suggest that polar glacier drainage systems may export meltwater with nutrient compositions that vary within and between melt seasons and watersheds. Because nutrient dynamics identified in this study differ between glacier watersheds with broadly similar hydrology, climate, and geology, we emphasize the need to develop conceptual models of nutrient export that thoroughly integrate the biogeochemical and hydrological processes that control the sources, fate, and export of nutrients from each system.     

Subglacial drainage system structure and morphology of Brewster Glacier,New Zealand

A global positioning system and ground penetrating radar surveys is used to produce digital elevation models of the surface and bed of Brewster Glacier. These are used to derive maps of subglacial hydraulic potential and drainage system structure using three different assumptions about the subglacial water pressure (P_w): (i) P_w = ice overburden; (ii) P_w = half ice overburden; (iii) P_w = atmospheric. Additionally, 16 dye‐tracing experiments at 12 locations were performed through a summer melt season. Dye return curve shape, together with calculations of transit velocity, dispersivity and storage, are used to infer the likely morphology of the subglacial drainage system. Taken together, the data indicate that the glacier is underlain by a channelised but hydraulically inefficient drainage system in the early summer in which water pressures are close to ice overburden. By mid‐summer, water pressures are closer to half‐ice overburden and the channelised drainage system is more hydraulically efficient. Surface streams that enter the glacier close to the location of major subglacial drainage pathways are routed quickly to the channels and then to the glacier snout. Streams that enter the glacier further away from the drainage pathways are routed slowly to the channels and then to the snout because they first flow through a distributed drainage system. Copyright © 2008 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.     

ISOTOPIC AND IONIC VARIATIONS IN GLACIER RIVER WATER DURING THREE CONTRASTING ABLATION SEASONS

The significance of the baseflow component of glacier river discharge in summer varies with geographical location, altitude, glacier geometry and glacier size. Baseflow is maintained by meltwater generated above the transient equilibrium line and by water released from temporary storage on, in or beneath the glacier. At the Norwegian glacier Austre Okstindbreen, where precipitation is generally high throughout the year and the summers are cool and wet, observations in three successive, but contrasting, years have shown that Na⁺ ion concentrations in the glacier river water are influenced strongly by the amount of snowmelt. This itself depends on the preceding winter conditions, which determine the amount of accumulation, and on the current summer's weather. The efficiency of the glacier's drainage systems depends on the general progress of summer ablation. The speed with which the systems develop influences ion provision from subglacial sources. Ca²⁺ ion concentrations are largely determined by subglacial conditions. Oxygen isotope variations in glacier river water reflect the relative contributions made to total discharge by snow meltwater and other sources; the composition of the snow cover, which is a function of winter temperatures, has a strong influence. Ice meltwater has low isotopic variability, but the isotopic composition of rainfall varies markedly. A simple model of mixing of englacial and subglacial waters, each of a constant composition, cannot be applied to a high-latitude glacier of the size and altitudinal range of Austre Okstindbreen.     

Effect of snowpack removal on energy balance,melt and runoff in a small supraglacial catchment

Modelling melt and runoff from snow‐ and ice‐covered catchments is important for water resource and hazard management and for the scientific study of glacier hydrology, dynamics and hydrochemistry. In this paper, a distributed, physically based model is used to determine the effects of the up‐glacier retreat of the snowline on spatial and temporal patterns of melt and water routing across a small (0·11 km²) supraglacial catchment on Haut Glacier d'Arolla, Switzerland. The melt model uses energy‐balance theory and accounts for the effects of slope angle, slope aspect and shading on the net radiation fluxes, and the effects of atmospheric stability on the turbulent fluxes. The water routing model uses simplified snow and open‐channel hydrology theory and accounts for the delaying effects of vertical and horizontal water flow through snow and across ice. The performance of the melt model is tested against hourly measurements of ablation in the catchment. Calculated and measured ablation rates show a high correlation (r² = 0·74) but some minor systematic discrepancies in the short term (hours). These probably result from the freezing of surface water at night, the melting of the frozen layer in the morning, and subsurface melting during the afternoon. The performance of the coupled melt/routing model is tested against hourly discharge variations measured in the supraglacial stream at the catchment outlet. Calculated and measured runoff variations show a high correlation (r² = 0·62). Five periods of anomalously high measured discharge that were not predicted by the model were associated with moulin overflow events. The radiation and turbulent fluxes contribute c. 86% and c. 14% of the total melt energy respectively. These proportions do not change significantly as the surface turns from snow to ice, because increases in the outgoing shortwave radiation flux (owing to lower albedo) happen to be accompanied by decreases in the incoming shortwave radiation flux (owing to lower solar incidence angles) and increases in the turbulent fluxes (owing to higher air temperatures and vapour pressures). Model sensitivity experiments reveal that the net effect of snow pack removal is to increase daily mean discharges by c. 50%, increase daily maximum discharges by >300%, decrease daily minimum discharges by c. 100%, increase daily discharge amplitudes by >1000%, and decrease the lag between peak melt rates and peak discharges from c. 3 h to c. 50 min. These changes have important implications for the development of subglacial drainage systems. Copyright © 2002 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.     

Seasonal sediment fluxes forcing supraglacial melting on the Wright Lower Glacier,McMurdo Dry Valleys,Antarctica

Cold‐based glaciers exist in low temperature and low humidity environments in which shortwave radiation is the largest source of energy to the glacier surface and the energy budget is very sensitive to the surface albedo. Consequently, the presence of relatively low volumes of debris on glacier surfaces has a significant impact on the timing, magnitude and rate of ablation at the surface. The aim of this study is to understand how the presence of sediment on the glacier surface at the start of the melt season can affect meltwater generation and delivery on a cold‐based glacier. A combination of field measurements, energy balance modelling and chemical mixing modelling were used on the Wright Lower Glacier, McMurdo Dry Valleys, Antarctica, between October 2005 and January 2006 to address this aim. In this system, sediment was transported onto the glacier surface during the winter months (March–October) by foehn winds, which reduced surface albedo at the start of the summer melt season. The areas of the glacier on which sediment accumulated began to melt earlier than other parts of the glacier and experienced a longer melt season. Over the study period, the total ablation on the dirty surfaces was nine times greater than for clean ice. Ablation on the dirty surfaces is dominated by melting, whereas sublimation dominates the clean ice. As the sediment was unevenly distributed over the glacier surface, the variation in melt amount and timing drove the development of a cryoconite hole system. Copyright © 2012 John Wiley & Sons, Ltd.     

Subglacial sediment production and snout marginal ice uplift during the late ablation season of a temperate valley glacier

Sediment export from glaciated basins involves complex interactions between ice flow, basal erosion and sediment transfer in subglacial and proglacial streams. In particular, we know very little about the processes associated with sediment transfer by subglacial streams. The Haut Glacier d'Arolla (VS, Switzerland) was investigated during the summer melt season of 2015. LiDAR survey revealed positive surface changes in the ablation zone, indicating glacier uplift, at the end of the morning during the period of peak ablation. Instream measures of sediment transport showed that suspended load and bedload responded differently to diurnal flow variability. Suspended load depended on the availability of fine material whereas bedload depended mainly on the competence of the flow. Interpretation of these results allowed development of a conceptual model of subglacial sediment transport dynamics. It is based upon the mechanisms of clogging (deposition) and flushing (transport/erosion) in sub-glacial channels as forced by diurnal flow variability. Through the melt season, the glacier hydrological response evolves from being buffered by glacier snow cover with a poorly developed subglacial drainage system to being dominated by more rapid ice melt with a more hydraulically efficient subglacial channel system. The resultant changes in the shape of diurnal discharge hydrographs, and notably higher peak flows and lower base flows, causes sediment transport to become discontinuous, with overnight clogging and late morning flushing of subglacial channels. Overnight clogging may be sufficient to reduce subglacial channel size, creating temporarily pressurized flow and lateral transfer of water away from the subglacial channels, leading to the late morning glacier surface uplift. However, without further data, we cannot exclude other hypotheses for the uplift. © 2018 John Wiley & Sons, Ltd.