Regularization uncertainty in density models estimated from normal mode data

Abstract

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 imposinga prioriassumptions 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 plausiblea prioriassumptions 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.