Pyrite-type FeO2Hx (P phase) has recently been suggested as a possible explanation for ultra-low velocity zones (ULVZs) due to its low seismic velocity and high density. Here we report the results on the congruent melting temperature and melt properties of P phase at high pressures from first-principles molecular dynamics simulations. The results show that P phase would likely be melted near the core-mantle boundary. Liquid FeO2Hx has smaller density and smaller bulk sound velocity compared to the isochemical P phase. As such, small amounts of liquid FeO2Hx could account for the observed seismic anomaly of ULVZs. However, to maintain the liquid FeO2Hx within the ULVZs against compaction requires special physical conditions, such as relatively high viscosity of the solid matrix and/or vigorous convection of the overlying mantle.
more »
« less
First‐Principles Study of FeO2Hx Solid and Melt System at High Pressures: Implications for Ultralow‐Velocity Zones
Pyrite‐type FeO2Hx (P phase) has recently been suggested as a possible alternative to explain ultralow‐velocity zones due to its low seismic velocity and high density. Here we report the results on the congruent melting temperature and melt properties of P phase at high pressures from first‐principles molecular dynamics simulations. The results show that P phase would likely be melted near the core–mantle boundary. Liquid FeO2Hx has smaller density and smaller bulk sound velocity compared to the isochemical P phase. As such, relatively small amounts of liquid FeO2Hx could account for the observed seismic anomaly of ultralow‐velocity zones. However, to maintain the liquid FeO2Hx within the ultralow‐velocity zones against compaction requires special physical conditions, such as relatively high viscosity of the solid matrix and/or vigorous convection of the overlying mantle.
more »
« less
- Award ID(s):
- 1764140
- PAR ID:
- 10100304
- Date Published:
- Journal Name:
- Journal of geophysical research
- ISSN:
- 2156-2202
- Page Range / eLocation ID:
- 1-10
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Qualitative and quantitative analysis of seismic waveforms sensitive to the core–mantle boundary (CMB) region reveal the presence of ultralow-velocity zones (ULVZs) that have a strong decrease in compressional (P) and shear (S) wave velocity, and an increase in density within thin structures. However, understanding their physical origin and relation to the other large-scale structures in the lowermost mantle are limited due to an incomplete mapping of ULVZs at the CMB. The SKS and SPdKS seismic waveforms is routinely used to infer ULVZ presence, but has thus far only been used in a limited epicentral distance range. As the SKS/SPdKS wavefield interacts with a ULVZ it generates additional seismic arrivals, thus increasing the complexity of the recorded wavefield. Here, we explore utilization of the multi-scale sample entropy method to search for ULVZ structures. We investigate the feasibility of this approach through analysis of synthetic seismograms computed for PREM, 1-, 2.5-, and 3-D ULVZs as well as heterogeneous structures with a strong increase in velocity in the lowermost mantle in 1- and 2.5-D. We find that the sample entropy technique may be useful across a wide range of epicentral distances from 100° to 130°. Such an analysis, when applied to real waveforms, could provide coverage of roughly 85% by surface area of the CMB.more » « less
-
Abstract We develop a 3‐D isotropic shear velocity model for the Alaska subduction zone using data from seafloor and land‐based seismographs to investigate along‐strike variations in structure. By applying ambient noise and teleseismic Helmholtz tomography, we derive Rayleigh wave group and phase velocity dispersion maps, then invert them for shear velocity structure using a Bayesian Monte Carlo algorithm. For land‐based stations, we perform a joint inversion of receiver functions and dispersion curves. The forearc crust is relatively thick (35–42 km) and has reduced lower crustal velocities beneath the Kodiak and Semidi segments, which may promote higher seismic coupling. Bristol Bay Basin crust is relatively thin and has a high‐velocity lower layer, suggesting a dense mafic lower crust emplaced by the rifting processes. The incoming plate shows low uppermost mantle velocities, indicating serpentinization. This hydration is more pronounced in the Shumagin segment, with greater velocity reduction extending to 18 ± 3 km depth, compared to the Semidi segment, showing smaller reductions extending to 14 ± 3 km depth. Our estimates of percent serpentinization from VSreduction and VP/VSare larger than those determined using VPreduction in prior studies, likely due to water in cracks affecting VSmore than VP. Revised estimates of serpentinization show that more water subducts than previous studies, and that twice as much mantle water is subducted in the Shumagin segment compared to the Semidi segment. Together with estimates from other subduction zones, the results indicate a wide variation in subducted mantle water between different subduction segments.more » « less
-
Abstract Mid‐ocean ridges generate basalt and harzburgite, which are introduced into the mantle through subduction as a mechanical mixture, contributing to both lateral and radial compositional heterogeneity. The possible accumulation of basalt in the mantle transition zone has been examined, but details of the mantle composition below the 660‐km discontinuity (hereafter d660) remain poorly constrained. In this study, we utilize the subtle waveform details ofS660S, the underside shear‐wave reflection off the d660, to interpret the seismic velocity, density, and compositional structure near, and particularly below, the d660. We identify a significant difference inS660Swaveform shape in subduction zones compared to other regions. The inversion results reveal globally enriched basalt at the d660, with a notably higher content in subduction zones, consistent with the smaller impedance jump andS660Speak amplitude. The basalt fraction decreases significantly to less than 10% near 800‐km depth, forming a global harzburgite‐enriched layer and resulting in a steep seismic velocity gradient just below the d660, in agreement with 1D global reference models. The striking compositional radial variations near the d660 verify geodynamic predictions and challenge the applicability of homogeneous radial compositional models in the mantle. These variations may also affect the viscosity profile and, consequently, the dynamics at the boundary between the upper and lower mantle.more » « less
-
Abstract We analyzed new recordings ofSPdKSseismic waveforms from a global set of broadband seismograms and horizontal tiltmeters from the Hi‐net array in Japan from 26 earthquakes in the Central American region. The anomalous waveforms are consistent with the presence of at least three ultralow‐velocity zones (ULVZs), on the core‐mantle boundary beneath northern Mexico and the southeastern United States. These ULVZs ring an area of high seismic wave speeds observed in tomographic models that has long been associated with past subduction. Waveform modeling using the PSVaxi method suggests that the ULVZs haveSandPwave velocity decreases of 40% and 10%, respectively. These velocity decreases are likely best explained by a partially molten origin where the melt is generated through melting of mid‐ocean ridge basalt atop the subducted slab.more » « less
An official website of the United States government

