New particle formation, growth and apparent shrinkage at a rural background site in western Saudi Arabia

Abstract

Atmospheric aerosols have significant effects on human health and the climate. A large fraction of these aerosols originates from secondary new particle formation (NPF), where atmospheric vapors form small particles that subsequently grow into larger sizes. In this study, we characterize NPF events observed at a rural background site of Hada Al Sham (21.802${}^{\circ }$ N, 39.729${}^{\circ }$ E), located in western Saudi Arabia, during the years 2013–2015. Our analysis shows that NPF events occur very frequently at the site, as 73 % of all the 454 classified days were NPF days. The high NPF frequency is likely explained by the typically prevailing conditions of clear skies and high solar radiation, in combination with sufficient amounts of precursor vapors for particle formation and growth. Several factors suggest that in Hada Al Sham these precursor vapors are related to the transport of anthropogenic emissions from the coastal urban and industrial areas. The median particle formation and growth rates for the NPF days were 8.7 cm${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$ ($J_{\mathrm{7}\mathrm{nm}}$) and 7.4 nm h${}^{-\mathrm{1}}$ (GR${}_{\mathrm{7}-\mathrm{12}\mathrm{nm}}$), respectively, both showing highest values during late summer. Interestingly, the formation and growth rates increase as a function of the condensation sink, likely reflecting the common anthropogenic sources of NPF precursor vapors and primary particles affecting the condensation sink. A total of 76 % of the NPF days showed an unusual progression, where the observed diameter of the newly formed particle mode started to decrease after the growth phase. In comparison to most long-term measurements, the NPF events in Hada Al Sham are exceptionally frequent and strong both in terms of formation and growth rates. In addition, the frequency of the decreasing mode diameter events is higher than anywhere else in the world.

1 Introduction

The effect of atmospheric aerosols on the Earth's radiative balance, via scattering, absorption and cloud interactions, is the single largest factor limiting our understanding of future and past climate changes (Stocker et al., 2013). In addition to the climate effects, aerosols are known to be detrimental to human health, with outdoor particulate pollution being the cause of more than 3 million premature deaths in the year 2010 (Lelieveld et al., 2015). These effects include the contribution of both primary and secondary aerosol particles. Primary particles are emitted into the atmosphere directly as particles, while secondary particles are formed from atmospheric vapors in new particle formation (NPF) events. Measurements of submicron particle number size distributions (PNSDs) have shown that NPF events are a global phenomenon (Kulmala et al., 2004), and they are estimated to produce around half (31 %–66 %) of the cloud condensation nuclei (CCN) in the lower troposphere (Yu and Luo, 2009; Merikanto et al., 2009; Gordon et al., 2017). Even in polluted regions, where the primary emissions are high, NPF events are found to occur frequently (Yu et al., 2017), and they are estimated to be a significant contributor to the particle number concentrations (Yue et al., 2011; Kulmala et al., 2016; Yao et al., 2018). Despite the importance of NPF, many aspects related to the initial formation and subsequent growth of secondary aerosol particles remain unknown. While sulfuric acid is widely regarded as the most important precursor for NPF, it is clear that other compounds are needed to explain particle formation and growth rates in ambient measurements, especially in the boundary layer (Kirkby et al., 2011; Ehrhart et al., 2016). Stabilizing bases, such as ammonia and dimethylamine, low-volatility oxidation products of VOCs (volatile organic compounds) and ions have been shown to enhance particle formation rates and bridge some of the gaps between theoretical evaluations, laboratory studies and ambient measurements (H. Yu et al., 2012; F. Yu et al., 2018; Almeida et al., 2013; Kürten et al., 2014, 2016, 2018; Zhang et al., 2004; Riccobono et al., 2014). The initial particle-forming compounds will also participate in growing the particles, but, in order to explain the observed growth rates, the presence of more abundant condensing (or otherwise particle mass forming, for example, via heterogeneous oligomer formation) species is required (Nieminen et al., 2010; Riccobono et al., 2012). These species are very likely to be low-volatility or semivolatile organic compounds, formed by the oxidation of either biogenic or anthropogenic VOCs (Smith et al., 2008; Tröstl et al., 2016; Dall'Osto et al., 2018). Overall, the mixture of compounds and the relative importance of different species participating in aerosol formation and growth are expected to vary depending on the ambient conditions; in some coastal environments NPF can be driven by iodine compounds (Sipilä et al., 2016), and, for example, in urban areas the uptake of nitrate can contribute significantly to aerosol mass (Li et al., 2018). Predicting all the occurring interactions in the atmosphere is impossible without observations from several different environments. While PNSD measurements have been conducted in a wide range of environments (Kerminen et al., 2018), especially continuous long-term measurements are still fairly uncommon and largely focused on Europe and the midlatitudes (Nieminen et al., 2018). Long-term measurements are needed for obtaining reliable estimates on the average properties and seasonal tendencies of atmospheric NPF. Such data are essential in evaluating the performance of global models, which currently have large uncertainties in simulating atmospheric NPF, as well as its contribution to CCN budgets and aerosol radiative effects in different environments (Pierce and Adams, 2009; Makkonen et al., 2012; Gordon et al., 2016; Semeniuk and Dastoor, 2018).

