1. Introduction
Over large areas in the subtropics and Tropics, the vertical redistribution of boundary layer air in shallow cumulus clouds represents an important source of heat and moisture in the lower troposphere. Nitta and Esbensen (1974) and Ciesielski et al. (1999) attribute apparent heat sinks on the order of 4 K day−1 and apparent moisture sources on the order of about 3 to 6 g kg−1 day−1 in the Atlantic trade wind regime primarily to detrainment of shallow clouds near the trade wind inversion layer. Episodic trade wind regimes over the western Pacific warm pool are frequent enough to impact heating and moistening budgets and may be important in the formation of deep convective clouds (Johnson and Lin 1997). Apparently there is a need to take into account the effects of shallow cumulus clouds in atmospheric general circulation models (AGCMs) used in climate change studies and for weather prediction. Additionally, small cumulus clouds can be expected to be relevant to air chemistry and aerosol research owing to efficient entrainment and mixing of chemical species from the lower troposphere and negligible scavenging of chemical reaction products by precipitation.
As is well known, AGCMs are typically designed to resolve processes with horizontal scales of a few hundred kilometers, vertical scales of tens to hundreds of meters and temporal scales of hours and longer. While it is true that the actual time step length in many such models is only a few minutes, this is mainly for computational reasons and it is not uncommon to impose periodic or continuous temporal smoothing of prognostic variables using filtering operators to prevent development of undesirable fluctuations at the highest resolved frequencies (Asselin 1972). Thus the major problem associated with the treatment of convective clouds in AGCMs is that individual convective elements, and even organized convective systems, have characteristic temporal and spatial scales that are much smaller than those usually resolved by these models. Thus as in most parameterization problems, the most important net effects of moist convection must be represented in AGCMs in an appropriately averaged way.
Research on simplified representations of the basic processes in cumulus clouds that are relevant to large-scale modeling has been active for several decades. Various, often substantially different, parameterizations of the effects of convective clouds have been suggested and used to account for the effects of shallow convection in large-scale models. For example, adjustment schemes are based on the principle that convection acts to adjust the larger-scale atmospheric state toward one that is neutrally stable and/or consistent with observed structures in the Tropics (Betts 1986).
Many of the parameterizations in current use are based on the mass flux approach wherein vertical fluxes of such quantities as mass, heat, moisture, and other prognostic scalars associated with convective clouds are represented using simple cloud models, often based on the entraining plume concept (e.g., Tiedtke 1989). However, both observations and modeling of properties of cumulus clouds within the last 3 decades have raised doubts about the validity of the entraining plume concept as a model of shallow cumulus clouds. Observational studies from a wide range of field experiments suggest that the distribution of cloud properties at any given level in individual cumulus clouds is typically nonhomogeneous and may be comprised of mixtures of air parcels that have entered the cloud from many different levels in the convective layer, including air from environmental levels up to a few hundred meters above the level of observation (Paluch 1979; Blyth et al. 1988). For the most part these observations are inconsistent with the entraining plume concept and this has motivated introduction of the concepts of episodic entrainment (Raymond and Blyth 1986) to explain the observations. Associated conceptual models such as the shedding thermal model (Blyth et al. 1988; Blyth 1993) have been proposed as a mechanism for such entrainment and mixing processes to occur in clouds. In this model entrainment occurs in and around thermals in the rising tops of cumulus clouds. These thermals shed cloudy air in their wake. Updrafts and downdrafts and turbulence would typically also occur in the cloudy wake region and this could in principle result in further entrainment and detrainment. However, most of the available observations suggest that the properties of the air throughout the cloud remains nonhomogeneous and often gives the appearance of being comprised of mixtures of air from near or above the ambient level and almost undiluted air from the vicinity of the cloud-base region. Results of Blyth et al. (1988) for small Montana cumuli and other observations (e.g., Stith 1992) are consistent with this schematic model. Recent simulations with cloud resolving models such as those of Carpenter et al. (1998) lend some support to the shedding thermal concept, although there does not appear to be substantial support for this concept from laboratory studies (Morton 1997).
Whether in association with the circulation around thermals or in some other way, episodic entrainment and associated mixing can explain how cloudy air at a given level would appear to be comprised of mixtures of environmental air from above or near the observation level and air from lower levels, including the cloud-base region. However, there is also evidence that observed departures from such mixing line structures may be associated with buoyancy sorting (Taylor and Baker 1991) such that individual cloudy parcels may, after the initial mixing associated with episodic entrainment, remain sufficiently distinct to ascend or descend to levels in the cloud where they are neutrally buoyant with respect the the cloud-free environment. The episodic mixing and buoyancy sorting hypotheses are central to the parameterization schemes of Raymond and Blyth (1986) and Emanuel and Zivković-Rothman (1999). These authors assume that environmental air is mixed into a cumulus cloud in such a way as to give rise to a spectrum of air parcels with different mixing fractions within the cloud. These parcels then ascend or descend toward levels of neutral buoyancy. While these parameterizations are conceptually in better agreement with the aforementioned observations of cumulus clouds than are entraining plume models, the plume models do permit realistic representation of some key features of cumulus clouds in a simple way. For example Raga et al. (1990) found that the mean thermodynamic properties and mass flux in much of the cloud layer for Hawaiian trade wind cumuli could be represented with reasonable accuracy using an entraining plume model. However, significant departures from the entraining plume model predictions occur in the upper portion of the cloud layer.
While it is difficult to summarize the details of the life cycle of small cumulus clouds, there are observations and model results that suggest a relatively short and pronounced collapse phase near the end of the life cycle of small cumuli. For example, vertical mass fluxes of Hawaiian trade wind cumuli derived from dual-Doppler radar observations are characterized by rapid decay within a small fraction of the total cloud lifetime (Grinnell et al. 1996). Distinct decay of cloud properties may be attributable to penetrative downdrafts (Squires 1958) that form in the upper portions of clouds as a result of collapsing turrets and evaporative cooling. According to the model results of Carpenter et al. (1998), the narrowest and strongest downdrafts in a cloud develop in the later stages of the life cycle along the edges of ascending turrets that overshoot their level of neutral buoyancy (LNB). Raga et al. (1990) came to similar conclusions in an earlier and different study.
As noted above, parameterizations of shallow moist convection that are based on the steady entraining plume concept are currently in wide use in AGCMs. In these parameterizations, the in-cloud properties in individual plumes are generally assumed to be homogeneously distributed at each level. Thus they are conceptually in basic disagreement with many of the aforementioned observations within convective clouds. Another shortcoming of most current entraining plume parameterizations is the assumption of steady-state clouds. Betts (1975), Fraedrich (1974), and Cho (1977) have discussed the effects of accounting for recycling of cloudy air into the environment as cumulus clouds decay. The work of Fraedrich (1974) in particular generalized to some extent the earlier work of Kuo (1965). Including such effects can in principle account for important aspects of the observed large-scale moistening of the convective layer that occurs in association with shallow cumulus activity and the basic features of the observed cooling in the upper parts of the convective layers. However, these studies do not take into explicit account the effects of nonhomogeneity in cloud properties and make ad hoc assumptions about the recycling rate and/or the properties of the recycled air in order to obtain a closed form parameterization or deduce it from observed heat and moisture budgets. In order to account for cumulus life cycle effects in the context of steady entraining plume models it is also necessary to introduce ad hoc assumptions concerning detrainment processes (Siebesma and Holtslag 1996) or to employ a spectral cumulus ensemble model with component plumes having different entrainment characteristics as in Arakawa and Schubert (1974). Although an appropriate choice of the elements of the ensemble results in realistic simulations of the effects of convective clouds (Lin and Arakawa 1997b) it does not address the inherent limitations of the steady entraining plume model and requires introduction of additional closure conditions.
The motivation for the development of the parameterization presented in this paper is to account for the combined effects of inhomogeneities in the horizontal distributions of in-cloud properties and of finite cumulus lifetimes by generalizing the earlier approach of Fraedrich (1974). As in that approach and the earlier one of Kuo (1965) it is assumed that individual cumulus clouds grow until they reach their maximum heights and then rapidly decay with an associated lateral mixing of cloud properties into the environment. However, as important generalizations, we will account for nonhomogeneity in cloud properties arising from entrainment of environmental air in narrow regions near the rising tops of cumulus clouds and we will also derive a self-consistent expression for the detrainment associated with recycling of cloudy air at the ends of cumulus life cycles. Since our goal is to represent the time-averaged (over cumulus life cycles) effects of shallow cumulus clouds we will employ an idealized representation of cumulus clouds in which cloud properties are assumed to be steady except near the rising tops of the cumuli.
First tests and applications of the new parameterization are discussed in this paper. The parameterization is calibrated and characteristic features are investigated using prescribed vertical profiles of environmental temperature and water vapor mixing ratios from phase 3 (22 to 30 June 1969) of the Barbados Oceanographic and Meteorological Experiment (BOMEX; Holland and Rasmusson 1973). Comparisons are done for results from large eddy simulations (LES) and large-scale heat and moisture budgets derived from observations (Nitta and Esbensen 1974). In order to put the results of the new parameterization into broader perspective, we also discuss comparisons with the results of using the modified version of the Tiedtke (1989) parameterization proposed by Siebesma and Holtslag (1996).
2. Basic model of transient shallow cumulus clouds
The basic concept of the proposed parameterization is illustrated in Fig. 1. A cumulus cloud is initiated by turbulent processes in a conditionally unstable planetary boundary layer by lifting of air parcels from the mixed layer into the cloud layer where, under favorable conditions, the water vapor in the parcels becomes saturated and condensation occurs with continued ascent. Upon reaching the level of free convection (LFC), the cloud which forms ascends under the action of buoyancy forces. The buoyancy and other in-cloud properties are strongly affected by entrainment and detrainment processes at the top and sides of the rising cloud. To account for this we will introduce separate representations for mass entrainment rates at the rising top (Ei and Di) and those at the sides of the cloud below the rising top (Es and Ds).
The cloud top finally arrives at its maximum height near the level of neutral buoyancy. The ascent of the cloud top above the LNB is inhibited by the decrease in buoyancy, leading to the formation of penetrative downdrafts as illustrated by model results by Carpenter et al. (1998). We assume that these processes lead to a rapid decay of the cloud as suggested by observed vertical cloud mass fluxes by Grinnell et al. (1996). This decay process is simply approximated as an abrupt lateral detrainment of cloudy air into the environment (Df in Fig. 1) at the end of the growth cycle of the cloud. In reality the mixing of cloudy and environmental air in a decaying cloud would likely result in the formation of negatively buoyant mixtures that may then sink within the cloud before entering the environment. While doing so it is also possible that descending parcels may undergo further mixing with other cloudy parcels before finally entering the environment. These processes are ignored for the sake of simplicity in the parameterization of shallow convection that is developed here. The observations of Blyth et al. (1988) and the modeling results of Carpenter et al. (1998) and Lin and Arakawa (1997a) suggest that negatively buoyant cloud air parcels usually descend less than a few hundred meters during the cloud life cycle. However, ignoring the final buoyancy sorting phase that may take place before cloudy air is recycled into the environment results in important differences between the simple parameterization that is developed here and those that attempt to account for the main effects of buoyancy sorting (Raymond and Blyth 1986; Emanuel 1991; Emanuel and Zivković-Rothman 1999).
As noted in the introduction, the intent of the parameterization that is developed here is to account for the effects of shallow cumulus clouds on the larger-scale environment in terms of averages over the life cycles of individual clouds. To begin with we will introduce a mathematical description of the simplified cloud life cycle model that is the basis of the parameterization. As will be seen, the simplifying assumptions introduced below give rise to an easily interpreted and implemented mathematical representation of the time-averaged effects of shallow cumulus clouds.
a. Basic model equations


