Large Eddy Simulation of Microphysics and Inﬂuencing Factors in Shallow Convective Clouds

: A ﬂight of shallow convective clouds during the SCMS95 (Small Cumulus Microphysics Study 1995) observation project is simulated by the large eddy simulation (LES) version of the Weather Research and Forecasting Model (WRF-LES) with spectral bin microphysics (SBM). This study focuses on relative dispersion of cloud droplet size distributions, since its inﬂuencing factors are still unclear. After validation of the simulation by aircraft observations, the factors affecting relative dispersion are analyzed. It is found that the relationships between relative dispersion and vertical velocity, and between relative dispersion and adiabatic fraction are both negative. Furthermore, the negative relationships are relatively weak near the cloud base, strengthen with the increasing height ﬁrst and then weaken again, which is related to the interplays among activation, condensation and evaporation for different vertical velocity and entrainment conditions. The results will be helpful to improve parameterizations related to relative dispersion (e.g., autoconversion and effective radius) in large-scale models.


Introduction
Shallow convective clouds are widely distributed in every corner of the Earth from the tropics to the polar regions, and from the oceans to the continents. Shallow convective clouds play a vital role in the atmospheric dynamic and thermodynamic processes, affecting radiation, precipitation, et cetera [1][2][3][4][5][6]. Therefore, shallow convective clouds have always been an important part of cumulus parameterization in both weather forecast models and general circulation models [7][8][9][10][11].
To study the microphysical properties of warm shallow convective clouds, cloud droplet size distributions are particularly important. Relative dispersion of cloud droplet size distribution, the ratio of standard deviation to mean radius, is a key parameter, affecting parameterizations of effective radius [17][18][19][20][21][22][23], and cloud-rain autoconversion in models [24][25][26]. Understanding the variation of relative dispersion in cloud is critical to improve the microphysical process parameterizations in numerical models [27]. Prabha et al. [28] found that it was reasonable to assume that constant standard deviation for premonsoon In the case of RF12, the small cumulus clouds form in the afternoon. These small cumulus clouds appear in the southeast of the area and move to the northwest. These clouds pass through the bay and the peninsula, then eventually move to the mainland. There is no effective surface precipitation. The temperature in clouds is always higher than 0 °C. Generally speaking, these clouds can be considered as warm-phase and nonprecipitation clouds. In the case of RF12, the small cumulus clouds form in the afternoon. These small cumulus clouds appear in the southeast of the area and move to the northwest. These clouds pass through the bay and the peninsula, then eventually move to the mainland. There is no effective surface precipitation. The temperature in clouds is always higher than 0 • C. Generally speaking, these clouds can be considered as warm-phase and nonprecipitation clouds.

