Articles | Volume 21, issue 13
Atmos. Chem. Phys., 21, 10499–10526, 2021

Special issue: Marine aerosols, trace gases, and clouds over the North Atlantic...

Atmos. Chem. Phys., 21, 10499–10526, 2021

Research article 13 Jul 2021

Research article | 13 Jul 2021

Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors

Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors
Hossein Dadashazar1, David Painemal2,3, Majid Alipanah4, Michael Brunke5, Seethala Chellappan6, Andrea F. Corral1, Ewan Crosbie2,3, Simon Kirschler7, Hongyu Liu8, Richard H. Moore2, Claire Robinson2,3, Amy Jo Scarino2,3, Michael Shook2, Kenneth Sinclair9,10, K. Lee Thornhill2, Christiane Voigt7, Hailong Wang11, Edward Winstead2,3, Xubin Zeng5, Luke Ziemba2, Paquita Zuidema6, and Armin Sorooshian1,5 Hossein Dadashazar et al.
  • 1Department of Chemical and Environmental Engineering, University of Arizona, Tucson, AZ, USA
  • 2NASA Langley Research Center, Hampton, VA, USA
  • 3Science Systems and Applications, Inc., Hampton, VA, USA
  • 4Department of Systems and Industrial Engineering, University of Arizona, Tucson, AZ, USA
  • 5Department of Hydrology and Atmospheric Sciences, University of Arizona, Tucson, AZ, USA
  • 6Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL, USA
  • 7Institute of Atmospheric Physics, German Aerospace Center, Oberpfaffenhofen, Germany​​​​​​​
  • 8National Institute of Aerospace, Hampton, VA, USA
  • 9NASA Goddard Institute for Space Studies, New York, NY, USA
  • 10Universities Space Research Association, Columbia, MD, USA
  • 11Atmospheric Sciences and Global Change Division, Pacific Northwest National Laboratory, Richland, WA, USA

Correspondence: Hossein Dadashazar (


Cloud drop number concentrations (Nd) over the western North Atlantic Ocean (WNAO) are generally highest during the winter (DJF) and lowest in summer (JJA), in contrast to aerosol proxy variables (aerosol optical depth, aerosol index, surface aerosol mass concentrations, surface cloud condensation nuclei (CCN) concentrations) that generally peak in spring (MAM) and JJA with minima in DJF. Using aircraft, satellite remote sensing, ground-based in situ measurement data, and reanalysis data, we characterize factors explaining the divergent seasonal cycles and furthermore probe into factors influencing Nd on seasonal timescales. The results can be summarized well by features most pronounced in DJF, including features associated with cold-air outbreak (CAO) conditions such as enhanced values of CAO index, planetary boundary layer height (PBLH), low-level liquid cloud fraction, and cloud-top height, in addition to winds aligned with continental outflow. Data sorted into high- and low-Nd days in each season, especially in DJF, revealed that all of these conditions were enhanced on the high-Nd days, including reduced sea level pressure and stronger wind speeds. Although aerosols may be more abundant in MAM and JJA, the conditions needed to activate those particles into cloud droplets are weaker than in colder months, which is demonstrated by calculations of the strongest (weakest) aerosol indirect effects in DJF (JJA) based on comparing Nd to perturbations in four different aerosol proxy variables (total and sulfate aerosol optical depth, aerosol index, surface mass concentration of sulfate). We used three machine learning models and up to 14 input variables to infer about most influential factors related to Nd for DJF and JJA, with the best performance obtained with gradient-boosted regression tree (GBRT) analysis. The model results indicated that cloud fraction was the most important input variable, followed by some combination (depending on season) of CAO index and surface mass concentrations of sulfate and organic carbon. Future work is recommended to further understand aspects uncovered here such as impacts of free tropospheric aerosol entrainment on clouds, degree of boundary layer coupling, wet scavenging, and giant CCN effects on aerosol–Nd relationships, updraft velocity, and vertical structure of cloud properties such as adiabaticity that impact the satellite estimation of Nd.

1 Introduction

Aerosol indirect effects remain the dominant source of uncertainty in estimates of total anthropogenic radiative forcing (Boucher et al., 2013; Myhre et al., 2013). Central to these effects is knowledge about cloud drop number concentration (Nd), as it is the connection between the subset of particles that activate into drops (cloud condensation nuclei, CCN) and cloud properties. It is widely accepted that warm clouds influenced by higher number concentrations of aerosol particles have elevated Nd and smaller drops (all else held fixed), resulting in enhanced cloud albedo at fixed liquid water path (Twomey, 1977) and potentially suppressed precipitation (Albrecht, 1989) and increased vulnerability to overlying air resulting from enhanced cloud-top entrainment (Ackerman et al., 2004).

Reducing uncertainty in how aerosols and clouds interact within a given meteorological context requires accurate estimates of Nd and aerosol concentrations and properties. Since intensive field studies struggle to obtain broad spatial and temporal coverage of such data, satellite remote sensing and reanalysis datasets are relied on for studies examining intra- and interannual features over large spatial areas. Limitations of satellite retrievals are important to recognize. Nd is not directly retrieved but derived using other parameters (e.g., cloud optical depth, cloud drop effective radius, cloud-top temperature) and with assumptions about cloud adiabatic growth and Nd being vertically constant (Grosvenor et al., 2018). Aerosol number concentrations are usually represented by a columnar parameter such as aerosol optical depth (AOD) and thus not directly below clouds, which is the aerosol layer most likely to interact with the clouds. Furthermore, aerosol data are difficult to retrieve in cloudy columns. Reanalysis datasets circumvent issues for the aerosol parameters as they provide vertically resolved data (e.g., surface layer and thus below clouds) and are available for cloudy columns.

Of special interest in this work is the western North Atlantic Ocean (WNAO) where decades of extensive research have been conducted for topics largely unrelated to aerosol–cloud interactions (Sorooshian et al., 2020), thereby providing an opportunity for closing knowledge gaps for this area in a region with a wide range of aerosol and meteorological conditions (Corral et al., 2021; Painemal et al., 2021). Past work showed different seasonal cycles of AOD and Nd in this region (Grosvenor et al., 2018; Sorooshian et al., 2019), which partly motivates this study to unravel why Nd behaves differently on seasonal timescales. A previous study investigating seasonal cycles of Nd in the North Atlantic region found that cloud microphysical properties were primarily dependent on CCN concentrations while cloud macrophysical properties were more dependent on meteorological conditions (e.g., Sinclair et al., 2020). However, due to the complexity of interactions involved and the co-variability between individual components, the magnitude and sign of these feedbacks remain uncertain.

This study uses a multitude of datasets to characterize the Nd seasonal cycle and factors related to Nd variability. The structure of the results and discussion is as follows: (i) case study flight highlighting the wide range of Nd in wintertime and factors potentially affecting that variability; (ii) seasonal cycle of Nd and aerosol concentrations based on different proxy variables; (iii) seasonal cycles of factors potentially influential for Nd such as aerosol size distribution, vertical distribution of aerosol, humidity effects, and aerosol–cloud interactions; (iv) composite analysis of influential factors on high- and low-Nd days in each season; (v) modeling analysis to probe more deeply into Nd relationships with other parameters for winter and summer seasons; and (vi) discussion of other factors relevant to Nd unexplored in this work.

2 Methods

2.1 Study region

We focus on the WNAO, defined here as being bounded by 25–50 N and 60–85 W. A subset of the results focuses on six individual sub-domains representative of different parts of the WNAO (shown later), with five just off the East Coast extending from south to north (south: S; central-south: C-S; central: C; central-north: C-N; north: N) and one over Bermuda.

2.2 Datasets

2.2.1 Satellite observations (CERES-MODIS/CALIPSO)

Relevant cloud parameters were obtained from the Clouds and the Earth's Radiant Energy System (CERES) edition 4 products (Minnis et al., 2011, 2020), which are based on the application of CERES's retrieval algorithms on the radiances measured by the MODerate resolution Imaging Spectroradiometer (MODIS) instrument aboard the Aqua satellite. Aqua observations used to estimate Nd were from the daytime overpasses of the satellite around 13:30 local time (LT). Level 3 daily cloud properties at 1× 1 spatial resolution (listed in Table 1) were used for the period between January 2013 and December 2017 from CERES-MODIS edition 4 Single Scanning Footprint (SSF) products (Loeb et al., 2016). The CERES-MODIS SSF Level 3 product includes 1× 1 averaged data according to the cloud-top pressure of individual pixels: low (heights below 700 hPa), mid-low (heights within 700–500 hPa), mid-high (heights within 500–300 hPa), and high (heights above 300 hPa) level clouds. For this study, we only use low-cloud averages.

Table 1Summary of various data products used in this study.​​​​​​​

Download Print Version | Download XLSX

Nd is estimated based on an adiabatic cloud model (Grosvenor et al., 2018):

(1) N d = 5 2 π k f ad C w τ Q ext ρ w r e 5 1 / 2 ,

where τ is cloud optical depth and re is cloud drop effective radius, both of which are obtained from CERES-MODIS for low-level (i.e., surface to 700 hPa) liquid clouds. Qext is the unitless extinction efficiency factor, assumed to be 2 for liquid cloud droplets, and ρw is the density of water (1 g cm−3). Methods described in Painemal (2018) were used to estimate parameters in Eq. (1) as follows. (i) Adiabatic water lapse rate (Cw) was determined using cloud-top pressure and temperature provided by CERES-MODIS. (ii) The Nd estimation is often corrected for the sub-adiabatic profile by applying the adiabatic value (fad), but in this work, a value of fad=1 was assumed due to both lack of consensus on its value and its relatively minor impact on Nd estimation (Grosvenor et al., 2018). (iii) k is the parameter representing the width of the droplet spectrum and was assumed to be 0.8 over the ocean. Statistics of Nd are often estimated after screening daily observations based on cloud fractions (Wood, 2012; Grosvenor et al., 2018). The purpose of such filters is to reduce the uncertainties associated with the estimation of Nd (Eq. 1) driven by the errors in the retrieval of re and τ from MODIS's observed reflectance in a highly heterogeneous cloud field. However, this may inadvertently mask the effects of cloud regime on aerosol–cloud interactions by only including certain low-level cloud types in the analyses (e.g., closed-cell stratocumulus). Therefore, we use all Nd data regardless of cloud fraction with exceptions being Sects. 3.5 and 4.2 where a filter of low-level liquid cloud fraction (i.e., CFlow-liq. 0.1) was applied.

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) provides data on the vertical distribution of aerosols (Winker et al., 2009). Nighttime extinction profiles were acquired from Level 2 version 4.20 products (i.e., 5 km aerosol profile data), between January 2013 and December 2017. We averaged the Level 2 daily extinctions in different 4× 5 sub-domains (shown later) to obtain the seasonal profiles after applying the screening scheme outlined in Tackett et al. (2018).