In the continuity equation for scalar properties χ, Eq. (2), Eχe refers to entrainment from the environment into the cloud, Dχc refers to detrainment of the cloud into the environment, and Sχ denotes in-cloud sources and sinks of χ. The scalar variables governed by Eq. (2) that are considered are the total water mixing ratio rt and the moist static energy h = cpT + Lυr + gz, where cp is the specific heat at constant pressure for dry air, T the temperature of the air, Lυ the latent heat of vaporization, r the water vapor mixing ratio, g the gravitational acceleration, and z the height. Effects of drizzle and radiation on small cumuli are usually small (e.g., Jiang and Cotton 2000) and the last term in Eq. (2) is omitted for rt and h in the following.
The vertical momentum budget in Eq. (3) includes terms which refer to detrainment of horizontally averaged vertical velocity wc, production of vertical momentum by buoyancy Bc, and effects of pressure perturbations Pc inside the cloud. No entrainment of vertical momentum is considered since it is assumed that vertical velocities are small in the environment compared to velocities in the cloud.
The mean perturbation pressure quantity Pc is retained to account for the dynamical response of the pressure field to horizontal inhomogeneity in the cloud vertical velocity field. This quantity is discussed further below but it should be noted that it does not account for all of the components that may be important in the pressure gradient force. As discussed by Yau (1979), the pressure gradient force in simple cumulus cloud models should typically have a component that partially offsets the buoyancy force and a dynamical component that is associated with the velocity field within the cloud. These two components are not distinguished here except for the part of the dynamical component that arises in response to the occurrence of inhomogeneity in the vertical velocity as indicated above. It is assumed that this dynamical component vanishes near the edge of the cloud. The reduction of the buoyancy force due to the pressure gradient force is accounted for in terms of a virtual mass parameter γ in the buoyancy term in Eq. (3), as is typical in one-dimensional cumulus models (e.g., Simpson and Wiggert 1969). The value of this parameter will be derived from results of LES models in section 3.


b. Cloud life cycle


For simplicity and without loss of generality we assume that the cloud begins its life cycle at time t = 0. The first Heaviside function in the square brackets in Eq. (5) represents the expansion of the cloud volume due to an ascending cloud top, the second represents the collapse of the cloud at t = τ at the time the ascending cloud top reaches its maximum height at z = zt. In Eq. (5), t∗ is the time elapsed between inception and the time at which the cloud reaches a certain height. Thus, because of the definition of the Heaviside function (zero when its argument is negative, unity otherwise), the fractional area and hence the mean vertical mass flux (Mc = ρawc) is nonzero only for t∗ ≤ t ≤ τ. Within this time period, however, it is assumed to be steady. As will be seen later, this fact has an important bearing on the form of the time-averaged mass flux.
The entrainment and detrainment rates at cloud top and an equation for t∗ are given below. In order to determine these quantities it is necessary to account for the role of spatial inhomogeneity within the rising cloud.
c. Effects of subcloud variability on vertical fluxes












While p(f) in Eq. (11) can in principle be derived from observations or cloud resolving model simulations, there is to our knowledge no published data on probability distributions of mixing fractions for shallow cumulus clouds. In order to obtain a sufficiently general parameterization of α for use in an AGCM, it would be necessary to specify the statistical quantities at different stages of the life cycles for many individual clouds. The alternative that is used in this study is to introduce physical assumptions and to test the sensitivity of ensemble-mean cloud properties and fluxes to these assumptions in the context of the BOMEX experiment. The probability density distributions that are used in associated sensitivity tests are described in section 3.
d. Solutions of basic equations
In the following, the approaches presented above are used to solve the fundamental set of equations of the parameterization [Eqs. (1) through (3)].














In order to determine τ = t∗(zt), it is simply assumed that the collapse of the cloud is initiated at that time at which the ascending cloud top reaches its LNB at z = zt. This assumption is consistent with the commonly used definition of the cloud top in the context of other parameterizations. However, this is not a necessary choice and other assumptions are in principle possible with regard to the application of the equations given above as long as wc > 0 at z = zt.
e. Organized entrainment and detrainment at lateral boundaries
Betts (1982) interpreted datasets by using convective pressure scales and saturation point coordinates. His results give evidence that the buoyancies of the air mixtures produced by cloud-top entrainment in typical trade wind cumulus regimes are too strong to produce downdrafts in clouds. To overcome this he suggests that a combination of cloud-top and lateral entrainment is required to explain up- and downdrafts in shallow clouds. The wide spread inclusion of lateral entrainment in entraining plume models stems from laboratory and field studies indicating the importance of lateral entrainment in penetrative convection.
Recent modeling studies give considerable insights into the nature of lateral entrainment and detrainment processes in shallow clouds (Bretherton and Smolarkiewicz 1989; Lin and Arakawa 1997a). Bretherton and Smolarkiewicz (1989) used a two-dimensional model to simulate an isolated long-lasting nonprecipitating cloud in a calm atmosphere with no surface forcing. Their results and an additional theoretical analysis show that lateral detrainment occurs preferably in regions with vertically decreasing buoyancy and entrainment in regions with increasing buoyancy in the cloud (Bretherton 1993). According to their results, buoyancy differences between cloud and the environment further away from the cloud produce circulations in the environment that tend to produce organized lateral flows toward the cloud in regions where the buoyancy increases with height and away from it where the buoyancy decreases with height.




The mixing rate scaling factor μ in Eq. (23) is a scalar that needs to be determined from closure assumptions. In our study, results from LES simulations and observations are used to determine μ.
f. Large-scale tendencies






