The Description of Mixing in Stratified Layers without Shear in Large-Scale Ocean Models

Large-scale geophysical flows often exhibit layers with negligible vertical shear and infinite gradient Richardson number Ri. It is well known that these layers may be regions of active mixing, even in the absence of local shear production of turbulence because, among other processes, turbulence may be supplied by vertical turbulent transport from neighboring regions. This observation is contrasted by the behavior of most turbulence parameterizations used in ocean climate modeling, predicting the collapse of mixing of mass and matter if the Richardson number exceeds a critical threshold. Here, the performance of a simple model without critical Richardson number is evaluated, taking into account the diffusion of turbulence into layers without shear production and therefore avoiding the suppression of mixing at large values of Ri. The model is based on the framework of second-moment turbulence closures, focusing on the consistent modeling of the turbulent length scale for strongly stratified turbulence. Results are compared to eddy-resolving simulations of stratified shear flows that have recently become available. The model is simple enough for inclusion in ocean climate models.