Open Access
How to translate text using browser tools
1 November 2017 The Influence of Seabirds on the Concentration of Selected Heavy Metals in Organic Soil on the Bellsund Coast, Western Spitsbergen
Marta Ziółek, Piotr Bartmiński, Alfred Stach
Author Affiliations +

Seabirds are an important factor affecting the chemical properties of the arctic soil environment. The objective of this work was to study the differences in content and distribution of heavy metals in organic soils resulting from the differential influence of seabirds. Studies were conducted in two stations in the southwest part of Spitsbergen—the first without the influence of seabirds and the second directly affected by an avian colony. Basic properties of soils as well as total content of Ca, Na, Mg, K, Fe, Al, Cu, Zn, Ni, Mn, Pb, Cd, Co, and Cr were determined. Reference was made to the previously published contents of different forms of phosphorus (P) for these locations. The studies showed that Zn, Mn, Cu, and Cd contents were higher in the soils that had been in the vicinity of the seabird colony. High statistically positive correlations of Cd, Cu, and Zn were noted with particular P forms. In the case of Pb, Cr, Co, and Ni content, the seabird influence was not dominant; probably other factors were more relevant (the processes of weathering, denudation, leaching, and atmospheric pollutants transported in the form of dust or gases). Clear segregation was observed of individual positions conditioned by selected soil features and by heavy metals content.


The arctic soil environment is subjected to changes, wherein natural processes play a fundamental role. Cryogenic processes influenced the current state of the Spitsbergen soil cover to the greatest extent, as well as humidity and weathering (Forman and Miller, 1984; Mann et al., 1986; Blümel et al., 1993; Melke and Chodorowski, 2006; Witkowska-Walczak et al., 2014; Bockheim, 2015). In the Bellsund region (southwestern Spitsbergen), organic and mineral-organic soils are relatively common, but they only cover small areas. Peat soils form under very moist tundra habitats, in areas of water stagnation, such as in local depressions of coastal plains, on hillsides, and at the bottom of slopes. Peat soil profiles are generally shallow and very shallow (Eurola, 1971; Zelikson, 1971; Låg, 1980; Göttlich and Hornburg, 1982; Klimowicz et al., 1997). They are characterized by organic matter that is poorly decomposed, and the peat layer is deposited on rock debris or heavy clay. The chemical composition of organic deposits in Spitsbergen has been poorly examined; only a few publications have reported on it (Dobrovolsky, 1990; Jóźwik and Magierski, 1992; Klimowicz et al., 1997; Melke, 1999; Melke and Chodorowski, 2006; Chmiel et al., 2009).

Polar land ecosystems, including soil cover, greatly depend on nutrients originating from the sea. Seabirds are the key organisms in the Arctic acting as biovectors and influencing energy and nutrient transportation from marine to land ecosystems (Tatur, 1989; Godzik, 1991; Headley, 1996; Ligęza and Smal, 2003; Michelutti et al., 2009, 2010; Mallory et al., 2015). From marine areas of foraging, birds carry nutrients to breeding areas, where they release them through defecation, regurgitation, dropping food, or mortality (Blais et al., 2005; Ellis et al., 2006; Mulder et al., 2011; Bauer and Hoye, 2014). In polar regions, the amount of nutrients delivered with avian excrement is much higher than from alternate sources and on a local scale it enriches the ecosystem significantly (Bildstein et al., 1992; Bokhorst et al., 2007). Introducing nutrients to the land food chain influences the chemical properties of soil and directly influences the flora in the proximity of avian colonies (Breuning-Madsen et al., 2008; Zwolicki et al., 2013, 2015; Ziółek and Melke, 2014).

In the bio-transportation of organic materials, transportation of pollutants and their concentration in the proximity of breeding areas can also be observed (Blais et al., 2005; Brimble et al., 2009; Mallory et al., 2015). Seabirds are the main source of the terrestrial environment's pollution because their food consists of organisms that are high in the food chain. As a result of bioaccumulation and biomagnification, these organisms contain a high concentration of pollutants, including heavy metals originating from the marine environment (Buckman et al., 2005; Campbell et al., 2005; Mallory and Braune, 2012). Enrichment of the terrestrial environment with heavy metals in the proximity of seabird breeding areas is documented for the polar (Godzik, 1991; Headley, 1996; Bargagli et al., 1998; Sun et al., 2004; Zhu et al., 2005; Evenset et al., 2007; Mallory et al., 2015) and temperate zones, as well as the tropical zone (Otero Perez, 1998; Hawke et al., 1999; Garcia et al., 2002; Liu et al., 2006). The efficiency of pollutant introduction to the terrestrial environment by seabirds is also visible in terms of the heavy metal content in flora. Plants in the proximity of avian colonies show higher concentrations than the same species growing in similar habitats distant from avian colonies (Godzik, 1991; Blais et al., 2005; Evenset, 2007).

Increased soil pollution is most often observed in areas subjected to direct pressure from human activity (Makuch, 2014). Nevertheless, the surveyed region's significant distance from any industrial areas does not exclude the influence of anthropogenic factors on the arctic environment. It is mostly manifested by deposition of gas and dust pollutants transferred over great distances (Pacyna, 1995; Law and Stohl, 2007; Ruman et al., 2013; Kozak et al., 2015). Local sources, such as coal mines, also have a small influence on the pollution of Svalbard's soil environment (Headley, 1996; Gulińska et al., 2003; Lewińska-Preis et al., 2009).

