A new gauge-invariant method for diagnosing eddy diffusivities

TL;DR

This paper introduces a gauge-invariant method for diagnosing eddy diffusivities in ocean models, addressing rotational flux contamination and improving parameterisation accuracy through a non-local force function approach.

Contribution

A novel methodology using eddy force functions to diagnose eddy diffusivities directly from flux divergences, avoiding rotational flux contamination.

Findings

Strictly down-gradient schemes have limited success.

Negative diffusivities are common near jets.

Diffusivity magnitude correlates with eddy energy.

Abstract

Coarse resolution numerical ocean models must typically include a parameterisation for mesoscale turbulence. A common recipe for such parameterisations is to invoke down-gradient mixing, or diffusion, of some tracer quantity, such as potential vorticity or buoyancy. However, it is well known that eddy fluxes include large rotational components which necessarily do not lead to any mixing; eddy diffusivities diagnosed from unfiltered fluxes are thus contaminated by the presence of these rotational components. Here a new methodology is applied whereby eddy diffusivities are diagnosed directly from the eddy force function. The eddy force function depends only upon flux divergences, is independent of any rotational flux components, and is inherently non-local and smooth. A one-shot inversion procedure is applied, minimising the mis-match between parameterised force functions and force functions derived from eddy resolving calculations. This enables diffusivities associated with the eddy potential vorticity and buoyancy fluxes to be diagnosed. The methodology is applied to multi-layer quasi-geostrophic ocean gyre simulations. It is found that: (i) a strictly down-gradient mixing scheme has limited success in reducing the mis-match compared to one with no sign constraint on the diffusivity; (ii) negative signals of diffusivities are prevalent around the time-mean jet; (iii) there is some indication that the magnitude of the diffusivity correlates well the eddy energy. Implications for parameterisation are discussed in light of these diagnostic results.