Recently, several NPF studies have pointed out an interesting phenomenon, where the average diameter of the particle mode formed in an NPF event begins to decrease after the growth phase. This is often referred to as aerosol shrinkage, but we will use the term DMD (decreasing mode diameter) event, since aerosol shrinkage quite directly implies a reduction in the size of individual particles, which is not necessarily the case. Such DMD events have been observed especially in subtropical regions (Yao et al., 2010; Backman et al., 2012; Cusack et al., 2013; Young et al., 2013; Zhang et al., 2016; Alonso-Blanco et al., 2017) but also in the temperate climate (Skrabalova, 2015; Salma et al., 2016). Typically, these DMD events are suggested to be caused by the evaporation of semivolatile compounds, due to changes in environmental conditions. However, the reduction in the mean diameter of a particle mode may also occur without evaporation, if smaller particles are transported to the measurement site (Kivekäs et al., 2016).

In this paper, we study the aerosol particle number size distribution dynamics at Hada Al Sham, Saudi Arabia, during February 2013–February 2015. The environment is quite unique due to high anthropogenic and low biogenic emissions, as well as a distinct segregation between the surroundings in different directions. To our knowledge, these are the first comprehensive long-term aerosol measurements conducted in the Arabian Peninsula. Two articles from the same measurement campaign are already published, describing the aerosol physical (Lihavainen et al., 2016) and optical properties (Lihavainen et al., 2017). The former article also contains analysis of particle number concentrations, which will not be presented in this study. This work focuses on identifying and characterizing NPF in the study region. We show that the NPF events in Hada al Sham are, in comparison to the locations analyzed in the exhaustive study by Nieminen et al. (2018), exceptionally frequent and intense, in terms of both particle formation and growth rates. We also make a detailed investigation of the typical diurnal cycle related to the NPF events and determine how the events are impacted by different environmental variables, including meteorological conditions and the concentrations and sources of primary aerosol particles and aerosol precursor compounds.

2 Measurements and methods

2.1 Measurement site and instrumentation

Hada Al Sham (21.802${}^{\circ }$ N, 39.729${}^{\circ }$ E; 254 m a.s.l.) is a small city in western Saudi Arabia (see Fig. 3). There are no major sources of anthropogenic emissions in the immediate vicinity, and the site can be described as a rural background site. Biogenic emissions are also presumably minor due to the arid desert climate and lack of vegetation. Sparsely inhabited desert-like areas cover the inland in the N–SE direction from the measurement site, but the coastal regions in the western sector are densely populated. Jeddah, the second largest city in Saudi Arabia, is located by the Red Sea $\sim \mathrm{60}$ km to the west from Hada Al Sham. In Jeddah there are several major emission sources, including power-generation plants, a seawater desalination plant, and an airport and a harbor, both of which experience heavy traffic due to the combined effect of economic growth and the vicinity of Mecca (also known as Makkah) located $\sim \mathrm{40}$ km to the south from Hada Al Sham. Other major emission sources in the region include a petrochemical refinery and steel industry in Rabigh ($\sim \mathrm{130}$ km NW; see Fig. 3) and a large oil refinery in Yanbu ($\sim \mathrm{290}$ km NW). The densely populated coastal region and the sparsely populated inland are separated by mountains running along the coast of the Red Sea. Hada Al Sham is located by the western slopes of these mountains and is thus topographically more connected to the coastal region.

The measurements were conducted at the Agricultural Research Station of King Abdulaziz University from February 2013 to February 2015. The instruments were placed inside a container, located on a sand field, with a distance of $\sim \mathrm{100}$ m to the nearest trees and other obstacles. The temperature was kept stable inside the container at $\sim \mathrm{25}$ ${}^{\circ }$C. The sample air inlets were located at a height of 4–4.5 m and the sample air flow rate was 16.7 L min${}^{-\mathrm{1}}$. A more thorough description of the measurement setup and the used instruments can be found in Lihavainen et al. (2016). In this study, we focus on the PNSD measurements in the mobility diameter range of 7–850 nm, obtained using a twin DMPS (differential mobility particle sizer; Wiedensohler et al., 2012), and the meteorological parameters (temperature, relative humidity, wind speed and wind direction), which were measured with a Vaisala WXT weather station. The twin DMPS used here consists of a short and a medium Hauke-type DMA (differential mobility analyzer, custom-made) and two CPCs (condensation particle counter, TSI 3772). To study the effects of mineral dust, we utilized the PNSDs in the aerodynamic diameter range of 0.5–10 $\mathrm{\mu }$m, measured with an APS (aerodynamic particle sizer, TSI 3321), and the mass of particles smaller than 10 $\mathrm{\mu }$m (PM${}_{\mathrm{10}}$), measured with a beta hybrid mass monitor (Thermo Scientific 5030). The inlet leading to the DMPS and the APS was equipped with a PM${}_{\mathrm{10}}$ filter and a twin diffusion dryer, which kept the sample air relative humidity mainly below 50 %. PNSD measurements were corrected for diffusional losses in the inlet line and they were at maximum around 15 % at the lower size limit of the DMPS. The PM${}_{\mathrm{10}}$ measurements were made from a separate inlet, equipped with a standard heater for sample drying. Since no gas-phase measurements were conducted during this campaign, we used data from the Ozone Monitoring Instrument (OMI) on board NASA's Aura satellite (Levelt et al., 2006) to estimate the ${\mathrm{SO}}_{\mathrm{2}}$ concentrations (OMI Level 2 ${\mathrm{SO}}_{\mathrm{2}}$ planetary boundary layer product; Li et al., 2013) at Hada Al Sham and its surroundings.

