Introduction

The Antarctic region is still the least polluted area on our planet. It is isolated by ocean, cyclonic storm belts, and hard weather conditions (Shaw 1988). However, since the arrival of first explorers in the nineteenth century and increasing human activity (scientific activity and tourism) and threaten environment of this region (Osyczka et al. 2007; Chwedorzewska and Korczak 2010), as already reported, the main sources of pollution in Antarctic area are burning of fuels and waste storing (Kabata-Pendias 2000; Bargagli 2008). According to increasing human activity, it is necessary to monitor pollution concentration and to control this unique environment. The King George Island is an example of Antarctic area with high-human activity. Ten polar stations (eight all-years and two summers only) are located on the island (Osyczka et al. 2007). They belong to Argentina, Chile, Brazil, China, Ecuador, South Korea, Peru, Poland, Russia, Uruguay, and USA. The King George Island is one of the South Shetlands Island group (Fig. 1). The mean annual temperatures are in range from − 1.7 to 2.4 °C, and the ice-free area of the island is about 8% (Kejna 1999; Bӧtler 2011). Soils on King George Island are mainly derived from volcanic rocks like basalts and andesite basalts, but in few sites also sedimentary rocks are present, and the soils formed on these rocks are cryosols, leptosols, regosols, and fluvisols. In this environment, mosses and lichens are two of the most important groups of organisms. Mosses and lichens are protecting initial soils from weather conditions and thus forming basic environments, which can be later populated by microbes and lower organisms (Tatur and Myrcha 1984; Bӧtler 2011). Since 1970s mosses and lichens are used as indicators of environmental pollution, especially air pollution (Sloof 1993; Steiness 1995; Szczepaniak and Biziuk 2003; Wu et al. 2014). Mosses and lichens have no roots unlike higher plants. Therefore, their basic nutrient source is atmospheric elements deposition. Because these groups of organisms have no protective cuticle, they can absorb ions from the air, rainwater, or snow. The growth rate of mosses and lichens is slow, so they can accumulate pollutants very effectively (Nash 1996; Turetsky et al. 2012). In this study, lichens Usnea antarctica (Du Rietz), U. aurantiaco-atra (Jacq.), and moss Sanionia uncinata (Hedw.) were used as bioindicators of concentration of heavy metals, rare earth elements, and trace and major elements in the air. Analysis of samples was carried out using Instrumental Neutron Activation Analysis (INAA). INAA is a non-destructive, multi-elemental method that allows simultaneous determination of about 50 elements in a single sample with mass about 300 mg (Bode 1996; Frontasyeva 2011). The advantages of INAA are low-detection limits, very precise results, and fast-sample preparation without complicated chemical treatment, thus INAA is widely applied in environmental pollution studies, for example, in the Intercontinental Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops (ICP Vegetation) research (Kłos et al. 2013; Thinova et al. 2014; Allajbeu et al. 2016; Harmens et al. 2016). The goal of this study was to estimate local environmental pollution levels and identify possible pollution sources including not only anthropogenic but also natural like, e.g., volcanic or biogenic activity, and analyze differences between trace element accumulation by mosses and lichens including investigation of rare earth elements distribution. Additionally, scanning electron microscopy (SEM) was applied to determine microparticles and presence of cosmic dust in moss and lichen samples. Cosmic dust is very abundant on earth (Brownlee 1985). The total mass of dust inside earth orbit is about 1015 tons, and every year 4∙103–104 tons is deposited on earth surface (Yada et al. 2004; Grachev et al. 2008) that means that cosmic dust deposition should be taken account during environmental trace analysis (Yada et al. 2004). Investigations of cosmic dust were started by the HMS Challenger expedition, during which traces of cosmic matter were discovered in sea sediments (Murray and Renard 1891). Later, metallic microspherules have been found also in Antarctic and Greenland glaciers and on deserts (Maurette et al. 1986; Yada et al. 2004). Cosmic dust can be subdivided into interplanetary cosmic dust particles (IDP) or micrometeorites depending to size of particles (IDP < 30 μm and micrometeorites > 50 μm), but there are no existing rigorous criteria for discrimination between cosmic dust and micrometeorites yet (Genge et al. 1997; Grachev et al. 2008). Cosmic dust particles may have very different shapes, so not only morphological criteria (spherical shape, textures, or metallic luster) but also their chemical composition is very important for confirmation of their cosmic origin. For example, volcanic-origin Fe microspheres have high-Ti content (more than 10%) and rarely have perfect spherical shape (Szӧőr et al. 2001; Grachev et al. 2008). There is still no clear origin of cosmic dust. It can be remained by primordial matter from protoplanetary cloud, or it can be produced by comets and asteroids destruction (Grachev et al. 2008).

