Despite progress in tomographic imaging of Earth's interior, a number of critical questions regarding the large-scale structure and dynamics of the mantle remain outstanding. One of those questions is the impact of phase-boundary undulations on global imaging of mantle heterogeneity and on geodynamic (i.e. convection-related) observables. To address this issue, we developed a joint seismic-geodynamic-mineral physical tomographic inversion procedure that incorporates lateral variations in the depths of the 410- and 660-km discontinuities. This inversion includes S-wave traveltimes, SS precursors that are sensitive to transition-zone topography, geodynamic observables/data (free-air gravity, dynamic surface topography, horizontal divergence of tectonic plates and excess core-mantle boundary ellipticity) and mineral physical constraints on thermal heterogeneity. Compared to joint tomography models that do not include data sensitivity to phase-boundary undulations in the transition zone, the inclusion of 410- and 660-km topography strongly influences the inference of volumetric anomalies in a depth interval that encompasses the transition zone and mid-mantle. It is notable that joint tomography inversions, which include constraints on transition-zone discontinuity topography by seismic and geodynamic data, yield more pronounced density anomalies associated with subduction zones and hotspots. We also find that the inclusion of 410- and 660-km topography may improve the fit to the geodynamic observables, depending on the weights applied to seismic and geodynamic data in the inversions. As a consequence, we find that the amplitude of non-thermal density anomalies required to explain the geodynamic data decreases in most of the mantle. These findings underline the sensitivity of the joint inversions to the inclusion of transition-zone complexity (e.g. phase-boundary topography) and the implications for the inferred non-thermal density anomalies in these depth regions. Finally, we underline that our inferences of 410- and 660-km topography avoid a commonly employed approximation that represents the contribution of volumetric heterogeneity to SS-wave precursor data. Our results suggest that this previously employed correction, based on a priori estimates of upper-mantle heterogeneity, might be a significant source of error in estimating the 410- and 660-km topography.
Inferring Earth’s discontinuous chemical layering from the 660-kilometer boundary topography
Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth’s upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.
more »
« less
- Award ID(s):
- 1736046
- NSF-PAR ID:
- 10108768
- Date Published:
- Journal Name:
- Science
- Volume:
- 363
- Issue:
- 6428
- ISSN:
- 0036-8075
- Page Range / eLocation ID:
- 736 to 740
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
Despite progress in tomographic imaging of Earth’s interior, a number of critical questions regarding the large-scale structure and dynamics of the mantle remain outstanding. One of those questions is the impact of phase-boundary undulations on global imaging of mantle heterogeneity and on geodynamic (i.e. convection-related) observables. To address this issue, we developed a joint seismic-geodynamic-mineral physical tomographic inversion procedure that incorporates lateral variations in the depths of the 410- and 660-km discontinuities. This inversion includes S-wave traveltimes, SS precursors that are sensitive to transition-zone topography, geodynamic observables/data (free-air gravity, dynamic surface topography, horizontal divergence of tectonic plates and excess core-mantle boundary ellipticity) and mineral physical constraints on thermal heterogeneity. Compared to joint tomography models that do not include data sensitivity to phase-boundary undulations in the transition zone, the inclusion of 410- and 660-km topography strongly influences the inference of volumetric anomalies in a depth interval that encompasses the transition zone and mid-mantle. It is notable that joint tomography inversions, which include constraints on transition-zone discontinuity topography by seismic and geodynamic data, yield more pronounced density anomalies associated with subduction zones and hotspots. We also find that the inclusion of 410- and 660-km topography may improve the fit to the geodynamic observables, depending on the weights applied to seismic and geodynamic data in the inversions. As a consequence, we find that the amplitude of non-thermal density anomalies required to explain the geodynamic data decreases in most of the mantle. These findings underline the sensitivity of the joint inversions to the inclusion of transition-zone complexity (e.g. phase-boundary topography) and the implications for the inferred non-thermal density anomalies in these depth regions. Finally, we underline that our inferences of 410- and 660-km topography avoid a commonly employed approximation that represents the contribution of volumetric heterogeneity to SS-wave precursor data. Our results suggest that this previously employed correction, based on a priori estimates of uppermantle heterogeneity, might be a significant source of error in estimating the 410- and 660-km topography.more » « less
-
more » « less
Abstract Topside reverberations off mantle discontinuities are commonly observed at long periods, but their interpretation is complicated because they include both near‐source and near‐receiver reflections. We have developed a method to isolate the stationside reflectors in large data sets with many sources and receivers. Analysis of USArray transverse‐component data from 3,200 earthquakes, using direct
S as a reference phase, shows clear reflections off the 410‐ and 660‐km discontinuities, which can be used to map the depth and brightness of these features. Because our results are sensitive to the impedance contrast (velocity and density), they provide a useful complement to receiver‐function studies, which are primarily sensitive to theS velocity jump alone. In addition, reflectors in our images are more spread out in time than in receiver functions, providing good depth resolution. Our images show strong discontinuities near 410 and 660 km across the entire USArray footprint, with intriguing reflectors at shallower depths in many regions. Overall, the discontinuities in the east appear simpler and more monotonous with a uniform transition zone thickness of 250 km compared to the western United States. In the west, we observe more complex discontinuity topography and small‐scale changes below the Great Basin and the Rocky Mountains, and a decrease in transition‐zone thickness along the western coast. We also observe a dipping reflector in the west that aligns with the top of the high‐velocity Farallon slab anomaly seen in some tomography models, but which also may be an artifact caused by near‐surface scattering of incomingS waves. -
SUMMARY We report finite-frequency imaging of the global 410- and 660-km discontinuities using boundary sensitivity kernels for traveltime measurements made on SS precursors. The application of finite-frequency sensitivity kernels overcomes resolution limits in previous studies associated with large Fresnel zones of SS precursors and their interferences with other seismic phases. In this study, we calculate the finite-frequency sensitivities of SS waves and their precursors based on a single-scattering (Born) approximation in the framework of travelling-wave mode summation. The global discontinuity surface is parametrized using a set of triangular gridpoints with a lateral spacing of about 4°, and we solve the linear finite-frequency inverse problem (2-D tomography) based on singular value decomposition (SVD). The new global models start to show a number of features that were absent (or weak) in ray-theoretical back-projection models at spherical harmonic degree l > 6. The thickness of the mantle transition zone correlates well with wave speed perturbations at a global scale, suggesting dominantly thermal origins for the lateral variations in the mantle transition zone. However, an anticorrelation between the topography of the 410-km discontinuity and wave speed variations is not observed at a global scale. Overall, the mantle transition zone is about 2–3 km thicker beneath the continents than in oceanic regions. The new models of the 410- and 660-km discontinuities show better agreement with the finite-frequency study by Lawrence & Shearer than other global models obtained using SS precursors. However, significant discrepancies between the two models exist in the Pacific Ocean and major subduction zones at spherical harmonic degree >6. This indicates the importance of accounting for wave interactions in the calculations of sensitivity kernels as well as the use of finite-frequency sensitivities in data quality control.more » « less
-
null (Ed.)SUMMARY The method of ScS reverberation migration is based on a ‘common reflection point’ analysis of multiple ScS reflections in the mantle transition zone (MTZ). We examine whether ray-theoretical traveltimes, slownesses and reflection points are sufficiently accurate for estimating the thickness H of the MTZ, defined by the distance between the 410- and 660-km phase transitions. First, we analyse ScS reverberations generated by 35 earthquakes and recorded at hundreds of seismic stations from the combined Arrays in China, Hi-NET in Japan and the Global Seismic Network. This analysis suggests that H varies by about 30 km and therefore that dynamic processes have modified the large-scale structure of the MTZ in eastern Asia and the western Pacific region. Second, we apply the same procedure to spectral-element synthetics for PREM and two 3-D models. One 3-D model incorporates degree-20 topography on the 410 and 660 discontinuities, otherwise preserving the PREM velocity model. The other model incorporates the degree-20 velocity heterogeneity of S20RTS and leaves the 410 and 660 flat. To optimize reflection point coverage, our synthetics were computed assuming a homogeneous grid of stations using 16 events, four of which are fictional. The resolved image using PREM synthetics resembles the PREM structure and indicates that the migration approach is correct. However, ScS reverberations are not as strongly sensitive to H as predicted ray-theoretically because the migration of synthetics for a model with degree-20 topography on the 410 and 660: H varies by less than 5 km in the resolved image but 10 km in the original model. In addition, the relatively strong influence of whole-mantle shear-velocity heterogeneity is evident from the migration of synthetics for the S20RTS velocity model and the broad sensitivity kernels of ScS reverberations at a period of 15 s. A ray-theoretical approach to modelling long-period ScS traveltimes appears inaccurate, at least for continental-scale regions with relatively sparse earthquake coverage. Additional modelling and comparisons with SS precursor and receiver function results should rely on 3-D waveform simulations for a variety of structures and ultimately the implementation of full wave theory.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1126/science.aav0822
