Articles | Volume 22, issue 11
Atmos. Chem. Phys., 22, 7353–7372, 2022
Atmos. Chem. Phys., 22, 7353–7372, 2022
Research article
08 Jun 2022
Research article | 08 Jun 2022

Addressing the difficulties in quantifying droplet number response to aerosol from satellite observations

Addressing the difficulties in quantifying droplet number response to aerosol from satellite observations
Hailing Jia1, Johannes Quaas1, Edward Gryspeerdt2,3, Christoph Böhm4, and Odran Sourdeval5 Hailing Jia et al.
  • 1Leipzig Institute for Meteorology, Universität Leipzig, Leipzig, Germany
  • 2Space and Atmospheric Physics Group, Imperial College London, London, UK
  • 3Grantham Institute for Climate Change and the Environment, Imperial College London, London, UK
  • 4Institute for Geophysics and Meteorology, University of Cologne, Cologne, Germany
  • 5Laboratoire d'Optique Atmosphérique, Université de Lille, CNRS, Lille, France

Correspondence: Hailing Jia (


Aerosol–cloud interaction is the most uncertain component of the overall anthropogenic forcing of the climate, in which cloud droplet number concentration (Nd) sensitivity to aerosol (S) is a key term for the overall estimation. However, satellite-based estimates of S are especially challenging, mainly due to the difficulty in disentangling aerosol effects on Nd from possible confounders. By combining multiple satellite observations and reanalysis, this study investigates the impacts of (a) updraft, (b) precipitation, (c) retrieval errors, and (d) vertical co-location between aerosol and cloud on the assessment of S in the context of marine warm (liquid) clouds. Our analysis suggests that S increases remarkably with both cloud-base height and cloud geometric thickness (proxies for vertical velocity at cloud base), consistent with stronger aerosol–cloud interactions at larger updraft velocity for midlatitude and low-latitude clouds. In turn, introducing the confounding effect of aerosol–precipitation interaction can artificially amplify S by an estimated 21 %, highlighting the necessity of removing precipitating clouds from analyses of S. It is noted that the retrieval biases in aerosol and cloud appear to underestimate S, in which cloud fraction acts as a key modulator, making it practically difficult to balance the accuracies of aerosol–cloud retrievals at aggregate scales (e.g., 1×1 grid). Moreover, we show that using column-integrated sulfate mass concentration (SO4C) to approximate sulfate concentration at cloud base (SO4B) can result in a degradation of correlation with Nd, along with a nearly twofold enhancement of S, mostly attributed to the inability of SO4C to capture the full spatiotemporal variability of SO4B. These findings point to several potential ways forward to practically account for the major influential factors by means of satellite observations and reanalysis, aiming at optimal observational estimates of global radiative forcings due to the Twomey effect and also cloud adjustments.

1 Introduction

Aerosol particles, by acting as cloud condensation nuclei (CCN), can modify cloud properties and precipitation formation, altering the radiative flux at the top of the atmosphere, which is known as effective radiative forcing from aerosol–cloud interactions (ERFaci) (Forster et al.2021). Additionally, absorbing aerosols can also alter the cloud distribution by perturbing the atmospheric temperature structure, known as semi-direct effects (Allen et al.2019). ERFaci may be further subdivided into (i) the radiative forcing due to aerosol–cloud interactions (RFaci), also known as the Twomey effect describing the increased cloud albedo resulting from enhancement in cloud droplet number concentration (Nd) due to an increase in anthropogenic aerosol emissions (Twomey1974), and (ii) rapid adjustments, which are essentially the consequent responses of liquid water path (LWP) and cloud horizontal extent to changed Nd via the Twomey effect (Albrecht1989; Ackerman et al.2004; Zhao and Garrett2015; Bellouin et al.2020). Although extensive investigations have been made to quantify the Twomey effect, significant uncertainties remain regarding its magnitude. This study will discuss the Twomey effect with a focus on the sensitivity of Nd to CCN perturbations due to its fundamental role in aerosol–cloud interactions. Note that the related radiative forcing will not be estimated here, as the anthropogenic perturbation to CCN concentrations is highly uncertain and not easily accessible from observational data.

Current climate models suggest diverse magnitudes of the Twomey effect even with identical anthropogenic aerosol emission perturbation (Gryspeerdt et al.2020; Smith et al.2020). Thus, observational data at the climate-relevant scale, i.e., satellite retrievals, are required to quantify and constrain the Twomey effect globally, which is basically the sensitivity of Nd to CCN perturbations (Seinfeld et al.2016). As reviewed recently by Quaas et al. (2020), there are, however, several uncertainties in inferring the Nd-to-CCN sensitivity (S=dlnNddlnNCCN, where NCCN means proxies for CCN number concentration) from satellite observations, hindering its applicability to further evaluate climate models or quantify RFaci from data. Most of them have been reported to bias S toward a lower value, in turn leading to an overall underestimated ERFaci, including (i) the instrument detectability limitations on aerosol loading in pristine environments (P.-L. Ma et al.2018), (ii) inadequate proxies (such as aerosol optical depth – AOD, or a variant thereof) for CCN owing to the lack of information on the aerosol size and chemical composition (Stier2016; Hasekamp et al.2019), (iii) the limited usability of the AOD–Nd relationship in the present day (PD) to determine the change in Nd caused by anthropogenic aerosol emission due to differing preindustrial (PI) and PD aerosol environments (Penner et al.2011; Gryspeerdt et al.2017), and (iv) satellite sampling biases, which tend to discard clouds with a high cloud fraction due to the inability to retrieve aerosol under cloudy conditions, thereby resulting in an artificial cloud regime selection (i.e., omitting more retrieval-reliable stratiform clouds; Gryspeerdt and Stier2012; Jia et al.2021). Additionally, meteorological conditions, e.g., lower tropospheric stability (P.-L. Ma et al.2018), relative humidity (Quaas et al.2010), availability of water vapor (Qiu et al.2017), and wind shear (Fan et al.2009), and vertical overlapping status of aerosol and cloud layers (Costantino and Bréon2013; Zhao et al.2019) also play roles in regulating aerosol–cloud interactions. It is worth noting that most of these studies calculated S based on cloud effective radius rather than Nd, so they are subject to even more errors from the problem of stratification by liquid water path. Currently, a key difficulty in interpreting satellite-observed aerosol–Nd relationships is to isolate the causal impact of aerosols on Nd from other confounding factors modifying the variations of aerosol and cloud simultaneously, specifically (i) updraft, determining cloud development as well as the maximum supersaturation at cloud base and thus aerosol population that can be activated, (ii) precipitation processes, depleting cloud droplets via coagulation and scavenging sub-cloud aerosol particles, and (iii) retrieval errors, biasing retrieved aerosol and cloud properties concurrently. However, a clear understanding of how they affect the estimates of S quantitatively is lacking from the perspective of satellite observations (Quaas et al.2020).

In terms of the updraft, in situ aircraft measurements (Berg et al.2011; Jia et al.2019b), ground-based remote sensing (Schmidt et al.2015; McComiskey et al.2009), and detailed parcel model simulations (Reutter et al.2009; Chen et al.2016) clearly showed the dependency of S on updraft, with generally larger S at stronger updraft. In particular, covariability of updrafts and aerosol concentrations has been found to result in a stronger S than keeping vertical velocity (w) constant (Bougiatioti et al.2020; Kacarab et al.2020). As noted by Gryspeerdt et al. (2017), the updraft may roughly explain 20 % of the variability in ΔNd from its PI–PD difference, adding to the uncertainty of the ERFaci estimate. Despite of the importance of dynamical constraint, it is not easily applicable to the analysis of satellite data due to the lack of updraft observations near cloud base at a global scale. As an alternative, cloud-base height (CBH) may potentially serve as a practical proxy for the updraft at the base of liquid cloud because of the tightly linear correlation illustrated by in situ observations of cumuliform clouds (Zheng and Rosenfeld2015). Although data used to draw this conclusion by Zheng and Rosenfeld (2015) were collected from only three locations, they covered various boundary conditions over both continents and oceans. Moreover, a theoretical framework has also been established to support the observed empirical relationship (Zheng2019), lending credibility to applying CBH as a proxy for the updraft. Building on this, recently developed CBH retrievals (Mülmenstädt et al.2018; Böhm et al.2019) offer an opportunity to gain some insight into the potential role updraft variability may play in the global ERFaci assessment.

In addition to the updraft, precipitation formation further complicates the derivation of the strength of S, since it can efficiently deplete cloud droplets and scavenge aerosols from clouds (Gryspeerdt et al.2015). In such cases, the change in Nd is not necessarily related to actual aerosol perturbations (Chen et al.2014) but rather to the intensity of cloud sink, and thus in principle, it should not be directly applied to infer ΔNd driven by anthropogenic emissions. However, due to the lack of simultaneous observations of precipitation and aerosol–cloud properties from passive satellite remote sensing alone, most aerosol–cloud interaction (ACI) estimates do not consider the influence of precipitation (Quaas et al.2008; Ma et al.2014; Gryspeerdt et al.2017; Jia et al.2021) or just roughly identify the occurrence of rain by relying on some simplified metrics, such as the threshold of 14 µm cloud effective radius (CER) for rain initiation (Gerber1996; Rosenfeld et al.2019; Yang et al.2021; Zhang et al.2022) or the difference of CER between retrievals employing the bands of 2.1 and 3.7 µm (Saponaro et al.2017; Jia et al.2019a). Even though few studies have explicitly accounted for this by combining simultaneous precipitation observations from active remote sensing (Chen et al.2014), how different treatments could influence the assessment of S remains unclear. Solving this problem is helpful to reconcile the current diverse ACI estimates in order to achieve a more confident observational constraint.

For satellite-based investigations, it is crucial but difficult to disentangle any physically meaningful attributable factors from artificial aerosol–cloud linkage induced by retrieval biases. In terms of Nd, retrievals for 3D-shaped clouds and partially cloudy pixels deviate from the retrieval assumptions of overcast homogenous cloud and 1D plane-parallel radiative transfer, thereby appearing to lead to an overestimation of CER (Coakley et al.2005; Matheson et al.2006; Zhang and Platnick2011; Zhao et al.2012), and in turn, an underestimated Nd (Grosvenor et al.2018). This issue was reported to be more pronounced for broken cloud regimes and could to some extent be addressed by only sampling Nd for pixels with either a high cloud fraction (Painemal et al.2020) or large cloud optical depth (COT; Zhu et al.2018). In addition to the assumptions for clouds, the existence of aerosols above clouds can also affect the retrieval of cloud optical depth (Haywood et al.2004; Li et al.2014), in turn bias Nd calculation. Meanwhile, the retrieved AOD or aerosol index (AI) can be biased to a larger value due to inability to detect thin clouds in an aerosol retrieval scene (Kaufman et al.2005) or due to enhanced reflectance from neighboring clouds (Várnai and Marshak2009). It is noteworthy that the overestimation of AOD tends to be enhanced with increasing cloud fraction (Zhang et al.2005) and COT (Várnai and Marshak2021) as a result of both retrieval problems and aerosol swelling (Quaas et al.2010). Therefore, the potential covariations between biases in Nd and AOD (AI) modulated by cloud macrophysical properties could incur a spurious correlation between the two variables, obscuring the causal interpretation. While a few studies pointed out that the AOD(AI)–Nd correlation is substantially enhanced when analyzing reliable Nd retrievals (Jia et al.2019a; Painemal et al.2020), how and to what extent the satellite-diagnosed S varies with the retrieval biases in terms of both aerosol and Nd, respectively, has not been fully understood. Such understanding is quite important for reconciling the previous estimates and proposing a meaningful method applicable to satellite-based investigations.

While the problem of vertical co-location between retrieved CCN proxies and clouds has been noticed in many previous studies, most of them placed focus on its influence on the correlation between aerosol and cloud (Stier2016; Painemal et al.2020), i.e., a much higher correlation between Nd and aerosol extinction coefficients near cloud base compared to Nd vs. column-integrated aerosol quantity (AOD or AI), rather than the influence on S. The latter is usually quantified as a regression coefficient (regression slope in log–log space) between Nd and the CCN proxy and is a key determinant of radiative forcing estimates. Using AI as a CCN proxy, Costantino and Bréon (2010) demonstrated a weaker cloud susceptibility for the case with separated aerosol–cloud layers than well-mixed ones. However, it is unclear how the S would change when switching commonly used column aerosol quantities to aerosol measures at cloud base. This understanding is particularly important for the intercomparison and further reconciliation between current ACI metrics relying on diverse CCN proxies, including column-integrated, near-surface, and cloud-level aerosol quantities.

In this study, we focus on the quantification of the impacts of three major confounders mentioned above, namely updraft, precipitation, and retrieval errors, as well as the problem of vertical co-location between aerosol and cloud, on the assessment of S in the context of marine warm clouds by combining multiple active and passive satellite sensors and reanalysis products. On the basis of current findings, this study further suggests several potential ways forward to practically account for, to the extent possible, the major influencing factors for the satellite-based quantification of S and hence the ERFaci.

2 Data and methods

This work is based on observational data from multiple instruments on board Terra, Aqua, and CloudSat platforms as well as reanalysis data from the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) (Randles et al.2017) and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) (Hersbach et al.2020). Table 1 summarizes the aerosol, cloud, and precipitation parameters and their corresponding sources, temporal–spatial resolutions, and time periods analyzed in the present study. Note that due to the requirement for co-located aerosol–cloud–precipitation observations, the data used in Sect. 3.2 are obtained from the A-Train constellation of satellites (Aqua and CloudSat), which are then interpolated to 5×5 km2 resolution for analysis, while the remaining parts are based on the observations from Terra, for which all data are interpolated to 1×1 resolution. The combination of datasets used in each section is summarized in Table 2. It is worth mentioning that, as S was found to vary with the spatial resolution of data (Sekiguchi et al.2003; McComiskey and Feingold2012), the different data resolutions between Sect. 3.2 and other sections can lead to a difference in S, but this is not the focus here. This study is restricted to the global ocean with latitude between 60 S and 60 N because of limited quality of retrievals of aerosol size parameters (Levy et al.2013) and Nd (Gryspeerdt et al.2021) over land and polar regions.