3. Calibration and testing for the undisturbed BOMEX period
For calibration and testing of the new parameterization, results for the undisturbed period from 22 to 26 June 1969 during phase 3 of BOMEX are used. This period is characterized by little active deep cumulus convection (Nitta and Esbensen 1974). The BOMEX rawindsonde data and derived vertical heat fluxes indicate a transition to relatively active and organized convection after that period. The occurrence of a clear trade wind inversion, nearly steady-state conditions, and extensive observations are some of the features of this period that have been documented and studied in numerous scientific papers on the characteristics and effects of shallow clouds (e.g., Betts 1975; Soong and Ogura 1976; Esbensen 1978; Siebesma and Cuijpers 1995; Jiang and Cotton 2000).
The parameterization introduced in the previous section is applied to two independent data sets for the BOMEX experiment. The first dataset refers to results of 6 LES models. The second dataset refers to the analysis of large-scale heating and moistening budgets of Nitta and Esbensen (1974).
The results of the LES models provide more detailed information about convective processes than is available from current analyses of the BOMEX observations such as that of Nitta and Esbensen (1974). Hence the time- and space-averaged results of the LES models are used to calibrate the parameterization. In detail, the calibration is performed for the treatment of cloud-top mixing [i.e., to determine p(f), introduced in section 2c], and organized lateral mixing (i.e., to determine μ, introduced in section 2e). The results of the calibrated parameterization are then compared to additional LES results in order to evaluate its ability to consistently reproduce fluxes and cloud properties that are associated with active shallow cumulus in the LES simulations. Based on these results, the sensitivity of the results to the value of γ (section 2a) and to the treatment of cloud-top detrainment (section 2d) are studied.
The second important application of the parameterization involves direct comparison of the results of the parameterization and the analyzed large-scale heat and moisture budgets as determined by Nitta and Esbensen (1974).
For both datasets, the relationships given in the previous sections are used to diagnose the ensemble-mean shallow cumulus cloud fields and budgets that occurred during the undisturbed period of BOMEX. Other processes, such as turbulent diffusion in the boundary layer or cloud–radiation interactions, are not considered here.
a. Statistical distribution of cloud properties
Three versions of the parameterization are tested in order to study different limiting cases with regard to the treatment of cloud inhomogeneities.
The first and simplest of the cases is that the conserved in-cloud properties are top-hat distributed; that is, there is no horizontal variability in the cloud. Accordingly, there is no cloud-top entrainment and p(f) [Eq. (11)] is simply a delta function, p(f) = δ(f). The version of the parameterization that corresponds to this assumption is labeled TOPHAT in the following in order to distinguish the results.
In the second case it is assumed that all mixtures have equal probabilities in the cloud at each level, that is the mixing fraction is uniformly distributed, which is expressed by p(f) =




b. Alternative parameterization of shallow cumulus






c. Application to LES results
In this first study, data from an intercomparison project for LES models are used. This project was performed by members of the GEWEX Cloud System Studies (GCSS) Boundary Layer Cloud working group, one of the activities of the Global Energy and Water Cycle Experiment (GEWEX). In total 11 LES models from various international groups were applied to the period 22 to 23 June in BOMEX. Extensive results1 for mass fluxes and cloud properties of clouds in the simulations have been presented at a group workshop from 24 to 25 July 1997 in Seattle, Washington (Siebesma et al. 2001, manuscript submitted to J. Atmos. Sci.; Siebesma and Cuijpers 1995). From the participating models, results of those 6 models are selected that are based on the Boussinesq approximation in order to be consistent with the current parameterizations. These 6 models are from the Royal Netherlands Meteorological Institute (KNMI; Cuijpers and Duynkerke 1993), the Max Planck Institute for Meteorology (MPI; Chlond 1992, 1994), the National Center for Atmospheric Research (NCAR; Moeng 1984, 1986), the Met Office (UKMO; MacVean 1993; Brown et al. 1994), the University of Oklahoma (UOK; Kogan et al. 1995), and the West Virginia University (WVU; Lewellen et al. 1996).
The size of the three-dimensional LES model domain is 6.4 km in both horizontal and 3 km in the vertical directions. The grid sizes are Δz = 100 m for the horizontal grid and Δz = 40 m for the vertical grid. Periodic boundary conditions are used at the lateral boundaries. A sponge layer above the inversion layer damps out any perturbations in the simulations. At the bottom, prescribed momentum, heat, and moisture fluxes are applied. Additional forcing in the simulations is applied with regard to large-scale subsidence, radiative cooling, and large-scale horizontal advection of moisture. The simulations are initialized with time-averaged vertical profiles of temperature and specific humidity that have been obtained by the Oceanographer, the ship at the northeast corner of the BOMEX domain. Random fluctuations of the virtual potential temperature and total water mixing ratio are used as initial perturbations at the lowest 40 levels. Hourly simulations results are available for a time period of 4 h after an initialization phase with a duration of 2 h. In this study we are only referring to results for visible clouds (i.e., nonzero LWC) in the cloud layer.
The mean vertical profiles of potential temperature and total water mixing ratio in the simulations, averaged over the horizontal domain and the simulation period, are shown in Fig. 2. The results show a well-mixed layer with a depth of about 500 m and an inversion above about 1500 m. The spread between the results of different models is largest in the inversion layer. The individual profiles are averaged in order to obtain vertical profiles for the applications of the parameterizations (Fig. 2). For this, the profiles of mean potential temperature and total water mixing ratio are provided on a grid with an increasing grid size with altitude. The grid size increases from Δz = 48 m near the surface to Δz = 121 m for the level centered at z = 592 m near the top of the boundary layer and to Δz = 340 m for the level centered at z = 2186 m. This resolution is close to that used in the third generation AGCM of the Canadian Centre for Climate Modelling and Analysis (CCCma), which uses a hybrid pressure-based vertical coordinate. The grid size is also typical of other such models.
The cloud base zb is determined as that level at which the mass fluxes and cloud cover are at maximum in the LES results, which is at about z = 600 m. In the applications of the parameterization, zb = 606 m is used, corresponding to the lifting condensation level (LCL) of a parcel that is lifted from the 5th level above ground at z = 470 m.
1) Calibration of subcloud vertical transport processes
The first issue that is addressed in the context of the LES results is how well cloud inhomogeneities and associated vertical fluxes can be represented by the parameterization. In particular, the goal is to find an approach for p(f) that produces a reasonably close match between the values of α from the LES results and the parameterization. For this purpose, two profiles of α are calculated for the mean vertical fluxes of liquid water potential temperature θl and the total water mixing ratio rt and from the mean vertical mass flux in the simulated LES cloud ensembles, using Eqs. (6) and (7). In the following, the versions DECORE and EQPROB of the parameterization, as described in section 3a, are applied and the results compared to the results of the LES models for α.
Figure 3 shows results for α from the LES models KNMI, MPI, NCAR, and UOK (thin lines). The results of the models UKMO and WVU are considerably different from those displayed in the figure. Additionally, the vertical profiles of
The small values of α obtained from the LES models near cloud base are consistent with top-hat distributed in-cloud properties at these levels. This suggests that downdrafts are not likely to penetrate to levels near cloud base. Above cloud base, α increases monotonically with height in all four simulations. Although the spread in the results also increases with height it is apparent that α has its maximum somewhere near the middle of the cloud layer. The cause of the scatter between the results of the individual models may at least in part be caused by different vertical velocities in the cloud ensembles in each simulation. The lowest values of wc are simulated in the MPI model (results are later shown in Fig. 7). A lower value of wc may be indicative for a stronger erosion of the cloud tops and hence larger values of α. This appears to be likely since the ascending cloud tops are longer exposed to mixing processes if the cloud top ascends slowly. Differences in the results near the base and the top of the cloud layer may also be caused by too short sampling times since the results at these levels are dominated by sporadic events.
Results of the parameterization in Fig. 3 are shown for version DECORE of the parameterization for φ = 1 × 10−3 s−1. There is good agreement between the results of the LES models and DECORE, especially near cloud base with this particular fit of φ to the data. Results near the inversion are highly uncertain due to large differences between the results of the LES models. In contrast, EQPROB tends to overestimate the cloud inhomogeneities throughout the cloud layer.
The results of EQPROB, DECORE, and TOPHAT (i.e., α = 0) cover roughly the range of uncertainty for α. Therefore, results of these versions of the parameterization can be used to estimate the sensitivity of the approach to this parameter.
2) Calibration of organized lateral mixing




For top-hat distributed and conserved cloud properties, σ represents a unique measure of the degree of dilution by organized lateral mixing processes above cloud base. For all other distributions, σ can only be uniquely related to the efficiency of the organized lateral mixing if d = 0 according to Eq. (28). It can also be shown that σ has the same value for all conserved in-cloud properties that fall along a single mixing line. Further, σ is unity if no organized lateral entrainment or detrainment occurs; that is, e = d = 0. It can be easily shown that σ equally measures the efficiency of organized lateral entrainment in the approach of SH1996 with
The application of Eq. (30) to the LES results and to the parameterizations leads to the vertical profiles of the excess flux fractions σ for liquid water potential temperature and total water mixing ratio shown in Fig. 4. The results for the models KNMI, MPI, NCAR, and UOK are very similar. Decreases of the values with height in the lower portion of the cloud layer correspond to reductions in
The corresponding results of the models UKMO and WVU are in considerable disagreement with the other LES models. These results are not shown in Fig. 4.
Values of σ, calculated with Tiedtke's parameterization and the modifications suggested by SH1996, are in good agreement with the LES results near cloud base. However, the results give evidence for much too strong lateral entrainment in the upper portion of the cloud layer (Fig. 4). Additionally, the parameterization does not reproduce the characteristic minimum at about z = 1200 m.
In contrast to the results of Tiedtke's parameterization, the parameterization of organized lateral entrainment and detrainment in terms of buoyancy gradients according to Eq. (23) gives very good agreement with the LES results at all levels for a fit with μ = 14 s2 m−1 to the data. Organized lateral entrainment occurs below z = 1473 m, and organized lateral detrainment above this level in this application of the parameterizations. Hence results for σ in Fig. 4 above this level are dependent upon the values of α at these levels [cf. Eq. (28)]. However, owing to relatively small values of α at these levels, results with d = 0 are very similar to the results shown in Fig. 4.
3) Comparisons of in-cloud properties and fluxes
With the values of μ and α given by calibration, the different versions of the parameterization can be tested with regard to a consistent representation of the mass fluxes and distributions of conserved properties in the simulated clouds.