2.2 NPF event classification

NPF event classification was done for the measurement days based on the visual interpretation of PNSDs, as described by Dal Maso et al. (2005). Each day was classified as either (1) an NPF day, (2) a nonevent day or (3) undefined. In short, a day is classified as an NPF day if a new growing mode of particles appears in the nucleation mode ($d_{\mathrm{p}}$ < 25 nm) and the growing mode is observed to persist for several hours (see Fig. 1 for an example of an NPF day and Fig. A1 for time series of PNSDs from selected periods). If some, but not all, of the abovementioned criteria are fulfilled, the day is classified as undefined. In addition, a day is classified as undefined if the time development of the newly formed particle mode is highly erratic, or if the mode is not continuous due to significant breaks. Only the days when there is clearly no indication of NPF are classified as nonevent days.

As an addition to the traditional classification, each NPF day was further classified based on whether the mean diameter of the mode formed in the NPF event clearly starts to decrease after the growth phase (see Fig. 1) or not. These days are referred to as DMD days and non-DMD days, respectively.

2.2.1 Event times

To describe the progression of the NPF events, we determined the points in time when (1) NPF is first observed in the smallest size bins of the DMPS measurements, (2) NPF is no longer observed, (3) the mode diameter of the newly formed particles starts to decrease and (4) the mode formed in the NPF event is no longer distinguishable from the background aerosols due to decreased number concentrations or, for example, changes in air masses. These times were determined visually from the PNSDs, and they are referred to as NPF start, NPF end, DMD start and NPF event end, respectively (Fig. 1). NPF start times were only determined for the days when NPF was observed all the way from the smallest size bins, and the NPF event end times were only determined for the days when the event ended during the same day as NPF had started.

2.3 Determination of particle formation and growth rates

The formation and growth rates of newly formed particles are important quantities in describing NPF (Kulmala et al., 2012). They provide information on the strength of the NPF events and are closely connected to the atmospheric factors driving the process, such as the concentrations of condensable vapors. The growth rates of freshly formed particles have a critical role in their probability to survive into climate-relevant sizes, as particles that grow too slowly are removed by coagulation with larger preexisting particles (Kuang et al., 2009).

In this study, the particle growth rates were determined by following the time development of the geometric mean diameters (GMDs) obtained from lognormal fits to the PNSD at each measurement time. The fitting was done using an automatic algorithm developed by Hussein et al. (2005), which analyses the measured PNSD, fits 2–3 lognormal modes and returns the fitting parameters. In practice, the growth rates are determined by plotting the GMDs of fitted modes together with the PNSD and making a linear fit to the points selected to represent the mode formed in the NPF event (Fig. 1). Now, the diameter growth rate (GR) in the size range $\mathrm{\Delta }{d}_{\mathrm{m}}$ can be calculated simply as the slope of the fitted line:

1

Results of the NPF event classification as fractions of classified days, separately for each month and all of the classified days. The numbers above the bars show the number of classified days. Some of the months contain days from more than a year, while in some months (June, July, September and December) data were available only for a single year.

[Figure omitted. See PDF]

The monthly fractions in Fig. 2 show that the NPF frequency is high (55 %–85 %) throughout the year and that no clear seasonal pattern is observed. This implies that the NPF events are not limited by any factor with a strong seasonal variability. The most notable deviations from the average frequency are found in June, November and December, which all have a higher-than-average fraction of nonevent days. Although no seasonal cycle is seen in the total NPF frequency, the DMD events are more frequent during the summer (and autumn) months and less frequent during winter. The fraction of DMD events from all NPF days is highly variable, ranging from 33 % in November to 95 % in September.