The objective of this research was to study the role of birds in the transport of pollutants, visible in the content and distribution of heavy metals in Spitsbergen organic soils. The presented results show a new perspective on the diverse accumulation of contaminants in organic soils and are useful for determining changes in the environment of the Arctic.

Study Area

The study was conducted at two stations in the northwest part of Wedel Jarlsberg Land, southwestern Spitsbergen (Fig. 1).

Within the area of the first station (Skilvika profile), there were organic soils not directly subjected to the influence of seabirds. The station is located on Calypsostranda plain in the direct neighborhood of a cliff coast on Sklivika Bay. In geological terms, it is a part of a system of elevated marine terraces with average height of 40–65 m a.s.l. where tundra polygons formed under the influence of frost weathering. The station area is covered with mossy flora with sparsely occurring clumps of grass and the Salix polaris willow.

The second station (Dunderdalen transect) is located on the southern slope of Dundrabeisen (435 m a.s.l.), at a height of around 30–90 m a.s.l. (Fig. 2.). There is a shallow, sloping peatland of less than 1 ha located here, formed on an alluvial cone, with a tilt of around 30°. Above the peatland, there is a very steep and plantless rock wall with bird nests, mostly blacklegged kittiwake, Rissa tridactyla. The number of nests was estimated at several hundred, and the share of other species was negligible. The influence of the avian colony can be seen in the whole area of the station (lush vegetation), with an obvious intensification directly under the nesting zone (organic remains: dry plants, bird droppings, and feathers). The research station's area is covered with mossy flora

The average air temperature in the area of study is around 5 °C in the summer and between -8 and -16 °C in the winter. Temperatures over 15 °C are recorded sporadically. The average annual temperature is around -5 °C. There is very little precipitation; the annual sum is below 400 mm. The highest amount of precipitation is observed in the fall season.

Materials and Methods

The Dunderdalen transect was designated along the fall of the alluvial cone and the peatland that formed on it. The material collection points were placed in an approximately 25-m sequence; the first one was located directly under the avian nesting zone and the last one on the flattening at the bottom of the cone. The thickness of the organic material was diversified, and one to three samples were collected from each individual point (15 in total).


Location of study area and survey points. Topographic map was used with the courtesy of the Norwegian Polar Institute.



The location of sample collection points at the Dunderdalen transect.


The Sklivika profile was described in the wedge between tundra polygons filled with peat moss. The depth of the wedges reached over 50 cm at the permafrost boundary. The station has an ombryogenic character. Peat, along with mineral substrate, was collected to a depth of 67 cm and divided into samples of 1.5–2.5 cm thickness (27 in total).

A detailed morphological description and results of basic properties (organic matter [OM] and pH) for the Skilvika profile and the Dunderdalen transect can be found in the work by Ziółek and Melke (2014).

The material was dried and ground in an agate mortar for further studies. The following basic analyses were carried out: soil organic matter content by burning at 550 °C, carbonate content using the volumetric method, and pH in 0.01 M CaCl2. Samples incinerated at 550 °C were dissolved in a solution of concentrated HF and HClO4 acids, evaporated to dryness, and then treated with 6 M HCl. Assessment of the total content of Fe, Al, Cu, Zn, Ni, Mn, Pb, Cd, Co, and Cr was carried out using the AAS method (Sparks et al., 1996). Analytical data were supplemented with data concerning the content of individual P forms in the studied soils (Ziółek and Melke, 2014).

The average values for the respective Skilvika profile depths were calculated for comparison with each other, and those depths refer directly to the depths from which samples along the Dunderdalen transect were collected, creating a summary of data referring to each other in terms of depth. Subsequently, the data were subjected to statistical analysis. Their purpose was to answer the following questions: (1) Do the concentrations of the individual components analyzed differ among sites? If so, what is the difference? To solve this problem, both single-variable and multi-variable statistical tests were used. Their more detailed description is provided below. (2) The difference between the posts may be because of not only the concentrations of the individual components, but also their specific combination. So, in the analyzed data, what specific classes are there and what are their characteristics? The multi-variable grouping methodology (described below) was used to answer this question. (3) Is there a relationship between the concentrations of the heavy metals tested and the concentrations of the different forms of P that may be seen as an indicator of the impact of the bird colony on the properties of the soils tested? Pearson correlations between individual parameters were calculated, and significance was adjusted by Holm's method (Holm, 1979). The results were confirmed by Spearman's rank correlation analysis.

Because of the small size of both study groups separated on the basis of location (n < 30), nonparametric methods were used to assess the significance of their differences. The following methods were used: U Mann-Whitney test (evaluation of the lack of differences in average), Brown-Forsythe test (assessment of the significance of differences in data variability), Wald-Wolfowitz runs test (conformity assessment of data distribution in populations). Typically, testing the significance of the differences among samples is limited to the mean values. It happens, however, that although the samples do not differ in average values, they may vary in terms of the variability (scattering of values around the mean) and/or the shape of the distribution (for example, its symmetry/skewness). Wald-Wolfowitz runs test is a nonparametric test that examines whether two populations differ in central tendency, variances, skewness, or any other distribution pattern. The above statistical tests allowed us to verify the assumptions about differences among particular (individual) variables in the populations from which the samples were derived. However, they did not allow us to determine the degree of similarity of the population resulting from the combination of many features. For multivariate assessment of similarity/dissimilarity of positions complementary methods: ANOSIM (one-way ANOSIM-ANalysis Of Similarities) and Simper (Similarity Percentage) were used (Clarke, 1993). The calculations were performed using PAST software (Hammer et al., 2001, Hammer and Harper, 2006) using Euclidean distance to the standardized data.

