# Source Apportionment of Ambient Particulate Matter (PM) in Two Western African Urban Sites (Dakar in Senegal and Bamako in Mali)

Thierno Doumbia, Catherine Liousse, Marie-Roumy Ouafo-Leumbe, Seydi Ababacar Ndiaye, Eric Gardrat, Corinne Galy-Lacaux, Cyril Zouiten,

Véronique Yoboué, Claire Granier

► To cite this version:

Thierno Doumbia, Catherine Liousse, Marie-Roumy Ouafo-Leumbe, Seydi Ababacar Ndiaye, Eric Gardrat, et al.. Source Apportionment of Ambient Particulate Matter (PM) in Two Western African Urban Sites (Dakar in Senegal and Bamako in Mali). Atmosphere, 2023, 14 (4), pp.684. 10.3390/atmos14040684. [https://arxiv.org/abs/2304.04684](https://arxiv.org/abs/2304.04684)

HAL Id: hal-04235213

https://cnrs.hal.science/hal-04235213v1

Submitted on 10 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC0 1.0 - Universal - International License
<figure><img src="image_2.png" /></figure>

Type of the Paper (Article)

1

Source Apportionment of Ambient Particulate Matter (PM) in two Western Africa Urban Sites (Dakar in Senegal and Bamako in Mali)

2

3

Thierno Doumbia¹*, Catherine Liousse¹, Marie-Roumy Ouaf-Leumbe², Seydi Ababacar Ndiaye³, Eric Gardrat¹, Corinne Galy-Lacaux¹, Cyril Zouiten⁴, Véronique Yoboué⁵, and Claire Granier¹,⁶

4

5

¹ Laboratoire d'Aérologie (LAERO), Université de Toulouse, CNRS/UPS, Toulouse, France; thiernodoumbia@yahoo.fr / thierno.doumbia@aero.obs-mip.fr

6

7

² Université de Douala, Douala, Cameroun

8

³ Laboratoire de Physique de l'Atmosphère et de l'Océan-Simeon Fongang (LPAO-SF), Université de Dakar, Sénégal

9

10

⁴ Géosciences Environnement Toulouse (GET), Université de Toulouse, Toulouse, France

11

⁵ Laboratoire des Sciences de la Matière, de l'Environnement et de l'énergie Solaire, Université Félix Houphouët-Boigny, Abidjan, Côte d'Ivoire

12

13

⁶ NOAA Chemical Sciences Laboratory and CIRES/University of Colorado, Boulder, CO, USA

14

* Correspondence: thiernodoumbia@yahoo.fr; Tel.: +33 5 61 33 27 66

15

**Abstract:** Following population growth and rapid urbanization, West African cities become major sources of anthropogenic pollution. Additionally, Saharan dust has a significant impact, representing a potentially toxic mix of sources for the population. This study characterizes atmospheric composition and its sources in two African capitals, Bamako, Mali and Dakar, Senegal. TSP, PM₁₀ and PM₂.₅ samples were collected during the dry season in 2009 when pollution levels were high: chemical analysis included organic carbon (OC), elemental carbon (EC), ions, and metals. PM₂.₅ and PM₁₀ concentrations were 5-10 times and 3-8 times higher, respectively, than the 2005 WHO 24-hour standards. Using PCA and PMF methodologies, five sources were identified in each city. In Bamako, traffic (motor vehicles and resuspended road dust) was the prevailing source of PM₂.₅ and PM₁₀, accounting for 47% and 45%, respectively. Crustal dust was the second most important source (24-30%), followed by solid fuel combustion (16-13%) and secondary aerosols (10-16%). In Dakar, the following sources of PM₂.₅ and PM₁₀ are identified: traffic (49%), mineral dust (16-25%), sea salts (15-20%) and industries (10-11%). Our study provides crucial information about the historical change of source characteristics in these two African cities, which can help for future mitigation strategies.

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

**Citation:** To be added by editorial staff during production.

32

Academic Editor: Firstname Last-name

**Keywords:** Air quality; particulate matter; source receptor models; PCA; PMF; West Africa.

33

34

Received: date

1. Introduction

35

Revised: date

Accepted: date

Published: date

The role of atmospheric aerosols in climate change has long been recognized through their interactions with solar radiation and influence on cloud formation. However, as reported in the IPCC reports [1], significant uncertainties about the effect of aerosols on climate persist. On the other hand, several epidemiological and toxicological studies have demonstrated the impact of particulate matter (PM) on human health [2] [3] [4]. Exposure to PM can adversely affect human health via inhalation according to their size and chemical composition, especially in urban environments. For example, aerosol particles with aerodynamic diameters of less than 2.5 µm (PM₂.₅) are recognized as being more likely to cause toxicological effects due to their ability to easily penetrate into human lungs [5] [6]. In

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Submitted for possible open access
publication under the terms and
conditions of the Creative Commons
Attribution (CC BY) license
(https://creativecommons.org/licenses/by/4.0/).
2010, ambient air pollution, primarily PM<sub>2.5</sub>, was responsible for 3.3 million premature deaths worldwide [7]. The chemical composition of PM, such as organic compounds and metals, contributes to their high toxicity [8] [9]. The study conducted in African urban areas by Val et al. [9] on the pro-inflammatory response of different particle sizes (ultrafine, fine, and coarse fractions) highlighted the role of ultrafine particles primarily from combustion sources. The identification and quantitative evaluation of the contributions of different aerosol sources is important for understanding their implications in health and environmental effects and for defining mitigation policies. In developed countries, source characterization of atmospheric pollution has been investigated for several years. This is not the case in developing countries with densely populated regions, such as the West Africa capitals, where very few studies have been conducted. This region is heavily influenced by Saharan dust aerosols, which contribute significantly to total aerosol loads, especially during the dry season, when North-Easterly trade winds (Harmattan) are predominant [10] [11] [12]. Deboudt et al. [13] demonstrated that in February 2006, desert dust dominated the coarse fraction and organic matter was the major constituent of the fine particles performed in Mbour, a coastal site in Senegal. In the urban areas, a large proportion of atmospheric particles is attributable to anthropogenic sources such as traffic (vehicular, two-wheel two-strokes and resuspended road particles), biomass and residential combustion (wood, charcoal and fuel), waste burning, and industries [14] [15] [16] [17] [18] [19] [20]. Indeed, mechanical processes such as resuspension of road particles remain a large source due to tire abrasion, brake wear and the high number of unpaved roads in West African cities. All of these emission sources combined result in West African PM pollution levels comparable to those observed in European, Asian, and United States megacities, and well above the World Health Organization (WHO) guidelines [19] [21]. Due to a lack of long-term PM concentration data, Doumbia et al. [19] estimated PM<sub>2.5</sub> from a one-year measurement of black carbon in Dakar (Senegal) in 2008 and found an annual average PM<sub>2.5</sub> concentration that was four times higher than the WHO recommendation of 10 µg/m<sup>3</sup>.

Consequently, rapid deterioration of air quality has now become a heavy challenge for the majority of capital cities in West Africa: in order to define policies to reduce particulate matter pollution at source level, there is an urgent need for a better quantification of the sources of this atmospheric pollution. Receptor-oriented source apportionment models are often used to identify the number of aerosol sources and to estimate their contributions to measured concentrations [22] [23] [24]. Principal component analysis (PCA) [25] [26], edge analysis (UNMIX) [27], chemical mass balance (CMB) [28], and positive matrix factorization (PMF) [29] are some popular multivariate models for determining the number of sources and elemental relative contributions. Although UNMIX and CMB methods have similar goals to PMF, they operate in different ways, with UNMIX not always able to resolve as many source factors as PMF [30]. The user must provide source profiles for mass apportionment when using CMB. The comparison of these different approaches shows that higher measurement uncertainty requirements in PMF produce better results than the other two models [31]. In our study, PCA is performed in conjunction with PMF to investigate and quantify the aerosol sources in West Africa urban sites. These two statistic models are considered here due to their complementarity, accessibility and relatively easy use. The PCA approach requires minimal inputs from the characteristics of the different sources, with an aim at quantifying information both on source profiles and magnitude of absolute concentrations. The PCA is used as an explanatory tool to identify major sources of air pollutants, whereas the PMF is used to quantify the contributions of identified sources to all measured pollutants. The PMF is a non-data sensitive technique that can manage and resolve inhomogeneous datasets without any previous univariate analysis. It does not require to know historical emission source profiles and introduces estimated error associated to each measured component. We also used enrichment factor
(EF) analyses in addition to PCA and PMF to categorize air pollution sources in our sites. 97
EF compares the ratios for various elements in measured atmospheric concentrations to 98
the corresponding ratios in geological material [28] [32]. 99

A dataset from the POLCA (POLlution des Capitales Africaines: African Capital Pollu-
tion), a jointly program between African and French universities, is used in this study [9]
[19] [33]. The goal was to characterize atmospheric particulate pollution and its health
impact at traffic sites in Bamako (Mali) and Dakar (Senegal), to handle different specifici-
ties typical of west-African megacities. Indeed, the measurement sites differ in terms of
vehicle fleet, fuel type, road quality, residential combustions, dust episodes, and biomass
burning impacts. Also, to study the upper limit of air quality impacts, samples were col-
lected during intensive campaigns occurring in the selected sites during the dry season
(January 2009 in Bamako and December 2009 in Dakar). The dry season is a time period
that is well known for exhibiting high levels of pollution in such areas. Samples were an-
alyzed for a comprehensive description of the physico-chemical characteristics of the aer-
osols, as well as their relative contributions to the aerosol particles with aerodynamic di-
ameters of less than 2.5 µm (PM2.5), aerosol particles with aerodynamic diameters of less
than 10 µm (PM10) and Total Suspended Particles (TSP). 101
102
103
104
105
106
107
108
109
110
111
112
113
114
The focus of this paper is on the characterization of aerosol chemical composition and
quantification of urban sources in Bamako and Dakar using PMF statistic model in asso-
ciation with PCA analysis. We describe the methodology in Section 2 and discuss the re-
sults in Section 3. 115
116
117
118

## 2. Materials and Methods 119

2.1. Description of monitoring sites 120

PM samples were collected during the dry season in Bamako (Mali) and Dakar (Sen-
egal), two West African cities. Bamako, the largest city of Mali, is located on the Niger
River in the southwestern part of the country (Figure 1) and had a population of approx-
imately 1.8 million inhabitants during the campaign in 2009 (12.5% of the national popu-
lation) [34]. Measurements were taken at a downtown crossroad (12°39′N, 8°04′W) be-
tween January 20th and February 1st 2009; a dust event episode occurred between January
22nd and January 24th, 2009 (Figure 2). In Dakar, a coastal city in west Senegal (Figure 1)
with a population of about 2.5 million people in 2009 (21% of the national population)
[35], measurements were performed between December 1st and 13th 2009, near a busy
crossroad (14°40′N, 17°25′W). Because of their proximity to markets, both sampling sites
are heavily influenced by intensive human activities, but with very different characteris-
tics important for air quality. 121
122
123
124
125
126
127
128
129
130
131
132
133
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Map of measurement sites with superimposed wind roses during the two field campaigns in Dakar (top) and Bamako (down). 134</figcaption></figure>

## 2.2. PM sample collection 137

Daily 24-hour integrated aerosol samples were collected by filtration at the two west Africa urban sites using automatic aerosol samplers designed and developed for the IN-DAAF (International Network to study Deposition and Atmospheric chemistry in Africa, https://indaafobs-mip.fr/) program [36]. The system is designed to collect samples simultaneously using six lines: each particle mode (Total Suspended Particle TSP, Particle Matters with diameters less than 10 µm PM₁₀ and Particle Matters with diameters less than 2.5 µm PM₂.₅) was sampled separately with an individual pump operated at a flow rate of 18±7 l.min⁻¹. 138 139 140 141 142 143 144 145 146

## 2.3. Environmental parameters 147

At each sampling site, an automated meteorological station collected wind speed and direction, temperature (T), and relative humidity (RH). The coefficients of light scattering (σscat) and aerosol absorption (σabs) were also measured using a nephelometer (520 nm wavelength, model M903, Ecotech) and aethalometers (model AE-42 with seven wavelengths (370 to 950 nm) and 880 nm wavelength, model AE-16, Magee Scientific). Figures 2-3 display the diurnal variations of all of these parameters for the duration of the campaigns. There are distinct diurnal patterns in σscat and σabs for both sites over the entire measurement period, with two daily peaks: a morning peak between 7 and 9 a.m. and an evening peak between 5 and 7 p.m., with a second evening peak at 9 p.m. in Bamako. This 148 149 150 151 152 153 154 155 156
diurnal variability is primarily caused by traffic and the evolution of the boundary layer during the day [19] [37], which is common in urban areas.

In Bamako (Figure 2), $\sigma_{scat}$ and $\sigma_{abs}$ generally show evening peaks (9 p.m.) more pronounced than in Dakar, coinciding with a thin boundary layer resulting from low ambient temperature and low westerly wind speed, except for the January 22-23 period, during which the diurnal variation is associated with high light scattering, low aerosol absorption, and high easterly wind speed (8 m.s⁻¹). The large difference observed between $\sigma_{scat}$ and $\sigma_{abs}$ is due to the dust episode that occurred on January 22-23, 2009.

In Dakar (Figure 3), $\sigma_{scat}$ and $\sigma_{abs}$ display large morning and evening peaks coinciding with morning and evening traffic peaks [19]. The lower values were observed late at night when traffic volume and air temperatures were low. Dakar had a lower $\sigma_{scat}$ and a higher $\sigma_{abs}$ than Bamako, indicating that the characteristics of the two cities differed. To support this analysis, a $\sigma_{abs}$ to $\sigma_{scat}$ ratio was calculated and plotted for Bamako and Dakar to characterize traffic or dust influence [38] (Figure S1, Supplement). The mean ratio obtained in Dakar is 2.55±0.99, which is four times higher than the values obtained in Bamako (0.70±0.33), indicating the importance of anthropogenic sources (e.g. traffic emissions) in Dakar, whereas the values in Bamako indicate an area influenced more by dust or resuspended dust particles.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Light scattering and absorption coefficients, as well as meteorological data collected in Bamako during the dry season of 2009 (January 19 to February 2).</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Light scattering and absorption coefficients, as well as meteorological data collected in Dakar during the dry season of 2009 (November 30 to December 14). 179</figcaption></figure>

## 2.4. Chemical analyses 183

Aerosol particles were collected on 47 mm quartz filters for carbonaceous analyses 184 and Teflon filters for water soluble ions and metal elements analyses. A calibrated micro-balance (Model Mettler MC21S) was used to weigh Teflon filters before and after aerosol sampling for determining the PM gravimetric masses. The samples were conditioned for 24-hour at relative ambient humidity of 30±15% prior to weighing and chemical analysis. 185 186 187 188 TSP, PM₁₀ and PM₂.₅ were analyzed for carbonaceous (elemental carbon EC, organic carbon OC and total carbon TC), water-soluble inorganic compounds and metal elements. 189 190 It should be noted that chemical analyses were performed a few days after sample collections. 191 192

### Water soluble inorganic compounds 194

Half of each Teflon filter was used in the analysis of inorganic ions (Na⁺, NH₄⁺, K⁺, Mg²⁺, Ca²⁺, SO₄²⁻, NO₃⁻, and Cl⁻), using ion chromatographic (IC) analyzer. These measurements 195 196 were conducted following the analytical protocol described in Adon et al. [39]. Briefly, cations are analyzed using Dionex ion chromatography (DX-100) and anions with a Dionex DX-500. Water extraction of the filter samples was performed by a 10 min-long sonication in plastic vials with 30 ml deionised ultra-pure water. All analyses were conducted immediately after extraction. 197 198 199 200

The quality control and quality assurance of the ion chemical analysis were ensured by measuring the blanks, the detection limit (DL), testing the reproducibility, and determining the measurement uncertainties. The DL is defined as three times the standard deviation (SD) of the procedure blank. The filter blank values and calculated DL for water soluble inorganic species are provided in Table S1 (Supplement). Only Mg²⁺ blank is lower than the DL, with a value of 1 ppb found in the 24 blanks tested. 201 202 203 204 205 206

The overall uncertainty of the measurements in this work is primarily due to three factors: sampling uncertainties, including possible gas-particle interactions and particle bounce-207 208 209
off; the relationship between blank values, detection limits and concentrations measured, and random contamination issues. The uncertainty associated with each ion was calculated using the following equation (1) from Alleman et al. [40]:

$$Unc = C\sqrt{u_r^2(m) + u_r^2(V) + \frac{u_c^2}{c^2}} \quad (1)$$

where C is the concentration of the element in the sample, u<sub>r</sub>(m) is the relative uncertainty of the mass m, u<sub>r</sub>(V) is the relative uncertainty of the air volume V, u<sub>c</sub> is the absolute uncertainty due to the contamination. Using the equation (1), the estimated uncertainties were 4% for SO<sub>4</sub><sup>2-</sup>, 18% for NH<sub>4</sub><sup>+</sup>, 44% for NO<sub>3</sub><sup>-</sup>, 14% for Cl<sup>-</sup>, 25% for Ca<sup>2+</sup>, and 1% for Mg<sup>2+</sup> and K<sup>+</sup>. NO<sub>3</sub><sup>-</sup> concentrations have a high degree of uncertainty due to evaporation issues [41]. When NO<sub>3</sub><sup>-</sup> is excluded, the average analytical uncertainty is estimated to be around 10%. A few NH<sub>4</sub><sup>+</sup> data were obtained in Bamako, most likely due to the same volatilization issue as with nitrate. As a result, NH<sub>4</sub><sup>+</sup> concentrations will not be shown here.

* Metal elements
The metal elements were determined by coupling microwave assisted digestion with Induced Coupled Plasma - Mass Spectrometry (ICP-MS). In this study, a special pre-treatment for loaded aerosols was set up following the Celo et al. [42] methodology. First, the samples, placed in 100 ml Teflon bombs, were subject to a microwave-assisted digestion within a mixture of ultra-high purity acids (10 ml 16N HNO<sub>3</sub>, 0.5 ml 28N HF). A closed vessel microwave assisted reaction system (MARS 5, CEM Corporation, Matthews, NC) was used, with a two-step digestion program. During the first step, the temperature was ramped within 15 min to 160 °C. After a 10 min dwell time, the temperature was ramped to 180 °C and samples were digested at this temperature for 30 min. The digests were then transferred to clean polypropylene tubes and diluted to 10 ml with HNO<sub>3</sub> 0.37N. Samples were then analyzed using ultra-sensitive ICP-MS (7500ce Element Technologies). The measurement quality of metal elements was ensured by controlling the blank and detection limit values. The concentrations of metals were determined in two types of blanks: vessel blanks (VB) and filter blanks (FB). The concentrations of trace metals in the vessel blank were used to evaluate the cleanliness of the microwave digestion and sample treatment processes. The VB was composed of acid solutions resembling those used in sample digestions. Metal concentrations in the FB were measured to establish a baseline from which field samples element concentrations were subtracted. It was made of VB and the same amount of filter paper that was used to collect aerosols in the field. The average concentrations in VB and FB, as well as the detection limits (DL), are shown in Table S2 (Supplement). The background concentrations of all elements are found to be higher in the filter blank than in the vessel blank. Most minor elements have DLs of less than 1 µg.l<sup>-1</sup>, whereas major elements (Ca, Na, K, Fe, Al, and Mg) have DLs ranging from 30 to 221 µg.l<sup>-1</sup>. Because our samples contain high concentration levels, the high sensitivity of the ICP-MS instrument accounts for the high variability. The high DL values do not pose a problem for atmospheric samples, because the concentrations of these elements are typically much higher. Field blank values for elements such as Na, K, Mg, Fe, Al, Ti, Zn, Mn, Pb, Be, and Ni are slightly higher than filter blanks (Figure S2, Supplement), indicating a small field contamination of field blanks for these elements. Standard reference materials (NIST SRM 1648 Urban Particulate Matter) prepared in the same manner as the digested samples, spikes, and duplicates were also analyzed to ensure the precision and accuracy of the analysis. Nine replicate NIST standard reference materials were measured. These materials were used directly without any treatment, except for acids added to match their matrix against samples analyzed in the same sequence. The reported recovery values for certified and non-certified elements are shown in Table S3 (Supplement). Except for Cr and Rb, which have recoveries of 60 and 69%, the digestion
method provides good recoveries (84-105%, with standard deviation SD range 2.7-6.5%) for SMRs referenced elements. Several authors reported comparable low Cr recovery percentages [43] [44]. This could be due to the amount of HF used during mineralization, with more HF increasing recovery.

Over 40 elements were determined, only a part of them being retained for the sake of interpretation after subtraction of blank filter concentrations. Uncertainties associated with each element, as required by the PMF statistical model, are determined through the analysis of blank filters and the intercomparison of analytical techniques. The estimated uncertainties in Dakar are less than 15%, with the highest values for Ni and Pb, and do not exceed 20% in Bamako. These uncertainties are acceptable for use in source-receptor models such as positive matrix factorization (PMF).

Finally, the concentrations of calcium, potassium, sodium, and magnesium in PM<sub>2.5</sub>, PM<sub>10</sub>, and TSP measured by ICP-MS (Inductively Coupled Plasma-Mass Spectrometry) and IC (Ion Chromatography) (see previous paragraph) were compared to ensure the quality of our data and to assess the ratio of soluble to non-soluble components. ICP-MS measures total elemental concentrations (Ca, K, Na, and Mg), whereas IC measures water-soluble fraction (Ca<sup>2+</sup>, K<sup>+</sup>, Na<sup>2+</sup>, and Mg<sup>2+</sup>). Linear regressions for elements such as calcium and sodium show slopes ranging from 0.9 to 1.3 in Dakar and from 1.9 to 2.1 in Bamako (Figure S3, Supplement), indicating that the majority of inorganic Ca and Na are soluble in Dakar (by 10-30%). Calcium and sodium concentrations measured by IC in Bamako are twice as low as those measured by ICP-MS, indicating that inorganic aerosols are less soluble in Bamako than in Dakar. This could be attributed to the strong influence of natural sources in Bamako, such as resuspended dust and Saharan Desert dust. Total potassium (K) concentrations in both sites are approximately 3 to 4 times that of soluble potassium (K<sup>+</sup>), and this value reached 10 during the dust episode in Bamako. This is significant because soluble potassium (K<sup>+</sup>) is a key indicator for wood combustion [45]. The increasing K/K<sup>+</sup> ratio during the dust event suggests that a fraction of the non-soluble potassium came from crustal dust. The same results are obtained for Mg/Mg<sup>2+</sup> ratios (about 3-4). Mg<sup>2+</sup> can be used as a tracer for marine aerosols. The increasing Mg/Mg<sup>2+</sup> ratio (from 3 to 16) during the dust event suggests that some of the non-soluble magnesium is also derived from crustal material [77].

* Carbonaceous elements 294

Quartz filters used for carbonaceous analysis were preheated for 24-hour at 600°C to lower carbon blanks before being exposed. 24-hour samples (PM<sub>2.5</sub>, PM<sub>10</sub> and TSP) were analyzed using a two-step combustion technique: total carbon (TC) and elemental carbon (EC) concentrations were determined using THERMAL (carbon analyzer model G4 ICARUS CS TF, Bruker Axes, Germany) method developed by Cachier et al. [46] and Thermal-Optical Reflectance (TOR, Desert Research Institute analyzer) method following the IMPROVE temperature protocol [47], as indicated in Table S4 (Supplement). Organic carbon (OC) concentrations were finally calculated as the difference between TC and EC. In the THERMAL technique, carbonates were removed (decarbonated or decarbonation) under HCl fumes prior to analyses due to carbonates interfering with carbon measurements. Sensitivity studies on the impact of decarbonation on Bamako and Dakar samples analyzed using TOR method were carried out. Decarbonation procedure is responsible of a decrease of OC concentration (range 40-60%) and an increase of EC (range 10-23%), depending on the site and size particle fraction (Tables S5 and S6, Supplement).

For analysis, lower concentrations of species in aerosol samples necessitate low and stable blank values. The DL was calculated using the mean value of blanks and their standard deviation (SD). The blank value must be small because the magnitude of SD normally increases with the magnitude of the blank value. The quality control measurements on EC and OC determinations performed in 2009 at the Bamako and Dakar sites are summarized in Table S7 (Supplement). The average DRI instrument DL for TC, EC, and OC was 0.45, 0.13 and 0.33 µg.cm⁻², respectively. These values are consistent with what the
manufacturer was reported. The obtained field blanks DL, on the other hand, are approximately 20 times higher for OC and 10 times higher for EC, indicating high field contamination.

Because different analytical methods can result in significant differences in OC and EC concentrations [48] [49], an intercomparison of the THERMAL and TOR methods with no-decarbonatation was also performed as a data quality test. TC values obtained by the TOR method were close (difference less than 10%) to those obtained by the THERMAL method in both Bamako and Dakar samples, with a slope close to unity (Figure S4, Supplement). These findings corroborate those reported in the literature. When it comes to EC and OC, there are differences between the two methods. OC TOR was approximately 0.94 times the value of OC THERMAL in Dakar samples, compared to 0.85 in Bamako samples. This means that, despite differences in charring correction and temperature protocol in the TOR method, the two methods are comparable in terms of OC analysis. The analytical method used influences EC values, which may be due to variable source contributions. The EC concentrations determined by the TOR method are higher than EC obtained by the THERMAL method. The EC TOR/EC THERMAL ratio in Dakar samples was 1.03, while it was 1.29 in Bamako samples. According to these findings, the differences between these two analytical methods are more pronounced in Bamako than in Dakar. This can be explained by three factors:

- The first one is susceptible to the presence of components such as brown carbon in carbonaceous aerosol [50], as brown carbon has the potential to influence splitting between OC and EC due to light absorption and medium thermal reactivity. Indeed, some brown carbon can be classified as OC, while the remainder is classified as EC [51]. Brown carbon was discovered to be emitted primarily by incomplete combustion, such as domestic fires, with this source heavily influencing the Bamako site [9]. Hitzenberger et al. [49] discovered comparable EC concentrations from diesel traffic sources using different methods. This finding is consistent with the fact that diesel emissions have a significant influence on Dakar samples (EC TOR to EC THERMAL ratio of 1.03 was found).

- The second factor is the mixing of aerosols. From Bamako, which is heavily influenced by dust, to Dakar, which is heavily influenced by traffic diesel, the EC TOR to EC THERMAL ratio decreases. The presence of dust particles in Bamako samples was assumed to interfere with thermal optical measurements, contributing to the observed large EC concentration variations.

- The third and final hypothesis is that the EC THERMAL is more sensitive to combustion aerosols originating from fossil fuels and biofuels [52], because these sources are more abundant in Bamako than in Dakar.

The estimated uncertainty is calculated using the following equation given by Chow et al. [53]:

$$Unc = \sqrt{(CV \times \overline{TC})^2 + DL^2} \quad (2)$$

where CV is the coefficient of variance for the replicate analysis, DL is the measured detection limit and $\overline{TC}$ is the average filter blanks. The uncertainty on measured TC with the TOR method was less than 10%.

As part of the THERMAL method quality control, filter blanks and instrument calibration were analyzed for two groups of samples: One for low carbon mass (< 25 µgC) and one for high carbon level (> 25 µgC). Before each analysis cycle, the instrument was calibrated by manually injecting sucrose into the analyzer. Table S7 (Supplement) shows good agreement between measured and target sucrose masses for both carbon
calibrations, with a ratio close to unity indicating excellent instrument calibration. When compared to the DRI analyzer (TOR method), the Carbon analyzer (THERMAL method) has lower detection limits. The THERMAL method produced an uncertainty of 5% using equation (2). The uncertainty in EC and OC varies by site and ranges from 5 to 30%. Based on all these results, the data measured using the THERMAL technique are considered in the following analysis.

2.5. Methods for identifying and quantifying air pollution sources 376

2.5.1. Principal Component Analysis (PCA) 377

The PCA technique was used to identify potential sources of air pollution in Dakar and Bamako. The goal of PCA in this case is to reduce the number of variables to a smaller set of factors or sources while retaining as much information as possible from the original dataset [26]. This method helps identify groups of correlated chemical compounds that come from the same air pollution sources. Due to the high variability of elemental concentrations, the chemical composition data are first transformed into a dimensionless standardized form. The new set of variables, known as principal components (PCs), are uncorrelated and ordered so that the first retains the majority of the variation present in all of the original variables. Assuming a linear relationship between observed variables (species concentrations) and a number of p principal components (sources), PCA is expressed in its simplified form as in equation (3):

$$Z_{ij} = \sum_{k=1}^{p} g_{ik} \times h_{kj} \qquad (3) \qquad 389$$

where Zij is the standardized concentration matrix, k=1,...,p is the number of principal components (sources), while gik is the factor loading matrix of correlations between compounds and PCs (sources) and hkj is the matrix of factor scores.

To better identify groups of correlated compounds, the statistical method called "VARI-MAX rotation" was applied to the standardized concentration matrix to redistribute variances and provide a more interpretable structure to PCs (sources). After rotation, species from the same source are linked to the same PC with high weights (factor loadings), and this PC is then associated with a specific air pollution source before using PMF to quantify this source.

2.5.2. Positive Matrix Factorization (PMF) 401

Another statistical method, called Positive Matrix Factorization (PMF version v5.0) was used to quantify the contribution of the various sources to the particulate emissions. The principles and usage of this approach are detailed in the EPA (Environment Protection Agency) PMF user guide (https://www.epa.gov/air-research/epa-positive-matrix-factorization-50-fundamentals-and-user-guide). This method involves a minimization of an objective cost function Q [29] [54] [55], given by equation (4):

$$Q = \sum_{i=1}^{n} \sum_{j=1}^{m} \frac{e_{ij}^2}{s_{ij}^2} \qquad (4) \qquad 409$$

where eij is the residual associated with i^th^ species concentrations measured in the j^th^ sample, sij is an uncertainty estimate for i^th^ species measured in the j^th^ sample, n and m are respectively the numbers of samples and of species.

PMF uses both measured concentrations and uncertainty estimates to generate calculated chemical profiles and time series associated with each of the chemical profiles. The selection of species to be included in the analysis data depends on the goal of the study and on the quality of available species measurements. PMF analysis, like the PCA method, included the selection of 27 chemical variables (Al, Be, Ti, V, Mn, Co, Fe, Tl, K, As, Rb, Mg,
Na, Ca, Ni, Se, Cr, Cu, Zn, Cd, Sb, EC, OC, Cl⁻, SO₄²⁻, and NO₃⁻) that were detected in the 39 and 36 samples in Bamako and Dakar, respectively. It should be noted that, total suspended particles (TSP) concentrations were included in the factor analysis in order to increase the robustness of the results. Receptor models do not generally provide strict guidelines upon the minimal number of required samples. This is because the number of required samples also depends on the number of analyzed chemical species (source tracers in particular), quality of the data and number of contributing sources [56]. Moreover, it has been seen that the number of samples in the dataset employed in the published literature ranges from ten [57] to tens of thousands [58]. Following these results and due to the number of sources impacting our two sites, 39 and 36 samples were considered as a satisfactory number of observations to conduct PMF apportionment.

To assess the contribution of the identified factors to the collected data, the analysis includes PM mass as an explicit species. In the species concentrations table, missing values were replaced with the median value of species with assigned uncertainty values multiplied by four times the standard deviation. Species with a percentage of missing data over 70% were excluded from the data set, which is the case for Pb and NH₄⁺. The basic criterion retained for PCA and PMF was that aerosol concentrations must be greater than the detection limit (DL). In this study, the DL values were calculated as in Alleman et al. [40]. When the measured concentrations are less than the DL, an adjustment is required; typically, data below the DL are rejected or replaced by DL/2, with an uncertainty of 5 x DL/6 [55]. This late methodology is applied here. Furthermore, a data quality classification based on the signal-to-noise ratio has been implemented, indicating whether the variability in measurements is real or due to data noise. The species were classified as "strong" if the signal to noise ratio was greater than 2, "weak" if it was between 0.2 and 2, and "bad" if it was less than 0.2 [59]. We multiplied the concentration by 3 for "weak" data. Only PM mass and Sb were categorized as "weak" in our measurements. As a result, the EPA PMF user guide recommended setting an extra modeling uncertainty. It was set at 20% in this study.

Before applying the PMF, isolated event such as dust episode was screened, as the contributions due to dust episodes could not be distinguished by the PMF technique [30]. The performance of PMF model was assessed by comparing reconstructed daily mean PM concentrations from all sources to measured values for each species (Tables S8 and S9, Supplement). The largest uncertainties in Bamako (15% for Ni, 16% for EC and 14% for OC) and in Dakar (10% for PM), can be attributed to data quality because missing values in some of these variables (EC and OC) were replaced by their median values, and associated uncertainties were assigned to be four times the mean values. PM is included as a variable, with high uncertainty determined as for missing values [60]. However, with significant correlation coefficients, there was generally good agreement between measured and calculated concentrations, suggesting that the identified factors effectively reproduce measured masses.

# 3. Results and discussion 461

## 3.1. Aerosol chemical mass concentrations 462

Figure 4 shows time series of PM₂.₅, PM₁₀, and TSP mass concentration measured during the POLCA intensive campaigns in Bamako (January 20, 2009 – February 01, 2009) and Dakar (December 01 to December 13, 2009), as well as the 48-hour back trajectories analysis from the Hybrid Single-Particle Lagrangian Intergrated Trajectory (HYSPLIT) modeling system [61]. The concentration of PM₂.₅, PM₁₀, and TSP in Bamako ranged from 131.6 to 1194.3 µg.m⁻³, 151.1 to 1513.1 µg.m⁻³ and 356.7 to 1768.9 µg.m⁻³, with average concentration of 276.8, 503.6 and 705.3 µg.m⁻³, respectively. The extreme daily values were measured on January 22nd, 2009 (Figure 4a): the HYSPLIT back-trajectories suggesting that these high values were caused by a Saharan dust event that occurred that day. According
to Figure 4b, the air mass that was stagnating over Bamako at the beginning of the measurements (January 22-23) was coming from the Southwest part of Algeria and had passed over the Great Western Erg. The dust episode has an impact on all PM modes. The average daily PM concentrations remain very high without the dust event (583 µg.m⁻³ for TSP, 390 µg.m⁻³ for PM₁₀, and 219 µg.m⁻³ for PM₂.₅). In Dakar, the average daily concentrations ranged from 152 to 479 µg.m⁻³ for TSP, 82 to 293 µg.m⁻³ for PM₁₀, and 81 to 258 µg.m⁻³ for PM₂.₅, with average values of 274.9, 155.9 and 138.2 µg.m⁻³, respectively. Both PM₂.₅ and PM₁₀ show peaks on December 3 and 10, corresponding to Thursdays when activity on the Sandaga market, where measuring instruments are located, is higher (Figure 4c). Only TSP shows a peak on December 8 that coincides with the origin of the air mass to the south-east (Figure 4d), implying a contribution from long-range transport in the high concentrations measured. This late hypothesis is in line with the observed dust episode. In addition, most air mass back trajectories were from continental origin, implying a potential mix of sources (dust, biomass burning, and anthropogenic) that may explain the high measured PM concentration during the measurement period (Figure 4d). Lower concentrations were measured in Dakar over the weekend of December 12-13, when cleaner air masses passed over the Atlantic Ocean before arriving at the measurement site. PM₁₀ and PM₂.₅ concentrations are up to twice and 3.4 times the Senegalese daily national standard of 150 and 75 µg.m⁻³, respectively. These limit values are two to three times higher than the other international standards. Mali currently lacks a national air quality standard. In both Bamako (Pearson correlation coefficient r = 0.89-0.94) and Dakar (r = 0.78-0.97), the temporal evolutions of the three aerosol fractions show significant correlations between the different particle sizes, indicating that these aerosols share sources. The concentrations in both Dakar and Bamako (excluding dust episode) are comparable to those measured in Conakry (Guinea) during the Harmattan period (winds blowing northeast from January to February) in 2004, namely 80-358 µg.m⁻³ for PM₁₀ and 38-177 µg.m⁻³ for PM₂.₅ [16]. Furthermore, the average PM₂.₅ concentrations in our study are of the same order of magnitude as the average values of 141.3, 175.3, and 269.7 µg.m⁻³ obtained in a traffic and waste burning site in Abidjan, and traffic site in Cotonou in January 2016 and 2017 [62]. More recent measurements of PM₂.₅ in Dakar show that, for the years 2018-2019, average daily concentrations ranged from 280.6 to 302.7 µg.m⁻³ [63], suggesting an increase in concentration over the previous ten years. This can be explained by rapid urbanization and population growth, as the population of Dakar increase of about 1 million residents between 2007 and 2017. According to these comparisons, the PM concentrations in Dakar and Bamako are typical of the urban pollution in western Africa.

472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Map Daily evolution of PM (TSP, PM₁₀, and PM₂.₅) mass concentrations (a) in Bamako and (c) Dakar, as well as HYSPLIT 48-hour air mass back trajectories arriving (b) in Bamako for the period from 20 January to 1 February 2009, and (d) in Dakar for the period from 1 to 13 December 2009.</figcaption></figure>

Table 1 summarizes the statistics (arithmetic mean, standard deviation and percentage in total mass) of measured chemical composition mass concentration in Bamako and Dakar, while figures S5 and S6 (Supplement) show the temporal evolution of OC, EC, ions and twenty-two metals measured in PM₂.₅, PM₁₀, and TSP.

In general, OC, EC, metals components, and ions account for 15-37%, 3-10%, 7-12%, and 3-6% of PM₂.₅, PM₁₀ and TSP in Bamako, respectively, whereas these percentages are 22-26%, 7-11%, 6-9%, and 13-25% in Dakar (Table 1). Total carbon (OC+EC) was the most abundant component in the atmospheric particle composition, accounting for 18-47% in Bamako and 32-37% in Dakar. These results are consistent with the values reported for some European and Asian urban aerosols by Putaud et al. [64], Querol et al. [65], and Yu et al. [66], who estimated contributions of 20-40% and 25-50% to the ambient PM₁₀ and PM₂.₅ mass, respectively. The average EC(OC) concentrations in Bamako for the entire period were 23.3(102.4) µg.m⁻³ for TSP, 18.8(82.4) µg.m⁻³ for PM₁₀, and 27.1(102.4) µg.m⁻³ for PM₂.₅. In Dakar, the average EC(OC) concentrations were 20.1(69.2) µg.m⁻³ for TSP, 16.4(34.9) µg.m⁻³ for PM₁₀, and 15.9(35.8) µg.m⁻³ for PM₂.₅. Results indicated that, the EC to OC ratio is higher in Dakar (0.29-0.47) than in Bamako (0.23-0.26). Figure S5 shows that during the dust event in Bamako, carbonaceous particles significantly decrease while elements associated with mineral dust significantly increase (especially in TSP and PM₁₀ fractions), highlighting the impact of high dust concentrations on carbonaceous measurements.
<table><caption>**Table 1:** Average 24-hour concentrations (in µg.m⁻³) for TSP, PM₁₀, PM₂.₅ and associated elemental compositions in Bamako (including dust episode) and Dakar. Values in parentheses are percentages in total masses.</caption><thead><tr><th rowspan="2">Species</th><th colspan="3">Bamako</th><th colspan="3">Dakar</th></tr><tr><th>TSP</th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th><th>TSP</th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th></tr></thead><tbody><tr><td><strong>Mass PM</strong></td><td>705.3 ± 99.3</td><td>503.6 ± 112.4</td><td>276.8 ± 94.7</td><td>274.9 ± 27.43</td><td>155.9 ± 15.7</td><td>138.2 ± 12.7</td></tr><tr><td>Cl<sup>-</sup></td><td>3.43 (&lt;1)</td><td>3.36 (&lt;1)</td><td>2.12 (&lt;1)</td><td>8.45 (3)</td><td>8.13 (5)</td><td>1.77 (1)</td></tr><tr><td>NO<sub>3</sub><sup>-</sup></td><td>1.87 (&lt;1)</td><td>1.87 (&lt;1)</td><td>1.18 (&lt;1)</td><td>2.06 (&lt;1)</td><td>2.31 (1)</td><td>0.97 (&lt;1)</td></tr><tr><td>SO<sub>4</sub><sup>2-</sup></td><td>4.44 (&lt;1)</td><td>4.29 (&lt;1)</td><td>2.98 (1)</td><td>6.18 (2)</td><td>7.49 (5)</td><td>4.75 (3)</td></tr><tr><td>Na<sup>+</sup></td><td>2.21 (&lt;1)</td><td>1.70 (&lt;1)</td><td>1.02 (&lt;1)</td><td>5.1 (2)</td><td>5.04 (3)</td><td>1.21 (&lt;1)</td></tr><tr><td>K<sup>+</sup></td><td>3.5 (&lt;1)</td><td>3.17 (&lt;1)</td><td>2.02 (&lt;1)</td><td>0.99 (&lt;1)</td><td>0.99 (&lt;1)</td><td>0.59 (&lt;1)</td></tr><tr><td>Mg<sup>2+</sup></td><td>0.87 (&lt;1)</td><td>0.75 (&lt;1)</td><td>0.49 (&lt;1)</td><td>0.72 (&lt;1)</td><td>0.71 (&lt;1)</td><td>0.28 (&lt;1)</td></tr><tr><td>Ca<sup>2+</sup></td><td>7.14 (1)</td><td>7.09 (1)</td><td>5.29 (2)</td><td>12.17 (4)</td><td>13.65 (9)</td><td>11.43 (8)</td></tr><tr><td><strong>Inorganic ions</strong></td><td>24.06 (3)</td><td>22.99 (5)</td><td>15.69 (6)</td><td>36.09 (13)</td><td>38.95 (25)</td><td>21.50 (16)</td></tr><tr><td>EC</td><td>23.34 (3)</td><td>18.85 (4)</td><td>27.11 (10)</td><td>20.14 (7)</td><td>16.41 (11)</td><td>15.88 (11)</td></tr><tr><td>OC</td><td>102.4 (15)</td><td>82.38 (16)</td><td>102.4 (37)</td><td>69.19 (25)</td><td>34.92 (22)</td><td>35.83 (26)</td></tr><tr><td>TC</td><td>125.74 (18)</td><td>101.23 (20)</td><td>129.51 (47)</td><td>89.33 (32)</td><td>51.33 (33)</td><td>51.71 (37)</td></tr><tr><td>Al</td><td>21.23 (3)</td><td>22.24 (4)</td><td>17.98 (6)</td><td>9.66 (3)</td><td>7.07 (4)</td><td>4.11 (3)</td></tr><tr><td>Fe</td><td>21.91 (3)</td><td>20.41 (4)</td><td>14.28 (5)</td><td>9.00 (3)</td><td>5.54 (3)</td><td>3.49 (2)</td></tr><tr><td>Ti</td><td>2.35</td><td>2.14</td><td>1.47</td><td>0.89</td><td>0.54</td><td>0.34</td></tr><tr><td>Mn</td><td>0.35</td><td>0.32</td><td>0.22</td><td>0.15</td><td>0.11</td><td>0.06</td></tr><tr><td>Zn</td><td>0.18</td><td>0.15</td><td>0.12</td><td>0.21</td><td>0.16</td><td>0.11</td></tr><tr><td>Cr</td><td>0.10</td><td>0.09</td><td>0.06</td><td>0.04</td><td>0.02</td><td>0.02</td></tr><tr><td>V</td><td>0.06</td><td>0.05</td><td>0.04</td><td>0.07</td><td>0.07</td><td>0.05</td></tr><tr><td>Cu</td><td>0.03</td><td>0.04</td><td>0.02</td><td>0.08</td><td>0.19</td><td>0.13</td></tr><tr><td>Ni</td><td>0.02</td><td>0.02</td><td>0.03</td><td>0.03</td><td>0.03</td><td>0.03</td></tr><tr><td>Pb</td><td>0.02</td><td>0.02</td><td>0.04</td><td>0.03</td><td>0.03</td><td>0.02</td></tr><tr><td>Rb</td><td>0.02</td><td>0.02</td><td>0.015</td><td>0.009</td><td>0.007</td><td>0.004</td></tr><tr><td>Co</td><td>0.009</td><td>0.008</td><td>0.005</td><td>0.004</td><td>0.002</td><td>0.0016</td></tr><tr><td>Sb</td><td>0.004</td><td>0.005</td><td>0.004</td><td>0.005</td><td>0.0044</td><td>0.0028</td></tr><tr><td>As</td><td>0.005</td><td>0.0045</td><td>0.003</td><td>0.003</td><td>0.0020</td><td>0.0014</td></tr><tr><td>Be</td><td>0.0008</td><td>0.0008</td><td>0.0005</td><td>0.0003</td><td>0.0002</td><td>0.0001</td></tr><tr><td>Cd</td><td>0.0007</td><td>0.0007</td><td>0.0008</td><td>0.0007</td><td>0.0006</td><td>0.0004</td></tr><tr><td>Se</td><td>0.0003</td><td>0.0003</td><td>0.0002</td><td>0.0007</td><td>0.0008</td><td>0.0005</td></tr><tr><td>Tl</td><td>0.0003</td><td>0.0003</td><td>0.0002</td><td>0.00005</td><td>0.00005</td><td>0.00003</td></tr><tr><td><strong>Metals</strong></td><td>46.29 (7)</td><td>45.53 (9)</td><td>34.29 (12)</td><td>20.18 (7)</td><td>13.78 (9)</td><td>8.38 (6)</td></tr></tbody></table>

As in Bamako, carbonaceous aerosols (EC+OC) are the dominant elements in Dakar, accounting for 32-37% of total mass (Table 1). This is consistent with the findings of Val et al. [9], who discovered a high total carbon content at the same sites, using 13-stage impactors. This is typical of urban areas. The second major element, water-soluble inorganic ions, account for 13% of TSP, 25% of PM₁₀, and 16% of PM₂.₅. Metal mass concentrations range from 8.4 to 20.2 µg.m⁻³, representing 6 to 9% of total mass, regardless of size fraction. Aside from the location of the measurement sites in the Sahara and Sahel regions, Bamako and to a lesser extent Dakar, are characterized by a large number of unpaved roads, which favor particle resuspension from urban sources such as traffic.

The concentration levels for all components were more important in Bamako than in Dakar, except for the water-soluble compounds for which concentrations in Dakar were 2-6 times higher than in Bamako. This is related to differences in location and meteorology,
since Bamako is located in a basin under strong influence of continental air masses, while Dakar is a typical coastal site with high impact of sea breeze which facilitates pollutant dispersion.

## 3.2. Identification and apportionment of sources

### 3.2.1. Enrichment factor calculations

In this section, we calculate Enrichment factors (EF), which are indicators for interpreting the metal element contributions in each source category. Enrichment factors of atmospheric elements were calculated as the ratio of the concentration of an element X to the concentration of a reference element in an aerosol sample compared to the same ratio in average upper crust geological samples [67], as given in the equation (5):

$$ EF = \frac{\left(\frac{x}{E_{ref}}\right)_{\text{aerosol}}}{\left(\frac{x}{E_{ref}}\right)_{\text{crust}}} \qquad (5) $$

where X and E<sub>ref</sub> respectively refer to the concentration (in ppb) of element X and the reference element. Many elements (Al, Ti, Si, Fe, Mn, and Sr), mostly of soil origin are used as references for E<sub>ref</sub> [12] [67]. EF interpretations widely vary according to authors, due to the uncertainties associated with reference elements.

The results obtained for Bamako and Dakar are reported in Figure 5. EF calculations are made using Al, Fe and Ti as crustal references [12] [67]. Alleman et al. [40] proposed a classification of sources based on EF, which we used in our study. Elements with EF values less than 5 (Al, Co, Cr, Fe, K, Mg, Mn, Na, Rb, Ti, Be, and Tl), regardless of the reference element, are considered of natural origin, whereas elements with EF values greater than 50 (Sb, Cd, and Cu) are generally from anthropogenic activities. As shown in the Figure 5, the average EF for Cu, however, exhibits values less than 10 in Bamako and greater than 50 in Dakar, which must be taken into account in the interpretation of sources. Elements with intermediate EF values (between 5 and 50) could originate from both natural and anthropogenic emissions (As, Ca, Ni, V, and Zn). The significant difference in EF for Ni, V, Ca, and Na between Bamako and Dakar indicates that these elements are most likely derived from different sources (Figure 5).
<figure><img src="image_8.png" /><figcaption>**Figure 5:** Enrichment Factor (EF) values for PM elemental composition calculated in a) Bamako and b) Dakar using various reference elements.</figcaption></figure>

### 3.2.2. PCA for source identification

The factor loadings after rotation for Bamako and Dakar are reported in Table 2. The number of optimal Principal Components (PCs) for each measurement site is determined based on the Eigenvalue, which represents the amount of variance explained by each PC greater than one. The table also displays the variability percentage, which expresses the contribution of each principal component, as well as the cumulative percentage, which represents the progressive accumulation of the variability percentage up to 100%. PCA identified five PCs (PC1 to PC5) or sources in both Bamako and Dakar, accounting for 92.9% and 89.5% of the variance, respectively. The factor loadings represent correlations between rotated principal components (PCs) and compounds. Factor loadings greater than 0.600 (highlighted in bold) are considered strong, while those between 0.300 and 0.600 (underlined) are considered moderate [68]. Some compounds have a strong relationship with a specific PC. If these compounds are known to be characteristic of a specific air pollution source, then that source is associated to this PC. Table 3 depicts the possible associations between PCs and pollution sources based on the characteristics of the compounds. However, it is not always possible to link clearly a source to a PC. In this case, a prior knowledge of the sampling sites and their possible sources is crucial.

584

585

586

587

588

589

590

591

592

593

594

595

596

597

598

599

600

601

602

603

604

605

606

607

608
<table><caption>**Table 2:** Correlations between chemical compounds and principal components (PCs) in Bamako and Dakar (factor loadings ≥ 0.600 are highlighted in bold, while moderate factor loadings (0.300-0.600) being underlined, categorization based on Nyanganyura et al. [68]).</caption><thead><tr><th></th><th>PC1</th><th>PC2</th><th>PC3</th><th>PC4</th><th>PC5</th><th>PC1</th><th>PC2</th><th>PC3</th><th>PC4</th><th>PC5</th></tr></thead><tbody><tr><td>Eigenvalue</td><td>15.386</td><td>3.180</td><td>1.284</td><td>1.228</td><td>1.012</td><td>17.126</td><td>3.231</td><td>1.560</td><td>1.137</td><td>1.107</td></tr><tr><td>Variability (%)</td><td>56.986</td><td>11.777</td><td>8.459</td><td>8.251</td><td>7.452</td><td>63.431</td><td>11.967</td><td>5.776</td><td>4.210</td><td>4.100</td></tr><tr><td>Cumulative (%)</td><td>56.986</td><td>68.763</td><td>77.222</td><td>85.472</td><td>92.924</td><td>63.431</td><td>75.398</td><td>81.174</td><td>85.384</td><td>89.484</td></tr><tr><td>Factor loadings</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Al</td><td><strong>0.982</strong></td><td>0.051</td><td>-0.041</td><td>0.035</td><td>-0.120</td><td><strong>0.688</strong></td><td><u>0.334</u></td><td><u>0.300</u></td><td><u>0.341</u></td><td><u>0.385</u></td></tr><tr><td>As</td><td><strong>0.938</strong></td><td>0.143</td><td>0.206</td><td>0.099</td><td>-0.111</td><td>0.154</td><td>0.177</td><td><u>0.470</u></td><td><u>0.577</u></td><td>0.245</td></tr><tr><td>Be</td><td><strong>0.976</strong></td><td>0.105</td><td>-0.086</td><td>0.032</td><td>-0.143</td><td><strong>0.629</strong></td><td><u>0.469</u></td><td>0.237</td><td><u>0.342</u></td><td><u>0.410</u></td></tr><tr><td>Ca</td><td><strong>0.620</strong></td><td><u>0.582</u></td><td>0.167</td><td><u>0.443</u></td><td>0.051</td><td>0.294</td><td><strong>0.746</strong></td><td>0.144</td><td><u>0.399</u></td><td><u>0.341</u></td></tr><tr><td>Cd</td><td>0.184</td><td>-0.218</td><td>0.251</td><td>0.233</td><td><strong>0.643</strong></td><td>0.277</td><td><strong>0.699</strong></td><td>0.262</td><td><u>0.441</u></td><td>0.252</td></tr><tr><td>Co</td><td><strong>0.972</strong></td><td>0.143</td><td>0.015</td><td>0.086</td><td>-0.121</td><td><u>0.542</u></td><td><u>0.347</u></td><td>0.225</td><td><u>0.407</u></td><td><u>0.596</u></td></tr><tr><td>Cr</td><td><strong>0.830</strong></td><td>0.291</td><td><u>0.358</u></td><td>0.232</td><td>0.074</td><td><u>0.509</u></td><td>0.158</td><td>0.259</td><td><u>0.428</u></td><td><strong>0.658</strong></td></tr><tr><td>Cu</td><td><strong>0.642</strong></td><td>0.291</td><td><u>0.512</u></td><td>0.276</td><td>0.089</td><td>0.097</td><td><u>0.360</u></td><td><u>0.312</u></td><td>-0.729</td><td>-0.107</td></tr><tr><td>Fe</td><td><strong>0.964</strong></td><td>0.181</td><td>0.085</td><td>0.129</td><td>-0.068</td><td><strong>0.653</strong></td><td>0.293</td><td><u>0.351</u></td><td><u>0.465</u></td><td><u>0.362</u></td></tr><tr><td>K</td><td><u>0.370</u></td><td><strong>0.879</strong></td><td>0.041</td><td>0.236</td><td>0.051</td><td><strong>0.702</strong></td><td><u>0.415</u></td><td><u>0.357</u></td><td>0.253</td><td>0.282</td></tr><tr><td>Mg</td><td><strong>0.920</strong></td><td>0.241</td><td>-0.069</td><td>0.181</td><td>-0.105</td><td><strong>0.620</strong></td><td>0.297</td><td><u>0.441</u></td><td><u>0.453</u></td><td><u>0.306</u></td></tr><tr><td>Mn</td><td><strong>0.974</strong></td><td>0.130</td><td>0.001</td><td>0.079</td><td>-0.123</td><td><u>0.594</u></td><td><u>0.425</u></td><td><u>0.319</u></td><td><u>0.481</u></td><td><u>0.318</u></td></tr><tr><td>Na</td><td>0.152</td><td><strong>0.930</strong></td><td>0.073</td><td>0.227</td><td>0.071</td><td><u>0.459</u></td><td>0.113</td><td><strong>0.782</strong></td><td>0.239</td><td>0.159</td></tr><tr><td>Ni</td><td><strong>0.870</strong></td><td>0.063</td><td>-0.134</td><td>0.012</td><td>0.119</td><td>-0.178</td><td>-0.004</td><td>-0.049</td><td>-0.062</td><td><strong>0.949</strong></td></tr><tr><td>Rb</td><td><strong>0.928</strong></td><td>0.059</td><td>0.178</td><td>0.128</td><td>-0.003</td><td><strong>0.707</strong></td><td><u>0.394</u></td><td>0.283</td><td><u>0.352</u></td><td><u>0.306</u></td></tr><tr><td>Sb</td><td>-0.202</td><td>-0.081</td><td><strong>0.892</strong></td><td>0.101</td><td><u>0.317</u></td><td><u>0.546</u></td><td><u>0.361</u></td><td><u>0.471</u></td><td>-0.100</td><td>0.267</td></tr><tr><td>Se</td><td><strong>0.841</strong></td><td>0.057</td><td>0.278</td><td>0.299</td><td>0.155</td><td>0.170</td><td><u>0.363</u></td><td><strong>0.657</strong></td><td>-0.319</td><td><u>0.376</u></td></tr><tr><td>Ti</td><td><strong>0.977</strong></td><td>0.107</td><td>-0.056</td><td>0.036</td><td>-0.134</td><td><strong>0.663</strong></td><td><u>0.352</u></td><td>0.282</td><td><u>0.424</u></td><td><u>0.385</u></td></tr><tr><td>Tl</td><td><strong>0.980</strong></td><td>0.038</td><td>0.095</td><td>0.078</td><td>-0.013</td><td><strong>0.612</strong></td><td><u>0.534</u></td><td>0.184</td><td>0.191</td><td><u>0.402</u></td></tr><tr><td>V</td><td><strong>0.980</strong></td><td>0.124</td><td>0.014</td><td>0.076</td><td>-0.104</td><td><u>0.346</u></td><td><strong>0.636</strong></td><td><u>0.342</u></td><td>-0.161</td><td><u>0.381</u></td></tr><tr><td>Zn</td><td>0.264</td><td><u>0.414</u></td><td><strong>0.698</strong></td><td>0.294</td><td><u>0.301</u></td><td><u>0.477</u></td><td><strong>0.659</strong></td><td><u>0.338</u></td><td>0.289</td><td>0.173</td></tr><tr><td>EC</td><td>-0.481</td><td>0.201</td><td>0.143</td><td>0.045</td><td><strong>0.753</strong></td><td>0.078</td><td><u>0.382</u></td><td>0.073</td><td><strong>0.807</strong></td><td>0.074</td></tr><tr><td>OC</td><td>-0.359</td><td><u>0.373</u></td><td>0.288</td><td>-0.013</td><td><strong>0.742</strong></td><td><u>0.472</u></td><td>0.078</td><td>0.117</td><td><strong>0.789</strong></td><td>0.091</td></tr><tr><td>Cl<sup>-</sup></td><td>0.187</td><td><u>0.349</u></td><td><u>0.443</u></td><td><strong>0.728</strong></td><td>0.213</td><td>0.197</td><td>0.097</td><td><strong>0.906</strong></td><td>0.152</td><td>-0.003</td></tr><tr><td>NO<sub>3</sub><sup>-</sup></td><td>0.133</td><td><u>0.312</u></td><td>0.095</td><td><strong>0.908</strong></td><td>0.080</td><td>-0.112</td><td><u>0.382</u></td><td><strong>0.832</strong></td><td>0.021</td><td>0.196</td></tr><tr><td>SO<sub>4</sub><sup>2-</sup></td><td><strong>0.889</strong></td><td>0.230</td><td>0.004</td><td>0.284</td><td>-0.058</td><td>-0.051</td><td><strong>0.778</strong></td><td><u>0.529</u></td><td>-0.042</td><td>0.193</td></tr></tbody></table>

609

610

611

612

613

614

615

616

617

618

619

620

621

622

623

624

625

626

627
<table><caption>**Table 3:** Principal component (PC) and potential source associations (chemical species with high factor loadings are in bold, while the other species have moderate factor loadings). 628 629</caption><thead><tr><th></th><th>Rotated PCs</th><th>Potential Sources</th><th>Characteristic Compounds</th></tr></thead><tbody><tr><td rowspan="5">Bamako</td><td>PC1</td><td>Dust</td><td>Al, As, Be, Cr, Co, Fe, Ti, Mg, Mn, V, Tl, Rb, Se, Ni, SO<sub>4</sub><sup>2-</sup>, Cu, Ca, K</td></tr><tr><td>PC2</td><td>Solid fuel combustion</td><td>Na, K, Ca, Zn, NO<sub>3</sub><sup>-</sup>, Cl<sup>-</sup>, OC</td></tr><tr><td>PC3</td><td>Resuspended road dust</td><td>Sb, Zn, Cu, Cl<sup>-</sup>, Cr</td></tr><tr><td>PC4</td><td>Secondary aerosols</td><td>NO<sub>3</sub><sup>-</sup>, Cl<sup>-</sup>, Ca</td></tr><tr><td>PC5</td><td>Vehicle</td><td>EC, OC, Cd, Zn</td></tr><tr><td rowspan="5">Dakar</td><td>PC1</td><td>Dust</td><td>Al, Be, Rb, K, Ti, Tl, Be, Fe, Mg, Co, Mn, Na, Sb, Cr, Zn, OC, V</td></tr><tr><td>PC2</td><td>Industries, Oil burning</td><td>Ca, Cd, V, Zn, SO<sub>4</sub><sup>2-</sup>, Tl, K, Mn, Be, NO<sub>3</sub><sup>-</sup>, EC, Ti, Al, Sb, Se, Rb, Cu</td></tr><tr><td>PC3</td><td>Salts</td><td>Na, Cl<sup>-</sup>, NO<sub>3</sub><sup>-</sup>, Se, SO<sub>4</sub><sup>2-</sup>, As, K, Mg, Al, Cu, Fe, Mn, Sb, V</td></tr><tr><td>PC4</td><td>Vehicle</td><td>EC, OC, Cd, Co, Cr, Fe, Mg, Mn, As, Al, Be, Ti, Ca, Rb</td></tr><tr><td>PC5</td><td>Resuspend road particles</td><td>Ni, Cr, Co, Be, Al, Tl, Ca, Fe, Mg, Mn, Rb, Se, Ti, V</td></tr></tbody></table>

### 3.2.3. *Description of sources* 630

The PMF was running in various configurations, testing a variety of factors (3 to 7 factors for each set of samples). The optimum solutions for the base run were obtained using 5 factors or sources in both Bamako and Dakar, matching the number of principal components (PCs) obtained using PCA. The correlations between the PCs obtained with PCA and the factors derived from the PMF receptor model were calculated to verify the coherence of these source identifications (Table 4). The two methods show a strong correspondence, which allows an investigation of the source profiles in Bamako and Dakar. The goal is to identify the nature of the sources by calculating their relative contribution to each of the chemical compounds.

<table><caption>**Table 4:** Pearson correlations between PCA principal components (PCs) and PMF model sources (Factors). P values ≤0.05 were used to determine significance. 642 643</caption><thead><tr><th>Bamako</th><th>Factor 1</th><th>Factor 2</th><th>Factor 3</th><th>Factor 4</th><th>Factor 5</th></tr></thead><tbody><tr><td>PC1</td><td>-0.72</td><td>0.00</td><td>-0.07</td><td><strong>0.84</strong></td><td>-0.06</td></tr><tr><td>PC2</td><td>-0.13</td><td>-0.46</td><td><strong>0.82</strong></td><td>-0.40</td><td>0.10</td></tr><tr><td>PC3</td><td>0.30</td><td><strong>0.71</strong></td><td>-0.13</td><td>-0.63</td><td>-0.07</td></tr><tr><td>PC4</td><td>-0.22</td><td>0.19</td><td>-0.01</td><td>-0.56</td><td><strong>0.78</strong></td></tr><tr><td>PC5</td><td><strong>0.81</strong></td><td>0.13</td><td>-0.17</td><td>-0.70</td><td>-0.05</td></tr></tbody></table>

<table><thead><tr><th>Dakar</th><th>Factor 1</th><th>Factor 2</th><th>Factor 3</th><th>Factor 4</th><th>Factor 5</th></tr></thead><tbody><tr><td>PC1</td><td>-0.24</td><td>0.14</td><td>-0.33</td><td><strong>0.84</strong></td><td>-0.34</td></tr><tr><td>PC2</td><td>-0.32</td><td>-0.13</td><td>-0.18</td><td>0.28</td><td><strong>0.57</strong></td></tr><tr><td>PC3</td><td><strong>0.89</strong></td><td>-0.38</td><td>-0.34</td><td>-0.33</td><td>0.05</td></tr><tr><td>PC4</td><td>-0.16</td><td><strong>0.78</strong></td><td>-0.15</td><td>0.28</td><td>-0.72</td></tr><tr><td>PC5</td><td>-0.35</td><td>-0.11</td><td><strong>0.57</strong></td><td>0.32</td><td>-0.38</td></tr></tbody></table>
<table>
  <tr>
    <td><strong>Bamako</strong></td>
    <td>644</td>
  </tr>
  <tr>
    <td>From PMF output, the comparison of reconstructed daily mean PM concentrations from all sources with measured PM concentrations demonstrates that the identified factors effectively reproduce measured masses. Table S8 (Supplement) compares calculated values with measured concentrations for each species to assess the performance of the PMF model. The largest uncertainties, such as 15% for Ni, 16% for EC and 14% for OC and the low correlation coefficients, can be attributed to data quality because EC and OC missing values were replaced by their median values and associated uncertainties were assigned to be four times the mean values. However, in general, significant correlation coefficients were obtained between measured and calculated concentrations.</td>
    <td>645<br>646<br>647<br>648<br>649<br>650<br>651<br>652<br>653<br>654</td>
  </tr>
  <tr>
    <td><em>Motor vehicles</em>: High factor loadings (0.742-0.753) were obtained for this factor in PCA (Table 2) and high contributions (83-94%) in PMF (Figure 6), both for OC and EC, which can be linked with fuel combustion emissions [69]. Several authors [70] [71] [72] have reported OC/EC mass concentration ratios of 1.1 for traffic emissions, 3.3 for secondary organic carbon formation, 6.6 for biofuel combustion and 12 for long-range transport. The OC/EC ratio for this factor ranges from 3.8 to 4.4 in Table 2, indicating probable mixtures of various combustion sources such as two-wheel two-stroke vehicles (using a mixture of oil and gasoline), personal and public car emissions. The high Cd (0.643) and moderate Sb (0.317) and Zn (0.301) loadings, as well as their presence in PMF (26%, 39% and 18%, respectively), support tire wear and brakes linings residuals [73] [74] as a contributing source.</td>
    <td>655<br>656<br>657<br>658<br>659<br>660<br>661<br>662<br>663<br>664<br>665</td>
  </tr>
  <tr>
    <td><em>Resuspended road dust</em>: In PMF, Sb (61%), Zn (51%) and Cd (48%) dominate the factor loadings. In the PCA results, Sb and Zn have also the highest factor loadings (0.892 and 0.698, respectively) (Table 2). Sb, Zn and Cd trends in Bamako did not show peaks during the dust episode (Figure S5), indicating that these elements are primarily emitted by anthropogenic activities such as brakes and tire wear, fuel combustion, or biomass burning [75] [76]. All of these point sources were present near the site, with open waste burning occurring in Bamako in the evenings. Furthermore, the high EF values reported in Figure 5 support the anthropogenic origin of these elements [40]. Sb, Zn, and Cd emissions from tire and brake wear have been well documented [73] [74] and the majority of roads surrounding the measurement site were unpaved. According to some studies [75] [76], particles from both natural and anthropogenic sources may aggregate in the soil before being resuspended as a result of vehicle traffic. This explains why this source could be classified as traffic.</td>
    <td>666<br>667<br>668<br>669<br>670<br>671<br>672<br>673<br>674<br>675<br>676<br>677<br>678<br>679</td>
  </tr>
  <tr>
    <td><em>Solid fuel combustion</em>: The factor was notable for its high Na and K contributions to total PM mass, which were around 99% for Na and 75% for K in the PMF analysis (Figure 6). Furthermore, in the PCA, these two species have the highest factor loadings (0.879-0.930) (Table 2, PC2). K is a tracer typical for emissions from wood combustion [77]. Zn, Fe, Cr and Co are also abundant in this factor, accounting for 20% of the total mass of each species, reflecting the importance of emissions from common local activities at the sampling site, such as residential combustions for cooking, street food preparation, traditional forges, smelters, open waste burning and unidentified combustion sources. Cr emissions from coal combustion for power generation and waste incineration are well established [16]. Because of the mix of anthropogenic combustion sources, we can conclude that this factor is most likely related to solid fuel combustions.</td>
    <td>680<br>681<br>682<br>683<br>684<br>685<br>686<br>687<br>688<br>689<br>690<br>691</td>
  </tr>
  <tr>
    <td><em>Crustal dust</em>: The factor is characterized by high contribution of Ti (67%), Be (70%), Mn (58%), Co (56%), V (56%), Al (63%), and Fe (46%) (Figure 6), with the majority of them exhibiting high loadings in the PCA results (Table 2, PC4). Metals were the second most abundant component in Bamako. Al was the most abundant metal found in all particle</td>
    <td>692<br>693<br>694<br>695</td>
  </tr>
</table>
modes (Table 1), with average daily concentrations of 17.9 µg.m⁻³ in PM₂.₅, 22.2 µg.m⁻³ in PM₁₀, and 21.2 µg.m⁻³ in TSP, followed by Fe, which had values of 14.3 µg.m⁻³ in PM₂.₅, 20.4 µg.m⁻³ in PM₁₀, and 21.9 µg.m⁻³ in TSP. During the dust episode, daily Al and Fe concentrations rise by nearly a factor of 13 (from 6.3 to 82.2 µg.m⁻³) and 12 (from 5.3 to 66.7 µg.m⁻³), respectively. The peak concentrations of Al and Fe during the dust episode (Figure S5) are very close to the values reported by Adepetu et al. [78] in typical continental crustal rocks of 81.3 and 50 µg.m⁻³, respectively. The average Fe/Al ratio varies from 0.79 (PM₂.₅) to 1.03 (TSP). These values are comparable to those reported in the literature for typical crustal dust elements [77] [79], which is consistent with their enrichment factors (EF) close to unity, regardless of the reference element (Figure 5). This source is most likely due to Saharan desert crustal dust predominating in PM concentrations during dust storms.

**Secondary aerosols:** The factor contains significant contributions from NO₃⁻ (65%), Cl⁻ and Mg (53%), SO₄²⁻ (45%), and Ca (43%) (Figure 6). The PCA results also show that some of these compounds have high factor loadings (0.908 for NO₃⁻ and 0.728 for Cl⁻) and Ca has moderate loading (0.443) (Table 2). Except for Mg and Ca, none of these water-soluble species showed a peak during the dust episode (Figure S5), indicating sources other than dust. SO₄²⁻ in this factor could result from heterogeneous processes in the presence of dust. Because the concentrations of non-soluble Mg and Ca peak during the dust episode, a fraction of these two compounds can be considered from crustal dust. This is consistent with the EF values of less than 1 for these two elements (Figure 5). Ca is also a common tracer typical for emissions from building construction and cement production [80]. NO₃⁻ and Cl⁻ are also linked to anthropogenic activities [77] [81]. Chemical processes in dust wells present during the field campaign may result in secondary aerosol formation because Bamako is located in a basin with low wind speeds and high temperatures that favor gas-to-particle conversion mechanisms [82]. This source results from the interaction of compounds primarily from anthropogenic (NO₃⁻, Cl⁻) and natural (Mg, Ca) sources.
<figure><img src="image_9.png" /><figcaption>**Figure 6:** Profiles of the five factors identified by the PMF in Bamako.</figcaption></figure>

## Dakar

The identified factors effectively reproduce measured masses when daily mean reconstructed PM concentrations from all sources are compared to measured values. Table S9 (Supplement) compares the measured and calculated concentrations for each species in Dakar. Uncertainties are generally low, with the highest being 10% for PM as a result of including PM as a variable and assigning an uncertainty value 4 times the PM concentration values [83] [84]. Significant correlation coefficients are also obtained between modeled and measured species concentrations.

**Sea salts:** This source is characterized by high contribution of Cl⁻ (97%), Na (80%), and moderate presence of Mg (29%) (Figure 7). PCA analysis also reveals that these elements have high factor loadings (0.906, 0.782, and 0.441, respectively) (Table 2, PC3). These three species account for 40, 24 and 3% of the total mass of this aerosol source, respectively, in close agreement with their corresponding percentages in seawater, which are 55, 30 and 3% [85]. The ratios Na/Mg and K/Mg were investigated to confirm the marine origin of these elements. The calculated ratio of 8.3 is very close to the 8.4 found in sea water [86]. Yuan et al. [87] and Alleman et al. [40] also found values of 8.4 and 8.5, respectively. The calculated K/Mg ratio for this factor (0.34) is also of the same order as in sea water (0.31)

725

726

727

728

729

730

731

732

733

734

735

736

737

738

739

740

741

742

743

744

745
[86]. This source also contains $NO_3^-$ (53%) and $SO_4^{2-}$ (31%) species with loading in the PCA of 0.832 and 0.529, respectively. As a result of condensation processes, sulfate and nitrate may be found in marine aerosols [80]. This is consistent with Dakar being surrounded on the three sides by the Atlantic Ocean (Figure 1).

746

747

748

749

750

*Motor Vehicles*: This source contains 86% of total OC, 62% of total EC, 68% of total As (Figure 7), and PCA results show the highest factor loadings of 0.807, 0.789 and 0.577 these species, respectively (Table 2, PC4). These species are common combustion tracers, such as motor vehicle exhaust fumes [80] [88]. The OC/EC ratio of 3.3 is comparable to the value found in Pio et al. [72] for traffic diesel and gasoline emissions. This value is of the same order as previously reported, namely 2.1±1.0 obtained from 13-stage impactors in Dakar [9]. PCA results also show that this factor contributes to the observed concentrations of Cd, Ca, Fe, Ti, Co, Cr, Mn, Sb, Rb, and Mg, with moderate factor loadings (0.340-0.480). The corresponding contributions from PMF analysis are ranged from 20 to 44%. Some of these elements are typically produced by natural sources (Ti, Fe, Cr, Co, and Mn), anthropogenic activities (Cd, Sb) or both sources (As, Ca). This implies the possibility of natural compounds aggregating with fresh aerosols from motor vehicles.

751

752

753

754

755

756

757

758

759

760

761

762

763

*Industries, Oil burning*: This source is characterized by high Ni (100%) and Zn (66%) contribution, as well as moderate Cr and Co contributions (Figure 7). PCA analysis also reveals high factor loadings for Ni (0.949), Cr (0.658) and moderate loadings for Co (0.596) and Se (0.376) (Table 2, PC5), indicating emissions from oil combustion [80] [89] related to hydrocarbon storage and transit at the port of Dakar located, which is close to the site. Zn is also produced by lubricating oil combustion and tire wear [75]. Additional species (Cr, Co, and Se) in this factor are most likely related to industrial activities such as manufacturing, local smelters, and activities like car repair. The measurement site was close to the two biggest car parks in Dakar. In addition, the Dakar industrial area is about one kilometer northeast of the measurement site, and winds from this direction prevailed throughout the campaign on average (Figure 1).

764

765

766

767

768

769

770

771

772

773

774

775

*Mineral dust*: The factor reveals high loadings of Al, Be, Ti, and Fe in PCA, generally greater than 0.6 (Table 2, PC1), and contains high contributions of these elements ranging from 40% to 60% in PMF (Figure 7). Additional elements such as Rb, Co, Cr and Mn are also abundant in this factor. EF calculations (Figure 5) show that all of these elements have values close to unity, confirming their natural origin [16] [12] [80] [90]. This result is consistent with mineral dust.

776

777

778

779

780

781

782

*Resuspended road particles*: The factor has high contribution of Cu (100%), V (52%), Se (51%), $SO_4^{2-}$ (46%), and moderate contribution of Cd, K, Sb, Tl, EC, and $NO_3^-$ (25-30%) (Figure 7). PCA also exhibits significant loadings of V, $SO_4^{2-}$, Cd, Zn (>0.600), and moderate loadings of Se, K, Sb, Tl, EC, and $NO_3^-$ (0.300-0.600) (Table 2, PC2). Cu, Sb, and Cd are thought to be anthropogenic due to their high EF values (>50) (Figure 5). Cu and Sb are components of brake linings [73] [74], while Se, K are good tracers for coal/wood combustion [91] [92]. $SO_4^{2-}$ in this factor could result from heterogeneous processes in the presence of dust. The presence of all these compounds indicates that this factor comes from a variety of sources. These compounds were resuspended by road traffic after being aggregated through complex processes.

783

784

785

786

787

788

789

790

791

792

793

794
<figure><img src="image_10.png" /><figcaption>**Figure 7:** Profiles of the five factors identified by the PMF in Dakar.</figcaption></figure>

### 3.2.4. Source apportionment via PMF

Only PM₁₀ and PM₂.₅ source contributions based on PCA/PMF in Bamako and Dakar are displayed due to their importance in health studies (Figure 8). Traffic emissions include motor vehicles (vehicular and two-wheel two-stroke) and resuspended road dust. At Bamako, this source comprised 47% of PM₂.₅ and 45% of PM₁₀ masses, due to the fact that a large proportion of Bamako's population driving two-wheel two-stroke vehicles, the frequent traffic congestion throughout the city, and the street qualities, which are mostly cobble or unpaved. Only a few major streets use asphalt. Crustal dust was the second dominant source of PM₁₀ mass (30%), but remains an important source of PM₂.₅ fraction (24%). This can be explained by the measurement period in the dry season, when the Har-mattan winds brought dust from the Sahara Desert. Two additional sources were identi-fied: one defined as secondary aerosols, which contributed 16% and 10% of the PM₁₀ and PM₂.₅ masses, respectively; and another identified as solid fuel combustion which account for 13% of the total PM₁₀ and 16% of the total PM₂.₅ (Figure 8a).

Traffic emissions, as expected, are the most significant source of PM₁₀ and PM₂.₅ at the Dakar site, accounting for approximately 49% of total PM₁₀ and PM₂.₅ concentrations (Fig-ure 8b). It should be noted that majority of the Senegal's vehicle fleet (more than half) operates in Dakar and is predominately powered by outdated diesel engines. Old diesel vehicles account for more than 60% of Dakar's public transportation. The proportion of marine aerosols or sea salts in PM₁₀ (20%) is slightly higher than in PM₂.₅ (15%). Mineral dust, which contributes approximately 25% and 16% of the aerosol mass on average for

795

796

797

798

799

800

801

802

803

804

805

806

807

808

809

810

811

812

813

814

815

816

817

818
PM₁₀ and PM₂.₅ mass concentrations, was the second most important source, as it was in Bamako. The industrial sector and oil-burning activities together contributed about 10-11% of the measured PM₁₀ and PM₂.₅ concentrations (Figure 8b). The identification of this source can be explained by location of the port industrial zone close to the sampling site (Figure 1), where oil refinery activities and the storage and transit of hydrocarbons were carried out. Using a year of measurements (2018-2019), Kebe et al. [63] identified four sources in Dakar: traffic emissions, mineral dust, industrial emissions and sea salt/secondary aerosols. However, they did not quantify their contribution to PM₂.₅ mass for a comparison with our study.

According to our findings in both Bamako and Dakar cities, anthropogenic activities are the primary sources air pollution in West Africa, accounting for 60-65% of total concentration regardless of PM mode.

<figure><img src="image_11.png" /><figcaption>**Figure 8:** PM₁₀ and PM₂.₅ source contributions estimated by PMF in a) Bamako and b) Dakar.</figcaption></figure>

