Stability Functions in the Stable Surface Layer Derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) Scheme

It is desirable that a surface layer scheme in an atmospheric numerical model is consistent with an atmospheric boundary layer scheme incorporated in the same model. In this study, stability functions based on Monin–Obukhov similarity theory for momentum and heat, φm and φh, in the stable surface layer are derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. The resulting stability functions are approximated by φm = 1 + 4.8z/L and φh = 0.74 + 6.0z/L, which can be analytically integrated with respect to height z to obtain momentum and heat fluxes, where L is the Obukhov length. The fluxes thus, obtained are compared with those acquired from stability functions in four previous studies: they are nearly the same as those from two of them, and show better agreement with observational data of the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) over sea ice than those from the other two studies. Detailed comparisons of the results of the MYNN scheme with the SHEBA data suggest that significant variations of the fluxes observed during “winter” when the ice was covered with dry snow may have been caused by those of the surface roughness around the observational site. The stability functions obtained from the MYNN scheme predict that the bulk and flux Richardson numbers approach critical values of 0.26 and 0.21, respectively, in the limit of z/L → ∞. These critical values result from the Kolmogorov hypothesis applied to the turbulent dissipation in the MYNN scheme and are considered to correspond to a transition from Kolmogorov to non-Kolmogorov turbulence.