Recent developments in modeling ice sheet deformation

Author(s):  
Felicity McCormack ◽  
Roland Warner ◽  
Adam Treverrow ◽  
Helene Seroussi

<p>Viscous deformation is the main process controlling ice flow in ice shelves and in slow-moving regions of polar ice sheets where ice is frozen to the bed. However, the role of deformation in flow in ice streams and fast-flowing regions is typically poorly represented in ice sheet models due to a major limitation in the current standard flow relation used in most large-scale ice sheet models – the Glen flow relation – which does not capture the steady-state flow of anisotropic ice that prevails in polar ice sheets. Here, we highlight recent advances in modeling deformation in the Ice Sheet System Model using the ESTAR (empirical, scalar, tertiary, anisotropic regime) flow relation – a new description of deformation that takes into account the impact of different types of stresses on the deformation rate. We contrast the influence of the ESTAR and Glen flow relations on the role of deformation in the dynamics of Thwaites Glacier, West Antarctica, using diagnostic simulations. We find key differences in: (1) the slow-flowing interior of the catchment where the unenhanced Glen flow relation simulates unphysical basal sliding; (2) over the floating Thwaites Glacier Tongue where the ESTAR flow relation outperforms the Glen flow relation in accounting for tertiary creep and the spatial differences in deformation rates inherent to ice anisotropy; and (3) in the grounded region within 80km of the grounding line where the ESTAR flow relation locally predicts up to three times more vertical shear deformation than the unenhanced Glen flow relation, from a combination of enhanced vertical shear flow and differences in the distribution of basal shear stresses. More broadly on grounded ice, the membrane stresses are found to play a key role in the patterns in basal shear stresses and the balance between basal shear stresses and gravitational forces simulated by each of the ESTAR and Glen flow relations. Our results have implications for the suitability of ice flow relations used to constrain uncertainty in reconstructions and projections of global sea levels, warranting further investigation into using the ESTAR flow relation in transient simulations of glacier and ice sheet dynamics. We conclude by discussing how geophysical data might be used to provide insight into the relationship between ice flow processes as captured by the ESTAR flow relation and ice fabric anisotropy.</p>

1999 ◽  
Vol 28 ◽  
pp. 23-32 ◽  
Author(s):  
Chris D. Clark

AbstractSubglacially-produced drift lineations provide spatially extensive evidence of ice flow that can be used to aid reconstructions of the evolution of former ice sheets. Such reconstructions, however, are highly sensitive to assumptions made about the glaciodynamic context of lineament generation; when during the glacial cycle and where within the ice sheet were they produced. A range of glaciodynamic contexts are explored which include: sheet-flow submarginally restricted; sheet-flow pervasive; sheet- flow patch; ice stream; and surge or re-advance. Examples of each are provided. The crux of deciphering the appropriate context is whether lineations were laid down time-trans-gressively or isochronously. It is proposed that spatial and morphometric characteristics of lineations, and their association with other landforms, can be used as objective criteria to help distinguish between these cases.A logically complete ice-sheet reconstruction must also account for the observed patches of older lineations and other relict surfaces and deposits that have survived erasure by subsequent ice flow. A range of potential preservation mechanisms are explored, including: cold- based ice; low basal-shear stresses; shallowing of the deforming layer; and basal uncoupling.


2015 ◽  
Vol 61 (227) ◽  
pp. 537-550 ◽  
Author(s):  
Joseph H. Kennedy ◽  
Erin C. Pettit

AbstractThe observable microstructures in ice are the result of many dynamic and competing processes. These processes are influenced by climate variables in the firn. Layers deposited in different climate regimes may show variations in fabric which can persist deep into the ice sheet; fabric may ‘remember’ these past climate regimes. We model the evolution of fabric variations below the firn–ice transition and show that the addition of shear to compressive-stress regimes preserves the modeled fabric variations longer than compression-only regimes, because shear drives a positive feedback between crystal rotation and deformation. Even without shear, the modeled ice retains memory of the fabric variation for 200 ka in typical polar ice-sheet conditions. Our model shows that temperature affects how long the fabric variation is preserved, but only affects the strain-integrated fabric evolution profile when comparing results straddling the thermal-activation-energy threshold (∼−10°C). Even at high temperatures, migration recrystallization does not eliminate the modeled fabric’s memory under most conditions. High levels of nearest-neighbor interactions will, however, eliminate the modeled fabric’s memory more quickly than low levels of nearest-neighbor interactions. Ultimately, our model predicts that fabrics will retain memory of past climatic variations when subject to a wide variety of conditions found in polar ice sheets.


1976 ◽  
Vol 16 (74) ◽  
pp. 41-71 ◽  
Author(s):  
T.J. Hughes

