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 (16) Citing Articles via Google Scholar Google Scholar Articles by Gottschämmer, E. Articles by Olsen, K. B. Search for Related Content GeoRef GeoRef Citation

Accuracy of the Explicit Planar Free-Surface Boundary Condition Implemented in a Fourth-Order Staggered-Grid Velocity-Stress Finite-Difference Scheme

Geophysikalisches Institut
Karlsruhe Universität
76187 Karlsruhe, Germany
(E.G.)

Institute for Crustal Studies
UC Santa Barbara
Santa Barbara, California 93106-1100
(K.B.O.)

We compute the accuracy of two implementations of the explicit planar free-surface boundary condition for 3D fourth-order velocity-stress staggered-grid finite differences, 1/2 grid apart vertically, in a uniform half-space. Due to the staggered grid, the closest distance between the free surface and some wave-field components for both implementations is 1/2-grid spacing. Overall, the differences in accuracy of the two implementations are small. When compared to a reflectivity solution computed at the staggered positions closest to the surface, the total misfit for all three components of the wave field is generally found to be larger for the free surface colocated with the normal stresses, compared to that for the free surface colocated with the xz and yz stresses. However, this trend is reversed when compared to the reflectivity solution exactly at the free surface (the misfit encountered in staggered-grid modeling). When the wave field is averaged across the free surface, thereby centering the staggered wave field exactly on the free surface, the free-surface condition colocated with the xz and yz stresses generates the smallest total misfit for increasing epicentral distance. For an epicentral distance/hypocentral depth of 10, the total misfit of this condition is about 15% smaller than that for the condition colocated with the normal stresses, mainly controlled by the misfit on the Rayleigh wave.

This article has been cited by other articles:

				K. B. Olsen, S. M. Day, J. B. Minster, Y. Cui, A. Chourasia, D. Okaya, P. Maechling, and T. Jordan TeraShake2: Spontaneous Rupture Simulations of Mw 7.7 Earthquakes on the Southern San Andreas Fault Bulletin of the Seismological Society of America, June 1, 2008; 98(3): 1162 - 1185. [Abstract] [Full Text] [PDF]

				S. M. Corish, C. R. Bradley, and K. B. Olsen Assessment of a Nonlinear Dynamic Rupture Inversion Technique Applied to a Synthetic Earthquake Bulletin of the Seismological Society of America, June 1, 2007; 97(3): 901 - 914. [Abstract] [Full Text] [PDF]

				L. Zhao, T. H. Jordan, K. B. Olsen, and P. Chen Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models Bulletin of the Seismological Society of America, December 1, 2005; 95(6): 2066 - 2080. [Abstract] [Full Text] [PDF]

				An Efficient Finite-Difference Method for Simulating 3D Seismic Response of Localized Basin Structures Bulletin of the Seismological Society of America, October 1, 2004; 94(5): 1690 - 1705.

				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]

				D.-J. Min, C. Shin, and H. S. Yoo Free-Surface Boundary Condition in Finite-Difference Elastic Wave Modeling Bulletin of the Seismological Society of America, February 1, 2004; 94(1): 237 - 250. [Abstract] [Full Text] [PDF]