@article{7667cedfe3f24b47a65469868dad6fdb,
title = "Anelastic sensitivity kernels with parsimonious storage for adjoint tomography and full waveform inversion",
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 timedomain 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 inwhichwe 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.",
keywords = "Computational seismology, Numerical solutions, Seismic attenuation, Seismic tomography, Tomography, Wave propagation",
author = "Dimitri Komatitsch and Zhinan Xie and Ebru Bozdağ and {de Andrade}, {Elliott Sales} and Daniel Peter and Qinya Liu and Jeroen Tromp",
note = "Funding Information: ACKNOWLEDGEMENTS We thank Mark Asch, Didier Auroux, C{\'e}dric Bellis, {\'E}lie Bretin, Andreas Fichtner, Josselin Garnier, Thomas Guillet, Ioannis G. Kevrekidis, Bruno Lombard, Vadim Monteiller and William W. Symes for fruitful discussion, and the Computational Infrastructure for Geodynamics (CIG) and Marie Cournille for support.We thank Heiner Igel and an anonymous reviewer for useful comments that improved the manuscript. Part of this work was funded by the Simone and Cino del Duca/Institut de France/French Academy of Sciences Foundation under grant no. 095164, by the European Union Horizon 2020 Marie Curie Action no. 641943 project 'WAVES' of call H2020-MSCA-ITN-2014, by U.S. NSF grant 1112906 and by China NSFC grant 51378479. ZX thanks the China Scholarship Council for financial support during his stay at LMA CNRS, and the continuous support from Prof Liao Zhenpeng. ES and QL were supported by the NSERC G8 Research Councils Initiative on Multilateral Research grant no. 490919 and Discovery grant no. 487237. This work was granted access to the European Partnership for Advanced Computing in Europe (PRACE) under allocation TGCC CURIE no. ra2410, to the French HPC resources of TGCC under allocation no. 2015-gen7165 made by GENCI and of the Aix-Marseille Supercomputing Mesocenter under allocations nos 14b013 and 15b034, to the Oak Ridge Leadership Computing Facility at Oak Ridge National Laboratory, USA, which is supported by the Office of Science of the U.S. Department of Energy under contract no. DE-AC05-00OR22725, and to the Sandybridge cluster at the SciNet HPC Consortium funded by the Canada Foundation for Innovation, the Ontario Research Fund and the University of Toronto Startup Fund. Part of this work waspresented at the GPU'2014 Conference in Roma, Italy, in September 2014. Publisher Copyright: {\textcopyright} The Authors 2016. Published by Oxford University Press on behalf of The Royal Astronomical Society.",
year = "2016",
month = sep,
day = "1",
doi = "10.1093/gji/ggw224",
language = "English (US)",
volume = "206",
pages = "1467--1478",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "3",
}