Shear-wave velocity structure beneath Alaska from a Bayesian joint inversion of Sp receiver functions and Rayleigh wave phase velocities

Abstract

The present-day lithospheric structure of Alaska is the result of a unique tectonic history of subduction and terrane accretion that controls upper plate thickness and rheology. To provide new constraints on the structure of the crust and upper mantle beneath Alaska, we jointly inverted Sp receiver functions and Rayleigh wave phase velocities to calculate shear-wave velocity profiles. Robust Sp receiver functions were obtained using a broad range of frequencies (2–100 s), time-domain deconvolution, and K-means cluster analysis. Tests of the Bayesian joint inversion with synthetic data illustrate that Sp receiver functions enhance the resolution of the velocity gradients at the Moho and the lithosphere-asthenosphere boundary, while Rayleigh surface waves provide information about absolute velocities. Our results show that in central Alaska, above the shallow slab, the continental lithosphere is thinnest (∼60 km) and the asthenosphere has its lowest velocities. This zone coincides with the Denali Volcanic Gap. The continental lithosphere thickens to the north beneath the Brooks Range and the northern Arctic Alaska terrane, reaching values of 110 to 130 km, with high lithospheric velocities that are comparable to Archean cratons. This pattern is consistent with a northward decrease in upper plate modification by melt and volatiles derived from the slab, in addition to intrinsic mantle velocity and viscosity differences between inherited lithospheric terranes. Lithospheric and asthenospheric velocities are not significantly different inside and outside of the Denali Volcanic Gap, but at the boundaries of this region lithospheric thickness increases rapidly to the north and gradually to the south. In the south, the subducting Yakutat plate is thicker (∼100–120 km) than the subducting Pacific plate (∼80–90 km), likely due to its thicker crust.