Marine boundary layer aerosol in Eastern North Atlantic: seasonal variations and key controlling processes

The response of marine low cloud systems to changes in aerosol concentration represents one of the largest uncertainties in 15 climate simulations. Major contributions to this uncertainty derive from poor understanding of aerosol under natural conditions and the perturbation by anthropogenic emissions. The Eastern North Atlantic (ENA) is a region of persistent but diverse marine boundary layer (MBL) clouds, whose albedo and precipitation are highly susceptible to perturbations in aerosol properties. In this study, we examine MBL aerosol properties, trace gas mixing ratios, and meteorological parameters measured at the Atmospheric Radiation Measurement Climate Research Facility’s ENA site on Graciosa Island, Azores, Portugal, during a three20 year period from 2015 to 2017. Measurements impacted by local pollutions on Graciosa Island and during occasional intense biomass burning and dust events are excluded from this study. Submicron aerosol size distribution typically consists of three modes: Aitken (At, diameter Dp < ~100 nm), Accumulation (Ac, Dp within ~100 to ~300 nm), and Larger Accumulation (LA, Dp > ~300 nm) modes, with average number concentrations (denoted as NAt, NAc and NLA below) of 330, 114, and 14 cm-3, respectively. NAt, NAc and NLA show contrasting seasonal variations, suggesting different sources and removal processes. NLA is 25 dominated by sea spray aerosol (SSA), and is higher in winter and lower in summer. This is due to the seasonal variations of SSA production, in-cloud coalescence scavenging, and dilution by entrained free troposphere (FT) air. In comparison, SSA typically contributes a relatively minor fraction to NAt (10 %) and NAc (21 %) on an annual basis. In addition to SSA, sources of Ac mode particles include entrained of FT aerosols and condensation growth of Aitken mode particles inside MBL, while incloud coalescence scavenging is the major sink of NAc. The observed seasonal variation of NAc, being higher in summer and 30 lower in winter, generally agrees with the steady-state concentration estimated from major sources and sinks. NAt is mainly controlled by entrainment of FT aerosol, coagulation loss, and growth of Aitken mode particles into the Ac mode size range. Our calculation suggests that besides the direct contribution from entrained FT Ac mode particles, growth of entrained FT Aitken mode particles in the MBL also represent a substantial source of cloud condensation nuclei (CCN), with the highest contribution potentially reaching 60 % during summer. The growth of Aitken mode particles to CCN size is an expected result of the 35 condensation of sulfuric acid, a product from dimethyl sulfide oxidation, suggesting that ocean ecosystems may have a substantial influence on MBL CCN populations in ENA.

Abstract. The response of marine low cloud systems to changes in aerosol concentration represents one of the largest uncertainties in climate simulations. Major contributions to this uncertainty are derived from poor understanding of aerosol under natural conditions and the perturbation by anthropogenic emissions. The eastern North Atlantic (ENA) is a region of persistent but diverse marine boundary layer (MBL) clouds, whose albedo and precipitation are highly susceptible to perturbations in aerosol properties. In this study, we examine MBL aerosol properties, trace gas mixing ratios, and meteorological parameters measured at the Atmospheric Radiation Measurement Climate Research Facility's ENA site on Graciosa Island, Azores, Portugal, during a 3year period from 2015 to 2017. Measurements impacted by local pollution on Graciosa Island and during occasional intense biomass burning and dust events are excluded from this study. Submicron aerosol size distribution typically consists of three modes: Aitken (At, diameter D p <∼ 100 nm), accumulation (Ac, D p within ∼ 100 to ∼ 300 nm), and larger accumulation (LA, D p >∼ 300 nm) modes, with average number concentrations (denoted as N At , N Ac , and N LA below) of 330, 114, and 14 cm −3 , respectively. N At , N Ac , and N LA show contrasting seasonal variations, suggesting different sources and removal processes. N LA is dominated by sea spray aerosol (SSA) and is higher in winter and lower in summer. This is due to the seasonal variations of SSA production, in-cloud coalescence scavenging, and dilution by entrained free troposphere (FT) air. In comparison, SSA typically contributes a relatively minor fraction to N At (10 %) and N Ac (21 %) on an annual basis. In addition to SSA, sources of Acmode particles include entrainment of FT aerosols and condensation growth of Aitken-mode particles inside the MBL, while in-cloud coalescence scavenging is the major sink of N Ac . The observed seasonal variation of N Ac , being higher in summer and lower in winter, generally agrees with the steady-state concentration estimated from major sources and sinks. N At is mainly controlled by entrainment of FT aerosol, coagulation loss, and growth of Aitken-mode particles into the Ac-mode size range. Our calculation suggests that besides the direct contribution from entrained FT Ac-mode particles, growth of entrained FT Aitken-mode particles in the MBL also represent a substantial source of cloud condensation nuclei (CCN), with the highest contribution potentially reaching 60 % during summer. The growth of Aitken-mode particles to CCN size is an expected result of the condensation of sulfuric acid, a product from dimethyl sulfide oxidation, suggesting that ocean ecosystems may have a substantial influence on MBL CCN populations in the ENA.

