Ismail Anil, Omar Alagha , Nawaf I. Blaisi, Iehab Abdelilah Mohamed, Mohammad Hisham Barghouthi, Mohammad Saood Manzar

Department of Environmental Engineering, College of Engineering, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia

Received: May 3, 2019
Revised: September 16, 2019
Accepted: October 23, 2019
Download Citation: ||  

  • Download: PDF

Cite this article:

Anil, I., Alagha, O., Blaisi, N.I., Mohamed, I.A., Barghouthi, M.H. and Manzar, M.S. (2019). Source Identification of Episodic Rain Pollutants by New Approach: Combining Satellite Observations and Backward Air Mass Trajectories. Aerosol Air Qual. Res. 19: 2827-2843.


  • Long-range transport from Kuwait, Iraq, Qatar contributed to chemical composition.
  • Source-identification by combining satellite-observations & back-trajectories.
  • Chemical composition was dominated by Ca2+, SO42–, Cl, HCO3, and Na+.
  • Ca2+ was responsible for 80% of total neutralization in rainwater, mean pH: 7.34.
  • SO42–, NO3, and NH4+ mainly originated from intense local anthropogenic activities.


In this study, various source apportionment methods, viz., enrichment factor, relative source contribution, and factor analysis in addition to cluster analysis coupled with satellite observations, were used to investigate the origin of rainwater pollutants in an arid urban area of Saudi Arabia. The rainwater samples were collected by an automatic wet-only sequential rain sampler and analyzed for their pH, electrical conductivity (EC), and major ions (Cl, NO3, SO42–, HCO3, Na+, K+, NH4+, Ca2+, and Mg2+). The results revealed that the chemical composition was dominated by Ca2+, SO42–, Cl, HCO3, and Na+, which represented 80% of the total ionic equivalent concentration. Ca2+, the most abundant ion, exhibited a mean concentration of 480 µeq L–1 and accounted for 30% of the total ionic equivalent concentration. The study area received a total ionic wet deposition flux (FWD) of 4.07 tons km–2 y–1, and SO42– and NO3 contributed 76% and 24%, respectively, of the rainwater acidity. Ca2+ was responsible for 80% of the rainwater’s neutralization. The source apportionment indicated that intense local human activity and medium- to long-range transport from Kuwait, Iraq, Qatar, Riyadh, and the Western Province produced the majority of the anthropogenic components (SO42–, NO3, and NH4+), which represented 47% of the total ionic FWD. However, moderate to high particulate matter loads over the southern part of Saudi Arabia and long-range transport from Iraq and Kuwait contributed the crustal components (Ca2+, HCO3, K+, and Mg2+), which accounted for 36% of the total ionic FWD, whereas the Arabian Gulf was the primary source of the marine components (Na+ and Cl). We found combining satellite observations and cluster analysis of backward air mass trajectories to be an effective new approach for assessing the source-receptor relationships for atmospheric air pollutants.

Keywords: Source apportionment; Backward trajectory; Cluster analysis; Satellite observation; Rainwater chemicalcomposition; Wet deposition flux.


The rapid development in urban areas could lead to serious air pollution concerns that affect millions of people around the world. Huge amounts of anthropogenic pollutants enter the atmosphere every instant via fossil fuel burning or from major power generation activities and industrial processes, which further lead to the development of unwanted gaseous pollutants such as oxides of sulfur and nitrogen in the atmosphere (Li et al., 2016; Wu et al., 2018). Thereafter, these gases are converted to strong acids, namely sulfuric and nitric acid, by gas-to-particle conversion, and they ultimately contribute to the formation of acid rain (Al-Momani et al., 2000; Başak and Alagha, 2004; Gaga and Ari, 2011; Burns et al., 2016).

Air pollutants are removed from the atmosphere by wet deposition and dry deposition processes (Likens, 1987; Budhavant et al., 2009; Tsai et al., 2014). Wet deposition is the main route and a powerful mechanism for the scavenging of the two phases, namely the gaseous and the particulate phase (Facchini Cerqueira et al., 2014). The deposition could be very acidic which may harm habitats and ecology systems by soil acidification and excess level of nitrogen leading to eutrophication (Driscoll, 2001; Bouwman et al., 2002; EMEP, 2013).

The rainwater chemical composition can reveal significant information on the degree of pollution and could help to assess the atmospheric deposition and transport of acidifying, neutralizing, and eutrophying species, heavy metals, and aerosols (Alagha and Tuncel, 2003; Tiwari et al., 2016; Anil et al., 2017). Furthermore, statistical evaluation of rainwater chemical composition data can provide insights into the necessary control and management practices to reduce and mitigate air pollution (Alagha and Tuncel, 2003; Başak and Alagha, 2004; Akkoyunlu et al., 2011). The study of rainwater chemical composition has been proven to be a powerful and reliable methodology to assess atmospheric pollution and its environmental impacts (Alagha and Tuncel, 2003; Al-Khashman, 2009; Akkoyunlu et al., 2011; Romero Orué et al., 2019; Vlastos et al., 2019).

Rain scarcity is a challenge to agricultural activities in the arid and semiarid regions (Whelpdale et al., 1997; DeNicola et al., 2015). Saudi Arabia is recognized by its arid and harsh weather with very scarce rainfall throughout the year (Al-Refeai and Al-Ghamdy, 1994; Mashat and Abdel Basset, 2011). There are very few studies in arid regions such as Saudi Arabia regarding rainwater chemistry and its association with the environment or atmosphere (Ahmed et al., 1990; Schemenauer and Cereceda, 1992; Alabdula’aly and Khan, 2000; Wedyan, 2013; Michelsen et al., 2015). Previous research points that the nature of the water collected from rain or wet deposition was mostly alkaline and partially acidic. Some of the studies in western regions of Saudi Arabia pinpointed the presence of terrestrial fungi in most of the rain samples (Nasser, 2005). In addition to these, other researchers investigated the occurrence and sources of polycyclic aromatic hydrocarbons in rainwater (Al-Daghri et al., 2013; Al-Saleh et al., 2013; Al-Daghri et al., 2014).

