Abstract:Abstract. Time dependent simulations of ice sheets require two equations to be solved: the mass transport equation, derived from the conservation of mass, and the stress balance equation, derived from the conservation of momentum. The mass transport equation controls the advection of ice from the interior of the ice sheet towards its periphery, thereby changing its geometry. Because it is based on a hyperbolic partial differential equation, a stabilization scheme needs to be employed when solved using the finite element method. Several stabilization schemes exist in the finite element method framework, but their respective accuracy and robustness have not yet been systematically assessed for glaciological applications. Here, we compare classical schemes used in the context of the finite element method: (i) Artificial Diffusion, (ii) Streamline Upwinding, (iii) Streamline Upwind Petrov-Galerkin, (iv) Discontinuous Galerkin, and (v) Flux Corrected Transport. We also look at the stress balance equation, which is responsible for computing the ice velocity that `advects' the ice dowstream. To improve the velocity computation accuracy, the ice sheet modeling community employs several sub-element parameterizations of physical processes at the grounding line, the point where the grounded ice starts to float onto the ocean. Here, we introduce a new sub-element parameterization for the driving stress, the force that drives the ice sheet flow. We analyze the response of each stabilization scheme by running transient simulations forced by ice shelf basal melt. The simulations are based on an idealized ice sheet geometry for which there is no influence of bedrock topography. We also perform transient simulations of the Amundsen Sea Sector, West Antarctica, where real bedrock and surface elevations are employed. In both idealized and real ice sheet experiments, stabilization schemes based on artificial diffusion lead systematically to a bias towards more mass loss in comparison to the other schemes, and therefore, should be avoided or employed with a sufficiently high mesh resolution in the vicinity of the grounding line. We also run diagnostic simulations to assess the accuracy of the driving stress parameterization, which in combination with an adequate parameterization for basal stress, provides improved numerical convergence in ice speed computations and more accurate results.

APPLICABILITY PISM FOR VELOCITY ANALYSIS OF THE AMERY ICE SHELF, EAST ANTARCTICA

The Spatiotemporal Surface Velocity Variations and Analysis of the Amery Ice Shelf from 2000 to 2022, East Antarctica

Surface Velocity Estimations of Ice Shelves in the Northern Antarctic Peninsula Derived from MODIS Data

RESEARCH PROGRESS ON ICE VELOCITY OF ANTARCTIC ICE SHEET

Simulation-Based Inference of Surface Accumulation and Basal Melt Rates of an Antarctic Ice Shelf from Isochronal Layers

Antarctic Ice Sheet GLIMMER Model Test and Its Simplified Model on 2-Dimensional Ice Flow

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

Mapping Complete Three-Dimensional Ice Velocities by Integrating Multi-Baseline and Multi-Aperture InSAR Measurements: A Case Study of the Grove Mountains Area, East Antarctic

A framework for time-dependent ice sheet uncertainty quantification, applied to three West Antarctic ice streams

Spatiotemporal variations in the surface velocities of Antarctic Peninsula glaciers

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

Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18

Impact of boundary conditions on the modeled thermal regime of the Antarctic ice sheet

Sensitivity of simulated water mass transformation on the Antarctic shelf to tides, topography and model resolution

The variation in basal channels and basal melt rates of Pine Island Ice Shelf

Quantifying the Effect of Ocean Bed Properties on Ice Sheet Geometry over 40,000 Years with a Full-Stokes Model

OBSERVATIONS OF CONTINUOUS INSTABILITY FOR SCAR INLET ICE SHELF, ANTARCTIC PENINSULA

Revisiting Ice Flux and Mass Balance of the Lambert Glacier–Amery Ice Shelf System Using Multi-Remote-Sensing Datasets, East Antarctica

Response of Onshore Oceanic Heat Supply to Yearly Changes in the Amundsen Sea Icescape (Antarctica)

Geometric amplification and suppression of ice-shelf basal melt in West Antarctica

Continent‐wide, interferometric SAR phase, mapping of Antarctic ice velocity