Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

被引:7
|
作者
Adourian, S. [1 ,2 ]
Lyu, C. [1 ]
Masson, Y. [3 ]
Munch, F. [1 ]
Romanowicz, B. [1 ,2 ,4 ]
机构
[1] Univ Calif Berkeley, Berkeley Seismol Lab, Berkeley, CA 94720 USA
[2] Inst Phys Globe, F-75005 Paris 05, France
[3] INRIA, F-64000 Pau, France
[4] Coll France, F-75005 Paris, France
基金
美国国家科学基金会; 瑞士国家科学基金会;
关键词
Numerical modelling; Tomography; Computational seismology; Wave propagation; Wave scattering and diffraction; SPECTRAL-ELEMENT METHOD; SHEAR VELOCITY STRUCTURE; TELESEISMIC BODY WAVES; UPPER-MANTLE; SYNTHETIC SEISMOGRAMS; FORM TOMOGRAPHY; FIELD INJECTION; ADJOINT METHODS; TIME-REVERSAL; ELASTIC-WAVES;
D O I
10.1093/gji/ggac394
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
In previous publications, we presented a general framework, which we called 'box tomography', that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or 'box'. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green's functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the 'box', and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.
引用
收藏
页码:1340 / 1356
页数:17
相关论文
共 35 条
  • [21] A new fast multi-domain BEM to model seismic wave propagation and amplification in 3-D geological structures
    Chaillat, S.
    Bonnet, M.
    Semblat, J. F.
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2009, 177 (02) : 509 - 531
  • [22] Efficient global wave propagation adapted to 3-D structural complexity: a pseudospectral/spectral-element approach
    Leng, Kuangdai
    Nissen-Meyer, Tarje
    van Driel, Martin
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2016, 207 (03) : 1700 - 1721
  • [23] 3-D lithospheric structure beneath southern Tibet from Rayleigh-wave tomography with a 2-D seismic array
    Jiang, Mingming
    Zhou, Shiyong
    Sandvol, Eric
    Chen, Xiaofei
    Liang, Xiaofeng
    Chen, Y. John
    Fan, Wenyuan
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2011, 185 (02) : 593 - 608
  • [24] Numerical simulation of seismic wave propagation in realistic 3-D geo-models with a Fourier pseudo-spectral method
    Klin, Peter
    Priolo, Enrico
    Seriani, Geza
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2010, 183 (02) : 905 - 922
  • [25] A New 3-D FDTD Model of Magnetized Cold Plasma for Modeling Electromagnetic Wave Propagation in the Earth-Ionosphere Waveguide
    Yu, Yaxin
    Simpson, Jamesina J.
    2009 IEEE ANTENNAS AND PROPAGATION SOCIETY INTERNATIONAL SYMPOSIUM AND USNC/URSI NATIONAL RADIO SCIENCE MEETING, VOLS 1-6, 2009, : 1480 - 1483
  • [26] 3-D Modeling of Electromagnetic Wave Propagation in the Uniform Earth-Ionosphere Cavity Using a Commercial FDTD Software Package
    Kwisanga, Christian
    Fourie, Coenrad J.
    IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2017, 65 (06) : 3275 - 3278
  • [27] 3-D body-wave tomography from the seismic ambient noise recorded by a dense array in the Dehdasht area, Iran
    Riahi, Ali
    Shomali, Zaher-Hossein
    Obermann, Anne
    Kamayestani, Ahmad
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2023, 234 (02) : 1332 - 1341
  • [28] Deep structure and tectonics of the northern-central Apennines as seen by regional-scale tomography and 3-D located earthquakes
    Carannante, S.
    Monachesi, G.
    Cattaneo, M.
    Amato, A.
    Chiarabba, C.
    JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2013, 118 (10) : 5391 - 5403
  • [29] High-resolution 3-D P wave attenuation structure of the New Madrid Seismic Zone using local earthquake tomography
    Bisrat, Shishay T.
    DeShon, Heather R.
    Pesicek, Jeremy
    Thurber, Clifford
    JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2014, 119 (01) : 409 - 424
  • [30] High-Frequency Simulations of Global Seismic Wave Propagation Using SPECFEM3D_GLOBE on 62K Processors
    Carrington, Laura
    Komatitsch, Dimitri
    Laurenzano, Michael
    Tikir, Mustafa M.
    Michea, David
    Le Goff, Nicolas
    Snavely, Allan
    Tromp, Jeroen
    INTERNATIONAL CONFERENCE FOR HIGH PERFORMANCE COMPUTING, NETWORKING, STORAGE AND ANALYSIS, 2008, : 1 - +