Small shallow cumulus clouds (< 1 km) over the tropical oceans appear to possess the ability to self-organise into mesoscale (10-100 km) patterns. To better understand the processes leading to such self-organized convection, we present Cloud Botany, an ensemble of 103 large-eddy simulations on domains of 150 km, produced by the Dutch Large Eddy Simulation (DALES) model on supercomputer Fugaku. Each simulation is run in an idealized, fixed, larger-scale environment, controlled by six free parameters. We vary these over characteristic ranges for the winter trades, including parameter combinations observed during the EUREC4A (Elucidating the role of clouds–circulation coupling in climate) field campaign. In contrast to simulation setups striving for maximum realism, Cloud Botany provides a platform for studying idealized, and therefore more clearly interpretable causal relationships between conditions in the larger-scale environment and patterns in mesoscale, self-organized shallow convection. We find that any simulation that supports cumulus clouds eventually develops mesoscale patterns in their cloud fields. We also find a rich variety in these patterns as our control parameters change, including cold pools lined by cloudy arcs, bands of cross-wind clouds and aggregated patches, sometimes topped by thin anvils. Many of these features are similar to cloud patterns found in nature. The published data set consists of raw simulation output on full 3D grids and 2D cross-sections, as well as post-processed quantities aggregated over the vertical (2D), horizontal (1D) and all spatial dimensions (time-series). The data set is directly accessible from Python through the use of the EUREC4A intake catalog.

In atmospheric modeling, superparameterization has gained popularity as a technique to improve cloud and convection representations in large scale models by coupling them locally to cloud-resolving models. We show how the different representations of cloud water in the local and the global models in superparameterization lead to a suppression of cloud advection and ultimately to a systematic underrepresentation of the cloud amount in the large scale model. We demonstrate this phenomenon in a regional superparameterization experiment with the global model OpenIFS coupled to the local model DALES (the Dutch Atmospheric Large Eddy Simulation), as well as in an idealized setup, where the large-scale model is replaced by a simple advection scheme. To mitigate the problem of suppressed cloud advection, we propose a scheme where the spatial variability of the local model’s total water content is enhanced in order to achieve the correct cloud condensate amount.

Motivated by an observed relationship between marine low cloud cover and surface wind speed, this study investigates how vertical wind shear affects trade-wind cumulus convection, including shallow cumulus and congestus with tops below the freezing level. We ran large-eddy simulations for an idealised case of trade-wind convection using different vertical shears in the zonal wind. Backward shear, whereby surface easterlies become upper westerlies, is effective at limiting vertical cloud development, which leads to a moister, shallower and cloudier trade-wind layer. Without shear or with forward shear, shallow convection tends to deepen more, but clouds tops are still limited under forward shear. A number of mechanisms explain the observed behaviour: First, shear leads to different surface wind speeds and, in turn, surface heat and moisture fluxes due to momentum transport, whereby the weakest surface wind speeds develop under backward shear. Second, a forward shear profile in the subcloud layer enhances moisture aggregation and leads to larger cloud clusters, but only on large domains that generally support cloud organization. Third, any absolute amount of shear across the cloud layer limits updraft speeds by enhancing the downward-oriented pressure perturbation force. Backward shear — the most typical shear found in the winter trades — can thus be argued a key ingredient at setting the typical structure of the trade-wind layer.

Momentum transport by boundary-layer turbulence causes a weak synoptic-scale vertical motion. The classical textbook solution for the strength of this Ekman pumping depends on the curl of the surface momentum flux. A new solution for Ekman pumping is derived in terms of the curl of the geostrophic wind and a term that depends in a non-trivial way on the vertical profile of the turbulent momentum flux. The solution is confined to a boundary-layer regime that is vertically well mixed and horizontally homogeneous. The momentum flux is computed from a commonly used bulk surface drag formula and a flux-jump relation to capture the entrainment flux of momentum at the top of the boundary layer. It is found that the strength of Ekman pumping is bounded. The weakening of Ekman pumping for enhanced turbulent surface friction can be explained from the fact that it will reduce the magnitude of the horizontal wind. It is demonstrated that entrainment of momentum across the top of the boundary layer tends to diminish the large-scale divergence of the wind. As momentum transport is parameterized in large-scale models, the analysis is relevant for the understanding and interpretation of the evolution of synoptic-scale vertical motions as predicted by such models.

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 large cloud structures. 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 (less than 100m). Therefore, high resolutions or improved unresolved scales 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 organise the tropical mesoscales.

Most cumulus parametrizations today make use of a simple conceptual model of convection, called the mass-flux approach. This approach depicts convection as an ensemble of updrafts and downdrafts occurring within a model grid-box. The aim of this study is to determine convective mass-fluxes and their constituents on the scale of a 100 km GCM grid-box from a C-band polarimetric radar and thereafter investigate the relative role of area fraction and vertical velocity in determining the shape and magnitude of bulk mass-flux profiles. We make use of observational estimates of these quantities spanning 13 wet seasons in the tropical region of Darwin. Following a bulk approach, the results show that the distribution of mass-flux is positively skewed and its mean profile peaks at 4 km. This is the result of constant area fractions and increasing vertical velocities below that level. Above 4 km, in-cloud vertical velocity plays a marginal role compared to the convective area fraction in controlling mass-flux profiles.

