Weak and vanishing upper mantle discontinuities generated by large-scale lithospheric delamination in the Longmenshan area, China

A large amount of high-quality teleseismic data is used for common conversion point (CCP) stacking of receiver functions in the Longmenshan area. The results show that a large-scale high-velocity anomaly or lithospheric delamination can completely destroy upper mantle discontinuities or erase the phase boundary of olivine, which is a very important finding and can be used to assess stagnant slabs in the mantle transition zone globally. The deepening region of the 660 km discontinuity beneath the Songpan-Ganzi terrane might indicate that the large-scale high-velocity anomaly in the mantle transition zone is a cold domain and can affect the topography of upper mantle discontinuities.

Imaging global mantle discontinuities: a test using full-waveforms and adjoint kernels

Summary We present a novel approach for imaging global mantle discontinuities based on full-waveform inversion (FWI). Over the past decades, extensive research has been done on imaging mantle discontinuities at approximately 400 km and 670 km depth. Accurate knowledge of their topography can put strong constraints on thermal and compositional variations and hence geodynamic modelling. So far, however, there is little consensus on their topography. We present an approach based on adjoint tomography, which has the advantage that Fréchet derivatives for discontinuities and measurements, to be inverted for, are fully consistent. Rather than working with real data, we focus on synthetic tests, where the answer is known in order to be able to evaluate the performance of the developed method. All calculations are based on the community code SPECFEM3D_GLOBE. We generate data in fixed 1-D or 3-D elastic background models of mantle velocity. Our ‘data’ to be inverted contain topography along the 400 km and 670 km mantle discontinuities. To investigate the approach, we perform several tests: (i) In a situation where we know the elastic background model 1-D or 3-D, we recover the target topography fast and accurately, (ii) The exact misfit is not of great importance here, except in terms of convergence speed, similar to a different inverse algorithm, (iii) In a situation where the background model is not known, the convergence is markedly slower, but there is reasonable convergence towards the correct target model of discontinuity topography. It has to be noted that our synthetic test is idealised and in a real data situation, the convergence to and uncertainty of the inferred model is bound to be larger. However, the use of data consistent with Fréchet kernels seems to pay off and might improve our consensus on the nature of mantle discontinuities. Our workflow could be incorporated in future FWI mantle models to adequately infer boundary interface topography.

Validity of Resolving the 785 km Discontinuity in the Lower Mantle with P′P′ Precursors?

P ′ P ′ precursors have been used to detect discontinuities in the lower mantle of the Earth, but some seismic phases propagating along asymmetric ray paths or scattered waves could be misinterpreted as reflections from mantle discontinuities. By forward modeling in standard 1D Earth models, we demonstrate that the frequency content, slowness, and decay with distance of precursors about 180 s before P′P′ arrival are consistent with those of the PKPPdiff phase (or PdiffPKP) at epicentral distances around 78° rather than a reflection from a lower mantle interface. Furthermore, a beamforming technique applied to waveform data recorded at the USArray demonstrates that PKPPdiff can be commonly observed from numerous earthquakes. Hence, a reference 1D Earth model without lower mantle discontinuities can explain many of the observed P′P′ precursors signals if they are interpreted as PKPPdiff, instead of P′785P′. However, this study does not exclude the possibility of 785 km interface beneath the Africa. If this interface indeed exists, P′P′ precursors at distances around 78° would better not be used for its detection to avoid interference from PKPPdiff. Indeed, it could be detected with P′P′ precursors at epicentral distances less than 76° or with other seismic phases such as backscattered PKP·PKP waves.

Imaging with pre-stack migration based on Sp scattering kernels

SUMMARY Sp receiver functions have been widely used to detect the lithosphere–asthenosphere boundary (LAB) and other mantle discontinuities. However, traditional common conversion point (CCP) stacking can be biased by the assumption of horizontal layers and this method typically underestimates scattering amplitudes from velocity boundaries with significant dips. A new pre-stack migration method based on recently developed Sp scattering kernels offers an alternative that more accurately captures the timing and amplitude of scattering. When calculating kernels, Sp-S times are estimated with the fast-marching method, and scattering amplitude versus direction, geometrical spreading and phase shifts are accounted for. To minimize imaging artefacts with larger station spacing, Sp receiver functions are interpolated to more closely spaced pseudo-stations using either compressive sampling or spatial averaging algorithms. To test the kernel-based stacking method, synthetic Sp phases were predicted using SPECFEM2D for velocity models with a flat Moho and a negative mantle velocity gradient with a ramp structure. The kernel-based stacking method resolves horizontal interfaces equally well as CCP stacking and outperforms CCP stacking when imaging boundaries with dips of more than 8°, although dip resolution is still limited. Use of more vertically incident phases such as SKSp improves retrieval of dipping discontinuity segments. A second approach is to down-weight the portions of the kernels that have the greatest positive interference among neighbouring stations, thus enhancing scattering from dipping structures where positive interference is lower. With this downweighting, the kernel-based stacking method applied to Sp data is able to continuously resolve LAB discontinuities with dips up to 15° and to partially resolve continuous LAB discontinuities with dips of ∼20°. The intrinsic properties of teleseismic Sp phase kernels limit their ability to resolve LAB structures with dips of ∼20–35°, but still larger dips of ∼40–50° are resolvable with dense and appropriately placed stations. Analysis of Sp scattering kernels also explains the effectiveness of CCP stacking for quasi-horizontal interfaces.