skip to main content


Title: Treatment of ice-shelf evolution combining flow and flexure
Abstract We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others.  more » « less
Award ID(s):
1841467 1841607
PAR ID:
10322206
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Glaciology
Volume:
67
Issue:
265
ISSN:
0022-1430
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Ice‐shelf basal channels form due to concentrated submarine melting. They are present in many Antarctic ice shelves and can reduce ice‐shelf structural integrity, potentially destabilizing ice shelves by full‐depth incision. Here, we describe the viscous ice response to a basal channel—secondary flow—which acts perpendicular to the channel axis and is induced by gradients in ice thickness. We use a full‐Stokes ice‐flow model to systematically assess the transient evolution of a basal channel in the presence of melting. Secondary flow increases with channel size and reduces the rate of channel incision, such that linear extrapolation or the Shallow‐Shelf Approximation cannot project future channel evolution. For thick ice shelves (m) secondary flow potentially stabilizes the channel, but is insufficient to significantly delay breakthrough for thinner ice (m). Using synthetic data, we assess the impact of secondary flow when inferring basal‐channel melt rates from satellite observations.

     
    more » « less
  2. Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner. 
    more » « less
  3. Flexure and extension of ice shelves in response to incident ocean surface gravity waves have been linked to iceberg calving, rift growth, and even disintegration of ice shelves. Most modeling studies utilize a plate bending model for the ice, focusing exclusively on flexural gravity waves. Ross Ice shelf seismic data shows not only flexural gravity waves, with dominantly vertical displacements, but also extensional Lamb waves, which propagate much faster with dominantly horizontal displacements. Our objective is to model the full-wave response of ice shelves, including ocean compressibility, ice elasticity, and gravity. Our model is a 2D vertical cross-section of the ice shelf and sub-shelf ocean cavity. We quantify the frequency-dependent excitation of flexural gravity and extensional Lamb waves and provide a quantitative theory for extensional Lamb wave generation by the horizontal force imparted by pressure changes on the vertical ice shelf edge exerted by gravity waves. Our model predicts a horizontal to vertical displacement ratio that increases with decreasing frequency, with ratio equal to unity at ~0.001 Hz. Furthermore, in the very long period band (<0.003 Hz), tilt from flexural gravity waves provides an order of magnitude larger contribution to seismometer horizontal components than horizontal displacements from extensional Lamb waves. 
    more » « less
  4. Abstract

    Surface melting and lakes are common to Antarctic ice shelves, and their existence and drainages have been invoked as a precursor for ice shelf collapse. Here, we present satellite observations over 2014–2020 of repeated, rapid drainages of a supraglacial lake at the grounding zone of Amery Ice Shelf, East Antarctica. Post‐drainage imagery in 2018 reveals lake bottom features characteristic of rapid, vertical lake drainage. Observed lake volumes indicate drainages are not associated with a threshold meltwater volume. Instead, drainages typically coincide with periods of high daily tidal amplitude, suggesting hydrofracture is assisted by tidally forced ice flexure inherent to the ice shelf grounding zone. Combined with observations of widespread grounding zone lake drainages on Amery, these findings indicate ice shelf meltwater accumulation may be inhibited by grounding zone drainage events, thus representing a potential stabilizing mechanism despite enhanced melting common to these regions.

     
    more » « less
  5. Abstract

    Rapid retreat of the Larsen A and B ice shelves has provided important clues about the ice shelf destabilization processes. The Larsen C Ice Shelf, the largest remaining ice shelf on the Antarctic Peninsula, may also be vulnerable to future collapse in a warming climate. Here, we utilize multisource satellite images collected over 1963–2020 to derive multidecadal time series of ice front, flow velocities, and critical rift features over Larsen C, with the aim of understanding the controls on its retreat. We complement these observations with modeling experiments using the Ice‐sheet and Sea‐level System Model to examine how front geometry conditions and mechanical weakening due to rifts affect ice shelf dynamics. Over the past six decades, Larsen C lost over 20% of its area, dominated by rift‐induced tabular iceberg calving. The Bawden Ice Rise and Gipps Ice Rise are critical areas for rift formation, through their impact on the longitudinal deviatoric stress field. Mechanical weakening around Gipps Ice Rise is found to be an important control on localized flow acceleration and the propagation of two rifts that caused a major calving event in 2017. Capturing the time‐varying effects of rifts on ice rigidity in ice shelf models is essential for making realistic predictions of ice shelf flow dynamics and instability. In the context of the Larsen A and Larsen B collapses, we infer a chronology of destabilization processes for embayment‐confined ice shelves, which provides a useful framework for understanding the historical and future destabilization of Antarctic ice shelves.

     
    more » « less