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

This Article Full Text (PDF) References Citation Map Services Email this article to a friend Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by DAINTY, A. M. Articles by HARRIS, D. B. Search for Related Content GeoRef GeoRef Citation

Phase velocity estimation of diffusely scattered waves

LAWRENCE LIVERMORE NATIONAL LABORATORY, L-205, LIVERMORE, CALIFORNIA 94550
EARTH SCIENCES DIVISION/LWH AIR FORCE GEOPHYSICS LABORATORY, HANSCOM AFB MASSACHUSETTS 01731-5000

Abstract

Two methods are investigated for estimating the phase velocity of diffusely scattered seismic waves simultaneously arriving from different azimuths and recorded by a two-dimensional array of seismometers. The Hankel spectrum is the average of the frequency-wavenumber (FK) spectrum over all azimuths, while the wavenumber spectrum is derived by integrating the FK spectrum around a contour of constant phase velocity, i.e., a circle centered on the origin in the wavenumber plane. If the conventional estimate of the FK spectrum using the covariance matrix of the seismometer signals is integrated, a closed form for both the Hankel spectrum and the wavenumber spectrum may be found; the two spectra are very similar, the wavenumber spectrum being equal to the Hankel spectrum times the wavenumber. In spite of this similarity, however, we find that the two formulations have significantly different behavior for small wavenumbers, i.e., high phase velocities. In both cases there is a highest (true) velocity that can be estimated from the spectral maximum for a given array aperture ("velocity cut-off"). The Hankel spectrum estimates too high a velocity; for true wavenumbers below a certain limit, infinite velocity is estimated. The wavenumber spectrum, on the other hand, estimates too low a velocity, and there is an upper limit on the estimated velocity. An example illustrating these difficulties for the two methods is given for teleseismic P coda of an event recorded at the NORESS array in southern Norway: in spite of the problems, the analysis is able to demonstrate that the coda consists of two components; a coherent P-wave component with a high phase velocity and a diffuse S-wave component of low phase velocity. The cut-off and bias problem are investigated by numerical simulation for the NORESS array using azimuthal averaging and synthetic signals. The results confirm and quantify the cut-off problem at low wavenumbers and indicate that wavenumbers estimated from the Hankel and wavenumber spectra maxima bracket the true wavenumber, with the Hankel spectrum estimate being low (phase velocity too high) and the wavenumber spectrum estimate high. The bias of both methods decreases with increasing wavenumber (decreasing phase velocity) and they are both asymptotically unbiased. The wavenumber spectrum has a superior performance at low wavenumbers (high phase velocity), but the Hankel spectrum gives superior results at high wavenumbers (low phase velocity). The product of the linear wavenumber (= 1/wavelength) and the array aperture define "high" and "low" wavenumbers; for low wavenumbers, the product is 1 or less. In an Appendix, we find absolute lower bounds on the cut-offs analytically. The problems could be mitigated by using high resolution methods.

Footnotes

^* Present address: Earth Resources Laboratory, Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139.

This article has been cited by other articles:

				Propagation of Crustal-Waveguide-Trapped Pg and Seismic Velocity Structure in the South Indian Shield Bulletin of the Seismological Society of America, October 1, 2000; 90(5): 1281 - 1294.

				B. Chouet, G. De Luca, G. Milana, P. Dawson, M. Martini, and R. Scarpa Shallow velocity structure of Stromboli volcano, Italy, derived from small-aperture array measurements of Strombolian tremor Bulletin of the Seismological Society of America, June 1, 1998; 88(3): 653 - 666. [Abstract] [PDF]

				A. M. Dainty and C. A. Schultz Crustal reflections and the nature of regional P coda Bulletin of the Seismological Society of America, June 1, 1995; 85(3): 851 - 858. [Abstract] [PDF]

				C. A. LANGSTON and C. J. AMMON Scattering of teleseismic body waves along the hayward-calaveras fault system Bulletin of the Seismological Society of America, April 1, 1991; 81(2): 576 - 591. [Abstract] [PDF]

				S. C. BANNISTER, E. S. HUSEBYE, and B. O. RUUD Teleseismic P coda analyzed by three-component and array techniques: Deterministic location of topographic P-to-Rg scattering near the NORESS array Bulletin of the Seismological Society of America, December 1, 1990; 80(6B): 1969 - 1986. [Abstract] [PDF]

				A. M. DAINTY and M. NAFI TOKSOZ Array analysis of seismic scattering Bulletin of the Seismological Society of America, December 1, 1990; 80(6B): 2242 - 2260. [Abstract] [PDF]