The average NPF event frequency of 73 % in Hada Al Sham is among the highest event frequencies obtained from long-term measurements. Nieminen et al. (2018) compared PNSD measurements, consisting of at least 1 full year of data, conducted at 36 different sites around the world. They observed that NPF events are most frequent in South Africa, where the NPF frequencies from three different sites were 69 %, 75 % and 86 % (Hirsikko et al., 2012; Vakkari et al., 2011, 2015). Thus, the NPF fraction of 73 % obtained from the measurements presented here would take the third place in this global comparison. The high NPF event frequency is a direct indication of typically favorable conditions for new particle formation and growth. NPF event frequency has been shown to be affected by at least solar radiation, ${\mathrm{SO}}_{\mathrm{2}}$ concentration, vapor and particle sinks, and air mass origins (Nieminen et al., 2015; Kerminen et al., 2018). The effects of these factors are discussed in the following paragraphs in order to explain the observed NPF frequency in Hada Al Sham.

The connection between solar radiation and NPF events is related to atmospheric photochemistry: the production of sulfuric acid, which is regarded as the driving compound of atmospheric new particle formation in most environments (Weber et al., 1997; Birmili et al., 2003; Kuang et al., 2008; Paasonen et al., 2010; Yao et al., 2018), occurs mainly via oxidation of ${\mathrm{SO}}_{\mathrm{2}}$ by OH, and the concentration of OH is dependent on the intensity of solar radiation. Baranizadeh et al. (2014) and Dada et al. (2017) studied the effect of cloudiness on the NPF frequency at the SMEAR II station in Hyytiälä, Finland, and showed that the NPF frequency increased from $\sim \mathrm{0}$ % to over 50 % with decreasing cloudiness. During some months, the NPF frequency on clear-sky days reached over 70 % also in Hyytiälä, which clearly demonstrates the importance of photochemistry. The exact number of clear-sky days in Hada Al Sham could not be determined as no radiation measurements were made, but in general the radiation conditions in this area (Alnaser et al., 2004) are highly favorable for the occurrence of NPF.

Solar radiation alone is, of course, not sufficient to cause NPF if no precursor vapors for the production of nucleating and condensing compounds are available. ${\mathrm{SO}}_{\mathrm{2}}$, required for the production of sulfuric acid, is emitted especially from traffic and industries that process or consume fossil fuels. In Hada Al Sham, there are no significant sources of ${\mathrm{SO}}_{\mathrm{2}}$ but there are plenty of sources within a 100 km radius of the site, which includes the large urban and industrial areas of Jeddah and Mecca. Figure 3 displays the average ${\mathrm{SO}}_{\mathrm{2}}$ concentration, retrieved by satellite measurements, in the surroundings of Hada Al Sham on NPF days. Even though these concentrations represent the amount of ${\mathrm{SO}}_{\mathrm{2}}$ in the whole vertical column, the values should reflect the concentrations in the boundary layer since the most significant ${\mathrm{SO}}_{\mathrm{2}}$ sources are located at ground level (for reference, 1 DU corresponds to 11 ppbv when distributed into a 1 km boundary layer at $T=\mathrm{300}$ K and $P=\mathrm{1}$ atm). The figure shows that the ${\mathrm{SO}}_{\mathrm{2}}$ concentrations on the coastal region are very high, close to those obtained from the most-${\mathrm{SO}}_{\mathrm{2}}$-polluted regions in the world (Krotkov et al., 2016). The horizontal transport of the emissions seems to be restricted by the mountains, which causes the accumulation of the concentrations on a relatively narrow area and creates a distinct boundary in the ${\mathrm{SO}}_{\mathrm{2}}$ levels compared to the inland. The observed high ${\mathrm{SO}}_{\mathrm{2}}$ concentrations, in combination with the radiation conditions, suggest that the production rate of sulfuric acid is high at this site.

Figure 3

Average concentration of ${\mathrm{SO}}_{\mathrm{2}}$ from the OMI Level 2 ${\mathrm{SO}}_{\mathrm{2}}$ planetary boundary layer product (Li et al., 2013) in the surroundings of Hada Al Sham during NPF days. The concentrations are shown in Dobson units (1 DU $=\mathrm{2.69}\times {\mathrm{10}}^{\mathrm{16}}$ molecules cm${}^{-\mathrm{2}}$).

[Figure omitted. See PDF]

In Fig. 4, we compare the air mass history during the mornings of NPF event and nonevent days. The shown emission sensitivities are calculated from the 24 h retroplumes, initiated at the time when NPF is typically taking place (10:00 LT). The comparison shows significant differences in the air mass origins between these cases. On NPF days (Fig. 4a), the air masses observed in Hada Al Sham originate mainly from a narrow strip that extends along the coast and includes the regions of significant ${\mathrm{SO}}_{\mathrm{2}}$ emissions (see Fig. 3). The formation of such an air mass source region can be explained by the typically prevailing large-scale winds that blow along the Red Sea, due to tunneling caused by the steep shores. Over the coastal cities, anthropogenic emissions are introduced into the air masses and then transported to Hada Al Sham with the sea breeze that usually starts to develop already early in the morning. On nonevent days (Fig. 4b) there are no signs of the sea breeze or the typically prevailing large-scale wind. Instead, the source regions initially point towards the east, after which they spread to cover larger areas in the inland. The inland is sparsely inhabited and the sources of anthropogenic emissions are few, which can also be seen as lower average ${\mathrm{SO}}_{\mathrm{2}}$ concentrations in the surroundings of Hada Al Sham on nonevent days than on NPF days (Fig. A2). The fact that all of the nonevent days are observed when the air masses are coming from the inland suggests that the NPF events observed in Hada Al Sham might be related to the emissions from the coastal anthropogenic activities. In addition, the influence of marine air could be beneficial for NPF occurrence, e.g., due to lower condensation sink or some specific precursors.