## 4. Conclusions

To identify the levels and the major sources of air pollution, this paper characterizes the concentrations of PM (TSP, PM₁₀ and PM₂.₅) and their chemical composition at two west African traffic sites (Bamako in Mali and Dakar in Senegal). Carbonaceous aerosol (OC, EC), major ions (Na⁺, K⁺, Mg²⁺, Ca²⁺, SO₄²⁻, NO₃⁻, and Cl⁻) and 21 metal elements were analyzed. Daily PM₂.₅ and PM₁₀ concentrations in the dry season exceeded the WHO's 24-hour air quality guidelines established in 2005 of 25 and 50 µg.m⁻³, respectively.

Using PCA and PMF methodologies, five sources were extracted and quantified in both Bamako (crustal dust, resuspended road dust, solid fuel combustion, motor vehicles and secondary aerosols) and Dakar (resuspended road particles, mineral dust, sea salts, motor vehicle emissions, and industries and oil combustion). Traffic sector, which includes motor vehicle emissions and resuspended road dust, was the dominant source of compounds in PM₂.₅ and PM₁₀ at both sites (45-49%). As a result, anthropogenic sources predominated, contributing up to 60-65% of the PM mass. These measurements were taken in 2009, and this detailed analysis could serve as a basis for future research on how the composition of particulate matter in African has changed since the late 2000s.

