Next Article in Journal
Evaluation of Utility of the Cement Solidification Process of Waste Ion Exchange Resin
Next Article in Special Issue
Seasonality of Aerosol Sources Calls for Distinct Air Quality Mitigation Strategies
Previous Article in Journal
A Rapid Method for Detecting Microplastics Based on Fluorescence Lifetime Imaging Technology (FLIM)
Previous Article in Special Issue
Personal Exposure to Black Carbon, Particulate Matter and Nitrogen Dioxide in the Paris Region Measured by Portable Sensors Worn by Volunteers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Source Apportionment of Atmospheric PM10 in Makkah Saudi Arabia by Modelling Its Ion and Trace Element Contents with Positive Matrix Factorization and Generalised Additive Model

1
Department of Environmental and Health Research, The Custodian of the Holy Two Mosques Institute for Hajj and Umrah Research, Umm Al Qura University, Makkah 24382, Saudi Arabia
2
Faculty of Environment, Institute for Transport Studies, University of Leeds, Leeds LS2 9JT, UK
*
Author to whom correspondence should be addressed.
Submission received: 29 January 2022 / Revised: 28 February 2022 / Accepted: 28 February 2022 / Published: 2 March 2022

Abstract

:
In this paper, the emission sources of PM10 are characterised by analysing its trace elements (TE) and ions contents. PM10 samples were collected for a year (2019–2020) at five sites and analysed. PM10 speciated data were analysed using graphical visualization, correlation analysis, generalised additive model (GAM), and positive matrix factorization (PMF). Annual average PM10 concentrations (µg/m3) were 304.68 ± 155.56 at Aziziyah, 219.59 ± 87.29 at Misfalah, 173.90 ± 103.08 at Abdeyah, 168.81 ± 82.50 at Askan, and 157.60 ± 80.10 at Sanaiyah in Makkah, which exceeded WHO (15 µg/m3), USEPA (50 µg/m3), and the Saudi Arabia national (80 µg/m3) annual air quality standards. A GAM model was developed using PM10 as a response and ions and TEs as predictors. Among the predictors Mg, Ca, Cr, Al, and Pb were highly significant (p < 0.01), Se, Cl, and NO2 were significant (p < 0.05), and PO4 and SO4 were significant (p < 0.1). The model showed R-squared (adj) 0.85 and deviance explained 88.1%. PMF identified four main emission sources of PM10 in Makkah: (1) Road traffic emissions (explained 51% variance); (2) Industrial emissions and mineral dust (explained 27.5% variance); (3) Restaurant and dwelling emissions (explained 13.6% variance); and (4) Fossil fuel combustion (explained 7.9% variance).

1. Introduction

Air pollution is one of the major public health concerns causing cardiovascular, pulmonary, respiratory, and cognitive problems. PM10 (particulate matter of aerodynamic diameter up to 10 μm) has been associated with millions of premature mortalities worldwide [1]. Air pollution is the cost of the rapid growth of urbanisation and industrialisation, and as a result, the air quality in large urban areas is getting worse day by day. Urban air quality has caught the interest of researchers worldwide, analysing emission sources, levels of various air pollutants, and predicting future concentrations and their impacts on human health [2,3,4,5]. In addition to the natural sources of air pollutants (e.g., dust storms, volcano eruptions, and forest fires), anthropogenic sources contribute a significant proportion of the pollutant loads to the atmosphere [6,7]. Therefore, numerous investigations have been made to quantify the contribution of each emission source (e.g., [8,9]).
Air pollution levels in the atmosphere are controlled by emission sources, meteorological conditions, and geographical characteristics. Previously, several studies have reported that particulate pollution is a serious environmental issue in Middle Eastern countries, including Saudi Arabia [10,11,12,13,14,15,16,17], which in addition to arid geographical conditions and frequent sand storms, is encouraged by rapid urbanisation and industrialisation. Saudi Arabia has large scale desert areas, arid or semi-arid climatic conditions, little rain, and strong winds, which play a positive role in enhancing the atmospheric burden of particulates. These natural conditions are supported by anthropogenic activities including busy roads, large and small scales industries, construction-and-demolition projects, and oil and coal burnings [11,18,19].
Air quality in Makkah has been analysed by several researchers in the past. Habeebullah et al. [20] carried out source apportionment of PM10 based on ionic analysis. In the study, SO42−, NO3, Ca2+, Na+, and Cl were the major contributors, contributing 37.81% to PM10 concentrations. Road traffic, mineral dust, industries, construction and demolition, sea spray, and marine aerosols were identified as the major emission sources of PM10 [20]. Farahat et al. [10] monitored air quality in Makkah, Madinah, and Jeddah during the Hajj seasons of 2019 and 2020 and compared the levels in the two years. They concluded that concentrations of PM10, NO2, and CO were reduced during the Hajj seasons in Makkah in 2020 compared to 2019 from 27.45 to 34.65 µg/m3, 9.95 to 4.50 ppbv, and 1.28 to 0.81 ppm, respectively. During the Hajj 2019 season, there were no restrictions, whereas in 2020, due to COVID-19, the number of pilgrims was restricted, therefore, the reductions seemed to be caused by the COVID-19 lockdown restriction [10]. Nayebare et al. [8] analysed PM2.5 samples to determine their emission sources in Makkah for a year. Average 24 h PM 2.5 concentration (35 µg/m3) was observed to be higher than the WHO standard (25 µg/m3). They reported that the Air Quality Index (AQI) was “unhealthy to hazardous.” According to this study, vehicular emission, industrial mixed dust, earth crust, and fossil fuel combustion were determined to be the major sources of PM2.5 [8].
The focus of this study is to measure the concentrations of several ions (Na+, Ca2+, Mg2+, F, Cl, Br, NO2, NO3, SO42−, and PO43−) and trace elements (Pb, As, Cd, Cr, Se, and Al) in PM10. This study builds on Habeebullah et al. [20], which analysed the concentrations of water-soluble ions with Principal Component Analysis (PCA) to identify the potential sources of PM10. In this study, the range of chemical species is increased, and advanced modelling approaches are employed, namely Generalised Additive Model (GAM), Positive Matrix Factorisation (PMF), and Enrichment Factor (EF), along with graphical visualization and correlation analysis to identify the main emission sources of PM10 and investigate the linear and nonlinear association between the various constituents of PM10.

2. Materials and Methods

2.1. Description of the Monitoring Sites

Makkah is the Holy City for all Muslims around the world. Makkah is situated in the western part of Saudi Arabia, among the Sarawat Mountains, about 80 km inward from the Red sea, and 277 m above sea level [10]. The population of Makkah city is around 1.96 million, and the population of the Makkah region is more than 8.5 million, with a growth rate of 1.8% [21,22]. Makkah is unique in the sense that it receives an additional 3 to 4 million pilgrims during the month of pilgrimage (Hajj) and Ramadhan (fasting) every year, except in 2020 and 2021, when Hajj was cancelled due to the COVID-19 pandemic [23]. The vast number of pilgrims leads to traffic congestion within and around the city, resulting in high levels of exhaust and non-exhaust emissions and dust resuspension [24]. Makkah experiences hot and dry climatic conditions, where the maximum temperature reaches 55 °C in the summer season [25].
The data analysed in this paper were collected at five sites around the city, namely Aziziyah, Sanaiyah, Abdeyah, Misfalah, and Askan (Figure 1). Aziziyah and Misfalah are considered urban traffic sites, Sanaiyah an industrial site, Askan a residential site, and Abdeyah a background site. The major emission sources in Makkah are related to road traffic, including both exhaust and non-exhaust emissions and the resuspension of dust particles on roadsides [20]. Makkah is not an industrial city like Jeddah; however, there is a power plant in Sanaiyah. Furthermore, emissions from restaurants and dwellings need to be quantified.

