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.     

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.     

Seasonal reorganization of subglacial drainage inferred from measurements in boreholes

The effect of the formation of a major subglacial drainage channel on the behaviour of the subglacial drainage system of Haut Glacier d'Arolla, Switzerland, was investigated using measurements of borehole water level and the electrical conductivity and turbidity of basal meltwaters. Electrical conductivity profiles were also measured within borehole water columns to identify the water sources driving water level changes, and to determine patterns of water circulation in boreholes. Prior to channel formation, boreholes showed idiosyncratic and poorly coordinated behaviour. Diurnal water level fluctuations were small and driven by supraglacial/englacial water inputs, even when boreholes were connected to a subglacial drainage system. This system appeared to consist of hydraulically impermeable patches interspersed with storage spaces, and transmitted a very low water flux. Drainage reorganization, which occurred around 31 July, 1993, in response to rapidly rising meltwater and rainfall inputs, seems to have involved the creation of a connection between an incipient channel and a well-established channelized system located further down-glacier. Once a major channel existed within the area of the borehole array, borehole water level fluctuations were forced by discharge-related changes in channel water pressure, although a diversity of responses was observed. These included (i) synchronous, (ii) damped and lagged, (iii) inverse, and (iv) alternating inverse/lagged responses. Synchronous responses occurred in boreholes connected directly to the channel, while damped and lagged responses occurred in boreholes connected to it by a more resistive drainage system. Pressure variations within the channel resulted in diurnal transfer of mechanical support for the ice overburden between connected and unconnected areas of the bed, producing inverse and alternating patterns of water level response. © 1998 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.     

Development and testing of a physically based,three-dimensional model of surface and subsurface hydrology

We present a numerical, catchment-scale model that solves flow equations of surface and subsurface flow in a three-dimensional domain. Surface flow is described by the two-dimensional parabolic approximation of the St. Venant equation, using Manning’s equation of motion; subsurface flow is described by the three-dimensional Richards’ equation for the unsaturated zone and by three-dimensional Darcy’s law for the saturated zone, using an integrated finite difference formulation. The hydrological component is a dynamic link library implemented within a comprehensive model which simulates surface energy, radiation budget, snow melt, potential evapotranspiration, plant development and plant water uptake. We tested the model by comparing distributed and integrated three-dimensional simulated and observed perched water depth (PWD), stream flow data, and soil water contents for a small catchment. Additional tests were performed for the snow melting algorithm as well as the different hydrological processes involved. The model successfully described the water balance and its components as evidenced by good agreement between measured and modelled data.     

Modelling Canadian prairie wetland hydrology using a semi‐distributed streamflow model

A hydrological model (SLURP) that was designed for simulating hydrological processes taking place in large river basins was, with minimal modification, used successfully to simulate water level variations over a 28‐year period (1969–1996) for a 3‐ha prairie wetland in Saskatchewan. The model calculates a water balance based on precipitation, snowmelt, evaporation, surface runoff and subsurface flow on a daily time‐step. The model was first calibrated for two periods (1969–1973 for cropland and 1987–1990 for grassland), then it was applied to records outside the calibration periods. The model reproduced the wetland water level variations during a 28‐year period with good accuracy. The wetland water levels were most sensitive to the infiltration coefficient of surface soil under frozen conditions and to maximum soil moisture storage. The applicability of the model and the calibrated parameters to a smaller wetland, with an area of 0·24 ha, was examined. This simulation indicated that scale effects are important, probably largely in relation to snow redistribution by wind. Copyright © 2000 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.     

HYDROCHEMISTRY AS AN INDICATOR OF SUBGLACIAL DRAINAGE SYSTEM STRUCTURE: A COMPARISON OF ALPINE AND SUB-POLAR ENVIRONMENTS