2.2.2 MERRA-2

Aerosol data were obtained from the Modern-Era Retrospective Analysis for Research and Applications-Version 2 (MERRA-2) (Gelaro et al., 2017). MERRA-2 is a multidecadal reanalysis where meteorological and aerosol observations are jointly assimilated into the Goddard Earth Observation System version 5 (GEOS-5) data assimilation system (Buchard et al., 2017; Randles et al., 2017). Aerosols in MERRA-2 are simulated with a radiatively coupled version of the Goddard Chemistry, Aerosol, Radiation, and Transport model (GOCART; Chin et al., 2002; Colarco et al., 2010). GOCART treats the sources, sinks, and chemistry of 15 externally mixed aerosol mass mixing ratio tracers, which include sulfate, hydrophobic and hydrophilic black and organic carbon, dust (five size bins), and sea salt (five size bins). MERRA-2 includes assimilation of bias-corrected Collection 5 MODIS AOD, bias-corrected AOD from the Advanced Very High Resolution Radiometer (AVHRR) instruments, AOD retrievals from the Multiangle Imaging SpectroRadiometer (MISR) over bright surfaces, and ground-based Aerosol Robotic Network (AERONET) direct measurements of AOD (Gelaro et al., 2017). In this study we used total and speciated (i.e., sea salt, dust, black carbon, organic carbon, and sulfate) AOD at 550 nm between January 2013 and December 2017 at times relevant to Aqua's overpass time (13:30 LT). Aerosol index was calculated as the product of AOD and the Ångström parameter. MERRA-2 also provides surface mass concentrations of aerosol species including sea salt, dust, black carbon, organic carbon, and sulfate, which were used as a measure of aerosol levels in the planetary boundary layer (PBL).

MERRA-2 data were also used for environmental variables including both thermodynamic (e.g., temperature and relative humidity) and dynamic parameters (e.g., sea level pressure (SLP) and geopotential heights) (Gelaro et al., 2017) listed in Table 1. Bilinear interpolation was applied to transfer all MERRA-2 variables (Table 1) from their original 0.5× 0.625 spatial resolution to the equivalent 1× 1 grid in CERES-MODIS Level 3 data.

2.2.3 Precipitation data

Daily precipitation data were obtained from the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks–Climate Data Record (PERSIANN-CDR) data product (Ashouri et al., 2015; Nguyen et al., 2018). Bilinear interpolation was applied to convert the PERSIANN-CDR data from their native spatial resolution (i.e., 0.25× 0.25) to equivalent 1× 1 grids in CERES-MODIS Level 3 data. It is important to note that we use daily averaged PERSIANN-CDR precipitation and, therefore, there is some temporal mismatch with the daily Nd value from MODIS-Aqua that comes at one time of the day. This can contribute to some level of uncertainty for the discussions based on analyses involving relationships between precipitation and Nd.

2.2.4 Surface-based CCN data

Cloud condensation nuclei (CCN) data were obtained from the U.S. Department of Energy's Two-Column Aerosol Project (TCAP) (Berg et al., 2016) to examine the seasonal variations in CCN number concentration at a representative site by Cape Cod, Massachusetts (41.67 N, 70.30 W), over the US East Coast. TCAP was a campaign conducted between June 2012 and June 2013 to investigate aerosol optical and physicochemical properties and interactions between aerosols and clouds (Berg et al., 2016; Liu and Li, 2019). CCN data were available between July 2012 and May 2013 at multiple supersaturations with some gaps in the data collection (i.e., November–December); for simplicity, we focused on CCN data measured at a single supersaturation of 1 % owing to relatively better data coverage compared to lower supersaturations. We note that this higher supersaturation is not necessarily representative of that relevant to the clouds of interest but is still insightful for understanding the seasonal cycle of CCN concentration. The qualitative seasonal cycle of CCN concentration at 1 % matches those at lower supersaturations (e.g., 0.15 %–0.8 %).

2.2.5 Airborne in situ data

We used airborne in situ data collected during the fifth research flight (RF05) of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) campaign. One flight is used both for simplicity and because it embodied conditions relevant to the discussion of other results. The mission concept involves joint flights between the NASA Langley UC-12 King Air and HU-25 Falcon such that the former flies around 8–10 km, and the latter flies in the boundary layer to simultaneously collect data on aerosol, cloud, gas, and meteorological parameters in the same column (Sorooshian et al., 2019). The Falcon flew in a systematic way to collect data at different vertical regions relative to cloud, including the following of relevance to this study: BCB – below cloud base; ACB – above cloud base; BCT – below cloud top; and Min. Alt – minimum altitude the plane flies at ( 150 m).

This study makes use of the HU-25 Falcon data from the following instruments: fast cloud droplet probe (FCDP; Dp 3–50 µm) (SPEC Inc.) aerosol and cloud droplet size distributions for quantification of cloud liquid water content (LWC), Nd, and aerosol number concentrations with Dp exceeding 3 µm in cloud-free air (termed FCDP-aerosol); two-dimensional stereo (2DS; Dp 28.5–1464.9 µm) (SPEC Inc.) probe for estimation of rain water content (RWC) by integrating raindrop (Dp 39.9 µm) size distributions; cloud condensation nuclei (CCN; DMT) counter for CCN number concentrations; laser aerosol spectrometer (LAS; TSI model 3340) and condensation particle counter (CPC; TSI model 3772) for aerosol number concentrations with Dp between 0.1–1 µm and above 10 nm, respectively; high-resolution time-of-flight aerosol mass spectrometer (AMS; Aerodyne) for submicrometer non-refractory aerosol composition (DeCarlo et al., 2008), operated in 1 Hz Fast-MS mode and averaged to 25 s time resolution; and turbulent air-motion measurement system (TAMMS) for winds and temperature (Thornhill et al., 2003).

CCN, LAS, CPC, and AMS data were collected downstream of an isokinetic double diffuser inlet (BMI, Inc.), whereas the AMS and LAS also sampled downstream of a counterflow virtual impactor (CVI) inlet (BMI, Inc.) when in cloud (Shingler et al., 2012). However, a filter was applied to remove LAS data when the CVI inlet was used. Measurements from the CCN counter, LAS, CPC, and FCDP-aerosol are only shown in cloud-free and rain-free conditions, distinguished by LWC < 0.05 g m−3 and RWC < 0.05 g m−3, respectively, and also excluding data collected 20 s before and after evidence of rain or cloud. Estimation of supermicrometer particles from FCDP measurements was performed after conducting the following additional screening steps to minimize cloud droplet artifacts: (i) only samples with RH < 98 % were included; (ii) data collected during ACB and BCT legs were excluded. CCN, LAS, CPC, and AMS measurements are reported at standard temperature and pressure (i.e., 273 K and 101.325 kPa) while FCDP and 2DS measurements correspond to ambient conditions.

2.3 Regression analyses

Regression modeling was conducted to investigate relationships between environmental variables and Nd. The gradient-boosted regression trees (GBRT) model, classified as a machine learning (ML) model, is used, consisting of several weak learners (i.e., regression trees with a fixed size) that are designed and subsequently trained to improve prediction accuracy by fitting the model's trees on residuals rather than response values (Hastie et al., 2009). Desirable characteristics of the GBRT model include both its capacity to capture non-linear relationships and being less vulnerable to overfitting (Persson et al., 2017; Fuchs et al., 2018; Dadashazar et al., 2020). Two separate GBRT models were trained using daily CERES-MODIS Nd data (1× 1) in winter (DJF) and summer (JJA) to reveal potential variables impacting Nd. Winter and summer are chosen as they exhibit the highest and lowest Nd concentrations, respectively, among all seasons over the WNAO.

Table 2List of input parameters used as predictor variables in the GBRT and linear models. Variables are grouped into three general categories.

 Skin potential temperature.

Download Print Version | Download XLSX

Many variables were picked as input parameters (Table 2) for the GBRT model, categorized as being aerosol, dynamic/thermodynamic, or cloud variables. Aerosol parameters included MERRA-2 surface mass concentrations for sulfate, sea salt, dust, and organic carbon. Black carbon concentration was removed from input parameters because of its high correlation (R2=0.6) with organic carbon. The following is the list of thermodynamic/dynamic input parameters derived from MERRA-2: vertical pressure velocity at 800 hPa (ω800), planetary boundary layer height (PBLH), cold-air outbreak (CAO) index, wind speed and wind direction at 2 m (wind2 m and wind-dir2 m), and relative humidity (RH) in the PBL and free troposphere represented by RH950 and RH800, respectively. CAO index is defined as the difference between skin potential temperature (θskt) and air potential temperature at 850 hPa (θ850) (Papritz et al., 2015). Updraft velocity plays a crucial role in the activation of aerosol into cloud droplets in warm clouds (Feingold, 2003; Reutter et al., 2009). Since the direct representation of updraft speed is not available from reanalysis data, near-surface wind speed (i.e., wind2 m) is used as a representative proxy parameter as an input parameter to the regression models. CERES-MODIS cloud parameters include liquid cloud fraction and cloud-top height for low-level clouds. In addition, PERSIANN-CDR daily precipitation (Rain) was included as a relevant cloud parameter.

Data were split into two sets: training/validation (70 %) and testing (30 %). Five-fold cross-validation was implemented to train the GBRT model using the training/validation data. Furthermore, both performance and generalizability of the trained models were tested via the aid of the test set, which was not used in the training process. Hyperparameters of the GBRT models were optimized through a combination of both random and grid search methods. Table S1 in the Supplement shows the list of important hyperparameters of the GBRT model and associated ranges tested via random and grid search methods. The optimized model hyperparameters can also be found in Table S1. The GBRT models were performed using the scikit-learn module designed in Python (Pedregosa et al., 2011).

The regression analyses were not performed solely to construct and provide a highly accurate model useful for prediction, but rather to disclose and examine the possible effects of the relevant input variables on Nd considering all the shortcomings of such analyses. For instance, there is some level of interdependency between input variables. To reduce unwanted consequences of correlated features, the interpretation of the results was done with the aid of accumulated local effect (ALE) plots, which are specifically designed to be unbiased to the correlated input variables (Apley and Zhu, 2020). ALE plots illustrate the influence of input variables on the response parameter in ML models. The ALE value for a particular variable s at a specific value of xs (i.e., fs,ALE (xs)) can be calculated as follows:

(2) f s , ALE x s = z 0 , 1 x s x c f s z s , x c P x c | z s d x c d z s - constant ,

