Search

Glacier basal motion is responsible for the majority of ice flux on fast-flowing glaciers, enables rapid changes in glacier motion and provides the means by which glaciers shape alpine landscapes. In an effort to enhance our understanding of basal motion, we investigate the evolution of glacier velocity and ice-marginal lake stage on Kennicott Glacier, Alaska, during the spring–summer transition, a time when subglacial drainage is undergoing rapid change. A complicated record of > 50 m fill-and-drain sequences on a hydraulically-connected ice-marginal lake likely reflects the punctuated establishment of efficient subglacial drainage as the melt season begins. The rate of change of lake stage generally correlates with diurnal velocity maxima, both in timing and magnitude. At the seasonal scale, the up-glacier progression of enhanced summer basal motion promotes uniformity of daily glacier velocity fluctuations throughout the 10 km study reach, and results in diurnal velocity patterns suggesting increasingly efficient meltwater delivery to and drainage from the subglacial channel system. Our findings suggest the potential of using an ice-marginal lake as a proxy for subglacial water pressure, and show how widespread basal motion affects bulk glacier behavior.

Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of glacial lakes connected by subglacial channels. Subglacial channel equations are supplied with effective pressure boundary conditions that are determined by a simple lake model. While the model can describe an arbitrary number of lakes, we solve it numerically with a finite element method for the case of two connected lakes. We examine the effect of relative lake size and spacing on the oscillations. Complex oscillations in the downstream lake are driven by discharge out of the upstream lake. These include multi-peaked and anti-phase filling–draining events. Similar filling–draining cycles have been observed on the Kennicott Glacier in Alaska and at the confluence of the Whillans and Mercer ice streams in West Antarctica. We further construct a simplified ordinary differential equation model that displays the same qualitative behavior as the full, spatially-dependent model. We analyze this model using dynamical systems theory to explain the appearance of filling–draining cycles as the meltwater supply varies.

Characterizing subglacial water flow is critical for understanding basal sliding and processes occurring under glaciers and ice sheets. Development of subglacial numerical models and acquisition of water pressure and tracer data have provided valuable insights into subglacial systems and their evolution. Despite these advances, numerical models, data conditioning, and uncertainty quantification are difficult, principally due to high number of unknown parameters and expensive forward computations. In this study, we aim to infer the properties of a subglacial drainage system in two dimensions using a framework that combines physical and geostatistical processes. The methodology is composed of three main components: (i) a channel generator to produce networks of the subglacial system, (ii) a physical model that computes water pressure and mass transport in steady state, and (iii) Bayesian inversion in which the outputs (pressure and tracer-transit times) are compared with synthetic data, thus allowing for parameter estimation and uncertainty quantification. We evaluate the ability of this framework to infer the subglacial characteristics of a synthetic ice sheet produced by a physically complex deterministic model, under different recharge scenarios. Results show that our methodology captures expected physical characteristics for each meltwater supply condition, while the precise locations of channels remain difficult to constrain. The framework enables uncertainty quantification, and the results highlight its potential to infer properties of real subglacial systems using observed water pressure and tracer-transit times.

Kuannersuit Glacier in west Greenland experienced a major surge in 1995–1998, where the glacier advanced 10.5 km down-valley. In this paper, we examine how the quiescent phase has progressed until 2015 with respect to ice dynamics, landform formations and glaciomorphological changes. In the initial quiescent phase (2001–2005), ice velocities along the center flowline were 37–55% higher than ice velocities (6.3–9.0 m yr −1 ) in 2005–2015, and a linear relationship between ice velocity and distance from the glacier front, which existed during the initial quiescent phase, had disappeared in 2005–2015. Between 2001 and 2015, the post-surge glacier tongue thinned by 84.0 ± 6.3 m, equal to 6.0 ± 0.5 m yr −1 . The 30-m high terminal moraine detached from the active glacier along the uppermost thrust band between 2005 and 2015, after which the glacier front receded 1.5 km. Between the terminal moraine and the glacier front, observations suggested that a proximal outwash plain was forming on top of glacier ice and that glacier naled (icing, aufeis) was incorporated into the outwash plain. In less dynamic areas with gently sloping topography along the lateral glacier margins, crevasse-squeeze ridges were melting out. The most prominent glaciomorphological feature in 2015 was a c. 1-km long gorge with c. 30–40 m high ice cliffs. This gorge had formed into the glacier front, following the collapse of a series of chasms above the main subglacial channel. Roof collapse processes above subglacial channels had a significant impact on the glacier's mass loss and on the glaciomorphology during the quiescent phase. Earlier observations from 1913 suggest that the formation of this gorge is a reoccurring phenomenon during the quiescent phase. A lower pothole field on the glacier tongue that existed during the initial quiescent phase had completely disappeared in 2015, while a larger upper pothole field in the reservoir area prevailed. During the initial quiescent phase, the subglacial drainage system developed into a stable channelized network, although a jökulhlaup from a large ice-dammed lake occurred. This single event was dated to have occurred between 11 and 13 August 2006, and the lake did not refill. The spatiotemporal association and evolution of glacial and glaciomorphological landforms at Kuannersuit Glacier are likely similar to processes at many other surging and rapidly receding glaciers of similar size.

Subglacial water drainage plays a significant role in glacier flow dynamics. Various forms of subglacial drainage have been observed beneath modern ice sheets, and the same forms of drainage likely occurred beneath paleo-ice sheets, such as the Laurentide Ice Sheet. Records of paleo drainage are occasionally preserved in the geomorphic record and can serve as a widespread and easily accessible means to investigate aspects of subglacial drainage that are difficult to directly study on modern-day glaciers. Linear surface depressions that extend tens of kilometers inward from the margin of Laurentide lobes are found throughout the Midwest of North America and are hypothesized to form from incision by channelized subglacial water flow. Here we estimate the subglacial hydropotential and its gradient along the western margin of the Marine Isotope Stage 2 Laurentide Ice Sheet's Green Bay Lobe to identify potential locations for subglacial water pooling. We find that linear surface depressions correlate well with areas where subglacial water likely pooled. We use a combination of active and passive seismic analyses to estimate the actual size of an incised subglacial channel in Waushara County, Wisconsin, that was subsequently infilled with sediment after the formation. We find that the channel incised ca. 65 m into the surrounding unlithified material, which is 6 times greater than the modern surface expression, and has a width of 450 m, which is nearly equal to the width of the surface expression. We estimate the channel would have taken ca. 60 days of water flow to form and that the flow of pore water from the surrounding till into the low-pressure channel could have increased the strength of the till by as much as 140 kPa following a single drainage by as much as 175 kPa in a smaller region following a series of 15-day-long drainage events. Till strengthening could have been a factor of 7 to 25× greater than predrainage strength values.