833

834

835

836

837

838

839

840

841

842

843

844

845

846

847

848

849

850

851

852

853

854
Supplementary Materials 855

The following supporting information can be downloaded at: www.mdpi.com/xxx/s1. 856

Author Contributions 857

Conceptualization, T.D. and C.L.; methodology, T.D.; C.L. and C.G-L.; validation, T.D.; E.G. and C.L.; formal analysis, T.D.; M-R.L-O.; E.G. and C. Z.; investigation, T.D.; C.L.; C.G-L.; S.A.N.; V.Y. and M-R.L-O.; writing—original draft preparation, T.D.; C.L. and C.G-L.; writing—review and editing, T.D.; C.L. and C.G.; visualization, T.D.; supervision, C.L. All authors have read and agreed to the published version of the manuscript. 858

Acknowledgments 864

This research was funded by the CORUS 2 program (Cooperation for Academic and Scientific Research), grant number 94-2009 and has been supported in part by the AQ-WATCH (Air Quality: Worldwide Analysis and Forecasting of Atmospheric Composition for Health) project funded by the European Commission under grant 870301. 865

Institutional Review Board Statement 871

Not applicable. 872

Informed Consent Statement 873

Not applicable. 874

Data Availability Statement 875

The data presented in this study are available on request from the corresponding author. 878

