INTERPRETATION OF SEISMIC SURFACE-WAVES IN REGIONAL NETWORKS - JOINT ESTIMATION OF WAVE-FIELD GEOMETRY AND LOCAL PHASE-VELOCITY - METHOD AND NUMERICAL TESTS

Citation
W. Friederich et E. Wielandt, INTERPRETATION OF SEISMIC SURFACE-WAVES IN REGIONAL NETWORKS - JOINT ESTIMATION OF WAVE-FIELD GEOMETRY AND LOCAL PHASE-VELOCITY - METHOD AND NUMERICAL TESTS, Geophysical journal international, 120(3), 1995, pp. 731-744
Citations number
15
Categorie Soggetti
Geosciences, Interdisciplinary
ISSN journal
0956540X
Volume
120
Issue
3
Year of publication
1995
Pages
731 - 744
Database
ISI
SICI code
0956-540X(1995)120:3<731:IOSSIR>2.0.ZU;2-U
Abstract
Teleseismic surface waves in general strongly deviate from plane waves as is evident from the analysis of surface-wave data recorded with de nse networks. This causes conventional, ray-based tomographic techniqu es to break down if applied to network surface-wave data. We present a new inversion method based on the acoustic-wave equation and applicab le to vertical-component surface-wave data which successfully deals wi th non-plane wavefield geometries. The basic idea of the method is a j oint estimation of the incoming wavefield and heterogeneous structure within the network region. Crucial to the success of the method is an adequate parametrization of the incoming wavefield which is realized u sing Hermite-Gaussian basis functions. Additionally, we apply a constr aint on the wavefield parameters that expresses the idea that the samp les of the wavefield taken at the station locations should be represen tative for the wavefield in the whole network region. In this way, wav efields that show stronger fluctuations in spectral amplitude than obs erved at the stations are rejected. To represent heterogeneous structu re within the network region we use an expansion into 2-D Hermite-Gaus sian functions. Provided that the density of stations is sufficient, t he proposed method retrieves heterogeneous structure in the network re gion very well. It is not sensitive to noise or non-uniform azimuthal coverage of earthquakes. Moreover, it yields smoothed versions of the true model if the roughness of the latter has been underestimated in t he inversion. Conversely, if the true model is much smoother than anti cipated, inspection of the trade-off between model smoothness and data misfit allows us to find the correct model. The limiting factor for t he inversion is the density of stations, which must allow for a reliab le interpolation of the observed wavefield within the network. Therefo re, in order to perform regional surface-wave studies it is essential to deploy seismic stations in the region of interest itself.