This study seeks to quantitatively and qualitatively understand how stability affects transport in the continuously turbulent stably stratified atmospheric boundary layer, based on a suite of large-eddy simulations. The test cases are based on the one adopted by the Global Energy and Water Cycle Experiment (GEWEX) Atmospheric Boundary Layer Study (GABLS) project, but with a largely expanded stability range where the gradient Richardson number (Rig) reaches up to around 1. The analysis is mainly focused on understanding the modification of turbulent structures and dynamics with increasing stability in order to improve the modeling of the stable atmospheric boundary layer in weather and climate models, a topic addressed in Part II of this work. It is found that at quasi equilibrium, an increase in stability results in stronger vertical gradients of the mean temperature, a lowered low-level jet, a decrease in vertical momentum transport, an increase in vertical buoyancy flux, and a shallower boundary layer. Analysis of coherent turbulent structures using two-point autocorrelation reveals that the autocorrelation of the streamwise velocity is horizontally anisotropic while the autocorrelation of the vertical velocity is relatively isotropic in the horizontal plane and its integral length scale decreases as stability increases. The effects of stability on the overall turbulent kinetic energy (TKE) and its budget terms are also investigated, and it is shown that the authors' large-eddy simulation results are in good agreement with previous experimental findings across varied stabilities. Finally, Nieuwstadt's local-scaling theory is reexamined and it is concluded that the height z is not a relevant scaling parameter and should be replaced by a constant length scale away from the surface, indicating that the z-less range starts lower than previously assumed.

No abstract

In polar regions, where the boundary layer is often stably stratified, atmospheric models produce large biases depending on the boundary-layer parametrizations and the parametrization of the exchange of energy at the surface. This model intercomparison focuses on the very stable stratification encountered over the Antarctic Plateau in 2009. Here, we analyze results from 10 large-eddy-simulation (LES) codes for different spatial resolutions over 24 consecutive hours, and compare them with observations acquired at the Concordia Research Station during summer. This is a challenging exercise for such simulations since they need to reproduce both the 300-m-deep convective boundary layer and the very thin stable boundary layer characterized by a strong vertical temperature gradient (10 K difference over the lowest 20 m) when the sun is low over the horizon. A large variability in surface fluxes among the different models is highlighted. The LES models correctly reproduce the convective boundary layer in terms of mean profiles and turbulent characteristics but display more spread during stable conditions, which is largely reduced by increasing the horizontal and vertical resolutions in additional simulations focusing only on the stable period. This highlights the fact that very fine resolution is needed to represent such conditions. Complementary sensitivity studies are conducted regarding the roughness length, the subgrid-scale turbulence closure as well as the resolution and domain size. While we find little dependence on the surface-flux parametrization, the results indicate a pronounced sensitivity to both the roughness length and the turbulence closure.

Terrestrial ecosystems are characterized by a wide range of canopy vegetation density, which is known to affect turbulent transport processes across the canopy-atmosphere interface. In the presence of a dense and horizontally homogeneous canopy, the canopy sublayer has been described as resembling a plane mixing layer. At the other extreme, where the canopy is essentially absent, the canopy sublayer is typically assumed to be similar to a turbulent boundary layer over a rough surface. However, it remains unclear how the canopy turbulence changes from boundary-layer-like to mixing-layer-like as the vegetation density increases. We use large-eddy simulation to study five different vegetation densities varying from an extremely sparse canopy to an extremely dense canopy. This investigation draws on the study of flow statistics as well as large-scale coherent turbulent structures within the canopy sublayer. The coherent structures are identified through the use of proper orthogonal decomposition. The results of skewness of velocity components and characteristic length scales suggest that, as the vegetation density increases, the canopy turbulence gradually undergoes a transition from resembling a rough-wall boundary layer to being similar to a mixing layer. As demonstrated by others, we found that the coherent structures within the canopy sublayer consist of a strong sweep/ejection motion framed by a counter-rotating vortex pair with elliptical cross-sections. As the canopy becomes denser, these important structures are shown here to be more elevated. Vegetation density does not appear to have a significant effect on the percentage of the total turbulent kinetic energy that is represented by the coherent structures.

No abstract