Acknowledgments 879

The authors would like to thank the CORUS2 POLCA and AMMA2 programs. We also want to thank the Paul Sabatier University for their financial support, the University of Bamako and the Laboratory of Atmospheric Physics and Oceanographic - Simeon Fongang at the University Cheikh Anta Diop of Dakar for their collaboration. Cyril Zouiten and the GET laboratory (*Géosciences Environnement Toulouse*) are acknowledged for element trace analysis. 880

Conflicts of Interest 881

The authors declare no conflict of interest. 887

References 888

[1] IPCC 2007, "AR4 Climate Change 2007: Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change," [Core Writing Team, Pachauri, R.K and Reisinger, A. (eds.)], Geneva, Switzerland, IPCC, 2007. Accessed: Mar. 28, 2023. [Online]. Available: https://www.ipcc.ch/report/ar4/syr/

[2] D. W. Dockery and C. A. Pope, "Acute Respiratory Effects of Particulate Air Pollution," *Annu. Rev. Public Health*, vol. 15, no. 1, pp. 107-132, 1994, doi: 10.1146/annurev.pu.15.050194.000543.
[3] C. A. Pope 3rd et al., "Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution," JAMA J. Am. Med. Assoc., vol. 287, no. 9, pp. 1132–1141, Mar. 2002, Accessed: Sep. 10, 2012. [Online]. Available: http://www.ncbi.nlm.nih.gov/pubmed/11879110

