Next Article in Journal
Parameterization of Radiation Fog-Top Height and Methods Evaluation in Tianjin
Next Article in Special Issue
Source Apportionment of PM2.5 in Florence (Italy) by PMF Analysis of Aerosol Composition Records
Previous Article in Journal
One-Year Measurements of Equivalent Black Carbon, Optical Properties, and Sources in the Urumqi River Valley, Tien Shan, China
Previous Article in Special Issue
Application of Miniaturized Sensors to Unmanned Aerial Systems, A New Pathway for the Survey of Polluted Areas: Preliminary Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chemical Composition of PM10 in 16 Urban, Industrial and Background Sites in Italy †

1
Institute of Atmospheric Pollution Research, National Research Council of Italy, Monterotondo st., 00015 Rome, Italy
2
Chemistry Department, Sapienza University of Rome, 00185 Rome, Italy
*
Author to whom correspondence should be addressed.
Presented by Cinzia Perrino, Maria Catrambone and Silvia Canepari on behalf of the IIA-SapCD study group.
Submission received: 16 March 2020 / Revised: 23 April 2020 / Accepted: 5 May 2020 / Published: 8 May 2020
(This article belongs to the Special Issue Recent Advances of Air Pollution Studies in Italy)

Abstract

:
Italy is characterized by a very variable configuration in terms of altitude, proximity to the sea, latitude and the presence of industrial plants. This paper summarizes the chemical characterization of PM10 obtained from 38 sampling campaigns carried out in 16 sites in Italy during the years 2008–2018. Chemical determinations include all macro-components (six macro-elements, eight ions, elemental carbon and organic carbon). The sum of the individual components agrees well with the PM10 mass. The chemical composition of the atmospheric aerosol clearly reflects the variety in the Italian territory and the pronounced seasonal variations in the meteoclimatic conditions that characterize the country. Macro-sources reconstruction allowed us to identify and evaluate the strength of the main PM10 sources in different areas. On 10 sampling sites, the soluble and insoluble fractions of 23 minor and trace elements were also determined. Principal Component Analysis was applied to these data to highlight the relationship between the elemental composition of PM10 and the characteristics of the sampling sites.

Graphical Abstract

1. Introduction

