Next Article in Journal
Smartphone-Based Cancer and Obesity Prevention Education Program for Chinese Women (SCOPE): A Pilot RCT
Next Article in Special Issue
Effect of Microplastic Particles on the Rheological Properties of Human Saliva and Mucus
Previous Article in Journal
Environmental Exposure Science and Human Health
Previous Article in Special Issue
Comment on Abdul Jabbar et al. Air Quality, Pollution and Sustainability Trends in South Asia: A Population-Based Study. Int. J. Environ. Res. Public Health 2022, 19, 7534
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment and Source Apportionment of Metals on Atmospheric Particulate Matter in a Suburban Background Area of Gran Canaria (Spain)

by
Yumara Martín-Cruz
1,* and
Álvaro Gómez-Losada
2
1
Independent Researcher, 35480 Agaete, Spain
2
Departament of Quantitative Methods, Universidad Loyola Andalucía, Avda de las Universidades s/n, 41704 Dos Hermanas, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(10), 5763; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20105763
Submission received: 28 February 2023 / Revised: 1 May 2023 / Accepted: 4 May 2023 / Published: 9 May 2023
(This article belongs to the Special Issue Impact of Particulate Matter on the Environment and Health)

Abstract

:
Concentration levels of 11 heavy metals were analyzed in PM 10 and PM 2.5 samples from a suburban area frequently affected by Saharan dust in which is located a school. The heavy metals risk assessment was carried out by the 2011 U.S. Environmental Protection Agency method, estimating the chronic and carcinogenic hazard levels both in adults and children. The highest level of chronic hazard was reached for Cr, with values of approximately 8 (PM 10 , adulthood), 2 (PM 10 , childhood) and 1.5 (PM 2.5 , adult age), significantly exceeding the limit value (equal to 1). Regarding the carcinogenic risk level, it was also high for Cr, with values between 10 3 and 10 1 for both study populations and particle size. For the rest of the studied metals, no health risk levels of concern were obtained. The positive matrix factorization method was used for the estimation of heavy metal emission sources apportionment. Non-exhaust vehicle emissions were the main source of Cr emissions under PM 2.5 , while industrial processes were the main source for PM 10 . Mineral dust and marine aerosol were common emission sources for both particles sizes—with different contributions. Vehicle emissions, construction and agricultural activities were the main emission sources for PM 10 , and fossil fuel combustion, road dust re-suspension and ammonium sulfate were for PM 2.5 . The results obtained in this study support the need to continue applying mitigation measures in suburban areas which are affected by nearby anthropogenic emissions, causing the consequent emission of materials hazardous to human health.

Graphical Abstract

1. Introduction

Particulate matter (PM) is a harmful atmospheric pollutant with critical effects on human health, causing respiratory and cardiovascular diseases, sleep disorders [1], dermal cancer [2] and premature death. Furthermore, PM can impact the climate [3,4], interfering in the Earth’s radiation balance by scattering and absorbing solar radiation and acting as cloud condensation nuclei, affecting the biogeochemical cycles of marine organisms [5], acidifying oceans and lakes, and contributing to acid rain. Regarding the latter effect, PM can cause deterioration and blackening of materials, with associated economic costs.
PM is a complex composite of a wide variety of chemical species, such as heavy metals, which are not metabolized by the human body and accumulate in the soft tissues [6], causing serious respiratory and cardiovascular diseases, cancer and premature death [7]. Additionally, these metals can be absorbed through dermal contact and oral ingestion [8,9]. Children’s exposure to heavy metals leads to neurodevelopmental deficits, skeletal damage and adverse effects on sexual function and fertility [10].
Heavy metals can be emitted from several sources. Natural sources include mineral dust (Al, Fe, Si, Ti, Mn) and volcanic eruptions [11,12]. Among anthropogenic sources are non-exhaust (Ba, Fe, Cu, Cd, Zn) and exhaust vehicle emissions (Ca, Mg, Ni, Ba, Cu, Cd, Zn), biomass burning and incineration processes (K, Cd) [13], fossil fuel combustion (Ni, V) [14], metallurgical industries (Cr, Cu, Mn, Pb, Sn, Zn, Co) [15], mining activities (As, Cd, Cu, Ni, Pb, Zn) [16] and agricultural activities (Cu, As, Cd, Mn, Pb, Zn) [17,18]. Apportionment of these emission sources is crucial for the development and implementation of air quality policies; it is also to report the effects of mitigation measures [19].
The assessment of the risks of air pollutants on human health has been extensively studied in urban environments, being much scarcer in peripheral areas to urban centers. This study aims to contribute to fill in this gap by assessing the health risks and estimating the contribution of emission sources from 11 atmospheric heavy metals in a suburban area in the southeast of Gran Canaria island (Spain).

2. Materials and Methods

2.1. Area of Study

The Canary Islands (Spain) is an archipelago located in the Atlantic Ocean which experiences intense Saharan outbreaks due to its proximity to the northeastern coast of Africa. This archipelago has eight islands, with Gran Canaria as the second most populated island. The area of study was Taliarte, an area with a suburban background in the southeast of Gran Canaria. The main sources of anthropogenic emissions in this area are vehicle (local and highway) and maritime traffic, together with agricultural activities. As for natural PM emission sources, the area of study is affected by Saharan dust and marine aerosols due to its proximity to the coast (Figure 1).
The area of study has a mild climate, with an average annual temperature of around 21 C and relative humidity of approximately 69%. Figure 2 shows the temporal evolution of the temperature and relative humidity, as well as the frequency of wind direction (%) according to the wind octants between 2010 and 2021.

2.2. Sampling

The sampling point was at the University of Las Palmas de Gran Canaria’s Parque Científico Tecnológico building.In total, 120 samples were collected during the second semester of 2017, using high volume samplers, according to the UNE-EN 1234 normalized method. Before and after sampling, quartz filters (Whatman) were conserved at constant relative humidity (45–50%) and temperature (20–25 C) conditions for at least 24 h. Samples were collected on random days, except days with Saharan dust outbreaks. In these cases, the sampling was continuous from the beginning to the end of the dust event. For the forecast of these outbreaks the Barcelona Supercomputing Center (BSC-DREAM-8b and NMMB/BSC-Dust) and Skiron University numeric prediction models were consulted.

2.3. Sample Treatment and Chemical Analysis

