A benchmark study of numerical implementations of the sea level equation in GIA modelling

The ocean load in glacial isostatic adjustment (GIA) modelling is represented by the so-called sea level equation (SLE). The SLE describes the mass redistribution of water between ice sheets and oceans on a deforming Earth. Despite various teams independently investigating GIA, there has been no systematic intercomparison among the numerical solvers of the SLE through which the methods may be validated. The goal of this paper is to present a series of synthetic examples designed for testing and comparing the numerical implementations of the SLE in GIA modelling. The 10 numerical codes tested combine various temporal and spatial parametrizations. The time-domain or Laplace-domain discretizations are used to solve the SLE through time, while spherical harmonics, finite differences or finite elements parametrize the GIA-related field variables spatially. The surface ice-water load and solid Earth's topography are represented spatially either on an equiangular grid, a Gauss-Legendre or an equiarea grid with icosahedron-shaped spherical pixels. Comparisons aremade in a series of five benchmark examples with an increasing degree of complexity. Due to the complexity of the SLE, there is no analytical solution to it. The accuracy of the numerical implementations is therefore assessed by the differences of the individual solutions with respect to a reference solution. While the benchmark study does not result in GIA predictions for a realistic loading scenario, we establish a set of agreed-upon results that can be extended in the future by including more complex case studies, such as solutions with realistic loading scenarios, the rotational feedback in the linear-momentum equation, and by considering a 3-D viscosity structure of the Earth's mantle. The test computations performed so far show very good agreement between the individual results and their ability to capture the main features of sea-surface variation and the surface vertical displacement. The differences found can often be attributed to the different approximations inherent in the various algorithms. This shows the accuracy that can be expected from different implementations of the SLE, which helps to assess differences noted in the literature between predictions for realistic loading cases.

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

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

[3]  On Postglacial Sea Level , 2007 .

[4]  Z. Martinec,et al.  ISMIP-HEINO experiment revisited: effect of higher-order approximation and sensitivity study , 2011, Journal of Glaciology.

[5]  J. Bamber,et al.  Sea‐level fingerprint of continental water and ice mass change from GRACE , 2010 .

[6]  Z. Martinec,et al.  Time-domain approach to linearized rotational response of a three-dimensional viscoelastic earth model induced by glacial-isostatic adjustment: I. Inertia-tensor perturbations , 2005 .

[7]  Z. Martinec,et al.  An Estimate of Global Mean Sea-level Rise Inferred from Tide-gauge Measurements Using Glacial-isostatic Models Consistent with the Relative Sea-level Record , 2007 .

[8]  Z. Martinec,et al.  Glacial isostasy and plate motion , 2008 .

[9]  Max Tegmark,et al.  An Icosahedron-based Method for Pixelizing the Celestial Sphere , 1996, The Astrophysical Journal.

[10]  Z. Martinec.  Spectral–finite element approach to three‐dimensional viscoelastic relaxation in a spherical earth , 2000 .

[11]  W. Wal,et al.  The effect of sediment loading in Fennoscandia and the Barents Sea during the last glacial cycle on glacial isostatic adjustment observations , 2017 .

[12]  T. James,et al.  A glacial isostatic adjustment model for the central and northern Laurentide Ice Sheet based on relative sea level and GPS measurements , 2016 .

[13]  W. Peltier,et al.  On postglacial geoid subsidence over the equatorial oceans , 1991 .

[14]  Carling C. Hay,et al.  On the robustness of predictions of sea level fingerprints , 2011 .

[15]  N. Dael.  Sensitivity of glacial isostaticadjustment induced rotational changesto the ice-load history , 2006 .

[16]  Sybren S. Drijfhout,et al.  Onderzoek naar bovengrensscenario's voor klimaatverandering voor overstromingsbescherming van Nederland : internationale wetenschappelijke beoordeling , 2008 .

[17]  L. Vermeersen,et al.  Sensitivity of glacial isostatic adjustment models with shallow low-viscosity earth layers to the ice-load history in relation to the performance of GOCE and GRACE , 2005 .

[18]  Giorgio Spada,et al.  SELEN: A Fortran 90 program for solving the "sea-level equation" , 2007, Comput. Geosci..

[19]  J. Wahr,et al.  Post‐Glacial Rebound Analysis for a Rotating Earth , 2013 .

[20]  K. Lambeck,et al.  Late Pleistocene and Holocene sea-level change along the Australian coast , 1990 .

[21]  D. Komatitsch,et al.  Glacial isostatic adjustment on 3‐D Earth models: a finite‐volume formulation , 2005 .

[22]  V. Barletta,et al.  Effect of different implementations of the same ice history in GIA modeling , 2013 .

[23]  E. Ivins,et al.  Ocean loading effects on the prediction of Antarctic glacial isostatic uplift and gravity rates , 2010 .

[24]  Z. Martinec,et al.  Contribution of glacial-isostatic adjustment to the geocenter motion , 2011 .

[25]  J. Mitrovica,et al.  Postglacial sea-level change on a rotating Earth , 1998 .

[26]  Lapo Boschi,et al.  MODELING EARTH'S POST-GLACIAL REBOUND , 2004 .

[27]  Patrick Wu,et al.  Lateral asthenospheric viscosity variations and postglacial rebound: A case study for the Barents Sea , 1998 .

[28]  K. Lambeck,et al.  Glacial rebound and relative sea-level variations: a new appraisal , 1987 .

[29]  I. Shapiro,et al.  A spectral formalism for computing three‐dimensional deformations due to surface loads: 1. Theory , 1994 .

[30]  K. Lambeck,et al.  Water-load definition in the glacio-hydro-isostatic sea-level equation , 2003 .

[31]  J. Mitrovica,et al.  Near-field hydro-isostasy: the implementation of a revised sea-level equation , 1999 .

[32]  Michael G. Sideris,et al.  Sea levels and uplift rate from composite rheology in glacial isostatic adjustment modeling , 2010 .

[33]  I. Sasgen,et al.  Sensitivity of grounding-line dynamics to viscoelastic deformation of the solid-earth in an idealized scenario , 2016 .

[34]  W. Peltier,et al.  Viscous gravitational relaxation , 1982 .

[35]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[36]  J. Mitrovica,et al.  On post-glacial sea level: I. General theory , 2003 .

[37]  E. Ivins,et al.  Predictions of Antarctic crustal motions driven by present-day ice sheet evolution and by isostatic memory of the Last Glacial Maximum , 1998 .

[38]  Riccardo E. M. Riva,et al.  A benchmark study for glacial isostatic adjustment codes , 2011 .

[39]  W. Wal,et al.  Postglacial sealevels on a spherical, self-gravitating viscoelastic earth: effects of lateral viscosity variations in the upper mantle on the inference of viscosity contrasts in the lower mantle , 2003 .

[40]  W. R. Peltier,et al.  Postglacial variations in the level of the sea: Implications for climate dynamics and solid‐Earth geophysics , 1998 .

[41]  J. Mitrovica,et al.  On post-glacial sea level – II. Numerical formulation and comparative results on spherically symmetric models , 2005 .

[42]  Z. Martinec,et al.  The rotational feedback on linear-momentum balance in glacial isostatic adjustment , 2014 .

[43]  P. Gasperini,et al.  Lateral heterogeneities in mantle viscosity and post‐glacial rebound , 1989 .

[44]  I. Sasgen,et al.  The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model , 2014, Surveys in Geophysics.

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

[46]  H. Schotman.  Shallow-earth rheology from glacial isostasy and satellite gravity : A sensitivity analysis for GOCE , 2008 .

[47]  J. Wahr,et al.  Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada , 2012 .