The anion compositions (SO²₄, HCO⁻₃ and Cl⁻) of runoff from the Haut Glacier d'Arolla, Switzerland and Austre Brøggerbreen, Svalbard are compared to assess whether or not variations in water chemistry with discharge are consistent with current understanding of the subglacial drainage structure of warm- and polythermal-based glaciers. These glacial catchments have very different bedrocks and the subglacial drainage structures are also believed to be different, yet the range of anion concentrations show considerable overlap for SO²⁻₄ and HCO⁻₃. Concentrations of Cl⁻ are higher at Austre Brøggerbreen because of the maritime location of the glacier. Correcting SO²⁻₄ for the snowpack component reveals that the variation in non-snowpack SO²⁻₄ with discharge and with HCO⁻₃ is similar to that observed at the Haut Glacier d'Arolla. Hence, if we assume that the provenance of the non-snowpack SO²⁻₄ is the same in both glacial drainage systems, a distributed drainage system also contributes to runoff at Austre Brøggerbreen. We have no independent means of testing the assumption at present. The lower concentrations of non-snowpack SO²⁻₄ at Austre Brøggerbreen may suggest that a smaller proportion of runoff originates from a distributed drainage system than at the Haut Glacier d'Arolla.     

VELOCITY–DISCHARGE RELATIONSHIPS DERIVED FROM DYE TRACER EXPERIMENTS IN GLACIAL MELTWATERS: IMPLICATIONS FOR SUBGLACIAL FLOW CONDITIONS

Repeated dye tracer tests were undertaken from individual moulins at Haut Glacier d'Arolla, Switzerland, over a number of diurnal discharge cycles during the summers of 1989–1991. It was hoped to use the concepts of at-a-station hydraulic geometry to infer flow conditions in subglacial channels from the form of the velocity–discharge relationships derived from these tests. The results, however, displayed both clockwise and anticlockwise velocity–discharge hysteresis, in addition to the simple power function relationship assumed in the hydraulic geometry approach. Clockwise hysteresis seems to indicate that a moulin drains into a small tributary channel rather than directly into an arterial channel, and that discharges in the two channels vary out of phase with each other. Anticlockwise hysteresis is accompanied by strong diurnal variations in the value of dispersivity derived from the dye breakthrough curve, and is best explained by hydraulic damming of moulins or sub/englacial passageways. Despite the complex velocity–discharge relationships observed, some indication of subglacial flow conditions may be obtained if tributary channels comprise only a small fraction of the drainage path and power function velocity–discharge relationships are derived from dye injections conducted during periods when the supraglacial discharge entering the moulin and the bulk discharge vary in phase. Analyses based on this premise suggest that both open and closed channel flow occur beneath Haut Glacier d'Arolla, and that flow conditions are highly variable at and between sites.     

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

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

Streamflow response to the rapid retreat of a lake‐calving glacier

There has been increasing attention over the last decade to the potential effects of glacier retreat on downstream discharge and aquatic habitat. This study focused on streamflow variability downstream of Bridge Glacier in the southern Coast Mountains of BC between 1979 and 2014, prior to and during a period in which the glacier experienced enhanced calving and rapid retreat across a lake‐filled basin. Here we combined empirical trend detection and a conceptual‐parametric hydrological model to address the following hypotheses: (1) streamflow trends in late summer and early autumn should reflect the opposing influences of climatic warming (which would tend to increase unit‐area meltwater production) and the reduction in glacier area (which would tend to reduce the total volume of meltwater generated), and (2) winter streamflow should increase because of displacement of lake water as ice flows past the grounding line and calves into the lake basin. In relation to the first hypothesis, we found no significant trends in monthly discharge during summer. However, applying regression analysis to account for air temperature and precipitation variations, weak but statistically significant negative trends were detected for August and melt season discharge. The HBV‐EC model was applied using time‐varying glacier cover, as derived from Landsat imagery. Relative to simulations based on constant glacier extent, model results indicated that glacier recession caused a decline in mean monthly streamflow of 9% in August and 11% in September. These declines in late‐summer streamflow are consistent with the results from our empirical analysis. The second hypothesis is supported by the finding of positive trends for December, January, and February discharge. Despite the modelled declines in late‐summer mean monthly streamflow, recorded discharge data exhibited neither positive nor negative trends during the melt season, suggesting that Bridge Glacier may currently be at or close to the point of peak water. Further analysis of the impact of lake‐terminating glaciers on downstream discharge is needed to refine the peak water model. Copyright © 2016 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.     

