REVEAL: A Global Full-Waveform Inversion Model

被引:19
作者
Thrastarson, Solvi [1 ]
van Herwaarden, Dirk-Philip [1 ]
Noe, Sebastian [1 ]
Schiller, Carl Josef [1 ]
Fichtner, Andreas [1 ]
机构
[1] Swiss Fed Inst Technol, Dept Earth Sci, Zurich, Switzerland
基金
欧盟地平线“2020”; 欧洲研究理事会;
关键词
SPECTRAL ELEMENT METHOD; SHEAR-VELOCITY MODEL; DISCONTINUOUS GALERKIN METHOD; UPPER-MANTLE STRUCTURE; EARTH STRUCTURE; MOMENT TENSORS; UNSTRUCTURED MESHES; ADJOINT TOMOGRAPHY; ELASTIC-WAVES; TRAVEL-TIME;
D O I
10.1785/0120230273
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
We present REVEAL, a global-scale, transversely isotropic full-waveform inversion model. REVEAL builds upon the earlier construction of the long-wavelength Earth (LOWE) model by lowering the minimum period from 100 to 33 s and by more than doubling the number of included earthquakes to 2366. In the course of 305 quasi-Newton iterations, REVEAL assimilated a total of 6,005,727 unique three-component waveforms. The inversion method rests on the combination of a stochastic mini -batch optimization and wavefield-adapted spectral-element meshes. Although the former naturally exploits redundancies in the data set, the latter reduces the cost of wavefield simulations by reducing the effective dimension of the numerical mesh. As a consequence, the average cost of an iteration in this inversion is only around 0.62% of an iteration that uses the complete data set with a standard cubed-sphere-type mesh. We calculated 3D synthetic seismograms using a graphics processing unit-accelerated spectral-element wave propagation solver, accommodating the effects of anelasticity, topography, bathymetry, ocean loading, and ellipticity. For a diverse range of global wavepaths, REVEAL predicts complete three-component seismograms at 33 s period that have not been included in the inversion. This generalization to unseen data suggests applications of REVEAL in event location and characterization, as well as in ground-motion modeling.
引用
收藏
页码:1392 / 1406
页数:15
相关论文
共 91 条
[1]   Modular and flexible spectral-element waveform modelling in two and three dimensions [J].
Afanasiev, Michael ;
Boehm, Christian ;
van Driel, Martin ;
Krischer, Lion ;
Rietmann, Max ;
May, Dave A. ;
Knepley, Matthew G. ;
Fichtner, Andreas .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2019, 216 (03) :1675-1692
[2]   Foundations for a multiscale collaborative global Earth model [J].
Afanasiev, Michael ;
Peter, Daniel ;
Sager, Korbinian ;
Simute, Saule ;
Ermert, Laura ;
Krischer, Lion ;
Fichtner, Andreas .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2016, 204 (01) :39-58
[3]   DETERMINATION OF 3-DIMENSIONAL SEISMIC STRUCTURE OF LITHOSPHERE [J].
AKI, K ;
CHRISTOFFERSSON, A ;
HUSEBYE, ES .
JOURNAL OF GEOPHYSICAL RESEARCH, 1977, 82 (02) :277-296
[4]  
AKI K, 1974, EOS T AM GEOPHYS UN, V55, P1145
[5]   Estimation of Full Moment Tensors, Including Uncertainties, for Nuclear Explosions, Volcanic Events, and Earthquakes [J].
Alvizuri, Celso ;
Silwal, Vipul ;
Krischer, Lion ;
Tape, Carl .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2018, 123 (06) :5099-5119
[6]  
Balay Satish, 2019, Petsc users manual
[7]  
Biswas R., 2017, SEG Technical Program Expanded Abstracts 2017, P1280
[8]   Parallel 3-D viscoelastic finite difference seismic modelling [J].
Bohlen, T .
COMPUTERS & GEOSCIENCES, 2002, 28 (08) :887-899
[9]   Global adjoint tomography: first-generation model [J].
Bozdag, Ebru ;
Peter, Daniel ;
Lefebvre, Matthieu ;
Komatitsch, Dimitri ;
Tromp, Jeroen ;
Hill, Judith ;
Podhorszki, Norbert ;
Pugmire, David .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2016, 207 (03) :1739-1766
[10]   Towards global earth tomography using the spectral element method: a technique based on source stacking [J].
Capdeville, Y ;
Gung, Y ;
Romanowicz, B .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2005, 162 (02) :541-554