Abstract. In this paper we present the current version of the Parallelized Large-Eddy Simulation Model (PALM) whose core has been developed at the Institute of Meteorology and Climatology at Leibniz Universität Hannover (Germany). PALM is a Fortran 95-based code with some Fortran 2003 extensions and has been applied for the simulation of a variety of atmospheric and oceanic boundary layers for more than 15 years. PALM is optimized for use on massively parallel computer architectures and was recently ported to general-purpose graphics processing units. In the present paper we give a detailed description of the current version of the model and its features, such as an embedded Lagrangian cloud model and the possibility to use Cartesian topography. Moreover, we discuss recent model developments and future perspectives for LES applications.

High‐resolution large‐eddy simulations (LESs) are used to investigate the effect of turbulence as well as the interaction between atmosphere and soil on the life cycle of nocturnal radiation fog. The first part of the article focuses on the validation of the LES model system for simulating radiation fog against measurement data from the meteorological super‐site at Cabauw (Netherlands). As in former studies, differences found in the fog life cycle, depth, and liquid water content between LES and observations can be largely ascribed to the presence of local advection processes in the observational data and uncertainty in the measurement data used to initialize the model. Moreover it is found that the choice of the droplet number concentration within the cloud microphysical representation has a high impact on the liquid water content within the fog layer, but a rather small effect on its life cycle. In the second part of the article, a set of LES runs is analysed with idealized initial conditions in order to study the effect of turbulent mixing as well as the initial state of the soil on the development of radiation fog. As expected, the results show that turbulent mixing has a strong impact on the time of fog formation, which is complicated by the interaction with both radiative cooling and water vapour removal by dew deposition. Furthermore, it is found that the near‐surface soil temperature plays a key role for the exact time of fog formation, whereas near‐surface soil moisture is decisive for the lifting and dissipation time of the fog layer.

Abstract. In this paper, we describe the PALM model system 6.0. PALM (formerly an abbreviation for Parallelized Large-eddy Simulation Model and now an independent name) is a Fortran-based code and has been applied for studying a variety of atmospheric and oceanic boundary layers for about 20 years. The model is optimized for use on massively parallel computer architectures. This is a follow-up paper to the PALM 4.0 model description in Maronga et al. (2015). During the last years, PALM has been significantly improved and now offers a variety of new components. In particular, much effort was made to enhance the model with components needed for applications in urban environments, like fully interactive land surface and radiation schemes, chemistry, and an indoor model. This paper serves as an overview paper of the PALM 6.0 model system and we describe its current model core. The individual components for urban applications, case studies, validation runs, and issues with suitable input data are presented and discussed in a series of companion papers in this special issue.

We investigate the impact of observed surface heterogeneities during the LITFASS-2003 experiment on the convective boundary layer (CBL). Large-eddy simulations (LES), driven by observed near-surface sensible and latent heat fluxes, were performed for the diurnal cycle and compare well with observations. As in former studies of idealized one-and two-dimensional heterogeneities, secondary circulations developed that are superimposed on the turbulent field and that partly take over the vertical transport of heat and moisture. The secondary circulation patterns vary between local and roll-like structures, depending on the background wind conditions. For higher background wind speeds, the flow feels an effective surface heat-flux pattern that derives from the original pattern by streamwise averaging. This effective pattern generates a roll-like secondary circulation with roll axes along the mean boundary-layer wind direction. Mainly the upstream surface conditions control the secondary circulation pattern, where the fetch increases with increasing background wind speed. Unlike the entrainment flux that appears to be slightly decreased compared to the homogeneously-heated CBL, the vertical flux of sensible heat appears not to be modified in the mixed layer, while the vertical flux of latent heat shows different responses to secondary circulations. The study illustrates that sufficient time averaging and ensemble averaging is required to separate the heterogeneity-induced signals from the raw LES turbulence data. This might be an important reason why experiments over heterogeneous terrain in the past did not give any clear evidence of heterogeneity-induced effects.

