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.