Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Rifts are full-thickness fractures that propagate laterally across an ice shelf. They cause ice-shelf weakening and calving of tabular icebergs, and control the initial size of calved icebergs. Here, we present a joint inverse and forward computational modeling framework to capture rifting by combining the vertically integrated momentum balance and anisotropic continuum damage mechanics formulations. We incorporate rift–flank boundary processes to investigate how the rift path is influenced by the pressure on rift–flank walls from seawater, contact between flanks, and ice mélange that may also transmit stress between flanks. To illustrate the viability of the framework, we simulate the final 2 years of rift propagation associated with the calving of tabular iceberg A68 in 2017. We find that the rift path can change with varying ice mélange conditions and the extent of contact between rift flanks. Combinations of parameters associated with slower rift widening rates yield simulated rift paths that best match observations. Our modeling framework lays the foundation for robust simulation of rifting and tabular calving processes, which can enable future studies on ice-sheet–climate interactions, and the effects of ice-shelf buttressing on land ice flow.more » « less
-
Abstract We develop a generalized interpolation material point method (GIMPM) for the shallow shelf approximation (SSA) of ice flow. The GIMPM, which can be viewed as a particle version of the finite element method, is used here to solve the shallow shelf approximations of the momentum balance and ice thickness evolution equations. We introduce novel numerical schemes for particle splitting and integration at domain boundaries to accurately simulate the spreading of an ice shelf. The advantages of the proposed GIMPM‐SSA framework include efficient advection of history or internal state variables without diffusion errors, automated tracking of the ice front and grounding line at sub‐element scales, and a weak formulation based on well‐established conventions of the finite element method with minimal additional computational cost. We demonstrate the numerical accuracy and stability of the GIMPM using 1‐D and 2‐D benchmark examples. We also compare the accuracy of the GIMPM with the standard material point method (sMPM) and a reweighted form of the sMPM. We find that the grid‐crossing error is very severe for SSA simulations with the sMPM, whereas the GIMPM successfully mitigates this error. While the grid‐crossing error can be reasonably reduced in the sMPM by implementing a simple material point reweighting scheme, this approach it not as accurate as the GIMPM. Thus, we illustrate that the GIMPM‐SSA framework is viable for the simulation of ice sheet‐shelf evolution and enables boundary tracking and error‐free advection of history or state variables, such as ice thickness or damage.
-
Abstract 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.