Radial Anisotropy in Receiver Function H−κ Stacks

Abstract Receiver functions can be used to estimate the Moho depth (H) and ratio of P to S wavespeed (α/β or κ) in the crust. This is commonly done by grid search, forward modeling travel times to produce so-called “H-κ” stacks of receiver function amplitude. However, radial anisotropy in the crust, which can be significant, is almost never considered in this process. Here, we show that radial anisotropy changes the H-κ stack, biasing interpretations of crustal structure by introducing errors up to ∼3% in H and ∼1% in κ for commonly observed anisotropy magnitudes. We propose a simple method to correct H-κ stacks by incorporating radial anisotropy in the forward calculation. Synthetic tests show that this approach almost completely removes error caused by radial anisotropy. We show examples of this procedure with stations in the eastern United States. We provide readers with code to construct radially anisotropic H-κ stacks.