Pressure stabilizes ferrous iron in bridgmanite under hydrous deep lower mantle conditions

Earth's lower mantle is a potential water reservoir. The physical and chemical properties of the region are in part controlled by the Fe3+/ΣFe ratio and total iron content in bridgmanite. However, the water effect on the chemistry of bridgmanite remains unclear. We carry out laser-heated diamond anvil cell experiments under hydrous conditions and observe dominant Fe2+ in bridgmanite (Mg, Fe)SiO3 above 105 GPa under the normal geotherm conditions corresponding to depth > 2300 km, whereas Fe3+-rich bridgmanite is obtained at lower pressures. We further observe FeO in coexistence with hydrous NiAs-type SiO2 under similar conditions, indicating that the stability of ferrous iron is a combined result of H2O effect and high pressure. The stability of ferrous iron in bridgmanite under hydrous conditions would provide an explanation for the nature of the low-shear-velocity anomalies in the deep lower mantle. In addition, entrainment from a hydrous dense layer may influence mantle plume dynamics and contribute to variations in the redox conditions of the mantle.

Imaging deep-mantle plumbing beneath La Réunion and Comores hot spots: Vertical plume conduits and horizontal ponding zones

Whether the two large low-shear velocity provinces (LLSVPs) at the base of Earth's mantle are wide compact structures extending thousands of kilometers upward or bundles of distinct mantle plumes is the subject of debate. Full waveform shear wave tomography of the deep mantle beneath the Indian Ocean highlights the presence of several separate broad low-velocity conduits anchored at the core-mantle boundary in the eastern part of the African LLSVP, most clearly beneath La Réunion and Comores hot spots. The deep plumbing system beneath these hot spots may also include alternating vertical conduits and horizontal ponding zones, from 1000-km depth to the top of the asthenosphere, reminiscent of dyke and sills in crustal volcanic systems, albeit at a whole-mantle scale.

Insights into the nature of the transition zone from physically constrained inversion of long-period seismic data

Imposing a thermal and compositional significance to the outcome of the inversion of seismic data facilitates their interpretation. Using long-period seismic waveforms and an inversion approach that includes constraints from mineral physics, we find that lateral variations of temperature can explain a large part of the data in the upper mantle. The additional compositional signature of cratons emerges in the global model as well. Above 300 km, we obtain seismic geotherms that span the range of expected temperatures in various tectonic regions. Absolute velocities and gradients with depth are well constrained by the seismic data throughout the upper mantle, except near discontinuities. The seismic data are consistent with a slower transition zone and an overall faster shallow upper mantle, which is not compatible with a homogenous dry pyrolite composition. A gradual enrichment with depth in a garnet-rich component helps to reduce the observed discrepancies. A hydrated transition zone would help to lower the velocities in the transition zone, but it does not explain the seismic structure above it.

Robust normal mode constraints on inner-core anisotropy from model space search

A technique for searching full model space that was applied to measurements of anomalously split normal modes showed a robust pattern of P-wave and S-wave anisotropy in the inner core. The parameter describing P-wave anisotropy changes sign around a radius of 400 kilometers, whereas S-wave anisotropy is small in the upper two-thirds of the inner core and becomes negative at greater depths. Our results agree with observed travel-time anomalies of rays traveling at epicentral distances varying from 150 degrees to 180 degrees. The models may be explained by progressively tilted hexagonal close-packed iron in the upper half of the inner core and could suggest a different iron phase in the center.

Mantle-circulation models with sequential data assimilation: inferring present-day mantle structure from plate-motion histories

Data assimilation is an approach to studying geodynamic models consistent simultaneously with observables and the governing equations of mantle flow. Such an approach is essential in mantle circulation models, where we seek to constrain an unknown initial condition some time in the past, and thus cannot hope to use first-principles convection calculations to infer the flow history of the mantle. One of the most important observables for mantle-flow history comes from models of Mesozoic and Cenozoic plate motion that provide constraints not only on the surface velocity of the mantle but also on the evolution of internal mantle-buoyancy forces due to subducted oceanic slabs. Here we present five mantle circulation models with an assimilated plate-motion history spanning the past 120 Myr, a time period for which reliable plate-motion reconstructions are available. All models agree well with upper- and mid-mantle heterogeneity imaged by seismic tomography. A simple standard model of whole-mantle convection, including a factor 40 viscosity increase from the upper to the lower mantle and predominantly internal heat generation, reveals downwellings related to Farallon and Tethys subduction. Adding 35% bottom heating from the core has the predictable effect of producing prominent high-temperature anomalies and a strong thermal boundary layer at the base of the mantle. Significantly delaying mantle flow through the transition zone either by modelling the dynamic effects of an endothermic phase reaction or by including a steep, factor 100, viscosity rise from the upper to the lower mantle results in substantial transition-zone heterogeneity, enhanced by the effects of trench migration implicit in the assimilated plate-motion history. An expected result is the failure to account for heterogeneity structure in the deepest mantle below 1500 km, which is influenced by Jurassic plate motions and thus cannot be modelled from sequential assimilation of plate motion histories limited in age to the Cretaceous. This result implies that sequential assimilation of past plate-motion models is ineffective in studying the temporal evolution of core-mantle-boundary heterogeneity, and that a method for extrapolating present-day information backwards in time is required. For short time periods (of the order of perhaps a few tens of Myr) such a method exists in the form of crude 'backward' convection calculations. For longer time periods (of the order of a mantle overturn), a rigorous approach to extrapolating information back in time exists in the form of iterative nonlinear optimization methods that carry assimilated information into the past through the use of an adjoint mantle convection model.

