Constraining velocity and density contrasts across the crust-mantle boundary with receiver function amplitudes

This paper investigates the ability of P-wave receiver functions to constrain both the velocity and density contrasts across the Moho. Expressions for receiver function amplitudes corresponding to a layer over a half-space are extended to media with depth-dependent properties by explicitly accounting for free-surface reflection coefficients. Forward numerical computations show that receiver function amplitudes become frequency-dependent for depth-dependent structures and that, for a given frequency, wavelengths associated to Ps converted phases are larger than wavelengths associated to multiply reverberated phases. When consistent measurements of the amplitudes of Ps phase and multiples can be obtained, the extended expressions suffice to account for the observed amplitudes. Receiver function amplitudes are sensitive to both velocity and density contrasts across the Moho, and a two-step grid-search procedure is proposed to recover the contrasts from their measurement. The near-surface velocity is recovered in the first step from the amplitude of the direct P wave and then used in the second step to work out the relative density and S-wave velocity contrast from the amplitude of the secondary phases. Examples from central Spain (PAB) and the Indian Shield (HYB) are discussed and demonstrate that receiver function amplitudes can constrain contrasts across the Moho accurately enough to be utilized in geological interpretation.