It is well known that air quality is strongly influenced by atmospheric conditions and that atmospheric conditions partly depend on the geographical location and orography of the site. From this point of view, Italy is a very complex environment, characterized by a wide geographic variability. It extends between 47°04′ and 35°29′ degrees north latitude and includes two great mountain ranges, the Alps, which encompass the northern part of the country, and the Apennine Mountains, which extend for about 1300 km along the peninsular region. It also includes a very wide plain to the north, the Po valley (46,000 km2), and about 7500 km of coastline. This varied landscape reflects in a wide variety of climate systems, which range from the cold, wet winters of the northern mountain areas and the strong winter stability of the Po valley to the mild conditions of the coastal regions, temperated by the breeze, and the long, hot summer in the south of the peninsula.
These different features influence air quality to varying extents. In general, air quality shows a wide variability not only according to the intensity of local emission sources but also as a consequence of the mixing properties of the lower atmosphere, the medium- and long-range transport of the air masses, the frequency and intensity of rainfall. Atmospheric particulate matter (PM), characterized by a very complex chemical composition, size distribution and source variety, is mainly influenced by the various geographical and climatic characteristics of the receptor site. It can thus be intriguing to observe similarities and differences in PM chemical composition among different Italian locations.
In the literature, many papers concern the chemical composition of PM in individual cities or geographic areas of Italy [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. However, to our knowledge, no studies take into consideration many locations in different areas of the country. Some papers focus on one city only: Milan [2,13], Genoa [3], Venice [8], Bologna [11], Ferrara [9,10], Lecce [16]. Some others compare a small number of sites in the same region (Lazio [4], Apulia [5], Emilia Romagna [14], Veneto [15]) or along the country [1,6]. Some peculiarities of the Italian region are discussed in the paper of Amato et al. [12], where the chemical composition of PM10 and PM2.5 in Milan and Florence is compared with those measured in Barcelona, Athens and Porto. Finally, the paper of Putaud et al. [7], discusses the chemical composition of PM in more than 60 sites across Europe, including four Italian cities, and highlights regional differences in PM characteristics among northwestern, central and southern Europe.
In this paper, we gather an extensive dataset about the mass concentration and chemical composition of PM10 (particulate matter with aerodynamic diameter less than 10 μm) collected during air quality studies carried out in different Italian sites during the period 2008–2018. The sites range from vast urban areas to small mountain villages. We also included some residential sites close to industrial areas, which are not infrequent in Italy.
We obtained the data during sampling campaigns carried out by using identical instrumental equipment and carried out the analyses by following the same analytical procedure. These features make the dataset very homogeneous and allow us to observe the variability of PM10 composition in different geographical and seasonal conditions.
Chemical characterization includes the determination of macro-components (seven macro-elements, eight ionic species, organic carbon, elemental carbon), which allowed us to obtain the mass closure, and 23 micro-elements, able to trace some sources of particular interest. In most cases, we collected the data during two different periods of the year (warm and cold season), for a total of about 1300 samples. To extract more information from the micro-elements dataset, we applied Principal Component Analysis (PCA); in this receptor model, the mathematical functions that best describe the relative distances between the samples are given by the direction with the maximum variance of the data [17].

2. Experiments

2.1. Sampling Sites and Periods

The sixteen sampling sites, shown on the map in Figure 1, are briefly described below, ordered from north to south.
Chamonix (CH) is a mountain village located in a small valley on the Mont Blanc, in a French territory about 8 km from the Italian border. We choose this area as it is representative of the alpine valleys. Monfalcone (MF) is a small industrial city overlooking the northern Adriatic Sea. Its territory includes a coal cogeneration plant, one of the leading shipyards of the Mediterranean and some engineering industries. Ferrara (FE) is a medium-sized town located in the Po valley. We carried out the samplings at about 5 km from Ferrara, in a site close to the industrial area that includes a waste-to-energy plant and a vast petrochemical complex. Rubbiano (RB) is a small village located in a valley surrounded by hills at the southern border of the Po valley. Its industrial area includes food and expanded clay production plants. La Spezia (LS) is one of the most important Italian military and commercial harbours, in the north-western part of the Mediterranean basin. A power plant is located in the area.
Civitavecchia (CV), on the Tyrrhenian Sea, is a major cruise and ferry port; the city is also the seat of two coal-fired thermal power stations. We carried out the samplings on the coastline at about 5 km from the city centre. Monterotondo scalo (MR) is a village located at about 20 Km NE of Rome, in central Italy. Close to the Tiber River, it is surrounded by agricultural areas. We carried out the samples in the Research Area of the National Research Council of Italy, at about 7 km NE from the village. Guidonia (GU) is a medium-sized town at about 20 km for the city of Rome. Its area is characterized by large travertine caves and cement industry. Rome (RM), the most extensive city in the European Union, about 25 km from the Tyrrhenian Sea, is characterized by the absence of heavy industry and by very high vehicular traffic, with about 2,700,000 private cars. The sampling site was inside the Sapienza University, not directly exposed to vehicular traffic. Fontechiari (FC) is a village in central Italy considered by the local Environmental Protection Agency as a regional background. Frosinone (FR) is an industrial and commercial city located on a hill facing the Sacco river valley and surrounded by mountains. Due to its position, the city experiences extreme atmospheric pollution events. We carried out the samplings in the industrial area, close to the river, at about 5 km from the city centre.
Brindisi (BR), on the eastern coast of Italy, is the leader area in electricity production in Italy. Our sampling site was close to a coal cogeneration plant, at about 10 km from the city centre. Viggiano (VG) is a village located in the Val d’Agri valley, a few km from the most extensive oil platform in continental Europe. Strongoli marina (ST) is a sea village in the south of Italy. A biomass power plant is located nearby (2 km). Palermo (PA), in Sicily, is one of the warmest cities in Italy and Europe, experiencing very long, hot and dry summers. It is one of the major passenger traffic harbours in the Mediterranean; it is also characterized by heavy vehicular traffic (about 750,000 private cars). Gela (GE), is located on the southern coast of Sicily. It is surrounded by crops, but it also borders with a large refinery that was still active during the sampling periods. We carried out the samplings at about 3 km from the city centre.
Table 1 reports the main characteristics of the 16 sampling locations, the number of sampling sites in each area, the years when we carried out the samplings and the total duration of the samplings during the warm (W) and the cold (C) periods. Warm and cold periods are defined as 1 April–30 September and 1 October–31 March, respectively.
The Latitude and longitude of the sampling areas and range of the mean daily temperature during the sampling periods considered in this work are reported in the Supplementary Material (Table S1).

2.2. Sampling and Analytical Procedure

PM10 was collected by using a dual-channel automatic monitor (SWAM 5a Dual Channel Monitor, FAI Instruments, Fonte Nuova, Rome, IT) providing the concentration of PM10 by the beta attenuation method. The instrument has been certified as a beta monitor by both the German TÜV (Technischer Überwachungsverein) and the British MCERTS (Environment Agency of England & Wales Monitoring Certification Scheme); it complies with EU equivalence criteria for PM measurements against EN 12341 and EN 14907 reference methods [18]. The two channels were equipped with Teflon membrane filters (TEFLO, 47 mm, 2.0-micron pore size, PALL Life Sciences, New York, NY, USA.) and quartz fibre filters (TISSUQUARTZ 2500QAT, 47 mm, PALL Life Sciences, New York, NY, USA.). Samples were collected daily, from midnight to midnight, at a flow rate of 2.3 m3/h.
After sampling, Teflon filters were analyzed by energy-dispersion X-ray fluorescence (XRF) (X-Lab 2000 and XEPOS, Spectro Analytical Instruments, Kleve, Germany) for Si, Al, Fe, Na, K, Mg, Ca and minor elements. Then, the filters were extracted for 20 min under sonication in 10 mL of deionized water. A total of 1 mL of the solution was analyzed for ions (chloride, nitrate, sulfate, sodium, potassium, ammonium, magnesium, calcium) by ion chromatography (IC) (ICS1000, Dionex Co., CA, USA). The remaining solution was added with 1 mL acetate buffer (CH3COOH/CH3COOK 0.1 M; pH 4.3). The rationale for the use of this extraction solution is widely discussed in Canepari et al. [19]. Briefly, the use of a pH-buffered solution makes the element solubility percentage independent from the concentration of the acidic species in the sample and directly attributable to the chemical form of elements released by each source. The solution was filtered with nitrocellulose filter (0.45 μm pore size) (Merck Millipore Ltd., MA, USA) and analyzed for the extracted fraction of elements (As, Ba, Bi, Cd, Ce, Co, Cs, Cu, Fe, La, Li, Mg, Mn, Mo, Ni, Pb, Rb, Sb, Sn, Ti, Tl, U, V) by inductively coupled plasma mass spectroscopy (ICP-MS) (Bruker 820 MS); the solid residual on both the sampling filter and the filtration filter was subjected to microwave-assisted acid digestion (HNO3:H2O2 2:1), filtered and analyzed by ICP-MS for the same elements (residual fraction). The results obtained for the two chemical fractions were added to get the total elemental concentration. A full description of the analytical procedure is reported in Canepari et al. [19,20].
Quartz filters were devoted to the analysis of elemental and organic carbon by thermo-optical analysis (OCEC Carbon Aerosol Analyser, Sunset Laboratory, OR, USA), applying the NIOSH–QUARTZ temperature protocol. In the case of site CH, a 1.5 cm2 punch of the filter was extracted in de-ionized water and analyzed for levoglucosan by High-Performance Anion-Exchange Chromatography with Pulsed Amperometric Detection (HPAEC-PAD) [21].
The analytical performances of the overall method were widely checked in its development stage and are reported in previous papers [22,23]. In the case of Mg and Fe, the measurements were carried out by both XRF and ICP. The results of the comparison between the two datasets considered in this work show good correlation for both elements and agree with those previously reported in Canepari et al. [23]: Mg:R2 = 0.93, slope = 1.18; Fe:R2 = 0.82, slope = 0.90.
We carried out explorative PCA using the statistical software CAT (Chemometric Agile Tool [24]) based on the R-project for statistical computing, Ver. 3.0, 32-bit. The data matrix was transformed by performing row and column autoscaling to correct for variations in the different scaling of the examined variables.

3. Results and Discussion

3.1. Macro-Components

For each PM10 sample, we checked the mass closure, that is the correspondence between the concentration of PM10 measured by beta attenuation and the concentration reconstructed by adding the main PM components. To account for non-measured elements (typically H, O), in this calculation, we applied some conversion factors. The rationale for the choice of the conversion factors is widely discussed in the review of Chow et al. [25] and is reported in the Supplementary Material (SM1). Briefly, we considered elements as oxides [26,27]), calculated carbonate as the sum of calcium ion multiplied by 1.5 and magnesium ion multiplied by 2.5, and multiplied organic carbon (OC) by a conversion factor to take into account non-C atoms in organic molecules and obtain organic matter (OM). The conversion factor of OC to OM was set as 2.1 for the agricultural background site FC, 1.6 for urban sites (LS, GU, PA) and 1.8 for all the other sites (residential and peri-urban) [28,29]. In the case of Rome, we applied factor 1.8 because the sampling sites were inside the university area, not directly exposed to traffic emission.
The comparison between the concentrations measured by beta attenuation and reconstructed from the chemical analysis is reported in the Appendix A (Figure A1). The agreement between the two data series is very satisfactory (slope: 0.95; intercept: 0.09; R2 = 0.96; N= 1293) and this ensures that the analytical results were of good quality and all the macro-components were included and confirms that the conversion coefficients for elements and organic carbon were appropriate.
The mean, median, 10th and 90th percentile of the analytical values are reported for all sites in Appendix B, Table A1, for the cold period of the year and Table A2 for the warm period. The data show that, in most cases, the variability of the concentrations, estimated as the ratio between 90th and 10th percentile, is in the range 2–5. Moreover, the difference between the mean and the median value is small and almost centred between the two percentiles, indicating that the concentration distribution is nearly symmetrical. As expected, the concentration variability is higher in the case of natural PM components, such as Al, Si, Na+, Cl. This variability is the consequence of the desert dust and sea-spray transport events that often occur in central-south areas and coastal areas of the country, respectively [5,16].
Considerable variability is also observed for nitrate and ammonium concentration at site MF during the cold season. This difference, which is also reflected in PM10 mass, is due to very different meteorological conditions occurring at this site during the two winter sampling periods. During the first one (12–25 February 2015), meteorological conditions in the area were characterized by moderate air mixing during both daytime and nighttime hours (mean wind velocity during the period: 1.0 m s−1; range of the daily mean: 0.1–2.8 m s−1). During the second one (22 January–4 February 2016), instead, conditions of windless calm and extreme atmospheric stability occurred during the whole period (mean wind velocity during the period: 0.25 m s−1; range of the daily mean; <0.1–1.1 m s−1), being responsible for the build-up of all atmospheric pollutants and favouring the formation and accumulation of secondary PM components. For this reason, in the following discussion, the two sampling periods at MF are considered separately.
For a more straightforward discussion and interpretation of the data, in Figure 2, we report the mean concentration recorded at the 16 sites during the cold period (Figure 2a) and the warm period (Figure 2b), ordered from north to south. For interpreting macro-components’ concentration, we decided to combine components in the case of natural sources (soil and sea) and to keep the rest of the parameters as individual components. Soil contribution was calculated by summing the concentration of elements (as metal oxides) that are generally associated with mineral dust: Al, Si, Fe plus the insoluble fractions of K, Mg and Ca (calculated as the difference between the XRF and the IC determinations), plus calcium and magnesium carbonate. Sea-salt contribution was calculated from the concentration of soluble sodium and chloride, determined by IC, multiplied by 1.176 to take into account minor seawater components (sulphate, magnesium, calcium, potassium) [28,29,30]. It is conceivable that this procedure neglects some minor sources (e.g., brake dust in the case of Fe), but it is suitable to get a comprehensive estimation of the relative weight of the main, ubiquitary PM sources. In the case of organic species, we report OM concentration, obtained by applying the same multiplication coefficients used for getting the reconstructed mass.
Although some variability in the results is certainly due to the different atmospheric conditions during the various sampling periods, which are not simultaneous, the data allow for the evaluation of some general trends. During the cold period, the data show a wide concentration variability among sites (up to 7–8 fold). This depends on the wide variability of meteorological conditions during the cold seasons, which span from windy and rainy days, which cause pollution reduction, to atmospheric stability periods, when the concentration of PM components, mainly secondary species, quickly increases [2,14,31]. This second condition is typical of plains surrounded by hills/mountains that limit air circulation. From a general point of view, the data show a concentration decrease with increasing latitude. Different to the northern areas, the higher temperatures and the milder conditions that occur in the south, due to the proximity to the sea, promote atmospheric mixing and impair the occurrence of long periods of atmospheric stability.
The individual differences in the composition of PM10 among areas depend on the peculiar geographical position of the sampling site and the local sources. Concerning soil, the concentration is particularly high at GU, due to the presence of travertine caves. When moving from north to south, the percent contribution of soil components tends to increase as a result of increasing soil dryness. In the urban area of Rome (RM), instead, the relevance of soil components (about one-quarter of the total PM10 concentration) is the result of traffic-induced soil re-suspension.
The per cent contribution of sea-salt ranges from a few percentage points in areas far from the sea, to about 25% in Sicily (GE). However, the contribution of this source is strongly dependent on the occurrence of transport events during the sampling campaign, and it is thus very variable as a function of the observation period [14].
Concerning inorganic secondary components, the widest variation is observed for nitrate, whose mean concentration ranges from below 1 μg/m3 at VG to more than 20 μg/m3 at MF during the second sampling period. The highest values are recorded in the Po valley and surrounding areas (MF, FE, RB), where a high concentration of precursors (NH3 and NOx) adds to high atmospheric stability and humid climate in favouring the production of ammonium nitrate [14,32]. The contribution of sulphate (in the form of ammonium salts) is generally low due to the low sulphur content of fossil fuels used in the country. Sulphate contribution above 10% are recorded at MF (second period), BR and VG. Here, the role of the industrial plants, whose SO2 emissions can be quickly transformed into sulphate, is probably not negligible (coal cogeneration plants at MF and BR and oil platform at VG). In Sicily, SO2 emission from the Etna volcano and from ships may also give a contribution.
The concentration of organic components shows wide variability, but their contribution to PM10 is high everywhere (on average, about 50%). One of the main winter sources of organics can be identified in biomass burning [12,33,34,35,36]. The highest values (30.2 μg/m3) and highest per cent OM contributions (78%) were recorded at CH, an Alpine village made of small houses heated by fireplaces and stoves. During this campaign, we also determined levoglucosan, a reliable tracer of biomass burning, and found an average concentration as high as 2 μg/m3, (range 0.75–3.2 μg/m3), indicating that this source was responsible, on average, for more than 50% of the organic mass. We recorded very high OM contributions also at FC (70%) and FR (58%), both located in areas where biomass burning is frequently used for domestic heating.
EC concentration and percent contribution are high at RM (9%) and FR (8%), where traffic emission is particularly intense, but also at CH (7%), as a consequence of both the proximity to the Mont Blanc Tunnel and the very high contribution of biomass burning.
During the warm period (Figure 2b), the concentrations in northern Italy are much lower than during the cold season, while they are only slightly lower in the south. Differently from the cold period, during the warm season, a general increase in concentration is observed when moving from north to south. This behaviour is mainly due to the higher temperature and solar radiation occurring at the lower latitudes, which result in the soil being dryer and more intense photochemical activity. During the cold season, besides, atmospheric stability periods generally occurring in northern regions are much less frequent and severe.
In comparison with the cold season, increased soil contribution to PM10 is observed all over the country, with an average value of about one third. In this case, as well, the highest concentration and contribution are recorded at GU. The contribution of sea-salts is very variable, as discussed for the cold period, but it is generally higher at the three harbours (LS, CV, PA). Nitrate is the component that experiences the widest concentration reduction from the cold to the warm semester, as the equilibrium between particulate ammonium nitrate and ammonia and nitric acid in the gaseous phase is sensitive to the temperature. Moreover, the reduction in NOx emissions from traffic and domestic heating may play a role [16]. Conversely, ammonium sulphate concentration increases during the warm period, due to photochemical conversion of SO2, particularly in areas (LS, PA, GE) where SO2 concentration is influenced by local emission sources (fossil fuel combustion from ship and power plants and, in Sicily, the Etna volcano) [5,16,37]. Organic matter also shows a significant concentration decrease during the warm period, despite the increase in photochemical bioVOCs oxidation. This confirms, once again, the relevance of the contribution of biomass burning for domestic heating, which switches off during the warm months. The decrease in nitrate and OM concentration is the main factor responsible for the general reduction in PM10 concentration during the warm period. EC concentration values are very similar to those recorded during the cold period, and percent contribution is slightly higher due to the reduction in the concentration of the other components.
Our findings are in good agreement with those of the AIRUSE project [12]: the average composition in Milan (OM 38%, EC 5%, sulphate 7%, nitrate 14%, ammonium 8%, sea-salt 2%, soil dust 9%) are in the range of our results for northern Italy. In the same paper, the chemical composition in Milan (in the Po valley) and in Florence (located in a closed basin in Central Italy) showed higher PM10 levels (by a factor of more than 2) in autumn and winter with respect to the warm season and the cities of Athens, Barcelona and Porto. These high levels were attributed to the stagnant meteorological conditions that prevail in Milan and Florence during the cold period. Very high PM10 values in comparison with the 60 European examined locations were also recorded in the Po valley (Bologna) by Putaud et al. [7].