Table 1The list of the parameters, sources, and their corresponding temporal–spatial resolutions applied in the present study.

Download Print Version | Download XLSX

Table 2The combination of datasets used in each subsection of the Results section.

Download Print Version | Download XLSX

Aerosol properties (Levy et al.2013) are obtained from the level 3 Moderate Resolution Imaging Spectroradiometer (MODIS) Dark Target product (MOD08 and MYD08; Platnick et al.2017c, d). In order to collect co-located (adjacent) aerosol and cloud retrievals for analysis, aerosol retrievals on a coarsely resolved grid (1× 1 on a latitude–longitude grid) are used to match cloud pixels (1×1 km2), assuming that aerosol properties in adjacent clear areas are homogeneous enough to represent those under cloudy conditions (Anderson et al.2003; Quaas et al.2008). Note that this assumption would be questionable, especially when aerosol is scavenged by precipitation (Gryspeerdt et al.2015). In addition to commonly used AOD, the aerosol index (AI = AOD × Ångström exponent) containing information on aerosol size is also employed since it is considered a better proxy for CCN (Nakajima et al.2001). The Ångström exponent is calculated from AOD at wavelengths of 460 and 660 nm. To eliminate 1 by 1 scenes in which the aerosol distribution is heterogeneous, retrievals with a standard deviation higher than the mean values are discarded (Saponaro et al.2017). As suggested by Hasekamp et al. (2019), the lowest 15 % of data for AOD (AI) at a global scale are excluded to avoid large retrieval uncertainty at low aerosol concentrations (P.-L. Ma et al.2018). Note that leaving out the low AOD (AI) yields a larger S compared to using all data (Hasekamp et al.2019).

Cloud optical properties, including CER and COT at 3.7 µm (Platnick et al.2017e), are obtained from the MODIS level 2 cloud products (MOD06 and MYD06; Platnick et al.2017a, b) and then applied to compute Nd based on the adiabatic approximation (Quaas et al.2006). It was found that the filtering of cloud adiabaticity only has a negligible impact on the estimate of S, but in turn results in a reduction of up to 63 % in the data volume (Gryspeerdt et al.2021). For this reason, we do not apply such filtering here. Note that Nd is calculated on the level of the satellite pixel (order 1 km) before being aggregated to larger scales. Thus, the aggregation bias caused by the derivation of Nd from the highly nonlinear function of CER and COT as shown by Feingold et al. (2022) does not affect the results presented here. To ensure confident retrievals, the Nd is filtered to include only single-layer liquid clouds with top temperature higher than 268 K. Pixels for which CER<4µm and COT<4 are discarded due to the large uncertainty of retrievals (Sourdeval et al.2016). In addition, only pixels with a cloud fraction at 5 km resolution (CF5×5km2) > 0.9 and with a sub-pixel inhomogeneity index (cloud_mask_SPI) < 30 are used to reduce the retrieval errors induced by cloud edges and broken clouds (Zhang and Platnick2011). Further, we only consider pixels with a solar zenith angle of less than 65 and a sensor zenith angle of less than 41.4 to minimize the influence of known biases as detailed in Grosvenor et al. (2018). With the above sampling strategy, the random uncertainty in Nd was reported at 78 % on a pixel level, and this dropped substantially when averaged to a 1 by 1 region (Grosvenor et al.2018). However, as stated in Gryspeerdt et al. (2021), the systematic bias in the Nd retrievals for in situ measurements is low, with determination coefficients of 0.48 for all cloud types and 0.5–0.8 for stratocumulus clouds.

To overcome the lack of global updraft observations, we utilize satellite-based retrievals for CBH as a proxy for cloud-base updraft for cumuliform clouds based on the finding that these two quantities exhibit an approximately linear correlation for convective clouds (Zheng and Rosenfeld2015). Here, clouds are considered convective for lower troposphere stability (LTS) less than 16 K (Rosenfeld et al.2019). Additionally, cloud geometrical thickness (CGT; the difference between cloud-top height and CBH) is used as an alternative proxy for the updraft regardless cloud types, since it has been observed to be associated with the cloud-base updraft for shallow cumuliform clouds (Lareau et al.2018) and also correlated with cloud-base updraft for stratiform clouds via modulating cloud-top cooling (Zheng et al.2016). To obtain CBH and CGT, we apply a recently developed retrieval algorithm (0.25× 0.25 resolution, Böhm et al.2019) based on Multi-angle Imaging SpectroRadiometer (MISR)/Terra observations, i.e., the MISR Level 2 Cloud Product (MIL2TCSP; NASA/LARC/SD/ASDC2012). The best performance of this algorithm is achieved for clouds with CBH around 1 km and CGT below 1 km. For such heights, which are characteristic for oceanic clouds considered in this analysis, the root mean square error ranges 300–350 m. It is important to note that the MISR cloud-base height retrieval is limited to CBH > 560 m (Böhm et al.2019). At this lower end of the detection range, a slight underestimation of the CBH is expected (Böhm et al.2019). The ERA5 reanalysis is employed here to calculate LTS as the difference in potential temperature between 700 and 1000 hPa (Klein and Hartmann1993). The hourly LTS is then matched to 10:30 local solar time to approximate the overpass time of the Terra satellite.


