Statistics and frequency-domain moveout for multiple-taper receiver functions

The multiple-taper correlation (MTC) algorithm for the estimation of teleseismic receiver functions (RFs) has desirable statistical properties. This paper presents several adaptations to the MTC algorithm that exploit its frequency-domain uncertainty estimates to generate stable RFs that include moveout corrections for deeper interfaces. Narrow-band frequency averaging implicit in spectral cross-correlation restricts the MTC-based RF estimates to resolve Ps converted phases only at short delay times, appropriate to the upper 100 km of Earth's lithosphere. The Ps conversions from deeper interfaces can be reconstructed by the MTC algorithm in two ways. Event cross-correlation computes a cross-correlation of single-taper spectrum estimates for a cluster of events rather than for a set of eigenspectrum estimates of a single P coda. To extend the reach of the algorithm, pre-stack moveout corrections in the frequency domain preserves the formal uncertainties of the RF estimates, which are used to weight RF stacks. Moving-window migration retains the multiple-taper approach, but cross-correlates the P-polarized motion with time-delayed SH and SV motion to focus on a Ps phase of interest. The frequency-domain uncertainties of bin-averaged RFs do not translate directly into the time domain. A jackknife over data records in each bin stack offers uncertainty estimates in the time domain while preserving uncertainty weighting in the frequency-domain RF stack.