Articles | Volume 19, issue 12
Atmos. Chem. Phys., 19, 7955–7971, 2019

Special issue: VAMOS Ocean-Cloud-Atmosphere-Land Study (VOCALS) (ACP/OS inter-journal...

Atmos. Chem. Phys., 19, 7955–7971, 2019

Research article 18 Jun 2019

Research article | 18 Jun 2019

Exploring aerosol–cloud interaction using VOCALS-REx aircraft measurements

Exploring aerosol–cloud interaction using VOCALS-REx aircraft measurements
Hailing Jia1,2, Xiaoyan Ma1,2, and Yangang Liu3 Hailing Jia et al.
  • 1Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China
  • 2Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science & Technology, Nanjing 210044, China
  • 3Environmental and Climate Sciences Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA

Correspondence: Xiaoyan Ma (


In situ aircraft measurements obtained during the VAMOS (Variability of the American Monsoons) Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx) field campaign are analyzed to study the aerosol–cloud interactions in the stratocumulus clouds over the southeastern Pacific Ocean (SEP), with a focus on three understudied topics (separation of aerosol effects from dynamic effects, dispersion effects, and turbulent entrainment-mixing processes). Our analysis suggests that an increase in aerosol concentration tends to simultaneously increase both cloud droplet number concentration (Nd) and relative dispersion (ε), while an increase in vertical velocity (w) often increases Nd but decreases ε. After constraining the differences of cloud dynamics, the positive correlation between ε and Nd becomes stronger, implying that perturbations of w could weaken the aerosol influence on ε and hence result in an underestimation of dispersion effect. A comparative analysis of the difference of cloud microphysical properties between the entrainment and non-entrainment zones suggests that the entrainment-mixing mechanism is predominantly extremely inhomogeneous in the stratocumulus that capped by a sharp inversion, whereby the variation in liquid water content (25 %) is similar to that of Nd (29 %) and the droplet size remains approximately constant. In entrainment zone, drier air entrained from the top induces fewer cloud droplets with respect to total in-cloud particles (0.56±0.22) than the case in the non-entrainment zone (0.73±0.13) by promoting cloud droplet evaporation. This study is helpful in reducing uncertainties in dispersion effects and entrainment mixing for stratocumulus, and the results of this study may benefit cloud parameterizations in global climate models to more accurately assess aerosol indirect effects.

1 Introduction

Stratocumulus clouds play a key role in the radiative energy budget of the Earth by reflecting incoming shortwave radiation and thus cooling the planet surface and offsetting the warming by greenhouse gases (Hartmann et al., 1992). Stratocumulus clouds are susceptible to aerosols, i.e., the aerosol indirect effect (Twomey, 1974; Albrecht, 1989), which currently remain with large uncertainties (Lohmann and Feichter, 2005; Chen and Penner, 2005; Carslaw et al., 2013; McCoy et al., 2017).

Globally, marine stratocumulus clouds overlaying the southeastern Pacific Ocean (SEP) are the largest and most persistent clouds (Klein and Hartmann, 1993; Bretherton et al., 2004). Sources of anthropogenic aerosols from the Chilean and Peruvian coasts, in contrast with the relatively clean air masses from the Pacific Ocean, make the SEP an ideal region for exploring the interaction between aerosols and stratocumulus cloud-topped boundary layers. The cloud properties from satellite retrievals exhibit a gradient off the northern Chile shore. For example, the cloud droplet number concentration decreased from 160 to 40 cm−3 (George and Wood, 2010), and the cloud droplet effective radius increased from 8 to 14 µm from the coast to approximately 1000 km offshore (Wood et al., 2007). This gradient is plausibly attributable to anthropogenic aerosols near the coast. Huneeus et al. (2006) found that during easterly wind events, sulfate increased by 1 order of magnitude over the SEP, which resulted in a 1.6-fold to 2-fold increase in cloud droplet number concentration. Based on observations from satellites and cruises, Wood et al. (2008) suggested that open cellular convection within an overcast stratocumulus is associated with reduced aerosol concentration and air masses not passing through the Chilean coast, which further confirms the impact of aerosols on stratocumulus over the SEP. However, it is difficult to establish the generality of previous studies based on satellite remote sensing due to the absence of in situ observations that provide vertical profiles of cloud and aerosol and detailed in-cloud processes.

The VAMOS (Variability of the American Monsoons) Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx), which includes multiple aircraft missions, ship, and land-based measurements, took place in the region extending from the near-coast of northern Chile and southern Peru to the remote ocean in the SEP during October–November 2008 (Wood et al., 2011). The data collected during this campaign were examined to investigate the properties of aerosols, clouds, and the marine boundary layer over the SEP. For instance, Bretherton et al. (2010) found that the boundary layer was shallow and fairly well-mixed near shore but deeper and decoupled offshore. Twohy et al. (2013) found that the clouds near the shore exhibited higher aerosol concentrations, greater droplet concentrations, smaller droplet sizes, and a smaller liquid water path (LWP), and suggested a combination of anthropogenic aerosols and physically thinner clouds near the shore. However, Zheng et al. (2010) found an increase in the LWP with cloud condensation nuclei (CCN) concentrations under the similar meteorological conditions. Additionally, chemical components and sources of aerosols during the VOCALS-REx campaign have been discussed in several studies (Chand et al., 2010; Hawkins et al., 2010; Allen et al., 2011; Twohy et al., 2013; Lee et al., 2014).

Although these studies have improved our understanding of some aspects related to aerosol, cloud, and boundary layer properties over SEP, several important factors remain understudied or unexplored. First, the aerosol effect on clouds is often intertwined with the effects of other factors, especially meteorological conditions (Fan et al., 2009; Koren et al., 2010). Currently, the impact of aerosols on the shape of the cloud droplet size spectrum (i.e., dispersion effect) is reported to remain large uncertainty. The observed correlations between relative dispersion (ε) and Nd can be positive, negative, or not evident (Table 1), which could be largely attributable to the coincidentally changing cloud dynamics. Thus, it is necessary to isolate the response of ε to aerosol perturbations from meteorological effects, which, to our knowledge, has not received adequate attention in many previous studies. Second, applying different assumptions to the entrainment-mixing mechanism can have a significant impact on the cloud albedo (Grabowski, 2006; Chosson et al., 2007; Slawinska et al., 2008). Additionally, more recent studies suggested that entrainment mixing may be a possible physical interpretation for the observed anti-Twomey effect (Ma et al., 2018; Jia et al., 2019). However, it remains unclear whether the entrainment-mixing mechanism is predominantly homogeneous, inhomogeneous, or in between (Andrejczuk et al., 2009; Lehmann et al., 2009). By using cloud observations obtained from G-1 aircraft during VOCALS-REx, Yum et al. (2015) found both homogeneous and inhomogeneous mixing in their analysis and attributed the mixing to the uncertainty in the methods they used. Uncertainty in the entrainment-mixing mechanism could lead to the inaccurate assessment of aerosol indirect effects. Thus, more attention should be paid to this topic.

Table 1Correlations between ε and Nd (Na) from observation studies.

Download Print Version | Download XLSX

Based on the useful information on the microphysical properties of aerosols and clouds provided by previous studies, in this study, we conduct additional explorations regarding aerosol–cloud interactions over the SEP by employing in situ aircraft data collected by the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft during VOCALS-REx, which include the following: (a) investigating the controlling factors of cloud droplet formation (e.g., cloud dynamics and aerosols), (b) evaluating the dispersion effect under relatively constant cloud dynamical conditions, and (c) re-examining the entrainment-mixing mechanism by using a different approach to that of Yum et al. (2015).

2 Data and method

2.1 Aircraft data

The Twin Otter aircraft operated by the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) was aimed to observe aerosol and cloud microphysics and turbulence near Point Alpha (20 S, 72 W) off the coast of northern Chile from 16 October to 13 November 2008. A total of 19 flights were carried out, each of which conducting about 3 h of sampling at Point Alpha and including several soundings and horizontal legs near the ocean surface, below the cloud, near the cloud base, within the cloud, near the cloud top, and above the cloud (Fig. 1). Since all flight tracks are similar, only one track (18 October) is shown in Fig. 1. As cloud and aerosol probe measurements failed during the flight on 5 November and drizzle processes occurred on the flights on 1 and 2 November, only the observations from other 16 non-drizzling flights are included in this paper.

Figure 1The flight track in 18 October; the colors represent flight time in hours (UTC).


