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.

Advancements in sea ice dynamics modeling based on a mixed least‐squares finite element study with nonconforming stress approximation

Least-Squares finite element method for the simulation of sea-ice motion

Least-Squares Finite Element Method for Sea-Ice Dynamics

Ocean circulation and sea ice distribution in a finite element global sea ice–ocean model

Dynamically Coupling Full Stokes and Shallow Shelf Approximation for Marine Ice Sheet Flow Using Elmer/Ice (V8.3)

A Computational Fluid Dynamics Model for the Small-Scale Dynamics of Wave, Ice Floe and Interstitial Grease Ice Interaction

Bonded discrete element simulations of sea ice with non-local failure: Applications to Nares Strait

Numerical Simulation of Sea Ice and Structure Interaction Using Common Node DEM-SPH Model

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

Numerical Simulation Study on the Ice Load and Ship Motion Response in Formation Navigation in Brash Ice Channels

A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling

An unstructured CD-grid variational formulation for sea ice dynamics

Understanding the influence of ocean waves on Arctic sea ice simulation: a modeling study with an atmosphere–ocean–wave–sea ice coupled model

Thermodynamic growth of sea ice: assessing the role of salinity using a quasi-static modelling framework

The Finite Element Sea ice-Ocean Model (FESOM): formulation of an unstructured-mesh ocean general circulation model

Increasing stable time-step sizes of the free-surface problem arising in ice-sheet simulations

Particle-Continuum Multiscale Modeling of Sea Ice Floes

SPH Modelling of Sea-ice Pack Dynamics

Improving numerical accuracy for the viscous-plastic formulation of sea ice

Wave-ice dynamical interaction: a numerical model and its application