Moss biomonitoring technique was used for a heavy-metal pollution study in Macedonia in the framework of the International Cooperative Program on Effects of Air Pollution on Natural Vegetation and Crops (UNECE IPC Vegetation). Moss samples (n = 72) were collected during the summers of 2002, 2005, and 2010. The contents of 41 elements were determined by neutron activation analysis, atomic absorption spectrometry, and inductively coupled plasma atomic emission spectrometry. Using factor and cluster analyses, three geogenic factors were determined (Factor 1, including Al, As, Co, Cs, Fe, Hf, Na, Rb, Sc, Ta, Th, Ti, U, V, Zr, and rare-earth elements–RE; Factor 4 with Ba, K, and Sr; and Factor 5 with Br and I), one anthropogenic factor (Factor 2, including Cd, Pb, Sb, and Zn), and one geogenic-anthropogenic factor (Factor 3, including Cr and Ni). The highest anthropogenic impact of heavy metal to the air pollution in the country was from the ferronickel smelter near Kavadraci (Ni and Cr), the lead and zinc mines in the vicinity of Makedonska Kamenica, Probištip, and Kriva Palanka in the eastern part of the country (Cd, Pb, and Zn), and the former lead and zinc smelter plant in Veles. Beside the anthropogenic influences, the lithology and the composition of the soil also play an important role in the distribution of the elements.

This paper aims to investigate the correlations between the concentrations of nine heavy metals in moss and atmospheric deposition within ecological land classes covering Europe. Additionally, it is examined to what extent the statistical relations are affected by the land use around the moss sampling sites. Based on moss data collected in 2010/2011 throughout Europe and data on total atmospheric deposition modelled by two chemical transport models (EMEP MSC-E, LOTOS-EUROS), correlation coefficients between concentrations of heavy metals in moss and in modelled atmospheric deposition were specified for spatial subsamples defined by ecological land classes of Europe (ELCE) as a spatial reference system. Linear discriminant analysis (LDA) and logistic regression (LR) were then used to separate moss sampling sites regarding their contribution to the strength of correlation considering the areal percentage of urban, agricultural and forestry land use around the sampling location. After verification LDA models by LR, LDA models were used to transform spatial information on the land use to maps of potential correlation levels, applicable for future network planning in the European Moss Survey. Correlations between concentrations of heavy metals in moss and in modelled atmospheric deposition were found to be specific for elements and ELCE units. Land use around the sampling sites mainly influences the correlation level. Small radiuses around the sampling sites examined (5 km) are more relevant for Cd, Cu, Ni, and Zn, while the areal percentage of urban and agricultural land use within large radiuses (75–100 km) is more relevant for As, Cr, Hg, Pb, and V. Most valid LDA models pattern with error rates of < 40% were found for As, Cr, Cu, Hg, Pb, and V. Land use-dependent predictions of spatial patterns split up Europe into investigation areas revealing potentially high (= above-average) or low (= below-average) correlation coefficients. LDA is an eligible method identifying and ranking boundary conditions of correlations between atmospheric deposition and respective concentrations of heavy metals in moss and related mapping considering the influence of the land use around moss sampling sites.

In the summer of 2005 and 2010, moss samples were collected from 72 sampling sites evenly distributed all over the territory of Republic of North Macedonia. Kjeldahl method was used to determine the nitrogen content in the samples. Descriptive statistics and distribution maps were prepared. Data obtained from these two surveys were compared, and additional comparison was done with data obtained from similar studies in the South-Eastern European countries and Finland as a clean area. The median value of N content in the samples collected in 2005 is 1.21%, varies from 0.70% to 1.54%, while the content of N in samples collected in 2010 ranges between 0.68% and 1.75% with the median value of 1.06%. High contents of N were found in the Northern and Central parts of the country mainly as a result of agricultural activities, industry and traffic.

The moss biomonitoring technique was used in 2002, 2005, 2010 and 2015 in a potentially toxic elements study (PTEs) in Macedonia. For that purpose, more than 70 moss samples from two dominant species (Hypnum cupressiforme and Homalothecium lutescens) were collected during the summers of the mentioned years. Total digestion of the samples was done using a microwave digestion system, whilst mercury was analyzed by cold vapour atomic absorption spectrometry (CV–AAS). The content of mercury ranged from 0.018 mg/kg to 0.26 mg/kg in 2002, from 0.010 mg/kg to 0.42 mg/kg in 2005, from 0.010 mg/kg to 0.60 mg/kg in 2010 and from 0.020 mg/kg to 0.25 mg/kg in 2015. Analysis of the median values shows the increase of the content in the period 2002–2010 and a slight reduction of the air pollution with Hg in the period 2010–2015. Mercury distribution maps show that sites with increased concentrations of mercury in moss are likely impacted by anthropogenic pollution. The results were compared to similar studies done during the same years in neighboring countries and in Norway—which is a pristine area and serves as a reference, and it was concluded that mercury air pollution in Macedonia is significant primarily in industrialized regions.