[4] R. B. Hamanaka and G. M. Mutlu, "Particulate Matter Air Pollution: Effects on the Cardiovascular System," Front. Endocrinol., vol. 9, p. 680, Nov. 2018, doi: 10.3389/fendo.2018.00680.

[5] S.-L. Huang, M.-K. Hsu, and C.-C. Chan, "Effects of submicrometer particle compositions on cytokine production and lipid peroxidation of human bronchial epithelial cells," Environ. Health Perspect., vol. 111, no. 4, pp. 478–482, Apr. 2003, Accessed: Mar. 28, 2023. [Online]. Available: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1241431/

[6] J. Seagrave et al., "Lung toxicity of ambient particulate matter from southeastern U.S. sites with different contributing sources: relationships between composition and effects," Environ. Health Perspect., vol. 114, no. 9, pp. 1387–1393, Sep. 2006, doi: 10.1289/ehp.9234.

[7] J. Lelieveld, J. S. Evans, M. Fnais, D. Giannadaki, and A. Pozzer, "The contribution of outdoor air pollution sources to premature mortality on a global scale," Nature, vol. 525, no. 7569, pp. 367–371, Sep. 2015, doi: 10.1038/nature15371.

[8] M. S. Happo et al., "Chemical compositions responsible for inflammation and tissue damage in the mouse lung by coarse and fine particulate samples from contrasting air pollution in Europe," Inhal. Toxicol., vol. 20, no. 14, pp. 1215–1231, Nov. 2008, doi: 10.1080/08958370802147282.