Air mass back trajectories could provide insights into the source regions contributing to concentrations of monitored air pollutants at a receptor site if the emission inventories of the potential source regions are well known (Kant et al., 2019). In practical application, presence or access to those kinds of datasets may not be straightforward. In that case, an alternative approach combining satellite observations and backward air mass trajectories could procure a better understanding of the source-receptor relationship for atmospheric air pollutants. Notaro et al. (2013) investigated the trajectory analysis of Saudi Arabian dust storms by utilizing weather observations and remotely sensed aerosol optical depth (AOD). The current study deals with remotely sensed air pollutants data such as AOD, NO2, SO2, and sea salt. The new addition in the current study is that it used the real-time measurements of pollutants from collected rainwater along with remotely sensed satellite data. The main objective of this study is to identify the sources affecting the elevated rainwater pollutants using statistical models including enrichment factor, relative source contribution, factor analysis, and cluster analysis. The approach of combining satellite observations and cluster analysis of backward air mass trajectories of rainy days was used as a tool in investigating the source-receptor relationship for rainwater pollutants.


Study Area and Sampling Methodology

The study area Dammam is the largest city in the Eastern Province Region of Saudi Arabia with an area of 800 km2 and a population of about 1.3 million inhabitants. It is also considered as the major center for oil, commerce, and the petrochemical industry. The urbanized part of the Dammam metropolitan area is approximately 912 km2 and the land use is distributed as follows: 36% residential and commercial, 32% zoning area, 19% industrial, and 13% seaports and airports. The study area receives 7–10 severe dust storms per year. During these episodic events, atmosphere of the study area is highly enriched with mineral particles when strong winds lift these particles from the regional soil which is reported to be composed of alluvium of gravel, sand, silt, clay deposits, and limestone subsoil (Alsaaran, 2008; El-Mubarak et al., 2014). The temperature is rarely below 7°C or above 46°C and usually ranges between 11°C and 44°C with an annual average of 28°C. The rainy season is typically between mid-October and mid-May with an average annual precipitation total of 100 ± 8 mm (Mashat and Abdel Basset, 2011). Accordingly, the precipitation regime of the region is categorized as “arid”. Wind systems from northern directions dominate over the area with an average ground-level wind speed of 4.2 m s–1.

The rain sampling was conducted at the air quality monitoring station located in the main campus of Imam Abdulrahman Bin Faisal University (26.390°N, 50.182°E). The station is 6.7 m above the sea level and 2 km away from the coastline. The location of the air quality monitoring station and the land use map of the study area are shown in Fig. 1. An Eigenbrodt® automatic wet-only sequential rain sampler was used. The total height and collection surface area of the sampler are 1.6 m and 500 cm2, respectively. The sampler had been operated continuously during the rainy season between October 2015 and May 2016. The start time of the sampling was set to 09:00 a.m. as it is recommended (GAW, 2004). A total number of 21 rainwater samples were collected between November 11th and April 17th with a total precipitation depth of 74.4 mm.

Fig. 1. Location of the sampling site and the land use map of the study area.Fig. 1. Location of the sampling site and the land use map of the study area.

Chemical Analysis

The pH and electrical conductivity (EC) measurements of the samples were performed immediately after the sample collection by using Consort™ C3010 Benchtop Multiparameter Analyzer. Following this, each sample was filtered through a 0.45-µm cellulose acetate membrane filter using a Nalgene® filtration system. The filtrate part (soluble fraction) of the sample was transferred into a Nalgene® bottle, sealed tightly, and stored in a refrigerator at 4°C prior to the time of analysis. Major anions (Cl, NO3, and SO42) and cations (Na+, K+, NH4+, Ca2+, and Mg2+) were quantified within a week after sample collection by using a high-performance liquid chromatograph (Prominence; Shimadzu). Bicarbonate ion (HCO3) in rainwater samples were determined using Standard Methods: 2320 B. Titration Method.

Quality Assurance Measures

The field and laboratory blanks were prepared to correct the possible contaminations of the samples during the sampling, handling, preparation, storage, and measurement steps. Field blanks were collected bimonthly by routine washing of the rainwater collection bucket with 100-mL ultrapure water. Laboratory blanks were prepared and carried through the same procedure as the samples along with each filtration set. The blanks were treated as a real sample and included in all analyses. The ultrapure water with a resistivity of 18.2 MΩcm used in rinsing, blank preparation, dilution, and calibration steps was produced by an ELGA® water purification system.

The detection limit (DL) of each analyte measured in rainwater samples was calculated to be mean of the blank values plus three times the standard deviation of the blank values (de Moraes Dias, et al., 2012). The analytical protocols were verified before measurements of the samples by analyzing Standard Reference Materials (SRMs; simulated rainwater SR-1 and SR-2; High-Purity Standards®). The precision of each method was determined by calculating the relative standard deviation (RSD; %) of 10 replicate measurements of the standards and samples. The possible contaminations of the samples were checked by blank to sample ratios. The calibration range, linearity (the coefficient of determination (R2) of the calibration), DL, precision (%), and sample to blank ratio (S/B) of each analyte are indicated in Table 1.

Table 1. Quality control performances of measured ions.

The calculated linearities ranging from 0.9989 to 0.9998 were found to be acceptable. The analyte concentration per sample was substantially higher than its DL. The precision results varying between 1.23% and 9.29% were satisfactory. S/B ratios were above 3, suggesting that the contaminations of the samples were negligible.

The equivalent ratio of the total anions to total cations (∑anions/∑cations) has been used as an indication of the analytical data quality and completeness of measured parameters (Al-Momani et al., 1995a, b). If this ratio is within the interval of 1 ± 0.25, data are commonly assumed acceptable (Keene et al., 1986). The mean ∑anions/∑cations ratio for the collected samples is 0.94 ± 0.34, suggesting that most of the important ions were included in measurements.

Data Analysis

Volume weighted mean (VWM) values are commonly used to avoid such variability in measured parameters due to precipitation amount. Wet deposition flux (FWD) could be defined as the deposited amount of atmospheric chemicals onto unit area of Earth’s surface via hydrometeorological processes. The VWM concentration and FWD of pollutants in rainwater were calculated by using Eqs. (1) and (2), respectively (Alagha and Tuncel, 2003):

where CVWM is the VWM concentration, [Xi] is the concentration of an ion in the ith rainwater sample, Pi is the rainfall depth (mm) during the ith rain event, and:

where FWD was calculated as mg m–2 or kg km–2 and A is the collection surface area of the sampler.

