General purpose graphic processing unit(GPU) calculation technology is gradually widely used in various fields.Its mode of single instruction,multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps.In this study,we introduce a GPU-based parallel calculation method of a precise integration method(PIM) for seismic forward modeling.Compared with CPU single-core calculation,GPU parallel calculating perfectly keeps the features of PIM,which has small bandwidth,high accuracy and capability of modeling complex substructures,and GPU calculation brings high computational efficiency,which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.
Some recent publications presented a result suggesting that Zipingpu reservoir hastened the occurrence of the 2008 MS8.0 Wenchuan earthquake by tens to hundreds of years.Their researches calculated the Coulomb stress change induced by Zipingpu reservoir on the rupturing fault of Wenchuan earthquake.Their results,however,are critically dependent upon the 3-D event location,reservoir location,and the fault plane orientation.We repeated Ge et al.'s work in this paper and found that an improper dip angle parameter of their 2-D fault model might lead to a wrong conclusion.Both the modeling results based on the 2-D model and 3-D model with proper fault parameters will show Coulomb stress changes alone were neither large enough nor had the correct orientation to affect the occurrence of Wenchuan earthquake,which supports our recent argument based on the local seismicity analysis and the induced Coulomb stress change calculation with a 3-D model.
We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.
