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.

Dust devils are convective vortices with a vertical axis of rotation made visible by lifted soil particles. Currently, there is great uncertainty about the extent to which dust devils contribute to the atmospheric aerosol input and thereby influence Earth's radiation budget. Past efforts to quantify the aerosol transport and study their formation, maintenance, and statistics using large‐eddy simulation (LES) have been of limited success. Therefore, some important features of dust devil‐like vortices simulated with LES still do not compare well with those of observed ones. One major difference is the simulated value of the core pressure drop, which is almost 1 order of magnitude smaller compared to the observed range of 250 to 450 Pa. However, most of the existing numerical simulations are based on highly idealized setups and coarse grid spacings. In this study, we investigate the effects of various factors on the simulated vortex strength with high‐resolution LES. For the fist time, we are able to reproduce observed core pressures by using a high spatial resolution of 2 m, a model setup with moderate background wind and a spatially heterogeneous surface heat flux. It is found that vortices mainly appear at the lines of horizontal flow convergence above the centers of the strongly heated patches, which is in contrast to some older observations in which vortices seemed to be created along the patch edges.

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.

First direct numerical simulations of dust devil-like structures for Rayleigh numbers up to 10 11 are presented• Dust devils are strongly affected by Rayleigh number and surface friction but less affected by the aspect ratio of the model domain• The 3-D structure of simulated dust devils is very similar to the one in large-eddy simulations of the atmospheric boundary layer

Dust devils are organized convective vortices with pressure drops of hundreds of pascals that spirally lift surface material into the air. This material modifies the radiation budget by contributing to the atmospheric aerosol concentration. Quantification of this contribution requires good knowledge of the dust devil statistics and dynamics. The latter can also help to understand vortex genesis, evolution and decay, in general. Dust devil-like vortices are numerically investigated mainly by large-eddy simulation (LES). A critical parameter in these simulations is the grid spacing, which has a great influence on the dust devil statistics. So far, it is unknown which grid size is sufficient to capture dust devils accurately. We investigate the convergence of simulated convective vertical vortices that resemble dust devils by using the LES model PALM. We use the nesting capabilities of PALM to explore grid spacings from 10 to 0.625 m. Grid spacings of 1 m or less have never been used for the analysis of dust devil-like vortices that develop in a horizontal domain of more than 10 km$$^2$$ 2 . Our results demonstrate that a minimum resolution of 1.25 m is necessary to achieve a convergence for sample-averaged quantities like the core pressure drop. This grid spacing or smaller should be used for future quantifications of dust devil sediment fluxes. However, sample maxima of the investigated dust devil population and peak velocity values of the general flow show no convergence. If a qualitative description of the dust devil flow pattern is sufficient, we recommend a grid spacing of 2.5 m or smaller.