Both the aerosols below and above clouds and the interstitial aerosols in-cloud were obtained by the Passive Cavity Aerosol Spectrometer Probe (PCASP), which counted and sized particles with a 0.1–2.0 µm diameter with 20 bins (Zheng et al., 2011; Cai et al., 2013; Twohy et al., 2013). The CCN number concentration was observed by the CCN spectrometer at a supersaturation of 0.2 % and 0.5 %. The cloud data include cloud droplet number concentration (Nd; size range: 2.07–40.2 µm with 20 bins) from the Cloud, Aerosol and Precipitation probe (CAS); effective radius of cloud droplets (Re); and liquid water content (LWC) from the PVM-100 probe (Gerber et al., 1994). All data sets used in this study are at a frequency of 1 Hz. Zheng et al. (2011) pointed out that uncertainties of aerosols and clouds measured by these probes are within 15 %. More detailed information about the instruments and measurements aboard the CIRPAS Twin Otter aircraft during VOCALS-REx can be found in Zheng et al. (2010) and Wood et al. (2011).

2.2 Data processing

In this study, the data collected near the land during both take-off and landing are removed to ensure that only the measurements close to Point Alpha (20 S, 72 W) are analyzed. Only the clouds with LWC >0.05 g m−3 and Nd>15 cm−3 are selected for analysis. We averaged the CCN number concentrations during the legs within 200 m above the cloud top to obtain the average above-cloud CCN and within 200 m below the cloud base to obtain the mean sub-cloud CCN. During the study period, the CCN spectrometer constantly measured CCN at a supersaturation of 0.2 %, except on the first four flights, where it measured at a supersaturation of 0.5 %. In order to have a consistent comparison between all flights, we adopted the method by Zheng et al. (2011) to adjust the CCN concentration from supersaturation of 0.5 % to 0.2 % on the first four flights. The effective radius of the aerosol particle is calculated from the PCASP-measured aerosol size distributions based on the following equation:

(1) R a = n i r i 3 / n i r i 2 ,

where ni is the aerosol number concentration in the ith bin of PCASP, and ri represents the arithmetic mean radius of the ith bin.

To investigate the impact of the entrainment-mixing processes on cloud properties and aerosol–cloud interactions, we defined the entrainment zone and non-entrainment zone, respectively. Gerber et al. (2005) showed that, in the marine stratocumulus, entrainment occurs when LWC begins to decrease from the bottom of the cloud. In this paper, the entrainment and non-entrainment zone are thus defined as the regions within 20 m above and below the height of maximal LWC, respectively. Given that the two zones are both thin layers, there is little difference in the dynamical and thermos-dynamical conditions. It is therefore assumed that the difference of cloud microphysical characteristics between the two zones is only caused by entrainment.

3 Results

3.1 Vertical profiles of aerosol, cloud, and meteorological variables

Figure 2 shows the vertical profiles of temperature (Fig. 2a), relative humidity (Fig. 2b), liquid water content (Fig. 2c), cloud droplet effective radius (Fig. 2d), cloud droplet number concentration (Fig. 2e), aerosol effective radius (Fig. 2f), and the ratio of CCN to condensation nuclei (Fig. 2g) during the 16 flights. Note that the vertical altitude is normalized by the inversion height (zi), defined as the height of the maximum vertical gradient of liquid water potential temperature (Zheng et al., 2011). The normalization minimizes the effect of the variation in zi between flights, allowing for better exploration of the average boundary layer (BL) structure during VOCALS-REx.

Figure 2Vertical profiles scaled by the inversion height; (a) temperature (K), (b) relative humidity (%), (c) liquid water content (g m−3), (d)  cloud droplet effective radius (µm), (e) cloud droplet number concentration (cm−3), (f) aerosol effective radius (µm), and (g) the number concentration ratio of CCN to aerosols for all 16 non-drizzling flights. The gray lines show all individual flights, and the orange lines indicate the average profiles. The red and green lines represent the polluted (18 October) and clean (9 November) cases, respectively.


As shown in Fig. 2a, temperature (T) decreased sharply with height within the BL, which is close to the dry adiabatic lapse rate. A strong inversion occurred at the top of the BL, with an average temperature increase of approximately 10 C. Due to the decrease in T with height and the nearly constant water vapor mixing ratio within the strong mixing BL, the relative humidity (RH) increased rapidly with increasing height (Fig. 2b). T and RH reached the minimum and maximum, respectively, when zzi was close to 0.9. Near the top of the BL (0.9<z/zi<1.0), the entrainment of the dry and warm air from the free atmosphere aloft resulted in a slight increase in T and a slight decrease in RH. As zzi varied from 1 to 1.1, T increased from 11 to approximately 18 C, and RH rapidly decreased to approximately 16 % (Fig. 2a, b). The vertical profiles of T and RH are overall consistent with the observations of other marine stratocumulus clouds (Martin et al., 1994; Keil and Haywood, 2003). For the cloud properties, an average of all profiles that are normalized by zi only may be insufficient for indicating the vertical variation in clouds due to different cloud base heights of each profile. Thus, the average profiles are not shown in Fig. 2c, d, and e, and the vertical variation in cloud properties can easily be seen from the single profile. Figure 2c shows that the LWC first increased with height from the cloud base, reached the maximum at z/zi=0.9, and then decreased with further increasing height when 0.9<z/zi<1.0. The profile of Re is similar to that of LWC (Fig. 2d). The profile of Nd remains relatively constant, with a slight increase and decrease near the base and top, respectively (Fig. S1 in the Supplement), which is consistent with the results from other VOCALS-REx observations (Painemal and Zuidema, 2011). It is interesting to note that the effective radius of aerosol particles (Ra) below cloud is larger than that above cloud, which is probably attributable to the differences in aerosol sources and aerosol properties (e.g., chemical composition; Fig. 2f). The profile of CCN∕CN is similar to that of Ra (Fig. 2g), suggesting that aerosols with large sizes are more likely to become CCN (Dusek et al., 2006; Zhang et al., 2011). Larger Ra and CCN∕CN values are also found in polluted cases than in clean cases.

3.2 Relationships between aerosol and cloud properties

The relationships between aerosol and cloud properties are essential for understanding and evaluating aerosol–cloud interactions. Most studies based on satellite data have employed aerosol optical depth or the aerosol index as a proxy for CCN number concentration to investigate aerosol–cloud interactions (Koren et al., 2005, 2010; Su et al., 2010; Tang et al., 2014; Ma et al., 2014, 2018; Wang et al., 2014, 2015; Saponaro et al., 2017). However, not all aerosols in the vertical column are actually involved in cloud formation; thus, this assumption is questionable, especially when the cloud layer is decoupled from the aerosol layer. For example, a few studies have shown that aerosols have little effect on cloud properties when aerosol and cloud layers are clearly separated (Costantino and Bréon, 2010, 2013; Liu et al., 2017). To further investigate this issue, the CCN number concentrations both below cloud (sub-CCN) and above cloud (abv-CCN) are examined for their impacts on the cloud properties.

Figure 3 shows the relationships between sub-CCN and cloud properties during all 16 non-drizzling flights. The red dots denote the 10 flights with a typical well-mixed boundary layer (BL). These flights also shared similar meteorological conditions, such as inversion heights and the jump of potential temperature and total water mixing ratio across the inversion (Zheng et al., 2010), and thus can be used to isolate the response of cloud properties to aerosol perturbations. The blue dots represent the other cases in which the conditions except the typical well-mixed BL, such as strong wind shear within the BL, moist layers above clouds, a strong decoupled BL, and so on, are involved (Table 2). For the cases with a typical well-mixed BL, both LWC (Fig. 3a) and Nd (Fig. 3b) exhibited positive correlations with sub-CCN, with correlation coefficients of 0.60 and 0.79, respectively. It is worth highlighting that the similar increases in Nd and LWC led to Re having no evident correlation with sub-CCN (Fig. 3c), as expected from the conventional first aerosol indirect effect whereby a constant LWC is assumed. For the other cases (blue dots), the sub-CCN impacts on the cloud properties were not evident due to the large differences in the meteorological conditions and the BL structure.

Figure 3(a) LWC (g cm−3), (b)  Nd (cm−3), and (c) Re (µm) as a function of sub-cloud CCN concentrations (SS =0.2 %) for all 16 non-drizzling flights. The error bars through these symbols indicate the standard deviation. Red symbols are the cases with typical well-mixed BL discussed in Zheng et al. (2011), and blue symbols are for other cases. Red (black) texts are the correlation coefficient for typical well-mixed cases (all cases).