Nearly all meltwater from glaciers and ice sheets is routed englacially through moulins, which collectively comprise approximately 10–14 % of the efficient englacial–subglacial hydrologic system. Therefore, the geometry and evolution of moulins has the potential to influence subglacial water pressure variations, ice motion, and the runoff hydrograph delivered to the ocean. We develop the Moulin Shape (MouSh) model, a time-evolving model of moulin geometry. MouSh models ice deformation around a moulin using both viscous and elastic rheologies and models melting within the moulin through heat dissipation from turbulent water flow, both above and below the water line. We force MouSh with idealized and realistic surface melt inputs. Our results show that variations in surface melt change the geometry of a moulin by approximately 30 % daily and by over 100 % seasonally. These size variations cause observable differences in moulin water storage capacity, moulin water levels, and subglacial channel size compared to a static, cylindrical moulin. Our results suggest that moulins are significant storage reservoirs for meltwater, with storage capacity and water levels varying over multiple timescales. Representing moulin geometry within subglacial hydrologic models would therefore improve their accuracy, especially over seasonal periods or in regions where overburden pressures are high.

Subglacial meltwater channels (N-channels) are attributed to erosion by meltwater in subglacial conduits. They exert a major control on meltwater accumulation at the base of ice sheets, serving as drainage pathways and modifying ice flow rates. The study of exposed relict subglacial channels offers a unique opportunity to characterize the geomorphologic fingerprint of subglacial erosion as well as study the structure and characteristics of ice sheet drainage systems. In this study we present detailed field and remote sensing observations of exposed subglacial meltwater channels in excellent preservation state on Devon Island (Canadian Arctic Archipelago). We characterize channel cross section, longitudinal profiles, and network morphologies and establish the spatial extent and distinctive characteristics of subglacial drainage systems. We use field-based GPS measurements of subglacial channel longitudinal profiles, along with stereo imagery-derived digital surface models (DSMs), and novel kinematic portable lidar data to establish a detailed characterization of subglacial channels in our field study area, including their distinction from rivers and other meltwater drainage systems. Subglacial channels typically cluster in groups of ∼10 channels and are oriented perpendicular to active or former ice margins. Although their overall direction generally follows topographic gradients, channels can be oblique to topographic gradients and have undulating longitudinal profiles. We also observe that the width of first-order tributaries is 1 to 2 orders of magnitude larger than in Devon Island river systems and approximately constant. Furthermore, our findings are consistent with theoretical expectations drawn from analyses of flow driven by gradients in effective water pressure related to variations in ice thickness. Our field and remote sensing observations represent the first high-resolution study of the subglacial geomorphology of the high Arctic, and provide quantitative and qualitative descriptions of subglacial channels that revisit well-established field identification guidelines. Distinguishing subglacial channels in topographic data is critical for understanding the emergence, geometry, and extent of channelized meltwater systems and their role in ice sheet drainage. The final aim of this study is to facilitate the identification of subglacial channel networks throughout the globe by using remote sensing techniques, which will improve the detection of these systems and help to build understanding of the underlying mechanics of subglacial channelized drainage.

This research observed variation of flow velocity on the David Glacier, East Antarctica observed by using Sentinel-1A satellite. Because of the 12-day revisit cycle, relatively slow-flowing area was observed. Using double-differential InSAR (DDInSAR) technique, we could find the circular shaped velocity anomaly. We combined the ascending and descending DDInSAR images to figure out the direction of the true displacement over the anomaly. As a result, the circular anomaly appears mainly to be a vertical movement, not in the flow direction. This vertical movement in the David Glacier, where subglacial lakes are known to exist, can be thought of as a signal of water level change of a subglacial lake forming a subglacial channel system.

Water that pressurizes the base of glaciers and ice sheets enhances glacier velocities and modulates glacial erosion. Predicting ice flow and erosion therefore requires knowledge of subglacial channel evolution, which remains observationally limited. Here we demonstrate that detailed analysis of seismic ground motion caused by subglacial water flow at Mendenhall Glacier (Alaska) allows for continuous measurement of daily to subseasonal changes in basal water pressure gradient, channel size, and sediment transport. We observe intermittent subglacial water pressure gradient changes during the melt season, at odds with common assumptions of slowly varying, low-pressure channels. These observations indicate that changes in channel size do not keep pace with changes in discharge. This behavior strongly affects glacier dynamics and subglacial channel erosion at Mendenhall Glacier, where episodic periods of high water pressure gradients enhance glacier surface velocity and channel sediment transport by up to 30% and 50%, respectively. We expect the application of this framework to future seismic observations acquired at glaciers worldwide to improve our understanding of subglacial processes.

The Rospuda subglacial channel hosts not only deep ribbon lakes and the Rospuda River itself, but also an outwash plain. The paper deals with a lithofacies analysis of deposits from exposures within the Rospuda outwash plain. With the support of geomorphological and geological analyses, it has been found that the outwash plain formed along the Rospuda subglacial channel that was being exposed from under the ice cover. The lithological record of the Rospuda outwash plain shows a distinct division into a proximal and a distal zone. The deposits under study reveal a record of a break in basial sedimentation, when ice-wedge casts developed under conditions of periglacial climate. The development of the Rospuda outwash plain is linked with the ice-marginal zone of the Pomeranian phase, similarly to the Olecko-Rajgród outwash plain that is located to the west.

The Neutral Hills Uplands of southern Alberta, Canada is an area of complex and varied glacial landforms dominated by glacitectonic compressional structures but also containing expansive areas of hummocky terrain and kame and kettle topography. It lies between the strongly streamlined trunks of the former Central Alberta (CAIS) and Maskwa palaeo-ice streams of the SW Laurentide Ice Sheet (LIS) and hence comprises an inter-ice stream regional moraine zone, constructed at around 15.5 cal ka BP. This study aimed to compile a regional map of the glacial geomorphology of central southeast Alberta in order to decipher the landform-sediment signatures of overprinted ice stream margins in terrestrial continental environments, and to refine the palaeoglaciological reconstructions for the southwest LIS. Detailed mapping from LiDAR and aerial imagery identifies distinctive glacial landsystems diagnostic of the partial overprinting of cross-cutting ice stream trunks and fast flow lobes. Widespread evidence of surge-diagnostic features indicates that the ice streams experienced repeated flow instabilities, consistent with the broader scenario of a highly dynamic and unstable SW LIS, characterised by markedly transitory and cross-cutting palaeo-ice streams. The inter-ice stream moraine zone is characterised by spectacular glacitectonic compression of bedrock, cupola hill construction and mega raft displacement but also displays evidence of multi-phase stagnant ice melt-out, where partially overprinted surge lobes advanced into large areas of buried glacier ice. Contemporaneous ice melting led to the widespread development of glacier karst and the production of eskers at a range of scales, the largest of which record deranged drainage patterns indicative of ice-walled channel sedimentation controlled by the regional bedrock slope towards the northeast. These process-form regimes have created a significant local relief that is a product of not only glacitectonic compression of bedrock but also the creation and melting of a melange of ice and bedrock/sediment blocks of variable ice volume, which are representative of former buried snout ice with a glacier karst system that was repeatedly proglacially thrust due to surging. Widespread evidence for subglacial channel cutting is likely strongly linked to the transitory, surging and cross-cutting nature of the palaeo-ice streams in the region, whereby ice streams switched on and surged in response to the build-up, migration and marginal outbursts of subglacial water reservoirs. In addition to the reduced basal friction caused by the low permeability of the Cretaceous bedrock, pressurized groundwater and potentially also shallow biogenic gas deposits were likely important to the process-form regimes of surging lobes of soft-bedded ice streams in a region where ice flow was against an adverse bed slope; a scenario that gave rise to a variety of enigmatic landforms such as doughnuts, doughnut chains, apparent blow-out features and possible till eskers, as well as glacitectonic mega-rafts.

