Estimation of Elements’ Concentration in Air in Kosovo through Mosses as Biomonitors

Elements atmospheric deposition was studied by mosses method in Kosovo. Mosses are a very useful tool for atmospheric deposition of elements monitoring, owing to their physiological and morphological characteristics. Moss samples were collected from 45 locations, they were cleaned from leaves, twigs and other materials, they were grinded, digested in a microwave system and analyzed by spectroscopy Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-AES) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). A total of 25 elements concentration in collected moss samples was determined: Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Sb, Sr, Ti, Tl, V and Zn. Extremely high concentrations of Pb and Zn were found in moss samples from the areas of Zveçan and Stanterg where the ore processing facilities are located, as well as the tailings landfill in Mitrovica which is nearby in the south, which reflects eminent air pollution. High concentrations of As and Cd were also measured in the same sampling points. The obtained results were processed by multivariate statistical analysis/factor analysis, which resulted in four factors, one anthropogenic and three geogenic or mixed geogenic-anthropogenic. The factor analysis revealed that the area with the highest atmospheric deposition of potentially toxic elements is the basin of Kosovo, from the north in Leposaviç to the south in Hani i Elezit because of the industrial plants laying in this line.


Introduction
Air pollution in general is a very eminent issue globally and particularly for developing countries which do not dispose with appropriate technology to reduce gasses, dust emissions and other waste products into the atmosphere. Among chemical pollutants, heavy metals are of great concern because they are known to inflict several health disorders in humans. For illustration a few health damaging effects of some of the heavy metals are mentioned as follows. Lead is largely known for its toxicity effects in humans, such as nervous system damage in children [1], children recovered from lead encephalopathy can suffer from mental retardation, optic neuropathy and epilepsy [2], adults may develop hypertension [3]. Mercury has neurotoxic effects and even deaths are reported as a result of mercury intoxication [4]. Cadmium shows adverse effects in respiratory system, bone, cardiovascular system, liver and other health complications [5]. Exposure to arsenic can result in several health disorders including cancer [6].
Heavy metals, not only are harmful to the environment, but they also are refractory so they do not ever destroy like some of the organic pollutants do. All elements occur naturally but human activities cause them to accumulate in different areas, where as concentrated

Sampling
The collection of moss samples was performed during the dry season of summer, August to September in 2019. The homogeneity of sampling points was to the maximum the presence of mosses allowed, as can be seen in Figure 2. The procedure of sampling

Sampling
The collection of moss samples was performed during the dry season of summer, August to September in 2019. The homogeneity of sampling points was to the maximum the presence of mosses allowed, as can be seen in Figure 2. The procedure of sampling was in accordance with the Monitoring Manual of "International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops, 2020". There were collected three types of mosses (Figure 3), Homalothecium sericeum (Hedw.) Schimp. 1851, Hypnum cupressiforme Hedwig, 1801 and Pseudoscleropodium purum (Hedw.) M.Fleisch. Three to ten (most of them over five) samples were collected in each sampling site, in an area of 50 × 50 m 2 and they were put in paper bags of 1 L which were stored for 5-7 days at ambient temperature (20-25 • C) prior to cleaning. There was no mixing of the types of mosses, in each sampling site only one type of mosses was collected. Samples were taken at least 100 m away from small roads and over 300 m away from main roads, villages and industries. All samples were collected from ground and rarely from rotten branches and trees in open fields with low shrubs or no shrubs at all. In case of forest covered areas mosses were taken from forest gaps of at least 10 m diameter and over 3 m from the nearest tree canopy drip.
Atmosphere 2021, 12, x FOR PEER REVIEW 4 of 17 was in accordance with the Monitoring Manual of "International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops, 2020". There were collected three types of mosses (Figure 3), Homalothecium sericeum (Hedw.) Schimp. 1851, Hypnum cupressiforme Hedwig, 1801 and Pseudoscleropodium purum (Hedw.) M.Fleisch. Three to ten (most of them over five) samples were collected in each sampling site, in an area of 50 × 50 m² and they were put in paper bags of 1 L which were stored for 5-7 days at ambient temperature (20-25 °C) prior to cleaning. There was no mixing of the types of mosses, in each sampling site only one type of mosses was collected. Samples were taken at least 100 m away from small roads and over 300 m away from main roads, villages and industries. All samples were collected from ground and rarely from rotten branches and trees in open fields with low shrubs or no shrubs at all. In case of forest covered areas mosses were taken from forest gaps of at least 10 m diameter and over 3 m from the nearest tree canopy drip.