[9] S. Val et al., "Physico-chemical characterization of African urban aerosols (Bamako in Mali and Dakar in Senegal) and their toxic effects in human bronchial epithelial cells: description of a worrying situation," Part. Fibre Toxicol., vol. 10, p. 10, Apr. 2013, doi: 10.1186/1743-8977-10-10.

[10] J. Bertrand, J. Baudet, and A. Drochon, "Importance des aérosols naturels en Afrique de l'Ouest," Journal de Recherche Atmosphérique, pp. 845–860, 1974. [Online]. Available: http://bibliotheque.meteo.fr/exl.php/cadcgp.php?CMD=CHERCHE&MO-DELE=vues/mf_-_internet_recherche_avancee_anonyme/tpl-r.html&WHERE_IS_DOC_REF_LIT=QUE00025379&&TABLE=ILS_DOC

[11] O. I. Asubiojo et al., "Elemental characterization of airborne particulates at two Nigerian locations during the Harmattan season," J. Radioanal. Nucl. Chem., vol. 167, no. 2, pp. 283–293, Jan. 1993, doi: 10.1007/BF02037187.

[12] M. A. H. Eltayeb, J. Injuk, W. Maenhaut, and R. E. Van Grieken, "Elemental Composition of Mineral Aerosol Generated from Sudan Sahara Sand," J. Atmospheric Chem., vol. 40, no. 3, pp. 247–273, Nov. 2001, doi: 10.1023/A:1012272208129.

