Articles | Volume 15, issue 9
Review and perspective paper
04 May 2022
Review and perspective paper |  | 04 May 2022

Empirical values and assumptions in the convection schemes of numerical models

Anahí Villalba-Pradas and Francisco J. Tapiador

Convection influences climate and weather events over a wide range of spatial and temporal scales. Therefore, accurate predictions of the time and location of convection and its development into severe weather are of great importance. Convection has to be parameterized in global climate models and Earth system models as the key physical processes occur at scales much lower than the model grid size. This parameterization is also used in some numerical weather prediction (NWP) models when convection is not explicitly resolved. The convection schemes described in the literature represent the physics by simplified models that require assumptions about the processes and the use of a number of parameters based on empirical values. These empirical values and assumptions are rarely discussed in the literature. The present paper examines these choices and their impacts on model outputs and emphasizes the importance of observations to improve our current understanding of the physics of convection. The focus is mainly on the empirical values and assumptions used in the activation of convection (trigger), the transport and microphysics (commonly referred to as the cloud model), and the intensity of convection (closure). Such information can assist satellite missions focused on elucidating convective processes (e.g., the INCUS mission) and the evaluation of model output uncertainties due to spatial and temporal variability of the empirical values embedded into the parameterizations.

1 Introduction

Numerical weather prediction models, global climate models, and Earth system models (NWP, GCMs, and ESMs) generate precipitation mainly through two parameterizations: microphysics of precipitation (MP hereafter) and cumulus parameterization (CP) schemes. They produce what is known as large-scale precipitation and convective precipitation, respectively. While other schemes, such as the planetary boundary layer (PBL) parameterization used to parameterize turbulence within the PBL without accounting for moist convection, also affect precipitation occurrence, the especially intricate processes by which water vapor becomes cloud droplets or ice crystals and then liquid or solid precipitation are mainly modeled by the two former modules.

The empirical values and assumptions embedded in the MP were explored in Tapiador et al. (2019b). The goal of the present paper is to provide a comprehensive account of the empirical choices and assumptions behind the representation of convective precipitation in models. There are indeed several reviews thoroughly discussing the empirical values and assumptions in convective models (e.g., de Roode et al., 2012), but they are generally focused on a particular parameter. To the best of our knowledge, there is no such extensive review of the empirical values and assumptions in the convection schemes available in the literature. Also, excellent recent reviews describing convection schemes already exist, namely Arakawa (2004) or Plant (2010), but the empiricisms in their physics have rarely been discussed. This paper aims to fill that void.

The scientific interest of our endeavor is twofold. First, it can assist dedicated satellite missions such as the Investigation of Convective Updrafts (INCUS) mission, a new Earth Venture Mission-3 (EVM-3) of three SmallSats expected to be launched in 2027 that aims to increase our knowledge of precipitation processes, specifically of the many nuances behind convection (Stephens et al., 2020). Indeed, INCUS aims to advance our present understanding and modeling of convection in the directions identified in the “decadal survey” (see Jakob, 2010; National Academies of Sciences, Engineering and Medicine, 2018, hereafter “decadal survey”). The precise description and rationale behind the empirical parameters in the parameterization of convection can help INCUS and similar missions to focus on the key parameters and to analyze their impacts on weather and climate models.

Another science goal of our review is to pinpoint the more relevant empirical values so systematic sensitivity studies can be readily carried out. We exemplify the latest goal, showing that the spread of a perturbed ensemble of just a few parameters can be substantial. Thus, we have used the European Centre for Medium-Range Forecasts (ECMWF) Integrated Forecasting System (IFS) to perform a sensitivity experiment with seven parameters (organized entrainment, entrainment for shallow convection, turbulent detrainment, adjustment time, rain conversion, momentum transport, and shallow vs. deep cloud thickness). While this is a small subset of the many parameters we have identified in this review and the experiment is intended as an illustration of the spread in the simulations for two tropical storms, the case invites more systematic runs in both space (global coverage) and time (decadal simulations) over the whole empirical set of parameters of any given model. The spread of the results will help to gauge the uncertainties due to the empiricisms embedded in the convection modules and constrain those through dedicated campaigns and targeted observations.

Precipitation is arguably the most important component of the water cycle. Extreme hydrological events in the form of floods are responsible for the loss of thousands of lives every year and great damage to property, while droughts affect water resources, livestock, and crop production. Both extremes represent important threats for human life and developing economies (e.g., Trenberth, 2011; Pham-Duc et al., 2020). Changes in the hydrological cycle also affect human activities such as the production of electricity in hydropower plants, where a better optimization of electricity production depends on water input (García-Morales and Dubus, 2007; Tapiador et al., 2011). Precipitation is also a key environmental parameter for biota. The types of vegetation and animal life that exist in a certain area are conditioned by temperature but even more by precipitation. Changes in the precipitation regime alter plant growth and survival and consequently impact the food chain (McLaughlin et al., 2002; Choat et al., 2012; Barros et al., 2014; Deguines et al., 2017). Prolonged droughts may increase the risk of wildfires, with the associated loss of local species (Holden et al., 2018). Therefore, it is not surprising that providing an accurate representation of precipitation in models is an active research topic. Specifically, in the climate realm it is already known that the effects of climate change will strongly modify the distribution and variability of precipitation around the world (Easterling et al., 2000; Dore, 2005; Giorgi and Lionello, 2008; Trenberth, 2011), posing many risks to life and human activities (Patz et al., 2005; McGranahan et al., 2007; IPCC, 2014; Woetzel et al., 2020). Thus, it is important to provide an explicit account of how models produce rain and snow in order to fully understand the outputs of the simulations.

The paper is organized as follows. A brief note on model parameterization, tuning, and the importance of convection follows (Sect. 1.1 and 1.2). Then, the main strategies to model cumulus convection are briefly presented to provide the framework to the rest of the paper (Sect. 2). The core of the review is in the following three sections, which present the assumptions and empirical values in the trigger (Sect. 3), the cloud model (Sect. 4), and the closure of the scheme (Sect. 5). The paper concludes with notes and considerations on the topic, bringing together the most important results. The acronyms used through the paper may be found in Appendix A.

1.1 Model parameterizations

Parameterizations in numerical models address the fact that some significant physical processes in nature occur at scales much lower than the grid size used in models (Arakawa and Schubert, 1974; Stensrud, 2007; McFarlane, 2011). That is the case of convection, for which spatial resolutions of at least 100 m are required to realistically solve its dynamics (Bryan et al., 2003). However, typical horizontal grid resolutions in current models range from a kilometer scale for high-resolution NWP applied to a particular area to dozens of kilometers in global NWPs, GCMs, and ESMs. With these model grids, convection is a subgrid-scale process not explicitly resolved. The physics are then represented by a simplified model that requires assumptions about the processes and the use of several parameters based on empirical values. These are used as thresholds, constraints, or mean values of a number of processes, whereas the former simplification requires a compromise between reducing complexity and a fair representation of the atmosphere.

While sometimes neglected and seldom explicit, tuning is an integral procedure of modeling (Hourdin et al., 2017; Schmidt et al., 2017; Tapiador et al., 2019a, b). It consists of estimating sensible values for the empirical parameters to reduce the discrepancies between model outputs and observations. An example of these discrepancies is shown in Figs. 1 and 2. Hence, tuning may have a significant influence on model results and can help identify the parts of the model that need further attention. However, blind tuning can mask fundamental problems within the parameterization, leading to nonrealistic physical states of the system and compensating for errors that translate into an inappropriate budget equilibrium or affect other metrics (Tapiador et al., 2019b). This is particularly important for climate models, since projections and simulations of future climates always include the ceteris paribus assumption (Smith, 2002), i.e., the tenet that in the future the multiple feedbacks between the many processes will operate in the same way as in the present.

As stated in Couvreux et al. (2021), different approaches have been proposed to avoid tuning, including the use of convection-permitting models or machine-learning approaches that replace some parameterizations by neural networks. In the former approach, the high spatial and temporal resolutions of the model allow simulating convection directly without resorting to parameterization. Couvreux et al. (2021) proposed a new method that performs a multi-case comparison between single cloud models (SCMs) and large eddy simulation (LES) to calibrate parameterizations. The method uses machine learning without replacing parameterizations due to their important role in the production of reliable climate projections. Indeed, the computing power required to perform global, centennial ensemble simulations below kilometer resolution and under several anthropogenic forcings would be enormous, so improving the parameterization of convection schemes still is a thriving research field, as described below.

Figure 1Comparison between simulated 6 h accumulated surface liquid precipitation with the new Tiedtke convection parameterization in the WRF model using GFS initial and boundary conditions (cumulus option 16 in WRF, a) and the GPM IMERG final run (b) for Typhoon Megi on 25 September 2016 from 18:00 UTC. The accumulated precipitation includes cumulus, shallow cumulus, and grid-scale rain. The domain is located over the Philippine Sea with a horizontal grid size of 10 km. Radiation scheme: RRTMG shortwave and longwave schemes, boundary layer scheme: Mellor–Yamada–Janjić scheme, microphysics scheme: NSSL two-moment scheme, land surface option: unified Noah land surface model, surface layer option: eta similarity scheme. Spinning time: 24 h. The typhoon was not seeded.


1.2 Convection: a key process in models