Results for all 6 LES models and for the parameterizations are displayed in Fig. 5. The vertical profiles are similar in shape to those in Fig. 4 but clearly have lower values at any level above cloud base. This means that inhomogeneities in the distributions of in-cloud properties cannot be neglected and that wcχc <
The overall best approximation to the LES results according to Fig. 4 is achieved by the results of version DECORE. The values are within in the range obtained from the LES results at all levels.
Excessive values of ν are computed in the version TOPHAT and in the approach of SH1996 near cloud base. As already apparent in the results for σ, SH1996's approach also leads to large differences in the profiles of the cloud-averaged values of θl and rt in the upper portions of the cloud ensemble. These findings for SH1996 are generally consistent with the results of Raga et al. (1990). They applied a simple entraining plume approach with constant fractional entrainment rates to cloud bands observed off the coast of Hawaii, similar to the treatment of mixing in the approach of SH1996. Their results give evidence for good agreement between simulated and observed wet equivalent potential temperature and liquid water content near cloud base. However, the simulated values are considerably lower than observed in the upper portion of the cloud layer.
The LNB derived from the results shown in Fig. 5 is at zt = 1720 m for DECORE, zt = 1552 m for EQPROB, zt = 2229 m for TOPHAT, and zt = 1593 m for SH1996. As pointed out before, the diagnostic results for the LES models are not representative near the top of the cloud layer owing to low statistical significance and uncertainties in the interpolates. Although these results cannot be used in an interpretation they are included in the figure in an attempt to present the complete diagnostic dataset.
Results of EQPROB and DECORE for the mean liquid water content agree well with the LES results (Fig. 6). More vigorous mixing at cloud top in EQPROB compared to DECORE leads to lower liquid water mixing ratios in this version of the parameterization (cf. Fig. 3). Similar to the results of EQPROB and DECORE, LES models that produce larger values of α also tend to produce smaller liquid water contents (e.g., compare results of MPI and NCAR in Fig. 3 and Fig. 6). In contrast, the omission of cloud-top mixing processes in TOPHAT and SH1996 leads to overall much larger liquid water contents. Values produced by these parameterizations are in considerable disagreement with the results of the LES models near cloud base. However, in the upper portion of the cloud layer, SH1996 produces unrealistically low values. These results give support to the earlier notion that the assumption of a constant fractional entrainment rate in SH1996 is not very realistic. Additionally, results in Fig. 6 give clear evidence for the importance of cloud-top mixing processes for the treatment of liquid water contents in shallow cumuli.


Generally, fluxes and cloud properties at the lower boundary of the cloud, such as the mean vertical velocity at cloud base, need to be determined through introduction of an additional closure assumption (Nitta 1975). As mentioned earlier, the relevant processes in the boundary layer that determine the magnitude of the vertical velocity at cloud base are not considered in the present study. Here wb = 0.3 m s−1 is used as the lower boundary condition for wc at z = zb.
As discussed earlier, there is considerable uncertainty with regard to the magnitude of γ. As a first approach, γ = 5 is used in all parameterizations. This choice produces good agreement between most of the mean parameterized vertical velocities and mean results from the LES models near cloud base within the spread between the LES results (Fig. 7a), despite being much larger than the value typically used in one-dimensional cumulus models (Simpson and Wiggert 1969).
The large sensitivity of the mean vertical velocities to the value of γ is illustrated in Fig. 7. In Fig. 7a, γ = 5 is used in the applications of the parameterizations. Compared to the simulation with γ = 0 (Fig. 7b), the mean vertical velocities are typically smaller by more than 50% in the results of the parameterizations.
However, the parameterized vertical velocities in DECORE and EQPROB compare well to the mean LES results obtained for positively buoyant cloudy air for γ = 0 (Fig. 7b). A comparison of the LES results in Figs. 7a,b give evidence for large contributions of negatively or neutrally buoyant air to the mean vertical velocities in the simulated clouds. An analysis of LES results for the mean fractional area occupied by positively buoyant air in the cloud gives evidence that 55% to 75% of the air in the clouds is positively buoyant below z = 1250 m (not shown). In contrast, an integration of the probability density function over the positively buoyant region in the clouds simulated in DECORE yields that more than 97.5% of the air in the cloud is positively buoyant at all levels below z = 1252 m. Above that level, the fraction of positively buoyant air decreases monotonically with height to 63% at the top of the cloud layer.
One possible interpretation of the differences between the mean vertical velocities and the amounts of positively buoyant air in DECORE and the LES results appears to be that negatively and neutrally buoyant air parcels are produced outside the active shallow cumulus clouds or during the final stage of the life cycle of the cumulus clouds in the LES simulations (see also the introduction and section 2). The mean vertical velocities associated with air parcels that occur after the growth phase of the cumulus life cycle (t > τ) can be expected to be negative or at least much smaller than the vertical velocities shown in Fig. 7b. Owing to the large relative humidities in the LES simulations (e.g., between RH = 95% to 96% in a layer between z = 450 and 1000 m) these air parcels may remain saturated for a long time period relative to the cloud lifetime τ in DECORE, EQPROB, and TOPHAT. However, according to the simplified cloud life cycle in these parameterizations, effects of any saturated air parcels that exist after the collapse of the cloud at t = τ are omitted. Parameterization results shown in Fig. 7 only refer to mean vertical velocities within the cumulus clouds at t ≤ τ. Hence mean buoyancies and vertical velocities in the parameterized clouds may be large compared to the LES results in Fig. 7a. However, the good agreement with results for the mean cloud properties (Fig. 5) and LWC (Fig. 6) suggests that these important parameters are not strongly affected by processes occurring in the final stage of the cloud life cycles.
Unfortunately, there are no LES results available to directly test the importance of air parcels that exist in the later stages of the shallow cumulus life cycle. Therefore, the mean LES results can only be used to determine the range of γ in the parameterization that is consistent with the LES results.
Another interpretation of the differences in Fig. 7 is that drag and perturbation pressure gradients at the edges of the clouds may have considerable impacts on the vertical momentum budget and that these effects are not sufficiently accounted for.
Figure 8a shows the time-averaged fractional cloud cover 〈a〉, normalized by cloud base mass flux Mb. The normalization allows a direct comparison between results of the LES models and the parameterizations without the need to evaluate Mb. It is evident that the assumptions made in SH1996, DECORE, and TOPHAT lead to very accurate representations of the fractional cloud cover with the choice γ = 5. For γ = 0, substantially lower cloud amounts are calculated in DECORE (Fig. 8b). This figure also includes results of a simulation with DECORE in which cloud-top detrainment is included with Di = Ei (section 2d), showing that the cloud amount is not very sensitive to the treatment of cloud-top detrainment.
Relatively large differences exist between the results of the LES models for the mean vertical mass fluxes normalized by Mb (Fig. 9). The results of EQPROB and DECORE are well within the range of the predicted values from the LES results. Large differences occur for TOPHAT owing to a large vertical extent of the clouds simulated with this version of the parameterization. A major difference between SH1996 and the other results is that the vertical mass flux drops from a relatively large value at cloud top to zero directly above the cloud top in this parameterization whereas the other profiles are smooth near cloud top. The discontinuity can be attributed to the omission of life cycle effects in SH1996's approach. The assumption of a time-dependent fractional cloud cover results in more realistic profiles near the top of the cloud layer.
The mean vertical mass flux calculated in DECORE is insensitive to the value of γ (Fig. 9b). Only slightly smaller vertical mass fluxes are obtained with γ = 0 compared to the reference simulation with γ = 5. The reason for the low sensitivity is that t∗ and τ scale in a similar way with γ. Hence Mc is insensitive to differences in t∗ and τ (see section 2). Similar to the cloud amount, the introduction of cloud-top detrainment does not lead to considerably different values of Mc (Fig. 9b).


