  We propose a new turbulence closure model based on the budget equations for
the key second moments: turbulent kinetic and potential energies: TKE and TPE
(comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent
fluxes of momentum and buoyancy (proportional to potential temperature).
Besides the concept of TTE, we take into account the non-gradient correction to
the traditional buoyancy flux formulation. The proposed model grants the
existence of turbulence at any gradient Richardson number, Ri. Instead of its
critical value separating - as usually assumed - the turbulent and the laminar
regimes, it reveals a transition interval, 0.1< Ri <1, which separates two
regimes of essentially different nature but both turbulent: strong turbulence
at Ri<<1; and weak turbulence, capable of transporting momentum but much less
efficient in transporting heat, at Ri>1. Predictions from this model are
consistent with available data from atmospheric and lab experiments, direct
numerical simulation (DNS) and large-eddy simulation (LES).
