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.     

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.     

Borehole drainage and its implications for the investigation of glacier hydrology: experiences from Haut Glacier d'Arolla,Switzerland

Studies of glacier hydrology rely increasingly on measurements made in boreholes as a basis for reconstructing the character and behaviour of subglacial drainage systems. In temperate glaciers, in which boreholes remain open to the atmosphere following drilling, the interpretation of such data may be complicated by supraglacial or englacial water flows to and from boreholes. We report on a suite of techniques used to identify borehole water sources and to reconstruct patterns of water circulation within boreholes at Haut Glacier d'Arolla, Switzerland. Results are used to define a number of borehole ‘drainage’ types. Examples of each drainage type are presented, along with the manner in which they influence interpretations of borehole water‐levels, borehole water‐quality data, and borehole dye traces. The analysis indicates that a full understanding of possible borehole drainage modes is required for the correct interpretation of many borehole observations, and that those observations provide an accurate indication of subglacial conditions only under relatively restricted circumstances. Copyright © 2001 John Wiley & Sons, Ltd.     

Geometry,bed topography and drainage system structure of the haut glacier d'Arolla,Switzerland

As part of an integrated study of the hydrology, meltwater quality and dynamics of the Haut Glacier d'Arolla, Switzerland, the glacier's drainage network structure was determined from patterns of dye recovery in 342 injection experiments conducted from 47 moulins distributed widely across the glacier. This structure was compared with theoretical predictions based upon reconstructed patterns of water flow governed by (a) the subglacial hydraulic potential surface, and (b) the subglacial bedrock surface. These reconstructions were based on measurements of ice surface and bedrock topography obtained by a combination of ground survey and radio-echo sounding techniques. The two reconstructions simulate the drainage system structures expected for (a) closed channels, in which water is pressurized by the overlying ice, and (b) gravity-driven, open-channel flow. The closed-channel model provides the best fit to the observed structure, even though theoretical calculations suggest that, under summer discharge conditions, open-channel flow may be widespread beneath the glacier. Possible reasons for this apparent discrepancy are discussed.     

Sulphide oxidation under partially anoxic conditions at the bed of the Haut Glacier d'Arolla,Switzerland

Oxygen isotopic data are presented for bulk glacial meltwaters draining the Haut Glacier d'Arolla, Valais, Switzerland and for the sulphate contained within them in an attempt to assess the redox status of the subglacial chemical weathering environment. The sulphate derived from subglacial chemical weathering is so depleted in ¹⁸O that it must have formed, at least partially, in an anoxic environment. Under these conditions, Fe³⁺ can act as an oxidizing agent and oxygen atoms incorporated into sulphate are derived from ¹⁸O‐depleted water molecules (by contrast, dissolved O₂ is strongly enriched in ¹⁸O). These data therefore support the hypothesis that sectors of the glacier bed are anoxic and that Fe³⁺ may act as a significant oxidizing agent under these conditions. Copyright © 2002 John Wiley & Sons, Ltd.     

Recent advances in understanding ice sheet dynamics

