Bas J. H. van de Wiel

Shear capacity as prognostic for nocturnal boundary layer regimes

Field observations and theoretical analysis are used to investigate the appearance of different nocturnal boundary layer regimes. Recent theoretical findings predict the appearance of two different regimes: the continuously turbulent (weakly stable) boundary layer and the relatively ‘‘quiet’’ (very stable) boundary layer. A large number of nights (approximately 4500 in total) are analyzed using an ensemble averaging technique. The observations support the existence of these two fundamentally different regimes: weakly stable (turbulent) nights rapidly reach a steady state (within 2‐3h). In contrast, very stable nights reach a steady state much later after a transition period (2‐6h). During this period turbulence is weak and nonstationary. To characterize the regime, a new parameter is introduced: the shear capacity. This parameter compares the actual shear after sunset with the minimum shear needed to sustain continuous turbulence. In turn, the minimum shear is dictated by the heat flux demand at the surface (net radiative cooling), so that the shear capacity combines flow information with knowledge of the boundary condition. It is shown that the shear capacity enables prediction of the flow regimes. The prognostic strength of this nondimensional parameter appears to outperform the traditional ones like the similarity parameter z/L and the gradient Richardson number Ri as a regime indicator.

Regime Transitions in Near-Surface Temperature Inversions: A Conceptual Model

AbstractA conceptual model is used in combination with observational analysis to understand regime transitions of near-surface temperature inversions at night as well as in Arctic conditions. The model combines a surface energy budget with a bulk parameterization for turbulent heat transport. Energy fluxes or feedbacks due to soil and radiative heat transfer are accounted for by a “lumped parameter closure,” which represents the “coupling strength” of the system.Observations from Cabauw, Netherlands, and Dome C, Antarctica, are analyzed. As expected, inversions are weak for strong winds, whereas large inversions are found under weak-wind conditions. However, a sharp transition is found between those regimes, as it occurs within a narrow wind range. This results in a typical S-shaped dependency. The conceptual model explains why this characteristic must be a robust feature. Differences between the Cabauw and Dome C cases are explained from differences in coupling strength (being weaker in the Antarctic). F...

A Diagnosis of Excessive Mixing in Smagorinsky Subfilter-Scale Turbulent Kinetic Energy Models

AbstractLarge-eddy simulation (LES) models are widely used to study atmospheric turbulence. The effects of small-scale motions that cannot be resolved need to be modeled by a subfilter-scale (SFS) model. The SFS contribution to the turbulent fluxes is typically significant in the surface layer. This study presents analytical solutions of the classical Smagorinsky SFS turbulent kinetic energy (TKE) model including a buoyancy flux contribution. Both a constant length scale and a stability-dependent one as proposed by Deardorff are considered. Analytical expressions for the mixing functions are derived and Monin–Obukhov similarity relations that are implicitly imposed by the SFS TKE model are diagnosed. For neutral and weakly stable conditions, observations indicate that the turbulent Prandtl number (PrT) is close to unity. However, based on observations in the convective boundary layer, a lower value for PrT is often applied in LES models. As a lower Prandtl number promotes a stronger mixing of heat, this m...

Local characteristics of the nocturnal boundary layer in response to external pressure forcing

AbstractGeostrophic wind speed data, derived from pressure observations, are used in combination with tower measurements to investigate the nocturnal stable boundary layer at Cabauw, the Netherlands. Since the geostrophic wind speed is not directly influenced by local nocturnal stability, it may be regarded as an external forcing parameter of the nocturnal stable boundary layer. This is in contrast to local parameters such as in situ wind speed, the Monin–Obukhov stability parameter (z/L), or the local Richardson number. To characterize the stable boundary layer, ensemble averages of clear-sky nights with similar geostrophic wind speeds are formed. In this manner, the mean dynamical behavior of near-surface turbulent characteristics and composite profiles of wind and temperature are systematically investigated. The classification is found to result in a gradual ordering of the diagnosed variables in terms of the geostrophic wind speed. In an ensemble sense the transition from the weakly stable to very sta...

Near-Surface Temperature Inversion Growth Rate during the Onset of the Stable Boundary Layer

AbstractThis study aims to find the typical growth rate of the temperature inversion during the onset of the stable boundary layer around sunset. The sunset transition is a very challenging period for numerical weather prediction, since neither accepted theories for the convective boundary layer, nor those for the stable boundary layer appear to be applicable. To gain more insight in this period a systematic investigation of the temperature inversion growth rate is conducted. A statistical procedure is used to analyze almost 16 years of observations from the Cabauw observational tower, supported by observations from two additional sites (Dome C and Karlsruhe). The results show that, on average, the growth rate of the temperature inversion (normalized by the maximum inversion during the night) weakly declines with increasing wind speed. The observed growth rate is quantitatively consistent among the sites, and it appears insensitive to various other parameters. The results were also insensitive to the afte...

Towards Adaptive Grids for Atmospheric Boundary-Layer Simulations

We present a proof-of-concept for the adaptive mesh refinement method applied to atmospheric boundary-layer simulations. Such a method may form an attractive alternative to static grids for studies on atmospheric flows that have a high degree of scale separation in space and/or time. Examples include the diurnal cycle and a convective boundary layer capped by a strong inversion. For such cases, large-eddy simulations using regular grids often have to rely on a subgrid-scale closure for the most challenging regions in the spatial and/or temporal domain. Here we analyze a flow configuration that describes the growth and subsequent decay of a convective boundary layer using direct numerical simulation (DNS). We validate the obtained results and benchmark the performance of the adaptive solver against two runs using fixed regular grids. It appears that the adaptive-mesh algorithm is able to coarsen and refine the grid dynamically whilst maintaining an accurate solution. In particular, during the initial growth of the convective boundary layer a high resolution is required compared to the subsequent stage of decaying turbulence. More specifically, the number of grid cells varies by two orders of magnitude over the course of the simulation. For this specific DNS case, the adaptive solver was not yet more efficient than the more traditional solver that is dedicated to these types of flows. However, the overall analysis shows that the method has a clear potential for numerical investigations of the most challenging atmospheric cases.