where fs(zsxc) is the gradient of model's response with respect to variable s (i.e., local effect) and P(xc|zs) is the conditional distribution of xc, where c denotes the other input variables rather than s, and xc is the associated point in the variable space of c. z0,1 is chosen arbitrarily below the smallest observation of feature s (Apley and Zhu, 2020). The steps in Eq. (2) can be summarized as follows (Molnar, 2019; Apley and Zhu, 2020): (i) the average change in the model's prediction is calculated using the conditional distribution of features; (ii) the average change will then be accumulated by integrating it over feature s; and (iii) a constant will be subtracted to vertically center (i.e., the average of ALE becomes zero) the ALE plot. The aforementioned steps, although seemingly complex, assure the avoidance of undesired extrapolation (especially an issue for correlated variables) occurring in alternative approaches such as partial dependence (PD) plots. The value of fs,ALE(xs) can be viewed as the difference between the model's response at xs and the average prediction. We used the source code available in (Jumelle et al., 2021) for the calculation of ALE plots.

3 Results and discussion

3.1 Aircraft case study of Nd gradient

ACTIVATE Research Flight 5 (RF05) on 22 February 2020 demonstrates the wide range in Nd offshore in the PBL (1.6 km) over the WNAO (Fig. 1). On this day, the ACTIVATE study region was dominated by a surface high-pressure system centered over the southeastern US, with a significant ridge axis extending from the main high to the east-northeast off the Virginia–North Carolina coast and into the WNAO. Aloft, the flight region was located in northwesterly flow behind a trough offshore. This setup led to subsidence in the region and generally clear skies, except where scattered to broken marine boundary layer clouds formed along and east of the Gulf Stream. The 2 d NOAA HYSPLIT (Stein et al., 2015; Rolph et al., 2017) back trajectories using the “model vertical velocity” method and “REANALYSIS” meteorology data indicate air in the flight region (between 0–3 km) had wrapped around the surface high from the north and left the New England coast 12–24 h beforehand (with a descending profile). Along the flight segment shown, winds were approximately 6 m s−1, out of the north-northwest during the initial descent, Min. Alt. 1, and BCB1 legs and primarily from the northeast for the other sections of the flight. Sea surface temperatures were 6–9 C near the coast during the descent and Min. Alt. 1 leg (readers are referred to Fig. 1's caption for the definition of different legs); 21–25 C over the Gulf Stream during the BCB1, ACB1, and BCB2 legs; and 17–20 C for the remainder of the flight segment shown. The majority of the segment was in or below the boundary layer clouds, with cloud base around 900–1100 m and cloud top around 1750 m. Note that the initial BCB1 leg was much lower at around 460 m, likely reflecting a shallower marine boundary layer and cloud base near the much colder waters close to the coast. Static air temperature ranged between 0–10 C, except for the BCT1 leg where temperatures were around 2.3 C.

Figure 1Time series of selected parameters measured by the HU-25 Falcon aircraft during a selected segment of RF05 on 22 February 2020. (a) Overlayed flight track on GOES 16 visible imagery obtained at 14:55:04 UTC. (b) Altitude, cloud liquid water content (LWC), and Nd, with the latter two obtained from the FCDP. (c) Rain water content (RWC) measured by 2DS probe, AMS speciated mass concentration in cloud/rain-free air, and AMS mass fractions for droplet residual particles in cloud as measured downstream of a CVI inlet. (d) Number concentrations for CCN at 0.43 % supersaturation and particles for three diameter ranges: above 10 nm (CPC), 100–1000 nm (LAS), and above 3 µm (FCDP). Shaded gray areas in panels (b)(d) highlight cloudy periods identified as having LWC  0.05 g m−3. Locations of the cloudy regions are pointed to with red arrows in the satellite imagery. Level legs are defined as follows: BCB: below cloud base; ACB: above cloud base; Min. Alt.: minimum altitude the plane flies at ( 150 m); ACT: above cloud top; BCT: below cloud top.

Nd values from the FCDP ranged from a maximum value of 1298 cm−3 closer to the coast during the ACB1 leg (35.00 N, 74.55 W) to a minimum of 19 cm−3 farther away in the BCT1 leg (34.32 N, 72.73 W). The minimum Nd value in the ACB3 leg was 85 cm−3 (34.11 N, 72.80 W), which is a fairer comparison to the ACB1 leg compared to the BCT1 leg in terms of being closer to cloud base. The mean Nd values (cm−3) in the cloudy portions of the ACB1, BCT1, and ACB3 legs were as follows: 849, 77, and 143.

Based on the nearest BCB legs adjacent to the maximum and minimum Nd values (BCB1 = 35.31 N, 74.95 W; BCB3 = 34.41 N, 72.70 W), there was a significant offshore gradient in LAS submicrometer particle number concentration and AMS non-refractory aerosol mass, ranging from as high as 424 cm−3 and 5.60 µg m−3 (during BCB1) to as low as 21 cm−3 and 0.27 µg m−3 (during BCB3). The mean values of submicrometer particle number concentration and AMS non-refractory aerosol for the two BCB legs were as follows: 277 cm−3/3.64 µg m−3 (BCB1) and 48 cm−3/0.42 µg m−3 (BCB3). The higher Nd value (1298 cm−3) relative to LAS aerosol concentration (424 cm−3) at the near-shore point is suggestive of aerosol smaller than 0.1 µm activating into drops. This is supported by the fact that both CCN (supersaturation = 0.43 %) and CPC number concentrations with Dp> 10 nm exhibited mean values of 980 and 1723 cm−3 in the BCB1 leg, respectively, dropping to 98 and 260 cm−3 in the BCB3 leg. For the duration of the flight portion shown in Fig. 1, supermicrometer concentrations varied over 2 orders of magnitude (0.002–0.51 cm−3) and expectedly did not exhibit a pronounced offshore gradient as it is naturally emitted from the ocean.

Closer to shore during the Min. Alt. 1 leg, nitrate was the dominant aerosol species ( 70 % mass fraction). Farther offshore during both the BCB1 leg and cloud-free portion of the ACB1 leg, organics were the dominant constituent ( 46 % mass fraction), whereas farther during the BCB3 leg, the mean mass fraction of sulfate was the highest (75 %). Droplet residual particle data show a greater contribution of organics farther offshore, increasing from 46 % to 75 % between the ACB1 and ACB3 legs, respectively. These composition results, albeit limited to the non-refractory portion of submicrometer aerosol particles, reveal significant changes with distance offshore indicative of varying chemical properties of particles activating into droplets.

The cloudy portions of ACB1 are characterized as having little or no rain with a maximum RWC value of 0.02 g m−3 and mean value of 0.003 g m−3. There is a notable RWC peak at the beginning of the Min. Alt. 2 leg, reaching as high as 1.81 g m−3 associated with clouds aloft. The precipitation occurrence was also evident in a subsequent BCT1 leg where RWC reached as high as 0.18 g m−3. GOES satellite imagery of the study region (Fig. 1) also reflects the effect of precipitation on cloud morphology where clouds farther offshore resemble open-cell structures. Associated scavenging of particles through the washout process is presumed to contribute to the decline in aerosol concentrations with distance offshore.

Figure 1 shows changes in aerosol characteristics coincident with the large gradient in Nd. While ACTIVATE airborne data collection is ongoing to build flight statistics over multiple years, the wide changes in microphysical properties in RF05 motivate looking at other datasets with broader spatiotemporal coverage to learn about potential seasonally dependent drivers of Nd, including meteorological parameters that vary throughout the year. Furthermore, other datasets can provide insight into the source(s) of seasonal discrepancy between columnar aerosol remote sensing parameters and Nd.

3.2 Seasonal cycles of Nd and AOD

Figure 2 illustrates the seasonal differences in MERRA-2 AOD and CERES-MODIS Nd over the WNAO that partly motivate this study. Seasonal mean values (± standard deviation) of AOD/Nd (cm−3) were as follows for the entire WNAO: DJF = 0.11 ± 0.03/64.1 ± 18.0; MAM = 0.16 ± 0.03/60.4 ± 13.1; JJA = 0.15 ± 0.03/49.1 ± 10.1; and SON = 0.11 ± 0.03/50.3 ± 13.9. In contrast to AOD, Nd values and low-cloud fraction (Fig. 2c) were highest in DJF and lowest in JJA. DJF showed notably high Nd near the coast, qualitatively consistent with the airborne data. The seasons with the greatest AOD values, accompanied by the most pronounced spatial gradient offshore, were JJA and MAM. The offshore gradient owes to continental pollution outflow (Corral et al., 2021, and references therein). In contrast, DJF and SON exhibited lower AOD values with a distinct area of higher AOD values offshore between  35–40 N accounted for by sea salt. MERRA-2 speciated AOD data (Fig. 3) indicate that sea salt and sulfate dominate total AOD regardless of season and that sulfate, organic carbon, and black carbon most closely follow the offshore gradient pattern owing to continental sources. Dust and sea salt have different spatial distributions with the former derived from sources such as North Africa leading to enhanced AODs < 30 N in particular in JJA and sea salt being enhanced offshore in particular in JJA.

Figure 2Seasonal spatial maps for (a) MERRA-2 aerosol optical depth (AOD), (b) MERRA-2 aerosol index (AI), and (c) cloud drop number concentration (Nd) over the western North Atlantic Ocean (WNAO). Contours in panel (c) represent low-level (cloud-top pressure > 700 hPa) liquid cloud fraction (CFlow-liq.). Cloud data are based on daily Level 3 data from CERES-MODIS. The maps are based on data between January 2013 and December 2017. The boxes in the top left panel represent sub-domains examined in more detail throughout the study, with the blue star denoting Bermuda.

Figure 3Seasonal maps of MERRA-2 speciated AOD based on data between January 2013 and December 2017. The boxes in the top left panel represent sub-domains examined in more detail throughout the study, with the blue star denoting Bermuda.

Table 3Average drop number concentration (Nd), MERRA-2 AOD, and vertically resolved AOD characteristics from CALIOP for each season over the sub-domains shown in Fig. 2. Total CALIOP AOD is shown outside parentheses, and numbers inside are the percent AOD fraction in the planetary boundary layer followed by in the free troposphere. Also shown are PBLHs (shown in Fig. 4) and the relative humidity in the PBL and FT.

Download Print Version | Download XLSX

Table 3 probes deeper into individual WNAO sub-domains to compare seasonal AOD and Nd values. For the six sub-domains in Fig. 2, MERRA-2 AOD peaks in MAM and JJA, while Nd peaks in DJF. The Bermuda sub-domain was unique in that mean Nd was slightly higher in MAM (53 cm−3) compared to DJF (48 cm−3). We attribute the slightly different seasonal cycle over Bermuda to its remote nature, leading to differences in meteorology and aerosol sources between seasons.

One factor that could bias AOD towards higher values with disproportionately less impact on Nd is aerosol hygroscopic growth in humid conditions. Table 3 summarizes mean MERRA-2 RH values in the PBL and free troposphere (FT). Results show that while RH is highest in JJA (except for FT of DJF in sub-domain N), differences between seasons were not very large. The maximum difference among the four seasons when considering mean RH in the PBL and FT for all sub-domains ranged between 3 %–9 % and 7 %–25 %, respectively. Consequently, humidity effects on remotely sensed aerosol parameters are less likely to be the sole explanation of the dissimilar seasonal cycle of Nd and AOD, but can plausibly contribute to some extent.

One factor that could drive the seasonal variation in Nd is the unwanted effects of retrieval errors in the estimation of Nd at low-cloud-coverage conditions. Uncertainty associated with the estimation of Nd from MODIS observation increases as cloud fraction decreases (Grosvenor et al., 2018). This is mainly because of the overestimation of droplet effective radius (re) in the retrieval algorithm due to the interference of cloud-free pixels and also high spatial inhomogeneity in low-cloud-coverage conditions that violates horizontal homogeneity assumptions in the retrieval of re and τ from radiative transfer modeling (Zhang et al., 2012, 2016). To test whether retrieval errors in Nd are the main driver of seasonal trends, Fig. S1 shows the seasonal cycle of Nd at various low-level liquid cloud fractions. The results show that as cloud fraction increases the average Nd increases, regardless of season. Perhaps the more important result is that the seasonal trend in spatial maps of Nd remains similar regardless of cloud fraction. This finding is important as it confirms that the seasonal cycle in Nd cannot be solely explained by the uncertainties associated with the retrieval of Nd at low cloud fraction.

3.3 Contrasting AOD and aerosol index

While previous studies have pointed to the limitations of AOD as an aerosol proxy (e.g., Stier, 2016; Gryspeerdt et al., 2017; Painemal et al., 2020), the Nd–AOD anticorrelation at seasonal scale over the WNAO is at odds with findings for other regions supporting the relationship between these two parameters (Nakajima et al., 2001; Sekiguchi et al., 2003; Quaas et al., 2006, 2008; Grandey and Stier, 2010; Penner et al., 2011; Gryspeerdt et al., 2016) and also that between sulfate and Nd (Boucher and Lohmann, 1995; Lowenthal et al., 2004; Storelvmo et al., 2009; McCoy et al., 2017, 2018; MacDonald et al., 2020). Values of Nd are influenced by the number concentration of available CCN, which is determined by aerosol properties (size distribution and composition) and supersaturation level. AOD is an imperfect CCN proxy variable because it does not provide information about composition and size distribution and is sensitive to relative humidity. Aerosol index (AI) is more closely related to CCN as it partially accounts for the size distribution of aerosols (Deuze et al., 2001; Nakajima et al., 2001; Breon et al., 2002; Hasekamp et al., 2019). The sensitivity of AI to size is evident in spatial maps for each season showing more of an offshore gradient (like sulfate AOD in Fig. 3) in each season and lacking both the offshore peak in sea salt between  35–40 N and the maximum AOD for dust south of 30 N in JJA. However, when comparing absolute values between the four seasons in Fig. 2b, AI exhibits a similar seasonal cycle to AOD, thereby indicating that size distribution alone cannot explain diverging seasonal cycles for Nd and AOD. We next compare Nd to aerosol data in the PBL where CCN more relevant to droplet activation are confined. Size distribution effects in the PBL can instead be more of a factor, especially as sea salt is abundant.

3.4 Aerosol size distribution and vertical aerosol distribution

Vertical profiles of aerosol extinction coefficient estimated from CALIOP nighttime observations are shown in Fig. 4 for the six sub-domains. Also shown are the seasonally representative planetary boundary layer heights (PBLHs) from MERRA-2, with numerical values of both PBLH and fractional AOD contributions to the PBL and FT in Table 3. Although here we used nighttime observations from CALIOP because of having higher signal-to-noise ratio than daytime observations, we expect the general seasonal trends discussed here to remain the same regardless of the observation time. The CALIOP results indicate that aerosol extinction more closely follows the Nd seasonal cycle with the highest (lowest) values in the PBL during DJF (JJA). However, aerosol extinction coefficient is sensitive to aerosol size distribution, and a plausible scenario is that DJF extinction in the PBL is primarily contributed by coarse sea salt particles, which are especially hygroscopic but do not contribute significantly to number concentration as demonstrated clearly by airborne observations (i.e., FCDP>3µm time series shown in Fig. 1d). This is supported in part by how DJF is marked by the highest fractional AOD contribution from the PBL (59 %–72 %) where sea salt is concentrated. In contrast, JJA has the lowest fractional AOD contribution from the PBL (11.3 %–52.6 %). It is also possible that the higher fractional AOD contribution from the PBL in winter is partly owed to aerosol particles being more strongly confined to the PBL compared to the summer. Sub-domains C-N and N exhibit the greatest changes in AOD fraction in the PBL between seasons with a maximum in DJF (59 %–61 %) and a minimum in JJA (11 %–19 %), suggesting they are relatively more sensitive to the aerosol vertical distribution in leading to contrasting AOD and Nd seasonal cycles. Bermuda stands out as having the highest AOD fractional contributions in the PBL in DJF (72 %) and SON (69 %) and among the highest seasonal total AODs in those two seasons (0.14 in DJF and 0.10 in SON) assisted in large part by sea salt (Fig. 3) (Aldhaif et al., 2021), coincident with high seasonal wind speeds (Corral et al., 2021).

Figure 4Vertical profiles of CALIPSO aerosol extinction for different seasons in (a)(f) six different sub-domains of the WNAO. Average seasonal planetary boundary layer heights (PBLH) from MERRA-2 are denoted with dashed lines.

Figure 5Monthly statistics of CCN concentration (1 % supersaturation) measured at Cape Cod between July 2012 and May 2013. Red lines represent the median, whiskers are the monthly range, and the top and bottom of the boxes represent the 75th and 25th percentiles, respectively. The notches in the box plots demonstrate whether medians are different from each other with 95 % confidence. Boxes with notches that do not overlap with each other have different medians with 95 % confidence.

To explore aerosol number concentration characteristics in the PBL in different seasons, we next discuss results from an opportune dataset over the US East Coast (Cape Cod, MA) providing an annual profile of CCN concentration at 1 % supersaturation (Fig. 5). Cape Cod is a coastal location representative of the outflow, providing an important fraction of the CCN impacting offshore low-level clouds. As the supersaturation examined is relatively high (1 %), the measured CCN include smaller particles representing high number concentrations that would not appreciably contribute to the high aerosol extinctions from CALIOP in the PBL in direct contrast to sea salt (i.e., high extinction due to fewer but larger particles). Seasonal mean CCN values do not follow the seasonal cycle of Nd nor CALIOP extinction in the PBL, with values being as follows: DJF = 1436 cm−3; MAM = 1533 cm−3; JJA = 1895 cm−3; and SON = 1326 cm−3. These results suggest the following: (i) size distribution effects are significant in the PBL when comparing extinction to number concentration, and (ii) aerosol vertical distribution behavior cannot alone explain the divergent seasonal cycles of Nd and aerosol parameters (e.g., AOD, AI, surface number concentrations).

Figure 6Seasonal maps of MERRA-2 speciated aerosol concentrations at the surface based on data between January 2013 and December 2017. The boxes in the top left panel represent sub-domains examined in more detail throughout the study, with the blue star denoting Bermuda.

We next compare MERRA-2 speciated aerosol concentrations at the surface (Fig. 6) to those of speciated AOD (Fig. 3). Surface mass concentrations have the limitation of being biased by larger particles (similar to extinction). The seasonal cycle of mean values for speciated AOD and surface concentration for individual sub-domains generally agree with the exception that there was disagreement for sulfate in each sub-domain (see seasonal mean values in Table S2). Sulfate exhibited higher AODs in JJA but with surface concentrations usually being highest in DJF or MAM; although differences in seasonal mean mass concentrations were relatively small (< 1 µg m−3). A plausible explanation is enhanced secondary production of sulfate via oxidation of SO2 or DMS convectively lifted to the free troposphere in JJA. An important result confirmed by the surface mass concentrations is that sea salt is an order of magnitude higher than the other species, supporting the previous speculation that sea salt dominates the aerosol extinction in the PBL from CALIOP.

3.5 Aerosol–cloud interactions

Studies of China's east coast have shown that the aerosol indirect effect is especially strong in wintertime, whereby pollution outflow leads to high Nd and suppressed precipitation (Berg et al., 2008; Bennartz et al., 2011). It is hypothesized that a similar effect is taking place off of North America's east coast, which could in part explain enhanced Nd without a significant jump in aerosol parameter (e.g., AOD, AI) values necessarily. Grosvenor et al. (2018) suggested that high cloud fractions in wintertime off these east coasts relative to other seasons are coincident with strong temperature inversions usually associated with cold-air outbreaks that serve to concentrate and confine surface layer aerosols. We examine the relative seasonal strength of the aerosol indirect effect via spatial maps of the following metric commonly used in aerosol–cloud interaction (ACI) studies:

(3) ACI = d ln ( N d ) / d ln ( α ) ,

where α represents an aerosol proxy parameter that is represented here as AI, AOD, the speciated sulfate AOD (SulfateAOD), and sulfate surface mass concentration (Sulfatesf-mass). The expected range by common convention is 0–1, with higher values suggestive of greater enhancement in Nd for the same increase in the aerosol proxy parameter.

Table 4Estimated values of ACI calculated four ways (dlog(Nd)/dlog(AOD); dlog(Nd)/dlog(AI); dlog(Nd)/dlog(SulfateAOD); dlog(Nd)/dlog(Sulfatesf-mass)) for the sub-domains shown in Fig. 2. The ACI values were obtained from log–log regression on average daily values of Nd and each of the aerosol proxy variables including only the pixels with CFlow-liq. greater than 0.1. Numbers in parentheses, in order, are R2 and the number of points used for linear regression. Statistically insignificant ACI values with p values greater than 0.05 are marked by bold font​​​​​​​.

Download Print Version | Download XLSX

Figure 7Seasonal maps of the aerosol–cloud interaction (ACI) parameters over the WNAO using daily Nd and four different aerosol proxy parameters (AI, AOD, SulfateAOD, Sulfatesf-mass) from CERES-MODIS and MERRA-2, respectively. ACI statistics associated with the six sub-domains shown are summarized in Table 4.

Table 4 shows that DJF always exhibits the highest ACI values regardless of the aerosol proxy used, consistent with a stronger aerosol indirect effect in DJF over East Asia. The mean ACI values in DJF using AI, AOD, SulfateAOD, and Sulfatesf-mass ranged from 0.25 to 0.55, 0.28–0.59, 0.25–0.53, and 0.22–0.47, respectively, depending on the sub-domain. Spatial maps of ACI (Fig. 7) do not point to significant geographic features. Coefficients of determination (R2) for the linear regression between ln(Nd) and ln(α) when computing seasonal ACI values were generally low ( 0.30), with spatial maps of R2 and data point numbers in Fig. S2. Poor correlations are suggestive of the non-linear nature of aerosol–cloud interactions (e.g., Gryspeerdt et al., 2017) and the influence of other likely factors such as dynamical processes and turbulence, data spatial resolution and dataset size, cloud adiabaticity, wet scavenging effects, and aerosol size distribution (McComiskey et al., 2009). The results of this section suggest though that aerosol indirect effects could be strongest in DJF, meaning that Nd values increase more for the same increase in aerosol. Factors that can contribute to higher ACI values in winter than summer include seasonal differences in the following: (i) dynamical processes and turbulent structures of the marine boundary layer; (ii) aerosol size distributions and consequently varying particle number concentrations for a fixed mass concentration; and (iii) hygroscopicity of particles, especially as a result of changes in the composition of the carbonaceous aerosol fraction. Regarding dynamical processes and the effects of turbulence, Fig. 2 in Painemal et al. (2021) shows that heat fluxes (i.e., latent and sensible fluxes) are strongest (lowest) in the winter (summer) over the WNAO. The greater heat fluxes in DJF can contribute to more turbulent and coupled marine boundary layer conditions in winter than summer, presumably resulting in more efficient transport and activation of aerosol in the marine boundary layer, leading to higher ACI values. Forthcoming work will probe this issue in greater detail.

4 Discussion of potential influential factors

We probe deeper into factors related to the Nd seasonal cycle by using (Sect. 4.1) composite analyses based on high- and low-Nd days and (Sect. 4.2) advanced regression techniques tackling non-linear relationships. We focus the analyses on one sub-domain (C-N) for both simplicity and intriguing characteristics: (i) among the highest anthropogenic AOD values over the WNAO; (ii) significant seasonal changes in fractional AOD contribution to the PBL; (iii) close to the Cape Cod site where CCN data were shown; and (iv) the aerosol indirect effect (Table 4) strongest (weakest) in DJF (JJA).

4.1 Composite analysis

Discussion first addresses the behavior of different environmental parameters on days with the highest and lowest Nd values. Seasonal histograms of averaged daily Nd were generated for sub-domain C-N. The histograms are based on the natural logarithm of Nd to better resemble a normal distribution. We assign values as being low in each season if they are less than 1 standard deviation below the seasonal value; conversely, high values are those exceeding 1 standard deviation above the seasonal mean. Cut-off Nd values (cm−3) are as follows (low/high): 33/153 (DJF), 29/118 (MAM), 38/100 (JJA), and 31/115 (SON). Next, composite maps for these groups were created (Figs. 8–12) for sea level pressure, near-surface wind, low-level cloud fraction, cold-air outbreak index, and AOD. The figures contrast the low- and high-Nd maps with those showing mean seasonal values to investigate potential factors that contribute to seasonal Nd variability. Interested readers are referred to Figs. S3–S20 where similar composite map results are shown for Nd itself and other parameters including those in Table 2.

Figure 8Seasonal climatology of sea level pressure (SLP) (middle column) and anomalies from seasonal averages for low-Nd days (left column) and high-Nd days (right column). In the left and right columns, red and blue contours are associated with positive and negative anomalies from the climatology, respectively. The green box represents sub-domain C-N for which the analysis was conducted.

Figure 9Seasonal climatology of near-surface (2 m above ground) wind speed (middle column) and mean values for low-Nd days (left column) and high-Nd days (right column). The reference wind vector is shown in the top left panel. The red box represents sub-domain C-N for which the analysis was conducted.

Figure 10Seasonal averages of low-level liquid cloud fraction (middle column) and associated anomalies on low-Nd days (left column) and high-Nd days (right column). The red box represents sub-domain C-N for which the analysis was conducted.

Figure 11Seasonal averages of cold-air outbreak (CAO) index (middle column) and associated anomalies on low-Nd days (left column) and high-Nd days (right column). The red box represents sub-domain C-N for which the analysis was conducted.

Figure 12Seasonal averages of MERRA-2 AOD (middle column) and associated anomalies on low-Nd days (left column) and high-Nd days (right column). The red box represents sub-domain C-N for which the analysis was conducted.

The resulting composite maps indicate high-Nd days are characterized by (i) reduced SLP; (ii) more northerly-northwesterly flow for all seasons (except JJA) and especially stronger winds in DJF and SON; (iii) higher low-level liquid cloud fraction, especially in DJF; (iv) higher CAO index in the seasons when CAO events occur more frequently (DJF, SON, MAM); and (v) enhanced AOD. Low-Nd days generally exhibited opposite conditions when compared to seasonal mean values: (i) enhanced SLP; (ii) wind ranging from southerly to westerly without any significant wind speed enhancement; (iii) reduced low-level liquid cloud fraction, especially in DJF; (iv) lower CAO index in DJF, SON, and MAM; and (v) reduced AOD in DJF and MAM, enhanced AOD in JJA, and limited change in SON. Noteworthy results from Figs. S3–S20 included the enhancement/reduction of PBLH on high-/low-Nd days (least pronounced in JJA), higher/lower RH at 950 and 800 hPa on high-/low-Nd days, and higher/lower sulfate AOD and surface concentrations on high-/low-Nd days for DJF and MAM. Furthermore, there was a general reduction in rain on low-Nd days for most seasons except SON, with rain enhancement on high-Nd days except for DJF (Fig. S6); this was unexpected as wet removal was hypothesized to be a reason for reduced Nd for at least the low-Nd days. This may be attributed to the rain product being for surface precipitation (and thus not capturing all drizzle) and for all cloud types, including more heavily precipitating clouds deeper and higher than the low-level clouds examined for Nd. Another factor potentially contributing to the observed counterintuitive trends is the temporal offset between Nd estimations from MODIS-Aqua and precipitation data from PERSIANN-CDR.

The mean seasonal climatological values and anomalies suggest that high-Nd cases are marked by continental outflow, high cloud fractions, high PBLH, and low SLP, all of which occur most commonly in DJF and are associated with cold-air outbreaks. These events are marked by cold air over the warm ocean leading to strong surface heat fluxes, boundary layer deepening, weakened inversion strength, and high and deep clouds (Brummer, 1996; Kolstad et al., 2009; Fletcher et al., 2016; Abel et al., 2017; Naud et al., 2018). Coincident with these features is the Icelandic Low, which is a significant climatological feature of the North Atlantic whereby subpolar low pressure builds in extratropic areas beginning in the fall with westerly winds in the boundary layer that shift more to northerly in the winter (Sorooshian et al., 2020; Painemal et al., 2021). This low-pressure system seems to be stronger on high-Nd days, resulting in more continental outflow and high number concentrations of CCN; the greater CAO index values near the coast promote high cloud coverage, affording more opportunity for cloud processing of particles to ultimately enhance droplet activation. While there can be considerable enhancement in Nd as cold-air outbreak air masses evolve over warmer waters, precipitation scavenging farther downwind will be an efficient method of boundary layer aerosol (and Nd) removal (Abel et al., 2017; Lloyd et al., 2018), which contributes at least in part to the sharp Nd gradients offshore demonstrated in Fig. 1.

4.2 Multivariate regression analysis

Modeling analysis focuses on the two seasons (DJF and JJA) with the extremes in terms of seasonal mean values for Nd and aerosol parameters. Added motivation for examining those two seasons stems from spatial maps of R2 based on ACI analysis (Fig. S2). Using the surface sulfate concentration as the aerosol proxy generally yielded higher R2 values in three seasons (DJF = 0.13, MAM = 0.05, SON = 0.08) except JJA (0.02) for which the choice did not matter owing to low R2 ( 0.03) values for all four aerosol proxy variables tested. Although the R2 values are all generally low, DJF and JJA are the seasons when surface sulfate levels are the most and least capable of explaining Nd, with R2 among the four proxy variables exhibiting the widest (DJF values: 0.07–0.13) and narrowest range (JJA: 0.01–0.03) of values. We address here how much improvement is gained in modeling Nd by advancing from linear regressions based on one input variable to (i) adding more input variables and (ii) moving to a more sophisticated model (GBRT) that captures non-linear relationships.

Table 5Performance of different models in predicting Nd assessed based on average R2 scores on both validation and test sets. The models were fitted separately for DJF and JJA seasons. Table 2 has the complete list of variables used in the GBRT model.

Download Print Version | Download XLSX

We show in Table 5 the performance of two linear models based on a single linear regression (with sulfate mass concentration) and a multi-regression that uses 14 input variables listed in Table 2. In addition, Table 5 also lists the performance of the GBRT model that ingests 14 input variables, similar to the linear multi-regression model. The average R2 scores of the test set for predicting Nd based on a linear regression using only sulfate surface mass concentration were 0.17 and 0.09 in DJF and JJA, respectively. In contrast, R2 between the multi-regression linear model and the test dataset increased to 0.28 and 0.25 for DJF and JJA, respectively. This increase in predictive capability was helpful to reduce the gap between seasons by presumably accounting for factors more important in JJA aside from surface concentration of sulfate. The R2 scores increased even more to 0.47 and 0.43 for DJF and JJA, respectively, for the GBRT model. Therefore, accounting for non-linear relationships improved predictive capability in both seasons. It is important to note that the GBRT model was robust in terms of overfitting and especially generalizability as R2 values of the test and validation sets were similar for both seasons.

Figure 13Average permutation feature importance of input parameters for (a) DJF and (b) JJA based on GBRT models trained in each season. Feature importance values were calculated based on using the test set. Error bars exhibit the range of feature importance values stemming from the variability of the obtained models from the cross-validation resampling procedure.


We next discuss the importance ranking of different parameters from Table 2 in terms of influencing Nd for DJF and JJA (Fig. 13). Low-level liquid cloud fraction was the most important parameter in both seasons with some commonality in the next three parameters for both seasons. In DJF, sulfate surface mass concentrations were the second most important factor, followed by organic carbon surface concentrations and low-level liquid cloud-top effective height. As sulfate is secondarily produced via gas-to-particle conversion processes, this result is consistent with those from Fig. 1 showing the presumed strong impact of particles smaller than 100 nm in impacting Nd values close to shore. In JJA, the CAO index was the second most important, followed by organic carbon and sulfate surface concentrations. Also, our results throughout the study and supported by modeling are in agreement with Quinn et al. (2017) that sulfate particles contribute more to the CCN budget than sea salt particles. In DJF and JJA, the fifth most important factor was CAO index (second most important in JJA) and PBLH (11th most important in DJF), respectively.

Figures 14 and 15 show accumulated local effect (ALE) plots for the various parameters ranked in Fig. 13. In both seasons, but especially DJF, enhanced surface concentrations of sulfate and organic carbon coincide with higher Nd, whereas there was not any obvious positive association between Nd and either sea salt or dust (Fig. 14). Dust in JJA and sea salt in DJF, seasons of which each respective aerosol type is most predominant, exhibited negative relationships with Nd. Such a negative relationship is plausibly related to differences between ACI when calculated using AOD versus AI (Painemal et al., 2021); for instance, coarse sea salt can expedite collision–coalescence and thus reduce Nd, which has the effect of reducing ACI (Eq. 3) and even possibly yielding negative values (Table 4). Negative values of other ACI constructs coincident with poor R2 values have previously been attributed to potential effects of giant CCN (Terai et al., 2015; Dadashazar et al., 2017), but further research needs to examine this in more detail.

Figure 14Average local accumulated effect (ALE) profiles based on GBRT modeling for surface mass concentrations of the following parameters: (a) dust, (b) organic carbon, (c) sea salt, and (d) sulfate. Blue and red profiles represent ALEs of DJF and JJA, respectively. Shaded areas show the ALE ranges stemming from the variability of the obtained models from the cross-validation resampling procedure. Markers on the bottom and top x axes denote the values of the 5th, 25th, 50th, 75th, and 95th percentiles for each input variable.


Figure 15Same as Fig. 14 but for the following input parameters: (a) low-level liquid cloud fraction (CFlow-liq.), (b) cloud-top effective height of low-level liquid cloud (cloud-toplow-liq.), (c) cold-air outbreak (CAO) index, and (d) vertical pressure velocity at 800 hPa (ω800).


Figure 15 shows the similarity in the positive relationship between cloud fraction and Nd in both seasons. Only in DJF did cloud-top effective height exhibit a clear relationship with Nd (positive), likely linked to the common phenomenon of CAOs noted in Sect. 4.1 based on heightened CAO index values, deepening of the boundary layer, and weakened inversion strength. This is supported by enhanced Nd values coincident with negative values for ω800 (i.e., rising motion) and CAO index values above 0 in DJF without such relationships in JJA (Fig. 15). The six parameters in Fig. S21 (PBLH, RH950, RH800, Rain, Wind2 m, Wind-dir2 m) did not reveal very pronounced trends with Nd in either season consistent with how they did not rank highly in importance (Fig. 13). Of particular interest is Wind2 m, which is used here as a proxy variable for updraft speed in the marine boundary layer, which is expected to have a high impact on Nd via its effect on in-cloud supersaturation. Although the ALE plot of Wind2 m suggested a small increase of about  10 cm−3 in Nd as the wind speed increased, Wind2 m did not come out as a very important parameter in either season. This may be due to the fact that environmental conditions representing updraft speed were already included in parameters such as cloud fraction and CAO index. Another explanation can be the shortcomings and high uncertainties associated with the use of Wind2 m as a proxy for updraft speed.

The results of regression analysis highlight the high sensitivity of Nd to cloud fraction regardless of season. As discussed earlier, this can be attributed largely to two factors: (i) the relationship between cloud type (e.g., stratocumulus, shallow cumulus) and cloud fraction, which can, in turn, influence cloud microphysical properties like Nd, and (ii) uncertainties associated with Nd estimates from satellite observations that can result in negative biases in Nd for low-cloud-coverage conditions. To further test the relative influence of various variables at different cloud fractions, two sensitivity tests with GBRT modeling were conducted using subsets of data with varying cloud fraction (0.2  CFlow-liq. 0.4 and CFlow-liq. 0.7).

Beginning with results for CFlow-liq. 0.7 (Figs. S22–S25), the average R2 scores for validation and test sets for these runs were 0.47/0.39 (DJF/JJA) and 0.49/0.38 (DJF/JJA), respectively. A feature that stands out is that for both DJF and JJA, surface mass concentrations of sulfate became the most important factor. ALE plots presented in Fig. S23 suggest that Nd has a very similar sensitivity to sulfate concentration in high-cloud-coverage conditions regardless of season in contrast to the results of the original run where Nd was more sensitive to the changes in sulfate level in DJF than JJA. These results are in agreement with previous studies where Nd values for marine boundary layer clouds were highly sensitive to sulfate concentrations at the level close to cloud base (Boucher and Lohmann, 1995; Lowenthal et al., 2004; Storelvmo et al., 2009; McCoy et al., 2017, 2018; MacDonald et al., 2020). The second most important factor for DJF was the surface mass concentrations of organic carbon followed by CFlow-liq. and sea salt surface mass concentrations. On the other hand, the second most important factor in JJA was CAO index followed by CFlow-liq. and wind direction. ALE plots presented in Figs. S23–S25 showed similar relationships between Nd and input parameters as observed for the original runs where full datasets were used as the input.

Figure S26 shows the results of the GBRT model using input data with a cloud fraction between 0.2 and 0.4, the condition relatively more prevalent in JJA than DJF. The average R2 scores for validation and test sets for these runs were 0.30/0.30 (DJF/JJA) and 0.33/0.31 (DJF/JJA), respectively. It is interesting to see that for both seasons, an aerosol parameter emerged as the most important factor. Mass concentrations of OC appeared as the most important factor in JJA (the fourth most important factor in DJF) while in DJF, sulfate concentration came out as the most important factor (the fourth most important factor in JJA) consistent with the results of previously discussed models for DJF. It should be noted that ALE plots also suggested less sensitivity of Nd to sulfate in JJA than DJF, similar to the results observed in the original model run including all the data points. The second most important factor in DJF turned out to be the cloud-top effective height of low-level liquid clouds followed by CAO index. On the other hand, CAO index was the second most important factor in JJA followed by PBLH. ALE plots presented in Figs. S27–S29 also showed similar qualitative trends observed in original and high-cloud-coverage runs.

4.3 Unexplored factors

Additional factors impacting the relationship between aerosol and Nd seasonal cycles are discussed here that warrant additional research with more detailed data at finer scales such as with aircraft. We are cognizant that this list is not fully exhaustive. As low-level cloud fraction impacted model results of Sect. 4.2 so substantially, the dynamics of the studied clouds require further characterization. As cloud fraction and CAO index are well related, especially in DJF, aerosol–cloud interactions are likely stronger than other seasons (as implied by Sect. 3.5) due in part to enhanced surface fluxes and turbulence and thus more droplet activation with higher cloud supersaturations (Painemal et al., 2021); in contrast, the smaller shallow cumulus clouds in summertime may be less favorable for droplet activation due to factors such as reduced turbulence and more lateral entrainment.

Entrainment of free tropospheric aerosol can impact Nd values, with potentially varying degrees of influence between seasons. It is presumed that with summertime convection, the more broken cumulus scenes are less adiabatic through the cloudy column and more affected by entrainment and mixing; hence, Nd values derived using data that remote sensors retrieve near cloud top could be considerably lower than values lower by cloud base. Satellite remote sensing studies of aerosol–cloud interactions will presumably be more challenging in winter periods versus the summer with regard to the spatial and temporal mismatch between cloud and aerosol retrievals. More specifically, it is easier to get nearly coincidental sampling in summertime due to lower cloud fractions, while in winter the frontal regions with high cloud fractions make it challenging to get aerosol retrievals. There is complexity in understanding how aerosols relate to Nd due to how giant CCN can reduce Nd and also since wet scavenging can remove aerosols efficiently. As aircraft data are limited and difficult to use for assessing seasonal cycles, new techniques of retrieving CCN and Nd from space will greatly assist such types of studies in the future.

5 Conclusions

This work investigates the seasonal cycle of Nd over the WNAO region in terms of concentration statistics and with discussion of potential influential factors. The results of this work have implications for increased understanding of aerosol–cloud interactions and meteorological factors influencing concentration of cloud droplets in the marine boundary layer. The results and interpretations can be summarized as follows in the order of how they were presented.

  • An ACTIVATE case flight during the DJF season shows a sharp offshore Nd gradient ranging from > 1000 to < 50 cm−3 explained in part by particles smaller than 100 nm activating into drops during a cold-air outbreak with post-frontal clouds. There were significant changes in aerosol composition in cloud-free air and also in droplet residual particles as a function of offshore distance. These changes included a sharp decrease in aerosol number concentration, a decrease in mass fraction of sulfate in droplet residual particles, and an increase in mass fraction of organic and chloride of droplet residual particles moving offshore.

  • Nd is generally highest (lowest) in DJF (JJA) over the WNAO, but aerosol parameters such as AOD, AI, surface-based aerosol mass concentrations for most species, and CCN concentrations (1 % supersaturation) are generally highest in JJA and MAM and are at (or near) their lowest values in DJF. While aerosol extinction in the PBL is highest in DJF, it is driven largely by sea salt (large but few in number) and thus cannot explain the Nd peak in wintertime.

  • While relative humidity was generally highest in JJA across the WNAO, the differences between seasons in the PBL and FT were not sufficiently large to explain the divergent seasonal cycles of AOD and Nd.

  • The susceptibility of Nd to aerosols (Eq. 3) was strongest in DJF using four different proxy variables for aerosols, suggestive of at least one reason why Nd can be highest when aerosol proxy variables for concentration are typically near or at their lowest values.

  • Composite maps of high- versus low-Nd days across the WNAO reveal that conditions associated with the highest-Nd days, regardless of season (but especially DJF), are reduced sea level pressure, stronger winds aligned with continental outflow, high low-level liquid cloud fraction, higher CAO index and PBLH, and enhanced AOD. Cold-air outbreaks are coincident with all of these conditions, especially in the colder months of DJF in sharp contrast to JJA when Nd is lowest.

  • Gradient-boosted regression analysis shows that the most important predictors of Nd in DJF and JJA vary to some extent, but with cloud fraction being the most important parameter, followed by either (for DJF) surface mass concentrations of sulfate and organic carbon and CAO index or (for JJA) CAO index, surface mass concentrations of organic carbon, and sulfate concentrations. Accumulated local effect plots confirm that sulfate and organics help drive the high Nd values via continental outflow, which is assisted in large part by conditions associated with CAOs such as high cloud fraction and high CAO index.

Therefore, the combination of continental pollution outflow and turbulence changes contributed by surface fluxes (manifested in the strongest CAO index values in DJF and weakest in JJA) markedly influence the Nd cycle, leading to differing annual cycles in cloud microphysics and aerosols. More detailed data such as from aircraft and modeling can help extend this line of research to confirm these findings and speculations such as how (i) the aerosol indirect effect is strongest in DJF due to boundary layer dynamics such as with more turbulence and mixing than other seasons (Painemal et al., 2021); (ii) enhanced giant CCN in forms such as sea salt and dust can reduce Nd via expediting the collision–coalescence process; and (iii) substantial aerosol removal can occur far offshore as postfrontal clouds associated with CAOs build and then begin to precipitate. The latter hypothesis may help explain why Bermuda (> 1000 km offshore the US East Coast) was the only selected sub-domain in this study to not have a seasonal Nd peak in DJF.

Code availability

ALE plots were calculated using the source code provided by Jumelle et al. (2021) available at

Data availability

All the data used in this study are open access and can be found online as follows. CERES-MODIS data as described by Minnis et al. (2011) are available at; CALIPSO data are available at (CALIPSO, 2021); PERSIANN-CDR data as described by Nguyen et al. (2019) are available at; MERRA-2 data are available at (GMAO, 2021); TCAP CCN data as described by Berg et al. (2016) are available at; ACTIVATE Airborne Data are available at (NASA/LARC/SD/ASDC, 2020a), (NASA/LARC/SD/ASDC, 2020b), and (NASA/LARC/SD/ASDC, 2020c).​​​​​​​


The supplement related to this article is available online at:

Author contributions

DP assisted with CERES-MODIS data usage. MA performed the CALIPSO analysis, and HD carried out all other data analysis. EC, SK, RHM, CR, MS, KLT, CV, EW, and LZ participated in airborne data collection. AS and HD prepared the manuscript, with DP, MA, MB, SC, AFC, EC, SK, HL, RHM, CR, AJS, MS, KS, KLT, CV, HW, EW, XZ, LZ, and PZ contributing with editing support.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Marine aerosols, trace gases, and clouds over the North Atlantic (ACP/AMT inter-journal SI)”. It is not associated with a conference.


The work was funded by NASA grant 80NSSC19K0442 in support of ACTIVATE, a NASA Earth Venture Suborbital-3 (EVS-3) investigation funded by NASA's Earth Science Division and managed through the Earth System Science Pathfinder Program Office. The authors acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and READY website (, last access: 14 June 2021) used in this work.

Financial support

This research has been supported by the National Aeronautics and Space Administration (grant no. 80NSSC19K0442).

Review statement

This paper was edited by Stefania Gilardoni and reviewed by two anonymous referees.


Abel, S. J., Boutle, I. A., Waite, K., Fox, S., Brown, P. R. A., Cotton, R., Lloyd, G., Choularton, T. W., and Bower, K. N.: The Role of Precipitation in Controlling the Transition from Stratocumulus to Cumulus Clouds in a Northern Hemisphere Cold-Air Outbreak, J. Atmos. Sci., 74, 2293–2314,, 2017. 

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. 

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

Aldhaif, A. M., Lopez, D. H., Dadashazar, H., Painemal, D., Peters, A. J., and Sorooshian, A.: An Aerosol Climatology and Implications for Clouds at a Remote Marine Site: Case Study Over Bermuda, J. Geophys. Res.-Atmos., 126, e2020JD034038,, 2021. 

Apley, D. W. and Zhu, J. Y.: Visualizing the effects of predictor variables in black box supervised learning models, J. Roy. Stat. Soc. B, 82, 1059–1086,, 2020. 

Ashouri, H., Hsu, K.-L., Sorooshian, S., Braithwaite, D. K., Knapp, K. R., Cecil, L. D., Nelson, B. R., and Prat, O. P.: PERSIANN-CDR: Daily Precipitation Climate Data Record from Multisatellite Observations for Hydrological and Climate Studies, B. Am. Meteorol. Soc., 96, 69–83,, 2015. 

Bennartz, R., Fan, J. W., Rausch, J., Leung, L. R., and Heidinger, A. K.: Pollution from China increases cloud droplet number, suppresses rain over the East China Sea, Geophys. Res. Lett., 38, L09704,, 2011. 

Berg, L. K., Fast, J. D., Barnard, J. C., Burton, S. P., Cairns, B., Chand, D., Comstock, J. M., Dunagan, S., Ferrare, R. A., Flynn, C. J., Hair, J. W., Hostetler, C. A., Hubbe, J., Jefferson, A., Johnson, R., Kassianov, E. I., Kluzek, C. D., Kollias, P., Lamer, K., Lantz, K., Mei, F., Miller, M. A., Michalsky, J., Ortega, I., Pekour, M., Rogers, R. R., Russell, P. B., Redemann, J., Sedlacek, A. J., Segal-Rosenheimer, M., Schmid, B., Shilling, J. E., Shinozuka, Y., Springston, S. R., Tomlinson, J. M., Tyrrell, M., Wilson, J. M., Volkamer, R., Zelenyuk, A., and Berkowitz, C. M.: The Two-Column Aerosol Project: Phase IOverview and impact of elevated aerosol layers on aerosol optical depth, J. Geophys. Res.-Atmos., 121, 336–361,, 2016 (data availalbe at:, last access: 14 June 2021). 

Berg, W., L'Ecuyer, T., and van den Heever, S.: Evidence for the impact of aerosols on the onset and microphysical properties of rainfall from a combination of satellite observations and cloud-resolving model simulations, J. Geophys. Res.-Atmos., 113, D14S23,, 2008. 

Boucher, O. and Lohmann, U.: The Sulfate-Ccn-Cloud Albedo Effect – a Sensitivity Study with 2 General-Circulation Models, Tellus B, 47, 281–300,, 1995. 

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S. K., Sherwood, S., Stevens, B., and Zhang, X. Y.: Clouds and aerosols, in: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Doschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, UK, 571–657, 2013. 

Breon, F. M., Tanre, D., and Generoso, S.: Aerosol effect on cloud droplet size monitored from satellite, Science, 295, 834–838,, 2002. 

Brummer, B.: Boundary-layer modification in wintertime cold-air outbreaks from the arctic sea ice, Bound.-Lay. Meteorol., 80, 109–125,, 1996. 

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. 

CALIPSO: CALIPSO LID L2 Standard HDF File – Version 3.30, NASA Langley Research Center Atmospheric Science Data Center DAAC [data set], available at:, last access: 14 June 2021. 

Chin, M., Ginoux, P., Kinne, S., Torres, O., Holben, B. N., Duncan, B. N., Martin, R. V., Logan, J. A., Higurashi, A., and Nakajima, T.: Tropospheric Aerosol Optical Thickness from the GOCART Model and Comparisons with Satellite and Sun Photometer Measurements, J. Atmos. Sci., 59, 461–483,<0461:TAOTFT>2.0.CO;2, 2002. 

Colarco, P., da Silva, A., Chin, M., and Diehl, T.: Online simulations of global aerosol distributions in the NASA GEOS-4 model and comparisons to satellite and ground-based aerosol optical depth, J. Geophys. Res.-Atmos., 115, D14207,, 2010. 

Corral, A. F., Braun, R., Cairns, B., Gorooh, V. A., Liu, H., Ma, L., Mardi, A., Painemal, D., Stamnes, S., van Diedenhoven, B., Wang, H., Yang, Y., Zhang, B., and Sorooshian, A.: An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast – Part 1: Analysis of Aerosols, Gases, and Wet Deposition Chemistry, J. Geophys. Res.-Atmos., 126, e2020JD032592,, 2021. 

Dadashazar, H., Wang, Z., Crosbie, E., Brunke, M., Zeng, X. B., Jonsson, H., Woods, R. K., Flagan, R. C., Seinfeld, J. H., and Sorooshian, A.: Relationships between giant sea salt particles and clouds inferred from aircraft physicochemical data, J. Geophys. Res.-Atmos., 122, 3421–3434,, 2017. 

Dadashazar, H., Crosbie, E., Majdi, M. S., Panahi, M., Moghaddam, M. A., Behrangi, A., Brunke, M., Zeng, X., Jonsson, H. H., and Sorooshian, A.: Stratocumulus cloud clearings: statistics from satellites, reanalysis models, and airborne measurements, Atmos. Chem. Phys., 20, 4637–4665,, 2020. 

DeCarlo, P. F., Dunlea, E. J., Kimmel, J. R., Aiken, A. C., Sueper, D., Crounse, J., Wennberg, P. O., Emmons, L., Shinozuka, Y., Clarke, A., Zhou, J., Tomlinson, J., Collins, D. R., Knapp, D., Weinheimer, A. J., Montzka, D. D., Campos, T., and Jimenez, J. L.: Fast airborne aerosol size and chemistry measurements above Mexico City and Central Mexico during the MILAGRO campaign, Atmos. Chem. Phys., 8, 4027–4048,, 2008. 

Deuze, J. L., Breon, F. M., Devaux, C., Goloub, P., Herman, M., Lafrance, B., Maignan, F., Marchand, A., Nadal, F., Perry, G., and Tanre, D.: Remote sensing of aerosols over land surfaces from POLDER-ADEOS-1 polarized measurements, J. Geophys. Res.-Atmos., 106, 4913–4926,, 2001. 

Feingold, G.: Modeling of the first indirect effect: Analysis of measurement requirements, Geophys. Res. Lett., 30, 1997,, 2003. 

Fletcher, J., Mason, S., and Jakob, C.: The Climatology, Meteorology, and Boundary Layer Structure of Marine Cold Air Outbreaks in Both Hemispheres, J. Climate, 29, 1999–2014,, 2016. 

Fuchs, J., Cermak, J., and Andersen, H.: Building a cloud in the southeast Atlantic: understanding low-cloud controls based on satellite observations with machine learning, Atmos. Chem. Phys., 18, 16537–16552,, 2018. 

Gelaro, R., McCarty, W., Suarez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G. K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. 

Global Modeling and Assimilation Office (GMAO): MERRA-2, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC), available at:, last access: 14 June 2021. 

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. 

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., Hunerbein, 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. B., 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. 

