Long-period (100-250 sec.) Love and Rayleigh waves are used to map heterogeneity and azimuthal anisotropy in the upper mantle. Spherical harmonic descriptions of anisotropy up to l=m=3 and 2&thgr; are derived. Azimuthal anisotropy obtains values as high as 1 1/2%. There is good correlation of fast Rayleigh wave directions with upper mantle return flow directions derived from kinematic considerations. This is consistent with the a-axis of olivine being aligned in the flow direction. The main differences between the flow models and the Rayleigh wave azimuthal variation maps occur in the vicinity of hotspots. Hawaii, for example, appears to perturb the return flow. There is strong correlation of the geoid and surface wave velocity at l=4 and 5. Slow regions at this scale are associated with geoid highs and high heat flow, consistent with upwelling convective flow or with isostatically compensated regions of low density. The correlations of azimuthal anisotropy with upper mantle return flow directions, rather than with plate directions, suggests that part of the return flow is in the upper mantle and this, in turn, implies a low viscosity channel.