The influence of poorly resolving mixing processes in the nocturnal boundary layer (NBL) on the development of the convective boundary layer the following day is studied using large-eddy simulation (LES). Guided by measurement data from meteorological sites in Cabauw (Netherlands) and Hamburg (Germany), the typical summertime NBL conditions for Western Europe are characterized, and used to design idealized (absence of moisture and large-scale forcings) numerical experiments of the diel cycle. Using the UCLA-LES code with a traditional Smagorinsky-Lilly subgrid model and a simplified land-surface scheme, a sensitivity study to grid spacing is performed. At horizontal grid spacings ranging from 3.125 m in which we are capable of resolving most turbulence in the cases of interest to grid a spacing of 100 m which is clearly insufficient to resolve the NBL, the ability of LES to represent the NBL and the influence of NBL biases on the subsequent daytime development of the convective boundary layer are examined. Although the low-resolution experiments produce substantial biases in the NBL, the influence on daytime convection is shown to be small, with biases in the afternoon boundary layer depth and temperature of approximately 100 m and 0.5 K, which partially cancel each other in terms of the mixed-layer top relative humidity.

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.

Abstract. This paper describes MicroHH 1.0, a new and open source (www.microhh.org) computational fluid dynamics code for the simulation of turbulent flows in the atmosphere. It is primarily made for direct numerical simulation, but also supports large-eddy simulation (LES). The paper covers the description of the governing equations, their numerical implementation, and the parametrizations included in the code. Furthermore, the paper presents the validation of the dynamical core in the form of convergence and conservation tests, and comparison of simulations of channel flows and slope flows against well-established test cases. The full numerical model, including the associated parametrizations for LES, has been tested for a set of cases under stable and unstable conditions, under the Boussinesq and anelastic approximation, and with dry and moist convection under stationary and time-varying boundary conditions. The paper presents performance tests showing good scaling from 256 to 32,768 processes. The Graphical Processing Unit-enabled version of the code reaches speedups of more than an order of magnitude with respect to the conventional code for a variety of cases.

Abstract. We study the interactions between atmospheric boundary layer (ABL) dynamics and atmospheric chemistry using a mixed-layer model coupled to chemical reaction schemes. Guided by both atmospheric and chemical measurements obtained during the DOMINO (Diel Oxidant Mechanisms in relation to Nitrogen Oxides) campaign (2008), numerical experiments are performed to study the role of ABL dynamics and the accuracy of chemical schemes with different complexity: the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) and a reduced mechanism of this chemical system. Both schemes produce satisfactory results, indicating that the reduced scheme is capable of reproducing the O 3 -NO x -VOC-HO x diurnal cycle during conditions characterized by a low NO x regime and small O 3 tendencies (less than 1 ppb per hour). By focusing on the budget equations of chemical species in the mixedlayer model, we show that for species like O 3 , NO and NO 2 , the influence of entrainment and boundary layer growth is of the same order as chemical production/loss. This indicates that an accurate representation of ABL processes is crucial in understanding the diel cycle of chemical species. By comparing the time scales of chemical reactive species with the mixing time scale of turbulence, we propose a classification based on the Damköhler number to further determine the importance of dynamics on chemistry during field campaigns. Our findings advocate an integrated approach, simultaneously solving the ABL dynamics and chemical reactions, in order to obtain a better understanding of chemical pathways and processes and the interpretation of the results obtained during measurement campaigns.