Gryspeerdt, E., Quaas, J., and Bellouin, N.: Constraining the aerosol influence on cloud fraction, J. Geophys. Res.-Atmos., 121, 3566–3583,, 2016. 

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. L., Wang, M. H., and Zhang, K.: Constraining the instantaneous aerosol influence on cloud albedo, P. Natl. Acad. Sci. USA, 114, 4899–4904,, 2017. 

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. 

Hastie, T., Tibshirani, R., and Friedman, J.: The elements of statistical learning: data mining, inference and prediction, 2nd edn., Springer, New York City, USA, 744 pp., 2009. 

Jumelle, M., Kuhn-Regnier, A., Rajaratnam, S., and Mahima, A.: Python Accumulated Local Effects package, GitHub, available at:, last access: 14 June 2021. 

Kolstad, E. W., Bracegirdle, T. J., and Seierstad, I. A.: Marine cold-air outbreaks in the North Atlantic: temporal distribution and associations with large-scale atmospheric circulation, Clim. Dynam., 33, 187–197,, 2009. 

Liu, J. and Li, Z.: Aerosol properties and their influences on low warm clouds during the Two-Column Aerosol Project, Atmos. Chem. Phys., 19, 9515–9529,, 2019. 

Lloyd, G., Choularton, T. W., Bower, K. N., Gallagher, M. W., Crosier, J., O'Shea, S., Abel, S. J., Fox, S., Cotton, R., and Boutle, I. A.: In situ measurements of cloud microphysical and aerosol properties during the break-up of stratocumulus cloud layers in cold air outbreaks over the North Atlantic, Atmos. Chem. Phys., 18, 17191–17206,, 2018. 

