In a previous study (Pelegri and Csanady, 1991) we found upward entrainment and two-way exchange between the upper thermocline and surface layers of the Gulf Stream. Here we use isopycnic coordinates to analyze the mechanism of diapycnal shear-induced mixing in western boundary currents. A hydrographic section across the Gulf Stream off the Mid-Atlantic Bight is used for a case study. It shows that the upper thermocline layers of the Gulf Stream are characterized by relatively small (16 m averaged) gradient Richardson numbers, Ri. Analysis of the Jacobian, J = partial-derivativez/partial-derivativerho (with z the depth and rho the potential density), total vorticity and potential vorticity fields shows that this region has low Jacobian values associated with high values of potential vorticity, as required theoretically for the case of dominant diapycnal mixing. The density tendency, w(rho) = Drho/Dt, is calculated from the vertical gradient of the Reynolds density flux, the latter parameterized as density eddy diffusivity divided by the Jacobian. We use Munk and Anderson's (1948) expression for the eddy diffusivity coefficient, modified to conform to data of Ueda et al. (1981), maintaining a Ri-3/2 dependence for large Ri. From this we obtain an explicit expression for w(rho) in terms of the Ri and J distributions. The distributions of w(rho), its diapycnal gradient, partial-derivitivew(rho)/partial-derivitiverho, and the diapycnal velocity, w(d) = Jw(rho), are then calculated. Regions of high diapycnal convergence (large negative partial-derivitivew(rho)/partial-derivitiverho values) coincide with the location of anomalously low Jacobians, in agreement with the idea that diapycnal mixing is the consequence of anomalies produced during frontogenesis in some phase of the meanders. We suggest that this process, taking place intermittently in space and time, is responsible for two-way exchange in western boundary currents.