Investigation of droplet dynamics in a convective cloud using a Lagrangian cloud model

A precipitating convective cloud is simulated successfully using the Lagrangian cloud model, in which the flow field is simulated by large eddy simulation and the droplets are treated as Lagrangian particles, and the results are analyzed to investigate precipitation initiation and to examine the parameterization of cloud microphysics. It is found that raindrops appear initially near the cloud top, in which strong turbulence and broadened droplet spectrum are induced by the entrainment of dry air, but high liquid–water mixing ratio is maintained within cloud parts because of insufficient mixing. Statistical analysis of the downward vertical velocity of a droplet W reveals that the transition from cloud droplets to raindrops occurs in the range 20 μm < r < 100 μm, while the variation of W depends on turbulence as well as the droplet radius r. The general pattern of the raindrop size distribution is found to be consistent with the Marshall–Palmer distribution. The precipitation flux can be underestimated substantially, if the terminal velocity ws\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$w_{\text{s}}$$\end{document} is used instead of W, but it is not sensitive to the choice of the critical droplet radius dividing cloud drops and raindrops. It is also found that precipitation starts earlier and becomes stronger if the effect of turbulence is included in the collection kernel.


Introduction
Simulations of clouds have been carried out traditionally by Eulerian models; either bulk models, which calculate a few numbers of bulk variables, or spectral bin models, which calculate the evolution of droplet spectra (see, e.g., Houze 1993). Although these models have been successful to simulate the cloud motion, there are many cloud processes that can be better represented by the Lagrangian cloud model (LCM), in which the flow field is simulated by large eddy simulation (LES), and the droplets are treated as Lagrangian particles; for example, the initiation of precipitation and its downward flux, the time history of individual droplets, the condensational growth of a droplet through inhomogeneous environment, etc.
The Langrangian approach has often been taken to investigate the condensational growth of a droplet (Vaillancourt et al. 2002;Lasher-Trapp et al. 2005;Lanotte et al. 2009;Sidin et al. 2009;Cooper et al. 2013) or to simulate ice particles (Sölch and Kärcher 2010). The development of LCM has not been achieved until recently, however, because of the difficulties of handling an extremely large number of droplets and the collision process, and only a few attempts have been reported so far (Andrejczuk et al. 2008(Andrejczuk et al. , 2010Shima et al. 2009;Riechelmann et al. 2012;Arabas and Shima 2013). The concept of an ensemble of 'super-droplets' is used in all LCM, where a super-droplet, the terminology coined by Shima et al. (2009), represents a large number of real droplets of a given size.
Collision schemes are developed with an approach either to represent individual collision events (Shima et al. 2009;Arabas and Shima 2013) or to create the collision process equivalent to that in the spectral bin model (Andrejczuk et al. 2010;Riechelmann et al. 2012). Shima et al. (2009) and Arabas and Shima (2013) used a Monte Carlo scheme for the collision of simulated droplets, in which collisions between randomly sampled pairs of droplets are calculated, and Andrejczuk et al. (2010) used a scheme in which collisions are assumed to occur between superdroplets within grid boxes subdivided by the droplet size, and new super-droplets are created as the outcome of collisions. On the other hand, Riechelmann et al. (2012) introduced a statistical approach, in which the growth of a super-droplet by collision/coalescence is calculated in terms of the background droplet spectrum and turbulence. The outcome is then described in terms of the modification of both the radius and the weighting factor of the superdroplet. Here the weighting factor, which represents the number of contributing real droplets to a super-droplet, differs for each super-droplet, and changes with time. Riechelmann et al. (2012) applied the LCM to the simulations of an idealized single cloud and shallow cumulus convection, under the condition corresponding to the Barbados Oceanographic and Meteorological Experiment (BO-MEX) (Holland and Rasmusson 1973). They showed that the general features of shallow cumulus clouds could be reproduced successfully in agreement with traditional bulk model results. The droplet spectra simulated by the LCM, using the collection kernels with and without the effects of turbulence (Hall 1980;Ayala et al. 2008a, b), resemble the results of a spectral bin model (Xue et al. 2008). Furthermore, the LCM results of shallow cumulus convection are in good agreement with those from the LES intercomparison under the same condition (Siebesma et al. 2003). The sensitivity tests of LCM to two important parameters of the model, the time step and the number of super-droplets, have also been carried out to investigate the convergence of the LCM solution.
In the present work, using the LCM developed by Riechelmann et al. (2012), we show that LCM can reproduce realistically the evolution of a precipitating convective cloud, including the distributions of vertical velocity, dissipation rate, liquid water mixing ratio, and droplet spectrum within a cloud at various stages of cloud development and the formation of raindrops. The results from the LCM are used to clarify precipitation initiation and to examine the parameterization of cloud microphysics such as the vertical velocity of droplets, the raindrop size distribution, and the precipitation flux. With this work we hope to show that LCM has a potential to become a powerful tool for understanding the cloud process and for developing its parameterization.

Model description
We used the LCM developed by Riechelmann et al. (2012), to which one can refer for a complete description of the model and the numerical method. Here only the basic framework of the model is discussed.
In order to handle an extremely large number of droplets in a cloud, the concept of a super-droplet is introduced. Each super-droplet represents a large number of real droplets of a given size. The number of real droplets belonging to a super-droplet of radius r n is called the 'weighting factor' A n . In the present model A n differs for each super-droplet and changes with time. For a given grid box of volume DV, the mixing ratio q l is calculated by where q 0 and q 1 are the density of air and liquid water and Nis the number of super-droplets in the grid box.
Since q 0 =q 1 \\1, the motion of each super-droplet is simulated by where V i (=dX i =dt) is the velocity of the super-droplet, X i is its position, s p is the droplet relaxation time, and u i is the fluid velocity at the droplet position. Note that, contrary to the Eulerian model, each droplet moves differently from each other and from the surrounding fluid. The inertial response time s p is calculated using the nonlinear drag law (Clift et al. 1978): In the absence of the fluid velocity (u i ¼ 0) in (2), the downward velocity of a droplet W (=ÀV 3 ) approaches the terminal velocity w s . The flow is simulated by an LES model based on the nonhydrostatic incompressible Boussinesq equation, and the equations of the fluid phase, i.e., the conservation of momentum, mass, energy, and moisture, are essentially the same as other Eulerian LES models of a cloud. The sink/ source terms through condensation/evaporation in the equations of specific humidity q and potential temperature h are directly determined by the change of q l in (1), however. The droplet growth due to the condensation/evaporation of each super-droplet is calculated using the equation suggested by Mason (1971), in which the solution and curvature effects on the droplet's equilibrium vapor pressure are neglected. The radius of all super-droplets is initially given by r = 1 lm, which corresponds to the typical size of the smallest droplet/activated cloud condensation nuclei (CCN). It also corresponds to the typical size of the smallest bin in spectral bin models (Ziegler 1985;Khairoutdinov and Kogan 2000). Droplets outside of a cloud experience evaporation, but their minimum radii are maintained at r = 0.1 lm.
In order to calculate the droplet growth by collision/ coalescence, a statistical approach is taken in which the growth of a super-droplet, arising from a large number of collision events of droplets belonging to the super-droplet, is calculated based on the background droplet spectrum and turbulence, using the collection kernel. Here the background droplet spectrum is evaluated by counting all superdroplets existing in the grid box to which the corresponding super-droplet belongs. The outcome is then described in terms of the modification of the radius (r n ) and the weighting factor (A n ) of each droplet. For this purpose, the growth of r n is calculated by applying the collection kernel using the stochastic growth equation (e.g., Rogers and Yau 1989), and the decrease of A n , representing the loss of real droplets within a super-droplet resulting from the coalescence to other bigger super-droplets, is calculated in a similar way to a spectral bin model.
The collection kernel used in the present model is equivalent to that used for a spectral bin model. The collection kernel developed by Ayala et al. (2008b) and Wang and Grabowski (2009) is used (hereafter the AW kernel), in which the effects of turbulence from relative velocity, preferential concentration, and the enhanced collision efficiency are parameterized as a function of the dissipation rate e and the Taylor-microscale Reynolds number. The Hall kernel (Hall 1980), which considers only gravitational collision, is also used for comparison. In this paper the results are from the LCM with the AW kernel, unless it is stated otherwise explicitly. Only the droplet size distribution within a grid box, not the exact locations of droplets within a grid, is relevant to the droplet growth from collision in the present model, as in a spectral model. Momentum exchange between droplets and fluid is not considered.
Results from LES and direct numerical simulation (DNS) of particle-laden turbulent flows have shown that singleparticle statistics, such as the turbulent dispersion of particles, is not significantly affected by the subgrid-scale (SGS) motions of LES, except near the wall, as they are mainly controlled by large-scale eddies, but particle-pair statistics, such as particle collision and preferential concentration, are sensitive to small-scale eddies (Yeh and Lei 1991;Uijttewaal and Oliemans 1996;Wang and Squires 1996;Armenio et al. 1999;Yamamoto et al. 2001;Fede and Simonin 2006). Therefore, in the present LCM, the transport of a droplet is calculated by following the trajectory of a Lagrangian droplet using LES, while the collision/coalescence process and preferential concentration are parameterized.

Simulation
A two-dimensional rising warm air bubble is triggered by an initial potential temperature difference h Ã given by where y c = 1,920 m and z c = 170 m mark the center of the bubble, a_{y} = 200 m and a z = 150 m the radius of the bubble, and h Ã 0 = 0.4 K the maximum temperature difference. It is homogeneous in the x-direction. The model domain is 1,280 m 9 3,840 m 9 3,840 m along the x, y, and z directions with a grid spacing of 20 m in all directions. The periodic boundary condition is applied laterally, and no slip and free slip boundary conditions are applied at the bottom and the top, respectively.
The background condition of the simulations is based on the Rain in Cumulus over the Ocean (RICO) field study (Rauber et al. 2007). The initial profiles of potential temperature and specific humidity derived from RICO are shown in Fig. 1. These are the same profiles used in the LES intercomparison of the cloud layer in RICO (van Zanten et al. 2011) except that there is no background wind in the present simulation. RICO is a comprehensive field study of shallow cumulus convection that was located in the winter trade-winds of the northwestern Atlantic Ocean, just upwind of the Islands of Antigua and Barbuda, during December 16 2004 and January 8 2005. The mean value of area-averaged precipitation during the period was about 21 W m -2 , which is just less than 1 mm day -1 .
Super-droplets are released at the beginning of the simulation and uniformly distributed all over the model domain, up to a height of 2,800 m. The average distance between super-droplets is initially 4.5 m, yielding a total number of roughly 1.5 9 10 8 and about 87 super-droplets per grid box. Using an initial weighting factor of 9 9 10 9 , the droplet concentration of approximately 100 cm -3 is represented. Riechelmann et al. (2012) showed that a maximum time step size of 0.1 s and a minimum number of about 80 initial droplets per grid box are necessary for the convergence of solutions from the sensitivity tests of the present LCM. We presume that these criteria can still be applied to the present simulation, because the roles of condensation and collision are less important during precipitation.

Evolution of cloud motion
Evolution of a convective cloud simulated by LCM is illustrated by vertical cross sections of the fluid vertical Investigation of droplet dynamics in a convective cloud velocity w, liquid water mixing ratio q l , and dissipation rate e, averaged along the x-direction at t = 10, 17, 23, and 30 min, which correspond to the stages of the initial updraft, the fully developed convective cloud, the initiation of rain, and the active rainfall, respectively (Figs. 2,3,4). Here e was calculated from the parameterization of SGS turbulence in the LES model (Riechelmann et al. 2012). Vertical profiles of the corresponding variables at the center (1,810 m \ y \ 2,010 m) are shown together, while the results from Hall and AW kernels are compared.
The typical pattern of a rising single cloud appears at t = 10 min, such as the dipole pattern of w and the increase of its size through entrainment. The buoyancy force generates a well-developed turbulent thermal characterizing a convective cloud at t = 17 min, with the maximum values of q l and e appearing near the cloud top. The values of q l and e are the largest at this stage. The initiation of precipitation at t = 23 min induces downdraft, induced by the evaporation of droplets. Most properties of a convective cloud are not affected by the difference in the collection kernels, but the value of q l is slightly smaller within the cloud and slightly larger below the cloud at t = 23 min in the AW kernel, as a result of earlier and stronger precipitation. The maximum value of e is about 30 cm 2 s -3 , which is within the range 10-100 cm 2 s -3 obtained for small cumulus clouds (Mac-Pherson and Isaac 1977;Mazin et al. 1984;Siebert et al. 2006). Figure 5 shows the evolution of mass density distribution gðln rÞ for both cases of Hall and AW kernels. Here gðln rÞ is obtained using the mass of the super-droplets inside the cloud, multiplied by the respective weighting factor. The cloud is defined as the region where q l [ 1 Â 10 À2 g/kg. Until t = 10 min, the droplet growth is essentially due to condensation, and the difference between Hall and AW kernels does not appear yet. It shows the broadening of the droplet size distribution, indicating the inhomogeneous mixing due to the entrained dry air (Brenguier and Grabowski 1993;Kogan et al. 1995;Lasher-Trapp et al. 2005). Entrainment also provides small droplets from outside of the cloud.

Precipitation initiation
On the other hand, the bimodal pattern of the droplet spectrum, attributed to the collisional growth, starts to appear at t = 17 min, with the second peak of the spectrum near r = 100 lm. At this stage the droplet growth is significantly enhanced when the effect of turbulence is included in the collection kernel (i.e., the AW kernel) in agreement with the previous studies (Xue et al. 2008;Grabowski and Wang 2009;Seifert et al. 2010;Wyszogrodzki et al. 2013). At t = 30 min large droplets with r [ 200 lm disappear from the spectrum, because they fall out of the cloud as rain.
In order to clarify the droplet growth in detail, we investigate the distribution of the droplet spectra together with the distributions of liquid-water mixing ratio (q l ) and fluid velocity (u i ) (Fig. 6). During the condensation-dominated stage at t = 10 min, the droplet spectra are unimodal everywhere. The droplet size distribution tends to be broader near the cloud top and, to a less degree, near the sides, which are the regions affected by the entrained dry air (Baker et al. 1980;Brenguier and Grabowski 1993;    -Trapp et al. 2005). At t = 17 min, the droplet spectra in the lower part of the core still show the unimodal pattern, but the bimodal pattern is clearly identified near the cloud top and, to a less degree, near the sides, suggesting the formation of raindrops due to collision effects. With the start of rain at t = 23 min, the spectral peak at the larger size in the lower part of the cloud is mainly attributed by sedimentation. As the direction of vertical velocity is reversed to downward, the droplet spectra near the cloud top become unimodal again and have a smaller mean size. As the rain progresses at t = 30 min, the bimodal pattern within the cloud tends to disappear, since most raindrops fall out of the cloud.
The formation of raindrops near the cloud top, shown in Fig. 6, can be more clearly identified by locating the appearance of raindrops directly from the LCM results. In the present paper we regard a raindrop as a droplet with r [ 40 lm. Figure 7 shows the evolution of instantaneous distributions of super-droplets in the vertical cross section with the band thickness of 100 m (500 \ x \ 600 m) during t = 15-23 min. Initially raindrops only appear near the cloud top and in the downdraft region away from the center (t = 15, 17, 19 min). Subsequently, raindrops settle gravitationally and are dispersed within a cloud by turbulent mixing, before falling out of the cloud (t = 21, 23 min). A remarkable agreement is found between the pattern of precipitation initiation shown in Fig. 7 and the observation during RICO by Small and Chuang (2008). They found that raindrops are found preferentially at the cloud top and in the downdraft region and tend to cluster with each other with the length scale about 100 m. The appearance of raindrops within 20 min in the present simulation is also consistent with the typical time scale of warm rain formation (e.g., Rogers and Yau 1989).
The critical role of entrainment and mixing in the formation of raindrops in a cumulus cloud has been proposed in many previous studies (Cooper et al. 2013;Lasher-Trapp et al. 2005;Brenguier and Grabowski 1993;Baker et al. 1980). Entrainment of unsaturated air from the environment results in the additional activation of entrained CCN, providing a source of small cloud droplets. Entrainment also generates strong temperature and moisture fluctuations near the cloud edge and thus causes droplets to experience different integral supersaturations over time. Furthermore, the enhanced turbulence, caused by evaporative cooling due to the entrained unsaturated air, helps enhance the collision efficiency. These processes promote broadening the droplet size distribution and thus produce raindrops. On the other hand, entrainment can be detrimental to raindrop formation by reducing LWC.  Investigation of droplet dynamics in a convective cloud cross section at t = 17 min (z = 1,100, 1,300 m) reveals strong fluctuation of these values, demonstrating that the entrained air remains largely unmixed with the cloud air (Fig. 8). The values of q l within cloud parts in the entrainment zone away from the center at z = 1,300 m are actually larger than those in the cloud core at z = 1,100 m, contrary to the zonal mean values shown in Fig. 3. Figure 8 also supports the fact that the values of e are highly fluctuating within a cloud up to a factor of 50 (Siebert et al. 2006).  Note that the droplet size distribution is much broader in the entrainment zone than in the cloud core before and during the raindrop formation (t = 10 and 17 min) (Fig. 6). It implies that all the conditions favorable for the formation of raindrops, such as large q l and e and the broad droplet size distribution, are found in cloud parts in the entrainment zone near the cloud top.
The distribution of super-droplets at the same horizontal cross section, over the vertical thickness of 100 m, reveals that raindrops tend to appear in the regions with strong entrainment, indicated by large velocity, and often in the regions with large q l and e (Fig. 8). Furthermore, it is interesting to observe that raindrops tend to cluster with each other, in agreement with the observation by Small and Chuang (2008). Finally, it is necessary to mention that raindrops appear about the same time at the cloud top, regardless of the collection kernel, although they are generated in a larger amount and grow faster with the AW kernel, as reflected in Fig. 5.
Contrary to the Eulerian cloud model, each droplet moves differently from each other and from the surrounding fluid, and the condensational growth and evaporation of each droplet are computed using the values of supersaturation interpolated to the position of a droplet within a grid cell. Therefore, each droplet experiences different integral supersaturations over time naturally in the LCM, and it leads to the broadening of droplet spectrum. For example, Cooper et al. (2013) and Lasher-Trapp et al. (2005) showed the broadening of droplet size distribution by calculating a microphysical parcel model that is run along trajectories produced by a standard cloud model. Meanwhile, in the entrainment zone, turbulence is generated by evaporative cooling, thus enhancing the collision efficiency. On the other hand, the cloud air remains largely unmixed with the entrained air in the entrainment zone, thus preventing q l from decreasing substantially. b Fig. 8 Distributions of vertical velocity w, liquid water mixing ratio q l , and dissipation rate e at the horizontal cross section at t = 17 min. At bottom are the distributions of super-droplets within the vertical thickness 100 m: a z = 1,100 m, b z = 1,300 m

Examination of microphysics parameterization
In previous sections the LCM is shown to reproduce realistically the evolution of a precipitating cloud. Meanwhile, LCM can provide us the information on the motion of individual droplets directly, contrary to bulk and spectral bin models. Most bulk models are developed based on the assumptions of the bimodal distribution of droplets, dividing into cloud droplets and raindrops according to the size (Kessler 1969;Emanuel 1994;Stensrud 2007). Cloud drops, whose radius is smaller than a critical radius r c , are assumed to follow the fluid velocity, while the settling velocity is neglected. On the other hand, raindrops, whose radius is larger than r c , are assumed to fall with the terminal velocity while being advected by the fluid velocity; i.e., or W ¼ w s if w j j\\w s . Note that (5) is obtained, when s p is much smaller than the time scale of a turbulent flow in (2). In Sect. 4.2, r c = 40 lm is used, but a wide range of values is used in bulk models (r c = 40-100 lm) (Kessler 1969;Shiino 1983;Beheng and Doms 1986;Emanuel 1994;Seifert and Beheng 2001). Table 1 Values of the rms of ðW À w s Þ=w s (=r W=ws ) in the inner (1,800 m B y \ 2,000 m, at t = 17 min, and 1,600 m B y \ 2,200 m, at t = 23 min) and outer region of cloud droplets for r = 10, 40, 100 lm (t = 17, 23 min) Here the thickness of the spectral band for the calculation is given proportional to r, i.e., Dr ¼ 0:1r (=1, 4, 10 lm) Bulk models are then made of equations to calculate the mixing ratios of cloud water and rainwater, q C and q R (q l ¼ q C þ q R ). For example, as the simplest representation of warm cloud microphysics, Kessler (1969) suggested a bulk model as where U represents the contribution from condensation, P the precipitation flux, A the autoconversion from cloud water to rainwater, C the collection of cloud water by rainwater, E evaporation, and D C and D R the turbulent diffusion of cloud water and rainwater, respectively. In many bulk models, including Kessler (1969), the raindrop size distribution is assumed to follow the Marshall-Palmer distribution (Marshall and Palmer 1948) as where n 0 and k are constants.
It is important to recognize that LCM can provide information for the parameterization of most terms in the RHS of (6) and (7) in principle. For example, P, D C and D R can be evaluated by analyzing the transport of droplets, and U, C and E can be evaluated by analyzing the corresponding contributions to droplets. Meanwhile, A can be evaluated by analyzing the evolution of droplet spectrum as in spectral bin models.
With an aim to show the applicability of the LCM in the present work, however, we examine the following three hypotheses used in the Kessler's model (1969) by analyzing the LCM results. First, as to how the vertical velocity of a droplet varies with r during the transition from a cloud drop to a raindrop. Second, how the raindrop size is dis-tributed and whether the Marshall-Palmer distribution (8) is observed. Finally, how accurately the precipitation flux based on w s predicts the real one based on W. Figure 10 shows the scatter plots of W=w s vs. r at t = 17 and 23 min. The plots are made separately for the inner region (1,800 m B y \ 2,000 m, at t = 17 min, and 1,600 m B y \ 2,200 m, at t = 23 min) and for the outer region. The asymmetry at t = 23 min is due to the presence of strong downdraft. The range of W=w s , which is very large for small droplets (r \ 20 lm), decreases rapidly with r, converging to W=w s ¼ 1 (r [ 100 lm). This implies that W is basically determined by the fluid velocity w, independent of w s , for droplets with r \ 20 lm, but W is mainly controlled by w s for droplets with r [ 100 lm, although a large variation of W comparable to w s is still present. One can regard the former (r \ 20 lm) as cloud droplets (W ¼ Àw) and the latter (r [ 100 lm) as raindrops (W ¼ Àw þ w s ), according to Fig. 10. The transition from cloud droplets to raindrops occurs for smaller r in the outer region, in which the fluid velocities, both mean and turbulent, are smaller. Further evidence is provided by the evaluation of the root mean square (rms) values of ðW À w s Þ=w s (=r W=w s ) in the inner and outer regions, as shown in Table 1. Table 1 also reveals that the rms of W is comparable to w s , even for r = 100 lm.
We can also examine how accurately (5) represents the real droplet velocity W. The relative frequency distribution of DW ¼ w þ W ð¼ w À V 3 Þ for droplets with different sizes r = 10, 40, and 100 lm reveals that the variance of DW À w s , r DW , is not so sensitive to r, although it tends to increase weakly with r, suggesting the effect of particle inertia ( Fig. 11; Table 2). On the other hand, r DW =w s decreases with r. These suggest that the deviation of W from w À w s is small in raindrops (r = 100 lm), but not negligible. Although the transition from a cloud droplet (W ¼ Àw) to raindrops (W ¼ Àw þ w s ) with increasing r is expected from (5), no information has been available so far on how W varies with r (Fig. 10) and what the probability distribution of W is ( Fig. 11; Tables 1 and 2). The present results show that they are affected by the turbulent flow field in a cloud as well as r. Figure 12 reveals that the raindrop size distribution from LCM decreases exponentially with r, as in the given by (8). Here nðrÞ is calculated by the average over the region with rain for droplets with r [ 40 lm at t = 23 min. The variations of nðrÞ with r using the empirical constants, n 0 = 8 9 10 -3 m -3 mm -1 and k ¼ 4:1P À0:21 R mm -1 from Marshall and Palmer (1948) Fig. 13 The precipitation flux profile at t = 24, 26, 28 min, using the actual fall velocity W (left) and the terminal velocity w s (right): a r c = 40 lm, b r c = 100 lm Investigation of droplet dynamics in a convective cloud with P R = 0.1 mm h -1 are also shown. It is important, however, to notice that the direct quantitative comparison of nðrÞ is not appropriate here, because only the precipitation from a single cloud is considered. It may be also necessary to mention that the formula (8) is obtained mainly from the measurements of larger raindrops by Marshall and Palmer (1948) and that the breakup of droplets is not included in the present model. Nonetheless, it is worthwhile to mention that Fig. 12 is the first attempt to obtain the raindrop size distribution directly from simulation, taking advantage of LCM. The decrease of nðrÞ with r tends to be slightly slower in the case of the AW kernel, as expected from the higher precipitation rate.
In the LCM, P can be calculated directly using the downward velocities of individual raindrops (r [ r c ) as where nðrÞ is the number of raindrops of size r per unit volume. The precipitation rate P R is then calculated by the average of P over the thickness of 50 m from the surface. In bulk models w s is usually used, instead of W, to calculate P in (9). Therefore, it is interesting to investigate how different the value of P obtained in this way is from the real one calculated directly from (9). Figure 13 shows that, if w s is used for (9), instead of W, it becomes substantially smaller than the actual values of P. It indicates that the downdraft associated with precipitation, as shown in Figs. 3 and 6, enhances the sedimentation velocity itself. Note that precipitation initiates downdraft by the evaporative cooling of droplets and the downward dragging by settling droplets, although the latter is not considered in the present simulation (e.g., Houze 1993).
Comparison is also made for the profiles of P from two different values of the critical radius r c (r c = 40 and 100 lm) to investigate its sensitivity. It is found, however, that the critical value of r c does not affect P significantly (Fig. 13). It may be associated with the fact that the droplet mass density is lower and W is smaller for droplets in the range r = 40-100 lm than for larger droplets (r [ 100 lm). The use of w s instead of W can lead to the underestimation of the peak precipitation rate by as much as 28 %.
Finally, it is examined how the inclusion of the turbulence effect in the collision kernel affects precipitation by comparing the profiles of P (Fig. 14) and the time series of precipitation rate P R (Fig. 15). It is found that precipitation not only starts earlier but also becomes stronger under the AW kernel, which is consistent with the LES results with the bulk microphysics (Seifert et al. 2010) and with the spectral bin microphysics (Wyszogrodzki et al. 2013). Figure 14 also shows that raindrops fall faster and in a larger amount in the AW scheme than in the Hall scheme, leading to earlier and stronger precipitation, although they appear about the same time in both schemes.

Conclusion
In the present work, the physical process influencing the microphysical property in a precipitating convective cloud is investigated using the Lagrangian cloud model (LCM), in which the flow field is simulated by LES, and the droplets are treated as Lagrangian particles. For the LCM, the concept of an ensemble of super-droplets, with each super-droplet representing a large number of real droplets of the same size, is introduced, and a statistical method is developed to calculate the growth of a super-droplet by collision/coalescence based on the background droplet spectrum and turbulence, using the collection kernel.
The LCM is shown to reproduce the general features of the cloud and precipitation process successfully, including the evolutions of vertical velocity, liquid-water mixing ratio, dissipation rate, and droplet spectrum. It is then applied to investigate precipitation initiation and to examine the parameterization of cloud microphysics.
It is found that raindrops appear initially near the cloud top and in the downdraft region and tend to cluster with each other. The pattern of simulated precipitation initiation is in good agreement with the observation during RICO by Small and Chuang (2008). The mechanism for the raindrop formation is explained by the fact that, in the entrainment zone near the cloud top, strong turbulence and broadened droplet spectrum are induced by the entrainment of dry air, but high liquid-water mixing ratio is maintained within cloud parts because of insufficient mixing.
Statistical analysis of the vertical velocity of a droplet W reveals that the transition from cloud droplets to raindrops occurs in the range 20 lm \ r \ 100 lm, while the variation of W depends on turbulence as well as the droplet radius r. The general pattern of the raindrop size distribution is consistent with the Marshall-Palmer distribution. The precipitation flux can be underestimated substantially, if the terminal velocity w s is used instead of W, but it is not sensitive to the choice of the critical droplet radius dividing cloud drops and raindrops. It is found that precipitation starts earlier and becomes stronger if the effect of turbulence is included in the collection kernel.
The present work demonstrates that the LCM, which is a more natural approach to simulate particle-laden turbulent flows, can be used to reproduce the realistic dynamic feature of clouds, including precipitation initiation, while providing information of individual droplets that could not be obtained with other models. It thus illustrates that the LCM has a potential to become a powerful tool for understanding the cloud processes and for developing its parameterizations.
Nonetheless, further improvement of the model in various aspects is necessary to simulate more realistic clouds in future. For example, the collision/coalescence scheme and the collection kernel should be further elaborated, possibly by using information from DNS with higher Reynolds numbers, or by comparing with the spectral bin model results. More sophisticated two-phase fluid dynamical processes can be taken into account, such as the effect of the SGS turbulence on the transport of a super-droplet and the exchange of momentum between particle and fluid phases. It is also necessary to investigate further the numerical convergence with increasing resolution and particle numbers under various schemes of SGS turbulence and collision. Finally, there are various other physical processes that must be included for more realistic simulations, such as the activation of CCN and breakup of large droplets.