Parallel PSM/FDM Hybrid Simulation of Ground Motions from the 1999 Chi-Chi, Taiwan, Earthquake

Abstract— A new technique for the parallel computing of 3-D seismic wave propagation simulation is developed by hybridizing the Fourier pseudospectral method (PSM) and the finite-difference method (FDM). This PSM/FDM hybrid offers a good speed-up rate using a large number of processors. To show the feasibility of the hybrid, a numerical 3-D simulation of strong ground motion was conducted for the 1999 Chi-Chi, Taiwan earthquake (Mw 7.6). Comparisons between the simulation results and observed waveforms from a dense strong ground motion network in Taiwan clearly demonstrate that the variation of the subsurface structure and the complex fault slip distribution greatly affect the damage during the Chi-Chi earthquake. The directivity effect of the fault rupture produced large S-wave pulses along the direction of the rupture propagation. Slips in the shallow part of the fault generate significant surface waves in Coastal Plain along the western coast. A large velocity gradient in the upper crust can propagate seismic waves to longer distances with minimum attenuation. The S waves and surface waves were finally amplified further by the site effect of low-velocity sediments in basins, and caused the significant disasters.

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

[2]  W. M. Kohler,et al.  Reflected seismic waves and their effect on strong shaking during the 1989 Loma Prieta, California, earthquake , 1996, Bulletin of the Seismological Society of America.

[3]  K. Ma,et al.  Rupture Process of the 1999 Chi-Chi, Taiwan, Earthquake from the Inversion of Teleseismic Data , 2000 .

[4]  Jeen‐Hwa Wang,et al.  Three-Dimensional Seismic Velocity Structure of the Crust and Uppermost Mantle beneath Taiwan , 1996 .

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

[6]  Y. Tsai,et al.  A preliminary report on the 1999 Chi-Chi (Taiwan) earthquake , 2000 .

[7]  Tzay-Chyn Shin,et al.  A Free-Field Strong Motion Network in Taiwan:TSMIP , 1999 .

[8]  Shu-Huei Hung,et al.  Modelling anisotropic wave propagation in oceanic inhomogeneous structures using the parallel multidomain pseudo-spectral method , 1998 .

[9]  Yeong Tein Yeh,et al.  Shallow crustal structure from short-period Rayleigh-wave dispersion data in southwestern Taiwan , 1997, Bulletin of the Seismological Society of America.

[10]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[11]  R. Rau,et al.  Tomographic imaging of lithospheric structures under Taiwan , 1995 .

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

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

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

[15]  Robert B. Herrmann,et al.  The effect of crustal structure on strong ground motion attenuation relations in eastern North America , 1987 .

[16]  T. Furumura,et al.  Parallel 3-D Simulation of Ground Motion for the 1995 Kobe Earthquake: The Component Decomposition Approach , 2000 .

[17]  P. Somerville,et al.  The influence of critical Moho Reflections on strong ground motions recorded in San Francisco and Oakland during the 1989 Loma Prieta Earthquake , 1990 .

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

[19]  Cheng‐Horng Lin,et al.  A study of crustal structures beneath northeastern Taiwan: Possible evidence of the western extension of Okinawa Trough , 1989 .

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