Marine and crustal enrichment factors are commonly practiced in identifying marine and crustal sources of the ions in rainwater, respectively. Marine enrichment factor (EFm) is calculated using Na+ as a reference for marine source (Eq. (3)) (Kulshrestha et al., 1996), while crustal enrichment factor (EFc) is determined using Ca2+ as a reference for crustal source (Eq. (4)) (Al-Momani et al., 1997):

where [X/Na+]rain and [X/Na+]sea are the ratios of concentration of an ion (X) to Na+ concentration in rainwater sample and seawater, respectively, and [X/Ca2+]rain and [X/Ca2+]crust are the ratios of concentration of an ion (X) to Ca2+ concentration in rainwater sample and earth crust, respectively. [X/Na+]sea and [X/Ca2+]crust ratios were obtained from Keene et al. (1986) and Taylor (1964).

The sea-salt fraction (SSF), crustal fraction (CF), and anthropogenic fraction (AF) of each ion measured in rainwater samples could be estimated by using Eqs. (5), (6), and (7), respectively (Keene et al., 1986; Cao et al., 2009; Lu et al., 2011):


where [ss – X] and [c – X] are the sea-salt origin and crustal origin concentrations of an ion in rainwater, respectively; [X] is the concentration of an ion in rainwater; and [Na+]rain and [Ca2+]rain are the Na+ and Ca2+ concentrations in rainwater, respectively.

The strong acids such as H2SO4 and HNO3 mainly produced from anthropogenic emissions dominate the acidity in wet deposition. Considering that nss-SO42 and NO3 are two principal acidic components in rainwater, neutralization of acidity can be assessed by fractional acidity (FA). The acidity contributed by H2SO4 and HNO3 is not neutralized at all if FA is equal to unity (Balasubramanian et al., 2001). FA can be estimated using Eq. (8):

where [H+] is the H+ concentration calculated from pH and nss-SO42 is the non-sea-salt SO42 concentration in rainwater sample.

The free acidity of rainwater resulting from NO3 and nss-SO42 is buffered by NH4+, Ca2+, K+, and Mg2+. Neutralization factors (NFs) are used to quantitatively reveal the roles of these cations on the neutralization of free acidity in rainwater by using Eqs. (9) and (10) (Possanzini et al., 1988; Kulshrestha et al., 1996):

The ratio of neutralizing potential (NP) to acidifying potential (AP) is used to assess the balance between acidity and alkalinity, and it is calculated by Eq. (11) (Alagha and Tuncel, 2003):

Factor analysis (FA) is a multivariate statistical method that has been frequently used to identify the sources of chemical composition in a wet deposition (Başak and Alagha, 2004; Uygur et al., 2010; Anil et al., 2017). This technique simply separates the variables into groups called factors according to their variation and common variance. The obtained factors represent a specific source or mixed sources of the analyzed chemical species. Varimax rotation, the most widely used rotation method, was preferred in rotating the initial factor loadings orthogonally in order to simplify the factor structure and make the factor loadings easily interpretable.

Cluster Analysis of Backward Air Mass Trajectories

The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model (Version 4, 2018 release) was used to compute backward air mass trajectories reaching the receptor site (Calvo et al., 2012; Stein et al., 2015; Rolph et al., 2017). GDAS (Global Data Assimilation System) 1° HYSPLIT-compatible meteorological data produced by the National Weather Service’s National Centers for Environmental Prediction (NCEP) were added to “trajectory setup” (Shi et al., 2014; Xin et al., 2016). The endpoint files of 3-day backward trajectories arriving at 500 m above ground level (AGL) were calculated for every 6 h per day within the study period. 84 backward trajectories were calculated and then classified regarding their typical patterns by using “trajectory clustering analysis” feature of HYSPLIT 4. The software initially calculates 30 clusters, among which the optimum number of clusters is determined by using an objective percent change in total spatial variance (TSV) criterion of either 20% or 30%.

Satellite Observation Data

The satellite observation data of the following parameters over the study period between November 2015 and April 2016 were downloaded from the EARTHDATA GIOVANNI database (Acker and Leptoukh, 2007): i) aerosol optical depth at 0.55 microns for land and ocean (daily 1° [MODIS-Terra MOD08_D3 v6.1]), ii) sea-salt surface mass concentration (monthly 0.5 × 0.625° [MERRA-2 Model M2TMNXAER v5.12.4]), iii) SO2 surface mass concentration (monthly 0.5 × 0.625° [MERRA-2 Model M2TMNXAER v5.12.4]), and iv) NO2 tropospheric column (30% cloud screened, daily 0.25° [OMI OMNO2d v003]).


Variations of Rainfall, pH, and EC

Time series and correlation statistics of rainfall, duration of rainfall, pH, and EC are demonstrated in Fig. 2 in order to give a preliminary comprehension of the general features of the collected rainwater data. The minimum collected volume of rainfall was 0.36 mm (on Dec. 15) and the maximum was 19.0 mm (on Nov. 25). During the rainy season, the sampling area received less than 25% of the long-term average annual precipitation total of 100 mm. A positive strong correlation between rainfall and duration of rainfall was observed, indicating that rainfall intensity per each rain event was almost steady and has an average of 1.1 ± 1.0 mm h–1. The pH values fluctuating between 6.4 and 8.3 revealed the alkaline nature of rainwater over the study area as the pH of the natural rainwater in equilibrium with atmospheric CO2 is 5.6 (Granat, 1972). The time series trend of pH is not quite similar to that of rainfall. The correlation between pH and rainfall depth were negative and moderate, which signifies that the pH was partially influenced by rainfall. The EC values varied between 22.7 and 264 µS cm–1, with an average of 126 ± 66.6 µS cm–1. A positive and moderate correlation among pH and EC could be linked to the possible neutralization of acidic constituents by inputs of alkaline substances into rainwater, which will be thoroughly discussed in subsequent sections. The EC of precipitation mainly reflects the total soluble components of the atmosphere which are comparatively bigger in size and tend to be scavenged from the atmosphere more quickly than smaller particles. A major part of the suspended components in the atmosphere is washed out in the first few mm of the precipitation (Anil et al., 2017). For these reasons, the EC value is generally high in the initial stages of the rainfall and gradually drops as the rainfall continues, explaining the negative strong correlation observed among EC and rainfall during the study period.

