We present ab initio (LDA + U
The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al3+, Fe3+)‐, (Al3+)‐, (Fe2+)‐, and (Fe3+)‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe3+‐, Al3+‐, or (Al3+, Fe3+)‐alloying increases the transition pressure, while Fe2+‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO3, and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D region by analyses of tomographic images.
more » « less- NSF-PAR ID:
- 10401840
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract s c ) studies of high‐temperature and high‐pressure elastic properties of pure as well as iron‐bearing (ferrous, Fe2+, and ferric, Fe3+) and aluminum‐bearing MgSiO3postperovskite, the likely dominant phase in the deep lower mantle of the Earth. Thermal effects are addressed within the quasiharmonic approximation by combining vibrational density of states and static elastic coefficients. Aggregate elastic moduli and sound velocities for the Mg end members are successfully compared with scarce experimental data available. Effects of iron (Fe) and aluminum (Al) substitutions on elastic properties and their pressure and temperature dependence have been thoroughly investigated. At the observed perovskite to postperovskite transition (P = 125 GPa andT = 2,500 K), compressional and shear velocities increase by 0–1% and 1.5–3.75%, respectively. This observation is consistent with some seismic studies of the D′ ′ discontinuity beneath the Caribbean, which suggests that our robust estimates of elastic properties of the postperovskite phase will be very helpful to understand lateral velocity variations in the deep lower mantle region and to constrain its composition and thermal structure. -
more » « less
Abstract Large Low Shear Velocity Provinces (LLSVPs) in the lowermost mantle are key to understanding the chemical composition and thermal structure of the deep Earth, but their origins have long been debated. Bridgmanite, the most abundant lower-mantle mineral, can incorporate extensive amounts of iron (Fe) with effects on various geophysical properties. Here our high-pressure experiments and ab initio calculations reveal that a ferric-iron-rich bridgmanite coexists with an Fe-poor bridgmanite in the 90 mol% MgSiO3–10 mol% Fe2O3system, rather than forming a homogeneous single phase. The Fe3+-rich bridgmanite has substantially lower velocities and a higher
V P /V S ratio than MgSiO3bridgmanite under lowermost-mantle conditions. Our modeling shows that the enrichment of Fe3+-rich bridgmanite in a pyrolitic composition can explain the observed features of the LLSVPs. The presence of Fe3+-rich materials within LLSVPs may have profound effects on the deep reservoirs of redox-sensitive elements and their isotopes. -
more » « less
Abstract Fe‐Al‐bearing bridgmanite may be the dominant host for ferric iron in Earth's lower mantle. Here we report the synthesis of (Mg0.5Fe3+0.5)(Al0.5Si0.5)O3bridgmanite (FA50) with the highest Fe3+‐Al3+coupled substitution known to date. X‐ray diffraction measurements showed that at ambient conditions, the FA50 adopted the LiNbO3structure. Upon compression at room temperature to 18 GPa, it transformed back into the bridgmanite structure, which remained stable up to 102 GPa and 2,600 K. Fitting Birch‐Murnaghan equation of state of FA50 bridgmanite yields
V 0 = 172.1(4) Å3,K 0 = 229(4) GPa withK 0′ = 4(fixed). The calculated bulk sound velocity of the FA50 bridgmanite is ~7.7% lower than MgSiO3bridgmanite, mainly because the presence of ferric iron increases the unit‐cell mass by 15.5%. This difference likely represents the upper limit of sound velocity anomaly introduced by Fe3+‐Al3+substitution. X‐ray emission and synchrotron Mössbauer spectroscopy measurements showed that after laser annealing, ~6% of Fe3+cations exchanged with Al3+and underwent the high‐ to low‐spin transition at 59 GPa. The low‐spin proportion of Fe3+increased gradually with pressure and reached 17–31% at 80 GPa. Since the cation exchange and spin transition in this Fe3+‐Al3+‐enriched bridgmanite do not cause resolvable unit‐cell volume reduction, and the increase of low‐spin Fe3+fraction with pressure occurs gradually, the spin transition would not produce a distinct seismic signature in the lower mantle. However, it may influence iron partitioning and isotopic fractionation, thus introducing chemical heterogeneity in the lower mantle. -
more » « less
Abstract The amount of ferric iron Fe3+in the lower mantle is largely unknown and may be influenced by the disproportionation reaction of ferrous iron Fe2+into metallic Fe and Fe3+triggered by the formation of bridgmanite. Recent work has shown that Fe3+has a strong effect on the density and seismic wave speeds of bridgmanite and the incorporation of impurities such as aluminum. In order to further investigate the effects of ferric iron on mineral behavior at lower mantle conditions, we conducted laser‐heated diamond‐anvil cell (LHDAC) experiments on two sets of samples nearly identical in composition (an aluminum‐rich pyroxenite glass) except for the Fe3+content; with one sample with more Fe3+(“oxidized”: Fe3+/ΣFe ~ 55%) and the other with less Fe3+(“reduced”: Fe3+/ΣFe ~ 11%). We heated the samples to lower mantle conditions, and the resulting assemblages were drastically different between the two sets of samples. For the reduced composition, we observed a multiphase assemblage dominated by bridgmanite and calcium perovskite. In contrast, the oxidized material yielded a single phase of Ca‐bearing bridgmanite. These Al‐rich pyroxenite samples show a difference in density and seismic velocities for these two redox states, where the reduced assemblage is denser than the oxidized assemblage by ~1.5% at the bottom of the lower mantle and slower (bulk sound speed) by ~2%. Thus, heterogeneities of Fe3+content may lead to density and seismic wave speed heterogeneities in Earth's lower mantle.
-
In recent years there have been several attempts to make the link between mineral properties and seismic anisotropy in the D’’ region but have yet to reach consensus with regards to the dynamics in lower mantle minerals that could give rise to the observed seismic anisotropy. Here, we aim to provide further constraints on the observed long wavelength shear velocity patterns seen in seismic tomography studies. We introduce a forward model of deformation in a subducting slab as it impacts the core mantle boundary (D’’ layer) and proceeds to upwelling at the edge of a simulated LLSVP. By implementing the most recent results from atomistic modeling and high-pressure deformation experiments coupled with a 3-dimensional geodynamic model, we compare the microstructural evolution of an aggregate with a pyrolytic composition to the macroscopically observed seismic anisotropy of the lowermost mantle. We account for topotaxial relations in the forward and reverse phase transitions of MgSiO3-perovskite (Pv) to post-perovskite (pPv) within the slab as well as explore the effects introduced by partial melting near the CMB. Comparisons in the two leading candidate deformation mechanisms in the post-perovskite phase, (001) and (010), are compared. In this study we find that the reverse transition (pPv to Pv) occurs at a depth which is ~ 150 km deeper than that of the forward transition due to increasing temperature near the CMB providing a varying topography of the D’’ discontinuity. Our model also produces good fits with the isotropic velocities of PREM for the bulk lower mantle. When coupled with temperature and pressure dependent forward and reverse phase transitions, a pPv system with dominant (001) slip provides good correlation with the currently observed VSH fast horizontal (~ 1 – 6%) in D’’ and with VSV consistently fast in upwelling areas. Azimuthal variations along the streamline are also investigated showing a symmetry lower than that of the assumed VTI in D’’ introduced by ‘rolling’ effects near the slab’s edge. The addition of 1% partial melting at the CMB is shown to increase S and P wave anisotropy beneath the slab at the base of upwelling with up to ~2.5 & 4.0% P and S wave reductions respectively compared to the global reference.more » « less
