Investigation of aerosol–cloud interactions under different absorptive aerosol regimes using Atmospheric Radiation Measurement (ARM) southern Great Plains (SGP) ground-based measurements

The aerosol indirect effect on cloud microphysical and radiative properties is one of the largest uncertainties in climate simulations. In order to investigate the aerosol–cloud interactions, a total of 16 low-level stratus cloud cases under daytime coupled boundary-layer conditions are selected over the southern Great Plains (SGP) region of the United States. The physicochemical properties of aerosols and their impacts on cloud microphysical properties are examined using data collected from the Department of Energy Atmospheric Radiation Measurement (ARM) facility at the SGP site. The aerosol–cloud interaction index (ACIr) is used to quantify the aerosol impacts with respect to cloud-droplet effective radius. The mean value of ACIr calculated from all selected samples is 0.145± 0.05 and ranges from 0.09 to 0.24 at a range of cloud liquid water paths (LWPs; LWP= 20–300 g m−2). The magnitude of ACIr decreases with an increasing LWP, which suggests a diminished cloud microphysical response to aerosol loading, presumably due to enhanced condensational growth processes and enlarged particle sizes. The impact of aerosols with different light-absorbing abilities on the sensitivity of cloud microphysical responses is also investigated. In the presence of weak light-absorbing aerosols, the low-level clouds feature a higher number concentration of cloud condensation nuclei (NCCN) and smaller effective radii (re), while the opposite is true for strong light-absorbing aerosols. Furthermore, the mean activation ratio of aerosols to CCN (NCCN/Na) for weakly (strongly) absorbing aerosols is 0.54 (0.45), owing to the aerosol microphysical effects, particularly the different aerosol compositions inferred by their absorptive properties. In terms of the sensitivity of cloud-droplet number concentration (Nd) to NCCN, the fraction of CCN that converted to cloud droplets (Nd/NCCN) for the weakly (strongly) absorptive regime is 0.69 (0.54). The measured ACIr values in the weakly absorptive regime are relatively higher, indicating that clouds have greater microphysical responses to aerosols, owing to the favorable thermodynamic condition. The reduced ACIr values in the strongly absorptive regime are due to the cloud-layer heating effect induced by strong lightabsorbing aerosols. Consequently, we expect larger shortwave radiative cooling effects from clouds in the weakly absorptive regime than those in the strongly absorptive regime.

