The representation of the atmospheric boundary layer is an important part of weather and climate models and impacts many applications such as air quality and wind energy. Over the years, the performance in modeling 2-m temperature and 10-m wind speed has improved but errors are still significant. This is in particular the case under clear skies and low wind speed conditions at night as well as during winter in stably stratified conditions over land and ice. In this paper, the authors review these issues and provide an overview of the current understanding and model performance. Results from weather forecast and climate models are used to illustrate the state of the art as well as findings and recommendations from three intercomparison studies held within the Global Energy and Water Exchanges (GEWEX) Atmospheric Boundary Layer Study (GABLS). Within GABLS, the focus has been on the examination of the representation of the stable boundary layer and the diurnal cycle over land in clear-sky conditions. For this purpose, single-column versions of weather and climate models have been compared with observations, research models, and large-eddy simulations. The intercomparison cases are based on observations taken in the Arctic, Kansas, and Cabauw in the Netherlands. From these studies, we find that even for the noncloudy boundary layer important parameterization challenges remain.

The collapse of turbulence in the nocturnal boundary layer is studied by means of a simple bulk model that describes the basic physical interactions in the surface energy balance. It is shown that for a given mechanical forcing, the amount of turbulent heat that can be transported downward is limited to a certain maximum. In the case of weak winds and clear skies, this maximum can be significantly smaller than the net radiative loss minus soil heat transport. In the case when the surface has low heat capacity, this imbalance generates rapid surface cooling that further suppresses the turbulent heat transport, so that eventually turbulence largely ceases (positive feedback mechanism). The model predicts the minimum wind speed for sustainable turbulence for the so-called crossing level. At this level, some decameters above the surface, the wind is relatively stationary compared to lower and higher levels. The critical speed is predicted in the range of about 5–7 m s−1, depending on radiative forcing and surface properties, and is in agreement with observations at Cabauw. The critical value appears not very sensitive to model details or to the exact values of the input parameters. Finally, results are interpreted in terms of external forcings, such as geostrophic wind. As it is generally larger than the speed at crossing height, a 5 m s−1 geostrophic wind may be considered as the typical limit below which sustainable, continuous turbulence under clear-sky conditions is unlikely to exist. Below this threshold emergence of the very stable nocturnal boundary layer is anticipated.

A new tuning-free subgrid-scale model, termed locally averaged scale-dependent dynamic (LASDD) model, is developed and implemented in large-eddy simulations (LES) of stable boundary layers. The new model dynamically computes the Smagorinsky coefficient and the subgrid-scale Prandtl number based on the local dynamics of the resolved velocity and temperature fields. Overall, the agreement between the statistics of the LES-generated turbulence and some well-established empirical formulations and theoretical predictions (e.g., the local scaling hypothesis) is remarkable. Moreover, the simulated statistics obtained with the LASDD model show relatively little resolution dependence for the range of grid sizes considered here. In essence, it is shown here that the new LASDD model is a robust subgrid-scale parameterization for reliable, tuning-free simulations of stable boundary layers, even with relatively coarse resolutions.