Computers & Geosciences最新文献

{"title":"Fast forward modeling of grounded electrical-source transient electromagnetic based on inverse Laplace transform adaptive hybrid algorithm","authors":"Xiran You ,&nbsp;Jifeng Zhang ,&nbsp;Jiao Luo","doi":"10.1016/j.cageo.2024.105661","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105661","url":null,"abstract":"<div><p>Frequency–time conversion is a crucial step in grounded electrical-source transient electromagnetic response calculation, and the performance of the algorithm is directly related to the overall accuracy and speed of forward modeling. In mainstream algorithms, algorithms with high accuracy often have slow computation speed while algorithms with high efficiency have unsatisfactory accuracy, especially when facing inversion problems that are difficult to meet requirements. This paper introduces three inverse Laplace transform algorithms for this problem: the Gaver–Stehfest algorithm, the Euler algorithm, and the Talbot algorithm. The performance of each algorithm in forward modeling was analyzed using half-space and layered models, and the optimal selection schemes for algorithm weight coefficients were provided. The numerical calculation results show that the Gaver–Stehfest algorithm has a unique advantage in computational efficiency, while the Talbot algorithm and Euler algorithm meet the accuracy requirements. After considering both accuracy and efficiency, the Talbot algorithm is selected to replace conventional algorithms for calculation of grounded electrical-source transient electromagnetic forward modeling. In addition, this paper combines the characteristics of the Gaver–Stehfest algorithm and the Talbot algorithm to implement an adaptive hybrid algorithm. This algorithm uses the Gaver–Stehfest algorithm for forward modeling in the early times and the Talbot algorithm to compensate for the decrease in accuracy in the later times. Through the comparison of forward modeling calculations, it can be seen that the hybrid algorithm proposed in this paper fully utilizes the advantages of both algorithms. The hybrid algorithm greatly improves computational speed while meeting accuracy requirements, and has significant advantages over conventional algorithms.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"191 ","pages":"Article 105661"},"PeriodicalIF":4.2,"publicationDate":"2024-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141593893","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"DeLIA: A Dependability Library for Iterative Applications applied to parallel geophysical problems","authors":"Carla Santana ,&nbsp;Ramon C.F. Araújo ,&nbsp;Idalmis Milian Sardina ,&nbsp;Ítalo A.S. Assis ,&nbsp;Tiago Barros ,&nbsp;Calebe P. Bianchini ,&nbsp;Antonio D. de S. Oliveira ,&nbsp;João M. de Araújo ,&nbsp;Hervé Chauris ,&nbsp;Claude Tadonki ,&nbsp;Samuel Xavier-de-Souza","doi":"10.1016/j.cageo.2024.105662","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105662","url":null,"abstract":"<div><p>Many geophysical imaging applications, such as full-waveform inversion, often rely on high-performance computing to meet their demanding computational requirements. The failure of a subset of computer nodes during the execution of such applications can have a significant impact, as it may take several days or even weeks to recover the lost computation. To mitigate the consequences of these failures, it is crucial to employ effective fault tolerance techniques that do not introduce substantial overhead or hinder code optimization efforts. This paper addresses the primary research challenge of developing fault tolerance techniques with minimal impact on execution and optimization. To achieve this, we propose DeLIA, a Dependability Library for Iterative Applications designed for parallel programs that require data synchronization among all processes to maintain a globally consistent state after each iteration. DeLIA efficiently performs checkpointing and rollback of both the application’s global state and each process’s local state. Furthermore, DeLIA incorporates interruption detection mechanisms. One of the key advantages of DeLIA is its flexibility, allowing users to configure various parameters such as checkpointing frequency, selection of data to be saved, and the specific fault tolerance techniques to be applied. To validate the effectiveness of DeLIA, we applied it to a 3D full-waveform inversion code and conducted experiments to measure its overhead under different configurations using two workload schedulers. We also analyzed its behavior in preemptive circumstances. Our experiments revealed a maximum overhead of 8.8%, and DeLIA demonstrated its capability to detect termination signals and save the state of nodes in preemptive scenarios. Overall, the results of our study demonstrate the suitability of DeLIA to provide fault tolerance for iterative parallel applications.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"191 ","pages":"Article 105662"},"PeriodicalIF":4.2,"publicationDate":"2024-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0098300424001456/pdfft?md5=d3a34eb9baf8c143c8aae12bcda4ed57&pid=1-s2.0-S0098300424001456-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141485228","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Detection of the low-velocity layer using a convolutional neural network on passive surface-wave data: An application in Hangzhou, China","authors":"Xinhua Chen,&nbsp;Jianghai Xia,&nbsp;Jingyin Pang,&nbsp;Changjiang Zhou","doi":"10.1016/j.cageo.2024.105663","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105663","url":null,"abstract":"<div><p>Passive surface-wave methods using dense seismic arrays have gained growing attention in near-surface high-resolution imaging in urban environments. Deep learning (DL) in the extraction of dispersion curves and inversion can release a tremendous workload brought by dense seismic arrays. We presented a case study of imaging shear-wave velocity (Vs) structure and detecting low-velocity layer (LVL) in the Hangzhou urban area (eastern China). We used traffic-induced passive surface-wave data recorded by dense linear arrays. We extracted phase-velocity dispersion curves from noise recordings using seismic interferometry and multichannel analysis of surface waves. We adopted a convolutional neural network to estimate near-surface Vs models by inverting Rayleigh-wave fundamental-mode phase velocities. To improve the accuracy of the inversion, we utilized the sensitivities to weight the loss function. The average root mean square error from the weighted inversion is 46% lower than that from the unweighted DL inversion. The estimated pseudo-2D Vs profiles correspond to the velocities obtained from downhole seismic measurements. Compared with an investigation on the same survey area, our inversion results are more consistent with the Vs provided by downhole seismic measurements within 50–60 m where the LVL exists. The trained neural network successfully identified that the LVL is located at 50–60 m deep. To check the applicability of the trained neural network, we applied it to a nearby passive surface-wave survey and the inversion results agree with the existing investigation results. The two applications demonstrate the accuracy and efficiency of delineating near-surface Vs structures with the LVL from traffic-induced noise using the DL technique. The DL inversion has great potential for monitoring subsurface medium changes in urban areas.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"190 ","pages":"Article 105663"},"PeriodicalIF":4.2,"publicationDate":"2024-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141444683","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Sediment grain segmentation in thin-section images using dual-modal Vision Transformer","authors":"Dongyu Zheng ,&nbsp;Li Hou ,&nbsp;Xiumian Hu ,&nbsp;Mingcai Hou ,&nbsp;Kai Dong ,&nbsp;Sihai Hu ,&nbsp;Runlin Teng ,&nbsp;Chao Ma","doi":"10.1016/j.cageo.2024.105664","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105664","url":null,"abstract":"<div><p>Accurately identifying grain types in thin sections of sandy sediments or sandstones is crucial for understanding their provenance, depositional environments, and potential as natural resources. Although traditional computer vision methods and machine learning algorithms have been used for automatic grain identification, recent advancements in deep learning techniques have opened up new possibilities for achieving more reliable results with less manual labor. In this study, we present Trans-SedNet, a state-of-the-art dual-modal Vision-Transformer (ViT) model that uses both cross- (XPL) and plane-polarized light (PPL) images to achieve semantic segmentation of thin-section images. Our model classifies a total of ten grain types, including subtypes of quartz, feldspar, and lithic fragments, to emulate the manual identification process in sedimentary petrology. To optimize performance, we use SegFormer as the model backbone and add window- and mix-attention to the encoder to identify local information in the images and to best use XPL and PPL images. We also use a combination of focal and dice loss and a smoothing procedure to address imbalances and reduce over-segmentation. Our comparative analysis of several deep convolution neural networks and ViT models, including FCN, U-Net, DeepLabV3Plus, SegNeXT, and CMX, shows that Trans-SedNet outperforms the other models with a significant increase in evaluation metrics of mIoU and mPA. We also conduct an experiment to test the models' ability to handle dual-modal information, which reveals that the dual-modal models, including Trans-SedNet, achieve better results than single-modal models with the extra input of PPL images. Our study demonstrates the potential of ViT models in semantic segmentation of thin-section images and highlights the importance of dual-modal models for handling complex input in various geoscience disciplines. By improving data quality and quantity, our model has the potential to enhance the efficiency and reliability of grain identification in sedimentary petrology and relevant subjects.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"191 ","pages":"Article 105664"},"PeriodicalIF":4.2,"publicationDate":"2024-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141485227","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"FracAbut: A python toolbox for computing fracture stratigraphy using interface impedance","authors":"Paul Joseph Namongo Soro ,&nbsp;Juliette Lamarche ,&nbsp;Sophie Viseur ,&nbsp;Pascal Richard ,&nbsp;Fateh Messaadi","doi":"10.1016/j.cageo.2024.105656","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105656","url":null,"abstract":"<div><p>In Naturally Fractured Reservoirs (NFR) diffuse fractures arrangement results from mechanical stratigraphy and tectonic history during failure. Thus, modelling Discrete Fracture Network (DFN) requires to understand and to account for fracture relationships at bed-interface (abutment or crosscutting) in 3D through time (loading path). However, sampling fractures data meaningfully in subsurface has always been a challenge for geologist due to data scarcity.</p><p>To better understand and forecast fracture networks in stratified rocks, we study outcrops with a focus on geometric relationships between stratigraphic interfaces and fractures. This paper presents an original python toolbox called FracAbut. It is composed of 1 main and 2 auxiliary codes that quantify the geometric relation between fractures and stratigraphic interfaces from 1D (wells, scan-line) and 2D (digital image, photographs data). We calculate the Interface Impedance (<em>II</em>) that accounts for fracture abutment (crossing or not), persistence (single- or multi-bed) and propagation polarity (upward or downward). For each stratigraphic interface FracAbut provides information on fractures (type, number) and interface sensitivity (coupling strength).</p><p>First, we apply FracAbut on synthetic case studies, then, on naturally fractured and stratified carbonates in Berat, Albania. Using both 1D scan-line and 2D outcrop photograph, we show that i) a mechanical interface can have different coupling above and below based on propagation polarity, ii) FracAbut results can give useful insight on fracture transmissivity, iii) FracAbut is fast and efficient to quantify fracture patterns and classify mechanical interface impact; iv) they are no relation between bed thickness and fracture propagation.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"190 ","pages":"Article 105656"},"PeriodicalIF":4.4,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141424358","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A deep autoencoder network connected to geographical random forest for spatially aware geochemical anomaly detection","authors":"Zeinab Soltani ,&nbsp;Hossein Hassani ,&nbsp;Saeid Esmaeiloghli","doi":"10.1016/j.cageo.2024.105657","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105657","url":null,"abstract":"<div><p>Machine learning (ML) and deep learning (DL) techniques have recently shown encouraging performance in recognizing metal-vectoring geochemical anomalies within complex Earth systems. However, the generalization of these techniques to detect subtle anomalies may be precluded due to overlooking non-stationary spatial structures and intra-pattern local dependencies contained in geochemical exploration data. Motivated by this, we conceptualize in this paper an innovative algorithm connecting a DL architecture to a spatial ML processor to account for local neighborhood information and spatial non-stationarities in support of spatially aware anomaly detection. A deep autoencoder network (DAN) is trained to abstract deep feature codings (DFCs) of multi-element input data. The encoded DFCs represent the typical performance of a nonlinear Earth system, i.e., multi-element signatures of geochemical background populations developed by different geo-processes. A local version of the random forest algorithm, geographical random forest (GRF), is then connected to the input and code layers of the DAN processor to establish nonlinear and spatially aware regressions between original geochemical signals (dependent variables) and DFCs (independent variables). After contributions of the latter on the former are determined, residuals of GRF regressions are quantified and interpreted as spatially aware anomaly scores related to mineralization. The proposed algorithm (i.e., DAN‒GRF) is implemented in the R language environment and examined in a case study with stream sediment geochemical data pertaining to the Takht-e-Soleyman district, Iran. The high-scored anomalies mapped by DAN‒GRF, compared to those by the stand-alone DAN technique, indicated a stronger spatial correlation with locations of known metal occurrences, which was statistically confirmed by success-rate curves, Student's <span><math><mrow><mi>t</mi></mrow></math></span>‒statistic method, and prediction-area plots. The findings suggested that the proposed algorithm has an enhanced capability to recognize subtle multi-element geochemical anomalies and extract reliable insights into metal exploration targeting.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"190 ","pages":"Article 105657"},"PeriodicalIF":4.2,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141433891","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A dual watermarking algorithm for trajectory data based on robust watermarking and fragile watermarking","authors":"Yuchen Hu ,&nbsp;Xingxiang Jiang ,&nbsp;Changqing Zhu ,&nbsp;Na Ren ,&nbsp;Shuitao Guo ,&nbsp;Jia Duan ,&nbsp;Luanyun Hu","doi":"10.1016/j.cageo.2024.105655","DOIUrl":"10.1016/j.cageo.2024.105655","url":null,"abstract":"<div><p>Digital watermarking technology plays a crucial role in securing trajectory data. However, as trajectory data usage scenarios continue to expand, the security requirements for it have changed from a single copyright protection to one that takes into account data integrity. Existing digital watermarking algorithms for trajectory data can only choose between implementing copyright protection or ensuring integrity, unable to simultaneously achieve both functionalities. This limitation impedes the sharing and utilization of trajectory data. A dual watermarking algorithm that combines robust and fragile watermarking was innovatively proposed to solve this problem based on the geometric domain. Firstly, a set of feature points is extracted from the trajectory, and the farthest point pair of the minimum convex hull of the feature points is set as fixed points. The robust watermark is then embedded in the angles constructed by the feature points and the fixed points using quantization index modulation. Meanwhile, the trajectory points are grouped based on the angle and distance ratio constructed from the trajectory points to the fixed points. In each group, the spatiotemporal attributes of the trajectory points are mapped to the fragile watermark, which is then embedded into the distance ratios constructed by the trajectory points. Experimental results show that the proposed algorithm achieves both copyright protection and integrity verification for trajectory data and exhibits stronger robustness and tampering localization ability. This research can provide security and privacy protection for trajectory data and contribute positively to the application of trajectory data.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"191 ","pages":"Article 105655"},"PeriodicalIF":4.2,"publicationDate":"2024-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141411656","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"GPUPRSI: GPU implementation of seismic interferometry for retrieving reflection responses from passive source seismic recordings","authors":"Jun Zheng,&nbsp;Guofeng Liu","doi":"10.1016/j.cageo.2024.105654","DOIUrl":"https://doi.org/10.1016/j.cageo.2024.105654","url":null,"abstract":"<div><p>Passive seismic exploration is an environmentally friendly, economical, and highly accessible exploration method that is widely used in different-scale subsurface imaging. Retrieving reflections from ambient noise is a newly developed passive seismic method that can be used in many fields such as mineral exploration and near-surface imaging, but the interferometry calculation is time-consuming because it requires a longer period of data acquisition to improve the signal-to-noise ratio of the retrieved reflections. In this study, we introduced a graphical processing unit (GPU)- based implementation of seismic interferometry for retrieving reflection responses from passive source seismic recordings. Because all traces are involved in the computation process, but the size of passive source data often exceeds available memory, repeated disk reads lead to a decrease in computational efficiency. We design a strategy of grouping computations followed by stacking to minimize disk input and output, simultaneously keeping the memory requirements low. Passive source data is read and written only once, and there is no requirement for the memory size to be greater than the data size. Additionally, acceleration technologies such as asynchronous execution, asynchronous memory transfer, and GPU-accelerated libraries are used. We test the efficiency using short data where the number of sampling points per trace is on the order of 30,000, and long data, where the number of sampling points per trace is on the order of 30,000,000. For short data, the average speedup is 4 compared with a multi-core central processing unit (CPU); for long data, the speedup can reach 24. The GPU-based implementation of interferometry greatly reduces the calculation time for retrieving reflections from passive source seismic recordings, providing a solution to the problem of large calculation in three-dimensional (3D) passive source reflection exploration.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"190 ","pages":"Article 105654"},"PeriodicalIF":4.4,"publicationDate":"2024-06-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141325281","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"VSP wavefields separating method based on parallel local slope scanning","authors":"Wu Li,&nbsp;Yuyong Yang,&nbsp;Bocheng Tao,&nbsp;Zhengyang Wang,&nbsp;Huailai Zhou,&nbsp;Yuanjun Wang","doi":"10.1016/j.cageo.2024.105643","DOIUrl":"10.1016/j.cageo.2024.105643","url":null,"abstract":"<div><p>Due to the detectors being distributed in the target medium, vertical seismic profile (VSP) is a seismic observation method that has the advantages of high resolution and a high signal-to-noise ratio. Separating the mixed wavefields into upgoing and downgoing waves can obtain more obvious dynamic and kinematic characteristics of seismic waves, guiding subsequent imaging and interpretation. The traditional method is mainly based on the pickup of first breaks. Flattening the global seismic data by first breaks can enhance the seismic events through techniques like median filter and singular value decomposition (SVD). However, this method relies on high-precision pickup of first breaks and is limited by zero offset. To address this limitation, we introduce an improved median filtering separation method. This method employs separations of the local dip angles into positive and negative through multi-window scanning (MWS). Due to the high accuracy and robustness of this method in 2-D VSP data, we propose to use the local dip angle of the wavefields to median filter the wavefield through positive and negative angles to obtain upgoing and downgoing waves. This method for wavefield separation is optimized by iteratively identifying the directions of the seismic data. However, this optimization comes at the cost of increased computational requirements, especially under high-precision conditions. To help alleviate this problem, we use multi-thread parallel computing technique on a multi-core central processing unit (CPU) to improve computational efficiency. Finally, we validate the proposed method by testing it on synthetic seismic data and field VSP data, respectively. The results show that this wavefield separation method has advantages in terms of accuracy and robustness compared to the median filtering method based on the first break picking.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"191 ","pages":"Article 105643"},"PeriodicalIF":4.2,"publicationDate":"2024-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141416258","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Cube2sph : A toolkit enabling flexible and accurate continental-scale seismic wave simulations using the SPECFEM3D_Cartesian package","authors":"Tianshi Liu ,&nbsp;Kai Wang ,&nbsp;Yujiang Xie ,&nbsp;Bin He ,&nbsp;Ting Lei ,&nbsp;Nanqiao Du ,&nbsp;Ping Tong ,&nbsp;Yingjie Yang ,&nbsp;Catherine A. Rychert ,&nbsp;Nicholas Harmon ,&nbsp;Giovanni Grasselli ,&nbsp;Qinya Liu","doi":"10.1016/j.cageo.2024.105644","DOIUrl":"10.1016/j.cageo.2024.105644","url":null,"abstract":"<div><p>To enable flexible and accurate seismic wave simulations at continental scales (<span><math><mrow><mn>10</mn><mo>°</mo><mo>−</mo><mn>60</mn><mo>°</mo></mrow></math></span>) based on the spectral-element method using the open-source <span>SPECFEM3D_Cartesian</span> package, we develop a toolkit, <span>Cube2sph</span>, 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°.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"190 ","pages":"Article 105644"},"PeriodicalIF":4.4,"publicationDate":"2024-06-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141395207","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}