Spectral-element simulations of global seismic wave propagation - I. Validation

We use a spectral-element method to simulate seismic wave propagation throughout the entire globe. The method is based upon a weak formulation of the equations of motion and combines the flexibility of a finite-element method with the accuracy of a global pseudospectral method. The finite-element mesh honours all first-and second-order discontinuities in the earth model. To maintain a relatively constant resolution throughout the model in terms of the number of grid points per wavelength, the size of the elements is increased with depth in a conforming fashion, thus retaining a diagonal mass matrix. In the Earth's mantle and inner core we solve the wave equation in terms of displacement, whereas in the liquid outer core we use a formulation based upon a scalar potential. The three domains are matched at the inner core and core-mantle boundaries, honouring the continuity of traction and the normal component of velocity. The effects of attenuation and anisotropy are fully incorporated. The method is implemented on a parallel computer using a message passing technique. We benchmark spectral-element synthetic seismograms against normal-mode synthetics for a spherically symmetric reference model. The two methods are in excellent agreement for all body-and surface-wave arrivals with periods greater than about 20 s.