Glacial erosion beneath ice streams and ice-stream tributaries: constraints on temporal and spatial distribution of erosion from numerical simulations of a West Antarctic ice stream

Geologic evidence such as subglacial troughs and grounding zone wedges indicate that soft-bedded, West Antarctic ice streams are capable of eroding, transporting and depositing large volumes of debris at high rates (˜100 m³ yr^-1 per meter width). In order to understand the dynamics of ice streams and the geologic effects of their activity, it is important to understand the physical mechanisms that control these high rates of sub-ice-stream sediment generation and transport. Here, we use a numerical model of Ice Stream C run over c. 8500 model years to quantify the effects of a recently proposed, till-ploughing mechanism of till formation and redistribution beneath ice streams (Tulaczyk et al. 2001; Clark et al. in press). Our results show that this 'transport-limited' mechanism, in which till transport rates scale with ice velocity and erosion rates with spatial gradients of velocity, is consistent with existing constraints. For instance, our model results predict significantly higher (˜0.6 mm yr^-1) average erosion rates beneath ice-stream tributaries, which are underlain by deep subglacial troughs, than beneath ice-stream trunks (˜0.2mm yr^-1), whose subglacial troughs have a significantly smaller relief. We would not obtain this satisfactory result if subglacial erosion was parametrized in the model using the more traditional approach of scaling erosion rates with ice velocity (what we call the 'production-limited' parametrization). Because of the requirement of ice continuity, the magnitude of ice velocity generally increases downstream in polar ice streams, and so do the production-limited erosion rates. Pending further investigations, we propose that geologic and geomorphic effects of soft-bedded ice streams should be quantified using some form of a 'transport-limited' parametrization of subglacial erosion rates, e.g. the till-ploughing mechanism.     

The role of meltwater in glacial processes

Water plays a dominant role in many glacial processes and the erosional, depositional and climatic significance of meltwaters and associated fluvioglacial processes cannot be overemphasized. At its maximum extent c. 20,000 years ago, the volume of the Laurentide ice sheet was 33 × 10⁶ km³ (about the same as the volume of all ice present today on planet Earth). The bulk of this was released as water in little more than 10,000 years. Pulses of meltwater flowing to the Atlantic Ocean from large ice dammed lakes altered thermohaline circulation of the world's oceans and global climate. One such discharge event via Hudson Bay at 8200 years BP released 160,000 km³ of water in 12 months. Global sea levels recovered from glacial maximum low stands reached at about 20,000 years ago at an average rate of 15 m per thousand years but estimates of shorter term rates suggest as much as 20 m sea level rise in 1000 years and for short periods, rates as high as 4 m per hundred years. Meltwaters played a key role in lubricating ice sheet motion (and thus areal abrasion) across the inner portions of the ice sheet where it slid over rigid crystalline bedrock of the Canadian Shield. The recharge of meltwater into the ice sheets bed was instrumental in generating poorly sorted diamict sediments (till) by sliding-induced shearing and deformation of overpressured sediment and soft rock. The transformation of overpressured till into hyperconcentrated slurries in subglacial channels may have generated a highly effective erosional tool for selective overdeepening and sculpting of bedrock substrates. Some workers credit catastrophic subglacial ‘megafloods’ with the formation of drumlins and flutes on till surfaces. Subglacial melt river systems were instrumental in reworking large volumes of glaciclastic sediment to marine basins; it has been estimated that less than 6% of the total volume of glaciclastic sediment produced during the Pleistocene remains on land. Fluvioglacial and glaciolacustrine sediments and landforms dominate large tracts of the ‘glacial’ landscape in North America. The recharge of subglacial meltwater into underlying bedrock and sediment aquifers created transient reversals in the long-term equilibrium flow directions of basinal fluids. With regard to pre-Pleistocene glacial record, meltwaters moved enormous volumes of terrestrial ‘glaciclastic’ sediment to marine basins and thus played a key role in preserving a record of glaciation, a record otherwise almost entirely lost on land.     

Drumlin formation in Southern Anglesey and Arvon,Northwest Wales

