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 present results from laboratory and computational experiments on the

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. In this paper we describe the PALM model system 6.0. PALM 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 by 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.

Scalar dispersion from ground-level sources in arrays of buildings is investigated using wind-tunnel measurements and large-eddy simulation (LES). An array of uniformheight buildings of equal dimensions and an array with an additional single tall building (wind tunnel) or a periodically repeated tall building (LES) are considered. The buildings in the array are aligned and form long streets. The sensitivity of the dispersion pattern to small changes in wind direction is demonstrated. Vertical scalar fluxes are decomposed into the advective and turbulent parts and the influences of wind direction and of the presence of the tall building on the scalar flux components are evaluated. In the uniform-height array turbulent scalar fluxes are dominant, whereas the tall building produces an increase of the magnitude of advective scalar fluxes that yields the largest component. The presence of the tall building causes either an increase or a decrease to the total vertical scalar flux depending on the position of the source with respect to the tall building. The results of the simulations can be used to develop parametrizations for street-canyon dispersion models and enhance their capabilities in areas with tall buildings.

Abstract. The PALM 6.0 model system has been rapidly developed in the recent years with respect to its capability to simulate physical processes within urban environments. In this regard, it includes e.g. energy-balance solvers for building and land surfaces, a radiative transfer model to account for multiple reflections and shading, as well as a plant-canopy model to consider the effects of plants on the (thermo)dynamics of the flow. This study provides a thorough evaluation of modelled meteorological, air chemistry and wall-surface quantities against dedicated in-situ measurements taken in an urban environment in Prague, Dejvice, Czech Republic. Measurements included e.g. monitoring of air quality and meteorology in street canyons, surface temperature scanning with infrared camera and monitoring of wall heat fluxes. Large-eddy simulations (LES) for multiple days within two summer and three winter episodes that are characterized by different atmospheric conditions were performed with the PALM model driven by boundary conditions obtained from a mesoscale model. For the simulated episodes, the resulting temperature, wind speed and concentrations of chemical compounds within street canyons agreed well with the observations, except the LES did not adequately capture nighttime cooling near the surface at certain meteorological conditions. In some situations, less turbulent mixing was modelled resulting in higher near-surface concentrations. At most of the surface evaluation points the simulated wall-surface temperature agreed fairly well with the observed one regarding its absolute value as well as daily amplitude. However, especially for the winter episodes and for modern buildings with multi-layer walls, the heat transfer through the wall is partly not well captured leading to discrepancies between the modelled and observed wall-surface temperature. Furthermore, we show that model results depend on the accuracy of the input data, particularly the temperatures of surfaces affected by nearby trees strongly depend on the spatial distribution of the leaf area density, land-surface temperatures at grass surfaces strongly depend on the initial soil moisture, or wall-surface temperatures depend on the correct prescription of wall material parameters, though these parameters are often not available with sufficient accuracy. Moreover, we also point out current model limitations, here we particularly focus on implications with respect to the discrete representation of topography on a Cartesian grid, complex heterogeneous facades, as well as glass facades that are not well represented in terms of radiative processes. With these findings presented, we aim to validate the representation of physical processes in PALM as well as to point out specific shortcomings. This will help to build a baseline for future developments of the model and for improvements of simulations of physical processes in an urban environment.

One-dimensional urban models embedded in mesoscale numerical models may place several grid points within the urban canopy. This requires an accurate parametrization for shear stresses (i.e. vertical momentum fluxes) including the dispersive stress and momentum sinks at these points. We used a case study with a packing density of 33% and checked rigorously the vertical variation of spatially-averaged total shear stress, which can be used in a one-dimensional column urban model. We found that the intrinsic spatial average, in which the volume or area of the solid parts are not included in the average process, yield greater time-spatial average of total stress within the canopy and a more evident abrupt change at the top of the buildings than the comprehensive spatial average, in which the volume or area of the solid parts are included in the average.

Abstract. The Radiative Transfer Model (RTM) is an explicitly resolved three-dimensional multi-reflection radiation model integrated into the PALM modelling system. It is responsible for modelling complex radiative interactions within the urban canopy. It represents a key component in modelling energy transfer inside the urban layer and consequently PALM's ability to provide explicit simulations of the urban canopy at metre-scale resolution. This paper presents RTM version 3.0, which is integrated into the PALM modelling system version 6.0. This version of RTM has been substantially improved over previous versions. A more realistic representation is enabled by the newly simulated processes, e.g. the interaction of longwave radiation with the plant canopy, evapotranspiration and latent heat flux, calculation of mean radiant temperature, and bidirectional interaction with the radiation forcing model. The new version also features novel discretization schemes and algorithms, namely the angular discretization and the azimuthal ray tracing, which offer significantly improved scalability and computational efficiency, enabling larger parallel simulations. It has been successfully tested on a realistic urban scenario with a horizontal size of over 6 million grid points using 8192 parallel processes.