[13] K. Deboudt, P. Flament, M. Choël, A. Gloter, S. Sobanska, and C. Colliex, "Mixing state of aerosols and direct observation of carbonaceous and marine coatings on African dust by individual particle analysis: MIXING STATE AND COATING IN AFRICAN DUST," J. Geophys. Res. Atmospheres, vol. 115, no. D24, Dec. 2010, doi: 10.1029/2010JD013921.

[14] R. E. Arku et al., "Characterizing air pollution in two low-income neighborhoods in Accra, Ghana," Sci. Total Environ., vol. 402, no. 2–3, pp. 217–231, Sep. 2008, doi: 10.1016/j.scitotenv.2008.04.042.

[15] K. S. Kouassi et al., "Oxidative damage induced in A549 cells by physically and chemically characterized air particulate matter (PM2.5) collected in Abidjan, Côte d'Ivoire," J. Appl. Toxicol. JAT, vol. 30, no. 4, pp. 310–320, May 2010, doi: 10.1002/jat.1496.

[16] J. P. Weinstein, S. R. Hedges, and S. Kimbrough, "Characterization and aerosol mass balance of PM2.5 and PM10 collected in Conakry, Guinea during the 2004 Harmattan period," Chemosphere, vol. 78, no. 8, pp. 980–988, Feb. 2010, doi: 10.1016/j.chemosphere.2009.12.022.

[17] K. L. Dionisio et al., "Within-neighborhood patterns and sources of particle pollution: mobile monitoring and geographic information system analysis in four communities in Accra, Ghana," Environ. Health Perspect., vol. 118, no. 5, pp. 607–613, May 2010, doi: 10.1289/ehp.0901365.

[18] E.-M. Assamoi and C. Liousse, "A new inventory for two-wheel vehicle emissions in West Africa for 2002," Atmos. Environ., vol. 44, no. 32, pp. 3985–3996, Oct. 2010, doi: 10.1016/j.atmosenv.2010.06.048.
[19] E. H. T. Doumbia et al., "Real time black carbon measurements in West and Central Africa urban sites," Atmos. Environ., vol. 54, pp. 529-537, Jul. 2012, doi: 10.1016/j.atmosenv.2012.02.005.

941

[20] A. Mukherjee and M. Agrawal, "World air particulate matter: sources, distribution and health effects," Environ. Chem. Lett., vol. 15, no. 2, pp. 283-309, Jun. 2017, doi: 10.1007/s10311-017-0611-9.

943

[21] S. Gnamien et al., "Particulate Pollution in Korhogo and Abidjan (Cote d'Ivoire) during the Dry Season," Aerosol Air Qual. Res., vol. 21, no. 1, p. 200201, 2021, doi: 10.4209/aaqr.2020.05.0201.

945

[22] R. K. Larsen and J. E. Baker, "Source Apportionment of Polycyclic Aromatic Hydrocarbons in the Urban Atmosphere: A Comparison of Three Methods," Environ. Sci. Technol., vol. 37, no. 9, pp. 1873-1881, May 2003, doi: 10.1021/es0206184.

946

[23] H. Yang et al., "The chemical composition of inorganic and carbonaceous materials in PM2.5 in Nanjing, China," Atmos. Environ., vol. 39, no. 20, pp. 3735-3749, Jun. 2005, doi: 10.1016/j.atmosenv.2005.03.010.

949

[24] L. Xu et al., "Characterization of PM10 atmospheric aerosol at urban and urban background sites in Fuzhou city, China," Environ. Sci. Pollut. Res., vol. 19, no. 5, pp. 1443-1453, 2012, doi: 10.1007/s11356-012-0759-6.

951

[25] G. E. Gordon, "Receptor models," Environ. Sci. Technol., vol. 22, no. 10, pp. 1132-1142, Oct. 1988, doi: 10.1021/es00175a002.

952

[26] P. K. Hopke, "An introduction to receptor modeling," Chemom. Intell. Lab. Syst., vol. 10, no. 1-2, pp. 21-43, Feb. 1991, doi: 10.1016/0169-7439(91)80032-L.

954

[27] R. C. Henry, "HISTORY AND FUNDAMENTALS OF MULTIVARIATE AIR QUALITY RECEPTOR MODELS," Chemom. Intell. Lab. Syst., vol. 37, no. 1, pp. 37-42, 1997, Accessed: Sep. 25, 2012. [Online]. Available: http://cat.inist.fr/?aModele=afficheN&cpsidt=10777286

957

[28] J. G. Watson, T. Zhu, J. C. Chow, J. Engelbrecht, E. M. Fujita, and W. E. Wilson, "Receptor modeling application framework for particle source apportionment," Chemosphere, vol. 49, no. 9, pp. 1093-1136, Dec. 2002, doi: 10.1016/S0045-6535(02)00243-6.

959

[29] P. Paatero and U. Tapper, "Analysis of different modes of factor analysis as least squares fit problems," Chemom. Intell. Lab. Syst., vol. 18, no. 2, pp. 183-194, Feb. 1993, doi: 10.1016/0169-7439(93)80055-M.

961

[30] Y. Song, S. Xie, Y. Zhang, L. Zeng, L. G. Salmon, and M. Zheng, "Source apportionment of PM2.5 in Beijing using principal component analysis/absolute principal component scores and UNMIX," Sci. Total Environ., vol. 372, no. 1, pp. 278-286, Dec. 2006, doi: 10.1016/j.scitotenv.2006.08.041.

964

[31] M. S. Callén, M. T. de la Cruz, J. M. López, M. V. Navarro, and A. M. Mastral, "Comparison of receptor models for source apportionment of the PM10 in Zaragoza (Spain)," Chemosphere, vol. 76, no. 8, pp. 1120-1129, Aug. 2009, doi: 10.1016/j.chemosphere.2009.04.015.

966

967

[32] R. Zhang, Z. Han, Z. Shen, and J. Cao, "Continuous measurement of number concentrations and elemental composition of aerosol particles for a dust storm event in Beijing," Adv. Atmospheric Sci., vol. 25, no. 1, pp. 89-95, 2008, doi: 10.1007/s00376-008-0089-2.

969

970

[33] M. Adon et al., "Measurements of NO2, SO2, NH3, HNO3 and O3 in West African urban environments," Atmos. Environ., vol. 135, pp. 31-40, Jun. 2016, doi: 10.1016/j.atmosenv.2016.03.050.

