Abstract
A new method to obtain theoretical seismograms for a layered medium is presented, where an appropriate path of integration for the integral with respect to wave-number is chosen. In this method, upper limit of the integral can be taken as a finite value depending on the time interval in the time domain for which seismograms are computed and pole contributions are implicitly included without consideration of residues at the poles.
Comparison of the present method with those in previous investigations is made and problems of the generation of seismic surface waves are numerically evaluated.