Airborne radionuclides and heavy metals in high Arctic terrestrial environment as the indicators of sources and transfers of contamination

A survey of airborne radioactive isotopes (137Cs, 238Pu, 239+240Pu, 241Am, and 210Pb) and trace metals (Pb, Cu, Zn, Cd, Fe, Al) in tundra soils and cryoconite hole material sampled from several locations in the Kaffiøyra region of Spitsbergen revealed significant variability in spatial concentration. Lithogenic radionuclides (230Th, 232Th, 234U, 238U) show less variability than the airborne radionuclides because their activity concentrations are controlled by mixing of local material derived from different types of bedrock. Activity ratios of the artificial radionuclides in most cryoconite samples differ from global fallout signatures. The contribution of radionuclides from additional and more specific sources might be enhanced by non-continuous exposure of cryoconite to atmospheric deposition. We assumed that the main source of Pu, which was detected only in cryoconite samples, is derived from nuclear tests and nonexploded weapons-grade material. Approximately one-third of the total observed Pu activity concentration is 238Pu, most likely originating from the SNAP9A satellite re-entry, which was powered by a 238Pu thermoelectric generator. In samples from Waldemarbreen the influence of glacial local morphology on the capability of cryoconite for trapping and accumulating airborne radionuclides is apparent. Local glacial morphology plays an important role in determining the accumulation of airborne pollutants. Trace metal concentrations in soils were typical or slightly higher than concentrations characteristic for the natural background; the 206Pb/207Pb signature was also close to the natural ratio of the parent rocks. Conversely, trace metal concentrations in cryoconite samples (Pb and Cd) were higher than in soil samples and exceeded natural values. Cryoconite is an effective monitor of the spread of artificial radionuclides and heavy metals in their surrounding environment.


Introduction
The radioactive nuclides in the environment originate from both natural and anthropogenic sources.The present work is E. Łokas et al.: Airborne radionuclides and heavy metals in high Arctic terrestrial environment thropogenic radionuclides in the Northern Hemisphere are (i) nuclear weapon tests and atmospheric explosions (e.g.Novaya Zemlya, Semipalatinsk, Nevada), (ii) nuclear accidents (e.g.Kyshtym 1957, Lake Karachay 1968, Tomsk 1993, Chernobyl 1986, Fukushima 2011), and (iii) the disintegration of satellites powered by nuclear thermoelectric generators (e.g.SNAP9A-1964, Cosmos 958-1978).Relatively low but constant releases of artificial radioisotopes are also associated with the reprocessing of spent nuclear fuel and energy generation in nuclear power plants.Heavy metals are natural elements of the Earth's crust.Since the beginning of the industrial era (from ∼ 1850) the deposition of metals in the environment due to human activities has increased nearly 10-fold (Nriagu, 1996).The emission of heavy metals increased substantially after World War II and is still increasing in some areas.The primary sources of anthropogenic heavy metals include industry, mining, agriculture, fuel burning, waste disposal, and transportation (Pacyna and Pacyna, 2001).Artificial radionuclides and heavy metals are transported across the world mainly via the atmosphere, oceanic currents, and rivers.These materials can be found in remote Earth locations like the polar cryosphere.Elevated radionuclide and heavy metal levels have been reported in Asian and European glaciers (Barbante et al., 2004;Aizen et al., 2009;Eyrikh et al., 2017;Baccollo et al., 2017) and in remote Arctic and Antarctic areas (Hur et al., 2007;Singh et al., 2013).
The most significant quantities of these pollutants on glaciers are most likely stored in cryoconite granules and microfauna (Segawa et al., 2013;Łokas et al., 2016;Baccolo et al., 2017;Zawierucha et al., 2018).Cryoconite granules are aggregates of mineral and organic components, mainly composed of archaea, algae, cyanobacteria, fungi, and heterotrophic bacteria (Cook et al., 2015;Simon and Bouville, 2002;Takeuchi et al., 2001a, b;Hodson et al., 2008).Cyanobacteria play a crucial role in the formation of these aggregates as they produce extracellular polymeric substances (EPSs) whose adhesive properties enhance the accumulation of mineral dust and microorganisms (Gadd, 2004;Francis, 2007).Cryoconite holes are formed on the ablation zone of glaciers and constitute nutrient-rich pools with diverse biota within the supraglacial zone (Wharton et al., 1985;Cook et al., 2016).Dark-coloured cryoconite significantly increases the absorption of solar light by glaciers, enhancing and accelerating melting processes (Simon and Bouville, 2002;Takeuchi et al., 2001a, b;Hodson et al., 2008).Glaciers in Canada, the Russian Arctic, Alaska, Antarctica, and Greenland are projected to lose about 40 % of their volume by the end of 2100; central Europe, low-latitude South America, Caucasus, northern Asia, western Canada, and the US are expected to lose about 80 % of their volume or entirely disappear within decades at current climatic conditions (Radic et al., 2014).Cryoconite granules contain high concentrations of contaminants and can persist on glacier surfaces for several decades (Łokas et al., 2014).As glacier ice disappears, their contents are deposited on freshly ex-posed areas.The discovery of hotspots of radioactivity in the proglacial zones of Spitsbergen glaciers (Łokas et al., 2014, 2017b) provided evidence that cryoconite-derived contaminants were dispersed by meltwater channels and accumulated in depressions of the glacier fore field.The monitoring and fate of anthropogenic radionuclides and heavy metal distribution are important as both pose a potential threat to the environment.During intense cryosphere melting, contaminants accumulated over many decades on glacier surfaces are released and can pollute local vegetation, surface waters, and even drinking water sources (Dong et al., 2017).Additionally, the determination of Pu isotopes, 241 Am, and 137 Cs in cryoconite has great potential for monitoring the environmental impact of military nuclear activities, the rising number of nuclear events in civilian nuclear installations, and the growing risk of nuclear terrorism and undeclared nuclear activities.