Superplumes from the core-mantle boundary to the lithosphere: implications for heat flux

Three-dimensional modeling of upper-mantle anelastic structure reveals that thermal upwellings associated with the two superplumes, imaged by seismic elastic tomography at the base of the mantle, persist through the upper-mantle transition zone and are deflected horizontally beneath the lithosphere. This explains the unique transverse shear wave isotropy in the central Pacific. We infer that the two superplumes may play a major and stable role in supplying heat and horizontal flow to the low-viscosity asthenospheric channel, lubricating plate motions and feeding hot spots. We suggest that more heat may be carried through the core-mantle boundary than is accounted for by hot spot fluxes alone.

Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data

Surface geophysical data that are related to the process of thermal convection in the Earth's mantle provide constraints on the rheological properties and density structure of the mantle. We show that these convection-related data imply the existence of a region of very high effective viscosity near 2,000 km depth. This inference is obtained using a viscous-flow model based on recent high-resolution seismic models of three-dimensional structure in the mantle. The high-viscosity layer near 2,000 km depth results in a re-organization of flow from short to long horizontal length scales, which agrees with seismic tomographic observations of very long wavelength structures in the deep mantle. The high-viscosity region also strongly suppresses flow-induced deformation and convective mixing in the deep mantle. Here we predict compositional and thermal heterogeneity in this region, using viscous-flow calculations based on the new viscosity profile, together with independent mineral physics data. These maps are consistent with the anti-correlation of anomalies in seismic shear and bulk sound velocity in the deep mantle. The maps also show that mega-plumes in the lower mantle below the central Pacific and Africa are, despite the presence of compositional heterogeneity, buoyant and actively upwelling structures.

Sediments at the top of Earth's core

Unusual physical properties at the core-mantle boundary have been inferred from seismic and geodetic observations in recent years. We show how both types of observations can be explained by a layer of silicate sediments, which accumulate at the top of the core as Earth cools. Compaction of the sediments expels most of the liquid iron but leaves behind a small amount of core material, which is entrained in mantle convection and may account for the isotopic signatures of core material in some hot spot plumes. Extraction of light elements from the liquid core also enhances the vigor of convection in the core and may increase the power available to the geodynamo.

Normal-mode and free-Air gravity constraints on lateral variations in velocity and density of Earth's mantle

With the use of a large collection of free-oscillation data and additional constraints imposed by the free-air gravity anomaly, lateral variations in shear velocity, compressional velocity, and density within the mantle; dynamic topography on the free surface; and topography on the 660-km discontinuity and the core-mantle boundary were determined. The velocity models are consistent with existing models based on travel-time and waveform inversions. In the lowermost mantle, near the core-mantle boundary, denser than average material is found beneath regions of upwellings centered on the Pacific Ocean and Africa that are characterized by slow shear velocities. These anomalies suggest the existence of compositional heterogeneity near the core-mantle boundary.

Lateral variations in Compressional/Shear velocities at the base of the mantle

Observations of core-diffracted P (Pdiff) and SH (SHdiff) waves recorded by the Missouri-to-Massachusetts (MOMA) seismic array show that the ratio of compressional (P) seismic velocities to horizontal shear (SH) velocities at the base of the mantle changes abruptly from beneath the mid-Pacific (VP/VS = 1.88, also the value predicted by reference Earth models) to beneath Alaska (VP/VS = 1.83). This change signifies a sudden lateral variation in material properties that may have a mineralogical or textural origin. A textural change could be a result of shear stresses induced during the arrival at the core of ancient lithosphere from the northern Pacific paleotrench.