On the bulk Richardson number and flux-profile relations in an atmospheric surface layer under weak wind stable conditions

For describing the surface fluxes in the atmospheric dispersion models, the commonly used linear universal similarity functions phi(m) and phi(h) for nondimensional wind and temperature profiles (J. Atmos. Sci. 28 (1971) 181; Boundary-Layer Meteorol. 7 (1974) 363) are analyzed under weak wind stable conditions. The stability parameter z/L (Z is the height above the ground and L is the Obukhov length) is not known a priori. On the other hand, the bulk Richardson number (Ri(B)) is readily available from profile measurements. Thus, a systematic mathematical analysis is carried out to analyze the applicability of the linear functions in terms of Ri(B). These linear functions are found to be valid as long as RiB is smaller than Pr(t)gamma/beta(2) where Pr-t is the turbulent Prandtl number and beta and gamma are constants appearing in the linear functions of phi(m) and phi(h). The linear functional forms for phi(m) and phi(h) underestimate the Obukhov length resulting in the large value of the stability parameter when Ri(B)greater than or equal toPr(t) gamma/beta(2) occurring in weak wind stable conditions. Using the data from plume validation field experiment conducted by Electric Power Research Institute at Kincaid, it is found that in 70% of the weak wind cases, the Ri(B) is larger than Pr-t gamma/beta(2). The similarity functions proposed by Beljaars and Holtslag (J. Appl. Meteorol. 30 (1991) 327) are shown to perform well in weak wind stable conditions for all the values of Ri(B). The computed surface fluxes are also compared with those based on turbulence measurements taken in Cooperative Atmospheric Surface Exchange Study. (C) 2003 Elsevier Ltd. All rights reserved.