Quick Search:    advanced search

GSW Home   GeoRef Home   My GSW Alerts   Contact GSW   About GSW   Journals List   Help

This Article Figures Only Full Text Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (27) Citing Articles via Google Scholar Google Scholar Articles by Day, S. M. Articles by Bradley, C. R. Search for Related Content GeoRef GeoRef Citation

Memory-Efficient Simulation of Anelastic Wave Propagation

Department of Geological Sciences
San Diego State University
San Diego, California, 92182
(S.M.D.)

Los Alamos National Laboratory
P.O. Box 1663, MS F665
Los Alamos, New Mexico, 87545
(C.R.B.)

Realistic anelastic attenuation can be incorporated rigorously into finite difference and other numerical wave propagation methods using internal or memory variables. The main impediment to the realistic treatment of anelastic attenuation in 3D is the very large computational storage requirement imposed by the additional variables. We previously proposed an alternative to the conventional memory-variable formulation, the method of coarse-grain memory variables, and demonstrated its effectiveness in acoustic problems. We generalize this memory-efficient formulation to 3D anelasticity and describe a fourth-order, staggered-grid finite-difference implementation. The anelastic coarse-grain method applied to plane wave propagation successfully simulates frequency-independent Q_p and Q_s. Apparent Q values are constant to within 4% tolerance over approximately two decades in frequency and biased less than 4% from specified target values. This performance is comparable to that achieved previously for acoustic-wave propagation, and accuracy could be further improved by optimizing the memory-variable relaxation times and weights. For a given assignment of relaxation times and weights, the coarse-grain method provides an eight-fold reduction in the storage requirement for memory variables, relative to the conventional approach. The method closely approximates the wavenumber-integration solution for the response of an anelastic half-space to a shallow dislocation source, accurately calculating all phases including the surface-diffracted SP phase and the Rayleigh wave. The half-space test demonstrates that the wave field-averaging concept underlying the coarse-grain method is effective near boundaries and in the presence of evanescent waves. We anticipate that this method will also be applicable to unstructured grid methods, such as the finite-element method and the spectral-element method, although additional numerical testing will be required to establish accuracy in the presence of grid irregularity. The method is not effective at wavelengths equal to and shorter than 4 grid cell dimensions, where it produces anomalous scattering effects. This limitation could be significant for very high-order numerical schemes under some circumstances (i.e., whenever wave-lengths as short as 4 grids are otherwise within the usable bandwidth of the scheme), but it is of no practical importance in our fourth-order finite-difference implementation.

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]

				B. T. Aagaard, T. M. Brocher, D. Dolenc, D. Dreger, R. W. Graves, S. Harmsen, S. Hartzell, S. Larsen, and M. L. Zoback Ground-Motion Modeling of the 1906 San Francisco Earthquake, Part I: Validation Using the 1989 Loma Prieta Earthquake Bulletin of the Seismological Society of America, April 1, 2008; 98(2): 989 - 1011. [Abstract] [Full Text] [PDF]

				R. W. Graves The Seismic Response of the San Bernardino Basin Region during the 2001 Big Bear Lake Earthquake Bulletin of the Seismological Society of America, February 1, 2008; 98(1): 241 - 252. [Abstract] [Full Text] [PDF]

				A. Askan, V. Akcelik, J. Bielak, and O. Ghattas Full Waveform Inversion for Seismic Velocity and Anelastic Losses in Heterogeneous Structures Bulletin of the Seismological Society of America, December 1, 2007; 97(6): 1990 - 2008. [Abstract] [Full Text] [PDF]

				S. Ma, R. J. Archuleta, and M. T. Page Effects of Large-Scale Surface Topography on Ground Motions, as Demonstrated by a Study of the San Gabriel Mountains, Los Angeles, California Bulletin of the Seismological Society of America, December 1, 2007; 97(6): 2066 - 2079. [Abstract] [Full Text] [PDF]

				K. Yamada and M. Horike Inference of Q-Values below 1 Hz from Borehole and Surface Data in the Osaka Basin by Three-Component Waveform Fitting Bulletin of the Seismological Society of America, August 1, 2007; 97(4): 1267 - 1278. [Abstract] [Full Text] [PDF]

				P. Liu, R. J. Archuleta, and S. H. Hartzell Prediction of Broadband Ground-Motion Time Histories: Hybrid Low/High- Frequency Method with Correlated Random Source Parameters Bulletin of the Seismological Society of America, December 1, 2006; 96(6): 2118 - 2130. [Abstract] [Full Text] [PDF]

				S. Ma and P. Liu Modeling of the Perfectly Matched Layer Absorbing Boundaries and Intrinsic Attenuation in Explicit Finite-Element Methods Bulletin of the Seismological Society of America, October 1, 2006; 96(5): 1779 - 1794. [Abstract] [Full Text] [PDF]

				P. Liu and R. J. Archuleta Efficient Modeling of Q for 3D Numerical Simulation of Wave Propagation Bulletin of the Seismological Society of America, August 1, 2006; 96(4A): 1352 - 1358. [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.

				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.

				Domain Reduction Method for Three-Dimensional Earthquake Modeling in Localized Regions, Part I: Theory Bulletin of the Seismological Society of America, April 1, 2003; 93(2): 817 - 824.

				Stability and Accuracy Analysis of Coarse-Grain Viscoelastic Simulations Bulletin of the Seismological Society of America, February 1, 2003; 93(1): 283 - 300.