Mantle transition zone shear velocity gradients beneath USArray

Broadband P-to-s scattering isolated by teleseismic receiver function analysis is used to investigate shear velocity (VS) gradients in the mantle transition zone beneath USArray. Receiver functions from 2244 stations were filtered in multiple frequency bands and migrated to depth through P and S tomography models. The depth-migrated receiver functions were stacked along their local 410 and 660 km discontinuity depths to reduce stack incoherence and more accurately recover the frequency-dependent amplitudes of P410s and P660s. The stacked waveforms were inverted for one-dimensional VS between 320 and 840 km depth. First, a gradient-based inversion was used to find a least-squares solution and a subsequent Monte Carlo search about that solution constrained the range of VS profiles that provide an acceptable fit to the receiver function stacks. Relative to standard references models, all the acceptable models have diminished VS gradients surrounding the 410, a local VS gradient maximum at 490–500 km depth, and an enhanced VS gradient above the 660. The total 410 VS increase of 6.3% is greater than in reference models, and it occurs over a thickness of 20 km. However, 60% of this VS increase occurs over only 6 km. The 20 km total thickness of the 410 and diminished VS gradients surrounding the 410 are potential indications of high water content in the regional transition zone. An enhanced VS gradient overlying the 660 likely results from remnants of subduction lingering at the base of the transition zone. Cool temperatures from slabs subducted since the late Cretaceous and longer-term accumulation of former ocean crust both may contribute to the high gradient above the 660. The shallow depth of the 520 km discontinuity, 490–500 km, implies that the regional mean temperature in the transition zone is 110–160 K cooler than the global mean. A concentrated Vs gradient maximum centered near 660 km depth and a low VS gradient below 675 km confirms that the ringwoodite to perovskite+magnesiowüstite reaction is the dominant cause of the 660 in the region.