This is a study of Late Devensian drumlins formed in southern Anglesey and Arvon, northwest Wales. This area was affected by ice sheet coalescence when the Welsh ice sheet met with the lrish Sea ice sheet, and drumlins were formed once the two had coalesced. It is suggested that the drumlins were the result of net subglacial soft-bed erosion, and that they represent more resistant cores within the subglacial deforming layer. The drumlins have either gravel or till cores, and where the core was deformable, large-scale compressive glaciotectonic structures were seen (e.g. Dinas Dinlle) with local subglacial compression of –59%. Where the cores were more resistant (e.g. Lleiniog) these were not deformed but remained as more competent masses within the deforming layer. It is suggested that the less competent material flowed around the cores, some remaining as a thin carapace, but most of the material being removed down glacier, leaving the drumlins as erosional remnants. In northwest Wales there is a multi-till sequence that traditionally has been interpreted as having been deposited as the result of separate ice-sheet advances and retreats. However, in this study, it is suggested that the different tills were deposited as the result of ice-sheet coalescence, and that sites such as Dinas Dinlle do not show evidence of a major readvance in the retreat of the Devensian ice, but are indicative of continuously changing conditions within the subglacial deforming bed.     

Subglacial processes and drumlin formation in a confined bedrock valley,northwest Ireland

Knight, J. 2010: Subglacial processes and drumlin formation in a confined bedrock valley, northwest Ireland. Boreas, 10.1111/j.1502‐3885.2010.00182.x. ISSN 0300‐9483. Subglacial processes beneath the Late Weichselian ice sheet in northwest Ireland are deduced from sediments and structures within drumlins in a bedrock valley at Loughros Beg, County Donegal. Here, a glacially smoothed bedrock surface underlies the drumlins, which are composed on their up‐ice side of stacked, angular rafts of local bedrock. Overlying and down‐ice from these rafts are down‐ice‐dipping beds of massive to bedded diamicton that contain sand and gravel interbeds. In a down‐ice direction the diamicton matrix coarsens and the beds become laterally transitional to water‐sorted gravels. The down‐ice end of one drumlin shows a concentrically bedded stratified gravel core aligned parallel to ice flow and resembling the internal structure of an esker. With distance away from this core, the gravels become more poorly sorted with an increase in matrix content, and are transitional to massive to stratified diamicton. A four‐stage model describes the formation of drumlins in this sediment‐poor setting. The sediments that are located directly above the bedrock represent deposition in a semi‐enclosed subglacial cavity. A trigger for this process was the formation of subglacial relief by the thrusting up of bedrock rafts, which created the leeside cavity. Subsequent sediment deposition into this cavity represents a form of feedback (self‐regulation), which may be a typical characteristic of subglacial processes in sediment‐poor settings.     

Subglacial erosion forms in northwest Ireland

Subglacial erosional forms are commonly found on bedrock substrates inside the Late Weichselian ice margin in County Donegal, northwest Ireland, and can be used to provide detailed information on subglacial processes and environments. The erosional forms occur on spatial scales from whalebacks (tens of metres in scale), to asymmetric and channelized bedrock-cut scours (tens of cm in scale) and striations (mm scale). Processes responsible for development of subglacial erosional forms occur along a continuum, from free meltwater existing as a laterally extensive sheet at the ice-bed interface, to abrasion by basal ice. Channelized bedrock-cut scours are particularly common in County Donegal, and show asymmetric and meandering thalwegs, U-shaped cross-profiles and steep lateral margins. Innermost parts of the scours are highly polished and have striations that follow thalweg direction. In places, bedrock surfaces are overlain by a delicate polish and thin calcite cement, and are buried beneath glacial till. Based on their morphology, the bedrock scours are interpreted as s-forms caused by high-pressure subglacial meltwater erosion. Striations within the scoured channels reflect periods of ice-bed coupling and subglacial abrasion. The range of features observed here was used to consider relationships between subglacial topography, hydraulic processes and ice-bed coupling. Precipitation of calcite cement took place in depressions on the bedrock surface by CO₂ degassing. Infilling of depressions by glacial till formed a new type of 'sticky spot' related to spatial variations in subglacial water pressure. The temporal evolution of sticky spots reflects interactions within the subglacial environment between subglacial relief, hydraulic regime and ice-bed coupling.     