Being a determinant factor of the glacier’s dynamic, subglacial water behavior needs a special attention.Water flowing from the glacier’s surface is the principal source supplying the subglacial drainage system. Therefore, insight into the state and evolution of the supraglacial drainage system is crucial for recognition of recharge pattern of the englacial and subglacial drainage pathways. Climate warming causes increased ablation generating higher amount of meltwater and thinning of glacier. Decadal timescale evolution of the supraglacial drainage leads to some modifications of the system in opposition to its nearly stable state on an annual timescale. For two studied glaciers Hansbreen and Werenskioldbreen in southern Svalbard surface meltwater is the main runoff component. During the ablation season 2015, 72.5 % of the total amount was provided by meltwater and 27.5 % by precipitations. Supraglacial catchments were determined on the high resolution digital elevation model using standard watershed modelling tool in ArcGIS, for each water-input area (WIA). Spatialized water runoff calculations for all the main WIAs have been done. Having data on the water sources from catchments delimited on glacier’s surface, modelling of a theoretical pattern of subglacial conduits was done considering discrete water recharge via moulins, shear fractures or crevasses. Classical modelling with an assumption of homogeneous water supply was done for comparison. Several water pressure conditions have been taken into account as well. Results show that models of subglacial drainage system with homogeneous water recharge are more realistic for tidewater glaciers with rather broad permeable firn areas and creased frontal zones, while discrete water recharge models are better for land-terminating glaciers with almost continuous impermeable superficial cold ice layer. Subglacial channel models are assumed to be valid for a minimum period of two decades taking into account evolution of supraglacial drainage system and ice thickness changes of Svalbard polythermal glaciers.

We used a finite element model to interpret anti-correlated pressure variations at the base of a glacier to show the importance of stress redistribution in the basal ice. Two pairs of load cells are installed 20 m apart at the base of the 210 m thick Engabreen glacier in Northern Norway. Pressurisation of a subglacial channel located over one pair led to anti-correlation in pressure between them. A full Stokes 3D model of a 210 m thick and 25-200 m wide glacier with a pressurised subglacial channel represented as a pressure boundary condition was used to investigate the anti-correlated response at the bed. The model reproduced the anti-correlated pressure response at the glacier bed and variations in pressure of the same order of magnitude as the load cell observations. The anti-correlation pattern was shown to depend on the bed/surface slope such that the anti-correlated pressure variations were reproduced at a distance greater than 10-20 m from the channel when the bed slope was zero, whereas anti-correlation occurred within 10 m of the channel when the bed was inclined by 5 degrees. Pressurisation of the channel led to lateral or vertical ice flow away from the channel, support of the overlying ice and reduction of the normal stress on the bed. If the modelled cross-section was laterally constrained and the bed flat, the resulting bridging effect diverted some of the normal forces acting on the bed to the sides. In contrast, if the bed was inclined, then channel support was vertical only. The model showed that the effect of stress redistribution depends on the slope as well as the geometry of the subglacial channel and glacier, and can lead to an opposite response in pressure at the same distance from the channel.

The region of the Lubusz Lakeland in western Poland where there are a lot of subglacial channels provides opportunity for multi-proxy palaeoenvironmental reconstructions. None of them has not been the object of a specific study. The developmental history of the palaeolakes and their vicinity in the subglacial trough Jordanowo-Niesulice, spanning the Late Glacial and beginning of the Holocene, was investigated using geological research, lithological and geomorphological analysis, geochemical composition, palynological and archaeological research, OSL and AMS-radiocarbon dating. Geological research shows varied morphology of subglacial channel where at least two different reservoirs functioned in the end of the Last Glacial period and at the beginning of the Holocene. Mostly during the Bølling-Allerød interval and at the beginning of the Younger Dryas there took place melting of buried ice-blocks which preserved the analysied course of the Jordanowo-Niesulice trough. The level of water, and especially depth of reservoirs underwent also changes. Palynological analysis shows very diversified course of the Allerød interval.

The theoretical basis of subglacial channel dynamics can be traced back to the work of Röthlisberger (1972) and Nye (1953). Röthlisberger (1972) considered the channels’ behaviour to be governed by a mix between water friction melting back the channel walls and the viscous closure of the surrounding ice; Nye (1953) derived a viscous closure rate for the ice. While their modelling is evidently well constructed, two aspects of their work have gone undeveloped. The first is the consideration of a finite glacier depth within the viscous closure law, instead of the assumption of an infinite glacier depth. The second is the allowance of a region of open channel flow, so that a channel’s water may transition from a region of closed channel flow to one where the water is exposed to the atmosphere. This paper helps close these two gaps, showing how Nye’s equation for the rate of ice closure can be modified, and how the point of transition between closed and open channel flow may be determined.

We use a simple numerical model to test whether surface water influx to the bed of the interior Greenland ice sheet has the potential to cause significant subglacial channel growth similar to that observed near the ice-sheet margin and at alpine glaciers. We examine the effects on channel growth from (1) rapid supraglacial lake drainage events and (2) sustained water input into moulins. By assuming that all drainage occurs through subglacial channels and by prescribing favorable pressure conditions at the domain inlet, the model can provide upper bounds on channel growth. Our results indicate that R-channels do not grow significantly within the limited period of high pressure associated with lake drainage events. Subsequent channel growth only occurs with sustained pressures above overburden. Rapid closure of channels at low pressures suggests channels in the interior are unlikely to draw significant quantities of water from nearby distributed networks. These results indicate that other drainage mechanisms such as turbulent sheets or linked-cavity networks are likely to be of greater importance for interior subglacial drainage than the growth of channels.

Interactions between subglacial hydrology and the ocean make the existence of estuaries at the grounding zones of ice sheets likely. Here we present geophysical observations of an estuary at the downstream end of the hydrologic system that links the active subglacial lakes beneath Whillans Ice Stream to the ocean beneath the Ross Ice Shelf, Antarctica. This subglacial estuary consists of a hydropotential low upstream of the grounding zone, which is linked to the ocean by a hydropotential trough and a large subglacial channel. This subglacial channel, which is imaged using active source seismic methods, has an apparent width of 1 km and a maximum depth of 7 m. The hydropotential trough continues upstream of the grounding zone and results from an along-flow depression in surface elevations. Pressure differences along the trough axis are within a range that can be overcome by tidally induced processes, making the interaction of subglacial and ocean water likely.