Digestion and Chemical Analysis
After cleaning mosses form other herbs, leaves, twigs and other materials, they were grinded by hand using gloves and dried at 40 • C for 48 h. Samples were digested in a microwave system by the wet digestion method. 0.5 g of only one type of mosses sample were put in the teflone tube (there was no mixing of types of mosses for sampling points, although three types of them were collected over the whole study area), where 7 mL of nitric acid and 2.5 mL of hydrogen peroxide were added. The mixture was left to react for 10 min in order to give some time reactions to take place and some gases to release, so that less pressure would develop in the teflon tube during microwave digestion. Then,

Digestion and Chemical Analysis
After cleaning mosses form other herbs, leaves, twigs and other materials, they w grinded by hand using gloves and dried at 40 °C for 48 h. Samples were digested microwave system by the wet digestion method. 0.5 g of only one type of mosses sam were put in the teflone tube (there was no mixing of types of mosses for sampling po although three types of them were collected over the whole study area), where 7 m nitric acid and 2.5 mL of hydrogen peroxide were added. The mixture was left to reac 10 min in order to give some time reactions to take place and some gases to release, so less pressure would develop in the teflon tube during microwave digestion. Then, Te tubes were put in the microwave system (Analytic Jena TOPwave) with the follow program: 5 min up to 170 °C, hold time of 10 min at 170 °C, 1 min up to 200 °C, hold of 15 min at 200 °C, 1 min down to 50 °C, hold time of 23 min at 50 °C. After digestion obtained solutions were filtered, poured in 25 mL plastic flasks and the rest of the v was filled with redistilled water and sent for elements chemical analysis.

Instrumentation
For this study reagents with analytical grade or better were used: nitric acid, pure (Merck, Germany) and hydrogen peroxide, p.a. (Merck, Germany). Redistilled w was used for the preparation of all solutions. Standard solutions of metals were prep by dilution of 1000 mg/L solutions (11355-ICP multi Element Standard).
By the application of inductively coupled plasma-atomic emission spectrom (ICP-AES) and inductively coupled plasma-mass spectrometry (ICP-MS) the follow 25 elements were determined: Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, K, Li, Mg, Mn, Na, Ni, P, Pb, Sb, Sr, Ti, Tl, V and Zn. ICP-AES (Varian, model 715ES) was used fo analysis of Al, Ba Ca, Co, Cr, Cu, Fe, K, Mn, Na, P, Sr, V and Zn. The QC/QA of the app technique was performed by standard addition method and it was found that the reco for the investigated elements ranges for ICP-AES was 98.5-101.2%. Quality control was ensured by standard moss reference materials M2 and M3, which are prepared for the