Late Quaternary submarine bedforms and ice‐sheet flow in Gerlache Strait and on the adjacent continental shelf,Antarctic Peninsula

Geophysical data from Gerlache Strait, Croker Passage, Bismarck Strait and the adjacent continental shelf reveal streamlined subglacial bedforms that were produced at the bed of the Antarctic Peninsula Ice Sheet (APIS) during the last glaciation. The spatial arrangement and orientation of these bedforms record the former drainage pattern and flow dynamics of an APIS outlet up‐flow, and feeding into, a palaeo‐ice stream in the Western Bransfield Basin. Evidence suggests that together, they represent a single ice‐flow system that drained the APIS during the last glaciation. The ice‐sheet outlet flowed north/northeastwards through Gerlache Strait and Croker Passage and converged with a second, more easterly ice‐flow tributary on the middle shelf to form the main palaeo‐ice stream. The dominance of drumlins with low elongation ratios suggests that ice‐sheet outlet draining through Gerlache Strait was comparatively slower than the main palaeo‐ice stream in the Western Bransfield Basin, although the low elongation ratios may also partly reflect the lack of sediment. Progressive elongation of drumlins further down‐flow indicates that the ice sheet accelerated through Croker Passage and the western tributary trough, and fed into the main zone of streaming flow in the Western Bransfield Basin. Topography would have exerted a strong control on the development of the palaeo‐ice stream system but subglacial geology may also have been significant given the transition from crystalline bedrock to sedimentary strata on the inner–mid‐shelf. In the broader context, the APIS was drained by a number of major fast‐flowing outlets through cross‐shelf troughs to the outer continental shelf during the last glaciation. Copyright © 2004 John Wiley & Sons, Ltd.     

Quantitative study of the distributary braidplain of the Preboreal ice-contact Gardermoen delta complex, southeastern Norway

This raised delta structure is an ice-contact deltaic complex with a volume of c. 4.4.10⁹ m³, deposited c . 9500 yr BP in a shallow wide 'fjord' during the retreat of the Scandinavian ice cap. The delta plain lies at an altitude of 200–223 m. It aggraded c . 20 m above the contemporaneous sea level during a regional marine regression. The braidplain palaeochannel characteristics indicate a peak meltwater discharge of 7–9 10³ m³/s. Calculations based on a glacial ablation model indicate a mid-summer discharge of c . 5.5 10³ m³/s. However, the fluvial topset of the delta has an erosive base whose altitude decreases upstream and indicates stream incision by more the 6 m below the contemporaneous sea level. The deep scour is ascribed to episodic floods over the relatively short delta plain, which exceeded direct ablation-associated discharges. The depositional time-span of the delta is assessed to have been 70 years, calculated from coastal gradient and shoreline displacement curves. The average sedimentation rate of the delta is thereby inferred to have been extremely high, c . 6. 10⁷ m³/yr. The sedimentation is thought to reflect 'extreme' ice-margin conditions, where the sediment and water discharge was maximized by full-scale ablation, with simultaneous subglacial, englacial and supraglacial sediment and water supply. These conditions might further coincide with an abundant rainfall in the catchment area or the drainage of dammed waters, initiating episodic floods which eroded deep beneath sea level. As a whole, the study illustrates the hydrological conditions of proglacial sedimentation at the front of the rapidly retreating last Scandinavian ice cap.     

The meltwater hypothesis for subglacial bedforms

