Journal cover Journal topic
Solid Earth An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 2.083 IF 2.083
  • IF 5-year<br/> value: 2.138 IF 5-year
    2.138
  • CiteScore<br/> value: 2.14 CiteScore
    2.14
  • SNIP value: 0.783 SNIP 0.783
  • SJR value: 1.039 SJR 1.039
  • IPP value: 1.987 IPP 1.987
  • h5-index value: 20 h5-index 20
doi:10.5194/se-2017-38
© Author(s) 2017. This work is distributed
under the Creative Commons Attribution 3.0 License.
Method article
08 May 2017
Review status
This discussion paper is under review for the journal Solid Earth (SE).
Element-by-element parallel spectral-element methods for 3-D acoustic-wave-equation-based teleseismic wave modeling
Shaolin Liu1, Dinghui Yang1, Xingpeng Dong1, Qiancheng Liu2, and Yongchang Zheng3 1Department of Mathematical Sciences, Tsinghua University, Beijing, 100084, China
2King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
3Department of Liver Surgery, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, 100730, China
Abstract. The increasing demand for the high-resolution imaging of deep lithosphere structures requires the utilization of a teleseismic dataset for waveform inversion. The construction of an efficient algorithm for the teleseismic wavefield modeling is valuable for the calculation of misfit kernels or Fr├ęchet derivatives when the teleseismic waveform is used for adjoint tomography. Here, we introduce an element-by-element parallel spectral-element method (EBE-SEM) for the efficient modeling of teleseismic wavefield propagation in a localized geology model. Under the assumption of the plane wave, the frequency-wavenumber (FK) technique is implemented to compute the boundary wavefield used for constructing the boundary condition of the teleseismic wave incidence. To reduce the memory required for the storage of the boundary wavefield for the incidence boundary condition, an economical strategy is introduced to store the boundary wavefield on the model boundary. The perfectly matched layers absorbing boundary condition (PML ABC) is formulated by the EBE-SEM to absorb the scattered wavefield from the model interior. The misfit kernel (derivatives of the waveform misfit with respect to model parameters) can be easily constructed without extra computational effort for the calculation of the element stiffness matrix per time step during the calculation of the adjoint wavefield. Three synthetic examples demonstrate the validity of EBE-SEM for use in teleseismic wavefield modeling and the misfit kernel calculation.

Citation: Liu, S., Yang, D., Dong, X., Liu, Q., and Zheng, Y.: Element-by-element parallel spectral-element methods for 3-D acoustic-wave-equation-based teleseismic wave modeling, Solid Earth Discuss., doi:10.5194/se-2017-38, in review, 2017.
Shaolin Liu et al.
Shaolin Liu et al.
Shaolin Liu et al.

Viewed

Total article views: 277 (including HTML, PDF, and XML)

HTML PDF XML Total BibTeX EndNote
230 42 5 277 2 3

Views and downloads (calculated since 08 May 2017)

Cumulative views and downloads (calculated since 08 May 2017)

Viewed (geographical distribution)

Total article views: 277 (including HTML, PDF, and XML)

Thereof 275 with geography defined and 2 with unknown origin.

Country # Views %
  • 1

Saved

Discussed

Latest update: 28 May 2017
Publications Copernicus
Download
Short summary
Teleseismic wave adjoint tomography is able to imaging the deep structure of the lithosphere. Here we proposed the element-by-element spectral-element method (EBE-SEM) for efficient teleseismic wave modeling. A detailed discussion of the teleseismic seismic incident boundary condition and PML absorbing boundary condition for the scattered waves are presented. Beside the high-efficiency of EBE-SEM for forward modeling of teleseismic waves, it is able to easily construct misfit kernels.
Teleseismic wave adjoint tomography is able to imaging the deep structure of the lithosphere....
Share