Table 2Flight information and parameters that represent the properties of entrainment during all 16 non-drizzling flights.

a,b,c w ave, w SD, and w skew are the average, standard deviation, and skewness of in-cloud vertical velocities, respectively. d,e,f PLWC, PNd, and PRe are the percentages
of reduction in LWC, Nd, and Re within entrainment zone relative to non-entrainment zone (unit: %). g,h AFent and AFnon-ent are adiabatic fraction in entrainment zone and non-entrainment zone, respectively. Here, adiabatic fraction is defined as the ratio of the measured LWC to its adiabatic LWC that is calculated using pressure and temperature
near cloud base.

Download Print Version | Download XLSX

Compared to sub-cloud CCN, the influence of above-cloud CCN on cloud properties is very weak, even for the cases with a typical well-mixed BL. The absolute values of the correlation coefficient between the abv-CCN and cloud properties are all less than 0.4 (not shown), and none of them passed the significance test (α=0.05). In this study, the above-cloud aerosol number concentration is very low (129±60 cm−3), and the inversion capping the cloud top is extremely strong, which weakens the aerosol mixing with cloud layer and hence the aerosol effects on cloud properties. Some previous studies based on aircraft observations for stratocumulus clouds also found that Nd exhibits a significantly positive correlation with sub-CCN but no correlation with abv-CCN (Martin et al., 1994; Hudson et al., 2010; Hegg et al., 2012).

Figure 4 contrasts the relationships of Nd (a) and Re (b) as functions of LWC between the most polluted (19 October) and cleanest (9 November) cases with aerosol concentrations of 647±60 and 268±35 cm−3, respectively. Also shown are the corresponding power-law fits. Although Nd and Re both increased with first increasing LWC and then leveled off, there were significant detailed differences between the polluted and clean cases. The polluted case exhibits a steeper increase in Nd with increasing LWC than the clean case when LWC is small, whereas the opposite was true for Re. The low aerosol concentrations under the clean case inhibit the increase in Nd with LWC (Fig. 4a), which hence promotes the rapid increase in Re with LWC (Fig. 4b). In contrast, there are enough particles that may potentially activated into cloud droplets under the polluted case; thus, Nd increases rapidly with LWC. As a certain amount of water is shared by large amount particles, the increase in Re is limited. The result is consistent with the study in Beijing by Zhang et al. (2011), but the difference in cloud formations between the clean and polluted conditions is less evident, which is likely attributable to the much smaller difference in aerosol concentration in this study (approximately 400 cm−3) than that in Zhang et al. (2011; approximately 7000 cm−3).

Figure 4Correlations between (a) Nd (cm−3), (b) Re (µm), and LWC (g m−3) for clean (green) and polluted (red) cases, respectively.


3.3 Cloud droplet formation and its controlling factors

Sub-cloud CCN are considered to be a good proxy for aerosols entering a cloud. However, during the actual flight, it is difficult to simultaneously collect enough samples of sub-cloud CCN and cloud droplets, which may result in statistical analysis uncertainty. This limitation can be remedied by using the total particle concentration, which equals the aerosol concentration outside the clouds and the sum of the droplet concentration and interstitial aerosol concentration inside the clouds. Interstitial aerosols are particles observed inside clouds that either have never activated into cloud droplets or have been deactivated into aerosols after cloud droplet evaporation. Kleinman et al. (2012) pointed out that the number concentration of interstitial aerosols (Ni) can be obtained either directly from the observation of in-cloud aerosols or indirectly from a number balance between sub-cloud and in-cloud particles. In this study, the interstitial aerosol properties are derived from direct measurements in the cloud. By employing aircraft observations over both land and ocean, Gultepe and Isaac (1996) found that the difference in the number concentration between the total in-cloud particles (Nd+Ni) measured directly and sub-cloud aerosols is very small. Thus, the total in-cloud particles are assumed to characterize the overall level of in-cloud aerosol concentration before activation. Figure 5 shows an example of the 18 October flight to support this assumption. It is shown that the number concentrations of sub-cloud aerosols and total in-cloud particles are very close, with values of 583±55 and 567±59 cm−3, respectively. Similar results are also found in the other flights. The average ratio of Nd+Ni to the sub-cloud aerosol concentration during all flights is 0.94, which is smaller than the value (1.29) found by Kleinman et al. (2012) based on G-1 aircraft measurements during VOCALS-REx. Therefore, the interstitial aerosol observations in this study are unlikely to be significantly interfered with by factors such as cloud droplet shatter and cloud droplet evaporation due to instrument heating, as discussed by Kleinman et al. (2012), which has the potential to create spurious extra aerosols in cloud.

Figure 5Vertical profiles of number concentrations of aerosols (Na), cloud droplets (Nd), and total in-cloud particles (Nd+Ni) during the flight on 18 October.


The relations between Nd and Nd+Ni during the 16 non-drizzling flights are shown in Fig. 6, where the color represents the in-cloud vertical velocity. All flights exhibited positive correlations between Nd and Nd+Ni, representing the aerosol–cloud interaction (IPCC, 2001, 2007, 2013; Hegg et al., 2012). In addition, the effect of dynamical conditions on cloud droplet formation is evident. As shown in Fig. 6, the data are close to the 1:1 line when the vertical velocity is relatively large; namely, the aerosols were almost entirely activated into cloud droplets. However, the data deviate from the 1:1 line when the vertical velocity is small or negative. For example, for all flights, the average ratio of Nd to Nd+Ni with vertical velocity greater than 1 m s−1 is 0.84±0.12, which is much larger than that with vertical velocity less than −1 m s−1 (0.64±0.14). The regime-dependent behavior is likely due to the high supersaturation caused by the adiabatic uplift when the vertical velocity is large (Reutter et al., 2009; Chen et al., 2016).

Figure 6Relationships between Nd and Ni+Nd during all 16 non-drizzling flights. The color represents in-cloud vertical velocities (m s−1), and gray line is 1:1 line. The mean and standard deviation of Nd/(Nd+Ni) for vertical velocity greater than 1 m s−1 (red) and less than −1 m s−1 (blue) are shown.


In addition to the dynamical conditions, aerosol microphysical properties such as size distribution and chemical components can also significantly affect the activation process (Nenes et al., 2002; Lance et al., 2004; Ervens et al., 2005; Dusek et al., 2006; McFiggans et al., 2006; Zhang et al., 2011; Almeida et al., 2014; Leck and Svensson, 2015). Since part of the aerosol population was activated, it is difficult to obtain information about aerosol size before activation. According to the Köhler theory, larger aerosols have smaller critical supersaturations, and, thus, they activate preferentially, suggesting that the effective radius of interstitial aerosols (Ri) is smaller than that of the aerosols before activation. Li et al. (2011) compared the difference in size distribution between interstitial aerosols and aerosols that have been activated to cloud droplets and found that the peak diameter of the former (0.45 µm) was much smaller than that of the latter (0.8 µm). It can be thus inferred that the size of aerosols activated to cloud droplets, and thus the size of initial aerosols, would be larger with an increase in Ri, though the quantitative relationship depends on in-cloud dynamics. Therefore, it is assumed that when compared with the data measured at different sampling locations during flight, the size of the interstitial aerosols can still represent the size of the aerosols before activation to some extent. As indicated in Fig. 7, the larger Ri is, the closer the data are to the 1:1 line, i.e., the higher the proportion of cloud droplets in total in-cloud particles (Nd/(Nd+Ni)). The averaged Nd/(Nd+Ni) for all flights is 0.76±0.13 when Ri is larger than 0.5 µm but only 0.64±0.23 when Ri is less than 0.25 µm. It is because those aerosols with large sizes are more likely to be activated into cloud droplets. Additionally, as larger aerosol particles form into larger cloud droplets (Twohy et al., 1989, 2013) that are relatively difficult to evaporate, large particles can also inhibit cloud droplet evaporation to a certain extent.

Figure 7Same as Fig. 6, but the color represents the effective radius of interstitial aerosol (Ri; µm). The mean and standard deviation of Nd/(Nd+Ni) for Ri greater than 0.5 µm (red) and less than 0.25 µm (blue) are shown.


3.4 Dispersion effect

