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
Ice viscosity is more sensitive to stress than commonly assumed
Abstract Accurate representation of the viscous flow of ice is fundamental to understanding glacier dynamics and projecting sea-level rise. Ice viscosity is often described by a simple but largely untested and uncalibrated constitutive relation, Glen’s Flow Law, wherein the rate of deformation is proportional to stress raised to the power n . The value n = 3 is commonly prescribed in ice-flow models, though observations and experiments support a range of values across stresses and temperatures found on Earth. Here, we leverage recent remotely-sensed observations of Antarctic ice shelves to show that Glen’s Flow Law approximates the viscous flow of ice with n = 4.1 ± 0.4 in fast-flowing areas. The viscosity and flow rate of ice are therefore more sensitive to changes in stress than most ice-flow models allow. By calibrating the governing equation of ice deformation, our result is a pathway towards improving projections of future glacier change.
more »
« less
- Award ID(s):
- 1853918
- PAR ID:
- 10417126
- Date Published:
- Journal Name:
- Communications Earth & Environment
- Volume:
- 3
- Issue:
- 1
- ISSN:
- 2662-4435
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Viscous flow in ice is often described by the Glen flow law – anon-Newtonian, power-law relationship between stress and strain rate with astress exponent n ∼ 3. The Glen law is attributed tograin-size-insensitive dislocation creep; however, laboratory and fieldstudies demonstrate that deformation in ice can be strongly dependent ongrain size. This has led to the hypothesis that at sufficiently lowstresses, ice flow is controlled by grain boundary sliding, which explicitly incorporates the grain size dependence of ice rheology. Experimental studiesfind that neither dislocation creep (n ∼ 4) nor grain boundarysliding (n ∼ 1.8) have stress exponents that match the value ofn ∼ 3 in the Glen law. Thus, although the Glen law provides anapproximate description of ice flow in glaciers and ice sheets, itsfunctional form is not explained by a single deformation mechanism. Here weseek to understand the origin of the n ∼ 3 dependence of theGlen law by using the “wattmeter” to model grain size evolution in ice.The wattmeter posits that grain size is controlled by a balance between themechanical work required for grain growth and dynamic grain size reduction.Using the wattmeter, we calculate grain size evolution in two end-membercases: (1) a 1-D shear zone and (2) as a function of depth within anice sheet. Calculated grain sizes match both laboratory data and ice coreobservations for the interior of ice sheets. Finally, we show thatvariations in grain size with deformation conditions result in an effectivestress exponent intermediate between grain boundary sliding and dislocationcreep, which is consistent with a value of n = 3 ± 0.5 over the rangeof strain rates found in most natural systems.more » « less
-
Abstract. Full-thickness crevasses can transport water from the glacier surface to the bedrock where high water pressures can open kilometre-long cracks along the basal interface, which can accelerate glacier flow. We present a first computational modelling study that describes time-dependent fracture propagation in an idealised glacier causing rapid supraglacial lake drainage. A novel two-scale numerical method is developed to capture the elastic and viscoelastic deformations of ice along with crevasse propagation. The fluid-conserving thermo–hydro–mechanical model incorporates turbulent fluid flow and accounts for melting and refreezing in fractures. Applying this model to observational data from a 2008 rapid-lake-drainage event indicates that viscous deformation exerts a much stronger control on hydrofracture propagation compared to thermal effects. This finding contradicts the conventional assumption that elastic deformation is adequate to describe fracture propagation in glaciers over short timescales (minutes to several hours) and instead demonstrates that viscous deformation must be considered to reproduce observations of lake drainage rates and local ice surface elevation changes. As supraglacial lakes continue expanding inland and as Greenland Ice Sheet temperatures become warmer than −8 °C, our results suggest rapid lake drainage events are likely to occur without refreezing, which has implications for the rate of sea level rise.more » « less
-
A quasi-one-dimensional ice mélange flow model based on continuum descriptions of granular materialsAbstract. Field and remote sensing studies suggest that ice mélange influences glacier-fjord systems by exerting stresses on glacier termini and releasing large amounts of freshwater into fjords. The broader impacts of ice mélange over long time scales are unknown, in part due to a lack of suitable ice mélange flow models. Previous efforts have included modifying existing viscous ice shelf models, despite the fact that ice mélange is fundamentally a granular material, and running computationally expensive discrete element simulations. Here, we draw on laboratory studies of granular materials, which exhibit viscous flow when stresses greatly exceed the yield point, plug flow when the stresses approach the yield point, and stress transfer via force chains. By implementing the nonlocal granular fluidity rheology into a depth- and width-integrated stress balance equation, we produce a numerical model of ice mélange flow that is consistent with our understanding of well-packed granular materials and that is suitable for long time-scale simulations. For parallel-sided fjords, the model exhibits two possible steady state solutions. When there is no calving of new icebergs or melting of previously calved icebergs, the ice mélange is pushed down fjord by the advancing glacier terminus, the velocity is constant along the length of the fjord, and the thickness profile is exponential. When calving and melting are included, the ice mélange evolves to another steady state in which its location is fixed relative to the fjord walls, the thickness profile is relatively steep, and the flow is extensional. For the latter case, the model predicts that the steady-state ice mélange buttressing force depends on the surface and basal melt rates through an inverse power law relationship, decays roughly exponentially with both fjord width and gradient in fjord width, and increases with the iceberg calving flux. The increase in buttressing force with the calving flux, which depends on glacier thickness, appears to occur more rapidly than the force required to prevent the capsize of full-glacier-thickness icebergs, suggesting that glaciers with high calving fluxes may be more strongly influenced by ice mélange than those with small fluxes.more » « less
-
A quasi-one-dimensional ice mélange flow model based on continuum descriptions of granular materialsField and remote sensing studies suggest that ice mélange influences glacier–fjord systems by exerting stresses on glacier termini and releasing large amounts of freshwater into fjords. The broader impacts of ice mélange over long timescales are unknown, in part due to a lack of suitable ice mélange flow models. Previous efforts have included modifying existing viscous ice shelf models, despite the fact that ice mélange is fundamentally a granular material, and running computationally expensive discrete element simulations. Here, we draw on laboratory studies of granular materials, which exhibit viscous flow when stresses greatly exceed the yield point, plug flow when the stresses approach the yield point, and exhibit stress transfer via force chains. By implementing the nonlocal granular fluidity rheology into a depth- and width-integrated stress balance equation, we produce a numerical model of ice mélange flow that is consistent with our understanding of well-packed granular materials and that is suitable for long-timescale simulations. For parallel-sided fjords, the model exhibits two possible steady-state solutions. When there is no calving of icebergs or melting of previously calved icebergs, the ice mélange is pushed down-fjord by the advancing glacier terminus, the velocity is constant along the length of the fjord, and the thickness profile is exponential. When calving and melting are included and treated as constants, the ice mélange evolves into another steady state in which its location is fixed relative to the fjord walls, the thickness profile is relatively steep, and the flow is extensional. For the latter case, the model predicts that the steady-state ice mélange buttressing force depends on the surface and basal melt rates through an inverse power-law relationship, decays roughly exponentially with both fjord width and gradient in fjord width, and increases with the iceberg calving flux. The buttressing force appears to increase with calving flux (i.e., glacier thickness) more rapidly than the force required to prevent the capsizing of full-glacier-thickness icebergs, suggesting that glaciers with high calving fluxes may be more strongly influenced by ice mélange than those with small fluxes.more » « less
An official website of the United States government

