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

Handle URI:
http://hdl.handle.net/10754/625887
Title:
Element-by-element parallel spectral-element methods for 3-D teleseismic wave modeling
Authors:
Liu, Shaolin; Yang, Dinghui; Dong, Xingpeng; Liu, Qiancheng; Zheng, Yongchang
Abstract:
The development of an efficient algorithm for teleseismic wave field modeling is valuable for calculating the gradients of the misfit function (termed misfit gradients) 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 wave field propagation in a reduced geology model. Under the plane-wave assumption, the frequency-wavenumber (FK) technique is implemented to compute the boundary wave field used to construct the boundary condition of the teleseismic wave incidence. To reduce the memory required for the storage of the boundary wave field for the incidence boundary condition, a strategy is introduced to efficiently store the boundary wave field on the model boundary. The perfectly matched layers absorbing boundary condition (PML ABC) is formulated using the EBE-SEM to absorb the scattered wave field from the model interior. The misfit gradient can easily be constructed in each time step during the calculation of the adjoint wave field. Three synthetic examples demonstrate the validity of the EBE-SEM for use in teleseismic wave field modeling and the misfit gradient calculation.
KAUST Department:
Physical Sciences and Engineering (PSE) Division
Citation:
Liu S, Yang D, Dong X, Liu Q, Zheng Y (2017) Element-by-element parallel spectral-element methods for 3-D teleseismic wave modeling. Solid Earth 8: 969–986. Available: http://dx.doi.org/10.5194/se-8-969-2017.
Publisher:
Copernicus GmbH
Journal:
Solid Earth
Issue Date:
28-Sep-2017
DOI:
10.5194/se-8-969-2017
Type:
Article
ISSN:
1869-9529
Sponsors:
We greatly appreciate the detailed suggestions from Michal Afanasiev and the anonymous reviewer. Their valuable suggestions greatly improved the quality of the paper. This study was supported by the National Natural Science Foundation of China (grant nos. 41230210 and 41604034). Shaolin Liu was financially supported by the China Postdoctoral Science Foundation (grant no. 2015M580085).
Additional Links:
https://www.solid-earth.net/8/969/2017/
Appears in Collections:
Articles; Physical Sciences and Engineering (PSE) Division

Full metadata record

DC FieldValue Language
dc.contributor.authorLiu, Shaolinen
dc.contributor.authorYang, Dinghuien
dc.contributor.authorDong, Xingpengen
dc.contributor.authorLiu, Qianchengen
dc.contributor.authorZheng, Yongchangen
dc.date.accessioned2017-10-17T11:47:40Z-
dc.date.available2017-10-17T11:47:40Z-
dc.date.issued2017-09-28en
dc.identifier.citationLiu S, Yang D, Dong X, Liu Q, Zheng Y (2017) Element-by-element parallel spectral-element methods for 3-D teleseismic wave modeling. Solid Earth 8: 969–986. Available: http://dx.doi.org/10.5194/se-8-969-2017.en
dc.identifier.issn1869-9529en
dc.identifier.doi10.5194/se-8-969-2017en
dc.identifier.urihttp://hdl.handle.net/10754/625887-
dc.description.abstractThe development of an efficient algorithm for teleseismic wave field modeling is valuable for calculating the gradients of the misfit function (termed misfit gradients) 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 wave field propagation in a reduced geology model. Under the plane-wave assumption, the frequency-wavenumber (FK) technique is implemented to compute the boundary wave field used to construct the boundary condition of the teleseismic wave incidence. To reduce the memory required for the storage of the boundary wave field for the incidence boundary condition, a strategy is introduced to efficiently store the boundary wave field on the model boundary. The perfectly matched layers absorbing boundary condition (PML ABC) is formulated using the EBE-SEM to absorb the scattered wave field from the model interior. The misfit gradient can easily be constructed in each time step during the calculation of the adjoint wave field. Three synthetic examples demonstrate the validity of the EBE-SEM for use in teleseismic wave field modeling and the misfit gradient calculation.en
dc.description.sponsorshipWe greatly appreciate the detailed suggestions from Michal Afanasiev and the anonymous reviewer. Their valuable suggestions greatly improved the quality of the paper. This study was supported by the National Natural Science Foundation of China (grant nos. 41230210 and 41604034). Shaolin Liu was financially supported by the China Postdoctoral Science Foundation (grant no. 2015M580085).en
dc.publisherCopernicus GmbHen
dc.relation.urlhttps://www.solid-earth.net/8/969/2017/en
dc.rightsThis work is distributed under the Creative Commons Attribution 3.0 License.en
dc.rights.urihttps://creativecommons.org/licenses/by/3.0/en
dc.titleElement-by-element parallel spectral-element methods for 3-D teleseismic wave modelingen
dc.typeArticleen
dc.contributor.departmentPhysical Sciences and Engineering (PSE) Divisionen
dc.identifier.journalSolid Earthen
dc.eprint.versionPublisher's Version/PDFen
dc.contributor.institutionDepartment of Mathematical Sciences, Tsinghua University, Beijing, 100084, , Chinaen
dc.contributor.institutionDepartment of Liver Surgery, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Peking Union Medical College, Beijing, 100730, , , Chinaen
kaust.authorLiu, Qianchengen
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.