Modern Ice Sheets Research Articles

Lessons learned from shallow subglacial bedrock drilling campaigns in Antarctica

Abstract We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

De‐Tuning Albedo Parameters in a Coupled Climate Ice Sheet Model to Simulate the North American Ice Sheet at the Last Glacial Maximum

AbstractThe Last Glacial Maximum extent of the North American Ice Sheets is well constrained empirically but has proven to be challenging to simulate with coupled Climate‐Ice Sheet models. Coupled Climate‐Ice Sheet models are often too computationally expensive to sufficiently explore uncertainty in input parameters, and it is unlikely that values calibrated to reproduce modern ice sheets will reproduce the known extent of the ice at the Last Glacial Maximum. To address this, we run an ensemble with a coupled Climate‐Ice Sheet model (FAMOUS‐ice), simulating the final stages of growth of the last North American Ice Sheets' maximum extent. Using this large ensemble approach, we explore the influence of numerous uncertain ice sheet albedo, ice sheet dynamics, atmospheric, and oceanic parameters on the ice sheet extent. We find that ice sheet albedo parameters determine the majority of uncertainty when simulating the Last Glacial Maximum North American Ice Sheets. Importantly, different albedo parameters are needed to produce a good match to the Last Glacial Maximum North American Ice Sheets than have previously been used to model the contemporary Greenland Ice Sheet due to differences in cloud cover over ablation zones. Thus, calibrating coupled climate‐ice sheet models on one ice sheet may produce strong biases when the model is applied to a new domain.

Patterns of valley incision beneath the Greenland Ice Sheet revealed using automated mapping and classification

The Greenland Ice Sheet covers an area of 1.7 million km2, equivalent to ∼79 % of the surface of Kalaallit Nunaat (Greenland) and ∼1.2 % of the Earth's land surface. The macro-scale geomorphology beneath the ice can provide a valuable record of past ice sheet behaviour, particularly during warm periods that may serve as analogues for present and future climates. However, despite extensive mapping of the landscape by airborne radar surveying, Greenland's subglacial geomorphology remains comparatively understudied. Here we construct an automated workflow to identify, extract, and quantify the morphology of valley cross-sectional profiles across Greenland, as observed in NASA Operation IceBridge radar data. We identify 5335 cross-sectional profiles and apply a supervised machine learning method to classify valleys based on their morphological similarity to those formed by glacial or fluvial incision elsewhere in the Northern Hemisphere. Approximately two thirds of the valleys are classified as ‘glacial’, some of which reflect active incision at the modern ice sheet margin, whereas others are situated beneath cold-based, slow-moving ice, indicating that they were incised under a different ice configuration earlier in Greenland's glacial history. The presence of ‘fluvial’ valleys in the low-lying interior of northern Greenland and in mountainous southern Greenland suggests parts of the inherited landscape formed under ice-free conditions during pre- or inter-glacial times have been preserved due to negligible long-term subglacial erosion rates. Some low-lying catchments show hallmarks of a combination of fluvial, glacial, and glacio-fluvial incision, hinting at complex interplays between valley-forming processes over the history of the Greenland Ice Sheet.

Climate Variations in the Past 250 Million Years and Contributing Factors

