Open AccessFinite‐frequency sensitivity kernels for global seismic wave propagation based upon adjoint methods
Author(s) -
Liu Qinya,
Tromp Jeroen
Publication year - 2008
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2008.03798.x
Subject(s) - mathematical analysis , adjoint equation , hydrostatic equilibrium , sensitivity (control systems) , finite element method , mathematics , geometry , physics , partial differential equation , quantum mechanics , electronic engineering , engineering , thermodynamics
SUMMARY We determine adjoint equations and Fréchet kernels for global seismic wave propagation based upon a Lagrange multiplier method. We start from the equations of motion for a rotating, self‐gravitating earth model initially in hydrostatic equilibrium, and derive the corresponding adjoint equations that involve motions on an earth model that rotates in the opposite direction. Variations in the misfit function χ then may be expressed as , where δ ln m =δ m / m denotes relative model perturbations in the volume V , δ ln d denotes relative topographic variations on solid–solid or fluid–solid boundaries Σ, and ∇ Σ δ ln d denotes surface gradients in relative topographic variations on fluid–solid boundaries Σ FS . The 3‐D Fréchet kernel K m determines the sensitivity to model perturbations δ ln m , and the 2‐D kernels K d and K d determine the sensitivity to topographic variations δ ln d . We demonstrate also how anelasticity may be incorporated within the framework of adjoint methods. Finite‐frequency sensitivity kernels are calculated by simultaneously computing the adjoint wavefield forward in time and reconstructing the regular wavefield backward in time. Both the forward and adjoint simulations are based upon a spectral‐element method. We apply the adjoint technique to generate finite‐frequency traveltime kernels for global seismic phases ( P , P diff , PKP , S , SKS , depth phases, surface‐reflected phases, surface waves, etc.) in both 1‐D and 3‐D earth models. For 1‐D models these adjoint‐generated kernels generally agree well with results obtained from ray‐based methods. However, adjoint methods do not have the same theoretical limitations as ray‐based methods, and can produce sensitivity kernels for any given phase in any 3‐D earth model. The Fréchet kernels presented in this paper illustrate the sensitivity of seismic observations to structural parameters and topography on internal discontinuities. These kernels form the basis of future 3‐D tomographic inversions.