Study area
The Kaffiøyra region (Fig. 1), located in northwestern Spitsbergen, the adjoining Aavatsmarkbreen (75 km 2 ), and the Dahlbreen (132 km 2 ) together comprise an area of about 310 km 2 (Sobota et al., 2016).It is a coastal lowland at Forlandsundet (Forland Sound) and accounts for 12 % of the area of Oscar II Land.Along with Kaffiøyra, which is 14 km long and 4 km wide, seven land glaciers are located in this region.The main pedogenetic processes in the Kaffiøyra region consist of the accumulation of organic material, gleying, and decarbonation.Cryogenic processes, which overlap the pedogenetic processes, also play an important role in the genesis of the soils.The substantial moisture content of the surface layers and stagnant waters result in a shallow permafrost (Plichta, 2005).During 1996-2015, the greatest average thickness of active permafrost (215 cm) was recorded at a location in the moraine.A much lower mean value was recorded on the beach (126 cm).An active layer in the tundra reached on average a depth of 150 cm (Sobota et al., 2018).Climate research in the Kaffiøyra region has been conducted since 1975.The increase in air temperature noted during the study confirms the general trend observed in Svalbard and in many areas of the Arctic (Nørdli et al., 2014).In the years 1997-2016 the average air temperature in the region during the summer season was 5.4 • C.
Glaciers of this region are polythermal (Sobota, 2009(Sobota, , 2011)).Since their maximum extent in the late 19th and early 20th centuries to 2015, the total area of this region's valley glaciers has decreased on average by about 43 % (Sobota et al., 2016).Waldemarbreen is an alpine-type glacier running down the valley towards Kaffiøyra.It contains one accumulation zone and a glacier tongue in the valley.The area of the glacier is 2.40 km 2 (Sobota, 2017).It borders the Prins Heinrichfjella ridge in the north and east (500-770 m above sea level), and Gråfjellet (300-350 m) in the south.The glacier consists of two distinct parts separated by a medial moraine.The mean annual mass balance of the Waldemarbreen in 1996-2015 was −0.72 m w.e.(Sobota et al., 2016), with very significant ablation in 2014.Since 1909, the glacier has been receding by 8 m a −1 .Between 1995 and 2009 the recession rate accelerated to 10 m a −1 , and between 2000 and 2009 reached 11 m a −1 (Sobota and Lankauf, 2010).From the time of maximum advance to 2015 the Waldemarbreen area has decreased by 31.4 % with an average annual recession of 8 m a −1 (Sobota et al., 2016).

Field campaign
Sample locations are shown in Fig. 1.Soil samples were collected in 2009 (coded S01-S06).Sample sites represent accumulations of typical tundra.These areas are predominantly covered with dry moss and lichen tundra with communities of Saxifraga oppositifolia, Salix polaris, Sanionia uncinata, Luzula arcuata ssp.confusa, Cetrariella delisei, and Scorpidium revolvens (Sobota et al., 2018).Soil profiles S01-S06 were collected at increasing distances from the front of Waldemarbreen (soil profile S06 was located on the flat outwash fan at the foreland of the moraine ridge of Aavatsmarkbreen).The profile depths varied from 12 to 17 cm.All profiles were covered by a layer of organic material.Recovered soil cores were divided into 1-2 cm sub-samples and dried at 105 • C to a constant weight, passed through a 2 mm sieve, and then prepared for radionuclide analysis.Additionally, 12 cryoconite samples (coded 2-13) were collected at Waldemarbreen in August-September 2014.The samples were taken from cryoconite holes at different altitudes, depths, and surface areas (see Table S1 in the Supplement).Material was collected with disposable plastic Pasteur pipettes and transferred to 15 cm 3 plastic test tubes.