Model Setup: WRF-LES and SBM Scheme
Due to the small cross-sectional area of shallow convective clouds, a high-resolution numerical model is necessary to simulate the detailed structure of shallow convective clouds. For a grid resolution of 500 m, the common shallow convection covers only one or only a few grid points [37]. Therefore, large eddy simulation has become one of the best ways to simulate shallow convection [38,39]. Compared with other large-scale systems, the life cycle and scale of shallow convective clouds are shorter. Therefore, for a fixed simulation area, the entire shallow convective life cycle from formation to dissipation can be included, and detailed simulations can be obtained with large eddy simulations [40][41][42].
The Weather Research and Forecasting Model (WRF, version 3.9.1) with the built-in large eddy simulation (LES, [43]) module is used here. LES can finely analyze the energy transportation of the turbulent motion processes, and it has been widely used for turbulence research and physical process research. For the cloud microphysics scheme, single-moment or double-moment bulk parameterization schemes can capture the main features of cloud droplet size distributions, but with many assumptions [44][45][46]. The bin microphysics scheme is a better choice [47], which has been successfully applied to simulate the radiation characteristics of warm clouds and precipitation and has played a good role in simulating ice-phase clouds [48][49][50]. The fast version of spectral bin microphysics parameterization scheme (Fast-SBM, [51,52]) is used. Full-SBM has the species of cloud, rain, ice(columns), ice(plates), ice(dendrites), snow, graupel, and hail. Compared with Full-SBM, Fast-SBM combines ice(columns), ice(plates) and ice(dendrites) to ice and combines graupel and hail to graupel ( [53]). Because we focus on liquid-phase clouds, the difference between Full-SBM and Fast-SBM is expected to have little effect on the results in this study. Previous studies also found that the results of Fast-SBM were similar to those of Full-SBM ( [51]). For cloud droplets, spectral bin microphysics (SBM) includes 33 bins with the mass of one bin being twice the previous one. For the cloud droplet bins, the mass of the initial bin is 3.35 × 10 −14 kg.
A 5-level nested grid area (D01, D02, D03, D04, D05, Figure 2a) is used in the simulations. LES is active on D04 and D05, which means that the boundary conditions from the real situation can be passed to D05 through D01-D04. The grid resolutions of D01-D03 are 12.5 km × 12.5 km, 2.5 km × 2.5 km and 0.5 km × 0.5 km, respectively. For the regions where large eddy simulation is active, the resolutions are D04~∆x = ∆y = 100 m and D05~∆x = ∆y = 33.3 m. Corresponding to the high horizontal resolution of 33.3 m, we are able to set a high vertical resolution. There are a total of 122 vertical levels, and the layers are dense below 5 km where the shallow convection exists. With the high resolution, the model can capture the microphysical, dynamical and thermodynamic processes well in shallow cumulus clouds. Similar resolutions have been used in previous studies [54,55]. For D04 and D05, the areas of the domains are 25.6 km × 25.6 km and 8.5 km × 8.5 km, respectively. The cumulus convection parameterization scheme is only enabled in D01. The boundary layer parameterization scheme is turned off in D04 and D05, because it is the replacement of large eddy simulation. The large-scale forcing data come from the Climate Forecast System Reanalysis (CFSR) provided by the National Centers for Environmental Prediction (NCEP), with a resolution of 0.5 • × 0.5 • . Lateral boundary conditions are updated every 6 h. The simulation time is from 1700 (UTC) to 2300 (UTC) and the first 3 h are taken as the spin-up time. We use the data from 2000 (UTC) to 2300 (UTC) for analysis. The model outputs the results every 15 s. As shown later in Figure 3, liquid water potential temperature, water vapor mixing ratio and cloud water mixing ratio from the simulation are consistent with observations. Therefore, we believe that it is enough to discard the first 3 h as spin-up time. At the beginning of the simulation, the relative humidity at the low level in the simulation domain (D05) is high, the surface temperature is about 29 • C and the lifting condensation level (LCL) is about 922 hPa (Figure 2b). The lower atmosphere is dominated by the southeast wind, which turns into the northeast wind above 500 hPa. Below 700 hPa, the atmosphere is humid. The cloud condensation nuclei (CCN) concentration is calculated by CCN = AS B , where A = 500 cm −3 , B = 0.5 and S is supersaturation. At the beginning of the simulation, the domain D05's average sensible heat flux is 117.96 W·m −2 and the average latent heat flux is 180.02 W·m −2 .
The boundary layer parameterization scheme is turned off in D04 and D05, because it is the replacement of large eddy simulation. The large-scale forcing data come from the Climate Forecast System Reanalysis (CFSR) provided by the National Centers for Environmental Prediction (NCEP), with a resolution of 0.5° × 0.5°. Lateral boundary conditions are updated every 6 h. The simulation time is from 1700 (UTC) to 2300 (UTC) and the first 3 h are taken as the spin-up time. We use the data from 2000 (UTC) to 2300 (UTC) for analysis. The model outputs the results every 15 s. As shown later in Figure 3, liquid water potential temperature, water vapor mixing ratio and cloud water mixing ratio from the simulation are consistent with observations. Therefore, we believe that it is enough to discard the first 3 h as spin-up time. At the beginning of the simulation, the relative humidity at the low level in the simulation domain (D05) is high, the surface temperature is about 29 °C and the lifting condensation level (LCL) is about 922 hPa (Figure 2b). The lower atmosphere is dominated by the southeast wind, which turns into the northeast wind above 500 hPa. Below 700 hPa, the atmosphere is humid. The cloud condensation nuclei (CCN) concentration is calculated by CCN = , where A = 500 cm −3 , B = 0.5 and S is supersaturation. At the beginning of the simulation, the domain D05's average sensible heat flux is 117.96 W·m −2 and the average latent heat flux is 180.02 W·m −2 .

Methods for Calculating Cloud Physical and Microphysical Properties
As mentioned above, the shallow convective clouds in this case are all warm-phase, and cloud droplets can be assumed to be spherical particles. The mass bins obtained by the SBM can be converted into radius bins, and the radius ( ) of the k-th mass bin ( )

Methods for Calculating Cloud Physical and Microphysical Properties
As mentioned above, the shallow convective clouds in this case are all warm-phase, and cloud droplets can be assumed to be spherical particles. The mass bins obtained by the SBM can be converted into radius bins, and the radius (r k ) of the k-th mass bin M(k) is: where ρ w is the density of liquid water and the initial value of k is 1. Because the mass of one bin is twice the previous one, Equation (1) indicates that the radius increases at the rate of 3 √ 2, and r k can be expressed as: The bin width (∆r k ) corresponding to the center radius of each bins is: The number density n(r k ) of each bin is: where ρ d is dry air density under standard atmospheric conditions, q c (k) is the cloud water mixing ratio per unit mass of air in the k-th mass bin. Total number concentration (N), average radius (r m ), standard deviation of cloud droplet size distribution (σ) and the relative dispersion of cloud droplet size distribution (ε) can be obtained by the following equations: respectively. Buoyancy (B), virtual temperature (T v ), cloud liquid water potential temperature (θ l , [56]) and latent heat (L v ) are calculated by the following formulas: Atmosphere 2021, 12, 485 6 of 21 Among them, the subscripts c and e represent cloud and environment, respectively; q l is the liquid water mixing ratio, equal to cloud water mixing ratio here, since there is no precipitation; q v is the water vapor mixing ratio; c p is the dry air specific heat capacity at constant pressure; ∆m is the mass variation at grid points only affected by condensation or evaporation within the time step and l v is the specific latent heat coefficient.
Cloud adiabatic fraction is obtained by the following method: where q ca is the adiabatic cloud water mixing ratio in the cloud.

