An official website of the United States government
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
This content will become publicly available on December 1, 2025
Substantial Global Radial Variations of Basalt Content Near the 660‐km Discontinuity
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
- Award ID(s):
- 2123529
- PAR ID:
- 10590380
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 5
- Issue:
- 6
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Seismic tomography has demonstrated that the shear‐wave velocity is relatively high over a 3,000‐km wide region in the lowermost mantle beneath southern and eastern Asia. This seismic anomaly demarcates the current position of slab remnants that may have subducted in the Cretaceous. To further characterize the seismic structure at smaller scales, we measure 929 residual travel time differences (δt) between the phasesScSandSusing recordings of eight earthquakes beneath the Indian Ocean at stations from the Chinese Digital Seismic Network. We interpret variations of δtup to 10 s as due to horizontal shear‐velocity variations in D″ beneath northern India, Nepal, and southwestern China. The shear velocity can vary by as much as 7% over distances shorter than 300 km. Our observations provide additional observational evidence that compositional heterogeneity and possibly melt contribute to the seismic structure of the lower mantle characterized by long‐term subduction and mantle downwelling.more » « less
-
Abstract We investigate broadband SPdKS waveforms from earthquakes occurring beneath Myanmar. These paths sample the core–mantle boundary beneath northwestern China. Waveform modeling shows that two ∼250 × 250 km wide ultra-low velocity zones (ULVZs) with a thickness of roughly 10 km exist in the region. The ULVZ models fitting these data have large S-wave velocity drops of 55% but relatively small 14% P-wave velocity reductions. This is almost a 4:1 S- to P-wave velocity ratio and is suggestive of a partial melt origin. These ULVZs exist in a region of the Circum-Pacific with a long history of subduction and far from large low-velocity province (LLVP) boundaries where ULVZs are more commonly observed. It is possible that these ULVZs are generated by partial melting of mid-ocean ridge basalt.more » « less
-
SUMMARY The detailed structure near the 410-km discontinuity provides key constraints of the dynamic interactions between the upper mantle and the lower mantle through the mantle transition zone (MTZ) via mass and heat exchange. Meanwhile, the temperature of the subducting slab, which can be derived from its fast wave speed perturbation, is critical for understanding the mantle dynamics in subduction zones where the slab enters the MTZ. Multipathing, i.e. triplicated, body waves that bottom near the MTZ carry rich information of the 410-km discontinuity structure and can be used to constrain the discontinuity depth and radial variations of wave speeds across it. In this study, we systematically analysed the trade-off between model parameters in triplication studies using synthetic examples. Specifically, we illustrated the necessity of using array-normalized amplitude. Two 1-D depth profiles of the wave speed below the Tatar Strait of Russia in the Kuril subduction zone are obtained. We have observed triplications due to both the 410-km discontinuity and the slab upper surface. And, seismic structures for these two interfaces are simultaneously inverted. Our derived 410-km discontinuity depths for the northern and southern regions are at 420$$\pm $$15 and 425$$\pm $$15 km, respectively, with no observable uplift. The slab upper surface is inverted to be located about 50–70 km below the 410-km discontinuity. This location is between the depths of the 1 and 2 per cent P-wave speed perturbation contours of a regional 3-D full-waveform inversion (FWI) model, but we found twice the wave speed perturbation amplitude. A wave speed increase of 3.9–4.6 per cent within the slab, compared to 2.0–2.4 per cent from the 3-D FWI model, is necessary to fit the waveforms with the shortest period of 2 s, indicating that high-frequency waves are required to accurately resolve the detailed structures near the MTZ.more » « less
-
Abstract A total of 10,586P‐to‐Sradial receiver functions recorded by 64 broadband seismic stations were utilized to image the 410 and 660 km discontinuities (d410 and d660, respectively) bordering the mantle transition zone (MTZ) beneath the Sumatra Island, the Malay Peninsula, and the western margin of the South China Sea. The d410 and d660 were imaged by stacking receiver functions in successive circular bins with a radius of 1°, after moveout corrections based on the 1‐D IASP91 Earth model. The resulting apparent depths of the discontinuities exhibit significant and spatially systematic variations. The apparent depths of the d410 and d660 range from 382 to 459 km and 637 to 700 km with an average of 406 ± 13 and 670± 12 km, respectively, while the corresponding values for the MTZ thickness are 217 to 295 km and 261 ± 13 km. Underneath southern Sumatra and adjacent regions, the MTZ is characterized by an uplifted d410 and a depressed d660. While the former is probably caused by the low temperature anomaly, the latter is most likely related to a combination of the low temperature anomaly and dehydration associated with the subducted Australian Plate that has reached at least the d660. In contrast, an abnormally thin MTZ is imaged to the southwest of the Toba Caldera. This observation, when combined with results from previous seismic tomography studies, can be explained by advective thermal upwelling through a slab window.more » « less
