Pranvera Lazo This email address is being protected from spambots. You need JavaScript enabled to view it.1, Alda Kika2, Flora Qarri3, Lirim Bekteshi4, Shaniko Allajbeu1, Trajce Stafilov5

1 Department of Chemistry, Faculty of Natural Sciences, University of Tirana, Tirana, Albania
2 Department of Computer Science, Faculty of Natural Sciences, University of Tirana, Tirana, Albania
3 Department of Chemistry, University of Vlora, Vlora, Albania
4 Department of Chemistry, Faculty of Natural Sciences, University of Elbasan, Elbasan, Albania
5 Institute of Chemistry, Faculty of Science, Ss. Cyril and Methodius University, Skopje, North Macedonia

Received: January 4, 2022
Revised: April 1, 2022
Accepted: April 4, 2022

 Copyright The Author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.

Download Citation: ||  

Cite this article:

Lazo, P., Kika, A., Qarri, F., Bekteshi, L., Allajbeu, S., Stafilov, T. (2022). Air Quality Assessment by Moss Biomonitoring and Trace Metals Atmospheric Deposition. Aerosol Air Qual. Res. 22, 220008.


  • Atmospheric deposition of toxic metals was assessed by moss biomonitoring in Albania.
  • Natural and anthropogenic factors have affected the concentration of metals in current moss.
  • High anthropogenic contribution of metals was found in the eastern region.
  • Local anomalies of Cr, Ni, and Co reveal both natural and anthropogenic sources.
  • Cd, Cu, Pb, and Zn were higher in cross-border polluted areas and heavy traffic regions.


Carpet-forming moss species are useful tools for monitoring the atmospheric deposition of pollutants like metals, organic compounds, nitrogen, and more recently, microplastics. Even with a dense sampling scheme, moss biomonitoring could extend the environmental study to national and/or continental scales. Moss species (Hypnum cupressiforme) was used as a biomonitor to assess the atmospheric deposition of metals in the entire territory of Albania and to identify the probable anthropogenic factors affecting the levels of the most toxic metals/metalloids (As, Cd, Co, Cr, Cu, Fe, Hg, Ni, Pb, and Zn) in the current moss. Statistical analysis is used to examine the relationships between elements accumulated in moss samples, and to assess the air quality of Albania. Considering the geochemical normalized data using lithium as a normalizing element, Fe, Cr, Co, and Ni show obvious anthropogenic anomalies. GIS maps identified the local enrichment of these elements, which is stretched in the north-southeastern direction and affected mostly by anthropogenic sources. The elements Fe, Cd, Cu, Pb, and Zn look likely to be homogenously distributed by indicating a common natural origin on a national scale. It is confirmed by the association of these elements with iron which is a typical crustal element that originates mostly from lithogenic dust in sparsely vegetation regions. Fe, Cd, Cu, Pb, and Zn were mostly higher between a few regions with few obvious anomalies positioned in the center of the country with high industrial activity, in the areas with heavy traffic, and in the cross-borders of Albania with Monte Negro and Greece mostly emerging from the local contribution of anthropogenic activities, local lithology, and transboundary pollution sources. The elements Cd and Hg are linked mostly with the long-range transport of pollutants.

Keywords: Air quality, Moss biomonitoring, Toxic metals, Statistical analysis, Albania


Environmental pollution is strongly linked with the anthropogenic emissions of contaminants that show negative effects on the health of human beings and different ecosystems. For this reason, understanding and estimating the transmission, spatial distribution, and temporal trends of the pollutants in order to identify the problematic areas, their local sources, and temporal trends have been the focus of the scientists. Reviewed study done by Mamun et al. (2020) indicate that trace metals in atmospheric deposition revealed from diverse emission sources and environments including urban, industrial, suburban, rural, traffic, and coastal sites and are significantly affected by meteorological conditions and the size of dust particles. Generally, the crustal elements deposit at higher rates than those of the anthropogenic elements due to the larger particle sizes of crustal elements. Since metals are non-biodegradable and remain persistent in the environment for a very long time, they cannot be broken down, and remain present for an extended period in soils and sediments until they are eluted to other compartments, or can also react with other elements in the soil or sediment and form or degrade to more toxic compounds (Briffa et al., 2020).

Air pollution, and especially pollution from particles containing high content of toxic metals, has serious impacts on human health by reducing life expectancy and causing many diseases (Sokhi et al., 2022). Toxic and carcinogenic metals such as Cu, Fe, Cd, Cr, Ni, Hg and V, can produce reactive radicals which damage DNA, perform lipid peroxidation and depletion of protein sulphhydryls, accompanied by other health effects (Briffa et al., 2020).

Road traffic is considered as a potential source of metals coming from fuel and oil combustion products, tires, brake linings, bearing con+sumption, corrosion of vehicle components, road construction materials, etc. (Richards, 2020; U.S. EPA, 2004; Kummer et al., 2009; Harmens et al., 2015, 2013). Old vehicles and outdated industrial technologies used in developing countries, the smoke from accidental fires and/or waste incineration (Shaman and Wheeler, 2007; Grodzinska and Szarek-Lukaszewska, 2001), have led to the increased content of different pollutants in the environment which became a serious problem in many parts of the world. Besides, different interventions in the environment without criteria have brought serious consequences linked with the growth of natural and anthropogenic emissions of various contaminants. Among them, forest degradation and forest burning, intervention in river beds, and the increase in barren land areas have increased the rate of erosion, the flooding in certain areas, the content of soil dust, fine and ultra-fine particulate matter (PM) in the air, and metal content near the emission sources. The increased level of toxic metals, Pb, Cr, As, Cd, Hg, Zn, Cu, Co, and Ni, is an inorganic chemical hazard to humans and environmental compartments (Järup, 2003). Atmospheric deposition removes pollutants from the atmosphere to natural ecosystems by causing negative effects on human health, crop yields, land, and marine ecosystems (WMO, 2019). Nowadays studies suggested that PM content in the air presents adverse effects on pulmonary and cardiovascular health by making this category of people more vulnerable to the severity of COVID-19 symptoms (Bourdrel et al., 2021; Dales et al., 2021; Donaldson and Mac Nee, 2001; Fang et al., 2021). The data of the global COVID-19 pandemic indicate that air pollution has a direct impact on the increasing mortality among infected people (Brunekreef et al., 2021; Fang et al., 2021; Fattorini and Regoli 2020; Pozzer et al., 2020; Zhang et al., 2020). Such effects make it necessary to set strict parameters regarding air quality and develop a continuous monitoring program to control the pollution level and the effects of pollutants.

