Cube2sph : A toolkit enabling flexible and accurate continental-scale seismic wave simulations using the SPECFEM3D_Cartesian package

• Published in: Computers & geosciences Vol. 190; p. 105644
• Main Authors: Liu, Tianshi, Wang, Kai, Xie, Yujiang, He, Bin, Lei, Ting, Du, Nanqiao, Tong, Ping, Yang, Yingjie, Rychert, Catherine A., Harmon, Nicholas, Grasselli, Giovanni, Liu, Qinya
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.08.2024
• Subjects: High-performance computing; Perfectly matched layer; Seismology; Spectral-element method; Perfectly matched layer; Spectral-element method; Seismology; High-performance computing
• ISSN: 0098-3004; 1873-7803
• DOI: 10.1016/j.cageo.2024.105644

To enable flexible and accurate seismic wave simulations at continental scales (10°−60°) based on the spectral-element method using the open-source SPECFEM3D_Cartesian package, we develop a toolkit, Cube2sph, that allows the generation of customized spherical meshes that account for the Earth’s curvature. This toolkit enables the usage of the perfectly matched layer (PML) absorbing boundary condition even when the artificial boundaries do not align with the coordinate axes. A series of numerical experiments are presented to validate the effectiveness of this toolkit. From these numerical experiments, we conclude that (1) continental-scale seismic wave simulations, especially surface wave simulations, can be more efficiently performed without the loss of accuracy by truncating the mesh at an appropriate depth, (2) curvilinear-grid PML can be used to effectively suppress artificial reflections for seismic wave simulations at continental scales, and (3) the Earth’s spherical geometry needs to be accurately meshed in order to obtain accurate simulation results for study regions larger than 8°. •Flexible mesh generation for accurate continental-scale seismic wave simulations.•Curvilinear-grid implementation of perfectly matched layer to accommodate deformed mesh.•Numerical examples with practical guidance for continental-scale seismic wave simulations.