On the validity of the planar wave approximation to compute synthetic seismograms of teleseismic body waves in a 3-D regional model

Injection methods are a very efficient means to compute synthetic seismograms of short-period teleseismic body waves in 3-D regional models. The principle is to inject an incident teleseismic wavefield inside a regional 3-D Cartesian spectral-element grid. We have developed an opern-source package that allows us to inject either an incident plane wave, computed with a frequency–wavenumber method, or the complete wavefield, computed in a spherically symmetric reference earth model with AxiSEM. The computations inside the regional spectral-element grid are performed with SPECFEM3D_Cartesian. We compare the efficiency and reliability of the two injection methods for teleseismic P waves, considering a wide range of epicentral distance and hypocentral depths. Our simulations demonstrate that in practice the effects of wave front and Earth curvature are negligible for moderate size regional domains (several hundreds of kilometres) and for periods larger than 2 s. The main differences observed in synthetic seismograms are related to secondary phases that have a different slowness from the one of the reference P phase.