A new open-source viscoelastic solid earth deformation module implemented in Elmer (v8.4)

We present a new, open-source viscoelastic solid earth deformation model, Elmer/Earth. Using the multi-physics finite-element package Elmer, a model to compute viscoelastic material deformation has been implemented into the existing linear elasticity solver routine. Unlike approaches often implemented in engineering codes, our solver accounts for the restoring force of buoyancy within a system of layers with depth-varying density. It does this by directly integrating the solution of the system rather than by applying stress-jump conditions in the form of Winkler foundations on inter-layer boundaries, as is usually needed when solving the minimization problem given by the stress divergence in commercial codes. We benchmarked the new model with results from a commercial finite-element engineering package (ABAQUS, v2018) and another open-source code that uses viscoelastic normal mode theory, TABOO, using a flat-earth setup loaded by a cylindrical disc of 100 km in diameter and 100 m in height at the density of ice. Evaluating the differences in predicted surface deformation at the centre of the load and two distinctive distances (100 and 200 km), average deviations of 7 and 2.7 cm of Elmer/Earth results to ABAQUS and TABOO, respectively, were observed. In view of more than 100 cm maximum vertical deformation and the different numerical methods and parameters, these are very encouraging results. Elmer is set up as a highly scalable parallel code and distributed under the (L)GPL license, meaning that large-scale computations can be made without any licensing restrictions. Scaling figures presented in this paper show good parallel performance of the new model. Additionally, the high-fidelity ice-sheet code Elmer/Ice utilizes the same source base as Elmer and thereby the new model opens the way to undertaking high-resolution coupled ice-flow–solid-earth deformation simulations, which are required for robust projections of future sea-level rise and glacial isostatic adjustment.

[1]  P. Whitehouse,et al.  Glacial isostatic adjustment modelling: historical perspectives, recent advances, and future directions , 2018, Earth Surface Dynamics.

[2]  D. Pollard,et al.  A 3-D coupled ice sheet – sea level model applied to Antarctica through the last 40 ky , 2013 .

[3]  Brian Hutchison,et al.  Getting started. , 2019, Healthcare policy = Politiques de sante.

[4]  H. Steffen,et al.  Three-dimensional finite-element modeling of the glacial isostatic adjustment in Fennoscandia , 2006 .

[5]  C. Schoof.  Ice sheet grounding line dynamics: Steady states, stability, and hysteresis , 2007 .

[6]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[7]  Mika Malinen,et al.  Capabilities and performance of Elmer/Ice, a new-generation ice sheet model , 2013 .

[8]  Peter J. Clarke,et al.  Rapid bedrock uplift in the Antarctic Peninsula explained by viscoelastic response to recent ice unloading , 2014 .

[9]  S. Nowicki,et al.  Future Antarctic bed topography and its implications for ice sheet dynamics , 2014 .

[10]  On Postglacial Sea Level , 2007 .

[11]  D. Pollard,et al.  Sea-level feedback lowers projections of future Antarctic Ice-Sheet mass loss , 2015, Nature Communications.

[12]  George Karypis,et al.  Multilevel k-way Partitioning Scheme for Irregular Graphs , 1998, J. Parallel Distributed Comput..

[13]  A. Payne,et al.  The Glimmer community ice sheet model , 2009 .

[14]  E. Ivins,et al.  Slowdown in Antarctic mass loss from solid Earth and sea-level feedbacks , 2019, Science.

[15]  Mika Malinen,et al.  Elmer Models Manual , 2015 .

[16]  S. Eisenstat,et al.  Variational Iterative Methods for Nonsymmetric Systems of Linear Equations , 1983 .

[17]  Michael Klemm,et al.  OpenMP ^* SIMD Vectorization and Threading of the Elmer Finite Element Software , 2017, IWOMP.

[18]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[19]  R. S. W. van de Wal,et al.  A fully coupled 3-D ice-sheet–sea-level model: algorithm and applications , 2014 .

[20]  Jed Brown,et al.  Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations , 2007, Annals of Glaciology.

[21]  Patrick Wu,et al.  Dynamics of the Ice Age Earth: A Modern Perspective , 1998 .

[22]  P. Clark,et al.  Sea level as a stabilizing factor for marine-ice-sheet grounding lines , 2010 .

[23]  B. Smith,et al.  Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability , 2018, Science.

[24]  Patrick Wu,et al.  Using commercial finite element packages for the study of earth deformations, sea levels and the state of stress , 2004 .

[25]  W. Peltier.  The impulse response of a Maxwell Earth , 1974 .

[26]  P. Huybrechts,et al.  A comparison of different ways of dealing with isostasy: examples from modelling the Antarctic ice sheet during the last glacial cycle , 1996, Annals of Glaciology.