Full-waveform inversion reveals diverse origins of lower mantle positive wave speed anomalies

Determining Earth's structure is paramount to unravel its interior dynamics. Seismic tomography reveals positive wave speed anomalies throughout the mantle that spatially correlate with the expected locations of subducted slabs. This correlation has been widely applied in plate reconstructions and geodynamic modelling. However, global travel-time tomography typically incorporates only a limited number of easily identifiable body wave phases and is therefore strongly dependent on the source-receiver geometry. Here, we show how global full-waveform inversion is less sensitive to source-receiver geometry and reveals numerous previously undetected positive wave speed anomalies in the lower mantle. Many of these previously undetected anomalies are situated below major oceans and continental interiors, with no geologic record of subduction, such as beneath the western Pacific Ocean. Moreover, we find no statistically significant correlation positive anomalies as imaged using full-waveform inversion and past subduction. These findings suggest more diverse origins for these anomalies in Earth's lower mantle, unlocking full-waveform inversion as an indispensable tool for mantle exploration.

Remnant of the late Permian superplume that generated the Siberian Traps inferred from geomagnetic data

Mantle plumes have played a key role in tectonic events such as continental break-up and large magmatic events since at least the formation of Gondwana. However, as their signatures on Earth's surface, many of large igneous provinces have disappeared into the mantle during Earth's long-term evolution, meaning that plume remnants in the mantle are crucial in advancing mantle plume theory and accurately reconstructing Earth history. Here we present an electrical conductivity model for North Asia constructed from geomagnetic data. The model shows a large high-electrical-conductivity anomaly in the mantle transition zone beneath the Siberian Traps at the time of their eruption that we interpret to be a thermal anomaly with trace amounts of melt. This anomaly lies almost directly over an isolated low-seismic-wave-velocity anomaly known as the Perm anomaly. The spatial correlation of our anomaly with the Siberian Traps suggests that it represents a remnant of a superplume that was generated from the Perm anomaly. This plume was responsible for the late Permian Siberian large igneous province. The model strengthens the validity of the mantle plume hypothesis.

Glacial isostatic adjustment: physical models and observational constraints

By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume

We report the discovery of anomalously high ³He/⁴He in “cold” geothermal fluids of central Panama, far from any active volcanoes. Combined with independent constraints from lava geochemistry, mantle source geochemical anomalies in Central America require a Galápagos plume contribution that is not derived from hotspot track recycling. Instead, these signals likely originate from large-scale transport of Galápagos plume material at sublithospheric depths. Mantle flow modeling and geophysical observations further indicate these geochemical anomalies could result from a Galápagos plume-influenced asthenospheric “mantle wind” that is actively “blowing” through a slab window beneath central Panama. The lateral transport of plume material represents a potentially widespread yet underappreciated mechanism that scatters enriched geochemical signatures in mantle domains far from plumes.

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high ³He/⁴He (up to 8.9R_A) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source ³He/⁴He >10.3R_A (and potentially up to 26R_A, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1R_A). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high ³He/⁴He values in central Panama are likely derived from the infiltration of a Galápagos plume–like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a “mantle wind” toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Ab Initio Thermoelasticity of Liquid Iron-Nickel-Light Element Alloys

The earth’s core is thought to be composed of Fe-Ni alloy including substantially large amounts of light elements. Although oxygen, silicon, carbon, nitrogen, sulfur, and hydrogen have been proposed as candidates for the light elements, little is known about the amount and the species so far, primarily because of the difficulties in measurements of liquid properties under the outer core pressure and temperature condition. Here, we carry out massive ab initio computations of liquid Fe-Ni-light element alloys with various compositions under the whole outer core P, T condition in order to quantitatively evaluate their thermoelasticity. Calculated results indicate that Si and S have larger effects on the density of liquid iron than O and H, but the seismological reference values of the outer core can be reproduced simultaneously by any light elements except for C. In order to place further constraints on the outer core chemistry, other information, in particular melting phase relations of iron light elements alloys at the inner core-outer core boundary, are necessary. The optimized best-fit compositions demonstrate that the major element composition of the bulk earth is expected to be CI chondritic for the Si-rich core with the pyrolytic mantle or for the Si-poor core and the (Mg,Fe)SiO3-dominant mantle. But the H-rich core likely causes a distinct Fe depletion for the bulk Earth composition.