The conventional technique of monitoring requires some expensive equipment limited only to the monitoring stations. For more than 60 years, moss biomonitors have been developed and widely used as an alternative method to assess air pollution and define and characterize metal pollution sources in atmospheric deposition in European countries (Berg et al., 1995; Basile et al., 2008, 2009; Macedo-Miranda et al., 2016; Esposito et al., 2018), and more recently in Asia, Brazil, and North America (Stankovi et al., 2018; Harmens et al., 2015, 2013; Gjengedal and Steinnes, 1990; Steinnes, 1989). Bryophytes (particularly mosses) have been used as biomonitors for assessing metals in atmospheric deposition since the early ’60s of the last century (Tyler, 1990; Harmens et al., 2015; Frontasyeva et al., 2020). Biomonitors are effective species that accumulate and document the impact of different pollutants on living organisms and the impacts of different pollutants to different ecosystems (Wright et al., 2018). Based on the accumulation capacity and sensibility to various chemicals, bryophytes can be used as accumulation indicators or as sensitive indicators, mainly for polluting substances (Zechmeister et al., 2003). Pleurozium schreberi, Hylocomium splendens, and Hypnum cupressiforme mosses species were recommended in the experimental protocol of the European Moss Survey for assessing metals in atmospheric deposition (ICP Vegetation, 2010). The accumulation of metals in mosses varied greatly according to the type of moss species (Zechmeister et al., 2003; Schröder et al., 2016; Schröder and Nickel, 2019). It is important to use the same or similar moss species for passive biomonitoring of air quality because it makes it possible to compare results from different studies (Boquete et al., 2020). According to Bargagli et al. (2002), Hypnum cupressiforme is widely used in trace element biomonitoring and in bioaccumulating studies of trace element pollution due to its wide-ranging distribution and because of the abundant data existing in the literature on this species. In such conditions, it is ideal for making comparisons with other areas of Europe. The importance of such studies is related to the fact that it can be realized both in local and intercontinental ranges and in different time periods by providing information on the spatial distribution of pollutants and the changes over time (Harmens et al., 2004; Meyer et al., 2014). Understanding the spatial and temporal distribution of total (wet and dry) deposition of pollutants on a global scale is critical in identifying the most vulnerable populations, ecosystems, and farmlands (WMO, 2019) a process easy doing by biomonitoring.

This paper deals with moss biomonitoring surveys in Albania, which is part of the framework of the International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops (ICP Vegetation). ICP Vegetation in European scale has evaluated atmospheric deposition of metals since 1990 by using moss biomonitoring at five years intervals (Rühling, 1994; Rühling and Steinnes, 1998; Harmens et al., 2011, 2013, 2015; ICP Vegetation, 2010). The aim of this research is to assess the atmospheric deposition of metals in the entire territory of Albania using low-cost moss biomonitoring, to identify the most probable sources of toxic metals in the current moss, and to differentiate between their natural and anthropogenic through the use of concentration and normalized data.


2.1 Sampling