Figure 4

Average 24 h emission sensitivity for air masses arriving at Hada Al Sham at 10:00 LT for (a) NPF days and (b) nonevent days. The arrival time of the air masses is chosen to represent the time when new particle formation is typically taking place on NPF days (see Fig. 5).

[Figure omitted. See PDF]

The markedly different wind conditions between the NPF and nonevent days can also be seen from the local measurements at Hada Al Sham (Fig. 5a and b). On NPF days, the weak nocturnal easterly wind (land breeze) turns westerly and its speed starts to increase between 08:00 and 10:00 due to the development of the sea breeze. This shifts the air mass source regions to the coastal areas, as seen in Fig. 4a. On nonevent days, such change in the wind direction is not seen. During the night and early morning, the wind is easterly, similar to the NPF days, but, on nonevent days, the easterly wind is significantly stronger. This inhibits the development of the sea breeze circulation and the westerly wind associated with it.

Figure 5

Diurnal variation of (a) wind speed, (b) wind direction, (c) condensation sink, and (d) mass of particles smaller than 10 $\mathrm{\mu }$m in Hada Al Sham during NPF days (red lines) and nonevent days (blue lines). The solid lines show the median values and the shaded areas represent the 25th–75th percentile range.

[Figure omitted. See PDF]

The strong easterly winds on nonevent days seem to resuspend dust from the inland desert, which can be seen as the simultaneous increase in PM${}_{\mathrm{10}}$ with the wind speed (Fig. 5d). Both the wind speed and the PM${}_{\mathrm{10}}$ obtain their largest values around the same time when NPF typically starts (see Fig. 6). It is therefore possible that, in addition to the fewer emission sources in the inland, the NPF is inhibited by the windblown dust, which can reduce the concentrations of the clustering and condensing vapors by limiting solar radiation and by acting as a sink for reactive gases and oxidants (Hanisch and Crowley, 2003; Usher et al., 2003). The condensation sink does not, however, seem to be an inhibiting factor for NPF occurrence, as significantly larger CS values are observed on NPF days than on nonevent days (Fig. 5c). Here the presented CS values are calculated including only the particles measured by the DMPS ($d_{\mathrm{p}}$ < 850 nm), but including the larger particles measured by the APS ($d_{\mathrm{p}}$ up to 10 $\mathrm{\mu }$m) did not have a significant effect on these results (increase in CS on nonevent days was around 10 %). Furthermore, recent studies have shown indications that mineral dust, mixed with anthropogenic emissions, could actually enhance NPF due to heterogeneous production of sulfates and hydroxyl radicals (Nie et al., 2014; Xie et al., 2015). In Hada Al Sham, the concurrence of high PM${}_{\mathrm{10}}$ values and nonevent days would then mainly highlight the lack of anthropogenic emissions in the inland, although we cannot quantify the effect of PM${}_{\mathrm{10}}$ on radiation.

Higher CS values during NPF days (Fig. 5c) are also reported from a high-altitude site in the Swiss Alps (Boulon et al., 2010). Here this connection is speculated to stem from the coupled appearance of NPF precursor vapors and CS due to their common lower-altitude sources. Similar situation could apply to Hada Al Sham if both the CS and the NPF precursor vapors originate mainly from the same sources. The high CS during the calm nights and early mornings of NPF days could then suggest enhanced accumulation of precursor vapors, thus facilitating the occurrence of the NPF process after sunrise.

3.2 NPF event progression and characteristics

Figure 6 displays the frequency histograms of the points in time describing the progression of the NPF events (see Sect. 2.2.1), together with the diurnal variation of meteorological parameters and CS. In Fig. 7, the NPF progression times are plotted for a year-long period from June 2013 to June 2014 to show their seasonal variation. NPF events are typically observed to start slightly before 09:00 (Fig. 6a). On a seasonal scale, the starting times change according to the changes in the time of sunrise (Fig. 7). This observation highlights the importance of photochemistry for the NPF, especially since none of the NPF events start before sunrise. There is, however, quite a significant day-to-day variation in the starting times, which is not explained by differences in the times of sunrise. Some of this variation can be attributed to differing growth rates between the NPF events. This is because the NPF starting times are defined here as the times when new particles are observed at the size of 7 nm, even though the formation of new particles actually starts from the molecular scale. Therefore, the time that it takes for the small particles/molecular clusters of the size $\sim \mathrm{1}$–2 nm (Kulmala et al., 2013) to grow and reach the lower limit of the DMPS affects the starting times shown here. In some cases, the later NPF starting times are caused by a delayed shift in the wind direction from the inland to the coastal side of the measurement site (Fig. A3). This is in agreement with the interpretation that the NPF events observed in Hada Al Sham are related to the transport of emissions from the coastal regions. The connection between the shift in the wind direction and the onset of NPF can also be seen from the average values, as the typical onset time of the NPF events coincides with the onset of the westerly sea breeze (Fig. 6a and c). In addition to the radiation conditions and the wind direction, the starting time seems to be connected to a drop in the CS, which is likely caused by the increasing ABL height (Fig. 6a and b).