AbstractWe simulate climate variations in the past 250 million years (Myr), using the fully coupled Community Earth System Model version 1.2.2 (CESM1.2.2) with the Community Atmosphere Model version 4 (CAM4). Three groups of simulations are performed, each including 26 simulations, with a 10‐million‐year interval. The Control group is constrained by paleogeography, increasing solar radiation, and reconstructed global mean surface temperatures (GMSTs) by tuning CO2 concentrations. No ice sheets are prescribed for all simulations except for the pre‐industrial (PI) simulation in which modern geography, ice sheets and vegetation are used. Simulated zonal mean surface temperatures are always higher than those of proxy reconstructions in the tropics, but lower than those of proxy reconstructions at middle latitudes. The relative importance of individual contributing factors for surface temperature variations in the past 250 Myr is diagnosed, using the energy‐balance analysis. Results show that greenhouse gases are the major driver in regulating GMST variations, with a maximum contribution of 12.2°C. Varying surface albedo contributes to GMST variations by 3.3°C. Increasing solar radiation leads to GMST increases by 1.5°C. Cloud radiative effects have relatively weak impacts on GMST variations, less than ±0.8°C. For comparison, two groups of sensitivity simulations are performed. One group has the CO2 concentration fixed at 10 times the PI value, and the other group has fixed CO2 concentration of 10 times the PI value and fixed solar radiation at the present‐day value, showing that varying both paleogeography and solar constant and varying paleogeography alone result in GMST changes by 7.3°C and 5.6°C, respectively.

Evidence for an extensive ice shelf in northern Baffin Bay during the Last Glacial Maximum

The glaciological significance of ice shelves is relatively well established for the stability of modern ice sheets of Antarctica. Past ice shelves of the Arctic, however, are poorly documented while their role for the stability of former ice sheets remains mostly unknown. Here we present swath bathymetry data and seismostratigraphic profiles that reveal a large moraine system extending along the continental slope off Baffin Island, demonstrating that a 500-m thick ice shelf covered northern Baffin Bay during the last glacial episode. We suggest that this ice shelf had a profound impact on the stability of a series of major ice streams that drained the interior of the Laurentide, Innuitian and Greenland ice sheets. Climate warming and global sea-level rise in the early stage of deglaciation possibly contributed to a large-scale break-up of the ice shelf, which led to the destabilisation and reorganisation of tributary ice streams from these three ice sheets.

Mechanism for the subglacial formation of cryogenic brines

AbstractCryogenic brines are under-studied, despite the fact that they may contain information about past ice-sheet behavior. Cryogenic brines form through cryoconcentration of seawater, although the specific setting and mechanism of formation have been debated. Previous conceptual models of brine formation require seawater isolation from the ocean in a closed basin experiencing freezing. We propose instead that they may form in pore spaces of marine sediments subjected to repeat cycles of ice-sheet advance and retreat. During periods of basal freezing, cryoconcentration produces hypersaline brines which experience downward flow driven by unstable density stratification. Our advection-diffusion model of porewater chemistry evolution successfully recreates the porewater chemistry of two deep Antarctic cores containing cryogenic brines (AND-1B and AND-2A), suggesting that cryogenic brines can be formed through the repeated isolation and cryoconcentration of marine waters within subglacial sediment pore spaces of modern and past ice sheets.

Character of advance and retreat of the southwest sector of the British-Irish Ice Sheet during the last glaciation

Relict landforms and sediments across former glaciated settings provide information about ice-sheet dynamics and can contribute to the understanding of the behaviour of contemporary ice masses, for which observations are limited in spatial and temporal extent. In this study, we focus on the shelf offshore southwest Ireland, in the Celtic Sea, which was once occupied by the Irish Sea Ice Stream (ISIS), the largest ice stream draining the southern portion of the marine-terminating British-Irish Ice Sheet (BIIS). Newly acquired high-resolution multibeam echosounder, sub-bottom and core data enabled the investigation of the shelf geomorphology and of the sedimentology and chronology of glacial and glacimarine sediments. A suite of drumlins records ice sheet flow from the coastline towards the central part of the shelf in southwest Ireland. Pre-existing highs in the seafloor topography promoted the formation of arcuate and transverse landforms interpreted as a grounding-zone wedge and moraines and they document episodic retreat of the ISIS across this portion of the shelf. Observed lithofacies show consolidated subglacial till and laminated fine muds. The sediments provide evidence of ice grounded ca. 30 km off the south-west Irish coastline with subsequent deglaciation occurring under glacimarine conditions. These new data refine the current reconstructions of the dynamics of the southern BIIS. They reveal for the first time the interplay of marine- and land-based ice and the presence of grounded ice offshore SW Ireland. This study highlights the importance of high-resolution data in revealing palaeo-landscapes as valuable analogues to test possible scenarios of modern ice sheet changes.