3.2. Micro- and Trace- Elements

The analysis of micro- and trace- elements was carried out on a subset of sites: MF, FE, RB, RM, BR, VG, ST, GE (cold period) MF, FE, RB, LS, CV, RM, GE (warm period). The mean, median, 10th percentile and 90th percentile of the concentrations are reported, for all sites, in Appendix C, Table A3, for the cold period of the year and Table A4 for the warm period.
Also in this dataset, the mean and median concentration values are generally very close to each other and in the middle range to the 10th and 90th percentiles, suggesting a symmetrical distribution of the data. As observed for the macro-components, the two winter periods at MF show very different results, with elemental concentration 2–4 times higher during the second campaign. In general, a higher variability range of the elemental concentrations is observed with respect to macro-components (ratio between 90th and 10th percentile up to about 10). This variability is probably due to the higher sensitivity of micro-components to local or discontinuous minor PM sources, which may show a wide strength variability in time (sources that have a small impact on PM mass concentration but may significantly affect the concentration of micro-tracers).
In the case of sites MF, FE, RB, RM and GE, it is possible to discuss the seasonal variations. In agreement with the behaviour observed for macro-components, the data show that the occurrence of atmospheric stability conditions has a definite effect on the measured concentration, mainly in sites that experience severe episodes (Po valley: MF and FE). The ratios between the cold and the warm periods widely exceed unity at MF (second campaign), FE and RM, mainly for elements related to industrial activities (As, Bi, Cd, Ni, Pb, Sn), traffic emission (Ba, Cu, Fe, Mo, Sb) and biomass burning (Cs, Rb, Tl) [10,36,38,39]. The difference along the year is instead negligible in the other sites (RB and GE).
For a more straightforward interpretation of the results, we run the PCA on the elemental concentration obtained by ICP-MS; in this case, we considered the two MF winter sampling periods separately. In this analysis, we chose not to perform PCA on individual campaigns but to include all the available data. This choice does not allow a reliable source apportionment in individual sites because each studied area is characterized by different, peculiar sources. However, it permits the easy, visual identification of differences and similarities in the elemental composition of PM among the considered sites.
PCA analysis extracted six components, able to explain 81.1% of the total variance; the loadings are reported as Supplementary Material (Table S2).
The first component, which explains 44.5% of the variance, mainly reflects the overall variation in PM concentration; the score and loading plot of Components 1 and 2 are reported in the Supplementary Material (Figure S1).
Figure 3 shows the score plot (a) and loading plot (b) of Components 2 and 3, which explain 15.8% and 8.1% of the variance, respectively. The elemental composition clearly differentiates the sites: the samples taken at each site group in a different area of the score plot, independently from the period of the year when they were collected. In particular, the samples collected at MF, FE, CV and RM are well separated from each other. MF data are found in the upper-right quarter of the score plot; in the loading plot, this area corresponds to elements of industrial origin (Bi, As, Mo, Cd, Sn, Pb). The same elements, probably with different relative weight, also characterize FE results (upper side of the score plot). The RM samples (upper-left side of the score plot) correspond, as expected, to the area of the loading plot where traffic elements group (Cu, Sb). In the literature, Cu/Sb is widely considered as a diagnostic ratio of traffic sources, with values around 4–5 [40]. Our data show that, in sites with a high traffic intensity (LS, RM, but also MF and FE), this ratio is much higher than 4 (range 7–12). These higher Cu/Sb values are probably due to the phasing-out of Sb from brake pads, due to the growing concern about Sb adverse effects on health [41]. Finally, CV samples spread over the left side of the score plot and are characterized by high concentrations of La, Ce, Ba, U, Mg and V.
Many literature studies have shown that the chemical fractionation of elements is a robust tool to increase their performance as source tracers [10,42]. In particular, combustion/industrial sources were found to release elements as chemical species that are generally more soluble than those released by abrasive processes, such as non-exhaust traffic emission [11,42,43,44,45,46]. The descriptive statistics of the concentration of the elemental soluble fraction during the two seasonal periods is reported as Supplementary Material (Tables S3 and S4). In Figure 4, we report the results of the PCA, run on the soluble fraction of elements, considering all sites and periods. Considering that this fraction is mainly associated with combustion/industrial sources, it is reasonable to think that that it is more effective in separating these contributions.
In this case, as well, we considered six components, explaining 73.1% of the total variance; the loadings are reported as Supplementary Material (Table S5). Again, the first component, which explains 36.0% of the variance, reflects only the variation in element concentration; the score and loading plot of Components 1 and 2 are reported in Supplementary Material (Figure S2). Figure 4 shows that MF results group in the upper-left part of the score plot and are well separated from FE data. In the loading plot, they correspond to Sn, Fe, Bi and Li, elements that typically are of industrial origin. It is worth noting that, despite the very different concentration levels recorded at MF during the two winter sampling periods, all the MF data align on the same trajectory. This finding further confirms that the stability of the atmosphere can increase/decrease pollutants’ concentration, but, in the case of primary pollutants, it does not alter the elemental profile. FE samples, instead, are grouped in the upper-right quarter of the score plot and correspond, in the loading plot, to a different group of industrial elements: Cd, As and Pb. These elements are recognised as reliable tracers of waste-to-energy plants [47]. It is interesting to note, however, that in a previous paper we showed that in the area of FE the concentration of these elements does not seem to be directly related to local sources but is homogeneous over a vast territory [10].
On the lower side of the score plot, we find LS, CV, GE, all sited on the coastline and close to a coal power plant or an oil refinery. In the loading plot, this area corresponds to Ni and V, related to heavy oil combustion [48,49]. Our data (Appendix C) show that, in agreement with the literature, the Ni/V ratio is lower in coastal sites. At LS, CV, BR, GE, the value of Ni/V it is well below one, while, in all the other cases, it is higher than one and increases during the cold period. It has been suggested that this increase may be influenced by biomass burning [50,51].
Further information about the variability of the source strength can be obtained by considering the solubility percentage of the elements. As reported in previous papers [10,42], the solubility percentage of each element is strictly dependent on its chemical form, which, on its turn, is typical of the source that has emitted the aerosol in the atmosphere. Variation in the solubility of elements therefore indicates variation in the relative strength of the aerosol sources. The descriptive statistics of the percent solubility of each element during the two seasonal periods are reported in the Supplementary Material, Tables S6 and S7. The data show that the solubility percentage is largely associated with the considered element. The soluble fraction typically constitutes more than 75% in the case of Cd and Tl and less than 25% for Bi, Ce, Fe, La, Sn, Ti, U, while the remaining elements are more distributed between the two fractions. However, we observe a certain variability in the solubility that is due to the different strength of individual sources in the various sites and periods. To evaluate these issues in our dataset, we run the PCA on the solubility percentage (six components, explaining 74.9% of the total variance; loadings in Supplementary Material, Table S8).
Our results confirm that the chemical fractionation increases the ability of elements to trace PM sources, as the PCA run on the solubility percentage provides better separation among sites than the PCA run on concentrations. As an example, we report, in Appendix D (Figure A2), the 3D score plot of Components 3, 4 and 5, where all the sites appear well separated from each other.
Figure 5 shows the score plot and loading plot f the solubility percentage for Components 2 and 3 (13.8% and 8.1%, respectively, of the total variance). The score and loading plot of Components 1 and 2 are reported in the Supplementary Material (Figure S3).
For sites MF, FE, RB, RM (inland), the cold period is drawn in green and the warm period in red. For sites LS, CV, VG, ST, GE, (coastal/petrochemical), instead, the data are drawn in black, independently of the season. The three groups take up different areas of the graph. Samples collected during the cold semester in inland areas are grouped in the lower-right quarter of the plot and correspond, in the loading plot, to high values of the soluble fraction of Rb, Cs and Li, which is a robust tracer of biomass burning [36,39]. Samples collected during the warm semester in inland sites (upper-right quarter) correspond to a high solubility of As, Sn, Sb e Bi. As discussed above, the soluble fraction of these elements traces industrial sources. Coastal/petrochemical sites, instead, do not show any significant seasonal difference in solubility percentage, as they are mostly influenced by the solubility of elements released by sea-spray (Mg) or by vessel traffic (Ni and V). In facts, Mg in sea-spray is in the form of the soluble MgCl2, and Ni and V in particles emitted by heavy oil-fuelled sources are generally more soluble than elements released by other sources.

