A New Look at Richardson Number Mixing Schemes for Equatorial Ocean Modeling

Abstract

A reexamination of turbulence dissipation measurements from the equatorial Pacific shows that the turbulence diffusivities are not a simple function of the gradient Richardson number. A widely used mixing scheme, the K-profile parameterization, overpredicts the turbulent vertical heat flux by roughly a factor of 4 in the stably stratified region between the surface mixed layer and the Equatorial Undercurrent (EUC). Additionally, the heat flux divergence is of the incorrect sign in the upper 80 m. An alternative class of parameterizations is examined that expresses the mixing coefficients in terms of the large-scale kinetic energy, shear, and Richardson number. These representations collapse the turbulence diffusivities above and below the Equatorial Undercurrent, and a tuned version is able to reproduce the vertical turbulence heat flux within the 50–180-m depth range. Kinetic energy is not Galilean invariant, so the collapse of the data with the new parameterization suggests that oceanic turbulence responds to boundary forcing at depths well below the surface mixed layer.