This paper presents an historical and in places informal account of the meltwater hypothesis, which invokes enormous outburst floods for the formation of subglacial bedforms. It begins with a brief discussion of the difficulties of determining processes of formation for landforms, which are not seen in formation. Analogy provides a solution to these difficulties. Analogy between erosional marks at the bases of turbidites and drumlins, which were the starting point for this hypothesis, rests on the idea that inverted erosional marks at the ice bed are subsequently infilled to form drumlins. Field tests on the sedimentology, architecture, and landform associations of drumlins in the Livingstone Lake drumlin field are outlined before more extensive work on bedrock erosional forms and flood routes is introduced. Bedrock erosional forms played a central part in establishing the hypothesis since their form and ornamentation are confidently interpreted as fluvial. Their form and genesis are discussed mainly with reference to sites at French River and Wilton Creek, Ontario, though some remarkable bedrock erosional forms in Antarctica support their regional extent. Initially in the meltwater hypothesis, drumlins were thought to be cavity fills and erosional drumlins were recognized later. This development is shown to be central to the realization that drumlin composition may be inferred from drumlin form. The scale of drumlin fields, measured at about 10³ km², and the magnitude of the inferred floods require that the flood events were regional. Regional-scale flood tracts in Ontario, Quebec, Alberta and the Northwest Territories extending over 1000 km in length and several hundred kilometers in width, support this suggestion. Floods, had they occurred, would have caused rapid rates of sea level rise and may have changed climate through their effects on ocean stratification and sea surface temperatures. The meltwater hypothesis covers a range of bedforms besides drumlins and bedrock erosional marks—fluting, Rogen moraine, hummocky terrain, and transverse ridges. Recent work shows how these forms are best explained by the meltwater hypothesis. The roles of water storage and release, which underpin the theory of the meltwater hypothesis, remain poorly understood.     

Drumlin formation: a time transgressive model

Graphical and numerical reconstructions of the Rainy and Superior lobes of the Laurentide Ice Sheet suggest that drumlin formation was time transgressive. Suites of glacial landforms including drumlins, tunnel valleys, eskers, and ice-collapse features can be correlated with specific recessional ice margins and are used as boundary conditions in the modeling. A contour map of the ice surface is then drawn using a specified basal shear stress. The shear stress can be constant or allowed to vary with position on the bed and is chosen to be consistent with the subglacial regime indicated by field evidence. Assuming that ice flow is parallel to drumlin orientations and perpendicular to the ice surface contours and moraines, the trend of drumlin axes is best accommodated by time transgressive drumlin formation during minor stillstands in the overall ice recession. The alternative, that drumlins were formed while the ice was at the Late Wisconsin maximum limit, requires large spatial variations in the basal shear stress distribution and therefore implies large mass-balance gradients or large variations in basal sliding velocities over small distances, for which there is little evidence.     

The internal structure of glacial landforms: an example from the Halton till plain, Scarborough Bluffs, Ontario

Current views on the internal structure of many glacial landforms need further definition. For example, drumlinized Halton till plain near the Scarborough Bluffs, Ontario would traditionally be. viewed as a lodgement till sheet, but it was found to consist of complex sedimentary assemblages including sediment flows, melt-out, deformation and lodgement tills. These facies vary spatially depending on whether deposition occurred beneath grounded ice or within subglacial cavities. Proglacial sediments bury portions of the till plain. Surface Rutings and drumlins clearly indicate the action of subglacial processes on the surface of Halton drift. Sedimentary structures at the contact between stratified sediments and diamictons within the Late Wisconsinan Halton drift are similar to those in older beds exposed at Scarborough Bluffs. The demonstration of the role of grounded ice in Halton drift and the similarity of sedimentary structures to those of the underlying Thorncliffe and Sunnybrook sediments suggests that the action of grounded ice cannot be ruled out in the case of the lower beds, as has been done by Eyles & Eyles ( Geology 11 , 146–152, 1983). Thus, surface Halton drift may be a model for recognition of similar environments of deposition in older beds beneath Halton. This analysis indicates flaws in a recent re-evaluation of Scarborough Bluffs sediment interpreted as solely lacustrine and not directly affected by glaciers.     

The last deglaciation of the Franz Victoria Trough, northern Barents Sea

