Abstract:We present an extension of the method recently introduced by Zhao & Chevrot for calculating Frechet kernels from a precomputed database of strain Green's tensors by normal mode summation. The extension involves two aspects: (1) we compute the strain Green's tensors using the Direct Solution Method, which allows us to go up to frequencies as high as 1 Hz; and (2) we develop a spatial interpolation scheme so that the Green's tensors can be computed with a relatively coarse grid, thus improving the efficiency in the computation of the sensitivity kernels. The only requirement is that the Green's tensors be computed with a fine enough spatial sampling rate to avoid spatial aliasing. The Green's tensors can then be interpolated to any location inside the Earth, avoiding the need to store and retrieve strain Green's tensors for a fine sampling grid. The interpolation scheme not only significantly reduces the CPU time required to calculate the Green's tensor database and the disk space to store it, but also enhances the efficiency in computing the kernels by reducing the number of I/O operations needed to retrieve the Green's tensors. Our new implementation allows us to calculate sensitivity kernels for high-frequency teleseismic body waves with very modest computational resources such as a laptop. We illustrate the potential of our approach for seismic tomography by computing traveltime and amplitude sensitivity kernels for high frequency P, PKP and Pdiff phases. A comparison of our PKP kernels with those computed by asymptotic ray theory clearly shows the limits of the latter. With ray theory, it is not possible to model waves diffracted by internal discontinuities such as the coremantle boundary, and it is also difficult to compute amplitudes for paths close to the B-caustic of the PKP phase. We also compute waveform partial derivatives for different parts of the seismic wavefield, a key ingredient for high resolution imaging by waveform inversion. Our computations of partial derivatives in the time window where PcP precursors are commonly observed show that the distribution of sensitivity is complex and counter-intuitive, with a large contribution from the mid-mantle region. This clearly emphasizes the need to use accurate and complete partial derivatives in waveform inversion.

Three-dimensional sensitivity kernels for surface wave observables

Multimode Surface Wave Sensitivity Kernels in Radially Anisotropic Earth Media

Surface-wave group-delay and attenuation kernels

3D Sensitivity Kernels with Full Attenuation Computed by a Combination of the Strong Stability Preserving Runge-Kutta Method and the Scattering-Integral Method

Surface wave sensitivity: mode summation versus adjoint SEM

Three-dimensional Frechet Differential Kernels for Seismicdelay Times

Finite-Frequency Effects in Global Surface-Wave Tomography

Numerical Investigation of 3D multichannel analysis of surface wave method

3D Nearly Analytic Central Difference Method for Computation of Sensitivity Kernels of Wave‐Equation‐Based Seismic Tomography

Finite-Frequency Structural Sensitivities of Short-Period Compressional Body Waves

Structural Sensitivities of Finite-Frequency Seismic Waves: a Full-Wave Approach

Finite-Frequency Sensitivity Kernels for Head Waves

Travelling-wave Green Tensor and Near-Field Rayleigh-wave Sensitivity

Investigation and Use of Surface-wave Characteristics for Near-surface Applications

An Efficient and Flexible Approach to the Calculation of Three-Dimensional Full-Wave Fréchet Kernels for Seismic Tomography—i. Theory

Surface-wave sensitivity to 3-D anelasticity

Sensitivities of Phase-Velocity Dispersion Curves of Surface Waves Due to High-Velocity-Layer and Low-Velocity-Layer Models

Adjoint sensitivity kernels for free oscillation spectra

Sensitivity of frequency-dependent traveltimes to laterally heterogeneous, anisotropic Earth structure

Estimating Q Factor from Multi-mode Shallow-Seismic Surface Waves

Discussion on Some Practical Equations with Implications to High-Frequency Surface-Wave Techniques