Equations (32) and (33) can be directly applied to the LES results in order to determine 〈E〉
Results of the approaches DECORE and TOPHAT in Figs. 10a,b are generally in good agreement with the LES entrainment fluxes. The overall best agreement is found for version DECORE. EQPROB produces excessive entrainment and detrainment fluxes.
A striking feature of the diagnosed detrainment flux for SH1996's parameterization is that a large fraction of the detrainment is confined to a narrow layer near the top of the cloud. This characteristic feature becomes even more evident at higher vertical grid resolution. The reason for this is that the detrainment at the top of the cloud layer is a delta function in this approach that is determined by the magnitude of the vertical mass flux at cloud top.
Similar to the results shown in Fig. 9, the entrainment and detrainment fluxes in DECORE are not very sensitive to the value of γ (Fig. 10c) and the formulation of cloud-top detrainment (Fig. 10d).
The relative contribution of the cloud-top entrainment fluxes Ei to the time-averaged ensemble-mean entrainment mass flux, 〈E〉 is shown in Fig. 11. According to the results for parameterization versions EQPROB and DECORE, the contribution of cloud-top entrainment is negligible near cloud base. In contrast, cloud-top entrainment is the predominant entrainment process in the upper portion of the cloud layer. At cloud top, the dilution of the cloud occurs exclusively by cloud-top entrainment. Organized lateral entrainment does not contribute to the total entrainment because the vertical gradient of the undiluted parcel buoyancy is negative near cloud top.
If accounted for, the contribution of cloud-top detrainment to the total detrainment is much smaller than the contribution of final detrainment. Less than 26% of the total detrainment is associated with Di in the simulation with DECORE and Di = Ei (Figs. 8 to 10).
d. Large-scale heat and moisture budgets in BOMEX
In this section, the parameterization is applied to the mean observed meteorological situation in the BOMEX domain. The results are compared with the analyses of large-scale heat and moisture budgets of Nitta and Esbensen (1974).
Only results for 22 and 23 June are considered since the fluxes during these days are very similar whereas the situation changes considerably on 24 June (Soong and Ogura 1976). The quasi-steady character of the situation is illustrated by the mean vertical profiles of potential temperature and total water vapor mixing ratio in Fig. 12. The profiles shown in this figure refer to time- and station-averaged results of rawindsonde soundings that were taken 15 times per day at each of the four corner ships in the experiment. The measurements were checked for inconsistencies and daytime humidity errors due to radiative heating of the hygristor elements of the rawindsondes were additionally corrected (Nitta and Esbensen 1974). There is good agreement between the measurements and the results of the LES models for the total water mixing ratios. However, the temperature is about 1 to 2 K lower in the lower part of the LES simulation domain, which is associated with slightly lower potential temperatures in the LES results (Fig. 12a).