Fig. 2. Time series of rainfall, duration, pH, and EC.Fig. 2. Time series of rainfall, duration, pH, and EC.

The monthly variations of rainfall, duration, pH, and EC are depicted in Fig. 3. Even though total durations of rainfall on Nov. 2015 and Dec. 2015 were very close, the total rainfall of Nov. 2015 was twofold higher than that of Dec. 2015 due to the presence of more intense rain events on Nov. 2015. The highest rainfall was recorded on Nov. 2015 with a total depth of 30.3 mm followed by Dec. 2015 (16.3 mm) and Mar 2016 (14.2 mm). Rainfall within these three months accounted for 82% of the total precipitation depth during the study period. In this period, EC values were lower than those of rain events collected on other months (Jan., Feb., and April) due to the dilution effect of abundant rainfall on the total soluble components in the atmospheric column. The negative correlation between rainfall and EC parameters shown in Fig. 3 becomes clearer and stronger (r = –0.94) than that in Fig. 2. On the other hand, the highest EC values were measured as 216 µS cm–1 and 162 µS cm–1 on Jan. 2016 and Feb. 2016 where the lowest precipitation depths were received. The monthly VWM values of pH showed slight variation during the study period (SD = 0.49) and ranged between 6.73 and 8.34 with an overall VWM of 7.34. The lowest monthly VWM of pH was recorded on Nov. 2015 and accompanied with the highest monthly rainfall total and the lowest EC. 

Fig. 3. Monthly variations of rainfall, duration, pH, and EC.Fig. 3.
 Monthly variations of rainfall, duration, pH, and EC.

Chemical Composition

The descriptive statistics of the measured parameters in Dammam rainwater such as VWM, mean, median, standard deviation (std. dev.), minimum (min.), maximum (max.), percentiles, and percent distribution of ions in rainwater are listed in Table 2. From the descriptive statistics data, arithmetic means and medians are commonly higher than VWMs, implying that the higher concentrations of ions are usually coupled with lower rainfall amounts. The reason behind this inverse relation could probably be attributed to the dominant wash-out mechanism during the precipitation (Wu et al., 2016; Szép et al., 2018). The high standard deviations of all ions revealed a wide fluctuation of the concentrations in individual precipitation events. The coefficients of variation (Std.Dev/CVWM) of HCO3, Na+, and NO3 were more than the unity, indicating that these ions had larger variations during the study period. 

Table 2. Descriptive statistics of the measured parameters in Dammam’s rainwater.

The VWM concentrations of the ions were in the order of Ca2+ > SO42 > Cl > HCO3 > Na+ > NO3 > Mg2+ > NH4+ > K+ >> H+, revealing that 80% of the rainwater chemical composition was dominated by Ca2+, SO42, Cl, HCO3, and Na+. Following this, the contributions of NO3, Mg2+, NH4+, and K+ to the total ionic composition were 6.93%, 5.95%, 5.88%, and 0.942%, respectively. H+ was the least abundant ion in the precipitation with the contribution of 0.007%. Among the precipitation components, Ca2+ was the most abundant single ion with VWM concentration of 480 µeq L1, accounted for 30.1% of the total ions and 58.6% of the total cations. The VWM concentration of Ca2+ together with Na+ reached 615 µeq L1, which explained 75% of the total cation budget in precipitation. The VWM concentration of acidic components, SO42 plus NO3, was 466 µeq L1, occupying 60% of all anions. Marine-originated ions (Na+ + Cl) with VWM concentration of 306 µeq L1 represented 19.2% of the total measured components.

The VWM concentrations of the parameters measured in this study were compared to the literature data from nearby and remote sites in Table 3. Compared to the previous study conducted during 1987 and 1988 (Ahmed et al., 1990), the NH4+, SO42, K+, and Cl concentrations increased by 106%, 53%, 31%, and 19%, respectively, whereas NO3 indicated 23% reduction. The relative changes in other ions were found to be less than 7%. This is the most interesting to note that the pH value increased from 5.48 to 7.34 in response to decreasing NP/AP ratio from 1.35 to 1.23. In addition to this, the most recent data report demonstrated that SO2 and NOx emissions over Saudi Arabia increased 75.6% and 47.7%, respectively, through the years 1988 and 2016 (Crippa et al., 2018). In this regard, one could expect lower pH value and higher NO3 concentrations in this current study compared to the previous one. These temporal changes in the precipitation chemistry of study area could be attributed to the variations in: i) atmospheric and meteorological conditions, ii) atmospheric residence times of the ions, iii) mechanisms governing wet scavenging of ions, and iv) long-range transport of air pollutants.

The pH value found in this study was the second highest among those listed in Table 3. The former studies conducted at Saudi Arabia (Riyadh), Jordan, Iran, and Turkey reveal that the pH structure of the Middle East region is alkaline, which is commonly linked to the calcareous soil type of the region, rather than lack of acid-forming ions in the rain (Alabdula'aly and Khan, 2000; Al-Khashman, 2009; Anil et al., 2017; Bayramoğlu Karşı et al., 2018; Naimabadi et al., 2018). The Ca-rich suspended particulate matter over the region, originating from local soil resuspension and long-range particulate matter transport from the Fertile Crescent region and Saharan desert, is the main reason enabling Ca2+ to be the most abundant acid-neutralizing component in the regional wet deposition (Ahmed et al., 1990; Alabdula'aly and Khan, 2000). The VWM concentrations of Mg2+, NH4+, and K+, acid-neutralizing species rather than Ca2+, were ranked as the median among those listed in Table 3.