2.2. Collection of PM10 Samples

High-volume samplers were used to collect PM10 samples from 8 March 2020 to 9 March 2021 at the five sites in Makkah. These samplers used high–volume glass fibre filters (8 × 10 inches, grade G 653, Whatman, Maidstone, UK) with inlet collection efficiency of a cut–point of 9.7 microns and flow rate of 1.13 m3/min. Before deploying in the field, the filters were put in an oven (LDO–060E, Lab Tech, Daejeon, Korea) at 300 °C for 5 h to remove their moisture and organic contents [26,27,28,29]. Samples were collected for 24-h periods covering all days of the week. The filters of the samples were collected strictly from 9:00 to 10:00 am. Samples were not collected from 7 April 2020 to 31 May 2020 due to the COVID-19 lockdown in Makkah. After collection, the samples were transported to the laboratory at the Custodian of the Two Holy Mosques Institute for Hajj and Umrah Research, Umm Al-Qura University Makkah, where the filters were kept at room temperature for 24 h in a desiccator. An analytical balance was used to check the weight of the filters until a constant weight was obtained, which was considered the final weight of the filters. The filters were then stored in the laboratory, while sealed in a polythene bag until the time of analysis [26,27,28,29].

2.3. Analysis of Water-Soluble Ions and Metal Contents

Before analysis, the filters were weighed again. Then each filter was cut into four equal portions and stored at 4 °C in a refrigerator. Each piece was analysed for different constituents. PM10 concentrations were calculated using the mass of the collected samples on the filters (the difference in the mass of the filters before and after exposure in the air for 24 h), and the volume of the air passed through the filters [26].
The filters were analysed for various ions in the PM10 samples. One-quarter of the filter was analysed for water-soluble ions, namely fluoride (F), chloride (Cl), nitrite (NO2), bromide (Br), nitrate (NO3), sulphate (SO42−), phosphate (PO43−), sodium (Na+), calcium (Ca2+), and magnesium (Mg2+). For this purpose, the filter was shredded into 25 mL deionised distilled water with a resistivity of 18 Ωcm in a conical flask. The conical flask was ultrasonicated for 1 h in an ultrasonic bath (ATM40–28LCD, Ovan, Badalona, Spain). To remove undissolved particles, the solution was filtered through a 0.45 μm pore size membrane (CHROMAFIL, CA–45/25 (S), Macherey-Nagel, Düren, Germany) and stored in a refrigerator at 4 °C [30]. Ion Chromatography (850 Professional, Metrohm, Herisau, Switzerland) was used to determine the ion concentrations in the extract. Then, HNO3 (3.2 mM) and Na2CO3 (1.8 mM) were used to dissolve cations and anions, respectively, using a flow rate of 0.7 mL/min and injection volume of 10 µL [31]. For quality assurance, triplicate samples and blank calibrations were carried out. The detection limit (ppm) of the ion chromatography was 0.001 for fluoride and phosphate, 0.002 for sodium and magnesium, and 0.005 for chloride, nitrite, bromide, nitrate, and calcium.
For metal analysis, a portion of the filter paper was digested using the ‘hot acid extraction procedure’ (US EPA Method IO-3.1). The filter paper was shredded into strips in a 150 mL small flask, adding 10 mL of freshly prepared regal water (HCl:HNO3, 75 mL:25 mL). Filter strips were made to be set down in such a way that they were covered by the acid mixture entirely. The flask was placed on the hot plate at 60–70 °C for 60 min and was not allowed to dry. Once digestion was completed, the flask was removed from the hot plate and allowed to cool. The flask walls were rinsed with 20 mL deionized water and left for another 30 min, allowing acid to diffuse into the water from the sample filter. A PTFE syringe filter (Lichen Cottage Syringe filters, Marlborough, MA, USA) 0.45 µm in size was used to filter the digested material into a 50 mL volumetric flask and filled up to the mark [32].
The concentrations of heavy metals were determined with a graphite furnace using an atomic absorption spectrophotometer, Thermo Scientific (ice 3000 series, Waltham, MA, USA). A deuterium lamp was used by the instrument to make background corrections for every reading. Stock calibration standards for each analysed metal (Pb, As, Cd, Cr, Se, Al) were made by 1000 ppm (Sigma-Aldrich, Saint Louis, MO, USA) standard for respective metals. Calibration standards for each metal were made based on concentration to absorbance values.
For QA procedures, respective certified reference materials (CRM) were used the same way as the samples. NIST 1648a (National Institute of Standards and Technology, Gaithersburg, MD, USA), IC-6-1 (nsi lab solutions, Raleigh, NC, USA), and IC-7-2 (nsi lab solutions, USA) were used as the CRM for metals, cations, and anions, respectively. Recoveries show satisfying results for the analyzed concentrations [33,34].

2.4. General Statistical Analysis

The speciated data of PM10 was analysed using R programming language [35] and three of its packages: ‘openair’ (version 2.8-6) [36], ‘mgcv’ (version 1.8-39) [37], and ‘ggplot2’ (version 3.3.5) [38]. Correlation analysis was performed to investigate the linear relationship between different ions. A correlation plot was developed in the openair–package [36] using its function ‘corPlot’ (version 0.92). The ‘mgcv-package’ was used for the implementation of GAM, and PMF (version 5.0) was used to identify the main emission sources of PM10 in Makkah.

2.5. Positive Matrix Factorization

Positive matrix factorization (PMF) is one of the most widely employed models for the source apportionment of different air pollutants [39]. PMF (version 5.0) was employed for the source apportionment analysis to identify the major sources of PM10 in Makkah, which explained over 80% variation in PM10 concentrations. PMF has proven to be a powerful technique for particulate matter source apportionment, which was developed by [40]. For more details on PMF and receptor modelling, the readers are referred to [41]. PMF solves the chemical mass balance equation between measured pollutant concentrations and source profiles as presented in Equation (1) [42]:
X ij = k = 1 p g i k f k j + e i j
In Equation (1), ‘p’ is the number of factors, ‘f’ is the species profile of each source, ‘g’ is the amount of mass contributed by each factor to each individual sample, and ‘e’ is the error (residual) of each species.

2.6. Generalised Additive Model (GAM)

GAM is a type of supervised machine learning model initially developed by [43], which combines the properties of generalized linear models with additive models. GAM is an advanced nonlinear model, which does not assume linearity between the response and explanatory variables and does not require the modelled variable to be normally distributed, permitting the response probability distribution to be any member of the exponential family (e.g., normal, exponential, gamma, Poisson and many other) [44]. For a given response variable, Y regressed over ‘m’ explanatory variables X1, X2, ..., Xm, a GAM in a general form can be described as shown in Equation (2):
Y = s1 (X1) + s2 (X2) + … + sm (Xm)
where Y is the response (PM10 concentrations in Equation (3)) variable and ‘s’ is the smoothing term, which corresponds to an associated explanatory (independent or predictor) variable (X). Using the above general form, the model developed in this study can be presented as below (Equation (3)):
PM10 ~ s1 (Pb) + s2 (Cd) + s3 (Cr) + s4 (As) + s5 (Se) + s6 (Al) + s7 (F) + s8 (Cl) + s9 (NO2) + s10 (Br) + s11 (NO3) + s12 (PO4) + s13 (SO4) + s14 (Na) + s15 (Ca) + s16 (Mg)