A quarter of quartz sample filter was used in both PM 10 and PM 2.5 to analyze 11 heavy metals—Al, Ba, Cd, Cr, Cu, Fe, Mn, Ni, Ti, V and Zn. Acid digestion (2.5 mL HNO 3 :5 mL HF:2.5 mL HClO 4 ) was carried out, and the samples were analyzed using a graphite furnace atomic absorption spectrometer. Two complementary analyses were performed in addition to the analysis of the eleven heavy metals in order to facilitate the subsequent identification and contribution of emission sources using positive matrix factorization (PMF): (1) Ca, Mg and K analysis (acid extraction and later analysis using a graphite furnace atomic absorption spectrometer), and (2) analysis of soluble species. To the latter, two-quarters of the sample filter were extracted with deionized water (resistivity: 18.2 M Ω -cm ) and ultrasonicated for 60 min. Ca 2 + , K + , Mg 2 + and Na + were analyzed with graphite furnace atomic absorption spectrometry, Cl with the argentometric method. Molecular absorption UV–Visible spectrometry (APHA 4500-NO-3B method) was used for NO 3 , and the turbidimetric method (EPA 9038) for SO 4 2 . Finally, NH 4 + was analyzed using the nesslerization method (EPA 350.2).

2.4. Health Risk Assessment

The health risks were estimated using the method developed by the USEPA 2011 (U.S. Enviromental Protection Agency), analyzing the daily chemical intake through oral ingestion, exposure concentration through inhalation, and dermal absorption dose through dermal contact [20,21]. Then, non-carcinogenic risks (e.g., asthma, rhinitis, pulmonary fibrosis, among others) and carcinogenic risk were estimated for childhood and adulthood. For this, the hazardous index (HI) and total carcinogenic risk (TCR) for each exposure pathway were calculated. A detailed explanation of these parameters is given in Appendix A. The limit values for determining the chronic and carcinogenic risk level are shown in Table 1 [22].

2.5. Source Apportionment

To estimate the contribution of anthropogenic sources to the concentrations of the metals being studied both in PM 10 and PM 2.5 , enrichment factors and positive matrix factorization were used. These are described next.

2.5.1. Enrichment Factors (EFs)

Enrichment factors were calculated according to Equation (1) [23,24]:
E F = ( X / R ) s a m p l e ( X / R ) c r u s t a l
where X is the concentration of the metal being studied and R is that of the reference metal, this being the metal with the lowest anthropogenic contribution at sampling area. An EF lower than 10 indicates that the metal has a predominately natural origin, while if EF exceeds 100, the metal is emitted only by anthropogenic sources. EF values between 10 and 100 reveal a confluence of sources, with those of anthropogenic origin predominating [25].
As reference metal, Ti was chosen for crustal analysis. This metal was selected based on the study of the ratio ([Al]/[Ti]) s a m p l e vs. ([Al]/[Ti]) c r u s t a l . If the first term is higher than the second term, Al is enriched, taking the other metal as reference. Standard composition of basalt from the Geochemical Database for Reference Materials (GEOREM) [26] was used to determine the (X/R) c r u s t a l value.

2.5.2. Positive Matrix Factorization (PMF)

PMF is a multivariate factorial analysis used to identify and quantify the sources of chemical species’ emissions at a receptor point using the fingerprint of these sources. It is the most widely used receptor-oriented model because it does not require much previous information on emission sources, enabling it to be applied to numerous studies of source apportionment [19]. This analysis, like any receptor model, is based on the principle of mass conservation between concentrations of measured species and source profiles (Equation (2)) [27,28].
X i j = k = 1 p ( g i k · f k j ) + e i j
The matrix of concentrations X i j corresponds to a set of measured data on i samples of j chemical species, and is divided into two sub-matrices for p sources: the g i k matrix shows the contribution of k source at i sample, and the matrix f j k contains the concentration of j species for k source (source profiles). e i j is the residue corresponding to each datum value. A weighted least-squares method is used to obtain the g and f values that minimize the Q function (Equation (3)), including the data uncertainties ( u i j ) in the input matrix. n and m indicate the number of samples and chemical species analyzed, respectively [29].
Q = i = 1 n j = 1 m e i j u i j 2
The PMF was carried out with EPA PMF 5.0 software, taking into account the concentrations of species studied and associated uncertainty around input data. Among sources of error on any measurement, the error associated with the analytical procedure is the most important [30]. The calculation procedure applied in this study to determine these uncertainties as well as the procedure performed with this software, is explained in Appendix B, Appendix C and Appendix D.

3. Results and Discussion

3.1. Descriptive Statistical Analysis

The average PM 10 concentration during the sampling period was 64.95 ± 36.11 μ g·m 3 , exceeding the limit value permitted under Directive 2008/50/EC (50 μ g·m 3 ), and reaching maximum values that quadrupled this limit. These results were the consequence of intense Saharan dust outbreaks over the islands during the sampling period, together with the presence of intense marine aerosol in the zone. In the case of PM 2.5 , the average concentration was 21.43 ± 9.71 μ g·m 3 , slightly higher than the limit established by Directive 2008/50/EC (20 μ g·m 3 ), where the effect of Saharan dust was less important. As for heavy metals, the statistical summaries for both PM 10 and PM 2.5 at the sampling site are shown in Table 2.
As mentioned, the study area is strongly affected by Saharan dust outbreaks, impacting the concentration level of some heavy metals. This effect was observed for both crustal metals such as Al, with a coefficient of variation of 136% in contrast to 68% of variability in PM 2.5 , and anthropogenic metals such as V, with a coefficient of variation equal to 84% in PM 10 and 58% in PM 2.5 . Figure 3 and Figure 4 show the behavior of each metal in the three Saharan dust regimes considered (absence, medium and high) for PM 10 and PM 2.5 , respectively. The Wilcoxon pairwise test was used for the comparison of means, since a pairwise behavior of non-parametric groups was performed. A significance level of less than 0.05 was considered. The influence of Saharan dust on concentration levels of heavy metals in the fine fraction (PM 2.5 ) was less significant (even in crustal metals such as Al or Mn) repeating the pattern of a greater presence of coarse particles in the dust mineral.
The predominance of fine fraction of crustal metals such as Fe and Al, during the period with non-Saharan events, with PM 2.5 to PM 10 ratios equal to 70% and 62.5%, respectively, could indicate their emission by anthropogenic sources such as exhaust emission from vehicle traffic [13]. Unlike these two metals, Ti and Mn showed a percentage exceeding 50% of the coarse fraction. Possible causes could be the natural local sources and non-exhaust emissions from vehicle traffic such as road dust re-suspension [31]. Two other metals linked to the latter type of metal emission are Cu and Ba, with more than 50% of their composition distributed in the coarse fraction. Both metals could have been emitted by brake pad wear, since they are typically used as filler materials in the form of barium sulfate and copper fibers [32].
Figure 5 shows the mass percentages of heavy metals on PM 10 concentration. The highest contribution was obtained during intense Saharan dust episodes, due to the increase in the concentration of crustal metals such as Al, Fe and Ti. As for the anthropogenic metals such as Zn, Cr or Cd, a decrease in mass contribution was observed under these conditions as their concentration was not affected during dust outbreaks. Unlike PM 10 , the lowest percentage mass contribution for PM 2.5 (Figure 6) was observed during intense Saharan dust outbreaks. This seems to be due to the metal composition of PM 2.5 , which was unaffected by the intensity of these Saharan dust events. Due to the mineral origin of this dust, the significant increase in metal concentration occurred in both PM 10 2.5 fractions.