Figure 7

Seasonal variation of the different phases of the NPF events observed in Hada Al Sham. The colored lines show the 20-point moving average for each of the different phases. The solid black lines show the times of sunrise and sunset, while the dashed black line shows the time of maximum solar radiation calculated based on the latitude of the measurement site.

[Figure omitted. See PDF]

A vast majority ($\sim \mathrm{85}$ %) of the NPF events observed in Hada Al Sham end during the same day they started, approximately 3 h after DMD start or 9 h after NPF start (Fig. 7). Here the NPF event end times were defined as the points in time when the number concentrations, associated with the particle mode formed in the NPF event, drop significantly, making the mode indistinguishable from the background aerosols (see Fig. 1). Based on the observations presented earlier in this work, explaining the ending of the NPF events would be rather straightforward if they were related to the wind turning easterly at around 23:00 (Fig. 6c). However, the NPF events are typically observed to end hours before this (Fig. 6a), while the winds are still westerly. This indicates that also in the westerly direction an area where NPF is not occurring is reached. When the air masses that resided in this area during the active NPF hours (NPF start–NPF end) are transported to the measurement site, the particle mode related to NPF is no longer seen. An order-of-magnitude estimate for the westerly extent of the NPF area can be obtained by multiplying the average duration of the NPF events by the average wind speed: $\mathrm{10}\mathrm{h}\times \mathrm{10}\mathrm{km}{\mathrm{h}}^{-\mathrm{1}}=\mathrm{100}\mathrm{km}$. This is comparable with the distance from Hada Al Sham to the coast of the Red Sea where the concentration of ${\mathrm{SO}}_{\mathrm{2}}$ drops rapidly (Fig. 3). Therefore, it seems that no NPF is happening outside the region of strong contribution from anthropogenic emissions. It is reasonable to assume that there is no discrete boundary between the regions where NPF is and is not occurring. Instead, when moving away from the region of high emissions, particle formation and growth rates can be expected to decrease gradually due to the decrease in the concentrations of participating vapors. This provides another possible explanation, in addition to particle evaporation, for the observed DMD events: when the particles formed in less favorable conditions are transported to the measurement site, it is possible that they would have grown less than the previously observed particles. Continuous observations of particles that have grown less and less can produce a DMD event without particle evaporation or other shrinkage (Kivekäs et al., 2016). The cause of the DMD events observed in Hada Al Sham is further investigated in a future study.

3.3 Particle formation and growth rates and their dependence on environmental conditions

