Element-by-element parallel spectral-element methods for 3-D acoustic-wave-equation-based teleseismic wave modeling

Shaolin Liu,Dinghui Yang,Xingpeng Dong,Qiancheng Liu,Yongchang Zheng
DOI: https://doi.org/10.5194/se-2017-38
2017-01-01
Solid Earth Discussions
Abstract: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.
What problem does this paper attempt to address?