Online citations, reference lists, and bibliographies.
← Back to Search

Assessment Of Numerical Schemes For Transient, Finite-element Ice Flow Models Using ISSM V4.18

Thiago Dias dos Santos, M. Morlighem, H. Seroussi
Published 2021 · Geology

Save to my Library
Download PDF
Analyze on Scholarcy Visualize in Litmaps
Share
Reduce the time it takes to create your bibliography by a factor of 10 by using the world’s favourite reference manager
Time to take this seriously.
Get Citationsy
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 an advection 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 downstream. 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 Embayment, 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.
This paper references
Triangular mesh methods for the neutron transport equation
W. H. Reed (1973)
10.1016/0021-9991(73)90147-2
Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works
J. Boris (1973)
10.1017/S0022143000023327
Stability of the junction of an ice sheet and an ice shelf
J. Weertman (1974)
10.1002/NME.1620100617
Finite element methods for second order differential equations with significant first derivatives
I. Christie (1976)
10.1016/0045-7930(76)90010-4
An evaluation of upwind and central difference approximations by a study of recirculating flow
G. D. Davis (1976)
10.1002/NME.1620110113
An ‘upwind’ finite element scheme for two‐dimensional convective transport equation
J. Heinrich (1977)
10.1016/0045-7825(78)90012-9
An analysis of the petrov—galerkin finite element method
D. Griffiths (1978)
10.1016/0021-9991(79)90051-2
Fully multidimensional flux-corrected transport algorithms for fluids
S. Zalesak (1979)
10.1002/NME.1620151111
A note on upwinding and anisotropic balancing dissipation in finite element approximations to convective diffusion problems
D. Kelly (1980)
10.2307/1551118
The Physics of Glaciers
R. Armstrong (1981)
10.1002/NME.1620200108
A Taylor–Galerkin method for convective transport problems
J. Donea (1983)
10.1016/0045-7825(84)90153-1
Time-accurate solution of advection-diffusion problems by finite elements
J. Donea (1984)
10.1016/0029-5493(84)90163-8
Recent advances in computational methods for steady and transient transport problems
J. Donea (1984)
10.1002/FLD.1650071007
Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations
R. Löhner (1987)
10.1029/JB094IB04P04071
Large‐scale ice flow over a viscous basal sediment: Theory and application to ice stream B, Antarctica
D. Macayeal (1989)
10.1016/0045-7825(82)90071-8
Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations
A. Brooks (1990)
10.2307/2938728
Numerical methods for conservation laws
R. LeVeque (1990)
10.1051/M2AN/1991250303371
The Runge-Kutta local projection $P^1$-discontinuous-Galerkin finite element method for scalar conservation laws
Bernardo Cockburn (1991)
10.1007/978-3-0348-8629-1_3
Scalar Conservation Laws
R. LeVeque (1992)
10.1016/0045-7825(93)90143-L
The node-centred finite volume approach: bridge between finite differences and finite elements
V. Selmin (1993)
10.2307/2153181
Finite Element Analysis.
B. Szabó (1993)
10.1080/10618569508904524
ANALYSIS AND DESIGN OF NUMERICAL SCHEMES FOR GAS DYNAMICS, 1: ARTIFICIAL DIFFUSION, UPWIND BIASING, LIMITERS AND THEIR EFFECT ON ACCURACY AND MULTIGRID CONVERGENCE
A. Jameson (1995)
10.1017/S002214300001621X
Velocity and stress fields in grounded glaciers: a simple algorithm for including deviatoric stress gradients
H. Blatter (1995)
10.1007/BFb0014497
Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator
J. Shewchuk (1996)
10.1016/S0045-7825(96)01108-5
A stable Petrov-Galerkin method for convection-dominated problems
R. C. Almeida (1997)
10.1007/BFB0096353
An introduction to the Discontinuous Galerkin method for convection-dominated problems
Bernardo Cockburn (1998)
10.1016/S0045-7825(97)00206-5
Comparison of some finite element methods for solving the diffusion-convection-reaction equation
R. Codina (1998)
10.1016/S0898-1221(99)00117-0
A discontinuous hp finite element method for diffusion problems: 1-D analysis☆
I. Babuska (1999)
10.1007/978-3-662-03882-6_2
Discontinuous Galerkin Methods for Convection-Dominated Problems
Bernardo Cockburn (1999)
10.1016/S0045-7825(98)00279-5
Stabilised finite element methods for steady incompressible flow
Henrik Blank (1999)
10.1002/1098-2426(200007)16:4<365::AID-NUM2>3.0.CO;2-Y
Discontinuous Galerkin approximations for elliptic problems
F. Brezzi (2000)
10.1002/ZAMM.200310088
Discontinuous Galerkin methods
Bernardo Cockburn (2000)
10.1016/S0045-7825(00)00177-8
On stabilized finite element methods for linear systems of convection-diffusion-reaction equations
R. Codina (2000)
10.1006/JCPH.2001.6955
Flux correction tools for finite elements
D. Kuzmin (2002)
10.1016/S0045-7825(02)00222-0
Stabilized finite element methods with shock capturing for advection–diffusion problems
T. Knopp (2002)
10.1137/S0036142901384162
Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
D. Arnold (2002)
10.1002/FLD.493
Multidimensional FEM‐FCT schemes for arbitrary time stepping
D. Kuzmin (2003)
10.1016/J.CMA.2004.01.026
Stability of the SUPG finite element method for transient advection-diffusion problems
P. Bochev (2004)
10.1007/978-3-642-18560-1_20
Implicit Flux-Corrected Transport Algorithm for Finite Element Simulation of the Compressible Euler Equations
M. Möller (2004)
10.1016/J.CMA.2003.12.050
Calculation of the advective limit of the SUPG stabilization parameter for linear and higher-order elements
J. Akin (2004)
10.1142/S0218202504003866
DISCONTINUOUS GALERKIN METHODS FOR FIRST-ORDER HYPERBOLIC PROBLEMS
F. Brezzi (2004)
10.1029/2004GL022024
Assessment of the importance of ice‐shelf buttressing to ice‐sheet flow
T. Dupont (2005)
10.1016/J.CMA.2004.06.036
Stabilized discontinuous Galerkin method for hyperbolic equations
Jorge L. Díaz Calle (2005)
10.1016/J.CMA.2005.05.028
Revisiting stabilized finite element methods for the advective–diffusive equation
L. Franca (2006)
A combined continuous-discontinuous finite element method for convection-diffusion problems
P. Devloo (2007)
10.1017/S0022112006003570
Marine ice-sheet dynamics. Part 1. The case of rapid sliding
C. Schoof (2007)
10.1029/2006JF000664
Ice sheet grounding line dynamics: Steady states, stability, and hysteresis
C. Schoof (2007)
10.1016/J.CMA.2008.08.016
On Finite Element Methods for 3D Time-Dependent Convection-Diffusion-Reaction Equations with Small Diffusion
V. John (2008)
10.1029/2008JF001227
Numerical and theoretical treatment of grounding line movement and ice shelf buttressing in marine ice sheets
D. Goldberg (2009)
10.1007/978-3-642-03415-2
Dynamics of Ice Sheets and Glaciers
R. Greve (2009)
10.3189/172756409789624283
Full Stokes modeling of marine ice sheets: influence of the grid size
G. Durand (2009)
10.1016/j.jcp.2008.12.011
Explicit and implicit FEM-FCT algorithms with flux linearization
D. Kuzmin (2009)
10.1016/J.CMA.2009.11.023
Consistent SUPG-method for transient transport problems: Stability and convergence
E. Burman (2010)
10.1029/2010GL043853
Spatial patterns of basal drag inferred using control methods from a full-Stokes and simpler models for Pine Island Glacier,
West Antarctica (2010)
10.1029/2011GL048892
Impact of bedrock description on modeling ice sheet dynamics
G. Durand (2011)
10.5194/TC-6-573-2012
Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP
F. Pattyn (2012)
10.1029/2011JF002140
Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM)
E. Larour (2012)
10.5194/TC-6-1497-2012
The stability of grounding lines on retrograde slopes
G. Gudmundsson (2012)
10.1016/j.jcp.2012.08.037
Adaptive mesh, finite volume modeling of marine ice sheets
S. Cornford (2013)
10.5194/GMD-6-1299-2013
Capabilities and performance of Elmer/Ice, a new-generation ice sheet model
O. Gagliardini (2013)
10.1002/JGRF.20125
Inversion of basal friction in Antarctica using exact and incomplete adjoints of a higher‐order model
M. Morlighem (2013)
10.3189/2013JoG12J129
Grounding-line migration in plan-view marine ice-sheet models: results of the ice2sea MISMIP3d intercomparison
F. Pattyn (2013)
10.5194/TC-8-1699-2014
Sensitivity of the dynamics of Pine Island Glacier, West Antarctica, to climate forcing for the next 50 years
H. Seroussi (2014)
10.3189/2014JOG14J051
Improved representation of East Antarctic surface mass balance in a regional atmospheric climate model
J. M. V. Wessem (2014)
10.1126/science.1249055
Marine Ice Sheet Collapse Potentially Under Way for the Thwaites Glacier Basin, West Antarctica
I. Joughin (2014)
10.3189/2014JoG13J093
Resolution-dependent performance of grounding line motion in a shallow model compared with a full-Stokes model according to the MISMIP3d intercomparison
J. Feldmann (2014)
10.1016/J.CMA.2015.06.008
Numerical solution of steady-state groundwater flow and solute transport problems: Discontinuous Galerkin based methods compared to the Streamline Diffusion approach
A. Ngo (2014)
10.5194/TC-8-2075-2014
Hydrostatic grounding line parameterization in ice sheet models
H. Seroussi (2014)
10.3189/2014JoG13J117
Thwaites Glacier grounding-line retreat: influence of width and buttressing parameterizations
D. Docquier (2014)
10.1038/NCLIMATE2094
Retreat of Pine Island Glacier controlled by marine ice-sheet instability
L. Favier (2014)
10.1038/nature16147
Potential sea-level rise from Antarctic ice-sheet instability constrained by observations
C. Ritz (2015)
10.5194/TC-9-1579-2015
Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate
S. Cornford (2015)
10.1002/2015WR017164
Numerical stability analysis of two‐dimensional solute transport along a discrete fracture in a porous rock matrix
N. Watanabe (2015)
10.1017/aog.2016.13
Adaptive mesh refinement versus subgrid friction interpolation in simulations of Antarctic ice dynamics
S. Cornford (2016)
10.1038/nature17145
Contribution of Antarctica to past and future sea-level rise
R. DeConto (2016)
10.1016/j.jcp.2017.03.049
Detailed finite element method modeling of evaporating multi-component droplets
C. Diddens (2017)
10.1002/2017GL072910
Continued retreat of Thwaites Glacier, West Antarctica, controlled by bed topography and ocean circulation
H. Seroussi (2017)
10.1007/s00791-018-0290-5
Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?
V. John (2018)
10.5194/TC-12-3085-2018
Representation of basal melting at the grounding line in ice flow models
H. Seroussi (2018)
10.1038/s41598-019-53190-6
Pathways of ocean heat towards Pine Island and Thwaites grounding lines
Y. Nakayama (2019)
10.1038/s41561-019-0510-8
Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet
M. Morlighem (2019)
10.1029/2019GL083826
Continent‐Wide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity
J. Mouginot (2019)
10.1002/cnm.3148
Finite element modeling of near-wall mass transport in cardiovascular flows.
Kirk B Hansen (2019)
10.1029/2018GL081229
Millennial-Scale Vulnerability of the Antarctic Ice Sheet to Regional Ice Shelf Collapse
Daniel F. Martin (2019)
10.1126/sciadv.aav9396
Contribution of the Greenland Ice Sheet to sea level over the next millennium
A. Aschwanden (2019)
10.1029/2019gl085027
Instantaneous Antarctic ice sheet mass loss driven by thinning ice shelves
G. Gudmundsson (2019)
10.5194/TC-14-3071-2020
The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6
H. Goelzer (2020)
10.5194/tc-2020-235
The transferability of adjoint inversion products between different ice flow models
J. M. Barnes (2020)
10.5194/TC-14-3033-2020
ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century
H. Seroussi (2020)
10.5194/TC-15-1975-2021
The transferability of adjoint inversion products between different ice flow models
J. M. Barnes (2021)



Semantic Scholar Logo Some data provided by SemanticScholar