Modelling post-glacial rebound with lateral viscosity variations

SUMMARY Observations of isostatic adjustment of the Earth's surface due to glacial loading provide important constraints on mantle viscosity structure. We solve the forward problem of glacial isostatic adjustment in two complementary ways: a spectral method for strictly 1-D (spherically symmetric) earth models, and a finite element method that can accommodate 3-D viscosity structure. We discuss how each method may be augmented in three ways: to accommodate motion of the centre of mass, to implement a gravitationally self-consistent ocean load via the sea level equation and to include the influence of polar wander. With all these effects implemented, the two methods are benchmarked against each other. We also study the influence of lateral viscosity variations upon measurements of post-glacial rebound (PGR) in two ways: first by observing the effect of viscosity perturbations in an idealized model and second by developing a realistic 3-D viscosity model and comparing it with results of related 1-D (spherically symmetric) models. The 3-D viscosity structure is derived starting from seismic tomography models. We conclude from both approaches that PGR observations are sensitive to both the local viscosity structure and to the viscosity structure beneath the loaded region, even if it is removed from where the observations are made. In particular, PGR measurements made at Hudson Bay tend to reflect the local viscosity structure beneath Hudson Bay; PGR measurements made along the east coast of North America, being sensitive to both the local (east coast) viscosity structure as well as the loaded (continental) viscosity structure, are not reproducible with a 1-D viscosity model.

[1]  A. Paulson,et al.  Three-dimensional finite-element modelling of Earth's viscoelastic deformation: effects of lateral variations in lithospheric thickness , 2003 .

[2]  W. Peltier,et al.  A comparison of methods for the inversion of viscoelastic relaxation spectra , 1992 .

[3]  On Postglacial Sea Level , 2007 .

[4]  Michael H. Ritzwoller,et al.  Thermodynamic constraints on seismic inversions , 2004 .

[5]  S. Lee.  High-resolution estimates of lithospheric thickness from Missouri to Massachusetts, USA , 2002 .

[6]  Kurt Lambeck,et al.  The Earth's Variable Rotation: Geophysical Causes and Consequences , 1980 .

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

[8]  J. Wahr,et al.  Geodetic measurements in Greenland and their implications , 2001 .

[9]  S. Karato,et al.  Importance of anelasticity in the interpretation of seismic tomography , 1993 .

[10]  S. Karato,et al.  Development of finite strain in the convecting lower mantle and its implications for seismic anisotropy , 2003 .

[11]  W. Peltier,et al.  Pleistocene deglaciation and the Earth's rotation: a new analysis , 1984 .

[12]  D. Yamazaki,et al.  Some mineral physics constraints on the rheology and geothermal structure of Earth’s lower mantle , 2001 .

[13]  Philippe Huybrechts,et al.  Short term mass variability in Greenland, from GRACE , 2005 .

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

[15]  J. Woodhouse,et al.  Complex Shear Wave Velocity Structure Imaged Beneath Africa and Iceland. , 1999, Science.

[16]  Georg Kaufmann,et al.  Glacial isostatic adjustment in Fennoscandia for a laterally heterogeneous earth , 2000 .

[17]  J. Mitrovica.  Recent controversies in predicting post-glacial sea-level change , 2003 .

[18]  Z. Martinec,et al.  Can the 1D viscosity profiles inferred from postglacial rebound data be affected by lateral viscosity variations in the tectosphere? , 2001 .

[19]  A. Paulson,et al.  The rotational stability of an ice-age earth , 2005 .

[20]  Walter Munk,et al.  The rotation of the earth , 1960 .

[21]  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 .

[22]  W. Peltier,et al.  ICE-3G: A new global model of late Pleistocene deglaciation based upon geophysical predictions of po , 1991 .

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

[24]  M. Ritzwoller,et al.  Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle , 2002 .

[25]  W. Farrell.  Deformation of the Earth by surface loads , 1972 .

[26]  T. James,et al.  Polar motion of a viscoelastic Earth due to glacial cycle mass loading , 1997 .

[27]  Louis Moresi,et al.  Role of temperature‐dependent viscosity and surface plates in spherical shell models of mantle convection , 2000 .