Loeb, N. G., Manalo-Smith, N., Su, W. Y., Shankar, M., and Thomas, S.: CERES Top-of-Atmosphere Earth Radiation Budget Climate Data Record: Accounting for in-Orbit Changes in Instrument Calibration, Remote Sens.-Basel, 8, 182,, 2016. 

Lowenthal, D. H., Borys, R. D., Choularton, T. W., Bower, K. N., Flynn, M. J., and Gallagher, M. W.: Parameterization of the cloud droplet-sulfate relationship, Atmos. Environ., 38, 287–292,, 2004. 

MacDonald, A. B., Hossein Mardi, A., Dadashazar, H., Azadi Aghdam, M., Crosbie, E., Jonsson, H. H., Flagan, R. C., Seinfeld, J. H., and Sorooshian, A.: On the relationship between cloud water composition and cloud droplet number concentration, Atmos. Chem. Phys., 20, 7645–7665,, 2020. 

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

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

McCoy, D. T., Bender, F. A.-M., Grosvenor, D. P., Mohrmann, J. K., Hartmann, D. L., Wood, R., and Field, P. R.: Predicting decadal trends in cloud droplet number concentration using reanalysis and satellite data, Atmos. Chem. Phys., 18, 2035–2047,, 2018. 

Minnis, P., Sun-Mack, S., Young, D. F., Heck, P. W., Garber, D. P., Chen, Y., Spangenberg, D. A., Arduini, R. F., Trepte, Q. Z., Smith, W. L., Ayers, J. K., Gibson, S. C., Miller, W. F., Hong, G., Chakrapani, V., Takano, Y., Liou, K. N., Xie, Y., and Yang, P.: CERES Edition-2 Cloud Property Retrievals Using TRMM VIRS and Terra and Aqua MODIS Data – Part I: Algorithms, IEEE T. Geosci. Remote, 49, 4374–4400,, 2011 (data available at:, last access: 14 June 2021). 

