Stochastic Inversion of P‐to‐S Converted Waves for Mantle Composition and Thermal Structure: Methodology and Application

• Published in: Journal of geophysical research. Solid earth Vol. 123; no. 12; pp. 10,706 - 10,726
• Main Authors: Munch, F. D., Khan, A., Tauzin, B., Zunino, A., Giardini, D.
• Format: Journal Article
• Language: English
• Published: Washington Blackwell Publishing Ltd 01-12-2018; American Geophysical Union
• Subjects: Anchors; Composition; Crustal structure; Data processing; Earth mantle; Elastic properties; Geophysical methods; Geophysics; Inversions; Methodology; Mineralogy; Phase equilibria; Physics; Properties; Sciences of the Universe; Seismic properties; Seismic velocities; Stochasticity; Temperature; Thermal structure; Topography (geology); Transition zone; Upper mantle; Waveforms
• ISSN: 2169-9313; 2169-9356
• DOI: 10.1029/2018JB016032

We present a new methodology for inverting P‐to‐S receiver function (RF) waveforms directly for mantle temperature and composition. This is achieved by interfacing the geophysical inversion with self‐consistent mineral phase equilibria calculations from which rock mineralogy and its elastic properties are predicted as a function of pressure, temperature, and bulk composition. This approach anchors temperatures, composition, seismic properties, and discontinuities that are in mineral physics data, while permitting the simultaneous use of geophysical inverse methods to optimize models of seismic properties to match RF waveforms. Resultant estimates of transition zone (TZ) topography and volumetric seismic velocities are independent of tomographic models usually required for correcting for upper mantle structure. We considered two end‐member compositional models: the equilibrated equilibrium assemblage (EA) and the disequilibrated mechanical mixture (MM) models. Thermal variations were found to influence arrival times of computed RF waveforms, whereas compositional variations affected amplitudes of waves converted at the TZ discontinuities. The robustness of the inversion strategy was tested by performing a set of synthetic inversions in which crustal structure was assumed both fixed and variable. These tests indicate that unaccounted‐for crustal structure strongly affects the retrieval of mantle properties, calling for a two‐step strategy presented herein to simultaneously recover both crustal and mantle parameters. As a proof of concept, the methodology is applied to data from two stations located in the Siberian and East European continental platforms.
Key Points
We present a new tool for inversion of P‐to‐S receiver function waveforms to probe mantle temperature and composition
The method allows for self‐consistent determination of depth, sharpness, and velocity contrast at seismic discontinuities
We account for the effects of seismic processing, crustal structure, and seismic attenuation on receiver function waveforms