skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: On the Origin of Small‐Scale Seismic Scatters at 660‐km Depth
Abstract Strong small‐scale seismic scatters (<10 km) have been recently observed at 660 km depth, but their origin remains uncertain. We systematically conduct both high‐resolution 2‐D geodynamic computations that include realistic thermodynamic properties, synthetic seismic waveforms, and insight from shallow seismic observations to explore their origin. We demonstrate that neither short‐term subduction, nor long‐term mechanical mantle mixing processes can produce sufficiently strong heterogeneities to explain the origin of such small‐scale seismic scatters. Instead, the intrinsic heterogeneities inside the oceanic lithosphere which subducts into the mantle transition zone and the uppermost lower mantle can explain the observed short‐wavelength scatter waves.  more » « less
Award ID(s):
2009935
PAR ID:
10383261
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
23
Issue:
12
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Seismic observations indicate that the lowermost mantle above the core-mantle boundary is strongly heterogeneous. Body waves reveal a variety of ultra-low velocity zones (ULVZs), which extend not more than 100 km above the core-mantle boundary and have shear velocity reductions of up to 30 per cent. While the nature and origin of these ULVZs remain uncertain, some have suggested they are evidence of partial melting at the base of mantle plumes. Here we use coupled geodynamic/thermodynamic modelling to explore the hypothesis that present-day deep mantle melting creates ULVZs and introduces compositional heterogeneity in the mantle. Our models explore the generation and migration of melt in a deforming and compacting host rock at the base of a plume in the lowermost mantle. We test whether the balance of gravitational and viscous forces can generate partially molten zones that are consistent with the seismic observations. We find that for a wide range of plausible melt densities, permeabilities and viscosities, lower mantle melt is too dense to be stirred into convective flow and instead sinks down to form a completely molten layer, which is inconsistent with observations of ULVZs. Only if melt is less dense or at most ca. 1 per cent more dense than the solid, or if melt pockets are trapped within the solid, can melt remain suspended in the partial melt zone. In these cases, seismic velocities would be reduced in a cone at the base of the plume. Generally, we find partial melt alone does not explain the observed ULVZ morphologies and solid-state compositional variation is required to explain the anomalies. Our findings provide a framework for testing whether seismically observed ULVZ shapes are consistent with a partial melt origin, which is an important step towards constraining the nature of the heterogeneities in the lowermost mantle and their influence on the thermal, compositional, and dynamical evolution of the Earth. 
    more » « less
  2. Abstract Regularization of seismic inversions has a strong imprint on tomographic images. We analyze recorded and spectral‐element S, Sdiff, and SS waveforms to evaluate the benefit of body‐wave amplitudes in global tomography. L‐curve analysis for S40RTS models with recorded and synthetic waveforms show that SS‐S traveltimes and SS/S amplitude ratios have minima within the same damping parameter range. SS/S ratios for S40RTS and model GLAD‐M25 show the trade‐off between scale‐length and strength of lowermost‐mantle heterogeneities. The recorded SS/Sdiff ratios are lower than predicted by 3D mantle models which may be explained by a decrease in the mean shear velocity by at the lowermost 200 km of the mantle. Our results suggest that SS/S amplitude measurements made with 3D waveforms can be used to constrain damping in linearized inversions, and amplitudes are essential for studying the size of heterogeneities. 
    more » « less
  3. Abstract The Earth's lowermost mantle is characterized by two large low shear velocity provinces (LLSVPs). The regions outside the LLSVPs have been suggested to be strongly influenced by subducted slabs and, therefore, much colder than the LLSVPs. However, localized low‐velocity seismic anomalies have been detected in the subduction‐influenced regions, whose origin remains unclear. Here, three‐dimensional geodynamic calculations are performed, and they show that linear, ridge‐like hot thermal anomalies, or thermal ridges, form in the relatively cold, downwelling regions of the lowermost mantle. Like the formation of Richter rolls due to sublithosphere small‐scale convection (SSC), the thermal ridges form as a result of SSC from the basal thermal boundary layer and they extend in directions parallel to the surrounding mantle flow. The formation of thermal ridges in subduction regions of the lowermost mantle is very sensitive to the thermal structures of the subducted materials, and thermal heterogeneities brought to the bottom of the mantle by subducting slabs greatly promote the formation of thermal ridges. The formation of thermal ridges is also facilitated by the increase of core‐mantle boundary heat flux and vigor of lowermost mantle convection. The thermal ridges may explain the low‐velocity seismic anomalies outside of the LLSVPs in the lowermost mantle. The results suggest that the relatively cold, subduction‐influenced regions of the Earth's lowermost mantle may contain localized hot anomalies. 
    more » « less
  4. Abstract Small‐scale intraslab heterogeneity is well documented seismically in multiple subduction zones, but its nature remains elusive. Previous efforts have been mostly focusing on the scattering strength at intermediate depth (<350 km), without constraining its evolution as a function of depth. Here, we illustrate that the inter‐source interferometry method, which turns deep earthquakes into virtual receivers, can resolve small‐scale intraslab heterogeneity in the mantle transition zone. The interferometric waveform observations in the Japan subduction zone require weak scattering (<1.0%) within the slab below 410 km. Combining with previous studies that suggest high heterogeneity level (∼2.5%) at intermediate depth, we conclude that the small‐scale intraslab heterogeneity weakens as slabs subduct. We suggest that the heterogeneities are caused by intraslab hydrous minerals, and the decrease in their scattering strength with depth reveals processes associated with dehydration of subducting slabs. 
    more » « less
  5. SUMMARY Seismic interrogation of the upper mantle from the base of the crust to the top of the mantle transition zone has revealed discontinuities that are variable in space, depth, lateral extent, amplitude and lack a unified explanation for their origin. Improved constraints on the detectability and properties of mantle discontinuities can be obtained with P-to-S receiver function (Ps-RF) where energy scatters from P to S as seismic waves propagate across discontinuities of interest. However, due to the interference of crustal multiples, uppermost mantle discontinuities are more commonly imaged with lower resolution S-to-P receiver function (Sp-RF). In this study, a new method called CRISP-RF (Clean Receiver-function Imaging using SParse Radon Filters) is proposed, which incorporates ideas from compressive sensing and model-based image reconstruction. The central idea involves applying a sparse Radon transform to effectively decompose the Ps-RF into its underlying wavefield contributions, that is direct conversions, multiples, and noise, based on the phase moveout and coherence. A masking filter is then designed and applied to create a multiple-free and denoised Ps-RF. We demonstrate, using synthetic experiment, that our implementation of the Radon transform using a sparsity-promoting regularization outperforms the conventional least-squares methods and can effectively isolate direct Ps conversions. We further apply the CRISP-RF workflow on real data, including single station data on cratons, common-conversion-point stack at continental margins and seismic data from ocean islands. The application of CRISP-RF to global data sets will advance our understanding of the enigmatic origins of the upper mantle discontinuities like the ubiquitous mid-lithospheric discontinuity and the elusive X-discontinuity. 
    more » « less