Minnis, P., Sun-Mack, S., Chen, Y., Chang, F., Yost, C. R., Smith, W. L., Heck, P. W., Arduini, R. F., Bedka, S. T., Yi, Y., Hong, G., Jin, Z., Painemal, D., Palikonda, R., Scarino, B. R., Spangenberg, D. A., Smith, R. A., Trepte, Q. Z., Yang, P., and Xie, Y.: CERES MODIS Cloud Product Retrievals for Edition 4 – Part I: Algorithm Changes, IEEE T. Geosci. Remote, 59, 1–37,, 2020. 

Molnar, C.: Interpretable Machine Learning, A Guide for Making Black Box Models Explainable, available at: (last access: 14 June 2021), 2019. 

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and natural radiative forcing, in: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Doschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, UK, 659–740, 2013. 

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. 

NASA/LARC/SD/ASDC: ACTIVATE Falcon In Situ Aerosol Data [data set],, 2020a. 

NASA/LARC/SD/ASDC: ACTIVATE Falcon In Situ Cloud Data [data set],, 2020b. 

NASA/LARC/SD/ASDC: ACTIVATE Falcon In-Situ Meteorological and Navigational Data [data set],, 2020c. 

Naud, C. M., Booth, J. F., and Lamraoui, F.: Post Cold Frontal Clouds at the ARM Eastern North Atlantic Site: An Examination of the Relationship Between Large-Scale Environment and Low-Level Cloud Properties, J. Geophys. Res.-Atmos., 123, 12117–12132,, 2018. 