The sampling process significantly affects the uncertainty of the analytical results. Thus, to guarantee the representative samples relatively free from the interference of external factors, sampling was carried out in accordance with the LRTAP Convention-ICP Vegetation protocol and sampling strategy of the European Program on Biomonitoring of Heavy Metal Atmospheric Deposition (ICP Vegetation, 2010). The moss samples (Hypnum cupressiforme (Hedv.), widely spread in Albania, were collected from 47 sampling sites (approx. 1.5 moss samples per 1000 km2) during the relatively dry periods of October–November 2010, and June–July 2011. The locations of the sampling sites were situated at least 300 m away from main roads or buildings and 100 m from small roads and single houses, located mostly in open areas. Composite moss samples were formed by five to ten sub-samples collected within an area of 50 × 50 m2. Disposable polyethylene gloves were used during the sampling and sample cleaning to avoid contamination of the samples. The sampling site distribution is depicted on a map of Albania (centered at 41°00′ north of the equator and 20°00′ east of Greenwich) (Fig. 1).

Fig. 1. The location of the sampling sites on the map of Albania position.Fig. 1. The location of the sampling sites on the map of Albania position.

2.2 Moss Analysis

Moss samples were cleaned from the adhering materials and brown parts of the plant tissues. Only the green and green-brown parts of moss tissues that represent at last three years of moss growth were selected for chemical analysis. Samples were opened separately on big paper sheets and then were dried at room temperature for about 72 hours. The content of metals in moss samples was determined by inductively coupled plasma-atomic emission spectrometry (ICP-AES) (Varian, 715ES equipped with ultrasonic nebulizer CETAC (ICP/U5000AT+) for better sensitivity of the elements), and electro-thermal atomic absorption spectrometry (ETAAS, SpectrAA 640Z for Cu and Pb) after acid digestion. Around 0.5 g moss samples were digested in Teflon digestion vessels by adding 7 mL HNO3 (69%, m/V) and 2 mL H2O2 (30%, m/V). Analytical grade reagents were used in this study: nitric acid, trace pure (Merck, Germany); hydrogen peroxide, p.a. (Merck, Germany); and bi-distilled water. The dried moss samples were digested with Mars, CEM, USA Microwave digestion systems with the following programme: step 1—temperature 180°C, 10 min ramp time, with power of 500 W and 20 bar pressure; step 2—temperature 180°C, 20 min hold time, with power of 500 W and 20 bar pressure (Stafilov et al., 2018). Finally, the vessels were cooled and carefully opened. The digested samples were quantitatively transferred into 25-mL calibrated flasks. The ICP-AES and AAS analyses were conducted in the Institute of Chemistry, Faculty of Science, St. Cyril and Methodius University, Skopje, North Macedonia.

The digested samples with nitric acid (9:1) at half-pressure Teflon tubes were analyzed for mercury content by using the cold vapor atomic absorption spectroscopy (CVAAS) as prescribed by Lazo et al. (2018). Varian 10+ atomic absorber equipped with homemade Hg cold vapor system was used for Hg determination. The analysis of Hg has been performed at the University of Tirana, Faculty of Natural Sciences.

The detection limits of elements calculated as 3 SD of the lowest instrumental measurements of the blanks are as follows: 0.002 mg kg−1 for Zn; 0.005 mg kg−1 for Cd and Fe; 0.01 kg−1 for As, Cu; 0.05 mg kg−1 for Cr, Co and V; 0.25 mg kg−1 for Ni and Pb, and 0.002 mg kg−1 for Hg. The limits of quantification calculated as 10 SD of the lowest instrumental measurements of the blanks are as follows: 0.007 mg kg−1 for Zn; 0.02 mg kg−1 for Cd and Fe; 0.035 kg−1 for As and Cu; 0.2 mg kg−1 for Cr and V; 0.85 mg kg−1 for Ni and Pb, and 0.008 mg kg−1 for Hg. Three replicates per moss sample were digested and three run measurements per digest were performed.

2.3 Quality Control

The quality assurance of ICP-AES was checked by two moss reference materials, M2 and M3, prepared first for the 1995/6 European moss survey (Steinnes et al., 1997). Blank samples were simultaneously analyzed with moss samples. The recovery of the investigated elements was checked by the standard addition method. It ranged between 98.5% and 101.2% for ICP-AES, and 96.9% to 103.2% for AAS.

2.4 Data Processing and Statistical Analysis

Descriptive statistics and spatial analysis were used as the first step of statistical analyses to investigate the variability and spatial distribution of the data. The concentration data (CD) were divided by their medians to standardize those at the same digit numbers. To account for natural metal variability and to highlight the anthropogenic component, the data were normalized using Li as a normalizing element (Loring and Rantala, 1992). The relationships between metals present in moss provide an indication of the factors that have influenced their distribution and origin. It is tested by Pearson correlation analysis and confirmed by the statistical significance at p < 0.05. Factor analysis (FA), based on the correlation matrix of concentration (CD) and normalized concentration data (NCD), may explore the hidden multivariate structures of the data (Astel et al., 2008; Reimann et al., 2002) and help to clarify the relationship between metals that may have similar origins or develop similar associations in moss samples. Statistical analysis was performed using the MINITAB 19 software package. The spatial distribution of the elements was visualized from the GIS maps, plotted with the Arc-GIS 10.2 system by applying the local deterministic methods and the inverse distance weighting.


3.1 Trace Metal Concentrations in Moss Samples

The most important statistical parameters, such as the concentration range, coefficient of variation (CV%), skewness, and kurtosis are shown in Table 1.

Table 1. Descriptive statistic analysis of trace metal data (N = 47, * in mg kg–1).

The data of CDs and NCDs of As, Cd, Hg, and Ni follow the lognormal distribution model, are positively skewed, and show high variability (CV > 75%). It indicates high asymmetry in data distribution. The sequence of the CDs of elements in moss samples were similar with that of NCDs, Fe > Cr > Zn > Ni > Cu > Pb > Co > As > Hg > Cd, except Zn and Cr NCDs that changed the position between them, Zn(N) > Cr(N)). The statistical parameters that characterize the asymmetry of the elemental data differ significantly that indicate different behaviors of the elements in the study areas. Several activities were developed in different regions over the last century such as mining activity, iron, ferrochromium, nickel, and copper metallurgy, chemical and tannery industry since the early ‘50s till the beginning of the ‘90s of the last century. It could explain high variability and the enriched anomalies observed in different areas. High contents of Cr, Zn, Ni, Cu, and Pb evaluated by respective median values were followed by high variation (CV% > 75%) would suggest an effect of a polluted air mass derived from long-transport or regional air pollution. In similar way as other global regions, the data collected from this study show that the content of major elements (Fe, Ca, Mg, K, and Na), represented here by Fe, are orders of magnitude greater than those of routinely monitored anthropogenic elements (Zn, Pb, Cu, Ni, Cd, Cr) (Cheng et al., 2021; Feng et al., 2021; Sokhi et al., 2022).

The use of mosses as a bioindicator in biomonitoring of metal atmospheric deposition is a relative method for identifying problematic areas based on the geographical fluctuation of metal concentrations. As different moss species accumulate metals in different ways, the allowed limits of metal have yet to be defined. In such conditions, the use of the same species for monitoring is recommended. To assess the extent of pollution in Albania, we compared the median values of metals in the present moss with the respective European medians, as shown in Fig. 2.