The similarity of individual soil samples was also determined, regardless of the place of collection. Similarity was determined from three perspectives: for all analyzed components of the soil, the macronutrients and pH (OM, Ca, Mg, Na, K, Al, and H+) and heavy metals only (Fe, Mn, Cu, Cb, Cd, Pb, Zn, Ni, Cr). Grouping using the Ward method in the Euclidean space was performed, on the basis of standardized data. A distance for which inflection of the connecting (agglomeration) curve followed was taken as a criterion for grouping. Such a collapse of the agglomeration curve indicated the start of grouping of objects of small similarity (above this level objects of low similarity are combined). In each case, similarity of variables was analyzed in addition to the similarity among the analyzed samples, because high similarity among variables may indicate their similar origins in the soil.

Physical Description

The organic soils of the studied stations are diversified in terms of morphological structure, organic matter content, and acidity. In the peat soil of Skilvika, the organic part of the profile reaches a depth of 57.5 cm, below which there is a mineral clay substrate. The organic part of the profile consists of peat moss in varying degrees of decomposition. In the upper part, it is very weakly decomposed and with the increase of depth its decomposition rises to the level of well-decomposed in the lower part.

In Dunderdalen, all samples in the seven points of the transect consisted of organic material deposited on rock debris. The first point of the transect—the surface sample—consisted of dry floral remains, avian excrement, and feathers. Samples of subsequent points consisted of peat moss at different degrees of decomposition: from weakly to mid-decomposed in samples at depths of 0–10 cm and 0–5 cm, and mid- to well-decomposed in deeper samples. Generally, the degree of peat decomposition rises along the transect and with depth.


Box plot for particular elements in soil samples from surveyed locations. D—Dunderdalen, S—Sklivika. Indents of boxes represent the range of 95% confidence interval for the median. White and gray filling of boxes was used to distinguish the groups of parameters (organic matter [OM] and pH macronutrient/light metal [i.e., Ca, Mg, Na, K, and Al] and heavy metals [i.e., Fe, Mn, Cu, Co, Cd, Pb, Zn, Ni, and Cr]). The dots represent data outliers, defined as greater than 3 quartile or less than 1 quartile more than 1.5 IQR (interquartile interval). Units for particular components are listed in Table 1.



For the most of the parameters analyzed, the differences between the positions were very significant and applied to the averages (medians), dispersion of data, and the general character of the distribution. Charts (Fig. 3) also showed a significantly higher incidence of outlier values in Skilvika (38 cases to 11 in Dunderdalen). Most of these samples were observed by reference to the Al concentrations (6 cases), and Mg, Fe, and Zn (4 cases). Large differences in medians for OM, Ca, Al, Fe, Pb, and Cr are visible. Mg, pH, Na, Zn, and Cr were very differentiated in terms of variability.

In the studied organic soils of Spitsbergen, samples from Dunderdalen transect were characterized by a much greater amount of Ca, Na, Zn, Mn, and Cu. The reverse situation occurred in the case of Fe, Al, Cr, Pb, and Co, where a significantly higher level was recorded for the Skilvika profile. There were only slight differences in the content of Mg, K, Cd, and Ni between Dunderdalen and Skilvika; the determined amounts of these elements were in a similar range (Table 1, Fig. 4).

The heavy metal content was inconsistently correlated with particular forms of P (Table 4). Zn, Cd, and Cu showed positive correlations with individual forms of P. Fe, Co, Ni, Pb, and Cr were not correlated with P.

The vertical distribution of the most studied heavy metals within the Skilvika profile indicated some relation with the ash content. The lowest content was observed in the surface part of the profile in the very weakly decomposed peat moss, with the highest content of organic matter (over 80%). The exceptions were Pb, which saw an increase in content between 2 and 5 cm, and Ni, for which a slight accumulation was observed in the surface sample (0–2 cm; Fig. 4). From the depth of 7 cm, where the ash content rises to over 40%, there was an increase in the occurrence of all studied elements, with a distinct accumulation at the depth of 15–17.5 cm. However, the maximal content for all studied microelements (Cu, Co, Cd, Pb, Zn, and Ni) was found in the well-decomposed peat moss at a depth of 52.5–55 cm. Mn and Cr are the exceptions; the highest content of Mn was reported in the upper part of the profile at a depth of 7–8.5 cm, whereas Cr reached its maximum content in the mineral substrate (Fig. 4).


Basic descriptive statistics of individual components in soil samples of surveyed stations. OM = loss combustion of organic matter. Prior to the calculation of the statistics for soil reaction pH values were converted to a concentration of hydrogen ion (H+). The results are restated in the pH scale.


The positioning of microelements in the Dunderdalen transect was more diversified both along the transect and in the respective points. Their maximum content was characteristic for some elements (Mn, Zn, Cu, and Cd) at point 2, just below the surface sample under the avian colony in the weakly decomposed peat moss (Fig. 4). This point also had the lowest pH in the whole transect. Ni, Co, and Cr also had heightened concentrations here, but their maximum content occurred at the final seventh point of the transect (Co and Ni in the 20–25 cm sample; Cr in the 5–15 cm sample of that point). The content of Pb exhibited a different trend, and its maximum occurred in the surface sample of the transect's fifth point (Fig. 4).

