Ice shelf fracture is responsible for roughly half of Antarctic ice mass loss in the form of calving and can weaken buttressing of upstream ice flow. Large uncertainties associated with the ice sheet response to climate variations are due to a poor understanding of these fracture processes and how to model them. Here, we address these problems by implementing an anisotropic, nonlocal integral formulation of creep damage within a large‐scale shallow‐shelf ice flow model. This model can be used to study the full evolution of fracture from initiation of crevassing to rifting that eventually causes tabular calving. While previous ice shelf fracture models have largely relied on simple expressions to estimate crevasse depths, our model parameterizes fracture as a progressive damage evolution process in three‐dimensions (3‐D). We also implement an efficient numerical framework based on the material point method, which avoids advection errors. Using an idealized marine ice sheet, we test the creep damage model and a crevasse‐depth based damage model, including a modified version of the latter that accounts for damage evolution due to necking and mass balance. We demonstrate that the creep damage model is best suited for capturing weakening and rifting over shorter (monthly/yearly) timescales, and that anisotropic damage reproduces typically observed fracture patterns better than isotropic damage. Because necking and mass balance can significantly influence damage on longer (decadal) timescales, we discuss the potential for a combined approach between models to best represent mechanical weakening and tabular calving within long‐term simulations.
- Award ID(s):
- 1847173
- PAR ID:
- 10495675
- Publisher / Repository:
- Cambridge Core
- Date Published:
- Journal Name:
- Journal of Glaciology
- ISSN:
- 0022-1430
- Page Range / eLocation ID:
- 1 to 14
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
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.
-
A quasi-one-dimensional ice mélange flow model based on continuum descriptions of granular materialsmore » « less
Abstract. 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
Abstract Tabular calving events occur from Antarctica's large ice shelves only every few decades, and are preceded by rift propagation. We used high‐resolution imagery and ICESat‐2 data to determine the propagation rates for the three active rifts on Amery Ice Shelf (AIS; T1, T2, and E3) and observe the calving of D‐28 on September 25, 2019 along T1. AIS front advance accelerated downstream of T1 in the years before calving, possibly increasing stress at the rift tip. T1 experienced significant acceleration for 12 days before calving, coinciding with a jump in propagation of E3. ICESat‐2's high resolution and repeat acquisitions every 91 days allowed for analysis of the ice front before and after calving, and rift detection where it was not visible in imagery as a ∼1 m surface depression, suggesting that it propagates as a basal fracture. Our results show that ICESat‐2 can provide process‐scale information about iceberg calving.
-
more » « less
Abstract Iceberg calving is a major contributor to Greenland’s ice mass loss. Ice mélange, tightly packed sea ice and icebergs, has been hypothesized to buttress the calving fronts. However, quantifying the mélange buttressing force from field observations remains a challenge. Here we show that such quantification can be achieved with a single field measurement: thickness of mélange at the glacier terminus. We develop the first three-dimensional discrete element model of m´elange along with a simple analytical model to quantify the mélange buttressing using mélange thickness data from ArcticDEM over 32 Greenland glacier termini. We observed a strong seasonality in mélange thickness: thin mélange (averaged thickness 34+17−15m) in summertime when terminus retreats, and thick mélange (averaged thickness 119+31−37m) in wintertime when terminus advances. The observed seasonal changes of mélange thickness strongly coincide with observed Greenland calving dynamics and the modeled buttressing effects.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jog.2023.71
