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 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 Web of Science (13) Citing Articles via Google Scholar Google Scholar Articles by Saito, T. Articles by Ohtake, M. Search for Related Content GeoRef GeoRef Citation

Simulating the Envelope of Scalar Waves in 2D Random Media Having Power-Law Spectra of Velocity Fluctuation

Department of Geophysics
Graduate School of Science
Tohoku University
Aramaki-aza Aoba, Aoba-ku
Sendai 980-8578, Japan
saito{at}zisin.geophys.tohoku.ac.jp
(T.S., H.S., M.O.)
Los Alamos Seismic Research Center
Los Alamos National Laboratory
Los Alamos, New Mexico 87545
(M.F.)

Manuscript received 10 April 2002.

During propagation through random media, impulsive waves radiated from a point source decrease in amplitude and increase in duration with increasing travel distance. The excitation of coda waves is prominent at long lapse time. We use a finite-difference method to numerically simulate scalar waves that propagate through random media characterized by a von Kármán autocorrelation function. The power spectral density function of fractional velocity fluctuation for -th order von Kármán-type random media obeys a power law at large wavenumbers. Such media are considered to be appropriate models for the random component of the structure of the Earth's lithosphere. The average of the square of numerically calculated wave traces over an ensemble of random media gives the reference envelope for the evaluation of envelope simulation methods. The Markov approximation method gives reliable quantitative predictions of the entire envelope for random media that are poor in short wavelength components of heterogeneity ( = 1.0), while it fails to predict the coda envelope for random media that have rich short-wavelength components ( = 0.1). The radiative-transfer theory reliably predicts the coda excitation for = 0.1 when the momentum-transfer scattering coefficient is used as the effective isotropic scattering coefficient. Replacing the direct term of the radiative-transfer solution with the envelope of the Markov approximation, we propose a new method for simulating the entire envelope from the direct arrival through the coda. The method quantitatively explains the whole envelope for = 0.1. For the case of = 0.5, however, our method predicts too much coda excitation. In such a case, the method can explain whole envelopes by using the effective scattering coefficient estimated from coda excitation.