A Year-Round Measurement of Water-Soluble Trace and Rare Earth Elements in Arctic Aerosol: Possible Inorganic Tracers of Speciﬁc Events

: This study presents the year-round variability of the water-soluble fraction of trace elements (wsTE) and rare earth elements (wsREE) among size segregated airborne particulate matter samples collected at Ny-Ålesund in the Svalbard Archipelago from 26 February 2018 to 26 February 2019. Six different aerosol dimensional fractions were collected using a multi-stage Andersen impactor to better understand local and global circulation with the aim of disentangling the source of inorganic tracers from speciﬁc natural or anthropogenic sources. The wsTE and wsREE content, especially in the ﬁnest fractions in remote areas, is primarily related to long-range transport and it gives valuable information on (1) the global circulation, (2) the natural sources and (3) the contribution of human activities to aerosol composition. A Factor Analysis was applied to the dataset, including levoglucosan and methanesulfonic acid (MSA), to assess the possibility of using certain inorganic tracers as indicators of speciﬁc transport events or circulation regimes. We also investigate back-trajectories to determine potential source areas.


Introduction
Polar regions are widely acknowledged as important areas for studying the effect of climate change, the impact of anomaly natural events and human activities on the climate itself. Warming in the Arctic has occurred much more rapidly than the global average and in particular the Ny-Ålesund archipelago has experienced a significant rise in temperature over the last two decades [1]; this phenomenon is known as Arctic amplification [2].
Particulate matter can influence the chemical and radiative properties of the atmosphere. These particles can transport material through the atmosphere, and they can act as cloud condensation nuclei (CCN), affecting the optical properties of the atmosphere [3]. The transport of particulate matter (PM) from low and medium latitudes towards the polar regions may dramatically impact the quality of these fragile areas. To understand the role that aerosol could have in variation in the budget of solar radiation, reduction of visibility and possible pollution of the ecosystem [4,5], it is useful to know the geochemical characteristics of the PM that reach Arctic areas. In addition, the aerosol composition could significantly affect the snow pack chemical composition and the load of specific elements [6]. Another important property of the aerosols is their size distribution, giving valuable information on the sources of these aerosols and on the atmospheric processes modifying their properties during atmospheric transportation.
The investigation of the chemical composition of Arctic aerosol is essential to discriminate between local sources or long-range atmospheric transport processes in order to define Forty-six aerosol samples were collected from 26 February 2018 to 26 February 2019 in the Gruvebadet atmospheric laboratory (Svalbard Islands, 78 • 55 03" N, 11 • 53 39" E, 50 m a.s.l.), close to Ny-Ålesund and 1.4 km far from the Kongsfjorden.
As described in   [17], a multi-stage Andersen impactor (TE-6000 series, Tisch Environmental Inc., Cleves, OH, USA) was used to collect aerosol samples on six pre-combusted (4 h at 400 • C in a muffle furnace) quartz fibre filters. The sampler accumulated particles with cut-off diameters of 10.0 µm, 7.2 µm, 3.0 µm, 1.5 µm and 0.95 µm on slotted quartz fibre filters (QFF) and <0.49 µm on the backup filters (also reported in this paper as "stages" and "dimensional fractions"). The frequency of sampling was 6 days in spring and summer and 10 days in autumn and winter, following the results and the protocols of previous campaigns [11,18]. The mean total air volume was 9000 and 16,000 m 3 , respectively. The details about the collection period are reported in Table S1. Field blanks were obtained using filters installed on the sampler for 5 min with the air pump switched off. Samples and blanks were stored at −20 • C until analysis. Referring to the division in seasons, we adopted the following intervals: end-Winter 2018, from 26