Total isostatic response to the complete unloading of the Greenland and Antarctic Ice Sheets

The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.

Review article: Existing and potential evidence for Holocene grounding line retreat and readvance in Antarctica

Abstract. Widespread existing geological records from above the modern ice sheet surface and outboard of the current ice margin show that the Antarctic Ice Sheet (AIS) was much more extensive at the Last Glacial Maximum (∼ 20 ka) than at present. However, whether it was ever smaller than present during the last few millennia, and (if so) by how much, is known only for a few locations because direct evidence lies within or beneath the ice sheet, which is challenging to access. Here, we describe how retreat and readvance (henceforth “readvance”) of AIS grounding lines during the Holocene could be detected and quantified using subglacial bedrock, subglacial sediments, marine sediment cores, relative sea-level (RSL) records, geodetic observations, radar data, and ice cores. Of these, only subglacial bedrock and subglacial sediments can provide direct evidence for readvance. Marine archives are of limited utility because readvance commonly covers evidence of earlier retreat. Nevertheless, stratigraphic transitions documenting change in environment may provide support for direct evidence from subglacial records, as can the presence of transgressions in RSL records, and isostatic subsidence. With independent age control, ice structure revealed by radar can be used to infer past changes in ice flow and geometry, and therefore potential readvance. Since ice cores capture changes in surface mass balance, elevation, and atmospheric and oceanic circulation that are known to drive grounding line migration, they also have potential for identifying readvance. A multidisciplinary approach is likely to provide the strongest evidence for or against a smaller-than-present AIS in the Holocene.

Formation of murtoos by repeated flooding of ribbed bedforms along subglacial meltwater corridors

Fluctuations in meltwater discharge below modern glaciers and ice sheets due to diurnal, seasonal and long-term temperature variations are modulated by complex interactions between subglacial drainage, basal processes and bedform development. The bed of palaeo-ice sheets contains a variety of bedforms recording these modulations and provides an open window into the subglacial environment. Through the morphometric analysis of natural and experimental bedforms, respectively mapped along Scandinavian meltwater corridors and produced in a physical model simulating transitory subglacial water flow, we observe a morphological and genetic bedform continuum corresponding to the progressive transformation of ribbed bedforms into murtoos. Two alternating drainage configurations, related to repeated subglacial flooding events, are involved in this transformation: (i) significant meltwater discharge, high hydraulic connectivity and ice-bed decoupling during flooding events lead to hydraulic alteration of ribbed bedforms by erosion, sediment deposition and channel incision, while (ii) limited meltwater flow, low hydraulic connectivity and ice-bed recoupling that follow flooding events lead to their deformational reshaping into murtoos. The degree of transformation of ribbed bedforms into murtoos can be quantified by combining two dimensionless morphometric parameters (circularity and sinuosity) and provides a convenient proxy to constrain magnitudes, durations and/or frequencies of subglacial floods in palaeo-meltwater corridors.

Reviewing the elements of marine ice cliff instability

Antarctica’s ice sheets are the largest potential sea-level rise contributors, but projections of future sea-level rise yield wide ranges of estimates under different emission scenarios. An important factor in the variability of estimates is marine ice cliff instability (MICI). Inclusion of MICI yields the highest potential sea-level rise cases but also the largest uncertainty due to poor understanding of the factors that control it and the mechanisms of how it happens. Although evidence for MICI has been implied by paleo-ice sheet studies and observations of keel plough mark on sea-floor, recent statistical and modelling studies have suggested a lower magnitude of MICI effect on sea-level rise due to thinning of ice sheets and buttressing forces exerted on potentially failing cliffs. This paper reviews the factors that control MICI with the goal of identifying priorities for modern ice sheet studies to better bound the estimates.