Method for Normalizing Variables
The following method is used to normalize the physical variable x: where x i is the data needing to be normalized, x is the mean value of the data, s is the standard deviation of the data and the normalized data is x s . The normalization can transform variables to non-dimensional variables, and can be used to fairly compare two different variables.

Comparison of Simulations and Observations
Before comparing simulations with observations, clouds need to be selected. Consider the grids with q c > 0.01 g/kg as cloud and those with q c < 0.01 g/kg as environment. Figure 3 shows the comparison of liquid water potential temperature, water vapor mixing ratio and cloud water mixing ratio in the simulation and observation. The observational data are from the plots in Neggers et al. [36]. The liquid water potential temperature is a conserved variable in the reversible moist adiabatic motion [57], which is very convenient to analyze the wet process of the air mass during the lifting and condensation process. Whether in the cloud or in the environment, the results of the liquid water potential temperature simulated by the LES are very close to the results observed by aircraft. The simulation and observation of the water vapor mixing ratio and cloud water mixing ratio are also in good agreement. Therefore, the LES results are thought to be reliable. Figure 4 shows cloud fraction at different heights from the simulation. The cloud fraction reaches its maximum value at 1.0 km, and then gradually decreases as the height increases. Most clouds have their cloud tops below 3.8 km. Since the lowest cloud base is around 0.6 km, the average thickness of most clouds does not exceed 3.2 km. There are also some clouds with their tops reaching 6 km or even higher. These clouds only occupy a small part and they are not analyzed here, to focus on shallow cumulus clouds. Based on aircraft observations, it is hard to determine if there were deep clouds in the field campaign, because aircraft observations only flew along a line and could not provide the whole picture of the clouds simulated here. In the simulated domain, clouds in different locations may have different cloud base heights. If the lowest cloud base height is taken to be the cloud base height for the whole domain, there would be big biases to plot vertical profiles of properties in cloud (method 1). Instead, identifying the cloud base for each cloud would cause huge computational cost (method 2). Therefore, a method different from the previous two methods is chosen. D05 is divided into 81 (9 × 9) sub-areas with each sub-area of approximately 950 m × 950 m, equivalent to the shallow cumulus cloud size. A cloud base is determined for each sub-area, where the lowest cloud grid has a cloud water mixing ratio exceeding 0.01 g/kg. In the following analysis, sub-areas of 950 m × 950 m are selected with the criteria of cloud top heights being below 3.8 km, cloud thicknesses smaller than 3.2 km and cloud base heights lower than 1.0 km. Cloud thickness is the difference between cloud top and base heights with cloud water mixing ratio greater than 0.01 g/kg.

Cloud Physics Properties
Vertical profiles of properties above cloud base are needed in the following analysis. As mentioned above, D05 is divided into 81 sub-areas with each sub-area of approximately 950 m × 950 m. Vertical profiles of properties above cloud base are calculated by combining the properties in each sub-area with respect to its own cloud base.
Vertical velocity ( ), and in the cloud are three important cloud physics properties affecting cloud development ( Figure 5). The positive vertical velocity within the cloud gradually increases from the cloud bases and reaches the maximum value near 1.0 km above the cloud bases; then the positive vertical velocity remains stable and decreases slightly with the increasing height. Instead, the absolute value of negative vertical velocity gradually increases. The mean vertical velocity changes from positive to negative near 1.8 km. The buoyancy has a similar trend with the vertical velocity, but the height at which the mean buoyancy changes from positive to negative is ~1.4 km, lower than 1.8 km for the vertical velocity, which is consistent with Newton's law of motion. From the point of view of latent heat, the mean latent heat in the cloud is positive within about 1.0 km above the cloud bases, indicating that the lower layers are dominated by condensation. Above about 1.0 km, the mean latent heat is negative, which could be due to evaporation during entrainment.