Instrumental Analysis
The description of the sample treatment has been reported in five previous studies [10,13,14,17,18]. Briefly, the slotted and backup QFF were cut in half using stainlesssteel scissors that were previously washed with methanol. Each filter was extracted twice in an ultra-sonic bath at 10 • C for 15 min using ultrapure water. The slotted filters and the back-up filters were extracted in a first step with 9 and 25 mL of ultrapure water followed by a second step with 1 and 5 mL of ultrapure water, respectively. The extracts were unified, filtered, using a 0.45 µm polytetrafluoroethylene (PTFE) filter (Whatman, Maidstone, Kent, UK), and immediately frozen until the analysis. All samples and field blank filters were handled inside an ISO5 clean room under a laminar flow bench (class 100) to avoid any contamination from laboratory air particles.
All 300 samples (46 collected samples per 6 fractions plus 4 blanks per 6 fractions) were analysed at the Institute of Polar Sciences-CNR/University of Ca' Foscari laboratories in Venice to determine TE and REE by means of an iCAP-RQ ICP-MS (Thermo Scientific™, Bremen, Germany) equipped with an autosampler. Analyses were performed on acidified samples (2% v/v of ultrapure nitric acid, Romil™ UPA, Cambridge, UK) and quantification was done using the external calibration method. The multi-elemental stock solutions used for TE and REE were IMS102 and IMS101 by Ultra Scientific (Santa Clara, CA, USA), respectively. Standards at concentrations of 0, 10, 20, 50, 100, 200 and 500 ng·L −1 and 1, 5, 10, 50 and 100 µg L −1 for TE and 0, 2, 5, 10, 20, 50, 100 and 200 ng·L −1 for REE were prepared gravimetrically, adding the appropriate quantity of a multi-elemental stock solution (IMS102 and IMS101 by Ultra Scientific, respectively). During each session of analysis, also six blank solutions were repeatedly analysed to evaluate the detection limits, calculated as three times the standard deviation of six blanks.
Levoglucosan and methanesulfonate (MSA) were used as specific markers for biomass burning and phytoplankton blooms, in comparison with the concentration of wsTE and wsREE to confirm potential emission sources. The determination and quantification of these compounds were already described elsewhere [25]. We used an ion chromatograph (Thermo Scientific™ Dionex™ ICS-5000, Waltham, MA, USA) coupled with a single quadrupole mass spectrometer (MSQ Plus™, Thermo Scientific™, Bremen, Germany). Briefly, the determination of levoglucosan was performed using a CarboPac MA1™ analytical column (Thermo Scientific, 2 × 250 mm) equipped with an AminoTrap column (2 × 50 mm). The NaOH eluent gradient at a flow rate of 0.25 mL min −1 was as follows: 20 mM (0-23 min), 100 mM (23-43 min) and 20 mM (43-53 min). The injection volume = 50 µL. Quantification was performed using labelled 13 C 6 -levoglucosan as the internal standard. MSA determination was achieved using an anionic exchange column (Dionex Ion Pac AS11 2 × 250 mm) and a guard column (Dionex Ion Pac AG11 2 × 50 mm). The NaOH gradient at a flow rate of 0.25 mL min −1 was as follows: 0-3.5 min gradient from 0.5 to 5 mM; 3.5-5 min gradient from 5 to 10 mM; 5-25 min gradient from 10 to 38 mM; 25-30 min, column cleaning with 38 mM; 30-35 min; and equilibration at 0.5 mM. The injection volume = 100 µL. Quantification was achieved using an external calibration curve. All specific details about the instrumental and methodological parameters are reported in previous papers [10,26].

Enrichment Factor
The use of trace element enrichment factors (EFs) relative to the Earth's crust concentrations is useful in highlighting the contribution of non-natural sources to elemental levels. EF is calculated as follows: where EF i is the enrichment factor of element i, j is a reference element of crustal origin, (i/j) atmosphere is the ratio of element i to element j in the atmosphere and (i/j) upper crust is the ratio of element i to element j in the upper crust. The composition of the upper crust was described in Wedepohl et al. (1995) [27] while, following the suggestion of Gao et al. (1992) [28], we used Al as the reference element in the EF calculations. An EF near 1 indicates a crustal derived element (geogenic), while significantly higher EF values indicate a non-crustal-derived element (non-geogenic): an EF between 10 and 100 indicates a moderate enrichment while an EF above 100 indicates a prevalent non-crustal origin. Calculating the EFs allows us to make a starting estimate of the potential source of various elements but the indications obtained must be confirmed by other evaluations.
To estimate the influence of marine source we also calculated the marine EFs (MEFs) relative to sea water concentrations, using the abundance in sea water reported by Nozaki et al. (2010) [29]. MEFs were calculated as Mg was used as a reference element (j) of marine origin. Analogous to the EF, an MEF near 1 indicates a marine-derived element, while higher values indicate moderate (10-100) to high (>100) enrichment relative to sea water.