Abstract. Urban areas are an important part of the climate system and many aspects of urban climate have direct effects on human health and living conditions. This implies that reliable tools for local urban climate studies supporting sustainable urban planning are needed. However, a realistic implementation of urban canopy processes still poses a serious challenge for weather and climate modelling for the current generation of numerical models. To address this demand, a new urban surface model (USM), describing the surface energy processes for urban environments, was developed and integrated as a module into the PALM large-eddy simulation model. The development of the presented first version of the USM originated from modelling the urban heat island during summer heat wave episodes and thus implements primarily processes important in such conditions. The USM contains a multi-reflection radiation model for shortwave and longwave radiation with an integrated model of absorption of radiation by resolved plant canopy (i.e. trees, shrubs). Furthermore, it consists of an energy balance solver for horizontal and vertical impervious surfaces, and thermal diffusion in ground, wall, and roof materials, and it includes a simple model for the consideration of anthropogenic heat sources. The USM was parallelized using the standard Message Passing Interface and performance testing demonstrates that the computational costs of the USM are reasonable on typical clusters for the tested configurations. The module was fully integrated into PALM and is available via its online repository under the GNU General Public License (GPL). The USM was tested on a summer heat-wave episode for a selected Prague crossroads. The general representation of the urban boundary layer and patterns of surface temperatures of various surface types (walls, pavement) are in good agreement with in situ observations made in Prague. Additional simulations were performed in order to assess the sensitivity of the results to uncertainties in the material parameters, the domain size, and the general effect of the USM itself. The first version of the USM is limited to the processes most relevant to the study of summer heat waves and serves as a basis for ongoing development which will address additional processes of the urban environment and lead to improvements to extend the utilization of the USM to other environments and conditions.

Abstract:The aim of the research project "Innovative Strategies for Observations in the Arctic Atmospheric Boundary Layer (ISOBAR)" is to substantially increase the understanding of the stable atmospheric boundary layer (SBL) through a combination of well-established and innovative observation methods as well as by models of different complexity. During three weeks in February 2017, a first field campaign was carried out over the sea ice of the Bothnian Bay in the vicinity of the Finnish island of Hailuoto. Observations were based on ground-based eddy-covariance (EC), automatic weather stations (AWS) and remote-sensing instrumentation as well as more than 150 flight missions by several different Unmanned Aerial Vehicles (UAVs) during mostly stable and very stable boundary layer conditions. The structure of the atmospheric boundary layer (ABL) and above could be resolved at a very high vertical resolution, especially close to the ground, by combining surface-based measurements with UAV observations, i.e., multicopter and fixed-wing profiles up to 200 m agl and 1800 m agl, respectively. Repeated multicopter profiles provided detailed information on the evolution of the SBL, in addition to the continuous SODAR and LIDAR windAtmosphere 2018, 9, 268; doi:10.3390/atmos9070268www.mdpi.com/journal/atmosphereAtmosphere 2018, 9, 268 2 of 29 measurements. The paper describes the campaign and the potential of the collected data set for future SBL research and focuses on both the UAV operations and the benefits of complementing established measurement methods by UAV measurements to enable SBL observations at an unprecedented spatial and temporal resolution.

We used a set of large-eddy simulations to investigate the effect of one-dimensional stripe-like surface heat-flux heterogeneities on mixed-layer top entrainment. The profiles of sensible heat flux and the temporal evolution of the boundary-layer depth revealed decreased entrainment for small heat-flux amplitudes and increased entrainment for large heat-flux amplitudes, compared to the homogeneously-heated mixed layer. For large heat-flux amplitudes the largest entrainment was observed for patch sizes in the order of the boundary-layer depth, while for significantly smaller or larger patch sizes entrainment was similar as in the homogeneous case. In order to understand the underlying physics of this impact, a new approach was developed to infer local information on entrainment by means of the local flux divergence. We found an entrainment maximum over the centre of the stronger heated surface patch, where thermal energy is accumulated by the secondary circulation (SC) that was induced by the surface heterogeneity. Furthermore, we observed an entrainment maximum over the less heated patch as well, which we suppose is to be linked to the SC-induced horizontal flow convergence at the top of the convective boundary layer (CBL). For small heat-flux amplitudes a counteracting effect dominates that decreases entrainment, which we suppose is the horizontal advection of cold air in the lower, and warm air in the upper, CBL by the SC, stabilizing the CBL and thus weakening thermal convection. Moreover, we found that a mean wind can reduce the heterogeneity-induced impact on entrainment. If the flow is aligned perpendicular to the border between the differentially-heated patches, the SC and thus its impact on entrainment vanishes due to increased horizontal mixing, even for moderate wind speeds. However, if the flow is directed parallel to the border between the differentially-heated patches, the SC and thus its impact on entrainment persists.

No abstract