The adjoint equations for thermochemical compressible mantle convection: derivation and verification by twin experiments

The adjoint method is an efficient way to obtain gradient information in a mantle convection model relative to past flow structure, allowing one to retrodict mantle flow from observations of the present-day mantle state. While adjoint equations for isochemical mantle flow have been derived for both incompressible and compressible flows, here we extend the method to thermochemical mantle flow models, and present thermochemical adjoint equations in the elastic-liquid approximation. We verify the method with twin experiments, and retrodict the flow history of a thermochemical reference model (reference twin) assuming for the final state, either a consistent thermochemical interpretation, using the thermochemical adjoint equations, or an inconsistent purely thermal interpretation, using the isochemical adjoint equations. The consistent simulation correctly retrodicts the flow evolution of the reference twin. The inconsistent case, instead, restores a false flow history whereby internal buoyancy forces and convectively maintained topography are overestimated. Because the cost function is reduced in either case, our results suggest that the adjoint method can be used to link assumptions on the role of chemical mantle heterogeneity to geologic inferences of dynamic topography, thus providing additional means to test hypotheses on mantle composition and dynamics.

Tidal tomography constrains Earth's deep-mantle buoyancy

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

Density structure of Earth's lowermost mantle from Stoneley mode splitting observations

Advances in our understanding of Earth's thermal evolution and the style of mantle convection rely on robust seismological constraints on lateral variations of density. The large-low-shear-wave velocity provinces (LLSVPs) atop the core-mantle boundary beneath Africa and the Pacific are the largest structures in the lower mantle, and hence severely affect the convective flow. Here, we show that anomalous splitting of Stoneley modes, a unique class of free oscillations that are perturbed primarily by velocity and density variations at the core-mantle boundary, is explained best when the overall density of the LLSVPs is lower than the surrounding mantle. The resolved density variations can be explained by the presence of post-perovskite, chemical heterogeneity or a combination of the two. Although we cannot rule out the presence of a ∼100-km-thick denser-than-average basal structure, our results support the hypothesis that LLSVPs signify large-scale mantle upwelling in two antipodal regions of the mantle.

Elasticity of Ferropericlase across the Spin Crossover in the Earth's Lower Mantle

Knowing the elasticity of ferropericlase across the spin transition can help explain seismic and mineralogical models of the lower-mantle including the origin of seismic heterogeneities in the middle to lowermost parts of the lower mantle1234. However, the effects of spin transition on full elastic constants of ferropericlase remain experimentally controversial due to technical challenges in directly measuring sound velocities under lower-mantle conditions12345. Here we have reliably measured both V_P and V_S of a single-crystal ferropericlase ((Mg_0.92,Fe_0.08)O) using complementary Brillouin Light Scattering and Impulsive Stimulated Light Scattering coupled with a diamond anvil cell up to 96 GPa. The derived elastic constants show drastically softened C₁₁ and C₁₂ within the spin transition at 40–60 GPa while C₄₄ is not affected. The spin transition is associated with a significant reduction of the aggregate V_P/V_S via the aggregate V_P softening because V_S softening does not visibly occur within the transition. Based on thermoelastic modelling along an expected geotherm, the spin crossover in ferropericlase can contribute to 2% reduction in V_P/V_S in a pyrolite mineralogical model in mid lower-mantle. Our results imply that the middle to lowermost parts of the lower-mantle would exhibit enhanced seismic heterogeneities due to the occurrence of the mixed-spin and low-spin ferropericlase.

Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle

Deciphering the origin of seismic velocity heterogeneities in the mantle is crucial to understanding internal structures and processes at work in the Earth. The spin crossover in iron in ferropericlase (Fp), the second most abundant phase in the lower mantle, introduces unfamiliar effects on seismic velocities. First-principles calculations indicate that anticorrelation between shear velocity (VS) and bulk sound velocity (Vφ) in the mantle, usually interpreted as compositional heterogeneity, can also be produced in homogeneous aggregates containing Fp. The spin crossover also suppresses thermally induced heterogeneity in longitudinal velocity (VP) at certain depths but not in VS. This effect is observed in tomography models at conditions where the spin crossover in Fp is expected in the lower mantle. In addition, the one-of-a-kind signature of this spin crossover in the RS/P (∂ ln VS/∂ ln VP) heterogeneity ratio might be a useful fingerprint to detect the presence of Fp in the lower mantle.