Radionuclide analyses
Gamma-emitting radionuclides ( 137 Cs and 210 Pb) were measured by a planar HPGe (high-purity germanium) detector.The activities of 137 Cs were determined using its emission peak at 662 keV and for 210 Pb 46.6 keV was used.Correc-tions were made for the effect of self-absorption of lowenergy γ -rays (46.6 keV) within the sample, although these corrections were insignificant because of the very low sample weight (2 to 12 g).The activities of 238 Pu, 239+240 Pu, 241 Am, 234,238 U, and 230,232 Th were determined for 0.5-1.15g cryoconite samples and 5 g dried soil samples.Organic matter contents were determined by loss on ignition (LOI) at 600 • C for 6 h.The samples were then dissolved using concentrated acids (HF, HNO 3 , HCl).Details of the sequential radiochemical procedure used for determination of 238 Pu, 239+240 Pu, 241 Am, 234,238 U, and 230,232 Th are described in previous publications (Łokas et al., 2010, 2018).The 240 Pu/ 239 Pu atom ratio was determined using Thermo Fisher Scientific Neptune multicollector inductively coupled plasma mass spectrometry (MC-ICP-MS; analytical method based on Taylor et al., 2001;Łokas et al., 2018).The combined measurement uncertainties of activity concentrations, which include counting statistics and counting efficiency uncertainties, are reported as 1σ .The Environmental Radioactivity Laboratory of the Henryk Niewodniczanski Institute of Nuclear Physics holds an ISO17025 accreditation for gamma-ray spectrometric measurements and Pu analyses by alpha-particle spectrometry.

Heavy metal analyses
Subsamples (0.5 g) were taken for measurement of selected heavy metals considered to be important anthropogenic contaminants (Pb, Cu, Zn, Cd).These metals and others (Fe, Al, Ni, Cr, Mn) were measured by flame atomic absorption spectrometry (AAS Shimadzu 6800) using deuterium background correction.Cd and Pb stable isotope ratios were measured by ICP-MS (PerkinElmer SCIEX ELAN 9000).Measurement procedures are described in detail in Zaborska et al. (2017).Metal enrichment factors (EFs) were calculated according to the formula detailed in Zaborska et al. (2017).Mean background metal concentrations were taken from literature for the Hornsund region.The fraction of anthropogenic Pb (Pb ant ) was calculated according to the formula detailed in Zaborska et al. (2017).The excess 206 Pb/ 207 Pb ( 206 Pb/ 207 Pb excess ) ratio calculation was used to source Pb pollution (Farmer et al., 1996;Bindler et al., 2001).The quality control and assurance for this method are presented in Zaborska et al. (2017).

Statistical analyses
Regression analysis was performed to find correlations between particular radionuclides and metal concentrations.The correlations characterized by r > 0.5 are assumed to show moderate correlation, with r > 0.7 indicating a strong correlation (e.g.Mukaka, 2012).The nonparametric Kruskal-Wallis test was used to observe the difference between radionuclide and metal concentrations in soil and cryoconite samples.All www.the-cryosphere.net/13/2075/2019/The Cryosphere, 13, 2075-2086, 2019 3 Results

Radionuclide contents and activity ratios in tundra soil profiles and cryoconite samples
Activity concentrations and inventories of anthropogenic radionuclides ( 137 Cs, 238 Pu, 239+240 Pu, 241 Am) for tundra soil profiles and for cryoconite are shown in Table S2 and Fig. 2 and in Table S4 and Fig. 2, respectively.The radionuclide inventory is defined here as the activity of a given radionuclide contained in the soil column at unit surface area.Most of the profiles have concordant depth distributions of all artificial radionuclide activity with distinct peaks occurring in the surface or initial subsurface layers.The highest activity concentrations of anthropogenic radionuclides (65 ± 7 Bq kg −1 for 137 Cs, 0.09±0.01Bq kg −1 for 238 Pu, 2.13±0.16Bq kg −1 for 239+240 Pu, and 0.90 ± 0.06 Bq kg −1 for 241 Am) are observed in profile S04 at 1 cm depth where significant amounts of organic matter were observed (21 % LOI).Profile S02 shows no detectable 137 Cs (with a minimum detectable activity (MDC) range of 0.1 to 3 Bq kg −1 ).Activities of other radionuclides in this profile were not determined.In most samples activity concentrations of 238 Pu are close to the detection limit; consequently activity ratios of 238 Pu/ 239+240 Pu were calculated for only three samples.The maximum activity concentrations for cryoconite are 2000 ± 300 Bq kg −1 for 137 Cs, 0.08±0.02 to 2.1±0.2Bq kg −1 for 238 Pu, 0.09±0.02 to 43 ± 3 Bq kg −1 for 239+240 Pu, and 0.25 ± 0.06 to 25 ± 2 Bq kg −1 for 241 Am.All of these values are significantly higher than those observed in soil samples collected near this glacier.
The activity ratios of 238 Pu/ 239+240 Pu, 239+240 Pu/ 137 Cs, and 241 Am/ 239+240 Pu, as well as atom ratios of 240 Pu/ 239 Pu, were calculated to potentially distinguish the sources of these radionuclides in tundra soils and cryoconite samples.The mean 238 Pu/ 239+240 Pu activity ratio for tundra is 0.039 ± 0.010, lower than for cryoconite at 0.050 ± 0.009.The analysis of atom ratios of 240 Pu/ 239 Pu provides important information which potentially allows a more precise identification of the origin of Pu isotopes.Average 240 Pu/ 239 Pu atom ratios in the tundra profiles were 0.178 ± 0.003 and 0.147 ± 0.009 in cryoconite.The mean 239+240 Pu/ 137 Cs activity ratio in tundra (0.042 ± 0.010) is higher than in the cryoconite samples (0.020 ± 0.004).Average 241 Am/ 239+240 Pu activity ratios for tundra and cryoconite were comparable (0.46 ± 0.07 and 0.47 ± 0.07).
Concentrations of the natural radioisotopes ( 210 Pb, 234,238 U, 230,232 Th) measured in the tundra soil profiles and cryoconite are shown in Table S3 and Table S5, respectively.Soil profile activity concentrations range from 15±1 to 173± 4 Bq kg −1 for total 210 Pb, 13 ± 1 to 28 ± 2 Bq kg −1 for 234 U, 13 ± 1 to 30 ± 2 Bq kg −1 for 238 U, 9 ± 1 to 22 ± 2 Bq kg −1 for 230 Th, and 8 ± 1 to 42 ± 3 Bq kg −1 for 232 Th.The highest activity concentration for 210 Pb was found in the surface layer of profile S04 (concentration in cryoconite was 12 300 ± 600 Bq kg −1 ).Activity concentrations for U and Th isotopes in cryoconite samples are slightly higher than those in soils. 234U/ 238 U activity ratios in soils and cryoconite are similar.

