Transition zone velocity gradients and the 520-km discontinuity

Stacks of long-period Global Digital Seismograph Network (GDSN) seismograms at 110 degrees to 180 degrees epicentral distance reveal precursors to SS that result from underside reflections off upper mantle seismic discontinuities. The 410- and 660-km discontinuities are obvious in these stacks, but identification and modeling of other transition zone discontinuities are complicated by sidelobes from the 410- and 660-km reflections. These sidelobes result from the limited bandwidth of the GDSN instrument responses and the effect of crustal reverberations on the SS reference phase. The crustal effects can be minimized by restricting the records to oceanic bounce points where the NG-km-thick crust has little effect on the long-period waveforms. Over 2000 long-period, transverse-component seismograms with oceanic SS bounce points recorded by the GDSN from 1976 to 1991 are manually edited, aligned on SS, and then stacked using a new procedure that weights the records by data quality. The resulting image shows a clear reflection from a 520-km discontinuity that cannot be explained as a sidelobe artifact, confirming earlier results of Shearer [1990, 1991] and Revenaugh and Jordan. [1991]. By stacking along the expected travel time curves for discontinuity phases, the time versus range image of the precursor wave field is reduced to a single trace that measures upper mantle reflectivity versus time. The features in this reflectivity profile are sensitive to the brightness and depth of the transition zone discontinuities and to the steepness of the velocity gradients between the interfaces. Using geometrical ray theory and assuming a constant velocity versus density scaling relationship, I fit this reflectivity profile with velocity models of the upper mantle using both forward modeling and direct inversion. The inverse problem. is addressed by performing a deconvolution of the profile with the SS reference phase (after a correction for attenuation), followed by a direct mapping of reflectivity versus time into velocity versus depth. Velocity-depth profiles resulting from these procedures are roughly in agreement with standard upper mantle velocity models, except that the SS precursor data require a minor discontinuity near 520 km and a steeper gradient just below the 660-km discontinuity. Estimated discontinuity shear impedance changes are 6.7 +/- 0.1% at 420 km, 2.9 +/- 0.7% at 519 km, and 9.9 +/- 1.5% at 663 km. The impedance change near 520 km is consistent with current mineral physics results for the olivine beta to gamma phase change and places constraints on the fraction of olivine in the transition zone.