Did Holocene climate changes drive West Antarctic grounding line retreat and readvance?

Abstract. Knowledge of past ice sheet configurations is useful for informing projections of future ice sheet dynamics and for calibrating ice sheet models. The topology of grounding line retreat in the Ross Sea sector of Antarctica has been much debated, but it has generally been assumed that the modern ice sheet is as small as it has been for more than 100 000 years (Conway et al., 1999; Lee et al., 2017; Lowry et al., 2019; McKay et al., 2016; Scherer et al., 1998). Recent findings suggest that the West Antarctic Ice Sheet (WAIS) grounding line retreated beyond its current location earlier in the Holocene and subsequently readvanced to reach its modern position (Bradley et al., 2015; Kingslake et al., 2018). Here, we further constrain the post-LGM (Last Glacial Maximum) grounding line retreat and readvance in the Ross Sea sector using a two-phase model of radiocarbon input and decay in subglacial sediments from six sub-ice sampling locations. In addition, we reinterpret high basal temperature gradients, measured previously at three sites in this region (Engelhardt, 2004), which we explain as resulting from recent ice shelf re-grounding accompanying grounding line readvance. At one location – Whillans Subglacial Lake (SLW) – for which a sediment porewater chemistry profile is known, we estimate the grounding line readvance by simulating ionic diffusion. Collectively, our analyses indicate that the grounding line retreated over SLW 4300-2500+1500 years ago, and over sites on Whillans Ice Stream (WIS), Kamb Ice Stream (KIS), and Bindschadler Ice Stream (BIS) 4700-2300+1500, 1800-700+2700, and 1700-600+2800 years ago, respectively. The grounding line only recently readvanced back over those sites 1100-100+200, 1500-200+500, 1000-300+200, and 800±100 years ago for SLW, WIS, KIS, and BIS, respectively. The timing of grounding line retreat coincided with a warm period in the mid-Holocene to late Holocene. Conversely, grounding line readvance is coincident with cooling climate in the last 1000–2000 years. Our estimates for the timing of grounding line retreat and readvance are also consistent with relatively low carbon-to-nitrogen ratios measured in our subglacial sediment samples (suggesting a marine source of organic matter) and with the lack of grounding zone wedges in front of modern grounding lines. Based on these results, we propose that the Siple Coast grounding line motions in the mid-Holocene to late Holocene were primarily driven by relatively modest changes in regional climate, rather than by ice sheet dynamics and glacioisostatic rebound, as hypothesized previously (Kingslake et al., 2018).

Deformation mechanisms and kinematics of a soft sedimentary bed beneath the Scandinavian Ice Sheet, north-central Poland, revealed by magnetic fabrics

Characterising deformable substrates of modern and past ice sheets is a key factor in understanding ice mass behaviour and dynamics. This is challenging when field-scale signatures of deformation in soft sedimentary beds are weakly developed or absent. To overcome these difficulties, in this study, we combine micro- and macroscale studies and quantitative analyses to decipher the nature of subglacial sediment deformation under a Pleistocene fast-flowing (streaming) ice sheet. We investigated a complex sedimentary succession at Nieszawa in north-central Poland. A multi-proxy approach involving till micromorphology in thin sections, ripple-cross lamination orientation, till fabric, clast microfabric and the anisotropy of magnetic susceptibility (AMS) data in a closely sampled vertical profile throughout the fluvial sand and till was used to reveal subglacial bed deformation mechanisms and kinematics. This records subglacial deformation of the sedimentary succession at Nieszawa under low basal shear stresses, where simple shear was the dominant bed deformation mechanism. Intergranular advection, decoupling at the ice-bed interface related to high subglacial water pressure, clast ploughing, and basal sliding also took place. The main finding of this study is that AMS reveals the shear strain development and ice flow direction in the structurally and texturally complex deposits under weak coupling conditions at the ice-bed interface beneath a fast flowing/streaming ice sheet, where strain and kinematic indicators other than AMS fabrics are vague or absent.