2.7. Enrichment Factor

Enrichment factor (EF) has been widely used to calculate the extent of anthropogenic contribution in comparison to the earth crust. The main aim of EF is to determine the sources of particulate matter, here PM10. Generally, aluminium (Al) is used as a reference element due to its abundance in earth crust. Several researchers (e.g., [8,9,45,46]) have previously used EF to identify the two main sources (anthropogenic sources and crustal earth) of trace elements in PM10 and PM2.5.
To calculate EF the following formula was used:
Ratio 1 = Concentration of TE in PM10/Concentration of Al in PM10
Ratio 2 = Concentration of TE in earth-curst/Concentration of Al in earth crust
EF = Ratio 1/Ratio 2
In Equations (4)–(6), EF stands for enrichment factor, TE for trace elements, whose emission sources are to be determined, and Al (aluminium) is the reference element. The concentrations of TE and Al in PM10 were determined as a part of this project, whereas the concentrations of TE and Al in earth crust were taken from [47]. An EF value of a TE less than 10 shows that a significant proportion of the TE is emitted from the earth’s crust, whereas an EF value of greater than 10 shows that the TE is mainly emitted by anthropogenic sources [8,9].

3. Results and Discussion

3.1. Descriptive Analysis of PM10 and Its Constituents

The levels of PM10 (µg/m3) at different monitoring sites are depicted in Figure 2, which showed that annual average PM10 concentrations exceeded the Saudi Arabia national annual air quality standards of 80 µg/m3 [48]. Annual average air quality standards developed by WHO (15 µg/m3) and USEPA (50 µg/m3) are even lower [49]; therefore, the WHO and EPA standards were exceeded as well. The annual average of PM10 concentrations (µg/m3) were 304.68 ± 155.56 at Aziziyah, 219.59 ± 87.29 at Misfalah, 173.90 ± 103.08 at Abdeyah, 168.81 ± 82.50 at Askan, and 157.60 ± 80.10 at Sanaiyah, in descending order. Average concentrations of PM10 and its major and minor constituents are depicted in Figure 2, and the time series of the monthly average is depicted in Figure 3. Aziziyah experienced the highest concentrations of PM10, ion contents, and trace elements, followed by the Misfalah site. These two sites are among the busiest urban areas in Makkah and remain very busy during the month of Hajj (Pilgrimage) and Ramadhan (fasting). However, Makkah also remains busy during other times of the year. It should be noted that in 2020 Hajj was cancelled due to the COVID-19 pandemic lockdown. Figure 3 showed that the levels of PM10 and its constituents were higher in summer (June to September), which could be due to seasonal effects, as generally, summer experiences the most number of dust storms in Saudi Arabia [19]. Missing data in May was due to the COVID-19 pandemic lockdown in Makkah.
High PM10 levels may pose a serious risk to the health of residents and visitors. Previously, several authors have reported high levels of particulate matter in Makkah e.g., [23]. Several reasons have been reported for the high levels of particulate matter in Makkah and the surrounding areas. Farahat [50] reported that large-scale infrastructure activities, overusing governmental subsidized energy, water desalination, heavy traffic in urban areas, and cement plants were the main reasons for the high levels of pollutants. Saudi Arabia is situated within a semi-arid region, and a major part of the country is occupied by deserts. Strong wind in the summer season uplifts desert soil and sand into the atmosphere, which causes heavy aerosol loading [18]. Due to large desert areas, frequent sandstorms, and little rainfall, the desert areas and the urban areas of Saudi Arabia experience higher concentrations of particulate matter [9,50,51].

3.2. Correlation Analysis

Correlation analysis shows the linear relationship between two variables. The strength of the relationship is expressed in terms of correlation coefficients (r), which range from −1 to +1, where zero shows no linear relationship, and ±1 shows a hundred percent positive/negative linear relationship.
The R-values between PM10 and its major and minor constituents are presented in the shape of correlation plots in Figure 4. The strength of the R-values is defined by the shape of the circle, its colour, and the number written inside. The R-values between the different chemical species varied from site to site (Figure 4). To provide an overall picture of the correlation, we combined the data from all five sites and calculated the R-values (Figure 4). Combining the data from all sites, we increased the power of the statistical test by increasing the number of observations of each chemical species [8]. PM10 showed the strongest correlation with Mg (0.83), followed by Ca (0.77), and the weakest correlation with Al (0.03), followed by As (−0.07). The strongest correlation was found between Na and Ca (0.91), followed by Ca and PO4 (0.87), considering all chemical species. Most of the pollutants had a positive association with each other, except As, which was negatively associated with all pollutants, except Cd and Cr (Figure 4). However, the correlation varied from site to site (not shown for brevity), probably indicating variations in their emission sources and factors responsible for their dispersion, dependent on the local micro-level environmental and infrastructure conditions. Generally, if two chemical species are emitted by the same source, they will have a high positive correlation coefficient. Otherwise, they will show a weak or even negative correlation coefficient. Furthermore, those pollutants emitted by the dominant sources that strongly affect PM10 concentrations (e.g., road traffic), will have stronger correlation coefficients with PM10. Ion concentrations dominate the PM10 concentrations and have a stronger effect on PM10. Therefore, they have stronger correlation coefficients than metals. However, linear correlation simply shows a relationship between a pair of pollutants and does not take into account the interaction between various pollutants. In addition, the relationship between different pollutants is not always linear (as shown in Section 3.3). Therefore, further advanced analysis is provided in the coming sections.

3.3. Generalised Additive Model (GAM)

As discussed in Section 3.2, the strength of the correlation coefficient is based on the linear association between two chemical species. However, the relationship between the two species is not always linear. Therefore, in this section, a nonlinear machine learning approach was employed to provide further insight into the association between PM10 and its constituents and see how much variation in PM10 concentrations could be explained by the predictors analysed here.
PM10 was used as the modelled (response or dependent) variable, and all constituents (TEs and ions) were used as predictor (explanatory or independent) variables. Among the predictors, Mg, Ca, Cr, Al, and Pb were highly significant (p-value < 0.01, ‘**’), Se, Cl, and NO2 were significant (p < 0.05, ‘*’), and PO4 and SO4 were significant (p < 0.1, ‘.’). The model showed an adjusted r-squared of 0.85, and deviance explained 88.1%. Deviance is a measure of goodness-of-fit of a model, similar to the R2 of Gaussian data. Several other statistical metrics were calculated by comparing observed and predicted PM10 concentrations, including correlation coefficients (0.94), a factor of two (0.99), mean biased error (1.19e−9), and root mean squared error (8.89 µg/m3). Definitions and formulae for calculating these statistical metrics can be found in [52,53]. Predicted and observed PM10 concentrations are compared in Figure 5 in the form of a scatter plot, which showed a strong association between predicted and measured concentrations of PM10.
Figure 6 shows the nonlinear association between PM10 and some of the predictors. Their mutual association varied with the levels of the predictor variables. For example, in the case of Al, there was a positive association between Al and PM10; however, when Al level reached about 1.8 µg/m3, the nature of association changed from positive to negative, which meant PM10 levels decreased with increasing Al. When Al levels reached about 2.8, the association changed again, and the curve became parallel with the x-axis (horizontal), meaning PM10 levels did not change with further increasing Al levels approximately from 2.8–4.0 µg/m3. However, it should be noted that at higher levels of Al, there were fewer data points (shown by the vertical black ticks on the x-axis), and therefore the uncertainties levels increased (shown by the dotted lines on both sides of the solid line). Therefore, at higher levels, the association between the two variables became weaker and less reliable. For more detail on the nonlinear association between PM10 and the predictors, see Figure 6. The nonlinear relationship between PM10 and the predictor variables was probably due to changes in emission sources and meteorological conditions that controlled the levels of PM10 and the predictor variables.