3.2. Health Risk Assessment

Figure 7 shows the results obtained after the chronic risk assessment for each studied heavy metal. The main exposure pathway was inhalation, with HQ values for ingestion and dermal contact equal or very close to zero. Inhalation risk levels for Cr were the highest, significantly exceeding (up to 8 times) the USEPA limit equal to 1. The risk of causing chronic effects of Cr was reported by previous studies, such as the one conducted by [31] in a marginal area of an urban core in China. Reference [9] performed risk assessments in an urban area with high population density in China. In this study, an HQ value greater than 1 for Cr was obtained in childhood but not in adulthood. Unlike this present study, the main exposure pathway was ingestion. Likewise, in a study by [21] in the urban area of Kitahyushu (Japan), the HI value for Cr was less than 1. As for the other metals analyzed in this present research, the HQ i n h values were less than 1, indicating the absence of chronic risks, except for Mn (PM 10 , adulthood). In the study by [10], high values were also obtained for Mn, although at PM 2.5 .
Regarding carcinogenic risks, the main exposure pathway also was inhalation (Table 3 and Table 4). In this case, only Cd, Ni and Cr were considered, as they are classified as carcinogenic to humans (group I) by the International Agency for Research on Cancer. Previous studies [31] demonstrated the carcinogenic effect of the former metals, with TCR values higher than 10 6 (USEPA limit). Conversely, in this study, TCR i n h values for both Ni and Cd were around this limit, especially in adulthood, considered a low risk level. A similar situation was obtained in the study by [9]. In the case of Cr, high carcinogenic risk was reached both in the childhood and adulthood study in PM 2.5 and PM 10 . These carcinogenic risks were also reported by other studies such as [8].

3.3. Source Apportionment

3.3.1. Enrichment Factors

EFs c r u s t a l for metal species studied are shown in Figure 8. Cd, Cr, Cu, Ni and Zn were emitted by anthropogenic sources in PM 10 and PM 2.5 , with values higher than 100. The EF c r u s t a l for Ca was greater than 10, implying that this metal was also emitted by anthropogenic sources, such as construction activities and vehicle traffic. In the case of Na, the EF c r u s t a l was remarkable, with a value higher than 100. This situation showed the influence of anthropogenic sources, such as exhaust emissions from vehicle traffic, since these high values were also observed for PM 2.5 . The case of Al and V is worth highlighting. In the PM 10 , these metals showed an EF c r u s t a l very close to 10 (9.8 and 12.9, respectively), which corresponds to a predominance of crustal sources combined with anthropogenic sources. Analyzing the values achieved in PM 2.5 for these metals, the EFs c r u s t a l were higher (14.1 for Al and 37.8 for V), implying that both showed a greater contribution of anthropogenic sources in the fine fraction.

3.3.2. Positive Matrix Factorization

After application of EPA PMF 5.0, six and five factors were obtained as the optimal solution for PM 10 and PM 2.5 , respectively. The high values of the correlation factor for the scatter plots between the observed and predicted PM concentrations (0.94 for PM 10 and 0.80 for PM 2.5 ) and the lack of collinearity among factors verified the suitability of the solution adopted.
  • PM 10 emission sources
  • The percentage contribution of each factor to the concentration of the studied metal species is shown in Table 5. The factors are listed from highest to lowest contribution to the concentration of PM and whose sum is equal to 100.
  • The first factor showed high contributions of Al, Ti and Mn (>60%) and, to a lesser extent, contributions of Fe, V, Ba and K (between 20 and 40%), which have a predominantly mineral origin [33]. According to the change in the contribution of this factor, this crustal matter was dominated by the Saharan dust outbreaks, with important contribution values observed during these events. Likewise, high contributions of Cl and NO 3 , possibly due to the presence of halite and sodium nitrate (aged marine aerosol) were observed. These two ions were also explained by the second factor, which may be associated with construction activities carried out during the sampling period. Both Cl and NO 3 are typically used as Portland cement additives to accelerate setting times, among other benefits [34]. The presence of Ba, Cu, Mn and Ni in this factor could be due to exhaust emissions from vehicles used on construction sites, used as diesel additives [35] and emitted as combustion sub-products [31,36]. Finally, Cl was also explained by the third factor, as was the Na + and a slight contribution of SO 4 2 , which suggests that this factor could corresponds to marine aerosol emission [30], due to the sampling site’s proximity to the coast.
  • The fourth factor is entirely accounted for by vehicle emissions, with a predominance of non-exhaust emissions. The intense Fe contribution in this factor can be explained by wear on brake pads, since it is their main component [32]. Other elements explained by this factor and associated with these emissions are Ba (used as filler material and considered a tracer [37]), Cu (used in reinforcement fibers [38]) and Al and Cr (employed in abrasives [38]). The second type of non-exhaust emissions corresponds to road dust re-suspension from the circulation of vehicles, with high percentages of Ca [39], K and V. The third type is due to tire wear, which could explain the contribution of this factor to the concentration of Zn, used as a vulcanization agent and the main source of ambient Zn [39]. Exhaust emissions may also be the cause of emission of Cd and Zn due to the combustion of lubricating oil [40,41] and of Mn used as a catalyst (Mn 2 O 3 ) and petrol additive [42].
  • The fifth factor was referred as agricultural activities and traffic, and is the main source of NH 4 + emission from the use of fertilizers and manure [43]. These agricultural activities could also be responsible for high contributions of Cd and K. The proximity of the highway accounts for the contributions of metals such as Ba, Cu, Zn and Al.
  • The sixth factor explained a significant percentage of Cr and, to a lesser extent, Ni, which would indicate a possible industrial source [15]. Likewise, the mechanical abrasion and sanding works in the port [44] and emissions from the aircraft engines that continuously circulate in the area could also be considered a source of Cr emission.
  • PM 2.5 emission sources
  • The percentage contributions for each of the factors obtained in PM 2.5 are shown in Table 6. As in the case of PM 10 , the total sum of the contributions is equal to 100. The high Ni loading in addition to V and SO 4 2 in the first factor were indicative of emissions from fossil fuel combustion [45]. The low V/Ni ratio, equal to 0.17, showed an additional source of Ni, such as emissions from motor vehicles, particularly, the diesel-powered ones. This fact was confirmed by the presence of other metals such as Ba, Cu and Mn.
  • The second factor may correspond to the road dust re-suspension due to high Fe and Ca presence, while the third factor was attributed to mineral dust, as it revealed high percentages of Al, Mn and Ti, which have a predominantly crustal origin, as already commented in the case of PM 10 .
  • The fourth factor was characterized by a high NH 4 + loading, more than 80%. As already mentioned for PM 10 , agricultural and livestock activities around of the sampling site were two significant sources of this ion. The NH 4 + /SO 4 2 molar ratio was equal to 1.6, indicating that the total ammonium was neutralized by sulfates, since the molar ratio was 2:1. The slight excess of the latter anions could be due to emissions from vehicle traffic, since the presence of Ba, Ti and K was also observed. Following [46], traffic emissions are considered a major sources of these inorganic ions in urban areas, so that the factor was cataloged as “ammonium sulfate + traffic”.
  • Marine aerosol also contributed to the PM 2.5 concentration, but very slightly, constituting the fifth factor and the one with the lowest percentage. This factor was characterized by high Na, Mg and K loading. It is a marine aerosol polluted by anthropogenic sources such as traffic, as indicated by the high contributions of Cd and Al.

