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.

Earth's climate sensitivity depends on how shallow clouds in the trades respond to changes in the large-scale tropical circulation with warming. In canonical theory for this cloud-circulation coupling, it is assumed that the clouds are controlled by the field of vertical motion on horizontal scales larger than the convection's depth (~1 km). This assumption has been challenged both by recent in-situ observations, and idealized large-eddy simulations (LESs). Here, we therefore bring together the recent observations, new analysis from satellite data, and a forty-day, large-domain (1600 x 900 km2) LES of the North Atlantic from the 2020 EUREC4A field campaign, to study the interaction between shallow convection and vertical motions on scales between 10-1000 km (mesoscales), in settings that are as realistic as possible. Across all datasets, the shallow mesoscale vertical motions are consistently represented, ubiquitous, frequently organised into circulations, and formed without imprinting themselves on the mesoscale buoyancy field. Therefore, we use the weak-temperature gradient approximation to show that between at least 12.5-400 km scales, the vertical motion balances heating fluctuations in groups of precipitating shallow cumuli. That is, across the mesoscales, shallow convection controls the vertical motion in the trades, and does not simply adjust to it. In turn, the mesoscale convective heating patterns appear to consistently grow through moisture-convection feedback. Therefore, to represent and understand the cloud-circulation coupling of trade cumuli, the full range of scales between the synoptics and the hectometre must be included in our conceptual and numerical models.

Recent observations of the trade-wind regions highlight the covariability between cold-pool properties and mesoscale cloud organization. Given the covariability of organization with cloud cover and albedo, this suggests a potential impact of cold pools on the cloud radiative effect (CRE). To explore this, we use an ensemble of 103 large-domain, high-resolution, large-eddy simulations and investigate how the variability in cold pools is determined by large-scale external cloud-controlling factors (CCFs) and shaped by processes within the mesoscale. It is demonstrated that the size and frequency of occurrence of cold pools are strongly influenced by the near-surface horizontal wind speed and large-scale subsidence. The temporal evolution of cold pools is strongly correlated with the diurnality in radiation. Even without external variability, we find a strong intermittent behaviour in the evolution of cold pools, governed by a complex interplay between cold pools and clouds which expresses itself in the form of shallow squall lines. These squall lines result from precipitating downdrafts, cold pool outflows and the resulting gust fronts, reinforcing parent clouds. Cold pools influence the CRE of trade cumuli, but only when they exist during the day. This emphasizes the importance of the synchronization between cold-pool events and the diurnal cycle of insolation for the dependence of the CRE on cold pools.

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.

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.