3.4. Enrichment Factors (EF)

Enrichment factors (EF) were calculated for the TE analysed in this paper (Table 1). Concentrations of TE in PM10 and earth-crust are given in Table 1. The ratios of TE to Al in PM10 were divided by their ratios in earth-crust to get the values of EF (Table 1). As mentioned above, all TEs having EF higher than 10 (EF > 10), were considered to be predominantly emitted by anthropogenic sources, whereas TE having EF less than 10 (EF < 10) were considered to have originated from the earth-crust. According to the data presented in Table 1, only two TE had EF higher than 10, which were Na and Ca. All the rest of TE had EF less than 10. Therefore, EF analysis showed that most of the TE originated from the earth-crust, except Na and Ca, which were attributed to anthropogenic sources. The shortcoming of the EF analysis here was that the elemental composition of local earth-crust was not available, which might have affected the outcome of the analysis. Khodeir et al. [9] and Nayebare et al. [8] have also used EF in their analysis for sources identification in Jeddah and Makkah, respectively. According to [8], the anthropogenic TE were Cu, Zn, Eu, Cl, Pb, S, Br, and Lu, whereas the TE derived from the earth-crust were Al, Si, Na, Mg, Rb, K, Zr, Ti, Fe, Mn, Sr, Y, Cr, Ga, Ca, Ni, and Ce.

3.5. Positive Matrix Factorization (PMF)

The outcomes of PMF are presented in Figure 7 and Table 2. Four factors were identified and ascribed as the four major sources of PM10 in Makkah. Factor fingerprints, factor profiles, and factor contributions to PM10 are presented in Figure 7 and Table 2. Factor 1 (industrial emissions and mineral dust) explained 27.5%, factor 2 (fossil fuel combustion) explained 7.9%, factor 3 (Road traffic) explained 51.0%, and factor 4 (restaurants and dwellings) explained 13.6% of the variations in PM10 concentrations (Figure 7). Looking at the fingerprints of different factors, it can be observed that most of the elements contributed to more than one factor, meaning they were emitted by more than one emission source in different proportions, except As and NO2, which contributed to only a single factor each, i.e., factor 4 and 3, respectively. Se and Cd contributed to only two factors each, i.e., factors 1 and 4 and 2 and 3, respectively. This showed that different emission sources contributed to the emission of a chemical species in various proportions. Such analysis helped us identify the main sources of emission, as explained below. The four factors of PMF were identified as four emission sources of PM10 in Makkah based on the loadings of chemical species in each factor.
Road traffic (exhaust and non-exhaust emission and resuspension of dust): Factor 3 explained the highest proportion (51%) of the PM10 concentration. Factor 3 was identified as road traffic emission because the pollutants are related to vehicle exhaust, non-exhaust, or resuspension of dust on roadsides [8]. Factor 3 mainly consisted of NO2, F, Cl, Br, PO4, NO3, SO4, Na, Ca, and Mg (Figure 7). Minor contributors were Cd, Cr, and Pb. It was reported that Ca, Mg, Na, and Al were some of the major elements emitted by the vehicle’s tyres [8], whereas Pb and Cd were found in the emission from brake wear [54,55].
Industrial emissions and mineral dust: Factor 1 explained the second largest proportion (27.5%) in PM10 concentrations. The chemical species that contributed to factor 1 were Se, Al, Pb, Cr, Mg, Ca, SO4, Br, Na, NO3, PO4, Cr, and F. The three main contributors were Se, Pb, and Al. About 98% of Se falls in factor 1, whereas the loadings of Al and Pb were about 40% and 35%, respectively. Forty percent (40%) of Se is reported to be emitted by anthropogenic sources, including coal combustion, metal smelting, and biomass burning [56]. These sources also emitted Al and Pb. Therefore, Factor 1 was identified as industrial emission. However, Na, Ca, Mg, Al, and PO4 are the markers for mineral dust [57]. Part of the mineral dust may also come from construction and demolition activities. Regional transport of marine aerosols from the Red Sea also contributed to the concentrations of Na, Cl, PO4, and Mg [8].
Fossil fuel combustion: Factor 2 explained a 7.9% variation in PM10 concentrations and consisted of Cd, Pb, Al, SO4, F, Cr, Br, NO3, Na, Ca, and Mg. A large proportion of factor 2 was made of Cd, Pb, and Al. Airborne Cd is predominantly (85–90%) emitted by anthropogenic sources, including fossil fuel combustion, the smelting and refinement of nonferrous metals, and municipal waste incineration [58]. Naturally, Cd is emitted by volcanic eruptions [58]. In addition to road traffic, Pb is also emitted during the open burning of municipal waste [59]. The presence of Al, F, Na, Ca, and Mg also indicated the effect of windblown dust particles, which is common in Makkah due to high temperature, dry weather, and frequent dust storms.
Restaurants and dwellings: Factor 4 explained 13.6% of the variance in PM10 concentrations and showed high loadings of As and Cr. Other elements that made a minor contribution to factor 4 were Al, Mg, Ca, Na, SO4, NO3, Br, Cr, F, and Se. ‘As’ is naturally emitted by volcanic eruptions, vegetation, and windblown dust, whereas anthropogenic sources of ‘As’ include fuel combustion especially burning of low-grade brown coal and metal smelting [58]. Cr is mainly emitted to the atmosphere by burning coal, oil, and natural gas [60]. In Makkah, coal is mainly burnt in restaurants, whereas natural gas is burnt in restaurants and homes. Combustion of coal and natural gas also contributes to the combustion-related pollutants, which include SO2 and NOx, which in turn are transformed to NO3 and SO4 ions in the atmosphere.
According to the enrichment factor (EF) analysis, most of the TEs were originated from the earth’s crust, except Na and Ca. However, according to the PMF model, the dominant emission sources, in addition to windblown dust and resuspension of dust on roadsides, in Makkah were road traffic, various industries, combustion of fossil fuels, and combustion of low-grade coal, oil, and natural gas in restaurants and dwellings. Previously, the authors of [8] identified four main sources of PM2.5, namely, vehicular emissions (30.1%), industrial-mixed dust (28.9%), soil/earth crust (24.7%), and fossil fuels/oil combustion (16.3%), which generally were in agreement with the findings of the current study. However, in this study, the contribution of restaurants and dwellings was also identified as a significant source. Makkah, due to the presence of the Holy Mosque, attracts millions of visitors every year. These visitors use the local restaurants for their food that burn coal and natural gas for cooking and thus emit a significant amount of pollutants to the atmosphere. In Makkah, several restaurants on each street burn coals and natural gas, especially in the centre of the city around the Holy Mosque. It is important to mention that large-scale construction-and-demolition activities in Makkah also add to the atmospheric load of particulate matter [53]. Previously, Na and Cl have been identified as the markers of sea spray and marine aerosols; however, in this study, their contribution to each factor is minor, probably indicating that these pollutants are emitted by combustion processes and are generated from the earth-crust.

4. Conclusions and Recommendation

