EFFECTS OF TRENCH-ZONE SCATTERING ON RECEIVER FUNCTIONS OVER A SUBDUCTION ZONE: A 3-D FINITE-DIFFERENCE MODELLING STUDY

Abstract

As Morozov [Morozov, I. B. (2004). Crustal scattering and some artefacts in receiver function images. Bull. Seismol. Soc. Am., 94 (4), 1492–1499.] suggested, for a teleseismic array targeting subducting crust in a zone of active subduction, scattering from the strong horizontal velocity heterogeneity beneath the trench zone itself produces subhorizontally-propagating waves that should be observed as coherent dipping events in receiver functions (RF). Due to similar RF delay times and moveouts, these events could be difficult to distinguish from backscattered P- and S-wave modes. To further verify this suggestion, we performed a full-waveform, 3-D visco-elastic finite-difference modelling of teleseismic wave propagation within a simplified model of a subduction zone. The synthetics show strong scattering from the area beneath the trench, dominated by the mantle and crustal P-waves propagating at 6.2–8.1 km/s and slower. These scattered waves occupy the same time and moveout intervals as the backscattered converted modes, and also have similar amplitudes. Although their amplitude decay characters are different, the uncertainty in the knowledge of the velocity and density structure of the subduction zone could make distinguishing between these modes difficult. However, under minimal assumptions, recent observations of receiver function amplitudes decreasing away from the trench support the interpretation of (sub-) trench-zone scattering.