Finite-Frequency Kernels Based on Adjoint Methods

We derive the adjoint equations associated with the calculation of Frechet derivatives for tomographic inversions based upon a Lagrange multiplier method. The Frechet derivative of an objective function χ(m), where m denotes the Earth model, may be written in the generic form δχ = ∫ K_m(x) δ ln m(x) d^3x, where δ ln m = δm/m denotes the relative model perturbation and K_m the associated 3D sensitivity or Frechet kernel. Complications due to artificial absorbing boundaries for regional simulations as well as finite sources are accommodated. We construct the 3D finite-frequency “banana-doughnut” kernel K_m by simultaneously computing the so-called “adjoint” wave field forward in time and reconstructing the regular wave field backward in time. The adjoint wave field is produced by using time- reversed signals at the receivers as fictitious, simultaneous sources, while the regular wave field is reconstructed on the fly by propagating the last frame of the wave field, saved by a previous forward simulation, backward in time. The approach is based on the spectral-element method, and only two simulations are needed to produce the 3D finite-frequency sensitivity kernels. The method is applied to 1D and 3D regional models. Various 3D shear- and compressional-wave sensitivity kernels are presented for different regional body- and surface-wave arrivals in the seismograms. These kernels illustrate the sensitivity of the observations to the structural parameters and form the basis of fully 3D tomographic inversions.

[1]  B. Engquist,et al.  Absorbing boundary conditions for acoustic and elastic wave equations , 1977, Bulletin of the Seismological Society of America.

[2]  D. Helmberger,et al.  Modeling the long-period body waves from shallow earthquakes at regional ranges , 1980 .

[3]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[4]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[5]  Alfio Quarteroni,et al.  Generalized Galerkin approximations of elastic waves with absorbing boundary conditions , 1998 .

[6]  Chrysoula Tsogka,et al.  Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic hete , 1998 .

[7]  Guust Nolet,et al.  Three-dimensional sensitivity kernels for finite-frequency traveltimes: the banana–doughnut paradox , 1999 .

[8]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[9]  Egill Hauksson,et al.  Crustal structure and seismicity distribution adjacent to the Pacific and North America plate boundary in southern California , 2000 .

[10]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 2000 .

[11]  Thomas H. Jordan,et al.  Three‐dimensional Fréchet differential kernels for seismicdelay times , 2000 .

[12]  Guust Nolet,et al.  Fréchet kernels for finite‐frequency traveltimes—II. Examples , 2000 .

[13]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[14]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation—I. Validation , 2002 .

[15]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation: II. Three-dimensional models, oceans, rotation and self-gravitation , 2002 .

[16]  3-D Sensitivity Kernels for Surface Wave Observables , 2002 .

[17]  Chen Ji,et al.  A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[18]  Jeroen Tromp,et al.  A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation , 2003 .

[19]  Jean-Pierre Vilotte,et al.  Solving elastodynamics in a fluid-solid heterogeneous sphere: a parallel spectral element approximation on non-conforming grids , 2003 .

[20]  J. Shaw.  P-wave seismic velocity structure derived from sonic logs and industry reflection data in the Los An , 2003 .

[21]  Guust Nolet,et al.  Three‐dimensional sensitivity kernels for surface wave observables , 2004 .

[22]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[23]  E. Engdahl,et al.  Finite-Frequency Tomography Reveals a Variety of Plumes in the Mantle , 2004, Science.

[24]  Jeroen Tromp,et al.  Spectral-element moment tensor inversions for earthquakes in Southern California , 2004 .

[25]  D. Komatitsch,et al.  Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method , 2004 .

[26]  J. Vilotte,et al.  The Newmark scheme as velocity–stress time-staggering: an efficient PML implementation for spectral element simulations of elastodynamics , 2005 .

[27]  Kim B. Olsen,et al.  Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models , 2005 .

[28]  J. Tromp,et al.  A Structural VP Model of the Salton Trough, California, and Its Implications for Seismic Hazard , 2006 .

[29]  Carl Tape,et al.  Finite‐frequency tomography using adjoint methods—Methodology and examples using membrane surface waves , 2007 .