In this paper, the main aim was to identify the major emission sources of PM10 in Makkah, employing various statistical and machine learning approaches, which included: (a) Descriptive statistics and graphical visualization, (b) Correlation analysis, (c) Enrichment factor, (d) Generalised additive model, and (e) Positive matrix factorization. Makkah is one of the densely populated cities in Saudi Arabia, where atmospheric particulate matter is the pollutant of concern. PM10 concentrations exceeded WHO (15 µg/m3), USEPA (50 µg/m3), and the Saudi Arabia national (80 µg/m3) annual air quality standards. For the first time, a GAM model was developed using PM10 as a response variable and ions and TEs as predictors. GAM analysed the nonlinear association between PM10 and its predictors, which varied at different levels of the predictors, and showed that these predictors were able to explain about 88% variation in PM10.
PMF identified four main emission sources of PM10 in Makkah:
  • Road traffic including both exhaust and non-exhaust emission, and resuspension of dust (explained 51% variance in PM10);
  • Industrial emissions and mineral dust (explained 27.5% variance);
  • Restaurants and dwellings (explained 13.6% variance); and
  • Other fossil fuel combustion (explained 7.9%).
In addition to the natural sources of PM10 such as windblown dust, several anthropogenic sources were identified that contributed significantly to PM10 concentrations. Although it is not feasible to control emissions from natural sources, government policy can be implemented to control emissions from anthropogenic sources, including emissions from road traffic, industries, and restaurants and dwellings. It is interesting to see that burning coals and natural gas in restaurants, hotels, and dwellings were identified as a significant source of particulate matter. However, like other urban areas, in Makkah, the main emission source identified was road traffics, which emit exhaust and non-exhaust emissions and cause resuspension of dust particles on roadsides. The main weakness of the study is that a limited number of ions and TEs were analysed. However, GAM analysis showed that these species explained nearly 90% variations in PM10. It is recommended that a detailed emission inventory be compiled in Makkah characterising road traffic flow, fleet composition, point emissions (e.g., industries, power plants), and area emissions (e.g., emission from restaurants, hotels, and dwellings). Furthermore, the impacts of air pollution should be analysed, including its impact on health, economy, and sustainability.
To manage air quality and cut the emission of pollutants from various sources, the following actions are recommended in Makkah:
  • Road traffic is the main emission source of gaseous and particle pollution in Makkah. Therefore, it is important to effectively manage traffic flow to avoid congestion during the peak hours, particularly during the months of the Hajj and Ramadhan [61]. Also, the use of public transport should be encouraged;
  • Further intervention for reducing air pollution in Makkah may include banning old polluting vehicles and retrofitting them with new technology [61];
  • More electric vehicles should be included in the fleet and better charging facilities should be provided [62];
  • To improve driving behaviour, further training should be provided to discourage abrupt acceleration, deceleration and idling [61]. Idling is a serious issue due to extremely high temperatures. Drivers keep their vehicle’s engines on during summer to keep the interior of the vehicle cold using the vehicle’s air conditioners. Otherwise, within a few minutes, the vehicle becomes unbearably hot. Environmentally friendly vehicles and technologies must be used in such situations;
  • As previously reported by several authors [63], trees not only control air pollution but also help moderate air temperature. Therefore, more trees should be grown in the City of Makkah, particularly on roadsides;
  • Large-scale construction and demolition activities increase the loading of atmospheric dust, which could be reduced by an effective water spray programme [64].

Author Contributions

