The theory of finite frequency body wave synthetic seismograms in inhomogeneous elastic media

b
A new method is presented by means of which one can compute finite frequency synthetic seismograms for media whose velocity and density are continuous functions of two or three spatial variables. Basically, the method is a generalization of the familiar phase integral method, to which it reduces in a stratified medium. For a given source location the travel-time and distance functions needed to compute synthetics are obtained by numerically tracing rays through the model. This information is then used to evaluate a double integral over frequency and take-off angle at the source. The solution obtained reduces to the geometrical optics solution wherever that is valid but it also works in shadows and at caustics without knowing explicitly where these may be located. The method can be used as a spectral method, in which the integral over take-off angle is evaluated first, or as a slowness method, in which the frequency integral is evaluated first.