4. Conclusions

In this study, PM 10 and PM 2.5 samples were collected in the University of Las Palmas de Gran Canaria’s Parque Científico Tecnológico building (Gran Canaria, Spain). This is a coastal and suburban area with frequent Saharan dust outbreaks. Air quality of this area was acceptable, with PM 10 and PM 2.5 concentration levels below Directive 2008/50/EC limits under normal conditions. Chemical characterization of the samples was carried out and the health risk assessment and sources apportionment of eleven heavy metals were performed. The highest risk levels were obtained for Cr, both for its chronic and carcinogenic effects, exceeding the established limits by USEPA. Cr emissions could come from the industrial activities states near the sampling point and from the wharf activities. Likewise, the proximity of a school to the sampling point reveals an increased exposure of children to these emissions due to outdoor activities. Based on this, the results obtained in this research support the need to implement mitigation measures in areas that are highly industrialized and seriously affected by emissions of substances that pose a health hazard to the population.

Author Contributions

Conceptualization, Y.M.-C.; Data curation, Y.M.-C.; Investigation, Y.M.-C.; Methodology, Y.M.-C.; Supervision, Á.G.-L.; Writing—original draft, Y.M.-C.; Writing—review and editing, Y.M.-C. and Á.G.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Parameters for Health Risks Assessment

Appendix A.1. Calculation Equations

This section shows the equations used for the calculation coefficients needs for the health risk assessment, as well as the significance of each parameter and the adopted value for each them.
C D I i n g = C 95 % · I R · E F · E D · C F B W · A T
E C i n h = C 95 % · E T · E F · E D A T n
D A D d e r = C 95 % · S A · A F · A B S · E F · E D · C F B W · A T
C 95 % = e x p l n X ¯ + 0.5 · S l n x 2 + S l n X · H 0.95 n 1
H Q i n g = C D I i n g R f D o
C R i n g = C D I i n g · S F o
H Q i n h = E C i n h R f C i · 1000
C R i n h = E C i n h · I U R
H Q d e r = D A D d e r R f D o · G I A B S
C R d e r = D A D d e r · S F o G I A B S

Appendix A.2. Parameter’s Values

Table A1. Required parameters (name and value) for health risks assessment.
Table A1. Required parameters (name and value) for health risks assessment.
AcronymNameUnitValue
Common parameters
EFExposure frequencyday·year 1 350
ETExposure timeh·day 1 24
ABSDermal absorption factor-0.001
ATAveraging lifetimedays365 × 70
(for carcinogens)
ATnAverage lifetimehours365 × 70 × 24
(for carcinogens)
CFUnit conversion factorKg·g 1 10 × 10 6
Parameters dependent on life stage:ChildhoodAdulthood
IRIngestion ratemg·day 1 200100
EDExposure durationyear624
BWBody weightKg1770
ATAveraging lifetimedaysED × 70
(for non carcinogens)
ATnAverage lifetimehoursED × 70 × 24
(for non carcinogens)
SASkin surface areacm 2 28005700
AFSkin adherence factormg·cm 2 0.20.07
Parameters dependent of heavy metalsBaCdCrCuMnNiVZn
RfD 0 Oral reference dosesmg·(kg·day) 1 0.20.0010.0030.040.0240.020.0070.3
RfC i Inhalation reference concentrationmg·m 3 0.00050.000010.000010.040.000050.000090.00010.3
GIABSGastrointestinal absorption factor-0.070.0250.002510.040.040.0261
Sf 0 Oral slope factor(mg·(kg·day) 1 ) 1 -6.30.5-----
IURInhalation unit risk( μ g·m 3 ) 1 --0.084-----

Appendix B. Calculation Procedure with EPA PMF 5.0

Once the input data were entered, the signal to noise (S/N) was studied for each variable, indicating if the variability of the measurements is real or within data noise [47]. This parameter can be interpreted as the relationship between concentrations of studied chemical species and uncertainties, allowing for classification of studied elements in the follow three categories. Species with S/N ≥ 2 were considered as “strong variables”, which indicates good quality data. Species with S/N between 0.2 and 2 were categorized as “weak variables”, which can indicate insufficient variations in the total concentration by them, contributing to noise on obtained results. In these cases, threefold uncertainty was employed with them by the model. Finally, species with S/N < 0.2 were labeled as “bad variables” and were excluded from analysis. In the majority of species in this study, S/N > 2 were obtained, except in the following cases: Cu and Cd were considered weak both in PM 10 and PM 2.5 ; Zn was referred to as weak in PM 10 and bad in PM 2.5 ; finally, Ti and Ni were considered as weak only in PM 2.5 . Values of S/N for each variable are shown in Table A2.
As already mentioned, the PMF model is applied to identify the number of emission sources in a specific receptor area. To determine the optimal number of factors, the following criteria were applied. Firstly, the value of Q theoretical, calculated as nm − p(n + m),was compared with the values of Q true and Q robust, which were parameters provided by the model. A good fit is obtained when these three parameters are similar to each other [48]. Secondly, scaled residuals were assessed. They are defined as the ratio of model residuals, e i j and uncertainty associated by the input data. For good fit, these residuals should be symmetrically distributed within a range of −3 to +3. Thirdly, dispersion diagrams [49] between observed value and predicted value and G diagrams were evaluated. These last indicate the existence of collinearity among sources (or factors) [50]. Several proofs were carried out, considering from 5 to 8 factors and different initial additional uncertainties (0–15%).

Appendix C. Uncertainties Calculation for PMF

In this section, the calculation procedure of required uncertainties to PMF analysis development is explained.

Appendix C.1. Particulate Matter

Associated uncertainties to PM 10 and PM 2.5 were estimated with standard UNE-EN 12341:2015 using the following equation:
u T = u f d 2 + u m 2 V 2 + C 2 · u j 2 100 2
u f d : field uncertainty (0.8 μ g·m 3 and 2.3 μ g·m 3 for PM 10 and PM 2.5 , respectively).
u m : mass uncertainty. Considering the associated uncertainty to the balance and the buoyancy effect, a value of 2.5 μ g·m 3 was considered.
u f w : 6.3% air flow uncertainty.
C , V : particulate matter concentration ( μ g·m 3 ) and sampled air volume (m 3 ), respectively, for PM 10 and PM 2.5

