Seismostratigraphy and thermal structure of Earth's core-mantle boundary region

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

Velocity and density characteristics of subducted oceanic crust and the origin of lower-mantle heterogeneities

Seismic heterogeneities detected in the lower mantle were proposed to be related to subducted oceanic crust. However, the velocity and density of subducted oceanic crust at lower-mantle conditions remain unknown. Here, we report ab initio results for the elastic properties of calcium ferrite-type phases and determine the velocities and density of oceanic crust along different mantle geotherms. We find that the subducted oceanic crust shows a large negative shear velocity anomaly at the phase boundary between stishovite and CaCl₂-type silica, which is highly consistent with the feature of mid-mantle scatterers. After this phase transition in silica, subducted oceanic crust will be visible as high-velocity heterogeneities as imaged by seismic tomography. This study suggests that the presence of subducted oceanic crust could provide good explanations for some lower-mantle seismic heterogeneities with different length scales except large low shear velocity provinces (LLSVPs).

Terrestrial Single-Station Analog for Constraining the Martian Core and Deep Interior: Implications for InSight

We used a terrestrial single-station seismometer to quantify the uncertainty of InSight (INterior explorations using Seismic Investigations, Geodesy and Heat Transport) data for determining Martian core size. To mimic Martian seismicity, we formed a catalog using 917 terrestrial earthquakes, from which we randomly selected events. We stacked ScS amplitudes on modeled arrival times and searched for where ScS produced coherent seismic amplitudes. A core detection was defined by a coherent peak with small offset between predicted and user-selected arrival times. Iterating the detection algorithm with varying signal-to-noise (SNR) ranges and quantity of events determined the selection frequency of each model and quantified core depth uncertainty. Increasing the quantity of events reduced core depth uncertainty while increasing the recovery rate, while increasing event SNR had little effect. Including ScS2 multiples increased the recovery rate and reduced core depth uncertainty when we used low quantities of events. The most-frequent core depths varied by back azimuth, suggesting our method is sensitive to the presence of mantle heterogeneities. When we added 1° in source distance errors, core depth uncertainty increased by up to 11 km and recovery rates decreased by <5%. Altering epicentral distances by 25% added ~35 km of uncertainty and reduced recovery rates to <50% in some cases. From these experiments, we estimate that if InSight can detect five events with high location precision (<10 % epicentral distance errors), that there is at least an 88% chance of core depth recovery using ScS alone with uncertainty in core depth approaching 18 km and decreasing as more events are located.

Non-monotonic pressure dependence of the thermal conductivity of boron arsenide

Recent experiments demonstrate that boron arsenide (BAs) is a showcase material to study the role of higher-order four-phonon interactions in affecting heat conduction in semiconductors. Here we use first-principles calculations to identify a phenomenon in BAs and a related material - boron antimonide, that has never been predicted or experimentally observed for any other material: competing responses of three-phonon and four-phonon interactions to pressure rise cause a non-monotonic pressure dependence of thermal conductivity, κ, which first increases similar to most materials and then decreases. The resulting peak in κ shows a strong temperature dependence from rapid strengthening of four-phonon interactions relative to three-phonon processes with temperature. Our results reveal pressure as a knob to tune the interplay between the competing phonon scattering mechanisms in BAs and similar compounds, and provide clear experimental guidelines for observation in a readily accessible measurement regime.

Thermal properties of materials are driven by complex many-body interactions among thermal atomic vibrations called phonons. Here the authors show, from first-principles, that a full description of the unusual thermal behaviour in boron arsenide requires considering often-neglected four-phonon interactions.

On the enigmatic birth of the Pacific Plate within the Panthalassa Ocean

The oceanic Pacific Plate started forming in Early Jurassic time within the vast Panthalassa Ocean that surrounded the supercontinent Pangea, and contains the oldest lithosphere that can directly constrain the geodynamic history of the circum-Pangean Earth. We show that the geometry of the oldest marine magnetic anomalies of the Pacific Plate attests to a unique plate kinematic event that sparked the plate's birth at virtually a point location, surrounded by the Izanagi, Farallon, and Phoenix Plates. We reconstruct the unstable triple junction that caused the plate reorganization, which led to the birth of the Pacific Plate, and present a model of the plate tectonic configuration that preconditioned this event. We show that a stable but migrating triple junction involving the gradual cessation of intraoceanic Panthalassa subduction culminated in the formation of an unstable transform-transform-transform triple junction. The consequent plate boundary reorganization resulted in the formation of a stable triangular three-ridge system from which the nascent Pacific Plate expanded. We link the birth of the Pacific Plate to the regional termination of intra-Panthalassa subduction. Remnants thereof have been identified in the deep lower mantle of which the locations may provide paleolongitudinal control on the absolute location of the early Pacific Plate. Our results constitute an essential step in unraveling the plate tectonic evolution of "Thalassa Incognita" that comprises the comprehensive Panthalassa Ocean surrounding Pangea.

Multiple seismic reflectors in Earth's lowermost mantle