4. Conclusions

We investigated the variability of PM10 composition in 16 sampling areas in Italy using both mass closure and a PCA receptor model. Variability between cold and warm periods and among sites was observed, reflecting the geographical position of the site and the influence of the local sources.
A yearly PM10 concentration gradient is recorded all over the country, with a decreasing concentration from north to south during the cold season and an increasing concentration during the warm season. Soil components increase in the same direction. They also increase during the warm semester, as a consequence of soil dryness. Very high ammonium nitrate concentrations are recorded in the Po valley during the cold period, while ammonium sulphate increases during the warm period, mainly in southern and coastal sites. Traffic contribution is more evident in urban areas, particularly in Rome, in the form of both direct exhaust emission and re-suspension. The contribution of organics is very high everywhere, but a significant concentration decrease is observed during the warm season, indicating a relevant contribution of biomass burning during the winter.
As expected, the concentration of micro- and trace elements is mainly driven by the strength of local sources. The elemental profiles are enriched in elements of industrial origin (Bi, As, Mo, Cd, Sn, Pb) at MF (power plant, shipyards, engineering industries) and FE (waste-to-energy plant, petrochemical complex), in elements from traffic exhausts and re-suspension (Ba, Cu, Fe, Mo, Sb) at the urban sites (LS, RM), in Ni, V and Mg at the coastal sites (LS, CV, BR, GE). The study of the soluble fraction of elements, which trace combustion sources more efficiently, allowed for better differentiation between industrial sites (MF and FE). Solubility percentage has been shown to be a robust tool for the identification of PM sources by elemental analysis. We observed a significant increase in the elemental solubility of some tracers of biomass burning (Rb, Li, Cs) during the cold season, as the switch-on of this seasonal PM source causes a substantial increase in their solubility.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4433/11/5/479/s1, SM1: Rationale for the choice of conversion factors; Figure S1: PCA of elemental concentration (all data): score plot (left panel) and loading plot (right panel) of components 1 and 2. _C and _W after the site code indicate the cold and warm periods, respectively. In the case of MF, the two winter periods are reported as _C1 and _C2; Figure S2: PCA of the soluble fraction of elements (all data): score plot (left panel) and loading plot (right panel) of Components 1 and 2; Figure S3: PCA of the solubility percentage of elements (all data): the score plot (left panel).and loading plot (right panel) of Components 1 and 2; Table S1: Latitude, longitude of the sampling sites and mean temperature during the sampling periods; Table S2: Loadings of PCA run on the total elemental concentrations; Table S3: Concentration of the micro- and trace elements (ng/m3) in the soluble fraction during the cold period; Table S4: Concentration of the micro- and trace elements (ng/m3) in the soluble fraction during the warm period; Table S5: Loadings of PCA run on the soluble fraction of elements, Table S6: Percent solubility of micro- and trace elements during the cold period; Table S6: Percent solubility of micro- and trace elements during the cold period; Table S7: Percent solubility of micro- and trace elements during the warm period; Table S8: Loadings of PCA run on the solubility percentage of elements; Figure S3: PCA of the soluble fraction of elements (all data): score plot (left panel) and loading plot (right panel) of components 1 and 2.

Author Contributions

Project supervision, C.P.; conceptualizzation, C.P. and S.C.; writing, C.P. and S.C.; supervision (IIA study group), C.P.; supervision (SapCD study group), S.C.; sampling campaigns supervision, M.C.; analytical supervision, M.C.; sampling activity, IIA study group; analytical activity and data elaboration, IIA-SapCD study group. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The Authors are indebted to T. Davanzo, M. Palozzo and E, Zappaterreno for their administrative support. IIA-SapCD study group: S. Dalla Torre, M. Giusto, F. Marcovecchio, S. Pareti, E. Rantica, T. Sargolini, L. Tofful (IIA); M.L. Astolfi, D. Frasca, M.A. Frezzini, M. Marcoccia, L. Massimi, G. Simonetti (SapCD).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Scatter plot of PM10 concentration as measured by β attenuation and as reconstructed by the sum of the chemical analyses.
Figure A1. Scatter plot of PM10 concentration as measured by β attenuation and as reconstructed by the sum of the chemical analyses.
Atmosphere 11 00479 g0a1

Appendix B