Abstract Application of thermal convection theory to polar ice sheets (Hughes, 1970, 1971. 1972[a],[c]) is reviewed and expanded. If it occurs, thermal convection is mainly concentrated near the bed of the ice sheet; resulting in active and passive convective flow, respectively below and above the ice density inversion. Convection begins as transient creep when a stress-independent critical Rayleigh number is exceeded, and stabilizes as steady-state creep when a stress-dependent critical Rayleigh number is exceeded. Transient- creep convection begins as unstable ripples in isotherms near the bed, with some ripples becoming upward bulges of basal ice which rapidly shrink laterally and grow vertically to become ascending dikes of recrystallized basal ice during steady-state creep. Sills of basal ice are injected horizontally between weakly coupled layers in the strata of cold ice slowly sinking en masse between dikes. Convection begins under domes of thick ice toward the ice-sheet center and a stable polygonal array of dikes may form if frictional heat creates hot ice at the bed as rapidly as convection flow redistributes hot basal ice in dikes and sills, Advective flow transports the converting ice toward the margin of the ice sheet where dikes converge at the heads of ice streams. Dike—sill convection then becomes ice-stream convection in which the entire ice stream behaves like a dike, uncoupling from the bed, and rising en masse. This would help explain why ice streams flow at surge velocities.


2003 ◽  
Vol 37 ◽  
pp. 108-112 ◽  
Author(s):  
Tim H. Jacka ◽  
Shavawn Donoghue ◽  
Jun Li ◽  
William F. Budd ◽  
Ross M. Anderson

AbstractIce-sheet basal ice is warmer than that above because of the heat from the Earth’s interior. The stresses acting on the basal ice are greatest. In addition, the basal ice often contains debris consisting of silt and small stones picked up from the rock over which the ice flows. Because the base is the warmest part of an ice sheet and the stress there is greatest, flow rates in the basal ice are large and often contribute most of the ice movement. It is therefore important, for accurate modelling of the ice sheets, to know whether the debris within the basal ice enhances or retards the flow of the ice. In this paper, we describe laboratory deformation tests in uniaxial compression and in simple shear, on sand-laden ice. We find no significant dependence of flow rate on sand content (up to 15% volume) in the stress range 0.13–0.5 MPa and temperature range –0.02 to –18.0°C. Further work needs to include laboratory tests on debris-laden ice extracted from the polar ice sheets. This work is underway.


2016 ◽  
Vol 790 ◽  
Author(s):  
Katarzyna N. Kowal ◽  
Samuel S. Pegler ◽  
M. Grae Worster

We present an experimental and theoretical study of the dynamics of laterally confined marine ice sheets in the natural limit in which the long, narrow channel into which they flow is wider than the depth of the ice. A marine ice sheet comprises a grounded ice sheet in contact with bedrock that floats away from the bedrock at a ‘grounding line’ to form a floating ice shelf. We model the grounded ice sheet as a viscous gravity current resisted dominantly by vertical shear stresses owing to the no-slip boundary condition applied at the bedrock. We model the ice shelf as a floating viscous current resisted dominantly by horizontal shear stresses owing to no-slip boundary conditions applied at the sidewalls of the channel. The two shear-dominated regions are coupled by jump conditions relating force and fluid flux across a short transition region downstream of the grounding line. We find that the influence of the stresses within the transition region becomes negligible at long times and we model the transition region as a singular interface across which the ice thickness and mass flux can be discontinuous. The confined shelf buttresses the sheet, causing the grounding line to advance more than it would otherwise. In the case that the sheet flows on a base of uniform slope, we find asymptotically that the grounding line advances indefinitely as $t^{1/3}$, where $t$ is time. This contrasts with the two-dimensional counterpart, for which the shelf provides no buttressing and the grounding line reaches a steady state (Robison, J. Fluid Mech., vol. 648, 2010, pp. 363–380).


Science ◽  
2019 ◽  
Vol 364 (6444) ◽  
pp. eaav7908 ◽  
Author(s):  
E. Larour ◽  
H. Seroussi ◽  
S. Adhikari ◽  
E. Ivins ◽  
L. Caron ◽  
...  

Geodetic investigations of crustal motions in the Amundsen Sea sector of West Antarctica and models of ice-sheet evolution in the past 10,000 years have recently highlighted the stabilizing role of solid-Earth uplift on polar ice sheets. One critical aspect, however, that has not been assessed is the impact of short-wavelength uplift generated by the solid-Earth response to unloading over short time scales close to ice-sheet grounding lines (areas where the ice becomes afloat). Here, we present a new global simulation of Antarctic evolution at high spatiotemporal resolution that captures all solid Earth processes that affect ice sheets and show a projected negative feedback in grounding line migration of 38% for Thwaites Glacier 350 years in the future, or 26.8% reduction in corresponding sea-level contribution.