Heavy metal analyses in tundra soil profiles and in cryoconite samples
The concentrations of metals and stable Pb isotope ratios ( 206 Pb/ 207 Pb, 208 Pb/ 206 Pb) measured in soil and cryoconite samples are shown in Tables S6 and S8, respectively.Pb and Cd concentrations in soils and cryoconite are shown in Fig. 3. Pb concentrations varied from 9.6 to 25.9 mg kg −1 .Cd concentrations range from 0.05 to 0.6 mg kg −1 , Zn concentrations range from 45.1 to 76.9 mg kg −1 , and Cu concentrations range from 11.8 to 26.6 mg kg −1 .The concentrations in cryoconite were significantly higher than in soils.Pb concentrations ranged from 19.9 to 97.7 mg kg −1 .Cd concentrations varied from 0.2 to 0.6 mg kg −1 , Zn concentrations ranged from 59.6 to 97.5 mg kg −1 , and Cu concentrations  S7) and cryoconite (r = 0.78) (Table S9).Pb also correlated with naturally derived metals Ni and Co (r = 0.48-0.59) in soils and in cryoconite Ni (r = 0.67).Cu correlated with Zn (r = 0.77) only in soils.Interestingly, Cd did not correlate with any of these metals in soils or cryoconite.Correlations between contaminating heavy metals and soil composition markers Fe and Al were assessed to test the influence of the natural variability of soil composition on metal concentration.Generally, there was no significant correlation in soils or cryoconite.
The metal enrichment factors (normalized to Al) are generally comparable for soil and cryoconite samples and ranged from 1.1 to 1.8 for Pb, from 0.9 to 1.6 for Zn, from 0.8 to 2.1 for Cu, from 0.4 to 7.7 for Cd (Table S6), from 1.5 to 6.1 for Pb, from 1.0 to 1.5 for Zn, from 1.3 to 2.3 for Cu, and from 0.5 to 4.5 for Cd (Table S8).For most metals there were significant differences in concentrations observed in the soil and cryoconite samples (Table S10).The only exception is Cd, which was similar (and elevated) in both environments.There was no difference in 208 Pb/ 206 Pb in soil and cryoconite samples.This is not surprising as this measure is less specific than 206 Pb/ 207 Pb (the 208 Pb/ 206 Pb coal burning signature has a ratio similar to a natural rock).

