Seismogram synthesis using normal mode superposition: the locked mode approximation

Summary. A normal mode superposition approach is used to synthesize complete seismic codas for flat layered earth models and the P-SV phases. Only modes which have real eigenwavenumbers are used so that the search for eigenvalues in the complex wavenumber plane is confined to the real axis. In order to synthesize early P -wave arrivals by summing a number of'trapped'modes, an anomalously high velocity cap layer is added to the bottom of the structure so that most of the seismic energy is contained in the upper layers as high-order surface waves. Causality arguments are used to define time windows for which the resulting synthetic seismograms are close approximations to the exact solutions without the cap layer. The traditional Thomson—Haskell matrix approach to computing the normal modes is reformulated so that numerical problems encountered at high frequencies are avoided and numerical results of the locked mode approximation are given.