Contribution of Volcanic and Fumarolic Emission to the Aerosol in Marine Atmosphere in the Central Mediterranean Sea: Results from Med-Oceanor 2017 Cruise Campaign

: This work studied the contribution of the geogenic sources volcanoes and fumaroles to the aerosol in marine atmosphere in the central Mediterranean basin. For this purpose, in the framework of the Med-Oceanor measurement program, we carried out a cruise campaign in the summer of 2017 to investigate the impact to the aerosol of the most important Mediterranean volcanoes (Mount Etna, Stromboli Island, and Marsili Seamount) and solfatara areas (Phlegraean Fields complex, Volcano Islands, Ischia Island, and Panarea submarine fumarole). We collected PM 10 and PM 2.5 samples in 12 sites and performed chemical characterization to gather information about the concentration of major and trace elements, elemental carbon (EC), organic carbon (OC), and ionic species. The use of triangular plots and the calculation of enrichment factors conﬁrmed the interception of volcanic plume. We integrated the outcomes from chemical characterization with the use of factor analysis and SEM / EDX analysis for the source apportionment. Anthropogenic and natural sources including shipping emissions, volcanic and fumarolic load, as well as sea spray were identiﬁed as the main factors a ﬀ ecting aerosol levels in the study area. Furthermore, we performed pattern recognition analysis by stepwise linear discriminant analysis to seek di ﬀ erences in the composition of PM 10 and PM 2.5 samples according to their volcanic or solfatara origin.


Introduction
Advancements in the current state of knowledge of atmospheric aerosols are an important research topic because of their implications in environmental and health issues [1]. Nowadays, much of the air pollution results from human activities, but natural sources also contribute to the increase of the exposure level of the population and ecosystem to the polluted air. Typically, the major natural emissions in terms of mass include sea spray (84%) and mineral dust (13%), along with other sources such as biological primary organic aerosols (POA), volcanic emissions, biogenic secondary organic aerosols (SOA), and biogenic sulfate particles [2]. Anthropogenic aerosols contribute only 2% to global

The Mount Etna
Mount Etna, one of the most active volcanoes in the world, covers an area of about 1250 km 2 and reaches an altitude of about 3340 m a.s.l. (37.754193° N, 14.978929° E) ( Figure 1). It has been active during historical time, with frequent paroxysmal episodes separated by passive degassing periods. The edifice consists of a lower shield unit overlain by a stratovolcano [27] that is truncated at 2800-2900 m of altitude by the Ellittico caldera and Cratere del Piano, while the upper part of the volcano culminates with a large summit cone endowed with four active craters (Bocca Nuova, Voragine, Northeast, and Southeast Craters). The general atmospheric circulation in the Etnean area is

The Mount Etna
Mount Etna, one of the most active volcanoes in the world, covers an area of about 1250 km 2 and reaches an altitude of about 3340 m a.s.l. (37.754193 • N, 14.978929 • E) ( Figure 1). It has been active during historical time, with frequent paroxysmal episodes separated by passive degassing periods. The edifice consists of a lower shield unit overlain by a stratovolcano [27] that is truncated at 2800-2900 m of altitude by the Ellittico caldera and Cratere del Piano, while the upper part of the volcano culminates with a large summit cone endowed with four active craters (Bocca Nuova, Voragine, Northeast, and Southeast Craters). The general atmospheric circulation in the Etnean area is predominantly westerly to northwesterly. According to these considerations, sampling sites in this zone were located in front of the Acireale and Acicastello coasts.

Stromboli
Stromboli is an island volcano belonging to the Eolian volcanic arc (926 m a.s.l; 38.7891 • N, 15.2131 • E) ( Figure 1). It is an isolated volcano characterized by permanent activity and continuous gas emissions from potassium-rich basalt. Although Stromboli has a significant hydrothermal system [28], its magma also degasses directly into the atmosphere (i.e., at very high temperature, >1000 • C) from lava-filled vents. Such systems are often described as "open-conduit" volcanoes and can degas very large quantities of unerupted magma. Three different sampling sites around the island of Stromboli were monitored labeled as Stromboli (1), Stromboli (2), and Stromboli (3).

Marsili
Marsili Seamount is about 60 km long and 20 km wide, rising about 3000 m from the seafloor at over 3500 m water depth (39.236276 • N, 14.369875 • E) ( Figure 1). The possible existence of a very large underwater explosive volcano, together with some encouraging clues that point to Marsili as an important and potential long-lasting-renewable energy resource [29], has reinforced research and exploitation efforts in this area. The existence of an active magmatic chamber at about 2.5 km below the summit was proposed on the basis of petrologic studies of basalts [30], and it is compatible with gravimetric and magnetic data modeling [31], as well as geochemical observations on the summit [32]. We monitored this area by a sampling stop above the crater of the volcano.

