loading page

Derivation and implementation of a non-gradient term to improve the oceanic convection representation within the k-ɛ parameterization
  • Alexandre Legay,
  • Bruno Deremble,
  • Hans Burchard
Alexandre Legay
Institut des Géosciences de l'Environnement

Corresponding Author:[email protected]

Author Profile
Bruno Deremble
Institut des Géosciences de l'Environnement
Author Profile
Hans Burchard
Leibniz Institute for Baltic Sea Research Warnemünde
Author Profile


The representation of turbulent fluxes during oceanic convective events is important to capture the evolution of the oceanic mixed layer. To improve the accuracy of turbulent fluxes, we examine the possibility of adding a non-gradient component in their expression in addition to the usual downgradient part. To do so, we extend the $k-\varepsilon$ algebraic second-moment closure by relaxing the assumption on the equilibrium of the temperature variance $\overline{\theta’^2}$. With this additional transport equation for the temperature variance, we obtain a $k - \varepsilon - \overline{\theta’^2}$ model (the “$k \varepsilon t$” model) which includes a non-gradient term for the temperature flux. We validate this new model against Large Eddy Simulations (LES) in both wind-forced and buoyancy-driven regimes. In both cases, we find that the vertical profile of temperature is well captured by the $k \varepsilon t$ model. Particularly, for the buoyancy-driven regime, the non-gradient term increases the portion of the mixed layer that is stably stratified. This is an improvement since this portion is too small with the $k - \varepsilon$ parameterization. Finally, a comparison of the non-gradient term with the KPP non-local term gives insights for refining the KPP’s ad hoc shape polynomial.
16 Feb 2024Submitted to ESS Open Archive
16 Feb 2024Published in ESS Open Archive