Shear Velocity Model of Alaska via Joint Inversion of Rayleigh Wave Ellipticity, Phase Velocities, and Receiver Functions across the Alaska Transportable Array

Through the Alaska Transportable Array deployment of over 200 stations, we create a 3-D tomographic model of Alaska with sensitivity ranging from the near surface (<1 km) into the upper mantle (similar to 140 km). We perform a Markov chain Monte Carlo joint inversion of Rayleigh wave ellipticity and phase velocities, from both ambient noise and earthquake measurements, along with receiver functions to create a shear wave velocity model. We also use a follow-up phase velocity inversion to resolve interstation structure. By comparing our results to previous tomography, geology, and geophysical studies we are able to validate our findings and connect localized near-surface studies with deeper, regional models. Specifically, we are able to resolve shallow basins, including the Copper River, Cook Inlet, Yukon Flats, Nenana, and a variety of other shallower basins. Additionally, we gain insight on the interaction between the upper mantle wedge, asthenosphere, and active and nonactive volcanism along the Aleutians and Denali volcanic gap, respectively. We observe thicker crust beneath the Brooks Range and south of the Denali fault within the Wrangellia Composite Terrane and thinner crust in the Yukon Composite Terrane in interior Alaska. We also gain new perspective on the Wrangell Volcanic Field and its interaction between surrounding asthenosphere and the Yakutat Terrane.