Phlegraean Fields
Solfatara, a crater located within the Phlegraean Fields (458 m a.s.l.; 40.8271 • N, 14.1391 • E), a quaternary volcanic complex near Naples (Figure 1), is renowned for its vigorous, low-temperature degassing. Indeed, the term Solfatara is applied to many similar manifestations found in other volcanic regions of the world. The hottest fumaroles at Solfatara reach temperatures of 140-160 • C. Previous studies indicate that gas emissions are dominated by H 2 O, followed by CO 2 and H 2 S [33]. In this area, sampling was performed at two different sites (labeled as Phlegraean Fields (1) and Phlegraean Fields (2)) along the direction of the inland Solfatara crater at a distance of about 1 km.

Vulcano
The volcanic island of Vulcano is the southernmost of the seven islands forming the Eolian Archipelago in the Tyrrhenian Sea north of Sicily (500 m a.s.l.; 38.4041 • N, 14.9621 • E). Since the last eruption of Vulcano in 1888-1890, a vigorous fumarole field has been developed in the summit crater, known as La Fossa. Numerous geochemical investigations have focused on Vulcano, especially since increases in gas fluxes and temperatures were observed in the 1970s. Similarly to the Solfatara site, Vulcano does not have magma on the surface but its fumaroles reach high temperatures, well over 300 • C at some vents. In Vulcano, we sampled at two sites labeled as Vulcano (1) and Vulcano (2) (Figure 1).

Ischia
Ischia (40.7352 • N, 13.8517 • E) is the westernmost active volcanic complex of the Campania area and belongs to the Phlegraean volcanic district, which also includes Phlegraean Fields (Figure 1). In Ischia, the last eruption took place in 1301, and since that time, several earthquakes affected the island. Thermal manifestations characterize Ischia including thermal waters and fumaroles, which have been well known since Roman times and are the main economic resource of the island at present. The fumarolic activity is mainly concentrated along the faults affecting Mt. Epomeo (789 m a.s.l.), a resurgent block in the central part of the island. In this area, we performed aerosol sampling in a site labeled Ischia (Figure 1).

Panarea
Panarea Island (Eolian Archipelago) lies in the western sector of a submarine 120 m deep platform constituting the 56 km 2 wide summit of a 2000 m high seamount (38.6369 • N, 15.1026 • E) ( Figure 1). Close to Panarea, a group of islets (Dattilo, Panarelli, Lisca Bianca, Bottaro, Lisca Nera, and Le Formiche) forms an archipelago that surrounds a submerged fumarolic area, which is characterized by relatively shallow depth (<35 m), hosting gas, and thermal water discharges. Geochemical studies of fluids from the submarine fumarolic field of this area started in the early 1990s [34], and were then fostered by the gas burst event that occurred on 3 November 2002 [35]. One sampling site, labeled Panarea, was selected offshore as a monitoring point.