A study of two piston cores and a 3.5 kHz seismic profile from the Franz Victoria Trough provides new stratigraphic, stable isotopic and foraminiferal AMS ¹⁴C data that help constrain the timing of ice-sheet retreat in the northern Barents Sea and the nature of the deglacial marine environment. Silty diamicton at the base of each core, interpreted as till or ice-marginal debris flow, suggests that the Barents ice sheet was grounded at the core sites (470 m water depth). Eight AMS ¹⁴C dates on sediment overlying the diamicton indicate that the ice sheet retreated from both core sites by 12.9 ka and that postglacial sedimentation began 10 ka ago. These dates, combined with a recently published ¹⁴C date from a nearby core, suggest that the Franz Victoria Trough may not have been deglaciated until c . 13 ka, 2000 years later than modeled ice-sheet reconstructions indicate. In the trough, oxygen isotopic ratios in planktonic foraminifera ^N. pachyderma (sinistral) were 0.5–0.750, lower during deglaciation than after, probably as a result of ice-sheet and/or iceberg melting. Foraminiferal assemblages suggest that Atlantic-derived intermediate water may have begun to penetrate the trough c . 13 ka ago.     

Regional reconstruction of subglacial hydrology and glaciodynamic behaviour along the southern margin of the Cordilleran Ice Sheet in British Columbia,Canada and northern Washington State,USA

Subglacial landsystems in and around Okanagan Valley, British Columbia, Canada are investigated in order to evaluate landscape development, subglacial hydrology and Cordilleran Ice Sheet dynamics along its southern margin. Major landscape elements include drumlin swarms and tunnel valleys. Drumlins are composed of bedrock, diamicton and glaciofluvial sediments; their form truncates the substrate. Tunnel valleys of various scales (km to 100s km length), incised into bedrock and sediment, exhibit convex longitudinal profiles, and truncate drumlin swarms. Okanagan Valley is the largest tunnel valley in the area and is eroded >300 m below sea level. Over 600 m of Late Wisconsin-age sediments, consisting of a fining-up sequence of cobble gravel, sand and silt fill Okanagan Valley. Landform–substrate relationships, landform associations, and sedimentary sequences are incompatible with prevailing explanations of landsystem development centred mainly on deforming beds. They are best explained by meltwater erosion and deposition during ice sheet underbursts.During the Late-Wisconsin glaciation, Okanagan Valley functioned as part of a subglacial lake spanning multiple connected valleys (few 100s km) of southern British Columbia. Subglacial lake development started either as glaciers advanced over a pre-existing sub-aerial lake (catch lake) or by incremental production and storage of basal meltwater. High geothermal heat flux, geothermal springs and/or subglacial volcanic eruptions contributed to ice melt, and may have triggered, along with priming from supraglacial lakes, subglacial lake drainage. During the underburst(s), sheetflows eroded drumlins in corridors and channelized flows eroded tunnel valleys. Progressive flow channelization focused flows toward major bedrock valleys. In Okanagan Valley, most of the pre-glacial and early-glacial sediment fill was removed. A fining-up sequence of boulder gravel and sand was deposited during waning stages of the underburst(s) and bedrock drumlins in Okanagan Valley were enhanced or wholly formed by this underburst(s).Subglacial lake development and drainage had an impact on ice sheet geometry and ice volumes. The prevailing conceptual model for growth and decay of the CIS suggests significantly thicker ice in valleys compared to plateaus. Subglacial lake development created a reversal of this ice sheet geometry where grounded ice on plateaus thickened while floating valley ice remained thinner (due to melting and enhanced sliding, with significant transfer of ice toward the ice sheet margin). Subglacial lake drainage may have hastened deglaciation by melting ice, lowering ice-surface elevations, and causing lid fracture. This paper highlights the importance of ice sheet hydrology: its control on ice flow dynamics, distribution and volume in continental ice masses.     

黄河内蒙古河段河道冲淤演变与凌情响应机制

黄河上游内蒙古冲积性河道凌汛问题突出,研究河道冲淤演变与凌情响应机制可为该河段防凌减灾提供技术支持。根据内蒙古河段凌情、河道冲淤演变资料,分析凌情变化表征指标及与之密切相关的河道冲淤演变特征指标,研究河道冲淤演变特征指标与凌情变化表征指标的响应关系。结果表明:河道冲淤演变的特征指标平滩流量与凌情表征指标冰下过流能力、槽蓄水增量关系密切,冰下过流能力为平滩流量的1/5左右,随着平滩流量减小而减小,而槽蓄水增量随着平滩流量的减小而增大,有利内蒙古河段防凌的平滩流量宜不小于2 000 m3/s,槽蓄水增量宜不超过14亿m3。本研究成果可为内蒙古河段冰凌灾害防治提供参考。     