Instrumentation
For this study reagents with analytical grade or better were used: nitric acid, trace pure (Merck, Germany) and hydrogen peroxide, p.a. (Merck, Germany). Redistilled water was used for the preparation of all solutions. Standard solutions of metals were prepared by dilution of 1000 mg/L solutions (11355-ICP multi Element Standard).
By the application of inductively coupled plasma-atomic emission spectrometry (ICP-AES) and inductively coupled plasma-mass spectrometry (ICP-MS) the following 25 elements were determined: Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Sb, Sr, Ti, Tl, V and Zn. ICP-AES (Varian, model 715ES) was used for the analysis of Al, Ba Ca, Co, Cr, Cu, Fe, K, Mn, Na, P, Sr, V and Zn. The QC/QA of the applied technique was performed by standard addition method and it was found that the recovery for the investigated elements ranges for ICP-AES was 98.5-101.2%. Quality control was also ensured by standard moss reference materials M2 and M3, which are prepared for the European Moss Survey [39]. The measured concentrations were in good agreement with the recommended values. Fourteen isotopes 7Li, 49Ti, 59Co, 60Ni, 75As, 78Se, 98Mo, 107Ag, 114Cd, 118Sn, 121Sb, 202Hg, 205Tl and Pb including the isotopes 206Pb, 207Pb, 208Pb were determined using inductively coupled plasma mass spectrometry (Plasma Quant ICP-MS, Analytic Jena, Germany). External calibration was performed by measuring standard solutions containing these elements at 0.5, 1, 2, 5, 10 and 100 µg/L concentrations for following isotopes 7Li, 49Ti, 59Co, 65Cu, 75As, 98Mo, 107Ag, 114Cd, 118Sn, 121Sb, 205Tl and Pb (206Pb, 207Pb, 208Pb). For 202Hg, external calibration was performed by measuring standards 0.5, 1, 2, 5 and 10 µg/L. No internal standards were used. Detector attenuation mode was used to cover wide range of linearity without of collision gas. For the problematic isotopes such as 60Ni and 78Se, He was used as a collision gas. 60Ni was determinate by external calibration in range of 1, 5, 10, 50 100 and 500 µg/L. The 78Se was determinate in range of 1-50 µg/L. The limit of detection and quantification was automatically calculated by the instrument software (Aspect MS 4.3, 2017). Standard solutions were inserted into the sample sequence at every 20 samples to verify sensitivity and repeatability. The recoveries of these elements were 82.1-105.1%.

Statistical Analysis
The data obtained from chemical analysis of mosses samples were statistically processed. Statistical analysis was performed by Statistica 13 software package (StatSoft, Inc., Tulsa, OK, USA), whereas the date visualization was performed by means of several software packages: Statistica 13 (StatSoft, Inc., Tuls, OK, USA), QGIS and Surfer 17 (Golden Software, Inc., Golden, CO, USA). Some of the basic statistical quantities calculated, were the mean, Box-Cox transformed data mean, median, minimum, maximum, tenth percentile, ninetieth percentile, standard deviation, coefficient of variation, standard deviation (standard error), MAD (median absolute deviation), skewness, kurtosis, skewness for Box-Cox transformed data and kurtosis for Box-Cox transformed data. Pearson correlation coefficients were calculated for all the elements concentration in samples and values 0.30 and higher were considered as significant. In order to reveal the possible associations between elements throughout the sampling area and reduce the number of variables, multivariate cluster and factor analysis was performed. Then, from the associations that were identified, geogenic and anthropogenic origin of elements was established.