Idea initiation, T.M.H., S.M., J.Z. and E.A.M.; PM10 samples collection, T.M.H., J.Z. and E.A.M.; sample lab analysis, J.Z.; statistical analysis, T.M.H. and S.M.; writing the first draft, S.M., T.M.H. and J.Z.; visualisation, S.M. and T.M.H.; review, T.M.H., S.M., J.Z. and E.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the King Abdulaziz City of Science and Technology (KACTS) (Research project number 14–ENV2582–10).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are thankful to the King Abdulaziz City of Science and Technology (KACTS) for funding this research project (14–ENV2582–10). The authors are also thankful to the Custodian of the Holy Two Mosques Institute for Hajj and Umrah Research and the Scientific Research Institute at Umm Al–Qura University for their support and assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organisation (WHO). Ambient (Outdoor) Air Pollution. 2018. Available online: https://www.who.int/news-room/fact-sheets/detail/ambient-(outdoor)-air-quality-and-health (accessed on 27 January 2022).
  2. Ramírez, O.; de la Campa, A.S.; Amato, F.; Catacolí, R.A.; Rojas, N.Y.; de la Rosa, J. Chemical composition and source apportionment of PM10 at an urban background site in a high–altitude Latin American megacity (Bogota, Colombia). Environ. Pollut. 2018, 233, 142–155. [Google Scholar] [CrossRef]
  3. Du, Y.; Xu, X.; Chu, M.; Guo, Y.; Wang, J. Air particulate matter and cardiovascular disease: The epidemiological, biomedical and clinical evidence. J. Thorac. Dis. 2016, 8, E8–E19. [Google Scholar] [CrossRef] [PubMed]
  4. Cowie, G.; Lawson, W.; Kim, N.D. Australian dust causing respiratory disease admissions in some North Island, New Zealand hospitals. N. Z. Med. J. 2010, 123. [Google Scholar]
  5. Weuve, J.; Puett, R.C.; Schwartz, J.; Yanosky, J.D.; Laden, F.; Grodstein, F. Exposure to Particulate Air Pollution and Cognitive Decline in Older Women. Arch. Intern. Med. 2012, 172, 219–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Baklanov, A.; Molina, L.T.; Gauss, M. Megacities, air quality and climate. Atmos. Environ. 2016, 126, 235–249. [Google Scholar] [CrossRef]
  7. Ryou, H.G.; Heo, J.; Kim, S.-Y. Source apportionment of PM10 and PM2.5 air pollution, and possible impacts of study characteristics in South Korea. Environ. Pollut. 2018, 240, 963–972. [Google Scholar] [CrossRef]
  8. Nayebare, S.R.; Aburizaiza, O.S.; Siddique, A.; Carpenter, D.O.; Hussain, M.M.; Zeb, J.; Aburiziza, A.J.; Khwaja, H.A. Ambient air quality in the holy city of Makkah: A source apportionment with elemental enrichment factors (EFs) and factor analysis (PMF). Environ. Pollut. 2018, 243, 1791–1801. [Google Scholar] [CrossRef]
  9. Khodeir, M.; Shamy, M.; Alghamdi, M.; Zhong, M.; Sun, H.; Costa, M.; Chen, L.-C.; Maciejczyk, P. Source apportionment and elemental composition of PM2.5 and PM10 in Jeddah City, Saudi Arabia. Atmos. Pollut. Res. 2012, 3, 331–340. [Google Scholar] [CrossRef] [Green Version]
  10. Farahat, A.; Chauhan, A.; Al Otaibi, M.; Singh, R.P. Air Quality Over Major Cities of Saudi Arabia During Hajj Periods of 2019 and 2020. Earth Syst. Environ. 2021, 5, 101–114. [Google Scholar] [CrossRef]
  11. Lim, C.C.; Thurston, G.D.; Shamy, M.; Alghamdi, M.; Khoder, M.; Mohorjy, A.M.; Alkhalaf, A.K.; Brocato, J.; Chen, L.C.; Costa, M.; et al. Temporal variations of fine and coarse particulate matter sources in Jeddah, Saudi Arabia. J. Air Waste Manag. Assoc. 2018, 68, 123–138. [Google Scholar] [CrossRef] [Green Version]
  12. Abbas, H.; Ghanem, S.; Abahussain, A. Assessing PM2.5 in Bahrain from the DPSIR Framework Perspective. In Proceedings of the 2021 Third International Sustainability and Resilience Conference: Climate Change, online, 15–16 November 2021; pp. 53–60. [Google Scholar]
  13. Javed, W.; Guo, B. Chemical characterization and source apportionment of fine and coarse atmospheric particulate matter in Doha, Qatar. Atmos. Pollut. Res. 2020, 12, 122–136. [Google Scholar] [CrossRef]
  14. Phanikumar, D.; Basha, G.; Ratnam, M.V.; Kondapalli, N.K.; Ouarda, T.; Pangaluru, K.; Kumar, K.N.; Kishore, P. Assessment of particulate matter concentration and gaseous pollutants in urban and rural regions over the Emirate of Abu Dhabi, UAE. J. Atmos. Sol.—Terr. Phys. 2020, 199, 105217. [Google Scholar] [CrossRef]
  15. Al-Thani, H.; Koc, M.; Fountoukis, C.; Isaifan, R.J. Evaluation of particulate matter emissions from non-passenger diesel vehicles in Qatar. J. Air Waste Manag. Assoc. 2020, 70, 228–242. [Google Scholar] [CrossRef] [PubMed]
  16. Alattar, N.; Yousif, J. Evaluating Particulate Matter (PM2.5 and PM10) Impact on Human Health in Oman Based on a Hybrid Artificial Neural Network and Mathematical Models. In Proceedings of the 2019 International Conference on Control, Artificial Intelligence, Robotics & Optimization (ICCAIRO), Majorca Island, Spain, 3–5 May 2019; pp. 129–135. [Google Scholar]
  17. Coskuner, G.; Jassim, M.S.; Munir, S. Characterizing temporal variability of PM2.5/PM10 ratio and its relationship with meteorological parameters in Bahrain. Environ. Forensics 2018, 19, 315–326. [Google Scholar] [CrossRef]
  18. Alam, K.; Trautmann, T.; Blaschke, T.; Subhan, F. Changes in aerosol optical properties due to dust storms in the Middle East and Southwest Asia. Remote Sens. Environ. 2014, 143, 216–227. [Google Scholar] [CrossRef]
  19. Farahat, A.; El-Askary, H.; Al-Shaibani, A. Study of Aerosols’ Characteristics and Dynamics over the Kingdom of Saudi Arabia Using a Multisensor Approach Combined with Ground Observations. Adv. Meteorol. 2015, 2015, 247531. [Google Scholar] [CrossRef] [Green Version]
  20. Habeebullah, T.M.; Munir, S.; Zeb, J.; Morsy, E.A. Analysis and Sources Identification of Atmospheric PM10 and Its Cation and Anion Contents in Makkah, Saudi Arabia. Atmoshere 2022, 13, 87. [Google Scholar] [CrossRef]
  21. Division, U.P. The World’s Cities in 2018: Data Booklet; United Nation: New York, NY, USA, 2018; ISBN 9789211483062. [Google Scholar]
  22. CDSI Central Department of Statistics and Information of Saudi Arabia. 2020. Available online: https://www.stats.gov.sa/en (accessed on 19 September 2021).
  23. Seroji, A. Particulates in the Atmosphere of Makkah and Mina Valley During Ramadan and Hajj Seasons of 1424 and 1425 H (2004–2005). Arab Gulf J. Sci. Res. 2008, 26, 199–206. [Google Scholar] [CrossRef] [Green Version]
  24. Habeebullah, T.M. An Analysis of Air Pollution in Makkah—A View Point of Source Identification. EnvironmentAsia 2013, 6, 11–17. [Google Scholar] [CrossRef]
  25. Almazroui, M.; Islam, M.N.; Athar, H.; Jones, P.D.; Rahman, M.A. Recent climate change in the Arabian Peninsula: Annual rainfall and temperature analysis of Saudi Arabia for 1978–2009. Int. J. Clim. 2012, 32, 953–966. [Google Scholar] [CrossRef]
  26. CARB. C.A.R.B. Standard Operating Procedure for the Determination of PM10 Mass by Gravimetric Analysis; Northern Laboratory Branch, Monitoring and Laboratory Division: Documents, CA, USA, 2018; Available online: https://fdocuments.in/document/table-of-contents-2020-6-30-determination-of-pm10-mass-by-gravimetric-analysis.html?page=1 (accessed on 12 January 2022).
  27. Ferreira, R.; Nunes, C.; Souza, M.; Canela, M. Multivariate Optimization of Extraction Variables of PAH in Particulate Matter (PM10) in Indoor/Outdoor Air at Campos dos Goytacazes, Brazil. J. Braz. Chem. Soc. 2021, 32, 618–625. [Google Scholar] [CrossRef]
  28. Waked, A.; Favez, O.; Alleman, L.Y.; Piot, C.; Petit, J.-E.; Delaunay, T.; Verlinden, E.; Golly, B.; Besombes, J.-L.; Jaffrezo, J.-L.; et al. Source apportionment of PM10 in a north-western Europe regional urban background site (Lens, France) using positive matrix factorization and including primary biogenic emissions. Atmos. Chem. Phys. 2014, 14, 3325–3346. [Google Scholar] [CrossRef] [Green Version]
  29. Rezaei, S.; Naddafi, K.; Hassanvand, M.S.; Nabizadeh, R.; Yunesian, M.; Ghanbarian, M.; Atafar, Z.; Faraji, M.; Nazmara, S.; Mahmoudi, B.; et al. Physiochemical characteristics and oxidative potential of ambient air particulate matter (PM10) during dust and non-dust storm events: A case study in Tehran, Iran. J. Environ. Health Sci. Eng. 2018, 16, 147–158. [Google Scholar] [CrossRef]
  30. Liu, H.; Zheng, J.; Qu, C.; Zhang, J.; Wang, Y.; Zhan, C.; Yao, R.; Cao, J. Characteristics and Source Analysis of Water-Soluble Inorganic Ions in PM10 in a Typical Mining City, Central China. Atmosphere 2017, 8, 74. [Google Scholar] [CrossRef] [Green Version]
  31. Sowlat, M.H.; Naddafi, K.; Yunesian, M.; Jackson, P.L.; Lotfi, S.; Shahsavani, A. PM10 Source Apportionment in Ahvaz, Iran, Using Positive Matrix Factorization | Enhanced Reader. Celan—Soil Air Water 2013, 41, 1143–1151. [Google Scholar] [CrossRef]
  32. EPA, U. IO Compendium Method IO-3.1: Compendium of Methods for the Determination of Inorganic Compunds in Ambient Air: Selection, Preparation and Extraction of Filter Material; EPA/625/R-; Centre for Environmental Research Information, Office of Research and Development, US EPA: Cinncinati, OH, USA, 1999. [Google Scholar]
  33. Ali, N.; Eqani, S.A.M.A.S.; Nazar, E.; Alhakamy, N.A.; Rashid, M.I.; Shahzad, K.; Zeb, J.; Shen, H.; Ismail, I.M.I.; Albar, H.M.S.A. Arsenic and lead in the indoor residential settings of different socio-economic status; assessment of human health risk via dust exposure. Environ. Sci. Pollut. Res. 2020, 28, 13288–13299. [Google Scholar] [CrossRef] [PubMed]
  34. Kadi, M.W.; Ismail, I.; Ali, N.; Shaltout, A.A. Spectroscopic Assessment of Platinum Group Elements of PM10 Particles Sampled in Three Different Areas in Jeddah, Saudi Arabia. Int. J. Environ. Res. Public Health 2020, 17, 3339. [Google Scholar] [CrossRef]
  35. Team, R.C.; R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: http://www.r-project.org/index.html (accessed on 20 April 2021).
  36. Carslaw, D. The Openair Manual: Open-Source Tools for Analysing Air Pollution Data. Available online: https://davidcarslaw.github.io/openair/ (accessed on 12 November 2021).
  37. Wood, S. The Mgcv Package: Mixed GAM Computation Vehicle with Automatic Smoothness Estimation. 2020. Available online: https://cran.r-project.org/web/packages/mgcv/index.html (accessed on 23 January 2022).
  38. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis, 2nd ed.; Springer: Cham, Switzerland, 2016; ISBN 9783319242774. [Google Scholar]
  39. Sun, X.; Wang, H.; Guo, Z.; Lu, P.; Song, F.; Liu, L.; Liu, J.; Rose, N.; Wang, F. Positive matrix factorization on source apportionment for typical pollutants in different environmental media: A review. Environ. Sci. Process. Impacts 2020, 22, 239–255. [Google Scholar] [CrossRef]
  40. Paatero, P. Least squares formulation of robust non-negative factor analysis. Chemom. Intell. Lab. Syst. 1997, 37, 23–35. [Google Scholar] [CrossRef]
  41. Hwang, I. Introduction and Application of the PMF Model to Estimate the Source Apportionment of PM2.5 at Various Sites. J. Inst. Eng. 1970, 8, 25–31. [Google Scholar] [CrossRef]
  42. EPA. EPA Positive Matrix Factorization (PMF) 5.0 Fundamentals and User Guide; EPA/600/R-.; U.S. Environmental Protection Agency, Office of Research and Development: Washington, DC, USA, 2014. [Google Scholar]
  43. Hastie, T.J. Generalized Additive Models, 1st ed.; Chapman and Hall: London, UK; CRC: London, UK, 1990; ISBN 0-412-34390-8. [Google Scholar]
  44. Wood, S.N. Generalized Additive Models: An Introduction with R; Chapman & Hall: Boca Raton, FL, USA, 2006; ISBN 9781584884743. [Google Scholar]
  45. Fabretti, J.-F.; Sauret, N.; Gal, J.-F.; Maria, P.-C.; Schärer, U. Elemental characterization and source identification of PM2.5 using Positive Matrix Factorization: The Malraux road tunnel, Nice, France. Atmos. Res. 2009, 94, 320–329. [Google Scholar] [CrossRef]
  46. Aprile, F.M.; Bouvy, M. Distribution and enrichment of heavy metals at the Tapacurá River basin, northeastern Brazil. Braz. J. Aquat. Sci. Technol. 2008, 12, 1–8. [Google Scholar] [CrossRef] [Green Version]
  47. Taylor, S.R. Abundance of chemical elements in the continental crust: A new table. Geochim. Cosmochim. Acta 1964, 28, 1273–1285. [Google Scholar] [CrossRef]
  48. ArRiyadh Air Quality. Air Quality Standards, Legislations and Policy. Available online: http://aq.arriyadhenv.gov.sa/aq/about-us/riyadh-air-quality-standards/ (accessed on 8 January 2022).
  49. Al-Taani, A.A.; Nazzal, Y.; Howari, F.M.; Yousef, A. Long-term trends in ambient fine particulate matter from 1980 to 2016 in United Arab Emirates. Environ. Monit. Assess. 2019, 191, 143. [Google Scholar] [CrossRef] [PubMed]
  50. Farahat, A.; El-Askary, H.; Dogan, A.U. Aerosols Size Distribution Characteristics and Role of Precipitation during Dust Storm Formation over Saudi Arabia. Aerosol Air Qual. Res. 2016, 16, 2523–2534. [Google Scholar] [CrossRef] [Green Version]
  51. Farahat, A. Comparative analysis of MODIS, MISR, and AERONET climatology over the Middle East and North Africa. Ann. Geophys. 2019, 37, 49–64. [Google Scholar] [CrossRef] [Green Version]
  52. Carslaw, D. Defra Regional and Transboundary Model Evaluation Analysis—Phase 1; Version 15; UK AIR: London, UK, 2011; Available online: https://uk-air.defra.gov.uk/library/reports?report_id=653 (accessed on 24 December 2021).
  53. Munir, S.; Habeebullah, T.M.; Seroji, A.; Morsy, E.A.; Mohammed, A.M.; Abu Saud, W.; Abdou, A.E.; Awad, A.H. Modeling Particulate Matter Concentrations in Makkah, Applying a Statistical Modeling Approach. Aerosol Air Qual. Res. 2013, 13, 901–910. [Google Scholar] [CrossRef] [Green Version]
  54. Grigoratos, T.; Martini, G. Brake wear particle emissions: A review. Environ. Sci. Pollut. Res. Int. 2015, 22, 2491–2504. [Google Scholar] [CrossRef] [Green Version]
  55. Kennedy, P.; Gadd, J. Preliminary Examination of Trace Elements in Tyres, Brake Pads and Road Bitumen in New Zealand; Infrastructure Auckland: Auckland, New Zealand, 2003. [Google Scholar]
  56. Feinberg, A.; Maliki, M.; Stenke, A.; Sudret, B.; Peter, T.; Winkel, L.H.E. Mapping the drivers of uncertainty in atmospheric selenium deposition with global sensitivity analysis. Atmos. Chem. Phys. 2020, 20, 1363–1390. [Google Scholar] [CrossRef]
  57. Chen, L.; Peng, C.; Gu, W.; Fu, H.; Jian, X.; Zhang, H.; Zhang, G.; Zhu, J.; Wang, X.; Tang, M. On mineral dust aerosol hygroscopicity. Atmos. Chem. Phys. 2020, 20, 13611–13626. [Google Scholar] [CrossRef]
  58. WHO. WHO Air Quality Guidelines for Europe, 2nd ed.; WHO: Geneva, Switzerland, 2000; Available online: https://www.euro.who.int/en/health-topics/environment-and-health/air-quality/publications/pre2009/who-air-quality-guidelines-for-europe,-2nd-edition,-2000-cd-rom-version. (accessed on 1 October 2021).
  59. Emsley, J. Nature’s Building Blocks: Everything You Need to Know about the Elements, 2nd ed.; Oxford University Press: Oxford, UK, 2011; ISBN 9780199605637. [Google Scholar]
  60. Cheng, H.; Zhou, T.; Li, Q.; Lu, L.; Lin, C. Anthropogenic Chromium Emissions in China from 1990 to 2009. PLoS ONE 2014, 9, e87753. [Google Scholar] [CrossRef] [PubMed]
  61. Simpson, I.J.; Aburizaiza, O.S.; Siddique, A.; Barletta, B.; Blake, N.J.; Gartner, A.; Khwaja, H.; Meinardi, S.; Zeb, J.; Blake, D.R. Air Quality in Mecca and Surrounding Holy Places in Saudi Arabia During Hajj: Initial Survey. Environ. Sci. Technol. 2014, 48, 8529–8537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Barisione, M. Electric Vehicles and Air Pollution: The Claims and the Facts. Available online: https://epha.org/electric-vehicles-and-air-pollution-the-claims-and-the-facts/ (accessed on 25 January 2022).
  63. Munir, S.; Habeebullah, T.M.; Morsy, E.A.; Mohammed, M.F. The Effect of Tree Plantations on Moderating Air Temperature in Arafat, Makkah for the Period 2002-2016. J. King Abdulaziz Univ. Environ. Arid Land Agric. Sci. 2018, 27, 63–72. [Google Scholar] [CrossRef]
  64. Munir, S.; Habeebullah, T.M.; Seroji, A.R.; Gabr, S.S.; Mohammed, A.M.; Morsy, E.A. Quantifying temporal trends of atmospheric pollutants in Makkah (1997–2012). Atmos. Environ. 2013, 77, 647–655. [Google Scholar] [CrossRef]