Cloud Physics Properties
Vertical profiles of properties above cloud base are needed in the following analysis. As mentioned above, D05 is divided into 81 sub-areas with each sub-area of approximately 950 m × 950 m. Vertical profiles of properties above cloud base are calculated by combining the properties in each sub-area with respect to its own cloud base.
Vertical velocity (w), B and L v in the cloud are three important cloud physics properties affecting cloud development ( Figure 5). The positive vertical velocity within the cloud gradually increases from the cloud bases and reaches the maximum value near 1.0 km above the cloud bases; then the positive vertical velocity remains stable and decreases slightly with the increasing height. Instead, the absolute value of negative vertical velocity gradually increases. The mean vertical velocity changes from positive to negative near 1.8 km. The buoyancy has a similar trend with the vertical velocity, but the height at which the mean buoyancy changes from positive to negative is~1.4 km, lower than 1.8 km for the vertical velocity, which is consistent with Newton's law of motion. From the point of view of latent heat, the mean latent heat in the cloud is positive within about 1.0 km above the cloud bases, indicating that the lower layers are dominated by condensation. Above about 1.0 km, the mean latent heat is negative, which could be due to evaporation during entrainment. The profiles of vertical velocity, buoyancy and latent heat can explain why most clouds maintain shallow convection: the buoyancy within the cloud is suppressed and the vertical velocity is weak. The formation of shallow convection can also be explained by the large-scale forcing. The clouds form in the west of the subtropical anticyclone. Figures  6 and 7 are both basic information of cloud properties. We would like to present them before introducing mechanisms in the following subsections. Figure 6 shows the cloud water mixing ratio and the cloud adiabatic fraction. The cloud water mixing ratio increases significantly from the cloud base to 1.0 km, and then increases gradually above 1.0 km; the maximum cloud water mixing ratio is achieved at ~3.2 km. The cloud base has the highest adiabatic fraction. As the height increases, the adiabatic fraction gradually decreases. Below 1.4 km, the rate of adiabatic fraction decline is the fastest; above 1.4 km, the adiabatic fraction decreases slowly and stabilizes at around 0.2. Corresponding to cloud water mixing ratio and cloud adiabatic fraction, the vertical profiles of cloud droplet spectral properties are also analyzed (Figure 7), including relative dispersion, standard deviation, number concentration and mean radius. gradually decreases with the increasing height within 0.2 km above the cloud base, and reaches the minimum mean value (~0.3) at 0.2 km. further increases until 1.0 km and then is almost constant. increases steadily with height, and the minimum is at the cloud base (1.14 μm). fluctuates slightly near the cloud base, and the maximum occurs at 0.1 km above the cloud base (206.1 cm −3 ). N then gradually decreases as the height increases.
is similar to , and the minimum is at the cloud base (3.21 μm). The profiles of vertical velocity, buoyancy and latent heat can explain why most clouds maintain shallow convection: the buoyancy within the cloud is suppressed and the vertical velocity is weak. The formation of shallow convection can also be explained by the largescale forcing. The clouds form in the west of the subtropical anticyclone. Figures 6 and 7 are both basic information of cloud properties. We would like to present them before introducing mechanisms in the following subsections. Figure 6 shows the cloud water mixing ratio and the cloud adiabatic fraction. The cloud water mixing ratio increases significantly from the cloud base to 1.0 km, and then increases gradually above 1.0 km; the maximum cloud water mixing ratio is achieved at~3.2 km. The cloud base has the highest adiabatic fraction. As the height increases, the adiabatic fraction gradually decreases. Below 1.4 km, the rate of adiabatic fraction decline is the fastest; above 1.4 km, the adiabatic fraction decreases slowly and stabilizes at around 0.2. Corresponding to cloud water mixing ratio and cloud adiabatic fraction, the vertical profiles of cloud droplet spectral properties are also analyzed (Figure 7), including relative dispersion, standard deviation, number concentration and mean radius. ε gradually decreases with the increasing height within 0.2 km above the cloud base, and reaches the minimum mean value (~0.3) at 0.2 km. ε further increases until 1.0 km and then is almost constant. σ increases steadily with height, and the minimum is at the cloud base (1.14 µm). N fluctuates slightly near the cloud base, and the maximum occurs at 0.1 km above the cloud base (206.1 cm −3 ). N then gradually decreases as the height increases. r m is similar to σ, and the minimum is at the cloud base (3.21 µm).

