Kawakatsu and Abe (2016) have highlighted the potential complicating effect of sediment reverberations on the analysis and interpretation of crust and mantle phases inferred from receiver functions analyzed from ocean-bottom seismograms. In their comment, they identify resonant peaks in the power spectrum at one of the stations, T06, and demonstrate with synthetic modeling how sedimentinduced resonances can cause instability in the recovered receiver-function (RF) traces. They also request a detailed explanation of how LQT rotation is conducted, and why its use leads to stable receiver functions in the analysis of . We welcome this query as an opportunity to highlight certain technical aspects of the data-analysis procedures used in . Our methods derive partly from methods recommended by previous studies of receiver functions estimated from seismic seafloor data, particularly the use of the modal wavefield decomposition (which we approximated by the LQT rotation) to suppress reverberation signals in the overlying water column.
Approach Used in Calculating Receiver FunctionWe use multiple-taper correlation (MTC) to estimate receiver functions [Park and Levin, 2000], which estimates Ps converted-waves from the portions of the P, SV, and SH motions that are mutually coherent in narrow band-passes (J. Park and V. Levin, Statistics and frequency-domain moveout for multiple-taper receiver functions, submitted to Geophysical Journal International, 2016), discarding the incoherent portion of the wavefield. We apply moveout corrections to the MTC RFs before stacking in the frequency domain to enhance primary Ps conversions from deep interfaces, and suppress reverberations [Helffrich, 2006;Bianchi et al, 2010] (Park and Levin, submitted manuscript, 2016). Finally, synthetic modeling suggests that, even at our noisiest stations, defocusing of seismic energy within the sediment layer leads to reverberations of much smaller amplitude than predicted by simple 1-D models. Real seafloor seismic observations exhibit other evidence for resonant behavior that plausibly arise from ambient motions at the seafloor, because they can be detected in the seismic noise prior to P-wave arrival. Our identification of the deep LAB interface relies on the observed moveout of its negative-amplitude Ps with epicentral distance, when we analyze P-coda for receiver functions. This Ps conversion (with moveout) is absent when preevent data are analyzed.
Adjustable Empirical LQT RotationWe use an empirical LQT rotation as an approximate method of wave-field decomposition [e.g., Vinnik, 1977;Bostock, 1998;Reading et al., 2003;Svenningsen, 2004]. An explicit modal wave-field decomposition is advocated by Bostock and Trehu [2012] to reduce contamination of water and sediment reverberations. The modal decomposition relies on either a collocated seafloor pressure sensor or prior knowledge of the sediment layer to effect its correction. The LQT rotation is crude by comparison, but can be scaled to less than a full rotation.We used an empirical r...