Figure 1. Location map of the five sampling sites (Aziziyah, Misfalah, Sanaiyah, Abdeyah, and Askan) in Makkah [20].
Figure 1. Location map of the five sampling sites (Aziziyah, Misfalah, Sanaiyah, Abdeyah, and Askan) in Makkah [20].
Toxics 10 00119 g001
Figure 2. The concentrations of PM10 (a), major constituents (b), and minor constituents (c) at different monitoring sites in Makkah.
Figure 2. The concentrations of PM10 (a), major constituents (b), and minor constituents (c) at different monitoring sites in Makkah.
Toxics 10 00119 g002aToxics 10 00119 g002b
Figure 3. Time variations of monthly average PM10(upper-panel), major constituents (middle-panel), and minor constituents (lower-panel) during the study period.
Figure 3. Time variations of monthly average PM10(upper-panel), major constituents (middle-panel), and minor constituents (lower-panel) during the study period.
Toxics 10 00119 g003aToxics 10 00119 g003b
Figure 4. Correlation plots of PM10, trace elements, and ions for ‘all sites’ in Makkah (individual sites are not shown for brevity).
Figure 4. Correlation plots of PM10, trace elements, and ions for ‘all sites’ in Makkah (individual sites are not shown for brevity).
Toxics 10 00119 g004
Figure 5. Scatter plot between measured and GAM predicted PM10 concentrations (µg/m3). The 1:1 line is solid black and the 1:0.5 and 1:2 lines are dashed black. Together these lines help show how close a group of points are to a 1:1 relationship and show the points that are within a factor of two (FAC2). The solid blue line shows the linear model line with 95% confidence intervals. The equation of the line and R-square value is shown in the top-left of the panel.
Figure 5. Scatter plot between measured and GAM predicted PM10 concentrations (µg/m3). The 1:1 line is solid black and the 1:0.5 and 1:2 lines are dashed black. Together these lines help show how close a group of points are to a 1:1 relationship and show the points that are within a factor of two (FAC2). The solid blue line shows the linear model line with 95% confidence intervals. The equation of the line and R-square value is shown in the top-left of the panel.
Toxics 10 00119 g005
Figure 6. Showing the nonlinear association of PM10 with the predictors (only four predictors i.e. Mg, PO4, Ca and Al are shown for brevity). The dashed lines show the confident intervals and the vertical lines on the x-axis show the data presence.
Figure 6. Showing the nonlinear association of PM10 with the predictors (only four predictors i.e. Mg, PO4, Ca and Al are shown for brevity). The dashed lines show the confident intervals and the vertical lines on the x-axis show the data presence.
Toxics 10 00119 g006
Figure 7. Results of the positive matrix factorization (PMF) run 3 (the convergent run), showing profiles % of species for factor 1 (a), factor 2 (b), factor 3 (c) and factor 4 (d).
Figure 7. Results of the positive matrix factorization (PMF) run 3 (the convergent run), showing profiles % of species for factor 1 (a), factor 2 (b), factor 3 (c) and factor 4 (d).
Toxics 10 00119 g007aToxics 10 00119 g007b
Table 1. Enrichment factors (EF) of various chemical species.
Table 1. Enrichment factors (EF) of various chemical species.
ElementEarth-Crust CxPM10 Cx(Cx/CAl) PM10(Cx/CAl) Earth-CrustEF
Al8.2301.7841.0001.0001.000
Pb12.5000.1610.0901.5190.059
Cd0.2000.0850.0480.0241.961
Cr100.0000.0170.01012.1510.001
As1.8000.0850.0480.2190.218
Se0.0500.0040.0020.0060.369
F525.0000.3300.18563.7910.003
Cl130.0006.8373.83215.7960.243
Br2.5000.0420.0240.3040.078
Na2.3608.5964.8180.28716.803
Ca4.1509.3415.2360.50410.384
Mg2.3300.4440.2490.2830.879
Note: The earth-crust concentrations were taken from [47]. Cx is the concentration of an element compared to ‘Al’ (CAl).
Table 2. Results of the positive matrix factorization run 3 (the convergent run), showing profile % of species, % of factor total and concentrations of species in each factor (factor 1 to factor 4).
Table 2. Results of the positive matrix factorization run 3 (the convergent run), showing profile % of species, % of factor total and concentrations of species in each factor (factor 1 to factor 4).
Profile (% of Species)Profile (% of Factor Total)Profile (Concentration of Species)
Factor 1Factor 2Factor 3Factor 4Factor 1Factor 2Factor 3Factor 4Factor 1Factor 2Factor 3Factor 4
PM1027.487.8751.0113.6478.5954.5671.5282.8356.6616.23105.1628.11
Pb34.0459.226.74<0.010.080.330.01<0.010.060.100.01<0.01
Cd<0.0193.196.390.42<0.010.28<0.01<0.01<0.010.080.01<0.01
Cr23.235.814.8866.080.01<0.01<0.010.03<0.01<0.01<0.010.01
As<0.00<0.01<0.01100.00<0.01<0.01<0.010.25<0.01<0.01<0.010.08
Se97.28<0.01<0.012.72<0.01<0.01<0.01<0.01<0.01<0.01<0.01<0.01
Al36.7451.03<0.0112.220.913.060.000.640.660.91<0.010.22
F6.287.4982.643.590.030.090.190.040.020.030.280.01
Cl14.5112.3862.5410.561.342.782.842.080.970.834.170.70
NO2<0.01<0.01100.00<0.01<0.01<0.010.01<0.01<0.01<0.010.01<0.01
Br24.1415.2557.363.250.010.020.02<0.010.010.010.02<0.01
NO312.2515.7165.406.632.156.675.622.471.551.988.260.84
PO45.165.8888.96<0.010.020.050.16<0.010.010.020.24<0.01
SO426.0119.9645.618.4113.0824.3211.248.999.437.2416.533.05
Na9.1912.6173.424.781.093.634.271.210.791.086.290.41
Ca19.5913.1062.534.782.534.103.961.311.831.225.830.45
Mg27.028.0953.1111.770.160.120.160.150.120.040.230.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Habeebullah, T.M.; Munir, S.; Zeb, J.; Morsy, E.A. Source Apportionment of Atmospheric PM10 in Makkah Saudi Arabia by Modelling Its Ion and Trace Element Contents with Positive Matrix Factorization and Generalised Additive Model. Toxics 2022, 10, 119. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10030119

AMA Style

Habeebullah TM, Munir S, Zeb J, Morsy EA. Source Apportionment of Atmospheric PM10 in Makkah Saudi Arabia by Modelling Its Ion and Trace Element Contents with Positive Matrix Factorization and Generalised Additive Model. Toxics. 2022; 10(3):119. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10030119

Chicago/Turabian Style

Habeebullah, Turki M., Said Munir, Jahan Zeb, and Essam A. Morsy. 2022. "Source Apportionment of Atmospheric PM10 in Makkah Saudi Arabia by Modelling Its Ion and Trace Element Contents with Positive Matrix Factorization and Generalised Additive Model" Toxics 10, no. 3: 119. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10030119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop