Seismic Velocity and Density Jumps Across the 410- and 660-Kilometer Discontinuities

Seismic evidence for global basalt accumulation in the mantle transition zone

The mantle's compositional structure reflects the thermochemical evolution of Earth. Yet, even the radial average composition of the mantle remains debated. Here, we analyze a global dataset of shear and compressional waves reflecting off the 410- and 660-km discontinuities that is 10 times larger than any previous studies. Our array analysis retrieves globally averaged amplitude-distance trends in SS and PP precursor reflectivity from which we infer relative wavespeed and density contrasts and associated mantle composition. Our results are best matched by a basalt-enriched mantle transition zone, with higher basalt fractions near 660 (~40%) than 410 (~18-31%). These are consistent with mantle-convection/plate-recycling simulations, which predict that basaltic crust accumulates in the mantle transition zone, with basalt fractions peaking near the 660. Basalt segregation in the mantle transition zone also implies that the overall mantle is more silica enriched than the often-assumed pyrolitic mantle reference composition.

Seismic evidence for a 1000 km mantle discontinuity under the Pacific

Seismic discontinuities in the mantle are indicators of its thermo-chemical state and offer clues to its dynamics. Ray-based seismic methods, though limited by the approximations made, have mapped mantle transition zone discontinuities in detail, but have yet to offer definitive conclusions on the presence and nature of mid-mantle discontinuities. Here, we show how to use a wave-equation-based imaging method, reverse-time migration of precursors to surface-reflected seismic body waves, to uncover both mantle transition zone and mid-mantle discontinuities, and interpret their physical nature. We observe a thinned mantle transition zone southeast of Hawaii, and a reduction in impedance contrast around 410 km depth in the same area, suggesting a hotter-than-average mantle in the region. Here, we furthermore reveal a 4000-5000 km-wide reflector in new images of the mid mantle below the central Pacific, at 950-1050 km depth. This deep discontinuity exhibits strong topography and generates reflections with polarity opposite to those originating at the 660 km discontinuity, implying an impedance reversal near 1000 km. We link this mid-mantle discontinuity to the upper reaches of deflected mantle plumes upwelling in the region. Reverse-time migration full-waveform imaging is a powerful approach to imaging Earth's interior, capable of broadening our understanding of its structure and dynamics and shrinking modeling uncertainties.

Earth's rugged lower mantle

Seismic data reveal kilometer-scale topography of the lower-mantle boundary

Inferring Earth's discontinuous chemical layering from the 660-kilometer boundary topography

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.

Sensitivity Kernels of PP Precursor Traveltimes and Their Limitations for Imaging Topography of Discontinuities

We analyze the sensitivity of PP precursor traveltimes that are often used to infer lateral variation in the depths of the 410- and 660-km discontinuities in the mantle. Previous results were inconclusive due to complex wave phenomena, such as multiple energy conversions and focusing/defocusing, that hamper their interpretation. Using spectral-element synthetics and Fréchet derivatives calculated with adjoint methods, we compute sensitivity kernels for volumetric and boundary parameters in a 1-D model for representative epicentral distances of past studies, and a dominant period of 11-25 s. Our results indicate that the boundary sensitivity of PP precursors is low and that these phases are not coherently seen in exact synthetics. Our most important finding is the strong sensitivity to both shear and compressional wave speeds, indicating that wave interference and wave conversions are dominant. The PP precursor traveltimes appear more sensitive to structural parameters, that is, compressional and shear wave speed, than to the boundaries; therefore, they are unlikely sources for valuable insight into discontinuity topography.

Compositional heterogeneity near the base of the mantle transition zone beneath Hawaii

Global seismic discontinuities near 410 and 660 km depth in Earth's mantle are expressions of solid-state phase transitions. These transitions modulate thermal and material fluxes across the mantle and variations in their depth are often attributed to temperature anomalies. Here we use novel seismic array analysis of SS waves reflecting off the 410 and 660 below the Hawaiian hotspot. We find amplitude-distance trends in reflectivity that imply lateral variations in wavespeed and density contrasts across 660 for which thermodynamic modeling precludes a thermal origin. No such variations are found along the 410. The inferred 660 contrasts can be explained by mantle composition varying from average (pyrolitic) mantle beneath Hawaii to a mixture with more melt-depleted harzburgite southeast of the hotspot. Such compositional segregation was predicted, from petrological and numerical convection studies, to occur near hot deep mantle upwellings like the one often invoked to cause volcanic activity on Hawaii.

Global observations of reflectors in the mid-mantle with implications for mantle structure and dynamics

Seismic tomography indicates that flow is commonly deflected in the mid-mantle. However, without a candidate mineral phase change, causative mechanisms remain controversial. Deflection of flow has been linked to radial changes in viscosity and/or composition, but a lack of global observations precludes comprehensive tests by seismically detectable features. Here we perform a systematic global-scale interrogation of mid-mantle seismic reflectors with lateral size 500–2000 km and depths 800–1300 km. Reflectors are detected globally with variable depth, lateral extent and seismic polarity and identify three distinct seismic domains in the mid-mantle. Near-absence of reflectors in seismically fast regions may relate to dominantly subvertical heterogeneous slab material or small impedance contrasts. Seismically slow thermochemical piles beneath the Pacific generate numerous reflections. Large reflectors at multiple depths within neutral regions possibly signify a compositional or textural transition, potentially linked to long-term slab stagnation. This variety of reflector properties indicates widespread compositional heterogeneity at mid-mantle depths.

The Earth’s mantle undergoes changes as temperature and pressure increase with depth. Here, the authors present a global interrogation of reflectors in the Earth’s mid-mantle revealing a significant variation in their properties, with widespread compositional heterogeneity and seismic velocity in the mid-mantle, which signify contrasting styles of mantle flow.

Constraints on the anisotropic contributions to velocity discontinuities at ∼60 km depth beneath the Pacific

Strong, sharp, negative seismic discontinuities, velocity decreases with depth, are observed beneath the Pacific seafloor at ∼60 km depth. It has been suggested that these are caused by an increase in radial anisotropy with depth, which occurs in global surface wave models. Here we test this hypothesis in two ways. We evaluate whether an increase in surface wave radial anisotropy with depth is robust with synthetic resolution tests. We do this by fitting an example surface wave data set near the East Pacific Rise. We also estimate the apparent isotropic seismic velocity discontinuities that could be caused by changes in radial anisotropy in S-to-P and P-to-S receiver functions and SS precursors using synthetic seismograms. We test one model where radial anisotropy is caused by olivine alignment and one model where it is caused by compositional layering. The result of our surface wave inversion suggests strong shallow azimuthal anisotropy beneath 0-10 Ma seafloor, which would also have a radial anisotropy signature. An increase in radial anisotropy with depth at 60 km depth is not well-resolved in surface wave models, and could be artificially observed. Shallow isotropy underlain by strong radial anisotropy could explain moderate apparent velocity drops (<6%) in SS precursor imaging, but not receiver functions. The effect is diminished if strong anisotropy also exists at 0-60 km depth as suggested by surface waves. Overall, an increase in radial anisotropy with depth may not exist at 60 km beneath the oceans and does not explain the scattered wave observations.

Mineral phases of the Earth´s mantle

Our knowledge of the structure of the Earth´s interior has been obtained by analysing seismic waves that travel in the Earth, and the reference Earth global models used by geophysicists are essentially seismological. Depth profiles of the seismic waves velocities reveal that the deep Earth is divided in several shells, separated by velocity and density discontinuities. The main discontinuity located at a depth of 2900 km corresponds to the transition between the mantle and the core. The Earth´s mantle can be further divided into the upper mantle and the lower mantle, with a transition zone characterised by two prominent increases in velocities observed at 410- and 660-km depths. This article will be focused on the mineral phases of the Earth´s mantle. The interpretation of seismological models in terms of chemical composition and temperature relies on the knowledge of the nature, structure and elastic properties of the candidate materials. We will describe to what extent recent advances in experimental mineral physics and X-ray diffraction have yielded essential knowledge on the structure and high-pressure high-temperature behaviour of pertinent materials, and major improvements in our understanding of the chemical and mineralogical composition of the Earth´s mantle.

Indoor seismology by probing the Earth's interior by using sound velocity measurements at high pressures and temperatures

The adiabatic bulk (K(S)) and shear (G) moduli of mantle materials at high pressure and temperature can be obtained directly by measuring compressional and shear wave velocities in the laboratory with experimental techniques based on physical acoustics. We present the application of the current state-of-the-art experimental techniques by using ultrasonic interferometry in conjunction with synchrotron x radiation to study the elasticity of olivine and pyroxenes and their high-pressure phases. By using these updated thermoelasticity data for these phases, velocity and density profiles for a pyrolite model are constructed and compared with radial seismic models. We conclude that pyrolite provides an adequate explanation of the major seismic discontinuities at 410- and 660-km depths, the gradient in the transition zone, as well as the velocities in the lower mantle, if the uncertainties in the modeling and the variations in different seismic models are considered. The characteristics of the seismic scaling factors in response to thermal anomalies suggest that anticorrelations between bulk sound and shear wave velocities, as well as the large positive density anomalies observed in the lower mantle, cannot be explained fully without invoking chemical variations.

The Nature of the 660-Kilometer Discontinuity in Earth's Mantle from Global Seismic Observations of PP Precursors

The 660-kilometer discontinuity, which separates Earth's upper and lower mantle, has been detected routinely on a global scale in underside reflections of precursors to SS shear waves. Here, we report observations of this discontinuity in many different regions, using precursors to compressional PP waves. The apparent absence of such precursors in previous studies had posed major problems for models of mantle composition. We find a complicated structure, showing single and double reflections ranging in depth from 640 to 720 kilometers, that requires the existence of multiple phase transitions at the base of the transition zone. The results are consistent with a pyrolite mantle composition.

Multiple seismic discontinuities near the base of the transition zone in the Earth's mantle

The seismologically defined boundary between the transition zone in the Earth's mantle (410-660 km depth) and the underlying lower mantle is generally interpreted to result from the breakdown of the gamma-spinel phase of olivine to magnesium-perovskite and magnesiowustite. Laboratory measurements of these transformations of olivine have determined that the phase boundary has a negative Clapeyron slope and does indeed occur near pressures corresponding to the base of the transition zone. But a computational study has indicated that, because of the presence of garnet minerals, multiple seismic discontinuities might exist near a depth of 660 km (ref. 4), which would alter the simple negative correlation of changes in temperature with changes in the depth of the phase boundary. In particular, garnet minerals undergo exothermic transformations near this depth, acting to complicate the phase relations and possibly effecting mantle convection processes in some regions. Here we present seismic evidence that supports the existence of such multiple transitions near a depth of 660 km beneath southern California. The observations are consistent with having been generated by garnet transformations coupling with the dissociation of the gamma-spinel phase of olivine. Temperature anomalies calculated from the imaged discontinuity depths--using Clapeyron slopes determined for the various transformations--generally match those predicted from an independent P-wave velocity model of the region.