skip to main content


Title: Freeze-on limits bed strength beneath sliding glaciers
Abstract

Discharge from sliding outlet glaciers controls uncertainty in projections for future sea level. Remarkably, over 90% of glacial area is subject to gravitational driving stresses below 150 kPa (median ∼70 kPa). Longstanding explanations that appeal to the shear-thinning rheology of ice tend to overpredict driving stresses and are restricted to areas where ice sheets only deform (roughly 50%). Over the more dynamic portions that slide, driving stresses must be balanced by thermo-mechanical interactions that control basal strength. Here we show that median bed strength is comparable to a threshold effective stress set by ice–liquid surface energy and till pore size. Above this threshold, ice infiltrates sediment to produce basal layers of debris-rich ice, even where net melting takes place. We demonstrate that the narrow range of inferred bed strengths can be explained by the mechanical resistance to sliding where roughness is enhanced by heterogeneous freeze-on.

 
more » « less
NSF-PAR ID:
10153576
Author(s) / Creator(s):
; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
9
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Globally, glaciers are shrinking in response to climate change, with implications for global sea level rise as well as downstream ecosystems and water resources. Sliding at the ice‐bed interface (basal motion) provides a mechanism for glaciers to respond rapidly to climate change. While the short‐term dynamics of glacier basal motion (<10 years) have received substantial attention, little is known about how basal motion and its sensitivity to subglacial hydrology changes over long (>50 year) timescales—this knowledge is required for accurate prediction of future glacier change. We compare historical data with modern estimates from field and satellite data at Athabasca Glacier and show that the glacier thinned by 60 m (−21%) over 1961–2020. However, a concurrent increase in surface slope results in minimal change in the average driving stress (−6 kPa and −4%). These geometric changes coincide with relatively uniform slowing (−15 m a−1and −45%). Simplified ice modeling suggests that declining basal motion accounts for most of this slow down (91% on average and 46% at minimum). A decline in basal motion can be explained by increasing basal friction resulting from geometric change in addition to increasing meltwater flux through a more efficient subglacial hydrologic system. These results highlight the need to include time‐varying dynamics of basal motion in glacier models and analyses. If these findings are generalizable, they suggest that declining basal motion reduces the flux of ice to lower elevations, helping to mitigate glacier mass loss in a warming climate.

     
    more » « less
  2. Abstract

    Glacier sliding has major environmental consequences, but friction caused by debris in the basal ice of glaciers is seldom considered in sliding models. To include such friction, divergent hypotheses for clast‐bed contact forces require testing. In experiments we rotate an ice ring (outside diameter = 0.9 m), with and without isolated till clasts, over a smooth rock bed. Ice is kept at its pressure‐melting temperature, and meltwater drains along a film at the bed to atmospheric pressure at its edges. The ice pressure or bed‐normal component of ice velocity is controlled, while bed shear stress is measured. Results with debris‐free ice indicate friction coefficients < 0.01. Shear stresses caused by clasts in ice are independent of ice pressure. This independence indicates that with increases in ice pressure the water pressure in cavities observed beneath clasts increases commensurately to allow drainage of cavities into the melt film, leaving clast‐bed contact forces unaffected. Shear stresses, instead, are proportional to bed‐normal ice velocity. Cavities and the absence of regelation ice indicate that, unlike model formulations, regelation past clasts does not control contact forces. Alternatively, heat from the bed melts ice above clasts, creating pressure gradients in adjacent meltwater films that cause contact forces to depend on bed‐normal ice velocity. This model can account for observations if rock friction predicated on Hertzian clast‐bed contacts is assumed. Including debris‐bed friction in glacier sliding models will require coupling the ice velocity field near the bed to contact forces rather than imposing a pressure‐based friction rule.

     
    more » « less
  3. Slip of marine-terminating ice streams over beds of deformable till is responsible for most of the contribution of the West Antarctic Ice Sheet to sea level rise. Flow models of the ice sheet and till-bedded glaciers elsewhere require a law that relates slip resistance, slip velocity, and water pressure at the bed. We present results of experiments in which pressurized ice at its melting temperature is slid over a water-saturated till bed. Steady-state slip resistance increases with slip velocity owing to sliding of ice across the bed, but above a threshold velocity, till shears at its rate-independent Coulomb strength. These results motivate a generalized slip law for glacier-flow models that combines processes of hard-bedded sliding and bed deformation.

     
    more » « less
  4. Basal slip along glaciers and ice streams can be significantly modified by external time-dependent forcing, although it is not clear why some systems are more sensitive to tidal stresses. We have conducted a series of laboratory experiments to explore the effect of time varying load point velocity on ice-on-rock friction. Varying the load point velocity induces shear stress forcing, making this an analogous simulation of aspects of ice stream tidal modulation. Ambient pressure, double-direct shear experiments were conducted in a cryogenic servo-controlled biaxial deformation apparatus at temperatures between −2°C and −16°C. In addition to a background, median velocity (1 and 10 μm/s), a sinusoidal velocity was applied to the central sliding sample over a range of periods and amplitudes. Normal stress was held constant over each run (0.1, 0.5 or 1 MPa) and the shear stress was measured. Over the range of parameters studied, the full spectrum of slip behavior from creeping to slow-slip to stick-slip was observed, similar to the diversity of sliding styles observed in Antarctic and Greenland ice streams. Under conditions in which the amplitude of oscillation is equal to the median velocity, significant healing occurs as velocity approaches zero, causing a high-amplitude change in friction. The amplitude of the event increases with increasing period (i.e. hold time). At high normal stress, velocity oscillations force an otherwise stable system to behave unstably, with consistently-timed events during every cycle. Rate-state friction parameters determined from velocity steps show that the ice-rock interface is velocity strengthening. A companion paper describes a method of analyzing the oscillatory data directly. Forward modeling of a sinusoidally-driven slider block, using rate-and-state dependent friction formulation and experimentally derived parameters, successfully predicts the experimental output in all but a few cases. 
    more » « less
  5. Abstract

    Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica.

     
    more » « less