Sampling and Laboratory Analysis
The aerosol sampling equipment was fixed on the deck at about 30 m above the sea level on the front of the ship to avoid contamination from the ship exhaust. Daily 24-h ambient PM 10 and PM 2.5 samples were collected on quartz fiber filters of 47 mm in diameter (2500 QAT-UP, Pallflex ® , Putnam, CT, USA) using a low volume sampler Echo PM (, TCR Tecora, Cogliate, Italy) operating at a flow of 2.3 m 3 h −1 . In accordance with the European Standard UNI EN 12341:2014, filters were conditioned for 48 h before and after sampling, at a temperature and relative humidity of 20 ± 1 • C and 50 ± 5%, respectively. The conditioned filters were weighed using a microbalance Crystal Micro (Gibertini, Novate, Italy) in order to determine the mass of both the PM fractions by standard gravimetric procedures, whereupon the filters were analyzed for the determination of the chemical composition. The analysis of the carbonaceous fraction of aerosol (TC = OC + EC) was performed on a 1 cm 2 portion of each filter by a thermal-optical transmission (TOT) technique using an OC/EC Analyzer (Sunset Laboratory Company, Tigard, OR, USA) and following the EUSAAR-II temperature protocol [36]. Briefly, the small punch of the filter was subjected to two stepwise heating cycles, which favored the volatilization and combustion processes of the carbonaceous material. The evolved carbon was converted to methane and subsequently detected by a flame ionization detector (FID) during both temperature cycles, which occurred in a non-oxidizing and oxidizing atmosphere, respectively. Laser light of 670 nm continuously monitored the transmittance through the sample filter in order to correct for pyrolyzed OC conversion to EC.
Chemical characterization of PM samples was performed by Inductively Coupled Plasma Atomic Emission Spectrometry (iCAP 6500, Thermo Scientific) for major elements (Al, Ca, Na, Mg, K, Fe, P, S) and Inductively Coupled Plasma Mass Spectrometry (iCAP-RQ, Thermo Scientific) for trace elements (Li, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, As, Se, Rb, Sr, Cd, Sn, Sb, Ba, La, Ce, Hf, Pb, Bi, Th, U). Half of each filter was acid digested with a mixture of HF:HNO 3 :HClO 4 (2.5:1.25:1.25 mL) in a Teflon reactor for 9 h at a temperature of 90 • C. Each digested sample was driven to dryness, re-dissolved with 1.25 mL HNO 3 , and made up to a final volume of 25 mL with ultrapure water [37]. Relative analytical errors were estimated using a small amount (15 mg) of the NIST-1633b (fly ash) reference material loaded on a 1 4 quartz fiber filter. Individual uncertainties were calculated taking into account the analytical uncertainty and the standard deviations of the element concentrations in the blank filters [38]. Since the portions of the PM 2.5 filters for the sites Acireale, Panarea, and Vulcano (2) were contaminated during the storage, we discarded these samples before the chemical analysis for major and trace elements. For the quantification of ammonium by ion-selective electrode (Orion 9512HPBNWP of Thermo Scientific) and soluble ions, such as sulfate, nitrite, nitrate, fluoride, bromide, chloride, and phosphate, by ion chromatography (Aquion, Thermo Scientific), the leftover filter was water leached in 20 mL sealed PVC bottles (for 14 h at 60 • C), preceded by an ultrasonic bath of 15 min. Sulfur dioxide (SO 2 ) measurements were performed using a Teledyne UV Fluorescence SO 2 Analyzer (API Model 100EU) with a sampling flow rate of 0.6 L min −1 and a time resolution of 1 min. A particle sampling cassette system (SKC, inlet diameter 1/8 inch), which housed transmission electron microscope (TEM) grids attached to filter support pads (25 mm in diameter), was used to collect particles for offline characterization following a similar procedure described elsewhere [39,40]. The system utilized an SKC Leland pump that was operated at a flow of 5 L min −1 for a duration of 4 h. The grids (Agar Scientific Quantifoil 200 Mesh Au) were analyzed offline for morphological and physicochemical particle characterization using a Field Emission Scanning Electron Microscope (JEDL J-7100) coupled with an energy-dispersive X-ray spectroscopy (EDS; retroscattered electron detector) spectrometer.