Seasonal evolution of the subglacial drainage system in the lower ablation area of the Hailuogou glacier ice tongue were revealed by repeated dye tracer (Rhodamine WT) experiments during the 2009 ablation season. Between April and October, 18 dye tracer experiments were conducted by injecting the tracer at one location of the lower ablation area of the Hailuogou Glacier to diagnose the seasonal variation of the subglacial drainage system of this section of glacier ice tongue. Using a simple advection-dispersion model (ADM), the flow velocity, hydrodynamic dispersion coefficient, and degree of tracer spreading were deduced. Tracer transit velocity through the tested subglacial channel varied from 0.148 to 0.555 m s&minus, 1 during the 2009 ablation season. Dispersivity showed a relatively high value than that found at other glaciers, which varied between 27.05 and 287.49 m2 s&minus, 1. Seasonal changes of these indexes indicated that the subglacial drainage system of the lower ablation area of the Hailougou Glacier is a relatively stable existing system in the case of its longitudinal shape, whereas its hydraulic efficiency is low in the early and late ablation seasons and high during the middle of summer due to subglacial channel expansion.

The Brampton kame belt represents one of the largest glaciofluvial complexes within the UK. It is composed of an array of landform-sediment assemblages, associated with a suite of meltwater channels and situated within a palimpsest landscape of glacial features in the heart of one of the most dynamic parts of the British–Irish Ice Sheet. Glacial geomorphological mapping and sedimentological analysis have allowed a detailed reconstruction of both the morphological features and the temporal evolution of the Brampton kame belt, with processes informed by analogues from modern ice margins. The kame belt demonstrates the development of a complex glacier karst typified by the evolution of subglacial meltwater tunnels into an englacial and supraglacial meltwater system dominated by ice-walled lakes and migrating ice-contact drainage networks. Topographic inversion led to the extensive reworking of sediments, with vertical collapse and debris flows causing partial disintegration of the morphology. The resultant landform comprises a series of kettle holes, discontinuous ridges and flat-topped hills. The Pennine escarpment meltwater network, which fed the Brampton kame belt, is composed of an anastomosing subglacial channel system and flights of lateral channels. The Brampton kame belt is envisaged to have formed during the stagnation of ice in the lee of the Pennines as ice retreated westwards into the Solway Lowlands. The formation of the Brampton kame belt also has particular conceptual resonance in terms of constraining the nature of kame genesis, whereby an evolving glacier karst is a key mechanism in the spatial and temporal development of ice-contact sediment-landform associations.

The sedimentology of three sand-cored drumlins in central Ulster. Northern Ireland is described. The sand-cores of the drumlins consist of complex sequences of cohesive sediment gravity forms, grain-flow deposits, turbiditic sands and deltaic sands and gravels, showing complex fades relationships. The bulk of the sediments formed in water-filled cavities associated with a major subglacial meltwater escape route excavated in bedrock. The sand-cores are overlain by a streamlined carapace of basal melt-out till. A model is presented to account for spatial variations in faeies arrangement within the sand-cores. Three specific meltwater depositional subenvironments were recognised: (1) main channel axis, (2) main channel flank, and (3) minor tributary channel. The sand-cores formed prior to the main phase of basal melt-out till deposition and before drumlin streamlining.

This work is presented in two parts. Part I presents observations on the coupling between subglacial channel flow and groundwater flow in determining subglacial hydraulic regime and creating eskers from an Icelandic glacier that is suggested as an analogue for many parts of Pleistocene ice sheets. Part II develops a theory of perennial subglacial stream flow and the origin of esker systems, and models the evolution of the subglacial stream system and associated groundwater flow in a glacier of the type described in Part I. It is suggested that groundwater flow may be the predominant mechanism whereby meltwater at the glacier bed finds its way to the major subglacial streams that discharge water to glacier margins.

A model subglacial drainage system, coupled to an ice-dammed reservoir that receives a time-varying meltwater input, is described and analysed. In numerical experiments an ice-marginal lake drains through a subglacial channel, producing periodic floods, and fills with meltwater at a rate dependent on air temperature, which varies seasonally with a peak value of Tm. The analysis reveals regions of Tm parameter space corresponding to ‘mode locking’, where flood repeat time is independent of Tm; resonance, where decreasing Tm counter-intuitively increases flood size; and chaotic dynamics, where flood cycles are sensitive to initial conditions, never repeat and exhibit phase-space mixing. Bifurcations associated with abrupt changes in flood size and timing within the year separate these regions. This is the first time these complex dynamics have been displayed by a glaciohydraulic model and these findings have implications for our understanding of ice-marginal lakes, moulins and subglacial lakes.

Subglacial drainage systems exert a major control on basal-sliding rates and glacier dynamics. However, comparatively few studies have examined the sedimentary record of subglacial drainage. This is due to the paucity of modern analogues, the limited recognition and preservation of upper flow regime deposits within the geological record, and the difficulty of distinguishing subglacial meltwater deposits from other meltwater sediments (e.g. glacier outburst flood deposits). Within this study, the sedimentological and structural evolution of a subglacial to subaerial (ice-marginal/proglacial) drainage system is examined. Particular emphasis is placed upon the genetic development and preservation of upper flow regime bedforms and specifically recognising them within a subglacial meltwater context. Facies are attributed to subglacial meltwater activity and record sedimentation within a confined, but progressively enlargening, subglacial channel system produced under dune to upper flow regime conditions. Bedforms include rare large-scale sinusoidal bedding with syn-depositional deformation produced by current-induced traction and shearing within the channel margins. Subglacial sedimentation culminated with the abrupt change to a more ephemeral drainage regime indicating channel-abandonment or a seasonal drainage regime. Retreat of the ice margin, led to the establishment of subaerial drainage with phases of sheet-flow punctuated by channel incision and anastomosing channel development under diurnal, ablation-related, seasonal discharge. The presence of extensive hydrofracture networks demonstrate that proglacial groundwater-levels fluctuated markedly and this may have influenced later overriding of the site by an ice stream.

The Wda River valley in northern Poland is a polygenetic valley located in a former subglacial meltwater channel that after ice sheet retreat hosted separate and then interconnected lake basins, subsequently replaced by a river. In one part of the river system an abandoned Lateglacial valley is found. This dry valley is a unique feature in the area located within the limit of the Scandinavian Ice Sheet during the Last Glaciation in the Central European Lowland. We investigated lacustrine and fluvial sediments and landforms in the valley and applied palynological analysis combined with radiocarbon dating to reconstruct the valley evolution from the early period of river inception in the Lateglacial through the formation of lakes during the Pleistocene/Holocene transition up to the establishment of a modern river system in the early Holocene. Three coeval processes were identified: (1) formation of lake basins in the subglacial channel connected by a river, (2) erosion of the river bed between the lakes, and (3) sediment deposition at the mouths of the channels in the lake basins first generating delta fans and eventually filling the lakes entirely. The valley formation was associated with diachronous melting of dead ice blocks buried in the subglacial channel controlled by the capacity of the hydrological system to evacuate the meltwater. Most of the lake basins in the study area were formed during the Bolling–Allerod period, but one section of the subglacial channel not affected by thermokarst processes survived protected by dead ice blocks throughout the entire Lateglacial. Accelerated dead-ice decay at the beginning of the Holocene triggered widespread lowering of the base level, which caused lake disappearance and the formation of the modern river system. The processes reconstructed in the central section of the Wda River indicate a highly dynamic and diachronous river valley development during the Lateglacial and early Holocene whereby the valley formation was interwoven with the formation of lake basins within the valley triggered by melting of dead ice blocks in the substratum.

