Global adjoint tomography - Model GLAD-M25

Wenjie Lei; Youyi Ruan; Ebru Bozdaǧ; Daniel Peter; Matthieu Lefebvre; Dimitri Komatitsch; Jeroen Tromp; Judith Hill; Norbert Podhorszki; David Pugmire

doi:10.1093/gji/ggaa253

Global adjoint tomography - Model GLAD-M25

Wenjie Lei, Youyi Ruan, Ebru Bozdaǧ, Daniel Peter, Matthieu Lefebvre, Dimitri Komatitsch, Jeroen Tromp, Judith Hill, Norbert Podhorszki, David Pugmire

Research output: Contribution to journal › Article › peer-review

154 Scopus citations

Abstract

Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 'Titan' and the IBM Power 9 'Summit' at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

Original language	English
Pages (from-to)	1-21
Number of pages	21
Journal	Geophysical Journal International
Volume	223
Issue number	1
DOIs	https://doi.org/10.1093/gji/ggaa253
State	Published - Oct 1 2020

Funding

We thank Jeroen Ritsema, Steve Grand and two anonymous reviewers for detailed comments and suggestions which helped to improve an earlier version of the manuscript. This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under contract DE-AC05-00OR22725. Additional computational resources were provided by the Princeton Institute for Computational Science & Engineering (PICSciE). We acknowledge IRIS (iris.edu) and ORFEUS (orfeus-eu.org) for providing the data used in this study. We thank Ryan Modrak, Ridvan Örsvuran, Frederik J. Simons and James Smith for fruitful discussions, and Caio Ciardelli for implementing the spherical harmonic model expansion. The open source spectral-element software package SPECFEM3D GLOBE and the seismicmeasurement software package FLEXWIN used for this article are freely available via the Computational Infrastructure for Geodynamics (CIG; geodynamics.org). This research was supported by NSF grant 1644826.

Keywords

Computational seismology
Seismic anisotropy
Seismic tomography
Wave propagation
Waveform inversion

Access to Document

10.1093/gji/ggaa253

Cite this

@article{f1ea135e781344719f3ab74533a79944,

title = "Global adjoint tomography - Model GLAD-M25",

abstract = "Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fr{\'e}chet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 'Titan' and the IBM Power 9 'Summit' at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.",

keywords = "Computational seismology, Seismic anisotropy, Seismic tomography, Wave propagation, Waveform inversion",

author = "Wenjie Lei and Youyi Ruan and Ebru Bozdaǧ and Daniel Peter and Matthieu Lefebvre and Dimitri Komatitsch and Jeroen Tromp and Judith Hill and Norbert Podhorszki and David Pugmire",

note = "Publisher Copyright: {\textcopyright} 2020 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society.",

year = "2020",

month = oct,

day = "1",

doi = "10.1093/gji/ggaa253",

language = "English",

volume = "223",

pages = "1--21",

journal = "Geophysical Journal International",

issn = "0956-540X",

publisher = "Oxford University Press",

number = "1",

}

TY - JOUR

T1 - Global adjoint tomography - Model GLAD-M25

AU - Lei, Wenjie

AU - Ruan, Youyi

AU - Bozdaǧ, Ebru

AU - Peter, Daniel

AU - Lefebvre, Matthieu

AU - Komatitsch, Dimitri

AU - Tromp, Jeroen

AU - Hill, Judith

AU - Podhorszki, Norbert

AU - Pugmire, David

PY - 2020/10/1

Y1 - 2020/10/1

N2 - Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 'Titan' and the IBM Power 9 'Summit' at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

AB - Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 'Titan' and the IBM Power 9 'Summit' at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

KW - Computational seismology

KW - Seismic anisotropy

KW - Seismic tomography

KW - Wave propagation

KW - Waveform inversion

UR - http://www.scopus.com/inward/record.url?scp=85094129553&partnerID=8YFLogxK

U2 - 10.1093/gji/ggaa253

DO - 10.1093/gji/ggaa253

M3 - Article

AN - SCOPUS:85094129553

SN - 0956-540X

VL - 223

SP - 1

EP - 21

JO - Geophysical Journal International

JF - Geophysical Journal International

IS - 1

ER -

Global adjoint tomography - Model GLAD-M25

Abstract

Funding

Keywords

Access to Document

Other files and links

Fingerprint

Cite this