Effects of Standard Deviation and Mean Radius on Relative Dispersion
Since ε is the ratio of σ to r m , it is critical to analyze the variations of σ and r m . Figure 7 shows that, the variation of ε with height is not monotonic, different from the variations of σ and r m that increase with the increasing height. Since σ and r m are both dimensional variables, they need to be transformed into non-dimensional variables with the normalized method in Section 3.2 (Figure 8a,b). The variability of the two variables within each 50 m height is analyzed (Figure 8c,d). If we take the derivative of Equation (8)  From Equation (15), the sign of is determined by a cross multiplication difference: In Equation (16), the influencing factors of the variation of relative dispersion with height can be divided into two categories, namely, the change rate of the variable with height ( ) and the value of the variable at a certain height. Since each term is of the same form, it is possible to estimate the positive or negative difference of the multiplication by simply subtracting the factors of the multiplication. Figure 9 is the normalized vertical profiles of − and − . However, it should be noted that as the height increases, the values of and become larger, and the above method may not work well. Therefore, only the data with the largest changes below 2.0 km are analyzed. It can be seen from Figure 9 that from the cloud base to 0.2 km, both − and − are mainly negative, which together determines the decreasing trend of relative dispersion with the increasing height (Figure 7a). In the range of 0.2-1 km, the increasing trend of relative dispersion with height mainly comes from positive − . In the range of 1-2 km, both − and − are near 0, so the relative dispersion changes slowly with height.  From Equation (15), the sign of dε dz is determined by a cross multiplication difference: In Equation (16), the influencing factors of the variation of relative dispersion with height can be divided into two categories, namely, the change rate of the variable with height ( d dz ) and the value of the variable at a certain height. Since each term is of the same form, it is possible to estimate the positive or negative difference of the multiplication by simply subtracting the factors of the multiplication. Figure 9 is the normalized vertical profiles of r m − σ and dσ dz − dr m dz . However, it should be noted that as the height increases, the values of σ and r m become larger, and the above method may not work well. Therefore, only the data with the largest ε changes below 2.0 km are analyzed. It can be seen from Figure 9 that from the cloud base to 0.2 km, both r m − σ and dσ dz − dr m dz are mainly negative, which together determines the decreasing trend of relative dispersion with the increasing height (Figure 7a). In the range of 0.2-1 km, the increasing trend of relative dispersion with height mainly comes from positive r m − σ. In the range of 1-2 km, both r m − σ and dσ dz − dr m dz are near 0, so the relative dispersion changes slowly with height. Besides and , it is very important to analyze cloud droplet size distributions, since represents the width of cloud droplet size distribution. Figure 10 illustrates that, cloud droplets of 2-5 μm gradually decrease with the increasing height below 0.2 km, and cloud droplets of 5-10 μm begin to appear and increase rapidly. This is in line with the characteristics of the cloud base itself. In the cloud base, a large number of cloud condensation nuclei are activated into small cloud droplets, and the condensation rate is large for the small droplets according to the adiabatic condensation theory [58]; as a result, the small droplets grow rapidly. This explains why the d d ⁄ within 0.2 km is large ( Figure   8d). Because of the fast increase in , decreases with the increasing height. The increase in above 0.2 km is related to entrainment-mixing processes. The radius of a large number of cloud droplets are in the range of 10-15 μm, but collision-coalescence is still thought to be weak. However, because of the entrainment of dry air either through the lateral boundary, cloud top or both, some droplets evaporate into small droplets, which leads to the increase of . One evidence is that cloud droplets less than 10 μm always exist, and do not significantly reduce with height. Besides r m and σ, it is very important to analyze cloud droplet size distributions, since ε represents the width of cloud droplet size distribution. Figure 10 illustrates that, cloud droplets of 2-5 µm gradually decrease with the increasing height below 0.2 km, and cloud droplets of 5-10 µm begin to appear and increase rapidly. This is in line with the characteristics of the cloud base itself. In the cloud base, a large number of cloud condensation nuclei are activated into small cloud droplets, and the condensation rate is large for the small droplets according to the adiabatic condensation theory [58]; as a result, the small droplets grow rapidly. This explains why the dr m /dz within 0.2 km is large (Figure 8d). Because of the fast increase in r m , ε decreases with the increasing height. The increase in ε above 0.2 km is related to entrainment-mixing processes. The radius of a large number of cloud droplets are in the range of 10-15 µm, but collision-coalescence is still thought to be weak. However, because of the entrainment of dry air either through the lateral boundary, cloud top or both, some droplets evaporate into small droplets, which leads to the increase of σ. One evidence is that cloud droplets less than 10 µm always exist, and do not significantly reduce with height.

Effects of Vertical Velocity and Entrainment on Relative Dispersion
In addition to r m and σ, which can directly affect ε, w and AF in the cloud can affect ε through physical processes [20,22,58]. Figure 11 shows a negative correlation between the relative dispersion and vertical velocity, and the correlation coefficient is −0.54. The negative correlation is consistent with previous studies [59,60]. On the one hand, larger w increases supersaturation, and further enhances condensational growth [61]. Larger supersaturation corresponds to narrower cloud droplet size distributions, in other words, smaller ε [17]. On the other hand, for larger w, less dry air is entrained into cloud [58,62,63] and droplet evaporation is weaker. As a result, broadening of cloud droplet size distribution during entrainment-mixing is weakened, in other words, smaller ε. Oppositely, when w is smaller or even negative, entrainment becomes stronger. Some droplets could evaporate faster than others and cloud droplet size distributions broaden significantly, in other words, larger ε [64].