In the applications of the parameterizations, a grid with a vertical grid size of 10 hPa is used. The cloud-base mass flux is prescribed to Mb = 288 hPa day−1 according to the value obtained by Nitta (1975) for the period 22 to 24 June in BOMEX. The parameterizations are separately applied to each day but only the mean results for the period 22 to 23 June are used in the comparisons.
In contrast to the comparison with the LES results in the previous section, it is not possible to directly obtain zb from the analysis results. In order to address the uncertainties in zb and cloud-base properties, the properties of the air at cloud base are estimated for the undiluted ascent of air parcels that are launched at different levels between 1 and 9 hPa above the sea surface at each day. The vertical distance between the launching levels is 1 hPa.
Mean vertical profiles of Q1 − QR for 22 and 23 June are shown in Fig. 13 using a given profile of QR (Holland and Rasmusson 1973). The results of Nitta and Esbensen (1974) give evidence for a maximum of the large-scale cooling slightly below z = 2000 m with a value of about 5 K day−1. According to Nitta and Esbensen (1974), this maximum can be attributed to reevaporation of cloud droplets that detrain near the top of the trade inversion. The small values above about z = 2300 m indicate that there are very few clouds that penetrate into the layers above the inversion. Positive values in the layer from about z = 600 m to about z = 1400 m may be caused by adiabatic warming due to compensating downward motion in the cumulus environment. Below about z = 600 m, positive values are caused by turbulent motions of clear air in the mixed layer that are not considered in the parameterizations. Hence the following discussion refers to the layer above z = 600 m only.
All parameterizations introduced in the previous sections consistently produce much too low cloud tops even for air parcels launched very close to the surface (Fig. 13). The deepest cloud layers are produced in TOPHAT. But even in this simulation, the cloud tops are below about z = 1800 m, which is much lower than the representative cloud top according to the profiles of analyzed apparent heat source and moisture sink.
In order to achieve some understanding of the considerable disagreement between the parameterizations and the BOMEX analyses it is useful to note that the cloud tops of the bulk of the clouds in the cloud ensembles appear to be below z = 1800 m in all LES model results according to Figs. 9 and 10. This is much lower than the corresponding heights according to the analysis results in Fig. 13. This conclusion is consistent with the results of Siebesma and Cuijpers (1995) who obtain maximum cooling rates at about z = 1600 m in simulations with the KNMI model for BOMEX. Their results show substantial underestimations in the magnitudes of Q1 and Q2 above this level and slight overpredictions below. From the calibration of μ and α in the parameterization based on the LES results, it follows that biases in the predicted fluxes in the LES models are inherited by the parameterization through the calibration process.
It appears that the differences in the large-scale temperatures in the LES model simulations and in the analysis of Nitta and Esbensen do not explain the low cloud tops in the LES simulations (Fig. 12). In contrast, the lower temperature in the mean LES results and associated higher relative humidities most likely tend to produce too high tops in the LES models rather than too low cloud tops compared to the BOMEX observations. This conclusion is supported by results of sensitivity studies using the parameterization described in the previous sections for different large-scale temperature profiles. According to the results of the sensitivity studies, the height of the cloud tops is typically reduced by up to a few hundred meters if the large-scale temperature is increased by 2 K in simulations with the parameterizations if the total water mixing ratio is kept constant. The differences in the prescribed large-scale values of temperature explain to a large extend the differences between the results of the parameterizations in Fig. 13 and results presented in the previous sections.
While it is not clear from the results available in this study what causes the low cloud depths in the LES results for BOMEX, it needs to be emphasized here that entrainment processes simulated in LES models may be sensitive to the parameterizations of subgrid-scale processes in these models, to the conditions at the boundaries of the model domain, and to the numerical approaches. For example, results of a different LES model intercomparison project of the GCSS Boundary Layer Cloud working group give evidence that LES models tend to produce much too large entrainment rates if the variabilities in height and thickness of strong stratocumulus inversions are not adequately resolved by the model grid (Bretherton et al. 1999; Stevens and Bretherton 1999). Additionally, Brown (1999) found a strong sensitivity of the sizes of individual shallow clouds in LES model simulations for BOMEX but also an encouraging low sensitivity of ensemble-averaged statistics to grid sizes.
Since it is not possible to determine the causes for the discrepancies in Fig. 13 due to the lack of data on the properties of the clouds in BOMEX, an alternative is to recalibrate the parameterization by finding new values of the free parameters in the parameterization that lead to better agreement between the large-scale budgets. The approach pursued here is simply to match the levels at which the maxima in the heating and moistening profiles occur in the analysis results in order to obtain first approximations. This approach is not completely different from the approaches in other studies of the large-scale heat and moisture budgets in BOMEX (e.g., Soong and Ogura 1976).
Results of different parameterizations with new calibrations are shown for the apparent heat source (Fig. 14) and the apparent moisture sink (Fig. 15). The parameters used in the simulations are summarized in Table 1. Results of the version DECORE_A in Figs. 14a and 15a refer to simulations without any organized lateral entrainment or detrainment but otherwise with the same treatment of cloud-top entrainment as in DECORE. TOPHAT_A in Figs. 14b and 15b does not include any cloud inhomogeneities. The organized lateral mixing in TOPHAT_A is relatively weak compared to DECORE. DECORE_B in Figs. 14c and 15c includes cloud-top entrainment and organized lateral mixing. The magnitude of the entrainment and detrainment fluxes is small compared to DECORE. Finally, Figs. 14d and 15d show results of EQPROB_A for a parameterization that uses the same treatment of cloud inhomogeneities as in EQPROB but that does not include any organized lateral mixing.
For DECORE_A, excessive cooling is diagnosed at heights below about z = 1500 m and insufficient cooling above this level (Fig. 14a). Only results for parcels launched at the levels between 5 and 9 hPa are shown for DECORE_A. Parcels launched at the first two levels produce clouds that overshoot the inversion and cannot be considered as part of a shallow cumulus ensemble.
Very good agreement between analyzed and mean simulated heating rates are obtained for TOPHAT_A (Fig. 14b). Compared to the results of DECORE_A, TOPHAT_A produces very different cloud-top heights for different parcel launching levels. The mean results above about z = 2000 m are strongly affected by the results for the first parcel launching level above ground.
The approach used to determine φ and μ in DECORE_B is to fit the results of DECORE_B to DECORE with respect to the relative contributions of cloud-top entrainment to the total entrainment at all levels in the clouds. Hence the relative magnitudes of cloud-top and organized lateral entrainment are similar in DECORE_B and DECORE. The results of this procedure are shown in Fig. 16.
Good general agreement with the analysis is found for the results of DECORE_B (Fig. 14c). The results for each parcel launching level are similar to those for slightly higher parcel launching levels in TOPHAT_A (e.g., compare results for 1 hPa in DECORE_B to 3 hPa in TOPHAT_A). It appears that parcels lifted from 1 hPa produce the best agreement with the observations.
The results of EQPROB_A give evidence that the assumption of equal probabilities for the air mixtures in the cloud generally produces too low cloud tops (Fig. 14d). The results of this parameterization are similar to those of DECORE_A. These results suggest that cloud-top entrainment, without additional consideration of organized lateral entrainment, tends to produce excessive cooling near the bottom of the cloud layer and deficits near the trade inversion.
Differences between the results of the parameterization and the analysis may in part be attributable to variability in the large-scale thermodynamic profiles which are not accounted for in the parameterization. Also, uncertainties associated with MB may affect the magnitudes of Q1 − QR and Q2. However, the magnitude of MB does not effect the relative magnitudes of the results shown in Figs. 13 to 15 [from Eq. (24)].
Finally, results for the apparent moisture sink shown in Fig. 15 are consistent with the results in Fig. 14 in the upper portion of the cloud layer. However, systematic overpredictions occur for the moistening rates in a layer that is centered at about z = 1100 m for all parameterizations. Although this discrepancy may indicate problems in the approaches it is interesting to note that the results of Soong and Ogura (1976), which are derived using a time-dependent, two-dimensional, axisymmetric cloud model, also show a similar minimum in Q2 at about the same height for the same BOMEX period. Soong and Ogura (1976) speculate that the local maximum at this height in Nitta and Esbensen's data may be caused by the correction method they use for the humidity data. Also, Holland and Rasmusson (1973) report different results for Q2 that they obtain with a different correction method. In particular, Holland and Rasmusson (1973) also find a local minimum in Q2 at about the same height.
4. Discussion and conclusions
A new parameterization of the bulk effects of transient shallow cumulus clouds is proposed. In the parameterization, the physical properties of rising air parcels are modified by entrainment of environmental air into the ascending top of the cloud and also by organized entrainment at the lateral boundaries of the cloud. The cloud-top mixing is associated with horizontal inhomogeneities in the properties and fluxes in the cloud. Different physical hypotheses are introduced in order to obtain simple representations for the effects of cloud inhomogeneities in the parameterization. Finally, the initial growth phase of the cloud in the parameterization is terminated by an abrupt and complete collapse of the cloud during which the air remaining in the cloud is detrained into the environment. Based on these assumptions, budget equations for the vertical mass flux, conserved cloud properties, and vertical velocity are derived in order to represent the mean effects of a shallow cumulus clouds averaged over their life cycles.
The new parameterization is applied to observations and LES results from the undisturbed period in BOMEX from 22 to 23 June 1969. Application of the parameterization to LES results for this period produce consistent profiles of mean fluxes and properties in the cloud ensembles for version DECORE of the parameterization. It is found that entrainment at ascending cloud tops is the largest contribution to the total entrainment flux near the trade inversion in the LES simulations for BOMEX. Additionally, the results show that penetrative downdrafts do not affect the in-cloud properties and fluxes near the base of the clouds in the simulations.
Including the simple life cycle approach in the parameterization produces good agreement with LES results with respect to fluxes and properties near the top of the cloud layer. The contribution of the final detrainment flux near the end of the cloud life cycle is by far the largest contribution to the total detrainment flux. In contrast, the contribution of cloud-top detrainment is negligible.
Additional tests of large-scale heat and moisture budgets for the same period in BOMEX using results of an analysis by Nitta and Esbensen (1974) are also encouraging. It is shown that the proposed parameterization produces realistic cooling and moistening profiles below the inversion. However, these comparisons give no conclusive evidence regarding the nature of the entrainment process although it appears that mixing by organized lateral entrainment is important for the large-scale budgets. There is also evidence that the entrainment processes in observed clouds in BOMEX are less efficient than suggested by the LES results.
As noted in section 2 above, the effects of buoyancy sorting are ignored in the parameterization proposed herein. Consequently certain aspects of the parameterization of large-scale heating and moistening effects are qualitatively different from current parameterizations based on the buoyancy sorting concept as applied to shallow cumulus regimes. An important difference is that in parameterizations based on the buoyancy sorting hypothesis the air, which is detrained into the environment, is nearly neutrally buoyant and free of liquid water (Emanuel 1991; Emanuel and Zivković-Rothman 1999). A result of this is that detrainment will have a less important effect on the large-scale heat and moisture budgets than is the case in our parameterization. In particular, since the mean mass flux is always positive in our parameterization, all of the net cooling at top of the cloud layer is accomplished through evaporation of detrained liquid water. For the buoyancy sorting parameterizations this cooling is achieved predominantly through the action of negatively buoyant downdrafts that arise as a result of the buoyancy sorting mechanism. Cho (1977) also assumed that air that is finally detrained by shallow convective clouds at the end of their life cycles is free of condensed water and has the same temperature as the environment. Consequently in applying his representation to the BOMEX data he determined that a net downward convective mass flux is required in the upward part of the convective layer to account for the net cooling in that region.
It is noteworthy that the LES models for the most part produce small but positive values of the net convective mass flux in the region just below the trade inversion. The parameterization proposed herein is in good agreement with the LES results. However, the LES models may also suffer from deficiencies that limit their ability to adequately represent the true behavior of negatively buoyant downdrafts that typically occur in reality in the upper part of the convective layer. Thus the absence of a net downward mass flux in that region in the LES simulations does not negate the importance of the net effects of buoyancy sorting in that region. However, despite this possible shortcoming we have demonstrated that the parameterization herein is capable of reproducing the observed large-scale heat and moisture budgets in shallow convective regimes with considerable fidelity.
This study is entirely based on the BOMEX case, which also has been used in many other studies on processes associated with shallow convection. The reason for the decision to use BOMEX results in the present study was motivated by the fact that the datasets for this experiment appear to be relatively detailed and well documented in comparison to other studies. However, further investigations are required in order to validate the parameterization under different conditions.
A major problem in this study is that the calibration and validation of the parameterization is relatively difficult owing to the fact that insufficient data with respect to natural shallow convection is available from observations and cloud resolving model simulations. For example, more detailed observations of liquid water contents and total water mixing ratios at different heights and at different stages of cloud life cycles would be very useful in future development and testing of parameterizations. The main problem regarding the parameterization of cloud-top entrainment processes in the current approach is that no data are available that can be used to directly determine the probability density distributions of the air mixtures in the ascending cloud tops. Although the present study and numerous other studies strongly suggest that inhomogeneities in the distributions of air mixtures are important features in shallow clouds there are still large uncertainties with regard to the treatment of cloud inhomogeneities in the context of a parameterization.
The parameterization described in this paper is currently being tested in the fourth generation AGCM of the CCCma in order to study the effects of shallow convection on simulations of the global distributions of heat, moisture, and sulfur compounds. Results of these applications and the treatment of cloud base mass fluxes (i.e., the cloud base closure) will be reported in a following paper.
Acknowledgments
We are grateful to two anonymous reviewers for their very helpful comments and to Steven Esbensen who provided BOMEX large-scale budgets and profiles of environmental meteorological parameters. Also, we wish to thank Pier Siebesma for giving us permission to use results of the GCSS Boundary Layer Cloud working group model intercomparison project for the BOMEX shallow cumulus study. Thanks are further due to our colleagues John Scinocca, Virginie Lorant, and Richard Harvey for comments and discussions.
This research was supported by the Meteorological Service of Canada through the Climate Change Action Fund (CCAF).
REFERENCES
Arakawa, A., and W. H. Schubert, 1974: Interaction of a cumulus cloud ensemble with the large-scale environment, Part I. J. Atmos. Sci., 31 , 674–701.
Asselin, R. A., 1972: Frequency filter for time integrations. Mon. Wea. Rev., 100 , 487–490.
Betts, A. K., 1975: Parametric interpretation of trade-wind cumulus budget studies. J. Atmos. Sci., 32 , 1934–1945.
Betts, A. K., . 1982: Cloud thermodynamic models in saturation point coordinates. J. Atmos. Sci., 39 , 2182–2191.
Betts, A. K., . 1986: A new convective adjustment scheme. Part I: Observational and theoretical basis. Quart. J. Roy. Meteor. Soc., 112 , 693–709.
Blyth, A. M., 1993: Entrainment in cumulus clouds. J. Appl. Meteor., 32 , 626–641.
Blyth, A. M., W. A. Cooper, and J. B. Jensen, 1988: A study of entrained air in Montana cumuli. J. Atmos. Sci., 45 , 3944–3964.
Bretherton, C. S., 1993: Understanding Albrecht's model of trade cumulus cloud fields. J. Atmos. Sci., 50 , 2264–2283.
Bretherton, C. S., and P. K. Smolarkiewicz, 1989: Gravity waves, compensating subsidence and detrainment around cumulus clouds. J. Atmos. Sci., 46 , 740–759.
Bretherton, C. S., and Coauthors. 1999: An intercomparison of radiatively driven entrainment and turbulence in a smoke cloud, as simulated by different numerical models. Quart. J. Roy. Meteor. Soc., 125 , 391–423.
Brown, A. R., 1999: The sensitivity of large-eddy simulations of shallow cumulus convection to resolution and subgrid model. Quart. J. Roy. Meteor. Soc., 125 , 469–482.
Brown, A. R., S. H. Derbyshire, and P. J. Mason, 1994: Large-eddy simulation of stable atmospheric boundary layers with a revised stochastic subgrid model. Quart. J. Roy. Meteor. Soc., 120 , 1485–1512.
Carpenter, R. L., K. K. Droegemeier, and A. M. Blyth, 1998: Entrainment and detrainment in numerically simulated cumulus congestus clouds. Part III: Parcel analysis. J. Atmos. Sci., 55 , 3440–3455.
Chlond, A., 1992: Three-dimensional simulation of cloud street development during a cold air outbreak. Bound.-Layer Meteor., 58 , 161–200.
Chlond, A., . 1994: Locally modified version of Bott's advection scheme. Mon. Wea. Rev., 122 , 111–125.
Cho, H-R., 1977: Contributions of cumulus life-cycle effects to the large-scale heat and moisture budget equations. J. Atmos. Sci., 34 , 87–97.
Ciesielski, P. E., W. H. Schubert, and R. H. Johnson, 1999: Large-scale heat and moisture budgets over the ASTEX region. J. Atmos. Sci., 56 , 3241–3261.
Cuijpers, J. W. M., and P. G. Duynkerke, 1993: Large-eddy simulation of trade-wind cumulus clouds. J. Atmos. Sci., 50 , 3894–3908.
Emanuel, K. A., 1991: A scheme for representing cumulus convection in large-scale models. J. Atmos. Sci., 48 , 2313–2335.
Emanuel, K. A., and M. Zivković-Rothman, 1999: Development and evaluation of a convection scheme for use in climate models. J. Atmos. Sci., 56 , 1766–1782.
Esbensen, S., 1978: Bulk thermodynamic effects and bulk properties of small tropical cumuli. J. Atmos. Sci., 35 , 826–837.
Fraedrich, K., 1974: Dynamic and thermodynamic aspects of the parameterization of cumulus convection: Part II. J. Atmos. Sci., 31 , 1838–1849.
Grinnell, S. A., C. S. Bretherton, D. E. Stevens, and A. M. Fraser, 1996: Vertical mass flux calculations in Hawaiian trade cumulus from dual-Doppler radar. J. Atmos. Sci., 53 , 1870–1886.
Holland, J. Z., and E. M. Rasmusson, 1973: Measurements of atmospheric mass, energy, and momentum budgets over a 500-kilometer square of tropical ocean. Mon. Wea. Rev., 101 , 44–55.
Jiang, H., and W. R. Cotton, 2000: Large eddy simulation of shallow cumulus convection during BOMEX: Sensitivity to microphysics and radiation. J. Atmos. Sci., 57 , 582–594.
Johnson, R. H., and X. Lin, 1997: Episodic trade wind regimes over the western Pacific warm pool. J. Atmos. Sci., 54 , 2020–2034.
Kogan, Y. L., M. P. Khairoutdinov, D. K. Lilly, Z. N. Kogan, and Q. Liu, 1995: Modeling of stratocumulus cloud layers in a large eddy simulation model with explicit microphysics. J. Atmos. Sci., 52 , 2923–2940.
Kuo, H. L., 1965: On the formation and intensification of tropical cyclones through latent heat release by cumulus convection. J. Atmos. Sci., 22 , 40–63.
Lewellen, D. C., W. S. Lewellen, and S. Yoh, 1996: Influence of Bowen ratio on boundary-layer cloud structure. J. Atmos. Sci., 53 , 175–187.
Lin, C., and A. Arakawa, 1997a: The macroscopic entrainment processes of simulated cumulus ensemble. Part I: Entrainment source. J. Atmos. Sci., 54 , 1027–1043.
Lin, C., and A. Arakawa, . 1997b: The macroscopic entrainment processes of simulated cumulus ensemble. Part II: Testing the entraining-plume model. J. Atmos. Sci., 54 , 1044–1053.
MacVean, M. K., 1993: A numerical investigation of the criterion for cloud-top entrainment instability. J. Atmos. Sci., 50 , 2481–2495.
Moeng, C-H., 1984: A large-eddy simulation model for the study of boundary layer turbulence. J. Atmos. Sci., 41 , 2052–2062.
Moeng, C-H., . 1986: Large-eddy simulation of a stratus-topped boundary layer. Part I: Structure and budgets. J. Atmos. Sci., 43 , 2886–2900.
Morton, B. R., 1997: Discrete dry convective entities: I Review. The Physics and Parameterization of Moist Atmospheric Convection, R. K. Smith, Ed., Kluwer, 143–173.
Nitta, T., 1975: Observational determination of cloud mass flux distributions. J. Atmos. Sci., 32 , 73–91.
Nitta, T., and S. Esbensen, 1974: Heat and moisture budget analyses using BOMEX data. Mon. Wea. Rev., 102 , 17–28.
Ogura, Y., and H. R. Cho, 1973: Diagnostic determination of cumulus cloud populations from observed large-scale variables. J. Atmos. Sci., 30 , 1276–1286.
Paluch, I. R., 1979: The entrainment mechanism in Colorado cumuli. J. Atmos. Sci., 36 , 2462–2478.
Raga, G. B., J. B. Jensen, and M. B. Baker, 1990: Characteristics of cumulus band clouds off the coast of Hawaii. J. Atmos. Sci., 47 , 338–355.
Raymond, D. J., and A. M. Blyth, 1986: A stochastic model for nonprecipitating cumulus clouds. J. Atmos. Sci., 43 , 2708–2718.
Siebesma, A. P., 1998: Shallow cumulus convection. Buoyant Convection in Geophysical Flows, E. J. Plate et al., Eds., Kluwer Academic, 441–486.
Siebesma, A. P., and W. M. Cuijpers, 1995: Evaluation of parametric assumptions for shallow cumulus convection. J. Atmos. Sci., 52 , 650–666.
Siebesma, A. P., and A. A. M. Holtslag, 1996: Model impacts of entrainment and detrainment rates in shallow cumulus convection. J. Atmos. Sci., 53 , 2354–2364.
Simpson, J., and V. Wiggert, 1969: Models of precipitating cumulus towers. Mon. Wea. Rev., 97 , 471–489.
Soong, S-T., and Y. Ogura, 1976: A determination of the trade-wind cumuli population using BOMEX data and an axisymmetric cloud model. J. Atmos. Sci., 33 , 992–1007.
Squires, P., 1958: Penetrative downdraughts in cumuli. Tellus, 10 , 381–389.
Stevens, D. E., and C. S. Bretherton, 1999: Effects of resolution on the simulation of stratocumulus entrainment. Quart. J. Roy. Meteor. Soc., 125 , 425–439.
Stith, J. L., 1992: Observations of cloud-top entrainment in cumuli. J. Atmos. Sci., 49 , 1334–1347.
Taylor, G. R., and M. B. Baker, 1991: Entrainment and detrainment in cumulus clouds. J. Atmos. Sci., 48 , 112–121.
Tiedtke, M., 1989: A comprehensive mass flux scheme for cumulus parameterization in large-scale models. Mon. Wea. Rev., 117 , 1779–1800.
Yanai, M., S. Esbensen, and J. H. Cho, 1973: Determination of bulk properties of tropical cloud clusters from large-scale heat and moisture budgets. J. Atmos. Sci., 30 , 611–627.
Yau, M. K., 1979: Perturbation pressure and cumulus convection. J. Atmos. Sci., 36 , 690–694.
APPENDIX A
In-Cloud Horizontal Distributions of Thermodynamic Properties
In order to obtain fc, the mean liquid water mixing ratio, and the buoyancy in the cloud, linearized expressions are derived for the distributions of thermodynamic properties in the cloud.










