skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on January 10, 2026

Title: Linear-viscous flow of temperate ice
Accurately modeling the deformation of temperate glacier ice, which is at its pressure-melting temperature and contains liquid water at grain boundaries, is essential for predicting ice sheet discharge to the ocean and associated sea-level rise. Central to such modeling is Glen’s flow law, in which strain rate depends on stress raised to a power ofn= 3 to 4. In sharp contrast to this nonlinearity, we found by conducting large-scale, shear-deformation experiments that temperate ice is linear-viscous (n ≈1.0) over common ranges of liquid water content and stress expected near glacier beds and in ice-stream margins. This linearity is likely caused by diffusive pressure melting and refreezing at grain boundaries and could help to stabilize modeled responses of ice sheets to shrinkage-induced stress increases.  more » « less
Award ID(s):
1643120
PAR ID:
10598064
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
AAAS
Date Published:
Journal Name:
Science
Volume:
387
Issue:
6730
ISSN:
0036-8075
Page Range / eLocation ID:
182 to 185
Subject(s) / Keyword(s):
Temperate ice linear-viscous flow ice streams
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ice at depth in ice-stream shear margins is thought to commonly be temperate, with interstitial meltwater that softens ice. Models that include this softening extrapolate results of a single experimental study in which ice effective viscosity decreased by a factor of ∼3 over water contents of ∼0.01–0.8%. Modeling indicates this softening by water localizes strain in shear margins and through shear heating increases meltwater at the bed, enhancing basal slip. To extend data to higher water contents, we shear lab-made ice in confined compression with a large ring-shear device. Ice rings with initial mean grain sizes of 2–4 mm are kept at the pressure-melting temperature and sheared at controlled rates with peak stresses of ∼0.06–0.20 MPa, spanning most of the estimated shear-stress range in West Antarctic shear margins. Final mean grain sizes are 8–13 mm. Water content is measured by inducing a freezing front at the ice-ring edges, tracking its movement inward with thermistors, and fitting the data with solutions of the relevant Stefan problem. Results indicate two creep regimes, below and above a water content of ∼0.6%. Comparison of effective viscosity values in secondary creep with those of tertiary creep from the earlier experimental study indicate that for water contents of 0.2–0.6%, viscosity in secondary creep is about twice as sensitive to water content than for ice sheared to tertiary creep. Above water contents of 0.6%, viscosity values in secondary creep are within 25% of those of tertiary creep, suggesting a stress-limiting mechanism at water contents greater than 0.6% that is insensitive to ice fabric development in tertiary creep. At water contents of ∼0.6–1.7%, effective viscosity is independent of water content, and ice is nearly linear-viscous. Minimization of intercrystalline stress heterogeneity by grain-scale melting and refreezing at rates that approach an upper bound as grain-boundary water films thicken might account for the two regimes. 
    more » « less
  2. Abstract Small quantities of liquid water lining triple junctions in polycrystalline glacier ice form connected vein networks that enable material exchange with underlying basal environments. Diffuse debris concentrations commonly observed in ice marginal regions might be attributed to this mechanism. Following recent cryogenic ring-shear experiments, we observed emplacement along grain boundaries of loess particles several tens of microns in size. Here, we describe an idealized model of vein liquid flow to elucidate conditions favoring such particle transport. Gradients in liquid potential drive flow toward colder temperatures and lower solute concentrations, while deviations of the ice stress state from hydrostatic balance produce additional suction toward anomalously low ice pressures. Our model predicts particle entrainment following both modest warming along the basal interface resulting from anticipated natural changes in effective stress, and the interior relaxation of temperature and solute concentration imposed by our experimental protocols. Comparisons with experimental observations are encouraging, but suggest that liquid flow rates are somewhat higher and/or more effective at dragging larger particles than predicted by our idealized model with nominal parameter choices. Diffuse debris entrainment extending several meters above the glacier bed likely requires a more sophisticated treatment that incorporates effects of ice deformation or other processes. 
    more » « less
  3. Abstract Results of ice-stream models that treat temperate ice deformation as a two-phase flow are sensitive to the ice permeability. We have constructed and begun using a custom, falling-head permeameter for measuring the permeability of temperate, polycrystalline ice. Chilled water is passed through an ice disk that is kept at the pressure-melting temperature while the rate of head decrease indicates the permeability. Fluorescein dye in the water allows water-vein geometry to be studied using fluorescence microscopy. Water flow over durations of seconds to hours is Darcian, and for grain diameter d increasing from 1.7 to 8.9 mm, average permeability decreases from 2 × 10 −12 to 4 × 10−15 m 2. In tests with dye on fine ( d= 2 mm) and coarse (d = 7 mm) ice, average area-weighted vein radii are nearly equal, 41 and 34 μm, respectively. These average radii, if included in a theory slightly modified from Nye and Frank (1973), yield permeability values within a factor of 2.0 of best-fit values based on regression of the data. Permeability values depend on d −3.4, rather than d−2 as predicted by models if vein radii are considered independent of d. In future experiments, the dependence of permeability on liquid water content will be measured. 
    more » « less
  4. Abstract The importance of glacier sliding has motivated a rich literature describing the thermomechanical interactions between ice, liquid water and bed materials. Early recognition of the gradient in melting temperature across small bed obstacles led to focused studies of regelation. An appreciation for the limits on ice deformation rates downstream of larger obstacles highlighted a role for cavitation, which has subsequently gained prominence in descriptions of subglacial drainage. Here, we show that the changes in melting temperature that accompany changes in normal stress along a sliding ice interface near cavities and other macroscopic drainage elements cause appreciable supercooling and basal mass exchange. This provides the basis of a novel formation mechanism for widely observed laminated debris-rich basal ice layers. 
    more » « less
  5. Abstract For brittle friction and rock deformation, the coefficientαin the general effective stress relationσe = σ − αPpcan be approximated as unity with sufficient accuracy. However, it is uncertain ifαdeviates from unity for semibrittle flow when both brittle and intracrystalline‐plastic deformation is involved. We conducted triaxial and isostatic compression experiments on synthetic salt‐rocks (∼300 ppm water) at room temperature to test the effective stress relation in the semibrittle regime using silicone oil and argon gas as pore fluids. Confining and pore pressures were cycled while their difference (differential pressure) was kept constant, such that changes in the mechanical behavior would indicate deviation ofαfrom unity. Microstructural observations were used to determine the dependence ofαon true area of grain contact from asperity yielding. In triaxial compression experiments, semibrittle flow involves grain boundary cracking and sliding, and intragranular dislocation glide and cracking. Flow strength remains constant for changes in pore fluid pressure of more than two orders of magnitude. In isostatic compression experiments, samples show combined processes of microcracking, grain boundary sliding, dislocation glide, and fluid‐assisted grain boundary migration recrystallization. Volumetric strain depends directly on the differential pressures (i.e.,αequals one). Analysis of grain‐contact area in both experiments indicates thatαis independent of the true area of contact defined by plastic yielding at grain boundaries. The observation ofαeffectively equals one may be explained by operation of pressure‐independent intracrystalline‐plastic mechanisms and transmission of pore pressure at grain boundaries through thin fluid films. 
    more » « less