V_S and density structure beneath the Colorado Plateau constrained by gravity anomalies and joint inversions of receiver function and phase velocity data

The Colorado Plateau is a physiographic province in the western US with an average elevation of ∼1.9 km where, in contrast to neighboring provinces, there is little evidence of large scale tectonic deformation or magmatism. Recent availability of Earthscope/USArray seismic data allow us to better examine the crust and upper mantle structure beneath the region and test proposed explanations for the plateau's uplift and relative stability. Using phase velocities for fundamental mode Rayleigh waves and P receiver functions, we perform over 800 joint inversions for 1-D shear wave velocity V_S profiles sampling the plateau and surrounding regions down to 150 km depth. We image a sharp change in crustal thickness at the western edge of the Colorado Plateau with a more gradual increase eastward moving into the Rocky Mountains. A relatively thick (≳100 km) lithosphere beneath the plateau extends into the Rocky Mountains to the north. We use empirical scaling relations to estimate densities from our V_S results, and predict the associated gravity anomalies, which are inconsistent with the observed distribution of the Bouguer gravity anomalies. We somewhat reconcile the prediction and observations by assuming that lateral density variations below 50 km can be ignored and the lithospheric root is therefore neutrally buoyant. While there is some evidence for small scale convection and lithospheric removal at its edges, the shape of the lithospheric mantle anomaly is consistent with a large scale uplift of the plateau by heating since removal of the Farallon slab. We conclude that the lithospheric root is key to the long term stability of the Colorado Plateau, leading to a colder, stronger crust.