The modern view of Earth's lowermost mantle considers a D″ region of enhanced (seismologically inferred) heterogeneity bounded by the core-mantle boundary and an interface some 150-300 km above it, with the latter often attributed to the postperovskite phase transition (in MgSiO3). Seismic exploration of Earth's deep interior suggests, however, that this view needs modification. So-called ScS and SKKS waves, which probe the lowermost mantle from above and below, respectively, reveal multiple reflectors beneath Central America and East Asia, two areas known for subduction of oceanic plates deep into Earth's mantle. This observation is inconsistent with expectations from a thermal response of a single isochemical postperovskite transition, but some of the newly observed structures can be explained with postperovskite transitions in differentiated slab materials. Our results imply that the lowermost mantle is more complex than hitherto thought and that interfaces and compositional heterogeneity occur beyond the D″ region sensu stricto.

Mineralogical effects on the detectability of the postperovskite boundary

The discovery of a phase transition in Mg-silicate perovskite (Pv) to postperovskite (pPv) at lowermost mantle pressure-temperature (P - T) conditions may provide an explanation for the discontinuous increase in shear wave velocity found in some regions at a depth range of 200 to 400 km above the core-mantle boundary, hereafter the D('') discontinuity. However, recent studies on binary and ternary systems showed that reasonable contents of Fe(2+) and Al for pyrolite increase the thickness (width of the mixed phase region) of the Pv - pPv boundary (400-600 km) to much larger than the D('') discontinuity (≤ 70 km). These results challenge the assignment of the D('') discontinuity to the Pv - pPv boundary in pyrolite (homogenized mantle composition). Furthermore, the mineralogy and composition of rocks that can host a detectable Pv → pPv boundary are still unknown. Here we report in situ measurements of the depths and thicknesses of the Pv → pPv transition in multiphase systems (San Carlos olivine, pyrolitic, and midocean ridge basaltic compositions) at the P - T conditions of the lowermost mantle, searching for candidate rocks with a sharp Pv - pPv discontinuity. Whereas the pyrolitic mantle may not have a seismologically detectable Pv → pPv transition due to the effect of Al, harzburgitic compositions have detectable transitions due to low Al content. In contrast, Al-rich basaltic compositions may have a detectable Pv - pPv boundary due to their distinct mineralogy. Therefore, the observation of the D('') discontinuity may be related to the Pv → pPv transition in the differentiated oceanic lithosphere materials transported to the lowermost mantle by subducting slabs.

Lattice thermal conductivity of MgO at conditions of Earth's interior

Thermal conductivity of the Earth's lower mantle greatly impacts the mantle convection style and affects the heat conduction from the core to the mantle. Direct laboratory measurement of thermal conductivity of mantle minerals remains a technical challenge at the pressure-temperature (P-T) conditions relevant to the lower mantle, and previously estimated values are extrapolated from low P-T data based on simple empirical thermal transport models. By using a numerical technique that combines first-principles electronic structure theory and Peierls-Boltzmann transport theory, we predict the lattice thermal conductivity of MgO, previously used to estimate the thermal conductivity in the Earth, at conditions from ambient to the core-mantle boundary (CMB). We show that our first-principles technique provides a realistic model for the P-T dependence of lattice thermal conductivity of MgO at conditions from ambient to the CMB, and we propose thermal conductivity profiles of MgO in the lower mantle based on geotherm models. The calculated conductivity increases from 15 -20 W/K-m at the 670 km seismic discontinuity to 40 -50 W/K-m at the CMB. This large depth variation in calculated thermal conductivity should be included in models of mantle convection, which has been traditionally studied based on the assumption of constant conductivity.

Temperature profile in the lowermost mantle from seismological and mineral physics joint modeling

The internal structure of the core-mantle boundary (CMB) region of the Earth plays a crucial role in controlling the dynamics and evolution of our planet. We have developed a comprehensive model based on the radial variations of shear velocity in the D'' layer (the base of the lower mantle) and the high-P,T elastic properties of major candidate minerals, including the effects of post-perovskite phase changes. This modeling shows a temperature profile in the lowermost mantle with a CMB temperature of 3,800 +/- 200 K, which suggests that lateral temperature variations of 200-300 K can explain much of the large velocity heterogeneity observed in D''. A single-crossing phase transition model was found to be more favorable in reproducing the observed seismic wave velocity structure than a double-crossing phase transition model.

Between a rock and a hot place: the core-mantle boundary

The boundary between the rocky mantle and iron core, almost 2900 km below the surface, is physically the most significant in the Earth's interior. It may be the terminus for subducted surface material, the source of mantle plumes and a control on the Earth's magnetic field. Its properties also have profound significance for the thermochemical and dynamic evolution of the solid Earth. Evidence from seismology shows that D'' (the lowermost few hundred kilometres of the mantle) has a variety of anomalous features. Understanding the origin of these observations requires an understanding of the elastic and deformation properties of the deep Earth minerals. Core-mantle boundary pressures and temperatures are achievable in the laboratory using diamond anvil cell (DAC) apparatus. Such experiments have led to the recent discovery of a new phase, 'post-perovskite', which may explain many hitherto poorly understood properties of D''. Experimental work is also done using analogue minerals at lower pressures and temperatures; these circumvent some of the limits imposed by the small sample size allowed by the DAC. A considerable contribution also comes from theoretical methods that provide a wealth of otherwise unavailable information, as well as verification and refinement of experimental results. The future of the study of the lowermost mantle will involve the linking of the ever-improving seismic observations with predictions of material properties from theoretical and experimental mineral physics in a quantitative fashion, including simulations of the dynamics of the deep Earth. This has the potential to dispel much of the mystery that still surrounds this remote but important region.