In addition to modulating the cloud droplet number concentration, aerosols also affect the shape of cloud droplet size spectrum (referred to as the “dispersion effect”) and thereby affect the cloud albedo (Liu and Daum, 2002). When the dispersion effect is taken into account, the estimated aerosol indirect forcing could be either reduced (Liu and Daum, 2002; Peng and Lohmann, 2003; Anil Kumar et al., 2016; Pandithurai et al., 2012) or enhanced (Ma et al., 2010), i.e., the dispersion effect could act to either offset or enhance the well-known Twomey effect, which mainly depends on the sensitivity of the relative dispersion (ε; the ratio of the standard deviation to the mean radius of the cloud droplet size distribution) to the aerosol number concentration (Na). However, the dependence of ε on Na is much less studied and remains even more uncertain than that of Nd. Table 1 summarizes the observed correlations between ε and Nd (or Na), being positive, being negative, or having no obvious correlations. The different relationships are indicative of the fact that the effect of aerosol on ε is often intertwined with the effects of other factors, especially cloud dynamical conditions (Pawlowska et al., 2006; Lu et al., 2012). In this section, the relationship between ε and Nd based on in-flight and flight-averaged data is discussed to distinguish the influences of aerosol and cloud dynamics on ε.

Within an individual flight, the aerosol number concentration and chemical components can be assumed to be similar, providing an opportunity to focus on the effect of cloud dynamics to the extent possible. Here, we employ the vertical velocity (w; m s−1) as a proxy for cloud dynamical conditions. As shown in Fig. 8, the correlations between ε and Nd based on in-flight data are significantly negative during all 16 non-drizzling flights, which is mainly modulated by w; i.e., a larger w corresponds to a smaller ε but a larger Nd. High supersaturation leads to more cloud droplets to activate and grow to the same size (i.e., narrowing the droplet spectrum) when w is relatively large, but a portion of the cloud droplets may evaporate into smaller sizes and even deactivate into interstitial aerosols when w is small or even negative, resulting in a decrease in Nd and broadening of the droplet spectrum.

Figure 8Relationships between relative dispersion (ε) and Nd during all 16 non-drizzling flights, in which the color represents in-cloud vertical velocities (m s−1).


It is interesting to see from Table 1 that the correlations between ε and Nd based on in-flight data are generally negative, while the correlations based on the flight-averaged data could be either positive, negative, or even uncorrelated. The latter uncertain relationships may result from variations in the strength of cloud dynamics between flights, which would disrupt or even cancel the real influence of aerosols on relative dispersion (Liu et al., 2006; Peng et al., 2007; Lu et al., 2012). However, many previous studies did not consider the difference in cloud dynamics between flights when correlating ε and Nd, which could result in some degree of overestimation or underestimation of dispersion effect. In this study, the data of all flights were sampled over the same location, i.e., Point Alpha, which can reduce the difference in dynamical conditions caused by variations in horizontal sampling locations. In addition, we also distinguish between the flights of a typical mixed BL and others to ensure relatively similar meteorological conditions (see Sect. 3.2). Figure 9 further shows the probability distribution function of w with mean values and standard deviations for 16 non-drizzling flights. The related statistics are given in Table 2. Except for other cases (crosses; especially 24 October, 29 October, 8 November, and 13 November), the difference in the in-cloud dynamics between typical well-mixed BL flights is very small, which confirms the assumption of similar meteorological conditions. As indicated in Fig. 10a, ε and Nd were positively correlated (correlation coefficient of 0.29 and slope of 1.9×10-4) in the case of the typical well-mixed BL, indicating that increased aerosols concurrently increased ε and Nd. However, the correlation coefficient and slope decrease to 0.11 and 7.7×10-5, respectively, in the all cases (i.e., w is not constrained), implying that the influence of aerosols on the εNd relationship tends to be weaker after intertwining with the effects of cloud dynamics. Although the perturbations of cloud dynamics have been eliminated as much as possible, Nd is still likely determined by both aerosol number concentrations and updraft velocity together. Therefore, a similar statistical analysis is also conducted for sub-cloud CCN. Similar positive correlations between ε and sub-cloud CCN were found, with much-improved correlation coefficients (slopes). The correlation coefficients (slopes) were 0.67 (3.1×10-4) and 0.31 (2.1×10-4) for the cases with a typical well-mixed BL and all cases, respectively (Fig. 10b).

Figure 9Probability distribution function (unit: %) of vertical velocity (w) for 16 non-drizzling flights. Black symbols are mean values of w, and error bars through these symbols indicate the standard deviation. Circles are the cases with typical well-mixed BL, and crosses represents the other cases.


Figure 10Relative dispersion (ε) as a function of (a) Nd and (b) sub-cloud CCN concentrations (SS =0.2 %) for all flights. The error bars through these symbols indicate the standard deviation. Red symbols are the cases with typical well-mixed BL, and blue symbols are for other cases. Red (black) texts are the correlation coefficient and slope for typical well-mixed cases (all cases).


3.5 Entrainment in stratocumulus

Entrainment is a key process that occurs in the clouds and plays an important role in the formation and evolution of clouds and the change of droplet spectrum as well as the aerosol indirect effect (Chen et al., 2014, 2015; Andersen and Cermak, 2015). The nature of entrainment is related to the cloud type. Entrainment in cumulus is primarily lateral with strong dilution of the cloud, which induces LWC to decrease rapidly to approximately 20 % of its adiabatic value (Warner, 1955). Entrainment in stratocumulus is mainly determined by the strength of the gradients in buoyancy and horizontal winds (Wang and Albrecht 1994; Gerber et al., 2005; de Roode and Wang 2007; Wood, 2012) and proceeds from the top and mostly affects a thin layer (Gerber et al., 2005), whose dilution effect is much weaker than that in cumulus (Warner, 1955, 1969a, b; Blyth et al., 1988; Gerber et al., 2008; Burnet and Brenguier, 2007; Haman et al., 2007). Aircraft observations of marine stratocumulus clouds showed that the vertical profile of LWC is essentially the same as the adiabatic profile; i.e., the cloud is almost adiabatic (Keil and Haywood, 2003). Furthermore, it remains unclear whether the subsequent entrainment-mixing mechanism is predominantly homogeneous, inhomogeneous, or in between (Andrejczuk et al., 2009; Lehmann et al., 2009). Some previous studies have shown that stratocumulus is generally dominated by the inhomogeneous mechanism (Pawlowska et al., 2000; Burnet and Brenguier, 2007; Haman et al., 2007; Lu et al., 2011; Yum et al., 2015). By employing a different vertical description in characterizing the region near the cloud top (Malinowski et al., 2013), Gerber et al. (2016) noted that both extremely inhomogeneous mixing and homogenous mixing play a role in unbroken stratocumulus, but the reduction in cloud droplet effective radius appears to be secondary in comparison to the dilution process that preserves the relative shape of the droplet spectrum.

To explore entrainment in stratocumulus during VOCALS-REx, we first compare the differences in cloud microphysics between the entrainment and non-entrainment zones near the cloud top. Here, the entrainment and non-entrainment zones are defined as the regions within 20 m above and below the height of the maximal LWC, respectively. As anticipated, the adiabatic fraction (AF; the ratio of the measured LWC to its adiabatic value) in the entrainment zone (AFent) is generally lower than that in the non-entrainment zone (AFnon-ent), with mean values for all flights of 0.64 and 0.77, respectively (Table 2), which further confirms the rationality of dividing the two zones. Compared with the non-entrainment zone, the peak radius of cloud droplets in the entrainment zone has few changes (Fig. 11), and the effective radius of cloud droplets (Re) increases by only 2 % (Table 2). However, Nd and LWC decrease significantly on average, by 29 % and 25 %, respectively (Table 2), especially during the flights on 18 October, 4 November, 9 November, and 13 November, where Nd decreases by 60 %, 56 %, 56 %, and 59 % and LWC decreases by 56 %, 62 %, 56 %, and 59 %, respectively (Table 2). It is suggested that dry and warm air entrained from the cloud top dilutes Nd and LWC by similar amounts, while the size of droplets is relatively unaffected, which is thought of as being an extremely inhomogeneous entrainment-mixing process. Additionally, both PLWC andPNd are negatively correlated with AFent∕AFnon-ent, with correlation coefficients of −0.60 and −0.47, respectively, indicating the dependence of the LWC and Nd changes on the adiabatic fraction changes (Fig. S2), where PLWC and PNd are the reduction percentages in LWC and Nd within the entrainment zone relative to the non-entrainment zone.

Figure 11Number size distributions of cloud droplets in the entrainment (orange) and non-entrainment zones (blue) during all 16 non-drizzling flights.