Figure 1Schematic diagram of the procedure for calculating the sensitivity (linear regression coefficient in log–log space) of Nd to CCN; AI is taken as a example. The upper panel shows the global joint Nd–AI histogram, wherein each column is normalized so that it sums to 1. The blue line is a linear regression on the 20 paired Nd–AI (blue dots) that are the medians of each AI bin with an equal number of samples, and the yellow dashed line shows a linear regression on all data points. Note that the lowest 15 % of AI values have been left out according to occurrence (bottom) before binning the data. The cloud susceptibilities to AI (SAI) derived from both approaches are shown along with 95 % uncertainty estimates (according to a Student's t test). The data used here are the same as in Sect. 3.1.


To identify the role of precipitation, CloudSat radar precipitation observations co-located with AOD and AI as well as Nd from MODIS/Aqua are adopted as well. Here, we use the precipitation flag from the 2B-CLDCLASS product (Sassen and Wang2008) to distinguish precipitating (with the flags of “liquid precipitation” and “possible drizzle”) and non-precipitation clouds (with the flag of “no precipitation”). As a sink of Nd, drizzle could also affect the aerosol–cloud interactions even without rain falling on ground (Yang et al.2021), so we also include drizzling clouds in precipitating cases. The CloudSat data at a 1.4×2.5 km2 resolution are matched to the nearest MYD06 5×5 km2 pixels for further analyses.

Figure 2Dependence of the linear regression slopes of ln Nd versus ln AOD (blue) and ln AI (green) on (a) CBH and (b) CGT derived via the pre-binned approach. Data are grouped into 10 fixed CBH (CGT) intervals for regressions. Error bars indicate the 95 % confidence interval of the linear regression, and the gray bars denote the total number of samples for each CBH (CGT) bin. The corresponding regression slopes computed from the data over all CBH (CGT) bins are shown as horizontal dashed lines (green for AI and blue for AOD). The equivalent figure (Fig. S1) shows similar results based on the all-data approach.


The MERRA-2 product assimilates observations of the atmospheric state as well as remotely sensed AOD so that it can generate reasonable aerosol horizontal and vertical distributions (Buchard et al.2017). The use of aerosol reanalysis also largely avoids the spuriously high AOD near clouds caused by retrieval artifacts from the satellite (Jia et al.2021). Given that variability in sulfate aerosols contributes the most strongly to variability in Nd among all aerosol species (McCoy et al.2017), the sulfate concentration is considered be to the CCN proxy here. We utilize vertically resolved sulfate mass concentrations from MERRA-2 reanalysis in combination with the MISR CBH retrieval to obtain sulfate mass concentrations near cloud base (SO4B). In addition, sulfate surface mass concentrations (SO4S) and column mass density (SO4C) are also used to investigate if there will be different behaviors of Nd-to-CCN sensitivity when applying CCN proxies at different levels. The MERRA-2 3 h averaged fields are interpolated to 10:30 local solar time to approximate the overpass time of the Terra satellite.

Figure 1 illustrates the regression procedure for calculating the S. After excluding the lowest 15 % AOD (AI), the data are then divided into 20 bins of CCN proxy, with each bin having an equal number of samples. The same number of samples ensures the same statistical representativeness within each bin. The values of Nd and the CCN proxy for a certain bin are the medians of all values in that bin. The generated 20 paired values of Nd and CCN proxies are then used in linear regression to determine S unless otherwise stated. The uncertainties of estimated S are reflected by the 95 % confidence interval of the regression slope. We also tried 100 and 1000 bins and found that the derived susceptibilities do not change significantly with the number of bins. Additionally, the linear regression on all data points is also shown (yellow dashed line) in Fig. 1 for comparison with the pre-binned approach, since both approaches have been used extensively by previous studies (Quaas et al.2008; Gryspeerdt et al.2017; Hasekamp et al.2019; Rosenfeld et al.2019), but it is still unclear how large the difference in estimates between the two approaches could be. Figure 1 illustrates that the pre-binned approach has a larger slope than lumping together all data points by 18 %, suggesting that attention should be paid when comparing S derived from different approaches. In our study, both approaches lead to similar conclusions; as such, we will only focus on the results from the pre-binned approach in the main text. Meanwhile, we also put the results associated with the all-data approach in the Supplement.

3 Results

3.1 Dependence on updraft

In adiabatic clouds, Nd is essentially a function of both CCN and updraft (Feingold et al.2001). To quantify how Nd responds to CCN perturbations, the variation of updraft must be constrained. In practical terms, however, the observation of in-cloud vertical velocity is possible only from in situ aircraft measurements or ground-based remote sensing, limiting the estimations to individual locations and sites. In order to obtain S at a global scale, which is only possible from a satellite, meteorological parameters (X. Ma et al.2018) or cloud regimes (Gryspeerdt and Stier2012) were generally employed to roughly approximate cloud dynamics. However, it should be noted that even in similar meteorological backgrounds and cloud regimes, the vertical velocity within individual clouds can still vary significantly (Hudson and Noble2014). Instead, based on previous findings from in situ observations (see the section “Data and methods”), our study utilizes CBH as a proxy for cloud-base updraft for cumuliform clouds and CGT as a proxy for the updraft regardless of cloud type. Note that with similar cloud-top heights, the higher cloud base means thinner cloud layer. To avoid the potential interference by CGT, the analysis of the dependence of S on CBH (Fig. 2a) is conducted within a quasi-constant CGT bin of 650–750 m. This range is chosen because of its relatively strong S, low possibility of precipitation, and sufficient data points (Fig. 2b).

Figure 3Joint histograms between AI and Nd (CER) created for weak and strong updraft conditions, as defined by the lowest and highest CGT quartiles, respectively. The difference plots between strong and weak cases are shown at the end of each row. The histograms are normalized so each column sums to 1 such that the histograms show the probability of observing a specific Nd (CER) given a certain AI.


Figure 2 shows the dependence of linear regression slopes of ln Nd versus ln AOD (ln AI), i.e., SAOD (SAI), on CBH and CGT, respectively. To constrain the variation of cloud dynamics, the data are grouped over CBH and CGT bins with intervals of 80 and 100 m, respectively. It is seen that SAOD and SAI exhibit increases with both CBH and CGT, consistent with the expectation of stronger aerosol–cloud interactions under larger in-cloud vertical velocity conditions. The result is in accord with previous findings based on surface remote sensing under stratus (McComiskey et al.2009) and altocumulus clouds (Schmidt et al.2015). Also, using ground-based observations, Feingold et al. (2003) quantified this linkage and gave a correlation of 0.67 between S and column maximum updraft. In our study, the correlation coefficients are 0.83 (0.98) for CBH–SAOD (SAI) and 0.96 (0.95) for CGT–SAOD (SAI). The higher correlations likely stem from the large volume of data used to stratify CBH(CGT), which enhances the representability of samples from a statistical perspective compared to the more limited number of cases used in Feingold et al. (2003).

Figure 4Joint histograms between AI and Nd created for (a) all clouds, (b) non-raining, and (c) raining clouds, as well as (d) the difference of joint histograms between the raining and non-raining cases. Cloud susceptibilities to AI derived via the pre-binned approach are also shown along with 95 % uncertainty estimates (according to a Student's t test). The fitting lines for three cases are merged into one single plot (e), with clean and polluted zones marked as blue and red, and the corresponding sample distributions are also shown (f, g).


It is also noted that, unlike the monotonic increase with CBH, SAOD (SAI) increases with CGT at a small to moderate CGT range (<900 m) and then levels off (Fig. 2b). This is likely due to the tighter linkage between the occurrence of precipitation and CGT than CBH. Specifically, larger CGT is an indicator of strong updraft, tending to generate larger SAOD (SAI), whereas in the meantime it is also associated with the higher possibility of precipitation, which acts as an efficient sink of droplets (see Sect. 3.2) and thereby partly offsets the increase in Nd induced by CCN, i.e., smaller SAOD (SAI). In short, the situation of SAOD (SAI) at larger CGT (Fig. 2b) is a result of the competition between the effects of updraft and precipitation. Comparing the different CCN proxies, we see that, in agreement with previous results (Hasekamp et al.2019), SAI is consistently higher than SAOD for both the all-data cases (dashed lines) and almost all CBH (CGT) bins except for CGT > 900 m. For the remainder of the paper, only AI that is a better CCN proxy is used unless otherwise stated.

To gain insight into the mechanism underlying the apparent dependence of S on updraft, we contrast AI–Nd (CER) joint histograms for weak and strong updraft conditions (Fig. 3). As the data volume for the CBH case is too small to populate the joint histogram, only the CGT-related result is shown. Here, the subsets of data with CGT lower than the 25th percentile and higher than the 75th percentile are defined as weak and strong updrafts, respectively. Note that applying the 10th and 90th percentiles also yields similar results as shown in Fig. S2. It is known that the aerosol–Nd relationship is nonlinear and particularly regime-dependent. Reutter et al. (2009) proposed three distinct regimes according to the ratio of vertical velocity and aerosol concentration: (a) an aerosol-limited regime characterized by a high ratio value, nearly linear dependence of Nd on aerosol, and insensitivity of Nd to updraft, (b) an updraft-limited regime characterized by a low ratio value and weak dependence of Nd on aerosol but quite strong dependence on updraft, and (c) a transitional regime falling between the above two regimes. Since we have limited the proxy for updraft (CGT) to a certain range, AI is thus assumed to be an indicator of regime. Specifically, the low AI zone is more likely aerosol-limited, while the high AI zone is close to an updraft-limited regime. As illustrated in the difference plots in Fig. 3, under polluted conditions with AI>0.4, the samples of the strong updraft case tend to concentrate in the larger Nd bins compared to the weak updraft (Fig. 3c), reflecting the critical role of updraft in facilitating activation of cloud droplets. Nevertheless, the distributions of CER do not exhibit systematic differences, except for less scattering for the strong updraft (Fig. 3f). As for clean conditions, what should be expected is a similar distribution of Nd between different cloud dynamics as determined by the nature of the aerosol-limited regime, or at least a slightly higher Nd for the strong updraft case. However, looking at the clean zone (AI<0.15) in Fig. 3, it is clear that strong updraft is associated with much lower Nd as well as larger CER (generally larger than 14 µm, the threshold for drizzle initiation suggested by Freud and Rosenfeld2012) compared to the weak updraft, indicating a higher possibility of precipitation and/or drizzle. Consequently, the strong sink of droplets via precipitation at low AI and the enhanced activation of droplets at high AI will jointly create a much larger regression slope of ln Nd versus ln AI for the strong updraft compared to the weak updraft condition. Moreover, these results also imply that the interference of precipitation tends to amplify realistic dependence of SAI on the updraft, highlighting the need to remove the influence of precipitation on the Nd budget.

Figure 5Dependence of the linear regression slopes of ln Nd (ln NdAll) versus ln AOD (ln AI) on (a) ΔL and (b) ΔNd derived via the pre-binned approach. Data are grouped into 10 fixed ΔL (ΔNd) intervals for the calculation of slopes. Error bars indicate the 95 % confidence interval of the linear regression, and the gray bars denote the total number of samples for each bin. The change in AOD (AI) with ΔL is also shown in the panel (a). The equivalent figure (Fig. S3) shows similar results based on the all-data approach.


3.2 Dependence on precipitation

In this section, the role of precipitation in the quantification of S will be explicitly accounted for by using the simultaneous aerosol–cloud–precipitation observations from CloudSat–MODIS combined datasets (see Sect. 2). The hypothesis is that for precipitating clouds, a sink to Nd exists (via the coagulation) that does not reflect the Twomey effect, so the CCN–Nd relationship is biased low in cases of precipitation formation. Figure 4 shows the AI–Nd joint histograms for non-raining, raining, and all clouds as well as the difference between non-raining and raining cases. The raining clouds exhibit a lower Nd relative to non-raining clouds over all AI bins, caused by the intensive sink of cloud droplets by collision–coalescence when precipitation forms (Fig. 4b, c, d). In addition, as the droplet sink and aerosol removal by precipitation can act together to veil the actual effect of aerosol on Nd, and the Nd in raining clouds shows a weaker response to increasing AI than that in non-raining clouds, with the corresponding SAI of 0.45 versus 0.56, respectively. The result is in agreement with Chen et al. (2014), who reported a consistently smaller CER-to-AI sensitivity in the precipitating case than in the non-precipitating case throughout different environmental conditions.

Interestingly, the regression slope of ln AI versus ln Nd is enhanced after lumping all cloud scenes together regardless of whether it rains or not (Fig. 4a). The corresponding SAI (0.68) increases by 21 % relative to the non-raining case (0.56). This phenomenon was also noted by Painemal et al. (2020), and they speculated that drizzle appears to strengthen the aerosol–Nd relationship, which is, however, contrary to the weaker SAI for raining clouds as illustrated above. For a clearer comparison of the SAI for non-raining, raining, and all clouds, the fitting lines for these three cases are put into one single plot (Fig. 4e), with clean and polluted zones marked as blue and red, and the corresponding sample distributions are presented in Fig. 4f and g. It is shown that the fitting line for all clouds nearly coincides with that for the non-raining case under polluted conditions but is closer to the raining case under clean conditions (Fig. 4e), consequently leading to a much steeper slope. This behavior is further corroborated by the different distributions of Nd. As shown in Fig. 4g, the polluted clouds consist predominately of the non-raining clouds as a result of the suppression of precipitation by aerosols, thus maintaining a high value of Nd. Instead, the majority of the clean clouds are raining ones that are significantly subjected to the sink processes for Nd and/or aerosol scavenging (Boucher and Quaas2013) (Fig. 4f), hence corresponding to a lower Nd. The results presented here imply that introducing the dependence of the possibility of precipitation on aerosols (i.e., cloud lifetime effect) into the estimation of the Twomey effect, as commonly done in most previous studies, would perturb the statistical analysis and artificially bias the strength of the Twomey effect to a higher value. Moreover, it should be noted that a more extensive zone with Nd being insensitive to aerosol is evident under low aerosol conditions after raining clouds are included (Fig. 4a), which means that, in addition to the overestimation of regression slope, the interference of precipitation also gives rise to an apparent nonlinearity of the aerosol–Nd relationship, hence adding substantial complexity in quantifying S using a linear regression (Gryspeerdt et al.2017).

3.3 Dependence on retrieval biases in AOD (AI) and Nd

Aerosol retrieval errors due to 3D radiative effects as well as cloud contamination, aerosol swelling, and cloud retrieval errors for 3D-shaped and heterogeneous clouds have been shown to artificially introduce biases in the estimation of aerosol–cloud interactions (Quaas et al.2010; Christensen et al.2017; Neubauer et al.2017; Jia et al.2019a, 2021). Here, we dig deeper on S as a function of retrieval errors by defining two metrics that characterize the retrieval biases quantitatively. In order to obtain horizontally “co-located” aerosol–cloud retrievals for analysis, the often adopted choice is a 1 by 1 gridding scale, at which aerosol concentrations are considered homogeneous (Anderson et al.2003). Within a 1 by 1 grid box, sub-grid clear-sky and cloudy pixels co-exist (if clouds are not fully overcast) and are used for retrieving cloud and aerosol properties, respectively. However, in the case that most clear-sky pixels are close to clouds, the problems of 3D radiative effects, cloud contamination, and aerosol swelling arise. Thus, the metric of aerosol retrieval errors (including 3D radiative effects and cloud contamination) and aerosol swelling is defined as the average distance to nearest cloudy pixel from clear pixels for aerosol retrieval (ΔL), which is provided directly by the MODIS L3 aerosol product. As for the cloud retrieval, the metric is the difference between Nd retrieved from all cloudy sub-pixels (NdAll, without the cloud screening on CER, COT, CF5×5km2, and the sub-pixel inhomogeneity index) and that retrieved from sub-pixels only with favorable situations for reliable cloud retrieval (see “Data and methods” for details), which is tightly related to the degree of cloud heterogeneity. Note that NdAll and Nd are concurrently calculated for each 1 by 1 cloud scene, and thus ΔNd (NdAll-Nd) only reflects the role of retrieval errors, with other conditions held constant (e.g., cloud types and meteorology). Generally, a negative value of ΔNd is expected since a positive bias in CER and a negative bias in COT for spatially inhomogeneous scenes act together to generate negatively biased NdAll according to the Eq. (1) in Quaas et al. (2006). In this section, we also look at AOD in addition to AI, since AOD is a directly retrieved quantity and thus more closely related to retrieval problems.

Figure 5a shows the dependences of both AOD (AI) and linear regression slopes of ln Nd (NdAll) versus ln AOD (ln AI) on ΔL. We note that AOD (AI) is the largest for the first ΔL bin with a value of 0.24 (0.17) and then drops rapidly to around 0.16 (0.13) for the other distances from clouds, indicating a quite strong near-cloud enhancement of AOD (AI) induced by retrieval biases and/or aerosol swelling. The Ångström exponent (AE) is calculated from AOD at wavelengths of 460 and 660 nm. As the AE was found to increase with ΔL (Várnai and Marshak2015), the reduction of AI with ΔL is thus less strong than AOD. Based on published in situ aircraft measurements, we roughly isolate the contribution of aerosol swelling from retrieval issues (i.e., 3D radiative effects and cloud contamination). During the Indian Ocean Experiment (INDOEX), Twohy et al. (2009) measured a rise in relative humidity (RH) from about 70 % more than 20 km from cloud to 90 % with 1–4 km of cloud edge (equivalent to the distances of the third and first ΔL bins in Fig. 5a), which in turn results in about a 69 % increase in the aerosol scattering cross section (Twohy et al.2009). Considering that aerosol humidification only occurs near cloud level, i.e., one-quarter to one-third of the aerosol column could be affected according to lidar observations (Twohy et al.2009), the increase in AOD by aerosol swelling is estimated to be 17 %–23 %. This is up to about a third of the relative increase in AOD from the third to first ΔL bins (64 %) in Fig. 5a, implying that the retrieval errors in aerosol could contribute the majority of the S reduction in the first ΔL bin. It should be noted that the estimated AOD rise due to humidification relies on observed RH variability surrounding cloud and also the vertical profile and chemical composition of aerosol, which could vary with geographic location.

Figure 6Relationships between ΔL and ΔNd, with the data grouped as a function of CF and each CF bin containing the same number of samples. A joint histogram between CF and CGT is shown in the inner plot; the blue dot shows the median CGT at each CF bin.


Corresponding to the biased-high AOD (AI), SAOD and SAI for the first ΔL bin are very low relative to other bins, especially for AOD, suggesting that both retrieval biases and aerosol swelling near clouds could result in a severe underestimation in S. These results imply that screening out the aerosol retrievals within the first ΔL bin (i.e., the average distance to the nearest cloud pixel less than 10 km) could be an applicable approach to side-step the interference of aerosol retrieval biases. It is also noted that SAOD (SAI) shows an increase first and then a decrease from the third ΔL bin. However, the following decrease is unlikely linked to the aerosol retrieval bias since the AOD (AI) remains almost constant (the upper panel in Fig. 5a). One interpretation for this would be that AOD and/or AI is getting less representative for the aerosol concentrations near cloud with increasing ΔL, especially for grid boxes with precipitation wherein aerosol is not as homogeneous as assumed (Anderson et al.2003). Moreover, as ΔL is also negatively correlated with CF (Várnai and Marshak2015), the decreasing SAOD (SAI) is probably associated with other factors modulated by CF (such as retrieval error in Nd as demonstrated in the following analysis).

Figure 7Two-dimensional probability density functions of ln Nd versus (a) ln SO4B, (b) ln SO4S, and (c) ln SO4C for the period 2006–2009. Sample numbers (N), correlation coefficients, and regression slopes with 95 % uncertainty estimates (according to Student's t test) for pre-binned SO4-Nd pairs are displayed in the upper left of each plot.


Interestingly, Fig. 5a also depicts the SAOD (SAI) calculated from NdAll as consistently lower than that from Nd for each ΔL bin, indicating that the cloud retrieval biases for partly cloudy pixels appear to lead to an underestimation of S. The increase in the difference between them with ΔL reveals that more underestimation occurs for high ΔL (typically low CF) conditions in which clouds are more partially cloudy, thereby deviating from the retrieval assumptions of overcast homogeneous cloud. As previously mentioned, ΔNd can act as a measure of some of the retrieval errors in cloud; the more negative the ΔNd, the larger the retrieval error in Nd. As shown in Fig. 5b, the SAOD (SAI) calculated from NdAll increases with ΔNd and then reaches its maximum when ΔNd approaches 0, demonstrating that the satellite-diagnosed S highly depends on the retrieval bias in cloud. In terms of the quality-assured Nd, the corresponding SAOD (SAI) is not anticipated to be affected by retrieval issues and is thus independent of ΔNd, but this is obviously not the case; the SAOD (SAI) also significantly increases with ΔNd, which means that the criteria used for selecting homogeneous clouds within a 5 km × 5 km grid would not be as sufficient for an optimal performance of retrieval (Grosvenor et al.2018) as we thought.

Figure 6 depicts relationships between ΔL and ΔNd, with the data grouped as a function of CF for 50 cloud fraction bins containing the same number of samples. It is clearly illustrated that CF regulates the negative correlation between ΔL and ΔNd. Under the condition of large CF, clear pixels are very close to the nearest cloud pixel, corresponding to a lower ΔL; meanwhile, most sub-grid cloud pixels meet the criteria for confident cloud retrievals, leading to a higher (near-zero) ΔNd. The reverse is true in the case of low CF. This means that it is practically difficult to balance the accuracies of retrievals on both aerosol and cloud, since the aerosol retrieval should stay away from clouds, requiring low CF, whereas the Nd retrieval should be performed in more homogeneous clouds (high CF) in order to satisfy the retrieval assumption of 1D plane-parallel radiative transfer. To avoid the spuriously high AOD (AI) retrieval near clouds, the use of aerosol reanalysis would be a way forward (Jia et al.2021). In terms of Nd, however, the situation is more complicated. Given that CF also correlates closely with cloud dynamics (CGT; Fig. 6), it does not make sense to simply restrict the analysis to low ΔNd (thus high CF) to reduce the retrieval uncertainty of Nd; in doing so, a selection of cloud regime could be artificially applied.

3.4 Dependence on vertical co-location between aerosol and cloud

Currently, the use of reanalyzed and/or modeled aerosol vertical profiles seems to be the only feasible alternative to exploit the problem of vertical co-location since it is impossible to obtain aerosol retrievals below or within clouds from satellites (Stier2016; McCoy et al.2017). Thus, unlike the previous sections based on satellite-retrieved AOD and AI, vertically resolved SO4 from the MERRA-2 reanalysis is utilized here to obtain the CCN proxies for different altitudes. Although not as commonly adopted as AOD and AI, SO4C and SO4S were also used as CCN proxies by previous studies (McCoy et al.2017; Jia et al.2021). Here, the SO4C and SO4S are used, respectively, to mimic the behaviors of AOD (AI) and the surface aerosol extinction coefficient, which are two commonly used CCN proxies in satellite-based and ground-based methods (Quaas et al.2008; Liu and Li2018), respectively. As demonstrated by Stier (2016), the SO4B derived in combination with CBH is expected to be more relevant to the amount of CCN actually activated at cloud base than SO4C and SO4S. The comparison of susceptibilities inferred from these three proxies helps us to understand whether the uses of column-integrated and near-surface aerosol quantities make sense and, more importantly, to reconcile the large range of existing estimates of the Twomey effect from different observational methods.

Figure 7 shows the two-dimensional probability density functions of ln Nd and ln SO4 along with fitting lines. We note that the pre-binned method yields similar high correlation coefficients (R) for SO4B (0.96), SO4S (0.95), and SO4C (0.98) due to the data stratification. When moving to the regression on all data points (Table S1), we can see that the R for SO4B is the highest (0.6), followed by SO4S (0.57), and the R for SO4C is the lowest (0.54), consistent with the results reported by Stier (2016) and Painemal et al. (2020). In contrast, the regression slopes for SO4C (0.88) are nearly twice as large as that for SO4B (0.47) and SO4S (0.46) (Fig. 7), implying that the strength of S derived on the basis of column-integrated aerosol quantity, which is often the case for most previous satellite-based estimates, is overestimated by nearly a factor of 2. Note that to explain the same change in ln Nd, ln SO4B and ln SO4S increase by about 5, while ln SO4C only increases by 2 (Fig. 7). Translating to the linear scale, this means that SO4B (SO4S) increases by 148-fold, whereas only a 10-fold increase can be seen in SO4C, resulting in the much larger slope of ln Nd versus ln SO4C. The underlying reason would be that the variability of SO4C is insufficient to explain the variabilities of SO4B (SO4S) .

Figure 8Map of coefficients of variation (CV) of (a) SO4B, (b) SO4S, and (c) SO4C, (d) the ratio of column mass of SO4 below clouds (SO4BC) to SO4C (%) and Pearson's correlation coefficients of SO4B with (e) SO4C and (f) SO4S, which are calculated for each 1× 1 grid box over the period of 2006–2009.

In order to verify whether SO4C has the capability to capture the variability of SO4B quantitatively, the coefficient of variation (CV; calculated as the ratio of the standard deviation to the mean) is employed, which is a measure of relative variability that is particularly useful for the comparison among quantities with different magnitudes and units, e.g., SO4C (in units of µg m−2) versus SO4B or SO4S (in units of µg m−3) here. Since S is generally inferred from the spatiotemporal variability of aerosol and cloud properties, here we calculate the temporal and spatial CVs, respectively; the temporal CV is calculated from the daily time series for the period 2006–2009 for each 1× 1 grid box, and the spatial CV is derived from the multi-annual averaged global geographical distribution. As shown in Fig. 8a, b, and c, the temporal CVs of SO4C are smaller than those of SO4B and SO4S almost everywhere, with globally averaged CVs of 0.52 versus 1.02 and 1.03. Spatially, the larger CVs are generally located over aerosol outflow regions, such as the western North Pacific, the Atlantic, and the east coasts of South America and southern Africa, indicative of an impact of the strong variation of continental, specifically anthropogenic, emissions. Similarly, the spatial CV of SO4C exhibits a much smaller (0.88) value than those of SO4B and SO4S (1.84 and 1.79). In other words, the variability of SO4C is only able to reflect about half of the variability of SO4 near cloud base. This is mainly due to the important role of SO4 above cloud in total column SO4. However, above-cloud aerosol is much more homogeneous compared to SO4B and SO4S that are directly driven by rapid changes in anthropogenic emissions near the surface.

This is demonstrated in Fig. 8d, which shows that the ratio of SO4C below cloud (SO4BC) to SO4C is quite low, with a global average of 11.89 %. Spatially, the ratio can be up to 35 % over aerosol outflow regions but generally below 10 % over vast remote oceans. The low ratio confirms the comparatively small role that sub-cloud aerosols have in determining the aerosol loading within a column. Interestingly, there is also good consistency between the spatial patterns of the ratio of SO4BC to SO4C and the correlation coefficient of SO4C with SO4B (Fig. 8d, e); i.e., the high-ratio regions (the ratio >15 %) generally have strong correlations (R>0.7). Therefore, with regard to the vertical co-location, it is comparatively sensible to use column-integrated quantities such as AOD and AI to represent CCN near cloud base over polluted continents and the immediate outflow region, where the correlation coefficient of SO4C with SO4B is overall larger than 0.7, but this is obviously not the case over remote oceans. The loose correlation between cloud-base and column-integrated aerosols found here (R<0.4), in combination with the detectability limitations of satellite instruments for aerosol loading (P.-L. Ma et al.2018), makes it more challenging to detect any meaningful aerosol–cloud associations in pristine environments from retrieved AOD and AI. Nevertheless, unlike the SO4C, rather strong correlations between SO4S and SO4B (R>0.7) can generally be found with the only exception of high-latitude oceans (Fig. 8f), which, in combination with the highly similar aerosol–Nd slopes and CVs between SO4S and SO4B, hints at surface observations as promising in terms of the vertical co-location issue.

4 Future improvements

Although this study has demonstrated the significant impacts of major confounders on the estimation of Nd-to-CCN sensitivity, some caveats remain. In order to achieve an optimal estimate of radiative forcing from the remote sensing perspective, the following sources of uncertainty should be accounted for in future investigations.

The derivation of Nd from satellite observations relies on a number of assumptions (Grosvenor et al.2018), making it prone to systematic biases. While some sampling strategies have been applied to side-step the biases in Nd retrieval (see Sect. 2), the uncertainties remain. To further ensure the cloud adiabaticity, there are two practical methods for use, including comparing the CER at different wavelengths (Bennartz and Rausch2017) and locating the cloud “core” (Zhu et al.2018). Appropriate Nd sampling strategies are beneficial in future investigations, though they have relatively little impact on S (and the implied RFaci) (Gryspeerdt et al.2021).

The retrieved AOD (AI) as well as reanalyzed SO4 were treated as CCN proxies in this study. However, the usability is limited due to the lack of information on the aerosol size and/or hygroscopicity for AOD (AI) and also due to the fact that SO4 cannot fully explain the variability of CCN since organic aerosols also contribute significantly (Ruehl et al.2016), particularly in the remote marine boundary layer (Zheng et al.2020). Therefore, the application of direct CCN retrievals from polarimetric satellites (Hasekamp et al.2019) is promising in future investigations of aerosol–cloud interactions. However, it would need to be combined with an estimate of the contribution of above-cloud aerosol, especially in regions unaffected by continental outflow. More importantly, the PD CCN–Nd relationship has been shown to be a better approximation of the PI and hence the “actual” sensitivity of Nd to aerosol perturbations than the AOD (AI)–Nd relationship, as it is not affected by the differing PI and PD aerosol environments (Gryspeerdt et al.2017). This again highlights the importance of directly retrieved CCN in the assessment of the radiative forcing from the Twomey effect.

Notably, using a linear regression slope from an ordinary least-squares (OLS) line-fitting method to describe the actual nonlinear aerosol–Nd (Fig. 1) can introduce additional uncertainties related to the problem of regression dilution (Pitkänen et al.2016; Quaas et al.2020). The OLS method is also likely to overestimate the change in Nd from PI to PD over polluted continents, as a saturation effect will occur as aerosols keep rising under a polluted background. A joint histogram method proposed by Gryspeerdt et al. (2017) can be useful to account for the nonlinearity.

In addition to the precipitation, entrainment mixing is a crucial droplet sink process (Blyth et al.1988). However, given that it is practically difficult to infer a quantitative measure of the strength of entrainment mixing from satellite observations, its impacts were not explicitly considered here. It has been proven that entrainment mixing process is associated with dynamical and cloud regimes (Warner1969; de Roode and Wang2007), so the updraft constraint in this study would also incorporate the effect of entrainment mixing to some extent. Although there have been some attempts to characterize entrainment mixing via the combination of lower tropospheric stability and relative humidity near cloud top (Chen et al.2014; Jia et al.2019a) or the Nd–LWP relationship at a certain phase relaxation timescale describing evaporation–entrainment feedback (Zhang et al.2022), they are relatively rough approximations or qualitative differentiations. An updated approach for deriving measures of entrainment mixing at the global scale would be highly beneficial.

It was found that S can vary not only with the spatial resolution of data (Sekiguchi et al.2003; McComiskey and Feingold2012) but also with the spatial scale at which the regression is performed (Grandey and Stier2010). Grandey and Stier (2010) demonstrated that conducting analysis over large regions could induce spurious aerosol–cloud correlations, mainly owing to the spatial covariations in aerosol type, cloud regime, and meteorological conditions. Despite the global analyses employed in this study, the applied updraft constraint may make our results less susceptible to this issue. It is expected that, with joint use of an updraft constraint and CCN retrieval that greatly eliminate the spatial gradient effects, global analysis would be preferable compared to a regional or local method, since the latter could lead to a large bias in the aerosol–Nd slope over pristine oceans where either the instrument detectability limitations on aerosol (P.-L. Ma et al.2018) or the inability of column-integrated measure to represent aerosol near cloud base for low aerosol conditions (see Sect. 3.4) could play a major role.

Given the impossibility of combining all datasets used in different sections together (e.g., the CBH and CGT from Terra are observed at 10:30 but the precipitation from Aqua at 13:30 local solar time), this work evaluates the individual impact of each bias on the estimate of S separately. Nevertheless, the sources of bias could also be correlated with each other; thus, an optimal estimate of S with all biases constrained is desirable. Future studies are being planned to make use of CALIOP/CloudSat satellite observations, which provide simultaneous retrievals of aerosol extinction profiles, precipitation, and cloud-base height (Mülmenstädt et al.2018) such that an analysis accounting for all potential sources of bias can be performed.

Table 3Issues highlighted in this study and their impacts on the overall estimation of S.

Download Print Version | Download XLSX

5 Conclusions and discussion

By employing a statistically robust dataset from multiple active and passive satellite sensors as well as a reanalysis product, we systematically assessed the aerosol impact on marine warm clouds and found that the Nd-to-CCN sensitivity (S) shows a strong dependence on (a) updraft proxy, (b) precipitation, (c) satellite retrieval biases, and (d) vertical co-location between aerosol and cloud layer. The key results and the corresponding implications are summarized as follows, and the impacts of issues highlighted here on the overall estimation of S are listed in Table 3.

  1. SAOD and SAI are found to increase remarkably with both CBH and CGT (treated as proxies for vertical velocity at cloud base), suggesting that stronger aerosol–cloud interactions generally occur under larger updraft velocity conditions. Although a similar dependency has been reported by some previous studies utilizing in situ aircraft measurements or ground-based remote sensing, they were limited to certain time periods and regions. Instead, S here is characterized as a function of CBH (CGT) based on 4 years of global satellite observations, which can thus reflect the full variability of cloud dynamic conditions. This functional relationship, as a better alternative to large-scale meteorological condition constraints (less directly linked to cloud dynamics on a cloud scale), could be promising in application to the estimation of global aerosol–cloud radiative forcing, by which the change in Nd from the PI to the PD may be inferred based on CBH (CGT) climatology from satellite and anthropogenic aerosol emission perturbation assuming first-order unchanged CBH distributions.

  2. There is an intensive sink of cloud droplets by precipitation, thereby leading to a much lower Nd in raining clouds (55 cm−3) compared to non-raining clouds (125 cm−3). In turn, a weaker S was found in raining clouds than in non-raining clouds, with the corresponding SAI of 0.45 versus 0.56, respectively. Surprisingly, after lumping all cloud scenes together, the derived SAI (0.68) is amplified by 21 % (51 %) relative to the non-raining (raining) case, and also a more nonlinear aerosol–Nd relationship is diagnosed. We showed that this amplification is just an artifact governed by the joint impacts of the suppression of precipitation by aerosols and the aerosol removal by precipitation. That is, introducing the confounding effect of aerosol–precipitation interactions into the estimation of the Twomey effect can artificially bias the S to a higher value. The finding highlights the necessity of removing precipitating clouds from statistical analyses when quantifying S and assessing the Twomey effect. To achieve this, the only way would be simultaneous aerosol–cloud–precipitation retrievals (e.g., from the A-Train satellite constellation). However, due to the fact that most of existing estimates of S and its radiative forcing did not take this aspect into consideration, the relative change in SAI from all clouds to non-raining clouds presented here could serve as a useful reference for the intercomparison of cloud susceptibilities from different studies.

  3. Aerosol retrieval biases (3D radiative effects and cloud contamination), aerosol swelling, and cloud retrieval bias (heterogeneity effect) tend to lead to an underestimation of S. Although SAI (SAOD) for the first ΔL bin, in which evident AI(AOD) enhancement exists, is about 29 % (50 %) less than other unaffected bins, the overall underestimation is only ∼3 % because of the small data volume in the first bin (Fig. 5a). Nevertheless, for low-ΔL-dominated regions (e.g., stratocumulus regions), the underestimation can be more pronounced. By comparing SAI (SAOD) calculated by NdAll and Nd, the underestimation by cloud retrieval issues is roughly estimated to be ∼8 % (∼17 %). It is noted that the CF can act as a key modulator of these two kinds of retrieval issues; i.e., an increase in CF enhances the aerosol retrieval biases via intensifying near-cloud enhancement of AOD (AI) but reduces cloud retrieval errors via alleviating the cloud heterogeneity, making it practically difficult to balance the accuracies of both retrievals within the same grid. In terms of aerosol, the use of aerosol reanalysis is a potential way to avoid the near-cloud enhancement of AOD (AI), but note that the issue of aerosol swelling remains to some extent. As for Nd, the retrievals under high CF (over a 1× 1 grid) conditions would be preferable even though strict criteria for cloud screening (Grosvenor et al.2018) have been applied, which, however, could incur an artificial selection of cloud regime since CF also covaries with cloud dynamics. Therefore, applying a CF updraft constraint in the Nd screening would be a path forward.

  4. Use of vertically integrated SO4 (SO4C) as a proxy for CCN near cloud base results in a degradation of correlation with Nd, with an approximately twofold enhancement of S compared to using SO4 near cloud base (SO4B). This is mostly attributed to the inability of SO4C to capture the full variability of SO4B. Generally, SO4C is dominated by SO4 above cloud, which is relatively homogeneous compared to SO4B that is tightly linked to rapid changes in anthropogenic but also natural emissions near the surface. As a result, to explain the same change in Nd, the corresponding fractional change in SO4C is much smaller than SO4B, hence leading to a higher regression slope that, however, is not associated with physically meaningful enhancement of S. The similar aerosol–Nd slopes, correlation coefficients, and relative variability between SO4S (SO4 near the surface) and SO4B suggest that the use of near-surface aerosol measurements, such as particulate matter (Guo et al.2018) or aerosol extinction coefficients (Liu and Li2018), is an effective solution to the problem of vertical co-location in the case that observations of the vertical profile of aerosol and cloud-base height are unavailable, although its suitability would depend on the degree of coupling of the boundary layer (Painemal et al.2020). Moreover, the result further raises complications to compare and reconcile the diverse cloud susceptibilities from studies utilizing CCN proxies at different altitudes. It should be noted that the derivation of Nd change from PI to PD (thus radiative forcing) is expected to be less affected given that the vertical co-location issue also applies to fractional change in aerosol due to anthropogenic emissions, thus partly compensating for the enhancement of S; nevertheless, the net effect on radiative forcing still needs further exploration.

Data availability

The MODIS Aqua and Terra Level 3 products are available from (Platnick et al.2017d) and (Platnick et al.2017c), and Level 2 products are available from (Platnick et al.2017b) and (Platnick et al.2017a). The CloudSat data are available from (last access: 6 July 2021; CDPC2021). The MISR Level 2 Cloud Product (MIL2TCSP) data are from (last access: 6 July 2021; NASA/LARC/SD/ASDC2012). The MERRA-2 and ERA5 reanalysis products are collected from (last access: 16 November 2021; GMAO2021) and (last access: 2 March 2022; ECMWF2022).


The supplement related to this article is available online at:

Author contributions

HJ and JQ designed the research. HJ performed the research and prepared the paper, with comments from JQ, EG, CB, and OS.

Competing interests

At least one of the co-authors is an associate editor of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


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


The authors gratefully acknowledge funding by the German Research Foundation (joint call between the National Science Foundation of China and the Deutsche Forschungsgemeinschaft, DFG, GZ QU 311/28-1, project “CloudTrend”) and by the European Union Horizon 2020 project FORCES (grant no. 821205). Edward Gryspeerdt would like to thank the Royal Society University Research Fellowship (URF/R1/191602). The publication of this article was funded by the Open Access Publishing Fund of Leipzig University supported by the German Research Foundation within the program Open Access Publication Funding.

MODIS data were acquired from the Level-1 and Atmosphere Archive & Distribution System (LAADS) Distributed Active Archive Center (DAAC). MISR products were obtained from the NASA Langley Research Center Atmospheric Science Data Center. CloudSat data products were provided by the CloudSat Data Processing Center at the Cooperative Institute for Research in the Atmosphere, Colorado State University. MERRA-2 reanalysis products were provided by NASA's Global Monitoring and Assimilation Office (GMAO). ERA5 reanalysis datasets were retrieved from ECMWF's Meteorological Archival and Retrieval System (MARS). This work also used JASMIN, the UK's collaborative data analysis environment (, last access: 1 March 2021).

Financial support

This research has been supported by the German Research Foundation (grant no.GZ QU 311/28-1) and by the European Union Horizon 2020 project FORCES (grant no. 821205). Edward Gryspeerdt was supported by a Royal Society University Research Fellowship (grant no. URF/R1/191602).

Review statement

This paper was edited by Timothy Garrett and reviewed by two anonymous referees.


Ackerman, A. S., Kirkpatrick, M. P., Stevens, D. E., and Toon, O. B.: The impact of humidity above stratiform clouds on indirect aerosol climate forcing, Nature, 432, 1014–1017,, 2004. a

Albrecht, B. A.: Aerosols, Cloud Microphysics, and Fractional Cloudiness, Science, 80, 245, 1227–1230,, 1989. a

Allen, R. J., Amiri-Farahani, A., Lamarque, J.-F., Smith, C., Shindell, D., Hassan, T., and Chung, C. E.: Observationally constrained aerosol–cloud semi-direct effects, Clim. Atmos. Sci., 2, 16,, 2019. a

Anderson, T. L., Charlson, R. J., Winker, D. M., Ogren, J. A., and Holmén, K.: Mesoscale Variations of Tropospheric Aerosols, J. Atmos. Sci., 60, 119–136,<0119:MVOTA>2.0.CO;2, 2003. a, b, c

Bellouin, N., Quaas, J., Gryspeerdt, E., Kinne, S., Stier, P., Watson‐Parris, D., Boucher, O., Carslaw, K. S., Christensen, M., Daniau, A., Dufresne, J., Feingold, G., Fiedler, S., Forster, P., Gettelman, A., Haywood, J. M., Lohmann, U., Malavelle, F., Mauritsen, T., McCoy, D. T., Myhre, G., Mülmenstädt, J., Neubauer, D., Possner, A., Rugenstein, M., Sato, Y., Schulz, M., Schwartz, S. E., Sourdeval, O., Storelvmo, T., Toll, V., Winker, D., and Stevens, B.: Bounding Global Aerosol Radiative Forcing of Climate Change, Rev. Geophys., 58, e2019RG000660,, 2020. a

Bennartz, R. and Rausch, J.: Global and regional estimates of warm cloud droplet number concentration based on 13 years of AQUA-MODIS observations, Atmos. Chem. Phys., 17, 9815–9836,, 2017. a

Berg, L. K., Berkowitz, C. M., Barnard, J. C., Senum, G., and Springston, S. R.: Observations of the first aerosol indirect effect in shallow cumuli, Geophys. Res. Lett., 38, L03809,, 2011. a

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

Böhm, C., Sourdeval, O., Mülmenstädt, J., Quaas, J., and Crewell, S.: Cloud base height retrieval from multi-angle satellite data, Atmos. Meas. Tech., 12, 1841–1860,, 2019. a, b, c, d

Boucher, O. and Quaas, J.: Water vapour affects both rain and aerosol optical depth, Nat. Geosci., 6, 4–5,, 2013. a

Bougiatioti, A., Nenes, A., Lin, J. J., Brock, C. A., de Gouw, J. A., Liao, J., Middlebrook, A. M., and Welti, A.: Drivers of cloud droplet number variability in the summertime in the southeastern United States, Atmos. Chem. Phys., 20, 12163–12176,, 2020. a

Buchard, V., Randles, C. A., da Silva, A. M., Darmenov, A., Colarco, P. R., Govindaraju, R., Ferrare, R., Hair, J., Beyersdorf, A. J., Ziemba, L. D., and Yu, H.: The MERRA-2 Aerosol Reanalysis, 1980 Onward. Part II: Evaluation and Case Studies, J. Climate, 30, 6851–6872,, 2017. a

CDPC: 2B-CLDCLASS product, CloudSat Data Processing Center at the Cooperative Institute for Research in the Atmosphere, Colorado State University [data set],, last access: 6 July 2021. a

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. a

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. a, b, c, d

Christensen, M. W., Neubauer, D., Poulsen, C. A., Thomas, G. E., McGarragh, G. R., Povey, A. C., Proud, S. R., and Grainger, R. G.: Unveiling aerosol–cloud interactions – Part 1: Cloud contamination in satellite products enhances the aerosol indirect forcing estimate, Atmos. Chem. Phys., 17, 13151–13164,, 2017. a

Coakley, J. A., Friedman, M. A., and Tahnk, W. R.: Retrieval of Cloud Properties for Partly Cloudy Imager Pixels, J. Atmos. Ocean. Technol., 22, 3–17,, 2005. a

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

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. a

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

ECMWF: ECMWF Reanalysis v5 (ERA5), European Centre for Medium-Range Weather Forecasts [data set],, last access: 2 March 2022. a

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

Feingold, G., Remer, L. A., Ramaprasad, J., and Kaufman, Y. J.: Analysis of smoke impact on clouds in Brazilian biomass burning regions: An extension of Twomey's approach, J. Geophys. Res.-Atmos., 106, 22907–22922,, 2001. a

Feingold, G., Eberhard, W. L., Veron, D. E., and Previdi, M.: First measurements of the Twomey indirect effect using ground-based remote sensors, Geophys. Res. Lett., 30, 1287,, 2003. a, b

Feingold, G., Goren, T., and Yamaguchi, T.: Quantifying albedo susceptibility biases in shallow clouds, Atmos. Chem. Phys., 22, 3303–3319,, 2022. a

Forster, P. M., Storelvmo, T., Armour, K., Collins, W., Dufresne, J. L., Frame, D., Lunt, D. J., Mauritsen, T., Palmer, M. D., Watanabe, M., Wild, M., and Zhang, H.: Chapter 7: The Earth's Energy Budget, Climate Feedbacks, and Climate Sensitivity, Clim. Change 2021 Phys. Sci. Basis. Contrib. Work. Gr. I to Sixth Assess. Rep. Intergov. Panel Clim. Change, 2021. a

Freud, E. and Rosenfeld, D.: Linear relation between convective cloud drop number concentration and depth for rain initiation, J. Geophys. Res.-Atmos., 117, D02207,, 2012. a

Gerber, H.: Microphysics of marine stratocumulus clouds with two drizzle modes, J. Atmos. Sci., 53, 1649–1662,<1649:MOMSCW>2.0.CO;2, 1996. a

GMAO: Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), Global Modeling and Assimilation Office [data set],, last access: 16 November 2021. a

Grandey, B. S. and Stier, P.: A critical look at spatial scale choices in satellite-based aerosol indirect effect studies, Atmos. Chem. Phys., 10, 11459–11470,, 2010. a, b

Grosvenor, D. P., Sourdeval, O., Zuidema, P., Ackerman, A., Alexandrov, M. D., Bennartz, R., Boers, R., Cairns, B., Chiu, J. C., Christensen, M., Deneke, H., Diamond, M., Feingold, G., Fridlind, A., Hünerbein, A., Knist, C., Kollias, P., Marshak, A., McCoy, D., Merk, D., Painemal, D., Rausch, J., Rosenfeld, D., Russchenberg, H., Seifert, P., Sinclair, K., Stier, P., van Diedenhoven, B., Wendisch, M., Werner, F., Wood, R., Zhang, Z., and Quaas, J.: Remote Sensing of Droplet Number Concentration in Warm Clouds: A Review of the Current State of Knowledge and Perspectives, Rev. Geophys., 56, 409–453,, 2018. a, b, c, d, e, f

Gryspeerdt, E. and Stier, P.: Regime-based analysis of aerosol-cloud interactions, Geophys. Res. Lett., 39, L21802,, 2012. a, b

Gryspeerdt, E., Stier, P., White, B. A., and Kipling, Z.: Wet scavenging limits the detection of aerosol effects on precipitation, Atmos. Chem. Phys., 15, 7557–7570,, 2015. a, b

Gryspeerdt, E., Quaas, J., Ferrachat, S., Gettelman, A., Ghan, S., Lohmann, U., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., Wang, M., and Zhang, K.: Constraining the instantaneous aerosol influence on cloud albedo, P. Natl. Acad. Sci. USA, 114, 4899–4904,, 2017. a, b, c, d, e, f, g

Gryspeerdt, E., Mülmenstädt, J., Gettelman, A., Malavelle, F. F., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., Wang, M., and Zhang, K.: Surprising similarities in model and observational aerosol radiative forcing estimates, Atmos. Chem. Phys., 20, 613–623,, 2020. a

Gryspeerdt, E., McCoy, D. T., Crosbie, E., Moore, R. H., Nott, G. J., Painemal, D., Small-Griswold, J., Sorooshian, A., and Ziemba, L.: The impact of sampling strategy on the cloud droplet number concentration estimated from satellite data, Atmos. Meas. Tech. Discuss. [preprint],, in review, 2021. a, b, c, d

Guo, J., Liu, H., Li, Z., Rosenfeld, D., Jiang, M., Xu, W., Jiang, J. H., He, J., Chen, D., Min, M., and Zhai, P.: Aerosol-induced changes in the vertical structure of precipitation: a perspective of TRMM precipitation radar, Atmos. Chem. Phys., 18, 13329–13343,, 2018. a

Hasekamp, O. P., Gryspeerdt, E., and Quaas, J.: Analysis of polarimetric satellite measurements suggests stronger cooling due to aerosol-cloud interactions, Nat. Commun., 10, 5405,, 2019. a, b, c, d, e, f

Haywood, J. M., Osborne, S. R., and Abel, S. J.: The effect of overlying absorbing aerosol layers on remote sensing retrievals of cloud effective radius and cloud optical depth, Q. J. Roy. Meteorol. Soc., 130, 779–800,, 2004. a

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz‐Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.: The ERA5 global reanalysis, Q. J. Roy. Meteorol. Soc., 146, 1999–2049,, 2020. a

Hudson, J. G. and Noble, S.: CCN and Vertical Velocity Influences on Droplet Concentrations and Supersaturations in Clean and Polluted Stratus Clouds, J. Atmos. Sci., 71, 312–331,, 2014. a

Jia, H., Ma, X., Quaas, J., Yin, Y., and Qiu, T.: Is positive correlation between cloud droplet effective radius and aerosol optical depth over land due to retrieval artifacts or real physical processes?, Atmos. Chem. Phys., 19, 8879–8896,, 2019a. a, b, c, d

Jia, H., Ma, X., Yu, F., Liu, Y., and Yin, Y.: Distinct Impacts of Increased Aerosols on Cloud Droplet Number Concentration of Stratus/Stratocumulus and Cumulus, Geophys. Res. Lett., 46, 13517–13525,, 2019b. a

Jia, H., Ma, X., Yu, F., and Quaas, J.: Significant underestimation of radiative forcing by aerosol–cloud interactions derived from satellite-based methods, Nat. Commun., 12, 3649,, 2021. a, b, c, d, e, f

Kacarab, M., Thornhill, K. L., Dobracki, A., Howell, S. G., O'Brien, J. R., Freitag, S., Poellot, M. R., Wood, R., Zuidema, P., Redemann, J., and Nenes, A.: Biomass burning aerosol as a modulator of the droplet number in the southeast Atlantic region, Atmos. Chem. Phys., 20, 3029–3040,, 2020. a

Kaufman, Y. J., Koren, I., Remer, L. A., Rosenfeld, D., and Rudich, Y.: The effect of smoke, dust, and pollution aerosol on shallow cloud development over the Atlantic Ocean, P. Natl. Acad. Sci. USA, 102, 11207–11212,, 2005. a

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

Lareau, N. P., Zhang, Y., and Klein, S. A.: Observed Boundary Layer Controls on Shallow Cumulus at the ARM Southern Great Plains Site, J. Atmos. Sci., 75, 2235–2255,, 2018. a

Levy, R. C., Mattoo, S., Munchak, L. A., Remer, L. A., Sayer, A. M., Patadia, F., and Hsu, N. C.: The Collection 6 MODIS aerosol products over land and ocean, Atmos. Meas. Tech., 6, 2989–3034,, 2013. a, b

Li, Z., Zhao, F., Liu, J., Jiang, M., Zhao, C., and Cribb, M.: Opposite effects of absorbing aerosols on the retrievals of cloud optical depth from spaceborne and ground-based measurements, J. Geophys. Res.-Atmos., 119, 5104–5114,, 2014. a

Liu, J. and Li, Z.: Significant Underestimation in the Optically Based Estimation of the Aerosol First Indirect Effect Induced by the Aerosol Swelling Effect, Geophys. Res. Lett., 45, 5690–5699,, 2018. a, b

Ma, P.-L., Rasch, P. J., Chepfer, H., Winker, D. M., and Ghan, S. J.: Observational constraint on cloud susceptibility weakened by aerosol retrieval limitations, Nat. Commun., 9, 2640,, 2018. a, b, c, d, e

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

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. a

Matheson, M. A., Coakley, J. A., and Tahnk, W. R.: Multiyear Advanced Very High Resolution Radiometer observations of summertime stratocumulus collocated with aerosols in the northeastern Atlantic, J. Geophys. Res., 111, D15206,, 2006. a

McComiskey, A. and Feingold, G.: The scale problem in quantifying aerosol indirect effects, Atmos. Chem. Phys., 12, 1031–1049,, 2012. a, b

McComiskey, A., Feingold, G., Frisch, A. S., Turner, D. D., Miller, M. A., Chiu, J. C., Min, Q., and Ogren, J. A.: An assess- ment of aerosol-cloud interactions in marine stratus clouds based on surface remote sensing, J. Geophys. Res., 114, D09203,, 2009. a, b

McCoy, D. T., Bender, F. A., 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. a, b, c

Mülmenstädt, J., Sourdeval, O., Henderson, D. S., L'Ecuyer, T. S., Unglaub, C., Jungandreas, L., Böhm, C., Russell, L. M., and Quaas, J.: Using CALIOP to estimate cloud-field base height and its uncertainty: the Cloud Base Altitude Spatial Extrapolator (CBASE) algorithm and dataset, Earth Syst. Sci. Data, 10, 2279–2293,, 2018. a, b

Nakajima, T., Higurashi, A., Kawamoto, K., and Penner, J. E.: A possible correlation between satellite-derived cloud and aerosol microphysical parameters, Geophys. Res. Lett., 28, 1171–1174,, 2001. a

NASA/LARC/SD/ASDC: MISR Level 2 TOA/Cloud Height and Motion parameters V001, NASA Langley Atmospheric Science Data Center DAAC [data set], (last access: 6 July 2021), 2012. a, b

Neubauer, D., Christensen, M. W., Poulsen, C. A., and Lohmann, U.: Unveiling aerosol–cloud interactions – Part 2: Minimising the effects of aerosol swelling and wet scavenging in ECHAM6-HAM2 for comparison to satellite data, Atmos. Chem. Phys., 17, 13165–13185,, 2017. a

Painemal, D., Chang, F.-L., Ferrare, R., Burton, S., Li, Z., Smith Jr., W. L., Minnis, P., Feng, Y., and Clayton, M.: Reducing uncertainties in satellite estimates of aerosol–cloud interactions over the subtropical ocean by integrating vertically resolved aerosol observations, Atmos. Chem. Phys., 20, 7167–7177,, 2020. a, b, c, d, e, f

Penner, J. E., Xu, L., and Wang, M.: Satellite methods underestimate indirect climate forcing by aerosols, P. Natl. Acad. Sci. USA, 108, 13404–13408,, 2011. a

Pitkänen, M. R. A., Mikkonen, S., Lehtinen, K. E. J., Lipponen, A., and Arola, A.: Artificial bias typically neglected in comparisons of uncertain atmospheric data, Geophys. Res. Lett., 43, 10003–10011,, 2016. a

Platnick, S., Ackerman, S., King, M. D., Wind, G., Meyer, K., Menzel, P., Frey, R., Holz, R. E., Baum, B., and Yang, P.: MODIS atmosphere L2 cloud product (06_L2), (Terra), LAADS [data set], 2017a. a, b

Platnick, S., Ackerman, S., King, M. D., Wind, G., Meyer, K., Menzel, P., Frey, R., Holz, R. E., Baum, B., and Yang, P.: MODIS atmosphere L2 cloud product (06_L2), (Aqua), LAADS [data set], 2017b. a, b

Platnick, S., King, M. D., and Hubanks, P. A.: MODIS Atmosphere L3 Daily Product, (Terra), LAADS [data set], 2017c. a, b

Platnick, S., King, M. D., and Hubanks, P. A.: MODIS Atmosphere L3 Daily Product, (Aqua), LAADS [data set], 2017d. a, b

Platnick, S., Meyer, K. G., King, M. D., Wind, G., Amarasinghe, N., Marchant, B., Arnold, G. T., Zhang, Z., Hubanks, P. A., Holz, R. E., Yang, P., Ridgway, W. L., and Riedi, J.: The MODIS Cloud Optical and Microphysical Products: Collection 6 Updates and Examples From Terra and Aqua, IEEE Trans. Geosci. Remote Sens., 55, 502–525,, 2017e. a

Qiu, Y., Zhao, C., Guo, J., and Li, J.: 8-Year ground-based observational analysis about the seasonal variation of the aerosol-cloud droplet effective radius relationship at SGP site, Atmos. Environ., 164, 139–146,, 2017. a

Quaas, J., Boucher, O., and Lohmann, U.: Constraining the total aerosol indirect effect in the LMDZ and ECHAM4 GCMs using MODIS satellite data, Atmos. Chem. Phys., 6, 947–955,, 2006. a, b

Quaas, J., Boucher, O., Bellouin, N., and Kinne, S.: Satellite-based estimate of the direct and indirect aerosol climate forcing, J. Geophys. Res.-Atmos., 113, D05204,, 2008. a, b, c, d

Quaas, J., Stevens, B., Stier, P., and Lohmann, U.: Interpreting the cloud cover – aerosol optical depth relationship found in satellite data using a general circulation model, Atmos. Chem. Phys., 10, 6129–6135,, 2010. a, b, c

Quaas, J., Arola, A., Cairns, B., Christensen, M., Deneke, H., Ekman, A. M. L., Feingold, G., Fridlind, A., Gryspeerdt, E., Hasekamp, O., Li, Z., Lipponen, A., Ma, P.-L., Mülmenstädt, J., Nenes, A., Penner, J. E., Rosenfeld, D., Schrödner, R., Sinclair, K., Sourdeval, O., Stier, P., Tesche, M., van Diedenhoven, B., and Wendisch, M.: Constraining the Twomey effect from satellite observations: issues and perspectives, Atmos. Chem. Phys., 20, 15079–15099,, 2020. a, b, c

Randles, C. A., da Silva, A. M., Buchard, V., Colarco, P. R., Darmenov, A., Govindaraju, R., Smirnov, A., Holben, B., Ferrare, R., Hair, J., Shinozuka, Y., and Flynn, C. J.: The MERRA-2 Aerosol Reanalysis, 1980 Onward. Part I: System Description and Data Assimilation Evaluation, J. Climate, 30, 6823–6850,, 2017. a

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. a, b

Rosenfeld, D., Zhu, Y., Wang, M., Zheng, Y., Goren, T., and Yu, S.: Aerosol-driven droplet concentrations dominate coverage and water of oceanic low-level clouds, Science, 80, 363,, 2019. a, b, c

Ruehl, C. R., Davies, J. F., and Wilson, K. R.: An interfacial mechanism for cloud droplet formation on organic aerosols, Science, 80, 351, 1447–1450,, 2016. a

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. a, b

Sassen, K. and Wang, Z.: Classifying clouds around the globe with the CloudSat radar: 1-year of results, Geophys. Res. Lett., 35, L04805,, 2008. a

Schmidt, J., Ansmann, A., Bühl, J., and Wandinger, U.: Strong aerosol–cloud interaction in altocumulus during updraft periods: lidar observations over central Europe, Atmos. Chem. Phys., 15, 10687–10700,, 2015. a, b

Seinfeld, J. H., Bretherton, C., Carslaw, K. S., Coe, H., DeMott, P. J., Dunlea, E. J., Feingold, G., Ghan, S., Guenther, A. B., Kahn, R., Kraucunas, I., Kreidenweis, S. M., Molina, M. J., Nenes, A., Penner, J. E., Prather, K. A., Ramanathan, V., Ramaswamy, V., Rasch, P. J., Ravishankara, A. R., Rosenfeld, D., Stephens, G., and Wood, R.: Improving our fundamental understanding of the role of aerosol-cloud interactions in the climate system, P. Natl. Acad. Sci. USA, 113, 5781–5790,, 2016. a

Sekiguchi, M., Nakajima, T., Suzuki, K., Kawamoto, K., Higurashi, A., Rosenfeld, D., Sano, I., and Mukai, S.: A study of the direct and indirect effects of aerosols using global satellite data sets of aerosol and cloud parameters, J. Geophys. Res.-Atmos., 108, 2002JD003359,, 2003. a, b

Smith, C. J., Kramer, R. J., Myhre, G., Alterskjær, K., Collins, W., Sima, A., Boucher, O., Dufresne, J.-L., Nabat, P., Michou, M., Yukimoto, S., Cole, J., Paynter, D., Shiogama, H., O'Connor, F. M., Robertson, E., Wiltshire, A., Andrews, T., Hannay, C., Miller, R., Nazarenko, L., Kirkevåg, A., Olivié, D., Fiedler, S., Lewinschal, A., Mackallah, C., Dix, M., Pincus, R., and Forster, P. M.: Effective radiative forcing and adjustments in CMIP6 models, Atmos. Chem. Phys., 20, 9591–9618,, 2020. a

Sourdeval, O., C.‐Labonnote, L., Baran, A. J., Mülmenstädt, J., and Brogniez, G.: A methodology for simultaneous retrieval of ice and liquid water cloud properties. Part 2: Near‐global retrievals and evaluation against A‐Train products, Q. J. Roy. Meteorol. Soc., 142, 3063–3081,, 2016. a

Stier, P.: Limitations of passive remote sensing to constrain global cloud condensation nuclei, Atmos. Chem. Phys., 16, 6595–6607,, 2016. a, b, c, d, e

Twohy, C. H., Coakley, J. A., and Tahnk, W. R.: Effect of changes in relative humidity on aerosol scattering near clouds, J. Geophys. Res., 114, 5205,, 2009. a, b, c

Twomey, S.: Pollution and the planetary albedo, Atmos. Environ., 8, 1251–1256,, 1974. a

Várnai, T. and Marshak, A.: MODIS observations of enhanced clear sky reflectance near clouds, Geophys. Res. Lett., 36, L06807,, 2009. a

Várnai, T. and Marshak, A.: Effect of Cloud Fraction on Near-Cloud Aerosol Behavior in the MODIS Atmospheric Correction Ocean Color Product, Remote Sens., 7, 5283–5299,, 2015. a, b

Várnai, T. and Marshak, A.: Analysis of Near-Cloud Changes in Atmospheric Aerosols Using Satellite Observations and Global Model Simulations, Remote Sens., 13, 1151,, 2021. a

Warner, J.: The Microstructure of Cumulus Cloud. Part II. The Effect on Droplet Size Distribution of the Cloud Nucleus Spectrum and Updraft Velocity, J. Atmos. Sci., 26, 1272–1282,<1272:TMOCCP>2.0.CO;2, 1969. a

Yang, Y., Zhao, C., Wang, Y., Zhao, X., Sun, W., Yang, J., Ma, Z., and Fan, H.: Multi‐Source Data Based Investigation of Aerosol‐Cloud Interaction Over the North China Plain and North of the Yangtze Plain, J. Geophys. Res.-Atmos., 126, e2021JD035609,, 2021. a, b

Zhang, J., Reid, J. S., and Holben, B. N.: An analysis of potential cloud artifacts in MODIS over ocean aerosol optical thickness products, Geophys. Res. Lett., 32, L15803,, 2005. a

Zhang, J., Zhou, X., Goren, T., and Feingold, G.: Albedo susceptibility of northeastern Pacific stratocumulus: the role of covarying meteorological conditions, Atmos. Chem. Phys., 22, 861–880,, 2022. a, b

Zhang, Z. and Platnick, S.: An assessment of differences between cloud effective particle radius retrievals for marine water clouds from three MODIS spectral bands, J. Geophys. Res., 116, D20215,, 2011. a, b

Zhao, C. and Garrett, T. J.: Effects of Arctic haze on surface cloud radiative forcing, Geophys. Res. Lett., 42, 557–564,, 2015. a

Zhao, C., Xie, S., Klein, S. A., Protat, A., Shupe, M. D., McFarlane, S. A., Comstock, J. M., Delanoë, J., Deng, M., Dunn, M., Hogan, R. J., Huang, D., Jensen, M. P., Mace, G. G., McCoy, R., O'Connor, E. J., Turner, D. D., and Wang, Z.: Toward understanding of differences in current cloud retrievals of ARM ground-based measurements, J. Geophys. Res.-Atmos., 117, D10206,, 2012. a

Zhao, C., Zhao, L., and Dong, X.: A Case Study of Stratus Cloud Properties Using In Situ Aircraft Observations over Huanghua, China, Atmosphere (Basel), 10, 19,, 2019. a

Zheng, G., Kuang, C., Uin, J., Watson, T., and Wang, J.: Large contribution of organics to condensational growth and formation of cloud condensation nuclei (CCN) in the remote marine boundary layer, Atmos. Chem. Phys., 20, 12515–12525,, 2020. a

Zheng, Y.: Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part I: Maritime Conditions, J. Atmos. Sci., 76, 2539–2558,, 2019. a

Zheng, Y. and Rosenfeld, D.: Linear relation between convective cloud base height and updrafts and application to satellite retrievals, Geophys. Res. Lett., 42, 6485–6491,, 2015.  a, b, c

Zheng, Y., Rosenfeld, D., and Li, Z.: Quantifying cloud base updraft speeds of marine stratocumulus from cloud top radiative cooling, Geophys. Res. Lett., 43, 11407–11413,, 2016. a

Zhu, Y., Rosenfeld, D., and Li, Z.: Under What Conditions Can We Trust Retrieved Cloud Drop Concentrations in Broken Marine Stratocumulus?, J. Geophys. Res.-Atmos., 123, 8754–8767,, 2018. a, b

Short summary
Aerosol–cloud interaction is the most uncertain component of the anthropogenic forcing of the climate. By combining satellite and reanalysis data, we show that the strength of the Twomey effect (S) increases remarkably with vertical velocity. Both the confounding effect of aerosol–precipitation interaction and the lack of vertical co-location between aerosol and cloud are found to overestimate S, whereas the retrieval biases in aerosol and cloud appear to underestimate S.
Final-revised paper