Deep mantle structure and the postperovskite phase transition

Apparent Splitting of S Waves Propagating Through an Isotropic Lowermost Mantle

Observations of shear wave anisotropy are key for understanding the mineralogical structure and flow in the mantle. Several researchers have reported the presence of seismic anisotropy in the lowermost 150-250 km of the mantle (i.e., D '' layer), based on differences in the arrival times of vertically (SV) and horizontally (SH) polarized shear waves. By computing waveforms at a period > 6 s for a wide range of 1-D and 3-D Earth structures, we illustrate that a time shift (i.e., apparent splitting) between SV and SH may appear in purely isotropic simulations. This may be misinterpreted as shear wave anisotropy. For near-surface earthquakes, apparent shear wave splitting can result from the interference of S with the surface reflection sS. For deep earthquakes, apparent splitting can be due to the S wave triplication in D '' , reflections off discontinuities in the upper mantle, and 3-D heterogeneity. The wave effects due to anomalous isotropic structure may not be easily distinguished from purely anisotropic effects if the analysis does not involve full waveform simulations.

Melting of peridotite to 140 gigapascals

Interrogating physical processes that occur within the lowermost mantle is a key to understanding Earth's evolution and present-day inner composition. Among such processes, partial melting has been proposed to explain mantle regions with ultralow seismic velocities near the core-mantle boundary, but experimental validation at the appropriate temperature and pressure regimes remains challenging. Using laser-heated diamond anvil cells, we constructed the solidus curve of a natural fertile peridotite between 36 and 140 gigapascals. Melting at core-mantle boundary pressures occurs at 4180 ± 150 kelvin, which is a value that matches estimated mantle geotherms. Molten regions may therefore exist at the base of the present-day mantle. Melting phase relations and element partitioning data also show that these liquids could host many incompatible elements at the base of the mantle.

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.

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

We used three-dimensional inverse scattering of core-reflected shear waves for large-scale, high-resolution exploration of Earth's deep interior (D'') and detected multiple, piecewise continuous interfaces in the lowermost layer (D'') beneath Central and North America. With thermodynamic properties of phase transitions in mantle silicates, we interpret the images and estimate in situ temperatures. A widespread wave-speed increase at 150 to 300 kilometers above the coremantle boundary is consistent with a transition from perovskite to postperovskite. Internal D'' stratification may be due to multiple phase-boundary crossings, and a deep wave-speed reduction may mark the base of a postperovskite lens about 2300 kilometers wide and 250 kilometers thick. The core-mantle boundary temperature is estimated at 3950 +/- 200 kelvin. Beneath Central America, a site of deep subduction, the D'' is relatively cold (DeltaT = 700 +/- 100 kelvin). Accounting for a factor-of-two uncertainty in thermal conductivity, core heat flux is 80 to 160 milliwatts per square meter (mW m(-2)) into the coldest D'' region and 35 to 70 mW m(-2) away from it. Combined with estimates from the central Pacific, this suggests a global average of 50 to 100 mW m(-2) and a total heat loss of 7.5 to 15 terawatts.

Seismological support for the metastable superplume model, sharp features, and phase changes within the lower mantle

Recently, a metastable thermal-chemical convection model was proposed to explain the African Superplume. Its bulk tabular shape remains relatively stable while its interior undergoes significant stirring with low-velocity conduits along its edges and down-welling near the middle. Here, we perform a mapping of chemistry and temperature into P and S velocity variations and replace a seismically derived structure with this hybrid model. Synthetic seismogram sections generated for this 2D model are then compared directly with corresponding seismic observations of P (P, P(C)P, and PKP) and S (S, S(C)S, and SKS) phases. These results explain the anticorrelation between the bulk velocity and shear velocity and the sharpness and level of SKS travel time delays. In addition, we present evidence for the existence of a D" triplication (a putative phase change) beneath the down-welling structure.

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.