We recorded electrical resistivity data at the base of four boreholes drilled through Haut Glacier d’Arolla, Switzerland. The data were acquired repetitively every hour over two diurnal hydrological cycles in the late melt season, separated by 10 days. Constrained three-dimensional (3-D) data inversion allowed reconstruction of hourly variations in bulk resistivity in the subglacial sediment layer. Inverted resistivity models reflect the establishment of channelized subglacial drainage in the study area between the two hydrological cycles, in agreement with previous work. Daily variations in bulk and water resistivity are in phase, and bulk resistivity amplitudes decrease away from the subglacial channel. Using selected electrical–hydraulic relationships, we estimate metre-scale changes in the hydraulic conductivity and porosity of the subglacial sediment layer, accounting for increasing clay content and decreasing median grain radius with distance from the channel. Hydraulic conductivity and porosity were respectively calculated to decrease from (6.4 ± 2.1) × 10–2ms–1 and 0.34 ± 0.01 at the channel to (3.3 ± 2.2) × 10–2ms–1 and 0.26 ± 0.01 at a distance of 5m from it. The hydraulic conductivity estimates are in agreement with previously inferred values, and the porosity estimates fall within the expected range for unlithified subglacial sediments. We conclude that collection and inversion of repeat 3-D subglacial resistivity data is feasible and has the capacity to generate multidimensional images of subglacial hydraulic processes and properties.

We illustrate here spectacular meltwater features associated with outburst floods beneath an ice sheet that overrode the Transantarctic Mountains in southern Victoria Land. Because of long-term hyperarid polar climate, these features are part of an ancient landscape preserved for about 14 million years. Some channels are associated with areal scouring of basement rocks extending from sea level to as much as 1200–2100 m elevation in coastal regions. Scablands with scallops, potholes and plunge pools are cut in Beacon Super group sandstones and Ferrar Dolerite and cover wide areas of high western plateaus near the mountain crest. Subglacial channel systems commonly originate near divides and converge downhill toward the northeast. We argue that the landforms were created beneath a major Antarctic Ice Sheet that submerged the whole area, with the possible exception of the high peaks of the Royal Society Range, as it flowed northeastward toward the outer Antarctic continental shelf. Areal scouring, associated with warm-based regimes, is restricted to the lower slopes close to the coast. In the higher terra in, meltwater channels and scabland alongside preserved patches of regolith are best explained by the breaching of cold-based ice on the mountain rim by subglacial melt water outbursts. Melt from warm-based ice, along with subglacial lakes trapped upstream of the mountain rim, are possible sources of the meltwater necessary to form the channel systems and scablands.

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66° N and 66°30′N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Röthlisberger’s equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Data from the mouth of the decelerating Whillans Ice Stream (WIS), West Antarctica, spanning 42 years are reviewed. Deceleration has continued, with local areas of both thinning and thickening occurring. The mean thinning rate is 0.48 ± 0.77 ma–1. No consistent overall pattern is observed. Ice thickens immediately upstream of Crary Ice Rise where deceleration and divergence are strongest, suggesting expanded upstream influence of the ice rise. Thinning is prevalent on the Ross Ice Shelf. Grounding-line advance at a rate of 0.3 km a–1 is detected in a few locations. Basal stresses vary across an ice-stream transect with a zone of enhanced flow at the margin. Marginal shear is felt at the ice-stream center. Mass-balance values are less negative, but larger errors of earlier measurements mask any possible temporal pattern. Comparisons of the recent flow field with flow stripes suggest WIS contributes less ice to the deep subglacial channel carved by Mercer Ice Stream and now flows straighter. The general lack of geometric changes suggests that the regional velocity decrease is due to changing basal conditions.

In the empirical study of jökulhlaups, the peak discharge, Qmax, and water volume drained by the ice-dammed lake during the floods, Vt, appear to follow a power-law relation , where K are b are constants determined from field data. First identified by Clague and Mathews (1973), this relation is a useful reference for predicting flood magnitude, but its physical origin remains unclear. Here, we develop the theory that connects it to contemporary models for simulating the flood hydrograph. We explain how the function Qmax = f(Vt) arises from Nye’s (1976) theory of time-dependent water flow in a subglacial channel coupled to a lake, and we describe how discharge–volume data record the (monotonically increasing) form of this function so long as the lake is not emptied in the floods. The Grímsvötn jökulhlaups present an example where, because of partial draining of the lake, agreement between the model-derived f and data is excellent. It is documented that other lake systems drain completely, but we explain how the exponent b ≈ 2/3 observed for them collectively is due primarily to a scaling effect related to their size, modified by other factors such as the flood initiation process.

Tracer injections into a subglacial channel at Unteraargletscher, Switzerland, were repeated at intervals of about 2 hours over two diurnal discharge cycles in August and September 2000. Records of dye concentration reveal a pronounced hysteresis in the velocity–dispersion relationship, thereby indicating alterations in the drainage system. Theoretical considerations for Röthlisberger channels suggest an evolution of the conduit cross-section in response to a diurnally varying discharge. We studied the relation between conduit cross-section and tracer dispersion with numerical tracer experiments. The velocity field for steady flow through a given conduit geometry is calculated using a commercial flow solver. Tracer transport is represented by a scalar volume which is advected by the velocity field. Experiments were conducted for several scenarios by varying flow velocity, conduit geometry and conduit roughness. Results show only a weak dependence of dispersion on conduit size. In contrast, changes in roughness of the conduit walls reveal a strong effect on tracer dispersion. Therefore, to explain the observed hysteresis in the velocity–dispersion relationship, we suggest that the evolution of a subglacial flow path might involve changes in roughness.

We explore a mathematical model that couples together a thermomechanically evolving subglacial channel, distributed cavity drainage, and basal sliding along a subglacial flood path fed by a jökulhlaup lake. It allows water transfer between channel and cavities and a migrating subglacial water divide or ‘seal’ to form between floods. Notably, it accounts for full coupling between the lake and subglacial drainage in terms of both discharge and pressure, unlike models that neglect the pressure coupling by imposing a known history of lake discharge at the channel inlet. This means that flood hydrographic evolution and its impact on glacier motion are consistently determined by our model. Numerical simulations for a model alpine lake yield stable limit cycles simulating repeating jökulhlaups, with the channel drawing water from the cavities at a varying rate that modulates basal sliding during each flood. A wave of fast sliding propagates down-glacier at flood initiation, followed by deceleration as the growing channel sucks water from the cavities. These behaviours cannot be correctly simulated without the full coupling. We show that the flood’s peak discharge, its initiation threshold and the magnitude of the ‘fast sliding’ wave decrease with the background water supply to the cavities.

