In this paper we advance physical background of the energy- and flux-budget turbulence closure based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time-scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth rotation. In accordance to modern experimental evidence, the closure implies maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: "strong turbulence" at Ri << 1 typical of boundary-layer flows and characterised by the practically constant turbulent Prandtl number; and "weak turbulence" at Ri > 1 typical of the free atmosphere or deep ocean, where the turbulent Prandtl number asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy-viscosity and eddy-conductivity, and general non-gradient model based on prognostic equations for all basic parameters of turbulence including turbulent fluxes.

Traditionally, turbulence energetics is characterized by turbulent kinetic energy (TKE) and modelled using solely the TKE budget equation. In stable stratification, TKE is generated by the velocity shear and expended through viscous dissipation and work against buoyancy forces. The effect of stratification is characterized by the ratio of the buoyancy gradient to squared shear, called Richardson number, Ri. It is widely believed that at Ri exceeding a critical value, Ric, local shear cannot maintain turbulence, and the flow becomes laminar. We revise this concept by extending the energy analysis to turbulent potential and total energies (TPE and TTE = TKE + TPE), consider their budget equations, and conclude that TTE is a conservative parameter maintained by shear in any stratification. Hence there is no "energetics Ric", in contrast to the hydrodynamic-instability threshold, Ric-instability, whose typical values vary from 0.25 to 1. We demonstrate that this interval, 0.25<Ri<1, separates two different turbulent regimes: strong mixing and weak mixing rather than the turbulent and the laminar regimes, as the classical concept states. This explains persistent occurrence of turbulence in the free atmosphere and deep ocean at Ri>>1, clarify principal difference between turbulent boundary layers and free flows, and provide basis for improving operational turbulence closure models.