Premium
Regularization uncertainty in density models estimated from normal mode data
Author(s) -
Resovsky Joseph S.,
Ritzwoller Michael H.
Publication year - 1999
Publication title -
geophysical research letters
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.007
H-Index - 273
eISSN - 1944-8007
pISSN - 0094-8276
DOI - 10.1029/1999gl900540
Subject(s) - a priori and a posteriori , decorrelation , regularization (linguistics) , normal mode , inverse problem , mathematics , statistical physics , algorithm , computer science , mathematical analysis , physics , artificial intelligence , philosophy , epistemology , quantum mechanics , vibration
Normal mode structure coefficients provide important constraints on the long‐wavelength component of 3‐D mantle density (ρ) structure, but inversions for independent models of υ s , υ p , and ρ using normal mode data alone are ill‐posed even at long wavelengths. Ill‐posed inversions typically are regularized by imposing a priori assumptions on the set of estimated models, but such regularization can introduce important uncertainties in the models. We characterize these uncertainties for ρ models estimated from current normal mode data using a set of 512 different “regularization schemes”. These schemes sample a variety of plausible a priori assumptions about the nature and distribution of mantle heterogeneity by specifying allowable υ s , υ p , ρ, and boundary topography structures. The estimated ρ models are fairly robust with respect to prior constraints on υ s , υ p , and topography. However, the character and amplitude of the estimated ρ models depend strongly on how ρ is allowed to decorrelate from υ s , and we display several models in which ρ and υ s , decorrelate in very different depth intervals. Because these models all result from plausible prior constraints and fit the data equally and acceptably well, inversions of current normal mode data cannot robustly locate the decorrelation of ρ from υ s . It remains possible that reliable ρ models may be obtained in the future as more normal mode measurements are introduced to break the strong tradeoffs between upper and lower mantle ρ structures that characterize current inversions.