Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff

Two recently proposed turbulence closure schemes are compared against the conventional Mellor-Yamada (MY) model for stably stratified atmospheric flows. The Energy and Flux-Budget (EFB) approach solves the budgets of turbulent momentum and heat fluxes and turbulent kinetic and potential energies. The Cospectral Budget (CSB) approach is formulated in wavenumber space and integrated across all turbulent scales to obtain flow variables in physical space. Unlike the MY model, which is subject to a critical gradient Richardson number, both EFB and CSB models allow turbulence to exist at any gradient Richardson number R-t and predict a saturation of flux Richardson number (R-f -> R-fm) at sufficiently large R-i. The CSB approach further predicts the value of Rim and reveals a unique expression linking the Rotta and von Karman constants. Hence, all constants in the CSB model are nontunable and stability independent. All models agree that the dimensionless sensible heat flux decays with increasing R-i. However, the decay rate and subsequent cutoff in the MY model appear abrupt. The MY model further exhibits an abrupt cutoff in the turbulent stress normalized by vertical velocity variance, while the CSB and EFB models display increasing trends. The EFB model produces a rapid increase in the ratio of turbulent potential energy and vertical velocity variance as Rim is approached, suggesting a strong self-preservation mechanism. Vertical anisotropy in the turbulent kinetic energy is parameterized in different ways in MY and EFB, but this consideration is not required in CSB. Differences between EFB and CSB model predictions originate from how the vertical anisotropy is specified in the EFB model.