Fig. 2. Sequence plot of mean, median, and European median values (in mg kg–1).Fig. 2. Sequence plot of mean, median, and European median values (in mg kg–1).

In general, the median concentration of elements in moss samples of Albania, except Cd and Zn, are higher than the respective European median values by indicating higher contamination level of these elements in atmospheric deposition of Albania, particularly for As, Cr and Ni.

3.2 Multivariate Analysis

Pearson correlation analysis was carried out to CDs and NCDs to investigate the linear relationship and the association between them. The obtained ‘’p’’ values of Pearson correlation between CDs and NCDs (r > 0.4) ranged from 0.000 to 0.005 (Table S1 that meets perfectly the significance condition (p < 0.05) for a reliable correlation. Based on Pearson correlation data, the elements are separated into three groups: 1) the 1st group shows very strong and significant correlations (r > 0.8, p ≤ 0.001) between CD-NCD data of the same element (Cr, Co, and Ni) characterized by a narrow confidence interval (CI) values (Table S1). It reflects a common chemical behavior of these metals and their common origin as mineral oxides and geochemical factors. 2) the 2nd group show strong and significant correlation (0.6 < r < 0.8, p < 0.005) between the pairs of the CDs of Cu, Cd, Zn, Pb, Hg with their respective NCDs (Table S1); 3) the 3rd group show moderate and significant correlation (0.4 < r < 0.6, p < 0.01) between the pairs of CDs of Cd, Cu, Pb, Hg, and Zn. The presence of these metals in the atmospheric deposition is mostly derived by long-range transport of pollutants (Cd and Hg) from other parts of Europe (LRTP) and from local emitting sources (Cu, Cd, Zn, Pb) such as high-temperature metal processing, traffic emission, and windblown dust that represents historical waste deposition and geochemical properties of the area (Harmens et al., 2015). Factor analysis (FA) was carried out to better explain the associations of the elements, to split the data into different factors based on the correlation matrix of elements in current moss, and to evaluate their probable sources in moss samples. The most important factors were extracted for an eigenvalue > 1 after varimax rotation. The associations of the elements with high loadings in the same factor (Reimann et al., 2002) are assumed as similar effects or similar characteristics of the elements. Only factor loadings (FL) larger than 0.4 are used to interpret the associations of the elements at the same factor. The results of FA are shown in Table 2, and are illustrated in GIS maps of Fig. 3

Table 2. Factor Analysis of the standardized CDs and NCDs of 2010 moss survey (N = 47). Factor Analysis of the Correlation Matrix;  Varimax Rotation; Sorted Rotated Factor Loadings and Communalities

Factor analysis revealed five main factors that represent 80.9% of the total variance. The first factor, (F1), is the strongest factor representing 27.7% of the total variance. The next two strong factors (F2 and F3) represent 32.4% of the total variance, respectively 16.3 and 16.1 %. The last two weak factors (F4, and F5) represent only 20.8% of the total variance (respectively 10.4, and 10.4 %). GIS maps of each factor are shown in Fig. 3. The associations of metals within the same factor represent similar properties and similar behaviors in the environment and could be explained as follows:

Factor 1 (F1) contains high loadings of CDs and NCDs of Cr, Ni, Co, and Fe. Strong associations of Cr, Ni, Co, and Fe CDs with their respective NCDs indicate strong anthropogenic sources of these elements. Cr, and the associated elements, enter into various environmental matrices (air, water, and soil) from several natural and anthropogenic sources, and the highest releases are attributed to industrial emission sources, such as metal processing, stainless steel welding, ferrochrome metallurgical, refractory, and chemical industries (Tchounwou et al., 2012), geogenic factors, and mining industry. The highest level of these elements was found in the eastern mineral belt rich with Cr, Ni, and Fe minerals. It indicates the geogenic origin of these elements. Soil dust emission is responsible for high levels of historical deposition of Cr, Ni, and Co in the soil throughout Albania. It has led to high soil concentrations of these metals, which means that the soil dust represents a major emission source of Cr, Ni, and Co in the air (Allajbeu et al., 2021; Lazo et al., 2019, 2021). The association of these elements with Fe, a typical crustal element, highly supports the discussion above. Arc GIS map is created to show the distribution pattern of high loading parameters of F1 (Fig. 3(a)).

Fig. 3. Distribution pattern of the standardized average factor loadings (FL) at each factor.
Fig. 3. Distribution pattern of the standardized average factor loadings (FL) at each factor.

