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 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.