Fig. 1
figure 1

Map of King George Island with sampling sites (numbers)

Materials and methods

Moss and lichen samples were collected in summer seasons 2002–2003 by Polish Antarctic Expeditions of the Polish Academy of Sciences to the Henryk Arctowski Station. Sampling sites were free from direct influence of visible aerosols. Four sampling sites (1, 2, 5, and 8) were located close to polar stations. Sampling site 1 was located near Marsch Airfield, close to Chilean, Eduardo Frei base with average winter population of 80 people. Samples from sampling site number 2 were collected close to Russian Bellingshausen base with average winter population of 25 people. Sampling site number 5 was located near Brazilian station Comandante Ferraz with 12 people winter crew, and site number 8 was located close to Peruvian Macchu Picchu station (summer only station). S. uncinata (Hedw.) is one of the most common mosses in the Antarctic region (Putzke et al. 2015). It forms large area carpets, associated with other mosses and lichens. U. antarctica (Du Rietz) and U. aurantiacoatra (Jacq.) are lichens mostly associated with S. uncinata (Victoria et al. 2009; Victoria et al. 2013). Sampling sites (Table 1) were located within Antarctic stations operating in full-year or summer only cycle. Collected material was packed into plastic bags. After identification, samples were dried and shipped to Poland. In 2015, samples were analyzed using INAA in Frank Laboratory of Neutron Physics, Joint Institute for Nuclear Research in Dubna. Samples were once again dried in 30 °C and cleaned mechanically from possible soil impurities using scissors and tweezers. From each sample, two portions about 300 mg were taken to form two tablets for short and long irradiation. The mass of 300 mg is a sufficient mass to use moss and lichen samples without homogenization (Korzewka et al. 2007; Steinnes et al. 2007). Tablets were packed in polyethylene bags and plastic containers (for short irradiation) and in aluminum foil and aluminum container for long irradiation. Samples were irradiated in IBR-2M reactor in the Frank Laboratory of Neutron Physics. For determination using short-lived isotopes (Mg, Al, Cl, Ca, V, Mn, I, Cu, K, S, Ti, In, Dy, Si), samples were irradiated for 3 min. After short irradiation in channel 2 (without Cd screen, Table 2), samples were measured twice. First measurement was for 5 min after 3 min of decay and second one for 20 min after 10 min of decay. For determination of long-lived isotopes (Na, K, Sc, Cr, Fe, Co, Ni, Sb, Zn, As, Rb, Sr, Cs, Ba, La, Sm, W, Th, U, Cs, Au, Eu, Nd, Zr, Rb, Br, Se, Cd, Ag, Mo, Ce, Hg, W, Ta, Hf, Lu, Yb, Tm, Tb, Gd), samples were irradiated by 100 h in Cd-screened channel 1. After irradiation, samples were repacked and measured for 45 min and for 3 h after 5 and 20 days, respectively. Measurements of gamma rays were performed using High-Purity Germanium (HPGe) detector with resolution 1.9 keV on 1332-keV line from 60Co. Data acquisition process was controlled using Genie 2000 software (Canberra). The processing of data and determination of element concentration in samples were performed using software developed in FLNP JINR (Ostrovnaya et al. 1993; Ostrovnaya 2000). For quality assurance purposes, certified reference materials (IAEA 336) and neutron flux comparators were used.

Table 1 Sample description
Table 2 Neutron flux parameters of irradiation channels

SEM microscopy and MPA analysis were performed using Tescan Vega II (Tescan, Czech Republic) with energy dispersive X-ray analyzer (EDS). Samples were analyzed with an accelerating voltage of 20 kV and beam current 0.2 nA. Samples of mosses and lichens with mass about 300 mg were milled in the agate mortar and sonicated. After homogenization, magnetic particles were extracted by permanent magnet. Magnetic particles collected on the magnet were transferred to double layer conductive carbon adhesive tape and placed on objective table using glass rod. Data acquisition time was 1–2 min, and obtained results were normalized to 100% (Pechersky et al. 2015a, b). Total number of 53 objects found in 14 samples of mosses and lichens were analyzed. Depending on size of object, 3 to 10 points were analyzed by EDS spectrometer.

Principal component analysis (PCA) was used for identification possible sources of elements in analyzed samples. For PCA analysis, we used Varimax rotation with Kaiser normalization (eigenvalues > 1) and p = 0.05.

Normalization of REE concentrations was done by using chondrite values following by Taylor and McLennan (1985).

Results and discussion