972

[34] DNSI, "Mali - Recensement général de la population et de l'habitat (2009)," Direction nationale de la statistique et de l'information, (DNSI) - Ministère de l'économie et du plan, Mali, MLI-INSTAT-RGPH-2009, 2019. Accessed: Mar. 28, 2023. [Online]. Available: https://demostaf.web.ined.fr/index.php/catalog/348

973

974

975

[35] ANSD, "Rapport national sur la Situation économique et sociale du Sénégal, Démographie," Agence National de la Statistique et de la Démographie, Edition 2009, 2009. [Online]. Available: https://www.ansd.sn/ressources/ses/chapitres/1-Démographie_2009.pdf

976

977

978

[36] M.-R. Ouafo-Leumbe et al., "Chemical composition and sources of atmospheric aerosols at Djougou (Benin)," Meteorol. Atmospheric Phys., vol. 130, no. 5, pp. 591-609, Oct. 2018, doi: 10.1007/s00703-017-0538-5.

979

980
[53] J. C. Chow, J. G. Watson, L.-W. A. Chen, M.-C. O. Chang, and G. P. Miranda, "Comparison of the DRI/OGC and Model 2001 Thermal/Optical Carbon Analyzers," 2005.

1023

[54] P. Paatero and U. Tapper, "Positive matrix factorization: A non-negative factor model with optimal utilization of error estimates of data values," *Environmetrics*, vol. 5, no. 2, pp. 111–126, 1994, doi: 10.1002/env.3170050203.

1025

[55] A. V. Polissar, P. K. Hopke, and R. L. Poirot, "Atmospheric Aerosol over Vermont: Chemical Composition and Sources," *Environ. Sci. Technol.*, vol. 35, no. 23, pp. 4604–4621, Dec. 2001, doi: 10.1021/es0105865.

1027

[56] Y. Chan et al., "Using multiple type composition data and wind data in PMF analysis to apportion and locate sources of air pollutants," *Atmos. Environ.*, vol. 45, no. 2, pp. 439–449, Jan. 2011, doi: 10.1016/j.atmosenv.2010.09.060.

1029

[57] H. Laupsa, B. Denby, S. Larssen, and J. Schaug, "Source apportionment of particulate matter (PM2.5) in an urban area using dispersion, receptor and inverse modelling," *Atmos. Environ.*, vol. 43, no. 31, pp. 4733–4744, Oct. 2009, doi: 10.1016/j.atmosenv.2008.07.010.

1032

[58] W. Yue et al., "Source apportionment of ambient fine particle size distribution using positive matrix factorization in Erfurt, Germany," *Sci. Total Environ.*, vol. 398, no. 1–3, pp. 133–144, Jul. 2008, doi: 10.1016/j.scitotenv.2008.02.049.

1034

[59] P. Paatero and P. K. Hopke, "Discarding or downweighting high-noise variables in factor analytic models," *Anal. Chim. Acta*, vol. 490, no. 1–2, pp. 277–289, Aug. 2003, doi: 10.1016/S0003-2670(02)01643-4.

1036

[60] I. Gupta, A. Salunkhe, and R. Kumar, "Source Apportionment of PM10 by Positive Matrix Factorization in Urban Area of Mumbai, India," *Sci. World J.*, vol. 2012, May 2012, doi: 10.1100/2012/585791.

1038

[61] R. R. Draxler and G. D. Hess, "An Overview of the HYSPLIT_4 Modeling System of Trajectories, Dispersion, and Deposition," *Aust. Meteorol. Mag.*, vol. 47, pp. 295–308, 1998.

1040

[62] A. J. Adon et al., "Physico-chemical characterization of urban aerosols from specific combustion sources in West Africa at Abidjan in Côte d'Ivoire and Cotonou in Benin in the frame of DACCIWA program," *Atmospheric Chem. Phys. Discuss.*, pp. 1–69, Jun. 2019, doi: 10.5194/acp-2019-406.

1043

[63] M. Kebe et al., "Source Apportionment and Assessment of Air Quality Index of PM2.5–10 and PM2.5 in at Two Different Sites in Urban Background Area in Senegal," *Atmosphere*, vol. 12, no. 2, p. 182, Jan. 2021, doi: 10.3390/atmos12020182.

1045

[64] J.-P. Putaud et al., "A European aerosol phenomenology—2: chemical characteristics of particulate matter at kerbside, urban, rural and background sites in Europe," *Atmos. Environ.*, vol. 38, no. 16, pp. 2579–2595, May 2004, doi: 10.1016/j.atmosenv.2004.01.041.

1048

[65] X. Querol et al., "Speciation and origin of PM10 and PM2.5 in Spain," *J. Aerosol Sci.*, vol. 35, no. 9, pp. 1151–1172, Sep. 2004, doi: 10.1016/j.jaerosci.2004.04.002.

1049

[66] J. Z. Yu, J. W. T. Tung, A. W. M. Wu, A. K. H. Lau, P. K.-K. Louie, and J. C. H. Fung, "Abundance and seasonal characteristics of elemental and organic carbon in Hong Kong PM10," *Atmos. Environ.*, vol. 38, no. 10, pp. 1511–1521, Mar. 2004, doi: 10.1016/j.atmosenv.2003.11.035.

1053

[67] S. M. McLennan, "Relationships between the trace element composition of sedimentary rocks and upper continental crust," *Geochem. Geophys. Geosystems*, vol. 2, no. 4, pp. 1021–24, Apr. 2001, doi: 10.1029/2000GC000109.

1055

[68] D. Nyanganyura, W. Maenhaut, M. Mathuthu, A. Makarau, and F. X. Meixner, "The chemical composition of tropospheric aerosols and their contributing sources to a continental background site in northern Zimbabwe from 1994 to 2000," *Atmos. Environ.*, vol. 41, no. 12, pp. 2644–2659, 2007.

1058

[69] M. Zheng et al., "Seasonal trends in PM2.5 source contributions in Beijing, China," *Atmos. Environ.*, vol. 39, no. 22, pp. 3967–3976, Jul. 2005, doi: 10.1016/j.atmosenv.2005.03.036.

1060

[70] B. Guinot, H. Cachier, J. Sciare, Y. Tong, W. Xin, and Y. Jianhua, "Beijing aerosol: Atmospheric interactions and new trends," *J. Geophys. Res.*, vol. 112, no. D14, p. D14314, Jul. 2007, doi: 10.1029/2006JD008195.

1062
[71] J. Sandradewi, A. S. H. Prévôt, E. Weingartner, R. Schmidhauser, M. Gysel, and U. Baltensperger, "A study of wood burning and traffic aerosols in an Alpine valley using a multi-wavelength Aethalometer," *Atmos. Environ.*, vol. 42, no. 1, pp. 101–112, Jan. 2008, doi: 10.1016/j.atmosenv.2007.09.034.

[72] C. Pio et al., "OC/EC ratio observations in Europe: Re-thinking the approach for apportionment between primary and secondary organic carbon," *Atmos. Environ.*, vol. 45, no. 34, pp. 6121–6132, Nov. 2011, doi: 10.1016/j.atmosenv.2011.08.045.

[73] N. Bukowiecki et al., "PM10 emission factors for non-exhaust particles generated by road traffic in an urban street canyon and along a freeway in Switzerland," *Atmos. Environ.*, vol. 44, no. 19, pp. 2330–2340, Jun. 2010, doi: 10.1016/j.atmosenv.2010.03.039.

[74] F. Amato et al., "Sources and variability of inhalable road dust particles in three European cities," *Atmos. Environ.*, vol. 45, no. 37, pp. 6777–6787, Dec. 2011, doi: 10.1016/j.atmosenv.2011.06.003.

[75] J. Sternbeck, Å. Sjödin, and K. Andréasson, "Metal emissions from road traffic and the influence of resuspension—results from two tunnel studies," *Atmos. Environ.*, vol. 36, no. 30, pp. 4735–4744, Oct. 2002, doi: 10.1016/S1352-2310(02)00561-7.

[76] A. Thorpe and R. M. Harrison, "Sources and properties of non-exhaust particulate matter from road traffic: a review," *Sci. Total Environ.*, vol. 400, no. 1–3, pp. 270–282, Aug. 2008, doi: 10.1016/j.scitotenv.2008.06.007.

[77] J. G. Watson and J. C. Chow, "Source characterization of major emission sources in the Imperial and Mexicali Valleys along the US/Mexico border," *Sci. Total Environ.*, vol. 276, no. 1–3, pp. 33–47, Aug. 2001, doi: 10.1016/S0048-9697(01)00770-7.

[78] J. A. Adepetu, O. I. Asubiojo, F. Y. Iskander, and T. L. Bauer, "Elemental composition of Nigerian harmattan dust," *J. Radioanal. Nucl. Chem.*, vol. 121, no. 1, pp. 141–147, 1988, Accessed: Aug. 24, 2012. [Online]. Available: http://www.springerlink.com/index/T3631906PT6M2428.pdf

[79] S. M. Almeida, C. A. Pio, M. C. Freitas, M. A. Reis, and M. A. Trancoso, "Source apportionment of fine and coarse particulate matter in a sub-urban area at the Western European Coast," *Atmos. Environ.*, vol. 39, no. 17, pp. 3127–3138, Jun. 2005, doi: 10.1016/j.atmosenv.2005.01.048.

[80] J. G. Watson, "Desert Research Institute, Protocol for Applying and Validating the CMB Model for PM2.5 and VOC. US Environmental Protection Agency, Air Quality Modeling Group," 2004.

[81] C. Liousse and C. Galy-Lacaux, "Urban pollution in West Africa," *La Météorologie*, no. 71, 2010, [Online]. Available: http://hdl.handle.net/2042/37377

[82] J. H. Seinfeld and S. Pandis, "Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 2nd Edition," *Wiley.com*, 2006. https://www.wiley.com/en-ae/Atmospheric+Chemistry+and+Physics%3A+From+Air+Pollution+to+Climate+Change%2C+2nd+Edition-p-9780471720188 (accessed Jun. 06, 2018).

[83] N. J. Pekney et al., "Major Source Categories for PM2.5 in Pittsburgh using PMF and UNMIX," *Aerosol Sci. Technol.*, vol. 40, no. 10, pp. 910–924, 2006, doi: 10.1080/02786820500380271.

[84] S. G. Bhanuprasad, C. Venkataraman, and M. Bhushan, "Positive matrix factorization and trajectory modelling for source identification: A new look at Indian Ocean Experiment ship observations," *Atmos. Environ.*, vol. 42, no. 20, pp. 4836–4852, 2008, doi: 10.1016/j.atmosenv.2008.02.041.

[85] E. D. Goldberg, *H.M. Hill (Ed.), Composition of sea Water, Comparative Oceanography*. in Vol.2 of the seaWiley. New York, 1963.

[86] F. J. Millero, R. Feistel, D. G. Wright, and T. J. McDougall, "The composition of Standard Seawater and the definition of the Reference-Composition Salinity Scale," *Deep Sea Res. Part Oceanogr. Res. Pap.*, vol. 55, no. 1, pp. 50–72, Jan. 2008, doi: 10.1016/j.dsr.2007.10.001.

[87] Z. Yuan, A. K. H. Lau, H. Zhang, J. Z. Yu, P. K. K. Louie, and J. C. H. Fung, "Identification and spatiotemporal variations of dominant PM10 sources over Hong Kong," *Atmos. Environ.*, vol. 40, no. 10, pp. 1803–1815, Mar. 2006, doi: 10.1016/j.atmosenv.2005.11.030.

[88] J. O. Allen, P. R. Mayo, L. S. Hughes, L. G. Salmon, and G. R. Cass, "Emissions of size-segregated aerosols from on-road vehicles in the Caldecott tunnel," *Environ. Sci. Technol.*, vol. 35, no. 21, pp. 4189–4197, Nov. 2001.
[89] J. M. Pacyna and E. G. Pacyna, "An assessment of global and regional emissions of trace metals to the atmosphere from anthropogenic sources worldwide," *Environ. Rev.*, vol. 9, no. 4, pp. 269–298, Dec. 2001, doi: 10.1139/a01-012.

1105

1106

[90] C. Petaloti, A. Triantafyllou, T. Kouimtzis, and C. Samara, "Trace elements in atmospheric particulate matter over a coal burning power production area of western Macedonia, Greece," *Chemosphere*, vol. 65, no. 11, pp. 2233–2243, Dec. 2006, doi: 10.1016/j.chemosphere.2006.05.053.

1107

1108

1109

[91] J. H. Lee, Y. P. Kim, K.-C. Moon, H.-K. Kim, and C. B. Lee, "Fine particle measurements at two background sites in Korea between 1996 and 1997," *Atmos. Environ.*, vol. 35, no. 4, pp. 635–643, 2001, doi: 10.1016/S1352-2310(00)00378-2.

1110

1111

[92] D. D. Cohen, J. Crawford, E. Stelcer, and V. T. Bac, "Characterisation and source apportionment of fine particulate sources at Hanoi from 2001 to 2008," *Atmos. Environ.*, vol. 44, no. 3, pp. 320–328, Jan. 2010, doi: 10.1016/j.atmosenv.2009.10.037.

1112

1113

1114

1115

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

1116

1117

1118