Table A1. Concentration (mean, median, 10th and 90th percentile) of the main PM10 components, in μg/m3, during the cold period.
Table A1. Concentration (mean, median, 10th and 90th percentile) of the main PM10 components, in μg/m3, during the cold period.
Site
Samples
StatisticsAlSiFeKMgCaClNO3SO42−Na+NH4+K+Mg2+Ca2+OCEC
CH
N = 14
MEAN0.050.270.340.27< 0.020.350.481.30.940.310.450.240.0140.2516.82.8
MEDIAN0.040.180.320.27< 0.020.330.331.30.910.210.440.260.0140.2615.62.8
10th0.020.050.060.13< 0.020.100.210.720.670.130.250.110.0040.0709.61.8
90th0.120.650.580.41< 0.020.751.11.81.20.680.620.360.0260.5125.03.9
MF
N = 107
MEAN0.120.380.390.610.080.950.22133.40.244.10.490.0600.859.40.96
MEDIAN0.120.340.260.480.070.670.134.03.40.231.90.380.0500.647.60.90
10th0.070.180.090.240.040.230.0320.480.930.0860.360.190.0190.283.60.31
90th0.180.630.941.170.141.880.51376.60.4211.10.990.111.716.91.6
FE
N = 82
MEAN0.140.470.250.600.120.750.45112.30.282.90.240.0560.669.21.0
MEDIAN0.130.390.200.530.100.560.338.42.10.232.50.260.0530.538.40.91
10th0.070.200.080.320.050.220.0692.30.830.0930.720.0200.0130.204.70.61
90th0.210.830.491.040.221.400.85234.00.525.90.530.111.415.61.7
RB
N = 21
MEAN0.321.300.430.380.230.870.305.00.990.301.40.310.0470.334.50.58
MEDIAN0.190.800.300.340.160.770.234.40.800.260.900.270.0390.284.10.60
10th0.100.300.090.180.070.290.0630.620.270.0900.180.150.0150.172.50.33
90th0.632.630.880.590.411.550.72112.20.643.10.470.0900.567.50.73
MR
N = 55
MEAN0.160.530.261.130.160.830.522.91.80.600.460.870.0860.578.40.77
MEDIAN0.100.330.231.000.110.700.192.91.60.420.350.830.0640.498.40.81
10th0.060.150.140.310.040.170.0560.830.970.0930.130.150.0180.183.10.22
90th0.240.860.371.970.311.461.34.62.71.30.861.60.200.8915.11.2
GU
N = 28
MEAN0.251.030.560.580.134.440.563.53.00.660.320.520.113.48.61.7
MEDIAN0.231.020.580.570.134.130.123.32.20.380.300.460.0892.97.21.5
10th0.090.380.150.240.061.510.0731.41.00.170.130.210.0411.13.30.53
90th0.431.801.030.990.207.501.45.86.51.30.510.930.186.716.32.6
RM
N = 42
MEAN0.150.580.580.460.181.141.01.91.20.820.380.270.120.596.92.6
MEDIAN0.110.430.510.390.120.910.231.51.00.390.250.230.0520.456.71.9
10th0.070.230.140.210.060.450.0430.580.330.0910.0920.0780.0200.233.00.77
90th0.261.040.890.740.381.923.43.91.82.40.720.500.340.9910.85.6
FC
N = 61
MEAN0.140.500.070.440.200.850.271.31.60.370.450.370.0720.579.50.69
MEDIAN0.040.100.010.430.090.540.171.11.60.220.420.340.0610.399.60.69
10th0.020.030.010.200.030.240.0650.370.580.100.170.180.0220.144.80.26
90th0.431.630.270.690.441.500.602.52.60.750.810.600.141.213.81.1
FR
N = 120
MEAN0.220.670.260.810.081.430.633.72.20.651.20.740.111.013.73.4
MEDIAN0.160.470.230.750.051.170.483.20.890.430.920.640.0830.8311.73.1
10th0.060.180.070.380.020.410.131.50.460.130.380.340.0210.366.41.3
90th0.361.190.531.360.192.891.27.56.81.52.31.30.242.023.36.0
BR
N = 42
MEAN0.240.850.170.460.150.820.682.53.40.741.0n.d.0.120.636.20.81
MEDIAN0.210.650.170.430.140.770.272.43.40.521.0n.d.0.110.635.50.77
10th0.080.270.060.300.080.440.150.921.70.190.24n.d.0.0510.312.20.36
90th0.471.750.310.720.251.111.84.75.11.91.7n.d.0.260.8510.61.2
VG
N = 28
MEAN0.080.260.040.150.090.240.290.550.990.330.330.0920.0500.152.40.20
MEDIAN0.060.170.020.150.080.190.0870.500.990.220.300.0890.0450.122.50.20
10th0.030.080.010.090.050.090.0310.300.380.0790.100.0450.0240.0691.30.10
90th0.130.480.070.210.140.450.690.891.80.630.660.140.0790.263.50.32
ST
N = 27
MEAN0.230.940.190.500.200.791.21.22.01.00.660.290.180.595.10.96
MEDIAN0.180.720.130.470.160.600.521.01.90.640.640.260.130.485.20.95
10th0.110.370.060.310.070.380.0720.590.830.170.220.140.0290.332.60.43
90th0.442.120.450.740.321.622.12.13.11.71.10.460.351.17.31.6
GE
N = 28
MEAN0.160.760.200.300.282.023.02.02.21.90.31n.d.0.271.32.70.75
MEDIAN0.100.510.150.210.271.752.12.11.91.60.24n.d.0.231.12.60.75
10th0.060.270.070.080.080.710.341.21.10.700.11n.d.0.110.912.20.56
90th0.341.510.340.400.463.647.52.53.53.90.64n.d.0.521.83.60.96
Table A2. Concentration (mean, median, 10th and 90th percentile) of the main PM10 components, in μg/m3, during the warm period.
Table A2. Concentration (mean, median, 10th and 90th percentile) of the main PM10 components, in μg/m3, during the warm period.
Site
Samples
StatisticsAlSiFeKMgCaClNO3SO42−Na+NH4+K+Mg2+Ca2+OCEC
CH
N = 14
MEAN0.170.770.160.180.070.450.0390.371.30.170.340.0910.0390.453.81.2
MEDIAN0.160.770.170.180.070.460.0340.321.40.170.330.0780.0350.464.01.2
0.10.040.190.030.090.030.140.0230.120.770.0660.170.0620.0200.141.90.87
0.90.281.210.260.260.090.680.0550.661.90.260.530.140.0620.685.11.4
MF
N = 111
MEAN0.130.380.190.170.080.820.0490.232.00.140.700.140.0250.264.10.35
MEDIAN0.120.360.180.170.080.810.0280.211.70.130.590.140.0220.243.90.36
0.10.110.260.110.140.060.530.0140.110.870.0830.300.120.0110.112.80.26
0.90.160.520.270.200.111.120.110.353.40.231.20.160.0430.425.60.45
FE
N = 84
MEAN0.160.620.240.190.120.980.0701.22.20.190.93n.d.0.0440.423.70.45
MEDIAN0.140.560.250.170.110.860.0450.922.20.150.88n.d.0.0440.373.70.42
0.10.090.270.050.100.050.410.0280.140.640.0850.38n.d.0.0060.0712.50.28
0.90.251.120.420.290.191.770.122.34.00.371.7n.d.0.0850.855.00.71
RB
N = 21
MEAN0.301.310.400.320.210.820.171.02.00.320.540.190.0520.264.10.55
MEDIAN0.231.040.330.280.160.770.0820.782.00.240.450.170.0430.233.80.49
0.10.160.720.230.180.110.580.0280.270.840.120.200.110.0180.132.90.38
0.90.462.110.660.500.311.200.531.83.30.621.00.300.100.415.70.81
LS
N = 14
MEAN0.271.080.350.250.421.050.991.53.51.40.740.120.260.955.51.7
MEDIAN0.291.170.390.260.321.060.201.43.50.880.680.110.200.975.71.7
0.10.190.750.180.210.230.670.110.821.60.460.170.100.130.643.90.96
0.90.351.450.490.290.711.323.32.15.03.31.30.170.521.26.42.3
CV
N = 21
MEAN0.281.100.300.300.301.841.92.21.81.60.350.120.281.23.11.2
MEDIAN0.200.800.210.240.301.491.22.01.81.40.320.130.241.12.80.73
0.10.100.330.070.140.160.730.400.970.570.500.150.0480.120.591.90.43
0.90.672.690.620.570.463.684.63.82.73.00.570.180.492.04.81.5
MR
N = 56
MEAN0.361.300.390.460.201.90.120.863.40.360.750.170.0710.994.20.50
MEDIAN0.301.020.330.460.161.80.0420.683.30.180.700.160.0550.824.10.49
0.10.160.460.190.290.080.570.0230.321.30.0630.390.100.0290.313.20.26
0.90.612.330.630.640.263.80.171.75.50.881.20.260.131.95.40.72
GU
N = 14
MEAN0.311.210.520.450.243.650.892.02.80.930.240.320.192.94.92.0
MEDIAN0.311.190.450.440.173.460.131.83.10.520.250.330.132.85.02.0
0.10.090.340.190.200.121.480.0571.31.60.250.0690.130.101.02.41.1
0.90.551.991.000.720.465.833.02.73.52.10.450.510.414.97.82.9
RM
N = 14
MEAN0.270.670.390.160.160.820.110.651.70.290.550.0820.0530.504.11.0
MEDIAN0.250.620.350.160.130.770.0520.561.60.300.570.0860.0480.524.00.97
0.10.150.470.280.130.120.630.0370.300.710.0920.250.0470.0260.303.30.83
0.90.380.830.520.210.221.030.311.22.50.490.860.120.0820.685.01.3
FC
N = 102
MEAN0.411.500.240.500.341.450.241.43.60.630.900.370.110.915.40.61
MEDIAN0.210.780.110.370.201.050.0591.13.30.470.870.330.0820.654.70.52
0.10.070.230.010.160.100.590.0340.441.70.130.460.140.0400.332.80.30
0.91.234.540.731.060.922.930.562.55.61.31.40.620.201.88.50.95
FR
N = 41
MEAN0.160.570.340.400.191.780.273.33.60.580.920.400.101.17.43.3
MEDIAN0.160.550.260.370.131.590.283.03.20.330.900.380.0880.816.43.4
0.10.100.290.110.230.061.050.141.91.50.190.470.190.0370.483.82.0
0.90.240.900.720.640.372.610.435.07.00.971.40.650.152.113.25.2
VG
N = 28
MEAN0.200.660.130.210.130.570.0370.372.30.190.690.130.0520.433.30.34
MEDIAN0.180.530.120.200.110.530.0150.282.10.130.670.140.0410.402.90.32
0.10.130.420.090.150.090.420.0100.221.50.0790.420.0820.0310.312.30.23
0.90.290.990.170.290.170.760.0810.723.20.451.00.180.0880.654.40.48
ST
N = 28
MEAN0.180.860.180.350.090.770.550.613.10.570.780.240.0860.374.40.88
MEDIAN0.180.850.190.360.090.740.330.563.50.400.720.250.0760.304.80.88
0.10.070.440.070.210.070.430.210.470.590.270.100.140.0410.192.30.39
0.90.331.490.290.490.111.050.940.864.81.01.60.320.140.605.81.3
PA
N = 56
MEAN0.170.680.250.230.572.230.661.94.31.11.40.170.291.44.51.2
MEDIAN0.150.570.240.230.492.070.251.73.80.801.30.170.231.34.41.0
0.10.100.320.070.150.311.190.0710.781.90.430.460.100.120.902.60.80
0.90.291.170.470.300.973.841.62.96.31.72.60.220.512.16.11.7
GE
N = 27
MEAN0.180.840.300.570.162.050.291.27.70.981.6n.d.0.111.63.70.80
MEDIAN0.170.800.300.410.151.980.161.28.00.911.6n.d.0.101.73.70.76
0.10.150.700.260.320.091.570.0350.585.60.481.0n.d.0.0650.862.90.51
0.90.231.070.340.600.232.610.792.09.71.82.5n.d.0.172.54.51.2

Appendix C