Tables 3 and 4 show ranges and medians for 52 element concentrations determined in moss and lichen samples. The concentrations of elements in investigated samples were different depending on sampling site and species. In general, concentrations of elements in moss samples were greater than in lichens. This observation can be explained by fact that S. uncinata can use rhizoids to ensure some part of water supply. Elements dissolved in soil-water may then be transported from soil to moss (Osyczka et al. 2007). Only for Cl, Ag, I, Lu, and W concentrations were greater in lichens than in mosses. The excess of I and Cl in lichens as compared to mosses can be probably an effect of organohalogen synthesis by fungal part of lichen (Matschullat et al. 1999; Gribble 2010).

Table 3 Results of lichens samples analysis
Table 4 Results of moss samples analysis

Rare earth elements

In this group, concentrations of 12 elements (Sc, La, Ce, Nd, Sm, Eu, Gd, Tb, Dy, Tm, Yb, Lu) were measured. Rare earth elements (REE) are naturally occurring ingredients of lithosphere commonly used for tracing geochemical processes, but also are widely applied in industry (e.g., electronics or medicine), and therefore anthropogenic contamination with REE is possible (Brioschi et al. 2013; Allajbeu et al. 2016). For lichens, the highest concentrations of REE were usually found on Red Hill and the lowest on Varueal Peak. In case of mosses, the highest concentrations were measured in samples collected near Ferraz Station and the lowest in samples from Blue Dyke. For RRE, chondrite normalization (Taylor and McLennan 1985) was applied and geochemical parameters were calculated to obtain data about their fractionation, and results are showed in Table 5.

Table 5 Geochemical parameters for rare earth elements (RRE) and chondrite normalized concentrations in samples and upper continental crust (UCC)

In Table 5, upper continental crust (UCC) is a reference composition (Rudnick and Gao 2004). Values of all REE concentrations are lower than reference values for UCC. For geochemical parameters, there were no significant differences between mosses and lichens, except La/Th and Tm/Tb ratios, but in these cases, standard deviation in lichens was high. That means that both groups of organisms are accumulating REE similarly. Both groups show fractionation of REE. For example, δEu is more than two times higher than reference value for UCC. Eu anomalies were calculated using equation proposed by Taylor and McLennan (1995):

$$ \updelta \mathrm{Eu}=\frac{{\mathrm{Eu}}_{\mathrm{n}}}{\sqrt{\left({\left[\mathrm{Sm}\right]}_{\mathrm{n}}\times {\left[\mathrm{Gd}\right]}_{\mathrm{n}}\right)}} $$
(1)

where Eun, [Sm]n, and [Gd]n are chondrite normalized concentrations of europium, samarium, and gadolinium. Similar observations were reported in other studies (Ryghaug 1983; Aubert et al. 2006) and show positive anomaly for heavy REE (Allajbeu et al. 2016). All of calculated parameters show REE distribution patterns different than UCC. The La/Th ratio higher than UCC can be evidence for sedimentary rock-originated REE’s deposited on mosses and lichens (McLennan et al. 1993).

Heavy metals

To asses King George contamination, 15 heavy metals (V, Cr, Mn, Fe, Ni, Co, Cu, Zn, Se, As, Mo, Cd, Sb, W, and Hg) were measured. For lichen samples, most contaminated area was Marsh Airfield. Only three elements (Cd, Mo, and Se) had maximal concentrations on Red Hill. In moss samples, it was difficult to choose one most contaminated site. Highest concentrations of Mn, Fe, Co, Zn, Sb, and V were founded near Bellingshausen station. Marsh Airfield had the highest Cd concentration. Maximal concentrations of Mo, As, and Se were measured in samples collected near Machu Picchu base. Varueal Peak was a sampling site with highest Cu concentration and Lions Ramp with Cr and Ni. Median values of heavy metal concentration in lichen and moss samples from King George Island can be compared with values from Norway (Steinnes 2005; Barandovski et al. 2006), reference plant model (Markert 1991), and literature data for Antarctic region (Bargagli et al. 1999; Bargagli et al. 2000; Gonzáles et al. 2002; Smykla et al. 2005; Lim et al. 2009; Osyczka et al. 2007; Zvěřina et al. 2014; Bubach et al. 2016) (Tables 6 and 7).

Table 6 Concentrations of heavy metals and Se compared with literature data, reference plant (RP), and Norway moss survey (Bargagli et al. 1999; Bargagli et al. 2000; Gonzáles et al. 2002; Ganeva and Yurukova 2004; Smykla et al. 2005; Barandovski et al. 2006; Osycza et al. 2007; Lim et al. 2009; Mão de Ferro et al. 2013; Zvěřina et al. 2014; Amaro et al. 2015; Bubach et al. 2016)
Table 7 Factor analysis of NAA data on moss and lichen samples (Rotation method: Varimax with Kaiser normalization)

