Numerical simulations of the tropical mesoscales often exhibit a self-reinforcing feedback between cumulus convection and shallow circulations, which leads to the self-aggregation of clouds into large clusters. We investigate whether this basic feedback can be adequately captured by large-eddy simulations (LESs). To do so, we simulate the non-precipitating, cumulus-topped boundary layer of the canonical “BOMEX” case over a range of numerical settings in two models. Since the energetic convective scales underpinning the self-aggregation are only slightly larger than typical LES grid spacings, aggregation timescales do not converge even at rather high resolutions (<100 m). Therefore, high resolutions or improved sub-filter scale models may be required to faithfully represent certain forms of trade-wind mesoscale cloud patterns and self-aggregating deep convection in large-eddy and cloud-resolving models, and to understand their significance relative to other processes that organize the tropical mesoscales.@en

Condensation in cumulus clouds plays a key role in structuring the mean, nonprecipitating trade wind boundary layer. Here, we summarize how this role also explains the spontaneous growth of mesoscale [.O(10) km] fluctuations in clouds and moisture around the mean state in a minimal-physics, large-eddy simulation of the undisturbed period during BOMEX on a large [O(100) km] domain. Small, spatial anomalies in condensation in cumulus clouds, which form on top of small moisture fluctuations, power circulations that transport moisture, but not heat, from dry to moist regions, and thus reinforce the condensation anomaly. We frame this positive feedback as a linear instability in mesoscale moisture fluctuations, whose time scale depends only on (i) a vertical velocity scale and (ii) the mean environment's vertical structure. In our minimal-physics setting, we show both ingredients are provided by the shallow cumulus convection itself: it is intrinsically unstable to length scale growth. The upshot is that energy released by clouds at kilometer scales may play a more profound and direct role in shaping the mesoscale trade wind environment than is generally appreciated, motivating further research into the mechanism's relevance. @en

Assessment of the impact of climate change on water resources over land requires knowledge on the origin of the precipitation and changes therein toward the future. We determine the origin of precipitation over the Mississippi River basin (MRB) using high-resolution (~25 km) climate model simulations for present and future climate (RCP4.5). Moisture resulting in precipitation over the MRB is tracked back in time using Eulerian offline moisture tracking, in order to find out from where this water originally evaporated (i.e., the moisture sources). We find that the most important continental moisture sources are the MRB itself and the area southwest of the basin. The two most relevant oceanic sources are the Gulf of Mexico/Caribbean and the Pacific. The distribution of sources varies per season, with more recycling of moisture within the basin during summer and more transport of moisture from the ocean toward the basin in winter. In future winters, we find an increase in moisture source from the oceans (related to higher sea surface temperatures), resulting in more precipitation over the MRB. In future summers, we find an approximately 5% decrease in moisture source from the basin itself, while the decrease in precipitation is smaller (i.e., lower recycling ratios). The results here are based on one climate model, and we do not study low-frequency climate variability. We conclude that Mis-sissippi’s moisture sources will become less local in a future climate, with more water originating from the oceans.@en

Observations of two typical contrasting weakly stable and very stable boundary layers from the winter at Dome C station, Antarctica, are used as a benchmark for two centimetre-scale-resolution large-eddy simulations. By taking the Antarctic winter, the effects of the diurnal cycle are eliminated, enabling the study of the long-lived steady stable boundary layer. With its homogeneous, flat snow surface, and extreme stabilities, the location is a natural laboratory for studies on the long-lived stable boundary layer. The two simulations differ only in the imposed geostrophic wind speed, which is identified as the main deciding factor for the resulting regime. In general, a good correspondence is found between the observed and simulated profiles of mean wind speed and temperature. Discrepancies in the temperature profiles are likely due to the exclusion of radiative transfer in the current simulations. The extreme stabilities result in a considerable contrast between the stable boundary layer at the Dome C site and that found at typical mid-latitudes. The boundary-layer height is found to range from approximately 50m to just 5m in the most extreme case. Remarkably, heating of the boundary layer by subsidence may result in thermal equilibrium of the boundary layer in which the associated heating is balanced by the turbulent cooling towards the surface. Using centimetre-scale resolutions, accurate large-eddy simulations of the extreme stabilities encountered in Antarctica appear to be possible. However, future simulations should aim to include radiative transfer and sub-surface heat transport to increase the degree of realism of these types of simulations.@en

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.@en

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.@en