A new map of glacigenic features and glacial landsystems in central mainland Nunavut, Canada

A new digital map of glacial geomorphic features and interpreted glacial landsystems was produced for an area covering ~415 000 km2 in the Keewatin Sector of the Laurentide Ice Sheet (LIS) in Nunavut. The map integrates information from previous surficial geology maps and >14 000 field stations, and is significantly improved by the detailed inventory of ~152 000 glacigenic features using high‐resolution ArcticDEM data and Landsat 8 imagery. From this, we identify and map coherent patterns of landform development (landsystems) between the Manitoba border and the Arctic coast, many of which are entirely new and others that are significantly modified or updated. In particular, we recognize six separate ice streams, including one probable remnant ice stream, and we delineate numerous palimpsest streamlined landscapes with associated ice‐flow trends and relative ages. A continuum of relict terrains with varying basal ice thermal conditions is mapped for the first time in the ice divide migration zone between Baker Lake and Wager Bay. In addition, deglacial cold‐based retreat terrains and preserved warm‐based landscapes unaffected by younger glacial events have been identified. These new georeferenced, multi‐scale data sets and interpreted glacial landsystems provide a comprehensive framework to strengthen reconstructions of the glacial history and dynamics of one of the largest ice domes of the LIS, identify distinct glacial sediment transport paths for applications to mineral exploration, and test numerical modelling of the LIS in support of climate change studies and long‐term evolution of modern ice sheets.

A new glacial varve chronology along the southern Laurentide Ice Sheet that spans the Younger Dryas–Holocene boundary

Abstract Glacial varves can detail ice-margin positions and provide a proxy for meltwater discharge at resolutions comparable to those of the Greenland ice core archives, and thus they can be critical paleorecords for assessing the response of both ancient and modern ice sheets to climate change. Here we provide an ∼1500 yr varve chronology straddling the Younger Dryas (YD)–Holocene boundary (11.65 cal. kyr B.P.), the first such chronology in North America. The varves are from glacial Lake Agassiz (central North America). The chronology is pinned on accelerator mass spectrometry radiocarbon-dated terrestrial macrofossils at the base of a widespread red-clay bed deposited during flooding from the Lake Superior basin. We illustrate the utility of this record by examining ice-margin retreat and melting through the late Younger Dryas and across the Holocene boundary. The ice margin receded at a constant rate, not only during the late YD, but for at least 300 yr after the onset of the Holocene. In contrast, varve thicknesses increased at the boundary, and a moraine formed over a 50 yr period, perhaps in response to the warming climate. Our expectation is that this time series will continue to be developed, expanded, and refined because it promises to be the longest and most geographically extensive glacial varve data set in North America.

New constraints on the last deglaciation of the Cordilleran Ice Sheet in coastal Southeast Alaska

AbstractUnderstanding marine-terminating ice sheet response to past climate transitions provides valuable long-term context for observations of modern ice sheet change. Here, we reconstruct the last deglaciation of marine-terminating Cordilleran Ice Sheet (CIS) margins in Southeast Alaska and explore potential forcings of western CIS retreat. We combine 27 new cosmogenic 10Be exposure ages, 13 recently published 10Be ages, and 25 new 14C ages from raised marine sediments to constrain CIS recession. Retreat from the outer coast was underway by 17 ka, and the inner fjords and sounds were ice-free by 15 ka. After 15 ka, the western margin of the CIS became primarily land-terminating and alpine glaciers disappeared from the outer coast. Isolated alpine glaciers may have persisted in high inland peaks until the early Holocene. Our results suggest that the most rapid phase of CIS retreat along the Pacific coast occurred between ~17 and 15 ka. This retreat was likely driven by processes operating at the ice-ocean interface, including sea level rise and ocean warming. CIS recession after ~15 ka occurred during a time of climatic amelioration in this region, when both ocean and air temperatures increased. These data highlight the sensitivity of marine-terminating CIS regions to deglacial climate change.