Nguyen, P., Ombadi, M., Sorooshian, S., Hsu, K., AghaKouchak, A., Braithwaite, D., Ashouri, H., and Thorstensen, A. R.: The PERSIANN family of global satellite precipitation data: a review and evaluation of products, Hydrol. Earth Syst. Sci., 22, 5801–5816,, 2018. 

Nguyen, P., Shearer, E. J., Tran, H., Ombadi, M., Hayatbini, N., Palacios, T., Huynh, P., Braithwaite, D., Updegraff, G., Hsu, K., Kuligowski, B., Logan, W. S., and Sorooshian, S.: The CHRS Data Portal, an easily accessible public repository for PERSIANN global satellite precipitation data, Scientific Data, 6, 1–10,, 2019 (data available at:, last access: 14 June 2021). 

Painemal, D.: Global Estimates of Changes in Shortwave Low-Cloud Albedo and Fluxes Due to Variations in Cloud Droplet Number Concentration Derived From CERES-MODIS Satellite Sensors, Geophys. Res. Lett., 45, 9288–9296,, 2018. 

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. 

Painemal, D., Corral, A. F., Sorooshian, A., Brunke, M. A., Chellappan, S., Gorooh, V. A., Ham, S., O'Neill, L., Smith Jr., W. L., Tselioudis, G., Wang, H., Zeng, X., and Zuidema, P.: An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast – Part 2: Circulation, Boundary Layer, and Clouds, J. Geophys. Res.-Atmos., 126, e2020JD033423,, 2021. 