Obtained results are general in the range of literature data. In moss samples, however, concentrations of Cr and Ni are elevated. The same situation can be observed for Se in lichen samples. In comparison to values from Norway, concentrations of almost all heavy metals are elevated; only depositions of W and Zn are lower.

Sources of pollutants

There are many possible sources of contamination of investigated area. The possible sources can be elevated natural background, local contamination from human activity, or long-range atmospheric transport. It was reported that As from Chilean Cu mines was found in Antarctic ice cores. Also, important source of As and Se is combustion of coal with arsenopyrite. The natural background can be also elevated due to volcanic activity (Sieprawska et al. 2015; Schwanck et al. 2016). On the other hand, Amouroux et al. (2001) described mechanism of releasing Se from oceans. In effect of biogenic activity, volatile compounds as dimethyl sulfide (DMS) and dimethyl selenide (DMSe) can be released from the ocean and later be deposited on a land surface. Oceans can be also source of Cd (Bargagli et al. 1998). For some element deposition, also local microclimate can be important and must be taken into account (Zvěřina et al. 2014). For establishing potential sources of pollution, statistical methods can be used. In our work, factor analysis was used. Based on this approach, four factors can be distinguished as an origin of pollutants in lichen and moss samples.

Factor 1 represents typical crustal composition and is associated with resuspension of soil and rock particles. Factor 2 can be probably associated with volcanic activity. Antarctic Peninsula has extensive volcanic system, and Deception Island is one of active volcanoes. Soils of King George Island have high content of tephra and ash (Deheyn et al. 2005), which probably are resuspended similar to factor 1. Factors 3 and 4 are probably connected with human activity on island and long-range transport of pollutants. There is also a possibility that some of pollutants in these factors are delivered by birds eating contaminated fishes. Pollutants could be then deposited on island with feces. That situation can be observed near Admiralty Bay where penguin colonies and ornithogenic soils are present, and arsenic concentrations were in general higher than in other sampling sites, and arsenic enriches in penguin dropping sediment (Xie and Sun 2008). Factor 4 includes Mo and As. Both of them can be local pollutants or transported to the island (Schwanck et al. 2016), where Mo may be reduced to insoluble form by bacterial activity (Ahmad et al. 2013). Chromium and Ni can be connected with combustion of diesel oil or with construction materials like stainless steel (Barałkiewicz and Siepak 1999; Kotaś and Stasicka 1999).

SEM results

SEM microscopy images of investigated samples of moss and lichen showed presence of volcano and cosmic-originated particles in all samples (Fig. 2). First six samples (1–6) are probably volcanic-originated. Sample 1 is titanomagnetite monocrystal. Number 2 is also titanomagnetite but with detrital structure. Number 3 is a native Al, with 5% Mg impurity. Sample number 4 is a pyrite (Fe 45%, S 48%, and O 7% by mass). Next sample (number 5) is a pyroxene and titanomagnetite. Large object in center of image is a pyroxene (44% O, 9% Mg, 2% Al, 23% Si, 13% Ca, and 10% Fe by mass), and small bright spots are titanomagnetites (35% O, 2% Al, 1% Si, 6% Ti, and 56% Fe); sample 6 is native iron. Samples 7, 9, 10, and 12 are spherical object space origin classified due to their dimensions as IDP and with chemical composition of Fe and O, without Ti. Samples 8 and 11 have differed chemical composition. In sample 8, we have 6% O, 16% Cr, and 78% of Fe and its chemical composition can be evidence for cosmic origin (Korchagin et al. 2010). Finally, sample 11 is made from Ti (2%), Ni (4%), Cu (49%), Zn (17%), and W (28%) and is a rare example of cosmic particle with Cu-Zn alloy (Korchagin et al. 2010) where Ti can be contaminated from basalts (Pechersky et al. 2015a, b). The spherical IDPs found in our samples are type I microspherules (Brownlee et al. 1997) and are composed mainly from Fe3O4 and FeO (Engrand et al. 2005). As the average flux of cosmic dust is 1/m2/d for 10 μm particles and 1/m2/y for 100 μm (McDonnell 1978; McDonnell et al. 1984), total annual deposition of cosmic particles on King George Island (1150 km2) can be estimated as 4.2∙1011 10 μm particles and 1.15∙109 100 μm particles.

Fig. 2
figure 2

SEM photos of earth-origin (objects 1–6) and space-origin (objects 6–12) microparticles

Summary

In this study, 50 element concentrations were determined using instrumental neutron activation analysis in moss and lichen samples. Based on INAA and statistical results (PCA analysis), a few possible sources of pollution of King George Island were suggested. Also, presence of extraterrestrial material in investigated samples was confirmed by SEM microscopy analysis.