Spatially heterogeneous diapycnal mixing in the abyssal ocean: A comparison of two parameterizations to observations

The spatial distributions of the diapycnal diffusivity predicted by two abyssal mixing schemes are compared to each other and to observational estimates based on microstructure surveys and large-scale hydrographic inversions. The parameterizations considered are the tidal mixing scheme by Jayne, St. Laurent and co-authors (JSL01) and the Roughness Diffusivity Model (RDM) by Decloedt and Luther. Comparison to microstructure surveys shows that both parameterizations are conservative in estimating the vertical extent to which bottom-intensified mixing penetrates into the stratified water column. In particular, the JSL01 exponential vertical structure function with fixed scale height decays to background values much nearer topography than observed. JSL01 and RDM yield dramatically different horizontal spatial distributions of diapycnal diffusivity, which would lead to quite different circulations in OGCMs, yet they produce similar basin-averaged diffusivity profiles. Both parameterizations are shown to yield smaller basin-mean diffusivity profiles than hydrographic inverse estimates for the major ocean basins, by factors ranging from 3 up to over an order of magnitude. The canonical 10(-4) m(2) s(-1) abyssal diffusivity is reached by the parameterizations only at depths below 3 km. Power consumption by diapycnal mixing below 1 km of depth, between roughly 32 degrees S and 48 degrees N, for the RDM and JSL01 parameterizations is 0.40 TW & 0.28 TW, respectively. The results presented here suggest that present-day mixing parameterizations significantly underestimate abyssal mixing. In conjunction with other recently published studies, a plausible interpretation is that parameterizing the dissipation of bottom-generated internal waves is not sufficient to approximate the global spatial distribution of diapycnal mixing in the abyssal ocean.