Postperovskite phase equilibria in the MgSiO3-Al2O3 system

We investigate high-P,T phase equilibria of the MgSiO(3)-Al(2)O(3) system by means of the density functional ab initio computation methods with multiconfiguration sampling. Being different from earlier studies based on the static substitution properties with no consideration of Rh(2)O(3)(II) phase, present calculations demonstrate that (i) dissolving Al(2)O(3) tends to decrease the postperovskite transition pressure of MgSiO(3) but the effect is not significant ( approximately -0.2 GPa/mol% Al(2)O(3)); (ii) Al(2)O(3) produces the narrow perovskite+postperovskite coexisting P,T area (approximately 1 GPa) for the pyrolitic concentration (x(Al2O3) approximately 6 mol%), which is sufficiently responsible to the deep-mantle D'' seismic discontinuity; (iii) the transition would be smeared (approximately 4 GPa) for the basaltic Al-rich composition (x(Al2O3) approximately 20 mol%), which is still seismically visible unless iron has significant effects; and last (iv) the perovskite structure spontaneously changes to the Rh(2)O(3)(II) with increasing the Al concentration involving small displacements of the Mg-site cations.

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

Crystal structure and thermoelastic properties of (Mg0.91Fe0.09)SiO3 postperovskite up to 135 GPa and 2,700 K

Intriguing seismic observations have been made for the bottom 400 km of Earth's mantle (the D'' region) over the past few decades, yet the origin of these seismic structures has not been well understood. Recent theoretical calculations have predicted many unusual changes in physical properties across the postperovskite transition, perovskite (Pv) --> postperovskite (PPv), that may provide explanations for the seismic observations. Here, we report measurements of the crystal structure of (Mg(0.91)Fe(0.09))SiO(3)-PPv under quasi-hydrostatic conditions up to the pressure (P)-temperature (T) conditions expected for the core-mantle boundary (CMB). The measured crystal structure is in excellent agreement with the first-principles calculations. We found that bulk sound speed (V(Phi)) decreases by 2.4 +/- 1.4% across the PPv transition. Combined with the predicted shear-wave velocity (V(S)) increase, our measurements indicate that lateral variations in mineralogy between Pv and PPv may result in the anticorrelation between the V(Phi) and V(S) anomalies at the D'' region. Also, density increases by 1.6 +/- 0.4% and Grüneisen parameter decreases by 21 +/- 15% across the PPv transition, which will dynamically stabilize the PPv lenses observed in recent seismic studies.

Structure and dynamics of Earth's lower mantle

Processes within the lowest several hundred kilometers of Earth's rocky mantle play a critical role in the evolution of the planet. Understanding Earth's lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earth's deepest magma chamber.

Deformation of (Mg,Fe)SiO 3 Post-Perovskite and D'' Anisotropy

Polycrystalline (Mg(0.9),Fe(0.1))SiO3 post-perovskite was plastically deformed in the diamond anvil cell between 145 and 157 gigapascals. The lattice-preferred orientations obtained in the sample suggest that slip on planes near (100) and (110) dominate plastic deformation under these conditions. Assuming similar behavior at lower mantle conditions, we simulated plastic strains and the contribution of post-perovskite to anisotropy in the D'' region at the Earth core-mantle boundary using numerical convection and viscoplastic polycrystal plasticity models. We find a significant depth dependence of the anisotropy that only develops near and beyond the turning point of a downwelling slab. Our calculated anisotropies are strongly dependent on the choice of elastic moduli and remain hard to reconcile with seismic observations.

Topology of the postperovskite phase transition and mantle dynamics

The postperovskite (ppv) phase transition occurs in the deep mantle close to the core-mantle boundary (CMB). For this reason, we must include in the dynamical considerations both the Clapeyron slope and the temperature intercept, T(int), which is the temperature of the phase transition at the CMB pressure. For a CMB temperature greater than T(int), there is a double crossing of the phase boundary by the geotherms associated with the descending flow. We have found a great sensitivity of the shape of the ppv surface due to the CMB from variations of various parameters such as the amount of internal heating, the Clapeyron slope, and the temperature intercept. Three-dimensional spherical models of mantle convection that can satisfy the seismological constraints depend on the Clapeyron slope. At moderate value, 8 MPa/K, the best fit is found with a core heat flow amounting for 40% of the total heat budget (approximately equal to 15 TW), whereas for 10 MPa/K the agreement is for a lower core heat flow (20%, approximately equal to 7.5 TW). In all cases, these solutions correspond to a temperature intercept 200 K lower than the CMB temperature. These models have holes of perovskite adjacent to ppv in regions of hot upwellings.