## Abstract

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 δ_{χ} = ∫_{v} K_{m} δ ln md^{3}r + ∫_{Σ} K_{d} δ ln dd^{2}r + ∫_{ΣFS} K_{d} · ∇^{Σ}δ ln dd^{2}r, 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.

Original language | English (US) |
---|---|

Pages (from-to) | 265-286 |

Number of pages | 22 |

Journal | Geophysical Journal International |

Volume | 174 |

Issue number | 1 |

DOIs | |

State | Published - Jul 2008 |

## All Science Journal Classification (ASJC) codes

- Geophysics
- Geochemistry and Petrology

## Keywords

- Computational seismology
- Inverse theory
- Seismic tomography
- Theoretical seismology
- Wave propagation