Introduction
Low clouds, especially stratocumulus, are the dominant cloud type in terms of spatial coverage of the Earth's surface, and are of vital importance to the Earth's climate . Major climate effects of low clouds are derived from their reflection of solar radiation (Seinfeld and Pandis, 2016).
The key parameters for quantifying climate effects of low clouds are the albedo (i.e., cloud reflectivity) and the cloud coverage, both of which are particularly sensitive to perturbations of aerosols. The concentration of cloud condensation nuclei (CCN) strongly influences the number concentration and sizes of cloud droplets and therefore the effective albedo of low clouds (i.e., first indirect effect of aerosol) (Twomey, 1974;Seinfeld and Pandis, 2016;Dong et al., 2015), especially in clean environments such as the remote marine boundary layer (MBL) (Reutter et al., 2009). In addition, CCN concentration and aerosol size distribution also influence cloud amount by impacting drizzle formation and precipitation (i.e., second indirect effect of aerosol) (Albrecht, 1989;Liu and Daum, 2004;Liu et al., 2006;Wood, 2005;Rémillard et al., 2012;Dong et al., 2014).
Currently, the aerosol indirect effects of marine low cloud systems remain one of the major uncertainties in climate change simulations (Lohmann and Feichter, 2005;Bony and Dufresne, 2005;Bony et al., 2006;. This large uncertainty is to a large degree a result of the incomplete understanding, and therefore representations, of aerosol properties, and the response of marine low clouds to aerosol changes. Therefore, it is imperative to understand MBL aerosol properties under natural conditions, the perturbation due to anthropogenic emissions, and the underlying controlling processes. The properties of aerosols in the remote MBL can be influenced by a variety of processes, including entrainment from the free troposphere (FT), production of sea spray aerosol (SSA), processing of aerosol particles both inside clouds and in clear air, depositions, and horizontal advection (Quinn and Bates, 2011;. Previous studies (O'Dowd et al., 2004;Clarke et al., 2013;Quinn et al., 2017;Wood et al., 2017;Pierce et al., 2015;Prather et al., 2013;Russell et al., 2010;Sanchez et al., 2018;Phinney et al., 2006;Langley et al., 2010) have greatly advanced our understanding of MBL aerosols, especially in the relative contributions of SSA versus long-range transported pollution in terms of the CCN budget (Blot et al., 2013;Clarke and Kapustin, 2010;Clarke et al., 2013;Quinn et al., 2017), and the removal of CCN by in-cloud coalescence scavenging . However, we are still lacking a quantitative understanding of the controlling processes sufficient to serve as a reliable foundation for developing global climate model parameterizations and representations that will adequately simulate aerosol in past, current, and future climates. The relative importance, the influence on different particle size ranges, and spatiotemporal variations of these processes are still not well quantified.
The eastern North Atlantic (ENA) is a region of persistent but diverse subtropical MBL clouds (Wood et al., 2015). Aerosols arriving in the ENA are of diverse origins, varying from marine clean air masses to air masses that are strongly influenced by continental emissions from North America or northern Europe (O'Dowd and Smith, 1993;Wood et al., 2015). As a result, ENA is among the regions with strong but uncertain aerosol indirect forcing (Carslaw et al., 2013). Several field campaigns, including the North Atlantic Regional Experiment (NARE) campaign during 1991 to 2001 (Parrish et al., 1998), the Atlantic Stratocumulus Transition Experiment (ASTEX) during June 1992, the second Aerosol Characterization Experiment (ACE-2) during summer 1997, and the Clouds, Aerosol, and Precipitation in the Marine Boundary Layer (CAP-MBL) campaign (Wood et al., 2015) from May 2009 to December 2010 took place in the ENA. However, they are either more focused on other subjects (e.g., ozone chemistry for NARE (Parrish et al., 1998) and cloud properties for CAP-MBL (Wood et al., 2015) or are short-term studies (e.g., ACE-2; Raes et al., 2000, AS-TEX;Albrecht et al., 1995). To our knowledge, the variation of aerosol properties and their controlling processes have not been systematically studied using long-term observation in the ENA.
Recently, a permanent ENA site was established by the Department of Energy Atmospheric Radiation Measurement (ARM) (Mather and Voyles, 2013) Climate Research Facility on Graciosa Island in the Azores, Portugal, providing an invaluable opportunity to study MBL aerosol properties and their interactions with low clouds. In this study, we examine the long-term variation of aerosol properties, trace gas mixing ratios, and meteorological parameters measured at the ARM ENA site from 2015 to 2017 (Sect. 2). The characteristics of the aerosol properties and their seasonal variations are summarized (Sect. 3). The governing equations of number concentration are established for different modes of MBL aerosol at the ENA site (Sect. 4). Subsequently, the seasonal variations of aerosol properties for different particle size modes are explained using key processes identified (Sects. 5 and 6). Finally, we present an overall picture of the processes that drive MBL aerosol properties in the ENA, and the implications are discussed (Sect. 7).

Measurement overview
Measurements of trace gases, meteorological parameters, and aerosol and cloud properties are conducted at the ENA site, located on Graciosa Island in the Azores, Portugal (39 • 5 30 N, 28 • 1 32 W; 30.48 m above mean sea level). The ENA site was initially set up in late 2013, with additional measurements added subsequently. The primary measurements used in this study and the available time periods are listed in Table 1. The measurements of trace gases (e.g., CO) and aerosol properties were first screened for impact from local pollution sources (see Sect. S1 in the Supplement). All measurements are then averaged into 1 h intervals. Here, we use 3 years of data from January 2015 to December 2017 to show the long-term variations and correlations among different parameters. For evaluation of the contributions of different controlling processes (Sect. 4), 1 year of data from September 2016 to August 2017 is used, during which most of the measurements are available.

Optical properties
Aerosol absorbing (B abs ) and scattering (B sca ) coefficients are measured by a three-wavelength Particle Soot Absorption Photometer (PSAP) and a nephelometer, respectively (Table 1). These two instruments share a common inlet, and the 50 % cut size of the inlet switches between 1 and 10 µm every hour (Springston, 2016). The corresponding B sca and B abs are denoted by "PM 1 " and "PM 10 ", respectively. In addition, properties of coarse-mode (1 < D p < 10 µm) aerosols, PM c , were derived by the difference between PM 10 and PM 1 . For example, "PM c B sca " refers to the difference between PM 10 B sca and PM 1 B sca hereafter, and PM c B abs is defined similarly. The mass flow calibration and filter loading correction are already applied to the PSAP data in the ARM data archive (Springston, 2016). In this study, additional corrections of contribution due to scattering for B abs (Bond et al., 1999;Virkkula et al., 2005;Virkkula, 2010;Costabile et al., 2013), and truncation and angular illumination for B sca (Anderson and Ogren, 1998;Müller et al., 2011) are applied, and the procedure is detailed in Sect. S2. The corrected PM 1 B sca shows strong correlation (correlation coefficient being 0.84) with the volume of PM 1 derived from the Ultra-High Sensitivity Aerosol Spectrometer (UHSAS) size distribution (Fig. S2). Potential particle losses for large particles (i.e., in the diameter range of 5-10 µm) are not corrected. However, we do not expect the losses to affect the relative trends of PM c B sca presented here (Sect. 5), or the correlation among PM c B sca and V LA (Fig. 6c).

Cloud and MBL properties
Properties of cloud and MBL, including MBL height, cloud thickness, and cloud fraction, are needed to estimate some of the key controlling processes that drive aerosol properties (see Sect. 4 for more details). The MBL height, H MBL , is derived from the backscatter signal from the ceilometer CL31 (Morris, 2012). Briefly, it is determined from the gradient of an idealized backscatter profile, the parameters of which are derived from fitting of the observed profile (Eresmaa et al., 2006). As the first boundary layer height given in the ceilometer data product is usually the surface layer (Lewis and Schwartz, 2004) below 100 m Emeis et al., 2007Emeis et al., , 2008Haeffelin et al., 2012;Morris, 2012), H MBL is chosen as the highest boundary layer height below 3 km (Zhou et al., 2015;Rémillard and Tselioudis, 2015;Rémillard et al., 2012).
Cloud thickness h is derived by combining H MBL and the cloud-base height derived from ceilometer data. In the ENA, H MBL usually represents the top height of boundary layer clouds (Rémillard et al., 2012). When multiple layers of clouds are detected, the layers with cloud-base heights higher than H MBL are first excluded, after which the highest layer is chosen to exclude potential influence of near-ground thin clouds. The cloud thickness h is then defined as the difference between H MBL (cloud top) and the base height of the chosen cloud. The value of h derived using the above approach is in general agreement with previous observations (Rémillard et al., 2012).
The cloud fraction, p cloud , is determined by the detection status information from the ceilometer (Morris, 2012). It is equal to the fraction of time with a detected boundary layer cloud base, or a determined full obscuration. Precipitation rate at cloud base is retrieved from the vertically pointing K-band cloud radar (Atmospheric Radiation Measurement Climate Research Facility, 1990) and the ceilometer (Morris, 2012) following the method of O' Connor et al. (2005).
3 Seasonal variation in synoptic conditions, trace gas mixing ratios, and aerosol properties in the ENA

Air mass origin
One major source of MBL aerosol in the ENA is the entrainment of FT air, which contains both particles from longrange transport of continental pollution and those formed through new particle formation (NPF) in the FT (Quinn and Bates, 2011;Sanchez et al., 2018). To examine the contribution from continental emissions and its seasonal variation, we analyze the back trajectories of air masses arriving at the ENA site. The cluster analysis results of four representative months from September 2016 to August 2017 (i.e., the main study period; see Sect. 2.1) are shown in Fig. 1. Results from other periods from 2015 to 2017 are similar (not shown). Most of the air masses arriving at the ENA site can be classified as one of the four clusters originating from North America, northern Europe, the Arctic, and the recirculating flow around the Azores High, respectively (O'Dowd and Smith, 1993;Wood et al., 2015). Among these clusters, the Azores High air masses usually linger within the MBL, as indicated by their stable and low-level trajectories (e.g., blue and red trajectory clusters in Fig. 1b). In comparison, other air masses usually undergo long-range transport within the FT before descending into the MBL. In addition, some air masses originating in the continental boundary layer were Figure 1. Cluster analysis of 10-day back trajectories arriving at 100 m above the ENA site in different seasons. The analysis was conducted using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) version 4 model (Stein et al., 2015). The 10-day back trajectories were simulated with a time step of 6 h using National Centers for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS) meteorological data as input. A cluster analysis of these trajectories was then performed, and for each season, the solution that captures most of the variance (e.g., Abdalmogith and Harrison, 2005) and has less than five identified clusters is chosen. The average trajectories of the clusters are represented by different colors, and the associated numbers denoted the arbitrarily given cluster ID and the occurrence percentages of this cluster. For example, the number "1" (90 %) beside the red trajectories indicates that the no. 1 cluster has an average trajectory shown by the red lines, and 90 % of the time the air masses arriving at the ENA site belong to this cluster. lofted up and then subsided into the MBL within 10 days (e.g., blue trajectory cluster in Fig. 1c).
The percentage of occurrence for each cluster shows strong seasonal variations (Figs. 1, S3). During fall (Figs. 1c, S3c) and winter (Figs. 1d, S3d), air masses influenced by anthropogenic emissions from North American (red lines) dominate, with the influence of clean maritime flow and northern European flow. In spring (Figs. 1a, S3a), contributions from Arctic/northern Europe air masses are more pronounced than during other seasons. For the summer months (Figs. 1b, S3b), the ENA site is dominated by the clean maritime flow associated with the recirculating Azores High. As the recirculating Azores High clusters are usually not associated with long-range transport, a reduced contribution to MBL aerosol from continental pollution through FT entrainment is expected during the summertime in ENA.

Mixing ratios of CO, O 3 , and water vapor
The mixing ratios of CO, O 3 , and water vapor within the MBL are expected to be strongly influenced by entrainment of FT air in the ENA. CO is a long-lived species with a lifetime of approximately 1 month (Seinfeld and Pandis, 2016) and therefore is a good indicator of long-range transported continental emissions for remote sites. At the ENA site, the influence of local emissions on trace gases and aerosol measurements is expected to be minimal after filtering of the data (Sect. S1). The lifetime of O 3 varies from hours in polluted urban regions (due to the high rate of photochemical reactions) to several weeks in the FT (Monks et al., 2015). Given its long lifetime in the FT, O 3 may also serve as a tracer for long-range transported pollutants. The local photochemical activities can be inferred from the correlation between O 3 and CO. In regions with strong local sources and sinks, O 3 and CO show a strong positive correlation during summer daytime due to photochemical reactions but a negative correlation during winter nights due to the stronger dry deposition of O 3 than CO (Poulida et al., 1991;Chin et al., 1994). In contrast, at the ENA site, CO and O 3 are positively correlated all year round, even in winter in nighttime with low wind speed (WS) < 2 m s −1 (Mao and Talbot, 2004) (Fig. S4a). This suggests that the variation of O 3 concentration observed at the ENA site is mainly influenced by the entrainment of FT air, in agreement with findings from previous modeling studies (Cooper et al., 2002;Voulgarakis et al., 2011). In addition, the strong anticorrelation (correlation coefficient being −0.75) of CO and O 3 with water vapor (Fig. S4b) also confirms this picture, as water mixing ratio usually negatively correlates with the extent of FT entrainment at remote marine sites (Helmig et al., 2002). Furthermore, the seasonal variations of O 3 and CO in the ENA (Fig. 2a, b) differ much from those observed at anthropogenic-influenced urban or rural sites, where ozone usually exhibits a summer peak due to strong photochemical production, while CO usually shows no clear seasonal variation (Poulida et al., 1991). In contrast, both CO and O 3 in the ENA show a summer minimum and spring-winter maximum, which is consistent with the FT entrainment as the dominant source and the seasonal variation of tropopause height. This suggests minor contributions from local emissions and in situ photochemistry (Parrish et al., 1998;Fischer et al., 2003;Mao and Talbot, 2004). The seasonal variations of CO and O 3 concentrations are also consistent with the cluster analysis of back trajectories, which indicates more influence from long-range transported pollution in winter-spring than in summer.

Absorbing aerosols
In the ENA boundary layer, absorbing aerosols, including black carbon, brown carbon, and dust, are likely entrained from the FT following transport from continental sources. Occasionally, air masses with very strong influence from biomass burning or dust are observed at the ENA site. These episodes are excluded from the analyses presented here to focus on the long-term background variations. These episodes are identified using the aerosol optical properties (Logan et al., 2013(Logan et al., , 2014Cazorla et al., 2013), particle chemical compositions (Clarke et al., 2007), and trace gas mixing ratios (Honrath et al., 2004). Identification of these dust and biomass burning episodes and characterization of aerosol properties during them will be discussed elsewhere. With these episodes excluded, the equivalent black carbon (EBC, following the naming convection suggested by Petzold et al., 2013) mass concentrations were estimated from PM 1 B abs with an assumed mass absorbing cross section of 7.5 m 2 g −1 at 529 nm (Bond et al., 2013).
While absorbing particles are entrained from the FT, the seasonal variation of EBC mass concentration is different from those of CO and O 3 (Fig. 2). As evidenced from a decreasing EBC/CO ratio with increasing precipitation rate at cloud base, P CB (Fig. S5), such differences are due to in-cloud coalescence scavenging (Sect. 4.2) both during the long-range transport and/or after entrainment into the MBL, which removes EBC but not CO or O 3 . Therefore, EBC mass concentration can be indicative of the overall effect of FT contribution from continental emissions and in-cloud coalescence scavenging. As shown in Fig. 2c, EBC mass concentrations are similar in all seasons but show larger annual variations than CO or O 3 , which are attributed to the larger annual variations of precipitation.

Modes of aerosol size distributions
The aerosol size distribution from 70 nm to 1 µm at the ENA site typically consists of three modes (Fig. 3): an Aitken (At) mode below ∼ 100 nm, an accumulation mode (Ac) which resides mostly from 100 to 300 nm, and a larger accumulation mode (LA) above ∼ 300 nm. Note that, due to the lower size limit of UHSAS, the Aitken mode is often not fully characterized. Therefore, its number concentration is derived by deducting fitted number concentrations of the other two modes from the total number concentration CN measured by the condensation particle counter (CPC), namely N At = CN-N Ac -N LA . With this definition, the derived Aitken-mode concentration also includes nucleation-mode particles (i.e., D p < 20 nm). However, previous studies have shown that NPF events within remote MBLs like the ENA are infrequent (Raes, 1995;Bates et al., 2000); therefore, nucleationmode particles likely represent a small fraction of the derived Aitken-mode number concentration for long-term measurements . The Ac mode is absent in 15 % of cases (Table 2), likely due to in-cloud coalescence scavenging or lack of cloud processing (Sect. 4). Among these three modes, aerosol number concentration is dominated by At (72.0 %) and Ac (24.9 %) modes (Fig. 3b1), while the volume concentration is controlled by the LA (74.3 %) and Ac (25.1 %) modes (Fig. 3b2). Based on the average volume size distributions (Fig. 3b2) and results shown in Sect. 5, the LA mode is essentially the sea spray aerosol coarse mode under vast majorities of the conditions.

Seasonal variations of each mode
Different seasonal variations are observed for the three particle modes. While there is substantial variation within each season, on average, the Ac mode exhibits higher number concentration, larger mode D p , and higher occurrence in summer than in winter (Table 2). In contrast, the LA mode shows opposite seasonal trends, with the number and volume concentrations in winter 1.5 times greater than those in summer (Table 2). These seasonal trends are also evident in the seasonally averaged size distributions (Fig. 4a). The monthly average concentrations and the seasonal trends of the Ac and LA modes are very consistent from 2015 to 2017, showing little annual variation (Fig. 4b). Despite the higher N Ac in summer, CN usually peaks in spring as a result of elevated N At (Fig. 4b). In comparison, the monthly average N At and CN exhibit some minor difference among the 3 years, while their seasonal trends remain the same (Fig. 4b). The mode-dependent seasonal trends indicate that the variations of N At , N Ac , and N LA are driven by different processes. Processes that may influence aerosol number concentrations in remote MBL are entrainment of the particles from the FT, SSA production, NPF inside the MBL, condensational growth (COND), coagulation (COAG), in-cloud scavenging of interstitial particles by droplets (INT), aqueousphase chemistry (AQ_CHEM), wet deposition, dry deposition, and advection. Among these processes, NPF within the MBL was shown to be infrequent in previous studies (Raes, 1995;Bates et al., 2000) and is neglected in the calculations of the long-term budget terms . Also, at remote marine sites like the ENA, the influence of advection is "averaged" out for long-term trends of particle concentrations. In addition, dry deposition is usually much slower compared to wet deposition for submicron particles, even after taking into account the time and spatial discontinuity of the wet deposition processes (see the discussion of in-cloud scavenging in Sect. 4.2) (Lewis and Schwartz, 2004;Henzing et al., 2006;Mohrmann et al., 2018). Thus,  it is neglected in further analysis. Wet deposition includes both in-cloud coalescence scavenging of activated droplets and therefore effectively CCN inside clouds (COALES) and the collection of aerosol particles by falling hydrometeors below clouds (i.e., washout). For aerosols between 10 nm and 1 µm, below-cloud washout is usually much less efficient than in-cloud coalescence scavenging (Garrett et al., 2006;Seinfeld and Pandis, 2016; and is neglected here. Earlier study suggests that the Ac mode in the MBL is formed through aqueous-phase chemistry inside cloud droplets (Hoppel et al., 1990). Therefore, we treat both Ac-and LA-mode particles as CCN, and At-mode particles as non-CCN (i.e., they remain as interstitial particles inside clouds). This treatment is also supported by the strong correlation between N Ac plus N LA and CCN concentration at 0.2 % ss, representative of marine low clouds Clarke and Kapustin, 2010) during all seasons (Fig. S6). As N Ac is usually 1 order of magnitude higher than N LA (Fig. 4), the CCN concentration at the ENA site is well represented by N Ac alone (Fig. S6). Therefore, the overall governing equation for each mode of MBL aerosol can be written as as depicted in Fig. 5 and discussed in detail below.

SSA
The change rate of MBL aerosol concentration due to SSA production flux, ∂ t N| SSA , can be expressed as (de Leeuw et al., 2011;) where 3.84×10 −6 WS 3.41 is the white cap fraction on the sea surface (Monahan et al., 1986) with WS in units of m s −1 , F SSA is the total SSA number production flux per white cap area in units of m −2 s −1 , H MBL is the MBL height in meters, and f SSA (lnD p ) is the lognormal number size distribution of SSA production flux curve. Thus, WS is the most important parameter in estimating total SSA contributions, while the detailed size distribution could differ with the f SSA (lnD p ) used (Gong, 2003;Lewis and Schwartz, 2004;Clarke et al., 2006;Grythe et al., 2014).

In-cloud coalescence scavenging
The rate of in-cloud coalescence scavenging of cloud droplets is given by (Wood, 2006; where E| X represents (∂ t N| X )/N , namely the percentage processing efficiency of process X. N d is cloud droplet number concentration which is assumed to be the same as CCN, or N Ac plus N LA (Sect. 3.4.2), and K is a constant of 2.25 m 2 kg −1 , while h H −1 MBL represents the in-cloud volume fraction of MBL aerosols (Mårtensson et al., 2010). Note that by setting the precipitation rate at cloud base, P CB , as 0 when there is no precipitation, the precipitation time fraction is already included in Eq. (3).

In-cloud scavenging of interstitial particles by activated droplets
Inside clouds, interstitial particles are scavenged when coagulating with cloud droplets. This process directly reduces At-mode particle number concentration, while also indirectly reducing CCN (i.e., Ac-and LA-mode) number concentration by removing particles that could otherwise grow and become CCN later (Pierce et al., 2015). The rate of scavenging scales with the probability that the particles are inside clouds, f cloud . Here, f cloud is defined as where p cloud is the probability that MBL cloud is encountered and is approximated by the in-cloud time fraction (Table 1), while h H −1 MBL is again indicative of the volume fraction of MBL aerosol particles inside the clouds (Mårtensson et al., 2010).
As At-mode particles are treated as non-CCN and remain as interstitial particles inside the clouds, the rate of the scavenging can be estimated by (Pierce et al., 2015) where N d is number concentration of cloud droplets assumed to be the sum of N Ac and N LA (Sect. 3.4.2), K int,d is the coagulation coefficient between D p,int and D p,d , where D p,int and D p,d represent the diameter of interstitial particles and cloud droplets, respectively. D p,d is assumed to be 10 µm (Pierce et al., 2015), while D p,int is assumed to be the corresponding wet diameter of D pg,At under a supersaturation (ss) of 0.12 % (Korolev and Mazin, 2003), where D pg,At is the geometric mean dry diameter of the At mode. The maximum supersaturation near the cloud base where CCN activation occurs is typically 0.2 % for marine low clouds Clarke and Kapustin, 2010;Leaitch et al., 2010). However, the supersaturation is usually lower above the cloud base where most of the interstitial scavenging occurs. Here, we assume the in-cloud ss of 0.12 % based on the work of Korolev and Mazin (2003). Assuming the At mode has a minimum D p of 23 nm , the D pg,At is estimated as 48 nm, and the corresponding wet particle diameter inside clouds, D p,int , is around 190 nm. Sensitivity of the interstitial scavenging rate to these parameters is discussed in Sect. 6.3.

Aqueous-phase chemistry
The aqueous-phase reaction (i.e., in-cloud production of sulfate) rate is positively related to the liquid water content (Seinfeld and Pandis, 2016;Meng and Seinfeld, 1994;Pandis et al., 1990;Cheng et al., 2016). As the liquid water content of cloud droplets is orders of magnitude higher than that of interstitial aerosols, only aqueous-phase reactions inside the cloud droplets are considered here (Pandis et al., 1990). As a result, the aqueous-phase reactions only promote the growth of CCN (i.e., Ac-and LA-mode particles). The influence of aqueous-phase reactions on the Aitken-mode particles is neglected until they reach CCN sizes through condensational growth (Hoppel et al., 1994;Pandis et al., 1990). Aqueous-phase reactions inside cloud droplets do not change total particle number concentration. On the other hand, they efficiently add mass to CCN and grow them into larger diameters when cloud droplets evaporate following the reactions. Therefore, the only influence of AQ_CHEM on number size distribution considered here is the growth of Acmode particles into the LA-mode size ranges. The magnitude of the influence depends on f cloud , liquid water content, precursor concentrations, and radiation which influences oxidant concentrations (Mårtensson et al., 2010).

Condensation growth
While condensation does not change the total particle number concentration, it grows the particles and therefore changes the number distribution among different modes (Seinfeld and Pandis, 2016). In this aspect, it functions similarly to aqueous-phase reactions, with the difference being that condensation acts on particles of all sizes, while aqueousphase reactions influence only CCN. The rate of a smaller mode A growing into a larger mode B through condensation can be estimated as where V A (in µm 3 ) is the volume difference between a particle with the minimum D p of mode B and a particle with the volume average of mode A, namely where the integrals are calculated from the binned aerosol size distribution using the binned simplification described in Pandis et al. (1994). J V (A) is the volume condensation rate of mode A in µm 3 m −3 s −1 , which can be estimated as (Seinfeld and Pandis, 2016) where K COND (A) is the condensation rate constant of mode A in s −1 , R is the gas constant of 8.314 J mol −1 , T is temperature in K, P is the atmospheric pressure being 1.013 × 10 5 Pa, ρ p is the aerosol density assumed to be 1 × 10 −12 g µm −3 , v i and v eq are the volume mixing ratio of condensate in the bulk gas-phase and at the aerosol surface, and m i is the molar mass of condensate. Here, we assume that the condensate is H 2 SO 4 , and thus m i = 98 g mol −1 , and v eq is 0 . Annual mean v i is assumed to be 1.0 ppt , while being 1.4, 1.3, 1.1, and 0.2 ppt in spring, summer, fall, and winter, respectively. This seasonal variation in v i is based on the monthly dimethyl sulfide (DMS) fluxes (assumed to be 7.0, 5.4, 2.9, and 1.0 µmol m −2 day −1 in spring, summer, fall, and winter, respectively) given in previous studies in the North Atlantic Ocean (Tarrasón et al., 1995), and the proposed dependence of H 2 SO 4 on DMS flux at the observed fluxes ranges Russell et al., 1994). Here, we assume that H 2 SO 4 is the dominant condensate. However, recent studies suggest that organics may play an important role in the growth of particles inside the MBL, and this is discussed later in Sect. 6.2. K COND (A) can be estimated by (Seinfeld and Pandis, 2016) where 10 −4 is the unit converter of µm cm −1 , D is the gas diffusivity of condensate in air equaling 0.1 cm 2 s −l , g f is the aerosol hygroscopic growth factor at ambient RH, D p is the dry aerosol diameter in µm, n(D p ) is the number size distribution of mode A in µm −1 cm −3 , D p1 and D p2 are the diameter boundaries of mode A and defined as the corresponding mode gap D p in Table 2 here, and Kn is the Knudsen number given by 2λ mfp (g f D p ) −1 , where λ mfp is the air mean free path. At the ENA site, observed ambient RH shows a modest diurnal variation of 75 % ± 10 %. Accordingly, the hygroscopic growth factor, g f , is assumed to be 1.3, based on the Hygroscopic Tandem Differential Mobility Analyzer measurements of At-and Ac-mode particles at the ENA site. The increased particle surface area due to hygroscopic growth leads to a factor of ∼ 1.7 increase in the estimated K COND compared with that under dry conditions. The term f (Kn, α) is the correction due to non-continuum effects (scaled by Kn) and imperfect surface accommodation (scaled by the mass accommodation coefficient α) estimated by the Fuchs-Sutugin approach as (Seinfeld and Pandis, 2016) f (Kn, α) = 0.75α(1 + Kn) where α is assumed to be 0.02 for H 2 SO 4 .

Coagulation
Unlike condensation, coagulation does not change the total mass concentration but reduces aerosol number concentrations. The intra-modal coagulation of particles in a smaller mode A (e.g., At mode) serves as both a source of particles in a larger mode B (corresponding rate denoted as J AA→B hereinafter) and a sink for particles of mode A. Given the typical aerosol size distribution observed at the ENA site, intramodal coagulation of mode-A particles is usually negligible when compared to inter-modal coagulation between mode A and another mode (B) with a different size range (corresponding rate denoted as J AB hereinafter) (Dal Maso et al., 2002). Therefore, we focus on the intra-modal coagulation as a source of a larger mode particles and inter-modal coagulation as a particle sink. The corresponding rates, J AA→B and J AB , are, respectively, estimated as where K 12 is the coagulation coefficient between two particles with diameters of g f D p1 and g f D p2 , respectively, and is calculated using the Fuchs form (Seinfeld and Pandis, 2016).
Similarly to the estimation of K COND , the growth factor under ambient RH, g f , is assumed to be 1.3. This increase in particle diameter results in a ∼ 20 % decrease in estimated K 12 . D p,min,A and D p,max,A are the boundary diameter of mode A (defined as the corresponding mode gap D p in Table 2 here), while D p,min,B and D p,max,B are defined similarly. D pc is defined by The coagulation loss rate of N A is thus while the coagulation production rate of N B is

Estimated rate of the potential key processes
The terms in the governing Eq. (1a-c) due to condensation, coagulation, scavenging of interstitial aerosol, and incloud coalescence scavenging of CCN are estimated using the equations described above and the size distribution parameters listed in Table 2. The values are listed in Table 3. The discussions of these estimates follow in Sects. 5 and 6.

Controlling processes of the larger accumulation mode
Potential processes that influence the LA-mode number concentration include Among these processes, SSA is expected to be the dominant source of N LA in the MBL, as suggested by the strong correlation between N LA and WS (Fig. 6b), a key parameter of SSA production flux (Sect. 4). Aqueous-phase reactions have been reported to produce "droplet-mode" particles in the LA-mode size range (Pandis et al., 1990;Meng and Seinfeld, 1994). However, if aqueous-phase reactions present a major source, we would expect the volume size distribution to exhibit a mode D p of 0.6-0.8 µm, corresponding to the size ranges that have the largest access to cloud water (Pandis et al., 1990;Seinfeld and Pandis, 2016). In Table 3. Estimated terms of the governing equations for three modes using size distribution parameters in Table 2   ; see sensitivity analysis in Sect. 6.2.
contrast, the volume size distribution indicates that the LA mode is actually the leading edge of a larger mode with peak D p in the supermicron range (Fig. 3b2). This is also supported by the strong correlation between V LA and PM c B sca (Fig. 6c). The PM c B sca is a surrogate for the supermicron mode (PM c , D p 1-10 µm) volume concentration (Sect. 2.2), while supermicron particles are dominated by SSA in remote MBL (Campuzano-Jost et al., 2003). Therefore, the strong correlation suggests that LA particles are also dominated by SSA, and the LA mode is essentially the sea spray aerosol coarse mode under vast majorities of the conditions. The contribution of aqueous-phase reactions to the LA-mode number concentration is likely minor and is neglected from the governing equation (Eq. 1c) in following analysis. Given the large sizes of LA-mode particles and that dust and biomass burning episodes are already excluded, we do not expect any significant FT sources. The lack of correlation between N LA and EBC mass concentration also suggests a low concentration of LA-mode particles in long-range transported continental pollution plumes. Here, we assume that the concentration of LA-mode particles in the FT is negligible when compared to that in the boundary layer. In such a case, the entrainment of FT air dilutes the MBL LA particles, serving as a sink rather than a source. At a typical entrainment velocity, ω e , of 3.5 mm s −1 (Mohrmann et al., 2018;Wood and Bretherton, 2004), the maximum dilution rate, −E LA | FT , equaling ω e H −1 MBL (Mohrmann et al., 2018), reaches ∼ 20 % per day. That is comparable to in-cloud coalescence scavenging, making the FT dilution an important sink of N LA . The sensitivity of the rate to entrainment velocity is discussed at the end of this section.
The terms of intra-modal coagulation (∂ t N LA | COAG ) and condensation (∂ t N LA | COND ) from the Ac mode are estimated as 0.02 and 0.6 cm −3 day −1 , respectively (Table 3). Both processes are too slow to exert significant influence on N LA during the typical aerosol lifetime of 7-10 days. The governing equation of N LA (Eq. 1c) can therefore be simplified into The seasonal variation of N LA is a result of the balance among the three processes (Fig. 6). Production flux of SSA is proportional to WS 3.41 H −1 MBL (Eq. 2), in-cloud coalescence scavenging efficiency is KP CB hH −1 MBL (Eq. 3), and FT dilution efficiency is estimated as ω e H −1 MBL . Among these three terms, the FT entrainment term (Fig. 6e) shows little seasonal variation. In comparison, both the in-cloud coalescence scavenging (Fig. 6d) and the SSA production (Fig. 6f) terms are lower in summer and higher in winter, with the SSA production exhibiting a stronger seasonal variation. The value of N LA under the quasi-steady state (i.e., when ∂ t N LA = 0) can be scaled using the three terms for each season. The scaled steady-state N LA (red markers in Fig. 6g) successfully produces the observed seasonal trend of N LA (box plots and black lines; Fig. 6g). Varying the value of the assumed entrainment velocity within the typical range of 2-5 mm s −1 does not affect the overall seasonal trend of the scaled N LA . 6 Controlling processes of the Aitken-accumulation mode

Contributions of SSA to Aitken and accumulation modes
Unlike N LA , N Ac is independent of the WS, and N At decreases with increasing WS (Fig. 7), indicating relatively minor contributions from SSA to At and Ac modes. The negative correlation between N At and WS may be due to the enhanced N LA with increasing WS (Fig. 6) and thus enhanced coagulation loss for Aitken-mode particles (see Sects. 4.2 and 6.3). In comparison, both N At and N Ac increase monotonically with EBC mass concentration (Fig. 7), suggesting the long-range transported anthropogenic aerosol is a major source of At-and Ac-mode particles in ENA. A semi-quantitative estimation of SSA contribution also supports the above conclusion. Assuming all LA-mode particles are from SSA, by combining N LA and an established size distribution of SSA production flux, one can estimate the upper limit of the SSA contribution to At and Ac modes (Fig. 8). For simplification, here we use number concentration of particles with D p in the range from 400 to 1000 nm, N 400 , to represent the observed SSA number concentration in the same D p range. SSA larger than ∼ 100 nm is CCN under ss of 0.1 % (Petters and Kreidenweis, 2007), while the measured CN has a cutoff diameter of roughly 10 nm. The contribution of SSA to CCN (0.1 %) and CN SSA can therefore be estimated by  4.3 f SSA, Ac (%) 21 ± 18 19 ± 15 12 ± 13 24 ± 18 31 ± 22 f a SSA, At (%) 10 ± 10 9 ± 8 8 ± 6 7 ± 9 16 ± 12 a Here, only an upper limit of k INT without considering the condensation growth is estimated. Correspondingly, the f SSA, At is also expected to be an upper limit.
where k CCN and k INT are factors that account for the size dependence of removal rate (see derivations in Sect. S3). The estimated k CCN is around 1, while k INT can vary from 1.7 to 4.4 as the removal efficiency is higher for CCN than non-CCN (Table 4, Fig. 9).
Here, we used four published f SSA (lnD p ) schemes (Gong, 2003;Lewis and Schwartz, 2004;Clarke et al., 2006;Grythe et al., 2014) to calculate the contribution of SSA to observed CCN (0.1 %) and CN (Fig. 8a). The initial calculation neglects the size dependence of the particle removal rate; therefore, the results represent lower limits on the contributions (Fig. 8b). This approach essentially assumes that the shape of the SSA size distribution in the MBL is the same as that of the SSA flux. Even for these lower limit estimates, CCN (0.1 %) SSA and CN SSA calculated using f SSA (lnD p ) from Gong et al. (2003) and Clark et al. (2006) exceed the observed total CCN (0.1 %) and CN for a substantial fraction of the data, suggesting that these two f SSA (lnD p ) functions result in overestimation of SSA contributions over the 10- 400 nm size range at the ENA site. This may be partially due to the parameter dependencies of sea surface temperatures, etc. (Gantt and Meskhidze, 2013;Gantt et al., 2015;Quinn et al., 2015), which are not considered here. The value of CCN (0.1 %) SSA and CN SSA are therefore estimated as the averages of predictions based on flux size distributions reported by Grythe et al. (2014) and Lewis and Schwartz (2004), with k CCN and k INT taken into consideration. The corresponding mean fractions of CCN (0.1 %) SSA and CN SSA in observed CCN (0.1 %) and CN are 24 % and 11 %, respectively. The estimated CCN (0.1 %) SSA fraction is consistent with a recent study that shows that the SSA contribution to CCN is smaller than 30 % globally (Quinn et al., 2017). In that study, the size distribution of SSA was derived by fitting the aerosol size distribution. If we follow the same approach (Quinn et al., 2017), the estimated SSA number concentration is actually N LA shown in this study, which represents 19 % of CCN (0.1 %).
Based on the estimated SSA contribution to CN and CCN (Eq. 8a, b), we can further estimate the SSA contribution to N Ac and N At , f Ac,SSA and f At,SSA , as and the corresponding annual mean f Ac,SSA and f At,SSA are 21 % and 10 %, respectively (Table 4).

Controlling processes of the accumulation mode
As shown in Sect. 6.1, the contribution of SSA to the Ac mode is likely substantial (annual average ∼ 21 %; Table 4). For the Ac mode, both intra-modal and inter-modal coagulations are much slower than in-cloud coalescence scavenging (Table 3) and can therefore be neglected from the governing equation of N Ac . On the other hand, condensational growth of Aitken-mode particles may represent a substantial source of the Ac mode. AQ_CHEM reduces N Ac by growing particles into LA size range. As discussed in Sect. 5, this process only makes a minor contribution to N LA . Given that N Ac is about 1 order of magnitude higher than N LA , the impact of AQ_CHEM on N Ac is therefore expected to be negligible. The governing equation of N Ac (Eq. 1b) can be simplified into The estimated values of ∂ t N Ac | COND and ∂ t N Ac | COALES are listed in Table 3. In-cloud coalescence scavenging is the only sink of Ac-mode particles among the four main processes, while the other three are sources. Under steady-state conditions, (∂ t N Ac = 0), we have Namely, The contribution due to the entrainment of FT air, ∂ t N Ac | FT , can be estimated as a residual using the f SSA,Ac value derived in the previous section (Table 4): The normalized rates of different processes are compared in Fig. 9a. The derived ∂ t N Ac | FT is stronger in winter-spring and lower in summer-fall, in general agreement with the seasonal trends of observed CO and EBC (Fig. 9b), consistent with the picture that anthropogenic emissions represent the main source of entrained FT Ac-mode particles. This agreement also suggests that the above analysis captures the major seasonal variation of the contribution of FT entrainment to Ac-mode particles.
Based on the first-order estimates shown in Fig. 9, we can see that, on an annual basis, entrainment from the FT represents the major source of N Ac , followed by condensational growth of Aitken-mode particles and SSA production. However, the relative importance of these three sources shows substantial seasonal variations. Contributions from SSA production are the lowest in summer (12 %) and the highest during winter (31 %; Table 4), a result of strong seasonal variation of surface wind speed. In contrast, condensation is negligible in winter due to the substantially lower DMS emissions and thus H 2 SO 4 concentrations (Sect. 4.2). In summer and fall, however, the contribution from condensational growth (60 % in summer and 42 % in fall) exceeded that from FT entrainment and became the dominate source of the Ac mode. This suggests that ocean ecosystems may have a substantial influence on Ac-mode particles and therefore MBL CCN populations in the ENA. The variation in relative importance despite similar E Ac | COND in spring to fall (Fig. 9a) is mainly due to the large seasonal variation in FT entrainment efficiencies. The estimated contribution from condensation is consistent with observations of individual aerosol particles in the western Atlantic (Sanchez et al., 2018) but is substantially higher than that simulated over the remote Southern Hemisphere oceans during summertime (Hannele et al., 2008). This difference is likely due to the much higher DMS sea surface concentration in the ENA (∼ 7.5 nM) than that in the southern oceans (∼ 2.5 nM) (Kettle et al., 1999), or due to the difference between observed and model-simulated aerosol size distributions, etc.
Major uncertainty in the above estimates comes from the concentration of condensates. First, there can be condensates other than H 2 SO 4 , such as organics. Common biogenic volatile organic compounds (BVOCs) such as isoprene and monoterpenes typically have a very low mixing ratio, and secondary organic aerosol (SOA) formation from these BVOCs is generally minor in the remote marine environment (Kavouras and Stephanou, 2002;Arnold et al., 2009;Gantt et al., 2009;Myriokefalitakis et al., 2010). However, recent studies suggest that photochemistry or heterogeneous oxidation at the sea surface microlayer may represent a substantial source of oxygenated gas-phase organic compounds (OVOCs), which potentially plays an important role in SOA formation and particle growth in the Arctic MBL Willis et al., 2017;Mungall et al., 2017). It is possible that the SOA formation from these OVOCs can contribute to the growth of Aitken-mode particles in the ENA as well. If so, the contribution to CCN by the growth of Aitken-mode particles would be even higher than the estimate here, which is based on condensation of H 2 SO 4 only. Second, several studies (Langley et al., 2010;Corbett and Fischbeck, 1997;Capaldo et al., 1999;Corbett et al., 2007;Wang et al., 2008;Johansson et al., 2017) have shown ship emissions represent a significant source of SO 2 in the MBL. In this study, the concentrations of SO 2 and H 2 SO 4 are estimated using DMS-SO 2 -H 2 SO 4 yields based on an observation-based parameterization Pandis et al., 1994). Therefore, H 2 SO 4 formed from ship-emitted SO 2 , and its contribution to condensational particle growth is implicitly included.

Controlling processes of the Aitken mode
The governing equation of N At is given by Following the same approach in Sect. 6.2, we have The contribution of SSA to the Aitken mode is even smaller than it is to the Ac mode, which is estimated to be no larger than 10 % (Table 4). As a result, the entrainment of FT At-mode particles represents the dominant source (Fig. 9a). ∂ t N At | FT is higher in spring-summer and lower in fallwinter, and such seasonal variation is somewhat different from that of the CO mixing ratio and EBC mass concentrations (Fig. 9b). These differences may be partially due to stronger new particle formation from biogenic precursors in the FT during spring and summer seasons (Sanchez et al., 2018). The strength of new particle formation is not correlated with CO or EBC concentrations, which are tracers for anthropogenic emissions. The contribution of NPF versus anthropogenic emissions to FT Aitken-mode particles cannot be quantitatively determined using data presented here alone and will be a subject of future study.
On an annual basis, inter-modal coagulation is the major (55 %) sink of N At (Fig. 9), followed by condensation growth (28 %) and interstitial scavenging (16 %). While it is less important when compared to inter-modal coagulation, interstitial scavenging is substantial and cannot be neglected. This is consistent with the finding of Pierce et al. (2015). The overall removal efficiency of N At (∼ 10 % day −1 ) is substantially lower than those of N Ac and N LA , which corresponds to a longer lifetime (∼ 10 days) for At-mode particles in the MBL. The N At removal efficiency is higher in summer and lower in winter, which is opposite that of Ac and LA modes. This is partially due to the less efficient removal of At-mode particles by coagulation and interstitial scavenging in winter, as a result of lower N Ac and therefore droplet number concentrations. In addition, the low DMS fluxes during winter (Sect. 4.2) lead to substantially weakened condensation growth of At-mode particles into Ac-mode size ranges, which also contribute to the lower overall removal efficiency in winter. Relative importance of these three removal processes is quite consistent in spring to fall, with contributions from coagulation, condensation, and interstitial scavenging being around 51 %, 33 %, and 16 %, respectively. In winter, condensation becomes a negligible (7 %) removal processes, while the contribution of coagulation dominates (71 %), with the remaining 22 % due to interstitial scavenging.
We note that there may be some uncertainties in the above estimates, especially the rate of interstitial scavenging, which depends on Aitken-mode size distribution and the supersaturation inside clouds, as well as the effective cloud droplet diameters. Based on the assumed baseline conditions (effective cloud droplet diameters of 10 µm, average dry interstitial aerosols of 48 nm, and average ss of 0.12 %; see Sect. 4.1), relative sensitivities of the K int,d are 10 %/µm, −5 %/ nm, and −8 %/ % with respect to changes in droplet diameter, dry in-terstitial aerosol diameter, and average ss, respectively. For average cloud droplet diameter at 15 µm, geometrical mean Aitken-mode diameter of 45 nm, and average ss of 0.1 %, a condition that is more favorable for intestinal scavenging, corresponding interstitial scavenging rate would increase by a factor of 1.8.
Given the low contribution of SSA to the At-mode particles, the governing equation for Aitken mode (Eq. 1b) can be simplified into

Conclusions
We examine the seasonal variations of aerosol properties, trace gas mixing ratios, and meteorological parameters measured at the ARM ENA site on Graciosa Island over a 3year period from 2015 to 2017. Aerosol size distributions from 70 nm to 1 µm typically consist of three modes: At (< 100 nm), Ac (100 to ∼ 300 nm), and LA (> 300 nm) modes. Observed CCN number concentrations are in general agreement with the sum of N Ac and N LA . The particle number concentration and mode diameter of the three modes exhibit different seasonal variations, suggesting that they are controlled by different processes. Sources of LA-mode particles are dominated by SSA. The major sinks of N LA are in-cloud coalescence scavenging and dilution by entrained FT air. N LA is higher in winter and lower in summer. The higher N LA during winter is attributed to strong SSA production flux due to high wind speed, which prevails over an increase in in-cloud coalescence scavenging. The seasonal variation of steady-state N LA is derived from scaling the rates of the major processes, and the result agrees well with the observation.
In comparison, SSA represents a relatively minor fraction of N Ac and N At , with estimated annual mean contributions being 21 % and no larger than 10 %, respectively. For N Ac , the other sources are entrained FT Ac-mode particles and condensational growth of Aitken-mode particles inside the MBL, while the major sink is in-cloud coalescence scavenging. The derived FT contribution to N Ac generally follows the seasonal trends of CO and EBC, namely higher in spring-winter and lower in summer, consistent with the picture that anthropogenic emissions represent the main source of entrained FT Ac-mode particles. While entrainment from the FT is the major source on the annual basis, the relative importance of the different sources varies strongly with the season. In summer and fall, condensation growth of the At mode may become the dominant source, contributing 60 % and 42 % of the Ac-mode particles in the MBL. In winter, SSA contributes to ∼ 31 % of the Ac mode, surpassing the contribution due to condensational growth. This is due to a For N At , entrainment from the FT is expected to be the dominant source, and coagulation represents the major sink. The derived FT contribution to N At is higher in springsummer and lower in fall-winter, possibly due to stronger NPF from biogenic precursors in the FT during spring and summer seasons (Tarrasón et al., 1995). The relative importance of NPF and long-range transported continental emissions to FT Aitken-and nucleation-mode particles, and the subsequent contribution to the MBL CCN population will be examined in future studies. On an annual basis, 52 %, 32 %, and 16 % of N At are removed by inter-modal coagulation, condensational growth, and interstitial scavenging, respectively. Relative importance of these three removal processes is quite similar from spring to fall. In winter, condensation becomes a negligible (7 %) removal process due to the low DMS fluxes, while contribution of coagulation increases to 71 %.
Based on the above results, the processes that control the concentrations of the different particle modes are summarized in Fig. 10. These results suggest particles entrained from the free troposphere represent the major source of CCN in the marine boundary layer. Some of the entrained particles directly contribute to the Ac-mode population in the MBL and are sufficiently large to serve as CCN. In addition, Aitken-mode particles in the free troposphere, which are attributed to NPF and long-range transported continental emissions, can grow and form CCN after their entrainment into the MBL. Our calculation suggests that this represents a significant source of MBL CCN all year, with the highest contribution of nearly 60 % during summer seasons. As the growth of Aitken-mode particles to CCN size is to a large degree the result of the condensation of sulfuric acid, a product of DMS oxidation, this suggests that ocean ecosystems may have a substantial influence on MBL CCN population in the ENA through emission of DMS.
Data availability. Data used in this study are publicly accessible at the permanent archive of data collected at the eastern North Atlantic site of the Atmospheric Radiation Measurement (ARM) user facility (available at https://www.archive.arm.gov/discovery; last access: 8 December 2018).
Author contributions. JW and GZ designed the research. GZ led the analyses. GZ and JW wrote the manuscript with major input from YW and MPJ, and further input from all other authors.
Competing interests. The authors declare that they have no conflict of interest.