The flight on 18 October with strong entrainment is chosen to investigate the difference between the entrainment and non-entrainment zones. As shown in Fig. 12b, dry and warm air entrained from the top reduced the RH in the entrainment zone by 9 % on average and hence acted to accelerate cloud droplet evaporation. Consequently, Nd/(Nd+Ni) in the entrainment zone (0.56±0.22) is much lower than that in non-entrainment zone (0.73±0.13; Fig. 12c). Additionally, the relative dispersion in the entrainment zone is generally larger than that in the non-entrainment zone (Fig. 12d), implying that drier air entrained from the top could broaden the cloud droplet spectrum by promoting cloud droplet evaporation. Some previous observations also showed that ε with a low AF tends to be larger than that with a high AF and attributed it to the effect of entrainment mixing (Pawlowska et al., 2006; Lu et al., 2009). It is noted that the occurrence frequency of Ri in the entrainment zone is significantly higher than that in the non-entrainment zone when Ri<0.38µm, but the opposite is true when Ri>0.5µm (Fig. 12a). This result suggests that in addition to dry and warm air, small particles are also entrained into clouds from the top (Fig. 2f) and that large particles are detrained out of the clouds simultaneously. However, the inversion capping the typical stratocumulus is usually too strong to excite significant updrafts near the cloud top (Stevens, 2002; Wood, 2012; Malinowski et al., 2013). Ghate et al. (2010) found that vertical velocities near the top of stratocumulus tend towards zero overall, with only approximately 4 % of updrafts being stronger than 0.5 m s−1. Therefore, although smaller aerosols are entrained into the entrainment zone, these aerosols seem unlikely to influence droplet formation. The effect of entrainment mixing on stratocumulus is mainly governed by the entrained dry air rather than small aerosols.

Figure 12Probability density functions of (a) Ri (µm), (b) RH (%), (c) Nd/(Nd+Ni), and (d) ε in the entrainment (orange) and non-entrainment zones (blue) during the flight on 18 October.


4 Summary

By using in situ aircraft data collected by the CIRPAS Twin Otter aircraft at Point Alpha during VOCALS-REx from 16 October to 13 November 2008, aerosol–cloud interactions are investigated with a focus on understudied factors, including separation of aerosol effects from dynamic effects, dispersion effects, and turbulent entrainment-mixing processes.

Vertical profiles of aerosol, cloud, and meteorological variables indicated that the BL was capped by a sharp inversion during 16 non-drizzling flights. Cloud properties, such as LWC and Nd, are positively correlated with sub-cloud CCN number concentration, with correlation coefficients of 0.60 and 0.79, respectively. No evident correlation was found between cloud properties and above-cloud CCN number concentrations. This is mainly due to the low aerosol number concentrations above cloud (129±60 cm−3) and the extremely strong inversion capping the cloud top, which inhibits the mixing of the above-cloud aerosols with the cloud layer. Therefore, the influence of the above-cloud CCN on cloud properties is weaker than the sub-cloud CCN.

The results showed that both dynamical conditions and aerosol microphysical properties have significant effects on cloud droplet formation. In the case of large vertical velocity and aerosol size, the cloud droplet proportion of total in-cloud particles is relatively high (e.g., 0.84±0.12 and 0.76±0.13, respectively); i.e., cloud droplets form more easily. Although aerosol chemical components are also critical in cloud droplet formation (Nenes et al., 2002; Lance et al., 2004; Ervens et al., 2005; McFiggans et al., 2006; Wang et al., 2008; Almeida et al., 2014), this topic was not discussed in this study due to the unavailability of measurements.

The correlations between ε and Nd based on the in-flight data, representing the w-induced correlation, are significantly negative, while the correlations derived from flight-averaged data (i.e., aerosol-induced correlation) are positive. This finding implies that an increase in aerosol concentration tends to concurrently increase ε and Nd, while an increase in w often increases Nd but decreases ε, which agrees with the theoretical analysis (Liu et al., 2006). After constraining the differences in cloud dynamics between flights, positive εNd correlations become stronger, indicating that perturbations of w could weaken the influence of aerosols on ε and hence may result in an underestimation of aerosol dispersion effect. Thus, this finding highlights the necessity of isolating the relative dispersion response to aerosol perturbations from dynamic effects when investigating the aerosol dispersion effect and estimating aerosol indirect forcing.

Overall, the entrainment in stratocumulus is quite weak and close to being adiabatic in some cases. In this study, the difference in cloud microphysics between the entrainment and non-entrainment zones indicated that the entrainment in stratocumulus is mostly dominated by an extremely inhomogeneous entrainment-mixing mechanism. On average, the entrainment reduced Nd and LWC by 29 % and 25 %, respectively, while having little effect on Re (only increases by 1.8 %). During the flights on 18 October, 4 November, 9 November, and 13 November, the entrainment was relatively strong and diluted Nd and LWC by about 50 %. In the entrainment zone, the drier air entrained from the top resulted in a smaller Nd/(Nd+Ni) (0.56±0.22) than that in the non-entrainment zone (0.73±0.13). This implies that entrainment may significantly influence cloud droplet formation and therefore influence the cloud properties near the top by promoting cloud droplet evaporation. Furthermore, we found that the relative dispersion in the entrainment zone is larger than that in the non-entrainment zone. In addition to the dry and warm air, aerosols with smaller sizes are also entrained into the entrainment zone, but these aerosols seem unlikely to influence cloud droplet formation due to the negligible droplet nucleation near the stratocumulus top. That is, the effect of entrainment mixing on stratocumulus is mainly determined by the entrained dry air instead of the aerosols with properties that are different from those near the cloud base. These results seem at odds with some studies on cumulus clouds. Slawinska et al. (2012) found that in a shallow cumulus, a significant part (40 %) of aerosols is activated above the cloud base (secondary activation), which is dominated by entrained aerosols. Using large-eddy simulations (LES), Hoffmann et al. (2015) suggested that sub-cloud aerosols and laterally entrained aerosols contribute to all activated aerosols inside the cloud by fractions of 70 % and 30 %, respectively. Evidently, the topics of how and to what extent entrained aerosols with properties that are different from sub-cloud aerosols can affect the formation and evolution of clouds merit further exploration.

Data availability

The aircraft measurement data during VOCALS-REx were obtained from the public FTP at (last access: 30 May 2019). The navigation and state parameters are available from (Albrecht, 2011a). The cloud data are available from (Albrecht, 2011b). The CCN data are available from (Albrecht, 2011c). The aerosol data are available from (Albrecht, 2011d).


The supplement related to this article is available online at:

Author contributions

HJ and XM designed the study and the statistical analysis. HJ processed the data and drafted the paper. XM and YL validated and debugged the results. All authors contributed to revising the paper.

Competing interests

The authors declare that they have no conflict of interest.


We are grateful for the dedicated efforts of several members of the support staff and scientists in making the observations from the CIRPAS Twin Otter during VOCALS-REx. We also thank Bruce Albrecht at the University of Miami for kindly providing the aerosol, cloud, and meteorological variable observations, which are the basis of this paper.

Financial support

This research has been supported by the National Natural Science Foundation of China grants (NSFC; grant no. 41475005), the National Key R&D Program of China grants (grant no. 2016YFA0600404), the NSFC (grant no. 41675004), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (grant no. KYCX18_1005), and the funding of Jiangsu innovation & entrepreneurship team. Yangang Liu is supported by the US Department of Energy's Atmospheric System Research (ASR) program.

Review statement

This paper was edited by Hugh Coe and reviewed by three anonymous referees.


Albrecht, B. A.: Aerosols, cloud microphysics, and fractional cloudiness, Science, 245, 1227–1230, 1989. 

Albrecht, B.: University of Miami, CIRPAS Twin Otter Navigation and State Parameters, Version 1.0, UCAR/NCAR – Earth Observing Laboratory, available at: (last access: 30 May 2019), 2011a. 

Albrecht, B.: University of Miami, CIRPAS Twin Otter CAPS Probe Data, Version 1.0, UCAR/NCAR – Earth Observing Laboratory, available at: (last access: 30 May 2019), 2011b. 

Albrecht, B.: University of Miami, CIRPAS Twin Otter CCN Data, Version 1.0, UCAR/NCAR – Earth Observing Laboratory, available at: (last access: 30 May 2019), 2011c. 

Albrecht, B.: University of Miami, CIRPAS Twin Otter PCASP Data, Version 1.0, UCAR/NCAR – Earth Observing Laboratory, available at: (last access: 30 May 2019), 2011d. 