Reconstructing past glacier dynamics and erosion from glacial geomorphic evidence: Snowdon,North Wales

Bedrock surfaces exposed around Llyn Llydaw, North Wales demonstrate contrasting styles of erosion beneath a Late Devensian ice sheet and a Loch Lomond Stadial (LLS) valley glacier. Ice sheet erosion involved lee-side fracturing, surface fracture wear and abrasive wear, while LLS erosion was primarily by abrasive wear. Preservation of ice sheet erosional features indicates limited rates of erosion during the LLS. Analysis of the geometry and distribution of erosional markings suggests that the low erosional capacity of the LLS glacier was due to a low basal sliding velocity. This prevented the formation of lee-side cavities, reduced the debris flux over the bed and minimised particle-bed contact loads. Reconstructions of the mass balance and geometry of the LLS glacier indicate that most of its balance velocity could be achieved by internal deformation alone. A combination of low subglacial water pressures and an unusually rough substrate explain the low sliding velocities. High bed roughness is due to the absence of leeside cavities and a change in flow orientation between ice sheet and LLS times, which meant that the LLS glacier was in contact with roughness elements which were generated in cavities beneath the ice sheet.     

Submarine sediment and landform record of a palaeo‐ice stream within the British−Irish Ice Sheet

This paper examines marine geophysical and geological data, and new multibeam bathymetry data to describe the Pleistocene sediment and landform record of a large ice‐stream system that drained ～3% of the entire British?Irish Ice Sheet at its maximum extent. Starting on the outer continental shelf NW of Scotland we describe: the ice‐stream terminus environment and depocentre on the outer shelf and continental slope; sediment architecture and subglacial landforms on the mid‐shelf and in a large marine embayment (the Minch); moraines and grounding line features on the inner shelf and in the fjordic zone. We identify new soft‐bed (sediment) and hard‐bed (bedrock) subglacial landform assemblages in the central and inner parts of the Minch that confirm the spatial distribution, coherence and trajectory of a grounded fast‐flowing ice‐sheet corridor. These include strongly streamlined bedrock forms and megagrooves indicating a high degree of ice‐bed coupling in a zone of flow convergence associated with ice‐stream onset; and a downstream bedform evolution (short drumlins to km‐scale glacial lineations) suggesting an ice‐flow velocity transition associated with a bed substrate and roughness change in the ice‐stream trunk. Chronology is still lacking for the timing of ice‐stream demise; however, the seismic stratigraphy, absence of moraines or grounding‐line features, and presence of well‐preserved subglacial bedforms and iceberg scours, combined with the landward deepening bathymetry, all suggest that frontal retreat in the Minch was probably rapid, via widespread calving, before stabilization in the nearshore zone. Large moraine complexes recording a coherent, apparently long‐lived, ice‐sheet margin position only 5–15 km offshore strongly support this model. Reconstructed ice‐discharge values for the Minch ice stream (12–20 Gt a^?1) are comparable to high mass‐flux ice streams today, underlining it as an excellent palaeo‐analogue for recent rapid change at the margins of the Greenland and West Antarctic Ice Sheets.     

Late Weichselian ice-sheet sensitivity over Franz Josef Land, Russian High Arctic, from numerical modelling experiments

