The surface wave full ray theory (FRT) is an efficient tool to calculate synthetic waveforms of surface waves. It combines the concept of local modes with exact ray tracing as a function of frequency, providing a more complete description of surface wave propagation than the widely used great circle approximation (GCA). The purpose of this study is to evaluate the ability of the FRT approach to model teleseismic long-period surface waveforms (T ∼ 45–150 s) in the context of current 3-D Earth models to empirically assess its validity domain and its scope for future studies in seismic tomography. To achieve this goal, we compute vertical and horizontal component fundamental mode synthetic Rayleigh waveforms using the FRT, which are compared with calculations using the highly accurate spectral element method. We use 13 global earth models including 3-D crustal and mantle structure, which are derived by successively varying the strength and lengthscale of heterogeneity in current tomographic models. For completeness, GCA waveforms are also compared with the spectral element method. We find that the FRT accurately predicts the phase and amplitude of long-period Rayleigh waves (T ∼ 45–150 s) for almost all the models considered, with errors in the modelling of the phase (amplitude) of Rayleigh waves being smaller than 5 per cent (10 per cent) in most cases. The largest errors in phase and amplitude are observed for T ∼ 45 s and for the three roughest earth models considered that exhibit shear wave anomalies of up to ∼20 per cent, which is much larger than in current global tomographic models. In addition, we find that overall the GCA does not predict Rayleigh wave amplitudes well, except for the longest wave periods (T ∼ 150 s) and the smoothest models considered. Although the GCA accurately predicts Rayleigh wave phase for current earth models such as S20RTS and S40RTS, FRT's phase errors are smaller, notably for the shortest wave periods considered (T ∼ 45 s and T ∼ 60 s). This suggests that the FRT approach is a useful means to build the next generation of elastic and anelastic surface wave tomography models. Finally, we observe a clear correlation between the FRT amplitude and phase errors and the roughness of the models. This allows us to quantify the limits of validity of the FRT in terms of model roughness thresholds, which can serve as useful guides in future seismic tomographic studies.
|Original language||English (US)|
|Number of pages||14|
|Journal||Geophysical Journal International|
|State||Published - Feb 10 2016|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: We sincerely thank Eric Debayle, an anonymous reviewer and the editor Andrea Morelli for their detailed reviews, which helped improve our manuscript. This research was carried out on the High Performance Computing Cluster supported by the Research and Specialist Computing Support services at the University of East Anglia and on Archer, the UK's National Supercomputing Service. Some figures were built using Generic Mapping Tools (GMT; Wessel & Smith 1998). This work was supported by the European Commission's Initial Training Network project QUEST (contract FP7-PEOPLE-ITN-2008-238007, http://www.quest-itn.org) and the University of East Anglia. The Incorporated Research Institutions for Seismology (IRIS) was used to access the waveform data from the IRIS/IDA network II (http://dx.doi.org/doi:10.7914/SN/II) and IU (http://dx.doi.org/doi:10.7914/SN/IU). AMGF also thanks funding by the Fundacao para a Ciencia e Tecnologia (FCT) project AQUAREL (PTDC/CTE-GIX/116819/2010).