HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Full Text (PDF) Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (34) Citing Articles via Google Scholar Google Scholar Articles by Moczo, P. Articles by Halada, L. Search for Related Content GeoRef GeoRef Citation

3D Heterogeneous Staggered-Grid Finite-Difference Modeling of Seismic Motion with Volume Harmonic and Arithmetic Averaging of Elastic Moduli and Densities

Peter Moczo, Jozef Kristek, Václav Vavryuk, Ralph J. Archuleta and Ladislav Halada

Department of Geophysics
Faculty of Mathematics, Physics and Informatics
Comenius University
Mlynská dolina F1
842 48 Bratislava
Slovak Republic
(P.M.)
Geophysical Institute
Slovak Academy of Sciences
Dúbravská cesta 9
845 28 Bratislava
Slovak Republic
(P.M., J.K., L.H.)
Geophysical Institute
Academy of Sciences of the Czech Republic
Boní II
141 31 Praha
Czech Republic
(V.V.)
Institute for Crustal Studies
University of California at Santa Barbara
Santa Barbara, California 93106
(R.J.A.)

Manuscript received 21 May 2001.

We analyze the problem of a heterogeneous formulation of the equation of motion and propose a new 3D fourth-order staggered-grid finite-difference (FD) scheme for modeling seismic motion and seismic-wave propagation.

We first consider a 1D problem for a welded planar interface of two half-spaces. A simple physical model of the contact of two media and mathematical considerations are shown to give an averaged medium representing the contact of two media. An exact heterogeneous formulation of the equation of motion is a basis for constructing the corresponding heterogeneous FD scheme.

In a much more complicated 3D problem we analyze a planar-interface contact of two isotropic media (both with interface parallel to a coordinate plane and interface in general position in the Cartesian coordinate system) and a nonplanar-interface contact of two isotropic media. Because in the latter case 21 elastic coefficients at each point are necessary to describe the averaged medium, we consider simplified boundary conditions for which the averaged medium can be described by only two elastic coefficients.

Based on the simplified approach we construct the explicit heterogeneous 3D fourth-order displacement-stress FD scheme on a staggered grid with the volume harmonic averaging of the shear modulus in grid positions of the stress-tensor components, volume harmonic averaging of the bulk modulus in grid positions of the normal stress-tensor components, and volume arithmetic averaging of density in grid positions of the displacement components.

Our displacement-stress FD scheme can be easily modified into the velocity-stress or displacement-velocity-stress FD schemes.

The scheme allows for an arbitrary position of the material discontinuity in the spatial grid. Numerical tests for 12 configurations in four types of models show that our scheme is more accurate than the staggered-grid schemes used so far.

Numerical examples also show that differences in thickness of a soft surface or interior layer smaller than one grid spacing can cause considerable changes in seismic motion. The results thus underline the importance of having a FD scheme with sufficient sensitivity to heterogeneity of the medium.

This article has been cited by other articles:

				T. Ichimura, M. Hori, and H. Kuwamoto Earthquake Motion Simulation with Multiscale Finite-Element Analysis on Hybrid Grid Bulletin of the Seismological Society of America, August 1, 2007; 97(4): 1133 - 1143. [Abstract] [Full Text] [PDF]

				J. Kristek and P. Moczo On the Accuracy of the Finite-Difference Schemes: The 1D Elastic Problem Bulletin of the Seismological Society of America, December 1, 2006; 96(6): 2398 - 2414. [Abstract] [Full Text] [PDF]

				L. A. Dalguer and S. M. Day Comparison of Fault Representation Methods in Finite Difference Simulations of Dynamic Rupture Bulletin of the Seismological Society of America, October 1, 2006; 96(5): 1764 - 1778. [Abstract] [Full Text] [PDF]

				P. Dangla, J.-F. Semblat, H. Xiao, and N. Delepine A Simple and Efficient Regularization Method for 3D BEM: Application to Frequency-Domain Elastodynamics Bulletin of the Seismological Society of America, October 1, 2005; 95(5): 1916 - 1927. [Abstract] [Full Text] [PDF]

				Hybrid Modeling of Elastic P-SV Wave Motion: A Combined Finite-Element and Staggered-Grid Finite-Difference Approach Bulletin of the Seismological Society of America, August 1, 2004; 94(4): 1557 - 1563.

				P. Moczo, J. Kristek, and M. Galis Simulation of the Planar Free Surface with Near-Surface Lateral Discontinuities in the Finite-Difference Modeling of Seismic Motion Bulletin of the Seismological Society of America, April 1, 2004; 94(2): 760 - 768. [Abstract] [Full Text] [PDF]

				Seismic-Wave Propagation in Viscoelastic Media with Material Discontinuities: A 3D Fourth-Order Staggered-Grid Finite-Difference Modeling Bulletin of the Seismological Society of America, October 1, 2003; 93(5): 2273 - 2280.

				Reply to "Comment on `3D Site Effects: A Thorough Analysis of a High-Quality Dataset,' by F. J. Chavez-Garcia, J. Castillo, and W. R. Stephenson," by R. Paolucci and E. Faccioli Bulletin of the Seismological Society of America, October 1, 2003; 93(5): 2306 - 2316.

				Estimation of Q for Long-Period (>2 sec) Waves in the Los Angeles Basin Bulletin of the Seismological Society of America, April 1, 2003; 93(2): 627 - 638.