Some recent advances in understanding the mineralogy of Earth's deep mantle

Understanding planetary structure and evolution requires a detailed knowledge of the properties of geological materials under the conditions of deep planetary interiors. Experiments under the extreme pressure-temperature conditions of the deep mantle are challenging, and many fundamental properties remain poorly constrained or are inferred only through uncertain extrapolations from lower pressure-temperature states. Nevertheless, the last several years have witnessed a number of new developments in this area, and a broad overview of the current understanding of the Earth's lower mantle is presented here. Some recent experimental and theoretical advances related to the lowermost mantle are highlighted. Measurements of the equation of state and deformation behaviour of (Mg,Fe)SiO3 in the CaIrO3-type (post-perovskite) structure yield insights into the nature of the core-mantle boundary region. Theoretical studies of the behaviour of MgSiO3 liquids under high pressure-temperature conditions provide constraints on melt volumes, diffusivities and viscosities that are relevant to understanding both the early Earth (e.g. deep magma oceans) and seismic structure observed in the present Earth (e.g. ultra-low-velocity zones).

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.

Evidence for a chemical-thermal structure at base of mantle from sharp lateral P-wave variations beneath Central America

Compressional waves that sample the lowermost mantle west of Central America show a rapid change in travel times of up to 4 s over a sampling distance of 300 km and a change in waveforms. The differential travel times of the PKP waves (which traverse Earth's core) correlate remarkably well with predictions for S-wave tomography. Our modeling suggests a sharp transition in the lowermost mantle from a broad slow region to a broad fast region with a narrow zone of slowest anomaly next to the boundary beneath the Cocos Plate and the Caribbean Plate. The structure may be the result of ponding of ancient subducted Farallon slabs situated near the edge of a thermal and chemical upwelling.

MgSiO3 postperovskite at D'' conditions

The postperovskite transition in MgSiO(3) at conditions similar to those expected at the D'' discontinuity of Earth's lower mantle offers a paradigm for interpreting the properties of this region. Despite consistent experimental and theoretical predictions of this phase transformation, the complexity of the D'' region raises questions about its geophysical significance. Here we report the thermoelastic properties of Cmcm postperovskite at appropriate conditions and evidences of its presence in the lowermost mantle. These are (i) the jumps in shear and longitudinal velocities similar to those observed in certain places of the D'' discontinuity and (ii) the anticorrelation between shear and bulk velocity anomalies as detected within the D'' region. In addition, the increase in shear modulus across the phase transition provides a possible explanation for the reported discrepancy between the calculated shear modulus of postperovskite free aggregates and the seismological counterpart in the lowermost mantle.

The physical and chemical composition of the lower mantle

This article reviews some of the recent advances made within the field of mineral physics. In order to link the observed seismic and density structures of the lower mantle with a particular mineral composition, knowledge of the thermodynamic properties of the candidate materials is required. Determining which compositional model best matches the observed data is difficult because of the wide variety of possible mineral structures and compositions. State-of-the-art experimental and analytical techniques have pushed forward our knowledge of mineral physics, yet certain properties, such as the elastic properties of lower mantle minerals at high pressures and temperatures, are difficult to determine experimentally and remain elusive. Fortunately, computational techniques are now sufficiently advanced to enable the prediction of these properties in a self-consistent manner, but more results are required.A fundamental question is whether or not the upper and lower mantles are mixing. Traditional models that involve chemically separate upper and lower mantles cannot yet be ruled out despite recent conflicting seismological evidence showing that subducting slabs penetrate deep into the lower mantle and that chemically distinct layers are, therefore, unlikely.Recent seismic tomography studies giving three-dimensional models of the seismic wave velocities in the Earth also base their interpretations on the thermodynamic properties of minerals. These studies reveal heterogeneous velocity and density anomalies in the lower mantle, which are difficult to reconcile with mineral physics data.

Spin transition of iron in magnesiowüstite in the Earth's lower mantle

Iron is the most abundant transition-metal element in the mantle and therefore plays an important role in the geochemistry and geodynamics of the Earth's interior. Pressure-induced electronic spin transitions of iron occur in magnesiowüstite, silicate perovskite and post-perovskite. Here we have studied the spin states of iron in magnesiowüstite and the isolated effects of the electronic transitions on the elasticity of magnesiowüstite with in situ X-ray emission spectroscopy and X-ray diffraction to pressures of the lowermost mantle. An observed high-spin to low-spin transition of iron in magnesiowüstite results in an abnormal compressional behaviour between the high-spin and the low-spin states. The high-pressure, low-spin state exhibits a much higher bulk modulus and bulk sound velocity than the low-pressure, high-spin state; the bulk modulus jumps by approximately 35 percent and bulk sound velocity increases by approximately 15 percent across the transition in (Mg0.83,Fe0.17)O. Although no significant density change is observed across the electronic transition, the jump in the sound velocities and the bulk modulus across the transition provides an additional explanation for the seismic wave heterogeneity in the lowermost mantle. The transition also affects current interpretations of the geophysical and geochemical models using extrapolated or calculated thermal equation-of-state data without considering the effects of the electronic transition.

Shear waves in the diamond-anvil cell reveal pressure-induced instability in (Mg,Fe)O

The emerging picture of Earth's deep interior from seismic tomography indicates more complexity than previously thought. The presence of lateral anisotropy and heterogeneity in Earth's mantle highlights the need for fully anisotropic elasticity data from mineral physics. A breakthrough in high-frequency (gigahertz) ultrasound has resulted in transmission of pure-mode elastic shear waves into a high-pressure diamond-anvil cell using a P-to-S elastic-wave conversion. The full elastic tensor (c(ij)) of high-pressure minerals or metals can be measured at extreme conditions without optical constraints. Here we report the effects of pressure and composition on shear-wave velocities in the major lower-mantle oxide, magnesiowüstite-(Mg,Fe)O. Magnesiowüstite containing more than approximately 50% iron exhibits pressure-induced c(44) shear-mode softening, indicating an instability in the rocksalt structure. The oxide closer to expected lower-mantle compositions ( approximately 20% iron) shows increasing shear velocities more similar to MgO, indicating that it also should have a wide pressure-stability field. A complete sign reversal in the c(44) pressure derivative points to a change in the topology of the (Mg,Fe)O phase diagram at approximately 50-60% iron. The relative stability of Mg-rich (Mg,Fe)O and the strong compositional dependence of shear-wave velocities (and partial differential c(44)/ partial differential P) in (Mg,Fe)O implies that seismic heterogeneity in Earth's lower mantle may result from compositional variations rather than phase changes in (Mg,Fe)O.

Computational mineral physics and the physical properties of perovskite

The inherent uncertainties in modern first-principles calculations are reviewed using geophysically relevant examples. The elastic constants of perovskite at lower-mantle temperatures and pressures are calculated using ab initio molecular dynamics. These are used in conjunction with seismic tomographic models to estimate that the lateral temperature contrasts in the Earth's lower mantle are 800 K at a depth of 1000 km, and 1500 K at a depth of 2000 km. The effect of Al(3+) on the compressibility of MgSiO(3) perovskite is calculated using three different pseudopotentials. The results confirm earlier work and show that the compressibility of perovskites with Al(3+) substituted for both Si(4+) and Mg(2+) is very similar to the compressibility of Al(3+)-free perovskite. Even when 100% of the Si(4+) and Mg(2+) ions are replaced with Al(3+), the bulk modulus is only 7% less than that for Al(3+)-free perovskite. In contrast, perovskites where Al(3+) substitutes for Si(4+) only and that are charge balanced by oxygen vacancies do show higher compressibilities. When corrected to similar concentrations of Al(3+), the calculated compressibilities of the oxygen-vacancy-rich perovskites are in agreement with experimental results.