Abstract. The aerosol indirect effect on cloud microphysical and radiative properties is one of the largest uncertainties in climate simulations. In order to investigate the aerosol-cloud interactions, a total of 16 low-level stratus cloud cases under daytime coupled boundary-layer conditions are selected over the southern Great Plains (SGP) region of the United States. The physicochemical properties of aerosols and their impacts on cloud microphysical properties are examined using data collected from the Department of Energy Atmospheric Radiation Measurement (ARM) facility at the SGP site. The aerosol-cloud interaction index (ACI r ) is used to quantify the aerosol impacts with respect to cloud-droplet effective radius. The mean value of ACI r calculated from all selected samples is 0.145 ± 0.05 and ranges from 0.09 to 0.24 at a range of cloud liquid water paths (LWPs; LWP = 20-300 g m −2 ). The magnitude of ACI r decreases with an increasing LWP, which suggests a diminished cloud microphysical response to aerosol loading, presumably due to enhanced condensational growth processes and enlarged particle sizes. The impact of aerosols with different light-absorbing abilities on the sensitivity of cloud microphysical responses is also investigated. In the presence of weak light-absorbing aerosols, the low-level clouds feature a higher number concentration of cloud condensation nuclei (N CCN ) and smaller effective radii (r e ), while the opposite is true for strong light-absorbing aerosols. Furthermore, the mean activation ratio of aerosols to CCN (N CCN /N a ) for weakly (strongly) absorbing aerosols is 0.54 (0.45), owing to the aerosol microphysical effects, particularly the different aerosol compositions inferred by their absorptive properties. In terms of the sensitivity of cloud-droplet number concentration (N d ) to N CCN , the fraction of CCN that converted to cloud droplets (N d /N CCN ) for the weakly (strongly) absorptive regime is 0.69 (0.54). The measured ACI r values in the weakly absorptive regime are relatively higher, indicating that clouds have greater microphysical responses to aerosols, owing to the favorable thermodynamic condition. The reduced ACI r values in the strongly absorptive regime are due to the cloud-layer heating effect induced by strong lightabsorbing aerosols. Consequently, we expect larger shortwave radiative cooling effects from clouds in the weakly absorptive regime than those in the strongly absorptive regime. ment of cloud-droplet number concentration (N d ) and the reduction of cloud-droplet effective radii (r e ), which results in an increase in cloud albedo. This phenomenon is defined as the aerosol first indirect effect (Twomey, 1977), and it is denoted as a general cooling effect in terms of global radiation balance. More fundamentally, the aerosol effects on cloud reflectance result from the cloud microphysical response to aerosol concentration (e.g., aerosol-cloud interaction -ACI).
The magnitude and sensitivity of ACIs in low-level clouds have been investigated by numerous studies, using various observational datasets such as ground-based measurements (Garrett et al., 2004;Kim et al., 2008;McComiskey et al., 2009;Wang et al., 2013aWang et al., , 2018a, satellite retrieved products (Sekiguchi et al., 2003;Su et al., 2010), and airborne in situ measurements (Twohy et al., 2013;Painemal and Zuidema, 2013;Zhao et al., 2018). However, large variations exist among various assessments because of intrinsic instrument uncertainty, differing analysis methods, and, more physically, the inherent variation in aerosol properties. The physical mechanism underlying the aerosol effect on clouds is that aerosols activate as cloud condensation nuclei (CCN) and then influence the cloud microphysical features. The efficacy of the activation of aerosol has been widely known to be influenced by aerosol size distribution and chemical composition, which are the primary sources of uncertainty in assessing ACI (Dusek et al., 2006;McFiggans et al., 2006;Liu and Li, 2014;Che et al., 2016).
Previous studies have suggested that the composition of aerosols can be inferred by their optical properties such as aerosol optical depth, single-scattering albedo, and the Ångström exponent (Clarke et al., 2004(Clarke et al., , 2007Bergstrom et al., 2007;Russell et al., 2010;Cazorla et al., 2013;Cappa et al., 2016). For instance, fine-mode carbonaceous particles (e.g., black and organic carbon) have strong light-absorbing abilities in the ultraviolet and visible spectra (Logan et al., 2013). Urban pollution aerosols associated with sulfate and nitrate particles are considered to be weakly absorbing aerosols (Eck et al., 1999(Eck et al., , 2005Bergstrom et al., 2007;Mian Chin et al., 2009). Although studies have been done to classify aerosol types using the absorption Ångström exponent, which is associated with the absorptive spectral dependence of particles, the measurements of this parameter typically carry large uncertainty and can provide limited information when there are mixtures of different aerosol species that share similar spectral dependences (Bergstrom et al., 2007;Lack and Cappa, 2010). Alternatively, the single-scattering albedo (SSA) and co-albedo (1-SSA) can be used to better separate the aerosol types because they focus on the relative absorbing ability of aerosols at specific wavelengths (Logan et al., 2013;Tian et al., 2017). Given the wide availability of aerosol optical property measurements, the feasibility of inferring aerosol species from their optical properties is useful particularly in areas with no direct measurements of aerosol chemical composition (Logan et al., 2013;Schmeisser et al., 2017).
The Atmospheric Radiation Measurement (ARM) program initiated by the US Department of Energy (DOE) aims to improve the parameterization of clouds in global climate models (Stokes and Schwartz, 1994). Thus far, the ARM program has established over 20 years of long-term ground-based measurements of cloud properties and surfacemeasured aerosol properties at the southern Great Plains (SGP) site, which represents typical continental conditions (Ackerman and Stokes, 2003;Dong et al., 2005). The size and composition of aerosols have been found to have a considerable seasonal and regional dependence, and their impacts on clouds also vary with different aerosol regimes (Sorooshian et al., 2010;Logan et al., 2018). The prevailing fine-mode aerosols at the ARM SGP site typically contain organic and black carbon associated with biomass burning and inorganic aerosols composed of sulfate and nitrate species (Parworth et al., 2015;Logan et al., 2018). The differences in intrinsic hygroscopicity among those aerosol species play various roles in aerosol activation processes and consequently lead to various interactions with clouds. Thus, it is necessary to investigate the aerosol and cloud properties as well as the magnitude of the ACI index at the ARM SGP site in order to (a) enhance the understanding of the ACI and (b) reduce the uncertainty in quantifying the ACI and associated radiative effects when modeling aerosol influences on low-level continental clouds.
In this study, the aerosol and cloud properties at the ARM SGP site from 16 selected non-precipitating low-level stratiform cloud cases during the 2007-2012 period are examined. Details of the observational measurement platforms and methods are introduced in Sect. 2. The development and analysis of the ACI for the 16 selected cases, the aerosol activation and cloud microphysical responses, and consequent cloud radiative effects under different aerosol absorptive regimes are investigated in Sect. 3. Lastly, a summary of our findings and future work is presented in Sect. 4.

Cloud boundaries
The cloud boundaries at the ARM SGP site were primarily determined by the ARM Active Remotely-Sensed Cloud Locations (ARSCL) product, which is a combination of data detected by multiple active remote-sensing instruments, in particular, the millimeter-wavelength cloud radar (MMCR). The MMCR operates at a frequency of 35 GHz (and wavelength of 8.7 mm) with a zenith-pointing beamwidth of 0.2 • and provides a continuous time-height profile of radar reflectivity with temporal and spatial resolutions of 10 s and 45 m, respectively (Clothiaux et al., 2000). After 2011, the MMCR was replaced by the Ka-band ARM zenith radar (KAZR), which has the same operating frequency and shares similar capabilities to the MMCR, but with the major improvement of a new receiver that allows for more sensitivity in cloud detection (Widener et al., 2012). The temporal and vertical resolutions of KAZR-detected reflectivity are 4 s and 30 m, respectively. The cloudy condition, as well as cloud top height, is identified via cloud radar reflectivity. The uncertainties of cloud top height detected by MMCR and KAZR are 45 and 30 m, respectively.
The cloud radar is sensitive to the sixth moment of droplet size distribution and can be contaminated by insects below the cloud base . The laser ceilometer measurement, which is sensitive to the second moment, is used to provide an accurate cloud base estimation. The uncertainty of cloud base height is around 10 m (Morris, 2016). Hence, the lidar-radar pair provides the most precise determination of cloud boundaries from a point-based perspective, with combined uncertainties of cloud thickness for MMCR and KAZR periods being 55 and 40 m, respectively. Note that this will not cause a significant difference in determining the cloud boundaries between these two radar periods. In this study, the cloud base and top heights were averaged into 5 min bins, where the low-level stratus cloud is defined as a cloud-top height lower than 3 km with no overlying cloud layer .

Cloud microphysical properties
The cloud liquid water path (LWP), defined as the columnintegrated cloud liquid water, was retrieved based on the measured brightness temperatures from the microwave radiometer (MWR) at 23.8 and 31.4 GHz, using the statistical method described in Liljegren et al. (2001). The uncertainty of LWP retrieval is 20 g m −2 for an LWP less than 200 g m −2 and around 10 % for an LWP higher than 200 g m −2 . In this study, we exclude the data points with LWPs less than 20 g m −2 to eliminate optically thin clouds as well as exclude the samples with LWPs greater than 300 g m −2 to prevent potential precipitation contamination issues (Dong et al., 2008).
For microphysical properties of low-level stratus, following the methods developed by Dong et al. (1998), the daytime layer-mean cloud-droplet effective radius (r e ) can be calculated by where γ is the solar transmission, µ 0 is the cosine of solar zenith angle, and the units of r e and the LWP are micrometers (µm) and 100 g m −2 , respectively. N d is obtained after r e and is known by the following calculation: where N d is number per cubic centimeter (cm −3 ), Z is cloud thickness determined from cloud boundaries with units of meters, and σ x is the width of the lognormal size distribution of cloud droplets, which is assumed to be a constant value of 0.38 (Miles et al., 2000). The sensitivities of retrieved r e and N d to the uncertainties of the cloud LWP, σ x , and γ have been investigated in Dong et al. (1997Dong et al. ( , 1998. The uncertainties of retrieved r e and N d have been estimated against aircraft in situ measurements over the ARM SGP site (Dong et al., 2002;Dong and Mace, 2003) and other regions (Dong et al., 1998). As a result, the 10 % change in the cloud LWP and downward SW at the surface would cause a 10 % uncertainty in r e retrieval. In addition, the N d uncertainty is statistically estimated to be 25 %, compared with the aircraft in situ measurements at the Pennsylvania State University surface site during fall 1996 (Dong et al., 1998) and at the ARM SGP site during the March 2000 cloud intensive observational period (IOP; Dong et al., 2002;Dong and Mace, 2003).

Aerosol properties
Surface aerosol properties were collected from the aerosol observing system (AOS), a platform consisting of an array of instruments to monitor real-time aerosol information. The total condensation nuclei number concentration (N a ), which represents the overall loading of aerosol particles with diameters larger than 10 nm, was obtained by the TSI model 3010 condensation particle counter. The aerosol scattering coefficient (σ sp ) was measured by the TSI model 3653 nephelometer at three wavelengths: 450, 500, and 700 nm. The relative humidity inside the nephelometer was set to 40 % to maintain a dry condition and prevent potential aerosol hygroscopic effects (Jefferson, 2011). Moreover, the quality of retrievals was assured using the Anderson and Ogren (1998) method. The absorption coefficient (σ ap ) was measured by the Radiance Research particle soot absorption photometer (PSAP) at three slightly different wavelengths (470, 528, and 660 nm), with the calibration and quality control process done by the method developed in Anderson et al. (1999). Note that both the nephelometer and PSAP employ two impactors with size cuts of 1 and 10 µm. The measurements switch between total aerosol (< 10 µm) and submicron aerosol (< 1 µm) every hour. In this study, the sub-10 µm aerosol optical properties with original 1 min temporal resolution were averaged into 5 min bins to match the cloud microphysical properties. The optical particle counter developed by Droplet Measurement Technologies is used to measure the CCN number concentration (N CCN ). The supersaturation (SS) level inside the instrument cycles between 0.15 % and 1.15 % every hour. The CCN activity can be presented as a function of SS: N CCN = cSS k (Twomey, 1959), where c and k are calculated by using a power-law fit for each hour. In this study, 0.2 % is used, as this represents typical supersaturation conditions of low-level stratus clouds (Hudson and Noble, 2013;Logan et al., 2014Logan et al., , 2018.

Boundary-layer condition and lower-tropospheric stability
Given the fact that the aerosol properties were measured at the surface, there is a question of whether surface aerosols can be linked to what actually happens in clouds aloft. This study adopts the method presented in Dong et al. (2015), which sorts the boundary-layer condition into two categories: coupled and decoupled. The vertical sounding profiles at a 1 min temporal resolution were collected from the ARM Merged Sounding product, with a vertical resolution of 20 m below 3 km (Mace et al., 2006;Troyan, 2012). The vertical profiles of liquid water potential temperature (θ L ) and the total water mixing ratio (q t ) for coupled and decoupled boundary-layer conditions, as well as the criteria to differentiate between them, are illustrated in Fig. 1. The coupled condition was identified by the change of θ L and q t from the surface layer to cloud base of less than 0.5 K and 0.5 g kg −1 , respectively. These thresholds are the same as in Dong et al. (2015), originally suggested by Jones et al. (2011). In that case, the boundary layer is considered to be well-mixed and suggests that the surface aerosols are comparable to incloud aerosols. In the decoupled condition, θ L and q t vary more drastically from the surface to the cloud base under decoupled conditions, which denotes a stratification of the subcloud layer, thereby disconnecting the surface aerosols from the ones aloft.
A study that was conducted by Delle Monache et al. (2004) used in situ aerosol measurements from 59 flights from March 2000 to March 2001 to compare with the surface aerosol measurements. Their results showed that the daytime layer-mean aerosol extensive properties such as the total extinction by particles measured within the well-mixed boundary layer were well-correlated with surface measurements (R 2 value of 0.88). Therefore, selecting cloud cases under coupled conditions can better constrain the thermodynamic condition, since the measured surface aerosols are representative in terms of aerosol-cloud interaction.
The lower-tropospheric stability (LTS), which is defined as the potential temperature difference between the surface and 700 hPa, is used to represent the large-scale thermodynamic condition. The LTS is obtained from the ECMWF model output, which specifically provides for analysis at the ARM SGP site. The value is obtained by averaging over a grid box of 0.56 • ×0.56 • , which is centered at the SGP. The original temporal resolution of LTS is 1 h and is then interpolated to 5 min to match the other variables, assuming that the large-scale forcing would not have significant changes during every 1 h window.

Shortwave radiation fluxes at the surface
The surface-measured broadband downwelling shortwave (SW) radiation fluxes and estimated clear-sky SW fluxes were collected from radiative flux analysis value-added prod-ucts (Long and Ackerman, 2000;Long and Turner, 2008), with an uncertainty of 10 W m −2 . The combination of cloudy and clear-sky SW fluxes were used to calculate the cloud radiative effect. In order to minimize the influence of non-cloud factors, such as the solar zenith angle and surface albedo, a representation of the relative cloud radiative effect (rCRE) is defined as where SW dn cld and SW dn clr are cloudy and clear-sky downwelling shortwave radiation fluxes, respectively (Betts and Viterbo, 2005;Vavrus, 2006;Liu et al., 2011).

Selection of low-level stratus cloud cases
As previously discussed, the selection of cloud cases is limited by the following criteria: non-precipitating and cloudtop height less than 3 km with a lifetime of more than 3 h under the limitation of 20 g m −2 < LWP < 300 g m −2 and the coupled boundary-layer conditions. Only daytime cloudy periods were considered in this study because the r e retrieval required the information of solar transmission (Dong et al., 1998). Note that all the variables used in the study are averaged into 5 min temporal resolution bins. A total of 16 cases were selected during the 6-year period from 2007 to 2012, which represents a total of 693 samples (∼ 58 h) in this study. The detailed period and the number of sample points of each case are listed in Table 1. Most cases occurred during the winter and spring months, since low-level cloud occurrences are higher during those seasons . The 72 h NOAA HYSPLIT backward trajectories (Stein et al., 2015) for sub-cloud air parcels that advected over the ARM SGP site are used to identify the aerosol source regions (Logan et al., 2018). Aerosol plumes consisting of different species from local sources and long-range transport can impact the ARM SGP site because of different transport pathways and can induce different cloud responses, which are further investigated in this study.

Aerosol and cloud properties of selected cases
The probability density functions (PDFs) of aerosol and cloud properties from all 16 cases are shown in Fig. 2; note that the distributions include each of the 5 min data points. For the aerosol properties shown in Fig. 2a-d, the Ångström exponent (AE) was calculated based on the nephelometer-aerosol extensive properties observed spectral scattering coefficient (σ sp ) at 450 and 700 nm, using the equation AE 450−700 nm = − log σ sp450 /σ sp700 / log(450/700). The negative log-log slope denotes the relative wavelength dependence of particle optical properties due to differences in particle sizes (Schuster et al., 2006). Therefore, the AE can be a good indicator Figure 1. Vertical profiles of liquid water potential temperature (θ L ) and total water mixing ratio (q t ) for coupled (a) and decoupled (b) boundary-layer conditions. Blue lines denote cloud top and base heights, respectively. of aerosol particle sizes, since AE > 1 indicates the particle size distributions dominated by fine-mode aerosols (submicron), while AE < 1 denotes the dominance of coarse-mode aerosols (Gobbi et al., 2007;Logan et al., 2010). The aerosol fine-mode fraction (FMF) is given by the ratio σ sp1 /σ sp10 , where σ sp1 and σ sp10 are the nephelometer-measured scattering coefficients at 550 nm for fine-mode aerosols (1 µm size cut) and total aerosols (10 µm size cut), respectively. This ratio indicates the dominant influence of fine-mode aerosols, owing to the physical properties of the entire aerosol plume. For example, FMF values greater than 0.6 represent the dom-inance of fine-mode aerosol in the total population, and values less than 0.2 represent the dominance of coarse-mode aerosols in the total population (Anderson et al., 2003). As illustrated in Fig. 2b and c, fine-mode aerosols are dominant in the 16 selected cases. All AE values are higher than 1, with most of the values ranging from 1.5 to 2. In addition, the majority of the FMF values are greater than 0.6 and range from 0.7 to 0.9.
The variation in aerosol single-scattering albedo (SSA) suggests different roles of the fine-mode aerosol absorptive properties that influence total light extinction, which in turn is a result of different aerosol species in the plume. This is further explained in Sect. 3.3. The distributions of N a , N CCN , and N d represent typical continental aerosol conditions with mean values of 1060, 475, and 297 cm −3 , respectively, and r e values are more normally distributed, with the majority of values between 7 and 9 µm. Note that the variation in the PDF of the LWP is relatively small, which allows for a better investigation of the LWP dependence of cloud microphysical properties.

Measured aerosol-cloud interaction
To examine the microphysical response of cloud to aerosol loading, the quantitative ACI term can be expressed as where α denotes aerosol loading. ACI r represents the relative change of layer mean r e with respect to the relative change of aerosol loading, thereby emphasizing the sensitivity of the cloud microphysical response (Feingold et al., 2003;Garrett et al., 2004). Note that values of ACI r have theoretical boundaries of 0-0.33, where the lower bound means no change of cloud microphysical properties with aerosol loading, and the upper bound indicates a linear relationship. As suggested by previous studies, the ACI r should be calculated and compared at a constant LWP, owing to the dependence of r e on the LWP (Twomey, 1977;Feingold et al., 2003). Therefore, in this study, we use six LWP bins ranging from 0 to 300 g m −2 with a bin size of 50 g m −2 and then group the sample data accordingly. Note that the first bin is actually 20-50 g m −2 due to the elimination of LWPs less than 20 g m −2 . The r e -N CCN relationship is presented in Fig. 3a, where only the samples from three LWP bins are used to illustrate the r e -N CCN response. In general, r e decreases with increasing CCN number concentration as expected. The ACI r values range from 0.09 to 0.24, with a mean value of 0.145 ± 0.05; the uncertainty of ACI r corresponds to the 95 % confidence interval. Note that the ACI r values from six LWP bins show a generally decreasing trend of ACI r with an increasing LWP (Fig. 3b). Particularly, this decreasing trend is more obvious in a range of LWPs that are less than 150 g m −2 . The higher values of ACI r at lower LWPs indicate that the clouds are more susceptible to aerosol loading under lower liquid water availability. Given that the ACI r describes the response of r e to N CCN change, under low-LWP conditions with more CCN entering the cloud, the smaller particles compete against each other for the limited water supply and cannot efficiently grow into larger sizes. In that case, the higher CCN loading could result in smaller r e , and thus the variable range of r e is relatively broad, which is reflected by enhanced ACI r . Under high-LWP conditions, typically associated with sufficient water supply, the newly activated cloud droplet can grow larger quickly via condensation. However, the efficacy of condensational growth de-creases with enlarged particle size. The enhanced condensational growth under high-LWP conditions can shift the clouddroplet population to larger sizes. Therefore, for a similar CCN perturbation, the variable range of r e is narrower, which is reflected by reduced ACI r .
Previous studies have focused on the aerosol-cloud interaction in stratocumulus clouds at the ARM SGP site. Based on the analysis of seven selected stratocumulus cases during the period -2000, Feingold et al. (2003 reported the first ground-based measured ACI r values of 0.02 to 0.16 using the lidar measured aerosol extinction at a wavelength of 355 nm as the proxy for aerosol loading. A later study conducted by  assessed the ACI r using different aerosol measurements as CCN proxies in three selected stratus cases during the intensive operation period in May 2003. They found that the ACI r values were unrealistic when using N a to represent CCN loading, while using the surface aerosol scattering coefficient (σ sp ) and aerosol extinction at an altitude of 350 m as CCN proxies yields similar ACI r values, ranging from 0.14 to 0.39 . A recent study conducted by Sena et al. (2016) within the SGP region showed the different methodologies in calculating ACI r . In particular, different retrieval methods of r e could induce large differences. Moreover, the assessment of ACI r can be largely affected by the usage of different aerosol measurements that served as CCN proxies due to their own characteristics. Aerosol scattering and extinction coefficients are known to be relatively reliable CCN proxies, since they are more sensitive to aerosols that have larger particle sizes. As for N a , which represents the concentration of aerosol par- ticles with diameters larger than 10 nm, it is likely to pick up the very small aerosols generated by new particle formation events. This proportion of aerosols is presumably hard to activate as CCN, so it would not be counted in N CCN , especially under the 0.2 % supersaturation used in this study. Hence, it is less representative to use N a to accurately represent N CCN without the prior knowledge of the aerosol capacity to activate as CCN. Therefore, the usage of N CCN in this study is favorable for yielding a more straightforward assessment of ACI r , since the CCN measurement directly represents the number of aerosol droplets that already activated and have the potential of further growth.
In order to better understand the aerosol particle activation process in typical continental low-level stratus clouds, the ratios between N CCN and N a are examined and used to represent the aerosol activating capacities in the latter part of this study. Aerosol activating capacities greatly depend on size and composition. In order to further examine the role of aerosol species in the aerosol activation process and the potential impact on ACI r , the samples from the 16 selected cases are divided into two groups according to their absorptive regime, which is discussed in the following section.

Aerosol absorptive properties of the 16 selected cases
The measured absorptive properties of aerosols can aid in inferring the general information of different aerosol species, since different types of aerosols can demonstrate different absorptive behaviors at certain wavelengths. Aerosol plumes dominated by organic carbonaceous particles tend to represent strong absorptive capabilities in the visible spectrum but weakly absorb in the near-infrared spectrum (Dubovik et al., 2002;Lewis et al., 2008), while black carbon particles (e.g., soot) absorb across the entire solar spectrum with a weak dependence on wavelength (Schuster et al., 2005;Lack and Cappa, 2010). However, when the aerosol plume is dominated by anthropogenic inorganic pollution, the absorbing ability becomes even weaker (Clarke et al., 2007), partly due to sulfate chemical species (Mian Chin et al., 2009). Therefore, the general existence of carbonaceous and pollution particles can be inferred via absorptive properties.
In this study, we adopt the classification method involving the AE and the ratio of the aerosol absorption coefficient to the total extinction coefficient or single-scattering co-albedo, (ω abs = σ abs / (σ abs + σ scat )) (Logan et al., 2013(Logan et al., , 2014. This parameter is more sensitive to the capabilities of aerosol light absorption (rather than scattering) to total aerosol light extinction and therefore can better infer the aerosol composition (Logan et al., 2013). The ω abs values at a wavelength of 450 nm along with the AE 450−700 nm of all the samples are shown in Fig. 4. A ω abs value of 0.07 is used as a demarcation line of aerosols that are weakly and strongly absorbing. This value was determined using a frequency analysis performed at four AERONET sites that are dominated by single aerosol modes (Logan et al., 2013). Of the 16 cases, six cases are dominated by strongly absorbing aerosols, six cases are dominated by weakly absorbing aerosols, and four cases have samples which broadly scatter across the ω abs domain, which denotes a mixture of different absorbing aerosol species.
Within the 693 selected samples, 360 data points are classified in the weakly absorptive aerosol regime, while the remaining data points are in the strongly absorptive aerosol regime. It is interesting to note that the majority of the winter cases are dominated by weakly absorbing aerosols, while most of the spring cases exhibit a strongly absorbing aerosol dominance, which suggests that the aerosol plumes over the SGP site also have a seasonal dependence. In spring, owing to the upper-level ridge centered over the western Atlantic, the SGP is located at the northwestern edge of the sub-tropical high. Under this synoptic pattern, the SGP is under the influence of relatively frequent southerly transport of the air masses from Central America, which is characterized by strongly absorbing carbonaceous aerosols produced from biomass burning, as well as the moisture transported from the Gulf of Mexico. During the winter, the SGP site experiences the transported air masses from higher latitudes, with less intrusion of air masses from the south (Andrews et al., 2011;Parworth et al., 2015;Logan et al., 2018).

Aerosol and cloud properties under different absorptive regimes
Figure 5a-c show the PDFs of total N a , N CCN , and AE for the two absorptive regimes classified by ω abs . The distributions of N a from the two absorptive regimes are comparable to one another. The mean N CCN for the weakly absorptive regime (559 cm −3 ) is larger than that from the strongly absorptive regime (384 cm −3 ), and the occurrence of high N CCN values (larger than 1000 cm −3 ) is also higher in the weakly absorptive regime. This suggests different responses of CCN concentration to aerosols that have similar magnitudes but different absorptive properties. The AE distributions suggest dominant fine-mode aerosol contributions for both regimes. As for the cloud microphysical property distributions, cloud samples between the two regimes exhibit different characteristics ( Fig. 5d-f). The numbers above the bars in LWP distribution (Fig. 5d) for the two absorptive regimes denote the number of data points which will be used in the analysis with the binned LWP in the later sections. Cloud LWPs and r e values under the strongly absorptive regime have larger values, which contrasts with those under the weakly absorptive regime. On average, the weakly absorbing regime has higher N d and smaller r e (374 cm −3 and 6.9 µm, respectively) compared to the strongly absorbing regime (214 cm −3 and 8.2 µm). Note that the LWPs under the strongly absorptive regime are generally higher than those under the weakly ab-sorptive regime. This LWP difference might be associated with the seasonality of air mass transport over the SGP as discussed in Sect. 3.3.1. Although the seasonality of aerosol distribution and the LWP has similar trends, no clear causality has been found between them. Thus, the question behind these results is whether the differences in cloud microphysical properties between the two regimes are due to the difference in LWPs. As previously stated by Dong et al. (2015), cloud droplets generally grow larger at higher LWPs, which eventually leads to lower droplet number concentration.

Relationship of aerosol activating as CCN under different absorptive regimes
The measured N a and N CCN under the strongly and weakly absorbing aerosol regimes are plotted in Fig. 6. Note that N a samples from both regimes cover a broad range of values from 200 to 3500 cm −3 , suggesting a wide variety of aerosol loading conditions. These highly overlapping distributions allow a quantitative comparison between the ratios of N CCN to N a . For a broad range of N a , especially 200-500 and 1100-3500 cm −3 , the majority (∼ 74 %) of sample points from the strongly absorbing regime are located below the samples from the weakly absorbing regime. The linear regressions (95 % confidence level) between N CCN and N a for two regimes demonstrate the sensitivity of CCN 0.2 % SS to total aerosol loading. Note that the slope derived from the weak regime is slightly steeper than the strong regime, indicating that the N CCN values in the weakly absorptive regime increase faster than in the strongly absorptive regime with the same amount of aerosol increment. On average, 54 % of weakly absorbing aerosols can effectively activate as CCN compared to 45 % of the strongly absorbing aerosols. Note that those ratios are computed for an observed supersaturation level of 0.2 %. The fraction of aerosols that can activate as CCN increases with an increase in the supersaturation level, under the same aerosol size and composition condition (Dusek et al., 2006). A sensitivity test of how the aerosol activation ratio varies with different supersaturation levels is done by first interpolating the N CCN from 0.2 % to 1.15 % and then calculating N CCN /N a . As a result, the ratios of N CCN /N a for the weakly absorptive regime range from 0.54 to 0.38, while the ratios for the strongly absorptive regime range from 0.45 to 0.25. Considering that a supersaturation of 1.15 % in continental boundary-layer stratus is nearly impossible to reach, the supersaturation level of 0.2 % used in this study, which represents the most typical condition for continental low-level stratus, yields reasonable results.
Although it is generally considered that the role of aerosol particle size distribution is more important than the chemical component in terms of becoming CCN (Seinfeld and Pandis, 2012;Dusek et al., 2006), many studies have found that aerosol chemical composition can also have a non-negligible impact on the aerosol activating ability under polluted and low-supersaturation conditions (Rose et al., 2011;Che et   al., 2016). According to Köhler theory, the critical level of supersaturation for aerosol activation depends on the aerosol solubility, which decreases with increasing soluble particle number concentration. Hence, the role of aerosol chemical composition is more important at lower supersaturation and diminishes with increasing supersaturation levels (Zhang et al., 2012).
As discussed in Sect. 3.3.1, both weakly and strongly absorptive regimes are linked to aerosol plumes that are dominated by pollution and carbonaceous aerosols, respectively. Therefore, the difference in the ability of aerosol activation between the two regimes can be explained by the different hygroscopicity factors of the particle types. For example, anthropogenic pollution is associated with inorganic particles that are highly hygroscopic and have great ability in taking up water (Hersey et al., 2009;Massling et al., 2009;, while carbonaceous species (e.g., black and organic carbon) exhibit varying degrees of hygroscopicity, with species dominated by hydrophobic soot and black carbon being the least hygroscopic (Shinozuka et al., 2009;Rose et al., 2010). Thus, for the given amount of aerosol loading, aerosols in the weakly absorptive regime can better attract water vapor molecules and result in more aerosol particles activating as CCN.
As shown in Fig. 6, for three N a ranges (200-500, 500-1100, and 1100-3500 cm −3 ), the strongly absorbing aerosols show different relationships compared to weakly absorbing aerosols. The mean N CCN /N a values for those three N a ranges for weakly absorptive regimes are 0.77, 0.58, and 0.42, respectively, while the mean N CCN /N a values for the strongly absorptive regimes are 0.35, 0.51, and 0.32, respectively. This phenomenon is due to the mixed effect of aerosol composition (inferred by absorbing ability), aerosol size, and water availability on the aerosol activation. In the 200-500 cm −3 N a range, where the samples from the two absorptive regimes are most separated, the mean values of LWPs (158 g m −2 and 162 g m −2 for weakly and strongly absorptive regimes, respectively) indicate relatively sufficient water availability with a smaller aerosol concentration. In addition, X. Zheng et al.: Aerosol-cloud interactions under different absorptive aerosol regimes the weakly absorbing aerosol sizes are larger (AE = 1.59) than the strongly absorbing aerosol (AE = 1.73). It is known that larger aerosol particles easily activate under the same composition (Dusek et al., 2006), considering that the weakly absorbing aerosols are more hydrophilic, and thus the largest activation ratio difference among these three ranges is to be expected. The samples in the 500-1100 cm −3 N a range have AE values near 1.40 (1.53) for strongly (weakly) absorbing aerosols. The LWP for the strongly absorptive regime is 167 and 138 g m −2 in the weakly absorptive regime. Hence, the combined effect of larger particles and more water in the strongly absorptive regime leads to a N CCN /N a ratio close to the N CCN /N a ratio in the weakly absorptive regime. The samples in the 1100-3500 cm −3 N a range exhibit a smaller (AE = 1.67 and 1.57 for weakly and strongly absorptive regimes, respectively) aerosol particle size and less water availability (LWP = 95 g m −2 and 127 g m −2 for weakly and strongly absorptive regimes, respectively), which results in the lowest activation ratio (N CCN /N a = 0.42 and 0.32 for weakly and strongly absorptive regimes, respectively) among the three ranges for both regimes.
Due to the lack of detailed chemical observations for all the cloud sample periods, as well as the uncertainties among aerosol optical and microphysical properties induced by aerosol transformation processes such as aging and mixing (Wang et al., 2010(Wang et al., , 2018b, the bulk activation ratios revealed from this study cannot be significantly distinguished from each other. However, the effect of different aerosol species inferred by the absorptive properties with respect to aerosol activation is evident, especially at the 0.2 % supersaturation level. Furthermore, in the following section, the values of N CCN /N a and AE are sorted by the LWP for the two absorptive regimes in order to rule out the influence of the LWP and AE on aerosol activation.

LWP dependence of aerosol and CCN activation under different absorptive regimes
In order to better understand the role of aerosol activation ability in the microphysical process from aerosol to CCN and then to cloud droplets, comparisons must be considered under similar available moisture conditions due to the discrepancy in the LWP between the two regimes. Accordingly, the sorted N a values by the stratified LWP are presented in Fig. 7a along with the activation ratios of N CCN /N a , which are denoted by solid lines. For a range of LWPs from 20 to 300 g m −2 , the ratios of N CCN /N a under both regimes increase slightly with an increased LWP. In addition, all binned N CCN /N a values from the weakly absorptive regime (ranging from 0.4 to 0.6) are higher than those from the strongly absorptive regime (ranging from 0.3 to 0.5). A Student's t test is performed to test the ratio difference in each LWP bin at the 95 % significance level. The results indicate that the ratio differences between two absorptive regimes are statistically significant.
Taking the variation in N CCN into account, the activation ratios of N a to N CCN under low-LWP conditions (< 50 g m −2 ) in both regimes could be simply due to the linear combination of high aerosol concentration and insufficient moisture supply such that aerosols are competing against each other, thus resulting in a low activation ratio. However, as the LWP increases, the activation ratios tend to increase as well, especially at LWP values higher than 100 g m −2 . In fact, the values of N a in both regimes are relatively small, with little variation for LWP > 100 g m −2 , while the N CCN /N a ratio demonstrates a more noticeable increasing trend in the weakly absorptive regime. Despite higher aerosol loading in the strongly absorptive regime at large LWPs, there are still more weakly absorbing aerosols being activated, which corresponds to greater water uptake ability. Moreover, in every LWP bin, the AE value for the weakly absorptive regime is either higher than or very close to the AE value for the strongly absorptive regime (figure not shown). Even with relatively smaller particle sizes, under similar water availability, the weakly absorbing aerosol can better activate as CCN. In conclusion, a significant impact of aerosol composition on aerosol activation capacity, which is inferred by the aerosol absorbing capability, does exist.
As for the process from CCN to cloud droplets, a similar assessment is presented in Fig. 7b, which illustrates the N CCN values and the ratios of N d /N CCN in relation to the LWP. The ratios of N d /N CCN in the weakly absorptive regime range from 0.58 to 0.86, with a mean value of 0.69, and highly fluctuate with the LWP. In contrast, N d /N CCN in the strongly absorptive regime shows lower values and less variability (from 0.47 to 0.64), with a mean value of 0.54. It is interesting to note that the variation in N d /N CCN in the strongly absorptive regime mimics the variation in N CCN with the LWP, indicating a relatively lower aerosol-to-CCN activating capacity. Therefore, N d /N CCN shows no significant dependence on the LWP, which is consistent with previous studies, which suggest that the response of N d to the change in N CCN has no fundamental relationship with the LWP (e.g., McComiskey et al., 2009). In addition, the sensitivity and uncertainty of N d are examined in order to estimate the impact of N d uncertainty on the assessment of N d /N CCN . To assess the contributions of different input parameter uncertainties to N d retrieval, every input parameter was perturbed by its own uncertainty, with other parameters held fixed. The results are as follows: (a) an increase (decrease) in the LWP of 20 g m −2 leads to 27.9 % (27.6 %) change in N d while an increase (decrease) in σ x of 0.15 leads to a 50.8 % (23.9 %) change in N d , (b) an increase (decrease) in cloud thickness of 0.15 leads to a 14.5 % (23.2 %) change in N d , and (c) an increase (decrease) in r e of 10 % leads to a 14.5 % (23.2 %) change in N d . The percentage changes in N d due to different input uncertainties range from 14.5 % to 50.8 %, with the majority falling between 20 % and 30 %. Note that the largest uncertainty of N d happens when increasing σ x by 0.15. However, when considering that continental stratocumulus generally contains smaller droplets, one might expect their distribution width to be smaller than 0.38 (Dong et al., 1997). Therefore, the overall uncertainty of 25 % compared to the aircraft in situ measurement should be a reasonable estimation. In this case, the mean ratio of N d /N CCN for the weakly absorptive aerosol regime ranges from 52 % to 86 %, while the mean ratio of N d /N CCN for the strongly absorptive aerosol regime ranges from 41 % to 67 %.
The overall differences in CCN conversion fractions are likely a result of the combined effects of meteorological factors and the aerosol heating effect on the cloud environment. To examine the meteorological influence on CCN conversion, the LTS parameter is used to investigate the difference in the large-scale thermodynamic condition. By sorting the LTS by the LWP for the two absorptive regimes, the LWP dependence on LTS can be ruled out, which can provide a better understanding of the potential role of LTS in cloud-droplet development. For each given LWP bin, the weakly absorptive regime has higher LTS values than the strongly absorptive regime (figure not shown). The LTS is largely impacted by the potential temperature difference throughout the mixed layer. If a strong capping boundary-layer temperature inversion is present, it will result in high LTS values and, in turn, a well-mixed boundary layer (Wood and Bretherton, 2006). Such results indicate that even under similar available moisture conditions, the more sufficient turbulence can transport the below-cloud moisture as well as the CCN that activated from weakly absorbing aerosols into the cloud more efficiently, contributing to a higher ratio of N d /N CCN in the weakly absorptive regime. However, the LTS emphasizes a general thermodynamic condition in the lower troposphere with a wider domain as compared to the single-point measurement.
In addition, the vertical velocity in pressure coordinate (ω) values at the 925 hPa level, which represent the large-scale forcing on the vertical motion between surface and cloud layer, are also sorted by the LWP for the two absorptive regimes in order to check the potential influence of the environmental dynamic state (figure not shown). However, the ω values for both absorptive regimes share the same mean value of 0.031 Pa s −1 and show no dependence on the LWP, indicating that the large-scale environments over the SGP are generally dominated by sinking motion. The synoptic patterns of composite geopotential height for the two absorptive regimes show that the ARM SGP site is located ahead of the 700 hPa ridge and is located within the surface high pressure. The meteorological pattern is favorable for the generation of downward motion at the lower troposphere, and the sinking motion induces relatively stable environments in the lower troposphere, which is consistent with the LTS measurements. Considering the fact that the ω value is obtained from a relatively larger domain surrounding the SGP, it is difficult to reflect the true cloud-scale dynamics, especially the vertical velocity or turbulence strength at the cloud base. Therefore, the influence of cloud-scale dynamics, presumably cloudbase updraft, is non-negligible, since the sensitivity of cloud droplet to aerosol loading is enhanced with increasing updraft velocity as reported in previous studies (e.g., Feingold et al., 2003;McComiskey et al., 2009).

X. Zheng et al.: Aerosol-cloud interactions under different absorptive aerosol regimes
Furthermore, the heating effect of light-absorbing aerosols on the cloud environment cannot be neglected. Strongly light-absorbing aerosols can absorb solar radiation and heat the in-cloud atmosphere by emission, which results in the reduction of relative humidity (or supersaturation) in the cloud layer (Bond et al., 2013;Wang et al., 2013b). This effect is evident in the observation, as the values of in-cloud relative humidity in the strongly absorptive regime are slightly lower than those in the weakly absorptive regime. Additionally, this aerosol heating effect disrupts the boundary-layer temperature structure by enhanced warming aloft and, consequently, inhibits the vertical transport of sensible and latent heat between surface and cloud layer. The impacts of lightabsorbing aerosol on cloud-scale thermodynamics and the dynamics state might eventually dampen the conversion process from CCN to cloud droplets. Unfortunately, due to the lack of measurement of cloud-base vertical velocity throughout the studying period, this competing effect of cloud thermodynamic and dynamic cannot be fully untangled from the aerosol effect given the currently available dataset. The differences in N d /N CCN between the two regimes might be affected by the combined effects of LTS, updraft velocity, and aerosol heating on the cloud environment.

r e and N d dependence on LWP under different absorptive regimes
In the previous section, we examined the activation ratios of aerosol to CCN and then from CCN to cloud droplets between the two regimes as well as their dependencies on the LWP, which eventually led to the cloud-droplet variation for a given LWP range. Figure 7c and d demonstrate that r e increases while N d decreases with an increased LWP up to roughly 150 g m −2 in both regimes. Note that as LWPs are greater than 150 g m −2 , N d values in both regimes show less variation with the LWP, while r e values in the strongly absorptive regime also show little variation, which implies limited growth even with increasing water availability. However, the r e values in the weakly absorptive regime increase from 7.8 to 8.8 µm, which suggests that under a given number concentration, the cloud droplet can grow by continuing to collect moisture. As shown in each LWP bin, the r e values in the weakly absorptive regime are smaller than those in the strongly absorptive regime, while the N d values in the strongly absorptive regime are much lower than those in the weakly absorptive regime. The combination of cloud thermodynamic, dynamic, and light-absorbing aerosol heating effects impact the conversion process from CCN to cloud droplets. Under a given moisture availability, a higher number of CCN in the weakly absorptive regime can be converted to cloud droplets. This results in higher number concentrations of smaller cloud droplets, while the dampened CCN conversion process in the strongly absorptive regime leads to fewer and larger cloud droplets at a fixed LWP.

Aerosol-cloud interaction under different absorptive regimes
To examine the sensitivity of clouds to both weakly and strongly absorbing aerosol loading, the relationships between cloud r e and N CCN are shown in Fig. 8. Two LWP ranges (0-50 and 200-250 g m −2 ) are selected in order to better represent ACI r at low-and high-LWP conditions. For the examination of r e as a function of N CCN (Fig. 8a for low LWP range), the ACI r values in the weakly absorptive regime are higher than those in the strongly absorptive regime. This suggests that the cloud droplets are more sensitive to weakly absorbing aerosols than to strongly absorbing aerosols in clouds with low LWPs. In other words, if there is some increment in aerosol particles, clouds influenced by weakly absorbing aerosols will respond to this increment more effectively and decrease faster in droplet sizes relatively. Under high-LWP conditions (Fig. 8b), the ACI r values are lower and show less difference between the two regimes, which is in agreement with previous discussions on the sensitivity of cloud microphysical properties to aerosol loading. Based on the sensitivity study, the 10 % change of the cloud LWP and downward SW at the surface would result in the 10 % uncertainty in r e retrieval (Dong et al., 1997). When compared with aircraft in situ measurements, the differences between retrievals and in situ measurements are around 10 % (Dong et al., 1998(Dong et al., , 2002. In order to assess the impact of r e uncertainty on ACI r , we use the Monte Carlo method to propagate the r e uncertainty on ACI r by the following procedure. The r e value for each data point is randomly perturbed by ±10 %, and thus the corresponding ACI r can be recalculated based on the perturbed set of r e . After 100 000 iterations, we obtain a distribution of ACI r values. The uncertainty of ACI r is given by 1 standard deviation of those 100 000 values of ACI r ; since the distribution of ACI r values follows a normal distribution with a narrow peak, this uncertainty value represents the uncertainty in the computed ACI r due to errors in the r e retrieval. The uncertainties of ACI r for the two absorptive regimes are denoted as the dashed line in Fig. 8. In the lower LWP range (Fig. 8a), the ACI r uncertainty is 0.020 (0.030) for the weakly (strongly) absorptive regime; to account for the uncertainties, the difference in ACI r between the two absorptive regimes is preserved. In the higher LWP range (Fig. 8b), the ACI r uncertainty is 0.044 (0.023) for the weakly (strongly) absorptive regime, which is nonnegligible. Taking the uncertainties of ACI r into account, the ACI r in the two absorptive regimes cannot be separated well, owing to the enhanced condensational growth process accompanied by a higher LWP and the diminished cloud response to aerosols associated with different ω abs values. In general, the 10 % uncertainty in r e retrieval contributes to 0.02-0.04 in ACI r uncertainties.
Note that the LTS values from the weakly absorptive regime (22.91 and 19.78 K) are higher than those from the strongly absorptive regime (21.72 and 17.83 K) for the se- Figure 8. r e as a function of N CCN and the values of ACI r under the strongly absorptive (in red) and the weakly absorptive (in blue) aerosol regimes at two LWP bins: 0-50 g m −2 (a) and 200-250 g m −2 (b). Note that the dashed lines denote the uncertainties of ACI r due to 10 % error in r e retrieval. lected two LWP bins. As discussed in the previous section, owing to the stronger temperature inversion indicated by the higher LTS values, low clouds are more closely connected to weakly absorbing aerosols and moisture below cloud by sufficient turbulence. In order to quantify the impact of LTS on ACI r , we adapted the criteria described in Grysperdt et al. (2016) in which the LTS value of 18 K denotes the demarcation line between high-and low-LTS regimes and constrains the ACI r for the two regimes by their LTS values accordingly. Owing to the highly limited sample points that fall into the low-LTS category, the ACI r can only be constrained in the high-LTS condition. For the 0-50 g m −2 LWP range, the ACI r for the weakly absorptive regime increases from 0.26 to 0.31, and the ACI r for the strongly absorptive regime increases from 0.21 to 0.24. The enhancement effect of LTS on the ACI r is noticeable, which is in accordance with the previous discussion in which the high-LTS environment is associated with (a) sufficient turbulence in the boundary layer and (b) a closer connection between the surface and cloud layer, which enhances the cloud microphysical responses to the CCN. The result is consistent with the previous study by Kim et al. (2008), who found that ACI r is enhanced under adiabatic cloud conditions and higher LTS values are associated with higher cloud adiabaticity. Note that though ACI r values are increased for both regimes, the difference between the regimes becomes larger (from 0.05 to 0.07) because low-level stratus clouds are more susceptible to weakly absorptive aerosol. Furthermore, the enhancement of ACI r in the high-LTS environment is more evident in the weakly absorptive regime. In the case of the 200-250 g m −2 LWP range, the LTS effects on ACI r are less significant compared to the lower LWP range. No significant change in the weakly absorptive regime is evident, and the ACI r in the strongly absorptive regime decreased from 0.12 to 0.10, partly owing to the enhanced condensational growth process accompanied by a higher LWP, and thus inhibits the impact of LTS on ACI r . Overall, ACI r values are enhanced under high-LTS conditions, but the difference between the two regimes indicates that both ω abs and LTS can be the impact factor of the ACI r , but there is not necessarily causality between them.
Furthermore, with the presence of strongly light-absorbing aerosols, the cloud-layer heating induced by the aerosol absorptive effect can result in the reduction of in-cloud supersaturation and leads to the damping of cloud microphysical sensitivity to strongly absorbing aerosols. A previous modeling study conducted at the ARM SGP site by Lin et al. (2016) estimated the shortwave heating rates in cloud layers by contrasting the simulations with and without lightabsorbing aerosols. The inclusion of light-absorbing aerosols was represented by an internal aerosol mixture with a mass combination of 95 % ammonium sulfate and 5 % black carbon. The SSA of this mixture is calculated to be roughly 0.9, as documented in the previous study of Wang et al. (2014). The different values of SSA used in their study (0.9 for lightabsorbing and 1.0 for non-absorbing) are comparable to this study (0.89 for strongly absorbing and 0.97 for weakly absorbing). The induced increments in cloud-layer shortwave heating rates have a maximum value of 3 K d −1 , compared to the simulation without light-absorbing aerosols. Note that the aerosol number concentration in Lin et al. (2016) was set to 2800 cm −3 . To get a simple comparison with the aerosol number concentration in this study, one might expect that the light-absorbing-aerosol-induced cloud-layer shortwave heating rates can have a similar maximum increment, and the general increment should be about 1 K d −1 , which is non-negligible. The absorption of solar radiation by lightabsorbing aerosols warms the cloud layer as well as the boundary layer below it, which in turn stabilizes the lower troposphere and results in reduced cloud susceptibility. In general, the results indicate that the ACI r can be counteracted by the light-absorbing aerosol heating effect and be enhanced under a thermodynamic environment of high static stability, especially under lower-LWP conditions.

Cloud shortwave radiative effects under different absorptive regimes
Aerosols with different absorptive properties can alter the ability of clouds to reflect incoming shortwave radiation. Accordingly, cloud radiative effects on shortwave radiation for the two absorptive regimes are investigated. Both cloudy and clear-sky downwelling shortwave fluxes for samples in the weakly absorptive regimes are generally higher than those in the strongly absorptive regime (not shown in here), largely owing to the discrepancies in solar zenith angle, seasonal variation in insolation, and surface albedo. Therefore, to ensure that the comparison is under minimum influence of noncloud factors, the shortwave relative cloud radiative effects (rCREs) are introduced and their dependencies on the LWP between the two regimes are examined. With all else being equal, as shown in Fig. 9, rCREs in both regimes noticeably increase with the LWP, especially for LWPs less than 150 g m −2 . Using a fixed LWP, rCREs in the weakly absorptive regime are always higher than those in the strongly absorptive regime because the greater activating ability of the weakly absorbing aerosols leads to higher N d and smaller r e as opposed to the strongly absorbing aerosols. Thus, clouds with a larger amount of small cloud droplets contribute more to the extinction of incident solar radiation. The difference in mean rCRE between the two regimes is small but nonnegligible (∼ 0.04). Quantitatively speaking, taking the climatological downwelling solar flux of the winter season (∼ 150 W m −2 , Dong et al., 2006) as an example, the extinction of incident solar radiation by clouds that develop from weakly absorbing aerosols is 6.0 W m −2 more than those by clouds from strongly absorbing aerosols. From independent radiative measurements, the phenomenon that clouds are more susceptible to weakly absorbing aerosols is further evident.

Conclusions
A total of 16 non-precipitating overcast low-level stratiform cloud cases under daytime coupled boundary-layer conditions were selected in order to investigate the sensitivity of cloud microphysical properties to aerosol physicochemical properties. The Ångström exponent and fine-mode fraction distributions indicate that the aerosol plumes that advected to the SGP site during all the selected cases were dominated by fine-mode particles, while the variation in aerosol singlescattering albedo suggests different characteristics of optical properties among the aerosol plumes. In terms of the sensitivity of cloud droplets to aerosol number concentration, the values of ACI r range from 0.09 to 0.24, with a mean of 0.145 ± 0.05, which supports the finding of previous studies Figure 9. Relative cloud radiative effect (rCRE) as a function of liquid water path (LWP) under the strongly absorptive (in red) and weakly absorptive (in blue) aerosol regimes. Whiskers denote 1 standard deviation for each bin.
using ground-based measurements. The magnitude of ACI r shows a decreasing trend with an increasing LWP, partly owing to the enhanced condensational growth process accompanied by a higher LWP. Clouds that develop under lower-LWP conditions are more susceptible to aerosol loading, owing to the enhanced competition between small cloud droplets to grow larger with limited water supply. The analysis of the N CCN /N a ratio under the two regimes further demonstrates that weakly absorbing aerosols have statistically significant higher activation ratios (mean ratio of 0.54) than the strongly absorbing aerosols (mean ratio of 0.45). The fraction of weakly absorbing aerosols that activate as CCN shows a noticeable increase with an increased LWP, while the activation ratios for strongly absorbing aerosols tend to slightly increase with the LWP under comparable aerosol loading conditions. This is likely related to the hygroscopicity associated with the aerosol species. For example, weakly absorbing aerosols are typically dominated by pollution aerosols that have greater water uptake ability, while strongly absorbing aerosols are generally hydrophobic, such as freshly emitted black and organic carbon.
The ratios of N d /N CCN in the weakly absorptive regime (mean ratio of 0.67) are higher than for the strongly absorptive regime (mean ratio of 0.54). This is due to the higher-LTS environment for the weakly absorptive regime, enhancing the connection between cloud and the below-cloud moisture and CCN. In addition, the cloud-layer heating effect induced by the strongly light-absorbing aerosols results in the reduction of in-cloud supersaturation and leads to the damping of the CCN conversion process for the strongly absorptive regime. As a result, cloud droplets that form from weakly absorbing aerosols tend to have smaller sizes and higher concentrations than cloud droplets forming from strongly absorbing aerosols. Furthermore, the cloud droplets in the weakly absorptive regime exhibit a greater growing ability, as given by larger r e values that increase with the LWP under similar N d . The differences in cloud-droplet development between the two regimes are a likely result of the combination of thermodynamics, dynamics, and aerosol heating effects.
Under low-LWP conditions, the measured ACI r values in the weakly absorptive regime are relatively higher, indicating that clouds have greater microphysical responses to aerosols in the weakly absorptive regime than in the strongly absorptive regime. The favorable LTS condition in the weakly absorptive regime enhanced the cloud susceptibility. The cloudlayer heating effect of light-absorbing aerosol reduced the ACI r in the strongly absorptive regime. The observed ACI r is enhanced after being constrained by high LTS, particularly under lower-LWP conditions. Under higher-LWP conditions, the enhanced condensational growth process diminishes the LTS impact on ACI r , and the damping of ACI r is more evident, which is consistent with the results from all the cases. In general, the 10 % uncertainty in r e retrieval contributes to ACI r uncertainties ranging from 0.02 to 0.04 for the two absorptive regimes, with the ACI r difference between the two absorptive regimes still being well-preserved. As a result, clouds that develop from weakly absorbing aerosols exhibit a stronger shortwave cloud radiative effect than clouds originating from strongly absorbing aerosols. Additional future work will focus on investigating in detail the composition of different aerosol plumes with respect to their physicochemical properties. The aerosol-cloud-interaction processes under the influence of different aerosol types associated with air masses and the sensitivity to dynamic and thermodynamic factors will be further examined.
Author contributions. The original idea of this study was discussed by XZ, BX, and XD. XZ performed the analyses and wrote the paper. XD, TL, YW, and PW participated in further scientific discussion and provided substantial comments on and edits to the paper.
Competing interests. The authors declare that they have no conflict of interest.