The total concentration of the SO42 and NO3 in the study area (466 µeq L1) succeeded the concentrations reported in Ahvaz, Iran (858 µeq L1); Riyadh, Saudi Arabia (722 µeq L1); and Bolu, Turkey (487 µeq L1) (Alabdula'aly and Khan, 2000; Bayramoğlu Karşı et al., 2018; Naimabadi et al., 2018). The higher concentration of total acidic components (SO42 + NO3) experienced in Dammam rainwater than those reported for megacities in Turkey (Istanbul), Pakistan (Karachi), India (New Delhi), and China (Beijing) could be attributed to: i) intense burning of fossil fuel for energy production and transportation activities; ii) accumulated SO42 and NO3 particles on the atmosphere due to limited precipitation and frequent dusty weather conditions, and iii) transboundary pollution.

HCO3 was found to be the second most abundant anion in rainwater monitored in Iran, India, Pakistan, Saudi Arabia, Jordan, and Romania (Ahmed et al., 1990; Al-Khashman, 2009; Kumar et al., 2014; Masood et al., 2018; Naimabadi et al., 2018; Szép et al., 2018). HCO3 concentration measured in this study (141 µeq L1) was lower than those observed in Ahvaz, Iran (667 µeq L1); New Delhi, India (205 µeq L1); and Karachi, Pakistan (164 µeq L1).

The equivalent ratio of Na+/Cl must be within the range between 0.5 and 1.5 to indicate a marine origin (Martins et al., 2019). The Na+/Cl ratio found in this study was 0.79 and only 8.1% lower than the seawater ratio of 0.86. The Na+/Cl ratios reported from the sampling sites listed in Table 3 are between 0.5 and 1.5 since they are located close to the coast. The sea-salt aerosol scavenging enriched the Na+ and Cl concentrations in rainwater in those stations except Karachi, Pakistan, and Beijing, China.

The ratio of total cations to total anions (∑cations/∑anions) of this study was calculated as 1.05, inferring a balance between total cations and total anions. This ratio close to unity gives an indication on neutralization of the rainwater acidity since the most abundant ions of the collected rainwaters samples were Ca2+ and SO42. Table 3 reveals that ∑cations/∑anions ratio of this study was lower than those of many studies except those reported for Turkey, Istanbul (0.81), and Ahvaz, Iran (0.79).

Wet Deposition Fluxes

The monthly variations of the wet deposition fluxes (FWD) of measured ions are presented in Fig. 4. The FWD of acid-neutralizing and soil-derived species (Ca2+, HCO3, NH4+, and Mg2+), acidic components (SO42 and NO3), and marine-originated ions (Na+, Cl, and K+) indicated similar temporal trends between each other. The highest FWD of SO42, NO3, Cl, Na+, Mg2+, and K+ on Nov. 2015 could be linked to the frequent rain events and highest rainfall amount on this period. The FWD of Ca2+ and HCO3 reached their maximums on Mar 2016 probably due to the higher accumulation of crustal particles appeared during January and February, months in which less rainfall was observed. Besides this, the lowest FWD of all measured ions were recorded on January and February, respectively.

Fig. 4. Monthly variations of wet deposition fluxes.Fig. 4. Monthly variations of wet deposition fluxes.

The monthly FWD were strongly correlated (0.704 < r < 0.986, p < 0.01) with the rainfall amount. In addition to that, a significant linear relationship between the total FWD of each event and the associated rainfall amount was found to be r = 0.848 (Fig. 5). It can be concluded from Fig. 5 that rainfall amount significantly contributes to the FWD of ions since the total ionic FWD exhibited a similar fluctuation tendency with the rainfall amount as shown in Fig. 5. In the meantime, the FWD of major chemical components over the study area was also affected by the temporal variations of emissions from crustal, marine, and anthropogenic sources. 

Fig. 5. Time series of total ionic wet deposition fluxes.Fig. 5. Time series of total ionic wet deposition fluxes.

The total FWD of the ions during the study period were in the order of SO42 (1,269 mg m–2) > Ca2+ (714 mg m–2) > HCO3 (638 mg m–2) > NO3 (510 mg m–2) > Cl (450 mg m–2) > Na+ (231 mg m–2) > NH4+ (126 mg m–2) > Mg2+ (85.8 mg m–2) > K+ (43.7 mg m–2) > H+ (8 µg m–2). During the study period, the Dammam metropolitan area received a total ionic FWD of 4.07 tons per km2 which was lower than those reported in Po Valley, Italy (43.3 tons km–2 y1) (Tositti et al., 2017); Beijing, China (23.3 tons km–2 y1) (Xu et al., 2015); Karachi, Pakistan (16.8 tons km–2 y1) (Masood et al., 2018) Istanbul, Turkey (9.01 tons km–2 y1) (Anil et al., 2017); Carpathians, Romania (7.13 tons km–2 y1) (Szép et al., 2018); and Ahvaz, Iran (4.57 tons km–2 y1) (Naimabadi et al., 2018). Consequently, rare precipitation events with low-volume rainfalls experienced in arid and semiarid regions such as Saudi Arabia commonly cause higher concentrations of ions in rainwater and lower annual FWD of ions relative to the wetter regions receiving more frequent rain events with larger volumes (Keene et al., 2015).

Acid Neutralization

The pH value of the rainwater is mainly controlled by the neutralization between acidic and basic components. The fractional acidity (FA), correlations between acidic and alkaline species, SO42/NO3 ratio, and neutralization factors (NFs) are the most well-known approaches to clarify the factors governing the acid neutralization in rainwater. The mean FA value of all rainwater samples was calculated to be 0.0003 ± 0.001, revealing that 99.97% of the acidity had been neutralized by alkaline species. The positive and strong linear correlation between ([Ca2+] + [Mg2+] + [NH4+] + [K+]) and ([SO42] + [NO3]) (r = 0.894) also supports that the alkaline species measured in the study were responsible for the neutralization of the rainwater acidity and the contribution of any other alkaline components to the neutralization process was not significant. The relative contributions of SO42 and NO3 in the precipitation acidity can be evaluated by using the equivalent ratio of [SO42–]/([NO3] + [SO42–]). The mean value of this ratio during the study period was 0.756 ± 0.039, reflecting that about 75.6% of rainwater acidity was due to SO42 while the contribution of NO3 was 24.4%. Comparing to the former study, the equivalent [SO42]/[NO3] ratio has increased from 1.61 to 3.21 since 1990, indicating that H2SO4 has an increasing effect on the precipitation acidity in Dammam (Ahmed et al., 1990). The NFs of Ca2+, Mg2+, NH4+, and K+ were calculated in order to evaluate the rainwater neutralization by alkaline components in detail. The Ca2+ showed the highest NF value (1.09) and followed by Mg2+ (0.135), NH4+ (0.117), and K+ (0.02). These results suggest that Ca2+ was the principal neutralizer in rainwater and its contribution to the total neutralization was about 80%. Furthermore, Ca2+ plays a dominant role in the neutralization process of rainwater in the Middle East due to the Ca-rich soil structure of the region.

Source Apportionment

Enrichment Factors and Relative Source Contributions

The enrichment factor (EF) model, a double normalization technique, has been used to determine possible sources of ions in rainwater (Başak and Alagha, 2004). An EF value of an ion much less than 1 or much higher than 1 is assumed to be diluted or enriched in regard to the reference source (Al-Momani et al., 1997). In this study, Na+ was selected as a marine reference ion since average equivalent ratios of [Cl]/[Na+] and [Mg2+]/[Na+] in rainwater samples were higher than those in standard seawater (Zhang et al., 2012). Ca and Al are commonly used as a soil reference element to identify the crustal origin of elements in rainwater (Okay et al., 2002). Ca2+ was selected as a crustal reference since Al was not measured in the rainwater samples. The mean values of marine enrichment factors (EFm) and crustal enrichment factors (EFc) of ions are given in Table 4. The EFm of NO3 was not calculated since its abundance in standard seawater is very meager (Cao et al., 2009).

Table 4. Enrichment factors of ions in Dammam’s rainwater relative to standard seawater and earth crust.

The EFm and EFc values of Cl were found to be 1.20 and 389, respectively, indicating that the origin of Cl in rainwater was mainly sea-salt particles due to the proximity of the sampling location to the western coast of Arabian Gulf. The EFm of Cl was slightly higher than the recommended seawater ratio of 1.16, which could be proof of possible terrestrial contribution. The EFm and EFc values of K+ and Mg2+ reveal that these ions were slightly enriched compared to seawater and diluted to earth crust. These results indicated that K+ and Mg2+ mostly originated from crustal sources and partly from marine sources. The high EFm of Ca2+ (117) showed that marine contribution to Ca2+ was negligible and crustal-originated particles were the main source of Ca2+ in rainwater. The contributions of NO3 and SO42 to rainwater were not mainly from marine or crustal sources since the EFc value of NO3 and the EFm and EFc values of SO42 were all much higher than 10.

The approximate sea-salt fraction (SSF), crustal fraction (CF), and anthropogenic fraction (AF) of each ionic component in rainwater were estimated by using Eqs. (5), (6), and (7) and shown in Table 5. In this study, Na+ was assumed to be of purely marine origin. Table 5 shows that 90.1% of Cl originated from sea-salt particles, 9.19% from anthropogenic sources, and only 0.705% from earth crust. Ca2+, K+, and Mg2+ in rainwater were mainly of crustal origin whereas their CFs were calculated to be 98.7%, 82.0%, and 68.7%, respectively. The contribution from marine sources to Mg2+ accounting for 31.3% was also considerable in comparison with other region-specific tracer cations of the earth crust. The terrestrial contribution to these three cations could be ascribed to only earth crust because it is very difficult to differentiate their crustal and anthropogenic sources with the available rainwater data (Cao et al., 2009; Xiao, 2016). The anthropogenic sources indicated the highest effect on NO3 and SO42 in rainwater with a contribution of 99.6% and 93.0%, respectively. The SSF and CF of SO42 were estimated to be 5.67% and 1.37%, respectively. Apart from AF, NO3 in rainwater had a very slight contribution only from crustal sources. In Saudi Arabia, road transportation, public electricity and heat production, manufacturing industries, and construction sectors were reported to be the most significant emission sources of atmospheric NOx and SO2 which are precursors of anthropogenic NO3 and SO42 in regional rainwater (Ahmed et al., 1990; Alabdula'aly and Khan, 2000; Crippa et al., 2018). 

Table 5. Marine and terrestrial contributions to ions in Dammam’s rainwater.

Factor Analysis

In this study, “principal components method” and “eigenvalues greater than 1” options were selected during the extraction of factors using SPSS ver. 24. Three major factors with eigenvalues higher than unity were extracted and the results are given in Table 6. The extracted factors explained 92.6% of the total variance of the dataset. The communalities of all variables were higher than 0.855, indicating that the extracted factors are reasonable.

Table 6. Factor analysis results of ions in the precipitation of Dammam.

Factor 1 explaining 51% of the total variance is the most significant factor and characterized by high loadings of Ca2+, K+, Mg2+, HCO3, SO42, and NO3. This group of ions describes mixed sources since it includes crustal markers and anthropogenic components. The surface soil structure of the study area and the region typically contains limestone (CaCO3 and CaMg(CO3)2) and is partially interbedded with feldspar (Al2SiO5 of K or CaO), gypsum (CaSO4·2H2O), and halite (NaCl) (Alsaaran, 2008). Considering this information, it can be inferred that accumulated coarse particles on the atmosphere due to resuspension of local soil by surface wind and long-range particulate transport from nearby deserts on North Africa, Fertile Crescent region, and Arabian Peninsula are the main sources of the enriched crustal cations, particularly Ca2+ in Dammam’s rainwater. The calcium bicarbonate (Ca(HCO3)2) is formed as a result of the reaction between the dissolved CO2 in rainwater and CaCO3. This soluble compound is eventually washed out from the atmosphere by wet deposition (Kulshrestha et al., 2003; Demirak et al., 2006). The considerable abundance of HCO3 in Dammam’s rainwater, 8.81% of total chemical composition, and a strong correlation between Ca2+ and HCO3 (r = 0.889) unveil two different scenarios that complement each other: i) CaCO3-rich particles over the atmosphere during the rain events and ii) the presence of elevated CO2 emissions to the atmosphere due to the burning of fossil fuels during the study period. The coexistence of SO42 and NO3 in this factor and strong correlation between these acidic components (r = 0.951) could be linked to their similar chemical behaviors in the atmosphere and common sources of their precursors SO2 and NOx. It was reported that emissions of SO2 and NOx in Saudi Arabia in 2016 were 3.37 Mt and 1.84 Mt, respectively (Crippa et al., 2018). Public electricity, heat production, and other energy industries were responsible for 78% of SO2 emissions while the contributions of the energy sector and road transportation to NOx emissions were 62% and 27%, respectively. The significant correlations between the acidifying components and crustal cations in this factor could be attributed to neutralizing reactions between H2SO4 and HNO3, and alkaline compounds in rainwater. Therefore, SO42 and NO3 in the precipitation mostly existed as the forms of CaSO4 and Ca(NO3)2 in this study whereas Ca2+ was responsible for 80% of the total neutralization in collected rainwater samples (see Section: “Acid Neutralization”).