Late-summer subglacial water pressures have been measured in a dense array of boreholes in the ablation area of Haut Glacier d’Arolla, Switzerland. Interpolated surfaces of minimum diurnal water pressure and diurnal water-pressure variation suggest the presence of a subglacial channel within a more widespread, distributed drainage system. The channel flows along the centre of avariable pressure axis(VPA), some tens of metres wide, that is characterized by low minimum diurnal water pressures (frequently atmospheric) and high diurnal water-pressure variations. These characteristics are transitional over a lateral distance of c. 70 m to higher and more stable subglacial water pressures in the adjacent distributed system. Water-pressure variations recorded in boreholes located close to the centre of the VPA reflect the delivery of surface-derived meltwater to the glacier bed and result in a diurnally reversing, transverse hydraulic gradient that drives water out from the channel into the distributed system during the afternoon and back to the channel overnight. Subglacial observations suggest that such flow occurs through a vertically confined sediment layer. Borehole turbidity records indicate that the resulting diurnal water flows are responsible for the mobilization and transport of fine debris in suspension. Analysis of the propagation velocity and amplitude attenuation cf the diurnal pressure waves suggests that the hydraulic conductivity of the sediment layer decreases exponentially with distance from the channel, falling from c. 10−4m s−1at the channel boundary to c. 10−7m s−170 m away. These apparent hydraulic conductivities are consistent with Darcian flow through clean sand and typical glacial till, respectively.We suggest that fine material is systematically flushed from basal sediments located adjacent to large, melt-season drainage channels beneath warm-based glaciers. This process may have important implications for patterns of glacier erosion, hydro-chemistry and dynamics.

The Bornhoved tunnel valley consists of a relatively broad, elongated depression filled primarily with tills, and a narrow channel in the axial part filled with outwash, glaciolacustrine sediments and tills. The tunnel valley is at least 13 km long, 222 m deep (bottom at 191 m b.s.l.), and extends parallel to the axis of a peripheral sink between two Permian salt diapirs in the substratum. It is suggested that the initiation of the tunnel valley was controlled by a differential response of the bed materials to ice overriding on top of and between the salt structures. The Bornhoved tunnel valley is a polygenetic feature that resulted from different erosional processes throughout all three main glaciations in northern Germany. During the Elsterian Glaciation glaciotectonic squeezing occurred, to be later replaced by subglacial meltwater erosion. During the ice retreat and the Holsteinian Interglacial, the valley was filled with fine-grained glaciolacustrine and marine sediments, which were subsequently removed and redeposited further down-ice as a 200-m-thick glaciotectonic melange by the first Saalian advance. The second Saalian advance caused only minor changes and the tunnel valley was reactivated as a subglacial channel during the Weichselian Glaciation. It is suggested that at the time of the first Weichselian ice advance a large subglacial water reservoir developed in the area of the Baltic Sea basin and caused a rapid, surge-like ice movement. As the ice sheet advanced out of the Baltic Sea basin, drainage of the water reservoir was prevented by the ice toe overriding the permafrost on the Saalian highlands. During the ice retreat, frozen ground was left beyond the ice margin and subglacial meltwater catastrophically drained through the tunnel valley with a discharge roughly estimated at about 3.75 × 103 m3/s. It seems possible that drainage of the subglacial water collector through the Bornhoved tunnel valley and other similar channels south of the Baltic Sea led to increased ice-sheet basal friction, and was one reason for successively smaller extents of the subsequent Weichselian advances. Evidence against pervasively deforming bed materials in the vicinity of the Bornhoved tunnel valley in the Weichselian Glaciation indicates that the valley served as an efficient drainage path for porewater from the subglacial sediments.

[1] Tracer injections into a moulin connecting to a subglacial channel at Unteraargletscher were repeated at intervals of about two hours over two diurnal discharge cycles in August and September 2000. The pronounced hysteresis in the relationships between tracer transit velocity and proglacial discharge may be explained by adjustments of the channel cross section in response to discharge variations or alternatively by modulation of inflow from the tributary moulin into a main subglacial channel. To test these hypotheses, we employ a physically based, time-dependent model of channelized subglacial drainage. In a series of numerical experiments, we explore the general characteristics of each mechanism before the model is applied to diurnal discharge variations measured during the tracer tests at Unteraargletscher. We found that the ability of a Rothlisberger channel to adjust its size to the prevailing hydraulic conditions can contribute to the hysteresis in the velocity-discharge relationships. However, using plausible parameter values, R channel hydraulics alone is not satisfactory in reproducing the magnitude of observed transit velocities. The model further confirms that retardation of the tracer in the tributary moulin due to inflow modulation has major effects on transit velocities both in magnitude and timing. With appropriate assumptions on channel sinuosity and moulin size, a combination of both mechanisms can be found that is capable of reproducing the observed transit velocities.

Interpretation of sidescan sonar and high resolution sub-bottom seismic data from Hudson Bay has led to the development of a model for final disintegration of the Laurentide Ice Sheet. The well-preserved geomorphic features observed from the bay floor are similar to subaerially exposed glaciogenic features observed around the perimeter of Hudson Bay, and can be related to the following three paleoenvironments: subglacial erosion/deposition, ice margin deposition and reworking, and post-glacial reworking and burial. Subglacial features observed include fluted terrains with superimposed rogen moraines, relief attributed to ice surging, eskers, subglacial channels and dead ice topography. De Geer moraines and recessional moraine complexes indicate the mechanism and orientation of ice breakup/retreat. The orientations of these features are locally consistent with regional trends indicating the overall direction of ice sheet flow and retreat during the last phase of the Laurentide Ice Sheet. Subglacial channel complexes incised into the till and bedrock, and high energy bedforms deposited on the till surface occur in present day water depths down to 160 m. Due to isostatic depression these could have been as much as 315 m deeper at the time of formation. They indicate a wet-based glacier with high energy subglacial drainage during the late phases of ice retreat. The sidescan sonar data also reveal numerous zones of distinctive morphology attributed to ice breakup and post-glacial processes. Cross-cutting relationships allow the determination of the relative ages of these zones. From oldest to youngest these zones are interpreted to represent scouring of the bedrock and overburden by the remnants of the ice sheet, seafloor scouring by tabular icebergs which moved in a tidally induced(?) arcuate pattern or under the influence of massive glacial lake outflows, development of dead ice topography under residual ice masses, and random iceberg scouring by smaller icebergs restricted to present day water depths of less than 180 m. Recently published reconstructions of the collapse of the Laurentide Ice Sheet, based on dated Quaternary sections on land, suggest that the bay was rapidly deglaciated in approximately 400 years. The well-preserved and distinctive glaciogenic and stratigraphic indicators observed from the bay floor suggest that they were formed by a succession of mechanisms which were too short lived to obliterate the underlying evidence of deglaciation. This supports the interpretation of extremely rapid deglaciation of Hudson Bay and points out that subsequent marine processes have not been sufficient to obliterate the glacial features.

A sedimentary sequence of till overlying a gravel aquifer was instrumented with water-pressure transducers prior to a small, anticipated surge of the margin of the glacier Breiðamerkurjökull in Iceland. The records of water pressure at each transducer site show a well-defined temporal sequence of hydraulic regimes that reflect the changing recharge of surface-derived meltwater, the pressure drop along the drainage pathway and the pattern of ice loading. The poroelastic and water-pressure response of glacially overridden sediments to the recharge rate is determined in the frequency domain through an analytic solution. This permits the in situ conductivity, compressibility and consolidation states of subglacial sediments to be derived, and reveals aquifer-scale compressibility that produces an important water-pressure wave associated with the advancing glacier. The model is then used to explore how varying conductivity/compressibility, largely determined by granulometry, can determine drainage states and instabilities that may have a large impact on glacier/ice-sheet dynamics, and how the drainage time of surface water to the bed can determine the frequency response of subglacial groundwater regimes and their influence on subglacial sediment stability. Mismatches between model predictions and specific events in water-pressure records are used to infer processes that are not incorporated in the model: hydrofracturing that changes the hydraulic properties of subglacial sediments; the impact on groundwater pressure of subglacial channel formation; upwelling beyond the glacier margin; and rapid variations in the state of consolidation. The poroelastic model also suggests how seismic methods can be developed further to monitor hydraulic conditions at the base of an ice sheet or glacier.

[1] This paper uses a combination of field data and three-dimensional modeling to investigate the spatial variability in basal conditions required to induce observed fluctuations in diurnal ice velocity at Haut Glacier d'Arolla, Switzerland. A network of surface velocity markers was observed at intervals of as little as four hours over diurnal cycles in both winter and late summer. Winter motion showed limited diurnal variability, presumably due to the absence of supraglacial meltwater inputs. By contrast, diurnal fluctuations in ice motion were recorded in summer across the lower and upper glacier. In the lower glacier, surface velocities were intimately linked to hydrological forcing in the vicinity of a subglacial channel. Previously observed diurnal excursions of meltwater away from the channel should reduce areas of basal drag adjacent to the channel thereby impacting on ice dynamics. Using a first-order ice flow approximation, we investigated the distribution of basal shear traction adjacent to the channel necessary to replicate the observed surface velocity field during periods of rapid ice motion. The modeling suggests that the observed variations in diurnal velocity will only occur with extensive reductions in basal drag across a transverse zone of up to 560 m across, well beyond the immediate vicinity and previously observed extent of diurnal excursions of meltwater away from the subglacial channel.

Slug testing allows estimation of subglacial hydraulic properties by evaluating the response of a coupled borehole–subglacial flow system to an artificial displacement of its steady-state hydraulic head. However, existing models developed specifically for application to slug-test data collected in glacier boreholes are challenging to apply in practice. Here, we demonstrate that conventional linear methods, which are relatively readily applicable and widely used in groundwater studies, can also be used to estimate subglacial hydraulic properties. Overdamped, underdamped and critically damped slug-test data were recorded in fluctuating boreholes drilled to the bed of Haut Glacier d’Arolla, Valais, Switzerland. We find that non-linear effects in the data are negligible, supporting the application of conventional hydrogeological methods. Results suggest that the hydraulic conductivity of the unconsolidated sediments within the area influenced by a major subglacial drainage channel is enhanced (several 10–2 m s–1; typical of gravelly sand) compared to areas outside the zone of influence (~0.1 × 10–2ms–1; typical of silty sand). A distance to a flow boundary within the subglacial aquifer, inferred on the basis of slug-test analysis, was found to coincide with the actual location of the subglacial channel. Sensitivity analyses reveal that uncertainties in borehole and filter radii, as well as the storage coefficient of the subglacial aquifer, are more likely to affect the accuracy of the hydraulic and distance estimates than uncertainties regarding the repeatability of the slug-test responses and curve-fitting procedures involved in the conventional hydrogeological methods.

[1] The morphology of the drainage system of Unteraargletscher, Switzerland, and the diurnal variability of drainage conditions were investigated by conducting a series of tracer tests over a number of discharge cycles during the ablation season 2000. Dye injections into a moulin were repeated at intervals of a few hours and were accompanied by simultaneous measurements of discharge of supraglacial meltwater draining into the moulin and bulk runoff in the proglacial stream. Fast transit velocities in conjunction with low dispersion values suggest that the tracer was routed through a hydraulically efficient, channelized drainage system. However, detailed analyses reveal a large diurnal variability in terms of transit velocity and dispersion coefficient. This finding underlines the difficulty of detecting a possible drainage system evolution based on single tracer tests conducted at coarse temporal intervals. Furthermore, the obtained velocity-discharge relationships display pronounced hysteresis. We suggest that the evolution of the cross-sectional area of an ice-walled conduit and the modulation of inflow at the junction of a tributary moulin to a main subglacial channel are responsible for the observed behavior.

The continental shelf of Queen Charlotte Sound, British Columbia, was inundated by piedmont glaciers which extended to the shelf edge by way of the shelf transverse troughs [Josenhans et al., 1992; Luternauer et al., 1989]. Seismic reflection profiles indicate that glacial erosion overdeepened these troughs and developed a smooth glacial unconformity on bedrock. Glacial and deglacial sediments comprising basal till, multiple ice-contact sequences, and stratified glaciomarine sediments overlie the glacial unconformity and typically attain 50m in thickness. The offshore banks, which separate the troughs, appear to be devoid of glacial deposits and may have remained ice free during the last glaciation. Regional seismic profiles from the trough axes reveal multiple generations of channel fill deposits in channels cut into the Plio-Pleistocene bedrock on the inner shelf. The upper surface of the channel deposits is truncated by the glacial erosional unconformity and overlain by till. These channel deposits have not been sampled and may represent a preglacial subaerial drainage system or possibly subglacial channel sequences that were deposited prior to the last glacial advance.

In the classical drainage theory proposed by Rothlisberger [1] and later extended by Nye [2], ice deformation closure of a subglacial channel is counteracted by melt-back of the channel walls, due to the heat dissipated by water flowing in the channel. The mathematical description of this is by means of an evolution equation for the channel cross-section which has been formulated under the assumption that the cross-sectional shape is circular. For steady flow, the equation reduces to the equilibrium relation derived by Rothlisberger [1], linking water flux, water pressure and hydraulic gradient. For unsteady flow, it has been applied with success to models of glacier outburst floods, or jokulhlaups, in Iceland [e.g., 2] and elsewhere. However, it has recently been suggested that this theory can be generalized to the consideration of non-circular cross-sections. For instance, Hooke et al. [3] suggested broad and low channels as the cause of why measured subglacial water pressures are frequently underestimated by theory. Fowler & Ng [4] proposed a theory of Icelandic jokulhlaups and an improved model for simulating flood hydrographs, by considering a wide channel that has independently evolving depth and width. This configuration resembles the canals in Walder & Fowler’s extension of the classical theory for deformable sediment beds [5]. In this paper, the first steps of analyzing the mathematics of wide subglacial channels of this type are presented. In particular, we outline the formulation of the evolution equations, for the case of ‘hard bed’ drainage where sediment processes can be excluded. We then derive the corresponding equilibrium results and discuss their implications.

Medial moraines and lineations on the surface of composite glaciers enable the detection of structural and hydrological features. A study of the composite glacier Breidamerkurjökull, south Iceland, indicates relationships between subglacial waterways, ice structure in the junction area, and development of the glacier terminus. First, rivers are situated in or near medial moraines because melt water percolates to the bed and moves from there with the subglacial rivers mainly in the direction of ice flow. Secondly, the contact between two feeder glaciers sometimes forms an angle at the glacier terminus. Then the meltwater river escaping from the contact, generally in a radial direction (away from the glacier front), during the retreat will be transferred to the front of the first receding glacier. The drainage of the contact zone then changes from a radial to a tangential direction, destroying the terminal moraines of the recession stages. Similar relations are found in relics of Pleistocene ice sheets, two examples of which are compared. The huge subglacial channel of the Münsterländer Kiessandzug-Esker below the ice sheet of the Saalian glaciation in the Münsterland, north-west Germany, was formed in the ice-flow direction. It therefore gives details of the morphology and of the great ice-dammed lake east of the Teutoburger Wald ridge. At the contact between the Norwegian and Baltic Sea glaciers, the terminus formed an angle during the maximum extent of the Weichselian glaciation in northern Jutland. During retreat, the Norwegian glacier receded first. The large melt-water river escaping from the contact between both glaciers had formed the huge Karup Sandur during the maximum, but now, during the recession, it changed to the front of the Norwegian ice, destroying the recession moraines there.

The physics of water flow within and under glacier ice is examined with special reference to the periodic catastrophic outbursts of water (jökulhlaups) from the subglacial lake Grímsvötn, Vatnajökull, Iceland. The lake is sealed until it reaches a critical level which enables it to lift the glacier, helped by a hydrostatic cantilever effect. The differential equations for non-steady water flow in a subglacial tunnel are derived and applied to the 1972 Grímsvötn outburst. The discharge: time relation observed during the growth stage, and the abrupt ending of the flood, are both very well accounted for by a theory which is insensitive to the details of the subglacial tunnel system. The steady state, in which an intergranular vein or tunnel is simultaneously melted open by frictional heat and closed by plastic deformation, may be stable or unstable according to the conditions imposed at the ends. This explains why the flow of water in a vein does not normally increase unstably as in a jökulhlaup. An ice-dammed lake does not drain away through the vein system because the driving force on the vein-water is towards the lake rather than away from it.

An extensive Late Pleistocene delta body, related to a water plane at 100 m above sea level, is exposed in the Upper Carey Valley, Northern Ireland. Exposures along the modern river show a lowermost facies association, of massive diamict and poorly-sorted gravel deposited as sediment gravity flows, overlain by a gravelly foreset and topset facies association recording the progradation of a large glacier-fed delta. Diamict facies were emplaced as subaqueous debris flows; associated gravels are chaotically-bedded and were deposited by the downslope slumping of poorly-sorted boulder gravel; their subtle grading characteristics can be linked with a variety of mass flow processes. Sediment gravity flow facies originated upvalley from rapidly deposited sediments which accumulated near the exit of a narrow subglacial channel below an ice margin fronting in water. Subsequent ice retreat allowed the growth of a Gilbert-type delta. Regional evidence shows that the delta top records high relative sea levels in response to isostatic downwarping around the periphery of the last British ice sheet. The sedimentology and structure of the Carey Valley infill may be typical of high-level glaciomarine deltas preserved in topographic lows around the coastline of north Britain.

We discuss the suggestion that small underwater transmitters might be used to illuminate the interior of major englacial water channels with radio waves. Once launched, the radio waves would naturally tend to be guided along the channels until attenuated by absorption and by radiative loss. Receivers placed within the channels or at the glacier surface could be used to detect the signals. They would provide valuable information about the connectivity of the water system. The electrical conductivity of the water is of crucial importance. A surface stream on Storglaciären, in Sweden, was found, using a low-frequency technique, to have a conductivity of approximately 4 × 10−4 S m−1. Although this is several hundred times higher than the conductivity of the surrounding glacier ice, the contrast is not sufficient to permit us simply to use electrical conductivity measurements to establish the connectivity of englacial water channels. However, the water conductivity is sufficiently small that, under favourable circumstances, radio signals should be detectable after travelling as much as a few hundred metres along an englacial water channel. In a preliminary field experiment, we demonstrated semi quantitatively that radio waves do indeed propagate as expected, at least in surface streams. We conclude that under-water radio transmitters could be of real practical value in the study of the englacial water system, provided that sufficiently robust devices can be constructed. In a subglacial channel, however, we expect the radio range would be much smaller, the environment much harsher, and the technique of less practical value.

A correction to Rothlisberger's model of the pressure melting point effect in subglacial or englacial channels is evaluated. The correction term can be very important near the terminus. Elsewhere it predicts a decrease in piezometric pressure gradients between about ten and seventeen percent. It is shown that subglacial channel pressures can be much higher on a wavy bed than on a flat bed with the same average slope. An instability phenomenon is discussed which can occur when the flux decreases in overpressurized englacial channels. The effects of cold ice upon piezometric pressure gradients are investigated. Pressure gradients are reasonably small provided that the channel flux is sufficiently high. It is doubtful that channels can be maintained through cold ice during the winter.

The paper attempts to relate four recent sets of measurements on the flow laws of ice: (1) Glen’s laboratory compression tests, (2) Gerrard, Perutz & Roch’s pipe experiment on the Jungfraufirn, (3) measurements of contraction rate in a tunnel in Skauthöe Glacier by McCall, and (4) measurements of contraction rate in tunnels in the Z’Mutt and Arolla Glaciers reported by Haefeli & Kasser. When they are analyzed within the framework of a general flow theory, satisfactory agreement is found between all these results, with one exception: the reported rate of collapse of the Arolla tunnel is faster than would be given by the present theory. The paper includes the theory of the slow collapse of a cylindrical or a spherical hole in a medium showing non-Newtonian viscosity.

A network of passages situated along three-grain intersections enables water to percolate through temperate glacier ice. The deformability of the ice allows the passages to expand and contract in response to changes in pressure, and melting of the passage walls by heat generated by viscous dissipation and carried by above-freezing water causes the larger passages gradually to increase in size at the expense of the smaller ones. Thus, the behavior of the passages is primarily the result of three basic characteristics: (1) the capacity of the system continually adjusts, though not instantly, to fluctuations in the supply of melt water; (2) the direction of movement of the water is determined mainly by the ambient pressure in the ice, which in turn is governed primarily by the slope of the ice surface and secondarily by the local topography of the glacier bed; and, most important, (3) the network of passages tends in time to become arborescent, with a superglacial part much like an ordinary river system in a karst region, an englacial part comprised of tree-like systems of passages penetrating the ice from bed to surface, and a subglacial part consisting of tunnels in the ice carrying water and sediment along the glacier bed. These characteristics indicate that a sheet-like basal water layer under a glacier would normally be unstable, the stable form being tunnels; and they explain, among other things, why ice-marginal melt-water streams and lakes are so common, why eskers, which are generally considered to have formed in subglacial passages, trend in the general direction of ice flow with a tendency to follow valley floors and to cross divides at their lowest points, why they are typically discontinuous where they cross ridge crests, why they sometimes contain fragments from bedrock outcrops near the esker but not actually crossed by it, and why they seem to be formed mostly during the later stages of glaciation.