Evaluating Models of The Neutral, Barotropic Planetary Boundary Layer using Integral Measures: Part I. Overview

Data for the cross-isobaric angle ₀, the geostrophic drag coefficient C_g, and the functions A and B of Rossby number similarity theory, obtained from meteorological field experiments, are used to evaluate a range of models of the neutral, barotropic planetary boundary layer. The data give well-defined relationships for ₀, C_g, and the integrated dissipation rate over the boundary layer, as a function of the surface Rossby number. Lettau's first-order closure mixing-length model gives an excellent fit to the data; other simple models give reasonable agreement. However more sophisticated models, e.g., higher-order closure, large-eddy simulation, direct numerical simulation and laboratory models, give poor fits to the data. The simplemodels have (at least) one free parameter in their turbulence closure that is matched toatmospheric observations; the more sophisticated models either base their closure onmore general flows or have no free closure parameters. It is suggested that all of theatmospheric experiments that we could locate violate the strict simplifying assumptionsof steady, homogeneous, neutral, barotropic flow required by the sophisticated models.The angle ₀ is more sensitive to violations of the assumptions than is C_g.The behaviour of the data varies in three latitude regimes. In middle and high latitudes the observed values of A and B exhibit little latitudinal dependence; the best estimates are A = 1.3 and B = 4.4. In lower latitudes the neutral, barotropic Rossby number theory breaks down. The value of B increases towards the Equator; the determination of A is ambiguous – the trend can increase or decrease towards the Equator. Between approximately 5° and 30° latitude, the scatter in the data is thought to be primarily due to the inherent presence of baroclinicity. The presence of the trade-wind inversion, thermal instability and the horizontal component of the Earth's rotation _H also contribute.Marked changes in the values of A and B occur in the region between the Equator andapproximately 5° latitude, as the Coriolis parameter |f| approaches zero. Although the variation of A and B with latitude suggests some similarity to the results obtained from the direct numerical simulations, the presence of additional complexities in the real atmosphere that are not included in the numerical model, precludes a meaningful direct comparison.