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.

A general formulation for computing spherical helioseismic sensitivity kernels while incorporating systematical effects

Numerical evaluation of time-distance helioseismic sensitivity kernels in spherical geometry

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

Interpretation of Helioseismic Travel Times - Sensitivity to Sound Speed, Pressure, Density, and Flows

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

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

Testing volume integrals of travel-time sensitivity kernels for flows

Testing the volume integrals of travel-time sensitivity kernels for flows

Linear Sensitivity of Helioseismic Travel Times to Local Flows

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

Adjoint sensitivity kernels for free oscillation spectra

Finite-Frequency Sensitivity Kernels for Head Waves

Fast and Accurate Spherical Harmonic Transform for Spatio-Temporal Regular Grid Data

Multimode Surface Wave Sensitivity Kernels in Radially Anisotropic Earth Media

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

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

An efficient solver for volumetric scattering based on fast spherical harmonics transforms

Three-dimensional Frechet Differential Kernels for Seismicdelay Times

Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations

Helioseismic Travel-Time Definitions and Sensitivity to Horizontal Flows Obtained From Simulations of Solar Convection

On the Small-Defect Perturbation and Sampling of Heterogeneous Solids