Table A3. Concentration (mean, median, 10th and 90th percentile) of the micro- and trace- elements, in ng/m3, during the cold period.
Table A3. Concentration (mean, median, 10th and 90th percentile) of the micro- and trace- elements, in ng/m3, during the cold period.
Site
Samples
StatisticsAsBa Bi Cd Ce Co Cs Cu Fe La LiMgMn Mo Ni Pb Rb Sb Sn TiTl UV
MF1
N = 53
MEAN0.603.80.130.160.120.080.0287.63410.060.126215.90.581.94.51.20.91.63.30.0230.0061.6
MEDIAN0.543.00.090.150.100.070.0277.02760.050.09478.80.501.64.30.90.81.22.30.0210.0051.2
0.10.181.80.050.060.050.030.0123.01250.020.04343.80.190.81.60.60.40.51.00.0080.0030.5
0.91.197.20.250.280.220.150.04612.36890.110.2011929.91.033.28.42.21.73.27.80.0440.0113.5
MF2
N = 54
MEAN1.1811.30.450.430.220.240.05321.86140.140.378948.12.585.611.62.42.87.57.90.0480.0164.3
MEDIAN1.159.70.360.400.180.220.05118.85070.120.308434.72.264.711.22.32.26.75.70.0470.0144.0
0.10.384.30.160.130.080.100.0248.02090.040.14389.10.881.94.80.91.01.62.00.0210.0060.8
0.92.0419.30.780.730.360.430.08441.611020.290.65148111.85.0410.120.84.05.115.019.20.0760.0268.4
FE
N = 82
MEAN1.496.80.340.410.170.200.05416.13950.210.2311012.82.032.613.00.92.25.42.90.0420.0091.3
MEDIAN1.355.80.290.300.140.150.04013.33310.150.189711.31.752.110.90.81.84.62.40.0340.0081.2
0.10.712.20.150.140.070.080.0246.41490.060.14355.00.901.14.90.31.01.81.10.0230.0050.6
0.92.2413.10.540.750.290.350.08731.87210.480.3720921.03.484.225.71.94.010.75.10.0650.0162.1
RB
N = 21
MEAN0.484.90.130.100.260.800.0227.63360.140.321539.00.875.03.50.60.82.09.00.0160.0091.0
MEDIAN0.415.10.120.080.150.600.0177.92990.130.191048.40.644.53.00.60.61.94.10.0140.0070.6
0.10.201.90.070.040.060.230.0104.11390.030.08434.20.381.51.30.30.30.91.60.0070.0030.4
0.90.817.80.180.220.661.830.03910.36160.300.6726713.31.569.06.30.91.33.420.50.0240.0162.3
RM
N = 42
MEAN0.5115.10.250.220.670.150.08530.96190.350.171688.51.812.66.22.03.74.35.40.1400.0151.2
MEDIAN0.5314.00.230.220.660.130.07226.35300.310.141107.11.402.45.61.62.73.64.20.1110.0121.0
0.10.2511.00.090.080.200.060.0319.12490.110.08643.50.601.61.90.71.01.62.30.0330.0080.3
0.90.8420.20.440.351.090.240.15352.69690.580.2636614.93.044.39.93.68.47.79.50.2420.0252.1
BR
N = 42
MEAN0.322.60.130.120.220.080.0262.9148n.a.0.151243.3n.a.1.93.90.70.90.74.10.0360.0103.9
MEDIAN0.322.50.060.110.200.070.0262.8145n.a0.141123.4n.a1.93.50.70.70.64.20.0280.0103.5
0.10.201.60.020.060.100.050.0181.577n.a.0.07541.6 n.a.1.02.10.40.30.31.60.0150.0052.1
0.90.463.60.350.190.370.120.0344.3221n.a.0.262354.7n.a.2.66.00.91.61.77.30.0790.0165.6
VG
N = 28
MEAN0.194.50.030.04n.a.0.020.0151.148n.a.0.05611.60.050.60.80.20.30.20.04n.a.n.a.0.7
MEDIAN0.144.50.020.03n.a0.020.0101.039n.a0.04611.50.030.60.90.20.20.20.03n.an.a0.6
0.10.101.20.010.01n.a.0.010.0050.619n.a.0.02370.60.010.30.20.10.10.10.01n.a.n.a.0.2
0.90.247.90.050.06n.a.0.040.0251.885n.a.0.08843.20.091.01.30.30.50.30.10n.a.n.a.1.2
ST
N = 41
MEAN0.381.9n.a.0.090.120.010.0201.8670.080.161766.20.172.11.02.20.50.40.40.1560.0051.5
MEDIAN0.361.8n.a0.090.100.010.0161.8630.050.151656.20.171.50.82.20.50.50.30.1470.0050.9
0.10.220.9n.a.0.050.060.010.0131.1340.040.081131.80.110.80.41.90.20.20.20.0670.0020.5
0.90.583.2n.a.0.150.200.020.0332.91000.130.252409.10.224.11.82.50.80.70.70.2770.0092.9
GE
N = 28
MEAN0.247.30.050.070.280.140.0233.82050.200.202305.30.692.13.10.40.80.54.70.0150.0154.4
MEDIAN0.276.70.040.070.200.130.0203.71570.190.142244.80.391.92.30.30.80.53.30.0130.0134.7
0.10.153.70.020.050.120.060.0122.51080.080.101312.90.291.51.90.20.50.31.80.0080.0072.6
0.90.3611.50.060.080.540.260.0304.63160.320.343498.21.542.94.90.61.40.69.00.0230.0236.0
Table A4. Concentration (mean, median, 10th and 90th percentile) of the micro- and trace elements, in ng/m3, during the warm period.
Table A4. Concentration (mean, median, 10th and 90th percentile) of the micro- and trace elements, in ng/m3, during the warm period.
Site
Samples
StatisticsAsBa Bi Cd Ce Co Cs Cu Fe La LiMgMn Mo Ni Pb Rb Sb Sn TiTl UV
MF
N = 111
MEAN0.495.40.130.100.220.090.0307.13130.110.197615.50.622.03.20.470.71.25.30.0220.0101.8
MEDIAN0.474.10.110.080.130.080.0206.62900.080.125811.60.551.73.00.430.61.23.40.0210.0081.3
0.10.241.70.050.030.060.040.0073.91370.030.04264.00.311.01.30.130.40.71.10.0060.0030.4
0.90.8011.30.280.180.520.180.06810.85270.250.4314434.81.033.15.50.890.91.912.50.0400.0203.6
FE
N = 84
MEAN0.624.20.130.130.210.110.0308.52470.120.181447.40.492.03.50.411.01.33.10.0290.0092.2
MEDIAN0.573.90.120.100.210.100.0258.12270.120.171267.20.491.73.10.380.81.22.50.0240.0081.8
0.10.272.20.080.050.100.050.0124.91300.050.09754.30.280.91.70.210.40.61.30.0090.0040.7
0.90.976.70.200.250.400.190.04912.34030.200.3022210.80.763.25.80.661.62.15.40.0500.0164.5
RB
N = 21
MEAN0.576.10.140.070.270.470.0218.93830.160.2816010.60.753.02.50.430.81.96.90.1460.0122.2
MEDIAN0.425.90.120.060.260.390.0208.03540.150.231369.40.752.82.30.420.71.75.00.1160.0112.0
0.10.324.60.090.030.160.170.0126.02650.080.16996.00.382.11.50.310.41.13.20.0580.0060.9
0.90.938.50.210.140.400.790.03112.75970.210.4823916.11.164.13.30.561.22.912.10.2910.0183.3
LS
N = 14
MEAN0.439.20.350.110.320.250.04816.25220.210.3833410.30.926.13.00.901.41.510.60.0120.0149.3
MEDIAN0.439.30.170.100.350.260.05015.45610.230.3730811.91.006.43.30.941.41.710.40.0130.01410.3
0.10.283.90.070.040.200.140.0245.92350.130.182225.10.433.71.90.610.40.74.80.0030.0085.3
0.90.5914.40.730.150.430.350.06823.57420.260.5151313.41.328.04.01.212.02.215.30.0200.01712.1
CV
N = 21
MEAN0.6725.10.060.060.920.400.0755.53080.590.292708.50.315.90.81.890.70.98.70.0320.03811.2
MEDIAN0.5413.80.050.060.610.360.0704.72900.440.232428.30.275.80.81.780.70.76.40.0250.0268.6
0.10.256.50.030.020.170.100.0212.31000.120.071342.50.212.10.31.410.20.32.00.0120.0073.4
0.91.3776.40.100.092.100.660.1559.65521.460.5743216.90.497.61.52.611.31.617.30.0580.09024.8
RM
N = 14
MEAN0.3510.30.160.130.430.120.04813.82980.210.14925.50.731.72.30.641.81.87.30.0330.0122.2
MEDIAN0.3210.10.160.130.430.120.04412.72900.210.13915.50.681.82.30.561.81.87.10.0280.0112.6
0.10.198.60.120.090.290.080.0388.82300.140.10624.30.481.21.60.541.31.34.80.0220.0090.8
0.90.5411.70.220.180.620.140.05817.13540.270.161217.00.952.13.30.822.52.38.80.0470.0133.3
GE
N = 27
MEAN0.526.60.070.14n.a.0.080.0566.0231n.a.0.271306.10.604.34.70.521.1n.a.6.40.0610.01511.2
MEDIAN0.526.40.070.11n.a0.070.0546.1223n.a0.261335.90.624.44.40.441.0n.a5.80.0570.01510.7
0.10.364.80.060.08n.a.0.060.0414.9199n.a.0.231075.30.332.83.40.350.7n.a.5.40.0340.0137.4
0.90.688.70.080.25n.a.0.090.0777.1280n.a.0.321597.00.875.36.10.581.5n.a.7.50.0890.01815.0

Appendix D

Figure A2. Three-dimensional score plot obtained by PCA of the solubility percentage of elements (all data).
Figure A2. Three-dimensional score plot obtained by PCA of the solubility percentage of elements (all data).
Atmosphere 11 00479 g0a2

