The inference of mantle viscosity from an inversion of the Fennoscandian relaxation spectrum

A formal inverse theory for mantle viscosity is here applied to a relaxation spectrum derived from the post-glacial uplift of Fennoscandia. the spectrum represents the set of eigenfrequencies (or inverse decay times) for the fundamental mode of viscous gravitational relaxation between the spherical...

Full description

Bibliographic Details
Published in:Geophysical Journal International
Main Authors: Mitrovica, J. X., Peltier, W. R.
Format: Text
Language:English
Published: Oxford University Press 1993
Subjects:
Online Access:http://gji.oxfordjournals.org/cgi/content/short/114/1/45
https://doi.org/10.1111/j.1365-246X.1993.tb01465.x
Description
Summary:A formal inverse theory for mantle viscosity is here applied to a relaxation spectrum derived from the post-glacial uplift of Fennoscandia. the spectrum represents the set of eigenfrequencies (or inverse decay times) for the fundamental mode of viscous gravitational relaxation between the spherical harmonic degrees 14 to 45 and 65 to 80. Theoretical predictions of the eigenfrequencies are based upon the determination of the zeroes of the secular determinant function derived for a spherically symmetric, self-gravitating, visco-elastic planet. Differential kernels relating shifts in the eigenfrequencies to arbitrary perturbations in the radial viscosity profile (i.e. Fréchet kernels) are computed using the variational principle derived by <cross-ref type="bib" refid="bib25">Peltier (1976)</cross-ref>. the inversions are performed within the framework of non-linear Bayesian inference, and the problem has been parameterized in terms of the logarithm of viscosity. The inversions have yielded a set of robust constraints which all models for the radial viscosity profile below Fennoscandia must satisfy. the a posteriori estimates and variance reduction are found to be insensitive to the a priori variance ascribed to the model layers. the constraints have, furthermore, been summarized into a set of a posteriori estimates of the average model viscosity value in radial regions consistent with the resolving power of the data (which decreases from a radial length scale of approximately 120 km at the base of the lithosphere to 1200 km at 1000 km depth; the data provide essentially no information regarding the mantle rheology below 1200 km depth). For example, for Earth models with a lithospheric thickness (LT) of 100 km, the volumetric average logarithm of viscosity in regions in the depth ranges 1040–400 km, 670–210 km and 235–100 km is constrained to be, respectively, 21.03±0.09, 20.70±0.08 and 20.37±0.19. We have repeated the inversions for a number of assumed lithospheric thicknesses and have found that a ...