Surface melt-driven seasonal behaviour (englacial and subglacial) from a soft-bedded temperate glacier recorded by in situ wireless probes

We investigate the spatial and temporal englacial and subglacial processes associated with a temperate glacier resting on a deformable bed using the unique Glacsweb wireless in situ probes (embedded in the ice and the till) combined with other techniques [including ground penetrating radar (GPR) and borehole analysis]. During the melt season (spring, summer and autumn), high surface melt leads to high water pressures in the englacial and subglacial environment. Winter is characterized by no surface melting on most days (‘base’) apart from a series of positive degree days. Once winter begins, a diurnal water pressure cycle is established in the ice and at the ice/sediment interface, with direct meltwater inputs from the positive degree days and a secondary slower englacial pathway with a five day lag. This direct surface melt also drives water pressure changes in the till. Till deformation occurred throughout the year, with the winter rate approximately 60% that of the melt season. We were able to show the bed comprised patches of till with different strengths, and were able to estimate their size, relative percentage and temporal stability. We show that the melt season is characterized by a high pressure distributed system, and winter by a low pressure channelized system. We contrast this with studies from Greenland (overlying rigid bedrock), where the opposite was found. We argue our results are typical of soft bedded glaciers with low englacial water content, and suggest this type of glacier can rapidly respond to surface-driven melt. Based on theoretical and field results we suggest that the subglacial hydrology comprises a melt season distributed system dominated by wide anastomosing broad flat channels and thin water sheets, which may become more channelized in winter, and more responsive to changes in meltwater inputs. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Mixing law versus discharge and electrical conductivity relationships: application to an alpine proglacial stream

Mixing laws have been used in hydrogeology for decades. In glacial hydrology, they are used to determine the contributions of the different reservoirs supplying the proglacial stream. However, some assumptions about discharge–solute matter and discharge–electrical conductivity (EC) relationships have led to erroneous conclusions. Analysis of a theoretical example indicated a hyperbolic relationship between discharge and solute matter concentration/EC for two reservoirs. We applied this approach to experimental data from the Baounet Glacier (Savoie, France), assuming that the proglacial stream is provided by a combination of watershed groundwater flow (diurnally stable and with high EC) and quick drainage from the glacier (low EC and discharge that change very quickly). Applying mixing laws allowed us to characterise these two flows and to separate the quick flows (glacial ablation and rainfalls). Copyright © 2011 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.     

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.     

Supraglacial weathering crust dynamics inferred from cryoconite hole hydrology

Water levels in cryoconite holes were monitored at high resolution over a 3‐week period on Austre Brøggerbreen (Svalbard). These data were combined with melt and energy balance modelling, providing insights into the evolution of the glacier's near‐surface hydrology and confirming that the hydrology of the near‐surface, porous ice known as the ‘weathering crust’ is dynamic and analogous to a shallow‐perched aquifer. A positive correlation between radiative forcing of melt and drainage efficiency was found within the weathering crust. This likely resulted from diurnal contraction and dilation of interstitial pore spaces driven by variations in radiative and turbulent fluxes in the surface energy balance, occasionally causing ‘sudden drainage events’. A linear decrease in water levels in cryoconite holes was also observed and attributed to cumulative increases in near‐surface ice porosity over the measurement period. The transport of particulate matter and microbes between cryoconite holes through the porous weathering crust is shown to be dependent upon weathering crust hydraulics and particle size. Cryoconite holes therefore yield an indication of the hydrological dynamics of the weathering crust and provide long‐term storage loci for cryoconite at the glacier surface. This study highlights the importance of the weathering crust as a crucial component of the hydrology, ecology and biogeochemistry of the glacier ecosystem and glacierized regions and demonstrates the utility of cryoconite holes as natural piezometers on glacier surfaces. 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.