Parameterizing mesoscale eddies with residual and Eulerian schemes, and a comparison with eddy-permitting models

The residual scheme is compared both to a conventional parameterization that uses a skew diffusion (or equivalently advection by a skew velocity), and to eddy-permitting calculations. Although in principle almost equivalent to certain forms of skew flux schemes, the residual formulation is found to have certain practical advantages over the conventional scheme in some circumstances, and in particular near the upper boundary where conventional schemes are sensitive to the choice of tapering but the residual scheme is less so. The residual scheme also enables the horizontal viscosity - which is mainly applied to maintain model stability - to be reduced. Finally, the residual scheme is somewhat easier to implement, and the tracer transport is easier to interpret. On the other hand, the residual scheme gives, at least formally, a transformed velocity, not the Eulerian velocity and will not be appropriate in all circumstances.