Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles

Internal wave - wave interaction theories and observations support a parameterization for the turbulent dissipation rate epsilon and eddy diffusivity K that depends on internal wave shear [V-z(2)] and strain [xi(2)(z)] variances. Its latest incarnation is applied to about 3500 lowered ADCP/CTD profiles from the Indian, Pacific, North Atlantic, and Southern Oceans. Inferred diffusivities K are functions of latitude and depth, ranging from 0.03 x 10(-4) m(2) s(-1) within 2 degrees of the equator to (0.4 - 0.5) x 10(-4) m(2) s(-1) at 50 degrees-70 degrees. Diffusivities K also increase with depth in tropical and subtropical waters. Diffusivities below 4500-m depth exhibit a peak of 0.7 x 10(-4) m(2) s(-1) between 20 degrees and 30 degrees, latitudes where semidiurnal parametric subharmonic instability is expected to be active. Turbulence is highly heterogeneous. Though the bulk of the vertically integrated dissipation integral(epsilon) is contributed from the main pycnocline, hotspots in integral(epsilon) show some correlation with small-scale bottom roughness and near-bottom flow at sites where strong surface tidal dissipation resulting from tide - topography interactions has been implicated. Average vertically integrated dissipation rates are 1.0 mW m(-2), lying closer to the 0.8 mW m(-2) expected for a canonical (Garrett and Munk) internal wave spectrum than the global-averaged deep-ocean surface tide loss of 3.3 mW m(-2).