1979 ◽  
Vol 23 (89) ◽  
pp. 425-428
Author(s):  
T. J. Hughes

Abstract Most numerical models of present ice-sheet dynamics predict basal thermal conditions for an assumed geothermal heat flux and measured ice thickness, surface temperature, and snow precipitation. These models are not ideally suited for reconstructing former ice sheets because what is known for present ice sheets is unknown for former ones, and vice versa. In particular, geothermal heat fluxes are immeasurable at an ice-sheet bed but can be measured after the ice sheet is gone, and the thermal conditions predicted at an ice-sheet bed can be inferred from the glacial-geological–topographic record after the ice sheet is gone. The Maine CLIMAP ice-sheet reconstruction model uses these inferred basal thermal conditions to compute ice thicknesses from basal shear stresses. Basal shear stress is assumed to reflect the degree of ice–bed coupling which, in turn, is assumed to reflect the amount and distribution of basal water under the ice sheet. Under the ice-sheet interior, basal water exists in a thin film of constant thickness covering the low places on the bed. This film expands for a melting bed and contracts for a freezing bed. Along the ice-sheet margin, basal water exists in narrow channels of varying thickness corresponding to troughs on the bed. These water channels become deeper for a melting bed and shallower for a freezing bed. In areas covered by the Laurentide and Scandinavian ice sheets, myriads of interconnected lakes in regions of greatest postglacial rebound are interpreted as evidence suggesting the interior basal water distribution, whereas eskers pointed toward terminal moraines and troughs across continental shelves are interpreted as evidence suggesting the basal water distribution toward the margins. Continental-shelf troughs were assumed to correspond to former ice streams, by analogy with observations in Greenland and Antarctica. Three modes of glacial erosion are considered to be responsible for the lakes, eskers, troughs, and associated topography. Quarrying is by a freeze-thaw mechanism which occurs where the melting-point isotherm intersects bedrock, so it is important only for freezing or melting beds because high places on the bed are frozen, low places are melted, and minor basal temperature fluctuations shift the isotherm separating them. Crushing results when rocks at the ice-bed interface are ground against each other and the bed by glacial sliding, so it occurs where the bed is melted and is most important when the entire bed is melted. Abrasion of bedrock occurs when rock cutting tools imbedded in the ice at the ice–rock interface are moved across the interface by glacial sliding, so it is also most important when the entire bed is melted. If basal melting continued after the entire bed is melted, abrasion-rates drop because the basal water layer thickens and drowns bedrock projections otherwise subjected to abrasion. Basal freezing reduces both crushing and abrasion-rates by coating quarried rocks with a sheath of relatively soft ice and transporting them upward from the ice–rock interface. An initially flat subglacial topography will develop depressions where glacial erosion is greatest and deposition is least, and ridges where the opposite conditions prevail. We interpret the central depressions represented today by Hudson Bay and the Gulf of Bothnia as caused by erosion on a melting bed under the Laurentide and Scandinavian ice sheets, respectively. The arc of lakes, gulfs, and shallow seas surrounding these depressions are interpreted as resulting from a freezing bed under the former ice sheets. The present watershed separating the depressions from the arcs marks the approximate former basal equilibrium line where the bed was melted. The Canadian and Baltic continental shields beyond these arcs are blanketed by material eroded from within the arcs, and represent areas having a frozen bed where evidence for abrasion is missing and a second zone having a melting bed where evidence for abrasion is present. This basic pattern was assumed to be imprinted on the bed during the steady-state period of maximum ice-sheet extent, and maintained in varying degrees during growth and shrinkage of these ice sheets.


1979 ◽  
Vol 23 (89) ◽  
pp. 425-428
Author(s):  
T. J. Hughes