Effects of Vertical Velocity and Entrainment on Relative Dispersion
In addition to and , which can directly affect , and in the cloud can affect through physical processes [20,22,58]. Figure 11 shows a negative correlation between the relative dispersion and vertical velocity, and the correlation coefficient is −0.54. The negative correlation is consistent with previous studies [59,60]. On the one hand, larger w increases supersaturation, and further enhances condensational growth [61]. Larger supersaturation corresponds to narrower cloud droplet size distributions, in other words, smaller [17]. On the other hand, for larger w, less dry air is entrained into cloud [58,62,63] and droplet evaporation is weaker. As a result, broadening of cloud droplet size distribution during entrainment-mixing is weakened, in other words, smaller . Oppositely, when w is smaller or even negative, entrainment becomes stronger. Some droplets could evaporate faster than others and cloud droplet size distributions broaden significantly, in other words, larger [64].  Figure 11 further shows that more than 3 × 10 5 data points form a single-centered distribution with the center concentrated in of 0.3-0.4 and w of 0-2 m/s. According to the vertical profiles of relative dispersion (Figure 7a), vertical velocity ( Figure 5a) and cloud fraction (Figure 4), most of these center points come from 0 to 2.0 km near the cloud   (Figure 4), most of these center points come from 0 to 2.0 km near the cloud base. The data points with large relative dispersion (0.4-0.5) and downdraft (w~−3 to 0 m/s) are mainly from the cloud upper part or the cloud top. Therefore, it is necessary to separately discuss the relationship between ε and w at different heights from the cloud base. In addition, with the increase of w, AF increases and ε decreases, which means that the cloud tends to be more adiabatic and cloud droplet size distribution tends to be narrower. When w is smaller, AF is also smaller and ε is larger. Therefore, evaporation during entrainment broadens cloud droplet size distributions.
The relationship between ε and w is fitted every 50 m. Figure 12a shows part of the relationships for every 400 m. Figure 12b,c shows the correlation coefficient and the slope of the linear regression equation as a function of height. Near the cloud base, the correlation coefficient is around −0.45 and the slope is around −0.02. As the height increases, the correlation coefficient and slope are both negative and their absolute values increase to the maximum absolute values of −0.59 and −0.04, respectively. As the height further increases, the number of data points decreases, and the center of the data points also moves from the area with positive w and small ε to the area with negative w and big ε, which is consistent with Figure 11. The physical mechanisms are summarized below. Near the cloud base, on the one hand, larger vertical velocity (updraft) makes cloud droplet activation more active, and more small droplets are formed to increase relative dispersion; on the other hand, larger vertical velocity increases supersaturation, promotes condensation and reduces relative dispersion. The two factors contradict each other and the negative relationship between ε and w is weak near the cloud base [65,66]. The latter factor is expected to dominate over the former one with the increasing height, because activation becomes weaker at higher heights. Therefore, the degree of negative correlation between ε and w increases to its maximum with the increasing height. As the height continues to increase, the cloud droplet size distributions tend to stabilize (Figure 10), and the relative dispersion range decreases. As a result, both of absolute values of the correlation coefficient and the slope between relative dispersion and vertical velocity decrease.
The above analysis indicates that the adiabatic fraction of the cloud is also an important factor affecting the relative dispersion of cloud droplet size distribution [17,19]. Similar to Figure 11, Figure 13 shows a negative correlation between ε and AF with the correlation coefficient of −0.79. Larger AF means that the cloud is less affected by entrainment of dry air, which results in weaker evaporation and smaller ε. Similar to Figure 12, Figure 14 shows the relationship between ε and AF for different heights. ε and AF have the weakest correlation near the cloud base, with the correlation coefficient of −0.22. Since AF has the largest value near the cloud base, it is expected to have small ε. However, as mentioned above, the strong activation of droplets near the cloud base increases ε, which contradicts the effect of AF. Therefore, the relationship between ε and AF is weak. As the height increases, the activation intensity decreases, and the influence of entrainment on ε increases; as a result, the absolute value of the negative correlation coefficient increases rapidly, with the maximum of −0.84. As the height continues to increase, the cloud droplet size distribution tends to stabilize (Figure 10), and the relative dispersion range decreases; therefore, the correlation coefficient between ε and AF decreases. The variation of slope is generally similar to that of the correlation coefficient. In the upper part of the cloud, the absolute value of slope increases, which could be due to the small AF range (<0.4). Since ε is the ratio of σ to r m , Figures S1 and S3 in the Supplementary Materials show the relationship of σ, r m with AF for all heights, respectively. Figures S2 and S4 are for different heights. Generally, the relationship of σ with AF is similar to that of ε. As expected, r m is positively correlated with AF at each height ( Figure S4). For a larger AF, entrainment and evaporation are weaker, so r m is larger. With the increasing height, the slope of the relationship between r m and AF increases, which contributes to the increasing absolute slope (negative) of the relationship between ε and AF ( Figure 14). Interestingly, when the data at all heights are combined, r m and AF are negatively correlated. As shown in Figure S2, the black ellipses schematically show the positive relationships between r m and AF at different heights. So, the negative relationship between r m and AF at all heights is caused by the varying height, in other words, the combination of the negative relationships between r m and AF at different heights. Figure 15 is plotted to show the vertical variations of ε, σ and r m considering different AF. It is clear that for the same height, ε decreases, σ decreases and r m increases for the increasing AF. This is consistent with theoretical expectation that for a more diluted cloud (smaller AF), cloud droplet size distributions are wider (larger ε and σ) and droplet sizes are smaller (smaller r m ). Figure 16 shows direct relationships of r m , σ and ε with cloud drop number concentration. All of r m , σ and ε are negatively correlated with number concentration. When AF is larger, a cloud is less affected by entrainment and number concentration is larger. As explained above, σ and ε are expected to be smaller in more adiabatic clouds. Prabha et al. [28] found that the relationships of r m , σ and ε with cloud drop number concentration were different under premonsoon and monsoon clouds. The negative relationship between ε and number concentration is consistent with many previous studies [59,67,68], but different from others [23,30].
Atmosphere 2021, 12, x FOR PEER REVIEW 15 of 23 becomes weaker at higher heights. Therefore, the degree of negative correlation between and w increases to its maximum with the increasing height. As the height continues to increase, the cloud droplet size distributions tend to stabilize (Figure 10), and the relative dispersion range decreases. As a result, both of absolute values of the correlation coefficient and the slope between relative dispersion and vertical velocity decrease. The above analysis indicates that the adiabatic fraction of the cloud is also an important factor affecting the relative dispersion of cloud droplet size distribution [17,19]. Similar to Figure 11, Figure 13 shows a negative correlation between and with the correlation coefficient of −0.79. Larger means that the cloud is less affected by entrainment of dry air, which results in weaker evaporation and smaller . Similar to Figure 12, Figure 14 shows the relationship between and for different heights. and have the weakest correlation near the cloud base, with the correlation coefficient of −0.22. Since has the largest value near the cloud base, it is expected to have small . How- ber concentration. All of , and are negatively correlated with number concentration. When is larger, a cloud is less affected by entrainment and number concentration is larger. As explained above, and are expected to be smaller in more adiabatic clouds. Prabha et al. [28] found that the relationships of , and with cloud drop number concentration were different under premonsoon and monsoon clouds. The negative relationship between and number concentration is consistent with many previous studies [59,67,68], but different from others [23,30].      , and are negatively correlated with small droplet number concentration (radius ≤10 μm, Figure 17a-c), respectively, similar to the relationships of , and with the total droplet number concentration ( Figure 16). and have weak negative correlations with large droplet number concentration (radius >10 μm, Figure 17d,e) and has a weak positive correlation with large droplet number concentration. For larger , both small and large droplet number concentrations are larger due to weaker entrainment and dilution. Prabha et al. [28] analyzed the relationships of , and with small droplet number concentration and . They found that the small droplet number concentration decreased when increased; for a given , decreased when the small droplet number concentration increased; was larger in the more adiabatic clouds; was negatively correlated with the small droplet number concentration in two premonsoon cases and one monsoon case.
In addition to using , the cloud droplet size distributions at different distances from the cloud core can also be analyzed to explore the impact of entrainment; at different heights, divide into different ranges and plot the cloud droplet size distribution in 4 typical ranges at the 1 km height, as shown in Figure 18. At the same height, from the edge to the core with increasing , decreases, increases slightly, and the peak of the cloud droplet size distribution gradually moves towards the big droplets. The reason is that the degree of adiabaticity increases from the boundary to the core, and the effect of evaporation due to entrainment is weakened. r m , σ and ε are negatively correlated with small droplet number concentration (radius ≤10 µm, Figure 17a-c), respectively, similar to the relationships of r m , σ and ε with the total droplet number concentration ( Figure 16). σ and r m have weak negative correlations with large droplet number concentration (radius >10 µm, Figure 17d,e) and ε has a weak positive correlation with large droplet number concentration. For larger AF, both small and large droplet number concentrations are larger due to weaker entrainment and dilution. Prabha et al. [28] analyzed the relationships of r m , σ and ε with small droplet number concentration and AF. They found that the small droplet number concentration decreased when AF increased; for a given AF, r m decreased when the small droplet number concentration increased; σ was larger in the more adiabatic clouds; ε was negatively correlated with the small droplet number concentration in two premonsoon cases and one monsoon case.
In addition to using AF, the cloud droplet size distributions at different distances from the cloud core can also be analyzed to explore the impact of entrainment; at different heights, divide q c into different ranges and plot the cloud droplet size distribution in 4 typical ranges at the 1 km height, as shown in Figure 18. At the same height, from the edge to the core with increasing q c , ε decreases, r m increases slightly, and the peak of the cloud droplet size distribution gradually moves towards the big droplets. The reason is that the degree of adiabaticity increases from the boundary to the core, and the effect of evaporation due to entrainment is weakened.