Of the 16 parameters taken into consideration, in 10 cases a statistically significant difference was found in averages, in 5 in variances, and in 9 in the whole distribution of the variables (Table 2). Two data sets—concentrations of Ca and Cr—were the most differentiated, because they differed in each analyzed aspect, and the probability of error was very low (p ≤ 0.01). The content of OM and concentrations of Fe, Al, and Pb differed significantly among stations because of average values and data distribution. In particular, concentrations of Fe and Pb differentiated the two groups, because the probability of error was less than 0.001. Reaction (pH) and concentration of Na and Cu were different in variances and the distribution of values. Concentrations of Mn, Co, Zn, and Ni were different only in mean values, and in the case of K and Cd there was no reason to reject a hypothesis of the lack of differences. It should also be noted that the test results for Cu and Ni were close to the threshold of statistical significance (0.035 < p < 0.049), so they should be interpreted with caution.

The results of the tests described above allowed us to evaluate the significance of differences in the statistics of individual variables. However, they did not allow us to determine the degree of differentiation resulting from the combination of many features simultaneously. For this purpose, ANOSIM and SIMPER methods were used. These calculations confirmed the difference between the two locations. For the data of all components analyzed, average rank of the intra-group distance in ANOSIM test was 368.3 (rb), and intergroup was 501.6 (rw). R test statistic calculated on their basis had a value of 0.309, and the probability (p) that the differences between the groups are random was only 0.0002. If ANOSIM used only the concentrations of heavy metals, the value of rb was 495.8, rw = 373.4, R = 0.2842, and p = 0.0006. The differences between soil in the studied positions, connected with the content of heavy metals, were also highly statistically significant, although somewhat weaker than in the other analyzed parameters (differences in R values: 0.309 vs. 0.284).


The distribution of heavy metals in the profile Skilvika and along Dunderdalen transect.


Results of SIMPER indicated that all the parameters determined in the soil samples showed fairly equal conditions of dissimilarity of soils at both positions (Table 3). The ratio of difference between the Ca that most affects the dissimilarity, and Ni, the effect of which was the smallest, was only like 1.69:1. The heavy metals' share in total dissimilarity was 54.8% (21,296/38,891) and was roughly proportional to their number in relation to all analyzed variables (9 of 16). Among them, the most important were Zn, Cd, and Pb, and the least, Ni and Co.


Test statistics and significance levels for comparison of the components in soil samples taken from the Dunderdalen and Skilvika stations. Explanations: M.-W. U = U test (Mann-Whitney), Brn-Fors. = Brown-Forsythe test, W.-W. = Wald-Wolfowitz runs test. Calculations for soil reaction (pH) were done on the basis of the data converted to hydrogen ion concentration. Cases with p < 0.05 are marked by gray background.


A dendrite of sample similarity constructed with the use of all parameters divided soils into five groups (Fig. 5). The division was made at a bond distance of 9.5, for which the agglomeration curve showed a distinct shoulder. Distinguished groups had from 2 to 17 objects. In classes 3 and 1, soil samples from only one station (S or D) were connected; in the other classes there were combinations of various samples. Nonetheless, even in the latter groups, on the lowest level of the dendrite (single branches), samples from different stations varied—there was no situation in which the sample of “S” was the most similar to the sample of “D” and vice versa. The largest class 3 was the most homogeneous one; it was characterized by a relatively small variation of all the parameters and values lying in the central part of their total range. The minimum concentrations of Ca and Mn were outstanding among the other classes.

Class number 5 showed high variability among its samples, in the case of most of the analyzed parameters. Its diversity was the result of a minimum content of OM and the highest concentrations of up to eight analyzed components (Mg, K, Fe, Al, Co, Pb, Ni, and Cr).

Class 1 varied the most from the other classes, having highest concentrations of Ca, Na, Mn, Cu, Cd, and Zn, and the lowest pH. The other two classes, 2 and 4, contained the samples of the two stations and had a rather complex internal structure. Class 2 was characterized by the highest mean content of OM and lowest concentrations of most of the other components (Mg, Na, K, Fe, Al, Cu, Co, Cd, Pb, Zn, Ni, and Cr). The highest average pH and relatively high variability of the other soil parameters (with average level of concentration) were specific for class 4.

Using the same methods, analyzed variables were also grouped on the basis of similarity (Fig. 6); four groups were distinguished. In the smallest one, OM content and the concentration of Ca were combined. The dendrite branch including this group also included groups consisting of Mn, Na, and Mg, and Cd, Zn, Cu, and H+. The last group, consisting of Pb, Cr, Al, Ni, Co, Fe, and K, was very dissimilar to both previous groups.

Additionally, we created a dendrite of similarity of soil samples, using only data for heavy metals (Fig. 7). Five groups were also distinguished, which tended to almost concurrently demonstrate a similar homogeneity as in the analysis of all parameters (Fig. 5). Only four samples from Dunderdalen: 3 (0–10), 3 (10–20), 4 (15–25), and 5 (10–20) were included in the other branch than in the first case. A dendrite of similarity of variables (Fig. 8) showed clear differentiation of heavy metals, which were combined into two groups consisting essentially of the elements having a higher concentration in Dunderdalen (Zn, Cu, and Mn) or Skilvika (Pb, Co, Cr, and Fe).


