From the reprocessed data set of Engdahl and co-workers we have carefully selected matching P and S data for tomographic imaging. We assess data and model error and conclude that our S model uncertainty is twice that of the P model. We account for this in our comparison of the perturbations in P and S-wavespeed. In accord with previous studies we find that P and S perturbations are positively correlated at all depths. However, in the deep mantle systematic differences occur between regions that have undergone subduction in the last 120 million years and those that have not. In particular, below 1500 km depth ∂ ln Vs/∂ ln Vp is significantly larger in mantle regions away from subduction than in mantle beneath convergent margins. This inference is substantiated by wavespeed analyses with random realizations of the slab/non-slab distribution. Through much of the mantle there is no significant correlation between bulk sound and S-wave perturbations, but they appear to be negatively correlated between 1700 and 2100 km depth, which is also where the largest differences in ∂ ln Vs/∂ ln Vp occur. This finding supports convection models with compositional heterogeneity in the lowermost mantle. ¿ 2001 American Geophysical Union