Results and Discussion
A first evaluation of the results is performed by descriptive statistics shown in Table 1. The original data were transformed by the Box-Cox method to improve the normality of distribution. Clearly, after transformation the absolute values of skewness and kurtosis were significantly reduced. The highest mean and Box-Cox mean concentrations are those for Ca, K and Mg, which also are the most concentrated elements in mosses samples. The element with the highest concentration of all is K (14,000 mg/kg) at sampling point Siboc-Obiliq, although Ca concentration (maximum concentration 12,000 mg/kg) is the highest for all other mosses sampling points, except for Siboc-Obiliq. Then, Al, Fe and P follow with maximum concentrations of 2700, 2000 and 1500 mg/kg respectively and Box-Cox means of 910 mg/kg Al, 820 mg/kg Fe and 770 mg/g P. High concentration of these elements in mosses is expected as they are the most common elements of the earth's crust and, as such, they can easily be transported by wind as dust from the surroundings and even the ground where mosses live in into the mosses surface.
The lowest minimum concentration was found for Tl 0.008 mg/kg and the Box-Cox mean of 0.016 mg/kg. Mo, Hg, Sb and as follow Tl with minimum concentrations of 0.054, 0.069, 0.073 and 0.087 mg/kg, respectively. Compared to the rest of the samples, As appears in higher concentrations in mosses from three sampling points, Zveçan, Stanterg-Mitrovicë and Cerajë-Leposaviq (6.85, 6.46 and 2.61 mg/kg), which obviously is a result of industrial processes of Trepça mines and smelters (smelter not operating since the year 2000). Then, the small number of moss samples with high content of as resulted in low value of P 90 which is 0.47 mg/kg. Although the sampling point Cerajë-Leposaviq, which is 11 km to the north of Zveçan is far to the north from Trepça facilities, it still can be influenced by this facility's emissions which can be transported by wind from the south to the north through the valley of Ibër river. Apart from Ca, K, Mg, Al, Fe and P, the minimum concentrations of which were high, all other elements minimum concentrations were ≤51 mg/ kg and most of them under 1 mg/kg (Table 1). It is worth noting the outstanding atmospheric deposition of Pb in sampling point Zveçan, 38 mg/kg and Ni in Harilaq, 79 mg/kg. Those concentrations represent also the maximum of these two elements and they can be explained by the industrial activity around the area of mosses sampling, Pb/Zn mining and ore processing activities in Zveçan and ferronickel mining processes in Harilaq. Otherwise, the Box-Cox transformed data mean for Pb is 5.2 mg/kg whereas for Ni is 1.7 mg/kg. Pb P 90 is 10 mg/kg whereas that of Ni is 6.1 mg/kg, which are much lower compared to the maximum concentration, indicating much smaller concentration for the remaining samples. X-arithmetical mean, Md-median, Min-minimum, Max-maximum, P 10 -10 percentile, P 90 -90 percentile, S-standard deviation, S X -standard error of mean, CV-coefficient of variation, MAD-median absolute deviation, A-skewness, E-kurtosis, (BC)-Box-Cox transformed data, a -Determined by ICP-AES; b -Determined by ICP-AES.
In Table 2, comparatively median concentrations and ranges for Kosovo (2019 samples), two neighboring countries (2015 samples) and Norway (2015 samples) as a pristine area, are presented. Calcium, K, Mg, Al and Fe are present at highest median concentrations in all countries. In Kosovo Ca, K and Mg have the highest concentrations compared to other countries, although there is not a very big difference. Mn median concentration in Kosovo is very much lower than in Norway, lower than in North Macedonia and higher than in Albania, whereas Na median concentration in Kosovo is the lowest. Cr, Cu and Ni are found in lower median concentrations in Kosovo than in Albania and North Macedonia, but Zn and Pb median is higher in Kosovo. Those five elements (Cr, Cu, Ni, Zn and Pb) appear in this order in compared countries because of the differences in geochemistry of the countries and industrial processes involving mineral exploitation [37,40]. Arsenic has a median of 0.17 mg/kg in Kosovo, which is lower than that of Albania and North Macedonia but higher than that of Norway. However, as concentration around the polluted site of Trepça facilities in Mitrovica is the highest of all, 6.8 mg/kg. Compared to Norway, most of the elements in mosses samples in Kosovo are at greater median concentration. Only Cu, Mn, Na, Ti and Tl are found at higher medians in Norway than in Kosovo, with Mn being the most distinct element which median is more than four times higher in Norway than in Kosovo. However, although the median and minimum concentrations of most of the elements are higher in Kosovo, the maximum concentration of many of the pollution related elements such as Co, Cr, Cu, Hg, Ni, Sb, V and Zn, are higher in Norway. Pb is the most distinctive pollutant element in Kosovo with regard to concentration compared to Norway, with the minimum concentration of 0.58 mg/kg which is higher than its maximum concentration in Norway that is 0.4 mg/kg.
In comparison with the moss 2010 survey in Kosovo [38], the median values of, Cr (2.63 mg/kg), Ni (2.0 mg/kg), Pb (7.78 mg/kg), Zn (38.5 mg/kg) and Fe (288 mg/kg), are higher than those of the present study which are: Cr (2.6 mg/kg), Ni (1.7 mg/kg), Pb (7.3 mg/kg), Zn (31 mg/kg) and Fe (820 mg/kg); meanwhile Cd (3.053 mg/kg), Pb (47.78 mg/kg) and Fe (3082 mg/kg) have lower maximum concentration then the maximum found in the present survey, which is for Cd 2.1 mg/kg, Pb 38 mg/kg and for Fe 2000 mg/kg. However, neither reliable comparison nor any reliable trend can be established, because this is only the second survey and also none of the sampling points are the same in both surveys. It is though worth noting that the differences are not very big. The association of elements was performed by hierarchical cluster analysis and factor analysis, based on the Pearson's correlation coefficients between elements. The Pearson's correlation coefficients of 25 elements analyzed in 45 moss samples are given in Table 3. There are 24 significant correlations with values of Pearson's coefficient between 0.5 and 0.94, of which 18 have values ≥0.5 and ≤0.7, whereas six of them were ≥0.7. Since there are e lot of numbers, the significant Pearson's coefficients shown in the table are given in bold so that they can be easily noticed. The strongest association is that of Fe and Al with a value 0.94 and they also have been found in great correlation in other works [35,36,41,42]. High correlation was found also for K and P with a value of 0.79.
Pearson's correlation coefficients matrix was used for the calculations in multivariate factor analysis, which resulted in data reduction down to four factors, given in Table 4. The analysis was performed on Box-Cox transformed data and the factors obtained were rotated orthogonally by the varimax method. The factors produced by factor analysis, as expected, were in decreasing order from F1 to F4 by the amount of variability that each of them could explain, and the significant loadings are given in bold to stand out of the rest. The threshold of the variables loading to be taken in consideration as important for factor value construction was set as ≥0.5.  The variables that loaded mostly on factor 1 were Al, Fe, V, Cr, Ti, Ba, Sr. This factor is the strongest and represents 24.4% of the total variability of 69.1%. The second factor groups five elements, As, Cd, Sb, Tl, Zn and it explains 21.4% of the total variance of all measured elements. Factor 3 explains 14.2% of the total variance with highest loading scores of Cu, Ca, K and P. Finally, factor 4, the weakest, explains 9.1% grouping only two elements Co and Ni. In the factor analysis seven elements were excluded (Hg, Li, Mg, Mn, Mo, Na and Pb) because they did not show any tendency to be associated to any group or other element.
The hierarchical clustering analysis is represented in a dendrogram shown in Figure 4. The clusters that can be observed were generated by taking the percentage of D(link)/D(max). Five groups of elements could be identified and they are compatible with elements association in factor analysis (except Cu that is found in two groups, 3 and 5) from which five factors or associations were extracted. Then, Al and Cr were the elements which correlated significantly with a greater number of other elements, Al correlated with six and Cr with five elements.
Factor analysis revealed four groups of elements which associated mostly to the given factor (four factors). Factor 2 can be considered as an anthropogenic factor, whereas factors 1, 3 and 4 are mixed factors, because they do not seem to be loaded significantly only by elements which are expected to be mostly or completely of anthropogenic origin. In Figure 5 are shown the spatial distribution map of the first extracted factor (A), and the standard factor values according to five districts of Kosovo and geological formations (B). This factor groups (Al, Ba, Cr, Fe, Sr, Ti, Tl and V), is the strongest factor as it represents the highest value of variability, 24.4% of the total 69.1% explained variance by all four factors. The highest values of this factor score appears in the east of Kosovo. The elements of this factor are considered as naturally distributed as common earth crust's components. Their provenance can be explained by the settlement of local fine soil particles previously airborne as dust during winds. Heavy rains can send soil particles into the mosses also. The presence of Cr though, can partly be as a result of mining in Novobërdë, cement facility in Hani i Elezit, ferronickel facility in Drenas (Glogovc) and the lignite exploitation in Obiliq. Moreover, Cr was determined to be the most concentrated metal in the lignite coal (19 mg/kg), excavated in Obiliq, for the purpose of electricity generation [43]. Consequently, it can be present in fly ash and transported by winds to other areas. The highest standard values of this factor are those of the district of Prishtina and Gjilan ( Figure 5). The district of Prizren represents the third highest standard value of the first factor, where traffic can be considered its main anthropogenic source, apart natural occurrence [44]. As can be seen in Figure 5, according to geological formations, this factor has the highest standard values in the clastitic rocks formed during Paleogene and Neogene periods of the Cenozoic era. These formations occur mostly in the Kosovo and Dukagjini basins but also in the center and other parts of Kosovo. High standard factor values are also found for Paleogene and Neogene magmatic rocks and Mesozoic flysch. The hierarchical clustering analysis is represented in a dendrogram shown in Figure  4. The clusters that can be observed were generated by taking the percentage of D(link)/D(max). Five groups of elements could be identified and they are compatible with elements association in factor analysis (except Cu that is found in two groups, 3 and 5) from which five factors or associations were extracted. Then, Al and Cr were the elements which correlated significantly with a greater number of other elements, Al correlated with six and Cr with five elements. Elements which load most significantly in factor 2 are As, Cd, Sb, Tl and Zn. The distribution map and the standard factor values according to districts and geology are shown in Figure 6. The element loading with the highest value is Zn and Sb load is the smallest of all. This is the second most important factor as it explains 21.4% of the total variance throughout the variables. In this group As, Cd and Zn are of anthropogenic origin as they are always associated with industrial activities [37,40,45]. However, it was found that Sb also co-exists with ores PbS and ZnS in the area of Trepça mines [46]. In a study performed in 2006 on the distribution of heavy metals in floodplain of agricultural soil along the Ibër River in the district of Mitrovica, almost the same group of elements was revealed by multivariate analysis [46]. These elements can be transported from this agricultural area by wind, to mosses locations during dry seasons lodged in fine particles of dust and this way contributing to factor 2. As can be seen in the distribution map of factor values, the highest values of this factor are spread along the Kosovo basin, it extends from Leposaviç in the north to the south in the Hani i Elezit, with a branch in the east between the city of Prishtina and Gjilan. The presence of these elements in this group clearly is a result of the pollution introduced into the atmosphere from anthropogenic activities. Industrial sites are located all along the Kosovo basin, starting in the north with the Trepça mines, lignite coal quarries and electricity generation plant, Kishnica mineral processing facility, Magnezite mine in Golesh, down to cement production plant in Hani i Elezit. The highest standard value of the second factor ( Figure 6B) is for the district of Mitrovica, which is expected because of heavy mineral exploitations and processing in that area; hence, the presence of factor 2 elements in the moss samples. Prishtina and Gjilan districts have also high factor 2 standard scores, which indicates that these two regions are also exposed to air pollution from industry. Particularly, regarding Prishtina, not only is it easily accessible for polluted air coming from Obiliq (lignite powered power plants) and Mitrovica (Trepça facilities and tailings), but also there is heavy usage of lignite coal in the city for domestic heating during winter and also there is very dens traffic, which obviously increases the introduction of these elements as pollutant in the atmosphere [32,47]. According to geological formations, the highest standard value of the factor 2 is found in magmatic rocks of Neogene and Paleogene periods. These formations cover relatively small areas and are found mostly in the north and the east of Kosovo, in the districts of Mitrovica and Prishtina. Antimony and Tl are found at higher concentrations in the north in Stanterg and Zveqan, probably with both geogenic and anthropogenic origin. Paleozoic era carbonates and metamorphic rocks also show high standard values for factor 2, although much lower than magmatic rocks.
given factor (four factors). Factor 2 can be considered as an anthropogenic factor, whereas factors 1, 3 and 4 are mixed factors, because they do not seem to be loaded significantly only by elements which are expected to be mostly or completely of anthropogenic origin. In Figure 5 are shown the spatial distribution map of the first extracted factor (A), and the standard factor values according to five districts of Kosovo and geological formations (B). This factor groups (Al, Ba, Cr, Fe, Sr, Ti, Tl and V), is the strongest factor as it represents the highest value of variability, 24.4% of the total 69.1% explained variance by all four factors. The highest values of this factor score appears in the east of Kosovo. The elements of this factor are considered as naturally distributed as common earth crust's components. Their provenance can be explained by the settlement of local fine soil particles previously airborne as dust during winds. Heavy rains can send soil particles into the mosses also. The presence of Cr though, can partly be as a result of mining in Novobërdë, cement facility in Hani i Elezit, ferronickel facility in Drenas (Glogovc) and the lignite exploitation in Obiliq. Moreover, Cr was determined to be the most concentrated metal in the lignite coal (19 mg/kg), excavated in Obiliq, for the purpose of electricity generation [43]. Consequently, it can be present in fly ash and transported by winds to other areas. The highest standard values of this factor are those of the district of Prishtina and Gjilan ( Figure 5). The district of Prizren represents the third highest standard value of the first factor, where traffic can be considered its main anthropogenic source, apart natural occurrence [44]. As can be seen in Figure 5, according to geological formations, this factor has the highest standard values in the clastitic rocks formed during Paleogene and Neogene periods of the Cenozoic era. These formations occur mostly in the Kosovo and Dukagjini basins but also in the center and other parts of Kosovo. High standard factor values are also found for Paleogene and Neogene magmatic rocks and Mesozoic flysch. Elements which load most significantly in factor 2 are As, Cd, Sb, Tl and Zn. The distribution map and the standard factor values according to districts and geology are shown in Figure 6. The element loading with the highest value is Zn and Sb load is the smallest of all. This is the second most important factor as it explains 21.4% of the total Factor 3 groups four elements, Ca, Cu, K and P, and it explains 14.2% of the total variance ( Figure 7). The distribution map of the factor values, as well as its standard values according to districts and geology of the studied area, are given in Figure 8. The highest values of the third factor are concentrated in the center and the south of Kosovo, around the lignite quarry in Obiliq and the cement plant in Hani i Elezit. This factor's standard values are the highest for the district of Gjilani, Mitrovica and Prishtina. This group is composed mostly by naturally occurring elements; however, anthropogenic origin is to be considered. K and P are widely used in agriculture as fertilizers, thus, during dry seasons K and P enriched soil particles can be carried by winds to mosses habitats. Ca, K and P are also found in significant quantities in the fly ash of Obiliq power plant [48], which regularly is spread out by winds and can contribute to the atmospheric deposition on mosses. Cu is mostly of anthropogenic origin and it is probably related to mining processes. Moreover, Cu is also found in fly ash of Obiliq power plant [48], but it can also reach mosses habitats from traffic emissions, since road environment was found to contain among other metals Cu which originates from vehicular usage [44,49,50].
phere [32,47]. According to geological formations, the highest standard value of the factor 2 is found in magmatic rocks of Neogene and Paleogene periods. These formations cover relatively small areas and are found mostly in the north and the east of Kosovo, in the districts of Mitrovica and Prishtina. Antimony and Tl are found at higher concentrations in the north in Stanterg and Zveqan, probably with both geogenic and anthropogenic origin. Paleozoic era carbonates and metamorphic rocks also show high standard values for factor 2, although much lower than magmatic rocks.  Factor 3 groups four elements, Ca, Cu, K and P, and it explains 14.2% of the total variance (Figure 7). The distribution map of the factor values, as well as its standard values according to districts and geology of the studied area, are given in Figure 8. The highest values of the third factor are concentrated in the center and the south of Kosovo, around the lignite quarry in Obiliq and the cement plant in Hani i Elezit. This factor's standard values are the highest for the district of Gjilani, Mitrovica and Prishtina. This group is composed mostly by naturally occurring elements; however, anthropogenic origin is to be considered. K and P are widely used in agriculture as fertilizers, thus, during dry seasons K and P enriched soil particles can be carried by winds to mosses habitats. Ca, K and P are also found in significant quantities in the fly ash of Obiliq power plant [48], which regularly is spread out by winds and can contribute to the atmospheric deposition on mosses. Cu is mostly of anthropogenic origin and it is probably related to mining processes. Moreover, Cu is also found in fly ash of Obiliq power plant [48], but it can also reach mosses habitats from traffic emissions, since road environment was found to contain among other metals Cu which originates from vehicular usage [44,49,50]. Following the geology of the studied area, this group of elements is mostly related to Paleozoic carbonates (small areas in the east center and west, Figure 3) and metamorphic rocks (in the east). Then, somewhat smaller standard factor values are found for quaternary deluvium/proluvium deposits, Paleogene and Neogene clastites and Paleogene and Neogene magmatic rocks. These rocks are weathered under atmospheric conditions and ground to fine particles of soil. Consequently, they are carried by wind into the atmosphere and then deposited over other areas and this way onto mosses also, enriching them with elements they contain. Figure 8 shows the elements grouped by factor 4. This is the weakest factor constructed by only two elements, Co and Ni and explaining only 9.1% of the total variance. Its distribution map shows that the highest values appear in the center of Kosovo and in the west. The highest standard value is that of the district of Prishtina, then Peja and Priz- Following the geology of the studied area, this group of elements is mostly related to Paleozoic carbonates (small areas in the east center and west, Figure 3) and metamorphic rocks (in the east). Then, somewhat smaller standard factor values are found for quaternary deluvium/proluvium deposits, Paleogene and Neogene clastites and Paleogene and Neogene magmatic rocks. These rocks are weathered under atmospheric conditions and ground to fine particles of soil. Consequently, they are carried by wind into the atmosphere and then deposited over other areas and this way onto mosses also, enriching them with elements they contain. its origin is mostly natural, as it may be for Ni in the west of Kosovo. The geological formation with the highest standard value of the fourth factor is magmatic rock of Mesozoic era. Neogene and Paleogene clastites come second, then Mesozoic carbonates and Quaternary deluvial/proluvial deposits also show high standard factor values.