Chemical Mass Closure
For the purpose of mass closure, the chemical components were grouped as sea salt (ss), mineral dust, organic matter and elemental carbon (OM + EC), secondary inorganic aerosols (SIA), and trace elements. Each of these components can originate from a variety of sources, some of which play a significant role in their production.
As regards to the sea salt component, six major ions represent more than 99% of the mass of salts dissolved in seawater: four cations, i.e., sodium (ss-Na + ), magnesium (ss-Mg 2+ ), calcium (ss-Ca 2+ ), and potassium (ss-K + ), and two anions, chloride (ss-Cl − ) and sulfate (ss-SO 4 2− To calculate the sea salt sodium [ss-Na + ], the value of non-sea salt sodium (nss-Na) was subtracted from the total Na value ([ss-Na + ] = [Na] − [nss-Na]) resulting from the chemical analysis. The [nss-Na] value was obtained from the Al concentration multiplied by the Na/Al ratio in the average crustal composition [42]. Based on the seawater composition, the sea salt sulfate [ss−SO 4 2− ] is calculated as Mineral dust is generated from the suspension of the minerals constituting the soil and can be expressed as the sum of the most common oxide forms [37]. In our study, dust content was estimated using Equation (2) as the sum of the concentrations of Al, Ca, Mg, K, Fe, P, and Ti, each multiplied by factors to convert them to their common oxides and SiO 2 , calculated multiplying [Al 2 O 3 ] by 2.5.
EC was obtained from thermal-optical analysis (as described in Section 3.1), while the OM contribution, which accounts for the unmeasured H, O, N, and S in organic compounds, was calculated according to Equation (3), which considers also the influence of the water-soluble organic carbon (WSOC) on the mass, that for non-urban aerosol is about 40% of the OC mass.
In this equation, OC was experimentally determined by thermal-optical analysis while the f multiplier is a coefficient that is not site-or time-specific. Depending on the extent of OM oxidation and secondary organic aerosol (SOA) formation, values for f vary from 1.2 for fresh aerosol in urban areas [43] to 2.6 for aged aerosol [44][45][46]. The conversion factor f is usually taken in the range from 1.4 for urban aerosols to 1.8 for remote aerosols. Turpin et al. revisited these conversion factors and proposed values of 1.6 ± 0.2 and 2.1 ± 0.2 for urban and non-urban aerosols, respectively [46]. Based on the previous consideration, the f value was set at 2.1 for the calculations in our work.  [37,47,48].
The last contribution to the PM mass closure, which derives from trace elements, was obtained as the sum of the trace element concentrations. The sum of all the above determinations accounted for 81.9% and 83.0% of the PM 10 and PM 2.5 average mass, respectively. The remaining undetermined mass can be attributed to the structural and adsorbed water that is not removed during the sample conditioning.

Data Analysis and Receptor Modeling
In order to characterize the daily atmospheric scenarios influencing PM levels, we used several complementary tools to interpret the different sources of air masses in the studied area. We calculated air mass backward trajectories using the US NOAA HYSPLIT model [49] with GDAS meteorological data. The occurrence of African dust outbreaks was also confirmed by the same tools coupled with the information of the NAAPS aerosol maps from the Naval Research Laboratory, USA (NRL) (Navy Aerosol Analysis and Prediction System) and the BSC-DREAM dust maps [50]. Enrichment factors (EFs) and factor analysis (FA) were used for source identification and apportionment of particulate air pollutants at receptor sites, whereas linear discriminant analysis (LDA) was exploited to perform supervised pattern recognition analysis using the chemical descriptors of PM samples according to their volcanic (Vo) or solfatara (So) origin.
The calculation of EFs is useful to rank the enrichment degree of an element compared to a known source; in this case, the Al concentration was considered as the basis for comparison. EF is defined as the ratio between the concentration of an element i (E i sample ) and that of a reference element (Al sample ) in the sample, normalized to the same ratio in the Earth's crust [42]. Specifically, for each element i, the EF i was calculated according to Equation (4).
where volc refers to the volcanic rocks (basalt and lava in this case), which were considered as the geogenic source. In our study, the EF i in Etna and Stromboli plumes was calculated using the aluminum as normalizing element, while the average Etnean lava [23] and the Stromboli basalt [51] were used, respectively, as reference material. Factor analysis (FA) is widely used to determine the number of independent sources that contribute to the system and elemental source profiles. The association of the factors with the source type is possible by comparing the elements with the highest loading within each component with the elements commonly emitted by recognized source types. In our work, the whole dataset was subjected to factor analysis using the software Statistica (StatSoft,) and Varimax rotation was applied to assure the maximum differentiation between the factors. The number of potential factors was limited to eight, but only six factors were obtained.
Linear discriminant analysis (LDA) is a chemometric tool that defines a set of delimiters that divide the multivariate space of the samples into as many subspaces as the number of the classes. Discriminant functions are computed as a linear combination of variables that maximizes the ratio of variance between categories to variance within categories. An essential restriction in the application of LDA concerns the ratio between the number of samples and the number of variables. Indeed, it is recommended that the number of samples is at least three times the number of variables to obtain a robust model. Therefore, in our study, a stepwise approach (S-LDA) was used for the selection of the variables with a major discriminant power. The S-LDA was used to differentiate between the PM samples collected in the volcanic (Vo) and solfatara (So) sites on the basis of their chemical composition. The dataset was slightly different from that used for FA because, for LDA, the variables that can be associated exclusively with an anthropogenic origin (Mo, EC, and OC) were not included in the dataset. Furthermore, the Acireale sample was excluded because it was affected by the Saharan Dust Outbreak (SDO), as pointed out from chemical characterization (see Section 4.2). Partial Wilks' λ values were used to figure out the variables that actually contributed to differentiation, i.e., the closer to zero its value is, the higher is the discriminatory importance of this variable in the model.

Meteorological Conditions
The trends of temperature, pressure, and wind speed recorded during the cruise by the on-board instrumentation are illustrated in Figure S1. During all the sampling days, the absence of precipitations, low wind intensity (with the maximum speed being 13.5 m s −1 ), and a mean atmospheric temperature of 26.1 • C characterized the local weather. Information on wind direction and intensity are reported for each sampling site by the wind rose diagrams ( Figure S2). Based on the average daily wind speed values and the Beaufort wind scale, we encountered light breeze for eight sampling days and gentle breeze for four sampling days. The atmospheric pressure exhibited the lowest average values and the largest variability between 1 September and 4 September 2017, with a minimum on 2 September, associated with episodes of strong wind and rough sea that forced us to stop in the Harbor of Lipari Island to shelter from the bad weather. The NCEP/NCAR Reanalysis [52] tool shows the composite means and anomalies of the sea level pressure over the Mediterranean basin during the campaign, supporting the worsening of the weather ( Figure S3). Indeed, while in August the expansion of the Azores Anticyclone characterized the synoptic conditions in line with seasonal climatology in the Mediterranean, on 2 September the anticyclonic system was confined over the Atlantic Ocean and North of Europe, thus favoring the development of low-pressure systems across Eastern Europe and the Mediterranean basin.

Percentage Chemical Composition of PM
The mass concentration values for PM 10 and PM 2.5 recorded during the campaign ranged between 12.3 and 30.7 µg m −3 and 9.4 and 18.3 µg m −3 , respectively ( Figure 2). The PM 10 concentration in each site was lower than the daily limit value of 50 µg m -3 , as set by the Directive 2015/1480/EC on ambient air quality. It is known that natural episodes, although they occur less frequently than anthropic ones, can influence air quality especially in the surroundings of marine areas and zones with tectonic activity [2]. The percentage chemical composition of the PM samples collected during the cruise is summed up in Table S1 and illustrated in Figure 2.
For the PM 10 , the main components were carbonaceous (OM + EC) and mineral dust aerosols, with average contributions of 27.2% and 26.6%, respectively. The dust load reached a peak of 44.7% at Acireale, owing to a North African dust outbreak as confirmed by NAAPS, BSC-DREAM, and backward trajectories maps presented in Figure S4. The carbonaceous fraction was above the 20% for the whole campaign with the maximum amount of 38.3%, which was found in Phlegraean Fields (1) and may be due to the interception of air masses moving from the urbanized coast near the sampling area, which is characterized by high vehicular traffic. This assumption is supported by the levels of EC (0.95 µg m −3 ) and OC (3.64 µg m −3 ), which were found at their highest concentrations at this stop. The SIA component had an average mass contribution to PM 10 of 19.9% with a maximum of 26.3% at Stromboli (2). Trace elements and sea salt contributed marginally to the mass, with tracers that accounted for less than 1% and sea salt for about 10%. We detected the maximum value of trace elements in the Ischia sample, where it represented 1.6% of PM 10 mass. The relatively low contribution of sea salt compared to what can be expected in a marine atmosphere may be attributed to the stable weather encountered for most of the campaign, except for rough sea and strong blowing wind arose at Vulcano (2), where we recorded a higher percentage of the sea salt component (23.7%).
Indeed, the occurrence of substantial sea spray production needs wind with an intensity and duration strong enough to allow its formation [53]. For PM 2.5 , the component prevailing on average was the carbonaceous aerosols (35.3%), which reached the maximum value in the Phlegraean Fields (1) (56.2%), thus reinforcing the preceding hypothesis involving the vehicular traffic contribution in this area. For SIA, the average loading reached 24.7%, thereby showing that its components (i.e., SO 4 2− , NH 4 +, and NO 3 -) prevailed in the PM 2.5 [37,54]. In both fractions, the maximum value of SIA was obtained in the sample Stromboli (2), where the upwind sampling (wind roses presented in Figure S2) led to a better interception of volcanic emissions characterized by high SO 2 values, which affected the increase of SO 4 2− concentration because of its reaction with the atmospheric water and the condensation on particles. In support of this, as showed in the time series plots of SO 2 ( Figure S5), we recorded a peak of SO 2 at this site, which lasted for about 5 h during that sampling day. As expected, the average mineral dust contribution was smaller in PM 2.5 than in PM 10 (17.6% and 25.6%, respectively), with the highest value reached in Acireale because of the SDO ( Figure S4). Average sea salt contribution to PM 2.5 was 4.3%, with the maximum value in Marsili (10.7%), while the trace elements were averagely higher in PM 2.5 (1.2%) than PM 10 (0.9%), with the maximum values recorded in Phlegraean Fields (1) and Phlegraean Fields (2) (1.4% and 1.6%, respectively).
precipitations, low wind intensity (with the maximum speed being 13.5 m s −1 ), and a mean atmospheric temperature of 26.1 °C characterized the local weather. Information on wind direction and intensity are reported for each sampling site by the wind rose diagrams ( Figure S2). Based on the average daily wind speed values and the Beaufort wind scale, we encountered light breeze for eight sampling days and gentle breeze for four sampling days. The atmospheric pressure exhibited the lowest average values and the largest variability between 1 September and 4 September 2017, with a minimum on 2 September, associated with episodes of strong wind and rough sea that forced us to stop in the Harbor of Lipari Island to shelter from the bad weather. The NCEP/NCAR Reanalysis [52] tool shows the composite means and anomalies of the sea level pressure over the Mediterranean basin during the campaign, supporting the worsening of the weather ( Figure S3). Indeed, while in August the expansion of the Azores Anticyclone characterized the synoptic conditions in line with seasonal climatology in the Mediterranean, on 2 September the anticyclonic system was confined over the Atlantic Ocean and North of Europe, thus favoring the development of low-pressure systems across Eastern Europe and the Mediterranean basin.

Percentage Chemical Composition of PM
The mass concentration values for PM10 and PM2.5 recorded during the campaign ranged between 12.3 and 30.7 µg m −3 and 9.4 and 18.3 µg m −3 , respectively ( Figure 2). The PM10 concentration in each site was lower than the daily limit value of 50 µg m -3 , as set by the Directive 2015/1480/EC on ambient air quality. It is known that natural episodes, although they occur less frequently than anthropic ones, can influence air quality especially in the surroundings of marine areas and zones with tectonic activity [2]. The percentage chemical composition of the PM samples collected during the cruise is summed up in Table S1 and illustrated in Figure 2.

Volcanic Areas
There are several studies reporting the relevance of volcanic emissions as sources of metal-bearing gases and PM into the atmosphere [51,[54][55][56][57]. A previous study revealed that volcanogenic trace metals are injected in the atmosphere as volatile elements, emitted from the silicate melt during magma degassing, and later dispersed in atmospheric plumes as sub-micron sized volcanic aerosols and low volatile (refractory or lithophile) elements, transported via coarse volcanic ash fragments [23]. The volcanic area monitored during the campaign, as previously discussed, included Etna and Stromboli, which are sub-aerial volcanoes with a significant part of their structure under the sea, while the Marsili Seamount is submerged and its activity is still open to debate. In this regard, at the Marsili site, we observed low values for each component constitutive of PM, except for sea salt in PM 2.5 , which reached the highest percentage in this site (10.7%) due to the low contributions of the other pollution sources because the sampling site was over 100 km from inland.
In our study, we used Al, Ca, and Fe, as tracers of crustal emissions to find out the volcanic-derived contribution to PM by comparing their average concentrations in the samples with the literature reference data. In particular, for information about the chemical composition of Etna's and Stromboli's lava and ash, we referred to the work of Calabrese et al. [23] and Allard et al. [51], respectively, while for the data on the chemical composition of Sahara and Sahel dust, we used the values provided by Moreno et al. [58]. The relative proportions of Al/Ca/Fe measured at Acireale and Acicastello are similar to the typical Etna's values in both ash and lava, confirming that samples, though collected at a considerable distance from the source (almost 40 km), are affected by the volcanic plume that spreads throughout the area (Figure 3a). The value of Acireale shown in Figure 3a is shifted toward the values of Sahara and Sahel, thus confirming the intrusion of air masses coming from the African desert region during the sampling. The ternary plot for the samples collected around Stromboli shows a similar composition with the reference values of magma reported in the literature [51]. In particular, the samples identified as Stromboli (2) and Stromboli (3) have a greater similitude with the composition of the basaltic lava, whereas the composition for the sample Stromboli (1) is closer to that of the volcanic ash ( Figure 3b). volcanic plume that spreads throughout the area (Figure 3a). The value of Acireale shown in Figure  3a is shifted toward the values of Sahara and Sahel, thus confirming the intrusion of air masses coming from the African desert region during the sampling. The ternary plot for the samples collected around Stromboli shows a similar composition with the reference values of magma reported in the literature [51]. In particular, the samples identified as Stromboli (2) and Stromboli (3) have a greater similitude with the composition of the basaltic lava, whereas the composition for the sample Stromboli (1) is closer to that of the volcanic ash (Figure 3b). Besides, the load of the volcanic plume aerosol in the sampled PM was investigated through the calculation of the enrichment factors, which is a useful tool to highlight the enrichment or depletion degree of an element in a geochemical medium relative to reference material. Although the use of EFs could hide some flaws, such as the variable composition of the reference material, or inadequate Besides, the load of the volcanic plume aerosol in the sampled PM was investigated through the calculation of the enrichment factors, which is a useful tool to highlight the enrichment or depletion degree of an element in a geochemical medium relative to reference material. Although the use of EFs could hide some flaws, such as the variable composition of the reference material, or inadequate choice of the element used to normalize the data [59], it provides a first-order insight on the geogenic or anthropogenic origin of the elements. Based on the results of the ternary plot, we calculated the EFs for Acicastello and Stromboli (3), which are the samples identified as those with the greatest volcanic contribution. EFs were computed according to Equation (4) using Al as a reference element because it is a low volatile element at magmatic conditions, and also one of the most immobile elements during the weathering of basalts [60]. Values of EF close to one indicate that there is no enrichment occurring for an element in sampled PM compared to the expected abundance in volcanic basalt. EFs far exceeding one characterize volatile elements enriched in the plume because of their affinity for the magmatic vapor phase or might indicate an anthropogenic origin. Despite some peculiarities, the EFs for the samples Stromboli (3) and Acicastello are overall consistent. The EFs illustrated in Figure S6 suggest that, for both volcanoes, the elements can be grouped into four categories according to the EF values. The strongly enriched elements in PM (EFs > 100) were those with a marked volcanic-volatile character or anthropogenic origin species (i.e., Sn, Sb, Cr, Se). The enriched elements such as Ni, Zn, U, As, Pb, and Tl were those with EFs ranging from 10 to 100, while the elements with EF ≈ 1-10 such as Cu, K, Ba, Rb, and some rare earth elements including Th, Sm, Ce, and La were associated with a geogenic origin. The elements with EFs < 1, such as Co, Ti, Mg, and Sr, can be considered depleted elements compared to geogenic elements and indicate a local geological reduction relative to volcanic rocks. The volatile-anthropogenic species, including Sb, Cr, and Sn, were enriched by two to three orders of magnitude in the PM of both volcanic sites. Furthermore, another enriched element was Ni which may be related to anthropogenic sources and is commonly associated, together with V, to shipping emissions. Among the enriched elements, the presence of Tl only in the Etna sample pointed out the volcanic contribution to PM [23].

Fumarole Areas
Sulfur dioxide is a valuable tracer of degassing processes by geothermal sources such as fumaroles and volcanoes. In this work, we averaged on a daily basis the 1 min recorded data obtained from the SO 2 analyzer ( Figure S5 As concerns Panarea, although it is a submarine fumarole source, a noteworthy value of both SO2 and nss-SO4 2− was detected, confirming the interception of fumarole degassing process. On the contrary, the analysis of trace metals did not present a significant contribution in PM of the elements typically emitted by natural sources, since they possibly dissolved in the path from the emission site and the seawater surface.

Factor Analysis and Source Identification
We used factor analysis, a tool for source identification and apportionment of particulate air pollutants, to characterize the different sources of aerosols during the campaign. We subjected the dataset to Varimax-rotated factor analysis identifying six factors with Eigenvalues greater than one ( Table 1). The six factors together explain 96.9% of the total variance of our dataset. As concerns Panarea, although it is a submarine fumarole source, a noteworthy value of both SO 2 and nss-SO 4 2− was detected, confirming the interception of fumarole degassing process. On the contrary, the analysis of trace metals did not present a significant contribution in PM of the elements typically emitted by natural sources, since they possibly dissolved in the path from the emission site and the seawater surface.

Factor Analysis and Source Identification
We used factor analysis, a tool for source identification and apportionment of particulate air pollutants, to characterize the different sources of aerosols during the campaign. We subjected the dataset to Varimax-rotated factor analysis identifying six factors with Eigenvalues greater than one ( Table 1). The six factors together explain 96.9% of the total variance of our dataset.

SEM/EDX Analysis
SEM/EDX analysis is utilized as a complementary tool for the quantitative chemical analysis (ICP). The elemental characteristics of aerosols are apportioned based on the respective sources. SEM/EDX confirms that particles with the specified compositions were present in the sampling area and their morphology is defined. Aerosols with a composition common for volcanic particles (e.g., SiO 2 , Mg, and Fe) have been observed, indicating that they may have originated from the volcanic plume (Figure 5a,d).
The anthropogenic factor has influenced the sampling (Figure 5f,j) and consequently the experimental results both in terms of chemical and microscopy analysis interpretation (i.e., factor and SEM analysis). In agreement with FA, aged sea salt (g) and several minerals were also observed.
Atmosphere 2020, 11, x FOR PEER REVIEW 14 of 19 SiO2, Mg, and Fe) have been observed, indicating that they may have originated from the volcanic plume (Figure 5a,d). The anthropogenic factor has influenced the sampling (Figure 5f,j) and consequently the experimental results both in terms of chemical and microscopy analysis interpretation (i.e., factor and SEM analysis). In agreement with FA, aged sea salt (g) and several minerals were also observed.

Stepwise Linear Discriminant (S-LDA) Analysis
To seek which variables best discriminate among samples coming from volcanoes (Vo) and solfatara (So), we performed stepwise linear discriminant analysis (S-LDA). The forward S-LDA (F to enter = 5.00 and F to remove = 4.00) retained six elements, namely Zr, As, Hf, Bi, Ca, and Zn, that allowed the construction of a model with satisfactory statistical parameters (Wilks' λ = 0.09560; F (6, 13) = 20.497; p < 0.00001). The low Wilks' λ values suggest that the model is highly discriminating, while the low p-level value shows that the classification, which resulted correct for all the samples, occurred with high reliability. Based on the partial Wilks' λ, Zn (0.38) is the independent variable that contributes most to the discrimination, followed by Ca (0.50). The scatterplot of canonical scores ( Figure S8) on the discriminant function shows the separation between the groups of the volcanic samples (Vo), with negative root values ranging from −4.38 to −1.56 and the solfatara samples (So) with positive root values in the range between 1.39 and 4.69. We found the predominant presence of Zn in solfatara samples, in accordance to Dekov et al. [32], which detected this element in the aerosol of sulfur-rich sites such as Panarea or other fumarole areas because of its presence as sphalerite (ZnS).

Stepwise Linear Discriminant (S-LDA) Analysis
To seek which variables best discriminate among samples coming from volcanoes (Vo) and solfatara (So), we performed stepwise linear discriminant analysis (S-LDA). The forward S-LDA (F to enter = 5.00 and F to remove = 4.00) retained six elements, namely Zr, As, Hf, Bi, Ca, and Zn, that allowed the construction of a model with satisfactory statistical parameters (Wilks' λ = 0.09560; F (6, 13) = 20.497; p < 0.00001). The low Wilks' λ values suggest that the model is highly discriminating, while the low p-level value shows that the classification, which resulted correct for all the samples, occurred with high reliability. Based on the partial Wilks' λ, Zn (0.38) is the independent variable that contributes most to the discrimination, followed by Ca (0.50). The scatterplot of canonical scores ( Figure S8) on the discriminant function shows the separation between the groups of the volcanic samples (Vo), with negative root values ranging from −4.38 to −1.56 and the solfatara samples (So) with positive root values in the range between 1.39 and 4.69. We found the predominant presence of Zn in solfatara samples, in accordance to Dekov et al. [32], which detected this element in the aerosol of sulfur-rich sites such as Panarea or other fumarole areas because of its presence as sphalerite (ZnS). On the contrary, the elements Hf and Zr contribute most to characterize volcanic samples. The presence of Hf and Zr agrees with the study of metal emissions from the Stromboli volcano by Allard et al. [51]. Indeed, Hf is one of the elements degassed from magmas as halides, sulfates, sulfides, and/or metals, and is typically found in the aerosol phase of the airborne plume, while the presence of Zr in volcanic samples can be due to its transport through volcanic ash, supporting the evidence previously obtained from factor analysis.

Conclusions
In this study, we investigated the impact of volcanic and fumarole sources on PM level and its chemical characteristics in the central Mediterranean Sea. Aerosol samples in marine atmosphere were collected during the research cruise campaign Med-Oceanor 2017 off the coast of the main active volcanic and fumarolic areas in the Mediterranean Sea, which included Mount Etna, the Eolian volcanic arc, Marsili Seamount, and the Phlegraean Fields. The mass concentration of collected PM 10 resulted, for each site, lower than the daily limit set by the European legislation. The PM 10 and PM 2.5 samples were characterized by quantification of major and trace elements, ionic species, as well as elemental and organic carbon content. These chemical data were illustrated by descriptive statistics, histograms, and then analyzed by the computation of mass closure, enrichment factors, factor analysis, and correlation analysis with the support of the daily SO 2 values. The PM 10 mainly consisted of carbonaceous material and dust, which together made up over 50% of the mass, while carbonaceous fraction and SIA, which reached 60% of the mass, prevalently composed the PM 2.5 . The trace element was the less abundant contributor with an average amount of about 1% in both PM fractions. From the EFs, we found that the particulate from the volcanic areas resulted enriched by Tl, Cu, La, Se, As, and alkali metals, whereas the chemical characterization of PM in fumarole areas showed a good correlation between the sulfur dioxide gas, emitted from fumaroles and sulfate ion, coming from the SO 2 reaction with water and condensation on particles. FA allowed for the identification of six potential sources that contribute to the PM composition, some of which have a natural origin such as volcanic, marine, and geogenic emissions, while other sources are anthropogenic, such as combustion processes, vehicular traffic, and shipping emissions. The SEM/EDX analysis on the sampled PM confirmed the influence of these sources identified with the FA. Stepwise linear discriminant analysis (S-LDA), which was used to seek for a model capable to discriminate between the sample of volcanic and fumarolic origin, pointed out the elements that discriminate the most between these two classes, namely Hf and Zr for volcano samples and Zn and Ca for solfatara samples.
Further investigation on this topic could be addressed in future research for a more thorough assessment of the impact of the volcanic and fumarolic emission on the air quality of the Mediterranean Sea. In this regard, the use of a high-volume sampler would allow sampling the PM with sufficient sensitivity even with a few hours of collection, thus improving the temporal profile of the impact. Similarly, further advancement of the study would be given by the use of cascade impactors, whose data could provide new insight to strengthen the characterization of the marine atmosphere in the neighborhood of volcanoes and fumaroles.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4433/11/2/149/s1, Figure S1: Hourly mean meteorological data; Temperature and Pressure (up) and wind speed (down), Figure  S2: Wind roses for sampling points, Figure S3: NCEP/NCAR reanalysis maps (mean and anomaly) of the period from the start of the campaign to 31 August (up) and 2 September (down), Figure S4: Acireale maps NAAPS, air mass back trajectory and BSC-Dream, Figure S5: Time series plot of SO 2 for each of the sampling sites; Y-axis: concentration (ppb), X-axis: sampling time (hours), Figure S6: Enrichment Factors for Etna (Acicastello sample) and Stromboli (Stromboli (3) sample). Dashed lines highlight thresholds for EFs = 10 and EFs = 100, Figure  S7: Correlation between SO 2 and nss-SO 4 2− (as µg m −3 ) in PM 10 and PM 2.5 . The plots (a) and (b) comprise all the sampling points, whereas the Vulcano (2) sample is excluded in (c) and (d) plots the plots., Figure S8: LDA scatterplot for aerosol samples collected in volcanic (Vo) or solfatara (So) sites, Table S1: Percentage chemical composition of PM 10 and PM 2.5 in each sampling site.
Funding: This research was funded by the European Commission-H2020, the ERA-PLANET programme (www.era-planet.eu; contract no. 689443) within the IGOSP project (www.igosp.eu).