Limitations on the inversion for mantle viscosity from postglacial rebound

SUMMARY Observations of postglacial rebound (PGR) can provide important constraints on mantle viscosity structure. In this study, we investigate how well PGR observations are able to constrain the spherically symmetric (1-D) viscosity structure of the Earth. We generate synthetic PGR data by calculating the response of an earth model with realistic 3-D viscosity. The viscosity model is constructed starting from seismic tomography models. We generate synthetic PGR data from this model including relative sea levels, exponential relaxation times, ˙ J2, polar wander, and GRACE time-variable gravity measurements, where most of the data is concentrated in the Laurentide region. We then attempt an inversion for a 1-D (spherically symmetric) viscosity structure based on minimizing the misfit to these PGR data. Using a Monte Carlo algorithm to invert for two layers of viscosity [upper and lower mantle (UM and LM)], we obtain well-constrained values which correspond to the two-layer average of the logarithm of the 3-D viscosity structure in the vicinity of Laurentia. We then attempt to invert for four layers of viscosity. In this case we find a ‘trade-off effect’ in which neighbouring layers may have highly variable viscosities while maintaining a constant average between them. Since the PGR data are insensitive to this trade-off in neighbouring layers, the viscosity of any one of the four layers cannot be well constrained. By repeating the inversion with synthetic data derived from an Earth model which itself is 1-D with four viscosity layers, we demonstrate that it is the insensitivity of the PGR data, not complications due to 3-D structure, that allow the trade-off effect and the resulting failure of the inversion. We also perform a resolution test to study the extent of the insensitivity of PGR to these viscosity trade-offs, finding that the limit of resolution in the UM is about the size of the entire UM, and similarly for the LM. This reinforces our findings that only two layers are obtainable in an inversion for viscosity from PGR.

[1]  J. Mitrovica,et al.  Haskell [1935] revisited , 1996 .

[2]  B. Hager,et al.  Geoid Anomalies in a Dynamic Earth , 1984 .

[3]  W. Peltier,et al.  Mantle viscosity from the simultaneous inversion of multiple data sets pertaining to postglacial rebound , 1996 .

[4]  M. Cheng,et al.  Variations in the Earth's oblateness during the past 28 years , 2004 .

[5]  Erik R. Ivins,et al.  On lateral viscosity contrast in the mantle and the rheology of low-frequency geodynamics , 1995 .

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

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

[8]  M. Sambridge Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[9]  On Postglacial Sea Level , 2007 .

[10]  J. Wahr,et al.  Postglacial rebound and Earth's viscosity structure from GRACE , 2002 .

[11]  Victor Zlotnicki,et al.  Time‐variable gravity from GRACE: First results , 2004 .

[12]  H. Steffen,et al.  Glacial isostatic adjustment of Scandinavia and northwestern Europe and the radial viscosity structure of the Earth's mantle , 2005 .

[13]  W. Peltier Glacial-Isostatic Adjustment—II. The Inverse Problem , 2007 .

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

[15]  W. Peltier,et al.  Constraints on mantle viscosity based upon the inversion of post-glacial uplift data from the Hudson Bay region , 1995 .

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

[17]  Archie Paulson,et al.  Modelling post-glacial rebound with lateral viscosity variations , 2005 .

[18]  E. Ivins,et al.  Lateral viscosity variations beneath Antarctica and their implications on regional rebound motions and seismotectonics , 2004 .

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

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

[21]  Duncan J. Wingham,et al.  A method of combining ICESat and GRACE satellite data to constrain Antarctic mass balance , 2000 .

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

[23]  K. Lambeck Late Devensian and Holocene shorelines of the British Isles and North Sea from models of glacio-hydro-isostatic rebound , 1995, Journal of the Geological Society.

[24]  W. R. Peltier,et al.  Validation of the ICE‐3G Model of Würm‐Wisconsin Deglaciation using a global data base of relative sea level histories , 1992 .

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

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

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

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

[29]  Erik R. Ivins,et al.  Deep mantle viscous structure with prior estimate and satellite constraint , 1993 .

[30]  B. Chao,et al.  Detection of a Large-Scale Mass Redistribution in the Terrestrial System Since 1998 , 2002, Science.

[31]  J. Mitrovica,et al.  A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data , 2004 .

[32]  Jan M. Johansson,et al.  Continuous GPS measurements of postglacial adjustment in Fennoscandia: 2. Modeling results: FENNOSCANDIAN GPS MODELING RESULTS , 2004 .

[33]  J. Mitrovica,et al.  New inferences of mantle viscosity from joint inversion of long‐wavelength mantle convection and post‐glacial rebound data , 1996 .

[34]  W. Peltier,et al.  Ice Age Paleotopography , 1994, Science.

[35]  Bradford H. Hager,et al.  Localization of the gravity field and the signature of glacial rebound , 1997, Nature.

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

[37]  K. Lambeck,et al.  Glacial isostatic adjustment and the radial viscosity profile from inverse modeling , 2002 .

[38]  B. Hager Subducted slabs and the geoid: Constraints on mantle rheology and flow , 1983 .

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

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

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

[42]  L. Cathles,et al.  The Viscosity of the Earth's Mantle , 1975 .

[43]  B. Luzum,et al.  Path of the Mean Rotational Pole From 1899 to 1994 , 1996 .

[44]  G. Kaufmann,et al.  Glacial isostatic adjustment in Fennoscandia with a three-dimensional viscosity structure as an inverse problem , 2002 .

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

[46]  G. Spada,et al.  Mantle viscosity inference: a comparison between simulated annealing and neighbourhood algorithm inversion methods , 2004 .

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

[48]  B. Hager,et al.  Long-wavelength variations in Earth’s geoid: physical models and dynamical implications , 1989, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[49]  J. Johansson,et al.  Continuous GPS measurements of postglacial adjustment in Fennoscandia 1. Geodetic results , 2002 .