Appendix C.2. Chemical Species

Associated uncertainties to studied chemical species were estimated using one of the three following equations, according to the ratio of species concentration ( X i j ) to detection limit ( L D j ):
u i j ( X i j < L D j ) = X i j + 2 / 3 L D j
u i j ( L D j < X i j < 3 L D j ) = 0.2 X i j + 2 / 3 L D j
u i j ( X i j > 3 L D j ) = 0.1 X i j + 2 / 3 L D j
0.1 and 0.2 are empirical coefficients which were obtained by [51].

Appendix D. Additional Results from PMF Analysis

The following table shows the categorization of the studied variables according to the signal to noise ratio (S/N) for PM 10 and PM 2.5 .
Table A2. S/N ratio value for PM 10 and PM 2.5 .
Table A2. S/N ratio value for PM 10 and PM 2.5 .
SpecieCategoryS/NCategoryS/N
PM 10 PM 2 . 5
PMStrong10.06Strong7.1
AlStrong6.8Strong5.3
BaStrong2.8Strong2.0
CdWeak0.9Weak0.5
CrStrong4.1Weak1.7
CuWeak1.5Weak0.6
FeStrong4.4Strong2.3
MnStrong5.5Strong2.4
NiStrong3.2Weak1.9
TiStrong4.2Weak1.9
VStrong6.9Strong6.9
ZnWeak0.2Bad0.1
Mg 2 + Strong5.7
Na + Strong7.9Strong6.4
Cl Strong8.5
SO 4 2 Strong8.8Strong9.0
NO 3 Strong8.8Strong9.0
NH 4 + Strong8.3Strong8.9
CaStrong6.6Strong4.4
KStrong8.0Strong5.9
MgStrong4.5