Radionuclides and heavy metal contamination
All soil profiles show activity concentrations and inventories typical of raised marine terraces covered by dry tundra in the western part of Spitsbergen (Łokas et al., 2014, 2017a, b;AMAP, 2015;Dowdall et al., 2003;Gwynn et al., 2004;Hardy et al., 1973).Retention of airborne radionuclides in tundra soils is facilitated by the presence of a well-developed organic horizon.Cryoconite materials are highly enriched with atmospheric radionuclides compared to tundra soils.This enrichment can be explained to some extent by the focusing of fine-grained material, washed down into very small depressions where the cryoconite samples were collected, and also by the properties of cryoconite material.Cryoconite granules are mixtures of mineral particles, organic matter, and microorganisms, which can retain and concentrate airborne radionuclides due to the metal-binding properties of extracellular substances that are excreted by microorganisms (Gadd, 2004;Francis, 2007).Differences in activity concentrations of anthropogenic radionuclides ( 137 Cs, Pu isotopes, and 241 Am) found among the cryoconite samples reflect the different locations of sampling sites, morphological features, and organic material content.A significant correlation was observed between the amount of organic matter present and the concentration of all airborne radionuclides (r = 0.79-0.83).We also note a significant correlation between 210 Pb, 137 Cs, Pu isotopes, and 241 Am with altitude (r = 0.64-0.91)and the location of the cryoconite holes (r = 0.58-0.63).This is not the case for 210 Pb (r = 0.38).Cryoconite hole depth correlates with 210 Pb and Pu isotopes (r = 0.60-0.68)(Table S8).Cryoconite samples collected at the highest locations (4, 12, 13) with the largest hole surface area (Table S1) were characterized by the highest observed activity concentrations.These three points were located close to each other so the material deposited could have originated from the same source and migrated with supraglacial streams.In contrast, sampling points 5, 6, and 7, which were located close to the moraine at lower altitude, were characterized by the lowest activity concentrations of anthropogenic radionuclides.These materials could have been washed out or mixed with mineral particles from the moraine.The highest 210 Pb activity concentration was observed in sample 4 (> 12 000 Bq kg −1 ).More uniform and comparable values (∼ 5000 Bq kg −1 ) were found in cryoconite samples 3, 8, 11, 12, and 13.This 210 Pb distribution is very similar to the distribution of anthropogenic radionuclides.Cryoconite granules can be locally redistributed on the surface of ice, depending on glacier morphology and altitude.It was noted that activity concentrations increase in relation to the altitude.Material deposited in the highest locations can migrate with supraglacial streams and be transported along with the glacier.Samples 4 and 13 were collected from typical cryoconite holes with the largest surface area (476 and 1200 cm 2 ) (Table S1).Regardless of the lifespan of individual cryoconite holes, their collapse does not imply removal of cryoconite from glacier surface as the dispersed cryoconite granules initiate formation of new holes (Takeuchi et al., 2001a, b).
The constant delivery of 210 Pb from the atmosphere suggests that the concentration of this radionuclide in cryoconite material should be proportional to the exposure time, while high concentrations of the artificial radionuclides indicate a significant contribution of material that was exposed to stratospheric or tropospheric fallout.The results show that 210 Pb is being deposited at a more constant rate; lower values for samples 3, 8, 11, 12, and 13 may suggest that the majority of the material was already removed with meltwater or that material in these samples is younger than in sample 4 (thereby having had less contact time with the atmosphere).
The effect of glacial morphology on effective trapping and storing of radionuclides was also observed in the Georgia glacier (Łokas et al., 2018).The differences in radionuclide concentrations between sampling points and the lack of clear differences in the elevation gradient from terminus towards icefall may reflect the heterogeneous topography of the glacier tongue.In the Waldemarbreen we observed differences in the elevation gradient from the top of the ice.
Activity concentrations of 234,238 U and 230,232 Th show little variability within the profiles and do not differ from values reported for soils globally (UNSCEAR, 1993), in the tundra soils in northern Spitsbergen (Dowdall et al., 2003) and in southwestern Spitsbergen (Łokas et al., 2017b).Soil activity concentrations of lithogenic radionuclides are related to their source minerals and are not influenced by the action of cryoconite.The presence of U and Th is the main contributor to natural terrestrial radioactivity.U isotopes ( 234 U and 238 U) in terrestrial samples (rocks, soils, and sediments) are usually present in radioactive equilibrium.The main sources of U in the natural environment are atmospheric precipitation of terrigenous material, soil resuspension, and rock weathering.
Our data for U and Th in cryoconite samples are comparable to other reported data.Our previous paper (Łokas et al., 2018) showed similar results in cryoconite granules from the Georgia glacier, while on a Swiss glacier (Baccolo et al., 2017), 238 U and 230 Th are higher.The 234 U/ 238 U activity ratio measured in cryoconite granules suggests a state of radioactive equilibrium.The 230 Th/ 232 Th activity ratio in cryoconite varies from 0.4 to 0.7.In only a few samples the activity ratio of 230 Th/ 232 Th is about 1 (samples 2, 3, 4, 12, 13).Isotopes derived from U ( 230 Th) and Th series ( 232 Th) are not correlated in cryoconite samples (r = 0.36) but are better correlated (r = 0.80) in tundra soil profiles.This may suggest that Th from the U series is removed with meltwater from cryoconite samples located at the top of the glacier.We also observed differences in U and Th activity concentrations in tundra soils and cryoconite samples, with tundra soils having lower values than in cryoconite.
The metal concentrations found in soils were typical or slightly higher than natural background levels.This is confirmed by generally low metal enrichment factors (EFs), with the exception of Cd.Metal distributions measured in this study are very similar to those of parent rocks (natural crust) in Spitsbergen.Rocks in Hornsund contain very low concentrations of heavy metals (Pb ∼ 11 mg kg −1 , Cd ∼ 0.04 mg kg −1 , Zn ∼ 58 mg kg −1 , Cu ∼ 20 mg kg −1 ; Samecka- Cymerman et al., 2011).Enrichment in Cd and Pb was found in surface and subsurface soil profile layers (Fig. 3).Similar enrichments (Pb up to 25.3 mg kg −1 , Cd up to 0.4 mg kg −1 ) in soil profiles of Kaffiøyra were found in the 1980s by Plichta and Kuczyńska (1991).Generally low metal concentrations have been mapped by Ottesen et al. (2010) in alluvial soils in the Kaffiøyra region.Alluvial sediments are believed to represent natural metal levels, although recently deposited layers can be polluted by anthropogenic activity; Cd concentrations range from 0.01 to 0.13 mg kg −1 , Cu from 13.4 to 43.4 mg kg −1 , Pb from 8.9 to 24.3 mg kg −1 , and Zn from 58.9 to 79.2 mg kg −1 .Metal concentrations were studied in peat cores collected from nearby Kongsfjorden.Cu ranged from 7 to 25 mg kg −1 , Pb from 10 to 35 mg kg −1 , and Zn from 20 to 75 mg kg −1 (Headley, 1996).The Kongsfjorden region might have been impacted by former mining activity, however (Hao et al., 2013).More recently, the metal distribution in soils of other Spitsbergen regions (Kongsfjorden and Adventfjorden) was studied by Halbach et al. (2017).Like our study, they found elevated concentrations of some metals, e.g.Cd (0.13-1.00 mg kg −1 ), in surface soils (covered by moss), but low concentrations of metals (Cd = 0.04-0.5,Cu = 3.3-28.4mg kg −1 , Pb = 2.9-22.7 mg kg −1 , and Zn = 25-106 mg kg −1 ) in the sandy fraction of soils (sampled at 20 cm depth).Slightly higher levels of metals (Cu: 17-92 mg kg −1 and Cd: 0.05-1.20 mg kg −1 ) were measured in biologically rich soils in the Hornsund region with the highest values in areas of rich vegetation (Wojtuń et al., 2013).Ziółek et al. (2017) report much higher metal concentrations in peat soils in Bellsund.In their study Pb reached 100 mg kg −1 , Zn reached 140 mg kg −1 , Cu 57 reached mg kg −1 , and Cd reached 8 mg kg −1 .In these cases the contamination of peat soils was caused by depositions of organic matter from a bird colony.In conclusion, an anthropogenic influence on soils in the Kaffiøyra region was limited.
While soils can receive metals from the atmosphere and from the surface run-off, cryoconites are fed only by atmospheric particles.Thus, interestingly, metal concentrations in cryoconite were higher than in soil samples (Table S8) and significantly exceeded natural values.Particularly elevated were the concentrations of Pb and Cd.The enrichment of other metals was lower (Zn and Cu).The maximum concentration of Pb in cryoconite samples from Hansbreen in Hornsund was lower (82.7 mg kg −1 ), but the concentrations of Zn (108 mg kg −1 ) and Cd (1.53 mg kg −1 ) were higher, however (Łokas et al., 2016).Singh et al. (2013) found a similar range of metal concentrations in cryoconite samples from a glacier located in the Kongsfjorden area (Pb up to 85.1 mg kg −1 , Cu up to 44 mg kg −1 , Cd up to 0.14 mg kg −1 , Zn up to 150 mg kg −1 ).Chmiel et al. (2009) did not find elevated metal levels in glacier waters in the Scott Glacier in Bellsund, but they did report enrichment in some metals (Pb = 190 mg kg −1 ) in cryogenic sediments.