The next strong factor (F2) is linked with high loadings of Cd, Zn, Cu, and Hg normalized data, obtained after the compensation of their natural variations and indicating mostly the anthropogenic fractions (Loring and Rantala, 1992). Hg and Cd are typical anthropogenic elements mostly from long-range transport of the pollutants, as well as from fuel combustion, waste incineration, and medical chemicals treatment indicating the possibility of anthropogenic sources (Tchounwou et al., 2012). A high anomaly of F2 is found in the areas with heavy traffic, and the next in the cross-borders with Monte Negro in the North-West and Greece in the South-East. Cu and Zn are known as tracers of non-exhaust brake and tire wear, so traffic emission followed by automobile exhaust, lubricant consumption, tire wear, brake wear, and mechanical parts friction are supposed as important sources of Cu and Zn (Adamiec, 2017; Harrison et al., 2003). The anomalies of Cd, Cu, Hg, and Zn NCDs in the cross border areas may indicate the effects of transboundary pollution. The lowest concentrations of these elements were found in the south. On the other hand, their distribution patterns in most parts of the territory look likely homogeneous by indicating the effects of long-range transport, as well as the effects of soil geochemistry and soil dust fine mineral particles. They may also originate from the oil-gas industry, and shipping activity in the west. Hg and Cd point also the sources of fire industrial activities, waste incineration, while the association of Cu and Zn indicates the traffic emission and atmospheric deposition sources. Hg and Cd are typical anthropogenic elements probably entrapped to soil dust fine particles. On the other hand, the association of Cd, Cu, Hg, and Zn NCDs may also indicate the geogenic factors, mining, and nonferrous metallurgy as important sources of these elements (Pacyna and Pacyna, 2001; UNEP, 2013). Although the copper mining, smelting, and processing industry in Albania had stopped since the beginning of ’90 years, and the emissions of metals such as copper, lead, cadmium, zinc, are decreased significantly, their concentrations in some parts of the country are still very high due to the historical deposition and the effects of mineral damps in the vicinity of the ex-copper industry particularly in the North (Lazo et al., 2019). Arc GIS map is created to show the distribution pattern of high loadings parameters of F2 (Fig. 3(b)).

Factor 3 (F3) contains high CD loadings of Fe, Hg, Zn, Cd, and Cu. These elements have both natural and anthropogenic origin. The presence of Fe together with these elements under the same factor supports the natural origin of these elements. Fe is naturally distributed as a typical soil element (Rudnick and Gao, 2003) that may indicate its soil dust origin. We noticed that CDs of this group explained a homogenous distribution, and were mostly higher between a few regions (Fig. 3) by indicating a common natural origin on the country scale. The highest loadings of Fe, Hg, Zn, Cd, and Cu are positioned in the center of the country with high industrial activity. Arc GIS map of Fig. 3(c) shows the distribution pattern of loading parameters of this factor.

Factor 4 (F4) and Factor 5 (F5) are associated with high loadings of Pb and Pb(N) (F4), and As and As(N) (F5). The main source of Pb is counted to be the vehicle exhaust, coal combustion, and industrial emission sources (Tchounwou et al., 2012) probably explaining the Pb anomaly in the central part of Albania. The strongest Pb anomaly founded in the cross border with Monte Negro is probably affected by the transboundary pollution of metals erased from the Al and Fe processing industry in Monte Negro (Lazo et al., 2018; Peck, 2004). All these findings show a strong anthropogenic origin of Pb in the current moss samples. A wide area with relatively high Pb content (13–21 mg kg1) was found in the North. It is probably originated from geogenic factors. Besides, GIS maps were created by using the local deterministic interpolation techniques that create surfaces from the measured neighboring sampling points, based on the extent of similarity of the data. The stretching of the surfaces with a similar concentration range of the elements depends on the distances between sampling points. High-density sampling points can create more accurate surfaces. Due to the conditions of the relief, the distance between sampling points in the North is relatively high, so the extension of the anomalies are relatively wide.

The highest As content is found in the South. Arsenic was used in agriculture activity as inorganic fertilizers and herbicides. It may lead to the enrichment of As in the local soils which may enter the atmosphere as fine soil dust particles, and then be re-suspended as atmospheric deposition in different distances controlled by the size of dust particles. The next area with relatively high loadings of As and As(N) belongs to the mineralized belt in the Eastern part of the country that is probably linked with geogenic factors and soil geochemistry of the area. The same as Pb, the extension of the areas with relatively high As content (11–19 mg kg1) in N–E and S–E directions were caused by the low density of sampling sites. For more accurate distribution, it is necessary to increase the density of sampling sites in different parts of the country. The distribution pattern maps of F4 and F5 are shown in Figs. 3(e) and 3(f).


The atmospheric deposition of metals was assessed in the entire territory of Albania by using moss as a bioindicator. It deals with the evaluation of the pollution level on the atmospheric deposition of trace metals at a national scale. Results showed that metal contaminations were arise from local anthropogenic activities, in particular, industry, road traffic, transboundary pollution, sea spry emission, and long-range transport of the pollutants from other parts of Europe. It is a complete research study that identified the main factors contributing to air quality in Albania. The study could contribute to other studies documenting the deposition of various elements that have been conducted in other parts of Europe by using mosses as bioindicators of trace metals atmospheric deposition under the framework of the ICP Vegetation Programme.

It is shown that moss data could provide a complementary measure of elemental deposition from the atmosphere to the terrestrial ecosystems, which makes it possible to achieve a relative comparison of the deposition data achieved at a high sampling density in national/and or continental monitoring.

The variability in the distribution pattern of the elements may identify the areas within the country with high levels of trace metals that need further investigations to determine the health and environmental risks. High concentration levels of Cr, Ni, Fe, and Co were found in the eastern part of the country, higher than the neighboring and EU countries.

The present study represents the need for more rigorous measures regarding the emission of atmospheric pollutants originating from anthropogenic sources. A continuous monitory program is necessary to control the pollution level and the effects of pollutants.


The study was conducted by the authors under the framework of UNECE-ICP Vegetation Programme and the Doctoral Program of University of Tirana, Faculty of Natural Sciences. The survey complies the ethical rules of ICP Vegetation Programme and the Doctoral Program of University of Tirana. It does not contain any study or experiment done with human participants by the participants or the third part.