The second factor accounting for 28.0% of the total variance indicated high loadings of Na+ and Cl and moderate loadings of Mg2+ and K+. This factor clearly describes the marine source and the incorporation of sea-salt particles into rainwater. In this study, the source of the Na+ in rainwater was assumed to be totally from the sea and the correlation coefficients of this ion with Cl, Mg2+, and K+ were found to be 0.995, 0.661, and 0.575, respectively. The estimated SSF fractions of Cl (90.1%), Mg2+ (31.3%), and K+ (18.0) from the previous section also supports the dominance of Na+-Cl couple and partial effects of Mg2+ and K+ in this factor. Mg2+ and K+ have higher loadings in Factor 1, explaining that their origins in rainwater were mainly from crustal sources and partly from marine sources. The existence of Cl with a weak loading (0.267) in Factor 1 could explain its anthropogenic fraction (9.19%) possibly due to emissions from fossil fuel combustion, burning of waste materials, and automobile exhaust within the study area (Palmer, 1976; Wu et al., 2016; Szép et al., 2018).

Factor 3 represents 13.6% of the total variance with high loading of NH4+ and moderate loadings of SO42 and NO3. Atmospheric NH4+ mainly originates from anthropogenic sources such as agricultural activities, fertilizer factories, biomass burning, and emissions from vehicles. The atmospheric NH3 reacts rapidly with H2SO4 and HNO3 to form secondary fine particles of (NH4)2SO4 and NH4NO3. In this neutralization reaction, the reaction of NH3 with H2SO4 is favored over reaction with HNO3. Since the atmospheric settling velocity of SO42 particles is slower than NH3 and HNO3, fine (NH4)2SO4 particles dominate the distribution of atmospheric NH4+ over longer distances (NOAA, 2000). In this regard, the higher factor loading of SO42 compared to NO3 in this factor suggests that (NH4)2SO4 particles could be more abundant than NH4NO3 particles over the study area. The possible sources of NH4+ in Dammam’s rainwater might be resulting from: i) local NH3 emissions from vehicles but not from other local sources since agricultural activities and biomass burning are quite limited in the study area and ii) long-range transport of secondary fine particles of (NH4)2SO4 and NH4NO3 over the regions where SO2, NOx, and NH3 emissions are intense.