The source of radionuclides and heavy metal contamination
Activity ratios of 238 Pu/ 239+240 Pu, 239+240 Pu/ 137 Cs 241 Am/ 239+240 Pu, and 240 Pu/ 239 Pu atom are commonly used to identify and distinguish between global (stratospheric) and regional (tropospheric) sources of these radionuclides (Oughton et al., 2004;Hirose and Povinec, 2015;Łokas et al., 2013, 2014, 2016, 2017b). 238Pu/ 239+240 Pu activity ratios in cryoconite samples are consistently higher than in tundra soils, suggesting a contribution from sources other than the global Pu fallout in this region.The proportion of global fallout versus local sources was modelled assuming that the 238 Pu/ 239+240 Pu activity ratios for global fallout and Chernobyl plutonium are about 0.025 and 0.24, respectively.The following formula described by Mietelski and W ąs (1995) was used: where f G represents the fraction of plutonium isotopes from global fallout sources (%), A R represents the 238 Pu/ 239+240 Pu ratio of Chernobyl sources (0.24),A M represents the measured 238 Pu/ 239+240 Pu activity ratio in cryoconite samples, and A G represents the 238 Pu/ 239+240 Pu activity ratio of global fallout (0.025).
Calculations were performed for each of the studied cryoconite samples; the fraction of global fallout varied from 82 ± 1 % to 97 ± 1 % (Fig. 4).Samples 4, 8, and 12, with the highest activity ratios of 238 Pu/ 239+240 Pu (0.062), also contained the highest contribution of Pu from other sources (18 % and 17 %, respectively).The highest activity ratios of 238 Pu/ 239+240 Pu (0.075 to 0.09) were observed in the initial soils from the vicinity of the Werenskiold glacier terminus in the southwestern part of Spitsbergen.These values fall within the range found for cryoconite from the Hans glacier (with a mean of 0.064) and other Arctic glaciers (Łokas et al., 2016; unpublished data) (Fig. 5a).The 240 Pu/ 239 Pu atom ratios measured in tundra soils (average of 0.179) (Fig. 5b) are comparable to global fallout ratios (0.183 ± 0.009) reported by Efurd et al. (2005) at 70 • N latitude and 0.180 reported by others (Muramatsu et al., 1999;Warneke et al., 2002).Similar results were reported for tundra soils from Spitsbergen (80 samples) with an average value of 0.179 (Łokas et al., 2017a).Thus, the tundra sites exposed to the atmosphere over the whole period of anthropogenic radionuclide release acquired their radionuclide contents primarily from direct atmospheric deposition.In tundra soils, 240 Pu/ 239 Pu atom ratios and the 238 Pu/ 239+240 Pu activity ratios are similar to the global fallout signature (0.180 and 0.025, respectively); however, in the cryoconite samples measured in this study (and other unpublished data from Arctic cryoconite) they are lower (0.147 and 0.050, respectively) (see Fig. 5a, b).Due to the slow accumulation of tundra soils and their continuous exposure to the atmosphere, the vertical profiles of the fallout radionuclides represent their time-averaged deposition on the soil surface.In contrast, different cryoconite samples from the same glacier could be exposed to fallout for various periods resulting in preferential uptake of radionuclides from single deposition events and deviations from the global fallout signatures.In a previous article (Łokas et al., 2017a) we assumed that the main sources of Pu in the proglacial soils (and now also in cryoconite samples) were nuclear tests and nonexploded weapons-grade material.The lack of correlation between excess 238 Pu and 240 Pu/ 239 Pu atom ratios (Fig. 5c) suggests that the excessive 238 Pu is derived from sources other than 239 Pu, such as pure 238 Pu from the SNAP9A satellite re-entry (Corcho Alvarado et al., 2014).A similar 238 Pu enrichment was observed in other cryoconite samples from Arctic glaciers (unpublished data) (Fig. 5c) and in air filters in Finland (Salminen-Paatero et al., 2012) and NW Poland (Kierepko et al., 2016).
The average 239+240 Pu/ 137 Cs ratios in tundra soils (0.042) are similar to the global fallout value (0.039) expected for the year 2013 (UNSCEAR, 1982); in cryoconite samples (0.020) it is lower.The observed average of about 0.021 ± 0.03 in the proglacial zone of the glacier from west Spitsbergen is remarkable as there were no other values as low as at this site.Similar or slightly higher values (0.021-0.035) were observed on Spitsbergen by Gwynn et al. (2004) and Łokas et al. (2014).The enrichment of soils in 137 Cs relative to plutonium might be due to the contributions from sources other than global fallout or these ratios may actually reflect the de-position ratios at this latitude.The mean 241 Am/ 239+240 Pu ratios for tundra soils (0.47) and cryoconite samples (0.46) fall within the ranges reported for soils from western Spitsbergen (Łokas et al., 2014, 2013) and other Arctic sites (Holm et al., 1983;Smith et al., 1995;Gwynn et al., 2005).
The relatively high 206 Pb/ 207 Pb ratio (1.190-1.217)measured in Kaffiøyra soils is close to the natural ratio for parent rocks (> 1.215) and/or represents a mixture of natural and anthropogenic Pb sources (Zaborska et al., 2017).The metallic trace elements, which are transported in the atmosphere attached to airborne particulate matter, tend to be retained and concentrated in the cryoconite material.This finding is confirmed by relatively lower 206 Pb/ 207 Pb ratios in cryoconite samples.These ratios suggested that Pb in cryoconite is a mixture of natural and anthropogenic sources and that anthropogenic Pb origin prevails.Using the endmember method and based on stable Pb isotope ratios, the contribution of anthropogenically derived Pb was estimated in cryoconite.It was found that anthropogenic Pb constituted 10 %-100 % of the Pb accumulated in cryoconite samples (Fig. 4).In the case of the Horsund site (Hansbreen) the anthropogenic lead fraction in cryoconite was lower and ranged from 29 % to 95 % (Łokas et al., 2016).It was also calculated that the excess Pb in cryoconite material (Pb measured -Pb natural ), was characterized by 206 Pb/ 207 Pb of 1.156-1.189(Fig. 4).Side stations located closer to rocky glacier sides were characterized by higher (more natural) calculated anthropogenic Pb isotopic ratios (1.17-1.19)compared to stations located in the central transect of the glacier (1.16-1.17).Since glaciers are fed by airborne particles 206 Pb/ 207 Pb ratios in cryoconite can be compared to aerosol 206 Pb/ 207 Pb.The study of aerosols in the Kongsfjorden region recognized the dominance of two sources of Pb.In the spring, isotope signatures of 206 Pb/ 207 Pb of ∼ 1.158 suggest an eastern Eurasian source, while in the summer a signature of 206 Pb/ 207 Pb of ∼ 1.167 suggests a mixture of North Amer- ican and Eurasian influence (Bazzano et al., 2015).Thus, the cryoconite 206 Pb/ 207 Pb ratios compare well with aerosol 206 Pb/ 207 Pb, particularly in the centre part of the glacier.This implies that it is mostly anthropogenic and not natural Pb that accumulates in the cryoconite samples and it can be said that an anthropogenic environmental impact is clearly visible in the Spitsbergen cryosphere.
Both radionuclides and heavy metals are deposited from the atmosphere on land, ocean, and glacier surfaces by wet and dry precipitation.Since there are no new sources of radioactive pollution in the atmosphere, the patterns of atmospheric circulation do not influence radionuclide distribution on different glaciers on Svalbard.The amount of precipitation received by the glacier can definitely influence the glacier pollution since coastal areas (e.g.Kaffiøyra) receive twice the amount received in inland areas (Førland et al., 2011).Post-depositional processes, mainly those connected to cryoconite hole formation and eventual material transport within the glacier, appear to have a more significant influence on the levels of radionuclide pollution within the glacier.A different scenario concerns the distribution of heavy metals.Since they are still emitted to the atmosphere from different sources, the direction of atmospheric circulation is an important factor shaping their distribution within the glacier.The primary atmospheric metal discharge originates from Europe and Russia (Isaksen et al., 2016).In coastal regions local winds can enhance the introduction of contaminants due to transport of contaminated sea salt aerosols (Lüdke et al., 2005Kozak et al., 2015).

