Seismic wavefield calculation for laterally heterogeneous whole earth models using the pseudospectral method

SUMMARY A method for simulating seismic wave propagation in a laterally heterogeneous whole earth model is presented by solving the elastodynamic equations in 2-D cylindrical coordinates using the pseudospectral method (PSM). The PSM is an attractive timedomain technique that uses the fast Fourier transform for an accurate diierentiation of ¢eld variables in the equations. Since no dispersion error arises in Fourier diierentiation, even when using a large grid spacing, computer memoryand time are reduced by several orders of magnitude compared to traditional ¢nite-diierence methods. In order to examine body-wave phases with current computing resources, a slice through the sphere is approximated with a 2-D cylindrical model. An irregular grid spacing is used in the vertical coordinate to improve the treatment of the various structural boundaries appearing in the earth model by matching the heterogeneity in the model. Synthetic seismograms obtained by the PSM calculation are compared with those calculated from an exact simulation method for a spherically homogeneous (1-D) earth model and achieve good agreement. The PSM method is illustrated by constructing the seismic P^SV wave¢eld for strongly heterogeneous earth models including a shield structure near the free surface and velocity perturbations just above the core^mantle boundary. The visualization of the evolution of seismic P and S waves in time and space is tracked using a sequence of snapshots and synthetic seismograms. These displays allow direct insight into the nature of the complex seismic wave behaviour in the Earth’s interior.

[1]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[2]  Thorne Lay,et al.  The core-mantle boundary region , 1995 .

[3]  Robert J. Geller,et al.  DSM complete synthetic seismograms : SH, spherically symmetric, case , 1994 .

[4]  Robert J. Geller,et al.  Computation of synthetic seismograms and their partial derivatives for heterogeneous media with arbitrary natural boundary conditions using the Direct Solution Method , 1994 .

[5]  Emmanuel Chaljub,et al.  Sensitivity of SS precursors to topography on the upper‐mantle 660‐km discontinuity , 1997 .

[6]  P. Williamson,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part 1, 2.5-D modeling method , 1995 .

[7]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[9]  Robert J. Geller,et al.  Computation of complete synthetic seismograms for laterally heterogeneous models using the Direct Solution Method , 1997 .

[10]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[11]  George A. McMechan,et al.  Simulation of long‐period 3‐D elastic responses for whole earth models , 1995 .

[12]  José M. Carcione,et al.  An accurate and efficient scheme for wave propagation in linear viscoelastic media , 1990 .

[13]  Robert B. Herrmann,et al.  SH-wave generation by dislocation sources—A numerical study , 1979, Bulletin of the Seismological Society of America.

[14]  Bengt Fornberg,et al.  The pseudospectral method; accurate representation of interfaces in elastic wave calculations , 1988 .

[15]  D. Kosloff,et al.  Acoustic wave propagation in 2-D cylindrical coordinates , 1990 .

[16]  Heiner Igel,et al.  P‐SV wave propagation in the Earth's mantle using finite differences: Application to heterogeneous lowermost mantle structure , 1996 .

[17]  Hiroshi Takenaka,et al.  2.5-D modelling of elastic waves using the pseudospectral method , 1996 .

[18]  Alfred Behle,et al.  Elastic wave propagation simulation in the presence of surface topography , 1992 .

[19]  D. Kosloff,et al.  Solution of the equations of dynamic elasticity by a Chebychev spectral method , 1990 .

[20]  M. Wysession,et al.  Visualization of whole mantle propagation of seismic shear energy using normal mode summation , 1994 .

[21]  Comment on “A comparison of finite-difference and fourier method calculations of synthetic seismograms” , 1990, Bulletin of the Seismological Society of America.

[22]  David Kessler,et al.  Elastic wave propagation using cylindrical coordinates , 1991 .

[23]  Bertrand Meyer,et al.  Crustal thickening in Gansu‐Qinghai, lithospheric mantle subduction, and oblique, strike‐slip controlled growth of the Tibet plateau , 1998 .

[24]  B. Fornberg.  The pseudospectral method: Comparisons with finite differences for the elastic wave equation , 1987 .

[25]  Z. Alterman,et al.  Pulse Propagation in a Laterally Heterogeneous Solid Elastic Sphere , 1970 .

[26]  Y. Fung.  Foundations of solid mechanics , 1965 .

[27]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[28]  Moshe Reshef,et al.  Elastic wave calculations by the Fourier method , 1984 .

[29]  George A. McMechan,et al.  2-D pseudo-spectral viscoacoustic modeling in a distributed-memory multi-processor computer , 1993, Bulletin of the Seismological Society of America.

[30]  R. Hilst,et al.  Upper-mantle shear velocity beneath eastern Australia from inversion of waveforms from SKIPPY portable arrays , 1996 .

[31]  L. W. Braile,et al.  Reply to J. Vidale's “Comment on ‘A comparison of finite-difference and fourier method calculations of synthetic seismograms’” , 1989, Bulletin of the Seismological Society of America.

[32]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[33]  R. Clayton,et al.  Finite-difference seismograms for SH waves , 1985 .

[34]  Hiroshi Takenaka,et al.  Parallel 3-D pseudospectral simulation of seismic wave propagation , 1998 .

[35]  Wolfgang Friederich,et al.  COMPLETE SYNTHETIC SEISMOGRAMS FOR A SPHERICALLY SYMMETRIC EARTH BY A NUMERICAL COMPUTATION OF THE GREEN'S FUNCTION IN THE FREQUENCY DOMAIN , 1995 .

[36]  Heiner Igel,et al.  Frequency-dependent effects on travel times and waveforms of long-period S and SS waves , 1997 .

[37]  Moshe Reshef,et al.  Three-dimensional elastic modeling by the Fourier method , 1988 .

[38]  George A. McMechan,et al.  Two-dimensional elastic pseudo-spectral modeling of wide-aperture seismic array data with application to the Wichita Uplift-Anadarko Basin region of southwestern Oklahoma , 1990, Bulletin of the Seismological Society of America.

[39]  Heiner Igel,et al.  SH-wave propagation in the whole mantle using high-order finite differences , 1995 .