 Solution

Write $q={\partial }_z\bar{\rho }>0$ for the unstable, horizontally averaged density gradient. A turbulent displacement over a mixing length of order $W$ encounters a density contrast $\Delta \rho \sim Wq$. Balancing its buoyancy acceleration with the inertial acceleration givesA turbulent diffusivity for density is therefore $D_T=K{u}_{\ast }W=K{W}^2\sqrt{gq/\bar{\rho }}$, absorbing dimensionless factors into $K>0$. The upward density flux is $J_{\rho }=⟨w^{\prime }{\rho }^{\prime }⟩=-D_{T}q$: turbulent exchange carries denser water downwards and lighter water upwards. Defining $F_{\rho }=-J_{\rho }$ as the positive downward flux yields the buoyancy-driven turbulent density diffusion lawThis is a flux per unit horizontal area, with units of mass per area per time. The total density transport across a horizontal section is $W^2{F}_{\rho }$.

The assumptions are established high-Reynolds number turbulence, width-limited eddies, an instantaneous local balance between buoyancy and inertia, approximately constant $K$, negligible molecular diffusion, and zero mean vertical volume flux. The last condition follows from incompressibility and impermeable sidewalls in a closed tube. Density variations and the eddy statistics must also vary slowly enough for a local mixing-length closure to be meaningful.

Horizontal averaging of mass conservation, with upward $z$, gives ${\bar{\rho }}_t=-{\partial }_z{J}_{\rho }={\partial }_z{F}_{\rho }$. Differentiating once more supplies the requested evolution of the density gradient:Under the Boussinesq approximation, replace $\bar{\rho }$ in the inertial prefactor by the constant reference density ${\rho }_0$; retain the changing density gradient inside the flux.