Allen, G., Coe, H., Clarke, A., Bretherton, C., Wood, R., Abel, S. J., Barrett, P., Brown, P., George, R., Freitag, S., McNaughton, C., Howell, S., Shank, L., Kapustin, V., Brekhovskikh, V., Kleinman, L., Lee, Y.-N., Springston, S., Toniazzo, T., Krejci, R., Fochesatto, J., Shaw, G., Krecl, P., Brooks, B., McMeeking, G., Bower, K. N., Williams, P. I., Crosier, J., Crawford, I., Connolly, P., Allan, J. D., Covert, D., Bandy, A. R., Russell, L. M., Trembath, J., Bart, M., McQuaid, J. B., Wang, J., and Chand, D.: South East Pacific atmospheric composition and variability sampled along 20 S during VOCALS-REx, Atmos. Chem. Phys., 11, 5237–5262,, 2011. 

Almeida, G. P., Brito, J., Morales, C. A., Andrade, M. F., and Artaxo, P.: Measured and modelled cloud condensation nuclei (CCN) concentration in São Paulo, Brazil: the importance of aerosol size-resolved chemical composition on CCN concentration prediction, Atmos. Chem. Phys., 14, 7559–7572,, 2014. 

Andersen, H. and Cermak, J.: How thermodynamic environments control stratocumulus microphysics and interactions with aerosols, Environ. Res. Lett., 10, 24004,, 2015. 

Andrejczuk, M., Grabowski, W. W., Malinowski, S. P., and Smolarkiewicz, P. K.: Numerical simulation of cloud-clear air interfacial mixing: Homogeneous versus inhomogenous mixing, J. Atmos. Sci., 66, 2493–2500, 2009. 

Anil Kumar, V., Pandithurai, G., Leena, P. P., Dani, K. K., Murugavel, P., Sonbawne, S. M., Patil, R. D., and Maheskumar, R. S.: Investigation of aerosol indirect effects on monsoon clouds using ground-based measurements over a high-altitude site in Western Ghats, Atmos. Chem. Phys., 16, 8423–8430,, 2016. 

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

Bretherton, C. S., Uttal, T., Fairall, C. W., Yuter, S. E., Weller, R. A., Baumgardner, D., Comstock, K., Wood, R., and Raga, G. B.: The Epic 2001 Stratocumulus Study, B. Am. Meteorol. Soc., 85, 967–977, 2004. 

Bretherton, C. S., Wood, R., George, R. C., Leon, D., Allen, G., and Zheng, X.: Southeast Pacific stratocumulus clouds, precipitation and boundary layer structure sampled along 20 S during VOCALS-REx, Atmos. Chem. Phys., 10, 10639–10654,, 2010. 

Burnet, F. and Brenguier, J. L.: Observational study of the entrainment-mixing process in warm convective clouds, J. Atmos. Sci., 64, 1995–2011, 2007. 

Cai, Y., Snider, J. R., and Wechsler, P.: Calibration of the passive cavity aerosol spectrometer probe for airborne determination of the size distribution, Atmos. Meas. Tech., 6, 2349–2358,, 2013. 

Carslaw, K. S., Lee, L. A., Reddington, C. L., Pringle, K. J., Rap, A., Forster, P. M., Mann, G. W., Spracklen, D. V., Woodhouse, M. T., Regayre, L. A., and Pierce, J. R.: Large contribution of natural aerosols to uncertainty in indirect forcing, Nature, 503, 67–71, 2013. 

Chand, D., Hegg, D. A., Wood, R., Shaw, G. E., Wallace, D., and Covert, D. S.: Source attribution of climatically important aerosol properties measured at Paposo (Chile) during VOCALS, Atmos. Chem. Phys., 10, 10789–10801,, 2010. 

Chen, J., Liu, Y., Zhang, M., and Peng, Y.: New understanding and quantification of the regime dependence of aerosol-cloud interaction for studying aerosol indirect effects, Geophys. Res. Lett., 43, 1780–1787,, 2016. 

Chen, Y. and Penner, J. E.: Uncertainty analysis for estimates of the first indirect aerosol effect, Atmos. Chem. Phys., 5, 2935–2948,, 2005. 

Chen, Y.-C., Christensen, M. W., Stephens, G. L., and Seinfeld, J. H.: Satellite-based estimate of global aerosol-cloud radiative forcing by marine warm clouds, Nat. Geosci., 7, 643–646, 2014. 

Chen, Y.-C., Christensen, M. W., Diner, D. J., and Garay, M. J.: Aerosol-cloud interactions in ship tracks using TerraMODIS/MISR, J. Geophys. Res.-Atmos., 120, 2819–2833, 2015. 

Chosson, F., Brenguier, J.-L., and Schuller, L.: Entrainment-mixing and radiative transfer simulation in boundary layer clouds, J. Atmos. Sci., 64, 2670–2682, 2007. 

Costantino, L. and Bréon, F.-M.: Analysis of aerosol-cloud interaction from multi-sensor satellite observations, Geophys. Res. Lett., 37, L11801,, 2010. 

Costantino, L. and Bréon, F.-M.: Aerosol indirect effect on warm clouds over South-East Atlantic, from co-located MODIS and CALIPSO observations, Atmos. Chem. Phys., 13, 69–88,, 2013. 

de Roode, S. R. and Wang, Q.: Do stratocumulus clouds detrain? FIRE I data revisited, Bound.-Lay. Meteorol., 122, 479–491,, 2007. 

Dusek, U., Frank, G. P., Hildebrandt, L., Curtius, J., Schneider, J., Walter, S., Chand, D., Drewnick, F., Hings, S., Jung, D., Borrmann, S., and Andreae, M. O.: Size matters more than chemistry for cloud-nucleating ability of aerosol particles, Science, 312, 1375–1378, 2006. 

Ervens, B., Feingold, G., and Kreidenweis, S. M.: Influence of water-soluble organic carbon on cloud drop number concentration, J. Geophys. Res., 110, D18211,, 2005. 

Fan, J., Yuan, T., Comstock, J. M., Ghan, S., Khain, A., Leung, L. R., Li, Z., Martins, J. V., and Ovchinnikov, M.: Dominant role by vertical wind shear in regulating aerosol effects on deep convective clouds, J. Geophys. Res., 114, D22206,, 2009. 

George, R. C. and Wood, R.: Subseasonal variability of low cloud radiative properties over the southeast Pacific Ocean, Atmos. Chem. Phys., 10, 4047–4063,, 2010. 

Gerber, H., Arends, B. G., and Ackerman, A. S.: A new microphysics sensor for aircraft use, Atmos. Res., 31, 235–252, 1994. 

Gerber, H., Frick, G., Malinowski, S. P., Brenguier J.-L., and Burnet, F.: Holes and entrainment in stratocumulus, J. Atmos. Sci., 62, 443–459, 2005. 

Gerber, H., Frick, G., Jensen, J. B., and Hudson, J. G.: Entrainment, mixing and microphysics in trade-wind cumulus, J. Meteorol. Soc. Jpn., 86A, 87–106, 2008. 

Gerber, H., Malinowski, S. P., and Jonsson, H.: Evaporative and radiative cooling in POST stratocumulus, J. Atmos. Sci., 73, 3877– 3884,, 2016. 

Ghate, V. P., Albrecht, B. A., and Kollias, P.: Vertical velocity structure of nonprecipitating continental boundary layer stratocumulus clouds, J. Geophys. Res., 115, D13204,, 2010. 

Grabowski, W. W.: Indirect impact of atmospheric aerosols in idealized simulations of convective-radfiative equilibrium, J. Climate, 19, 4664–4682, 2006. 

Gultepe, I. and Isaac, G.: The relationship between cloud droplet and aerosol number concentrations for climate models, Int. J. Climatol., 16, 941–946,<941::AID-JOC57>3.0.CO;2-O, 1996. 

Haman, K. E., Malinowski, S. P., Kurowski, M. J., Gerber, H., and Brenguier, J.-L.: Small scale mixing processes at the top of a marine stratocumulus – A case study, Q. J. Roy. Meteor. Soc., 133, 213–226, 2007. 

Hartmann, D. L., Ockert-Bell, M. E., and Michelsen, M. L.: The effect of cloud type on Earth's energy balance: Global analysis, J. Climate, 5, 1281–1304, 1992. 

Hawkins, L. N., Russell, L. M., Covert, D. S., Quinn, P. K., and Bates, T. S.: Carboxylic acids, sulfates, and organosulfates in processed continental organic aerosol over the southeast Pacific Ocean during VOCALS-Rex 2008, J. Geophys. Res., 115, D13201,, 2010. 

