Imaging of structure at and near the core-mantle boundary using a generalized radon transform: 2. Statistical inference of singularities

[1] We present the second part of our approach to high-resolution imaging of deep Earth's interfaces with large volumes of broadband, three-component seismograms. We focus on the lowermost mantle, also referred to as D″ region, but the methodology can be applied more generally. The first part describes the generalized radon transform (GRT) of broadband ScS data (comprising main arrival, precursors, and coda). The GRT produces “image gathers,” which represent multiple images of medium contrasts at the same image point near the base of the mantle. With a method for statistical inference we use this redundancy (1) to enhance the GRT images through improved recovery of weak contrasts and through suppression of spurious oscillations in the GRT image gathers and (2) to provide uncertainty estimates that can be used to identify the robust features in the images. Using the image gathers from paper 1 as input, we use mixed effects statistical modeling to produce the best estimates of reflectivity along with their uncertainty. In this framework, random noise in the signal is separated into white and coherent components using the geometry of the (GRT) imaging operators and a generalized cross-validation method. With synthetic data we show that conventional GRT images deteriorate substantially, in some cases to the point at which weak reflectors can no longer be detected, due to effects of uneven sampling, wave phenomena that are not accounted for in the underlying single scattering approximation, or errors in the assumed background wave speed model. We demonstrate that even in these circumstances, statistical analysis can yield adequate estimates of the true model. GRT imaging produces robust images of the core-mantle boundary (CMB) beneath Central America and suggests the presence of several structures in the D″ region, in particular between 100 and 200 and between 270 and 320 km above the CMB proper. Most of these structures are significant at the 1σ (that is, 68%) level, but at 2σ (95%) confidence the images show, at various depths above the CMB, intermittent instead of laterally contiguous features.

[1]  Chong Gu,et al.  Smoothing spline Gaussian regression: more scalable computation via efficient approximation , 2004 .

[2]  L. Tenorio,et al.  Imaging of structure at and near the core mantle boundary using a generalized radon transform: 1. Construction of image gathers , 2006 .

[3]  M. Gurnis,et al.  Evidence for a ubiquitous seismic discontinuity at the base of the mantle , 1999, Science.

[4]  Barry G. Quinn,et al.  A fast efficient technique for the estimation of frequency , 1991 .

[5]  G. Robinson.  That BLUP is a Good Thing: The Estimation of Random Effects , 1991 .

[6]  Yuedong Wang.  Mixed effects smoothing spline analysis of variance , 1998 .

[7]  P. Ma,et al.  Seismostratigraphy and Thermal Structure of Earth's Core-Mantle Boundary Region , 2007, Science.

[8]  E. Vonesh,et al.  Linear and Nonlinear Models for the Analysis of Repeated Measurements , 1996 .

[9]  G. Wahba.  Bayesian "Confidence Intervals" for the Cross-validated Smoothing Spline , 1983 .

[10]  T. Lay,et al.  Deep mantle structure and the postperovskite phase transition. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Yuedong Wang,et al.  Mixed-Effects Smoothing Spline ANOVA , 1998 .

[12]  S. Shim.  Stability of MgSiO3 perovskite in the lower mantle , 2013 .

[13]  J. Bass,et al.  Earth's deep mantle : structure, composition, and evolution , 2005 .

[14]  T. Lay,et al.  A shear velocity discontinuity in the lower mantle , 1983 .

[15]  Y. Ohishi,et al.  Post-Perovskite Phase Transition , 2004 .

[16]  G. Wahba.  Spline models for observational data , 1990 .

[17]  Wensheng Guo,et al.  Functional mixed effects models , 2012, Biometrics.

[18]  Douglas Nychka,et al.  Bayesian Confidence Intervals for Smoothing Splines , 1988 .

[19]  Post-perovskite Phase Transition in MgSiO , 2022 .

[20]  J. Tromp,et al.  Two Views of the Deep Mantle , 1998, Science.

[21]  E. Garnero.  Heterogeneity of the Lowermost Mantle , 2000 .

[22]  Maarten V. de Hoop,et al.  Microlocal analysis of seismic inverse scattering in anisotropic elastic media , 2002 .

[23]  Wensheng Guo.  Functional Mixed Effects Models , 2002 .

[24]  Martin Schimmel,et al.  Noise reduction and detection of weak, coherent signals through phase-weighted stacks , 1997 .

[25]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[26]  A. Oganov,et al.  Theoretical and experimental evidence for a post-perovskite phase of MgSiO3 in Earth's D″ layer , 2004, Nature.

[27]  R. Jeanloz,et al.  Stability and crystal structure of MgSiO3 perovskite to the core‐mantle boundary , 2004 .

[28]  Chong Gu.  Smoothing Spline Anova Models , 2002 .

[29]  Chong Gu,et al.  Optimal smoothing in nonparametric mixed-effect models , 2005 .