Widespread early Holocene deglaciation, Washington Land, northwest Greenland

Constraining the Greenland Ice Sheet (GIS) response to past climate change is imperative for understanding modern ice sheet sensitivity to rising temperatures. In this study, we investigate the response of the GIS to warming during the last deglaciation across the currently ice-free area of Washington Land in northwest Greenland. We apply 10Be cosmogenic surface exposure dating to glacially-transported boulders deposited across the landscape, the interpretation of the dates offering insight into GIS retreat after the Last Glacial Maximum. We find that the northwestern sector of the GIS sustained full ice cover across Washington Land until ∼11 ka, then began thinning, with widespread ice retreat culminating by 8.5 ± 1.1 ka. The record of our new exposure ages indicates that deglaciation occurred 2–3 ka after abrupt warming at the end of the Younger Dryas cold period, which suggests a delayed response of the GIS to this regional warming.

Characteristics of the deforming bed: till properties on the deglaciated Antarctic continental shelf

ABSTRACTContemporary ice stream flow is directly linked to conditions at the ice/bed interface, yet this environment is logistically difficult to access. Instead, we investigate subglacial processes important for ice stream flow by studying tills on the deglaciated Antarctic continental shelf. We test currently-accepted hypotheses surrounding subglacial processes and till properties with a Ross Sea dataset. Till shear strengths indicate a continuum of simultaneous processes acting at the bed, rather than discrete ‘deformation’ and ‘lodgement’ end-members. We identify a threshold water content representing saturated pore spaces, leading to basal sliding and meltwater channelization. Based on observations of till properties relative to glacial landforms, we challenge the assumption that low shear strength is linked to intense deformation. Spatial variability in landform morphology reflects variability in deforming processes at the sub-ice stream scale and suggests a maximum deforming bed thickness of 2 m at the grounding line. Regional till properties generally correlate with seafloor geology and deglacial history; the western Ross Sea is characterized by higher and more variable shear strengths and water contents, while lower-shear strength till was preserved in the Eastern Basin. These observations inform till interpretation and provide context for deforming beds beneath the modern ice sheet and on glaciated continental shelves.

The effects of tunnel channel formation on the Green Bay Lobe, Wisconsin, USA

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.

Microscale evidence of liquefaction and its potential triggers during soft-bed deformation within subglacial traction tills

Published conceptual models argue that much of the forward motion of modern and ancient glaciers is accommodated by deformation of soft-sediments within the underlying bed. At a microscale this deformation results in the development of a range of ductile and brittle structures in water-saturated sediments as they accommodate the stresses being applied by the overriding glacier. Detailed micromorphological studies of subglacial traction tills reveal that these polydeformed sediments may also contain evidence of having undergone repeated phases of liquefaction followed by solid-state shear deformation. This spatially and temporally restricted liquefaction of subglacial traction tills lowers the shear strength of the sediment and promotes the formation of “transient mobile zones” within the bed, which accommodate the shear imposed by the overriding ice. This process of soft-bed sliding, alternating with bed deformation, facilitates glacier movement by way of ‘stick-slip’ events. The various controls on the slip events have previously been identified as: (i) the introduction of pressurised meltwater into the bed, a process limited by the porosity and permeability of the till; and (ii) pressurisation of porewater as a result of subglacial deformation; to which we include (iii) episodic liquefaction of water-saturated subglacial traction tills in response to glacier seismic activity (icequakes), which are increasingly being recognized as significant processes in modern glaciers and ice sheets. As liquefaction operates only in materials already at very low values of effective stress, its process-form signatures are likely indicative of glacier sub-marginal tills.