Normal modes of a laterally heterogeneous body: A one-dimensional example

Various methods, including first- and second-order perturbation theory and variational methods have been proposed for calculating the normal modes of a laterally heterogeneous earth. In this paper, we test all three of these methods for a simple one-dimensional example for which the exact solution is available: an initially homogeneous "string" in which the density and stiffness are increased in one half and decreased in the other by equal amounts. It is found that first-order perturbation theory (as commonly applied in seismology) yields only the eigenvalues and eigenfunctions for a string with the average elastic properties; second-order perturbation theory is worse, because the eigenfunction is assumed to be the original eigenfunction plus small correction terms, but actually may be almost completely different. The variational method (Rayleigh-Ritz), using the unperturbed modes as trial functions, succeeds in giving correct eigenvalues and eigenfunctions even for modes of high-order number. For the example problem only the variational solution correctly yields the transient solution for excitation by a point force, including correct amplitudes for waves reflected by and transmitted through the discontinuity. Our result suggests but does not demonstrate, that the variational method may be the most appropriate method for finding the normal modes of a laterally heterogeneous earth model, particularly if the transient solution is desired.