Results of SIMPER analysis for all soil parameters analyzed and only for heavy metals.



Dendrite of similarities among soil samples, determined on the basis of all the components analyzed. Numbers (1–5) specify classes described in the text.



Similarity of all variables (components in soil samples).



Dendrite of similarities among soil samples, determined on the basis of only heavy metals concentrations (Fe, Mn, Cu, Co, Cd, Pb, Zn, Ni, and Cr).



Similarity in the group of variables of heavy metals.



Pearson correlation coefficients between heavy metals concentration and particular phosphorus forms. Marked values are significant at p < 0.05 (p values were adjusted using Holm's method). Explanations: PT = total P, OP = organic P, IP = inorganic P, Pi-L = labile inorganic P, Pi-FeAl = inorganic P associated with Fe and Al, Pi-CaMg = inorganic P associated with Ca and Mg, Po-HuAc = organic P associated with humic acids, Po-Res = residual organic P.



Our study has shown that seabird colonies, as a consequence of bio-transportation processes, play an important role in increasing the content of and changing the distribution of heavy metals in arctic soils, and that this is a widely described phenomenon (Godzik, 1991; Headley, 1996; Blais et al., 2005; Michelutti et al., 2010; Mallory et al., 2015). This is indicated by the highest scores of Zn, Mn, Cu, and Cd contents in the peat under the avian colony, which is the direct receiver of those inputs (the second point of the transect), and is in agreement with other studies (e.g., Brimble et al., 2009; Michelutti et al., 2009). Moreover, the surface part of the peatland, not subjected to the direct influence of birds, exhibited lower contents of all heavy metals in relation to the substrate under the avian colony (with the exception of Pb).

The process of bio-transportation does not, however, increase the contents of all heavy metals. One of the surprising results is the higher contents of Pb, Cr, and Co found in the soil profile not subjected to the direct influence of birds. Apart from the influence of seabirds, organic soils of Spitsbergen are exposed to the accumulation of heavy metals from other sources. The supply and loss of the particular analyzed elements in the sediment can be (a) the processes of weathering/type and composition of the in situ bedrock; (b) denudation (surface runoff erosion, mass movements, deflation), leaching, accumulation, and other agencies (Gulińska et al., 2003; Mallory et al., 2015); (c) accumulation of atmospheric pollutants transported in the form of dust or gases (Plichta and Kuczyńska, 1991; Steinnes and Friedland, 2006; Rumam et al., 2013; Kozak et al., 2015); (d) biogenic accumulation (in the case of Spitsbergen, mainly related to birds). Each type of delivery is related to specific relationships among the elements. The dominant system of similarities arising from the elemental composition of the minerals that build bedrock and the relative mobility of individual elements can be modified when there is a significant supply or loss from other sources listed above. In the case of a comparison between the positions of Dunderdalen and Skilvika, significant differences in geomorphology of both objects should be noted—in the first case the slope, in the second, uplifted marine terrace. However, the material in which the designations were performed is peat, characterized by a permanent bonding of heavy metals, notwithstanding some environmental differences. Sampling along the slope profile enables the balancing of local (point) effects of erosion and accumulation processes. In addition, the relatively small distance between the positions allows the conclusion that the differences cannot be more the result of modern atmospheric input of pollutants (dusts and gases).

Results indicate the evident influence of the avian colony on heightened soil pollution concentration in the Dunderdalen region, though the influence is not the same for all cases of heavy metals. In the context of seabirds' influence, attention is especially drawn to the low correlation of most metals with individual forms of P—generally only Cu and Zn, and Cd to a smaller extent, correlate with P. The dependence of these elements on seabirds' influence was already indicated by researchers (Godzik, 1991; Brimble et al., 2009; Mallory et al., 2015). Seabirds' activity increases the amount of P in soils (Ziółek and Melke, 2014), so there is a clear impact of seabirds on the concentration of elements mentioned above.

Significantly lower concentrations of Pb in soil samples from Dunderdalen were determined, which excludes the influence of avian colonies on the content of this element in soil, what is also indicated by not significant correlation (Table 4). A similar observation was presented by Mallory et al. (2015) and was explained as a high natural content of lead. In the discussed case, probably the age of Skilvika peat (2040 ± 80 yr B.P., at a depth of 55 cm, Ziółek, not published) was of great importance, in connection with long-distant transport of pollutants in the past, presented in literature (Weiss et al., 1999; Liang and Mao, 2015).

The dendrite of similarity based on all analyzed environmental variables (Figs. 5, 6) depicts the relatively distinctive segregation of the research stations subjected to the pressure of seabirds and those without it. On the basis of test statistics (Table 2), it should be stated that parameters associated with macro-elements had a high input into dissimilarity, and there is an evident occurrence of organic matter as well. To verify to what extent the dissimilarity of both localizations is associated with heavy metals, an additional analysis was conducted, with the exclusion of the remaining features (Figs. 7, 8). The segregation, although less evident, is still clearly visible.


Studies show that the concentrations of Zn, Mn, Cu, and Cd in the organic soils of Spitsbergen are higher within the range of influence of seabird colonies. However, seabirds are not always the greatest vectors of heavy metals—for the concentrations of Pb, Cr, Co, and Ni, the influence of seabirds is not dominant; other factors are more relevant, for example, atmospheric pollution (current and historical) and geological substrate.

High correlation of certain heavy metal contents (Cu, Cd, and Zn) with various forms of P is a strong evidence for birds' influence. In the case of Co, Ni, Cr, and Pb, no significant correlation was observed.

Significant differences between the locations subjected to the influence of birds and those not influenced by seabirds were determined. Individual points show a clear dissimilarity, conditioned more or less equally by all analyzed parameters. Segregation of positions also remains clear, taking into account only heavy metal concentrations as variables determining the dissimilarity of two groups of soils.


The research was carried out on the basis of material collected by Prof. Jerzy Melke.

References Cited


Bargagli, R., Sanchez-Hernandez, J. C., Martella, L., and Monaci, F., 1998: Mercury, cadmium and lead accumulation in Antarctic mosses growing along nutrient and moisture gradients. Polar Biology , 19: 316–322. Google Scholar


Bauer, S., and Hoye, B. J., 2014: Migratory animals couple biodiversity and ecosystem functioning worldwide. Science, 344(6179): 1242552-1–1242552-8. Google Scholar


Bildstein, K., Blood, E., and Frederick, P., 1992: The relative importance of biotic and abiotic vectors in nutrient transport. Estuaries ,15: 147–157. Google Scholar


Blais, J. M., Kimpe, L. E., McMahon, D., Keatley, B. E., Mallory, M. L., Douglas, M. S. V., and Smol, J. P., 2005: Arctic seabirds transport marine-derived contaminants. Science , 309(5733): 445. Google Scholar


Blümel, W. D., Eberle, J., and Weber, L., 1993: Verwitterung, Genese und Bodenverbreitung im Liefdefjord/Bockfjordgebiet (NW-Spitzbergen)—Untersuchungsmethoden und erste Ergebnisse. (Weathering, genesis and soil distribution in Liefdefjord/Bockfjord region [NW Spitsbergen]—Research methods and first results). Basler Beiträge zur Physiogeographie , 15: 65–70. Google Scholar


Bockheim, J. H., 2015: Cryopedology. Cham-Heidelberg-New York-Dordrecht-London: Springer, 177 pp. Google Scholar


Bokhorst, S., Huiskes, A., Convey, P., and Aerts, R., 2007: External nutrient inputs into terrestrial ecosystems of the Falkland Islands and the Maritime Antarctic region. Polar Biology , 30: 1315–1321. Google Scholar


Breuning-Madsen, H., Ehlers, C. B., and Borggaard, O. K., 2008: The impact of perennial cormorant colonies on soil phosphorus status. Geoderma , 148: 51–54. Google Scholar


Brimble, S. K., Foster, K. L., Mallory, M. L., Macdonald, R. W., Smol, J. P., and Blais, J. M., 2009: High arctic ponds receiving biotransported nutrients from a nearby seabird colony are also subject to potentially toxic loadings of arsenic, cadmium, and zinc. Environmental Toxicology and Chemistry , 28(11): 2426–2433. Google Scholar


Buckman, A. H., Norstrom, R. J., Hobson, K. A., Karnovsky, N. J., Duffe, J., and Fisk, A. T., 2005: Organochlorine contaminants in seven species of Arctic seabirds from northern Baffin Bay. Environmental Pollution , 128: 327–338. Google Scholar


Campbell, L. M., Norstrom, R. J., Hobson, K. A., Muir, D. C. G., Backus, S., and Fisk, A. T., 2005: Mercury and other trace elements in a pelagic Arctic marine food web (Northwater Polynya, Baffin Bay). Science of the Total Environment , 351-352: 247–263. Google Scholar


Chmiel, S., Reszka, M., and Rysiak, A., 2009: Heavy metals and radioactivity in environmental samples of the Scott Glacier Region on Spitsbergen in summer 2005. Quaestiones Geographicae , 28A: 23–29. Google Scholar


Clarke, K. R., 1993: Non-parametric multivariate analysis of changes in community structure. Australian Journal of Ecology , 18: 117–143.[cite this reference] Google Scholar


Dobrovolsky, V. V., 1990: Geohimija počv Špicbergena. (Geochemistry of Spitsbergen soils). Počvovedenie , 2: 5–20. Google Scholar


Ellis, J. C., Farina, J. M., and Witman, J. D., 2006: Nutrient transfer from sea to land: the case of gulls and cormorants in the Gulf of Maine. Journal of Animal Ecology , 75: 565–574. Google Scholar


Eurola, S., 1971: The middle Arctic mire vegetation in Spitsbergen. Acta Agralia Fennica , 123: 87–107. Google Scholar


Evenset, A., Christensen, G. N., Carroll, J., Zaborska, A., Berger, U., Herzke, D., and Gregor, D., 2007: Historical changes in persistent organic pollutants and metals recorded in sediment from Lake Ellasjøen, Bjørnøya, Norwegian Arctic. Environmental Pollution , 146: 196–205. Google Scholar


Forman, S. L., and Miller, G. H., 1984: Time-dependent soil morphologies and pedogenic processes on raised beaches, Bröggerhalvöya, Spitsbergen, Svalbard Archipelago. Arctic and Alpine Research , 16: 381–394. Google Scholar


Garcia, L. V., Maranon, T., Ojeda, F., Clemente, L., and Redondo, R., 2002: Seagull influence on soil properties, chenopod shrub distribution, and leaf nutrient status in semi-arid Mediterranean islands. Oikos , 98: 75–86. Google Scholar


Godzik, B., 1991: Heavy metals and macroelements in the tundra of southern Spitsbergen: the effect of little auk Alle alle (L.) colonies. Polar Research , 9: 121–131. Google Scholar


Göttlich, K., and Hornburg, P., 1982: Ein Zeuge wärmzeitlicher Moore im Adventdalen auf Spitzbergen (Svalbard-Archipel.). (The peatbog in Adventdalen on Spitsbergen (Svalbard-Archipel.) as a witness of warm period). Telma , 12: 253–260. Google Scholar


Gulińska, J., Rachlewicz, G., Szczuciński, W., Barałkiewicz, D., Kózka, M., Bulska, E., and Burzyk, M., 2003: Soil contamination in High Arctic areas of human impact, central Spitsbergen, Svalbard. Polish Journal of Environmental Studies , 12(6): 701–707. Google Scholar


Hammer, Ø., and Harper, D. A. T., 2006: Paleontological Data Analysis. Malden, Massachusetts: Blackwell Publishing Ltd, 351 pp. Google Scholar


Hammer, Ø., Harper, D. A. T., and Ryan, P. D., 2001: PAST: paleontological statistics software package for education and data analysis. Palaeontologia Electronica , 4(1): 9 pp. Google Scholar


Hawke, D. J., Holdaway, R. N., Causer, J. E., and Ogden, S., 1999: Soil indicators of pre-European seabird breeding in New Zealand at sites identified by predator deposits. Australian Journal of Soil Research , 37: 103–113. Google Scholar


Headley, A. D., 1996: Heavy metal concentrations in peat profiles from the High Arctic. Science of the Total Environment , 177: 105–111. Google Scholar


Holm, S., 1979: A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics , 6(2): 65–70. Google Scholar


Jóźwik, Z., and Magierski, J., 1992: Trace elements in plants and soils of coastal plains of south Bellsund (western Spitsbergen). In Repelewska-Pękalowa, J., and Pękala, K. (eds.), Polar Session. Arctic Natural Environment Problems. Lublin: Maria Curie-Skłodowska University Press, 161–169. Google Scholar


Klimowicz, Z., Melke, J., and Uziak, S., 1997: Peat soils in the Bellsund region, Spitsbergen. Polish Polar Research , 18: 25–39. Google Scholar


Kozak, K., Kozioł, K., Luks, B., Chmiel, S., Ruman, M., Marć, M., Namieśnik, J., and Polkowska, Ż., 2015: The role of atmospheric precipitation in introducing contaminants to the surface waters of the Fuglebekken catchment, Spitsbergen. Polar Research , 34(1): 24207. Google Scholar


Låg, J., 1980: Special peat formation in Svalbard. Acta Agriculturae Scandinavica , 30: 205–210. Google Scholar


Law, K. S., and Stohl, A., 2007: Arctic air pollution: origins and impacts. Science , 315(5818): 1537–1540. Google Scholar


Lewińska-Preis, L., Fabiańska, M. J., Ćmiel, S., and Kita, A., 2009: Geochemical distribution of trace elements in Kaffioyra and Longyearbyen coals, Spitsbergen, Norway. International Journal of Coal Geology , 80: 211–223. Google Scholar


Liang, J., and Mao, J., 2015: Source analysis of global anthropogenic lead emissions: their quantities and species. Environmental Science and Pollution Research , 22: 7129–7138. Google Scholar


Ligęza, S., and Smal, H., 2003: Accumulation of nutrients in soils affected by perennial colonies of piscivorous birds with reference to biogeochemical cycles of elements. Chemosphere , 52: 595–602. Google Scholar


Liu, X., Zhao, S., Sun, L., Yin, X., Xie, Z., Honghao, L., and Wang, Y., 2006: P and trace metal contents in biomaterials, soils, sediments and plants in colony of red-footed booby (Sula sula) in the Dongdao Island of South China Sea. Chemosphere , 65: 707–715. Google Scholar


Makuch, I., 2014: Content of chosen heavy metals in garden and arable soils. Polish Journal of Soil Science , 47: 89–99. Google Scholar


Mallory, M. L., and Braune, B. M., 2012: Tracking contaminants in seabirds of Arctic Canada: temporal and spatial insights. Marine Pollution Bulletin , 64: 1475–1484. Google Scholar


Mallory, M. L., Mahon, L., Tomlik, M. D., White, Ch., Milton, G. R., and Spooner, I., 2015: Colonial marine birds influence island soil chemistry through biotransport of trace elements. Water, Air, & Soil Pollution , 226: 31. Google Scholar


Mann, O. H., Sletten, R. S., and Ugolini, F. C., 1986: Soil development at Kongsfjorden, Spitsbergen. Polar Research , 4: 1–16. Google Scholar


Melke, J., 1999: Preliminary studies of fissure and slope peats in the Calypsostranda and Lyellstranda Region, West Spitsbergen. In Repelewska-Pękalowa, J. (ed.), Polish Polar Studies. 26th International Polar Symposium. Lublin: Maria Curie-Skłodowska University Press, 171–175. Google Scholar


Melke, J., and Chodorowski, J., 2006: Formation of Arctic soils in Chamberlindalen, Bellsund, Spitsbergen. Polish Polar Research , 27: 119–132. Google Scholar


Michelutti, N., Keatley, B. E., Brimble, S., Blais, J. M., Liu, H., Douglas, M. S. V., Mallory, M. L., Macdonald, R. W., and Smol, J. P., 2009: Seabird-driven shifts in Arctic pond ecosystems. Proceedings of the Royal Society, B: Biological Sciences , 276: 591–596. Google Scholar


Michelutti, N., Blais, J. M., Mallory, M. L., Brash, J., Thienpont, J., Kimpe, L. E., Douglas, M. S. V., and Smol, J. P., 2010: Trophic position influences the efficacy of seabirds as metal biovectors. Proceedings of the National Academy Sciences of the USA , 107(23): 10543–10548. Google Scholar


Mulder, C. P. H., Jones, H., Kameda, K., Palmborg, C., Schmidt, S., Ellis, J. C., Orrock, J. L., Wait, D. A., Wardle, D. A., Yang, L., Young, H., Croll, D. A., and Vidal, E., 2011: Impacts of seabirds on plant and soil properties. In Mulder, C. P. H., Anderson, W. B., Towns, D. R., and Bellingham, P. J. (eds.), Seabird Islands: Ecology, Invasion, and Restoration. New York: Oxford University Press, 135–176. Google Scholar


Otero Perez, X. L., 1998: Effects of nesting yellow-legged gulls (Larus cachinnans Pallas) on the heavy metal content of soils in the Cies Islands (Galicia, North-west Spain). Marine Pollution Bulletin , 36: 267–272. Google Scholar


Pacyna, J. M., 1995: The origin of Arctic air pollutants: lessons learned and future research. Science of the Total Environment , 160: 39–53. Google Scholar


Plichta, W., and Kuczyńska, I., 1991: Metal contents in soils of Kaffiöyra, Spitsbergen. Polish Polar Research , 12: 183–193. Google Scholar


Ruman, M., Kozak, K., Lehmann, S., Kozioł, K., and Polkowska, Ż., 2013: Pollutants present in different components of the Svalbard Archipelago environment. Ecological Chemistry and Engineering , 19: 571–584. Google Scholar


Sparks, D. L., Page, A. L., Helmke, P. A., Loeppert, R. H., Soltanpour, P. N., Tabatabai, M. A., Johnston, C. T., and Sumner, M. E. (eds.), 1996: Methods of Soil Analysis. Part 3. Chemical Methods. Madison, Wisconsin: Soil Science Society of America, Book Series 5: 1390 pp. Google Scholar


Steinnes, E., and Friedland, A. J., 2006: Metal contamination of natural surface soils from long-range atmospheric transport: existing and missing knowledge. Environmental Reviews , 14: 169–186. Google Scholar


Sun, L. G., Renbin, Z., Xuebin, Y., Xiaodong, L., Zhouqing, X., and Yuhong, W. A., 2004: Geochemical method for the reconstruction of the occupation history of a penguin colony in the maritime Antarctic. Polar Biology , 27: 670–678. Google Scholar


Tatur, A., 1989: Ornithogenic soils of the maritime Antarctic. Polish Polar Research , 10: 481–532. Google Scholar


Weiss, D., Shotyk, W., and Kempf, O., 1999: Archives of atmospheric lead pollution. Naturwissenschaften , 86: 262–275. Google Scholar


Witkowska-Walczak, B., Sławiński, C., Bartmiński, P., Melke, J., and Cymerman, J., 2014: Water conductivity of arctic zone soils (Spitsbergen). International Agrophysics , 28: 529–535. Google Scholar


Zelikson, E. M., 1971: Palinologičeskoe issledovanie golocenovogo torfjanika na Špicbergene. (Palynological investigation of a Holocene peat deposit from Spitsbergen.) In Neustadt, M. I. (ed.), Palinologija golocena. (Holocene palynology). Moskwa: Academy of Science SSSR. 199–212. Google Scholar


Zhu, R., Sun, L., Yin, X., Xie, Z., and Liu, X., 2005: Geochemical evidence for rapid enlargement of a gentoo penguin colony on Barton Peninsula in the maritime Antarctic. Antarctic Science , 17: 11–16. Google Scholar


Ziółek, M., and Melke, J., 2014: The impact of seabirds on the content of various forms of phosphorus in organic soils of the Bellsund coast, western Spitsbergen. Polar Research , 33: 19986. Google Scholar


Zwolicki, A., Zmudczyńska-Skarbek, K. M., Iliszko, L., and Stempniewicz, L., 2013: Guano deposition and nutrient enrichment in the vicinity of planktivorous and piscivorous seabird colonies in Spitsbergen. Polar Biology , 36: 363–372. Google Scholar


Zwolicki, A., Barcikowski, M., Barcikowski, A., Cymerski, M., Stempniewicz, L., and Convey, P., 2015: Seabird colony effects on soil properties and vegetation zonation patterns on King George Island, Maritime Antarctic. Polar Biology , 38: 1645–1655. Google Scholar
© 2017 Regents of the University of Colorado
Marta Ziółek, Piotr Bartmiński, and Alfred Stach "The Influence of Seabirds on the Concentration of Selected Heavy Metals in Organic Soil on the Bellsund Coast, Western Spitsbergen," Arctic, Antarctic, and Alpine Research 49(4), 507-520, (1 November 2017).
Received: 15 March 2016; Accepted: 1 July 2017; Published: 1 November 2017
Back to Top