Sensitivity analysis of seismic waveforms to upper-mantle discontinuities using the adjoint method

Abstract

Using spectral-element simulations of wave propagation, we investigated the sensitivity of seismic waveforms, recorded on transverse components, to upper-mantle discontinuities in 1-D and 3-D background models. These sensitivity kernels, or Fréchet derivatives, illustrate the spatial sensitivity to model parameters, of which those for shear wave speed and the surface topography of internal boundaries are discussed in this paper. We focus on the boundaries at 400 and 670 km depth of the mantle transition zone. SS precursors have frequently been used to infer the topography of upper-mantle discontinuities. These seismic phases are underside reflections off these boundaries and are usually analysed in the distance range of 110°–160°. This distance range is chosen to minimize the interference from other waves. We show sensitivity kernels for consecutive time windows at three characteristic epicentral distances within the 110°–160° range. The sensitivity kernels are computed with the adjoint method using synthetic data. From our simulations we can draw three main conclusions: (i) The exact Fréchet derivatives show that in all time windows, and also in those centred on the SS precursors, there is interference from other waves. This explains the difficulty reported in the literature to correct for 3-D shear wave speed perturbations, even if the 3-D structure is perfectly known. (ii) All studies attempting to map the topography of the 400 and 670 km discontinuities to date assume that the traveltimes of SS precursors can be linearly decomposed into a 3-D elastic structure and a topography part. We recently showed that such a linear decomposition is not possible for SS precursors, and the sensitivity kernels presented in this paper explain why. (iii) In agreement with previous work, we show that other parts of the seismograms have greater sensitivity to upper-mantle discontinuities than SS precursors, especially multiply bouncing S waves exploiting the S-wave triplications due to the mantle transition zone. These phases can potentially improve the inference of global topographic variations of the upper-mantle discontinuities in the context of full waveform inversion in a joint inversion for (an)elastic parameters and topography.