Cluster Analysis Coupled with Satellite Observations

The cluster analysis of backward air mass trajectories arriving at the sampling area for the rainy days of the study period resulted in a 7-cluster solution. This was determined as the optimum number of clusters regarding its TSV changes among other cluster solutions. The percentages of the classified 3-day backward air mass trajectories, clustered rain event codes, and total FWD values of ions in each cluster were calculated and shown in Table 7. The obtained satellite observation data of aerosol optical depth at 0.55 microns, sea-salt surface mass concentration, SO2 surface mass concentration, NO2 tropospheric column (30% cloud screened), and the average backward air mass trajectories in 7-cluster solution were visualized for the same domain (Region 9E, 15N, 60E, 45N) by using ArcGIS ver. 10.3 and depicted in Fig. 6.

Table 7. Cluster statistics and total FWD of ions (kg km–2) in each cluster.

During the study period, 53% of the observed precipitation events were mainly controlled by the air masses in Clusters 2, 4, 5, and 6 arriving at the sampling site from north-northwest (NNW) sector. Even though Clusters 2, 4, and 6 arose from different regions, they followed a similar pathway with Cluster 5 for about 24 h before their arrival to the sampling area. Their common trajectory started from southern Iraq, traversed Kuwait and Arabian Gulf, respectively, and they reached Dammam from NNW direction. According to their longer pathways among others, Clusters 4 and 6 originating from central Tunisia and Black Sea region of Turkey, respectively, were ascribed to fast-moving trajectories with the mean altitudes of 2,412 m and 2,019 m, respectively. On the other hand, trajectories with a mean altitude of 1,333 m in Cluster 5 were the slowest ones among northwesterly clusters. The air masses deriving from the west-southwest (WSW) sector and arriving at Dammam from northeast (NE) direction were attributed to Cluster 1. This group of trajectories with a mean altitude of 1,650 m representing 18% of all rain systems of the study could be characterized as marine- and terrestrial-impacted nationwide cluster as they originated from the Red Sea coast of Jeddah, passed over Western Province, Riyadh, and Arabian Gulf coast of the study region. The rain systems in Clusters 3 and 7 reaching the sampling site from southerly directions accounted for 29% of all precipitation events. The trajectories in Cluster 7 deriving from Oman transported over southeast (SE) part of Saudi Arabia, the western section of United Arab Emirates, Arabian Gulf, Qatar, Bahrain, and eventually reached over Dammam from southeast (SE) direction. The rain systems in this cluster could be enriched by marine and terrestrial sources throughout their trajectories. Cluster 3 representing 16% of the total back trajectories originated from 50 km off the Dammam shoreline (on the Arabian Gulf side), made a circular motion in a clockwise direction around Qatar, entered the land border of Saudi Arabia between Qatar and United Arab Emirates, and completed its circular motion by reaching the sampling point from south direction. Among all clusters, this distinctive group of trajectories was identified as the slowest- and lowest-altitude-moving ones due to their shortest range and lowest mean altitude (673 m AGL).

Clusters 1, 3, and 5 describing 56% of the total back trajectories dominated 75% of the total FWD during the study period. The rain systems in Cluster 1 represented 32% of the total FWD, including the highest depositions of anthropogenic ions (SO42, NO3, and NH4+) and high depositions of crustal and sea-salt components on the study area (Table 7). The potential long-range and medium-range sources of the highest depositions of these anthropogenic ions were highly populated and heavily industrialized cities throughout the route of Cluster 1 such as Western Province and Riyadh where high emissions of SO2 and NO2 were observed (Figs. 6(c) and 6(d)). The southern part of Saudi Arabia indicated moderate to high aerosol optical depth values within the map domain after hot spots over Chad, Niger, and Oman (Fig. 6(a)), which could be attributed to the high FWD of crustal ions brought by air masses in this cluster. The back trajectories deriving over the Red Sea have potentials to transport high quantities of Na+ and Cl to the study area due to the highest sea-salt concentrations observed over the Red Sea throughout the map domain (Fig. 6(b)). 