Chemometric Approach and Back Trajectories Analysis
Factor Analysis (FA) is a chemometric approach that allows us to reduce the size of our data matrix without losing the information it contains. FA with varimax rotation was performed on the auto-scaled data matrix using Statistica 10.0 (StatSoft, Inc., 2011, Hamburg, Germany) [30]. The varimax rotation is the most commonly rotation strategy used with the aim to better clarify the pattern of loading, highlighting the variables that mark the extracted factors. The concentrations below the method detection limit (MDL) were substituted with a value of 1 2 MDL in the FA. Data are presented as a plot of two factors; the X,Y values for each sample represent the importance of the variables that characterize the factor, for that sample itself, giving an estimation if the sample is related to a certain factor or not. High values of a factor indicate its importance for the sample.
To understand the general transport pattern of air masses recorded at the sampling site, backward air trajectories ending at the Gruvebadet aerosol laboratory were computed using Hybrid Single-Particle Lagrangian Integrated Trajectories (HYSPLIT) transport and dispersion models [31,32]. Back trajectories were calculated using a vertical velocity model with an endpoint at 500 m above ground level and starting a trajectory every 6 h. Knowing that airborne particles in the troposphere have typical lifetimes of 3-10 days on average, also considering dry and wet deposition, we chose to compute 7-days back trajectories [33]. The meteorological data used for computing all the backward trajectories were the NCEP/GDAS (Global Data Assimilation System) at a 1 × 1 • resolution. Based on the backward particle release simulation, the cluster aggregation was displayed for each month ( Figure S2).

