Anelastic sensitivity kernels with parsimonious storage for adjoint tomography and full waveform inversion. Issue 3 (13th June 2016)

• Record Type: Journal Article
• Authors: Komatitsch, Dimitri; Xie, Zhinan; Bozdağ, Ebru; Sales de Andrade, Elliott; Peter, Daniel; Liu, Qinya; Tromp, Jeroen
• Is Part Of: Geophysical journal international. Volume 206:Issue 3(2016:Sep.)
• Page Start: 1467
• Page End: 1478
• Publication Date: 2016-06-13
• Subjects: Numerical solutions -- Tomography -- Seismic attenuation -- Seismic tomography -- Computational seismology -- Wave propagation; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggw224
• Languages: English
• ISSNs: 0956-540X

Abstract

Abstract: We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required time steps is unchanged compared to usual acoustic or elastic approaches. The cost is reduced by a factor of 4/3 compared to the case in which anelasticity is partially accounted for by accommodating the effects of physical dispersion. We validate our technique by performing a test in which we compare the K α sensitivity kernel to the exact kernel obtained by saving the entire forward calculation. This benchmark confirms that our approach is also exact. We illustrate the importance of including full attenuation in the calculation of sensitivity kernels by showing significant differences with physical-dispersion-only kernels.