Conclusions
The atmospheric deposition of elements was studied by moss biomonitoring method in order to estimate air pollution in Kosovo. After the element's determination in 45 moss samples, variations of concentrations throughout the country were observed for each element. Peaks of concentrations in moss samples were detected for Pb, Zn, As and Cd, in the district of Mitrovica around the Trepça Pb-Zn ore processing unit and tailings dump as well as smelter plant with a landfill of slag deposit. In comparison to Norway and neighboring countries Albania and North Macedonia, mosses have higher Ca, K and Mg concentrations in Kosovo. Cr, Cu and Ni are found in lower median concentrations in Kosovo than in Albania and North Macedonia, but Zn and Pb median is higher in Kosovo. In fact, compared to all countries Pb, Zn, Hg and Cd are the heavy metals mostly polluting the atmosphere as they are found with the highest median concentration in mosses samples. Factor analysis revealed four factors, of them three were geogenic or mixed and one was anthropogenic. The anthropogenic factor includes As, Cd, Sb, Tl and Zn and is mostly present along the Basin of Kosovo. The most polluted districts of Kosovo are Mitrovica, Prishtina, Gjilan and Ferizaj, as a result of industry, traffic and domestic routines.    Figure 8 shows the elements grouped by factor 4. This is the weakest factor constructed by only two elements, Co and Ni and explaining only 9.1% of the total variance. Its distribution map shows that the highest values appear in the center of Kosovo and in the west. The highest standard value is that of the district of Prishtina, then Peja and Prizren have slightly smaller values. Ni presence in this group can be explained by the ferronickel smelter in Drenas where the ore concentration takes place and its mine in Golesh (Harilaq sampling point), both of these located around the center of Kosovo. However, no Co involving industry or processing of ore containing Co takes place in Kosovo, probably its origin is mostly natural, as it may be for Ni in the west of Kosovo. The geological formation with the highest standard value of the fourth factor is magmatic rock of Mesozoic era. Neogene and Paleogene clastites come second, then Mesozoic carbonates and Quaternary deluvial/proluvial deposits also show high standard factor values.

Conclusions
The atmospheric deposition of elements was studied by moss biomonitoring method in order to estimate air pollution in Kosovo. After the element's determination in 45 moss samples, variations of concentrations throughout the country were observed for each element. Peaks of concentrations in moss samples were detected for Pb, Zn, As and Cd, in the district of Mitrovica around the Trepça Pb-Zn ore processing unit and tailings dump as well as smelter plant with a landfill of slag deposit. In comparison to Norway and neighboring countries Albania and North Macedonia, mosses have higher Ca, K and Mg concentrations in Kosovo. Cr, Cu and Ni are found in lower median concentrations in Kosovo than in Albania and North Macedonia, but Zn and Pb median is higher in Kosovo. In fact, compared to all countries Pb, Zn, Hg and Cd are the heavy metals mostly polluting the atmosphere as they are found with the highest median concentration in mosses samples. Factor analysis revealed four factors, of them three were geogenic or mixed and one was anthropogenic. The anthropogenic factor includes As, Cd, Sb, Tl and Zn and is mostly present along the Basin of Kosovo. The most polluted districts of Kosovo are Mitrovica, Prishtina, Gjilan and Ferizaj, as a result of industry, traffic and domestic routines.