Three-dimensional models of Earth’s seismic structure can be used to identify temperature-dependent phenomena, including mineralogical phase and spin transformations, that are obscured in 1-D spherical averages. Full-waveform tomography maps seismic wave-speeds inside the Earth in three dimensions, at a higher resolution than classical methods. By providing absolute wave speeds (rather than perturbations) and simultaneously constraining bulk and shear wave speeds over the same frequency range, it becomes feasible to distinguish variations in temperature from changes in composition or spin state. We present a quantitative joint interpretation of bulk and shear wave speeds in the lower mantle, using a recently published full-waveform tomography model. At all depths the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in ferropericlase provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by an enrichment in silica towards the core-mantle boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean.
The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the ~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below ~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.
more » « less- Award ID(s):
- 1918126
- NSF-PAR ID:
- 10308202
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Abstract We present LDA + U sc calculations of high-spin (HS) and low-spin (LS) states in ferropericlase (fp) with an iron concentration of 18.75%. The Hubbard parameter U is determined self-consistently with structures optimized at arbitrary pressures. We confirm a strong dependence of U on the pressure and spin state. Static calculations confirm that the antiferromagnetic configuration is more stable than the ferromagnetic one in the HS state, consistent with low-temperature measurements. Phonon calculations guarantee the dynamical stability of HS and LS states throughout the pressure range of the Earth mantle. Compression curves for HS and LS states agree well with experiments. Using a non-ideal mixing model for the HS to LS states solid solution, we obtain a crossover starting at ∼45 GPa at room temperature and considerably broader than previous results. The spin-crossover phase diagram is calculated, including vibrational, magnetic, electronic, and non-ideal HS–LS entropic contributions. Our results suggest the mixed-spin state predominates in fp in most of the lower mantle.more » « less
-
null (Ed.)Abstract The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth’s lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe 3+ spin crossover (~ 45 GPa), the Debye sound velocity ( v D ) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant v D . Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities ( v S , v P , and v Φ ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low- v Φ and low- v P anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between v S and v Φ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle.more » « less
-
more » « less
SUMMARY We present a new, 3-D model of seismic velocity and anisotropy in the Pacific upper mantle, PAC13E. We invert a data set of single-station surface-wave phase-anomaly measurements sensitive only to Pacific structure for the full set of 13 anisotropic parameters that describe surface-wave anisotropy. Realistic scaling relationships for surface-wave azimuthal anisotropy are calculated from petrological information about the oceanic upper mantle and are used to help constrain the model. The strong age dependence in the oceanic velocities associated with plate cooling is also used as a priori information to constrain the model. We find strong radial anisotropy with vSH > vSV in the upper mantle; the signal peaks at depths of 100–160 km. We observe an age dependence in the depth of peak anisotropy and the thickness of the anisotropic layer, which both increase with seafloor age, but see little age dependence in the depth to the top of the radially anisotropic layer. We also find strong azimuthal anisotropy, which typically peaks in the asthenosphere. The azimuthal anisotropy at asthenospheric depths aligns better with absolute-plate-motion directions while the anisotropy within the lithosphere aligns better with palaeospreading directions. The relative strengths of radial and azimuthal anisotropy are consistent with A-type olivine fabric. Our findings are generally consistent with an explanation in which corner flow at the ridge leads to the development and freezing-in of anisotropy in the lithosphere, and shear between the lithosphere and underlying asthenosphere leads to anisotropy beneath the plate. We also observe large regions within the Pacific basin where the orientation of anisotropy and the absolute-plate-motion direction differ; this disagreement suggests the presence of shear in the asthenosphere that is not aligned with absolute-plate-motion directions. Azimuthal-anisotropy orientation rotates with depth; the depth of the maximum vertical gradient in the fast-axis orientation tends to be age dependent and agrees well with a thermally controlled lithosphere–asthenosphere boundary. We observe that azimuthal-anisotropy strength at shallow depths depends on half-spreading rate, with higher spreading rates associated with stronger anisotropy. Our model implies that corner flow is more efficient at aligning olivine to form lattice-preferred orientation anisotropy fabrics in the asthenosphere when the spreading rate at the ridge is higher.
-
more » « less
Abstract The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases
δ ‐AlOOH (delta phase), MgSiO2(OH)2(phase H), andε ‐FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core‐mantle boundary. We have constrained the equation of state, including the effects of a spin crossover in the Fe3+atoms, of (Al, Fe)‐phase H: Al0.84Fe3+0.07Mg0.02Si0.06OOH, using powder X‐ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)‐phase H andδ ‐(Al, Fe)OOH. The changes in spin state of Fe3+in (Al, Fe)‐phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared withδ ‐(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)‐phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)‐phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid‐lower mantle (∼1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)‐phase H could reproduce some of the seismic signatures of large, low seismic‐velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ ‐AlOOH)‐(MgSiO2(OH)2)‐(ε ‐FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.
