A Generalized Interpolation Material Point Method for Shallow Ice Shelves. 1: Shallow Shelf Approximation and Ice Thickness Evolution

Huth, Alex; Duddu, Ravindra; Smith, Ben

doi:10.1029/2020ms002277

Cited by 4 publications

(9 citation statements)

References 71 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…This diffusion is inherent to purely Eulerian advection schemes, where the mesh is not moved with the computed velocity field, and can also arise when working in a Lagrangian frame (moving-mesh) due to frequent remeshing that may be required when modeling large-deformation processes like ice shelf flow. While our creep damage model may be adopted for any flow-modeling framework, we implement it here within our GIMPM-SSA code to avoid these diffusion errors (Huth et al, 2021). The GIMPM (Bardenhagen & Kober, 2004) is one of several material point methods, which all share the same basic procedure.…”

Section: Generalized Interpolation Materials Point Methods (Gimpm)mentioning

confidence: 99%

“…Consequently, horizontal velocities and the corresponding strain-rates can be assumed constant with depth. Assuming that vertical normal stress is equal to the overburden pressure, excluding vertical shear components from the incompressible Stokes equations, and vertically integrating yields the 2-D SSA, or SSA (Huth et al, 2021;MacAyeal, 1989;Weis, 2001…”

Section: Shallow Shelf Approximationmentioning

confidence: 99%

“…where a  (m a −1 ) is the sum of the basal and surface accumulation rates (Huth et al, 2021). In our GIMPM-SSA framework (Huth et al, 2021), ice velocity is determined using the Eulerian description, whereas ice thickness is explicitly updated at material points or particles using the Lagrangian description. Thus, it can be classified as a hybrid Eulerian-Lagrangian method.…”

Section: Shallow Shelf Approximationmentioning

confidence: 99%

“…Damage is incorporated into the SSA solution by replacing  ε in 2 with   ε, which is calculated from Equation 7using D as the damage variable. This substitution modifies the original SSA-GIMPM discretization (see Huth et al, 2021), yielding the following element submatrices of the tangent matrix, K, that are computed by summing over material points:…”

Section: Ssa Solution and Materials Point Updatesmentioning

confidence: 99%

“…We adapt several schemes within this framework to improve model performance and physical consistency, including extension of the damage variable to nonlocal form, adaptive time-stepping based on damage accumulation, brittle rupture criteria, and numerical treatment once maximum damage is reached. The damage model is implemented within our generalized interpolation material point method (GIMPM) code, which is a hybrid Eulerian-Lagrangian particle variation of the finite element method (Huth et al, 2021). Traditional Eulerian ice flow models are subject to artificial diffusion when advecting the damage field (e.g., Albrecht & Levermann, 2014;Borstad et al, 2016), whereas this error is avoided when using our GIMPM-SSA model, thereby allowing sharpness of cracks to be preserved regardless of flow.…”

mentioning

confidence: 99%

See 4 more Smart Citations

A Generalized Interpolation Material Point Method for Shallow Ice Shelves. 2: Anisotropic Nonlocal Damage Mechanics and Rift Propagation

Huth

Duddu

Smith

2021

J Adv Model Earth Syst

Self Cite

View full text Add to dashboard Cite

Our shallow-shelf creep damage model can represent the full evolution of ice shelf fracture from crevasse initiation to tabular calving  Strongly anisotropic creep damage produces sharp, arcuate rift patterns more consistent with observations than isotropic creep damage  The zero-stress damage model poorly captures rifting, but a modified form accounts for damage evolution due to mass balance and necking

show abstract

Section: Generalized Interpolation Materials Point Methods (Gimpm)mentioning

confidence: 99%

Section: Shallow Shelf Approximationmentioning

confidence: 99%

Section: Shallow Shelf Approximationmentioning

confidence: 99%

Section: Ssa Solution and Materials Point Updatesmentioning

confidence: 99%

mentioning

confidence: 99%

See 3 more Smart Citations

A Generalized Interpolation Material Point Method for Shallow Ice Shelves. 2: Anisotropic Nonlocal Damage Mechanics and Rift Propagation

Huth

Duddu

Smith

2021

J Adv Model Earth Syst

Self Cite

View full text Add to dashboard Cite

show abstract

Simulating the processes controlling ice-shelf rift paths using damage mechanics

Huth,

Duddu,

Smith

et al. 2023

J. Glaciol.

Self Cite

View full text Add to dashboard Cite

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.

show abstract