The representation of snow processes in most large-scale hydrological and climate models is known to introduce considerable uncertainty into the predictions and projections of water availability. During the critical snowmelt period, the main challenge in snow modeling is that net radiation is spatially highly variable for a patchy snow cover, resulting in large horizontal differences in temperatures and heat fluxes. When a wind blows over such a system, these differences can drive advection of sensible and latent heat from the snow-free areas to the snow patches, potentially enhancing the melt rates at the leading edge and increasing the variability of subgrid melt rates. To get more insight into these processes, we examine the melt along the upwind and downwind edges of a 50 m long snow patch in the Finseelvi catchment, Norway, and try to explain the observed behavior with idealized simulations of heat fluxes and air movement over patchy snow. The melt of the snow patch was monitored from 11 June until 15 June 2019 by making use of height maps obtained through the photogrammetric structure-from-motion principle. A vertical melt of 23 ± 2.0 cm was observed at the upwind edge over the course of the field campaign, whereas the downwind edge melted only 3 ± 0.4 cm. When comparing this with meteorological measurements, we estimate the turbulent heat fluxes to be responsible for 60 % to 80 % of the upwind melt, of which a significant part is caused by the latent heat flux. The melt at the downwind edge approximately matches the melt occurring due to net radiation. To better understand the dominant processes, we represented this behavior in idealized direct numerical simulations, which are based on the measurements on a single snow patch by and resemble a flat, patchy snow cover with typical snow patch sizes of 15, 30, and 60 m. Using these simulations, we found that the reduction of the vertical temperature gradient over the snow patch was the main cause of the reductions in vertical sensible heat flux over distance from the leading edge, independent of the typical snow patch size. Moreover, we observed that the sensible heat fluxes at the leading edge and the decay over distance were independent of snow patch size as well, which resulted in a 15 % and 25 % reduction in average snowmelt for, respectively, a doubling and quadrupling of the typical snow patch size. These findings lay out pathways to include the effect of highly variable turbulent heat fluxes based on the typical snow patch size in large-scale hydrological and climate models to improve snowmelt modeling.@en

High-resolution large-eddy simulations of the Antarctic very stable boundary layer reveal a mechanism for systematic and periodic intermittent bursting. A nonbursting state with a boundary layer height of just 3 m is alternated by a bursting state with a height of ≈5 m. The bursts result from unstable wave growth triggered by a shear-generated Kelvin–Helmholtz instability, as confirmed by linear stability analysis. The shear at the top of the boundary layer is built up by two processes. The upper, quasi-laminar layer accelerates due to the combined effect of the pressure force and rotation by the Coriolis force, while the lower layer decelerates by turbulent friction. During the burst, this shear is eroded and the initial cause of the instability is removed. Subsequently, the interfacial shear builds up again, causing the entire sequence to repeat itself with a time scale of ≈10 min. Despite the clear intermittent bursting, the overall change of the mean wind profile is remarkably small during the cycle. This enables such a fast erosion and recovery of the shear. This mechanism for cyclic bursting is remarkably similar to the mechanism hypothesized by Businger in 1973, with one key difference. Whereas Businger proposes that the flow acceleration in the upper layer results from downward turbulent transfer of high-momentum flow, the current results indicate no turbulent activity in the upper layer, hence requiring another source of momentum. Finally, it would be interesting to construct a climatology of shear-generated intermittency in relation to large-scale conditions to assess the generality of this Businger mechanism.@en

Adaptive grids hold an important advantage over the usage of static and pre-tuned grids, as the most challenging atmospheric scenarios are typically localized in the spatial and temporal domain. Examples include the radiative cooling of cumulus cloud interfaces and the diurnal cycle. We present results obtained with applying an adaptivegrid approach to a variety of modeling fidelities, ranging from turbulence resolving direct numerical simulations and large-eddy simulations to Reynolds-averaged techniques. Corresponding to existing literature, the overall analysis shows that the adaptive-grid approach may play an important role in the future of atmospheric modeling as it enables the required computational effort for running a model at an increased resolution to scale correspondingly with the complexity of the resolved physics. Which is necessarily favorable compared to the scaling of fixed-grid approaches. @en