3D Fourth-Order Staggered-Grid Finite-Difference Schemes: Stability and Grid Dispersion

We investigated stability and grid dispersion in the 3D fourth-order in space, second-order in time, displacement-stress staggered-grid finite-difference scheme. Though only displacement-stress scheme is explicitly treated, all results also apply to the velocity-stress and displacement-velocity-stress finite-difference schemes. We derived independent stability conditions for the P and S waves by exact separation of equations for the two types of waves. Since the S -wave group velocity can differ from the actual velocity as much as 5% for the sampling ratio 1/5 (that is usually used in modeling), we recommend to sample a minimum S wavelength by six grid spacings. Grid dispersion is strongest for a wave propagating in the direction of a coordinate axis and weakest for a wave propagating along a body diagonal. Grid dispersion in the fourth-order scheme for the sampling ratios s = 1/5 and s = 1/6 is smaller than grid dispersion in the second-order scheme for s = 1/10 and s = 1/12, respectively.

[1]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[2]  Jean Virieux,et al.  SH-wave propagation in heterogeneous media; velocity-stress finite-difference method , 1984 .

[3]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[4]  Bent O. Ruud,et al.  3-D finite-difference elastic wave modeling including surface topography , 1998 .

[5]  Hiroyuki Fujiwara,et al.  3D finite-difference method using discontinuous grids , 1999, Bulletin of the Seismological Society of America.

[6]  Peter Moczo,et al.  EFFICIENCY AND OPTIMIZATION OF THE 3-D FINITE-DIFFERENCE MODELING OF SEISMIC GROUND MOTION , 2001 .

[7]  Robert W. Graves,et al.  Ground-motion amplification in the Santa Monica area: Effects of shallow basin-edge structure , 1998, Bulletin of the Seismological Society of America.

[8]  Arben Pitarka,et al.  3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing , 1999 .

[9]  Kojiro Irikura,et al.  Three-dimensional simulation of the near-fault ground motion for the 1995 Hyogo-Ken Nanbu (Kobe), Japan, earthquake , 1998, Bulletin of the Seismological Society of America.

[10]  Gerard T. Schuster,et al.  Parsimonious staggered grid finite‐differencing of the wave equation , 1990 .

[11]  J. Neumann.  Mathematische grundlagen der Quantenmechanik , 1935 .

[12]  John Etgen,et al.  3-D wave propagation in the Los Angeles basin for the Whittier-Narrows earthquake , 1993 .

[13]  R. Graves.  Modeling three-dimensional site response effects in the Marina District Basin, San Francisco, California , 1993, Bulletin of the Seismological Society of America.

[14]  Bernard A. Chouet,et al.  A free-surface boundary condition for including 3D topography in the finite-difference method , 1997, Bulletin of the Seismological Society of America.

[15]  Kojiro Irikura,et al.  Modelling of ground motion in the Higashinada (Kobe) area for an aftershock of the 1995 January 17 Hyogo-ken Nanbu, Japan, earthquake , 1997 .

[16]  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.

[17]  Robert W. Graves,et al.  The seismic response of the Los Angeles basin, California , 1998, Bulletin of the Seismological Society of America.

[18]  Kim B. Olsen,et al.  Three-Dimensional Simulation of a Magnitude 7.75 Earthquake on the San Andreas Fault , 1995, Science.

[19]  Mark Noble,et al.  Nonlinear elastic waveform inversion of land seismic reflection data , 1992 .