References

  1. D’Alessandro, A.; Lucarelli, F.; Mandò, P.A.; Marcazzan, G.; Nava, S.; Prati, P.; Valli, G.; Vecchi, R.; Zucchiatti, A. Hourly elemental composition and sources identification of fine and coarse PM10 particulate matter in four Italian towns. J. Aerosol. Sci. 2003, 34, 243–259. [Google Scholar] [CrossRef]
  2. Vecchi, R.; Marcazzan, G.; Valli, G.; Ceriani, M.; Antoniazzi, C. The role of atmospheric dispersion in the seasonal variation of PM1 and PM2. 5 concentration and composition in the urban area of Milan (Italy). Atmos. Environ. 2004, 38, 4437–4446. [Google Scholar] [CrossRef]
  3. Ariola, V.; D’Alessandro, A.; Lucarelli, F.; Marcazzan, G.; Mazzei, F.; Nava, S.; Garcia-Orellana, I.; Prati, P.; Valli, G.; Vecchi, R.; et al. Elemental characterization of PM10, PM2. 5 and PM1 in the town of Genoa (Italy). Chemosphere 2006, 62, 226–232. [Google Scholar] [CrossRef]
  4. Perrino, C.; Canepari, S.; Catrambone, M.; Dalla Torre, S.; Rantica, E.; Sargolini, T. Influence of natural events on the concentration and composition of atmospheric particulate matter. Atmos. Environ. 2009, 43, 4766–4779. [Google Scholar] [CrossRef]
  5. Amodio, M.; Andriani, E.; Cafagna, I.; Caselli, M.; Daresta, B.E.; de Gennaro, G.; Di Gilio, A.; Placentino, C.M.; Tutino, M. A statistical investigation about sources of PM in South Italy. Atmos. Res. 2010, 98, 207–218. [Google Scholar] [CrossRef]
  6. Carbone, C.; Decesari, S.; Mircea, M.; Giulianelli, L.; Finessi, E.; Rinaldi, M.; Fuzzi, S.; Marinoni, A.; Duchi, R.; Perrino, C. Size-resolved aerosol chemical composition over the Italian Peninsula during typical summer and winter conditions. Atmos. Environ. 2010, 44, 5269–5278. [Google Scholar] [CrossRef]
  7. Putaud, J.P.; Van Dingenen, R.; Alastuey, A.; Bauer, H.; Birmili, W.; Cyrys, J.; Flentje, H.; Fuzzi, S.; Gehrig, R.; Hansson, H.-C. A European aerosol phenomenology–3: Physical and chemical characteristics of particulate matter from 60 rural, urban, and kerbside sites across Europe. Atmos. Environ. 2010, 44, 1308–1320. [Google Scholar] [CrossRef]
  8. Contini, D.; Belosi, F.; Gambaro, A.; Cesari, D.; Stortini, A.M.; Bove, M.C. Comparison of PM10 concentrations and metal content in three different sites of the Venice Lagoon: An analysis of possible aerosol sources. J. Environ. Sci. 2012, 24, 1954–1965. [Google Scholar] [CrossRef] [Green Version]
  9. Perrino, C.; Catrambone, M.; Dalla Torre, S.; Rantica, E.; Sargolini, T.; Canepari, S. Seasonal variations in the chemical composition of particulate matter: A case study in the Po Valley. Part I: Macro-components and mass closure. Environ. Sci. Pollut. Res. 2014, 21, 3999–4009. [Google Scholar] [CrossRef] [PubMed]
  10. Canepari, S.; Astolfi, M.L.; Farao, C.; Maretto, M.; Frasca, D.; Marcoccia, M.; Perrino, C. Seasonal variations in the chemical composition of particulate matter: A case study in the Po Valley. Part II: Concentration and solubility of micro- and trace-elements”. Environ. Sci. Pollut. Res. 2014, 21, 4010–4022. [Google Scholar] [CrossRef]
  11. Sarti, E.; Pasti, L.; Rossi, M.; Ascanelli, M.; Pagnoni, A.; Trombini, M.; Remelli, M. The composition of PM1 and PM2. 5 samples, metals and their water soluble fractions in the Bologna area (Italy). Atmos. Pollut. Res. 2015, 6, 708–718. [Google Scholar] [CrossRef] [Green Version]
  12. Amato, F.; Alastuey, A.; Karanasiou, A.; Lucarelli, F.; Nava, S.; Calzolai, G.; Severi, M.; Becagli, S.; Gianelle, V.L.; Colombi, C. AIRUSE-LIFE+: A harmonized PM speciation and source apportionment in 5 Southern European cities. Atmos. Chem. Phys. 2016, 16, 3289–3309. [Google Scholar] [CrossRef] [Green Version]
  13. Perrone, M.G.; Zhou, J.; Malandrino, M.; Sangiorgi, G.; Rizzi, C.; Ferrero, L.; Dommen, J.; Bolzacchini, E. PM chemical composition and oxidative potential of the soluble fraction of particles at two sites in the urban area of Milan, Northern Italy. Atmos. Environ. 2016, 128, 104–113. [Google Scholar] [CrossRef] [Green Version]
  14. Ricciardelli, I.; Bacco, D.; Rinaldi, M.; Bonafè, G.; Scotto, F.; Trentini, A.; Bertacci, G.; Ugolini, P.; Zigola, C.; Rovere, F.; et al. A three-year investigation of daily PM2. 5 main chemical components in four sites: The routine measurement program of the Supersito Project (Po Valley, Italy). Atmos. Environ. 2017, 152, 418–430. [Google Scholar] [CrossRef]
  15. Benetello, F.; Squizzato, S.; Masiol, M.; Khan, M.B.; Visin, F.; Formenton, G.; Pavoni, B. A procedure to evaluate the factors determining the elemental composition of PM 2.5. Case study: The Veneto region (northeastern Italy). Environ. Sci. Pollut. Res. 2018, 25, 3823–3839. [Google Scholar] [CrossRef] [PubMed]
  16. Cesari, D.; De Benedetto, G.E.; Bonasoni, P.; Busetto, M.; Dinoi, A.; Merico, E.; Chirizzi, D.; Cristofanelli, P.; Donateo, A.; Grasso, F.M.; et al. Seasonal variability of PM2. 5 and PM10 composition and sources in an urban background site in Southern Italy. Sci. Total Environ. 2018, 612, 202–213. [Google Scholar] [CrossRef]
  17. Thusrton, G.D.; Spengler, J.D. A quantitative assessment of source contributions to inhalable particulate matter pollution in metropolitan Boston. Atmos. Environ. 1985, 19, 9–25. [Google Scholar]
  18. TÜV Certificate 0000028733_02; TÜV Rheinland Energy GmbH: Köln, Germany, 12 April 2018.
  19. Canepari, S.; Astolfi, M.L.; Moretti, S.; Curini, R. Comparison of extracting solutions for elemental fractionation in airborne particulate matter. Talanta 2010, 82, 834–844. [Google Scholar] [CrossRef]
  20. Canepari, S.; Cardarelli, E.; Pietrodangelo, A.; Strincone, M. Determination of metals, metalloids and non-volatile ions in airborne particulate matter by a new two-step sequential leaching procedure: Part B: Validation on equivalent real samples. Talanta 2006, 69, 588–595. [Google Scholar] [CrossRef]
  21. Engling, G.; Carrico, C.M.; Kreidenweis, S.M.; Collett, J.L., Jr.; Day, D.E.; Malm, W.C.; Lincoln, E.; Hao, W.M.; Iinuma, Y.; Herrmann, H. Determination of levoglucosan in biomass combustion aerosol by high-performance anion-exchange chromatography with pulsed amperometric detection. Atmos. Environ. 2006, 40, 299–311. [Google Scholar] [CrossRef]
  22. Perrino, C.; Canepari, S.; Cardarelli, E.; Catrambone, M.; Sargolini, T. Inorganic constituents of urban air pollution in the Lazio region (Central Italy). Environ. Monit. Assess. 2008, 136, 69–86. [Google Scholar] [CrossRef] [PubMed]
  23. Canepari, S.; Perrino, C.; Astolfi, M.L.; Catrambone, M.; Perret, D. Determination of soluble ions and elements in ambient air suspended particulate matter: Inter-technique comparison of XRF, IC and ICP for sample-by-sample quality control. Talanta 2009, 77, 1821–1829. [Google Scholar] [CrossRef] [PubMed]
  24. Leardi, R.; Melzi, C.; Polotti, G. CAT (Chemometric Agile Tool). Available online: http://gruppochemiometria.it/index.php/software (accessed on 1 February 2020).
  25. Chow, J.C.; Lowenthal, D.H.; Chen, L.W.; Wand, X.; Watson, J.G. Mass reconstruction methods for PM2.5: A review. Air Qual. Atmos. Health 2015, 8, 243–263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chan, Y.C.; Simpson, R.W.; McTainsh, G.H.; Vowles, P.D. Characterisation of Chemical Species in PM2.5 and PM10 Aerosols in Brisbane, Australia. Atmos. Environ. 1997, 31, 3773–3785. [Google Scholar] [CrossRef]
  27. Marcazzan, G.M.; Vaccaro, S.; Valli, G.; Vecchi, R. Characterisation of PM10 and PM2.5 Particulate Matter in the Ambient Air of Milan (Italy). Atmos. Environ. 2001, 35, 4639–4650. [Google Scholar] [CrossRef]
  28. Turpin, B.J.; Lim, H. Species Contribution to PM2.5 Mass Concentration: Revisiting Common Assumptions for Estimating Organic Mass. Aerosol. Sci. Technol. 2001, 35, 602–610. [Google Scholar] [CrossRef]
  29. Viidanoja, J.; Sillanpää, M.; Laakia, J.; Kerminen, V.M.; Hillamo, R.; Aarnio, P.; Koskentalo, T. Organic and black carbon in PM2.5 and PM10: 1 year of data from an urban site in Helsinki, Finland. Atmos. Environ. 2002, 36, 3183–3193. [Google Scholar] [CrossRef]
  30. Farao, C.; Canepari, S.; Perrino, C.; Harrison, R.M. Sources of PM in an Industrial Area: Comparison between receptor model Results and semiempirical calculations of source contributions. Aerosol Air Qual. Res. 2014, 14, 1558–1572. [Google Scholar] [CrossRef] [Green Version]
  31. Perrino, C.; Catrambone, M.; Pietrodangelo, A. Influence of atmospheric stability on the mass concentration and chemical composition of atmospheric particles: A case study in Rome, Italy. Environ. Int. 2008, 34, 621–628. [Google Scholar] [CrossRef]
  32. Squizzato, S.; Masiol, M.; Brunelli, A.; Pistollato, S.; Tarabotti, E.; Rampazzo, G. Factors determining the formation of secondary inorganic aerosol: A case study in the Po Valley (Italy). Atmos. Chem. Phys. 2013, 13, 1927. [Google Scholar] [CrossRef] [Green Version]
  33. Perrone, M.G.; Larsen, B.R.; Ferrero, L.; Sangiorgi, G.; De Gennaro, G.; Udisti, R.; Zangrando, R.; Gambaro, A.; Bolzacchini, E. Sources of high PM2.5 concentrations in Milan, Northern Italy: Molecular marker data and CMB modelling. Sci. Total Environ. 2012, 414, 343–355. [Google Scholar] [CrossRef] [PubMed]
  34. Nava, S.; Lucarelli, F.; Amato, F.; Becagli, S.; Calzolai, G.; Chiari, M.; Giannoni, M.; Traversi, R.; Udisti, R. Biomass burning contributions estimated by synergistic coupling of daily and hourly aerosol compo sition records. Sci. Total Environ. 2015, 511, 11–20. [Google Scholar] [CrossRef] [PubMed]
  35. Pietrogrande, M.C.; Bacco, D.; Ferrari, S.; Ricciardelli, I.; Scotto, F.; Trentini, A.; Visentin, M. Characteristics and major sources of carbonaceous aerosols in PM2. 5 in Emilia Romagna Region (Northern Italy) from four-year observations. Sci. Total Environ. 2016, 553, 172–183. [Google Scholar] [CrossRef] [PubMed]
  36. Perrino, C.; Tofful, L.; Dalla Torre, S.; Sargolini, T.; Canepari, S. Biomass burning contribution to PM10 concentration in Rome (Italy): Seasonal, daily and two-hourly variations. Chemosphere 2019, 222, 839–848. [Google Scholar] [CrossRef] [PubMed]
  37. Querol, X.; Alastuey, A.; Moreno, T.; Viana, M.M.; Castillo, S.; Pey, J.; Rodríguez, S.; Artiñano, B.; Salvador, P.; Sánchez, M.; et al. Spatial and temporal variations in airborne particulate matter (PM10 and PM2. 5) across Spain 1999–2005. Atmos. Environ. 2008, 42, 3964–3979. [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. Massimi, L.; Simonetti, G.; Buiarelli, F.; Di Filippo, P.; Pomata, D.; Riccardi, C.; Ristorini, M.; Astolfi, M.L.; Canepari, S. Spatial distribution of levoglucosan and alternative biomass burning tracers in atmospheric aerosols, in an urban and industrial hot-spot of Central Italy. Atmos. Res. 2020, 104904. [Google Scholar] [CrossRef]
  40. Wåhlin, P.; Berkowicz, R.; Palmgren, F. Characterisation of traffic-generated particulate matter in Copenhagen. Atmos. Environ. 2006, 40, 2151–2159. [Google Scholar] [CrossRef]
  41. Mahale, V.; Bijwe, J.; Sinha, S. Efforts towards green friction materials. Tribol. Int. 2019, 136, 196–206. [Google Scholar] [CrossRef]
  42. Canepari, S.; Pietrodangelo, A.; Perrino, C.; Astolfi, M.L.; Marzo, M.L. Enhancement of source traceability of atmospheric PM by elemental chemical fractionation. Atmos. Environ. 2009, 43, 4754–4765. [Google Scholar] [CrossRef]
  43. Sato, K.; Tamura, T.; Furuta, N. Partitioning between soluble and insoluble fractions of major and trace elements in size-classified airborne particulate matter collected in Tokyo. J. Environ. Monit. 2008, 10, 211–218. [Google Scholar] [CrossRef] [PubMed]
  44. Canepari, S.; Perrino, C.; Olivieri, F.; Astolfi, M.L. Characterisation of the traffic sources of PM through size-segregated sampling, sequential leaching and ICP analysis. Atmos. Environ. 2008, 42, 8161–8175. [Google Scholar] [CrossRef]
  45. Canepari, S.; Astolfi, M.L.; Catrambone, M.; Frasca, D.; Marcoccia, M.; Marcovecchio, F.; Massimi, L.; Rantica, E.; Perrino, C. A combined chemical/size fractionation approach to study winter/summer variations, ageing and source strength of atmospheric particles. Environ. Poll. 2019, 253, 19–28. [Google Scholar] [CrossRef] [PubMed]
  46. Dos Santos, M.; Gómez, D.; Dawidowski, L.; Gautier, E.; Smichowski, P. Determination of water-soluble and insoluble compounds in size classified airborne particulate matter. Microchem. J. 2009, 91, 133–139. [Google Scholar] [CrossRef]
  47. Font, A.; de Hoogh, K.; Leal-Sanchez, M.; Ashworth, D.C.; Brown, R.J.; Hansell, A.L.; Fuller, G.W. Using metal ratios to detect emissions from municipal waste incinerators in ambient air pollution data. Atmos. Environ. 2015, 113, 177–186. [Google Scholar] [CrossRef] [Green Version]
  48. Arditsoglou, A.; Samara, C. Levels of total suspended particulate matter and major trace elements in Kosovo: A source identification and apportionment study. Chemosphere 2005, 59, 669–678. [Google Scholar] [CrossRef]
  49. Mazzei, F.; D’alessandro, A.; Lucarelli, F.; Nava, S.; Prati, P.; Valli, G.; Vecchi, R. Characterization of particulate matter sources in an urban environment. Sci. Total Environ. 2008, 401, 81–89. [Google Scholar] [CrossRef]
  50. Nzihou, A.; Stanmore, B. The fate of heavy metals during combustion and gasification of contaminated biomass—A brief review. J. Hazard. Mater. 2013, 256, 56–66. [Google Scholar] [CrossRef] [Green Version]
  51. Cesari, D.; Genga, A.; Ielpo, P.; Siciliano, M.; Mascolo, G.; Grasso, F.M.; Contini, D. Source apportionment of PM2. 5 in the harbour–industrial area of Brindisi (Italy): Identification and estimation of the contribution of in-port ship emissions. Sci. Total Environ. 2014, 497, 392–400. [Google Scholar] [CrossRef]
Figure 1. Location of the sixteen sampling sites.
Figure 1. Location of the sixteen sampling sites.
Atmosphere 11 00479 g001
Figure 2. Chemical composition of PM10 at the sixteen sites during the cold (a) and the warm (b) period.
Figure 2. Chemical composition of PM10 at the sixteen sites during the cold (a) and the warm (b) period.
Atmosphere 11 00479 g002
Figure 3. Principal Component Analysis (PCA) of elemental concentration (all data): score plot (a) and loading plot (b) of Components 2 and 3. _C and _W after the site code indicate the cold and warm period, respectively. In the case of MF, the two winter periods are reported as _C1 and _C2.
Figure 3. Principal Component Analysis (PCA) of elemental concentration (all data): score plot (a) and loading plot (b) of Components 2 and 3. _C and _W after the site code indicate the cold and warm period, respectively. In the case of MF, the two winter periods are reported as _C1 and _C2.
Atmosphere 11 00479 g003
Figure 4. PCA of the soluble fraction of elements (all data): the score plot (a) and loading plot (b) of Components 2 and 3. The site codes are as seen in Figure 3.
Figure 4. PCA of the soluble fraction of elements (all data): the score plot (a) and loading plot (b) of Components 2 and 3. The site codes are as seen in Figure 3.
Atmosphere 11 00479 g004
Figure 5. PCA of the solubility percentage of elements (all data). (a) Score plot and (b) Loading plot; For MF, FE, RB, RM (inland sites), the cold season is shown in green, and the warm season is shown in red. For LS, CV, VG, ST, GE (coastal/ petrochemical), all data are shown in black.
Figure 5. PCA of the solubility percentage of elements (all data). (a) Score plot and (b) Loading plot; For MF, FE, RB, RM (inland sites), the cold season is shown in green, and the warm season is shown in red. For LS, CV, VG, ST, GE (coastal/ petrochemical), all data are shown in black.
Atmosphere 11 00479 g005
Table 1. Overview of the sampling sites and sampling periods. W and C indicate the warm and cold semesters, respectively.
Table 1. Overview of the sampling sites and sampling periods. W and C indicate the warm and cold semesters, respectively.
LocationSampling SiteYearSampling Duration (Weeks) and Season (W, C)
City/VillageClassificationResidentsDistance1 (km) Altitude (m a.s.l.)N° of Sites
Chamonix (F) CHresidential9400-1035120082 W + 2 C
Monfalcone MFperi-urban/ind.28,500-<1042014–20164 W + 4 C
Ferrara FEperi-urban/ind.131,8005<1022013–20156 W + 6 C
Rubbiano RBresidential /ind.500-16612015–20163 W + 3 C
La Spezia LSurban/ind.93,000-<10120122 W
Civitavecchia CVperi-urban/ind.52,5005<10120163 W
M.rotondo sc. MRperi-urban/agric.850072012014–20158 W + 8 C
Guidonia GU urban90,000-100120082 W + 4 C
Rome RMurban2,900,000-2112017–20182 W + 6 C
Fontechiari FCagricultural1300-37512007–200815 W + 9 C
Frosinone FRperi-urban/ind.46,000516132007–20082 W + 6 C
Brindisi BRperi-urban/ind.86,000101522009–20103 C
Viggiano VGresidential/ind.3400-97512017–20184 W + 4 C
Strongoli m. STresidential/ind.2500-<1022011–20122 W + 3 C
Palermo PAurban660,000-14120118 W
Gela GEperi-urban/ind.46,50034622013–20142 W + 2 C
1 Distance of the sampling site to the main center.

Share and Cite

MDPI and ACS Style

Perrino, C.; Catrambone, M.; Canepari, S. Chemical Composition of PM10 in 16 Urban, Industrial and Background Sites in Italy. Atmosphere 2020, 11, 479. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11050479

AMA Style

Perrino C, Catrambone M, Canepari S. Chemical Composition of PM10 in 16 Urban, Industrial and Background Sites in Italy. Atmosphere. 2020; 11(5):479. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11050479

Chicago/Turabian Style

Perrino, Cinzia, Maria Catrambone, and Silvia Canepari. 2020. "Chemical Composition of PM10 in 16 Urban, Industrial and Background Sites in Italy" Atmosphere 11, no. 5: 479. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11050479

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