References

  1. Li, L.; Zhang, W.; Xie, L.; Jia, S.; Feng, T.; Yu, H.; Huang, J.; Qian, B. Effects of atmospheric particulate matter pollution on sleep disorders and sleep duration: A cross-sectional study in the UK biobank. Sleep Med. 2020, 74, 152–164. [Google Scholar] [CrossRef] [PubMed]
  2. Kim, K.E.; Cho, D.; Park, H.J. Air pollution and skin diseases: Adverse effects of airborne particulate matter on various skin diseases. Life Sci. 2016, 152, 126–134. [Google Scholar] [CrossRef] [PubMed]
  3. Kaufman, Y.J.; Tanré, D.; Boucher, O. A satellite view of aerosols in the climate system. Nature 2002, 419, 215–223. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, H.; Hu, J.; Kleeman, M.; Ying, Q. Source apportionment of sulphate and nitrate particulate matter in the Eastern United States and effectiveness of emission control programs. Sci. Total Environ. 2014, 490, 171–181. [Google Scholar] [CrossRef]
  5. Li, H.; Chen, Y.; Zhou, S.; Wang, F.; Yang, T.; Zhu, Y.; Qingwie, M. Change of dominant phytoplankton groups in the eutrophic coastal sea due to atmospheric deposition. Sci. Total Environ. 2021, 753, 141961. [Google Scholar] [CrossRef]
  6. Moryani, H.T.; Kong, S.; Du, J.; Bao, J. Health risk assessment of heavy metals accumulated on PM2.5 fractioned road dust from two cities of Pakistan. Int. J. Environ. Res. Public Health 2020, 17, 7124. [Google Scholar] [CrossRef]
  7. Kastury, F.; Smith, E.; Juhasz, A.L. A critical review of approaches and limitations of inhalation bioavailability and bioaccessibility of metal(loid)s from ambient particulate matter or dust. Sci. Total Environ. 2017, 574, 1054–1074. [Google Scholar] [CrossRef]
  8. Aminiyan, M.M.; Baalousha, M.; Aminiyan, F.M. Evolution of human health risk based on EPA modeling for adults and children and pollution level of potentially toxic metals in Rafsanhan road dust: A case study in a semi-arid region, Iran. Environ. Sci. Pollut. Res. 2018, 25, 19767–19778. [Google Scholar] [CrossRef]
  9. Liu, P.; Ren, H.; Xu, H.; Lei, Y.; Shen, Z. Assessment of heavy metals characteristics and health risks associated with PM2.5 en Xi’an, the largest city in northwestern, China. Air Qual. Atmos. Health 2018, 11, 1037–1047. [Google Scholar] [CrossRef]
  10. Zhang, H.; Zhenxing, M.; Huang, K.; Wang, X.; Cheng, L.; Zeng, L.; Zhou, Y.; Jing, T. Multiple exposure pathways and health risk assessment of heavy metal(loid)s for children living in fourth-tier cities in Hubei Province. Environ. Int. 2019, 129, 517–524. [Google Scholar] [CrossRef]
  11. Morillas, H.; Maguregui, M.; Paris, C.; Bellot-Gurlet, L.; Colomban, P.; Madariaga, J.M. The role of marine aerosol in the formation on double sulfate/nitrate salts in plasters. Microchem. J. 2015, 123, 148–157. [Google Scholar] [CrossRef]
  12. Rodriguez, S.; Calzolai, G.; Chiari, M.; Nava, S.; García, M.I.; López-Solano, J.; Marrero, C.; López-Darias, J.; Cuevas, E.; Alonso-Pérez, S.; et al. Rapid changes of dust geochemistry in the Saharan Air Layer linked to sources and meteorology. Atmos. Environ. 2020, 223, 117186. [Google Scholar] [CrossRef]
  13. Saggu, G.S.; Mittal, S.K. Source apportionment of PM10 by positive matrix factorization model at a source region of biomass burning. J. Environ. Manag. 2020, 266, 110545. [Google Scholar] [CrossRef] [PubMed]
  14. Manousakas, M.; Diapouli, E.; Papaefthymiou, H.; Migliori, A.; Karydas, A.G.; Padilla-Alvarez, R.; Bogovac, M.; Kaiser, R.B.; Jaksic, M.; Bogdanovic-Radovic, I.; et al. Source apportionment by PMF on elemental concentrations obtained by PIXE analysis of PM10 samples collected at the vicinity of lignite power plants and mines in Megalopolis, Greece. Nucl. Instrum. Methods Phys. Res. Sect. Beam Interact. Mater. Atoms 2015, 349, 114–124. [Google Scholar] [CrossRef]
  15. Sadeghi, B.; Choi, Y.; Yoon, S.; Flynn, J.; Kotsakis, A.; Lee, S. The characterization of fine particulate matter downwind of Houston: Using integrated factor analysis to identify anthropogenic and natural sources. Environ. Pollut. 2020, 262, 114345. [Google Scholar] [CrossRef] [PubMed]
  16. Blondet, I.; Schreck, E.; Viers, J.; Casas, S.; Jubany, I.; Bahí, N.; Zouiten, C.; Dufréchou, G.; Freydier, R.; Galy-Lacaux, C.; et al. Atmospheric dust characterization in the mining distric of Cartagena-La Unión, Spain: Air quality and health risks assessment. Sci. Total Environ. 2019, 693, 133496. [Google Scholar] [CrossRef]
  17. Kumar, A.; Chauchan, A.; Arora, S.; Tripathi, A.; Mohammed, S.; Alghanem, S.; Khan, K.A.; Ghramh, H.A.; Özdemir, A.; Ansari, M.J. Chemical analysis of trace metal contamination in the air of industrial area of Gajraula (U.P.), India. J. King Saud Univ.-Sci. 2020, 32, 1106–1110. [Google Scholar] [CrossRef]
  18. Tyagi, V.; Gurjar, B.; Joshi, N.; Kumar, P. PM10 and heavy metals in suburban and rural atmospheric environments of Norhthern India. Pract. Period. Hazard. Toxic Radioact. Waste Manag. 2012, 16, 175–182. [Google Scholar] [CrossRef]
  19. Sun, X.; Wang, H.; Guo, Z.; Lu, P.; Song, F.; Liu, L.; Liu, J.; Rose, N.L.; Wang, F. Positive matrix factorization on source apportionment for typical pollutants in different environmental media: A review. Environ. Sci. Process. Impacts 2020, 22, 239–255. [Google Scholar] [CrossRef] [PubMed]
  20. MohseniBandpi, A.; Eslami, A.; Ghaderpoori, M.; Shahsavani, A.; Jeihooni, A.K.; Ghaderpoury, A.; Alinejad, A. Health risk assessment of heavy metals on PM2.5 in Tehran air, Iran. Data Brief 2018, 17, 347–355. [Google Scholar] [CrossRef]
  21. Zhang, X.; Eto, Y.; Aikawa, M. Risk assessment and management of PM2.5-bound heavy metals in the urban area of Kitakyushu, Japan. Sci. Total Environ. 2021, 795, 148748. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, R.; Jia, B.; Tian, Y.; Feng, Y. Source-specific health risk assessment of PM2.5-bound heavy metals based on high time-resolved measurement in a chinese megacity: Insights into seasonal and diurnal variations. Ecotoxicol. Environ. Saf. 2021, 216, 112167. [Google Scholar] [CrossRef] [PubMed]
  23. Shelley, R.U.; Morton, P.L.; Langing, W.M. Elemental ratios and enrichment factors in aerosols from the US-GEOTRACES North Atlantic transects. Deep-Sea Res. II 2015, 116, 262–272. [Google Scholar] [CrossRef]
  24. Buck, C.S.; Aguilar-Islas, A.; Marsay, C.; Kadko, D.; Landing, W.M. Trace element concentrations, elemental ratios and enrichment factors observed in aerosol samples collected during the US GEOTRACES eastern Pacific Ocean transect (GP16). Chem. Geol. 2019, 511, 212–224. [Google Scholar] [CrossRef]
  25. Heldari-Farsani, M.; Shirmardi, M.; Goudarzi, G.; Alavi, N.; Ahmadi-Ankali, K.; Zallaghi, E.; Naimabadi, A.; Hashemzadeh, B. The evaluation of heavy metals concentration related to PM10 in ambient air of Ahvaz City, Irane. J. Adv. Environ. Health Res. 2014, 1, 120–128. [Google Scholar]
  26. Jochum, K.P.; Nohl, U.; Herwing, K.; Lammel, E.; Stoll, B.; Hofmann, A.W. GeOReM: A new geochemical database for reference materials and isotopic standards. Geostand. Geonalytical Res. 2007, 29, 333–338. [Google Scholar] [CrossRef]
  27. Hopke, P.K. Review of receptor modeling methods for source apportionment. J. Air Waste Manag. Assoc. 2016, 66, 237–259. [Google Scholar] [CrossRef]
  28. Mircea, M.; Calori, G.; Pirovano, G.; Bellis, C. European guide on air pollution source apportionment for particulate matter with source oriented models and their combined use with receptor models. EUR 30082 EN Publ. Off. Eur. Union Luxemb. 2008, 10, 142–149. [Google Scholar]
  29. Bellis, C.A.; Larsen, B.R.; Amato, F.; El Haddad, I.; Favez, O.; Harrison, R.M.; Hopke, P.K.; Nava, S.; Paatero, P.; Prévôt, A.; et al. European Guide on Air Pollution Source Apportionment with Receptor Models; European Commission, EUR 26080 EN; Publications Office of the European Union: Luxembourg, 2013. [Google Scholar]
  30. Millan-Martinez, M.; Sánchez-Rodas, D.; Sánchez de la Campa, A.; de la Rosa, J. Contribution of anthropogenic and natural sources in PM10 during North African dust events in Southern Europe. Environ. Pollut. 2021, 290, 118065. [Google Scholar] [CrossRef]
  31. Wang, J.M.; Jeong, C.H.; Healy, R.M.; Sofowote, U.; Debosz, J.; Su, Y.; Munoz, A.; Evans, G.J. Quantifying metal emission from vehicle traffic using real world emissions factors. Environ. Pollut. 2021, 268, 115805. [Google Scholar] [CrossRef]
  32. Cai, R.; Zhang, J.; Nie, X.; Tjong, J.; Matthews, D.T.A. Wear mechanism evolution on brake discs for reduced wear and particulate emissions. Wear 2020, 452–453, 203283. [Google Scholar] [CrossRef]
  33. Lokorai, K.; Ali-Khodja, H.; Khardi, S.; Bencharif-Madani, F.; Naidja, L.; Bouziane, M. Influence of mineral dust on the concentration and composition of PM10 in the city of Constantine. Aeolian Res. 2021, 50, 100677. [Google Scholar] [CrossRef]
  34. Dorn, T.; Blask, O.; Stephan, D. Acceleration of cement hydration—A review of the working mechanisms, effects on setting time, and compressive strength development of accelerating admixtures. Contruction Build. Mater. 2022, 323, 126554. [Google Scholar] [CrossRef]
  35. Hosseinzadeh-Bandbafha, H.; Tabatebaei, M.; Aghbashlo, M.; Khanali, M.; Demirbas, A. A comprehensive review on the environmental impacts of diesel/biodiesel additives. Energy Convers. Manag. 2018, 174, 579–614. [Google Scholar] [CrossRef]
  36. Tian, H.Z.; Lu, L.; Cheng, K.; Hao, J.M.; Zhao, D.; Wang, Y.; Jia, W.X.; Qiu, P.P. Anthropogenic atmospheric nickel emissions and its distribution characteristics in China. Sci. Total Environ. 2012, 417–418, 148–157. [Google Scholar] [CrossRef]
  37. Beddows, D.C.S.; Dall’Osto, M.; Olatunbosum, O.A.; Harrison, R.M. Detection of brake wear aerosols by aerosol time-of-flight mass spectrometry. Atmos. Environ. 2016, 129, 167–175. [Google Scholar] [CrossRef]
  38. Thorpe, A.; Harrison, R.M. Sources and properties of non-exhaust particulate matter from road traffic: A review. Sci. Total Environ. 2008, 400, 270–282. [Google Scholar] [CrossRef]
  39. Amato, F.; Alastuey, A.; Karanasiou, A.; Lucarelli, F.; Nava, S.; Calzolai, G.; Severi, M.; Becagli, M.; Gianelle, V.L.; Colombi, C.; et al. AIRUSE-LIFE+: A harmonized PM speciation and source apportionment in five southern European cities. Atmos. Chem. Phys. 2016, 16, 3289–3309. [Google Scholar] [CrossRef]
  40. Hao, Y.; Meng, X.; Yu, X.; Lei, M.; Li, W.; Shi, F.; Yang, W.; Zhang, S.; Xie, S. Characteristics of trace elements in PM2.5 and PM10 of Chifeng, northeast China: Insights into spatiotemporal variations and sources. Atmos. Res. 2018, 213, 550–561. [Google Scholar] [CrossRef]
  41. Li, W.; Dryfhout-Clark, H.; Hung, H. PM10-bound trace elements in the Great Lakes Basin (1988–2017) indicates effectiveness of regulatory actions, variations in sources and reduction in human health risks. Environ. Int. 2020, 143, 106008. [Google Scholar] [CrossRef]
  42. Ghosh, S.; Rabba, R.; Chowdhury, M.; Padhy, P.K. Source and chemical species characterization of PM10 and human health risk assessment of semi-urban, urban and industrial areas of West Bengal, India. Chemosphere 2018, 207, 626–636. [Google Scholar] [CrossRef] [PubMed]
  43. Kabelitz, T.; Ammon, C.; Funk, R.; Münch, S.; Biniasch, O.; Nübel, U.; Thiel, N.; Rösler, U.; Siller, P.; Amon, B.; et al. Functional relationship of particulate matter (PM) emissions, animal species, and moisture content during manure application. Environ. Int. 2020, 143, 105577. [Google Scholar] [CrossRef] [PubMed]
  44. Wu, S.H.; Cai, M.J.; Xu, C.; Zhang, N.; Zhou, J.B.; Yan, J.P.; Schwab, J.J.; Yuan, C.S. Chemical nature of PM2.5 and PM10 in the coastal urban Xiamen, China: Insights into the impacts of shipping emissions and health risk. Atmos. Environ. 2020, 227, 117383. [Google Scholar] [CrossRef]
  45. Moreno, T.; Querol, X.; Alastuey, A.; de la Rosa, J.; Sánchez de la Campa, A.M.; Minguillón, M.; Pandolfi, M.; González-Castanedo, Y.; Monfort, E.; Gibbons, W. Variations in Vanadium, nickel and lanthnoid element concentrations in urban air. Sci. Total Environ. 2010, 408, 4569–4579. [Google Scholar] [CrossRef] [PubMed]
  46. Xing, W.; Yang, L.; Zhang, H.; Zhang, X.; Wang, Y.; Bai, P.; Zhang, L.; Hayakawa, K.; Nagao, S.; Tang, N. Variations in traffic-related water soluble inorganic ions in PM2.5 in Kanazawa, Japan, after the implementation of a new vehicle emissions regulation. Atmos. Pollut. Res. 2021, 12, 101233. [Google Scholar] [CrossRef]
  47. Paatero, P.; Hopke, P.K. Discarding or downweighting high-noise variables in factor analytic models. Anal. Chim. Acta 2003, 490, 277–289. [Google Scholar] [CrossRef]
  48. Paatero, P.; Hopke, P.K. Rotational tools for factor analytic models. Chemometrics 2009, 23, 91–100. [Google Scholar] [CrossRef]
  49. Alvi, M.U.; Kistler, M.; Mahmud, T.; Shahid, I.; Alam, K.; Chisthie, F.; Hussain, R.; Kasper-Giebl, A. The composition and sources of water soluble ions in PM10 at an urban site in the Indo-Ganfetic Plain. J. Atmos. Sol.-Terr. Phys. 2019, 196, 105142. [Google Scholar] [CrossRef]
  50. Paatero, P.; Hopke, P.K.; Begum, B.A.; Biswas, S.K. A graphical diagnostic method for assessing the rotation in factor analytical models of atmospheric pollution. Atmos. Environ. 2005, 39, 193–201. [Google Scholar] [CrossRef]
  51. Zabalza, J.; Ogulei, D.; Hopke, P.; Lee, J.; Hwang, I.; Querol, X.; Alastuey, A.; Santamaria, J. Concentration and sources of PM10 and its constituents in Alsasua, Spain. Water Air Soil Pollut. 2006, 174, 385–404. [Google Scholar] [CrossRef]
