skip to main content


Title: A Generalized Interpolation Material Point Method for Shallow Ice Shelves. 2: Anisotropic Nonlocal Damage Mechanics and Rift Propagation
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.

 
more » « less
Award ID(s):
1847173
PAR ID:
10374688
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
13
Issue:
8
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    Frontal ablation processes at marine‐terminating glaciers are challenging to observe and difficult to represent in numerical ice flow models, yet play critical roles in modulating ice sheet mass balance. Current ice sheet models typically rely on simple iceberg calving models to prescribe either terminus positions or iceberg calving rates, but the relative accuracies and uncertainties of these calving models remain largely unconstrained at the ice sheet scale. Here, we evaluate six published iceberg calving models against spatially and temporally diverse observations from 50 marine‐terminating outlet glaciers in Greenland. We seek the single model that best reproduces observed conditions across all glaciers, at all observation times, and with low sensitivity to calibration uncertainty. Five of six calving models can produce unbiased estimates of calving position or calving rate at the ice sheet scale. However, time series analysis reveals that, when using a single, optimized model parameter, rate‐predicting calving models frequently yield calving rate errors in excess of 10 m d−1. In comparison, terminus position‐predicting calving models more accurately track observed changes in terminus position (remaining within ~1 km of the observed grounded terminus position). Overall, our results indicate that the crevasse depth calving model provides the best balance of high accuracy and low sensitivity to imperfect parameter calibration. While the crevasse depth model appears unlikely to capture the true controls on crevasse penetration, numerically, it reproduces observed terminus dynamics with high fidelity and should be considered a leading candidate for use in models of the Greenland Ice Sheet.

     
    more » « less
  3. Irina Tezaur, Josefin Ahlkrona (Ed.)
    We present a finite-element-based cohesive zone model for simulating the nonlinear fracture process driving the propagation of water-filled surface crevasses in floating ice tongues. The fracture process is captured using an interface element whose constitutive behavior is described by a bilinear cohesive law, and the bulk rheology of ice is described by a nonlinear elasto-viscoplastic model. The additional loading due to meltwater pressure within the crevasse is incorporated by combining the ideas of poromechanics and damage mechanics.We performed several numerical studies to explore the parametric sensitivity of surface crevasse depth to ice rheology, cohesive strength, density, and temperature for different levels of meltwater depth.We find that viscous (creep) strain accumulation promotes crevasse propagation and that surface crevasses propagate deeper in ice shelves/tongues if we consider depth-varying ice density and temperature profiles. Therefore, ice flow models must account for depth-varying density and temperature-dependent viscosity to appropriately describe calving outcomes. 
    more » « less
  4. Abstract. Many floating ice shelves in Antarctica buttress the ice streams feeding them, thereby reducing the discharge of icebergs into the ocean. The rate at which ice shelves calve icebergs and how fast they flow determines whether they advance, retreat, or remain stable, exerting a first-order control on ice discharge. To parameterize calving within ice sheet models, several empirical and physical calving “laws” have been proposed in the past few decades. Such laws emphasize dissimilar features, including along- and across-flow strain rates (the eigencalving law), a fracture yield criterion (the von Mises law), longitudinal stretching (the crevasse depth law), and a simple ice thickness threshold (the minimum thickness law), among others. Despite the multitude of established calving laws, these laws remain largely unvalidated for the Antarctic Ice Sheet, rendering it difficult to assess the broad applicability of any given law in Antarctica. We address this shortcoming through a set of numerical experiments that evaluate existing calving laws for ten ice shelves around the Antarctic Ice Sheet. We utilize the Ice-sheet and Sea-level System Model (ISSM) and implement four calving laws under constant external forcing, calibrating the free parameter of each of these calving laws by assuming that the current position of the ice front is in steady state and finding the set of parameters that best achieves this position over a simulation of 200 years. We find that, in general, the eigencalving and von Mises laws best reproduce observed calving front positions under the steady state position assumption. These results will streamline future modeling efforts of Antarctic ice shelves by better informing the relevant physics of Antarctic-style calving on a shelf-by-shelf basis.

     
    more » « less
  5. Abstract

    The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

     
    more » « less