AbstractMost numerical models of present ice-sheet dynamics predict basal thermal conditions for an assumed geothermal heat flux and measured ice thickness, surface temperature, and snow precipitation. These models are not ideally suited for reconstructing former ice sheets because what is known for present ice sheets is unknown for former ones, and vice versa. In particular, geothermal heat fluxes are immeasurable at an ice-sheet bed but can be measured after the ice sheet is gone, and the thermal conditions predicted at an ice-sheet bed can be inferred from the glacial-geological–topographic record after the ice sheet is gone. The Maine CLIMAP ice-sheet reconstruction model uses these inferred basal thermal conditions to compute ice thicknesses from basal shear stresses.Basal shear stress is assumed to reflect the degree of ice–bed coupling which, in turn, is assumed to reflect the amount and distribution of basal water under the ice sheet. Under the ice-sheet interior, basal water exists in a thin film of constant thickness covering the low places on the bed. This film expands for a melting bed and contracts for a freezing bed. Along the ice-sheet margin, basal water exists in narrow channels of varying thickness corresponding to troughs on the bed. These water channels become deeper for a melting bed and shallower for a freezing bed. In areas covered by the Laurentide and Scandinavian ice sheets, myriads of interconnected lakes in regions of greatest postglacial rebound are interpreted as evidence suggesting the interior basal water distribution, whereas eskers pointed toward terminal moraines and troughs across continental shelves are interpreted as evidence suggesting the basal water distribution toward the margins. Continental-shelf troughs were assumed to correspond to former ice streams, by analogy with observations in Greenland and Antarctica.Three modes of glacial erosion are considered to be responsible for the lakes, eskers, troughs, and associated topography. Quarrying is by a freeze-thaw mechanism which occurs where the melting-point isotherm intersects bedrock, so it is important only for freezing or melting beds because high places on the bed are frozen, low places are melted, and minor basal temperature fluctuations shift the isotherm separating them. Crushing results when rocks at the ice-bed interface are ground against each other and the bed by glacial sliding, so it occurs where the bed is melted and is most important when the entire bed is melted. Abrasion of bedrock occurs when rock cutting tools imbedded in the ice at the ice–rock interface are moved across the interface by glacial sliding, so it is also most important when the entire bed is melted. If basal melting continued after the entire bed is melted, abrasion-rates drop because the basal water layer thickens and drowns bedrock projections otherwise subjected to abrasion. Basal freezing reduces both crushing and abrasion-rates by coating quarried rocks with a sheath of relatively soft ice and transporting them upward from the ice–rock interface.An initially flat subglacial topography will develop depressions where glacial erosion is greatest and deposition is least, and ridges where the opposite conditions prevail. We interpret the central depressions represented today by Hudson Bay and the Gulf of Bothnia as caused by erosion on a melting bed under the Laurentide and Scandinavian ice sheets, respectively. The arc of lakes, gulfs, and shallow seas surrounding these depressions are interpreted as resulting from a freezing bed under the former ice sheets. The present watershed separating the depressions from the arcs marks the approximate former basal equilibrium line where the bed was melted. The Canadian and Baltic continental shields beyond these arcs are blanketed by material eroded from within the arcs, and represent areas having a frozen bed where evidence for abrasion is missing and a second zone having a melting bed where evidence for abrasion is present. This basic pattern was assumed to be imprinted on the bed during the steady-state period of maximum ice-sheet extent, and maintained in varying degrees during growth and shrinkage of these ice sheets.


2020 ◽  
Vol 12 (24) ◽  
pp. 4178
Author(s):  
Xiaoli Su ◽  
Junyi Guo ◽  
C. K. Shum ◽  
Zhicai Luo ◽  
Yu Zhang

Replacing estimates of C20 from the Gravity Recovery and Climate Experiment (GRACE) monthly gravity field solutions by those from satellite laser ranging (SLR) data and including degree one terms has become a standard procedure for proper science applications in the satellite gravimetry community. Here, we assess the impact of degree one terms, SLR-based C20 and C30 estimates on GRACE-derived polar ice sheet mass variations. We report that degree one terms recommended for GRACE Release 06 (RL06) data have an impact of 2.5 times more than those for GRACE RL05 data on the mass trend estimates over the Greenland and the Antarctic ice sheets. The latest recommended C20 solutions in GRACE Technical Note 14 (TN14) affect the mass trend estimates of ice sheets in absolute value by more than 50%, as compared to those in TN11 and TN07. The SLR-based C30 replacement has some impact on the Antarctic ice sheet mass variations, mainly depending on the length of the study period. This study emphasizes that reliable solutions of low degree spherical harmonics are crucial for accurately deriving ice sheet mass balance from satellite gravimetry.


1997 ◽  
Vol 25 ◽  
pp. 418-422 ◽  
Author(s):  
Eric J. Steig

An important component of models of the cryosphere is the calculation of accumulation rates over polar ice sheets. As a first-order approximation, many models rely on the assumption that temperature is the main controlling factor for precipitation. However, compilation of available ice-core data, including a new core from Taylor Dome, East Antarctica, suggests that precipitation is significantly decoupled from temperature for a large proportion of both the Greenland and Antarctic ice sheets. While the estimated glacial-to-interglacial change in temperature does not differ greatly among ice cores from each ice sheet, the estimated change in accumulation rate varies by more than a factor of 2. A simple vapor-pressure parameterization gives reasonable estimates of accumulation in the ice-sheet interior, but this is not necessarily the case close to the ice-sheet margin, where synoptic weather systems are important.


Sign in / Sign up

Export Citation Format

Share Document