engleski

Generalizing "z-less" mixing length for stable boundary layers

Current research shows, recent and former suggests, that the nature and evolution of the stable atmospheric boundary layer(s) (SABL) is still understood and modelled inadequately. The "classical" SABL, almost always stratified weakly (i.e. gradient Richardson number, Ri << infinity, typically, 0 < Ri <= 1), has been modelled reasonably well during the last few decades or so, but very stable case, i.e. the VSABL (Ri >> 1), is generally not well understood. Excessively diffusive and much too deep VSABL flows, as often appearing in numerical models, were recently addressed ; the over-diffusion was alleviated by improving the local turbulent mixing length. This demands an explicit inclusion of the vertical shear of horizontal wind, S, in the mixing length, besides the previously known role of buoyancy frequency, N. A generalization of this recent work is given here by a simplified turbulent kinetic energy (TKE) equation and a set of subsequent parameterizations for the eddy diffusivity and conductivity, i.e. K- parameterizations, in terms of a generalized "z-less" mixing length, LAMBDA. The aim is to produce a parameterization that is uniformly valid for all Ri >= 0. It is shown that LAMBDA ~ (TKE)^1/2 / |S| f(Ri, Pr), uniformly valid for 0 <= Ri, regardless of the other parameterization details (the details appear as corrections) ; Pr is turbulent Prandtl number and f(Ri, Pr) is a simple set of derived functions depending on the parameterization properties. A couple of important shortcomings of the current turbulence parameterizations for the SABL, as modelled in numerical weather prediction, air-chemistry and climate models, will be remedied by using this new generalized "z-less" mixing length. This approach also recommends that it should be better if various mixing length scales were derived from simplified main principles, instead of only being guessed from plausible reasoning or dimensional analysis. In particular, it has often been assumed that K and LAMBDA profiles could be chosen for modelling purposes more or less independently from each other. It is shown, based on a simple renormalization for LAMBDA, that this is not the case if one wishes to use a more consistent parameterization scheme.

diffusion; parameterization; Prandtl number; Richardson gradient number; stably-stratified turbulence; very stable boundary layers

Rad je potpomognut i iz sredstava međunarodnog projekta EMEP4HR, 175183/S30 (Research Council of Norway).