“…In order to provide a model with resolution of Vs and Vp/Vs in the upper few km, we combine the complementary sensitivities of Rayleigh‐wave phase velocities (upper crust), ellipticity (upper few km), and the initial pulse of teleseismic receiver functions (shallow Vp/Vs ratio and shallow interfaces) to create a self‐consistent model at the regional scale across southern California. The idea to combine receiver functions and surface wave data in a Bayesian joint inversion to determine Vs and Vp/Vs is relatively new (Dreiling et al., 2020; Ojo et al., 2019), and only recently shown to be promising in resolving near‐surface Vs and Vp/Vs in sediments (Li et al., 2019). By including Vp/Vs as a parameter we are able to fit receiver functions on a regional scale for the first time across 231 Southern California stations, including in basins where receiver functions have long been discarded as nuisance signals or “corrected” with ad‐hoc models, as reverberations overprint Moho and other crustal signatures (e.g., Yeck et al., 2013).…”