There is a wide range of recent research topics in convection. These topics include machine learning to parameterize moist convection (e.g., Gentine et al., 2018; O'Gorman and Dwyer, 2018; Rasp et al., 2018), stochastic parameterizations of deep convection (e.g., Buizza et al., 1999; Majda et al., 1999, 2001; Majda and Khouider, 2002; Khouider et al., 2003; Majda et al., 2003; Shutts, 2005; Plant and Craig, 2008; Dorrestijn et al., 2013a, b; Khouider, 2014; Wang et al., 2016), the use of convective parameterization on “gray zones” (e.g., Wyngaard, 2004; Kuell et al., 2007; Mironov, 2009; Gerard et al., 2009; Yano et al., 2010; Mahoney, 2016; Honnert et al., 2020), aerosols and their influence on convection (e.g., van den Heever and Cotton, 2007; Storer et al., 2010; van den Heever et al., 2011; Morrison and Grabowski, 2013; Grell and Freitas, 2014; Kawecki et al., 2016; Peng et al., 2016; Han et al., 2017; Grabowski, 2018), microphysics impacts (e.g., Grabowski, 2015), the impact of new cumulus entrainment (e.g., Chikira and Sugiyama, 2010; Lu and Ren, 2016); orographic effects on convection (e.g., Panosetti et al., 2016), new mass flux formulations (e.g., Gerard and Geleyn, 2005; Piriou et al., 2007; Guérémy, 2011; Arakawa and Wu, 2013; Park, 2014a, b; Grell and Freitas, 2014; Yano, 2014; Gerard, 2015; Kwon and Hong, 2017; Han et al., 2017), large eddy simulations (LESs) (e.g., Siebesma and Cuijpers, 1995; Gerard et al., 2002; De Rooy and Siebesma, 2008; Heus and Jonker, 2008; Neggers et al., 2009; Dawe and Austin, 2013), and scale-aware cumulus parameterization (e.g., Kuell et al., 2007; Arakawa et al., 2011; Arakawa and Wu, 2013; Grell and Freitas, 2014; Zheng et al., 2016; Kwon and Hong, 2017; Wagner et al., 2018).

Such a wealth of papers illustrates the strength of this research topic in a vast number of fields. Of these, developing parameterization schemes for models is a thriving subfield, with several teams advancing the field (see Sect. 2 below). Difficulties persist, however. Convective processes have been identified in the latest decadal survey as a major source of uncertainty, and dedicated efforts are needed to fill the gaps in our present knowledge of the processes involved. Owing to the influence of convection on climate and weather events over a large range of spatial and temporal scales, one of the most important objectives of the decadal survey is to improve the predictions of the timing and location of convective storms, as well as their evolution into severe weather. Besides the drawbacks associated with the spatial resolution, the multiscale interactions leading to the organization and evolution of convective systems are difficult to observe and represent. Improving the observed and modeled representation of natural, low-frequency modes of weather and climate variability was also identified in the decadal survey as one of the most important challenges of the coming decade. Including interactions between large-scale circulation and organization of convection, such as the Madden–Julian oscillation (MJO) and the El Niño–Southern Oscillation (ENSO), aims to improve predictions by 50 % at lead times of 1 week to 2 months, which will have a high societal impact. It is therefore essential to further understand the physics and dynamics of the underlying processes, which are currently described with simple parameterizations in many models. Advanced observations of atmospheric convection and high-resolution models are also needed. While models will likely increase their nominal resolution in the next decade, it is also likely that global, century-long simulations from multi-ensembles under different assumptions will need to resort to parameterizing convection to reduce the computational burden.

Figure 2Simulated 6 h accumulated surface liquid precipitation for Typhoon Megi without using a CP (a) and using five different CPs in the WRF model. The accumulated precipitation includes cumulus, shallow cumulus, and grid-scale rain. The simulations start on 25 September 2016 at 18:00 UTC. The domain is located over the Philippine Sea with a horizontal grid size of 10 km. Radiation scheme: RRTMG shortwave and longwave schemes, boundary layer scheme: Mellor–Yamada–Janjić scheme, microphysics scheme: NSSL two-moment scheme, land surface option: unified Noah land surface model, surface layer option: eta similarity scheme. Spinning time: 24 h. GFS data were used to perform these simulations. The typhoon was not seeded.


2 Overview of the main schemes in cumulus convection modeling

Soon after Charney and Eliassen (1964) and Ooyama (1964) introduced the idea of cumulus parameterization, two approaches emerged: the convergence and the adjustment schemes (Arakawa, 2004). Later, a new scheme was introduced by Ooyama (1971): mass flux parameterization. Despite all these schemes attempting to explain the interaction between cumulus clouds and the large-scale environment, the choice of empirical values for certain parameters and the simplifications in the physics yield different convective parameterizations and strategies. Indeed, as shown in Fig. 2 for the 6 h total accumulated precipitation for Typhoon Megi, even today model outputs look different depending on the cumulus parameterization used. Many operational weather models and most climate models still use updated version of schemes described in the 1980s and 1990s. However, in recent years, new developments have emerged such as parameterizations including stochastic elements in the cumulus scheme, scale-aware approaches, or the addition of processes such as cold pools, among others (Rio et al., 2019). Many of these new schemes have been developed to simulate convection across the so-called gray zones, i.e., zones where traditional convective parameterizations are no longer valid but convection cannot be yet resolved explicitly (Wyngaard, 2004). Different treatments for shallow and deep convection have been traditionally used in convection parameterizations. However, this trend has changed towards a unified treatment in recent years based on the seamless transition between shallow and deep convection observed in nature (e.g., Park, 2014a, b).

As of 2021, the main cumulus convection schemes publicly available for NWPs are convergence schemes, adjustment schemes, mass flux schemes, cloud-system-resolving models (CSRM), super-parameterization (SP), PDF-based schemes (PDF: probability density function), unified models, scale-aware and scale-adaptive models, and models that account for convective memory and spatial organization The purpose of this paper is not to compare the performances of the schemes but to make explicit and investigate their empirical values and assumptions, so the focus of the following section is on these. The other drive of the paper, the assumptions in convective parameterizations, concerns the trigger model, the transport and microphysics, which are commonly referred to as the cloud model in classical convection schemes, and the closure of the scheme (Fig. 3 right). These are also described in the sections below.

Figure 3Schematic ensemble of cumulus cloud (a, c) and bulk convection scheme (b, d) showing the main components of a bulk convection scheme: trigger, updraft, downdraft, entrainment, detrainment, closure, conversion of cloud water to rainwater, precipitation and evaporation, and subsidence. Schemes based on Arakawa and Schubert (1974, a, c) and Bechtold (2019, b, d).


Figure 4Simulated 24 h position and pressure for Typhoon Megi (up) and Typhoon Chaba (down) using 15 ensembles in the ECMWF IFS model at 18 km horizontal grid size. Each marker represents one ensemble member. Square markers indicate observations. The simulations start on 26 September 2016 at 06:00 UTC for Typhoon Megi and on 3 October 2016 at 00:00 UTC for Typhoon Chaba. Figures on the left depict observations (obs) and perturbed initial conditions (pert init conds), while figures on the right show seven perturbed convection parameters (pert conv param) using the ECMWF stochastically perturbed parameterization (SPP). The perturbed parameters are organized entrainment, entrainment for shallow convection, turbulent detrainment, adjustment time, rain conversion, momentum transport, and shallow vs. deep cloud thickness.

2.1 Convergence schemes: the key role of the total moisture convergence parameter

Convergence schemes consider synoptic-scale convergence to destabilize the atmosphere, while the heat released through condensation in cumulus clouds stabilizes it. Typical examples of this approach are Charney and Eliassen (1964), Ooyama (1964), and Kuo (1974). Charney and Eliassen (1964) did not use cloud models to explain these interactions. Instead, the concept of conditional instability of the second kind (CISK) was introduced. In the tropical cyclone (TC) case, CISK states that cyclones provide moisture that maintains cumulus clouds, and cumulus clouds provide the heat that cyclones need. Ooyama (1964) used a similar formulation, but represented the heating released through condensation in cumulus clouds in terms of a mass flux and considered the entrainment of ambient air. Kuo (1965, 1974) used a simple cloud model scheme to describe the interaction between a large-scale environment and cumulus clouds. One of the key assumptions in this scheme is that the total moisture convergence can be divided into a fraction b, which is stored in the atmosphere, and the remaining fraction (1−b), which precipitates and heats the atmosphere. This parameter was further modified by Anthes (1977), who proposed a relationship between b and the mean relative humidity (RH) in the troposphere, with b≤1. In the evaluation of rainfall rates using the Global Atmospheric Research Program Atlantic Tropical Experiment (GATE) scale phase III, Krishnamurti et al. (1980) obtained the most realistic precipitation rates for b≈0 for the Kuo scheme (Kuo, 1974). This value of b is not realistic as it implies that no moisture is stored in the atmosphere. In a later paper, Krishnamurti et al. (1983) introduced an additional subgrid-scale moisture supply to account for the observed vertical distributions of heat and moisture that the Kuo scheme failed to reproduce, as well as to address the major limitation of b=0 reported in Krishnamurti et al. (1980). The total moisture supply was expressed as I=1+ηIL, with IL the large-scale moisture supply. The authors used a multiple regression approach to find the values of b and η. Another approach consists of using the wet-bulb characteristics to locally determine the partition between precipitation and moistening (Geleyn, 1985).

Due to its formulation, the Kuo scheme cannot produce a realistic moistening of the atmosphere and cannot represent shallow convection. Moreover, it assumes that convection consumes water and not energy, which violates causality (Raymond and Emanuel, 1993; Emanuel, 1994). Despite these drawbacks, it can produce acceptable results in various applications (e.g., Kuo and Anthes, 1984; Molinari, 1985; Pezzi et al., 2008), such as in GCMs and NWP models (e.g., Rocha and Caetano, 2010; Mbienda et al., 2017). This convective parameterization scheme demands the least computational power and is thus sometimes used for large, centennial simulations.

2.2 Adjustment schemes: two strategies to remove instability

In adjustment schemes, the atmospheric instability is removed through an adjustment towards a reference state. Therefore, the physical properties of clouds are implicit and no cloud model has to be explicitly specified. The first proposed adjustment scheme was the moist convective adjustment by Manabe et al. (1965), also known as the hard adjustment. In this parameterization, moist convection occurs if the air is supersaturated and conditionally unstable. The instability is removed through an instantaneous adjustment of the temperature to a moist adiabatic lapse rate and of water vapor mixing ratio to saturation. Moreover, all the condensed water in this process precipitates immediately. The main problems of this scheme are the production of very large precipitation rates and its saturated final state after convection, which is rarely observed in nature (Emanuel and Raymond, 1993).

The so-called soft or relaxed adjustment schemes attempt to alleviate these problems by assuming that the hard adjustment occurs only over a fraction a of the grid area or by specifying the final mean RH (Cotton and Anthes, 1992). For example, Miyakoda et al. (1969) defined saturation as 80 % RH, while Kurihara (1973) performed the adjustment based on the buoyancy condition of a hypothetical cloud element instead of the saturation criterion.

Further improvements to the adjustment schemes were introduced by Betts and Miller (1986), whose scheme is also known as a penetrative adjustment scheme. The authors proposed an adjustment of large-scale atmospheric temperature T and moisture q to reference profiles over a specified timescale τ (adjustment timescale):

(1) ( T / t ) cu = ( T ref - T ) / τ ( q / t ) cu = ( q ref - q ) / τ ,

where subscript “cu” refers to cumulus convection and “ref” to the reference profile for each field. The reference profiles, which are different for shallow and deep convection, are quasi-equilibrium states based on observational data from GATE, Barbados Oceanographic and Meteorological Experiment (BOMEX), and Atlantic Trade-Wind EXperiment (ATEX). For the construction of the temperature reference profile, Betts (1986) used a mixing line model (Betts, 1982, 1985). Then, the moisture reference profile was calculated from the temperature profile by specifying the pressure difference between air parcel saturation level and pressure level at cloud base, freezing level, and cloud top. Therefore, the three adjustment parameters used in this scheme are the adjustment timescale τ, the stability weight Ws, and the saturation pressure departure, Sp.

The sensitivity of the scheme to the adjustment parameters has been evaluated by numerous authors. For instance, Baik et al. (1990) analyzed the influence of different values of each adjustment parameter on the simulation of a tropical cyclone, while Vaidya and Singh (1997) did the same for the simulation of a monsoon depression using four sets of values, including those from Betts and Miller (1986) and Slingo et al. (1994). In all cases, the adjustment parameters had to be modified depending on the different climate regimes. While Baik et al. (1990) set Ws=0.95 and Sp=(-30,-37.5,-38) hPa as the optimal parameters to simulate a tropical cyclone, Vaidya and Singh (1997) obtained the best forecast for a monsoon depression with Ws=1.0 and Sp=(-60,-70,-50) hPa. Despite the improvements achieved through adjusting the parameters for different climate conditions, the original Betts–Miller scheme occasionally produced heavy spurious rainfall over warm water and light precipitation over oceanic regions (Janjić, 1994). To overcome this problem, Janjić (1994) proposed considering a range of reference equilibrium states and characterizing the convective regimes by a parameter called “cloud efficiency”, which is related to precipitation production and depends on cloud entropy. This parameter is the sort of empirical value that requires attention when future climates are to be simulated. The modified scheme, known as the Betts–Miller–Janjić (BMJ) scheme, is one of the most widely used adjustment schemes in NWP models (e.g., Vaidya and Singh, 2000; Evans et al., 2012; Fiori et al., 2014; Fonseca et al., 2015; García-Ortega et al., 2017), despite its large bias for light rainfall (e.g., Gallus and Segal, 2001; Jankov and Gallus, 2004; Jankov et al., 2005). Convective adjustment schemes are computationally efficient, which makes them suitable for large-scale simulations.

2.3 Mass flux schemes: assuming the rates of mass detrainment and entrainment

Because of the nature of both convergence and adjustment schemes, a cloud model does not have to be explicitly specified to describe the interaction between cumulus clouds and the large-scale environment. This is not the case for the mass flux schemes, wherein convective instability is removed through the vertical eddy transport of heat, moisture, and momentum. The main objective of mass flux schemes is to describe this convective vertical eddy transport in terms of convective mass flux (Plant and Yano, 2015). To do so, the total flux is defined as ωψ, where ω is the vertical velocity and ψ a physical variable, e.g., the total specific humidity q. Then, the total flux is expressed as the sum of a large-scale mean ωψ and an unresolved eddy contribution ωψ (Reynolds averaging). Decomposing the total flux into flux contributions from cumulus cover areas and environmental regions, defining an active cloud fractional area a, and again using Reynolds averaging, the turbulent flux is expressed as

(2) ω ψ = a ω ψ c + ( 1 - a ) ω ψ e + a ( 1 - a ) ( ω c - ω e ) ( ψ c - ψ e ) ,

where the overbar indexes “c” and “e” denote cloud (environmental) average of the fluctuations with respect to the cloud (environmental) average, and the superscripts “c” and “e” denote active cloud and passive environmental averages (Siebesma and Cuijpers, 1995). Commonly, the so-called “top-hat” approximation is used in convective schemes. This approximation implies neglecting the first two terms on the right-hand side in Eq. (2) in favor of the third one (the organized turbulent term due to organized updraft and compensating subsidence), which is considered dominant. Classical convective parameterizations further assumed that a is small compared to the large-scale system, i.e., a≪1 (e.g., Yanai et al., 1973; Arakawa and Schubert, 1974, hereafter AS). Then, the mass flux formulation using the definition of the convective mass flux is


where wc represents the in-cloud vertical velocity. The reader is referred to Bechtold (2019) and Siebesma and Cuijpers (1995) for detailed derivation of these equations. Using a simple entraining plume model and setting ρ to unity, the continuity equations for the mass, updraft properties, and vertical momentum are


where E and D refer to entrainment and detrainment rates, respectively, Sψ represents sources and sinks of ψ, ζ is a virtual mass parameter that reduces buoyancy due to the pressure gradient force, Pc includes pressure perturbations within the cloud, and the overbar denotes average values. The first formulation of this type was introduced by Ooyama (1971). The author assumed that cumulus clouds of different sizes coexist and that they could be represented by an ensemble of independent non-interacting buoyant elements. The definition of the so-called dispatcher function would close the parameterization. However, the author left this question open. Numerous schemes have been proposed since then, mostly using the steady-state assumption, i.e., /t=0 (e.g., Yanai et al., 1973; Arakawa and Schubert, 1974; Kain and Fritsch, 1990). As mentioned in de Roode et al. (2012), early mass flux schemes did not apply a vertical velocity equation for convective updrafts (Eq. 5c) and used an ad hoc assumption to specify the cloud top that depended on the vertical resolution. To alleviate this issue, recent mass flux parameterizations include a vertical velocity equation for updrafts in their formulation inspired by Simpson and Wiggert (1969):

(6) 1 2 w c 2 z = a w B - b w ε w c 2 ,

where ε is the fractional entrainment (E=εM), and aw and bw are tunable parameters related to pressure perturbation and sub-plume contributions, respectively (see Table 1). Since then, numerous convection schemes have applied equations similar to Eq. (6) for the in-cloud vertical velocity (e.g., Bechtold et al., 2001; Gregory, 2001; von Salzen and McFarlane, 2002; Jakob and Siebesma, 2003; Bretherton et al., 2004; Cheinet, 2004; Soares et al., 2004; Rio and Hourdin, 2008; Neggers et al., 2009; Pergaud et al., 2009; Rio et al., 2010; De Rooy and Siebesma, 2010; Kim and Kang, 2012; de Roode et al., 2012; Sušelj et al., 2012, 2013; Wang and Zhang, 2014; Morrison, 2016a, b; Peters, 2016; Suselj et al., 2019a). The reader is referred to de Roode et al. (2012) for a detail derivation of Eq. (6) from Eq. (5c) and a discussion about the values of the tunable parameters aw and bw.

Table 1A sample of values aw and bw used in Eq. (6). Based on de Roode et al. (2012).

Download Print Version | Download XLSX

To overcome the gray-zone issue, schemes should be scale-aware, which requires dropping the traditional assumption of a≪1 in convective parameterizations (Arakawa et al., 2011). Numerous cumulus schemes no longer use this assumption (e.g., Neggers et al., 2009; Arakawa and Wu, 2013; Grell and Freitas, 2014).

Mass flux convective parameterization schemes are still the most common convective parameterizations used in ESMs, regional climate models (RCMs), and NWP models.

2.4 Cloud-system-resolving models (CSRMs)

The performances of the previous schemes prompted the search for new strategies to model convection. Krueger (1988) put forward the CSRM idea (also known as explicit convection, convection-permitting, or cloud ensemble models) to explicitly simulate convective processes over a kilometer scale instead of using parameterizations. Most convective parameterizations tend to produce too little heavy rain and too much light rain (e.g., Dai and Trenberth, 2004; Sun et al., 2006; Dai, 2006; Allan and Soden, 2008; Stephens et al., 2010), though these results depend on the model used for the simulations, and have problems representing diurnal precipitation cycles over land (e.g., Yang and Slingo, 2001; Guichard et al., 2004). The use of convection-permitting models can solve errors associated with other convective parameterizations (e.g., Kendon et al., 2012; Prein et al., 2013; Brisson et al., 2016) but entails higher computational costs, which limits their application in climate modeling (e.g., Wagner et al., 2018; Randall et al., 2019). They are also increasingly used in NWP, however (e.g., Kain et al., 2006; Gebhardt et al., 2011). Recently, Prein et al. (2015) reviewed prospects and challenges in regional convection-permitting climate modeling.

2.5 Super-parameterization (SP)

Hybrid approaches also exist. SP (also known as cloud-resolving convective parameterization – CRCP – or multiscale model framework – MMF) is an approach between parameterized and explicit convection, which consists of replacing the convective parameterizations by 2D cloud-resolving models (CRMs), or even a 3D LES model, at each grid cell of a GCM (Grabowski and Smolarkiewicz, 1999; Grabowski, 2016). Randall et al. (2003) proposed SP as “the only way to break the cloud parameterization deadlock.” SP is mostly applied in GCMs (e.g., Grabowski, 2001; Khairoutdinov and Randall, 2003; Khairoutdinov et al., 2005; Zhu et al., 2009; Jung and Arakawa, 2014; Sun and Pritchard, 2016). Several studies have compared the performance of SP with convective parameterizations, in particular using the Community Atmosphere Model (CAM).

Among the most notable improvements achieved by SP in CAM are simulations of heavy rainfall events that are much more similar to observations, a better diurnal precipitation cycle over land (e.g., (Khairoutdinov et al., 2005; DeMott et al., 2007; Zhu et al., 2009; Holloway et al., 2012; Rosa and Collins, 2013), and the production of a realistic MJO (e.g., Thayer-Calder and Randall, 2009; Holloway et al., 2013). However, simulations with SP also have problems that need solving, such as the failure to simulate light rainfall rates reported by Zhu et al. (2009). The computational cost of this approach is also higher than the one for convective parameterizations (Krishnamurthy and Stan, 2015) but smaller than the computational cost for global CSRMs performing climate simulations (Randall et al., 2003).

2.6 PDF-based schemes

Numerous cloud and stochastic parameterizations are based on probability density functions (PDFs) of moist conserved thermodynamic variables. The so-called statistical schemes use PDFs to improve the simulations of cloud clover so important in the planetary energy budget (e.g., Cahalan et al., 1994; Bony and Dufresne, 2005; Neggers and Siebesma, 2013; Bony et al., 2015). To our knowledge, the first scheme suggesting a joint PDF to compute cloud cover was that of Sommeria and Deardorff (1977), followed by Mellor (1977). These schemes used a single-Gaussian PDF. Various PDF distributions have been proposed since the formulation of the first statistical scheme, including gamma (Bougeault, 1982), Gaussian (Sommeria and Deardorff, 1977; Mellor, 1977; Bechtold et al., 1992), triangular (Smith, 1990), uniform (Le Trent and Li, 1991), lognormal (Bony and Emanuel, 2001), beta (Tompkins, 2002), and double-Gaussian (Lewellen and Yoh, 1993; Larson et al., 2002; Golaz et al., 2002a; Naumann et al., 2013). Studies such as those of Tompkins (2002) and Watanabe et al. (2009) included prognostic equations for the shape parameters of the PDF, which reduced cloud cover bias when tested in ECHAM5 (Tompkins, 2002) and MIROC (Model for Interdisciplinary Research on Climate, Watanabe et al., 2009), respectively.

In the stochastic parameterization context, Craig and Cohen (2006) used statistical mechanics to describe fluctuations about a large-scale equilibrium to provide a theoretical basis for stochastic parameterizations. A PDF in the form of an exponential law provides random values of the mass flux per cloud. Plant and Craig (2008) followed this scheme and used a PDF in their formulation together with a plume model and closure assumption adapted from the Kain–Fritsch scheme (Kain and Fritsch, 1990, KF hereafter), while Teixeira and Reynolds (2008) obtained a stochastic component from a normal PDF to perturb the tendencies related to the convective parameterization. Tompkins and Berner (2008) used a similar approach to perturb the initial humidity of the convective parcel and/or the humidity of the air entrained during ascent. More recently, Sakradzija et al. (2015) extended the deep convective formulation in Plant and Craig (2008) to shallow convection.

PDFs are also used to unify the representation of moist convection and boundary layer turbulence into one single scheme (see Sect. 2.7). Randall et al. (1992) and Lappen and Randall (2001a, b) used double-delta PDF to model the subgrid-scale variability of vertical velocity, temperature, and moisture. The scheme is called assumed-distribution higher-order closure (ADHOC), and it is a combination of assumed distributions of higher-order closure and mass flux closure. Bechtold et al. (1995) used a positively skewed distribution function to account for shallow clouds. Later, Chaboureau and Bechtold (2002, 2005) extended this approach to include all types of clouds. Based on results from Larson et al. (2002) and the binormal model of Lewellen and Yoh (1993), Golaz et al. (2002a, b) proposed the Cloud Layers Unified By Binomials (CLUBB) approach that uses a double-Gaussian PDF instead of a double-delta PDF. More recently, Jam et al. (2013), Hourdin et al. (2013), and Qin et al. (2018) represented shallow cumulus clouds with the PDF variances diagnosed from the turbulent and shallow convective processes. In the context of the eddy diffusivity mass flux (EDMF) framework, Cheinet (2003, 2004) used a Gaussian distribution of the thermodynamic variables, Soares et al. (2004) parameterized cloudiness with a PDF, and Sušelj et al. (2012) and further modifications of the scheme (Sušelj et al., 2013, 2014; Suselj et al., 2019b, a) use a PDF to describe the moist updraft characteristics. Sakradzija et al. (2016) coupled the extension of the Plant and Craig (2008) described in Sakradzija et al. (2015) to the EDMF parameterization in ICON.

A number of studies that attempt to unify the representation of shallow and deep convection also use PDFs (e.g., Park, 2014a, b; see Sect. 2.8).

2.7 Unified models

Traditionally, models have used separate parameterizations for the boundary layer as well as shallow and deep convection. Deficiencies associated to deep convection schemes, such as the representation of the MJO, the diurnal cycle of precipitation, and the double Intertropical Convergence Zone (ITCZ), have been addressed by introducing different modifications in existing models. However, Guichard et al. (2004) showed that these modifications are not sufficient to resolve deficiencies of convection parameterization and stressed the necessity of using an ensemble of parameterizations that represents a succession of convective regimes. Numerous attempts to merge shallow and deep convection parameterizations into a single framework can be found in the literature (e.g., Bechtold et al., 2001; Kain, 2004; Kuang and Bretherton, 2006; Hohenegger and Bretherton, 2011; Mapes and Neale, 2011; D'Andrea et al., 2014; Park, 2014a, b). Hohenegger and Bretherton (2011) proposed a unified parameterization by modifying the University of Washington (UW) shallow convection scheme (Bretherton et al., 2004; Park and Bretherton, 2009) to make it more suitable for deep convection. The authors kept the assumption that mass flux at cloud base is proportional to the convection inhibition factor (CIN) / turbulent kinetic energy (TKE) but modified the proportionality factor following Fletcher and Bretherton (2010), who set it to 0.06. Besides, the increase of the average TKE over the depth of the boundary layer due to cold pools is included in the calculations of TKE and therefore in the closure. Mapes and Neale (2011) also modified the UW shallow convection scheme by making entrainment dependent on a prognostic variable called organization (see Sect. 2.9). Guérémy (2011) proposed a new mass flux scheme based on continuous buoyancy, and D'Andrea et al. (2014) extended the shallow convection of Gentine et al. (2013a, b) to deep convection. Park (2014a, b) described a unified convection scheme (UNICON) for both shallow and deep convection without relying on an equilibrium closure. The scheme diagnoses the dynamics, macrophysics, and microphysics of multiple plumes. Also, it includes a prognostic cold pool parameterization and mesoscale organized flow within the PBL, thus accounting for convective memory. Later, Park et al. (2017) modified UNICON to diagnose additional detrainment following Tiedtke (1993) and Teixeira and Kim (2008). More recently, Shin and Park (2020) developed a stochastic UNICON model wherein the correlated multivariate Gaussian distribution for updraft vertical velocity and thermodynamic scalars is used to randomly sample convective updraft plumes.

In general, models split the turbulence parameterization among the PBL and moist convection (usually based on different conceptual models), simplifying the treatment of turbulence but requiring the addition of an artificial closure to match both schemes (Sušelj et al., 2014). Examples of PBL schemes that produce precipitation include the IFS EDMF, the EDMF developed by Neggers (2009), and the CLUBB scheme implemented in CAM (Thayer-Calder et al., 2015), among others. To our knowledge, the first scheme proposing a unified scheme in this way was that of Chatfield and Brost (1987), which was further evaluated by Petersen et al. (1999) and extended by Lappen and Randall (2001a, b) (see Sect. 2.6 for further details). Golaz et al. (2002a, b) and Larson et al. (2002) proposed an approach to combine the representation of shallow convection and turbulence, the so-called Cloud Layers Unified By Binomials (CLUBB, Sect. 2.6 for more details). Efforts to apply CLUBB to deep convection include those of Cheng and Xu (2006) and Bogenschutz and Krueger (2013) in CRMs or Davies et al. (2013a) in an SCM. To improve deep convective simulations, Storer et al. (2015) and Thayer-Calder et al. (2015) used a subgrid importance Latin hypercube sampler (SILHS; Larson et al., 2005; Larson and Schanen, 2013) that draws samples from the joint PDF to drive microphysical processes. More recently, Larson (2020) described the unified configuration of CLUBB-SILHS, wherein no separated deep parameterization is used (the reader is referred to this paper for a detailed explanation of CLUBB-SILHS).

The EDMF approach was proposed by Siebesma and Teixeira (2000) and Siebesma et al. (2007) to overcome the commonly ad hoc matching between the mass flux approach for convective transport within the clouds and the eddy diffusivity approach to parameterize turbulent transport in the atmospheric boundary layer. Starting from Eq. (2), assuming a≪1, and identifying the third term in the equation with the convective mass flux,

(7) ω ψ = ω ψ e + M ( ψ c - ψ ) .

Then, the first term in Eq. (7) is approximated by an eddy diffusivity approach (Siebesma et al., 2007).

(8) w ψ - K ψ z + M ( ψ u - ψ )

Thus, the transport in the atmospheric boundary layer is determined as the sum of an eddy diffusivity component, defined as the product of a diffusivity coefficient K and the local gradient of a thermodynamic state variable ψ, and a mass flux part, defined as the product of a mass flux and the difference between ψ in the updraft and its horizontal mean value. The authors used a K profile (Holtslag, 1998) for the eddy diffusivity coefficient, took the updraft fractional area as a constant, and scaled the mass flux with the standard deviation of the vertical velocity σw. Despite being originally used for dry convective boundary layers (Siebesma and Teixeira, 2000; Siebesma et al., 2007; Witek et al., 2011), numerous versions of the scheme extended it to moist convection (e.g., Soares et al., 2004; Angevine, 2005; Rio and Hourdin, 2008; Neggers et al., 2009; Neggers, 2009; Pergaud et al., 2009; Angevine et al., 2010; Köhler et al., 2011; Sušelj et al., 2012, 2013, 2014; Suselj et al., 2019a, b).

Besides extending the EDMF model to moist convection, a number of versions included a multiple-plume formulation. For example, Cheinet (2003) combined the EDMF model with the multi-parcel model described in Neggers et al. (2002). With the goal of finding the least complex mass flux framework that can reproduce the smoothly varying coupling between the sub-cloud mixed layer and the shallow convective cloud layer, Neggers et al. (2009) and Neggers (2009) proposed a new formulation combining the EDMF concept with a dual mass flux (DualM) framework. There, two different updrafts are considered: a dry updraft and a moist updraft. Each of the updrafts are characterized by an area fraction (see Table 15) that varies in time, with a continuous area partitioning between moist and dry updraft. In order to realistically represent not only convectively driven boundary layers but also the transition between shallow and deep convection, Sušelj et al. (2013) further developed the scheme described in Sušelj et al. (2012). One of the main innovations included the use of a Monte Carlo sampling of the PDF of updraft properties at cloud base. Sušelj et al. (2014) described a simplified version of the Sušelj et al. (2013) stochastic model wherein the eddy diffusivity parameterization is based on Louis (1979), among other modifications. Later, Tan et al. (2018) extended the EDMF approach by using prognostic plumes and adding downdrafts, among other changes.

Neggers (2015) reformulated the EDMF approach in terms of discretized size densities with a limited number n of bins. This new version, referred as to ED(MF)n, was studied in an SCM. J. Han et al. (2016) proposed a hybrid EDMF parameterization whereby EDMF is used only for the strongly unstable PBL. For weakly unstable PBL, the scheme uses a nonlocal PBL scheme with an eddy diffusivity counter-gradient approach (Deardorff, 1966; Troen and Mahrt, 1986; Hong and Pan, 1996; Han and Pan, 2011). Han and Bretherton (2019) replaced the ED parameterization in this scheme with a new TKE-based moist EDMF parameterization for vertical turbulence mixing, included downdrafts, and assumed a decrease in the updraft mass flux with decreasing grid size, which makes the scheme scale-aware. More recently, Wu et al. (2020) implemented a new downdraft parameterization in EDMF through a Mellor–Yamada–Nakanishi–Niino (MYNN) ED component. Kurowski et al. (2019) implemented a stochastic multi-plume EDMF scheme into CAM5, and Sakradzija et al. (2016) coupled Sakradzija et al. (2015) to EDMF in ICON. Several NWP models have included EDMF approaches, i.e., ECMWF (Köhler, 2005; Köhler et al., 2011), AROME (Pergaud et al., 2009), NCEP GFS (J. Han et al., 2016), the Navy Global Environmental Model (NAVGEM) (Sušelj et al., 2014), and the Laboratoire de Météorologie Dynamique Zoom (LMDZ; Hourdin et al., 2013) model. Recently, Bhattacharya et al. (2018) and Wu et al. (2020) implemented different versions of the EDMF scheme in WRF.

2.8 Scale-aware and scale-adaptive models

Wyngaard (2004) coined the terms terra incognita or “gray zone” to refer to zones where traditional convective parameterizations are no longer valid but convection cannot be resolved explicitly yet. To palliate the gray zone parameterizations should become scale-aware and scale-adaptive. This means that the scheme is aware of the processes that need to be parameterized and parameterizes only those processes. Recently, Honnert et al. (2020) reviewed schemes that have been proposed for the convective boundary layer in the gray zone.

In the context of mass flux representations, the quasi-equilibrium (QE) assumption on a negligible small cloud area fraction a has to be eliminated to make parameterizations scale-aware (Arakawa et al., 2011). Arakawa et al. (2011) and Arakawa and Wu (2013) described a seamless approach in their unified parameterization wherein the assumption about a is eliminated, the vertical eddy transport is rederived, and the parameterization is forced to converge to an explicit simulation as a→1. Following this approach, Grell and Freitas (2014) extended the Grell and Dévényi (2002) scheme based on Grell (1993) by specifying a as a function of the convective updraft radius R obtained from the traditional definition of entrainment ε (Siebesma and Cuijpers, 1995; Simpson and Wiggert, 1969; Simpson, 1971), i.e., ε=0.2/R. Later, Freitas et al. (2017) tested this scheme in the Brazilian developments of the Regional Atmospheric Modeling System (BRAMS) version 5.2, obtaining a smooth transition between convective and grid-scale precipitation even at gray-zone scales.

Lim et al. (2014) modified the simplified Arakawa–Schubert scheme (SAS; e.g., Grell, 1993; Pan and Wu, 1995; Hong and Pan, 1998; Han and Pan, 2011) in NCEP GFS by introducing a grid-scale dependency in the trigger. More recently, Kwon and Hong (2017) extended this grid-scale dependency to the convective inhibition, mass flux, and detrainment of hydrometeors, and Han et al. (2017) updated the SAS scheme with a cloud mass flux that decreases with increasing grid resolution to include scale dependency.

Zheng et al. (2016) modified the adjustment timescale in the KF scheme following Bechtold et al. (2008) and included a scale-aware entrainment equation, among other modifications.

Other approaches to overcome the gray-zone issue include spreading subsidence to neighboring cells in the Grell3D scheme (Grell and Freitas, 2014) and a hybrid parameterization for nonhydrostatic weather prediction models as described in Kuell et al. (2007). This scheme uses a traditional cumulus parameterization for mass and energy transport in the updraft and downdraft and treats environmental subsidence by grid-scale equations. More recently, Freitas et al. (2018) implemented and tested a new version of the Grell and Freitas (2014) scheme in the NASA Goddard Earth Observing System (GEOS). The new scheme uses a trimodal formulation with different entrainment rates that depend on the normalized mass flux profile, which is prescribed by a beta PDF, among other modifications. Gao et al. (2017) compared the performance of the traditional KF scheme with the Grell and Freitas (2014) scheme in the simulation of summer precipitation across gray-zone resolutions. Better results were reported with the scale-aware scheme. An integrated package of subgrid- and grid-scale parameterizations in the range 2–10 km, also known as the Modular Multiscale Microphysics and Transport (3MT), was proposed by Gerard (2007). Zheng et al. (2016) added scale awareness to the KF scheme (Kain and Fritsch, 1990, 1993; Kain, 2004) by introducing scale dependency in in-cloud properties, such as entrainment and grid-scale vertical velocity.

Another way to introduce scale awareness and adaptivity consists of using multiple plumes instead of a single one. The first scheme using multiple plumes is that of Arakawa and Schubert (1974). Different schemes have been proposed based on multiple plumes for deep (e.g., Donner, 1993; Donner et al., 2001; Nober and Graf, 2005; Wagner and Graf, 2010) and shallow convection (e.g., Neggers et al., 2002; Sušelj et al., 2012; Neggers, 2015). Due to the lack of observations on cloud entrainment, Neggers et al. (2002) used LES results to formulate an expression for the lateral entrainment rate as a function of the vertical velocity of each parcel, while Sušelj et al. (2012) described moist updraft characteristic through a PDF. Other parameterizations, such as those of Wagner and Graf (2010), Nober and Graf (2005), and Neggers and Siebesma (2013), make use of active population dynamics such as those in the Lotka–Volterra equations (Lotka, 1910, 1920; Volterra, 1926), wherein two species interact with a predator–prey behavior. Neggers (2015) also introduced population dynamics in a new EDMF called the ED(MF)n. The author used bin macrophysics, wherein plumes are described in terms of discrete size densities formed by a limited number n of bins. The scale adaptivity of this scheme was further evaluated in Brast et al. (2018). Population dynamics were also used by Park (2014a, b) in his multi-cloud model in UNICON and by Hagos et al. (2018) in the STOchastic framework for Modeling Population dynamics of convective clouds (STOMP), among others. Khouider et al. (2010) described a stochastic multi-cloud model based on the deterministic multi-cloud model of Khouider and Majda (2006) but using a Markov chain lattice model. In this scheme, four possible convective states in each lattice are considered, namely clear sky, deep, congestus, or stratiform clouds, that randomly evolve in time as a birth–death process (Gillespie, 1975, 1977). Dorrestijn et al. (2013a, b, 2015) also used this approach but estimate transition probabilities from one state to another using LES results and observations, respectively. Further works followed, such as those of Deng et al. (2015) for representing the MJO, the coupling of Khouider et al. (2010) to simplified primitive equations of Frenkel et al. (2012), the use of observations to estimate transition probabilities in Peters et al. (2013), and the implementation of a stochastic multi-cloud scheme in ECHAM6.3 by Peters et al. (2017), among others. Later, Khouider (2014) improved Khouider et al. (2010) by using a coarse-grained Markov chain lattice model. Examples of stochastic parameterizations based on concepts from statistical mechanics include Plant and Craig (2008) for deep convection or Sakradzija et al. (2015, 2016) and Sakradzija and Klocke (2018) for shallow convection. Recently, Keane et al. (2014) evaluated the scale adaptivity of Plant and Craig (2008) in the ICON model. Rochetin et al. (2014a, b) added a stochastic component to the trigger function in LMDZ5B, and Sakradzija et al. (2016) introduced scale awareness in the ICON model by coupling the stochastic scheme described in Sakradzija et al. (2015) to the EDMF scheme. Other scale-aware schemes include CLUBB due to its limitation of the turbulent length scale to the horizontal grid spacing (Larson et al., 2012).

Other studies have included a scale-dependent entrainment and/or convective timescale (e.g., Bechtold et al., 2008; Zheng et al., 2016; Han et al., 2017; Gao et al., 2020) based on results obtained in entrainment and mixing studies (e.g., Burnet and Brenguier, 2007; Lu et al., 2011, 2014; Kumar et al., 2018; Kooperman et al., 2018).

The best way to achieve scale-aware and scale-adaptive cumulus schemes is still unknown, but the field is rapidly evolving.

2.9 Models accounting for convective memory and spatial organization

As pointed out in Davies et al. (2009), the QE hypothesis does not account for convective memory, which can be defined as the dependence of convection on past states. Different strategies have been proposed to include it in convective parameterizations, such as the use of prognostic variables or cold pools, among others. The first scheme to include convective memory was that of Pan and Randall (1998). The authors chose a cumulus kinetic energy prognostic closure in their formulation. Later, Gerard and Geleyn (2005) also accounted for convective memory. Based on Bougeault (1985), the authors defined cloud-base mass flux as the product of a prognostic vertical updraft velocity and a prognostic updraft fraction area, obtained by a moist static energy closure. Gerard (2007) and Gerard et al. (2009) also used this approach and even applied it for downdrafts (Gerard et al. 2009). Piriou et al. (2007) used precipitation evaporation as the source of convective memory and related entrainment to the probability of undiluted updrafts. Mapes and Neale (2011) also chose precipitation evaporation as the source of convective memory and introduced a prognostic variable called organization that links precipitation evaporation with the entrainment rate. Other authors selected the precipitation at convective cloud base as the source of convective memory and made entrainment a function of it (e.g., Hohenegger and Bretherton, 2011; Willett and Whitall, 2017, in the UK Met Office model). Another way to introduce convective memory consists of using a master equation or Markov chains, such as the schemes of Hagos et al. (2018) and Khouider et al. (2010). In their extended EDMF, Tan et al. (2018) included convective memory using prognostic equations for updrafts and downdrafts as well as for the area fraction (see Table 15).

Evaporation of precipitation from deep convective clouds gives rise to cold pools that, when spread at the surface, are able to initiate further convective events, therefore adding memory to the system (e.g., Khairoutdinov and Randall, 2006; Rio et al., 2009; Böing et al., 2012; Schlemmer and Hohenegger, 2014). Based on this, recent studies include convective memory through cold pools (e.g., Grandpeix and Lafore, 2010; Park, 2014a, b; Del Genio et al., 2015). The prognostic variables are the cold pool thermodynamic properties and fractional area (Grandpeix and Lafore, 2010) as well as the cold pool depth (Del Genio et al., 2015) or the mesoscale organized flow (Park, 2014a, b). More recently, Colin et al. (2019) performed numerical experiments to identify the source of convective memory using CRMs. The results showed that memory comes from low-level thermodynamic process such as rain evaporation, cold pools, or hot thermals, among others.

Based on the “game of life” (Chopard, 2009), Bengtsson et al. (2011) used a cellular automaton (CA) in their subgrid scheme. The authors introduced convective memory by assigning a prescribed lifetime to each active cell. Bengtsson et al. (2013) also included memory in their stochastic parameterization for deep convection using this approach in Aire Limitée Adaptation/Application de la Recherche à l'Opérationnel (ALARO). The definition of the area fraction in the cumulus scheme (Gerard et al., 2009) now includes the contribution from CA. Sakradzija et al. (2015) accounted for convective memory by considering that the cloud rate distribution in shallow convection comes from the superposition of two modes. These two modes consider passive and active clouds, respectively. In their work, the authors considered convective memory due to the finite lifetime of individual clouds. Later, Sakradzija et al. (2016) used this scheme in the calculation of the moist convective area fraction in EDMF in ICON.

Results from Davies et al. (2013b) suggested that spatial organization could strongly affect convective memory more than the microphysics parameterizations. Later, Colin (2020) confirmed this hypothesis.

Understanding spatial organization of convection is not only important for developing stochastic and scale-aware parameterizations but also due to its impact in the radiative–convective equilibrium (Neggers and Griewank, 2021) . Few studies have proposed parameterizations to represent convective organization in GCMs (e.g., Donner, 1993; Donner et al., 2001; Mapes and Neale, 2011; Donner et al., 2011; Khouider and Moncrieff, 2015; Moncrieff et al., 2017). Donner (1993), Alexander and Cotton (1998), and Donner et al. (2001) represented the effects of mesoscale circulations and downdrafts based on the Leary and Houze (1980) water budget model. A similar model was developed by Gray (2000), who also considered momentum fluxes and related the strength of mesoscale circulation to detrainment of the convective mass flux. As mentioned before, Mapes and Neale (2011) introduced a prognostic variable called organization into the UW shallow convection scheme (Bretherton et al., 2004; Park and Bretherton, 2009). This variable, which represents the degree of subgrid organization, could affect plume calculations in terms of plume-base vertical velocity, convective inhibition, preferential rising of warmer air in updrafts, area fraction and closure, and a shift in the spectrum toward wider plumes with lower lateral mixing and a preferential growth in preconditioned local environments. All this would lead to more and deeper convection and therefore more organization.

Other studies accounted for convective organization by including surface cold pools in their convective parameterizations (e.g., Rio et al., 2009; Grandpeix and Lafore, 2010; Rochetin et al., 2014a, b; Park, 2014a, b; Böing, 2016). Grandpeix and Lafore (2010) proposed a density current parameterization based on the first convective wake parameterization described by Qian et al. (1998). The impact of the cold pools on convection is implemented through two variables: the available lifting energy (ALE) provided by the density current and the available lifting power (ALP; see Sect. 5.1.1). In the UNICON model, Park (2014a) parameterized subgrid mesoscale convective organization in terms of the evaporation of convective precipitation and downdrafts. Later, Böing (2016) described an object-based model of the organization of moist convection by cold pools inspired by Abelian sandpile models (Bak et al., 1987). The model is a two-way feedback between instability and convection, whereby convection and instability are represented as particles coupled to a lattice grid. The authors suggested that an object-based model could capture properties of convective organization. Stratton and Stirling (2012) used the height of the lifting condensation level as a variable to introduce convective organization into their parameterization, while Folkins et al. (2014) introduced a dependency on the local precipitation generated by the convective scheme over the past 2 h. Khouider and Majda (2006) developed a multi-cloud parameterization wherein three cloud types control the heating fields of organized convection in the tropics. It was later refined by Khouider and Majda (2008) and applied by Khouider and Moncrieff (2015) in their parameterization of organized convection in the ITCZ. Moncrieff et al. (2017) proposed a new method referred to as multiscale coherent structure parameterization (MCSP) to parameterize physical and dynamical effects of organized convection. This new approach consists of using a slantwise overturning model with a special focus on top-heavy heating and upgradient momentum transport. Despite all these proposals, the model of Donner et al. (2011) is the only operational GCM representing all aspects of mesoscale convective systems (Rio et al., 2019).

In Shutts (2005) the spatial and temporal correlations of the atmospheric mesoscale are represented by a CA. Bengtsson et al. (2011) extended the implemented CA in the ECMWF Ensemble Prediction System to be able to interact with the numerical model. Later, Bengtsson et al. (2013) introduced this CA approach in ALARO and analyzed it in a regional gray-zone resolution model over Europe. This approach produced a precipitation intensity and convective organization in better agreement with OPERA observations than results obtained from the reference model. In Bengtsson et al. (2019), CA is conditioned by a prescribed stochastically generated skewed distribution with the goal of introducing subgrid-scale organization.

Other attempts to represent convective organization include the use of a damped-driven oscillator (Davies et al., 2009), spatially coupled oscillators (Feingold and Koren, 2013), or a Markov chain lattice model (e.g., Khouider et al., 2010). Moncrieff and Liu (2006) proposed a hybrid approach to represent convective organization. Mesoscale organization is represented by explicit convectively driven circulations using a CSRM and transient cumulus by the BMJ convective parameterization (Betts, 1986; Betts and Miller, 1986; Janjić, 1994). PDF-based or spectral schemes based on a discretized distribution (e.g., Neggers et al., 2003; Wagner and Graf, 2010; Neggers, 2012; Park, 2014a, b; Neggers, 2015) include size information into the system, which allows representing impacts of spatial organization (Neggers et al., 2019; van Laar, 2019). More recently, Neggers and Griewank (2021) developed a binomial stochastic framework referred to as the Binomial Objects on Microgrids (BiOMi) model, which proved to capture convective memory and simple forms of spatial organization, among other important convective behaviors, at a cheap computational cost.

This paper considers all the aforementioned convective parameterizations with an emphasis on the mass flux schemes.

3 Trigger function: assumptions and empiricisms

In a CP, the accurate simulation of convection greatly depends on the trigger function. The trigger function determines whether convectively unstable air at the boundary layer leads to the onset of convection and, if so, activates the CP.

There are as many strategies to initiate convection as there are convection schemes. This section focuses on the assumptions and empirical values of the most important trigger functions, the starting levels, and the impacts of the trigger formulations on the simulation of convective processes. Table 2 lists the most common choices used in the main trigger function types.

3.1 Trigger function types

According to the physical variable used as the main trigger condition, the most used trigger functions in CPs may be classified into (1) moisture convergence, (2) cloud work function (CWF), (3) cloud-base stability and convective available potential energy (CAPE) triggers, and (4) large-scale vertical velocity. Other triggers used are (5) stochastic and heated condensation framework (HCF) triggers. Table 2 lists the assumptions and empirical values used in the main trigger function types, which are discussed below.

3.1.1 Moisture convergence trigger

The main condition to activate convection, together with the existence of a deep layer of conditional instability, is exceeding a minimum threshold value of the vertically integrated moisture convergence. This is the case in the Anthes–Kuo scheme (Kuo, 1965; Anthes, 1977) and in the original Tiedtke scheme (Tiedtke, 1989). The latter has undergone several modifications since its publication. For instance, Gregory et al. (2000) substituted the condition of positive moisture convergence to activate deep convection by a minimum cloud-depth threshold in the European Centre for Medium-Range Forecasts (ECMWF) convective parameterization. Other authors replaced the moisture convergence trigger in the Tiedtke scheme with triggers based on positive buoyancy (Zhang et al., 2011) or the existence of an unstable parcel within some height above the ground (Bechtold et al., 2004). Therefore, these schemes are no longer classified as moisture convergence triggers.

Table 2A sample of empirical values and assumptions used in the main trigger function types.

Download Print Version | Download XLSX

3.1.2 CWF trigger

The first CWF trigger was introduced by AS, who proposed that convection activation depends on a threshold value of the CWF, which is defined as the integral buoyancy force of each entraining cloud between cloud base and cloud top. Several variations of the original CWF trigger function have been suggested. Tokioka et al. (1988) included a modification in the AS to suppress deep convection in areas where the depth of the PBL is not sufficiently thick. This modification is defined on a critical value of the entrainment rate below which deep convection is suppressed and moist air can accumulate in the large-scale low-level convergence zone. For example, the GFDL global atmosphere and land model (AM2–LM2; Anderson et al., 2004) includes this modification. In the relaxed Arakawa–Schubert scheme (RAS) (Moorthi and Suarez, 1992), the activation of convection depends on a critical value of the CWF, while the SAS scheme (Grell, 1993; Pan and Wu, 1995) triggers convection if the CWF is positive, as shown in Table 2. Another condition to activate convection in SAS is based on the pressure difference between the starting level, i.e., the level of maximum moist static energy between the surface and the 700 hPa level, and the level of free convection (LFC), which defines a threshold value for the convection inhibition (CIN) factor. With the aim of decreasing convection in large-scale subsidence regions and increasing it in large-scale convergent regions, Han and Pan (2011) modified the limit to reach the LFC, which is now proportional to large-scale vertical velocity ω. Further improvements to the SAS activation criteria include a grid spacing dependency in the convective trigger function (Lim et al., 2014), considering the spatial resolution dependency, and a new definition of the CIN threshold value applying a scale-aware factor (Kwon and Hong, 2017). Different versions of the AS scheme are currently used in the Global Forecast System (GFS) of the National Centers for Environmental Prediction (NCEP), the Mesoscale Model 5 (MM5), the Goddard Earth Observing System model version 5 (GEOS-5), the Geophysical Fluid Dynamics Laboratory (GFDL) model, and the WRF model.

To improve the representation of the diurnal cycle, Rio et al. (2009) proposed a new trigger for deep convection: the so-called available lifting energy (ALE). This trigger is defined as the kinetic energy of the parcel inside thermals and activates deep convection when it overcomes CIN. In this case, convection activation is controlled by lifting processes in the sub-cloud layer, e.g., gust fronts. The authors obtained a better representation of the diurnal cycle with their new formulation. Grandpeix and Lafore (2010) also used the ALE trigger in their coupled wake–convection scheme. Together with a closure based on the flux of kinetic energy associated with thermals and the splitting of convective heating and drying, a more realistic representation of moist convection was possible. More recently, Hourdin et al. (2013) confirmed these results in the implementation of the ALE trigger into a new version of the LMDZ atmospheric general circulation (LMDZ5B).

3.1.3 Cloud-base stability and CAPE triggers

Many CPs have been proposed to simplify the formulation and implementation of the AS scheme. Among other assumptions, some CPs substitute the convection trigger based on CWF by CAPE, defined in a similar way as CWF but without including dilution of an ascending parcel by entrainment. For instance, BMJ developed a new parameterization based on empirical results, in which the activation of convection requires the existence of CAPE. In this scheme, cloud base is the lifting condensation level (LCL) of a lifted parcel with the largest CAPE in the lowest 130 hPa of the model. From there, the parcel is lifted moist adiabatically until the equilibrium level (EL) is reached. In general, the cloud top is at the level immediately beneath EL. Moreover, deep convection continues if the cloud depth is greater than a certain value and covers at least two model layers (Baldwin et al., 2002). Finally, deep convection activates if the adjustment using reference profiles of temperature (based on a moist adiabat) and moisture (based on imposed subsaturation at the cloud base) results in the column drying. The reference profiles computed in the BMJ scheme are different for shallow and deep convection.The scheme is currently used in the NCEP North American Mesoscale model (NAM), MM5, and the WRF model. Another important convective parameterization also using a CAPE trigger is the Zhang–McFarlane scheme (Zhang and McFarlane, 1995, hereafter ZM). To improve climate simulations in the Canadian Climate Center GCM, the authors proposed a simplified version of the AS scheme that includes a positive CAPE trigger. However, it initiates convection too often during the day, which led Xie and Zhang (2000) to modify the scheme. They kept the positive CAPE condition and added a second condition based on the change in CAPE due to large-scale forcing (dCAPE). This new trigger improved the simulations of the ITCZ and MJO (Zhang, 2002; Song and Zhang, 2009; Zhang and Song, 2010). Alternative formulations of convection triggers include the addition of an RH threshold of 80 % in the convection trigger (Zhang and Mu, 2005a, b) to suppress convection if the boundary layer air is too dry. Another modification is the inclusion of dilution in CAPE calculation due to entrainment (dilute CAPE) by Neale et al. (2008) to reduce excessive precipitation over land in the simulations of ENSO.

Unlike some of the trigger criteria already discussed, a more recent trigger function by Tawfik and Dirmeyer (2014), the HCF, is not based on the lifting parcel method but uses vertical profiles of temperature and humidity. First, it finds the buoyant condensation level (BCL) and determines several variables such as the buoyant mixing potential temperature, θBM, defined as the 2 m potential temperature needs to reach the BCL, and the potential temperature deficit, θdef, defined as the difference between the θBM and the 2 m potential temperature, or the sum of all the temperature increments needed to attain the BCL. In HCF, convection will activate when θdef≤0. The HCF trigger reduces the number of false positives compared to the parcel-based trigger. When the HCF trigger is implemented in the NCEP Climate Forecast System version 2 (CFSv2), the representation of the Indian monsoon and tropical cyclone intensity improves (Bombardi et al., 2016). In the Community Earth System Model (CESM), the strategy improves the frequency of heavy precipitation events and reduces the overactivation of convection in the model (Tawfik et al., 2017).

3.1.4 Large-scale vertical velocity trigger

Drawing on the observations in Fritsch and Chappell (1980) suggesting a positive impact of background vertical motion on convective development, Kain and Fritsch (1990) (KF) proposed a trigger based on large-scale vertical velocity. In this scheme, the first potential source layer for convection, also known as the updraft source layer (USL), is a layer of at least 60 hPa thickness that is constructed by mixing vertically adjacent layers, beginning at the surface. The temperature and pressure of the parcel at its LCL are calculated, as is temperature perturbation δT, which is proportional to ω (see Table 3). If the sum of the parcel temperature and the temperature perturbation is higher than the environmental temperature, the parcel is released from its LCL. Above the LCL, the parcel is lifted upwards with entrainment, detrainment, water loading, and a vertical velocity determined by the Lagrangian parcel method (Bechtold et al., 2001). Convection is activated if the vertical velocity remains positive for a minimum depth of 3–4 km. Otherwise, the USL is moved up one model level and the procedure starts again. This process continues until a suitable USL is found or the search has moved up above the lowest 300 hPa of the atmosphere, where the search is terminated. The lake-effect snow observations of Niziol et al. (1995) forced a reduction of the minimum cloud-depth threshold in Kain and Fritsch (1993) from 3–4 to 2 km as they showed that clouds with this depth can produce significant snowfall. In Plant and Craig (2008), the temperature perturbation to find the USL is set to 0.2 as in Gregory and Rowntree (1990). If no buoyant source layer can be found, then the process (like in KF) is repeated with a temperature perturbation of 0.1 K. The plume radii are determined with an exponential PDF.

Other authors, such as Ma and Tan (2009), included moisture advection in the temperature perturbation to improve the KF scheme for the case of weak synoptic forcing. Berg et al. (2013) defined a PDF that generates a range of virtual potential temperature and water vapor mixing ratio to substitute δT in the trigger function. With this new trigger, the scheme more realistically accounts for subgrid variability within the convective boundary layer in a way. Both the modified version of the KF scheme and the KF itself are used in the WRF model.

As for the trigger of shallow convection, Bechtold et al. (2001) proposed a deep convective scheme based on Kain and Fritsch (1990, 1993) but also included a shallow parameterization. In this regard, the triggering criterion is only based on a cloud-depth condition without using the temperature perturbation included in the deep scheme. Also, cloud-depth condition and cloud radius take smaller values than those use for deep convection (see Table 3). Jakob and Siebesma (2003) also used a cloud-depth condition to decide whether deep or shallow convection is triggered. In this case, the maximum value of the cloud depth to activate shallow convection is set to 200 hPa. The procedure of finding cloud base is the same for both parameterizations.

In the shallow convection parameterization for mesoscale models described in Deng et al. (2003) based on Kain and Fritsch (1990, 1993), maximum cloud depth is set to 4 km and cloud radius is allowed to increase smoothly with time from a minimum value of 0.15 km to a maximum value of 1.50 km. Moreover, shallow convection trigger is a function of boundary layer TKE. In Han and Pan (2011), the USL is set to the level of maximum moist static energy within the PBL, and the maximum cloud top for shallow convection is restricted by the ratio between the layer pressure and surface pressure that cannot be higher than 0.7. A cloud-depth criterion to activate shallow or deep convection is also used in this case. Han et al. (2017) developed a scale-aware parameterization for NCEP GFS, wherein the cloud-depth criterion is increased to 200 hPa compared to the 150 hPa used in Han and Pan (2011).

In Kain (2004) the conditions to trigger shallow convection are the same as for deep convection except for the cloud depth, which must be smaller than the one for deep convection (see Table 3). In this parameterization, the values of cloud radius are the same for both shallow and deep convection for computational reasons. Bretherton et al. (2004) trigger convection if the vertical velocity of the parcel is equal to or higher than a critical value derived from the vertical velocity equation (Eq. 6). This critical velocity takes the form wcrit,sh=2aw(CIN), where aw is the virtual mass coefficient used in the updraft vertical velocity equation (Eq. 6; see de Roode et al., 2012). Park and Bretherton (2009) used the same triggering conditions as Bretherton et al. (2004).

Table 3A sample of empirical values and assumptions used in the trigger (note: subscript “sh” refers to shallow convection).

Download XLSX

3.1.5 Stochastic trigger

The traditional convective triggers lead to deficiencies in the simulation of different atmospheric events, as stated in Sect. 2. A promising strategy to reduce these deficiencies is the use of stochastic triggering (Rochetin et al. 2014a, b). Instead of using a deterministic parameterization in which the subgrid-scale response is fixed to a certain resolved-scale state, the response is sampled from a suitable probability distribution (Dorrestijn et al., 2013b). For example, Majda and Khouider (2002) and Khouider et al. (2003) used a stochastic model based on CIN using a Markov process. Stechmann and Neelin (2011) used a two-state Markov jump process as their stochastic trigger. Bright and Mullen (2002) modified the KF trigger function by applying stochastic perturbation to w, while Song et al. (2007) included several random parameters in the trigger criteria using a Bayesian learning procedure. Zhang et al. (2014) added a stochastic term to the SAS trigger function in the Hurricane Weather Research and Forecasting model (HWRF), and Rochetin et al. (2014a, b) used LES to introduce a stochastic trigger in the Emanuel parameterization (Emanuel, 1991).

3.2 Starting levels

The LFC is located at, or near, the cloud base or at the top of the PBL. Different methods are applied for calculating the LFC in the literature. For instance, KF used the potential source layers for clouds (USL) in their procedure to find LFC, while Pan and Wu (1995) first determined the convection starting level and then imposed a critical depth to find the LFC (see Sect. 3.1). In their stochastic parameterization, Plant and Craig (2008) set the depth of potential source layers to 50 hPa, with the base of each being 5 hPa higher than the potential layer previously tested. To trigger convection, both deep and shallow, Han and Pan (2011) set a threshold value for the pressure difference between LFC with and without sub-cloud-layer entrainment. Differences higher than this threshold value, set to 25 hPa, will activate convection. The authors also assumed that the convection starting level for deep convection is at the level of maximum moist static energy h between the surface and the level of 700 hPa, while for shallow convection it starts at the level of maximum h within the PBL. Table 4 lists a sample of the main assumptions and empirical values used to determine the starting levels.

Table 4A sample of empirical values and assumptions used in the starting levels (note: subscript “sh” refers to shallow convection).

Download Print Version | Download XLSX

While the starting level for the ascending currents (updrafts) is reasonably evident, the starting level for the descending currents (downdrafts), usually called the level of free sinking (LFS), may start at any vertical level no lower than the cloud base. Several convective parameterizations, such as those proposed by Tiedtke (1989) or Bechtold et al. (2001), follow the definition suggested by Fritsch and Chappell (1980), who assumed that LFS is the level at which the temperature of a saturated mixture of equal amounts of updraft and environmental air becomes smaller than the environmental temperature. In contrast, Grell et al. (1991) determined LFS as the minimum value of h, and Zhang and McFarlane (1995) matched LFS with the lowest updraft detrainment level. However, if the minimum value of h is lower than the bottom level of updraft detrainment, LFS is determined as in Grell (1993).

3.3 Impact of trigger functions on convective models

Differences between trigger functions depend on the identification of the source layer of convective air and on how this layer of unstable air can give rise to convection. While near-surface air is selected as the source layer in some CPs (Tiedtke, 1989; Donner, 1993; Bechtold et al., 2001; Tawfik and Dirmeyer, 2014), in others, the choice is the layer of maximum moist static energy, h (Arakawa and Schubert, 1974; Grell, 1993; Zhang and McFarlane, 1995; Wu, 2012). On the other hand, different convection triggers are used to determine whether unstable air turns into convection, as mentioned in the previous section. However, the best way to construct a trigger function is still unknown and, in many cases, an ad hoc formulation leads to poor performance in the activation of convection at the right location and time (Suhas and Zhang, 2014; Song and Zhang, 2017). Comparison between the performance of different trigger functions and observations from different climates leads to improvements in the formulation of the activation criteria for convection. Suhas and Zhang (2014) used three intensive observation period (IOP) datasets from the Atmospheric Radiation Measurement (ARM) program and long-term single-column models (SCMs) to evaluate the performance of different trigger functions (AS scheme, Bechtold scheme, Donner scheme, KF scheme, Tiedtke scheme, and four variants of the ZM scheme). The dilute dCAPE trigger function showed the best performance in both the tropics and midlatitudes, while the undilute dCAPE was as good as the dilute dCAPE only for the tropics. Furthermore, the Bechtold and dilute CAPE trigger functions were among the best-performing schemes. As a follow-up, Song and Zhang (2017) used observations from the Green Ocean Amazon (GOAmazon) field campaign to evaluate and improve the trigger functions selected in Suhas and Zhang (2014), with the addition of the HCF. In their study, the dCAPE-type triggers also ranked first, followed by the Bechtold and HCF triggers. The undilute dCAPE trigger performed better with the inclusion of a 700 hPa upward motion, while the dCAPE trigger improved with an optimization of the entrainment rate and dCAPE threshold. Using GOAmazon, the authors set the values for the dCAPE threshold and entrainment rate. The new values are 55 J kg−1 s−1 for the dCAPE threshold and 2.5×10-4m-1 for the entrainment rate.

The convection trigger criterion plays a crucial role in the simulation of a wide number of atmospheric events. The impact of the trigger function on the correct simulation of the diurnal cycle of convection and precipitation in atmospheric models has been widely studied, especially over land (Bechtold et al., 2004; Knievel et al., 2004; Lee et al., 2007a, b, 2008; Hara et al., 2009; Evans and Westra, 2012). The common problem in the simulation of the diurnal cycle is that it peaks too early and its amplitude is too high (Yang and Slingo, 2001; Collier and Bowman, 2004). Moreover, the diurnal cycle of precipitation peaks too early over land (in general, 2 to 4 h before the observed maxima) (Dai, 2006), which is related to the formulation of the trigger function (Betts and Jakob, 2002; Bechtold et al., 2004). Lee et al. (2008) performed a sensitivity analysis with four different trigger functions implemented in the RAS scheme and found significant differences in the diurnal cycle of precipitation over the Great Plains in the United States. Several studies have performed sensitivity analyses and found possible ways to improve the simulation of the diurnal cycle. Models with finer resolution provided a better simulation in the amplitude, variability, and timing of the diurnal cycle (Wang et al., 2007; Sato et al., 2009). The inclusion of the effect of moisture advection in the trigger function improved the distribution and intensity of convective precipitation in the MM5 (Ma and Tan, 2009). The use of different initiation and termination conditions in the SAS scheme led to a better diurnal variation of precipitation (Han et al., 2019), although it increased the excessive precipitation and did not alleviate the bias in the phase of precipitation intensity. The modification of both the trigger and closure criteria by considering cold pools could minimize the bias in the diurnal cycle of convection (Rio et al., 2009, 2013). Another important case involves the deficiencies in the simulation of the MJO (Lin et al., 2006), which are often improved by the modification of the trigger function. For example, Wang and Schlesinger (1999) found that a better representation of the MJO was possible by adding a moisture trigger to the convective parameterization used in the atmospheric general circulation model at the University of Illinois, Urban–Champaign (UIUC). Zhang and Mu (2005b) used the same approach in the National Center for Atmospheric Research (NCAR) Community Climate Model version 3 (CCM3) as did Lin et al. (2008) in the Seoul National University (SNU) atmospheric general circulation model. Another example is a better representation of the Indian summer monsoon rainfall by the addition of HCF to the trigger function in the Climate Forecast System version 2 (CFSv2) (Bombardi et al., 2015).

The lack of “convective memory” effects in the models based on the QE assumption causes a convective parameterization to be triggered, regardless of the convection stage, as long as the convection criteria are met. Different ways to include the memory effect have been proposed, such as using prognostic cumulus kinetic energy (Pan and Randall, 1998), or an ensemble of cold pools (Grandpeix and Lafore, 2010; Del Genio et al., 2015) (see Sect. 2.9).

4 Cloud model: types and choices

The cloud model represents the interaction between cumulus clouds and the large-scale environment. Thus, it determines the vertical distribution of convective heat and moisture through the parameterization of the mass flux profile, the entrainment–detrainment, and the microphysics. This section discusses the main types of mass flux and entrainment–detrainment schemes adopted in the literature, as well as the main assumptions and empirical values employed in the formulation of the cloud model.

4.1 Mass flux scheme types

According to the approach used to estimate the unknown quantities in Eq. (5a), (5b), and (5c), mass flux schemes are classified into spectral, bulk, and episodic mixing models.

4.1.1 Spectral models

Spectral models represent the ensemble of clouds within a grid box with a spectrum of clouds, each of them with a cloud model. Therefore, multiple types of convection are considered in these models in contrast to the bulk ones, wherein the use of only one cloud model for each grid box makes it necessary to decide a priori the type of convection and to characterize the cloud model by averages over the ensemble of clouds.

In spectral models, clouds within a grid box are grouped into different cloud models according to a certain parameter. The majority of spectral schemes generate an ensemble of plumes based on a distribution of entrainment rates (Arakawa and Schubert, 1974; Hack et al., 1984; Nober and Graf, 2005; Chikira and Sugiyama, 2010), although care has to be taken such that the results (convective regime) are not dominated by the least entraining parcels. Each cloud type contributes a different amount to the ensemble mean depending on their cloud-base mass flux. This type of model was original proposed by AS. Since then, the scheme has undergone several modifications, some of them make the scheme no longer a spectral model but a bulk mass flux scheme (e.g., Grell, 1993; Pan and Wu, 1995). For example, Moorthi and Suarez (1992) modified the closure in the AS scheme by replacing the QE assumption with a relaxation towards equilibrium. This scheme is also known as the relaxed Arakawa–Schubert (RAS). Numerous studies described models based on the spectral representation (e.g., Wagner and Graf, 2010; Donner, 1993; Sušelj et al., 2012, 2013; Hong et al., 2013; Neggers, 2015; Olson et al., 2019; Brast et al., 2018; Hagos et al., 2018).

4.1.2 Bulk models

The ensemble of clouds within a grid box is represented by a single cloud model, in contrast to spectral models. Yanai et al. (1973) describe the main representatives of this type of scheme. In their diagnostic study, clouds are classified according to their cloud tops, and the steady plume hypothesis (Morton et al., 1956) is applied. It is assumed that all clouds have a common cloud-base height and that the values of detrainment are identical to the values inside the plume. In mesoscale models, Fritsch and Chappell (1980) and Kain and Fritsch (1992) also applied the steady hypothesis, as did Singh et al. (2019) in their study of the relationship between humidity, instability, and precipitation in the tropics. Tiedtke (1989) and Gregory and Rowntree (1990) applied the same approach as Yanai et al. (1973) in their schemes at the ECMWF and at the UK Meteorological Office. The scheme used at ECMWF has undergone several modifications since then (e.g., Nordeng, 1994; Gregory et al., 2000; Li et al., 2007; Zhang et al., 2011; Kim and Kang, 2012; Stevens et al., 2013). Other studies, such as Grell (1993), changed the spectrum of cloud sizes in AS for a simple non-entraining cloud within a single grid box. Pan and Wu (1995) developed the so-called simplified Arakawa–Schubert model (SAS), which is a modified version of the model proposed by Grell (1993). The cloud ensemble is also represented by a single non-entraining cloud, and the downdraft starting level is modified to avoid excessive cooling below cloud base. Han and Pan (2011) further modified entrainment, detrainment, and cloud-base mass flux in SAS to overcome unrealistic grid-scale precipitation and developed a bulk mass flux parameterization for shallow convection. Many mass flux parameterizations use the bulk-cloud approach (e.g., Siebesma and Holtslag, 1996; Bechtold et al., 2001; Neggers et al., 2009; Yano and Baizig, 2012; Loriaux et al., 2013) with different formulations of their cloud models (i.e., formulation of the mass flux at cloud base, entrainment, detrainment, microphysics).

4.1.3 Episodic mixing models

Drawing on the continuous entrainment and average buoyancy used in entraining–detraining plume models in both bulk and spectral formulations, Emanuel (1991, 1994) proposed the so-called episodic mixing model, which is based on the stochastic mixing model of Raymond and Blyth (1986) and the observations of Taylor and Baker (1991), among others. Thus, Emanuel assumed that mixing is highly inhomogeneous and episodic and applied the buoyancy sorting hypothesis (Telford, 1975; Taylor and Baker, 1991), which is the basis of a number of cumulus parameterizations (e.g., James and Markowski, 2010; Park, 2014a), especially those focused on shallow convection (e.g., Bretherton et al., 2004; De Rooy and Siebesma, 2008; Neggers et al., 2009; Pergaud et al., 2009). The Emanuel scheme and its modified versions (Emanuel and Živković-Rothman, 1999; Grandpeix et al., 2004; Peng et al., 2004) are widely used in RCMs (e.g., Zou et al., 2014; Raju et al., 2015; Bhatla et al., 2016; Gao et al., 2016; Kumar and Dimri, 2020).

The aforementioned mass flux scheme types are explained from the point of view of the ascending currents. However, convective downdrafts, i.e., descendent currents caused by evaporation of condensate and rainwater loading, should be taken into account. Simply put, they may be considered bottom-up updrafts. Downdrafts are of great importance in atmospheric convection. As Plant and Yano (2015) highlighted, they have opposite effects on the organization and evolution of convective systems. The transport of cooler and drier air into the sub-cloud layer may stabilize it and therefore inhibit convection or may lead to the development of new convective elements if downdrafts cause an increase in low-level convergence. The majority of convective parameterizations include downdrafts with assumptions about their starting level, entrained and detrained air, and the amount of condensate available for evaporation. However, many schemes, such as Grell (1993), the ZM scheme used in CESM, and the Tiedtke scheme in the ECHAM model, have described downdrafts as simple saturated plumes, i.e., “inverse plume”, with a mass flux proportional to the updraft mass flux (Thayer-Calder, 2012). Other authors have proposed a more complex parameterization including unsaturated downdrafts in their formulations and a downdraft mass flux based on Eq. (5a), (5b), and (5c) (e.g., Emanuel, 1991; Xu et al., 2002).

4.2 Entrainment and detrainment

The mixing of air masses due to entrainment of environmental air into clouds and detrainment of cloudy air into the environment is a key process in convective parameterizations (Blyth, 1993; Luo et al., 2010; Donner et al., 2016) as it modifies the vertical profiles of heat and moisture within cloudy air. Sanderson et al. (2008) identified the entrainment rate as one of the dominant parameters affecting climate sensitivity after evaluating thousands of GCM simulations. Other authors, such as Rougier et al. (2009), Klocke et al. (2011), and Zhao (2014), have obtained similar conclusions in their analyses. In addition, the influence of convective detrainment of water vapor and hydrometeors from cumulus clouds is an important source of water that strongly impacts climate simulations (e.g., Ramanathan and Collins, 1991; Lindzen et al., 2001).

In this section, attention is drawn to the most important model types of entrainment and detrainment, the main assumptions and empirical values used in the literature, and the impact that the different formulations have in convective models. The main assumptions and empirical values used in the formulation of entrainment and detrainment are listed in Tables 5 and 6 and in Tables 7 and 8, respectively.

4.2.1 The choice of lateral vs. cloud-top entrainment

Since Stommel (1947) provided the first description of cumulus cloud dilution by entrainment of environmental air, two conceptual models are still competing: the lateral entrainment model and the cloud-top entrainment model.

In the lateral entrainment model, Stommel (1947) considered environmental air to enter the cloud through the lateral cloud edges and continuously dilute cloudy air during its ascent, regardless of whether it is considered a plume or a bubble. Several aircraft observations and experiments in water tanks (Turner, 1962; Morton, 1965) contributed to the formulation of the lateral entrainment theory. However, authors such as Warner (1970) pointed out the deficiencies of this theory in predicting the right profile of liquid water content (LWC).

In order to address these deficiencies, Squires (1958) proposed another entrainment model, the cloud-top entrainment. This author suggested that environmental air enters the cloud predominantly at or near the cloud top, descends through penetrative downdrafts created by evaporative cooling, and dilutes the cloud by turbulent mixing. Paluch (1979) provided more evidence for cloud-top entrainment in her study on cumulus clouds over Colorado. The author found that the cloud water mixing ratio and the wet equivalent potential temperature follow a line at a single level, the so-called “mixing line”, which connects cloud base and cloud top. Paluch interpreted it as evidence for a two-point mixing scenario. Further studies (Boatman and Auer, 1983; Lamontagne and Telford, 1983; Jensen et al., 1985; Reuter and Yau, 1987) confirmed Paluch's results. However, several authors have criticized the mixing line source levels (e.g., Blyth et al., 1988; Malinowski and Pawlowska-Mankiewicz, 1989; Raga et al., 1990; Grabowski and Pawlowska, 1993; Neggers et al., 2002; Zhao and Austin, 2005) and the interpretation of the mixing line (e.g., Betts and Albrecht, 1987; Taylor and Baker, 1991; Grabowski and Pawlowska, 1993; Siebesma, 1998; Böing et al., 2014).

Which of the two models predominates in cumulus convection remained unclear for many years. The increase in computational power in recent decades has promoted the use of LES to study entrainment and detrainment mainly in shallow cumulus clouds. Several authors, such as Heus et al. (2008) and Böing et al. (2014), have applied LES to identify the dominant process in mixing in cumulus clouds, concluding that cloud-top entrainment is insignificant compared to lateral entrainment.

4.2.2 Main empirical values in entrainment and detrainment formulations

Aircraft observations and experiments in water tanks (Turner, 1962; Morton, 1965) led to the formulation of the lateral entrainment theory, which anticipates that the fractional entrainment rate (hereafter entrainment rate) changes with the cloud radius (Malkus, 1959; Squires and Turner, 1962; Simpson and Wiggert, 1969; Simpson, 1971):

(9) 1 M M z = ε C R ,

where M is the mass flux, z is the height, ε denotes the entrainment rate, C is a constant, and R is the radius of the rising plume. These first parameterizations set C=0.2 based on laboratory results. As De Rooy et al. (2013) pointed out in their review article on entrainment and detrainment in cumulus convection, many cloud models still use this formulation (e.g., Arakawa and Schubert, 1974; Kain and Fritsch, 1990; Donner, 1993), sometimes assuming a constant entrainment rate.

Houghton and Cramer (1951) improved this theory by taking into account the increase in vertical velocity due to buoyancy. Thus, the authors distinguish between dynamical entrainment due to larger-scale organized inflow, εdyn, and turbulent entrainment caused by turbulent mixing, εturb. The turbulent entrainment rate is related to the flux across the updraft boundary, which is often described with an eddy diffusivity approach (Kuo, 1962; Asai and Kasahara, 1967; De Rooy et al., 2013; Cohen et al., 2020). Under the eddy diffusivity approach, the eddy flux is modeled by a downgradient and an eddy diffusivity that for the case of the turbulent entrainment is proportional to the radial scale of a plume (used as a mixing length) and the turbulent velocity scale of the environment. The change in mass flux with height, including the detrainment δ of negatively buoyant mixtures, is given by

(10) 1 M M z = ε dyn + ε turb - δ dyn - δ turb ,

Tiedtke (1989) and Nordeng (1994) assumed that turbulent entrainment is inversely proportional to cloud radii, as in Simpson and Wiggert (1969) and Simpson (1971). They used typical cloud sizes based on observations for different types of convection to fix the values of entrainment rates. For penetrative and midlevel convection, the entrainment rate was fixed to εturb=1×10-4m-1. This is a typical value for tropical clouds as shown in the analysis of aircraft observations in Simpson (1971). For shallow convection, the entrainment rate was based on typical values for large trade cumuli: εturb=3×10-4m-1 (Nitta, 1975). Gregory and Rowntree (1990) also assumed a turbulent entrainment rate, but inversely proportional to the height, while in Bechtold et al. (2008), εturb is O(1×10-3m-1) in better agreement with CRM results and also relative-humidity-dependent, which turned out to be important to represent realistic tropical variability (Table 5). Dynamical entrainment εdyn is proportional to moisture convergence and occurs only in the lower part of the cloud layer up to the level of strongest vertical ascent in Tiedtke (1989). In Nordeng (1994), it is based on momentum convergence. Gregory and Rowntree (1990) did not include it in their parameterization, whereas in Bechtold et al. (2008), it depends on RH and is only applied to deep convection. For downdraft, Bechtold et al. (2014) set εturb=3×10-4m-1 and εdyn as a function of B. A common practice in the definition of entrainment rates for downdraft consists of assuming a similar parameterization as for updrafts (Table 6).

Kain and Fritsch (1990) introduced another type of parameterization based on the buoyancy sorting. In their parameterization, homogeneous mixing of cloudy and environmental air was assumed, leading to mixtures with different buoyancy properties that have the same probability of occurrence. Moreover, the authors modified Eq. (9) to make it pressure-dependent. The fraction of environmental air that makes the mixture neutrally buoyant is the so-called critical mixing fraction χc, which determines whether a mixture entrains or detrains after mixing. Thus, entrainment of positively buoyant mixtures occurs if χ<χc, while χ>χc leads to immediate detrainment of negatively buoyant mixtures. Therefore, detrainment can occur at any level at which χ>χc, unlike in the AS scheme, wherein only the cloud-top detrainment is considered. Moreover, the maximum entrainment rate is proportional to pressure and inversely proportional to updraft radius. However, the KF scheme had deficiencies, such as excessive detrainment and the production of unrealistic deep saturated layers. In newer versions of the KF scheme, a mitigation of unrealistic deep saturated layers is achieved by assuming that the entrainment of environmental air cannot be lower than 50 % of the total environmental air involved in the mixing process in the updraft and that cloud radius depends on the convergence of the sub-cloud layer (Kain, 2004). Recently, Zheng et al. (2016) modified the minimum entrainment equation in Kain (2004) to include both organized and turbulent entrainment. The authors made the equation scale-dependent and expressed it in terms of sub-cloud-layer depth instead of cloud radius. Another scheme based on the buoyancy sorting hypothesis, but assuming episodic mixing, is the Emanuel scheme (Emanuel, 1991), wherein, in contrast to the KF scheme, the resulting mixtures just ascend or descend to their level of neutral buoyancy to detrain.

Other approaches use in-cloud quantities instead of only the environmental quantities to estimate the entrainment rate. For instance, Gregory (2001) proposed an entrainment rate that depends on B and inversely on the square of the updraft speed w calculated using Eq. (6). The value of aw also comes from the equation and is selected by comparing SCM simulations against LES and CRM studies as well as available observations. This parameterization deals with both shallow and deep convection. What distinguishes one type of convection from another is the value of a constant Ce, whose values were specified by using an SCM in ECMWF model.

Apart from buoyancy, another environmental quantity that might influence entrainment, and therefore convection, is RH. A number of studies have analyzed the effect of RH in parameterization of entrainment–detrainment rates, drawing different conclusions. For instance, Jensen and Del Genio (2006) found a positive correlation between entrainment rate and RH in their analysis of remote sensing observations and soundings at Nauru Island, while Bechtold et al. (2008) and Zhao et al. (2018) found a negative correlation using the Atmospheric Model version 4 (AM4.0). The same conclusion was achieved by Stirling and Stratton (2012) using a CRM formulation and the Met Office Unified Model (Met Office UM).

Mapes and Neale (2011) addressed the so-called “entrainment dilemma”, in which excessive entrainment values tend to excessively restrain convection, while insufficient entrainment values abundantly ease its activation. To overcome this, they proposed a new formulation of the entrainment rate dependent on a prognostic variable called organization, which expresses the interaction between the environment and convection. In their formulation, the rain evaporation rate controls the organization and produces more deep convection for lower values of the entrainment rate.

The previous discussion about entrainment and detrainment rates was focused on deep convective schemes with some references to unified schemes. However, parameterizations of these processes are also important in shallow convection. Tiedtke (1989) fixed the entrainment and detrainment rates for shallow convection to ε=δ=3×10-4m-1 based on typical values for large trade cumuli (Nitta, 1975). Using LES based on BOMEX, Siebesma and Cuijpers (1995) found typical values of entrainment for the core between 1.5×10-3 and 2×10-3m-1 and around 3×10-3m-1 for the updraft. Siebesma (1998) found typical values for entrainment in shallow convection in the range 1.5–2.5×10-3m-1. In their revision and performance analysis of the ECMWF IFS, Gregory et al. (2000) found values of ε=1.2×10-3m-1 at cloud base and ε=3×10-3m-1 150 hPa above it by employing a control physics package that included a cloud scheme based on Tiedtke (1989, 1993).

Grant and Brown (1999) and Grant and Lock (2004) described a similarity theory for shallow convective transport. In this theory, buoyancy production and turbulent dissipation are assumed to nearly balance within QE shallow convective fields. As for the entrainment formulation, it is scaled based on observable quantities such as CAPE or mass flux at cloud base with a constant Aε that represents the fraction of TKE available for entrainment. The value of this constant is derived from LES results. Kirshbaum and Grant (2012) used this formulation with Aε=0.06. Drueke et al. (2019) also used this TKE similarity theory for cloud ensembles to retrieve values of entrainment rates based on sub-cloud and environmental conditions. The authors also compared this method with the parcel model of Jensen and Del Genio (2006), which coupled surface remote sensing observations and soundings at Nauru Island to a parcel model, and the Entrainment Rate In Cumulus Algorithm (ERICA) proposed by Wagner et al. (2013), which uses an algorithm to retrieve values of entrainment from ground-based remote sensing observations. The analysis was performed using LESs of a range of shallow cumulus over ocean and land, showing a strong contrast in entrainment between them, as well as a lower dilution for wider clouds. The parcel method and TKE similarity theory better capture the sensitivity within continental cumuli and showed a lower mean error compared to ERICA. The diurnal variations of entrainment within continental shallow cumulus were only reproduced by the TKE method. With this method, the authors found values of Aε in the range 0.037–0.035. More recently, Kirshbaum and Lamer (2021) performed a climatological sensitivity analysis of shallow cumulus entrainment in oceanic and continental locations using the parcel method and the TKE as in Drueke et al. (2019). A total of 4 years of observations at two ARM observatories were used. The analysis confirmed the results obtained by Drueke et al. (2019) and identified other sources of entrainment variability such as sub-cloud wind speed in oceanic flows and cloud-base mass flux in individual cumuli. Median values of entrainment at a continental site range between 0.5 and 0.6 km−1 and between 1.0 and 1.1 km−1 at the oceanic site.

Neggers et al. (2002) developed a new formulation using LES. The authors proposed an entrainment rate inversely proportional to a turnover timescale that seems to be independent of cloud depth and the vertical velocity of the parcel. Thus, each parcel will have its own entrainment rate depending on its vertical velocity. For the ensemble of parcels, the fractional entrainment rate is of the order of the values shown in Siebesma and Cuijpers (1995). Sušelj et al. (2012) followed Neggers et al. (2002) but with a different value of the turnover timescale (see Table 5). Model results using an SCM proved to be sensitive to the choices of this parameter.

In their EDMF model, Soares et al. (2004) used a constant entrainment rate within the cloud layer following the entrainment rate in Siebesma (1998), while in the sub-cloud layer the entrainment is inversely proportional to height.

Bretherton et al. (2004) proposed an entrainment formulation similar to that of KF but modified χc by defining a critical eddy mixing distance dc based on observations and LES results that revealed fractions of negatively buoyant air in the updrafts (Taylor and Baker, 1991; Siebesma and Cuijpers, 1995). The so-called fractional mixing rate ε0 is defined as inversely proportional to the top of the cumulus layer H. In their unified scheme, Hohenegger and Bretherton (2011) applied the buoyancy sorting idea to compute entrainment and detrainment rates as in Bretherton et al. (2004), defining ε0 in a different way. Taking into account LESs performed with the System for Atmospheric Modelling (SAM), this value is link to the convective precipitation at cloud base here (see Table 5).

Based on the results obtained from using tracers in LESs of shallow convection during BOMEX that pointed to a description of entrainment through a stochastic Poisson process, Romps and Kuang (2010b) developed a parcel model with stochastic entrainment similar to the one proposed in Romps and Kuang (2010a). The authors used a Monte Carlo method to model entrainment rate. The parameterization uses two probability functions characterized by two parameters, i.e., the mean ratio of the entrained mass ment and the distance that parcel travels between entrainment events dent. The mean fractional entrainment per distance is given by the ratio of these two parameters. The values that best fit the CRM results were dent=226 mm and ment=0.91, i.e., ε=4.0×10-3m-1. Nie and Kuang (2012) specified ment=0.32 and dent=125 m for their LESs of BOMEX to reduce the number of undilute updrafts to a number comparable to their 25 m resolution run. For the sub-cloud layer, the parameters were set to dent=30 m and ment=0.06. Sušelj et al. (2013) replaced the entrainment parameterization in Sušelj et al. (2012) with a stochastic formulation. The authors considered a constant entrainment rate for dry updrafts below the condensation level and an entrainment formulation similar to the one proposed by Romps and Kuang (2010b). In this case, the authors found a typical distance of 100 m between entrainment events for BOMEX phase-3 experiment. Sušelj et al. (2014) parameterized the entrainment rate as in Sušelj et al. (2013) although with different values for the constant entrainment rate and dent.

Recently, in their shallow cumulus study, Lu et al. (2018) identified deficiencies in previous studies about the impact of RH on entrainment that could lead to erroneous conclusions regarding the effects of RH on entrainment, such as the use of conserved quantities related to RH to estimate entrainment rates and that no observations had thus far been used to determine the relationship between RH and entrainment. To address these deficiencies, the authors analyzed aircraft observations from the Routine AAF (ARM Aerial Facility) CLOWD (Clouds with Low Optical Water Depths) Optical Radiative Observations (RACORO) (Vogelmann et al., 2012) and Rain In Cumulus over the Ocean (RICO) field campaigns (Rauber et al., 2007) for shallow cumulus and concluded that ε and RH are positively correlated. Nonetheless, there is no general consensus on the effects of environmental RH on entrainment rates (Lu et al., 2018).

Table 5A sample of empirical values and assumptions used in the parameterization of entrainment in the updraft (note: subscript “sh” refers to shallow convection).

Download XLSX

Table 6A sample of empirical values and assumptions used in the parameterization of entrainment in the downdraft.

Download Print Version | Download XLSX

Less attention has been paid to the parameterizations of the detrainment process. Many convection schemes set it as a constant value (see Tables 7 and 8), while others consider detrainment to be negligible (Lu et al., 2012). Tiedtke (1989) and Nordeng (1994) assumed a turbulent detrainment inversely proportional to cloud radii and fixed its value to δturb=1×10-4m-1 for penetrative and midlevel convection (see Table 7). On the other hand, Gregory and Rowntree (1990) assumed a turbulent detrainment rate inversely proportional to the height and smaller than εturb, while Bechtold et al. (2008) set δturb to a constant value. Dynamical detrainment δdyn is defined to occur in Tiedtke (1989), Bechtold et al. (2008), and Gregory and Rowntree (1990) when the updraft buoyancy becomes negative. In the former two schemes it is then set proportional to the decrease in updraft kinetic energy, while in the latter it is computed implicitly. For downdraft, Bechtold et al. (2014) set δturb=εturb and enforced δdyn over the lowest 50 hPa. As in the case of entrainment rates in downdrafts, a common practice in the definition of detrainment rates for downdraft consists of assuming a similar parameterization as for updrafts (Table 8).

Table 7A sample of empirical values and assumptions used in the parameterization of detrainment in the updraft (note: subscript “sh” refers to shallow convection).

Download Print Version | Download XLSX

Table 8A sample of empirical values and assumptions used in the parameterization of detrainment in the downdraft.

Download Print Version | Download XLSX

In the parameterization of detrainment in shallow convection schemes, De Rooy and Siebesma (2008) treated the mass flux and the entrainment formulation separately based on LES results, which suggests that variations in the mass flux profile are mostly related to the fractional detrainment (Jonker et al., 2006; De Rooy and Siebesma, 2008). De Rooy and Siebesma (2008) kept ε fixed as an inverse function of height and developed a dynamical formulation for δ dependent on the average of χc from cloud base up to the middle of the cloud layer χc (the reader is referred to Eq. A11 in De Rooy and Siebesma, 2008, for a detailed calculation of χc) and on the cloud layer depth. For shallow convection, Siebesma and Cuijpers (1995) found vales of detrainment rates that were rather constant at around 3×10-3m-1 for the core and 4×10-3m-1 for the updraft. Using LES output from BOMEX, Siebesma (1998) found typical values of detrainment in the range 2.5–3×10-3m-1. Other studies, such as Soares et al. (2004), used a constant detrainment rate following Siebesma (1998) and set it to the value of entrainment at cloud base (e.g., Han and Pan, 2011) or proportional to the entrainment rates (e.g., Bechtold et al., 2014), among others.

4.2.3 Impact of entrainment and detrainment on convective models

The discussion above illustrates the many nuances in the modeling of convection, the importance of empirical values in the final results, and the need for further research to disentangle the many details involved. It is accepted that the parameterizations of entrainment and detrainment still have great uncertainties (e.g., Romps, 2010; Becker and Hohenegger, 2018) and problems in producing a realistic representation of convection (e.g., Mapes and Neale, 2011). For example, Stratton and Stirling (2012) improved the timing and amplitude of the diurnal cycle of tropical convection in the Met Office climate model by setting the entrainment for deep convection as a function of the height of LCL.

Perhaps not surprisingly, MJO simulations are also sensitive to entrainment (e.g., Hannah and Maloney, 2011; Del Genio et al., 2012; Kim et al., 2012; Hirons et al., 2013; Klingaman and Woolnough, 2014). Hannah and Maloney (2011) applied the RAS scheme in a GCM and analyzed the influence of minimum entrainment rate and rain evaporation fraction in the simulation of MJO. Larger values of either of the two parameters led to a better representation of the MJO and interseasonal variability, although higher values of minimum entrainment produced a drier and cooler atmosphere in contrast to the effect of higher values of rain precipitation fraction. Klingaman and Woolnough (2014) evaluated the effects of 22 model configurations and subgrid parameterizations on the simulation of MJO in the Hadley Centre Global Environmental model Global Atmosphere version 2 (HadGEM3 GA2.0) and tested the changes in 14 hindcast cases. A better representation of the MJO for both hindcast and climate simulations was achieved by increasing entrainment and detrainment rates for midlevel and deep convection. A better representation of MJO was also achieved by Kim et al. (2012) using a GCM to evaluate the tropical subseasonal variability. However, this improvement was at the expense of an increased bias in the mean state, which is typical for other GCMs with stronger MJO (Kim et al., 2011).

The entrainment parameterization proposed by Gregory (2001) for both deep and shallow convection achieved satisfactory results in various analyses (e.g., Chikira and Sugiyama, 2010; Del Genio and Wu, 2010) but proved to be cloud- and altitude-dependent. Recently, Baba (2019) modified Gregory's parameterization of the entrainment rate by relating it to the detrainment rate and B. This new parameterization led to improvements in the positive bias of precipitation in western Pacific region, in the positive bias of outgoing shortwave radiation over the ocean, and in the simulation of MJO, equatorial waves, and precipitation over the western Pacific region. Using an RCM over the Maritime Continent region, Wang et al. (2007) demonstrated that changes in the values of the fractional entrainment–detrainment rates in Tiedtke scheme, including both shallow and deep convection, affect the simulation of the tropical precipitation diurnal cycle. Over land, Del Genio and Wu (2010) used a CRM to study the transition from shallow to deep convection in diurnal cycles and inferred entrainment rates. Subsequently, the authors compared results from three different entrainment parameterizations to the results obtained with the CRM and concluded that the best results were achieved by the entrainment parameterization of Gregory (2001). Through a version of the Goddard Institute for Space Studies Global Climate Model (GISS GCM) with the entrainment rate proposed by Gregory (2001), Del Genio et al. (2012) efficiently reproduced the MJO transition from shallow to deep convection.

The advantage of the formulation of entrainment and detrainment rates in the unified scheme of Hohenegger and Bretherton (2011) is that it does not require an explicit distinction between deep and shallow convection. This formulation linking the fractional mixing rate ε0 to the convective precipitation at cloud base improved the simulation of the precipitation diurnal cycle compared to CAM, as well as relative humidity, cloud cover, and mass flux profiles, and could realistically simulate the transition between shallow and deep convection. Willet and Whitall (2017) also achieved a more realistic representation of the diurnal cycle in the tropics with this fractional mixing rate in their parameterization of entrainment in the UK Met Office model.

Other studies have evaluated the impact of entrainment–detrainment formulation on large-scale features, such as the double ITCZ (e.g., Chikira, 2010; Chikira and Sugiyama, 2010; Möbis and Stevens, 2012; Oueslati and Bellon, 2013). Möbis and Stevens (2012) used both the Tiedtke and Nordeng schemes in an aquaplanet GCM to evaluate the sensitivity of ITCZ to the choice of the convective parameterization. The Tiedtke scheme produced a double ITCZ, while the Nordeng scheme, with a higher lateral entrainment rate, led to a single ITCZ. In the works by Chikira (2010) and Chikira and Sugiyama (2010), the entrainment rate from AS was replaced by a formulation that depends on the surrounding environment following Gregory (2001) and Neggers et al. (2002). With this new formulation, variability and climatology improved, including the double ITCZ and the South Pacific Convergence Zone (SPCZ). Oueslati and Bellon (2013) obtained similar improvements in their study of the effects of entrainment on ITCZ by increasing entrainment in a hierarchy of models (coupled ocean–atmosphere GCM, atmospheric GCM, and aquaplanet GCM) at the cost of an overestimation of precipitation in the center of convergence zones. The role of entrainment in large-scale features was also underlined by Hirota et al. (2014) in their comparison of four atmospheric models with different entrainment formulations over tropical oceans.

Based on Zhang (2002) and using sounding data from the Coupled Ocean–Atmosphere Response Experiment (COARE), the South Pacific Convergence Zone (SGP97), and the Tropical Warm Pool – International Cloud Experiment (TWP-ICE), Zhang (2009) concluded that the entrainment of environmental air also affects CAPE and closure assumptions in CPs. The drier the entrained air, the stronger the dilution effect that acts to reduce CAPE. Moreover, dilute CAPE shows a better correlation with the consumption of CAPE than undilute CAPE.

As for the impact of entrainment and detrainment formulations for shallow convection, Siebesma and Holtslag (1996) evaluated a mass flux shallow cumulus based on BOMEX results and found that lateral entrainment and detrainment rates were 1 order of magnitude larger than those used in the Tiedtke scheme. Neggers et al. (2002) evaluated their multi-parcel model with LES results based on BOMEX and the Small Cumulus Microphysics Study (SCMS). The model reproduced the features of the buoyant part of the clouds and the variability of temperature, moisture, and velocity observed in cumulus clouds. Romps and Kuang (2010b) found that their stochastic formulation of entrainment reproduces the variability observed in the CRM well even when the cloud-base variability is turned off. While the convective updrafts simulated with the approach proposed by Sušelj et al. (2012) did not reach high enough compared to LES results and observations, the stochastic entrainment formulation described in Sušelj et al. (2013) properly simulated shallow cumulus, including the height of the updrafts and their reduction of horizontal area with height.

As mentioned in Sect. 4.2.2, less attention has been paid to the parameterizations of the detrainment process. Based on LES results for shallow convection, De Rooy and Siebesma (2008) proposed a new detrainment parameterization that led to improvements for ARM, BOMEX, and RICO shallow convection cases compared to the standard parameterizations of entrainment and detrainment (Siebesma and Cuijpers, 1995; Siebesma et al., 2003). Moreover, the authors revealed a greater variation in the detrainment rates from hour to hour and case to case than the variation in the entrainment rates. Derbyshire et al. (2011) confirmed this finding using a CRM and an adaptive detrainment proportional to the environmental relative humidity. Later, De Rooy and Siebesma (2010) showed that detrainment strongly influences the vertical structure of the mass flux.

4.3 Microphysics in convective clouds

The representation of microphysical processes in cumulus parameterizations is key to simulations of climate change (e.g., Ramanathan and Collins, 1991; Rennó et al., 1994; Lindzen et al., 2001). Convective microphysics greatly affect the representation of convective clouds due to their influence on detrainment of water vapor and hydrometeors, as well as the interaction between clouds and aerosols (e.g., Khain et al., 2005; Koren et al., 2005; Rosenfeld et al., 2008; Song and Zhang, 2011; Song et al., 2012; Tao et al., 2012). However, many convective parameterization schemes treat microphysical processes crudely, specifying an empirically determined conversion rate from cloud water to rainwater (e.g., Arakawa and Schubert, 1974; Tiedtke, 1989; Zhang and McFarlane, 1995; Han and Pan, 2011) or a certain precipitation efficiency, defined as the fraction of condensed cloud water converted to precipitation (Emanuel, 1991). The reader should keep in mind that other authors also take into account the effect of precipitation evaporation, and thus precipitation efficiency is defined as the fraction of condensate that reaches the surface (see Table 9). This is used in the calculations of the initial downdraft mass flux like in Bechtold et al. (2001). A brief description of the main assumptions and empirical values used in the representation of microphysics in CPs is presented here for the sake of completeness. For a detailed review of microphysics parameterizations, the reader is referred to Zhang and Song (2016) for convection and Tapiador et al. (2019b) for a full account.

Table 9A sample of empirical values and assumptions used in precipitation efficiency accounting for evaporation.

Download Print Version | Download XLSX

4.3.1 Conversion of cloud water to precipitation

Despite the importance of microphysical processes in the simulation of surface precipitation, radiation, or cloud cover, only a few convection schemes attempted to realistically represent these processes. A common approach is to assume that a specified fraction of the condensate is instantaneously removed as rain. In Yanai et al. (1973) and Tiedtke (1989), the conversion rate from cloud water to rainwater is assumed to be proportional to cloud water mixing ratio lw with an empirical function K(z) conversion coefficient that depends on height, as shown in Table 10. Other assumptions include a constant conversion coefficient Cc (Arakawa and Schubert, 1974; Grell, 1993; Zhang and McFarlane, 1995) or define a temperature-dependent threshold water content lwc, above which all cloud water is converted to precipitation (Emanuel and Živković-Rothman, 1999). Park and Bretherton (2009) modified the shallow cumulus parameterization described in Hack (1994) and used it in the UW scheme based on the shallow convective parameterization of Bretherton et al. (2004). Among the modifications introduced, cloud condensate exceeding a certain threshold value of the cloud condensate mixing ratio is converted into precipitation and includes the evaporation of convective precipitation above cloud base. In general, shallow convective schemes do not include a parameterization of conversion to precipitation.

Few schemes with a more realistic treatment of the conversion of cloud water to rainwater can be found in the literature on convection. Autoconversion of cloud water in the convection scheme is considered in Sud and Walker (1999), following Sundqvist (1978), as well as in Zhang et al. (2005). The latter included the autoconversion of cloud water and other microphysical processes for both cloud water and ice in the Tiedtke scheme. However, neither the size nor the number concentration of both hydrometeors is considered explicitly. This makes it impossible to account for aerosol–convection interaction, which is of great importance in climate simulations. To overcome this shortcoming, Song and Zhang (2011) and Song et al. (2012) added the mass mixing ratio and number concentration of each hydrometeor in their parameterization. Another more realistic treatment of condensation is that proposed by Bony and Emanuel (2001). In this scheme, the condensed water produced at the subgrid scale is predicted by the convection scheme, while its spatial distribution is predicted by a statistical cloud scheme through a probability distribution function of the total water. Indeed, the parameterization of the microphysics is more comprehensively devoted to this specific problem.

Table 10A sample of empirical values and assumptions used in the conversion of cloud water to precipitation (note: subscript “sh” refers to shallow convection).

Download Print Version | Download XLSX

4.3.2 Evaporation in downdrafts

Downdrafts are greatly affected by evaporation of hydrometeors and detrained cloud droplets due to latent cooling. Therefore, a realistic representation of this microphysical process is needed. However, only a limited number of convective parameterizations, such as Emanuel (1991), include an explicit calculation of this process, as shown in Table 11. Instead, crude assumptions can be found in the literature. The evaporation in downdrafts is often implicitly computed by assuming that the evaporation maintains a saturated or quasi-saturated downdraft while the equivalent potential temperature is conserved (e.g., Fritsch and Chappell, 1980; Zhang and McFarlane, 1995). More sophisticated formulations include those of Kreitzberg and Perkey (1976) based on Kessler (1969) and Song and Zhang (2011) based on Sundqvist (1988).

Table 11A sample of empirical values and assumptions used in the evaporation in the downdraft.

Download Print Version | Download XLSX

4.3.3 Aerosols

Aerosols play a key role in the climate system due to their influence on the Earth's energy budget through absorption and scattering of solar radiation. Focused on microphysical processes, aerosols serve as cloud condensation nuclei (CCN) and ice nuclei (IN) and thus affect cloud properties, dynamics, and precipitation. However, aerosol–convection interactions are very complex processes seldom included in convection microphysics. Zhang et al. (2005) developed a new parameterization accounting for the effects of aerosols in stratiform and convective clouds. This was later modified by Lohmann (2008) to include droplet activation by aerosols in terms of the updraft velocity w, temperature, aerosol number concentration, and size distribution, while ice nucleation is a function of w, aerosol properties, and air temperature. More recently, Grell and Freitas (2014) developed a new convective parameterization that includes an interaction with aerosols through an autoconversion of cloud water to rainwater dependent on CCN, parameterized in terms of the aerosol optical thickness (AOT) at 550 nm, as well as an aerosol-dependent evaporation of cloud drops. The authors also included tracer transport and wet scavenging in their parameterization. This convection scheme is currently available in WRF.

5 Closure: strategies to close the budget equation

Closure consists of defining the intensity or strength of convection, i.e., the amount of convection regulated by large-scale variables. Therefore, it is essential to close the budget equations (Eq. 5a, 5b and 5c). Despite the number of hypotheses proposed in the literature, it is still considered an unresolved problem (Yano et al., 2013). The following subsections discuss the main closure types, as well as their main assumptions and empirical values. The impact of the closure formulation in convective models concludes the section.

5.1 Closure types

Existing convective closures can be classified into diagnostic, prognostic, and stochastic. While diagnostic closures relate cumulus effects to the large-scale dynamics at a particular timescale, prognostic closures perform a time integration of explicitly formulated transient processes. Stochastic closures include randomness elements in closure schemes.

5.1.1 Diagnostic closures

Diagnostic closures include different types of closures based on a certain physical variable that expresses the intensity of convection. Table 12 shows a sample of empirical values and assumptions used in the closure in the updraft. In moisture convergence schemes, moisture convergence or vertical advection of moisture is selected as the closure variable (e.g., Kuo, 1974; Anthes, 1977; Krishnamurti et al., 1980, 1983; Kuo and Anthes, 1984; Molinari and Corsetti, 1985; Tiedtke, 1989), therefore assuming that convection consumes the moisture supplied by the large-scale processes.

The first parameterizations based on moisture convergence were too crude to produce results similar to those observed in nature, which led to the formulation of mass flux schemes. Early parameterizations lacked a theoretical framework to explain the interactions between the large-scale dynamics and convection or were incomplete, such as in Ooyama (1971). In an attempt to overcome this drawback, Arakawa and Schubert (1974) proposed a closed theory based on the QE of the CWF, which is similar to CAPE. Since then, many CPs have used CAPE-like closures, generally assuming that the adjustment occurs at a relaxed timescale in contrast to the instantaneous adjustment proposed in Arakawa (1969), among others. Table 13 lists the most important choices made for the relaxation timescale.

Following Lin et al. (2015), CAPE-like closures can be classified into two types according to the decomposition and constraints applied to the closure variable: the flux type and the state type. In the flux type, the change in the CAPE-like variable is decomposed into its large-scale and convective components. Of these types of closures, CAPE is the most commonly used closure variable in CPs (Fritsch and Chappell, 1980; Kain and Fritsch, 1993; Zhang and McFarlane, 1995; Gregory et al., 2000; Bechtold et al., 2001) with adjustment timescales varying from constant values to functional forms (Bechtold et al., 2008). Other schemes with CAPE closure include the KF scheme in WRF (Kain, 2004), as well as in CAM (Neale et al., 2008; Wang and Zhang, 2013), CAM6, and the Met Office Unified Model Global Atmosphere 7.0 (GA7.0) (Walters et al., 2019) for deep convection schemes. While the preceding schemes applied convective closure to the full troposphere, Emanuel (1995) and Raymond (1995) proposed the so-called boundary layer QE, wherein only the boundary layer component of the CAPE closure is considered. On the other hand, Zhang (2002) introduced a modified version of the QE assumption, in which only dCAPE is employed as the closure variable, without considering the effect of boundary layer forcing. This type of closure, known as the free tropospheric QE or the parcel-environment QE, provides a better simulation of the diurnal cycle of precipitation than the boundary layer QE (Zhang, 2003a), as well as a better representation of MJO and ITCZ than the QE assumption used in the Zhang–McFarlane scheme (Zhang and Mu, 2005b; Zhang and Wang, 2006; Song and Zhang, 2009; Zhang and Song, 2010). Donner and Phillips (2003) confirmed these results in their analysis over oceanic tropical areas and the midlatitude continental location of ARM. More recently, Bechtold et al. (2014) used the QE assumption to formulate a closure for the free troposphere based on boundary layer forcing. The dCAPE closure variable was replaced by PCAPE, defined as the integral over pressure of the buoyancy of an entraining ascending parcel with density scaling. The authors defined a convective adjustment timescale following Bechtold et al. (2008). This adjustment time is defined as the product of a convective turnover timescale τc and empirical scaling function f(n) that decreases with increasing spectral truncation. At the same time, τc is given by the ratio of the convective cloud depth and the vertical averaged updraft velocity. The authors stressed the dependency of τc with PCAPE through the velocity, which agrees with the observations in Zimmer et al. (2011). The implementation of this closure in the ECMWF IFS led to a better representation of the diurnal cycle of precipitation.

Table 12A sample of empirical values and assumptions used in the closure in the updraft.

Download Print Version | Download XLSX

Table 13A sample of the empirical values and assumptions in the relaxation timescale (note: subscript “sh” refers to shallow convection).

Download Print Version | Download XLSX

In contrast to the previous flux-type closures, state-type closures decompose the change in the CAPE-like variable into its boundary layer component and free troposphere component instead of its large-scale and convective component. The main representatives of state-type closures are the convective adjustment schemes of Betts (1986), wherein mesoscale and subgrid-scale cloud processes maintain QE, and Emanuel (1994), wherein QE is related to fluctuations of entropy in the sub-cloud layer. Differences between these adjustment schemes are in the adjustment timescale and reference profiles selected for the adjustment. For example, Emanuel (1994) included an adjustment timescale for the sub-cloud layer of the order of half a day, while Betts and Miller (1986) found good results for values between 1 and 2 h based on GATE wave data. More recently, authors such as Khouider and Majda (2006, 2008) and Kuang (2008) applied a state-type scheme only to the lower troposphere.

An alternative principle to QE is the so-called activation control proposed by Mapes (1997), in which the intensity of deep convection is controlled by inhibition and initiation processes at low levels, and closure is formulated in terms of CIN and the turbulent kinetic energy (TKE) (Mapes, 2000; Fletcher and Bretherton, 2010). However, as highlighted in Yano and Plant (2012b), this formulation is not self-consistent, which is a must, as models are intended to test physical hypotheses (the reader is referred to Yano et al., 2013, for a detailed explanation). In Rio et al. (2009) the intensity of convection is controlled by sub-cloud processes, such as boundary layer thermals. The authors defined the closure in terms of the so-called available lifting power (ALP), which is the flux of kinetic energy associated with thermals. Grandpeix and Lafore (2010) also used an ALP closure in their wake parameterization for GCMs coupled with Emanuel's scheme (Emanuel, 1991), as did Hourdin et al. (2013) in the development of the LMDZ5B. While in Grandpeix and Lafore (2010) the source of ALP comes from the collapse of the wakes, in Hourdin et al. (2013) the thermal plumes and the spread of cold pools are the ones providing the power.

This section presented the assumptions and empirical values used in the formulation of the closure for updrafts. However, the magnitude of the downdrafts should also be addressed. In the schemes in which it is included, it is commonly expressed as a fraction γd of the closure of the corresponding updraft, setting γd as a certain value (Johnson, 1976; Tiedtke, 1989; Baba, 2019). Alternatively, other authors have related γd to precipitation efficiency (Emanuel, 1995; Bechtold et al., 2001), the RH in the LFS (Kain, 2004), or proposed a formula for γd in terms of the total precipitation rate within the updraft (Zhang and McFarlane, 1995). Table 14 lists some of the empirical values and assumptions used in closure in the downdraft.

Table 14A sample of empirical values and assumptions used in the closure in the downdraft.

Download Print Version | Download XLSX

The discussion above focused on closure in deep convective and unified schemes. As for shallow convection closures, different approaches have been proposed since the publication of the first convection schemes. In this paper, we present a framework for the main empirical values and assumptions for shallow convection following the classification in Neggers et al. (2004). The authors classified the main shallow convection closures into moist static energy convergence, CAPE adjustment, and sub-cloud convective velocity scaling.

In the moist static energy closures, the QE budget for moist static energy controls shallow convection activity. Based on the results obtained by LeMone and Pennell (1976) from trade wind cumuli and the moisture convergence hypothesis from Kuo (1965, 1974) and (Lindzen, 1988), Tiedtke (1989) proposed a shallow convection closure based on the moist static energy closure. Later, Raymond (1995) and Emanuel (1995) used it in the boundary layer quasi-equilibrium for shallow convection, and Gregory et al. (2000) included it in a revised version of the ECMWF scheme. More recently, Bechtold et al. (2014) parameterized the mass flux for shallow convection in terms of the vertically integrated moist static energy tendency.

Other authors proposed shallow convection closures based on the relaxation of the system towards a certain reference state within a relaxation timescale, i.e., adjustment scheme. For example, Albrecht et al. (1979) used this closure in their study of the trade wind boundary layer, specifying a constant adjustment time set to 1/3 d according to the observation results obtained by Betts (1975) for BOMEX. Later, based on observations from BOMEX and ATEX, Betts (1986) used an adjustment scheme for shallow convection in which the thermodynamic structure tends towards a mixing line with an adjustment time set to 3 h. Bechtold et al. (2001) used the same value for the relaxation time in their CAPE closure formulation for shallow convection.

One of the main representatives of TKE budget closures is in Grant (2001), who assumed that mass flux at cloud base is proportional to the convective velocity scale proposed by Deardorff et al. (1969), w. The proportionality constant is the area fraction of cumulus updrafts and was determined by plotting the cloud-base mass flux versus the sub-cloud-layer velocity scale in LES (see Table 15). This shallow closure was further used by other authors, such as Soares et al. (2004), Siebesma et al. (2007), and Pergaud et al. (2009) in an EDMF, and Han and Pan (2011) and Han et al. (2017) in their revision of the NCEP GFS, among others. While Soares et al. (2004) defined the mass flux as the product of the updraft vertical velocity and a constant updraft fraction, Siebesma et al. (2007) scaled the mass flux with the standard vertical velocity deviation and set the proportionality constant to 0.3. In Pergaud et al. (2009) the closure is based on the mass flux near the surface instead of at the LCL. The authors set the proportionality constant to 0.065 based on LES results. Han et al. (2017) modified the closure by making the cloud-base mass flux a function of the mean updraft velocity. This way, shallow convection can be triggered in the stable boundary layer. Another closure based on the relationship between mass flux and TKE is that described in Kain (2004), wherein the mass flux is scaled with the maximum TKE in the sub-cloud layer. The convective time period in this parameterization ranges from 1800 to 3600 s.

Similar to these parameterizations, Hourdin et al. (2002) developed a new mass parameterization of vertical transport in the convective boundary layer, known as the thermal plume model, wherein the closure depends on the maximum vertical velocity and an area fraction. As stated in Rio and Hourdin (2008) the area fraction is predicted according to the entrainment and detrainment in contrast to the constant values used in Soares et al. (2004) and Siebesma et al. (2007), among others.

Using LESs and observations, Grant and Lock (2004) proposed a shallow convective closure proportional to CAPE and the convective velocity scale w. More recently, Zheng et al. (2016) extended the shallow convection study of Grant and Lock (2004) and expressed the closure in terms of CAPE and cloud-depth-averaged vertical velocity.

In the DualM framework, Neggers et al. (2009) defined the vertical structure of the updraft mass flux as the product of the updraft vertical velocity and updraft fraction. Based on results from De Rooy and Siebesma (2008) and the statistical distribution type in Sommeria and Deardorff (1977), the authors used a moist zero buoyancy deficit to estimate the updraft area fraction and, through it, the vertical velocity and mass flux.

A different shallow convection closure was suggested by Mapes (2000). The author expressed the mass flux in terms of CIN and TKE. Later, Bretherton et al. (2004) developed a new parameterization consisting of coupling a PBL turbulence model based on Grenier and Bretherton (2001) with a shallow convective mass flux scheme based on an entraining–detraining single-plume model. The closure assumes that a buoyant cumulus cloud can form if the vertical velocity of source air is high enough to penetrate the inversion layer in the sub-cloud layer and reach its LFC. The critical velocity is a function of CIN and the distribution of velocities is assumed to be Gaussian. The mass flux closure has a form similar to that proposed by Mapes (2000). In this case, it is an exponential function of the ratio between CIN and the average TKE in the sub-cloud layer calculated by the PBL scheme. In their simulations of the transition from shallow to deep convection, Kuang and Bretherton (2006) applied the CIN-based closure proposed by Mapes (2000) with the updraft velocity at cloud base set to the sub-cloud-layer TKE as in Bretherton et al. (2004). In the unified scheme of Hohenegger and Bretherton (2011), the shallow closure is a function of the ratio between CIN and mean planetary boundary layer TKE. Despite its use in several convection schemes, this parameterization is not self-consistent as already mentioned in Sect. 5.1.1.

In the MM5, Deng et al. (2003) proposed three different shallow convection closures depending on the values of the cloud depth (CD), cloud-top height zt, and LFC height zLFC and assumed a uniform updraft geometry. The closures include a TKE-based closure, a CAPE closure, and a hybrid closure. TKE-based closure is used when ztzLFC. In this closure, the cloud-base mass flux in the sub-cloud layer scales with the maximum diagnosed TKE in the sub-cloud mass-source layer over a relaxation timescale. If CD≥4 km, the CAPE closure applies, while for CD<4 km and zt>zLFC a hybrid closure between TKE and CAPE closures is used. The transition is done through a simple linear averaging. More recently, Freitas et al. (2020) proposed a trimodal formulation instead of the unimodal deep plume used in Grell and Freitas (2014) to represent shallow, congestus, and deep convection. Closures for shallow convection include the boundary layer quasi-equilibrium from Raymond (1995), the closure proposed in Grant (2001), and a closure based on the heat engine treatment of convection applied in Rennó et al. (1994). This closure relates the updraft cloud-base mass flux to the buoyancy surface flux, a certain thermodynamic efficiency, and the total CAPE that is equivalent to the standard CAPE.

5.1.2 Prognostic closures

Compared to the QE assumption used in the majority of the diagnostic closures mentioned above, prognostic closures do not distinguish between large-scale and convective processes and substitute the QE assumption with time integration of prognostic equations. These equations explicitly account for the time changes of different physical variables, i.e., convective kinetic energy or h, which are related to the cloud-base mass flux through a dimensional parameter. Energy dissipation rate is also included in this type of closure through a dissipation term, either determined by a second dimensional parameter called dissipation time (e.g., Randall and Pan, 1993; Pan and Randall, 1998; Yano and Plant, 2012a) or expressed in terms of the entrainment rate and an aerodynamic friction coefficient (e.g., Gerard and Geleyn, 2005). Gerard and Geleyn (2005) defined cloud-base mass flux as Mu=-auwu, where au is a prognostic updraft fraction area obtained by a moist static energy closure, and wu is a prognostic vertical updraft velocity. Gerard (2007) and Gerard et al. (2009) also used this approach and even applied it for downdrafts (Gerard et al., 2009). Other schemes using prognostic updraft fractional areas include those of Grandpeix and Lafore (2010), Mapes and Neale (2011), and Tan et al. (2018), among others (see Table 15).

Table 15A sample of empirical values and assumptions used in the cloud fraction (note: subscript “sh” refers to shallow convection).

Download Print Version | Download XLSX

5.1.3 Stochastic closures

Usually subgrid-scale processes are considered in an ensemble mean sense in CPs (Lin and Neelin, 2000, 2002). Stochastic closures include randomness elements in convective scheme closures to represent these subgrid-scale processes in a more realistic way. Numerous stochastic convective parameterizations have been proposed (e.g., Lin and Neelin, 2000, 2002; Majda and Khouider, 2002; Lin and Neelin, 2003; Khouider et al., 2003; Khouider, 2014). However, as Stechmann and Neelin (2011) stated, sometimes the distinction between stochastic triggers and stochastic closures is not clear. Differences between the proposed closures are in the type of stochastic process employed. For instance, Stechmann and Neelin (2011) proposed a stochastic closure for precipitation using Gaussian white noise, while Majda and Khouider (2002) and Khouider et al. (2003) used a Markov jump process.

For deep convection, Lin and Neelin (2000) include a first-order autoregressive random noise component in the convective parameterization of Betts and Miller (1986), keeping the convective relaxation timescale. This random noise is expressed as ξt=cξξt-1+zt, where cξ is an autoregressive coefficient that yields an autocorrelation time τξ for the process, and zt is white noise with zero mean and standard deviation σz. The authors evaluated three values for τξ, i.e., 20 min, 2 h, and 1 d, with three different σz, i.e., 4.5, 0.8, and 0.1 K, respectively. Longer τξ produced better results compared to observations. Lin and Neelin (2003) introduced this stochastic component in the ZM closure with τξ= 1 d and σz=1000Jkg-1. This scheme increased precipitation variance toward observations. Based on the variability around the equilibrium state, Plant and Craig (2008) and Groenemeijer and Craig (2012) used a PDF to obtain random values for the cloud-base mass flux. This PDF expresses the chance of launching a cloud with a certain radius between two calls of the convective scheme. The radius is assumed to be related to the mass flux. It is defined as pmdm=1mexp-mmdm, where m is the mass flux per cloud and m is its ensemble average, both related through the definition of updraft radius m=mR2R2. Moreover, the closure timescale in Plant and Craig (2008) is defined as τc=kL=kmM, where M is the ensemble mean total cloud-base mass flux calculated as in Kain and Fritsch (1990), and k is a constant that depends on the definition of adjustment. The default parameter choices in Plant and Craig (2008) are m=2×10-7kgs-1, a root mean squared cloud radius of R21/2=450 m, and k=0.3sm-1. In Groenemeijer and Craig (2012) these values did not produce enough convective, so they were changed to m=1×10-7kgs-1 and R21/2=1200 m, with fixed τc=600 s. Bengtsson et al. (2013) introduced a CA in the parameterization of the updraft mesh fraction au used in the Gerard et al. (2009) cumulus convective scheme closure. Using observational data, Dorrestijn et al. (2015) determined the au for various cloud types using Markov chains. The one for deep convection was later implemented in the Tiedtke cumulus scheme in Simplified Parameterizations, Primitive Equation Dynamics (SPEEDY).

For shallow convection, Sakradzija et al. (2015) developed a stochastic shallow parameterization following the studies of Craig and Cohen (2006) and Plant and Craig (2008) for deep convection. In this scheme, the number of new clouds is a sample from a Poisson distribution, while the lifetime average mass flux for each new cloud is randomly sampled from a Weibull distribution with two modes, namely forced and passive clouds on one hand and active clouds on the other. This Weibull distribution is defined through a scale λ and a shape k parameter. The cloud lifetime is defined as τclt=αimβi, where the coefficients are obtained from the nonlinear least square fitting of the joint distribution of cloud mass flux and cloud lifetime. The total cloud-base mass flux is then calculated by integrating the instantaneous mass flux distribution, i.e., M=0mτclt(m)Gp(m)dm or M=Gαλk+1Γ2+1k, where G is the cloud-generating rate. The following values were used for this parameterization: k=0.7, λ1=7269.08kgs-1, λ1=29868.48kgs-1, α1=0.02kg-1, α2=0.33kg-1, and G=4.55s-1 (subscript 1 refers to forced and passive clouds, and subscript 2 is for active clouds). The reader is referred to Sakradzija et al. (2015) for values of other parameters. This scheme was later implemented in EDMF (Sakradzija et al., 2016) and ICON (Sakradzija and Klocke, 2018) with variations in the values of the aforementioned parameters.

5.2 Impact of closure on convective models

The closure problem is one of the major challenges in CPs. As well as being essential to close the budget equations (Eq. 5a, 5b and 5c), it plays an important role in the performance of CPs. For instance, Bechtold et al. (2008) obtained a better representation of the rainfall pattern and tropical wave activity with their modifications of the entrainment and convective adjustment time in the deep convection scheme in IFS. In Rio et al. (2009), the representation of the diurnal cycle of precipitation is greatly improved using the ALP deep closure in a 1D model. In their formulation, the convective mass flux scheme is coupled with cold pools and the thermal plume model through the ALP. Using a dilute CAPE closure together with convective momentum transport, Neale et al. (2008) improved the representation of ENSO in CAM3. Adding a stochastic component to the deep convection closure in BMJ, Lin and Neelin (2000) obtained a better representation of the intraseasonal variability. Later, Lin and Neelin (2003) included a stochastic component in the deep closure of the ZM scheme. The daily variance was much closer to observations than without the stochastic component. Moreover, the SPCZ was better placed.

Replacing the CAPE closure used in the ZM scheme with a dCAPE closure, Zhang (2002) improved the simulation of precipitation, moisture, and temperature for midlatitude continental convection. This closure also improved the diurnal cycle of precipitation over the Southern Great Plains in the US (Zhang, 2003b). The replacement of the ZM closure by dCAPE provided a better representation of the tropical precipitation in NCAR CCM in Zhang and Mu (2005a). With this closure, the precipitation was enhanced over the western Pacific monsoon region during June, July, and August, as well as the SPCZ during December, January, and February. In the representation of the MJO, Zhang and Mu (2005b) used the closure and convection trigger proposed in Zhang and Mu (2005a) and removed the restriction in the convection originating level. The simulated MJO was more consistent with the observations in terms of variability in precipitation, outgoing longwave radiation, and zonal wind, and it exhibited a clear eastward propagation. However, the precipitation signal and the time period of the MJO differ from the observations. This revision of the ZM scheme used in the NCAR Community Climate System Model (CCSM3) also alleviates the biases related to the double ITCZ in precipitation and the cold tongue in sea surface temperature (SST) over the Equator, among other benefits (Zhang and Wang, 2006; Song and Zhang, 2009; Zhang and Song, 2010). Wang and Zhang (2013) evaluated three different trigger and closure assumptions in CAM4 and CAM5 and highlighted the need to use multiple independent observations simultaneously to constrain models to reduce the degrees of freedom as well as the need to avoid the individual treatment of model physical parameterizations. Wang et al. (2016) obtained a better representation of the precipitation intensity, especially over the tropical belt, as well as improved simulations of the eastward-propagating intraseasonal signals of precipitation and zonal wind by coupling the Plant and Craig (2008) stochastic parameterization with the ZM scheme in CAM5. More recently, Becker et al. (2021) showed a better representation of the propagation and organization of mesoscale convective systems, such as African squall lines, when adding a term for the integrated and scaled total advective moisture tendency to the CAPE closure.

Using CRM simulations, Kuang and Bretherton (2006) tested the viability of representing the transition from shallow to deep convection using a CIN-based closure similar to the shallow closure in Bretherton et al. (2004). Results from an idealized numerical experiment of shallow-to-deep convection transition are in agreement with the CIN-based closure and do not support a closure based solely on CAPE. Later, Fletcher and Bretherton (2010) extended the Bretherton et al. (2004) shallow closure to deep convection with the goal of finding a closure that works well for both shallow and deep convection without changing any parameter. Three CRM simulations forced with observations from ARM Great Plains, Kwajalein Experiment (KWAJEX), and BOMEX were used to test this closure as were a CAPE and a Grant closure (Grant, 2001). The CIN-based closure was more skillful in the prediction of the cloud-base mass flux and performed well for both deep and shallow convection. Hohenegger and Bretherton (2011) modified the UW shallow convection scheme to develop a unified scheme for shallow and deep convection. The closure introduced also relates the cloud-base mass flux to TKE and CIN, taking into account the contribution of cold pools to the increase in TKE. LESs as well as BOMEX, KWAJEX, and ARM were used to formulate and improve this parameterization. Tested in the Single-column Community Atmosphere Model (SCAM) single-column modeling framework, this parameterization was able to represent both shallow and deep convection as well as midlatitude continental convection. Han and Pan (2011) modified the deep scheme in SAS (Pan and Wu, 1995) by increasing the allowable cloud-base mass flux, which was originally set to 0.1kg(m2s-1)-1, with a Courant–Friedrichs–Lewy (CFL) criterion to make cumulus deeper and stronger. This scheme effectively eliminated the remaining instability in the atmospheric column that was producing excessive grid-scale precipitation in the original formulation. Using a PCAPE closure with boundary layer forcing, the scheme for shallow and deep convection described in Bechtold et al. (2014) represented the observed daytime evolution of convection over land fairly well when compared with observations such as satellite data. Moreover, the evolution of shallow and deep convection agreed with CRM results. Over Europe, it better represented the mainly surface-driven convection over the Balkans and the Atlas Mountains, as well as forced convection over central Europe, and reduced unrealistic rates of snowfall along the coast of the British Isles and near the European continent for a particular winter case. Han et al. (2020) obtained similar results using this closure in KIM (the Netherlands Institute for Transport Policy Analysis). The afternoon peak was delayed, and the biases of the overestimated precipitation over land in the morning and late afternoon were reduced.

Focused on closures for shallow convection, different authors have analyzed the impact that shallow convection closures have on the simulation of the diurnal cycle. For instance, Neggers et al. (2004) evaluated moist static energy closure, CAPE adjustment, and sub-cloud convective velocity scaling closure against LESs and analyzed the impact of each closure on the simulation of the diurnal cycle. Among those, the sub-cloud convective velocity scaling closure showed the best results. The onset, dissipation time, and cloud cover of cumulus clouds were well captured by the EDMF scheme in Soares et al. (2004). Scaling the mass flux with the standard vertical velocity deviation in the EDMF, Siebesma et al. (2007) obtained realistic representation of the main properties of dry convective boundary layers. Using a similar closure, Pergaud et al. (2009) showed the ability of the EDMF scheme to represent mixing in the counter-gradient zone and to handle the diurnal cycle of boundary layer cumulus clouds. Similar results were obtained by Rio and Hourdin (2008) in terms of the diurnal cycle of the boundary layer. The shallow cumulus parameterization developed by Bretherton et al. (2004) reproduced LES results obtained by Siebesma and Cuijpers (1995) and Siebesma et al. (2003) well for a subperiod of BOMEX, as well as those by Wyant et al. (1997) for the transition from stratocumulus to trade. However, this transition was slightly abrupt in the simulations with the shallow parameterization. McCaa and Bretherton (2004) further analyzed the performance of this scheme in a regional climate simulation of the subtropical northeastern Pacific Ocean in MM5. The regional mean shortwave cloud radiative forcing and vertical structure were better represented by this scheme compared to other parameterizations of cloud-topped boundary layer processes. In the DualM framework, Neggers et al. (2009) defined the cloud-base mass flux as the product of updraft fraction and updraft vertical velocity. Examined for ATEX, this closure produced steeper gradients closer to LES results than the ones obtained with a fixed structure of the mass flux, and it was concluded that this result is an indicator of the interaction between the mass flux and environmental humidity introduced by the closure. Han and Pan (2011) replaced the shallow convection in SAS with a new formulation using the shallow closure described in Grant (2001). Compared to the original formulation, this new scheme did not destroy stratocumulus clouds off the west coasts of South America and Africa.

6 Conclusions

Numerical models need simplifications in order to cope with the complexity of the physical processes actually occurring in the atmosphere. The degree of simplification in the physics is evolving at a pace inverse to the availability of computational power. Thus, early convective parameterizations (as well as parameterizations of radiation, turbulence, microphysics, etc.) were based on very simple assumptions, such as the conditional instability of the second kind (CISK) first presented by Charney and Eliassen (1964) and Ooyama (1964) in tropical cyclone modeling. Manabe et al. (1965) proposed a different parameterization, the so-called adjustment scheme, wherein atmospheric instability is removed through an adjustment towards a reference state. The instability was removed instantaneously, and condensed water precipitated immediately. However, the scheme produced very large precipitation rates and a saturated final state after convection, which is rarely observed in nature (Emanuel and Raymond, 1993). To alleviate these issues, relaxed adjustment schemes and penetrative adjustment schemes (Betts, 1986; Betts and Miller, 1986) were proposed. Such improvements were only possible when more powerful computers became available. However, novel theoretical approaches ahead of the technological capabilities of the time have also greatly impacted the field. Thus, the first parameterizations based on moisture convergence were too crude to produce results similar to those observed in nature, which led to the formulation of mass flux schemes. Simulations improved with further refinements of the interaction of cumulus clouds with the large-scale environment by, for instance, Ooyama (1971) (a statistical ensemble of bubbles represents cumulus convection) or Yanai et al. (1973) (detrainment and cumulus-induced subsidence). Early parameterizations lacked a theoretical framework to explain the interactions between the large-scale dynamics and convection or were incomplete, such as in Ooyama (1971). In an attempt to overcome this drawback, Arakawa and Schubert (1974) proposed a closed theory based on the cloud work function and adjustment towards QE. A few years after, thanks to the increase in computational power, more complex parameterizations and new variables based on observations were implemented to achieve better spatial and temporal resolutions. Krueger (1988) put forward the cloud-system-resolving model (CSRM) idea to explicitly simulate convective processes over a kilometer scale instead of using parameterizations. However, this approach entails an extremely high computational cost. As an alternative with a lower computational cost, a multiscale model framework (MMF) or super-parameterizations (SPs) emerged. In this case, convective parameterizations are replaced by 2D cloud-resolving models (CRMs), or even a 3D LES model, at each grid cell of a GCM (Grabowski and Smolarkiewicz, 1999).

To alleviate problems associated with traditional convective parameterizations, e.g., the representation of the diurnal cycle of convection (e.g., Yang and Slingo, 2001; Guichard et al., 2004), several studies introduced modifications to existing models. Challenges remain for convective parameterizations. As highlighted in Rio et al. (2019), three of these major challenges include (a) improving the representation of convective cloud ensembles, (b) improving the representation of convective memory and organization, and (c) improving the representation of convection to large-scale interactions. The reader is referred to Rio et al. (2019) for a comprehensive review. Here, only the main representatives of each challenge are mentioned.

Regarding the first challenge, current approaches to improve the representation of convective cloud ensembles include unified and multi-object framework parameterizations that account for the coexistence of more numerous cloud types within a model grid cell and different methods to compute the vertical profile of cloud properties. Traditionally, models have used separate parameterizations for shallow and deep convection. Guichard et al. (2004) stressed the necessity of using an ensemble of parameterizations that represents a succession of convective regimes. Some modelers proposed keeping shallow and deep convection parameterizations separate due to their different nature and then using a parameterization to couple them (e.g., Rio et al., 2013), while others proposed unified schemes that attempt to merge shallow and deep convection into one parameterization (e.g., Guérémy, 2011; Arakawa and Wu, 2013; Wu and Arakawa, 2014; Park, 2014a, b; D'Andrea et al., 2014; Kwon and Hong, 2017; Zhao et al., 2018). Also, models traditionally split the turbulence parameterization among the PBL and moist convection, simplifying the treatment of turbulence but requiring the addition of an artificial closure to match both schemes (Sušelj et al., 2014). Unified models have been also used to merge these parameterizations, such as the so-called Cloud Layers Unified By Binomials (CLUBB) (Golaz et al., 2002a, b; Larson et al., 2002). Two different approaches have been proposed that unify the PBL as well as shallow and deep convection. Those approaches are the so-called EDMF framework (e.g., Hourdin et al., 2002; Köhler et al., 2011; Hourdin et al., 2013; Bhattacharya et al., 2018) and third-order turbulent schemes (e.g., H. Guo et al., 2014, 2015). Parameterizations that account for the coexistence of more numerous cloud types within a model grid cell include the use of Markov chains considering a certain number of cloud types (Khouider et al., 2010; Dorrestijn et al., 2013b; Peters et al., 2013) or the use of a probability density function (PDF) (e.g., Plant and Craig, 2008; Sakradzija et al., 2016) , among others. As for the methods to compute the vertical profile of cloud properties, numerous studies apply a deterministic entrainment to different cloud types; others use stochastic entrainment parameterizations (e.g., Raymond and Blyth, 1986; Emanuel and Živković-Rothman, 1999; Grandpeix et al., 2004; Romps and Kuang, 2010b; Sušelj et al., 2013; Romps, 2016). The vertical profile of vertical velocity also needs further attention as many schemes do not solve an equation for the vertical velocity, and the ones that do it are mostly based on the equation proposed by Simpson and Wiggert (1969) as highlighted in de Roode et al. (2012).

For the second challenge, improving the representation of convective memory and organization, there are at least two outstanding issues. On the one hand, as pointed out in Davies et al. (2009), the QE hypothesis does not account for convective memory. Different strategies have been proposed to include it in convective parameterizations, such as the use of prognostic variables (e.g., Pan and Randall, 1998; Gerard and Geleyn, 2005; Piriou et al., 2007; Mapes and Neale, 2011; Hohenegger and Bretherton, 2011; Willet and Whitall, 2017; Tan et al., 2018), Markov chains (e.g., Khouider et al., 2010; Hagos et al., 2018), cellular automaton (CA) assigning a prescribed lifetime to each active cell (e.g., Bengtsson et al., 2011, 2013) or cold pools (e.g., Grandpeix and Lafore, 2010; Park, 2014a, b; Del Genio et al., 2015; Colin et al., 2019). On the other hand, as for the representation of convective organization, Donner (1993), Alexander and Cotton (1998), and Donner et al. (2001) represented the effects of mesoscale circulations and Mapes and Neale (2011) introduced a prognostic variable called organization that represents the degree of subgrid organization. Other studies accounting for convective organization use surface cold pools (e.g., Rio et al., 2009; Grandpeix and Lafore, 2010; Rochetin et al., 2014a, b; Park, 2014a, b; Böing, 2016), a slantwise overturning model (e.g., Moncrieff et al., 2017), CA (e.g., Shutts, 2005; Bengtsson et al., 2011, 2013, 2019, 2021), or PDF-based or spectral schemes based on a discretized distribution (e.g., Neggers et al., 2003; Wagner and Graf, 2010; Neggers, 2012; Park, 2014a, b; Neggers, 2015). Accurate representations of precipitation and cloud cover are important for the spatial organization and time evolution of convective systems. Parameterizations accounting for the microphysics of precipitation include those of Feingold (2003), Genio et al. (2005), McFiggans et al. (2006), and Heymsfield et al. (2013), among others. Several studies have also attempted to improve convective cloud radiative effects using PDFs (e.g., Bogenschutz et al., 2010; Perraud et al., 2011; Hourdin et al., 2013; Storer et al., 2015; Qin et al., 2018).

The third main challenge is to achieve better representations of convection to large-scale interactions, i.e., shallow convection and transitions from shallow to deep and from deep to organized convection. For transitions from shallow to deep, various approaches have been proposed, especially focused on the representation of the diurnal cycle of precipitation (e.g., Rio et al., 2009; Stratton and Stirling, 2012; Rio et al., 2013; Bechtold et al., 2014; Rochetin et al., 2014a; Peters et al., 2017). Other aspects that deserve more attention, among others, are the representation of the impact of sea breeze in deep convection initiation over islands and the tendency to show strong positive tropical rain biases for models with strong intraseasonal variability due to the sensitivity of convection to free tropospheric humidity through entrainment (Rio et al., 2019). Transitions from deep to organized convection also deserve more attention due to the role that mesoscale convective systems play in weather and climate.

The field of modeling convection is full of details and intricacies. As already mentioned, mass flux convective parameterization schemes are still the most common convective parameterizations used in ESMs, RCMs, and NWP models. Models have traditionally used separate parameterizations for shallow and deep convection Therefore, we mainly focused our attention to the assumptions and empirical values used in shallow and deep mass flux schemes for their three main elements, i.e., trigger, cloud model, and closure. In the activation of convection, the main differences between shallow and deep convection are in the cloud-depth criterion, the updraft radius, and the buoyancy threshold. Both cloud depth and radius are always set to smaller values compared to deep convection. As for the temperature perturbation that some deep convective parameterizations include in the buoyancy threshold, it is absent in shallow convection triggers. Commonly, the procedure followed to find cloud base and trigger convection is the same for both schemes, though some studies set different conditions for the USL (Han and Pan, 2011) or use a vertical velocity criterion to trigger shallow convection (Bretherton et al., 2004; Park and Bretherton, 2009). The cloud-depth criterion is what decides which type of convection activates.

Numerous parameterizations of entrainment and detrainment have been proposed for shallow and deep convection including turbulent and dynamical components (e.g., Tiedtke, 1989, and Nordeng, 1994, for deep and shallow convection), constant values (e.g., Song and Zhang, 2017, for deep and Siebesma, 1998, for shallow convection), inverse proportionality to height (e.g., Siebesma and Cuijpers, 1995, for deep and Jakob and Siebesma, 2003, for shallow convection) or to the vertical velocity of the parcel (e.g., Gregory, 2001, for both deep and shallow convection), or dependence on a critical mixing fraction (e.g., Kain and Fritsch, 1990, for deep and Bretherton et al., 2004, for shallow convection), among others. For schemes using the same parameterization for shallow and deep convection, the main difference between the two types is in the values, which are higher for shallow than for deep convection. Entrainment and detrainment formulations for downdrafts usually use similar parameterization as for updrafts. In terms of the microphysics, shallow convective schemes usually do not include a parameterization of conversion to precipitation.

As for the closure formulation, numerous deep convective schemes use CAPE-based closures, although formulations based on convective adjustment in terms of CIN and TKE or using stochastic closure have also been proposed. For shallow convection, the most often used are TKE-based closures. Other closures such as moist static energy convergence (Tiedtke, 1989) and CAPE adjustment closures (Betts, 1986) are also used in shallow convection. For the latter, the adjustment time is usually higher for shallow than for deep convection. In the parameterizations in which it is included, downdraft closure is commonly expressed as a fraction of the closure of the corresponding updraft.

Convective parameters require fine tuning, but there is no explicit methodology to do so. In some cases, authors use the variables that are easiest to measure. In others, mean values describe processes that cannot be modeled in sufficient detail, or the values represent particular conditions for certain locations and atmospheric events (Mauritsen et al., 2012). For instance, Bony and Emanuel (2001) adjusted their water vapor and temperature prediction using the TOGA-COARE data measured in the western Pacific Ocean in 1993, while Betts and Miller (1986) used GATE datasets measured over the tropical Atlantic Ocean in 1974 to develop their deep convection scheme. Hence, empirical values and assumptions selected this way might yield good results when compared to observations from certain locations and less good results for others. Commonly, manual tuning of convective parameters is used, although various automatic methods have recently been used to estimate parameters, including the variational method (Emanuel and Živković-Rothman, 1999), Bayesian calibration (e.g., Hararuk et al., 2014; Wu et al., 2018), simulated annealing method (e.g., Jackson et al., 2004, 2008; Liang et al., 2014), genetic algorithm (e.g., Lee et al., 2006), ensemble data assimilation (e.g., Ruiz et al., 2013; Li et al., 2018), and machine learning (e.g., Schneider et al., 2017), among others. Recently, Couvreux et al. (2021) proposed a new method that performs a multi-case comparison between SCM and LES results to calibrate parameterizations. The method uses machine learning without replacing parameterizations.

Comparisons with observations were, and still are, crucial to the development of convective parameterizations. For instance, the underprediction of large-scale precipitation by dry adiabatic models compared to observations led to the inclusion of moist adiabatic processes in NWP models (Smagorinsky, 1956), and lake-effect snow observations (Niziol et al., 1995) forced a reduction of the minimum cloud-depth threshold in Kain and Fritsch (1993) to 2 km. However, observations suffer from data gaps and the instruments used are not able to sample key variables in parametric equations. The use of observations by the convective modeling community has not been sufficient so far. The reasons are twofold. Basic convective quantities like mass flux and important parameters like adjustment timescale, entrainment, and microphysical parameters can often be only indirectly inferred from observations like infrared and microwave satellite data, radar data, rainfall rates, radiosonde networks, and reanalysis data. When we say that they are indirectly inferred, we mean that these quantities are adjusted to optimize the model fit to the observed radiative and surface fluxes as well as the observed temperature and wind field. On the other hand, long-term instrumentation deployment at meteorological supersites (e.g., Neggers et al., 2012; Song et al., 2013; Gustafson et al., 2020; Zheng et al., 2021) or dedicated convection field campaigns like GATE, TOGA-COARE, DYNAMO, PECAN (Geerts et al., 2017), and EUREC4A (Bony et al., 2017), to mention a few, have been conducted to quantify convection and its effect on the large-scale flow, and powerful LES data are available with statistical samples of the convective updraft and downdraft properties. However, the dilemma is that these data are only available locally or for specific setups; LES data also need to be constrained by observations and an accurate convection parameterization in a global model needs to be constrained globally.

Modern extensive big datasets such as those derived from COPERNICUS data are very relevant to constrain assumptions and calibrate parameterizations. Recently, Neggers et al. (2012) and Gustafson et al. (2020), among others, have provided a successful attempt to reconcile observations and LES data. This new approach consists of combining LES outputs with observations. Indeed, high-resolution models provide additional information in 4D that is not possible to obtain from point-based measurements (Gustafson et al., 2020). The complementary approach consists of new dedicated satellite missions such as INCUS or the follow-on to CloudSat and CALIPSO, which can provide global, homogeneous, and time-extended observations. Satellite estimates of the convective mass flux are becoming available (Jeyaratnam et al., 2021) and new missions are in the planning to fill the gap in global, multiple-regime observations of convection. Although observations have long been used to tune parameters in convective schemes to reduce errors, it is still unclear whether these tuned parameters based on particular datasets can improve model skills across different locations, model resolutions, or atmospheric events. As described above, it is known that model results are sensitive to the empirical values in convection. To summarize the numerous sensitivity studies here, some have reported that the location and intensity of precipitation are extremely sensitive to cumulus parameterization (e.g., Bechtold et al., 2008; Ma and Tan, 2009; Chikira and Sugiyama, 2010). For instance, Wang et al. (2007) improved the simulated diurnal cycle over land and ocean by increasing the entrainment–detrainment rates for deep and shallow convection used in the Tiedtke scheme, which tends to simulate convective precipitation too early in the day and with an unrealistic amplitude over land. Thus, the choice of a convective scheme impacts the diurnal cycle (e.g., Bechtold et al., 2004; Wang et al., 2007), as well as the simulation of monsoon precipitation in climate models (e.g., Mukhopadhyay et al., 2010), the MJO (e.g., Lin et al., 2006), the ENSO (e.g., Wu et al., 2007; Neale et al., 2008), the ITCZ configuration (e.g., Liu et al., 2019), and cloud cover and precipitation over urban areas (e.g., Karlický et al., 2020), among others. This topic has profound practical effects: it has been shown that choices in the convective parameterization affect the prediction of the track, intensity, and associated rainfall of tropical cyclones (e.g., Mohandas and Ashrit, 2014). However, the impacts of the empirical values in convection are extremely code-specific, and often errors in calibration of one parameter are hidden by errors in another. Examples of these include masking errors in vertical structure due to errors in cloud overlap (Neggers and Siebesma, 2013) or the “too few, too bright” problem (e.g., Nam et al., 2014). Therefore, results obtained in one GCM with a particular set of empirical values might differ from results obtained in a different GCM with the same set of empirical values.

Timely provision of the correct amount of precipitation at the right location is still a challenge for models. In the weather realm, Fig. 2 is an example of how different the precipitation field may look depending on the cumulus parameterization used. All a priori sensible methods locate the maximum and minima in different parts of Typhoon Megi and predict different areas and total accumulations. Figure 4 shows differences in the location and pressure of Typhoons Megi and Chaba with initial perturbations, and when seven different convection parameters are perturbed using SPP. Compared to the initial perturbations, changes in convection parameters show a bigger dispersion and yield a wider range of pressure values for each of the cyclones. In the climate model realm, validation exercises focusing on precipitation (Tapiador et al., 2012, 2017, 2018) have shown the importance and challenges of comparing model outputs with precipitation measurements in order to improve model performance. Indeed, the difficulties of quantitative precipitation estimation suggest precipitation as a privileged metric to gauge model performance (Tapiador et al., 2019a). The “ultimate test”, as it has been described, makes precipitation science an active field of research. As discussed in such papers, there is no complete agreement even in the reference data, with datasets differing even on aggregated values such as the global mean value of precipitation on Earth. Advances in satellite precipitation estimation (Kummerow et al., 1998; Joyce et al., 2004; Okamoto et al., 2005; Ushio and Kachi, 2010; Watanabe et al., 2010, 2011; Kucera et al., 2013; Hou et al., 2014; Huffman et al., 2015; Xie et al., 2017; Levizzani and Cattani, 2019; Skofronick-Jackson et al., 2019) are indispensable for further advancement, since direct estimates of precipitation (pluviometers, disdrometers) and ground radars are limited to land areas. In the near future, it is likely that satellites will continue to play a vital role in validating models and therefore in opening new directions in the way key physical processes are modeled. These advances need to be parallel with an explicit account of what is empirical in models in order to benefit fields, observations, and models. Algorithm developers in the satellite realm are perhaps more used to specifying their assumptions through the Algorithm Theoretical Basis Documents (ATBD), but a full comparison between the physics and empirical values behind both algorithms and parameterizations is much needed to advance the field. On that note, it is clear that better access to climate model code would contribute to addressing scientific gaps in climate models and improving their reliability (Añel et al., 2021). It would also be highly desirable that scientists not only specify the parameterizations they have used, but also the assumptions and empirical values they have actually selected within these. Tables 1–15 can be used to easily identify and pinpoint their choices. The benefit will be immense as some discrepancies could be readily attributed to known issues (i.e., heavy spurious rainfall over warm water in adjustment schemes) or identified as confounding variables. As in the case of the microphysics, making transparent the codes, the assumptions, and the empiricisms can only benefit the community and dispel any potential concerns.

As a final comment, it is important to note that the focus of this paper is not comparing the publicly available convection schemes or steering users towards one or another but to explore the physics behind the modules and to do that from an objective and independent point of view. Neither is the paper about criticizing the simplifications that are inherent to modeling the atmosphere or the limitations of current methods. On the contrary, the research arises from the conviction that models are the way forward to advance climate research. Being aware of the potential misuse of the results shown here to discredit models, it is important to discourage uninformed critics and futile attempts: neither this paper nor Tapiador et al. (2019b) cast any shadow of doubt on model outputs. On the contrary, they display and celebrate the delicate intricacies, nuances, precise measurements, and careful choices made by the community to craft complex tools to forecast, simulate, and predict precipitation.

Appendix A: List of acronyms
Acronym Meaning Acronym Meaning
ADHOC Assumed-distribution higher-order closure HWRF Hurricane Weather Research and Forecasting model
ALARO Aire Limitée Adaptation/Application de la Recherche à l'Opérationnel (ALARO). ICON Icosahedral nonhydrostatic model
ALE Available lifting energy IFS Integrated Forecasting System
ALP Available lifting power IN INCUS Ice Nuclei Investigation of Convective Updrafts Mission
AM4.0 Atmospheric Model version 4 IOP Intensive observation period
AOT Aerosol optical thickness ITCZ Intertropical Convergence Zone
ARM Atmospheric radiation measurement KF Kain–Fritsch scheme
ARW Advanced Research WRF KIM Koel isolatie maatschappij (The Netherlands Institute for Transport Policy Analysis)
AS Arakawa–Schubert scheme KWAJEX Kwajalein Experiment
ATBD Algorithm Theoretical Basis Documents LBN Level of neutral buoyancy
ATEX Atlantic Trade-Wind Experiment LCL Lifting condensation level
BCL Buoyant condensation level LFC Level of free convection
BMJ Betts–Miller–Janjić LFS Level of free sinking
BRAMS Brazilian developments on the Regional Atmospheric Modeling System LMDZ Laboratoire de Météorologie Dynamique Zoom
BOMEX Barbados Oceanographic and Meteorological EXperiment LWC Liquid water content
CA Cellular automaton MIROC Model for Interdisciplinary Research on Climate
CAM Community Atmosphere Model MJO Madden–Julian oscillation
CAPE Convective available potential energy MM5 Mesoscale Model version 5
CCM3 Community Climate Model version 3 MMF Multiscale model framework
CCN Cloud condensation nuclei MP Microphysics parameterization
CCSM Community Climate System Model NAM North American Mesoscale model
CDNC Cloud droplet number concentration NAVGEM Navy Global Environmental Model
CESM Community Earth System Model NCAR National Center for Atmospheric Research
CFSv2 Climate Forecast System version 2 NCEP National Centers for Environmental Prediction
CIN Convective inhibition NWP Numerical weather prediction
CISK Conditional instability of the second kind PBL Planetary boundary layer
CLUBB Cloud Layers Unified By Binomials PCAPE Integral over pressure of the buoyancy of an entraining ascending parcel with density scaling
COARE Coupled Ocean–Atmosphere Response Experiment PDF Probability density function
CP Cumulus parameterization PECAN Plains Elevated Convection at Night
CRCP Cloud-resolving convective parameterization PML Potential mixed layer
CRM Cloud-resolving model QE Quasi-equilibrium
CSRM Cloud-system-resolving model RACORO Routine AAF (ARM Aerial Facility) CLOWD (Clouds with Low Optical Water Depths) Optical Radiative Observations
CWF Cloud work function RAS Relaxed Arakawa–Schubert scheme
DBL Downdraft base layer RCM Regional climate model
dCAPE Dynamic convective available potential energy RH Relative humidity
Acronym Meaning Acronym Meaning
DDL Downdraft detrainment level RICO Rain In Cumulus over the Ocean field campaign
DualM Dual mass flux framework SAS Simplified Arakawa–Schubert scheme
DYNAMO Dynamics of the Madden–Julian Oscillation SCAM Single-column Community Atmosphere Model
ECHAM General circulation model developed by the Max Planck Institute for Meteorology SCM Single cloud model
ECMWF European Centre for Medium-Range Forecasts SGP97 Southern Great Plains 97
EDMF Eddy diffusivity mass flux SILHS Subgrid importance Latin hypercube sampler
EL Equilibrium level SNU Seoul National University
ENSO El Niño–Southern Oscillation SP Super-parameterization
EPS Ensemble prediction system
ESM Earth system model SPCZ South Pacific Convergence Zone
EUREC4A Elucidating the Role of Clouds–Circulation Coupling in Climate SST Sea surface temperature
GARP Global Atmospheric Research Program STOMP STOchastic framework for Modeling Population dynamics of convective clouds
GATE GARP Atlantic Tropical Experiment TC Tropical cyclone
GCM Global circulation/climate model TKE Turbulent kinetic energy
GEOS-5 Goddard Earth Observing System version 5 model TOGA Tropical Ocean–Global Atmosphere
GFDL Geophysical Fluid Dynamics Laboratory TWP-ICE Tropical Warm Pool – International Cloud Experiment
GFS Global Forecast System UIUC University of Illinois, Urban–Champaign
GISS GCM Goddard Institute for Space Studies Global Climate Model UM Unified Model
GOAmazon Green Ocean Amazon field campaign UNICON Unified convection scheme
HadGEM3 GA2.0 Hadley Centre Global Environmental model Global Atmosphere version 2 USL Updraft source layer
HCF Heated condensation framework WRF Weather Research and Forecasting model
Data availability

GFS data have been used to run the WRF simulations in Figs. 1 and 2. These datasets are publicly available online at (National Centers for Environmental Prediction/National Weather Service/NOAA/US Department of Commerce, 2015). Also, the GPM IMERG final run daily precipitation dataset was used in the comparison in Fig. 1. This dataset is publicly available online at (Huffman et al., 2019).

Author contributions

Conceptualization was done by FJT and AVP. FJT acquired funding. Investigation was handled by FJT and AVP. Methodology was conducted by FJT and AVP. Supervision was provided by FJT. AVP contributed to writing the original draft. FJT and AVP contributed to writing, review, and editing.

Competing interests

The contact author has declared that neither they nor their co-author has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Funding from projects PID2019-108470RB-C21 (AEI/FEDER, UE) and CGL2016-80609-R is gratefully acknowledged. Anahí Villalba-Pradas acknowledges support from grant FPI BES-2017-079685 for conducting her PhD. We are grateful to two anonymous referees for their valuable comments. Special thanks are due to Peter Bechtold for kindly performing the sensitivity experiments depicted in Fig. 3 with the IFS model during the research stay of Anahí Villalba-Pradas at ECMWF in February 2022 and for making some observations and suggestions that certainly improved the revised version of the paper.

Financial support

This research has been supported by the Agencia Estatal de Investigación (grant no. PID2019-108470RB-C21), the Ministerio de Economía y Competitividad (grant no. CGL2016-80609-R), and the Ministerio de Ciencia, Innovación y Universidades (grant no. BES-2017-079685).

Review statement

This paper was edited by Chiel van Heerwaarden and reviewed by two anonymous referees.


Albrecht, B. A., Betts, A. K., Schubert, W. H., and Cox, S. K.: Model of the Thermodynamic Structure of the Trade-Wind Boundary Layer: Part I. Theoretical Formulation and Sensitivity Tests, J. Atmos. Sci., 36, 73–89,<0073:MOTTSO>2.0.CO;2, 1979. 

Alexander, G. D. and Cotton, W. R.: The Use of Cloud-Resolving Simulations of Mesoscale Convective Systems to Build a Mesoscale Parameterization Scheme, J. Atmos. Sci., 55, 2137–2161,<2137:TUOCRS>2.0.CO;2, 1998. 

Allan, R. P. and Soden, B. J.: Atmospheric Warming and the Amplification of Precipitation Extremes, Science, 321, 1481–1484,,2008. 

Anderson, J. L., Balaji, V., Broccoli, A. J., Cooke, W. F., Delworth, T. L., Dixon, K. W., Donner, L. J., Dunne, K. A., Freidenreich, S. M., Garner, S. T., and Gudgel, R. G.: The New GFDL Global Atmosphere and Land Model AM2–LM2: Evaluation with Prescribed SST Simulations, J. Climate, 17, 4641–4673,, 2004. 

Añel, J. A., García-Rodríguez, M., and Rodeiro, J.: Current status on the need for improved accessibility to climate models code, Geosci. Model Dev., 14, 923–934,, 2021. 

Angevine, W. M.: An Integrated Turbulence Scheme for Boundary Layers with Shallow Cumulus Applied to Pollutant Transport, J. Appl. Meteorol., 44, 1436–1452,, 2005. 

Angevine, W. M., Jiang, H., and Mauritsen, T.: Performance of an Eddy Diffusivity–Mass Flux Scheme for Shallow Cumulus Boundary Layers, Mon. Weather Rev., 138, 2895–2912,, 2010. 

Anthes, R. A.: A Cumulus Parameterization Scheme Utilizing a One-Dimensional Cloud Model, Mon. Weather Rev., 138, 2895–2912,<0270:ACPSUA>2.0.CO;2, 1977. 

Arakawa, A.: Parameterization of cumulus convection, Proceedings of WMO/IUGG Symposium, Numerical Weather Prediction, Japan Meteorological Agency, IV,8, 1–6, 1969. 

Arakawa, A.: The Cumulus Parameterization Problem: Past, Present, and Future, J. Climate, 17, 2493–2525,<2493:RATCPP>2.0.CO;2, 2004. 

Arakawa, A. and Schubert, W. H.: Interaction of a Cumulus Cloud Ensemble with the Large-Scale Environment, Part I., J. Atmos. Sci., 31, 674–701,<0674:IOACCE>2.0.CO;2, 1974. 

Arakawa, A. and Wu, C.-M.: A Unified Representation of Deep Moist Convection in Numerical Modeling of the Atmosphere. Part I, J. Atmos. Sci., 70, 1977–1992,, 2013. 

Arakawa, A., Jung, J.-H., and Wu, C.-M.: Toward unification of the multiscale modeling of the atmosphere, Atmos. Chem. Phys., 11, 3731–3742,, 2011. 

Asai, T. and Kasahara, A.: A Theoretical Study of the Compensating Downward Motions Associated with Cumulus Clouds, J. Atmos. Sci., 24, 487–496,<0487:ATSOTC>2.0.CO;2, 1967. 

Baba, Y.: Spectral cumulus parameterization based on cloud-resolving model, Clim. Dynam., 52, 309–334,, 2019. 

Baik, J.-J., DeMaria, M., and Raman, S.: Tropical Cyclone Simulations with the Betts Convective Adjustment Scheme. Part II: Sensitivity Experiments, Mon. Weather Rev., 118, 529–541,<0529:TCSWTB>2.0.CO;2, 1990. 

Bak, P., Tang, C., and Wiesenfeld, K.: Self-organized criticality: An explanation of the 1/f noise, Phys. Rev. Lett., 59, 381,, 1987. 

Baldwin, M. E., Kain, J. S., and Kay, M. P.: Properties of the Convection Scheme in NCEP's Eta Model that Affect Forecast Sounding Interpretation, Weather Forecast., 17, 1063–1079,<1063:POTCSI>2.0.CO;2, 2002. 

Barros, D. F., Albernaz, A. L. M., Barros, D. F., and Albernaz, A. L. M.: Possible impacts of climate change on wetlands and its biota in the Brazilian Amazon, Braz. J. Biol., 74, 810–820,, 2014. 

Bechtold, P. (Ed.): Atmospheric moist convection, Meteorological Training Course Lecture Series, ECMWF, (last access: 10 September 2021), 2019. 

Bechtold, P., Pinty, J. P., and Fravalo, C.: A Model of Marine Boundary-Layer Cloudiness for Mesoscale Applications, J. Atmos. Sci., 49, 1723–1744,<1723:AMOMBL>2.0.CO;2, 1992. 

Bechtold, P., Cuijpers, J. W. M., Mascart, P., and Trouilhet, P.: Modeling of Trade Wind Cumuli with a Low-Order Turbulence Model: Toward a Unified Description of Cu and Se Clouds in Meteorological Models, J. Atmos. Sci., 52, 455–463,<0455:MOTWCW>2.0.CO;2, 1995. 

Bechtold, P., Bazile, E., Guichard, F., Mascart, P., and Richard, E.: A mass-flux convection scheme for regional and global models, Q. J. Roy. Meteor. Soc., 127, 869–886,, 2001. 

Bechtold, P., Chaboureau, J.-P., Beljaars, A., Betts, A. K., Köhler, M., Miller, M., and Redelsperger, J.-L.: The simulation of the diurnal cycle of convective precipitation over land in a global model, Q. J. Roy. Meteor. Soc., 130, 3119–3137,, 2004. 

Bechtold, P., Köhler, M., Jung, T., Doblas-Reyes, F., Leutbecher, M., Rodwell, M. J., Vitart, F., and Balsamo, G.: Advances in simulating atmospheric variability with the ECMWF model: From synoptic to decadal time-scales, Q. J. Roy. Meteor. Soc., 134, 1337–1351,, 2008. 

Bechtold, P., Semane, N., Lopez, P., Chaboureau, J.-P., Beljaars, A., and Bormann, N.: Representing Equilibrium and Nonequilibrium Convection in Large-Scale Models, J. Atmos. Sci., 71, 734–753,, 2014. 

Becker, T. and Hohenegger, C.: Estimating Bulk Entrainment for Deep Convection – from Idealized to Realistic Simulations, American Geophysical Union, Fall Meeting 2018, Washington, D.C., abstract #A21K-2864, 21, 2018. 

Becker, T., Bechtold, P., and Sandu, I.: Characteristics of convective precipitation over tropical Africa in storm‐resolving global simulations, Q. J. Roy. Meteor. Soc., 147, 4388–4407,, 2021. 

Bengtsson, L., Körnich, H., Källén, E., and Svensson, G.: Large-Scale Dynamical Response to Subgrid-Scale Organization Provided by Cellular Automata, J. Atmos. Sci., 68, 3132–3144,, 2011. 

Bengtsson, L., Steinheimer, M., Bechtold, P., and Geleyn, J.-F.: A stochastic parametrization for deep convection using cellular automata, Q. J. Roy. Meteor. Soc., 139, 1533–1543,, 2013. 

Bengtsson, L., Bao, J.-W., Pegion, P., Penland, C., Michelson, S., and Whitaker, J.: A Model Framework for Stochastic Representation of Uncertainties Associated with Physical Processes in NOAA's Next Generation Global Prediction System (NGGPS), Mon. Weather Rev., 147, 893–911,, 2019. 

Bengtsson, L., Dias, J., Tulich, S., Gehne, M., and Bao, J.-W.: A Stochastic Parameterization of Organized Tropical Convection Using Cellular Automata for Global Forecasts in NOAA's Unified Forecast System, J. Adv. Model Earth Sy., 13, e2020MS002260,, 2021. 

Berg, L. K., Gustafson, W. I., Kassianov, E. I., and Deng, L.: Evaluation of a Modified Scheme for Shallow Convection: Implementation of CuP and Case Studies, Mon. Weather Rev., 141, 134–147,, 2013. 

Betts, A. K.: Parametric Interpretation of Trade-Wind Cumulus Budget Studies, J. Atmos. Sci., 32, 1934–1945,<1934:PIOTWC>2.0.CO;2, 1975. 

Betts, A. K.: Saturation Point Analysis of Moist Convective Overturning, J. Atmos. Sci., 39, 1484–1505,<1484:SPAOMC>2.0.CO;2, 1982. 

Betts, A. K.: Mixing Line Analysis of Clouds and Cloudy Boundary Layers, J. Atmos. Sci., 42, 2751–2763,<2751:MLAOCA>2.0.CO;2, 1985. 

Betts, A. K.: A new convective adjustment scheme. Part I: Observational and theoretical basis, Q. J. Roy. Meteor. Soc., 112, 677–691,, 1986. 

Betts, A. K. and Albrecht, B. A.: Conserved Variable Analysis of the Convective Boundary Layer Thermodynamic Structure over the Tropical Oceans, J. Atmos. Sci., 44, 83–99,<0083:CVAOTC>2.0.CO;2, 1987. 

Betts, A. K. and Jakob, C.: Evaluation of the diurnal cycle of precipitation, surface thermodynamics, and surface fluxes in the ECMWF model using LBA data, J. Geophys. Res.-Atmos., 107, LBA 12-1–LBA 12-8,, 2002. 

Betts, A. K. and Miller, M. J.: A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, ATEX and arctic air-mass data sets, Q. J. Roy. Meteor. Soc., 112, 693–709,, 1986. 

Bhatla, R., Ghosh, S., Mandal, B., Mall, R. K., and Sharma, K.: Simulation of Indian summer monsoon onset with different parameterization convection schemes of RegCM-4.3, Atmos. Res., 176–177, 10–18,, 2016. 

Bhattacharya, R., Bordoni, S., Suselj, K., and Teixeira, J.: Parameterization Interactions in Global Aquaplanet Simulations, J. Adv. Model Earth Sy., 10, 403–420,, 2018. 

Blyth, A. M.: Entrainment in Cumulus Clouds, J. Appl. Meteorol. Clim., 32, 626–641,<0626:EICC>2.0.CO;2, 1993. 

Blyth, A. M., Cooper, W. A., and Jensen, J. B.: A Study of the Source of Entrained Air in Montana Cumuli, J. Atmos. Sci., 45, 3944–3964,<3944:ASOTSO>2.0.CO;2, 1988. 

Boatman, J. F. and Auer, A. H.: The Role of Cloud Top Entrainment in Cumulus Clouds, J. Atmos. Sci., 40, 1517–1534,<1517:TROCTE>2.0.CO;2, 1983. 

Bogenschutz, P. A. and Krueger, S. K.: A simplified PDF parameterization of subgrid-scale clouds and turbulence for cloud-resolving models, J. Adv. Model Earth Sy., 5, 195–211,, 2013. 

Bogenschutz, P. A., Krueger, S. K., and Khairoutdinov, M.: Assumed Probability Density Functions for Shallow and Deep Convection, J. Adv. Model Earth Sy., 2, 10,, 2010. 

Böing, S. J.: An object-based model for convective cold pool dynamics, Mathematics of Climate and Weather Forecasting, 2, 43–60,, 2016. 

Böing, S. J., Jonker, H. J. J., Siebesma, A. P., and Grabowski, W. W.: Influence of the Subcloud Layer on the Development of a Deep Convective Ensemble, J. Atmos. Sci., 69, 2682–2698,, 2012. 

Böing, S. J., Jonker, H. J. J., Nawara, W. A., and Siebesma, A. P.: On the Deceiving Aspects of Mixing Diagrams of Deep Cumulus Convection, J. Atmos. Sci., 71, 56–68,, 2014. 

Bombardi, R. J., Schneider, E. K., Marx, L., Halder, S., Singh, B., Tawfik, A. B., Dirmeyer, P. A., and Kinter, J. L.: Improvements in the representation of the Indian summer monsoon in the NCEP climate forecast system version 2, Clim. Dynam., 45, 2485–2498,, 2015. 

Bombardi, R. J., Tawfik, A. B., Manganello, J. V., Marx, L., Shin, C.-S., Halder, S., Schneider, E. K., Dirmeyer, P. A., and Kinter, J. L.: The heated condensation framework as a convective trigger in the NCEP Climate Forecast System version 2, J. Adv. Model Earth Sy., 8, 1310–1329,, 2016. 

Bony, S. and Dufresne, J.-L.: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models, Geophys. Res. Lett., 32, L20806,, 2005. 

Bony, S. and Emanuel, K. A.: A Parameterization of the Cloudiness Associated with Cumulus Convection; Evaluation Using TOGA COARE Data, J. Atmos. Sci., 58, 3158–3183,<3158:APOTCA>2.0.CO;2, 2001. 

Bony, S., Stevens, B., Frierson, D. M. W., Jakob, C., Kageyama, M., Pincus, R., Shepherd, T. G., Sherwood, S. C., Siebesma, A. P., Sobel, A. H., Watanabe, M., and Webb, M. J.: Clouds, circulation and climate sensitivity, Nat. Geosci, 8, 261–268,, 2015. 

Bony, S., Stevens, B., Ament, F., Bigorre, S., Chazette, P., Crewell, S., Delanoë, J., Emanuel, K., Farrell, D., Flamant, C., Gross, S., Hirsch, L., Karstensen, J., Mayer, B., Nuijens, L., Ruppert, J. H., Sandu, I., Siebesma, P., Speich, S., Szczap, F., Totems, J., Vogel, R., Wendisch, M., and Wirth, M.: EUREC4A: A Field Campaign to Elucidate the Couplings Between Clouds, Convection and Circulation, Surv. Geophys., 38, 1529–1568,, 2017. 

Bougeault, P.: Cloud-Ensemble Relations Based on the Gamma Probability Distribution for the Higher-Order Models of the Planetary Boundary Layer, J. Atmos. Sci., 39, 2691–2700,<2691:CERBOT>2.0.CO;2, 1982. 

Bougeault, P.: A Simple Parameterization of the Large-Scale Effects of Cumulus Convection, Mon. Weather Rev., 113, 2108–2121,<2108:ASPOTL>2.0.CO;2, 1985. 

Brast, M., Schemann, V., and Neggers, R. A. J.: Investigating the Scale Adaptivity of a Size-Filtered Mass Flux Parameterization in the Gray Zone of Shallow Cumulus Convection, J. Atmos. Sci., 75, 1195–1214,, 2018. 

Bretherton, C. S., McCaa, J. R., and Grenier, H.: A New Parameterization for Shallow Cumulus Convection and Its Application to Marine Subtropical Cloud-Topped Boundary Layers. Part I: Description and 1D Results, Mon. Weather Rev., 132, 864–882,<0864:ANPFSC>2.0.CO;2, 2004. 

Bright, D. R. and Mullen, S. L.: Short-Range Ensemble Forecasts of Precipitation during the Southwest Monsoon, Weather Forecast., 17, 1080–1100,<1080:SREFOP>2.0.CO;2, 2002. 

Brisson, E., Van Weverberg, K., Demuzere, M., Devis, A., Saeed, S., Stengel, M., and van Lipzig, N. P. M.: How well can a convection-permitting climate model reproduce decadal statistics of precipitation, temperature and cloud characteristics?, Clim. Dynam., 47, 3043–3061,, 2016. 

Bryan, G. H., Wyngaard, J. C., and Fritsch, J. M.: Resolution Requirements for the Simulation of Deep Moist Convection, Mon. Weather Rev., 131, 2394–2416,<2394:RRFTSO>2.0.CO;2, 2003. 

Buizza, R., Milleer, M., and Palmer, T. N.: Stochastic representation of model uncertainties in the ECMWF ensemble prediction system, Q. J. Roy. Meteor. Soc., 125, 2887–2908,, 1999. 

Burnet, F. and Brenguier, J.-L.: Observational Study of the Entrainment-Mixing Process in Warm Convective Clouds, J. Atmos. Sci., 64, 1995–2011,, 2007. 

Cahalan, R. F., Ridgway, W., Wiscombe, W. J., Bell, T. L., and Snider, J. B.: The Albedo of Fractal Stratocumulus Clouds, J. Atmos. Sci., 51, 2434–2455,<2434:TAOFSC>2.0.CO;2, 1994. 

Chaboureau, J.-P. and Bechtold, P.: A Simple Cloud Parameterization Derived from Cloud Resolving Model Data: Diagnostic and Prognostic Applications, J. Atmos. Sci., 59, 2362–2372,<2362:ASCPDF>2.0.CO;2, 2002. 

Chaboureau, J.-P. and Bechtold, P.: Statistical representation of clouds in a regional model and the impact on the diurnal cycle of convection during Tropical Convection, Cirrus and Nitrogen Oxides (TROCCINOX), J. Geophys. Res.-Atmos., 110, D17103,, 2005. 

Charney, J. G. and Eliassen, A.: On the Growth of the Hurricane Depression, J. Atmos. Sci., 21, 68–75,<0068:OTGOTH>2.0.CO;2, 1964. 

Chatfield, R. B. and Brost, R. A.: A two-stream model of the vertical transport of trace species in the convective boundary layer, J. Geophys. Res.-Atmos., 92, 13263–13276,, 1987. 

Cheinet, S.: A Multiple Mass-Flux Parameterization for the Surface-Generated Convection. Part I: Dry Plumes, J. Atmos. Sci., 60, 2313–2327,<2313:AMMPFT>2.0.CO;2, 2003. 

Cheinet, S.: A Multiple Mass Flux Parameterization for the Surface-Generated Convection. Part II: Cloudy Cores, J. Atmos. Sci., 61, 1093–1113,<1093:AMMFPF>2.0.CO;2, 2004. 

Cheng, A. and Xu, K.-M.: Simulation of shallow cumuli and their transition to deep convective clouds by cloud-resolving models with different third-order turbulence closures, Q. J. Roy. Meteor. Soc., 132, 359–382,, 2006. 

Chikira, M.: A Cumulus Parameterization with State-Dependent Entrainment Rate. Part II: Impact on Climatology in a General Circulation Model, J. Atmos. Sci., 67, 2194–2211,, 2010. 

Chikira, M. and Sugiyama, M.: A Cumulus Parameterization with State-Dependent Entrainment Rate. Part I: Description and Sensitivity to Temperature and Humidity Profiles, J. Atmos. Sci., 67, 2171–2193,, 2010. 

Choat, B., Jansen, S., Brodribb, T. J., Cochard, H., Delzon, S., Bhaskar, R., Bucci, S. J., Feild, T. S., Gleason, S. M., Hacke, U. G., Jacobsen, A. L., Lens, F., Maherali, H., Martínez-Vilalta, J., Mayr, S., Mencuccini, M., Mitchell, P. J., Nardini, A., Pittermann, J., Pratt, R. B., Sperry, J. S., Westoby, M., Wright, I. J., and Zanne, A. E.: Global convergence in the vulnerability of forests to drought, Nature, 491, 752–755,, 2012. 

Chopard, B.: Cellular Automata Modeling of Physical Systems, in: Encyclopedia of Complexity and Systems Science, edited by: Meyers, R. A., Encyclopedia of Complexity and Systems Science Springer, New York, NY, 865–892,, 2009. 

Cohen, Y., Lopez-Gomez, I., Jaruga, A., He, J., Kaul, C. M., and Schneider, T.: Unified Entrainment and Detrainment Closures for Extended Eddy-Diffusivity Mass-Flux Schemes, J. Adv. Model Earth Sy., 12, e2020MS002162,, 2020. 

Colin, M.: Convective memory, and the role of cold pools, Meteorology, Sorbonne Université, HAL Id: tel-02864797, 2018. 

Colin, M., Sherwood, S., Geoffroy, O., Bony, S., and Fuchs, D.: Identifying the Sources of Convective Memory in Cloud-Resolving Simulations, J. Atmos. Sci., 76, 947–962,, 2019. 

Collier, J. C. and Bowman, K. P.: Diurnal cycle of tropical precipitation in a general circulation model, J. Geophys. Res.-Atmos., 109, D17105,, 2004. 

Couvreux, F., Hourdin, F., Williamson, D., Roehrig, R., Volodina, V., Villefranque, N., Rio, C., Audouin, O., Salter, J., Bazile, E., Brient, F., Favot, F., Honnert, R., Lefebvre, M.-P., Madeleine, J.-B., Rodier, Q., and Xu, W.: Process-Based Climate Model Development Harnessing Machine Learning: I. A Calibration Tool for Parameterization Improvement, J. Adv. Model Earth Sy., 13, e2020MS002217,, 2021. 

Cotton, W. and Anthes, R.: Storm and Cloud Dynamics, 1st edn., Academic Press, 1992. 

Craig, G. C. and Cohen, B. G.: Fluctuations in an Equilibrium Convective Ensemble. Part I: Theoretical Formulation, J. Atmos. Sci., 63, 1996–2004,, 2006. 

Dai, A.: Precipitation Characteristics in Eighteen Coupled Climate Models, J. Climate, 19, 4605–4630,, 2006. 

Dai, A. and Trenberth, K. E.: The Diurnal Cycle and Its Depiction in the Community Climate System Model, J. Climate, 17, 930–951,<0930:TDCAID>2.0.CO;2, 2004. 

D'Andrea, F., Gentine, P., Betts, A. K., and Lintner, B. R.: Triggering Deep Convection with a Probabilistic Plume Model, J. Atmos. Sci., 71, 3881–3901,, 2014. 

Davies, L., Plant, R. S., and Derbyshire, S. H.: A simple model of convection with memory, J. Geophys. Res.-Atmos., 114, D17202,, 2009. 

Davies, L., Jakob, C., Cheung, K., Genio, A. D., Hill, A., Hume, T., Keane, R. J., Komori, T., Larson, V. E., Lin, Y., Liu, X., Nielsen, B. J., Petch, J., Plant, R. S., Singh, M. S., Shi, X., Song, X., Wang, W., Whitall, M. A., Wolf, A., Xie, S., and Zhang, G.: A single-column model ensemble approach applied to the TWP-ICE experiment, J. Geophys. Res.-Atmos., 118, 6544–6563,, 2013a. 

Davies, L., Plant, R. S., and Derbyshire, S. H.: Departures from convective equilibrium with a rapidly varying surface forcing, Q. J. Roy. Meteor. Soc., 139, 1731–1746,, 2013b. 

Dawe, J. T. and Austin, P. H.: Direct entrainment and detrainment rate distributions of individual shallow cumulus clouds in an LES, Atmos. Chem. Phys., 13, 7795–7811,, 2013. 

Deardorff, J. W.: The Counter-Gradient Heat Flux in the Lower Atmosphere and in the Laboratory, J. Atmos. Sci., 23, 503–506,<0503:TCGHFI>2.0.CO;2, 1966. 

Deardorff, J. W., Willis, G. E., and Lilly, D. K.: Laboratory investigation of non-steady penetrative convection, J. Fliud Mech., 35, 7–31,, 1969. 

Deguines, N., Brashares, J. S., and Prugh, L. R.: Precipitation alters interactions in a grassland ecological community, J. Anim. Ecol., 86, 262–272,, 2017. 

Del Genio, A. D. and Wu, J.: The Role of Entrainment in the Diurnal Cycle of Continental Convection, J. Climate, 23, 2722–2738,, 2010. 

Del Genio, A. D., Chen, Y., Kim, D., and Yao, M.-S.: The MJO Transition from Shallow to Deep Convection in CloudSat/CALIPSO Data and GISS GCM Simulations, J. Climate, 25, 3755–3770,, 2012. 

Del Genio, A. D., Wu, J., Wolf, A. B., Chen, Y., Yao, M.-S., and Kim, D.: Constraints on Cumulus Parameterization from Simulations of Observed MJO Events, J. Climate, 28, 6419–6442,, 2015. 

DeMott, C. A., Randall, D. A., and Khairoutdinov, M.: Convective Precipitation Variability as a Tool for General Circulation Model Analysis, J. Climate, 20, 91–112,, 2007. 

Deng, A., Seaman, N. L., and Kain, J. S.: A Shallow-Convection Parameterization for Mesoscale Models. Part I: Submodel Description and Preliminary Applications, J. Atmos. Sci., 60, 34–56,<0034:ASCPFM>2.0.CO;2, 2003. 

Deng, Q., Khouider, B., and Majda, A. J.: The MJO in a Coarse-Resolution GCM with a Stochastic Multicloud Parameterization, J. Atmos. Sci., 72, 55–74,, 2015. 

Derbyshire, S. H., Maidens, A. V., Milton, S. F., Stratton, R. A., and Willett, M. R.: Adaptive detrainment in a convective parametrization, Q. J. Roy. Meteor. Soc., 137, 1856–1871,, 2011. 

de Roode, S. R., Siebesma, A. P., Jonker, H. J. J., and de Voogd, Y.: Parameterization of the Vertical Velocity Equation for Shallow Cumulus Clouds, Mon. Weather Rev., 140, 2424–2436,, 2012. 

De Rooy, W. C. and Siebesma, A. P.: A Simple Parameterization for Detrainment in Shallow Cumulus, Mon. Weather Rev., 136, 560–576,, 2008. 

De Rooy, W. C. and Siebesma, A. P.: Analytical expressions for entrainment and detrainment in cumulus convection, Q. J. Roy. Meteor. Soc., 136, 1216–1227,, 2010. 

De Rooy, W. C., Bechtold, P., Fröhlich, K., Hohenegger, C., Jonker, H., Mironov, D., Siebesma, A. P., Teixeira, J., and Yano, J.-I.: Entrainment and detrainment in cumulus convection: an overview, Q. J. Roy. Meteor. Soc., 139, 1–19,, 2013. 

Donner, L. J.: A Cumulus Parameterization Including Mass Fluxes, Vertical Momentum Dynamics, and Mesoscale Effects, J. Climate, 50, 889–906,<0889:ACPIMF>2.0.CO;2, 1993. 

Donner, L. J. and Phillips, V. T.: Boundary layer control on convective available potential energy: Implications for cumulus parameterization, J. Geophys. Res.-Atmos., 108, 4701,, 2003. 

Donner, L. J., Seman, C. J., Hemler, R. S., and Fan, S.: A Cumulus Parameterization Including Mass Fluxes, Convective Vertical Velocities, and Mesoscale Effects: Thermodynamic and Hydrological Aspects in a General Circulation Model, J. Atmos. Sci., 14, 3444–3463,<3444:ACPIMF>2.0.CO;2, 2001. 

Donner, L. J., Wyman, B. L., Hemler, R. S., Horowitz, L. W., Ming, Y., Zhao, M., Golaz, J.-C., Ginoux, P., Lin, S.-J., Schwarzkopf, M. D., Austin, J., Alaka, G., Cooke, W. F., Delworth, T. L., Freidenreich, S. M., Gordon, C. T., Griffies, S. M., Held, I. M., Hurlin, W. J., Klein, S. A., Knutson, T. R., Langenhorst, A. R., Lee, H.-C., Lin, Y., Magi, B. I., Malyshev, S. L., Milly, P. C. D., Naik, V., Nath, M. J., Pincus, R., Ploshay, J. J., Ramaswamy, V., Seman, C. J., Shevliakova, E., Sirutis, J. J., Stern, W. F., Stouffer, R. J., Wilson, R. J., Winton, M., Wittenberg, A. T., and Zeng, F.: The Dynamical Core, Physical Parameterizations, and Basic Simulation Characteristics of the Atmospheric Component AM3 of the GFDL Global Coupled Model CM3, J. Climate, 24, 3484–3519,, 2011. 

Donner, L. J., O'Brien, T. A., Rieger, D., Vogel, B., and Cooke, W. F.: Are atmospheric updrafts a key to unlocking climate forcing and sensitivity?, Atmos. Chem. Phys., 16, 12983–12992,, 2016. 

Dore, M. H. I.: Climate change and changes in global precipitation patterns: What do we know?, Environ. Int., 31, 1167–1181,, 2005. 

Dorrestijn, J., Crommelin, D. T., Biello, J. A., and Böing, S. J.: A data-driven multi-cloud model for stochastic parametrization of deep convection, Philos. T. Roy. Soc. A., 371, 20120374,, 2013a. 

Dorrestijn, J., Crommelin, D. T., Siebesma, A. Pier., and Jonker, H. J. J.: Stochastic parameterization of shallow cumulus convection estimated from high-resolution model data, Theor. Comp. Fluid Dyn., 27, 133–148,, 2013b. 

Dorrestijn, J., Crommelin, D. T., Siebesma, A. P., Jonker, H. J. J., and Jakob, C.: Stochastic Parameterization of Convective Area Fractions with a Multicloud Model Inferred from Observational Data, J. Atmos. Sci., 72, 854–869,, 2015. 

Drueke, S., Kirshbaum, D. J., and Kollias, P.: Evaluation of Shallow-Cumulus Entrainment Rate Retrievals Using Large-Eddy Simulation, J. Geophys. Res.-Atmos., 124, 9624–9643,, 2019. 

Easterling, D. R., Meehl, G. A., Parmesan, C., Changnon, S. A., Karl, T. R., and Mearns, L. O.: Climate Extremes: Observations, Modeling, and Impacts, Science, 289, 2068–2074,, 2000. 

Emanuel, K.: Atmospheric convection, Oxford University Press, 592 pp., ISBN 0-19-506630-8, 1994. 

Emanuel, K. and Raymond, D. J. (Eds.): The Representation of Cumulus Convection in Numerical Models of the Atmosphere, American Meteorological Society, 246 pp.,, 1993. 

Emanuel, K. A.: The Finite-Amplitude Nature of Tropical Cyclogenesis, J. Atmos. Sci., 46, 3431–3456,<3431:TFANOT>2.0.CO;2, 1989. 

Emanuel, K. A.: A Scheme for Representing Cumulus Convection in Large-Scale Models, J. Atmos. Sci., 48, 2313–2329,<2313:ASFRCC>2.0.CO;2, 1991. 

Emanuel, K. A.: The Behavior of a Simple Hurricane Model Using a Convective Scheme Based on Subcloud-Layer Entropy Equilibrium, J. Atmos. Sci., 52, 3960–3968,<3960:TBOASH>2.0.CO;2, 1995. 

Emanuel, K. A. and Živković-Rothman, M.: Development and Evaluation of a Convection Scheme for Use in Climate Models, J. Atmos. Sci., 56, 1766–1782,<1766:DAEOAC>2.0.CO;2, 1999. 

Evans, J. P. and Westra, S.: Investigating the Mechanisms of Diurnal Rainfall Variability Using a Regional Climate Model, J. Climate, 25, 7232–7247,, 2012. 

Evans, J. P., Ekström, M., and Ji, F.: Evaluating the performance of a WRF physics ensemble over South-East Australia, Clim. Dynam., 39, 1241–1258,, 2012. 

Feingold, G.: Modeling of the first indirect effect: Analysis of measurement requirements, Geophys. Res. Lett., 30, 19,, 2003. 

Feingold, G. and Koren, I.: A model of coupled oscillators applied to the aerosol–cloud–precipitation system, Nonlin. Processes Geophys., 20, 1011–1021,, 2013. 

Fiori, E., Comellas, A., Molini, L., Rebora, N., Siccardi, F., Gochis, D. J., Tanelli, S., and Parodi, A.: Analysis and hindcast simulations of an extreme rainfall event in the Mediterranean area: The Genoa 2011 case, Atmos. Res., 138, 13–29,, 2014. 

Fletcher, J. K. and Bretherton, C. S.: Evaluating Boundary Layer–Based Mass Flux Closures Using Cloud-Resolving Model Simulations of Deep Convection, J. Atmos. Sci., 67, 2212–2225,, 2010. 

Folkins, I., Mitovski, T., and Pierce, J. R.: A simple way to improve the diurnal cycle in convective rainfall over land in climate models, J. Geophys. Res.-Atmos., 119, 2113–2130,, 2014. 

Fonseca, R. M., Zhang, T., and Yong, K.-T.: Improved simulation of precipitation in the tropics using a modified BMJ scheme in the WRF model, Geosci. Model Dev., 8, 2915–2928,, 2015. 

Freitas, S. R., Panetta, J., Longo, K. M., Rodrigues, L. F., Moreira, D. S., Rosário, N. E., Silva Dias, P. L., Silva Dias, M. A. F., Souza, E. P., Freitas, E. D., Longo, M., Frassoni, A., Fazenda, A. L., Santos e Silva, C. M., Pavani, C. A. B., Eiras, D., França, D. A., Massaru, D., Silva, F. B., Santos, F. C., Pereira, G., Camponogara, G., Ferrada, G. A., Campos Velho, H. F., Menezes, I., Freire, J. L., Alonso, M. F., Gácita, M. S., Zarzur, M., Fonseca, R. M., Lima, R. S., Siqueira, R. A., Braz, R., Tomita, S., Oliveira, V., and Martins, L. D.: The Brazilian developments on the Regional Atmospheric Modeling System (BRAMS 5.2): an integrated environmental model tuned for tropical areas, Geosci. Model Dev., 10, 189–222,, 2017. 

Freitas, S. R., Grell, G. A., Molod, A., Thompson, M. A., Putman, W. M., Silva, C. M. S. e, and Souza, E. P.: Assessing the Grell-Freitas Convection Parameterization in the NASA GEOS Modeling System, J. Adv. Model Earth Sy., 10, 1266–1289,, 2018. 

Freitas, S. R., Grell, G. A., and Li, H.: The Grell–Freitas (GF) convection parameterization: recent developments, extensions, and applications, Geosci. Model Dev., 14, 5393–5411,, 2021. 

Frenkel, Y., Majda, A. J., and Khouider, B.: Using the Stochastic Multicloud Model to Improve Tropical Convective Parameterization: A Paradigm Example, J. Atmos. Sci., 69, 1080–1105,, 2012. 

Fritsch, J. M. and Chappell, C. F.: Numerical Prediction of Convectively Driven Mesoscale Pressure Systems. Part I: Convective Parameterization, J. Atmos. Sci., 37, 1722–1733,<1722:NPOCDM>2.0.CO;2, 1980. 

Gallus, W. and Segal, M.: Impact of improved initialization of mesoscale features on convective system rainfall in 10-km Eta simulations, Weather Forecast., 16, 680–696,<0680:IOIIOM>2.0.CO;2, 2001. 

Gao, S., Lu, C., Liu, Y., Mei, F., Wang, J., Zhu, L., and Yan, S.: Contrasting Scale Dependence of Entrainment-Mixing Mechanisms in Stratocumulus Clouds, Geophys. Res. Lett., 47, e2020GL086970,, 2020. 

Gao, X.-J., Shi, Y., and Giorgi, F.: Comparison of convective parameterizations in RegCM4 experiments over China with CLM as the land surface model, Atmos. Ocean. Sc. Lett., 9, 246–254,, 2016. 

Gao, Y., Leung, L. R., Zhao, C., and Hagos, S.: Sensitivity of U.S. summer precipitation to model resolution and convective parameterizations across gray zone resolutions, J. Geophys. Res.-Atmos., 122, 2714–2733,, 2017. 

García-Morales, M. B. and Dubus, L.: Forecasting precipitation for hydroelectric power management: how to exploit GCM's seasonal ensemble forecasts, Int. J. Climatol., 27, 1691–1705,, 2007. 

García-Ortega, E., Lorenzana, J., Merino, A., Fernández-González, S., López, L., and Sánchez, J. L.: Performance of multi-physics ensembles in convective precipitation events over northeastern Spain, Atmos. Res., 190, 55–67,, 2017. 

Gebhardt, C., Theis, S. E., Paulat, M., and Ben Bouallègue, Z.: Uncertainties in COSMO-DE precipitation forecasts introduced by model perturbations and variation of lateral boundaries, Atmos. Res., 100, 168–177,, 2011. 

Geerts, B., Parsons, D., Ziegler, C. L., Weckwerth, T. M., Biggerstaff, M. I., Clark, R. D., Coniglio, M. C., Demoz, B. B., Ferrare, R. A., Gallus, W. A., Haghi, K., Hanesiak, J. M., Klein, P. M., Knupp, K. R., Kosiba, K., McFarquhar, G. M., Moore, J. A., Nehrir, A. R., Parker, M. D., Pinto, J. O., Rauber, R. M., Schumacher, R. S., Turner, D. D., Wang, Q., Wang, X., Wang, Z., and Wurman, J.: The 2015 Plains Elevated Convection at Night Field Project, B. Am. Meteorol. Soc., 98, 767–786,, 2017. 

Geleyn, J.-F.: On a Simple, Parameter-Free Partition between Moistening and Precipitation in the Kuo Scheme, Mon. Weather Rev., 113, 405–407,<0405:OASPFP>2.0.CO;2, 1985. 

Genio, A. D. D., Kovari, W., Yao, M.-S., and Jonas, J.: Cumulus Microphysics and Climate Sensitivity, J. Climate, 18, 2376–2387,, 2005. 

Gentine, P., Betts, A. K., Lintner, B. R., Findell, K. L., van Heerwaarden, C. C., Tzella, A., and D'Andrea, F.: A Probabilistic Bulk Model of Coupled Mixed Layer and Convection. Part I: Clear-Sky Case, J. Atmos. Sci., 70, 1543–1556,, 2013a. 

Gentine, P., Betts, A. K., Lintner, B. R., Findell, K. L., van Heerwaarden, C. C., and D'Andrea, F.: A Probabilistic Bulk Model of Coupled Mixed Layer and Convection. Part II: Shallow Convection Case, J. Atmos., Sci., 70, 1557–1576,, 2013b. 

Gentine, P., Pritchard, M., Rasp, S., Reinaudi, G., and Yacalis, G.: Could Machine Learning Break the Convection Parameterization Deadlock?, Geophys. Res. Lett., 45, 5742–5751,, 2018. 

Gerard, L.: An integrated package for subgrid convection, clouds and precipitation compatible with meso-gamma scales, Q. J. Roy. Meteor. Soc., 133, 711–730,, 2007. 

Gerard, L.: Bulk mass-flux perturbation formulation for a unified approach of deep convection at high resolution, Mon. Weather Rev., 143, 4038–4063,, 2015. 

Gerard, L. and Geleyn, J.-F.: Evolution of a subgrid deep convection parametrization in a limited-area model with increasing resolution, Q. J. Roy. Meteor. Soc., 131, 2293–2312,, 2005. 

Gerard, L., Brown, A. R., Cederwall, R. T., Chlond, A., Duynkerke, P. G., Golaz, J.-C., Khairoutdinov, M., Lewellen, D. C., Lock, A. P., MacVean, M. K., Moeng, C.-H., Neggers, R. a. J., Siebesma, A. P., and Stevens, B.: Large-eddy simulation of the diurnal cycle of shallow cumulus convection over land, Q. J. Roy. Meteor. Soc., 128, 1075–1093,, 2002. 

Gerard, L., Piriou, J.-M., Brožková, R., Geleyn, J.-F., and Banciu, D.: Cloud and Precipitation Parameterization in a Meso-Gamma-Scale Operational Weather Prediction Model, Mon. Weather Rev., 137, 3960–3977,, 2009. 

Gillespie, D. T.: An Exact Method for Numerically Simulating the Stochastic Coalescence Process in a Cloud, J. Atmos. Sci., 32, 1977–1989,<1977:AEMFNS>2.0.CO;2, 1975. 

Gillespie, D. T.: Exact stochastic simulation of coupled chemical reactions, J. Phys. Chem., 81, 2340–2361,, 1977. 

Giorgi, F. and Lionello, P.: Climate change projections for the Mediterranean region, Global Planet.Change, 63, 90–104,, 2008. 

Golaz, J.-C., Larson, V. E., and Cotton, W. R.: A PDF-Based Model for Boundary Layer Clouds. Part I: Method and Model Description, J. Atmos. Sci., 59, 3540–3551,<3540:APBMFB>2.0.CO;2, 2002a. 

Golaz, J.-C., Larson, V. E., and Cotton, W. R.: A PDF-Based Model for Boundary Layer Clouds. Part II: Model Results, J. Atmos. Sci., 59, 3552–3571,<3552:APBMFB>2.0.CO;2, 2002b. 

Gottwald, G. A., Peters, K., and Davies, L.: A data-driven method for the stochastic parametrisation of subgrid-scale tropical convective area fraction, J. Roy. Meteor. Soc., 142, 349–359,, 2016. 

Grabowski, W. W.: Coupling Cloud Processes with the Large-Scale Dynamics Using the Cloud-Resolving Convection Parameterization (CRCP), J. Atmos. Sci., 58, 978–997,<0978:CCPWTL>2.0.CO;2, 2001. 

Grabowski, W. W.: Untangling Microphysical Impacts on Deep Convection Applying a Novel Modeling Methodology, J. Atmos. Sci., 72, 2446–2464,, 2015. 

Grabowski, W. W.: Towards Global Large Eddy Simulation: Super-Parameterization Revisited, J. Meteorol. Soc. Jpn., 94, 327–344,, 2016. 

Grabowski, W. W.: Can the Impact of Aerosols on Deep Convection be Isolated from Meteorological Effects in Atmospheric Observations?, J. Atmos. Sci., 75, 3347–3363,, 2018. 

Grabowski, W. W. and Pawlowska, H.: Entrainment and Mixing in Clouds: The Paluch Mixing Diagram Revisited, J. Appl. Meteorol. Clim., 32, 1767–1773,<1767:EAMICT>2.0.CO;2, 1993. 

Grabowski, W. W. and Smolarkiewicz, P. K.: CRCP: a Cloud Resolving Convection Parameterization for modeling the tropical convecting atmosphere, Physica D: Nonlinear Phenomena, 133, 171–178,, 1999. 

Grandpeix, J.-Y. and Lafore, J.-P.: A Density Current Parameterization Coupled with Emanuel's Convection Scheme. Part I: The Models, J. Atmos. Sci., 67, 881–897,, 2010. 

Grandpeix, J.-Y., Phillips, V., and Tailleux, R.: Improved mixing representation in Emanuel's convection scheme, Q. J. Roy. Meteor. Soc., 130, 3207–3222,, 2004. 

Grant, A. L. M.: Cloud-base fluxes in the cumulus-capped boundary layer, Q. J. Roy. Meteor. Soc., 127, 407–421,, 2001. 

Grant, A. L. M. and Brown, A. R.: A similarity hypothesis for shallow-cumulus transports, Q. J. Roy. Meteor. Soc., 125, 1913–1936,, 1999. 

Grant, A. L. M. and Lock, A. P.: The turbulent kinetic energy budget for shallow cumulus convection, Q. J. Roy. Meteor. Soc., 130, 401–422,, 2004. 

Gray, M. E. B.: Characteristics of Numerically Simulated Mesoscale Convective Systems and Their Application to Parameterization, J. Atmos. Sci., 57, 3953–3970,<3953:CONSMC>2.0.CO;2, 2000. 

Gregory, D.: Estimation of entrainment rate in simple models of convective clouds, Q. J. Roy. Meteor. Soc., 127, 53–72,, 2001. 

Gregory, D. and Rowntree, P. R.: A Mass Flux Convection Scheme with Representation of Cloud Ensemble Characteristics and Stability-Dependent Closure, Mon. Weather Rev., 118, 1483–1506,<1483:AMFCSW>2.0.CO;2, 1990. 

Gregory, D., Morcrette, J.-J., Jakob, C., Beljaars, A. C. M., and Stockdale, T.: Revision of convection, radiation and cloud schemes in the ECMWF integrated forecasting system, Q. J. Roy. Meteor. Soc., 126, 1685–1710,, 2000. 

Grell, A. G., Dudhia, J., and Stauffer, D.: A description of the fifthgeneration Penn State/NCAR Mesoscale Model (MM5), University Corporation for Atmospheric Research, No. NCAR/TN-398+STR,, 1994. 

Grell, G. A.: Prognostic Evaluation of Assumptions Used by Cumulus Parameterizations, Mon. Weather Rev., 121, 764–787,<0764:PEOAUB>2.0.CO;2, 1993. 

Grell, G. A. and Dévényi, D.: A generalized approach to parameterizing convection combining ensemble and data assimilation techniques, Geophys. Res. Lett., 29, 38-1–38–4,, 2002. 

Grell, G. A. and Freitas, S. R.: A scale and aerosol aware stochastic convective parameterization for weather and air quality modeling, Atmos. Chem. Phys., 14, 5233–5250,, 2014. 

Grell, G. A., Kuo, Y.-H., and Pasch, R. J.: Semiprognostic Tests of Cumulus Parameterization Schemes in the Middle Latitudes, Mon. Weather Rev., 119, 5–31,<0005:STOCPS>2.0.CO;2, 1991. 

Grenier, H. and Bretherton, C. S.: A Moist PBL Parameterization for Large-Scale Models and Its Application to Subtropical Cloud-Topped Marine Boundary Layers, Mon. Weather Rev., 129, 357–377,<0357:AMPPFL>2.0.CO;2, 2001. 

Groenemeijer, P. and Craig, G. C.: Ensemble forecasting with a stochastic convective parametrization based on equilibrium statistics, Atmos. Chem. Phys., 12, 4555–4565,, 2012. 

Guérémy, J. F.: A continuous buoyancy based convection scheme: one-and three-dimensional validation, Tellus A, 63, 687–706,, 2011. 

Guichard, F., Petch, J. C., Redelsperger, J.-L., Bechtold, P., Chaboureau, J.-P., Cheinet, S., Grabowski, W., Grenier, H., Jones, C. G., Köhler, M., Piriou, J.-M., Tailleux, R., and Tomasini, M.: Modelling the diurnal cycle of deep precipitating convection over land with cloud-resolving models and single-column models, Q. J. Roy. Meteor. Soc., 130, 3139–3172,, 2004. 

Guo, H., Golaz, J.-C., Donner, L. J., Ginoux, P., and Hemler, R. S.: Multivariate Probability Density Functions with Dynamics in the GFDL Atmospheric General Circulation Model: Global Tests, J. Climate, 27, 2087–2108,, 2014. 

Guo, H., Golaz, J.-C., Donner, L. J., Wyman, B., Zhao, M., and Ginoux, P.: CLUBB as a unified cloud parameterization: Opportunities and challenges, Geophys. Res. Lett., 42, 4540–4547,, 2015a. 

Guo, X., Lu, C., Zhao, T., Zhang, G. J., and Liu, Y.: An Observational Study of Entrainment Rate in Deep Convection, Atmosphere-Basel, 6, 1362–1376,, 2015b. 

Gustafson, W. I., Vogelmann, A. M., Li, Z., Cheng, X., Dumas, K. K., Endo, S., Johnson, K. L., Krishna, B., Fairless, T., and Xiao, H.: The Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) Activity for Continental Shallow Convection, B. Am. Meteorol. Soc., 101, E462–E479,, 2020. 

Hack, J. J.: Parameterization of moist convection in the National Center for Atmospheric Research community climate model (CCM2), J. Geophys. Res.-Atmos., 99, 5551–5568,, 1994. 

Hack, J. J., Schubert, W. H., and Dias, P. L. S.: A Spectral Cumulus Parameterization for Use in Numerical Models of the Tropical Atmosphere, Mon. Weather Rev., 112, 704–716, 1984. 

Hagos, S., Feng, Z., Plant, R. S., Houze, R. A., and Xiao, H.: A Stochastic Framework for Modeling the Population Dynamics of Convective Clouds, J. Adv. Model. Earth Sy., 10, 448–465,, 2018. 

Han, J. and Bretherton, C. S.: TKE-Based Moist Eddy-Diffusivity Mass-Flux (EDMF) Parameterization for Vertical Turbulent Mixing, Weather Forecast., 34, 869–886,, 2019. 

Han, J. and Pan, H.-L.: Revision of Convection and Vertical Diffusion Schemes in the NCEP Global Forecast System, Weather Forecast., 26, 520–533,, 2011. 

Han, J., Witek, M. L., Teixeira, J., Sun, R., Pan, H.-L., Fletcher, J. K., and Bretherton, C. S.: Implementation in the NCEP GFS of a Hybrid Eddy-Diffusivity Mass-Flux (EDMF) Boundary Layer Parameterization with Dissipative Heating and Modified Stable Boundary Layer Mixing, Weather Forecast., 31, 341–352,, 2016a. 

Han, J., Wang, W., Kwon, Y. C., Hong, S.-Y., Tallapragada, V., and Yang, F.: Updates in the NCEP GFS Cumulus Convection Schemes with Scale and Aerosol Awareness, Weather Forecast., 32, 2005–2017,, 2017. 

Han, J.-Y., Hong, S.-Y., Lim, K.-S. S., and Han, J.: Sensitivity of a Cumulus Parameterization Scheme to Precipitation Production Representation and Its Impact on a Heavy Rain Event over Korea, Mon. Weather Rev., 144, 2125–2135,, 2016b. 

Han, J.-Y., Kim, S.-Y., Choi, I.-J., and Jin, E. K.: Effects of the Convective Triggering Process in a Cumulus Parameterization Scheme on the Diurnal Variation of Precipitation over East Asia, Atmosphere-Basel, 10, 28,, 2019. 

Han, J.-Y., Hong, S.-Y., and Kwon, Y. C.: The Performance of a Revised Simplified Arakawa–Schubert (SAS) Convection Scheme in the Medium-Range Forecasts of the Korean Integrated Model (KIM), Weather Forecast., 35, 1113–1128,, 2020. 

Hannah, W. M. and Maloney, E. D.: The Role of Moisture–Convection Feedbacks in Simulating the Madden–Julian Oscillation, J. Climate, 24, 2754–2770,, 2011. 

Hara, M., Yoshikane, T., Takahashi, H. G., Kimura, F., Noda, A., and Tokioka, T.: Assessment of the Diurnal Cycle of Precipitation over the Maritime Continent Simulated by a 20 km Mesh GCM Using TRMM PR Data, J. Meteorol. Soc. Jpn., 87A, 413–424,, 2009. 

Hararuk, O., Xia, J., and Luo, Y.: Evaluation and improvement of a global land model against soil carbon data using a Bayesian Markov chain Monte Carlo method, J. Geophys. Res.-Biogeo., 119, 403–417,, 2014. 

Heus, T. and Jonker, H. J. J.: Subsiding Shells around Shallow Cumulus Clouds, J. Atmos. Sci., 65, 1003–1018,, 2008. 

Heus, T., van Dijk, G., Jonker, H. J. J., and Akker, H. E. A. V. den: Mixing in Shallow Cumulus Clouds Studied by Lagrangian Particle Tracking, J. Atmos. Sci., 65, 2581–2597,, 2008. 

Heymsfield, A. J., Schmitt, C., and Bansemer, A.: Ice Cloud Particle Size Distributions and Pressure-Dependent Terminal Velocities from In Situ Observations at Temperatures from 0 to −86C, J. Atmos. Sci., 70, 4123–4154,, 2013. 

Hirons, L. C., Inness, P., Vitart, F., and Bechtold, P.: Understanding advances in the simulation of intraseasonal variability in the ECMWF model. Part I: The representation of the MJO, Q. J. Roy. Meteor. Soc., 139675, 1417–1426,, 2013. 

Hirota, N., Takayabu, Y. N., Watanabe, M., Kimoto, M., and Chikira, M.: Role of Convective Entrainment in Spatial Distributions of and Temporal Variations in Precipitation over Tropical Oceans, J. Climate, 27, 8707–8723,, 2014. 

Hohenegger, C. and Bretherton, C. S.: Simulating deep convection with a shallow convection scheme, Atmos. Chem. Phys., 11, 10389–10406,, 2011. 

Holden, Z. A., Swanson, A., Luce, C. H., Jolly, W. M., Maneta, M., Oyler, J. W., Warren, D. A., Parsons, R., and Affleck, D.: Decreasing fire season precipitation increased recent western US forest wildfire activity, P. Natl. Acad. Sci. USA, 115, E8349–E8357,, 2018. 

Holloway, C. E., Woolnough, S. J., and Lister, G. M. S.: Precipitation distributions for explicit versus parametrized convection in a large-domain high-resolution tropical case study, Q. J. Roy. Meteor. Soc., 138, 1692–1708,, 2012. 

Holloway, C. E., Woolnough, S. J., and Lister, G. M. S.: The Effects of Explicit versus Parameterized Convection on the MJO in a Large-Domain High-Resolution Tropical Case Study. Part I: Characterization of Large-Scale Organization and Propagation, J. ATmos. Sci., 70, 1342–1369,, 2013. 

Holtslag, A. A. M.: Modelling of atmospheric boundary layers, Royal Netherlands Academy of Arts and Sciences, 85, 110, 1998. 

Hong, S.-Y. and Pan, H.-L.: Nonlocal Boundary Layer Vertical Diffusion in a Medium-Range Forecast Model, Mon. Weather Rev., 124, 2322–2339,<2322:NBLVDI>2.0.CO;2, 1996. 

Hong, S.-Y. and Pan, H.-L.: Convective Trigger Function for a Mass-Flux Cumulus Parameterization Scheme, Mon. Weather Rev., 126, 2599–2620,<2599:CTFFAM>2.0.CO;2, 1998. 

Hong, S.-Y., Park, H., Cheong, H.-B., Kim, J.-E. E., Koo, M.-S., Jang, J., Ham, S., Hwang, S.-O., Park, B.-K., Chang, E.-C., and Li, H.: The Global/Regional Integrated Model system (GRIMs), Asia-Pac. J. Atmos. Sci., 49, 219–243,, 2013. 

Honnert, R., Efstathiou, G. A., Beare, R. J., Ito, J., Lock, A., Neggers, R., Plant, R. S., Shin, H. H., Tomassini, L., and Zhou, B.: The Atmospheric Boundary Layer and the “Gray Zone” of Turbulence: A Critical Review, J. Geophys. Res.-Atmos., 125, e2019JD030317,, 2020. 

Hou, A. Y., Kakar, R. K., Neeck, S., Azarbarzin, A. A., Kummerow, C. D., Kojima, M., Oki, R., Nakamura, K., and Iguchi, T.: The Global Precipitation Measurement Mission, B. Am. Meteorol. Soc., 95, 701–722,, 2014. 

Houghton, H. G. and Cramer, H. E.: a Theory of Entrainment in Convective Currents, J. Atmos. Sci., 8, 95–102,<0095:ATOEIC>2.0.CO;2, 1951. 

Hourdin, F., Couvreux, F., and Menut, L.: Parameterization of the Dry Convective Boundary Layer Based on a Mass Flux Representation of Thermals, J. Atmos. Sci., 59, 1105–1123,<1105:POTDCB>2.0.CO;2, 2002. 

Hourdin, F., Grandpeix, J.-Y., Rio, C., Bony, S., Jam, A., Cheruy, F., Rochetin, N., Fairhead, L., Idelkadi, A., Musat, I., Dufresne, J.-L., Lahellec, A., Lefebvre, M.-P., and Roehrig, R.: LMDZ5B: the atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection, Clim. Dynam., 40, 2193–2222,, 2013. 

Hourdin, F., Mauritsen, T., Gettelman, A., Golaz, J.-C., Balaji, V., Duan, Q., Folini, D., Ji, D., Klocke, D., Qian, Y., Rauser, F., Rio, C., Tomassini, L., Watanabe, M., and Williamson, D.: The Art and Science of Climate Model Tuning, B. Am. Meteorol. Soc., 98, 589–602,, 2017. 

Huffman, G. J., Bolvin, D. T., Braithwaite, D., Hsu, K., Joyce, R., Kidd, C., Nelkin, E. J., and Xie, P.: NASA Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals for GPM (IMERG), Algorithm Theoretical Basis Document (ATBD), 4, 26, 2015. 

Huffman, G. J., Stocker, E. F., Bolvin, D. T., Nelkin, E. J., Tan, J.: GPM IMERG Final Precipitation L3 1 day 0.1 degree x 0.1 degree V06, edited by: Savtchenko, A., Greenbelt, MD, Goddard Earth Sciences Data and Information Services Center (GES DISC) [data set], (last access: 20 November 2020), 2019. 

IPCC: Climate Change 2014: synthesis report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, 151, 2014. 

Jackson, C., Sen, M. K., and Stoffa, P. L.: An Efficient Stochastic Bayesian Approach to Optimal Parameter and Uncertainty Estimation for Climate Model Predictions, J. Climate, 17, 2828–2841,<2828:AESBAT>2.0.CO;2, 2004. 

Jakob, C.: Accelerating progress in global atmospheric model development through improved parameterizations: Challenges, opportunities, and strategies, B. Am. Meterol. Soc., 91, 869–876,, 2010. 

Jackson, C. S., Sen, M. K., Huerta, G., Deng, Y., and Bowman, K. P.: Error Reduction and Convergence in Climate Prediction, J. Climate, 21, 6698–6709,, 2008. 

Jakob, C. and Siebesma, A. P.: A New Subcloud Model for Mass-Flux Convection Schemes: Influence on Triggering, Updraft Properties, and Model Climate, Mon. Weather Rev., 131, 2765–2778,<2765:ANSMFM>2.0.CO;2, 2003. 

Jam, A., Hourdin, F., Rio, C., and Couvreux, F.: Resolved Versus Parametrized Boundary-Layer Plumes. Part III: Derivation of a Statistical Scheme for Cumulus Clouds, Boundary-Layer Meteorol, Bound.-Lay. Meteorol., 147, 421–441,, 2013. 

James, R. P. and Markowski, P. M.: A Numerical Investigation of the Effects of Dry Air Aloft on Deep Convection, Mon. Weather Rev., 138, 140–161,, 2010. 

Janjić, Z. I.: The Step-Mountain Eta Coordinate Model: Further Developments of the Convection, Viscous Sublayer, and Turbulence Closure Schemes, 122, 927–945,<0927:TSMECM>2.0.CO;2, 1994. 

Jankov, I. and Gallus, W. A.: Some contrasts between good and bad forecasts of warm season MCS rainfall, Journal of Hydrology, Mon. Weather Rev., 288, 122–152,, 2004. 

Jankov, I., Gallus, W. A., Segal, M., Shaw, B., and Koch, S. E.: The Impact of Different WRF Model Physical Parameterizations and Their Interactions on Warm Season MCS Rainfall, Weather Forecast., 20, 1048–1060,, 2005. 

Jensen, J. B., Austin, P. H., Baker, M. B., and Blyth, A. M.: Turbulent Mixing, Spectral Evolution and Dynamics in a Warm Cumulus Cloud, J. Atmos. Sci., 42, 173–192,<0173:TMSEAD>2.0.CO;2, 1985. 

Jensen, M. P. and Del Genio, A. D.: Factors Limiting Convective Cloud-Top Height at the ARM Nauru Island Climate Research Facility, J. Climate, 19, 2105–2117,, 2006. 

Jeyaratnam, J., Luo, Z. J., Giangrande, S. E., Wang, D., and Masunaga, H.: A Satellite-Based Estimate of Convective Vertical Velocity and Convective Mass Flux: Global Survey and Comparison With Radar Wind Profiler Observations, Geophys. Res. Lett., 48, e2020GL090675,, 2021. 

Jiang, H., Feingold, G., and Sorooshian, A.: Effect of Aerosol on the Susceptibility and Efficiency of Precipitation in Warm Trade Cumulus Clouds, J. Atmos. Sci., 67, 3525–3540,, 2010. 

Johnson, R. H.: The Role of Convective-Scale Precipitation Downdrafts in Cumulus and Synoptic-Scale Interactions, J. Atmos. Sci., 33, 1890–1910,<1890:TROCSP>2.0.CO;2, 1976. 

Johnson, R. H.: Diagnosis of Convective and Mesoscale Motions During Phase IH of Gate, J. Atmos. Sci., 37, 733–753,<0733:DOCAMM>2.0.CO;2, 1980. 

Jonker, H. J. J., Verzijlbergh, R. A., Heus, T., and Siebesma, A. P.: The Influence of the Sub-Cloud Moisture Field on Cloud Size Distributions and the Consequences for Entrainment, Extended Abstracts, 17th Symp. on Boundary Layers and Turbulence, 23 May 2006, SanDiego, CA, 2006. 

Joyce, R. J., Janowiak, J. E., Arkin, P. A., and Xie, P.: CMORPH:: A Method that Produces Global Precipitation Estimates from Passive Microwave and Infrared Data at High Spatial and Temporal Resolution, J. Hydrometeorol., 5, 487–503, 2004. 

Jung, J.-H. and Arakawa, A.: Modeling the moist-convective atmosphere with a Quasi-3-D Multiscale Modeling Framework (Q3D MMF), J. Adv. Model. Earth Sy., 6, 185–205,, 2014. 

Kain, J. S.: The Kain–Fritsch Convective Parameterization: An Update, J. Appl. Meteorol. Clim., 43, 170–181,<0170:TKCPAU>2.0.CO;2, 2004. 

Kain, J. S. and Fritsch, J. M.: A One-Dimensional Entraining/Detraining Plume Model and Its Application in Convective Parameterization, J. Atmos. Sci., 47, 2784–2802,<2784:AODEPM>2.0.CO;2, 1990. 

Kain, J. S. and Fritsch, J. M.: The role of the convective “trigger function” in numerical forecasts of mesoscale convective systems, Meteorol. Atmos. Phys., 49, 93–106,, 1992. 

Kain, J. S. and Fritsch, J. M.: Convective Parameterization for Mesoscale Models: The Kain-Fritsch Scheme, in: The Representation of Cumulus Convection in Numerical Models, Meteorological Monographs, American Meteorological Society, 246 pp.,, 1993. 

Kain, J. S., Weiss, S. J., Levit, J. J., Baldwin, M. E., and Bright, D. R.: Examination of Convection-Allowing Configurations of the WRF Model for the Prediction of Severe Convective Weather: The SPC/NSSL Spring Program 2004, Weather Forecast., 21, 167–181,, 2006. 

Karlický, J., Huszár, P., Nováková, T., Belda, M., Švábik, F., Ďoubalová, J., and Halenka, T.: The “urban meteorology island”: a multi-model ensemble analysis, Atmos. Chem. Phys., 20, 15061–15077,, 2020. 

Kawecki, S., Henebry, G. M., and Steiner, A. L.: Effects of Urban Plume Aerosols on a Mesoscale Convective System, J. Atmos. Sci., 73, 4641–4660,, 2016. 

Keane, R. J., Craig, G. C., Keil, C., and Zängl, G.: The Plant–Craig Stochastic Convection Scheme in ICON and Its Scale Adaptivity, J. Atmos. Sci., 71, 3404–3415,, 2014. 

Kendon, E. J., Roberts, N. M., Senior, C. A., and Roberts, M. J.: Realism of Rainfall in a Very High-Resolution Regional Climate Model, J. Climate, 25, 5791–5806,, 2012. 

Kessler, E.: On the Distribution and Continuity of Water Substance in Atmospheric Circulations, in: On the Distribution and Continuity of Water Substance in Atmospheric Circulations, Meteorological Monographs, vol. 10, American Meteorological Society,, 1969. 

Khain, A., Rosenfeld, D., and Pokrovsky, A.: Aerosol impact on the dynamics and microphysics of deep convective clouds, Q. J. Roy. Meteor. Soc., 131, 2639–2663,, 2005. 

Khairoutdinov, M. and Randall, D.: High-Resolution Simulation of Shallow-to-Deep Convection Transition over Land, J. Atmos. Sci., 63, 3421–3436,, 2006. 

Khairoutdinov, M., Randall, D., and DeMott, C.: Simulations of the Atmospheric General Circulation Using a Cloud-Resolving Model as a Superparameterization of Physical Processes, J. Atmos. Sci., 62, 2136–2154,, 2005. 

Khairoutdinov, M. F. and Randall, D. A.: Cloud Resolving Modeling of the ARM Summer 1997 IOP: Model Formulation, Results, Uncertainties, and Sensitivities, J. Atmos. Sci., 60, 607–625,<0607:CRMOTA>2.0.CO;2, 2003. 

Khouider, B.: A coarse grained stochastic multi-type particle interacting model for tropical convection: Nearest neighbour interactions, Comm. Math. Sci., 12, 1379–1407,, 2014. 

Khouider, B. and Majda, A.: Multicloud Models for Organized Tropical Convection: Enhanced Congestus Heating, J. Atmos. Sci., 65, 895–914,, 2008. 

Khouider, B. and Majda, A. J.: A Simple Multicloud Parameterization for Convectively Coupled Tropical Waves. Part I: Linear Analysis, J. Atmos. Sci., 63, 1308–1323,, 2006. 

Khouider, B. and Moncrieff, M. W.: Organized Convection Parameterization for the ITCZ, J. Atmos. Sci., 72, 3073–3096,, 2015. 

Khouider, B., Majda, A. J., and Katsoulakis, M. A.: Coarse-grained stochastic models for tropical convection and climate, P. Natl. Acad. Sci. USA, 100, 11941–11946,, 2003. 

Khouider, B., Biello, J., and Majda, A. J.: A stochastic multicloud model for tropical convection, Comm. Math. Sci., 8, 187–216, 2010. 

Kim, D. and Kang, I.-S.: A bulk mass flux convection scheme for climate model: description and moisture sensitivity, Clim. Dynam., 38, 411–429,, 2012. 

Kim, D., Sobel, A. H., Maloney, E. D., Frierson, D. M. W., and Kang, I.-S.: A Systematic Relationship between Intraseasonal Variability and Mean State Bias in AGCM Simulations, J. Climate, 24, 5506–5520,, 2011. 

Kim, D., Sobel, A. H., Del Genio, A. D., Chen, Y., Camargo, S. J., Yao, M.-S., Kelley, M., and Nazarenko, L.: The Tropical Subseasonal Variability Simulated in the NASA GISS General Circulation Model, J. Climate, 25, 4641–4659,, 2012. 

Kim, D., Del Genio, A. D., and Yao, M.-S.: Moist convection scheme in Model E2, arXiv preprint arXiv:1312.7496, 2013. 

Kirshbaum, D. J. and Grant, A. L. M.: Invigoration of cumulus cloud fields by mesoscale ascent, Q. J. Roy. Meteor. Soc., 138, 2136–2150,, 2012. 

Kirshbaum, D. J. and Lamer, K.: Climatological Sensitivities of Shallow-Cumulus Bulk Entrainment in Continental and Oceanic Locations, J. Atmos. Sci., 78, 2429–2443,, 2021. 

Klein, S. A. and Hartmann, D. L.: The Seasonal Cycle of Low Stratiform Clouds, J. Climate, 6, 1587–1606,<1587:TSCOLS>2.0.CO;2, 1993. 

Klingaman, N. P. and Woolnough, S. J.: Using a case-study approach to improve the Madden–Julian oscillation in the Hadley Centre model, Q. J. Roy. Meteor. Soc., 140, 2491–2505,, 2014. 

Klocke, D., Pincus, R., and Quaas, J.: On Constraining Estimates of Climate Sensitivity with Present-Day Observations through Model Weighting, J. Climate, 24, 6092–6099,, 2011. 

Knievel, J. C., Ahijevych, D. A., and Manning, K. W.: Using Temporal Modes of Rainfall to Evaluate the Performance of a Numerical Weather Prediction Model, Mon. Weather Rev., 132, 2995–3009,, 2004. 

Köhler, M.: Improved prediction of boundary layer clouds, ECMWF Newsletter, 104, 18–22, 2005. 

Köhler, M., Ahlgrimm, M., and Beljaars, A.: Unified treatment of dry convective and stratocumulus-topped boundary layer in the ECMWF model, Q. J. Roy. Meteor. Soc.,, 137, 43–57,, 2011. 

Kooperman, G. J., Pritchard, M. S., O'Brien, T. A., and Timmermans, B. W.: Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present-Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model, J. Adv. Model. Earth Sy., 10, 971–988,, 2018. 

Koren, I., Kaufman, Y. J., Rosenfeld, D., Remer, L. A., and Rudich, Y.: Aerosol invigoration and restructuring of Atlantic convective clouds, Geophys. Res. Lett., 32,, 2005. 

Kreitzberg, C. W. and Perkey, D. J.: Release of Potential Instability: Part I. A Sequential Plume Model within a Hydrostatic Primitive Equation Model, J. Atmos. Sci., 33, 456–475,<0456:ROPIPI>2.0.CO;2, 1976. 

Krishnamurthy, V. and Stan, C.: Simulation of the South American climate by a coupled model with super-parameterized convection, Clim. Dynam., 44, 2369–2382,, 2015. 

Krishnamurti, T. N., Ramanathan, Y., Pan, H.-L., Pasch, R. J., and Molinari, J.: Cumulus Parameterization and Rainfall Rates I, Mon. Weather Rev., 108, 465–472,<0465:CPARRI>2.0.CO;2, 1980. 

Krishnamurti, T. N., Low-Nam, S., and Pasch, R.: Cumulus Parameterization and Rainfall Rates II, Mon. Weather Rev., 111, 815–828,<0815:CPARRI>2.0.CO;2, 1983. 

Krueger, S. K.: Numerical Simulation of Tropical Cumulus Clouds and Their Interaction with the Subcloud Layer, J. Atmos. Sci., 45, 2221–2250,<2221:NSOTCC>2.0.CO;2, 1988. 

Kuang, Z.: Modeling the Interaction between Cumulus Convection and Linear Gravity Waves Using a Limited-Domain Cloud System–Resolving Model, J. Atmos. Sci., 65, 576–591,, 2008. 

Kuang, Z. and Bretherton, C. S.: A Mass-Flux Scheme View of a High-Resolution Simulation of a Transition from Shallow to Deep Cumulus Convection, J. Atmos. Sci., 63, 1895–1909,, 2006. 

Kucera, P. A., Ebert, E. E., Turk, F. J., Levizzani, V., Kirschbaum, D., Tapiador, F. J., Loew, A., and Borsche, M.: Precipitation from Space: Advancing Earth System Science, B. Am. Meteorol. Soc., 94, 365–375,, 2013. 

Kuell, V., Gassmann, A., and Bott, A.: Towards a new hybrid cumulus parametrization scheme for use in non-hydrostatic weather prediction models, Q. J. Roy. Meteor. Soc., 133, 479–490,, 2007. 

Kumar, B., Götzfried, P., Suresh, N., Schumacher, J., and Shaw, R. A.: Scale Dependence of Cloud Microphysical Response to Turbulent Entrainment and Mixing, J. Adv. Model. Earth Sy.,10, 2777–2785,, 2018. 

Kumar, D. and Dimri, A. P.: Sensitivity of convective and land surface parameterization in the simulation of contrasting monsoons over CORDEX-South Asia domain using RegCM-, Theor. Appl. Climatol., 139, 297–322,, 2020. 

Kummerow, C., Barnes, W., Kozu, T., Shiue, J., and Simpson, J.: The Tropical Rainfall Measuring Mission (TRMM) Sensor Package, J. Atmos. Ocean. Tech., 15, 809–817,<0809:TTRMMT>2.0.CO;2, 1998. 

Kuo, H. L.: On the Controlling Influences of Eddy Diffusion on Thermal Convection, J. Atmos. Sci., 19, 236–243,<0236:OTCIOE>2.0.CO;2, 1962. 

Kuo, H. L.: On Formation and Intensification of Tropical Cyclones Through Latent Heat Release by Cumulus Convection, J. Atmos. Sci., 22, 40–63,<0040:OFAIOT>2.0.CO;2, 1965. 

Kuo, H. L.: Further Studies of the Parameterization of the Influence of Cumulus Convection on Large-Scale Flow, J. Atmos. Sci., 31, 1232–1240,<1232:FSOTPO>2.0.CO;2, 1974. 

Kuo, Y.-H. and Anthes, R. A.: Semiprognostic Tests of Kuo–Type Cumulus Parameterization Schemes in an Extratropical Convective System, Mon. Weather Rev., 112, 1498–1509,<1498:STOKCP>2.0.CO;2, 1984. 

Kurihara, Y.: A Scheme of Moist Convective Adjustment, Mon. Weather Rev., 101, 547–553,<0547:ASOMCA>2.3.CO;2, 1973. 

Kurowski, M. J., Thrastarson, H. T., Suselj, K., and Teixeira, J.: Towards unifying the planetary boundary layer and shallow convection in CAM5 with the eddy-diffusivity/mass-flux approach, Atmosphere-Basel, 10, 484,, 2019. 

Kwon, Y. C. and Hong, S.-Y.: A Mass-Flux Cumulus Parameterization Scheme across Gray-Zone Resolutions, Mon. Weather Rev., 145, 583–598,, 2017. 

Lamontagne, R. G. and Telford, J. W.: Cloud Top Mixing in Small Cumuli., Journal of Atmospheric Sciences, J. Atmos. Sci., 40, 2148–2156,<2148:CTMISC>2.0.CO;2, 1983. 

Lappen, C.-L. and Randall, D. A.: Toward a Unified Parameterization of the Boundary Layer and Moist Convection. Part I: A New Type of Mass-Flux Model, J. Atmos. Sci., 58, 2021–2036,<2021:TAUPOT>2.0.CO;2, 2001a. 

Lappen, C.-L. and Randall, D. A.: Toward a Unified Parameterization of the Boundary Layer and Moist Convection. Part II: Lateral Mass Exchanges and Subplume-Scale Fluxes, J. Atmos. Sci., 58, 2037–2051,<2037:TAUPOT>2.0.CO;2, 2001b. 

Larson, V. E.: CLUBB-SILHS: A parameterization of subgrid variability in the atmosphere, arXiv preprint arXiv:1711.03675, 2020. 

Larson, V. E. and Schanen, D. P.: The Subgrid Importance Latin Hypercube Sampler (SILHS): a multivariate subcolumn generator, Geosci. Model Dev., 6, 1813–1829,, 2013. 

Larson, V. E., Golaz, J.-C., and Cotton, W. R.: Small-Scale and Mesoscale Variability in Cloudy Boundary Layers: Joint Probability Density Functions, J. Atmos. Sci., 59, 3519–3539,<3519:SSAMVI>2.0.CO;2, 2002. 

Larson, V. E., Golaz, J.-C., Jiang, H., and Cotton, W. R.: Supplying Local Microphysics Parameterizations with Information about Subgrid Variability: Latin Hypercube Sampling, J. Atmos. Sci., 62, 4010–4026,, 2005. 

Larson, V. E., Schanen, D. P., Wang, M., Ovchinnikov, M., and Ghan, S.: PDF Parameterization of Boundary Layer Clouds in Models with Horizontal Grid Spacings from 2 to 16 km, Mon. Weather Rev., 140, 285–306,, 2012. 

Le Trent, H. and Li, Z.-X.: Sensitivity of an atmospheric general circulation model to prescribed SST changes: feedback effects associated with the simulation of cloud optical properties, Clim. Dynam., 5, 175–187,, 1991. 

Leary, C. A. and Houze, R. A.: The Contribution of Mesoscale Motions to the Mass and Heat Fluxes of an Intense Tropical Convective System, J. Atmos. Sci., 37, 784–796,<0784:TCOMMT>2.0.CO;2, 1980. 

Lee, M.-I., Schubert, S. D., Suarez, M. J., Held, I. M., Lau, N.-C., Ploshay, J. J., Kumar, A., Kim, H.-K., and Schemm, J.-K. E.: An Analysis of the Warm-Season Diurnal Cycle over the Continental United States and Northern Mexico in General Circulation Models, J. Hydormeteorol., 8, 344–366,, 2007a. 

Lee, M.-I., Schubert, S. D., Suarez, M. J., Held, I. M., Kumar, A., Bell, T. L., Schemm, J.-K. E., Lau, N.-C., Ploshay, J. J., Kim, H.-K., and Yoo, S.-H.: Sensitivity to Horizontal Resolution in the AGCM Simulations of Warm Season Diurnal Cycle of Precipitation over the United States and Northern Mexico, J. Climate, 20, 1862–1881,, 2007b. 

Lee, M.-I., Schubert, S. D., Suarez, M. J., Schemm, J.-K. E., Pan, H.-L., Han, J., and Yoo, S.-H.: Role of convection triggers in the simulation of the diurnal cycle of precipitation over the United States Great Plains in a general circulation model, J. Geophys. Res.-Atmos., 113, D02111,, 2008. 

Lee, Y. H., Park, S. K., and Chang, D.-E.: Parameter estimation using the genetic algorithm and its impact on quantitative precipitation forecast, Ann. Geophys., 24, 3185–3189,, 2006. 

LeMone, M. A. and Pennell, W. T.: The Relationship of Trade Wind Cumulus Distribution to Subcloud Layer Fluxes and Structure, Mon. Weather Rev., 104, 524–539,<0524:TROTWC>2.0.CO;2, 1976. 

Levizzani, V. and Cattani, E.: Satellite Remote Sensing of Precipitation and the Terrestrial Water Cycle in a Changing Climate, Remote Sens.-Basel, 11, 2301,, 2019. 

Lewellen, W. S. and Yoh, S.: Binormal Model of Ensemble Partial Cloudiness, J. Atmos. Sci., 50, 1228–1237,<1228:BMOEPC>2.0.CO;2, 1993. 

Li, L., Wang, B., Yuqing, W., and Hui, W.: Improvements in climate simulation with modifications to the Tiedtke convective parameterization in the grid-point atmospheric model of IAP LASG (GAMIL), Adv. Atmos. Sci., 24, 323–335,, 2007. 

Li, S., Zhang, S., Liu, Z., Lu, L., Zhu, J., Zhang, X., Wu, X., Zhao, M., Vecchi, G. A., Zhang, R.-H., and Lin, X.: Estimating Convection Parameters in the GFDL CM2.1 Model Using Ensemble Data Assimilation, J. Adv. Model. Earth Sy., 10, 989–1010,, 2018. 

Liang, F., Cheng, Y., and Lin, G.: Simulated Stochastic Approximation Annealing for Global Optimization With a Square-Root Cooling Schedule, J. Am. Stat. Assoc., 109, 847–863,, 2014. 

Lim, K.-S. S., Hong, S.-Y., Yoon, J.-H., and Han, J.: Simulation of the Summer Monsoon Rainfall over East Asia Using the NCEP GFS Cumulus Parameterization at Different Horizontal Resolutions, Weather Forecast., 29, 1143–1154,, 2014. 

Lin, J. W.-B. and Neelin, J. D.: Influence of a stochastic moist convective parameterization on tropical climate variability, Geophys. Res. Lett., 27, 3691–3694,, 2000. 

Lin, J. W.-B. and Neelin, J. D.: Considerations for Stochastic Convective Parameterization, J. Atmos. Sci., 59, 959–975,<0959:CFSCP>2.0.CO;2, 2002. 

Lin, J. W.-B. and Neelin, J. D.: Toward stochastic deep convective parameterization in general circulation models, Geophys. Res. Lett., 30, 1162,, 2003. 

Lin, J.-L., Kiladis, G. N., Mapes, B. E., Weickmann, K. M., Sperber, K. R., Lin, W., Wheeler, M. C., Schubert, S. D., Genio, A. D., Donner, L. J., Emori, S., Gueremy, J.-F., Hourdin, F., Rasch, P. J., Roeckner, E., and Scinocca, J. F.: Tropical Intraseasonal Variability in 14 IPCC AR4 Climate Models. Part I: Convective Signals, J. Climate, 19, 2665–2690,, 2006. 

Lin, J.-L., Lee, M.-I., Kim, D., Kang, I.-S., and Frierson, D. M. W.: The Impacts of Convective Parameterization and Moisture Triggering on AGCM-Simulated Convectively Coupled Equatorial Waves, J. Climate, 21, 883–909,, 2008. 

Lin, J.-L., Qian, T., Shinoda, T., and Li, S.: Is the Tropical Atmosphere in Convective Quasi-Equilibrium?, J. Climate, 28, 4357–4372,, 2015. 

Lindzen, R. S.: Some remarks on cumulus parameterization, Pure Appl. Geophys., 126, 123–135,, 1988. 

Lindzen, R. S., Chou, M.-D., and Hou, A. Y.: Does the Earth Have an Adaptive Infrared Iris?, B. Am. Meteorol. Soc., 82, 417–432,<0417:DTEHAA>2.3.CO;2, 2001. 

Liu, C., Fedorovich, E., Huang, J., Hu, X.-M., Wang, Y., and Lee, X.: Impact of Aerosol Shortwave Radiative Heating on Entrainment in the Atmospheric Convective Boundary Layer: A Large-Eddy Simulation Study, J. Atmos. Sci., 76, 785–799,, 2019. 

Lohmann, U.: Global anthropogenic aerosol effects on convective clouds in ECHAM5-HAM, Atmos. Chem. Phys., 8, 2115–2131,, 2008. 

Lord, S. J., Chao, W. C., and Arakawa, A.: Interaction of a Cumulus Cloud Ensemble with the Large-Scale Environment. Part IV: The Discrete Model, J. Atmos. Sci., 39, 104–113,<0104:IOACCE>2.0.CO;2, 1982. 

Loriaux, J. M., Lenderink, G., Roode, S. R. D., and Siebesma, A. P.: Understanding Convective Extreme Precipitation Scaling Using Observations and an Entraining Plume Model, J. Atmos. Sci., 70, 3641–3655,, 2013. 

Lotka, A. J.: Contribution to the Theory of Periodic Reactions, J. Phys. Chem., 14, 271-274,, 1910. 

Lotka, A. J.: Analytical Note on Certain Rhythmic Relations in Organic Systems, P. Natl. Acad. Sci. USA, 6, 410–415,, 1920. 

Louis, J.-F.: A parametric model of vertical eddy fluxes in the atmosphere, Bound.-Lay. Meteorol, 17, 187–202,, 1979. 

Lu, B. and Ren, H.-L.: Improving ENSO periodicity simulation by adjusting cumulus entrainment in BCC_CSMs, Dynam. Atmos. Oceans, 76, 127–140,, 2016. 

Lu, C., Liu, Y., and Niu, S.: Examination of turbulent entrainment-mixing mechanisms using a combined approach, J. Geophys. Res.-Atmos., 116, D20207,, 2011. 

Lu, C., Liu, Y., Yum, S. S., Niu, S., and Endo, S.: A new approach for estimating entrainment rate in cumulus clouds, Geophys. Res. Lett., 39, L04802,, 2012. 

Lu, C., Liu, Y., Niu, S., and Endo, S.: Scale dependence of entrainment-mixing mechanisms in cumulus clouds, J. Geophys. Res.-Atmos., 119, 13877–13890,, 2014. 

Lu, C., Sun, C., Liu, Y., Zhang, G. J., Lin, Y., Gao, W., Niu, S., Yin, Y., Qiu, Y., and Jin, L.: Observational Relationship Between Entrainment Rate and Environmental Relative Humidity and Implications for Convection Parameterization, Geophys. Res. Lett., 45, 13495–13504,, 2018. 

Luo, Z. J., Liu, G. Y., and Stephens, G. L.: Use of A-Train data to estimate convective buoyancy and entrainment rate, Geophys. Res. Lett., 37, L09804,, 2010. 

Ma, L.-M. and Tan, Z.-M.: Improving the behavior of the cumulus parameterization for tropical cyclone prediction: Convection trigger, Atmospheric Research, Atmos. Res., 92, 190–211,, 2009. 

Mahoney, K. M.: The representation of cumulus convection in high-resolution simulations of the 2013 Colorado Front Range flood, Mon. Weather Rev., 144, 4265–4278, 2016. 

Majda, A. J. and Khouider, B.: Stochastic and mesoscopic models for tropical convection, P. Natl. Acad. Sci. USA, 99, 1123–1128,, 2002. 

Majda, A. J., Timofeyev, I., and Eijnden, E. V.: Models for stochastic climate prediction, P. Natl. Acad. Sci. USA, 96, 14687–14691,, 1999. 

Majda, A. J., Timofeyev, I., and Eijnden, E. V.: A mathematical framework for stochastic climate models, Commun. Pur. Appl. Math., 54, 891–974,, 2001. 

Majda, A. J., Timofeyev, I., and Vanden-Eijnden, E.: Systematic Strategies for Stochastic Mode Reduction in Climate, J. Atmos. Sci., 60, 1705–1722,<1705:SSFSMR>2.0.CO;2, 2003. 

Malinowski, S. P. and Pawlowska-Mankiewicz, H.: On Estimating the Entraininent Level in Cumulus Clouds, J. Atmos. Sci., 46, 2463–2465,<2463:OETELI>2.0.CO;2, 1989. 

Malkus, J. S.: Recent developments in studies of penetrative convection and an application to hurricane cumulonimbus towers, Cumulus Dynamics: Proceedings First Conference on Cumulus Convection, 19–22 May 1959, Portsmouth, N.H., edited by: Anderson, C. E., Pergamon Press, London, New York, 65–84, 1960. 

Manabe, S., Smagorinsky, J., and Strickler, R. F.: Simulated Climatology of a General Circulation Model with Hydrologic Cycle, Mon. Weather Rev., 93, 769–798,<0769:SCOAGC>2.3.CO;2, 1965. 

Mapes, B. and Neale, R.: Parameterizing Convective Organization to Escape the Entrainment Dilemma, J. Adv. Model. Earth Sy., 3, M06004,, 2011. 

Mapes, B. E.: Equilibrium Vs. Activation Control of Large-Scale Variations of Tropical Deep Convection, in: The Physics and Parameterization of Moist Atmospheric Convection, edited by: Smith, R. K., Springer Netherlands, Dordrecht, 321–358,, 1997. 

Mapes, B. E.: Convective Inhibition, Subgrid-Scale Triggering Energy, and Stratiform Instability in a Toy Tropical Wave Model, J. Atmos. Sci., 57, 1515–1535,<1515:CISSTE>2.0.CO;2, 2000. 

Mauritsen, T., Stevens, B., Roeckner, E., Crueger, T., Esch, M., Giorgetta, M., Haak, H., Jungclaus, J., Klocke, D., Matei, D., Mikolajewicz, U., Notz, D., Pincus, R., Schmidt, H., and Tomassini, L.: Tuning the climate of a global model, J. Adv. Model. Earth Sy., 4, M00A01,, 2012. 

Mbienda, A. J. K., Tchawoua, C., Vondou, D. A., Choumbou, P., Sadem, C. K., and Dey, S.: Sensitivity experiments of RegCM4 simulations to different convective schemes over Central Africa, Int. J. Climatol., 37, 328–342,, 2017. 

McCaa, J. R. and Bretherton, C. S.: A New Parameterization for Shallow Cumulus Convection and Its Application to Marine Subtropical Cloud-Topped Boundary Layers. Part II: Regional Simulations of Marine Boundary Layer Clouds, Mon. Weather Rev., 132, 883–896,<0883:ANPFSC>2.0.CO;2, 2004. 

McFarlane, N.: Parameterizations: representing key processes in climate models without resolving them, WIRES CLim. Change, 2, 482–497,, 2011. 

McFiggans, G., Artaxo, P., Baltensperger, U., Coe, H., Facchini, M. C., Feingold, G., Fuzzi, S., Gysel, M., Laaksonen, A., Lohmann, U., Mentel, T. F., Murphy, D. M., O'Dowd, C. D., Snider, J. R., and Weingartner, E.: The effect of physical and chemical aerosol properties on warm cloud droplet activation, Atmos. Chem. Phys., 6, 2593–2649,, 2006. 

McGranahan, G., Balk, D., and Anderson, B.: The rising tide: assessing the risks of climate change and human settlements in low elevation coastal zones, Environ. Urban., 19, 17–37,, 2007. 

McLaughlin, J. F., Hellmann, J. J., Boggs, C. L., and Ehrlich, P. R.: Climate change hastens population extinctions, P. Natl. Acad. Sci. USA, 99, 6070–6074,, 2002. 

Mellor, G. L.: The Gaussian Cloud Model Relations, J. Atmos. Sci., 34, 356–358,<0356:TGCMR> 2.0.CO;2, 1977. 

Mironov, D. V.: Turbulence in the Lower Troposphere: Second-Order Closure and Mass–Flux Modelling Frameworks, in Interdisciplinary Aspects of Turbulence, Lect. Notes Phys., Springer, Berlin, Heidelberg, 161–221,, 2009. 

Miyakoda, K., Smagorinsky, J., Strickler, R. F., and Hembree, G. D.: Experimental predictions with a nine-level hemispheric model, Mon. Weather Rev., 97, 1–76,<0001:EEPWAN>2.3.CO;2, 1969. 

Möbis, B. and Stevens, B.: Factors controlling the position of the Intertropical Convergence Zone on an aquaplanet, J. Adv. Model. Earth Sy., 4, M00A04,, 2012. 

Mohandas, S. and Ashrit, R.: Sensitivity of different convective parameterization schemes on tropical cyclone prediction using a mesoscale model, Nat. Hazards, 73, 213–235,, 2014. 

Molinari, J.: A General Form of Kuo's Cumulus Parameterization, Mon. Weather Rev., 113, 1411–1416,<1411:AGFOKC>2.0.CO;2, 1985. 

Molinari, J. and Corsetti, T.: Incorporation of Cloud-Scale and Mesoscale Downdrafts into a Cumulus Parameterization: Results of One- and Three-Dimensional Integrations, Mon. Weather Rev., 113, 485–501,<0485:IOCSAM>2.0.CO;2, 1985. 

Moncrieff, M. W. and Liu, C.: Representing convective organization in prediction models by a hybrid strategy, J. Atmos. Sci., 63, 3404–3420,, 2006. 

Moncrieff, M. W., Liu, C., and Bogenschutz, P.: Simulation, Modeling, and Dynamically Based Parameterization of Organized Tropical Convection for Global Climate Models, J. Atmos. Sci., 74, 1363–1380,, 2017. 

Moorthi, S. and Suarez, M. J.: Relaxed Arakawa-Schubert. A Parameterization of Moist Convection for General Circulation Models, Mon. Weather Rev., 120, 978–1002,<0978:RASAPO>2.0.CO;2, 1992. 

Morrison, H.: Impacts of updraft size and dimensionality on the perturbation pressure and vertical velocity in cumulus convection. Part I: Simple, generalized analytic solutions, J. Atmos. Sci., 73, 1441–1454,, 2016a. 

Morrison, H.: Impacts of updraft size and dimensionality on the perturbation pressure and vertical velocity in cumulus convection. Part II: Comparison of theoretical and numerical solutions and fully dynamical simulations, J. Atmos. Sci., 73, 1455–1480,, 2016. 

Morrison, H. and Grabowski, W. W.: Response of Tropical Deep Convection to Localized Heating Perturbations: Implications for Aerosol-Induced Convective Invigoration, J. Atmos. Sci., 70, 3533–555,, 2013. 

Morton, B. R.: Modeling fire plumes, Symposium (International) on Combustion, 10, 973–982,, 1965. 

Morton, B. R., Taylor, G. I., and Turner, J. S.: Turbulent gravitational convection from maintained and instantaneous sources, P. Roy. Soc. Lond. A Mat., 234, 1–23,, 1956. 

Mukhopadhyay, P., Taraphdar, S., Goswami, B. N., and Krishnakumar, K.: Indian Summer Monsoon Precipitation Climatology in a High-Resolution Regional Climate Model: Impacts of Convective Parameterization on Systematic Biases, Weather Forecast., 25, 369–387,, 2010. 

Nam, C. C. W., Quaas, J., Neggers, R., Drian, C. S.-L., and Isotta, F.: Evaluation of boundary layer cloud parameterizations in the ECHAM5 general circulation model using CALIPSO and CloudSat satellite data, J. Adv. Model. Earth Sy., 6, 300–314,, 2014. 

National Academies of Sciences, Engineering and Medicine: Thriving on Our Changing Planet: A Decadal Strategy for Earth Observation from Space, The National Academies Press, Washington, D.C.,, 2018. 

National Centers for Environmental Prediction/National Weather Service/NOAA/U.S. Department of Commerce: NCEP GFS 0.25 Degree Global Forecast Grids Historical Archive, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory [data set], (last access: 22 November 2020), 2015, updated daily. 

Naumann, A. K., Seifert, A., and Mellado, J. P.: A refined statistical cloud closure using double-Gaussian probability density functions, Geosci. Model Dev., 6, 1641–1657,, 2013. 

Neale, R. B., Richter, J. H., and Jochum, M.: The Impact of Convection on ENSO: From a Delayed Oscillator to a Series of Events, J. Climate, 21, 5904–5924,, 2008. 

Neggers, R.: Humidity-convection feedbacks in a mass flux scheme based on resolved size densities, ECMWF Workshop on Parametrization of Clouds and Precipitation, 10, (last access: 5 September 2021), 2012. 

Neggers, R. A. J.: A Dual Mass Flux Framework for Boundary Layer Convection. Part II: Clouds, J. Atmos. Sci., 66, 1489–1506,, 2009. 

Neggers, R. A. J.: Exploring bin-macrophysics models for moist convective transport and clouds, J. Adv. Model. Earth Sy., 7, 2079–2104,, 2015. 

Neggers, R. A. J. and Griewank, P. J.: A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations, J. Adv. Model. Earth Sy., 13, e2020MS002229,, 2021. 

Neggers, R. A. J. and Siebesma, A. P.: Constraining a System of Interacting Parameterizations through Multiple-Parameter Evaluation: Tracing a Compensating Error between Cloud Vertical Structure and Cloud Overlap, J. Climate, 26, 6698–6715,, 2013. 

Neggers, R. A. J., Siebesma, A. P., and Jonker, H. J. J.: A Multiparcel Model for Shallow Cumulus Convection, J. Atmos. Sci., 59, 1655–1668,<1655:AMMFSC>2.0.CO;2, 2002. 

Neggers, R. A. J., Jonker, H. J. J., and Siebesma, A. P.: Size Statistics of Cumulus Cloud Populations in Large-Eddy Simulations, J. Atmos. Sci., 60, 1060–1074,<1060:SSOCCP>2.0.CO;2, 2003. 

Neggers, R. A. J., Siebesma, A. P., Lenderink, G., and Holtslag, A. A. M.: An Evaluation of Mass Flux Closures for Diurnal Cycles of Shallow Cumulus, Mon. Weather Rev., 132, 2525–2538,, 2004. 

Neggers, R. A. J., Stevens, B., and Neelin, J. D.: Variance scaling in shallow-cumulus-topped mixed layers, Q. J. Roy. Meteor. Soc., 133, 1629–1641,, 2007. 

Neggers, R. A. J., Köhler, M., and Beljaars, A. C. M.: A Dual Mass Flux Framework for Boundary Layer Convection. Part I: Transport, J. Atmos. Sci., 66, 1465–1487,, 2009. 

Neggers, R. A. J., Siebesma, A. P., and Heus, T.: Continuous Single-Column Model Evaluation at a Permanent Meteorological Supersite, B. Am. Meteorol. Soc., 93, 1389–1400,, 2012. 

Neggers, R. A. J., Griewank, P. J., and Heus, T.: Power-Law Scaling in the Internal Variability of Cumulus Cloud Size Distributions due to Subsampling and Spatial Organization, J. Atmos. Sci., 76, 1489–1503,, 2019. 

Nie, J. and Kuang, Z.: Responses of Shallow Cumulus Convection to Large-Scale Temperature and Moisture Perturbations: A Comparison of Large-Eddy Simulations and a Convective Parameterization Based on Stochastically Entraining Parcels, J. Atmos. Sci., 69, 1936–1956,, 2012. 

Nitta, T.: Observational Determination of Cloud Mass Flux Distributions, J. Atmos. Sci., 32, 73–91,<0073:ODOCMF>2.0.CO;2, 1975. 

Niziol, T. A., Snyder, W. R., and Waldstreicher, J. S.: Winter Weather Forecasting throughout the Eastern United States. Part IV: Lake Effect Snow, Weather Forecast., 10, 61–77,<0061:WWFTTE>2.0.CO;2, 1995. 

Nober, F. J. and Graf, H. F.: A new convective cloud field model based on principles of self-organisation, Atmos. Chem. Phys., 5, 2749–2759,, 2005. 

Nober, F. J., Graf, H.-F., and Rosenfeld, D.: Sensitivity of the global circulation to the suppression of precipitation by anthropogenic aerosols, Global Planet. Change, 37, 57–80,, 2003. 

Nordeng, T.-E.: Extended versions of the convective parametrization scheme at ECMWF and their impact on the mean and transient activity of the model in the tropics, Research Department Technical Memorandum, no. 206, 1994. 

O'Gorman, P. A. and Dwyer, J. G.: Using Machine Learning to Parameterize Moist Convection: Potential for Modeling of Climate, Climate Change, and Extreme Events, J. Adv. Model Earth Sy., 10, 2548–2563,, 2018. 

Okamoto, K. I., Ushio, T., Iguchi, T., Takahashi, N., and Iwanami, K.: The global satellite mapping of precipitation (GSMaP) project, Aqua (AMSR-E), 3414–3416,, 2005. 

Olson, J., Kenyon, J., Angevine, W. A., Brown, J. M., Pagowski, M., and Sušelj, K. (Eds.): A Description of the MYNN-EDMF Scheme and the Coupling to Other Components in WRF–ARW, NOAA Technical Memorandum OAR GSD, 61,, 2019. 

Ooyama, K.: A dynamical model for the study of tropical cyclone development., Geofis. Int., 4, 187–198, 1964. 

Ooyama, K.: A Theory on Parameterization of Cumulus Convection, J. Meteorol. Soc. Jpn., 49A, 744–756,, 1971. 

Oueslati, B. and Bellon, G.: Convective Entrainment and Large-Scale Organization of Tropical Precipitation: Sensitivity of the CNRM-CM5 Hierarchy of Models, J. Climate, 26, 2931–2946,, 2013. 

Paluch, I. R.: The Entrainment Mechanism in Colorado Cumuli, J. Atmos. Sci., 36, 2467–2478,<2467:TEMICC>2.0.CO;2, 1979. 

Pan, D.-M. and Randall, D. D. A.: A cumulus parameterization with a prognostic closure, Q. J. Roy. Meteor. Soc., 124, 949–981,, 1998. 

Pan, H.-L. and Wu, W.-S.: Implementing a mass flux convection parameterization package for the NMC medium-range forecast model, National Centers for Environmental Prediction (U.S.), Office note (National Centers for Environmental Prediction (U.S.)), 409, (last access: 5 September 2021), 1995. 

Panosetti, D., Böing, S., Schlemmer, L., and Schmidli, J.: Idealized Large-Eddy and Convection-Resolving Simulations of Moist Convection over Mountainous Terrain, J. Atmos. Sci., 73, 4021–4041,, 2016. 

Park, S.: A Unified Convection Scheme (UNICON). Part I: Formulation, J. Atmos. Sci., 71, 3902–3930,, 2014a. 

Park, S.: A Unified Convection Scheme (UNICON). Part II: Simulation, J. Atmos. Sci., 71, 3931–3973,, 2014b. 

Park, S. and Bretherton, C. S.: The University of Washington Shallow Convection and Moist Turbulence Schemes and Their Impact on Climate Simulations with the Community Atmosphere Model, J. Climate, 22, 3449–3469,, 2009. 

Park, S., Baek, E.-H., Kim, B.-M., and Kim, S.-J.: Impact of detrained cumulus on climate simulated by the Community Atmosphere Model Version 5 with a unified convection scheme, J. Adv. Model. Earth Sy., 9, 1399–1411,, 2017. 

Patz, J. A., Campbell-Lendrum, D., Holloway, T., and Foley, J. A.: Impact of regional climate change on human health, Nature, 438, 310–317,, 2005. 

Peng, J., Li, Z., Zhang, H., Liu, J., and Cribb, M.: Systematic Changes in Cloud Radiative Forcing with Aerosol Loading for Deep Clouds in the Tropics, J. Atmos. Sci., 73, 231–249,, 2016. 

Peng, M. S., Ridout, J. A., and Hogan, T. F.: Recent Modifications of the Emanuel Convective Scheme in the Navy Operational Global Atmospheric Prediction System, Mon. Weather Rev., 132, 1254–1268,<1254:RMOTEC>2.0.CO;2, 2004. 

Pergaud, J., Masson, V., Malardel, S., and Couvreux, F.: A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction, Bound.-Lay. Meteorol., 132, 83,, 2009. 

Perraud, E., Couvreux, F., Malardel, S., Lac, C., Masson, V., and Thouron, O.: Evaluation of Statistical Distributions for the Parametrization of Subgrid Boundary-Layer Clouds, Bound.-Lay. Meteorol., 140, 263–294,, 2011. 

Peters, J. M.: The impact of effective buoyancy and dynamic pressure forcing on vertical velocities within two-dimensional updrafts, J. Atmos. Sci., t3, 4531–4551,, 2016. 

Peters, K., Jakob, C., Davies, L., Khouider, B., and Majda, A. J.: Stochastic Behavior of Tropical Convection in Observations and a Multicloud Model, J. Atmos. Sci., 70, 3556–3575,, 2013. 

Peters, K., Crueger, T., Jakob, C., and Möbis, B.: Improved MJO-simulation in ECHAM6.3 by coupling a Stochastic Multicloud Model to the convection scheme, J. Adv. Model. Earth Sy., 9, 193–219,, 2017. 

Petersen, A. C., Beets, C., Dop, H. van, Duynkerke, P. G., and Siebesma, A. P.: Mass-Flux Characteristics of Reactive Scalars in the Convective Boundary Layer, J. Atmos. Sci., 56, 37–56,<0037:MFCORS>2.0.CO;2, 1999. 

Pezzi, L. P., Cavalcanti, I. F. A., and Mendonça, A. M.: A sensitivity study using two different convection schemes over south america, Revista Brasileira de Meteorologia, 23, 170–189,, 2008. 

Pham-Duc, B., Sylvestre, F., Papa, F., Frappart, F., Bouchez, C., and Crétaux, J.-F.: The Lake Chad hydrology under current climate change, Sci. Rep.-UK, 10, 5498,, 2020. 

Piriou, J.-M., Redelsperger, J.-L., Geleyn, J.-F., Lafore, J.-P., and Guichard, F.: An Approach for Convective Parameterization with Memory: Separating Microphysics and Transport in Grid-Scale Equations, J. Atmos. Sci., 64, 4127–4139,, 2007. 

Plant, R. S.: A review of the theoretical basis for bulk mass flux convective parameterization, Atmos. Chem. Phys., 10, 3529–3544,, 2010. 

Plant, R. S. and Craig, G. C.: A Stochastic Parameterization for Deep Convection Based on Equilibrium Statistics, J. Atmos. Sci., 65, 87–105,, 2008. 

Plant, R. S. and Yano, J.-I.: Parameterization of Atmospheric Convection: (In 2 Volumes) Volume 1: Theoretical Background and FormulationVolume 2: Current Issues and New Theories, Imperial College Press,, 2015. 

Prein, A. F., Gobiet, A., Suklitsch, M., Truhetz, H., Awan, N. K., Keuler, K., and Georgievski, G.: Added value of convection permitting seasonal simulations, Clim. Dynam., 41, 2655–2677,, 2013. 

Prein, A. F., Langhans, W., Fosser, G., Ferrone, A., Ban, N., Goergen, K., Keller, M., Tölle, M., Gutjahr, O., Feser, F., Brisson, E., Kollet, S., Schmidli, J., Lipzig, N. P. M. van, and Leung, R.: A review on regional convection-permitting climate modeling: Demonstrations, prospects, and challenges, Rev. Geophys., 53, 323–361,, 2015. 

Qian, L., Young, G. S., and Frank, W. M.: A Convective Wake Parameterization Scheme for Use in General Circulation Models, Mon. Weather Rev., 126, 456–469,<0456:ACWPSF>2.0.CO;2, 1998. 

Qin, Y., Lin, Y., Xu, S., Ma, H.-Y., and Xie, S.: A Diagnostic PDF Cloud Scheme to Improve Subtropical Low Clouds in NCAR Community Atmosphere Model (CAM5), J. Adv. Model. Earth Sy., 10, 320–341,, 2018. 

Raga, G. B., Jensen, J. B., and Baker, M. B.: Characteristics of Cumulus Band Clouds off the Coast of Hawaii, J. Atmos. Sci., 47, 338–356,<0338:COCBCO>2.0.CO;2, 1990. 

Raju, P. V. S., Bhatla, R., Almazroui, M., and Assiri, M.: Performance of convection schemes on the simulation of summer monsoon features over the South Asia CORDEX domain using RegCM-4.3, Int. J. Climatol., 35, 4695–4706,, 2015. 

Ramanathan, V. and Collins, W.: Thermodynamic regulation of ocean warming by cirrus clouds deduced from observations of the 1987 El Niño, Nature, 351, 27–32,, 1991. 

Randall, D., Khairoutdinov, M., Arakawa, A., and Grabowski, W.: Breaking the Cloud Parameterization Deadlock, B. Am. Meteorol. Soc., 84, 1547–1564,, 2003. 

Randall, D. A. and Pan, D.-M.: Implementation of the Arakawa-Schubert Cumulus Parameterization with a Prognostic Closure, in: The Representation of Cumulus Convection in Numerical Models, edited by: Emanuel, K. A. and Raymond, D. J., American Meteorological Society, Boston, MA, 137–144,, 1993. 

Randall, D. A., Shao, Q., and Moeng, C.-H.: A Second-Order Bulk Boundary-Layer Model, J. Atmos. Sci., 49, 1903–1923,<1903:ASOBBL>2.0.CO;2, 1992. 

Randall, D. A., Srinivasan, J., Nanjundiah, R. A., and Mukhopadhyay, P. (Eds.): Current Trends in the Representation of Physical Processes in Weather and Climate Models, Springer Singapore,, 2019. 

Rasp, S., Pritchard, M. S., and Gentine, P.: Deep learning to represent subgrid processes in climate models, P. Natl. Acad. Sci. USA, 115, 9684–9689,, 2018. 

Rauber, R. M., Stevens, B., Ochs, H. T., Knight, C., Albrecht, B. A., Blythe, A. M., Fairall, C. W., Jensen, J. B., Lasher-Trapp, S. G., Mayol-Bracero, O. L., Vali, G., Anderson, J. R., Baker, B. A., Bandy, A. R., Brunet, E., Brenguier, J. L., Brewer, W. A., Brown, P. R. A., Chuang, P., Cotton, W. R., Girolamo, L. D., Geerts, B., Gerber, H., Göke, S., Gomes, L., Heikes, B. G., Hudson, J. G., Kollias, P., Lawson, R. P., Krueger, S. K., Lenschow, D. H., Nuijens, L., O'Sullivan, D. W., Rilling, R. A., Rogers, D. C., Siebesma, A. P., Snodgrass, F., Stith, J. L., Thornton, D. C., Tucker, S., Twohy, C. H., and Zuidema, P.: Rain in shallow cumulus over the ocean: The RICO campaign, B. Am. Meteorol. Soc., 88, 1912–1928,, 2007. 

Raymond, D. J.: Regulation of Moist Convection over the West Pacific Warm Pool, J. Atmos. Sci., 52, 3945–3959,<3945:ROMCOT>2.0.CO;2, 1995. 

Raymond, D. J. and Blyth, A. M.: A Stochastic Mixing Model for Nonprecipitating Cumulus Clouds, J. Atmos. Sci., 43, 2708–2718,<2708:ASMMFN>2.0.CO;2, 1986. 

Raymond, D. J. and Emanuel, K. A.: The Kuo Cumulus Parameterization, in: The Representation of Cumulus Convection in Numerical Models, edited by: Emanuel, K. A. and Raymond, D. J., American Meteorological Society, Boston, MA, 145–147,, 1993. 

Rennó, N. O., Emanuel, K. A., and Stone, P. H.: Radiative-convective model with an explicit hydrologic cycle: 1. Formulation and sensitivity to model parameters, J. Geophys. Res.-Atmos., 99, 14429–14441,, 1994. 

Reuter, G. W. and Yau, M. K.: Mixing Mechanisms in Cumulus Congestus Clouds. Part II: Numerical Simulations, J. Atmos. Sci., 44, 798–827,<0798:MMICCC>2.0.CO;2, 1987. 

Rio, C. and Hourdin, F.: A Thermal Plume Model for the Convective Boundary Layer: Representation of Cumulus Clouds, J. Atmos. Sci., 65, 407–425,, 2008. 

Rio, C., Hourdin, F., Grandpeix, J.-Y., and Lafore, J.-P.: Shifting the diurnal cycle of parameterized deep convection over land, Geophys. Res. Lett., 36, L07809,, 2009. 

Rio, C., Hourdin, F., Couvreux, F., and Jam, A.: Resolved Versus Parametrized Boundary-Layer Plumes. Part II: Continuous Formulations of Mixing Rates for Mass-Flux Schemes, Bound.-Lay. Meteorol., 135, 469–483,, 2010. 

Rio, C., Grandpeix, J.-Y., Hourdin, F., Guichard, F., Couvreux, F., Lafore, J.-P., Fridlind, A., Mrowiec, A., Roehrig, R., Rochetin, N., Lefebvre, M.-P., and Idelkadi, A.: Control of deep convection by sub-cloud lifting processes: the ALP closure in the LMDZ5B general circulation model, Clim. Dynam., 40, 2271–2292,, 2013. 

Rio, C., Del Genio, A. D., and Hourdin, F.: Ongoing Breakthroughs in Convective Parameterization, Curr. Clim. Change Rep., 5, 95–111,, 2019. 

Rocha, R. P. D. and Caetano, E.: The role of convective parameterization in the simulation of a cyclone over the South Atlantic, Atmosfera. 23, 1–23, 2010. 

Rochetin, N., Couvreux, F., Grandpeix, J.-Y., and Rio, C.: Deep Convection Triggering by Boundary Layer Thermals. Part I: LES Analysis and Stochastic Triggering Formulation, J. Atmos. Sci., 71, 496–514,, 2014a. 

Rochetin, N., Grandpeix, J.-Y., Rio, C., and Couvreux, F.: Deep Convection Triggering by Boundary Layer Thermals. Part II: Stochastic Triggering Parameterization for the LMDZ GCM, J. Atmos. Sci., 71, 515–538,, 2014b. 

Romps, D. M.: A Direct Measure of Entrainment, J. Atmos. Sci., 67, 1908–1927,, 2010. 

Romps, D. M.: The Stochastic Parcel Model: A deterministic parameterization of stochastically entraining convection, J. Adv. Model. Earth Sy., 8, 319–344,, 2016. 

Romps, D. M. and Kuang, Z.: Do Undiluted Convective Plumes Exist in the Upper Tropical Troposphere?, J. Atmos. Sci., 67, 468–484,, 2010a. 

Romps, D. M. and Kuang, Z.: Nature versus Nurture in Shallow Convection, J. Atmos. Sci., 67, 1655–1666,, 2010b. 

Rosa, D. and Collins, W. D.: A case study of subdaily simulated and observed continental convective precipitation: CMIP5 and multiscale global climate models comparison, Geophys. Res. Lett., 40, 5999–6003,, 2013. 

Rosenfeld, D., Lohmann, U., Raga, G., O'Dowd, C., Kulmala, M., Sandro, F., Reissell, A., and Andreae, M.: Flood or drought: How do aerosols affect precipitation?, Science, 321, 1309–1313, 2008. 

Rougier, J., Sexton, D. M. H., Murphy, J. M., and Stainforth, D.: Analyzing the Climate Sensitivity of the HadSM3 Climate Model Using Ensembles from Different but Related Experiments, J. Climate, 22, 3540–3557,, 2009. 

Ruiz, J. J., Pulido, M., and Miyoshi, T.: Estimating Model Parameters with Ensemble-Based Data Assimilation: A Review, J. Meteorol. Soc. Jpn., 91, 79–99,, 2013. 

Sakradzija, M. and Klocke, D.: Physically Constrained Stochastic Shallow Convection in Realistic Kilometer-Scale Simulations, J. Adv. Model. Earth Sy., 10, 2755–2776,, 2018. 

Sakradzija, M., Seifert, A., and Heus, T.: Fluctuations in a quasi-stationary shallow cumulus cloud ensemble, Nonlin. Processes Geophys., 22, 65–85,, 2015. 

Sakradzija, M., Seifert, A., and Dipankar, A.: A stochastic scale-aware parameterization of shallow cumulus convection across the convective gray zone, J. Adv. Model. Earth Sy., 8, 786–812,, 2016. 

Sanderson, B. M., Piani, C., Ingram, W. J., Stone, D. A., and Allen, M. R.: Towards constraining climate sensitivity by linear analysis of feedback patterns in thousands of perturbed-physics GCM simulations, Clim. Dynam., 30, 175–190,, 2008. 

Sato, T., Miura, H., Satoh, M., Takayabu, Y. N., and Wang, Y.: Diurnal Cycle of Precipitation in the Tropics Simulated in a Global Cloud-Resolving Model, J. Climate, 22, 4809–4826,, 2009. 

Schlemmer, L. and Hohenegger, C.: The Formation of Wider and Deeper Clouds as a Result of Cold-Pool Dynamics, J. Atmos. Sci., 71, 2842–2858,, 2014. 

Schmidt, G. A., Bader, D., Donner, L. J., Elsaesser, G. S., Golaz, J.-C., Hannay, C., Molod, A., Neale, R. B., and Saha, S.: Practice and philosophy of climate model tuning across six US modeling centers, Geosci. Model Dev., 10, 3207–3223,, 2017. 

Schneider, T., Lan, S., Stuart, A., and Teixeira, J.: Earth System Modeling 2.0: A Blueprint for Models That Learn From Observations and Targeted High-Resolution Simulations, Geophys. Res. Lett., 44, 12396–12417,, 2017. 

Shin, J. and Park, S.: A Stochastic Unified Convection Scheme (UNICON). Part I: Formulation and Single-Column Simulation for Shallow Convection, J. Atmos. Sci., 77, 583–610,, 2020. 

Shutts, G.: A kinetic energy backscatter algorithm for use in ensemble prediction systems, Q. J. Roy. Meteor. Soc., 131, 3079–3102,, 2005. 

Siebesma, A. P.: Shallow Cumulus Convection, in: Buoyant Convection in Geophysical Flows, edited by: Plate, E. J., Fedorovich, E. E., Viegas, D. X., and Wyngaard, J. C., Springer Netherlands, Dordrecht, 441–486,, 1998. 

Siebesma, A. P. and Cuijpers, J. W. M.: Evaluation of Parametric Assumptions for Shallow Cumulus Convection, J. Atmos. Sci., 52, 650–666,<0650:EOPAFS>2.0.CO;2, 1995. 

Siebesma, A. P. and Holtslag, A. A. M.: Model Impacts of Entrainment and Detrainment Rates in Shallow Cumulus Convection, J. Atmos. Sci., 53, 2354–2364,<2354:MIOEAD>2.0.CO;2, 1996. 

Siebesma, A. P. and Teixeira, J.: An Advection-Diffusion scheme for the convective boundary layer: description and 1d-results, 14th Symp. on Boundary Layers and Turbulence, 9 August 2000, Aspen, CO, 133–136, 2000. 

Siebesma, A. P., Bretherton, C. S., Brown, A., Chlond, A., Cuxart, J., Duynkerke, P. G., Jiang, H., Khairoutdinov, M., Lewellen, D., Moeng, C.-H., Sanchez, E., Stevens, B., and Stevens, D. E.: A Large Eddy Simulation Intercomparison Study of Shallow Cumulus Convection, J. Atmos. Sci., 60, 1201–1219,<1201:ALESIS>2.0.CO;2, 2003. 

Siebesma, A. P., Soares, P. M. M., and Teixeira, J.: A Combined Eddy-Diffusivity Mass-Flux Approach for the Convective Boundary Layer, J. Atmos. Sci., 64, 1230–1248,, 2007. 

Simpson, J.: On Cumulus Entrainment and One-Dimensional Models, J. Atmos. Sci., 28, 449–455,<0449:OCEAOD>2.0.CO;2, 1971. 

Simpson, J. and Wiggert, V.: Modes of Precipitating Cumulus Towers, Mon. Weather Rev., 97, 471–489,<0471:MOPCT>2.3.CO;2, 1969. 

Singh, M. S., Warren, R. A., and Jakob, C.: A Steady-State Model for the Relationship Between Humidity, Instability, and Precipitation in the Tropics, J. Adv. Model. Earth Sy.,11, 3973–3994,, 2019. 

Skofronick-Jackson, G., Kulie, M., Milani, L., Munchak, S. J., Wood, N. B., and Levizzani, V.: Satellite Estimation of Falling Snow: A Global Precipitation Measurement (GPM) Core Observatory Perspective, J. Appl. Meteorol. Clim., 58, 1429–1448,, 2019. 

Slingo, J., Blackburn, M., Betts, A., Brugge, R., Hodges, K., Hoskins, B., Miller, M., Steenman-Clark, L., and Thuburn, J.: Mean climate and transience in the tropics of the UGAMP GCM: Sensitivity to convective parametrization, Q. J. Roy. Meteor. Soc., 120, 881–922,, 1994. 

Smagorinsky, J.: On the inclusion of moist adiabatic processes in numerical prediction models, Ber. Dtsch. Wetterdienstes, 38, 82–90, 1956. 

Smith, L. A.: What might we learn from climate forecasts?, P. Natl. Acad. Sci. USA, 99, 2487–2492,, 2002. 

Smith, R. N. B.: A scheme for predicting layer clouds and their water content in a general circulation model, Q. J. Roy. Meteor. Soc., 116, 435–460,, 1990. 

Soares, P. M. M., Miranda, P. M. A., Siebesma, A. P., and Teixeira, J.: An eddy-diffusivity/mass-flux parametrization for dry and shallow cumulus convection, Q. J. Roy. Meteor. Soc., 130, 3365–3383,, 2004. 

Sommeria, G. and Deardorff, J. W.: Subgrid-Scale Condensation in Models of Nonprecipitating Clouds, J. Atmos. Sci., 34, 344–355,<0344:SSCIMO>2.0.CO;2, 1977. 

Song, F. and Zhang, G. J.: Improving Trigger Functions for Convective Parameterization Schemes Using GOAmazon Observations, J. Climate, 30, 8711–8726,, 2017. 

Song, H., Lin, W., Lin, Y., Wolf, A. B., Neggers, R., Donner, L. J., Genio, A. D. D., and Liu, Y.: Evaluation of Precipitation Simulated by Seven SCMs against the ARM Observations at the SGP Site, J. Climate, 26, 5467–5492,, 2013. 

Song, X. and Zhang, G. J.: Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part I: Climatology and Atmospheric Feedback, J. Climate, 22, 4299–4315,, 2009. 

Song, X. and Zhang, G. J.: Microphysics parameterization for convective clouds in a global climate model: Description and single-column model tests, J. Geophys. Res.-Atmos., 116, D02201,, 2011. 

Song, X. and Zhang, G. J.: The Roles of Convection Parameterization in the Formation of Double ITCZ Syndrome in the NCAR CESM: I. Atmospheric Processes, J. Adv. Model. Earth Sy., 10, 842–866,, 2018. 

Song, X., Zhang, G. J., and Li, J.-L. F.: Evaluation of Microphysics Parameterization for Convective Clouds in the NCAR Community Atmosphere Model CAM5, J. Climate, 25, 8568–8590,, 2012. 

Song, Y., Wikle, C. K., Anderson, C. J., and Lack, S. A.: Bayesian Estimation of Stochastic Parameterizations in a Numerical Weather Forecasting Model, Mon. Weather Rev., 135, 4045–4059,, 2007. 

Squires, P.: Penetrative Downdraughts in Cumuli, Tellus, 10, 381–389,, 1958. 

Squires, P. and Turner, J. S.: An entraining jet model for cumulo-nimbus updraughts, Tellus, 14, 422–434,, 1962. 

Stechmann, S. N. and Neelin, J. D.: A Stochastic Model for the Transition to Strong Convection, J. Atmos. Sci., 68, 2955–2970,, 2011. 

Stensrud, D. J.: Parameterization Schemes: Keys to Understanding Numerical Weather Prediction Models, Cambridge University Press, Cambridge,, 2007. 

Stephens, G. L., L'Ecuyer, T., Forbes, R., Gettelmen, A., Golaz, J.-C., Bodas-Salcedo, A., Suzuki, K., Gabriel, P., and Haynes, J.: Dreary state of precipitation in global models, J. Geophys. Res.-Atmos., 115, D24211,, 2010. 

Stephens, G. L., van den Heever, S. C., Haddad, Z. S., Posselt, D. J., Storer, R. L., Grant, L. D., Sy, O. O., Rao, T. N., Tanelli, S., and Peral, E.: A distributed small satellite approach for measuring convective transports in the Earth's atmosphere, IEEE T. Geosci. Remote, 58, 4–13,, 2020. 

Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric component of the MPI-M Earth System Model: ECHAM6, J. Adv. Model. Earth Sy., 5, 146–172,, 2013. 

Stirling, A. J. and Stratton, R. A.: Entrainment processes in the diurnal cycle of deep convection over land, Q. J. Roy. Meteor. Soc., 138, 1135–1149,, 2012. 

Stommel, H.: ENTRAINMENT OF AIR INTO A CUMULUS CLOUD: (Paper presented 27 December 1946 at the Annual Meeting, A.M.S., Cambridge, Massachusetts), J. Atmos. Sci., 4, 91–94,<0091:EOAIAC>2.0.CO;2, 1947. 

Storer, R. L., van den Heever, S. C., and Stephens, G. L.: Modeling Aerosol Impacts on Convective Storms in Different Environments, J. Atmos. Sci., 67, 3904–3915,, 2010. 

Storer, R. L., Zhang, G. J., and Song, X.: Effects of Convective Microphysics Parameterization on Large-Scale Cloud Hydrological Cycle and Radiative Budget in Tropical and Midlatitude Convective Regions, J. Climate, 28, 9277–9297,, 2015. 

Stratton, R. A. and Stirling, A. J.: Improving the diurnal cycle of convection in GCMs, Q. J. Roy. Meteor. Soc., 138, 1121–1134,, 2012. 

Sud, Y. C. and Walker, G. K.: Microphysics of Clouds with the Relaxed Arakawa–Schubert Scheme (McRAS). Part I: Design and Evaluation with GATE Phase III Data, J. Atmos. Sci., 56, 3196–3220,<3196:MOCWTR>2.0.CO;2, 1999. 

Suhas, E. and Zhang, G. J.: Evaluation of Trigger Functions for Convective Parameterization Schemes Using Observations, J. Climate, 27, 7647–7666,, 2014. 

Sun, J. and Pritchard, M. S.: Effects of explicit convection on global land-atmosphere coupling in the superparameterized CAM, J. Adv. Model. Earth Sy., 8, 1248–1269,, 2016. 

Sun, Y., Solomon, S., Dai, A., and Portmann, R. W.: How Often Does It Rain?, J. Climate, 19, 916–934,, 2006. 

Sundqvist, H.: A parameterization scheme for non-convective condensation including prediction of cloud water content, Q. J. Roy. Meteor. Soc., 104, 677–690,, 1978. 

Sundqvist, H.: Parameterization of Condensation and Associated Clouds in Models for Weather Prediction and General Circulation Simulation, in: Physically-Based Modelling and Simulation of Climate and Climatic Change: Part 1, edited by: Schlesinger, M. E., Springer Netherlands, Dordrecht, 433–461,, 1988. 

Sušelj, K., Teixeira, J., and Matheou, G.: Eddy Diffusivity/Mass Flux and Shallow Cumulus Boundary Layer: An Updraft PDF Multiple Mass Flux Scheme, J. Atmos. Sci., 69, 1513–1533,, 2012. 

Sušelj, K., Teixeira, J., and Chung, D.: A Unified Model for Moist Convective Boundary Layers Based on a Stochastic Eddy-Diffusivity/Mass-Flux Parameterization, J. Atmos. Sci., 70, 1929–1953,, 2013. 

Sušelj, K., Hogan, T. F., and Teixeira, J.: Implementation of a Stochastic Eddy-Diffusivity/Mass-Flux Parameterization into the Navy Global Environmental Model, Weather Forecast., 29, 1374–1390,, 2014. 

Suselj, K., Kurowski, M. J., and Teixeira, J.: A Unified Eddy-Diffusivity/Mass-Flux Approach for Modeling Atmospheric Convection, J. Atmos. Sci., 76, 2505–2537,, 2019a. 

Suselj, K., Kurowski, M. J., and Teixeira, J.: On the Factors Controlling the Development of Shallow Convection in Eddy-Diffusivity/Mass-Flux Models, J. Atmos. Sci., 76, 433–456,, 2019b. 

Tan, Z., Kaul, C. M., Pressel, K. G., Cohen, Y., Schneider, T., and Teixeira, J.: An Extended Eddy-Diffusivity Mass-Flux Scheme for Unified Representation of Subgrid-Scale Turbulence and Convection, J. Adv. Model. Earth Sy., 10, 770–800,, 2018. 

Tao, W.-K., Chen, J.-P., Li, Z., Wang, C., and Zhang, C.: Impact of aerosols on convective clouds and precipitation, Rev. Geophys., 50, RG2001,, 2012. 

Tapiador, F. J., Hou, A. Y., de Castro, M., Checa, R., Cuartero, F., and Barros, A. P.: Precipitation estimates for hydroelectricity, Energ. Environ. Sci., 4, 4435–4448,, 2011. 

Tapiador, F. J., Turk, F. J., Petersen, W., Hou, A. Y., García-Ortega, E., Machado, L. A. T., Angelis, C. F., Salio, P., Kidd, C., Huffman, G. J., and de Castro, M.: Global precipitation measurement: Methods, datasets and applications, Atmos. Res., 104–105, 70–97,, 2012. 

Tapiador, F. J., Navarro, A., Levizzani, V., García-Ortega, E., Huffman, G. J., Kidd, C., Kucera, P. A., Kummerow, C. D., Masunaga, H., Petersen, W. A., Roca, R., Sánchez, J.-L., Tao, W.-K., and Turk, F. J.: Global precipitation measurements for validating climate models, Atmos. Res., 197, 1–20,, 2017. 

Tapiador, F. J., Navarro, A., Jiménez, A., Moreno, R., and García-Ortega, E.: Discrepancies with satellite observations in the spatial structure of global precipitation as derived from global climate models, . J. Roy. Meteor. Soc., 144, 419–435,, 2018. 

Tapiador, F. J., Roca, R., Del Genio, A., Dewitte, B., Petersen, W., and Zhang, F.: Is Precipitation a Good Metric for Model Performance?, B. Am. Meteorol. Soc., 100, 223–233,, 2019a. 

Tapiador, F. J., Sánchez, J.-L., and García-Ortega, E.: Empirical values and assumptions in the microphysics of numerical models, Atmos. Res., 215, 214–238,, 2019b. 

Tawfik, A. B. and Dirmeyer, P. A.: A process-based framework for quantifying the atmospheric preconditioning of surface-triggered convection, Geophys. Res. Lett., 41, 173–178,, 2014. 

Tawfik, A. B., Lawrence, D. M., and Dirmeyer, P. A.: Representing subgrid convective initiation in the Community Earth System Model, J. Adv. Model. Earth Sy., 9, 1740–1758,, 2017. 

Taylor, G. R. and Baker, M. B.: Entrainment and Detrainment in Cumulus Clouds, J. Atmos. Sci., 48, 112–121,<0112:EADICC>2.0.CO;2, 1991. 

Teixeira, J. and Kim, Y. J.: On a simple parameterization of convective cloud fraction, Asia-Pac. J. Atmos. Sci., 44, 191–199, 2008. 

Teixeira, J. and Reynolds, C. A.: Stochastic Nature of Physical Parameterizations in Ensemble Prediction: A Stochastic Convection Approach, Mon. Weather Rev., 136, 483–496,, 2008. 

Telford, J. W.: Turbulence, entrainment, and mixing in cloud dynamics, Pure Appl. Geophys., 113, 1067–1084,, 1975. 

Thayer-Calder, K.: Downdraft impacts on tropical convection, Colorado State University, Publication Number: AAT 3565466, ISBN 9781303152504, 2012. 

Thayer-Calder, K. and Randall, D. A.: The Role of Convective Moistening in the Madden–Julian Oscillation, J. Atmos. Sci., 66, 3297–3312,, 2009. 

Thayer-Calder, K., Gettelman, A., Craig, C., Goldhaber, S., Bogenschutz, P. A., Chen, C.-C., Morrison, H., Höft, J., Raut, E., Griffin, B. M., Weber, J. K., Larson, V. E., Wyant, M. C., Wang, M., Guo, Z., and Ghan, S. J.: A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model, Geosci. Model Dev., 8, 3801–3821,, 2015. 

Tiedtke, M.: A Comprehensive Mass Flux Scheme for Cumulus Parameterization in Large-Scale Models, Mon. Weather Rev., 117, 1779–1800,<1779:ACMFSF>2.0.CO;2, 1989. 

Tiedtke, M.: Representation of Clouds in Large-Scale Models, Mon. Weather Rev., 121, 3040–3061,<3040:ROCILS>2.0.CO;2, 1993. 

Tokioka, T., Yamazaki, K., Kitoh, A., and Ose, T.: The Equatorial 30-60 day Oscillation and the Arakawa-Schubert Penetrative Cumulus Parameterization, J. Meteorol. Soc. Jpn., 66, 883–901,, 1988. 

Tompkins, A., Bechtold, P., Beljaars, A., Benedetti, A., Cheinet, S., Janiskova, M., Köhler, M., Lopez, P., and Morcrette, J.-J.: Moist physical processes in the IFS: Progress and Plans, Technical memorandum,, 2004. 

Tompkins, A. M.: A Prognostic Parameterization for the Subgrid-Scale Variability of Water Vapor and Clouds in Large-Scale Models and Its Use to Diagnose Cloud Cover, J. Atmos. Sci., 59, 1917–1942,<1917:APPFTS>2.0.CO;2, 2002. 

Tompkins, A. M. and Berner, J.: A stochastic convective approach to account for model uncertainty due to unresolved humidity variability, J. Geophys. Res.-Atmos., 113, D18101,, 2008. 

Trenberth, K. E.: Changes in precipitation with climate change, Clim. Res., 47, 123–138, 2011. 

Troen, I. B. and Mahrt, L.: A simple model of the atmospheric boundary layer; sensitivity to surface evaporation, Bound.-Lay. Meteorol., 37, 129–148,, 1986. 

Turner, J. S.: The “starting plume” in neutral surroundings, J. Fluid Mech., 13, 356–368,, 1962. 

Ushio, T. and Kachi, M.: Kalman Filtering Applications for Global Satellite Mapping of Precipitation (GSMaP), in: Satellite Rainfall Applications for Surface Hydrology, edited by: Gebremichael, M. and Hossain, F., Springer Netherlands, Dordrecht, 105–123,, 2010. 

Vaidya, S. S. and Singh, S. S.: Thermodynamic Adjustment Parameters in the Betts–Miller Scheme of Convection, Weather Forecast., 12, 819–825,<0819:TAPITB>2.0.CO;2, 1997. 

Vaidya, S. S. and Singh, S. S.: Applying the Betts–Miller–Janjić Scheme of Convection in Prediction of the Indian Monsoon, Weather Forecast., 15, 349–356,<0349:ATBMJS>2.0.CO;2, 2000. 

van den Heever, S. C. and Cotton, W. R.: Urban Aerosol Impacts on Downwind Convective Storms, J. Appl. Meteorol. Clim., 46, 828–850,, 2007. 

van den Heever, S. C., Stephens, G. L., and Wood, N. B.: Aerosol Indirect Effects on Tropical Convection Characteristics under Conditions of Radiative–Convective Equilibrium, J. Atmos. Sci., 68, 699–718,, 2011. 

van Laar, T. W.: Spatial patterns in shallow cumulus cloud populations over a heterogeneous surface, text.thesis.doctoral, Universität zu Köln, (last access: 19 September 2021), 2019. 

Vogelmann, A. M., McFarquhar, G. M., Ogren, J. A., Turner, D. D., Comstock, J. M., Feingold, G., Long, C. N., Jonsson, H. H., Bucholtz, A., Collins, D. R., Diskin, G. S., Gerber, H., Lawson, R. P., Woods, R. K., Andrews, E., Yang, H.-J., Chiu, J. C., Hartsock, D., Hubbe, J. M., Lo, C., Marshak, A., Monroe, J. W., McFarlane, S. A., Schmid, B., Tomlinson, J. M., and Toto, T.: RACORO Extended-Term Aircraft Observations of Boundary Layer Clouds, B. Am. Meteorol. Soc., 93, 861–878,, 2012. 

Volterra, V.: Variazioni e fluttuazioni del numero d'individui in specie animali conviventi, Memoria della Reale Accademia Nazionale dei Lincei, 2, 209, 1926. 

von Salzen, K. and McFarlane, N. A.: Parameterization of the Bulk Effects of Lateral and Cloud-Top Entrainment in Transient Shallow Cumulus Clouds, Atmos. Sci., 59, 1405–1430,<1405:POTBEO>2.0.CO;2, 2002. 

Wagner, A., Heinzeller, D., Wagner, S., Rummler, T., and Kunstmann, H.: Explicit Convection and Scale-Aware Cumulus Parameterizations: High-Resolution Simulations over Areas of Different Topography in Germany, Mon. Weather Rev., 146, 1925–1944,, 2018. 

Wagner, T. J., Turner, D. D., Berg, L. K., and Krueger, S. K.: Ground-Based Remote Retrievals of Cumulus Entrainment Rates, J. Atmos. Ocean Tech., 30, 1460–1471,, 2013. 

Wagner, T. M. and Graf, H.-F.: An Ensemble Cumulus Convection Parameterization with Explicit Cloud Treatment, J. Atmos. Sci., 67, 3854–3869,, 2010. 

Walters, D., Baran, A. J., Boutle, I., Brooks, M., Earnshaw, P., Edwards, J., Furtado, K., Hill, P., Lock, A., Manners, J., Morcrette, C., Mulcahy, J., Sanchez, C., Smith, C., Stratton, R., Tennant, W., Tomassini, L., Van Weverberg, K., Vosper, S., Willett, M., Browse, J., Bushell, A., Carslaw, K., Dalvi, M., Essery, R., Gedney, N., Hardiman, S., Johnson, B., Johnson, C., Jones, A., Jones, C., Mann, G., Milton, S., Rumbold, H., Sellar, A., Ujiie, M., Whitall, M., Williams, K., and Zerroukat, M.: The Met Office Unified Model Global Atmosphere 7.0/7.1 and JULES Global Land 7.0 configurations, Geosci. Model Dev., 12, 1909–1963,, 2019. 

Wang, W. and Schlesinger, M. E.: The Dependence on Convection Parameterization of the Tropical Intraseasonal Oscillation Simulated by the UIUC 11-Layer Atmospheric GCM, J. Climate, 12, 1423–1457,<1423:TDOCPO>2.0.CO;2, 1999. 

Wang, X. and Zhang, M.: An analysis of parameterization interactions and sensitivity of single-column model simulations to convection schemes in CAM4 and CAM5, J. Geophys. Res.-Atmos., 118, 8869–8880,, 2013. 

Wang, X. and Zhang, M.: Vertical velocity in shallow convection for different plume types, J. Adv. Model. Earth Sy., 6, 478–489,, 2014. 

Wang, Y., Zhou, L., and Hamilton, K.: Effect of Convective Entrainment/Detrainment on the Simulation of the Tropical Precipitation Diurnal Cycle, Mon. Weather Rev., 135, 567–585,, 2007. 

Wang, Y., Zhang, G. J., and Craig, G. C.: Stochastic convective parameterization improving the simulation of tropical precipitation variability in the NCAR CAM5, Geophys. Res. Lett., 43, 6612–6619,, 2016. 

Warner, J.: The Microstructure of Cumulus Cloud. Part III. The Nature of the Updraft, J. Atmos. Sci., 27, 682–688,<0682:TMOCCP>2.0.CO;2, 1970. 

Watanabe, M., Emori, S., Satoh, M., and Miura, H.: A PDF-based hybrid prognostic cloud scheme for general circulation models, Clim. Dynam., 33, 795–816,, 2009. 

Watanabe, M., Suzuki, T., O'ishi, R., Komuro, Y., Watanabe, S., Emori, S., Takemura, T., Chikira, M., Ogura, T., Sekiguchi, M., Takata, K., Yamazaki, D., Yokohata, T., Nozawa, T., Hasumi, H., Tatebe, H., and Kimoto, M.: Improved Climate Simulation by MIROC5: Mean States, Variability, and Climate Sensitivity, J. Climate, 23, 6312–6335,, 2010. 

Watanabe, S., Hajima, T., Sudo, K., Nagashima, T., Takemura, T., Okajima, H., Nozawa, T., Kawase, H., Abe, M., Yokohata, T., Ise, T., Sato, H., Kato, E., Takata, K., Emori, S., and Kawamiya, M.: MIROC-ESM 2010: model description and basic results of CMIP5-20c3m experiments, Geosci. Model Dev., 4, 845–872,, 2011. 

Wilcox, E. M. and Donner, L. J.: The Frequency of Extreme Rain Events in Satellite Rain-Rate Estimates and an Atmospheric General Circulation Model, J. Climate, 20, 53–69,, 2007. 

Willet, M. R. and Whitall, M. A.: A simple prognostic based convective entrainment rate for the Unified Model: Description and tests, Met Office internal) Forecasting Research Technical Reports617, 2017. 

Witek, M. L., Teixeira, J., and Matheou, G.: An Integrated TKE-Based Eddy Diffusivity/Mass Flux Boundary Layer Closure for the Dry Convective Boundary Layer, J. Atmos. Sci., 68, 1526–1540,, 2011. 

Woetzel, J., Pinner, D., Samandari, H., Engel, H., Krishnan, M., Boland, B., and Powis, C.: Climate and risk response: Physical hazars and socioeconomic impacts, McKinsey Global Institute, 18, 164,, 2020. 

Wu, C.-M. and Arakawa, A.: A Unified Representation of Deep Moist Convection in Numerical Modeling of the Atmosphere. Part II, J. Atmos. Sci., 71, 2089–2103,, 2014. 

Wu, E., Yang, H., Kleissl, J., Suselj, K., Kurowski, M. J., and Teixeira, J.: On the Parameterization of Convective Downdrafts for Marine Stratocumulus Clouds, Mon. Weather Rev., 148, 1931–1950,, 2020. 

Wu, L., Wong, S., Wang, T., and Huffman, G. J.: Moist convection: a key to tropical wave–moisture interaction in Indian monsoon intraseasonal oscillation, Clim. Dynam., 51, 3673–3684,, 2018. 

Wu, T.: A mass-flux cumulus parameterization scheme for large-scale models: description and test with observations, Clim. Dynam., 38, 725–744,, 2012. 

Wu, X., Deng, L., Song, X., Vettoretti, G., Peltier, W. R., and Zhang, G. J.: Impact of a modified convective scheme on the Madden-Julian Oscillation and El Niño–Southern Oscillation in a coupled climate model, Geophys. Res. Lett., 34, L16823,, 2007. 

Wyant, M. C., Bretherton, C. S., Rand, H. A., and Stevens, D. E.: Numerical Simulations and a Conceptual Model of the Stratocumulus to Trade Cumulus Transition, J. Atmos. Sci., 54, 168–192,<0168:NSAACM>2.0.CO;2, 1997. 

Wyngaard, J. C.: Toward Numerical Modeling in the “Terra Incognita”, J. Atmos. Sci., 61, 1816–1826,<1816:TNMITT>2.0.CO;2, 2004. 

Xie, P., Joyce, R., Wu, S., Yoo, S.-H., Yarosh, Y., Sun, F., and Lin, R.: Reprocessed, Bias-Corrected CMORPH Global High-Resolution Precipitation Estimates from 1998, J. Hydrometeorol., 18, 1617–1641,, 2017. 

Xie, S. and Zhang, M.: Impact of the convection triggering function on single-column model simulations, J. Geophys. Res.-Atmos., 105, 14983–14996,, 2000. 

Xu, K.-M. and Randall, D. A.: A Semiempirical Cloudiness Parameterization for Use in Climate Models, J. Atmos. Sci., 53, 3084–3102,<3084:ASCPFU>2.0.CO;2, 1996. 

Xu, K.-M., Cederwall, R. T., Donner, L. J., Grabowski, W. W., Guichard, F., Johnson, D. E., Khairoutdinov, M., Krueger, S. K., Petch, J. C., Randall, D. A., Seman, C. J., Tao, W.-K., Wang, D., Xie, S. C., Yio, J. J., and Zhang, M.-H.: An intercomparison of cloud-resolving models with the atmospheric radiation measurement summer 1997 intensive observation period data, Q. J. Roy. Meteor. Soc., 128, 593–624,, 2002. 

Yanai, M., Esbensen, S., and Chu, J.-H.: Determination of Bulk Properties of Tropical Cloud Clusters from Large-Scale Heat and Moisture Budgets, J. Atmos. Sci., 30, 611–627,<0611:DOBPOT>2.0.CO;2, 1973. 

Yang, G.-Y. and Slingo, J.: The Diurnal Cycle in the Tropics, Mon. Weather Rev., 129, 784–801,<0784:TDCITT>2.0.CO;2, 2001. 

Yano, J., Bénard, P., Couvreux, F., and Lahellec, A.: NAM–SCA: A Nonhydrostatic Anelastic Model with Segmentally Constant Approximations, Mon. Weather Rev., 138, 1957–1974,, 2010. 

Yano, J. I.: Formulation structure of the mass-flux convection parameterization, Dynam. Atmos. Oceans, 67, 1–28,, 2014. 

Yano, J.-I. and Baizig, H.: Single SCA-plume dynamics, Dynam. Atmos. Oceans, 58, 62–94,, 2012. 

Yano, J.-I. and Plant, R.: Finite departure from convective quasi-equilibrium: periodic cycle and discharge–recharge mechanism, Q. J. Roy. Meteor. Soc., 138, 626–637,, 2012a. 

Yano, J.-I. and Plant, R. S.: Convective quasi-equilibrium, Rev. Geophys., 50, RG4004,, 2012b. 

Yano, J.-I., Bister, M., Fuchs, Ž., Gerard, L., Phillips, V. T. J., Barkidija, S., and Piriou, J.-M.: Phenomenology of convection-parameterization closure, Atmos. Chem. Phys., 13, 4111–4131,, 2013. 

Zhang, C., Wang, Y., and Hamilton, K.: Improved Representation of Boundary Layer Clouds over the Southeast Pacific in ARW-WRF Using a Modified Tiedtke Cumulus Parameterization Scheme, Mon. Weather Rev., 139, 3489–3513,, 2011. 

Zhang, D.-L. and Fritsch, J. M.: Numerical Simulation of the Meso-β Scale Structure and Evolution of the 1977 Johnstown Flood. Part I: Model Description and Verification, J. Atmos. Sci., 43, 1913–1944,<1913:NSOTMS>2.0.CO;2, 1986. 

Zhang, G. J.: Convective quasi-equilibrium in midlatitude continental environment and its effect on convective parameterization, J. Geophys. Res.-Atmos., 107, ACL 12-1–ACL 12-16,, 2002. 

Zhang, G. J.: Convective quasi-equilibrium in the tropical western Pacific: Comparison with midlatitude continental environment, J. Geophys. Res.-Atmos., 108,, 2003a. 

Zhang, G. J.: Roles of tropospheric and boundary layer forcing in the diurnal cycle of convection in the U.S. southern great plains, Geophys. Res. Lett., 30, 2281,, 2003b. 

Zhang, G. J.: Effects of entrainment on convective available potential energy and closure assumptions in convection parameterization, J. Geophys. Res.-Atmos., 114, D07109,, 2009. 

Zhang, G. J. and McFarlane, N. A.: Sensitivity of climate simulations to the parameterization of cumulus convection in the Canadian climate centre general circulation model, Atmosphere-Ocean, 33, 407–446,, 1995. 

Zhang, G. J. and Mu, M.: Effects of modifications to the Zhang-McFarlane convection parameterization on the simulation of the tropical precipitation in the National Center for Atmospheric Research Community Climate Model, version 3, J. Geophys. Res.-Atmos., 110, D09109,, 2005a. 

Zhang, G. J. and Mu, M.: Simulation of the Madden–Julian Oscillation in the NCAR CCM3 Using a Revised Zhang–McFarlane Convection Parameterization Scheme, J. Climate, 18, 4046–4064,, 2005b. 

Zhang, G. J. and Song, X.: Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part II: Coupled Feedback and the Role of Ocean Heat Transport, J. Climate, 23, 800–812,, 2010. 

Zhang, G. J. and Song, X.: Parameterization of Microphysical Processes in Convective Clouds in Global Climate Models, Meteor. Mon., 56, 12.1–12.18,, 2016. 

Zhang, G. J. and Wang, H.: Toward mitigating the double ITCZ problem in NCAR CCSM3, Geophys. Res. Lett., 33, L06709,, 2006. 

Zhang, J., Lohmann, U., and Stier, P.: A microphysical parameterization for convective clouds in the ECHAM5 climate model: Single-column model results evaluated at the Oklahoma Atmospheric Radiation Measurement Program site, J. Geophys. Res.-Atmos., 110, D15S07,, 2005. 

Zhang, Z., Tallapragada, V., Kieu, C., Trahan, S., and Wang, W.: HWRF Based Ensemble Prediction System Using Perturbations from GEFS and Stochastic Convective Trigger Function, Tropical Cyclone Research and Review, 3, 145–161,, 2014.  

Zhao, M.: An Investigation of the Connections among Convection, Clouds, and Climate Sensitivity in a Global Climate Model, J. Climate, 27, 1845–1862,, 2014. 

Zhao, M. and Austin, P. H.: Life Cycle of Numerically Simulated Shallow Cumulus Clouds. Part II: Mixing Dynamics, J. Atmos. Sci., 62, 1291–1310,, 2005. 

Zhao, M., Golaz, J.-C., Held, I. M., Guo, H., Balaji, V., Benson, R., Chen, J.-H., Chen, X., Donner, L. J., Dunne, J. P., Dunne, K., Durachta, J., Fan, S.-M., Freidenreich, S. M., Garner, S. T., Ginoux, P., Harris, L. M., Horowitz, L. W., Krasting, J. P., Langenhorst, A. R., Liang, Z., Lin, P., Lin, S.-J., Malyshev, S. L., Mason, E., Milly, P. C. D., Ming, Y., Naik, V., Paulot, F., Paynter, D., Phillipps, P., Radhakrishnan, A., Ramaswamy, V., Robinson, T., Schwarzkopf, D., Seman, C. J., Shevliakova, E., Shen, Z., Shin, H., Silvers, L. G., Wilson, J. R., Winton, M., Wittenberg, A. T., Wyman, B., and Xiang, B.: The GFDL Global Atmosphere and Land Model AM4.0/LM4.0: 2. Model Description, Sensitivity Studies, and Tuning Strategies, J. Adv. Model. Earth Sy., 10, 735–769,, 2018. 

Zheng, Y., Alapaty, K., Herwehe, J. A., Del Genio, A. D., and Niyogi, D.: Improving High-Resolution Weather Forecasts Using the Weather Research and Forecasting (WRF) Model with an Updated Kain–Fritsch Scheme, Mon. Weather Rev., 144, 833–860,, 2016. 

Zheng, Y., Rosenfeld, D., and Li, Z.: Sub-Cloud Turbulence Explains Cloud-Base Updrafts for Shallow Cumulus Ensembles: First Observational Evidence, Geophys. Res. Lett., 48, e2020GL091881,, 2021. 

Zhu, H., Hendon, H., and Jakob, C.: Convection in a Parameterized and Superparameterized Model and Its Role in the Representation of the MJO, J. Atmos. Sci., 66, 2796–2811,, 2009. 

Zimmer, M., Craig, G. C., Keil, C., and Wernli, H.: Classification of precipitation events with a convective response timescale and their forecasting characteristics, Geophys. Res. Lett., 38, L05802,, 2011. 

Zou, L., Qian, Y., Zhou, T., and Yang, B.: Parameter Tuning and Calibration of RegCM3 with MIT–Emanuel Cumulus Parameterization Scheme over CORDEX East Asia Domain, J. Climate, 27, 7687–7701,, 2014. 

Short summary
The paper provides a comprehensive review of the empirical values and assumptions used in the convection schemes of numerical models. The focus is on the values and assumptions used in the activation of convection (trigger), the transport and microphysics (commonly referred to as the cloud model), and the intensity of convection (closure). Such information can assist satellite missions focused on elucidating convective processes and the evaluation of model output uncertainties.