Figure 1. (A) Canary Islands. (B) Geographical location of the area of study (Taliarte) in Gran Canaria island. (C) Taliarte area.
Figure 1. (A) Canary Islands. (B) Geographical location of the area of study (Taliarte) in Gran Canaria island. (C) Taliarte area.
Ijerph 20 05763 g001
Figure 2. Meteorological variables of the study zone from 2010 to 2021. (A) Minimium and maximium temperature (T) in Canary Island. (B) Minimium and maximium relative humidity (RH). (C) Wind direction frequency by octant. Source: Aemet (Spanish Meteorology Agency).
Figure 2. Meteorological variables of the study zone from 2010 to 2021. (A) Minimium and maximium temperature (T) in Canary Island. (B) Minimium and maximium relative humidity (RH). (C) Wind direction frequency by octant. Source: Aemet (Spanish Meteorology Agency).
Ijerph 20 05763 g002
Figure 3. Box–whiskers plots for heavy metals studied in PM 10 for each Saharan dust regime. Red circles indicate the average concentration and numbers show p-values resulting from the Wilcoxon pairwise test.
Figure 3. Box–whiskers plots for heavy metals studied in PM 10 for each Saharan dust regime. Red circles indicate the average concentration and numbers show p-values resulting from the Wilcoxon pairwise test.
Ijerph 20 05763 g003
Figure 4. Box–whiskers plots for heavy metals studied in PM 2.5 for each Saharan dust regimes. Green circles indicate the average concentration and numbers show p-values resulting from the Wilcoxon pairwise test.
Figure 4. Box–whiskers plots for heavy metals studied in PM 2.5 for each Saharan dust regimes. Green circles indicate the average concentration and numbers show p-values resulting from the Wilcoxon pairwise test.
Ijerph 20 05763 g004
Figure 5. Mass percentages of each metal studied over the PM 10 composition. The contribution of the total sum of the metals analyzed is indicated as a percentage.
Figure 5. Mass percentages of each metal studied over the PM 10 composition. The contribution of the total sum of the metals analyzed is indicated as a percentage.
Ijerph 20 05763 g005
Figure 6. Mass percentages of each metal studied over the PM 2.5 composition. The text indicates the contribution of the total sum of the metals analyzed.
Figure 6. Mass percentages of each metal studied over the PM 2.5 composition. The text indicates the contribution of the total sum of the metals analyzed.
Ijerph 20 05763 g006
Figure 7. Chronic health risks at adulthood (Ad) and childhood (Ch) as measured by the values of hazardous quotient (HQ) and the hazard index by inhalation (HI i n h ), oral ingestion (HI i n g ) and dermal contact (HI d e r ). A dashed red line indicates the limit value for non-carcinogenic risk.
Figure 7. Chronic health risks at adulthood (Ad) and childhood (Ch) as measured by the values of hazardous quotient (HQ) and the hazard index by inhalation (HI i n h ), oral ingestion (HI i n g ) and dermal contact (HI d e r ). A dashed red line indicates the limit value for non-carcinogenic risk.
Ijerph 20 05763 g007
Figure 8. Crustal enrichment factors for the metal species studied at the sampling site. The y-axis is shown as a logarithmic scale for clarity. The enrichment scale is represented by different colors and point size. The black diamond indicates the reference metals considered.
Figure 8. Crustal enrichment factors for the metal species studied at the sampling site. The y-axis is shown as a logarithmic scale for clarity. The enrichment scale is represented by different colors and point size. The black diamond indicates the reference metals considered.
Ijerph 20 05763 g008
Table 1. Risk level according to chronic and carcinogenic risks and values for hazardous index (HI) and total carcinogenic risk (TCR).
Table 1. Risk level according to chronic and carcinogenic risks and values for hazardous index (HI) and total carcinogenic risk (TCR).
Type of RiskValueRisk Level
Chronic (non-carcinogenic)HI ≥ 1Causal effects
HI < 1Non causal effects
CarcinogenicTCR ≥ 10 1 Very high
10 3 ≤ TCR <10 1 High
10 4 ≤ TCR < 10 3 Moderate
10 6 ≤ TCR < 10 4 Low
TCR < 10 6 Very low
Table 2. Average concentrations, standard deviation ( μ ± σ ), coefficient of variation (CV) in % and maximum and minimum values (in ng·m 3 ), for each metal studied for PM 10 and PM 2.5 .
Table 2. Average concentrations, standard deviation ( μ ± σ ), coefficient of variation (CV) in % and maximum and minimum values (in ng·m 3 ), for each metal studied for PM 10 and PM 2.5 .
MetalSize μ  ±  σ CVMaxMin
AlPM 10 2430.93 ± 3306.5513615,897.7713.00
PM 2.5 645.00 ± 435.6968875.9959.78
BaPM 10 38.40 ± 20.3553115.2512.00
PM 2.5 26.39 ± 12.094664.6310.86
CdPM 10 0.91 ± 0.71782.740.02
PM 2.5 0.54 ± 0.551021.950.01
CrPM 10 114.41 ± 136.81120856.660.58
PM 2.5 31.59 ± 35.36112189.883.19
CuPM 10 12.01 ± 9.317849.340.34
PM 2.5 5.91 ± 5.799823.020.62
FePM 10 1062.48 ± 1617.881528169.592.73
PM 2.5 304.44 ± 375.771231869.681.58
MnPM 10 27.73 ± 31.79115166.761.05
PM 2.5 5.26 ± 4.548623.930.33
NiPM 10 18.81 ± 10.805763.532.35
PM 2.5 10.96 ± 9.788961.721.87
TiPM 10 89.85 ± 131.19146630.162.66
PM 2.5 17.41 ± 22.95132135.960.76
VPM 10 9.90 ± 8.288438.452.16
PM 2.5 5.81 ± 3.375816.781.76
ZnPM 10 74.81 ± 70.3994391.3210.90
PM 2.5 46.98 ± 30.3865148.417.65
Table 3. TCR i n h values for PM 10 .
Table 3. TCR i n h values for PM 10 .
Heavy MetalChildhoodAdulthood
Cd9.04 × 10 7 3.62 × 10 6
Cr5.33 × 10 3 2.13 × 10 2
Ni9.97 × 10 7 3.99 × 10 6
Table 4. TCR í n h values for PM 2.5 .
Table 4. TCR í n h values for PM 2.5 .
Heavy MetalChildhoodAdulthood
Cd5.83 × 10 7 2.33 × 10 6
Cr1.14 × 10 3 4.54 × 10 3
Ni1.61 × 10 6 6.44 × 10 6
Table 5. Contribution of each factor obtained by the PMF model to each chemical species on PM 10 study (in %).
Table 5. Contribution of each factor obtained by the PMF model to each chemical species on PM 10 study (in %).
SpeciesFactor 1Factor 2Factor 3Factor 4Factor 5Factor 6
PM 10 2523191697
Al72047161
Ba2420020306
Cd001154277
Cr002010070
Cu02119182617
Fe30006901
Mn65112777
Ni637091334
Ti8590024
V3716330112
Zn003219490
Ca035937190
K200237338
Mg 2 + 15268582
Na + 006914180
Cl 24451004
NH 4 + 95706414
NO 3 5423171914
SO 4 2 0141335390
Table 6. Contribution of each factor obtained by the PMF model to each chemical species on PM 2.5 (in %).
Table 6. Contribution of each factor obtained by the PMF model to each chemical species on PM 2.5 (in %).
SpeciesFactor 1Factor 2Factor 3Factor 4Factor 5
PM 2.5 442516114
Al0563132
Ba340322410
Cd0520245
Cr74810035
Cu185301811
Fe0821008
Mn38041813
Ni7023270
Ti12453310
V24369219
Ca9683200
K10033452
Mg1125478
Na + 33010057
NH 4 + 0058510
NO 3 14004143
SO 4 2 103402432
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.

Share and Cite

MDPI and ACS Style

Martín-Cruz, Y.; Gómez-Losada, Á. Risk Assessment and Source Apportionment of Metals on Atmospheric Particulate Matter in a Suburban Background Area of Gran Canaria (Spain). Int. J. Environ. Res. Public Health 2023, 20, 5763. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20105763

AMA Style

Martín-Cruz Y, Gómez-Losada Á. Risk Assessment and Source Apportionment of Metals on Atmospheric Particulate Matter in a Suburban Background Area of Gran Canaria (Spain). International Journal of Environmental Research and Public Health. 2023; 20(10):5763. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20105763

Chicago/Turabian Style

Martín-Cruz, Yumara, and Álvaro Gómez-Losada. 2023. "Risk Assessment and Source Apportionment of Metals on Atmospheric Particulate Matter in a Suburban Background Area of Gran Canaria (Spain)" International Journal of Environmental Research and Public Health 20, no. 10: 5763. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20105763

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

Article Metrics

Back to TopTop