Summary
To analyze the factors affecting relative dispersion of cloud droplet size distributions, shallow convective clouds are simulated by large eddy simulation. The simulation results are first verified by aircraft observations by comparing vertical profiles of liquid water

Summary
To analyze the factors affecting relative dispersion of cloud droplet size distributions, shallow convective clouds are simulated by large eddy simulation. The simulation results are first verified by aircraft observations by comparing vertical profiles of liquid water Figure 18. Cloud droplet size distributions in different cloud water mixing ratio (q c ) at the 1 km height. ε is relative dispersion and r m is mean radius. The simulation is from 2000 (UTC) to 2300 (UTC) on 5 August 1995.

Summary
To analyze the factors affecting relative dispersion of cloud droplet size distributions, shallow convective clouds are simulated by large eddy simulation. The simulation results are first verified by aircraft observations by comparing vertical profiles of liquid water potential temperature and water vapor mixing ratio in cloud and environment. After that, vertical profiles of vertical velocity, buoyancy, latent heat, cloud water mixing ratio and cloud adiabatic fraction are plotted to understand the general macro characteristics of the simulated shallow cumulus clouds.
Since relative dispersion is the ratio of standard deviation to mean radius, it is naturally important to analyze the effects of standard deviation and mean radius on relative dispersion. From the cloud base to 0.2 km, relative dispersion decreases with the increasing height, which is mainly due to negative r m − σ and dσ dz − dr m dz . The increase in relative dispersion with the increasing height in the range of 0.2-1 km is mainly due to positive r m − σ. Relative dispersion changes slowly with height above 1 km.
Besides standard deviation and mean radius, relative dispersion is significantly affected by vertical velocity in cloud and entrainment-mixing processes. Both of the relationships between relative dispersion and vertical velocity, and between relative dispersion and adiabatic fraction are negative. Larger positive vertical velocity increases supersaturation and corresponds to narrower cloud droplet size distributions. Larger adiabatic fraction means that the cloud is less affected by evaporation during entrainment of dry air, in other words, narrower cloud droplet size distributions. The impact of entrainment is also examined by analyzing cloud droplet size distributions at different distances from the cloud core. When entrainment weakens from the edge to the core, relative dispersion decreases.
It is interesting to note that the negative relationships are relatively weak near the cloud base, become the strongest in the middle of clouds, and then weaken with further increasing height. The weak relationship between relative dispersion and vertical velocity near the cloud base is due to the two contradicting factors: larger vertical velocity causes more active droplet activation to increase relative dispersion; as mentioned above, larger vertical velocity also reduces relative dispersion by increasing supersaturation. With the increasing height, the latter factor dominates the former one to produce the strongest relationship between relative dispersion and vertical velocity. The relationship weakens with further increasing height because the relative dispersion range decreases. Similarly, the weak relationship between relative dispersion and adiabatic fraction near the cloud base is related to the strong activation of droplets. With increasing height, the influence of entrainment dominates the activation intensity, thus the negative relationship strengthens. The relative dispersion range decreases with further increasing height, so the negative relationship weakens.
It is confirmed that when vertical velocity and entrainment dominate, relative dispersion is negatively correlated with number concentration. The mechanisms are that supersaturation is larger for larger vertical velocity, which further causes stronger activation (larger number concentration) and stronger condensational growth (smaller relative dispersion) [59,60]; for smaller adiabatic fractions, entrainment is stronger, which decreases number concentration and increases relative dispersion [69,70]. However, this study does not consider the effects of aerosol, which may cause a positive correlation between relative dispersion and number concentration [23]. In future studies, the effects of aerosol on relative dispersion, vertical velocity and adiabatic fraction will be included and the relationship between relative dispersion and number concentration will be analyzed more comprehensively.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/atmos12040485/s1, Figure S1: The relationship between standard deviation (σ) and cloud adiabatic fraction (AF) at all heights, Figure S2: (a) Relationships between standard deviation (σ) and adiabatic fraction (AF) at typical heights; the varia-tions of (b) correlation coefficients and (c) slopes of σ vs AF at every 50 m with height, Figure S3: The relationship between mean radius (rm) and cloud adiabatic fraction (AF) at all heights. The ellipses schematically show the positive relationships between rm and AF at different heights, Figure S4: (a) Relationships between mean radius (rm) and adiabatic fraction (AF) at typical heights; the variations of (b) correlation coefficients and (c) slopes of rm vs AF at every 50 m with height.