Along-strike variation in slab geometry at the southern Mariana subduction zone revealed by seismicity through ocean bottom seismic experiments

Author Posting. © The Authors, 2019. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 218(3), (2019): 2122-2135, doi: 10.1093/gji/ggz272.We have conducted the first passive Ocean Bottom Seismograph (OBS) experiment near the Challenger Deep at the southernmost Mariana subduction zone by deploying and recovering an array of 6 broad-band OBSs during December 2016–June 2017. The obtained passive-source seismic records provide the first-ever near-field seismic observations in the southernmost Mariana subduction zone. We first correct clock errors of the OBS recordings based on both teleseismic waveforms and ambient noise cross-correlation. We then perform matched filter earthquake detection using 53 template events in the catalogue of the US Geological Survey and find >7000 local earthquakes during the 6-month OBS deployment period. Results of the two independent approaches show that the maximum clock drifting was ∼2 s on one instrument (OBS PA01), while the rest of OBS waveforms had negligible time drifting. After timing correction, we locate the detected earthquakes using a newly refined local velocity model that was derived from a companion active source experiment in the same region. In total, 2004 earthquakes are located with relatively high resolution. Furthermore, we calibrate the magnitudes of the detected earthquakes by measuring the relative amplitudes to their nearest relocated templates on all OBSs and acquire a high-resolution local earthquake catalogue. The magnitudes of earthquakes in our new catalogue range from 1.1 to 5.6. The earthquakes span over the Southwest Mariana rift, the megathrust interface, forearc and outer-rise regions. While most earthquakes are shallow, depths of the slab earthquakes increase from ∼100 to ∼240 km from west to east towards Guam. We also delineate the subducting interface from seismicity distribution and find an increasing trend in dip angles from west to east. The observed along-strike variation in slab dip angles and its downdip extents provide new constraints on geodynamic processes of the southernmost Mariana subduction zone.We express our appreciation to the science parties and crew members of the R/V Shiyan 3 for deployment and collection of the OBS instruments during the Mariana expeditions. This study is supported by the Hong Kong Research Grant Council Grants (No. 14313816), Faculty of Science at CUHK, Chinese Academy of Sciences (No. Y4SL021001, QYZDY-SSW-DQC005, 133244KYSB20180029), the National Natural Science Foundation of China (No. 41890813, 91628301, 41676042, U1701641, 41576041, 91858207 and U1606401), the National Key R&D Program of China (2018YFC0309800 and 2018YFC0310100). Generic Mapping Tools (Wessel & Smith 1991) and PSSAC (developed by Prof Lupei Zhu) are used for data analysis and figure preparation in this study. Constructive comments from Dr Lidong Bie and two anonymous reviewers are helpful in improving the manuscript