Papritz, L., Pfahl, S., Sodemann, H., and Wernli, H.: A Climatology of Cold Air Outbreaks and Their Impact on Air-Sea Heat Fluxes in the High-Latitude South Pacific, J. Climate, 28, 342–364,, 2015. 

Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., and Duchesnay, E.: Scikit-learn: Machine Learning in Python, J. Mach. Learn. Res., 12, 2825–2830, 2011. 

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

Persson, C., Bacher, P., Shiga, T., and Madsen, H.: Multi-site solar power forecasting using gradient boosted regression trees, Sol. Energy, 150, 423–436,, 2017. 

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. 

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. 

Quinn, P. K., Coffman, D. J., Johnson, J. E., Upchurch, L. M., and Bates, T. S.: Small fraction of marine cloud condensation nuclei made up of sea spray aerosol, Nat. Geosci., 10, 674–679,, 2017. 

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. 

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. 

Rolph, G., Stein, A., and Stunder, B.: Real-time Environmental Applications and Display sYstem: READY, Environ. Model. Softw., 95, 210–228,, 2017. 

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, 4699,, 2003. 

Shingler, T., Dey, S., Sorooshian, A., Brechtel, F. J., Wang, Z., Metcalf, A., Coggon, M., Mülmenstädt, J., Russell, L. M., Jonsson, H. H., and Seinfeld, J. H.: Characterisation and airborne deployment of a new counterflow virtual impactor inlet, Atmos. Meas. Tech., 5, 1259–1269,, 2012. 

Sinclair, K., van Diedenhoven, B., Cairns, B., Alexandrov, M., Moore, R., Ziemba, L. D., and Crosbie, E.: Observations of Aerosol-Cloud Interactions During the North Atlantic Aerosol and Marine Ecosystem Study, Geophys. Res. Lett., 47, e2019GL085851,, 2020. 

Sorooshian, A., Anderson, B., Bauer, S. E., Braun, R. A., Cairns, B., Crosbie, E., Dadashazar, H., Diskin, G., Ferrare, R., Flagan, R. C., Hair, J., Hostetler, C., Jonsson, H. H., Kleb, M. M., Liu, H. Y., MacDonald, A. B., McComiskey, A., Moore, R., Painemal, D., Russell, L. M., Seinfeld, J. H., Shook, M., Smith, W. L., Thornhill, K., Tselioudis, G., Wang, H. L., Zeng, X. B., Zhang, B., Ziemba, L., and Zuidema, P.: Aerosol-Cloud-Meteorology Interaction Airborne Field Investigations: Using Lessons Learned from the US West Coast in the Design of ACTIVATE off the US East Coast, B. Am. Meteorol. Soc., 100, 1511–1528,, 2019. 

Sorooshian, A., Corral, A. F., Braun, R. A., Cairns, B., Crosbie, E., Ferrare, R., Hair, J., Kleb, M. M., Mardi, A. H., Maring, H., McComiskey, A., Moore, R., Painemal, D., Scarino, A. J., Schlosser, J., Shingler, T., Shook, M., Wang, H. L., Zeng, X. B., Ziemba, L., and Zuidema, P.: Atmospheric Research Over the Western North Atlantic Ocean Region and North American East Coast: A Review of Past Work and Challenges Ahead, J. Geophys. Res.-Atmos., 125, e2019JD031626,, 2020. 

Stein, A. F., Draxler, R. R., Rolph, G. D., Stunder, B. J. B., Cohen, M. D., and Ngan, F.: NOAA's HYSPLIT Atmospheric Transport and Dispersion Modeling System, B. Am. Meteorol. Soc., 96, 2059–2077,, 2015. 

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

Storelvmo, T., Lohmann, U., and Bennartz, R.: What governs the spread in shortwave forcings in the transient IPCC AR4 models?, Geophys. Res. Lett., 36, L01806,, 2009. 

Tackett, J. L., Winker, D. M., Getzewich, B. J., Vaughan, M. A., Young, S. A., and Kar, J.: CALIPSO lidar level 3 aerosol profile product: version 3 algorithm design, Atmos. Meas. Tech., 11, 4129–4152,, 2018. 

Terai, C. R., Wood, R., and Kubar, T. L.: Satellite estimates of precipitation susceptibility in low-level marine stratiform clouds, J. Geophys. Res.-Atmos., 120, 8878–8889,, 2015. 

Thornhill, K. L., Anderson, B. E., Barrick, J. D. W., Bagwell, D. R., Friesen, R., and Lenschow, D. H.: Air motion intercomparison flights during Transport and Chemical Evolution in the Pacific (TRACE-P)/ACE-ASIA, J. Geophys. Res.-Atmos., 108, 8783,, 2003.  

Twomey, S.: Influence of Pollution on Shortwave Albedo of Clouds, J. Atmos. Sci., 34, 1149–1152,<1149:Tiopot>2.0.Co;2, 1977. 

Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y. X., Powell, K. A., Liu, Z. Y., Hunt, W. H., and Young, S. A.: Overview of the CALIPSO Mission and CALIOP Data Processing Algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323,, 2009. 

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

Zhang, Z., Ackerman, A. S., Feingold, G., Platnick, S., Pincus, R., and Xue, H.: Effects of cloud horizontal inhomogeneity and drizzle on remote sensing of cloud droplet effective radius: Case studies based on large-eddy simulations, J. Geophys. Res.-Atmos., 117, D19208,, 2012. 

Zhang, Z., Werner, F., Cho, H.-M., Wind, G., Platnick, S., Ackerman, A. S., Di Girolamo, L., Marshak, A., and Meyer, K.: A framework based on 2-D Taylor expansion for quantifying the impacts of subpixel reflectance variance and covariance on cloud optical thickness and effective radius retrievals based on the bispectral method, J. Geophys. Res.-Atmos., 121, 7007–7025,, 2016. 

Short summary
This study investigates the seasonal cycle of cloud drop number concentration (Nd) over the western North Atlantic Ocean (WNAO) using multiple datasets. Reasons for the puzzling discrepancy between the seasonal cycles of Nd and aerosol concentration were identified. Results indicate that Nd is highest in winter (when aerosol proxy values are often lowest) due to conditions both linked to cold-air outbreaks and that promote greater droplet activation.
Final-revised paper