Fig. 6. (a) Aerosol optical depth at 0.55 micron for land and ocean, (b) sea-salt surface mass concentration (µg m–3), (c) SO2 surface mass concentration (µg m–3), (d) NO2 tropospheric column (1 × 10–14 cm–2), and (e) average 3-day backward trajectories arriving at 500 m AGL in 7-cluster solution.Fig. 6. (a) Aerosol optical depth at 0.55 micron for land and ocean, (b) sea-salt surface mass concentration (µg m3), (c) SO2 surface mass concentration (µg m3), (d) NO2 tropospheric column (1 × 1014 cm2), and (e) average 3-day backward trajectories arriving at 500 m AGL in 7-cluster solution.

Cluster 5 representing the highest fraction of the total back trajectories (23% of total) accounted for 27% of the total FWD over the study area. The rain events associated with the air masses from this cluster were characterized by the highest FWD of both crustal and sea-salt components and also high FWD of anthropogenic ions due to high particle loads and high emissions of SO2 and NO2 over southern Iraq, Kuwait, and the Dammam metropolitan area and also moderate sea-salt mass concentration over the Arabian Gulf throughout the pathway of their travel. In addition to this, the slow- and medium-altitude-moving characteristic of the air masses in this cluster resulted in efficient transport and better enrichment of the species in the atmospheric columns along the route from southern Iraq to the receptor point. Among the source regions, Kuwait was possibly the most important contributor to the FWD of anthropogenic components since it was the highest SO2 and NO2 emitter within this cluster.

Similar to Cluster 5, atmospheric species within the path of Cluster 3 were efficiently acquired by the rain clouds and deposited over the receptor area as this cluster grouped the slowest- and lowest-altitude-moving air masses and indicated more local characteristics as well. This cluster explaining 16% of the back trajectories represented the third highest FWD of both crustal and anthropogenic ions and medium sea-salt enrichment as well. Apart from the local effects on the total FWD, eastern part of Qatar emitting considerable amounts of NO2 partially contributed to FWD of NO3.

Even though Clusters 2, 4, 6, and 7 together described 44% of the total back trajectories, the regions that were passed over by all these clusters contributed only 25% to the total FWD over the study area. This could be attributed to: i) limited enrichment capacities of fast- and high-altitude moving systems in Clusters 2, 4, and 6 by atmospheric species on the atmospheric columns that they were passing through and ii) low trajectory frequencies of the individual clusters. Among these four clusters, air masses in Cluster 6 showed the highest contribution (11%) to the total FWD. About quarter of both Na+ and Cl, 14% of K+, 12% of Mg2+, 11% of Ca2+, 10% of NO3, and 8% of SO42 in total FWD corresponded to the rain systems in Cluster 6. The source locations of the sea-salt ions enriched together were substantially from the Arabian Gulf and partially from the central Black Sea region of Turkey where the air masses in this cluster originated. Like in Cluster 5, the anthropogenic emissions from Iraq and especially from Kuwait were the long-range origins of NO3 and SO42 deposited to the sampling area by rain systems in Cluster 6. The crustal components transported by Cluster 6 could be attributed to the medium particle loads over Iraq and Kuwait shown in Fig. 6(a).


21 samples of rainwater, which totaled 74 mm in rainfall depth, were collected by an automatic wet-only precipitation sampler during the rainy season between November 2015 and April 2016 in Dammam, Saudi Arabia. The VWM value of the pH was found to be 7.34 (± 0.49), indicating the alkaline nature of the rainwater in the study area. The chemical composition of the rainwater was dominated by Ca2+, SO42, Cl, HCO3, and Na+, which accounted for 80% of the total ionic equivalent concentration. Ca2+, with a VWM concentration of 480 µeq L1, was the most abundant ion and contributed 30.1% of the total ionic equivalent concentration. During the study period, the Dammam metropolitan area received a total ionic FWD of 4.07 tons per km2, which consisted of the following components in descending order: SO42 > Ca2+ > HCO3 > NO3 > Cl > Na+ > NH4+ > Mg2+ > K+.

The concentration of the combined SO42 and NO3 in the study area (466 µeq L1) was higher than those reported for megacities in Turkey, Pakistan, India, and China. Approximately 75.6% of the rainwater acidity was due to SO42 and 24.4%, to NO3. The extremely low fractional acidity of the rainwater indicated that almost all of the acidity was neutralized by alkaline species, with Ca2+ being responsible for 80% of the samples’ neutralization, revealing that SO42– and NO3 were mostly present as CaSO4 and Ca(NO3)2 in the precipitation.

Source apportionment based on enrichment factor, relative source contribution, and factor analysis as well as cluster analysis coupled with satellite observations showed that the anthropogenic components (SO42, NO3, and NH4+), which contributed 47% of the total ionic FWD, mainly arose from i) the intensive burning of fossil fuel for energy production and transportation activities in the vicinity of the study area and ii) the medium-range and long-range transport of anthropogenic pollutants from Kuwait, Iraq, Qatar, the Western Province, and Riyadh. Local soil resuspension, moderate to high particulate matter loads over the southern part of Saudi Arabia, and long-range transport from Iraq and Kuwait contributed the crustal components (Ca2+, HCO3, K+, and Mg2+), which accounted for 36% of the total ionic FWD. Finally, the marine components, Na+ and Cl, which primarily originated in the Arabian Gulf, with smaller contributions from the coast of the Red Sea in Jeddah and the central Black Sea region of Turkey, formed 17% of the total ionic FWD.

In this study, we developed a new method that combined satellite observations and cluster analysis of backward air mass trajectories to determine the potential source regions for the chemical components in Dammam’s rainwater. This technique was proven to be an effective tool for investigating the source receptor relationships of atmospheric air pollutants.


This study was supported by the Deanship of Scientific Research at Imam Abdulrahman Bin Faisal University with the grant numbers 2011-067-Eng and 2012-165-Eng. The authors gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and/or READY website ( used in this publication. The authors deeply acknowledge the Goddard Earth Sciences Data and Information Services Center (GES DISC) for the satellite data downloaded from EARTHDATA GIOVANNI website (


The authors declare no conflict of interest.

Don't forget to share this article 


Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

Latest coronavirus research from Aerosol and Air Quality Research

2018 Impact Factor: 2.735

5-Year Impact Factor: 2.827

SCImago Journal & Country Rank

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal, promotes submissions of high-quality research, and strives to be one of the leading aerosol and air quality open-access journals in the world.