Structure,morphology and water flux of a subglacial drainage system,Midtdalsbreen, Norway

Digital elevation models of the surface and bed of Midtdalsbreen, Norway are used to calculate subglacial hydraulic potential and infer drainage system structure for a series of subglacial water pressure assumptions ranging from atmospheric to ice overburden. A distributed degree‐day model is used to calculate the spatial distribution of melt on the glacier surface throughout a typical summer, which is accumulated along the various drainage system structures to calculate water fluxes beneath the glacier and exiting the portals for the different water pressure assumptions. In addition, 78 dye‐tracing tests were performed from 33 injection sites and numerous measurements of water discharge were made on the main proglacial streams over several summer melt seasons. Comparison of the calculated drainage system structures and water fluxes with dye tracing results and measured proglacial stream discharges suggests that the temporally and spatially averaged steady‐state water pressures beneath the glacier are ~70% of ice overburden. Analysis of the dye return curves, together with the calculated subglacial water fluxes shows that the main drainage network on the eastern half of the glacier consists of a hydraulically efficient system of broad, low channels (average width/height ratio ≈ 75). The smaller drainage network on the west consists of a hydraulically inefficient distributed system, dominated by channels that are exceptionally broad and very low (average width/height ratio ≈ 350). The even smaller central drainage network also consists of a hydraulically inefficient distributed system, dominated by channels that are very broad and exceptionally low (average width/height ratio ≈ 450). The channels beneath the western and central glacier must be so broad and low that they can essentially be thought of as a linked cavity system. Copyright © 2011 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.     

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.     

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.     

The hydrochemistry of runoff from a ‘cold-based’ glacier in the High Arctic (Scott Turnerbreen,Svalbard)

There are still relatively few hydrochemical studies of glacial runoff and meltwater routing from the high latitudes, where non-temperate glacier ice is frequently encountered. Representative samples of glacier meltwater were obtained from Scott Turnerbreen, a ‘cold-based’ glacier at 78° N in the Norwegian high Arctic archipelago of Svalbard, during the 1993 melt season and analysed for major ion chemistry. Laboratory dissolution experiments were also conducted, using suspended sediment from the runoff. Significant concentrations of crustal weathering derived SO²⁻₄ are present in the runoff, which is characterized by high ratios of SO²⁻₄: (SO²⁻₄+HCO⁻₃) and high p(CO₂). Meltwater is not routed subglacially, but flows to the glacier terminus through subaerial, ice marginal channels, and partly flows through a proglacial icing, containing highly concentrated interstitial waters, immediately afront the terminus. The hydrochemistry of the runoff is controlled by: (1) seasonal variations in the input of solutes from snow- and icemelt; (2) proglacial solute acquisition from the icing; and (3) subaerial chemical weathering within saturated, ice-cored lateral moraine adjoining drainage channels at the glacier margins, sediment and concentrated pore water from which is entrained by flowing meltwater. Diurnal variations in solute concentration arise from the net effects of variable sediment pore water entrainment and dilution in the ice marginal streams. Explanation of the hydrochemistry of Scott Turnerbreen requires only one major subaerial flow path, the ice marginal channel system, in which seasonally varying inputs of concentrated snowmelt and dilute icemelt are modified by seepage or entrainment of concentrated pore waters from sediment in lateral moraine, and by concentrated interstitial waters from the proglacial icing, supplied by leaching, slow drainage at grain intersections or simple melting of the icing itself. The ice marginal channels are analogous neither to dilute supra/englacial nor to concentrated subglacial flow components. © 1998 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.     

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.     

Till deposition by glacier submarginal,incremental thickening

Macro‐ and micro‐scale sedimentological analyses of recently deposited tills and complex push/squeeze moraines on the forelands of Icelandic glaciers and in a stacked till sequence at the former Younger Dryas margin of the Loch Lomond glacier lobe in Scotland are used to assess the depositional processes involved in glacier submarginal emplacement of sediment. Where subglacial meltwater is unable to flush out subglacial sediment or construct thick debris‐rich basal ice by cumulative freeze‐on processes, glacier submarginal processes are dictated by seasonal cycles of refreezing and melt‐out of tills advected from up‐ice by a combination of lodgement, deformation and ice keel and clast ploughing. Although individual till layers may display typical A and B horizon deformation characteristics, the spatially and temporally variable mosaic of subglacial processes will overprint sedimentary and structural signatures on till sequences to the extent that they would be almost impossible to classify genetically in the ancient sediment record. At the macro‐scale, Icelandic tills display moderately strong clast fabrics that conform to the ice flow directions documented by surface flutings; very strong fabrics typify unequivocally lodged clasts. Despite previous interpretations of these tills as subglacial deforming layers, micro‐morphological analysis reveals that shearing played only a partial role in the emplacement of till matrixes, and water escape and sediment flowage features are widespread. A model of submarginal incremental thickening is presented as an explanation of these data, involving till slab emplacement over several seasonal cycles. Each cycle involves: (1) late summer subglacial lodgement, bedrock and sediment plucking, subglacial deformation and ice keel ploughing; (2) early winter freeze‐on of subglacial sediment to the thin outer snout; (3) late winter readvance and failure along a decollement plane within the till, resulting in the carriage of till onto the proximal side of the previous year's push moraine; (4) early summer melt‐out of the till slab, initiating porewater migration, water escape and sediment flow and extrusion. Repeated reworking of the thin end of submarginal till wedges produces overprinted strain signatures and clast pavements. Copyright © 2005 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.     

Flow separation on Zongo Glacier,Cordillera Real,Bolivia

Meltwaters collected from the proglacial stream escaping from Zongo Glacier (2·1 km²), Bolivia (16°S), have been monitored in order to analyse the internal drainage system of an Andean glacier. Electrical conductivity has been measured sporadically between February 1995 and March 1996, during 16 one-day field surveys, under various meteorological conditions in summer and winter. The mixing-model technique based on the electrical conductivity is used for a quantitative separation of discharge which is derived from continuous water level registration. Tracer experiments (mainly uranine dye and NaCl salt) have been carried out from March to June 1997 to obtain information about the internal drainage system. In the tropical Andes, accumulation only occurs in austral summer, whereas ablation occurs throughout the year and is higher during the accumulation season, between November and March. The assumptions involved in the use of mixing models for analysis of glacial drainage structure are applicable for tropical glaciers because glacial conduits do not suffer complete closure, and are permanently supplied by meltwaters, even in wintertime. Two components of discharge are separated: an englacial flow originating from surface meltwater which is routed without chemical enrichment, and offering low electrical conductivity; and a subglacial one routed in contact with bedrock or sediments showing high ionic concentrations. Electrical conductivity of meltwater varies diurnally, inversely to discharge fluctuations. According to this behaviour, total discharge is mainly formed by the englacial component. The drainage structures for englacial and subglacial flow have to be widely interconnected, as indicated by diurnal variations of the subglacial discharge. Comparison of hydrograph separation based on conductivity and on ¹⁸O isotope confirms that the subglacial flow is influenced by surface melting. A hydrograph separation of the subglacial flow is proposed, between a diurnal variable component, composed of water coming from the englacial network, and a base flow, which may vary seasonally. The dye tracing experiments confirm the drainage complexity of Zongo Glacier and demonstrate the interest of identifying three main drainage components. © 1998 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.     

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.     

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.     

VARIABILITY IN THE CHEMICAL COMPOSITION OF IN SITU SUBGLACIAL MELTWATERS

Meltwaters collected from boreholes drilled to the base of the Haut Glacier d'Arolla, Switzerland have chemical compositions that can be classified into three main groups. The first group is dilute, whereas the second group is similar to, though generally less concentrated in major ions, than contemporaneous bulk glacial runoff. The third group is more concentrated than any observed bulk runoff, including periods of flow recession. Waters of the first group are believed to represent supraglacial meltwater and ice melted during drilling. Limited solutes may be derived from interactions with debris in the borehole. The spatial pattern of borehole water levels and borehole water column stratification, combined with the chemical composition of the different groups, suggest that the second group represent samples of subglacial waters that exchange with channel water on a diurnal basis, and that the third group represent samples of water draining through a ‘distributed’ subglacial hydraulic system. High NO⁻₃ concentrations in the third group suggest that snowmelt may provide a significant proportion of the waters and that the residence time of the waters at the bed in this particular section of the distributed system is of the order of a few months. The high NO⁻₃ concentrations also suggest that some snowmelt is routed along different subglacial flowpaths to those used by icemelt. The average SO²⁻₄: (HCO⁻₃ + SO²⁻₄) ratio of the third group of meltwaters is 0.3, suggesting that sulphide oxidation and carbonate dissolution (which gives rise to a ratio of 0.5) cannot provide all the HCO⁻₃ to solution. Hence, carbonate hydrolysis may be occurring before sulphide oxidation, or there may be subglacial sources of CO₂, perhaps arising from microbial oxidation of organic C in bedrock, air bubbles in glacier ice or pockets of air trapped in subglacial cavities. The channel marginal zone is identified as an area that may influence the composition of bulk meltwater during periods of recession flow and low diurnal discharge regimes. © 1997 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.     

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.     

Meltwater discharge through the subglacial bed and its land‐forming consequences from numerical experiments in the Polish lowland during the last glaciation

Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 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.     

Subglacial observations from øksfjordjøkelen,North Norway

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

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.