Hegg, D. A., Covert, D. S., Jonsson, H. H., and Woods, R. K.: A simple relationship between cloud drop number concentration and precursor aerosol concentration for the regions of Earth's large marine stratocumulus decks, Atmos. Chem. Phys., 12, 1229–1238,, 2012. 

Hoffmann, F., Raasch, S., and Noh, Y.: Entrainment of aerosols and their activation in a shallow cumulus cloud studied with a coupled LCM–LES approach, Atmos. Res., 156, 43–57, 2015. 

Hudson, J. G., Noble, S., and Jha, V.: Stratus cloud supersaturations, Geophys. Res. Lett., 37, 21813,, 2010. 

Hudson, J. G., Noble, S., and Jha, V.: Cloud droplet spectral width relationship to CCN spectra and vertical velocity, J. Geophys. Res., 117, D11211,, 2012. 

Huneeus, N., Gallardo, L., and Rutllant, J. A.: Offshore transport episodes of anthropogenic sulfur in northern Chile: Potential impact on the stratocumulus cloud deck, Geophys. Res. Lett., 33, L19819,, 2006. 

IPCC: Climate Change: The Scientific Basis: contributions of Working Group I to the Third Assessment Report of the IPCC, edited by: Houghton, J. T., Ding, Y., Griggs, D. J., Noguer, M., van der Linden, P. J., Dai, X., Maskell, K., and Johnson, C. A., Cambridge University Press, New York, 881 pp., 2001. 

IPCC: Climate Change: The Physical Science Basis: contribution of Working Group 1 to the Fourth Assessment Report of the IPCC, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, New York, 996 pp., 2007. 

IPCC: Climate change: The Physical Science Basis: contribution of Working group I to the Fifth Assessment Report of the IPCC, edited by: Stocker, T. F., Dahe, Q., Plattner, G. K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P. M., Cambridge University Press, Cambridge, 1535 pp., 2013. 

Jia, H., Ma, X., Quaas, J., Yin, Y., and Qiu, T.: Is Positive Correlation between Cloud Droplet Effective Radius and Aerosol Index over Land Due to Retrieval Artifacts or Real Physical Processes?, Atmos. Chem. Phys. Discuss.,, in review, 2019. 

Keil, A. and Haywood, J. M.: Solar radiative forcing by biomass burning aerosol particles during safari 2000: a case study based on measured aerosol and cloud properties, J. Geophys. Res., 108, D13,, 2003. 

Klein, S. A. and Hartmann, D. L.: The Seasonal Cycle of Low Stratiform Clouds, J. Climate, 6, 1587–1606, 1993. 

Kleinman, L. I., Daum, P. H., Lee, Y.-N., Lewis, E. R., Sedlacek III, A. J., Senum, G. I., Springston, S. R., Wang, J., Hubbe, J., Jayne, J., Min, Q., Yum, S. S., and Allen, G.: Aerosol concentration and size distribution measured below, in, and above cloud from the DOE G-1 during VOCALS-REx, Atmos. Chem. Phys., 12, 207–223,, 2012. 

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

Koren, I., Feingold, G., and Remer, L. A.: The invigoration of deep convective clouds over the Atlantic: aerosol effect, meteorology or retrieval artifact?, Atmos. Chem. Phys., 10, 8855–8872,, 2010. 

Lance, S., Nenes, A., and Rissman, T. A.: Chemical and dynamical effects on cloud droplet number: Implications for estimates of the aerosol indirect effect, J. Geophys. Res., 109, D22208,, 2004. 

Leck, C. and Svensson, E.: Importance of aerosol composition and mixing state for cloud droplet activation over the Arctic pack ice in summer, Atmos. Chem. Phys., 15, 2545–2568,, 2015. 

Lee, Y.-N., Springston, S., Jayne, J., Wang, J., Hubbe, J., Senum, G., Kleinman, L., and Daum, P. H.: Chemical composition and sources of coastal marine aerosol particles during the 2008 VOCALS-REx campaign, Atmos. Chem. Phys., 14, 5057–5072,, 2014. 

Lehmann, K., Siebert, H., and Shaw, R. A.: Homogeneous and Inhomogeneous Mixing in Cumulus Clouds: Dependence on Local Turbulence Structure, J. Atmos. Sci., 66, 3641–3659, 2009. 

Li, W., Li, P., Sun, G., Zhou, S., Yuan, Q., and Wang, W.: Cloud residues and interstitial aerosols from non-precipitating clouds over an industrial and urban area in northern China, Atmos. Environ., 45, 2488–2495, 2011. 

Liu, Y., Daum, P. H., and Yum, S. S.: Analytical expression for the relative dispersion of the cloud droplet size distribution, Geophys. Res. Lett., 33, L02810,, 2006. 

Liu, Y., de Leeuw, G., Kerminen, V.-M., Zhang, J., Zhou, P., Nie, W., Qi, X., Hong, J., Wang, Y., Ding, A., Guo, H., Krüger, O., Kulmala, M., and Petäjä, T.: Analysis of aerosol effects on warm clouds over the Yangtze River Delta from multi-sensor satellite observations, Atmos. Chem. Phys., 17, 5623–5641,, 2017. 

Liu, Y. G. and Daum, P. H.: Anthropogenic aerosols – Indirect warming effect from dispersion forcing, Nature, 419, 580–581, 2002. 

Lohmann, U. and Feichter, J.: Global indirect aerosol effects: a review, Atmos. Chem. Phys., 5, 715–737,, 2005. 

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

Lu, C., Liu, Y., Niu, S., and Vogelmann, A. M.: Observed impacts of vertical velocity on cloud microphysics and implications for aerosol indirect effects, Geophys. Res. Lett., 39, L21808,, 2012. 

Lu, M.-L., Conant, W. C., Jonsson, H. H., Varutbangkul, V., Flagan, R. C., and Seinfeld, J. H.: The Marine Stratus/ Stratocumulus Experiment (MASE): Aerosol-cloud relationships in marine stratocumulus, J. Geophys. Res., 112, D10209,, 2007. 

Lu, M.-L., Sorooshian, A., Jonsson, H. H., Feingold, G., Flagan, R. C., and Seinfeld, J. H.: Marine stratocumulus aerosol-cloud relationships in the MASE-II experiment: Precipitation susceptibility in eastern Pacific marine stratocumulus, J. Geophys. Res., 114, D24203,, 2009. 

Ma, J., Chen, Y., Wang, W., Yan, P., Liu, H., Yang, S., Hu, Z., and Lelieveld, J.: Strong air pollution causes widespread haze-clouds over China, J. Geophys. Res., 115, D18204,, 2010. 

Ma, X., Yu, F., and Quaas, J.: Reassessment of satellite-based estimate of aerosol climate forcing, J. Geophys. Res.-Atmos., 119, 10394–10409, 2014. 

Ma, X., Jia, H., Yu, F., and Quaas, J.: Opposite aerosol index-cloud droplet effective radius correlations over major industrial regions and their adjacent oceans, Geophys. Res. Lett., 45, 5771–5778, 2018. 

Malinowski, S. P., Gerber, H., Jen-La Plante, I., Kopec, M. K., Kumala, W., Nurowska, K., Chuang, P. Y., Khelif, D., and Haman, K. E.: Physics of Stratocumulus Top (POST): turbulent mixing across capping inversion, Atmos. Chem. Phys., 13, 12171–12186,, 2013. 

Martin, G. M., Johnson, D. W., and Spice, A.: The measurement and parameterization of effective radius of droplets in warm stratocumulus clouds, J. Atmos. Sci., 51, 1823–1842, 1994. 

McCoy, D. T., Bender, F. M., Mohrmann, J. K. C., Hartmann, D. L., Wood, R., and Grosvenor, D. P.: The global aerosol-cloud first indirect effect estimated using MODIS, MERRA, and AeroCom, J. Geophys. Res.-Atmos., 122, 1779–1796, 2017. 

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

Nenes, A., Charlson, R. J., Facchini, M. C., Kulmala, M., Laaksonen, A., and Seinfeld, J. H.: Can chemical effects on cloud droplet number rival the first indirect effect?, Geophys. Res. Lett., 29, 1848,, 2002. 

Painemal, D. and Zuidema, P.: Assessment of MODIS cloud effective radius and optical thickness retrievals over the Southeast Pacific with VOCALS-REx in situ measurements, J. Geophys. Res., 116, D24206,, 2011. 

