Abstract: Recently, we presented a parameterization of an arbitrary tracer 3D mesoscale flux that describes both diabatic and adiabatic regimes without using arbitrary tapering functions. However, we did not parameterize the mesoscale diffusivity, which is the subject of this work. A key difference between the present and previous diffusivity parameterizations is that in the latter, the two main ingredients, mesoscale drift velocity and eddy kinetic energy, were not parameterized but determined using p…