A numerical ice-sheet model was run in order to produce reconstructions of the Late Weichselian ice coverage of Franz Josef Land, Russian High Arctic. The model grid covers the archipelago and surrounding shelf, but does not include the whole Barents-Kara region or the extensive ice cover that may have built up there. One experiment, where rates of iceberg calving at the grounded margin were curtailed because of the assumed presence of permanent thick sea ice, yielded a single I.8 km-thick ice dome which covered the entire archipelago and surrounding sea. If, however, iceberg calving were included in the model's environmental input, the extent of the ice sheet would be limited to the periphery of the archipelago. If a large ice sheet existed over Franz Josef Land, the deglaciation of the islands may have been linked to the decay of the adjacent Barents-Kara Sea Ice Sheet, permitting iceberg calving (enhanced by relative sea-level rise) to occur. The introduction of a water-depth-related iceberg calving function at 15 000 yr ago forced an initial rapid rate of ice-sheet decay of 30 000 km³ 1000 yr'. However, as the ice sheet thinned, and isostatic rebound began, the calculated rate of iceberg calving was reduced such that ice remained over the archipelago at 8000 yr ago. The model's failure to simulate complete ice-sheet decay by 8000 yr ago is at variance with radiocarbon-dated raised terraces on Franz Josef Land, which indicates the complete deglaciation of the islands at this time.     

Thermo-mechanical facies representative of fast and slow flowing ice sheets: the Weichselian ice sheet,a central west Poland case study

This study deals with an issue of thermo-mechanical facies, reflecting specific thermal and mechanical properties of the subglacial environment. The main objective of this study was to develop a model of glacitectionic deformation and its sedimentary record beneath fast and slow flowing ice sheets, based on investigations conducted in Wielkopolska (west central Poland). Sedimentary structures, mainly at the contact between subglacial tills and glacifluvial sediments, were recognized to delineate typical facies associations in a Weichselian glacigenic succession. Each association was interpreted as a record of the different depositional environments related to different subglacial conditions. Those investigations suggest the substratum was composed of frozen and dry, and wet and mobile spots, and four thermo-mechanical facies were distinguished: A – is representative of slower ice flow, dry and cold subglacial conditions, where driving stresses and normal effective pressure were high; B – is also related to slow ice flow and occurrence of cold subglacial permafrost, but with little amount of unfrozen water (however, higher than in facies A), with similar physical characteristics of the ice sheet as facies A; thermo-mechanical facies C and D represent wet and warm ice sole, with low normal effective pressure and driving stresses, thus lowering sediments’ shear strength and enabling high ice-flow velocities. We suggest that these facies have specific and non-random location, thereby revealing the relationship between subglacial thermo-mechanical conditions and ice sheet dynamics. Slow moving, cold-based ice occurred along ice sheet margins and inter-stream areas, whereas fast-moving, warm-based, well-lubricated ice, was typical of the axial parts of ice streams.     

Drainage beneath ice sheets: groundwater–channel coupling,and the origin of esker systems from former ice sheets

The nature of the drainage system beneath ice sheets is crucial to their dynamic behaviour but remains problematic. An experimentally based theory of coupling between groundwater and major channel systems is applied to the esker systems in the area occupied the last ice sheet in Europe, which we regard as a fossil imprint of major longitudinal drainage channels. We conclude that the large-scale distribution and spacing of major eskers is consistent with the theory of groundwater control, in which esker spacing is partly controlled by the transmissivity of the bed.It is concluded that esker patterns reflect the large-scale organisation of the subglacial drainage pattern in which channel development is coupled to groundwater flow and to the ice sheet's dynamic regime. The theory is then used to deduce: basal meltwater recharge rates and their spatial variability from esker spacing in an area in which the ice sheet was actively streaming during its final retreat; patterns of palaeo-groundwater flow and head distribution; and the seasonally varying magnitude of discharge from stream tunnels at the retreating ice sheet margin. Major channel/esker systems appear to have been stable at least over several hundred of years during the retreat of the ice sheet, although major dynamic events are demonstrably associated with major shifts in the hydraulic regime.Modelling suggests: that glaciation can stimulate deep groundwater circulation cells that are spatially linked to channel locations, with groundwater flow predominantly transverse to ice flow; that the circulation pattern has the potential to create large-scale anomalies in groundwater chemistry; and that the spacing of channels will change through the glacial cycle, influencing water pressures in stream tunnels, subglacial hydraulic gradients and effective pressure. If the latter is reduced sufficiently, it could trigger enhanced bed deformation, thus coupling drainage to ice sheet movement. It suggests the possibility of distinctive phases of sediment deformation and drumlin mobilisation during a glacial cycle.