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 Google Scholar Google Scholar Articles by YAO, Z. X. Articles by HARKRIDER, D. G. Search for Related Content

A generalized reflection-transmission coefficient matrix and discrete wavenumber method for synthetic seismograms

INSTITUTE OF GEOPHYSICS ACADEMIA SINICA, BEIJING, People's Republic of China
SEISMOLOGICAL LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 91125

Abstract

Expressions for displacements on the surface of a layered half-space due to point force are given in terms of generalized reflection and transmission coefficient matrices (Kennett, 1980) and the discrete wavenumber summation method (Bouchon, 1981). The Bouchon method with complex frequencies yields accurate near-field dynamic and static solutions.

The algorithm is extended to include simultaneous evaluation of multiple sources at different depths. This feature is the same as in Olson's finite element discrete Fourier Bessel code (DWFE) (Olson, 1982).

As numerical examples, we calculate some layered half-space problems. The results agree with synthetics generated with the Cagniard-de Hoop technique, P-SV modes, and DWFE codes. For a 10-layered crust upper mantle model with a bandwidth of 0 to 10 Hz, this technique requires one-tenth the time of the DWFE calculation. In the presence of velocity gradients, where finer layering is required, the DWFE code is more efficient.