APPENDIX B
Cloud Ensembles









Simplified life cycle of a shallow cumulus cloud. The cloud begins to form at time t∗ = 0 and collapses at t∗ = τ. The air in the cloud is modified by entrainment and detrainment at the top of the cloud (Ei and Di) and at the lateral boundaries (Es and Ds). Final and complete detrainment (Df) occurs at t∗ = τ
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Simplified life cycle of a shallow cumulus cloud. The cloud begins to form at time t∗ = 0 and collapses at t∗ = τ. The air in the cloud is modified by entrainment and detrainment at the top of the cloud (Ei and Di) and at the lateral boundaries (Es and Ds). Final and complete detrainment (Df) occurs at t∗ = τ
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Simplified life cycle of a shallow cumulus cloud. The cloud begins to form at time t∗ = 0 and collapses at t∗ = τ. The air in the cloud is modified by entrainment and detrainment at the top of the cloud (Ei and Di) and at the lateral boundaries (Es and Ds). Final and complete detrainment (Df) occurs at t∗ = τ
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of (a) potential temperature and (b) total water mixing ratio. The thin lines refer to the mean domain-averaged results of the 6 LES models according to the chart in (a) for a simulation period of 4 h. The squares and thick line refer to the average result for all LES models and the vertical grid used in the current simulations
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of (a) potential temperature and (b) total water mixing ratio. The thin lines refer to the mean domain-averaged results of the 6 LES models according to the chart in (a) for a simulation period of 4 h. The squares and thick line refer to the average result for all LES models and the vertical grid used in the current simulations
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean profiles of (a) potential temperature and (b) total water mixing ratio. The thin lines refer to the mean domain-averaged results of the 6 LES models according to the chart in (a) for a simulation period of 4 h. The squares and thick line refer to the average result for all LES models and the vertical grid used in the current simulations
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Cloud-top entrainment parameter α for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the LES models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of the version DECORE and circles to results of EQPROB of the parameterization as described in the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Cloud-top entrainment parameter α for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the LES models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of the version DECORE and circles to results of EQPROB of the parameterization as described in the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Cloud-top entrainment parameter α for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the LES models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of the version DECORE and circles to results of EQPROB of the parameterization as described in the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Excess flux fractions σ for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of simulation DECORE in which the organized lateral mixing is parameterized in terms of the vertical gradient of buoyancy. Diamonds refer to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Excess flux fractions σ for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of simulation DECORE in which the organized lateral mixing is parameterized in terms of the vertical gradient of buoyancy. Diamonds refer to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Excess flux fractions σ for (a) liquid water potential temperature θl, (b) total water mixing ratio rt for the results of the models KNMI, MPI, NCAR, and UOK (thin lines), according to the chart included in Fig. 2. Triangles refer to results of simulation DECORE in which the organized lateral mixing is parameterized in terms of the vertical gradient of buoyancy. Diamonds refer to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mixing fractions ν for (a) liquid water potential temperature θl, (b) total water mixing ratio rt. Thin lines refer to LES model results according to Fig. 2. Squares refer to results of simulation TOPHAT, circles to EQPROB, triangles to DECORE, and diamonds to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mixing fractions ν for (a) liquid water potential temperature θl, (b) total water mixing ratio rt. Thin lines refer to LES model results according to Fig. 2. Squares refer to results of simulation TOPHAT, circles to EQPROB, triangles to DECORE, and diamonds to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mixing fractions ν for (a) liquid water potential temperature θl, (b) total water mixing ratio rt. Thin lines refer to LES model results according to Fig. 2. Squares refer to results of simulation TOPHAT, circles to EQPROB, triangles to DECORE, and diamonds to the results of SH1996
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean liquid water content in the clouds for results of the LES models and parameterizations according to Figs. 2 and 5. Solid downward pointing triangles refer to adiabatic values
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean liquid water content in the clouds for results of the LES models and parameterizations according to Figs. 2 and 5. Solid downward pointing triangles refer to adiabatic values
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean liquid water content in the clouds for results of the LES models and parameterizations according to Figs. 2 and 5. Solid downward pointing triangles refer to adiabatic values
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean vertical velocities in the clouds for the results of the LES models and parameterizations according to Figs. 2 and 5. LES results in (a) refer to air with nonzero LWC and (b) to positively buoyant air with nonzero LWC. Results of the parameterizations are shown for the choices (a) γ = 5 and (b) γ = 0
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean vertical velocities in the clouds for the results of the LES models and parameterizations according to Figs. 2 and 5. LES results in (a) refer to air with nonzero LWC and (b) to positively buoyant air with nonzero LWC. Results of the parameterizations are shown for the choices (a) γ = 5 and (b) γ = 0
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean vertical velocities in the clouds for the results of the LES models and parameterizations according to Figs. 2 and 5. LES results in (a) refer to air with nonzero LWC and (b) to positively buoyant air with nonzero LWC. Results of the parameterizations are shown for the choices (a) γ = 5 and (b) γ = 0
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of the normalized fractional cloud cover (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of the normalized fractional cloud cover (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean profiles of the normalized fractional cloud cover (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of the normalized vertical mass flux (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of the normalized vertical mass flux (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean profiles of the normalized vertical mass flux (a) for the results of the LES models and parameterizations according to Figs. 2 and 5. Results in (b) refer to DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of normalized (a) entrainment and (b) detrainment mass fluxes for the results of the LES models and parameterizations according to Figs. 2 and 5. Normalized (c) entrainment and (d) detrainment for DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of normalized (a) entrainment and (b) detrainment mass fluxes for the results of the LES models and parameterizations according to Figs. 2 and 5. Normalized (c) entrainment and (d) detrainment for DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean profiles of normalized (a) entrainment and (b) detrainment mass fluxes for the results of the LES models and parameterizations according to Figs. 2 and 5. Normalized (c) entrainment and (d) detrainment for DECORE for γ = 5 and without cloud-top detrainment (triangles), γ = 0 and without cloud-top detrainment (thin full line), and for γ = 5 and with cloud-top detrainment (thin dotted line)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Relative contributions of cloud-top entrainment to the total entrainment mass fluxes in the cloud ensembles for DECORE (triangles) and EQPROB (circles)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Relative contributions of cloud-top entrainment to the total entrainment mass fluxes in the cloud ensembles for DECORE (triangles) and EQPROB (circles)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Relative contributions of cloud-top entrainment to the total entrainment mass fluxes in the cloud ensembles for DECORE (triangles) and EQPROB (circles)
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of (a) potential temperature and (b) total water mixing ratio for 22 and 23 Jun 1969 during BOMEX from observations at the four corners of the experimental domain (Nitta and Esbensen 1974). Also shown are the mean results of the LES models used in Figs. 2 to 12 (thick line and squares) for comparison
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Mean profiles of (a) potential temperature and (b) total water mixing ratio for 22 and 23 Jun 1969 during BOMEX from observations at the four corners of the experimental domain (Nitta and Esbensen 1974). Also shown are the mean results of the LES models used in Figs. 2 to 12 (thick line and squares) for comparison
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Mean profiles of (a) potential temperature and (b) total water mixing ratio for 22 and 23 Jun 1969 during BOMEX from observations at the four corners of the experimental domain (Nitta and Esbensen 1974). Also shown are the mean results of the LES models used in Figs. 2 to 12 (thick line and squares) for comparison
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Differences between apparent heat source and radiational heating for 22 to 23 Jun in simulations (a) DECORE, (b) EQPROB, (c) TOPHAT, and (d) for the results of SH1996. The thick line and solid downward pointing triangles refer to the results of Nitta and Esbensen (1974) during this period. Thin lines refer to the results for different launching levels of the initiating parcel of air. The height of these levels in hPa above sea level is given in the chart in (c). Averages of the simulation results at the different parcel launching levels over the 2-day period are marked by thick lines and symbols in (a) to (d). Only results above z = 500 m are shown
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Differences between apparent heat source and radiational heating for 22 to 23 Jun in simulations (a) DECORE, (b) EQPROB, (c) TOPHAT, and (d) for the results of SH1996. The thick line and solid downward pointing triangles refer to the results of Nitta and Esbensen (1974) during this period. Thin lines refer to the results for different launching levels of the initiating parcel of air. The height of these levels in hPa above sea level is given in the chart in (c). Averages of the simulation results at the different parcel launching levels over the 2-day period are marked by thick lines and symbols in (a) to (d). Only results above z = 500 m are shown
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Differences between apparent heat source and radiational heating for 22 to 23 Jun in simulations (a) DECORE, (b) EQPROB, (c) TOPHAT, and (d) for the results of SH1996. The thick line and solid downward pointing triangles refer to the results of Nitta and Esbensen (1974) during this period. Thin lines refer to the results for different launching levels of the initiating parcel of air. The height of these levels in hPa above sea level is given in the chart in (c). Averages of the simulation results at the different parcel launching levels over the 2-day period are marked by thick lines and symbols in (a) to (d). Only results above z = 500 m are shown
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

As in Fig. 13, after refitting, for simulations (a) DECORE_A, (b) TOPHAT_A, (c) DECORE_B, and (d) EQPROB_A according to the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

As in Fig. 13, after refitting, for simulations (a) DECORE_A, (b) TOPHAT_A, (c) DECORE_B, and (d) EQPROB_A according to the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
As in Fig. 13, after refitting, for simulations (a) DECORE_A, (b) TOPHAT_A, (c) DECORE_B, and (d) EQPROB_A according to the text
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

As in Fig. 14, for the apparent moisture sink
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

As in Fig. 14, for the apparent moisture sink
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
As in Fig. 14, for the apparent moisture sink
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Relative contributions of cloud-top entrainment to the total entrainment in simulation DECORE (thick line and triangles) and in DECORE_B for all days and parcel launching levels
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2

Relative contributions of cloud-top entrainment to the total entrainment in simulation DECORE (thick line and triangles) and in DECORE_B for all days and parcel launching levels
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Relative contributions of cloud-top entrainment to the total entrainment in simulation DECORE (thick line and triangles) and in DECORE_B for all days and parcel launching levels
Citation: Journal of the Atmospheric Sciences 59, 8; 10.1175/1520-0469(2002)059<1405:POTBEO>2.0.CO;2
Parameters of the numerical experiments


Available at http://www.knmi.nl/~siebesma/gcss/bomex.html.