Characterizing the 410 km Discontinuity Low‐Velocity Layer Beneath the LA RISTRA Array in the North American Southwest

Receiver functions recorded by the 54-station 920 km long Program for Array Seismic Studies of the Continental Lithosphere–Incorporated Research Institutions for Seismology Colorado Plateau/Rio Grande Rift Seismic Transect Experiment (LA RISTRA) line array display a pervasive negative polarity P to S conversion (Pds) arrival preceding the positive polarity 410 km discontinuity arrival. These arrivals are modeled as a low-velocity layer atop the 410 km discontinuity (410-LVL) and are inverted for a velocity profile via a grid search using a five-parameter linear gradient velocity model. Model parameter likelihood and correlations are assessed via calculation of one- and two-dimensional marginal posterior probability distributions. The maximum likelihood model parameter values found are top velocity gradient thickness of 0.0 km with a 4.6% (−0.22 km/s) shear velocity reduction, a 19.8 km constant velocity layer, and bottom gradient thickness of 25.0 km with a 3.5% (+0.17 km/s) shear velocity increase. The estimated mean thickness of the 410-LVL is 32.3 km. The top gradient of the 410-LVL is sharp within vertical resolution limits of P to S conversion (km), and the diffuse 410 km velocity gradient is consistent with hydration of the olivine-wadsleyite phase transformation. The 410-LVL is interpreted as a melt layer created by the Transition Zone Water Filter model. Two secondary observations are found: (1) the 410-LVL is absent from the SE end of the array and (2) an intermittent negative polarity P525s arrival is observed. We speculate that upper mantle shear velocity anomalies above the 410 km discontinuity may manifest Rayleigh-Taylor instabilities nucleated from the 410-LVL melt layer that are being shed upward on time scales of tens of millions of years

Crustal structure and thickness along the Yellowstone hot spot track : evidence for lower crustal outflow from beneath the eastern Snake River Plain

Receiver functions from seismic stations about the Yellowstone hot spot track are migrated to depth using a Vp/Vs map constructed from stacking of the direct and free surface Moho reverberations (i.e., H‐K analysis) and a shear velocity tomogram constructed from surface wave measurements. The thickest crust (48–54 km) resides in the Wyoming province beneath the sampled Laramide age blocks, and the thinnest crust (32–37 km) resides in the Montana Basin and Range province. The eastern Snake River Plain (ESRP) crust is thickest (47 km) at its NE end beneath the young calderas and thinnest (40 km) at its SW end beneath the older Twin Falls caldera. Two ESRP crustal thickness domains are found: (1) at the older Twin Falls and Picabo calderas, the mean ESRP crust is 4 km thicker with respect to its margins and (2) adjacent to the Heise caldera field, the mean ESRP crust is 4 km thicker with respect to its SE margin crust but no thicker with respect to its NW margin crust. This lobe of anomalously thick crust is explained as resulting from lower crustal outflow from beneath the Heise caldera field. Confirmation of these crustal thickness variations is provided by inspection of common conversion point (CCP) stacks that delineate several secondary features: the top of a thick high‐velocity (3.9 km/s) lower crust layer within the Wyoming province up to 17 km thick and a paired negative and positive amplitude arrival at 12 km depth and 18 km depth beneath the Yellowstone Caldera. This paired arrival would be consistent with a low‐velocity zone perhaps associated with magma staging beneath the caldera. Our most important finding is that the magmatic loads injected into the ESRP crust over the last 4–12 Myr, in tandem with the ESRP crustal viscosity structure, have been sufficient to drive significant outflow of the ESRP lower crust