The science guiding the \EURECA campaign and its measurements are presented. \EURECA comprised roughly five weeks of measurements in the downstream winter trades of the North Atlantic — eastward and south-eastward of Barbados. Through its ability to characterize processes operating across a wide range of scales, \EURECA marked a turning point in our ability to observationally study factors influencing clouds in the trades, how they will respond to warming, and their link to other components of the earth system, such as upper-ocean processes or, or the life-cycle of particulate matter. This characterization was made possible by thousands (2500) of sondes distributed to measure circulations on meso (200 km) and larger (500 km) scales, roughly four hundred hours of flight time by four heavily instrumented research aircraft, four global-ocean class research vessels, an advanced ground-based cloud observatory, a flotilla of autonomous or tethered measurement devices operating in the upper ocean (nearly 10000 profiles), lower atmosphere (continuous profiling), and along the air-sea interface, a network of water stable isotopologue measurements, complemented by special programmes of satellite remote sensing and modeling with a new generation of weather/climate models. In addition to providing an outline of the novel measurements and their composition into a unified and coordinated campaign, the six distinct scientific facets that \EURECA explored — from Brazil Ring Current Eddies to turbulence induced clustering of cloud droplets and its influence on warm-rain formation — are presented along with an overview \EURECA’s outreach activities, environmental impact, and guidelines for scientific practice.

This study investigates how wind shear and momentum fluxes in the surface- and boundary layer vary across wind and cloud regimes. We use a nine-year-long data set from the Cabauw tower of the Ruisdael Observatory (NL) complemented by (8.2 x 8.2 km^2) daily LES hindcasts. An automated algorithm classifies observed and simulated days into different cloud regimes: 1) clear-sky days, 2) days with convective clouds (cumulus) rooted in the surface layer, with three ranges of cloud cover, and 3) days with clouds not rooted near the surface. Categorized days in observations and LES do not fully match, with a tendency of the LES to develop convective clouds on clear-sky days and less frequently produce non-rooted clouds, whose scales are far larger than the LES domain. Even so, the climatology and diurnal cycle of winds are for each regime very similar in LES and observations, strengthening our confidence in LES’ skill to reproduce certain clouds for an atmospheric state. Wind shear is smallest in clear-sky and cumulus regimes with limited cloud cover (CLCC), which also have the weakest 200 m wind speed and largest surface buoyancy flux. They have notably larger cross-wind fluxes, although along-wind momentum flux profiles are similar across all regimes. Cloudy days have larger momentum fluxes distributed over deeper layers, sustaining up to 20% of the surface flux value at cloud base. Compared to clear-sky, the CLCC regimes have stronger updrafts and deeper mixed-layers. At similar atmospheric stability, surface friction is larger and underestimated by Monin-Obukhov Similarity Theory.

Shallow cloud fields over the subtropical ocean exhibit many spatial patterns. The frequency of occurrence of these patterns can change under global warming. Hence, they may influence subtropical marine clouds’ climate feedback. While numerous metrics have been proposed to quantify cloud patterns, a systematic, widely accepted description is still missing. Therefore, this paper suggests one. We compute 21 metrics for 5000 satellite scenes of shallow clouds over the subtropical Atlantic Ocean and translate the resulting dataset to its principal components (PCs). This yields a unimodal, continuous distribution without distinct classes, whose first four PCs explain 82% of all 21 metrics’ variance. The PCs correspond to four interpretable dimensions: Characteristic length, void size, directional alignment and horizontal cloud-top height variance. These dimensions span a space in which an effective pattern description can be given, which may be used to better understand the patterns’ underlying physics and feedback on climate.

A new generation of operational atmospheric models operating at horizontal resolutions in the range 200 m ~ 2 km is becoming increasingly popular for operational use in numerical weather prediction and climate applications. Such grid spacings are becoming sufficiently fine to resolve a fraction of the turbulent transports. Here we analyze LES results of a convective boundary layer obtained by coarsening horizontal grid spacings up to 800 m. The aim is to explore the dependency of the mean state and turbulent fluxes on the grid resolution. Both isotropic and anisotropic eddy diffusion approaches are evaluated, where in the latter case the horizontal and vertical eddy diffusivities differ in accord with their horizontal and vertical grid spacings. For coarsening horizontal grid sizes entrainment at the top of the boundary layer tends to get slightly enhanced for isotropic diffusion. An analysis of the energy spectrum shows that anisotropic diffusion causes relatively more dissipation of variance at smaller length scales. This leads, in turn, to a shift of spectral energy towards larger length scales. This can also be clearly seen from the different kinds of spatial organization. The present study therefore suggests that details with regards to the representation of processes at small scales might impact the organization at length scales much larger than the smallest scales that can be resolved by the model.