Receiver function studies of Northern Sumatra

The Northern Sumatra subduction zone is distinguished by the occurrence of the 2004 Sumatra-Andaman megathrust earthquake and has a peculiar subduction of two major bathymetric structures; the Investigator fracture zone and the Wharton fossil ridge. Four stations in Northern Sumatra (BSI, PSI, PPI, GSI) and two stations in Malaysia (KUM and KOM) have been selected to construct migrated images based on receiver functions (RF) in order to study Earth structure and subduction processes in the region. Waveforms from 304 teleseismic earthquakes with Mb >5.5 and a distance range of 30o to 95o recorded from April 2006 to December 2008 were used for the analysis. The number of RF for each station varies from 20 to 192 depending on the signal/noise ratio. The computed RF clearly show pS conversions at major seismic velocity discontinuities associated with the subduction process where the Moho is visible at 5.5, 4, 3.5, and 2 sec for BSI, PSI, PPI, and GSI, respectively. RF for KUM and KOM have only conversions at the Moho near ~4 sec. The subducted slab is visible below Sumatra as a positive amplitude conversion preceded by a negative one, which we interpret as a low-velocity structure above the subducted slab. RF for PSI located at Toba supervolcano reveal pockets of low-velocity zones extending from a ~50 km depth down to the subducted slab. Forward modellings of RF suggest that seismic velocity contrasts can reach ~18% that is in accordance with previous local tomographic studies.