Results
All six dimensional fractions of each sample collected during the February 2018 to February 2019 campaign were analysed to determine the content of the water-soluble REE (La-Lu) and water-soluble TE (Li, Be, Mg, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, As, Y, Rb, Sr, Ag, Cd, Cs, Ba, Tl, Pb, Bi, U, Th) content ( Figure 1A

Discussion
The year-round study of chemical composition of PM particle-size distribution, gathered from February 2018 to February 2019 at Ny-Ålesund, a key region to understand air mass circulation, provides new information about the origin of the PM and the importance of their different dimensional size in discriminating possible source areas and specific events. The back-trajectory evaluation helps us to interpret the results ( Figure S2) because this method is widely used in recognizing the routes of pollutants [36]. In fact, the back trajectories distinctly indicate that the source areas were strictly related to seasonality. During summer, transport from the mid-latitudes decreased due to the contraction of the polar vortex, while in winter the long-range atmospheric transport covers a more extended area [23].
Different approaches were applied to interpret the results and to recognize the PM sources using wsTE and wsREE as tracers. Firstly, we considered the EFs calculated for all analysed elements related both to the upper crust and seawater composition. Then, the geogenic or non-geogenic origin of the samples was evaluated using ternary diagrams applied to REEs. Finally, a chemometric approach was used to reduce the number of variables and to link the inorganic and organic components of samples. In fact, Factor Analysis (FA) was performed, including in our dataset also levoglucosan and MSA. Levoglucosan is recognized as a tracer of biomass burning [37] while MSA is related to algal bloom [38], suggesting primary marine production as the source.

Enrichment Factors
The Figure 3 shows the EFs calculated for the measured TE from end-winter 2018 to winter 2019 as an average of the seasonal trend related to the upper crust composition [27]. Zn, Cd, Cu and As are the most enriched elements with values above 100, indicating a non-geogenic origin. Ag, Pb, Ni, Cr and Bi have values around 100, indicating that they are only weakly related to the crust. Th, Fe, Ga, REE, Rb, Ca, Be, Cs, Y and U, with values quite always below 10, show a geogenic origin. K, V, Li, Mn, Ba, Co, Tl, Sr and Mg have values that indicated a probably geogenic source; but, especially in Summer and Fall, these elements have values around 100 or higher, indicating probable mixed sources. Values strongly above 100 from Bi, Cr, Ni, Pb, Ag, As, Cu, Cd and Zn suggest that anthropogenic sources are responsible for the concentrations of these elements. These values are according to the previous results reported by Turetta et al. (2016) [18], where the investigation of the particle-size distribution of these species was performed.
It has to be noted that, during summer, the EF values of all stages are similar for each wsTE due to the presence of local sources as the main input. In contrast, during winters and spring, different values of EFs are found in the different stages because long-range transport plays an important role in the PM composition. Although size distribution is not directly correlated to EFs, the information about the EF values in the different stages can support the seasonal variation of input in the Arctic atmosphere. The most abundant wsREE was Ce (2 ± 1 pg m −3 ), followed by La (0.6 ± 0.5 pg m −3 ) and Nd (0.6 ± 0.4 pg m −3 ), following the normal order of abundance of most crustal materials (Ce > La > Nd > Sm > Yb). The wsREE were equally distributed between coarse and fine fractions ( Figure 1B). A slight shift toward the fine fraction in winter was detected ( Figure 2D).
The concentrations found in this campaign were comparable with the PM 10 values determined during spring and summer in some previous studies at Gruvebadet laboratory [2,9,34] and in the Chinese Arctic "Yellow River Station" at Ny-Ålesund [35].
The shift reported toward the fine fraction in spring and winter is frequently linked with an increase in concentration of wsTE and wsREE ( Figure 2). This suggest that the long-range atmospheric transport (LRAT) seems to be the dominant factor of the presence of metals in the Arctic atmosphere.

Discussion
The year-round study of chemical composition of PM particle-size distribution, gathered from February 2018 to February 2019 at Ny-Ålesund, a key region to understand air mass circulation, provides new information about the origin of the PM and the importance of their different dimensional size in discriminating possible source areas and specific events. The back-trajectory evaluation helps us to interpret the results ( Figure S2) because this method is widely used in recognizing the routes of pollutants [36]. In fact, the back trajectories distinctly indicate that the source areas were strictly related to seasonality. During summer, transport from the mid-latitudes decreased due to the contraction of the polar vortex, while in winter the long-range atmospheric transport covers a more extended area [23].
Different approaches were applied to interpret the results and to recognize the PM sources using wsTE and wsREE as tracers. Firstly, we considered the EFs calculated for all analysed elements related both to the upper crust and seawater composition. Then, the geogenic or non-geogenic origin of the samples was evaluated using ternary diagrams applied to REEs. Finally, a chemometric approach was used to reduce the number of variables and to link the inorganic and organic components of samples. In fact, Factor Analysis (FA) was performed, including in our dataset also levoglucosan and MSA. Levoglucosan is recognized as a tracer of biomass burning [37] while MSA is related to algal bloom [38], suggesting primary marine production as the source.

Enrichment Factors
The Figure 3 shows the EFs calculated for the measured TE from end-winter 2018 to winter 2019 as an average of the seasonal trend related to the upper crust composition [27]. Zn, Cd, Cu and As are the most enriched elements with values above 100, indicating a non-geogenic origin. Ag, Pb, Ni, Cr and Bi have values around 100, indicating that they are only weakly related to the crust. Th, Fe, Ga, REE, Rb, Ca, Be, Cs, Y and U, with values quite always below 10, show a geogenic origin. K, V, Li, Mn, Ba, Co, Tl, Sr and Mg have values that indicated a probably geogenic source; but, especially in Summer and Fall, these elements have values around 100 or higher, indicating probable mixed sources. Values strongly above 100 from Bi, Cr, Ni, Pb, Ag, As, Cu, Cd and Zn suggest that anthropogenic sources are responsible for the concentrations of these elements. These values are according to the previous results reported by Turetta et al. (2016) [18], where the investigation of the particle-size distribution of these species was performed. The MEFs ( Figure S3A,B) were calculated for all elements analysed as an average of the seasonal trend related to the seawater mean composition [29]. Ca and K have values close to 1, confirming their marine origin. Li, with values between 1 and 10, Rb, with values between 1 and 30, and U, with values between 2 and 50, seem to have a prevailing marine input but other sources must be considered. Figure S3B emphasizes the differences between seasons, highlighting the importance of the dimensional size of the elements of mixed origin. In spring, a higher variability is observed for those sizes with a possible long-range origin, while in summer the finer fractions do not seem to be relevant. Endwinter 2018 and Winter 2019 seems to be quite different, but it is to note that both seasons show the same trend of considered marine origin elements and a relevant contribution of long-range transport for those elements of mixed origin. All other elements, with values higher than 100, appear to have sources other than the marine one. From the comparison between the two calculated series of EFs, related to the Upper Crust and Sea Water composition, respectively, we can hypothesize that the mixed source detected by the EFs for Ca, K, Sr, Li, Rb and U could be related to marine inputs due to the fact that their MEF values range between 1 and 50.

Ternary Diagrams
With the aim to understand the potential sources of the analysed elements, we evaluate our data using ternary diagrams applied to REEs to evaluate the geogenic or nongeogenic origin of the samples. The diagram in Figure 4a compares La-Ce-Sm, whose values were adjusted to put the crustal concentration in the center of the ternary plot. The coarse fraction (  It has to be noted that, during summer, the EF values of all stages are similar for each wsTE due to the presence of local sources as the main input. In contrast, during winters and spring, different values of EFs are found in the different stages because long-range transport plays an important role in the PM composition. Although size distribution is not directly correlated to EFs, the information about the EF values in the different stages can support the seasonal variation of input in the Arctic atmosphere. The MEFs ( Figure S3A,B) were calculated for all elements analysed as an average of the seasonal trend related to the seawater mean composition [29]. Ca and K have values close to 1, confirming their marine origin. Li, with values between 1 and 10, Rb, with values between 1 and 30, and U, with values between 2 and 50, seem to have a prevailing marine input but other sources must be considered. Figure S3B emphasizes the differences between seasons, highlighting the importance of the dimensional size of the elements of mixed origin. In spring, a higher variability is observed for those sizes with a possible long-range origin, while in summer the finer fractions do not seem to be relevant. End-winter 2018 and Winter 2019 seems to be quite different, but it is to note that both seasons show the same trend of considered marine origin elements and a relevant contribution of long-range transport for those elements of mixed origin. All other elements, with values higher than 100, appear to have sources other than the marine one. From the comparison between the two calculated series of EFs, related to the Upper Crust and Sea Water composition, respectively, we can hypothesize that the mixed source detected by the EFs for Ca, K, Sr, Li, Rb and U could be related to marine inputs due to the fact that their MEF values range between 1 and 50.

Ternary Diagrams
With the aim to understand the potential sources of the analysed elements, we evaluate our data using ternary diagrams applied to REEs to evaluate the geogenic or non-geogenic origin of the samples. The diagram in Figure 4a compares La-Ce-Sm, whose values were adjusted to put the crustal concentration in the center of the ternary plot. The coarse fraction (10-3 µm) is prevailingly distributed around the crustal position while the medium fraction   In Figure 4b, we report the ternary diagram V-Ce-La to highlight the anthropogenic sources that can characterize our samples. It is known from the literature that V can help distinguish different potential anthropogenic sources when compared with some REEs. In particular, V may differentiate the contributions of emissions from ships or other oil combustion processes from those of refineries or motor vehicles [19,39,40].
Our samples are principally in the right corner, the V apex, indicating as prevailing In Figure 4b, we report the ternary diagram V-Ce-La to highlight the anthropogenic sources that can characterize our samples. It is known from the literature that V can help distinguish different potential anthropogenic sources when compared with some REEs. In particular, V may differentiate the contributions of emissions from ships or other oil combustion processes from those of refineries or motor vehicles [19,39,40].
Our samples are principally in the right corner, the V apex, indicating as prevailing source the oil combustion, but it is to note that the position of a relevant part of samples shifts from the V apex toward the crustal composition, especially samples of medium and fine dimensions. This shift could indicate a mixed source for these samples: oil combustion and weathering of the surface are the more probable, but, due to a displacement toward high values of Ce of some samples, we cannot exclude motor vehicles as a possible source.
V and Ni can be also used as typical markers for fossil fuel combustion from ships. Viana et al. (2008) [40] established V and Ni as valid tracers of shipping emissions, with a typical V/Ni ratio of 2.5-3.0. In order to better characterize this factor, the V/Ni ratio was calculated: the ratio was always below 1 during the whole sampling period. Furthermore, V and Ni showed only a slight concentration increase during the 20-25 June and 14-19 July samples (215 and 186 pg m −3 , respectively). The wind roses in the first period ( Figure S4) suggest a source of the particles from the fjord due to ship emissions. On the contrary, from 13 to 19 June, no favourable winds were detected. Considering that the marine traffic in the Kongsfjorden is generally high in June and July (3 to 10 landing ships per day in the Ny-Ålesund harbour), no valuable relation between local ship emission and V and Ni was reported in 2018.

Chemometric Evaluation
With the aim to better evaluate all samples and all variables together, and considering the dimension of our data-matrix, 276 samples (46 samples per 6 dimensional fractions) per 41 elements plus 2 organic compounds used to compare inorganic results with the information obtained from organic fraction, we have performed a statistical evaluation of all the data. In particular, Factor Analysis was performed, extracting four factors that account together for 75.6% of the explained variance.
In Figure 5 and Figure S5, we present the results of the FA as graphs of all the samples subdivided for the dimensional fractions in order to highlight the differences between the dimensions with respect to the possible sources. For clarity, we have plotted in three separate graphs the results obtained from the FA for each couple of factors on a dimensional basis. The particle-size distribution influences the results of the FA because the difference between wsTEs and wsREEs is clearly shown in Table 1, where the wsREEs are loaded in Factor 1.
The four extracted factors (Table 1) account together for 75.6% of the explained variance. The first factor brings together all REEs and Be, Mn, Co and U; the second levoglucosan, Be, V, Fe, Ga, As, Rb, Cd, Cs, Tl, Pb and Bi; the third Li, Mg, K, Ca and Sr; and the fourth MSA, Cr, Ni, Cu and Zn. Considering these clusters of elements, we can hypothesize Factor 1 (F1) as related to crustal source. In the score plot ( Figure 5A), on the basis of F1, the separation of the samples collected in spring, from March to April in coarser fractions and from March to June for finer ones, is evident compared to the others. This separation could be related to meteorological factors considering that springtime is very often characterized by strong winds. Factor 2 (F2) seems to be related to combustion, both biomass burning, as suggested from the presence in this factor of levoglucosan, and oil combustion, as suggested from V. In the finer fractions, this factor separates winter samples from all the others, while no-separation is evident in coarser fractions. The seasonal variation, showing elevated concentrations in winter (Tables S2 and S3), may be due to the transport from the midlatitudes to the Arctic intensifying during the winter and spring [22,23,41], as confirmed also by the back trajectories ( Figure S2) that well overlay the Northern Russian areas. Factor 3 (F3), mainly associated with elements related to marine origin, seems partially separate, both in the finer and coarser fractions, in the spring-summer samples compared to the others, with some exceptions, such as the sample from 26 February 2018 that show a high value for F3 ( Figure 5B). The importance of the marine source in the summer samples is related to the lack of sea-ice cover and to the input from the Arctic Ocean in those periods, as highlighted from back trajectory evaluations ( Figure S2). The medium fractions seem to be not affected by this factor. Factor 4 (F4), related to marine primary production on the basis of the presence of MSA, seems to have a strong effect only on the finest fraction (<0.49 µm), separating the summer and partially spring samples from the others ( Figure S5). The presence of micronutrients, such as Cu and Zn, support this interpretation of F4. Table 1. Factor loadings matrix after varimax rotation was applied. In bold are the loadings that are prevailing in a single factor, and in italic the loadings that are distributed in more than one factor.    (Table 1).

Conclusions
On the basis of all the presented data, we can conclude that the year-round sampling campaign lets us evaluate the evolution of PM through the different seasons in an Arctic area and gives an indication of the possible source areas for PM that reach the Arctic. In particular, the one-year sampling allowed us to recognize the different sources that contribute to the composition of wsTE and wsREE in particulate matter: crust source, combustion, marine inputs and marine primary productivity. On the other hand, the six different dimensional sizes of the samples separate, in the same season, the short-range from long-range contributions, highlighting the crustal source and different anthropogenic sources; for example, ternary diagrams help us to recognize both vehicular traffic and oil combustion inputs. The information obtained considering EFs, ternary diagrams and Factor Analysis, with the support of back-trajectory analyses, confirms the ability of TE and REE in tracing PM, both in terms of source area and seasonality when different dimensional sizes of the particles are considered.
From the comparison with previous research, in the same area, we can state that one year is probably not enough to fully understand the behaviour of TEs and REEs. However, this comparison makes it possible to hypothesize a similar behaviour for the same season in different years, obviously considering the ever-possible local variability and the occurrence of specific events, such as a volcanic eruption or wildfires.
The importance of recognizing sources of PM is highlighted here: the pollutants and, more generally, PM reaching the Arctic regions from the mid-latitudes may result in potential negative feedback regarding climate change. To understand if and how the PM is related to short-or long-range transport, and to different dimensional sizes, can help us in recognizing the role that it plays in the Earth's climate system.