Glaciers and ice sheets play a dynamic role in Earth's climate system, influencing regional- and global-scale climate and responding to climate change on time scales from years to millennia. They are also an integral part of Earth's landscape in alpine and polar regions, where they are an active agent in isostatic, tectonic, and Earth surface processes. This review paper summarizes recent progress in understanding and modelling ice sheet dynamics, from the microphysical processes of ice deformation in glaciers to continental-scale processes that influence ice dynamics. Based on recent insights and research directions, it can be expected that a new generation of ice sheet models will soon replace the current standard. Improvements that can be foreseen in the near future include: (i) the addition of internally-consistent evolutionary equations for ice crystal fabric (anisotropic flow laws), (ii) more generalized flow laws that include different deformation mechanisms under different stress regimes, (iii) explicit incorporation of the effects of chemical impurities and grain size (dynamic recrystallization) on ice deformation, (iv) higher-order stress solutions to the momentum balance (Stokes' equation) that governs ice sheet flow, and (v) the continued merger of ice sheet models with increasingly complex Earth systems models, which include fully-coupled subglacial hydrological and geological processes. Examples from the Greenland Ice Sheet and Vatnajökull Ice Cap, Iceland are used to illustrate several of these new directions and their importance to glacier dynamics.     

Gravity profiling as a technique for determining the thickness of glacier ice

Summary Detailed gravity measurements recently carried out on the Gorner glacier, Switzerland, are used to determine the variation of thickness across the glacier ice. The Gorner glacier was chosen as a test site because seismic control was available. The glacier ice at a profile near the Monte Rosa massif is associated with a relative gravity low of about –23 mgal. Model oalculations yield a corresponding ice thickness of about 400 m at the central part of the profile. A comparison of the derived residual gravity anomaly with the calculated effect of the 3-D ice model based on seismic information is made. It is shown that the regional field determined for the Gorner glacier is appropriate and gives the correct residual anomaly associated with the glacier ice. Therefore, the proposed gravity technique for determining variations of the thickness of glacier ice appears to be a valuable and rather inexpensive method for surveying glaciers.Institut für Geophysik, ETH-Zürich, Contribution No 145.     

A radar transparent layer in a temperate valley glacier: Bench Glacier,Alaska

Radar surveys of Bench Glacier, Alaska, collected over five field seasons between 2002 and 2006 reveal a surface layer of radar transparent ice in this temperate valley glacier. The transparent layer covers the up‐glacier half of the ablation zone and is defined by a distinct lack of the radar scattering events considered typical of temperate ice. Radar scattering ice underlies the transparent zone, and extends to the surface elsewhere on the glacier. We observed the layering in constant offset radar surveys conducted with characteristic frequencies ranging from 5 MHz to 100 MHz. The radar transparent layer extends from the surface to 20 m depth on average, but up to 50 m in some places. Bench Glacier's transparent layer appears similar to the cold surface layer of polythermal glaciers, however, observations in over 50 boreholes on Bench Glacier suggest there is no cold ice corresponding to the radar transparent layer. We conclude that spatially extensive radar‐transparent layers normally used to identify cold ice in polythermal glaciers are present in some temperate glaciers. Copyright © 2009 John Wiley & Sons, Ltd.     

The hydrochemistry of meltwaters draining a polythermal-based,high Arctic glacier,south Svalbard: I. The ablation season

Solute and runoff time-series at Finsterwalderbreen, Svalbard, provide evidence for considerable basal routing of water and the existence of at least two contrasting subglacial chemical weathering environments. The hydrochemistry of a subglacial upwelling provides evidence for a snowmelt-fed subglacial reservoir that dominates bulk runoff during recession flow. High concentrations of Cl⁻ and crustal ions, high pCO₂ and ratios of [*SO²⁻₄/(*SO²⁻₄+HCO⁻₃)] close to 0·5 indicate the passage of snowmelt through a subglacial weathering environment characterized by high rock:water ratios, prolonged residence times and restricted access to the atmosphere. At higher discharges, bulk runoff becomes dominated by icemelt from the lower part of the glacier that is conveyed through a chemical weathering environment characterized by low rock:water ratios, short residence times and free contact with atmospheric gases. These observations suggest that icemelt is routed via a hydrological system composed of basal/ice-marginal, englacial and supraglacial components and is directed to the glacier margins by the ice surface slope. Upwelling water flows relatively independently of icemelt to the terminus via a subglacial drainage system, possibly constituting flow through a sediment layer. Cold basal ice at the terminus forces it to take a subterranean routing in its latter stages. The existence of spatially discrete flow paths conveying icemelt and subglacial snowmelt to the terminus may be the norm for polythermal-based glaciers on Svalbard. Proglacial mixing of these components to form the bulk meltwaters gives rise to hydrochemical trends that resemble those of warm-based glaciers. These hydrochemical characteristics of bulk runoff have not been documented on any other glacier on Svalbard to date and have significance for understanding interactions between thermal regime and glacier hydrology. © 1998 John Wiley & Sons, Ltd.     

Melt water driven stream and groundwater stage fluctuations on a glacier forefield (Dammagletscher,Switzerland)

In many mountain regions, large land areas with heterogeneous soils have become ice‐free with the ongoing glacier retreat. On these recently formed proglacial fields, the melt of the remaining glaciers typically drives pronounced diurnal stream level fluctuations that propagate into the riparian zone. This behaviour was measured on the Damma glacier forefield in central Switzerland with stage recorders in the stream and groundwater monitoring wells along four transects. In spite of the large groundwater stage variations, radon measurements in the near‐stream riparian zone indicate that there is little mixing between stream water and groundwater on daily time scales. At all four transects, including both losing and gaining reaches, the groundwater level fluctuations lagged the stream stage variations and were often damped with distance from the stream. Similar behaviours have been modelled using the diffusion equation in coastal regions influenced by tidal sea level variations. We thus tested the ability of such a model to predict groundwater level fluctuations in proglacial fields. The model reproduced several key features of the observed fluctuations at three of four locations, although discrepancies also arise due to non representative input data and model simplifications. Nevertheless, calibration of the model for the individual transects yielded realistic estimates of hydraulic diffusivities between the stream and groundwater monitoring wells. We conclude that studying diurnal groundwater fluctuations can provide important information about the subsurface hydrology of alpine watersheds dominated by glacier melt. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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.     

Quantifying ice storage in upper Indus river basin using ground-penetrating radar measurements and glacier bed topography model version 2

Ice reserve estimates is a fundamental prerequisite for water resources management. We selected the UIB (upper Indus river basin) as study area because it contains the most abundant mid-latitude glaciers outside the polar region, however, the ice reserve estimates remain unclear due to the harsh topography. In this study, we validated the parameters of the GlabTop2 model (Glacier Bed Topography version 2) using the ice thickness measured by GPR (ground-penetrating radar) and compared the “GPR-measured ice thickness and ice bed elevation” versus “the estimated results obtained from GlabTop2.” Integrated with IDW (inverse distance weighted) interpolated results of glaciers of various sizes, a reasonable parametric scheme (τ = 120 kPa and f = 0.8) of GlabTop2 was applied on vast amounts of glaciers in the UIB region. The GlabTop2 estimates indicated that the ice thickness of the UIB varied from 0 to 736.0 ± 110.0 m, with an average value of 74.5 ± 11.2 m. A significant spatial heterogeneity exists in the sub-basins. The Shyok, Shigar and Hunza that contain the most abundant ice reserve. Lesser quantities are stored in the Western Himalaya and Hindu Kush ranges, which account for 11.3 % and 6.7 % of the total ice reserve in the UIB, respectively. A total volume of 1162.4 ± 175.1 km³ of glacier can be converted to 1046.2 ± 157.6 Gt ice reserve; this is 13.6 times the annual average discharge obtained from the outlet of the Besham hydrological station. We aim to present estimates that can provide the baseline information for glaciology study of the Indus river.     

Structural composition and sediment transfer in a composite cirque glacier: Glacier de St. Sorlin,France

This paper considers the links between structure, sediment transport and sediment delivery at Glacier de St. Sorlin, France. Sediment transported by the glacier is concentrated at flow‐unit boundaries as medial moraines, controlled by the position of bedrock outcrops in the accumulation area. Rockfall entrained within primary stratification is tightly folded at flow‐unit boundaries under high cumulative strains and laterally compressive stress. High cumulative strains and laterally compressive stresses lead to the development of longitudinal foliation from primary stratification. Folding elevates subglacial sediments into foliation‐parallel debris ridges, which are exposed in the ablation area. Crevasses and shear planes within the glacier have little control on sediment transport. Debris stripes in the proglacial area are morphologically similar to foliation‐parallel debris ridges; however, they are not structurally controlled, but formed by fluvial erosion. The conclusion of this study is that at Glacier de St. Sorlin proglacial sediment‐landform associations are subjected to intense syn‐ and post‐depositional modification by high melt‐water discharges, hence their composition does not reflect that of sediments melting out at the terminus. The action of melt water limits the potential of the sedimentary record to be used to constrain numerical models of past glacier dynamics in debris‐poor glacierized Alpine catchments. Copyright © 2008 John Wiley & Sons, Ltd.     

Debris entrainment by basal freeze‐on and thrusting during the 1995–1998 surge of Kuannersuit Glacier on Disko Island,west Greenland

Kuannersuit Glacier, a valley glacier on Disko Island in west Greenland, experienced a major surge from 1995 to 1998 where the glacier advanced 10·5 km and produced a ～65 m thick stacked sequence of debris‐rich basal ice and meteoric glacier ice. The aim of this study is to describe the tectonic evolution of large englacial thrusts and the processes of basal ice formation using a multiproxy approach including structural glaciology, stable isotope composition (δ¹⁸O and δD), sedimentology and ground‐penetrating radar. We argue that the major debris layers that can be traced in the terminal zone represent englacial thrusts that were formed early during the surge. Thrust overthrow was at least 200–300 m and this lead to a 30 m thick repetition of basal ice at the ice margin. It is assumed that the englacial thrusting was initiated at the transition between warm ice from the interior and the cold snout. The basal debris‐rich ice was mainly formed after the thrusting phase. Two sub‐facies of stratified basal ice have been identified; a lower massive ice facies (S_M) composed of frozen diamict enriched with heavy stable isotopes overlain by laminated ice facies (S_L) consisting of millimetre thick lamina of alternating debris‐poor and debris‐rich ice. We interpret the stratified basal ice as a continuum formed mainly by freeze‐on processes and localized regelation. First laminated basal ice is formed and as meltwater is depleted more sediment is entrained and finally the glacier freezes to the base and massive diamict is frozen‐on. The increased ability to entrain sediments may partly be associated with higher basal freezing rates enhanced by loss of frictional heat from cessation of fast flow and conductive cooling through a thin heavily crevassed ice during the final phase of the glacier surge. The dispersed basal ice facies (D) was mainly formed by secondary processes where fine‐grained sediment is mobilized in the vein system of ice. Our results have important implications for understanding the significance of basal ice formation and englacial thrusting beneath fast‐flowing glaciers and it provides new information about the development of landforms during a glacier surge. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Retreating ice: research in pro‐glacial areas matters

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

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.     

Stream temperature response to variable glacier coverage in coastal watersheds of Southeast Alaska

We measured stream temperature continuously during the 2011 summer run‐off season (May through October) in nine watersheds of Southeast Alaska that provide spawning habitat for Pacific salmon. The nine watersheds have glacier coverage ranging from 0% to 63%. Our goal was to determine how air temperature and watershed land cover, particularly glacier coverage, influence stream temperature across the seasonal glacial meltwater hydrograph. Multiple linear regression models identified mean watershed elevation (related to glacier extent) and watershed lake coverage (%) as the strongest landscape controls on mean monthly stream temperature, with the weakest (May) and strongest (July) models explaining 86% and 97% of the temperature variability, respectively. Mean weekly stream temperature was significantly correlated with mean weekly air temperature in seven streams; however, the relationships were weak to non‐significant in the streams influenced by glacial run‐off. Streams with >30% glacier coverage showed decreasing stream temperatures with rising summer air temperatures, whereas those with <30% glacier coverage exhibited summertime warming. Glaciers also had a cooling effect on monthly mean stream temperature during the summer (July through September) equivalent to a decrease of 1.1 °C for each 10% increase in glacier coverage. The maximum weekly average temperature (an index of thermal suitability for salmon) in the six glacial streams was substantially below the lower threshold for optimum salmon growth. This finding suggests that although glaciers are important for moderating summer stream temperatures, future reductions in glacier run‐off may actually improve the thermal suitability of some glacially dominated streams in Southeast Alaska for salmon. Copyright © 2013 John Wiley & Sons, Ltd.