Figure 8 shows the seasonal variation of the particle formation rates ($J_{\mathrm{7}\mathrm{nm}}$) and growth rates (GR${}_{\mathrm{7}-\mathrm{12}\mathrm{nm}}$) determined for the NPF events observed in Hada Al Sham. The formation rates (Fig. 8a) vary mostly between 1 and 50 cm${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$, having an annual median (determined as the median of all observations) of 8.7 cm${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$. The variation in the growth rates (Fig. 8b) is slightly smaller with the values ranging mainly from 2 to 20 nm h${}^{-\mathrm{1}}$ and having a median of 7.4 nm h${}^{-\mathrm{1}}$. The variation of these values is similar to the observations from several locations around the world (Kulmala et al., 2004), but the median values are very much in the high end of observations (Nieminen et al., 2018; Kerminen et al., 2018). Nieminen et al. (2018) reported, in their global study, the highest annual average particle formation rate from Beijing ($J_{\mathrm{10}\mathrm{nm}}\sim \mathrm{7}$ cm${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$) and found an overall increasing trend of formation rates with increasing degree of anthropogenic influence.

Figure 9

Particle formation rate ($J_{\mathrm{7}\mathrm{nm}}$) as a function of (a) temperature, (b) relative humidity, (c) wind speed, (d) atmospheric boundary layer height, (e) PM${}_{\mathrm{10}}$ and (f) the condensation sink. The black lines show the least squares fits to the log-linear (a, b, c, d) or log–log (e and f) data, and the $r$ and $p$ values denote the Pearson correlation coefficients and their significance levels, respectively. The values on the horizontal axis are event-time averages (from NPF start to NPF end; see Fig. 1). In the case of the CS, the averaging was done for 1 h before the NPF start to make sure that the correlation is not influenced by the particles formed in the NPF event itself. The data are colored according to the season, so that summertime is represented by yellow and wintertime by blue colors.

[Figure omitted. See PDF]

Figure 9b shows that the formation rate increases with increasing relative humidity (RH). This is the expected relationship between these two variables (Almeida et al., 2013; Duplissy et al., 2016; Kürten et al., 2016), as water vapor is known to participate in the cluster formation with sulfuric acid. However, in ambient measurements such a correlation could be caused by processes that are not necessarily related to the RH effect itself. For example, here the higher RH could be related to the coastal origin of the air masses and simultaneously to higher anthropogenic emissions from the coastal sector. To examine this, the correlation was calculated separately for winds coming from the S–W sector, where the ${\mathrm{SO}}_{\mathrm{2}}$ distribution seems most uniform (Fig. 3), but a similar relationship to the case with all data was found. Furthermore, the RH relationship does not seem to be related to the seasonal variation, as both high and low RH values are observed throughout the year.

Higher formation rates seem to be favored by low wind speed and low ABL height (Fig. 9c and d). During low ABL conditions, the near-surface anthropogenic emissions are distributed into a smaller volume, which could then lead to higher vapor concentrations and particle formation rates. Analogously, the accumulation of emissions, per unit volume of air, from a spatially limited emission source area is increased during low wind speed conditions. Interestingly, the lowest event-time ABL heights are observed during the summer months (Fig. 9d), meaning that, in addition to the summer NPF events occurring earlier in the absolute sense (see Fig. 7), they also take place earlier with respect to the boundary layer development. This observation could be caused by the higher emissions during summer (see Fig. 8a), since if the concentrations of NPF precursors are higher, the onset of NPF events is likely to be more sensitive to an increase in solar radiation. We note that the RH dependence might also arise partly from higher RH in low ABL and wind speed conditions. However, since the correlation of RH and either ABL height or wind speed is weaker than that between RH and the formation rate, we expect this connection to be of secondary importance.

In the end of Sect. 3.1, we discussed briefly the possible interactions between NPF and mineral dust, which is likely the major component of PM${}_{\mathrm{10}}$ in Hada Al Sham. We stated that mineral dust can possibly either weaken (increasing sink and decreased solar radiation) or enhance (increasing production of hydroxyl radicals and sulfates via heterogeneous reactions) the NPF events. The enhancing effect has been observed specifically in situations where mineral dust is mixed with anthropogenic pollution (Nie et al., 2014), which could correspond to the situation in Hada Al Sham. Despite this, no clear correlation was found between PM${}_{\mathrm{10}}$ and the formation (or growth) rates (Figs. 9e and A4e). It should be noted that in the case presented by Nie et al. (2014) the timescale of the process is significantly longer, as both the sources of the dust and the anthropogenic emissions are located further away from the measurement site. This allows for a longer interaction time between the dust and the emissions, which might be crucial for the enhancing effect to occur.

Out of the included variables, the strongest correlation is found between particle formation rates and CS (Fig. 9f). This positive correlation is quite interesting, since the concentrations of vapors participating in NPF are expected to decrease with increasing CS due to their faster loss rate. However, this is generally valid only if the sources of the CS and the condensing vapors are independent from one another. Here, this is likely not the case, but instead the increasing CS presumably represents increasing contribution from the (primary) anthropogenic emissions and is therefore simultaneously linked to higher concentrations of NPF precursor vapors. This is supported by the observation that both the CS and ${\mathrm{SO}}_{\mathrm{2}}$ concentration are higher during the NPF days than nonevent days (Figs. 5c and A2).

Figure A4 shows the particle growth rates as a function of the same variables as the particle formation rates in Fig. 9. Overall, the correlations are qualitatively similar but weaker in the case of GRs. Similarly to the formation rates, the strongest correlation is found with the CS ($r=\mathrm{0.37}$, $p=\mathrm{2.2}\times {\mathrm{10}}^{-\mathrm{5}}$). The most notable difference is that, in the case of GRs, a weak positive correlation with temperature ($r=\mathrm{0.19}$, $p=\mathrm{0.025}$) is observed.

4 Summary and conclusions

The analysis of the aerosol number size distribution measurements showed that NPF events are a highly frequent phenomenon in Hada Al Sham, with the fraction of NPF days accounting for 73 % of all the classified days. The high NPF frequency is likely explained by the high production of NPF precursor vapors, especially sulfuric acid, in the transported emission plumes from the coastal cities and industrial areas during the typically prevailing cloud-free and high-solar-intensity conditions. The fraction of nonevent days was only 4 %, and these days were shown to be linked to strong easterly winds that block the development of the sea breeze, which typically brings the polluted air masses to Hada Al Sham.

Most of the NPF events in Hada Al Sham displayed an unusual progression, where the diameter of the particle mode related to the NPF event started to decrease after the growth phase. Similar DMD events have been observed in other measurement sites as well, but in Hada Al Sham the frequency of these events was found to be exceptionally high (76 % of all NPF days). The DMD events were more frequent during the summer, and the average onset time of the DMD events was during the afternoon, approximately 6 h after NPF start.

The median particle formation and growth rates associated with the NPF events were 8.7 cm${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$ ($J_{\mathrm{7}\mathrm{nm}}$) and 7.4 nm h${}^{-\mathrm{1}}$ (GR${}_{\mathrm{7}-\mathrm{12}\mathrm{nm}}$), respectively. These values correspond to those typically obtained from polluted urban measurement sites. Both the formation and growth rates showed the highest values during summer and autumn months, presumably due to the increased emissions from energy production and the effect of stronger solar radiation on the rate of photochemical reactions. The formation rates were found to obtain higher values in calm conditions where both the wind speed and the ABL height were low and the relative humidity was high. Under such circumstances, anthropogenic emissions are likely to spread and accumulate throughout the coastal zone, including Hada Al Sham. Both the formation and growth rates obtained higher values in conditions of high CS, which is likely associated with the common anthropogenic sources of NPF precursor vapors and the large primary particles that control the CS.

Overall, the findings of this study highlight the importance of anthropogenic emissions and photochemistry for NPF. Due to the transportation of emissions from urban and industrial areas, NPF events were found to be very frequent in Hada Al Sham, located tens of kilometers away from the major sources. The frequency and strength of NPF observed here implies that NPF events might contribute significantly to the budget of both ultrafine and CCN particles, making their health and climate effects relevant topics for further studies in this region. The local conditions at Hada Al Sham, with high levels of regional anthropogenic emissions but presumably low concentrations of biogenic vapors, also allow us to research anthropogenic NPF events in detail. However, further experiments with a broader spectrum of instruments are required for determining the vapors responsible for new particle formation and growth, as well as the underlying reasons for the occurrence of the DMD events.

Data availability

Data used in this study are available from the corresponding author upon request ([email protected]).

Appendix A

Figure A1

Time series of particle number size distributions during February 2013, April 2014, July 2013 and December 2013 illustrating the high frequency and typical characteristics of NPF, as well as periods of nonevent and undefined days, e.g., during 23–27 December 2013. The shown time series were selected so that they would cover different seasons with continuous data availability.

[Figure omitted. See PDF]

Figure A2

The average ${\mathrm{SO}}_{\mathrm{2}}$ concentration within the radius of 50 km from Hada Al Sham during NPF days and nonevent days from the OMI Level 2 ${\mathrm{SO}}_{\mathrm{2}}$ planetary boundary layer product (Li et al., 2013). The concentrations are shown in Dobson units (1 DU $=\mathrm{2.69}\times {\mathrm{10}}^{\mathrm{16}}$ molecules cm${}^{-\mathrm{2}}$).

[Figure omitted. See PDF]

Figure A3

An example illustrating a case where the late starting time of an NPF event is related to a delayed shift in the wind direction in Hada Al Sham, 5 February 2013.

[Figure omitted. See PDF]

Figure A4

Particle growth rate (GR${}_{\mathrm{7}-\mathrm{12}\mathrm{nm}}$) as a function of (a) temperature, (b) relative humidity, (c) wind speed, (d) atmospheric boundary layer height, (e) PM${}_{\mathrm{10}}$ and (f) the condensation sink. The black lines show the least squares fits to the log-linear (a, b, c, d) or log–log (e and f) data, and the $r$ and $p$ values denote the Pearson correlation coefficients and their significance levels, respectively. The values on the horizontal axis are event-time averages (from NPF start to NPF end; see Fig. 1). In the case of the CS, the averaging was done for 1 h before the NPF start to make sure that the correlation is not influenced by the particles formed in the NPF event itself. The data are colored according to the season, so that summertime is represented by yellow and wintertime by blue colors.

[Figure omitted. See PDF]

Author contributions

MIK, HAJ, MAA, HL, APH and TH coordinated the measurement program and carried it out with KN, LD, ASA, IIS and FMA. ASA, IIS and FMA also provided resources and data curation. VV provided the essential means for calculating and analyzing the air mass trajectories. AMS produced the figures related to the OMI ${\mathrm{SO}}_{\mathrm{2}}$ data. SH and PP performed the data analysis, while VV, VMK, JK and MK contributed to the interpretation of the results. SH wrote the manuscript with contributions from all coauthors.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We would like to thank the editor, Fangqun Yu, and the two anonymous referees for their insightful comments that helped us improve the manuscript.

Financial support

This research was funded by the Deanship of Scientific Research, King Abdulaziz University (grant no. I-122-30); the Academy of Finland (grant nos. 307331 and 316114); and the European Commission (grant no. 742206).Open access funding provided by Helsinki University Library.

Review statement

This paper was edited by Fangqun Yu and reviewed by two anonymous referees.