Pandithurai, G., Dipu, S., Prabha, T. V., Maheskumar, R. S., Kulkarni, J. R., and Goswami, B. N.: Aerosol effect on droplet spectral dispersion in warm continental cumuli, J. Geophys. Res., 117, D16202,, 2012. 

Pawlowska, H., Brenguier, J. L., and Burnet, F.: Microphysical properties of stratocumulus clouds, Atmos. Res., 55, 15–33, 2000. 

Pawlowska, H., Grabowski, W. W., and Brenguier, J.-L.: Observations of the width of cloud droplet spectra in stratocumulus, Geophys. Res. Lett., 33, L19810,, 2006. 

Peng, Y. and Lohmann, U.: Sensitivity study of the spectral dispersion of the cloud droplet size distribution on the indirect aerosol effect, Geophys. Res. Lett., 30, 1507,, 2003. 

Peng, Y., Lohmann, U., Leaitch, R., and Kulmala, M.: An investigation into the aerosol dispersion effect through the activation process in marine stratus clouds, J. Geophys. Res., 112, D11117,, 2007. 

Reutter, P., Su, H., Trentmann, J., Simmel, M., Rose, D., Gunthe, S. S., Wernli, H., Andreae, M. O., and Pöschl, U.: Aerosol- and updraft-limited regimes of cloud droplet formation: influence of particle number, size and hygroscopicity on the activation of cloud condensation nuclei (CCN), Atmos. Chem. Phys., 9, 7067–7080,, 2009. 

Saponaro, G., Kolmonen, P., Sogacheva, L., Rodriguez, E., Virtanen, T., and de Leeuw, G.: Estimates of the aerosol indirect effect over the Baltic Sea region derived from 12 years of MODIS observations, Atmos. Chem. Phys., 17, 3133–3143,, 2017. 

Slawinska, J., Grabowski, W. W., Pawlowska, H., and Wyszogrodzki, A. A.: Optical properties of shallow convective clouds diagnosed from a bulk-microphysics large-eddy simulation, J. Climate, 21, 1639–1647, 2008. 

Slawinska, J., Grabowski, W. W., Pawlowska, H., and Morrison, H.: Droplet activation and mixing in large-eddy simulation of a shallow cumulus field, J. Atmos. Sci., 69, 444–462, 2012. 

Stevens, B.: Entrainment in stratocumulus-topped mixed layers, Q. J. Roy. Meteor. Soc., 128, 2663–2690,, 2002 

Su, W., Loeb, N. G., Xu, K.-M., Schuster, G. L., and Eitzen, Z. A. : An estimate of aerosol indirect effect from satellite measurements with concurrent meteorological analysis, J. Geophys. Res., 115, D18219,, 2010. 

Tang, J., Wang, P., Mickley, L. J., Xia, X., Liao, H., Yue, X., Sun, L., and Xia, J.: Positive relationship between liquid cloud droplet effective radius and aerosol optical depth over Eastern China from satellite data, Atmos. Environ., 84, 244–253, 2014. 

Twohy, C. H., Austin, P. H., and Charlson, R. J.: Chemical consequences of the initial diffusional growth of cloud droplets: a clean marine case, Tellus B, 41, 51–60, 1989. 

Twohy, C. H., Anderson, J. R., Toohey, D. W., Andrejczuk, M., Adams, A., Lytle, M., George, R. C., Wood, R., Saide, P., Spak, S., Zuidema, P., and Leon, D.: Impacts of aerosol particles on the microphysical and radiative properties of stratocumulus clouds over the southeast Pacific Ocean, Atmos. Chem. Phys., 13, 2541–2562,, 2013. 

Twomey, S.: Pollution and Planetary Albedo, Atmos. Environ., 8, 1251–1256, 1974. 

Wang, F., Guo, J., Wu, Y., Zhang, X., Deng, M., Li, X., Zhang, J., and Zhao, J.: Satellite observed aerosol-induced variability in warm cloud properties under different meteorological conditions over eastern China, Atmos. Environ., 84, 122–132, 2014. 

Wang, F., Guo, J., Zhang, J., Huang, J., Min, M., Chen, T., Liu, H., Deng, M., and Li, X.: Multi-sensor quantification of aerosol-induced variability in warm clouds over eastern China, Atmos. Environ., 113, 1–9, 2015. 

Wang, Q. and Albrecht, B. A.: Observations of cloud-top entrainment in marine stratocumulus clouds, J. Atmos. Sci., 51, 1530–1547, 1994. 

Warner, J.: The water content of cumuliform cloud, Tellus, 7, 449–457, 1955. 

Warner, J.: The Microstructure of Cumulus Cloud. Part I. General Features of the Droplet Spectrum, J. Atmos. Sci., 26, 1049–1059, 1969a. 

Warner, J.: The microstructure of cumulus clouds. Part II. The effect on droplet size distribution of cloud nucleus spectrum and updraft velocity, J. Atmos. Sci., 26, 1272–1282, 1969b. 

Wood, R.: Stratocumulus clouds, Mon. Weather Rev., 140, 2373–2423,, 2012. 

Wood, R., Bretherton, C., Huebert, B., Mechoso, C. R., and Weller, R.: VOCALS-SouthEast Pacific Regional Experiment (REx) Scientific Program Overview, available at: overview.pdf (last access: January 2011), June 2007. 

Wood, R., Comstock, K. K., Bretherton, C. S., Cornish, C., Tomlinson, J., Collins, D. R., and Fairall, C.: Open cellular structure in marine stratocumulus sheets, J. Geophys. Res., 113, D12207,, 2008. 

Wood, R., Mechoso, C. R., Bretherton, C. S., Weller, R. A., Huebert, B., Straneo, F., Albrecht, B. A., Coe, H., Allen, G., Vaughan, G., Daum, P., Fairall, C., Chand, D., Gallardo Klenner, L., Garreaud, R., Grados, C., Covert, D. S., Bates, T. S., Krejci, R., Russell, L. M., de Szoeke, S., Brewer, A., Yuter, S. E., Springston, S. R., Chaigneau, A., Toniazzo, T., Minnis, P., Palikonda, R., Abel, S. J., Brown, W. O. J., Williams, S., Fochesatto, J., Brioude, J., and Bower, K. N.: The VAMOS Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx): goals, platforms, and field operations, Atmos. Chem. Phys., 11, 627–654,, 2011.  

Yum, S. S., Wang, J., Liu, Y., Senum, G., Springston, S., McGraw, R., and Yeom, J. M.: Cloud microphysical relationships and their implication on entrainment and mixing mechanism for the stratocumulus clouds measured during the VOCALS project, J. Geophys. Res.-Atmos., 120, 5047–5069, 2015. 

Zhang, Q., Quan, J., Tie, X., Huang, M., and Ma, X.: Impact Aerosol Particles on Cloud Formation: Aircraft Measurements in Beijing, China, Atmos. Environ., 45, 665–672, 2011. 

Zhao, C., Tie, X., Brasseur, G., Noone, K. J., Nakajima, T., Zhang, Q., Zhang, R., Huang, M., Duan, Y., Li, G., and Ishizaka, Y.: Aircraft measurements of cloud droplet spectral dispersion and implications for indirect aerosol radiative forcing, Geophys. Res. Lett., 33, L16809,, 2006. 

Zheng, X., Albrecht, B. A., Minnis, P., Ayers, K., and Jonson, H. H.: Observed aerosol and liquid water path relationships in marine stratocumulus, Geophys. Res. Lett., 37, L17803,, 2010. 

Zheng, X., Albrecht, B., Jonsson, H. H., Khelif, D., Feingold, G., Minnis, P., Ayers, K., Chuang, P., Donaher, S., Rossiter, D., Ghate, V., Ruiz-Plancarte, J., and Sun-Mack, S.: Observations of the boundary layer, cloud, and aerosol variability in the southeast Pacific near-coastal marine stratocumulus during VOCALS-REx, Atmos. Chem. Phys., 11, 9943–9959,, 2011. 

Short summary
Aircraft measurements are used to study aerosol–cloud interaction, with a focus on three understudied topics (separation of aerosol effects from dynamic effects, dispersion effects, and entrainment-mixing processes). After constraining cloud dynamics, positive correlation between relative dispersion and CCN concentration became stronger, implying that perturbations of dynamics could underestimate dispersion effect. Entrainment mixing is predominantly extremely inhomogeneous in the stratocumulus.
Final-revised paper