Conclusions
The extent of variations in activity concentrations of airborne radionuclides and some trace metals (Pb and Cd) in cryoconite samples from Waldemarbreen exceeds ranges observed in soils surrounding this glacier.The mechanisms responsible for accumulation of atmospheric contaminates are not clear but we observed the effect of glacial morphology on trapping and storing of airborne radionuclides.The average atom ratios of 240 Pu/ 239 Pu for tundra soils are comparable to the characteristic ratio for global fallout, but for cryoconite samples (0.147) they are lower than for surrounding soils and also for other initial soils from the west coast of Spitsbergen.We observed an excess of 238 Pu in cryoconite samples similar to that observed in other initial soils in Spitsbergen.This enrichment of 238 Pu is likely derived from sources other than 239 Pu.This could be derived from pure 238 Pu from a satellite re-entry after injection from the stratosphere into the troposphere.Monitoring of cryoconite samples permits the identification of Pu release events in the form of a fallout associated with each separate event.Therefore, cryoconite is an excellent candidate for atmospheric contamination monitoring.The lithogenic radionuclide contents vary in much narrower ranges and their variability reflects changing proportions of weathering products derived from diverse parent rocks.
Data availability.All raw research data are presented in the Supplement.
Supplement.Full data about environmental radioactivity are available in the Supplement.The supplement related to this article is available online at: https://doi.org/10.5194/tc-13-2075-2019supplement.

Figure 1 .
Figure 1.Study area location.Profiles S01-S06 were collected at increasing distances from the front of Waldemarbreen and cryoconite samples (2-13) from different altitudes of this glacier.Pleiades 1 satellite scene composite acquired in 31 July 2017 ©CNES.

Figure 3 .
Figure 3. Concentrations of selected heavy metals (Pb, Cd) in cryoconite samples and tundra soil profiles.

Figure 4 .
Figure 4. (a) The contribution of global fallout Pu calculated using plutonium isotope ratios and (b) contribution of anthropogenic origin Pb calculated using Pb isotopes and the end-member method.

Figure 5 .
Figure 5.The relationship between 240 Pu/ 239 Pu atom ratios and 238 Pu/ 239+240 Pu activity ratios in cryoconite samples from this study and other Arctic cryoconite (unpublished data) (a); the relationship between 240 Pu/ 239 Pu atom ratio and sample number in tundra soils, cryoconite samples (this study), and other Arctic cryoconite (unpublished data) (b); the correlation between excess 238 Pu and 240 Pu/ 239 Pu atom ratios in cryoconite samples (this study) and other Arctic cryoconite (unpublished data) (c).