The authors express the gratitude to the staff of the Institute of Chemistry, Faculty of Science, Sts. Cyril and Methodius University, Skopje, Macedonia for analyzing the Albanian moss samples.


  1. Adamiec, E. (2017). Chemical fractionation and mobility of traffic-related elements in road environments. Environ. Geochem. Health 39, 1457–1468.

  2. Allajbeu, S., Bekteshi, L., Lazo, P., Qarri, F., Stafilov, T. (2021). Spatial Series and Multivariate Analysis in Assessing the Essential (Cu and Zn) and Toxic (As, Cd, Cr, Co, Hg, Ni and Pb) Metals Linked with Health Risk and Ecological Effects of Atmospheric Deposition by Using Bryophyte Moss as Bioindicator, in: Balabanova, B., Stafilov, T. (Eds.), Contaminant Levels and Ecological Effects: Understanding and Predicting with Chemometric Methods, Springer International Publishing, Cham, pp. 33–74.

  3. Astel, A., Astel, K., Biziuk, M. (2008). PCA and multidimensional visualization techniques united to aid in the bioindication of elements from transplanted Sphagnum palustre moss exposed in the Gdańsk City Area. Environ. Sci. Pollut. Res. Int. 15, 41–50.​espr2007.05.422

  4. Bargagli, F., Monaci, F., Borghini, F., Bravi, F., Agnorelli, C. (2002). Mosses and lichens as biomonitors of trace metals. A comparison study on Hypnum cupressiforme and Parmelia caperata in a former mining district in Italy. Environ. Pollut. 116, 279–287.​10.1016/S0269-7491(01)00125-7

  5. Basile, A., Sorbo, S., Aprile, G.G., Conte, B., Cobianchi, R.C. (2008). Comparison of the heavy metal bioaccumulation capacity of an epiphytic moss and an epiphytic lichen. Environ. Pollut. 151, 401-407.

  6. Basile, A., Sorbo, S., Aprile, G.G., Conte, B., Cobianchi, R.C., Pisani, T., Loppi, S. (2009). Heavy metal deposition in the Italian ‘‘triangle of death’’ determined with the moss Scorpiurum circinatum. Environ. Pollut. 157, 2255–2260.

  7. Berg, T., Røyset, O., Steinnes, E. (1995). Moss (Hylocomium splendens) used as biomonitor of atmospheric trace element deposition: estimation of uptake efficiencies. Atmos. Environ. 29, 353–360.

  8. Berg, T., Steinnes, E. (1997). Use of mosses (Hylocomium splendens and Pleurozium schreberi) as biomonitors of heavy metal deposition: From relative to absolute deposition values. Environ. Pollut. 98, 61–71.

  9. Boquete, M.T., Ares, A., Fernández, J.A., Aboal, J.R. (2020). Matching times: Trying to improve the correlation between heavy metal levels in mosses and bulk deposition. Sci. Total Environ. 715, 136955.

  10. Bourdrel, T., Annesi-Maesano, I., Alahmed, B., Maesano, C.N., Bind, M.A. (2021). The impact of outdoor air pollution on COVID-19: A review of evidence from in vitro, animal, and human studies. Eur. Respir. Rev. 30, 200242.

  11. Briffa, J., Sinagra, E., Blundell, R. (2020). Heavy metal pollution in the environment and their toxicological effects on humans. Heliyon 6, E04691.​e04691

  12. Brunekreef, B., Downward, G., Forastiere, F., Gehring, U., Heederik, D.J.J., Hoek, G., Koopmans, M.P.G., Lam, S., Roel, C. (2020). Air pollution and COVID-19: Including elements of air pollution in rural areas, indoor air pollution, vulnerability and resilience aspects of our society against respiratory disease, social inequality stemming from air pollution. Publications Office of the European Union, LU.

  13. Cheng, I., Mamun, A.A., Zhang, L. (2021). A synthesis review on atmospheric wet deposition of particulate elements: Scavenging ratios, solubility, and flux measurements. Environ. Rev. 29, 340–353

  14. Dales, R., Blanco-Vidal, C., Romero-Meza, R., Schoen, S., Lukina, A., Cakmak, S. (2021). The association between air pollution and COVID-19 related mortality in Santiago, Chile: A daily time series analysis. Environ. Res. 198, 111284.

  15. Donaldson, K., Mac Nee, W. (2001). Potential mechanisms of adverse pulmonary and cardiovascular effects of particulate air pollution (PM10). Int. J. Hyg. Environ. Health 203, 411–415.

  16. Esposito, S., Loppi, S., Monaci, F., Paoli, L., Vannini, A., Sorbo, S., Maresca, V., Fusaro, L., Karam, E.A., Lentini, M., De Lillo, A., Conte, B., Cianciullo, P., Basile, A. (2018). In-field and in-vitro study of the moss Leptodictyum riparium as bioindicator of toxic metal pollution in the aquatic environment: Ultrastructural damage, oxidative stress and HSP70 induction. PLoS One 13, e0195717.

  17. Fang, F., Mu, L., Zhu, Y., Rao, J., Heymann, J., Zhang, Z.F. (2021). Long-term exposure to PM2.5, facemask mandates, stay home orders and COVID-19 incidence in the United States. Int. J Environ. Res. Public Health 18, 6274.

  18. Fattorini, D., Regoli, F. (2020). Role of the chronic air pollution levels in the Covid-19 outbreak risk in Italy. Environ. Pollut. 264, 114732.

  19. Feng, J., Vet, R., Cole, A., Zhang, L., Cheng, I., O’Brien, J., Macdonald A.M. (2021). Inorganic chemical components in precipitation in the eastern U.S. and Eastern Canada during 1989–2016: Temporal and regional trends of wet concentration and wet deposition from the NADP and CAPMoN measurements. Atmos. Environ. 254. 118367.​2021.118367

  20. Frontasyeva, M., Harmens, H., Uzhinskiy, A., Chaligava, O., participants of the moss survey (2020). Mosses as biomonitors of air pollution: 2015/2016 survey on heavy metals, nitrogen and POPs in Europe and beyond. Report of the ICP Vegetation Moss Survey Coordination Centre, Joint Institute for Nuclear Research, Dubna, Russian Federation.

  21. Gjengedal, E., Steinnes, E. (1990). Uptake of metal ions in moss from artificial precipitation. Environ. Monit. Assess. 14, 77–87.

  22. Grodzinska, K., Szarek-Lukaszewska, G. (2001). Response of mosses to the heavy metal deposition in Poland – An overview. Environ. Pollut. 114, 443–451.

  23. Harmens, H., Buse, A., Büker, P., Norris, D., Mills, G., Williams, B., Reynolds, B., Ashenden, T.W., Rühling, Å., Steinnes, E. (2004). Heavy metal concentrations in European mosses: 2000/2001 survey. J. Atmos. Chem. 49, 425–436.

  24. Harmens, H, Mills, G, Hayes, F, Norris, D., the participants of the ICP Vegetation (2011). Air Pollution and Vegetation ICP Vegetation Annual Report 2010/2011. ICP Vegetation Programme Coordination Centre, Centre for Ecology and Hydrology, Bangor, UK.

  25. Harmens, H., Norris, D., Mills, G., the participants of the moss survey (2013). Heavy metals and nitrogen in mosses: Spatial patterns in 2010/2011 and long-term temporal trends in Europe. ICP Vegetation Programme Coordination Centre, Centre for Ecology and Hydrology, Bangor, UK. 

  26. Harmens, H., Norris, D.A., Sharps, K., Mills, G., Alber, R., Aleksiayenak, Y., Blum, O., Cucu-Man, S.M., Dam, M., De Temmerman, L., Ene, A., Fernández, J.A., Martinez-Abaigar, J., Frontasyeva, M., Godzik, B., Jeran, Z., Lazo, P., Leblond, S., Liiv, S., Magnússon, S.H., et al. (2015). Heavy metal and nitrogen concentrations in mosses are declining across Europe whilst some “hotspots” remain in 2010. Environ. Polluti. 200, 93–104.​2015.01.036

  27. Harrison, R.M., Tilling, R., Callén Romero, M.S, Harrad, S., Jarvis, K. (2003). A study of trace metals and polycyclic aromatic hydrocarbons in the roadside environment. Atmos. Environ. 37, 2391–2402.

  28. ICP Vegetation (2010). Heavy metals in European mosses: Survey. Monitoring manual. International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops. (accessed 4 October 2012).

  29. Järup, L. (2003). Hazards of heavy metal contamination. Br. Med. Bull. 68, 167–182.

  30. Kummer, U., Pacyna, J., Pacyna, E., Friedrich, R. (2009). Assessment of heavy metal releases from the use phase of road transport in Europe. Atmos. Environ. 43, 640–647.​10.1016/j.atmosenv.2008.10.007

  31. Lazo, P., Steinnes, E., Qarri, F., Allajbeu, S., Stafilov, T., Frontasyeva, M., Harmens, H. (2018). Origin and spatial distribution of metals in moss samples in Albania: A hotspot of heavy metal contamination in Europe. Chemosphere 190, 337–349.​ere.2017.09.132

  32. Lazo, P., Stafilov, T., Qarri, F., Allajbeu, Sh., Bekteshi, L., Fronasyeva, M., Harmens, H. (2019). Spatial and temporal trend of airborne metal deposition in Albania studied by moss biomonitoring. Ecol. Indic. 101, 1007–1017.

  33. Lazo, P., Qarri, F., Allajbeu, S., Kane, S., Bekteshi, L., Frontasyeva, M., Stafilov, T. (2021). The Evaluation of Air Quality in Albania by Moss Biomonitoring and Metals Atmospheric Deposition. Springer International Publishing, Cham.

  34. Loring, D.H., Rantala, R.T.T. (1992). Manual for the geochemical analyses of marine sediments and suspended particulate matter. Earth Sci. Rev. 32, 235–283.

  35. Macedo-Miranda, G., Avila‑Pérez, P., Gil‑Vargas, P., Zarazúa, G., Sánchez‑Meza, J.C., Zepeda-Gómez, C., Tejeda, S. (2016). Accumulation of heavy metals in mosses: A biomonitoring study. SpringerPlus 5, 715.

  36. Mamun, A.A., Cheng, I., Zhang, L., Dabek-Zlotorzynska, E., Charland, J.P. (2020). Overview of size distribution, concentration, and dry deposition of airborne particulate elements measured worldwide. Environ. Rev. 28, 77–88.

  37. Meyer, M., Pesch, R., Schröder, W., Steinnes, E., Uggerud, H.T. (2014). Spatial patterns and temporal trends of heavy metal concentrations in moss and surface soil specimens collected in Norway between 1990 and 2010. Environ. Sci. Eur. 26, 27.

  38. Pacyna, J.M., Pacyna, E.G. (2001). An assessment of global and regional emissions of trace metals to the atmosphere from anthropogenic sources worldwide. Environ. Rev. 9, 269–298.

  39. Peck, P. (2004). Reducing Environment & Security Risks from Mining in South Eastern Europe: Desk-assessment study for the Environment and Security Initiative Project. UNEP Division of Technology.

  40. Pozzer, A., Dominici, F., Haines, A., Witt, C., Münzel, T., Lelieveld., J. (2020). Regional and global contributions of air pollution to risk of death from COVID-19. Cardiovasc. Res. 116, 2247–2253.

  41. Reimann, C., Filzmoser, P., Garrett, R. (2002). Factor Analysis applied to regional geochemical data: Problems and possibilities. Appl. Geochem. 17, 185–206.​S0883-2927(01)00066-X

  42. Richards, J.R. (2020). Control of Particulate Matter Emissions. APTI Course 413 Third Edition. ICES Ltd. EPA Contract No. 68D99022.

  43. Rudnick, R.L., Gao, S. (2003). 3.01 - Composition of the Continental Crust, in: Holland, H.D., Turekian, K.K. (Eds.), Treatise on Geochemistry, Pergamon, Oxford, pp. 1–64.​10.1016/B0-08-043751-6/03016-4

  44. Rühling, A. (1994). Atmospheric Heavy Metal Deposition in Europe—Estimation Based on Moss Analysis; NORD 1994:9; Nordic Council of Ministers: Copenhagen, Denmark.

  45. Rühling, Å., Steinnes, E. (1998). Atmospheric heavy metal deposition in Europe 1995-1996. NORD 1998:15. Nordic Council of Ministers, Copenhagen.

  46. Schröder, W., Nickel, S., Schönrock, S., Meyer, M., Wosniok, W., Harmens, H., Frontasyeva, M.V., Alber, R., Aleksiayenak, J., Barandovski, L., Carballeira, A., Danielsson, H., de Temmermann, L., Godzik, B., Jeran, Z., Karlsson, G.P., Lazo, P., Leblond, S., Lindroos, A.J., Liiv, S., et al. (2016). Spatially valid data of atmospheric deposition of heavy metals and nitrogen derived by moss surveys for pollution risk assessments of ecosystems. Environ. Sci. Pollut. Res. 23, 10457–10476.

  47. Schröder, W., Nickel, S. (2019). Moss species-specific accumulation of atmospheric deposition? Environ. Sci. Eur. 31, 78.

  48. Shaman, D., Wheeler, D. (2007). Controlling industrial pollution in the developing world. Environ. Qual. Manage. 8, 69–75.

  49. Sokhi, R.S., Moussiopoulos, N., Baklanov, A., Bartzis, J., Coll, I., Finardi, S., Friedrich, R., Geels, C., Grönholm, T., Halenka, T., Ketzel, M., Maragkidou, A., Matthias, V., Moldanova, J., Ntziachristos, L., Schäfer, K., Suppan, P., Tsegas, G., Carmichael, G., Franco, V., et al. (2022). Advances in air quality research – current and emerging challenges. Atmos. Chem. Phys. 22, 4615–4703.

  50. Stafilov, T., Šajn, R., Barandovski, L., Andonovska, K.B., Malinovska, S. (2018). Moss biomonitoring of atmospheric deposition study of minor and trace elements in Macedonia. Air Qual. Atmos. Health 11, 137–152.

  51. Steinnes, E. (1989). Biomonitors of Air Pollution by Heavy Metals, in: Pacyna, J.M., Ottar, B. (Eds.), Control and Fate of Atmospheric Trace Metals, Springer Netherlands, Dordrecht, pp. 321–338.

  52. Steinnes, E., Rühling, Å., Lippo, H., Makinen, H. (1997). Reference materials for large-scale metal deposition survey. Accred. Qual. Assur. 2, 243–249.

  53. Tchounwou, P.B., Yedjou, C.G., Patlolla, A.K., Sutton, D.J. (2012). Heavy metals toxicity and the environment. NIH Public Access. 101, 133–164.

  54. Tyler, G. (1990). Bryophytes and heavy metals: A literature review. Bot. J. Limn. Soc. 104, 231–253.

  55. U.S. Environmental Protection Agency (U.S. EPA) (2004). Air Quality Criteria for Particulate Matter (Final Report, 2004). U.S. Environmental Protection Agency, Washington, DC, EPA 600/P-99/002aF-bF. Vol. I, pp. 2–11.

  56. United Nations Environment Programme (UNEP) (2013). Environmental Risks and Challenges of Anthropogenic Metals Flows and Cycles, A Report of the Working Group on the Global Metal Flows to the International Resource Panel. Van Der Voet, E., Salminen, R., Eckelman, M., Mudd, G., Norgate, T., Hischier, R. Copyright © United Nations Environment Programme.

  57. World Meteorological Organization (WMO) (2019). Atmospheric Deposition.

  58. Wright, L.P., Zhang, L., Cheng, I., Aherne, J., Wentworth, G.R. (2018). Impacts and effects indicators of atmospheric deposition of major pollutants to various ecosystems - A review. Aerosol Air Qual. Res. 18, 1953–1992.

  59. Zechmeister, H.G., Grodzińska, K., Szarek-Łukaszewska, G. (2003). Chapter 10 Bryophytes, in: Markert, B.A., Breure, A.M., Zechmeister, H.G. (Eds.), Bioindicators & Biomonitors: Principles, Concepts and Applications. Elsevier, Amsterdam.

  60. Zhang, Z., Xue, T., Jin, X. (2020). Effects of meteorological conditions and air pollution on COVID-19 transmission: Evidence from Chinese cities. Sci. Total Environ. 741, 140244.​10.1016/j.scitotEnviron.2020.140244

Share this article with your colleagues 


Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

81st percentile
Powered by

Aerosol and Air Quality Research partners with Publons

Special Call for Papers Air Pollution and its impact in South and Southeast Asia

2020 Impact Factor: 3.063
5-Year Impact Factor: 2.857

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal that promotes submissions of high-quality research and strives to be one of the leading aerosol and air quality open-access journals in the world. We use cookies on this website to personalize content to improve your user experience and analyze our traffic. By using this site you agree to its use of cookies.