The relationships between soil chemical properties and vegetation succession in the aspect of changes of distance from the glacier forehead and time elapsed after glacier retreat in the Irenebreen foreland (NW Svalbard)

The aim of this study was to determine relationships between soil chemical parameters (i.e. content of total organic carbon, total nitrogen, total sulphur, soil pH) and vegetation development in relation to distance from the current glacier forehead and time elapsed after glacier retreat in the Irenebreen foreland. Three transects were designated along the foreland. Species and vegetation cover were investigated in 1 m2 plots, placed every 50 m along each transect; corresponding soil samples were collected and chemically analysed. The total organic carbon and total nitrogen contents in soil change according to power and exponential functions, respectively, whereas soil pH decreases linearly with increased time elapsed after glacier retreat. The wide variation in total sulphur contents prevent the determination of clear relationships. Bryophytes and vascular plants dominate in the younger part of the foreland, whereas epigeic lichens prevail in the older part. Vegetation cover seems to be the main factor effects on soil properties; however, chemical soil properties and distance from the glacier forehead affect species distribution and vegetation cover. Considered so far as the first pioneers, the epigeic lichens need more time than bryophytes and vascular plants to colonise the foreland.


Introduction
The Arctic is one of the regions most threatened by climate change (Arndt et al. 2010;Hansen et al. 2012;IPCC 2007IPCC , 2014. Many studies suggest that climate changes may have an important impact on element cycles in the forelands of Arctic glaciers, and thus on the development of vegetation cover during the primary succession (e.g. Oechel and Vourlitis 1994;Wietrzyk et al. 2016;Yoshitake et al. 2011). Documentation of glacier range changes is very broad and includes both polar and mountain regions (Burga et al. 2010;Dai et al. 2002;Dümig et al. 2011;Elberling et al. 2004;Grieve 2000;Kabala and Zapart 2012;Nakatsubo et al. 2005). Studies clearly show a loss of glaciers on all continents, progressing especially in the second half of the twentieth century (Nuth et al. 2013;Zemp et al. 2008). In Svalbard, the progressive recession of glaciers has been observed since the end of the Little Ice Age (Błaszczyk et al. 2009). Glaciers have been reduced in terms of both thickness and area (e.g. James et al. 2012;Kohler et al. 2007;Nuth et al. 2013). Retreating glaciers uncover new land surfaces within postglacial valleys and these freshly deglaciated surfaces are characterised by a lack of biological life, including vegetation cover, whose development is connected with the availability of nutrients and water in soil (Burga et al. 2010). On the other hand, the supply of dead organic matter increases along with vegetation development, which results in increasing carbon accumulation in the soil. Nowadays, the phenomenon of glacier recession affords an opportunity to investigate changes in carbon and nitrogen accumulation in soils developed in forelands from the barren ground in the immediate vicinity of the actual glacier front to the oldest part of the foreland, as well as the interrelated nature of vegetation succession and chemical properties of soils.
In previous studies in the glacier forelands of the High Arctic (Elberling et al. 2004;Kabala and Zapart 2012;Nakatsubo et al. 2005) and in alpine regions (Burga et al. 2010;Dai et al. 2002;Dümig et al. 2011;Grieve 2000), great attention was paid to carbon accumulation and its cycle. Very few studies refer to the relationship between carbon storage and vegetation (e.g. Nakatsubo et al. 2005;Vilmundardóttir et al. 2015;Yoshitake et al. 2011). Apart from carbon, several studies consider nitrogen accumulation in the soil of Arctic glacier forelands (Crocker and Dickson 1957;Hodkinson et al. 2003;Jacobson and Birks 1980;Kabala and Zapart 2012;Yoshitake et al. 2007) and alpine glaciers (Noll and Wellinger 2008). Up to now, total sulphur content has not been studied in the soils of Arctic glacier forelands, despite the fact that sulphur is one of the most important elements taking part in the metabolic processes of plants such as photosynthesis and fixation of nitrogen (Brady and Weil 2004). Moreover, its content in the soil may be responsible for primary vegetation succession in the glacier foreland. Sulphate and sulphur contents have been examined only in the foreland of the Damma Glacier in Switzerland and were positively correlated with the organic carbon distribution (Noll and Wellinger 2008;Prietzel et al. 2013).
In terms of vegetation succession, the most detailed analyses in Svalbard's glacier forelands were associated with studies of vascular plants (e.g. Pirożnikow and Górniak 1992;Prach and Rachlewicz 2012;Yoshitake et al. 2011). Other components of glacier foreland vegetation, in particular cryptogams, have usually been overlooked and have not been identified at the species level. The high number of cryptogamic species colonising freshly deglaciated forelands, considerably exceeds the number of vascular plants species (Wietrzyk et al. 2016(Wietrzyk et al. , 2017. This indicates the need for comprehensive research including both vascular plants and cryptogamic species in analyses of vegetation development and its relationships with soil properties. In the early twentieth century, the theory of a linear and directional character of primary succession of vegetation had begun to be questioned because besides the time after glacier retreat, local environmental variables had started to be considered as factors affecting vegetation development as well (Walker and Del Moral 2003). Previous research conducted in Svalbard shows that in addition to the activity of glacier rivers (Moreau et al. 2008), distance from the glacier forehead (Wietrzyk et al. 2016) is also an important factor influencing the process of vegetation succession.
The main aim of the present study was to determine how selected chemical properties of the soil, i.e. content of total organic carbon (TOC), total nitrogen (TN), total sulphur (TS), and soil pH, change with (1) distance from the current glacier forehead; (2) time elapsed after glacier retreat; (3) vegetation cover development in the foreland of Irenebreen. The second aim, considering the same problem from another side, was to estimate which factor, i.e. distance from the current glacier forehead, time elapsed after glacier retreat or soil properties, is of stronger impact for succession of vegetation. In addition, relationships between particular species occurrence, distance from the glacier forehead, time elapsed after glacier retreat and soil chemical properties were analysed. We set the following hypotheses: (1) both increased distance from the current glacier forehead and time elapsed after deglaciation positively influence TOC, TN, and TS contents in soil and negatively influence soil pH; (2) with regard to the varying rate of the glacier's retreat over time, the relationship between time elapsed after glacier retreat and soil properties is close to being linear; in contrast, the relationship of distance to soil properties is non-linear in character; (3) spontaneous development of vegetation contributes to increased contents of TOC, TN and TS and reduced pH in the soil developed in the foreland; (4) development of vegetation cover is dependent mainly on the distance from the glacier forehead.

Study area
The fieldwork was conducted in the Irenebreen foreland (Fig. 1). The Irenebreen, one of the forelands bordering the Kaffiøyra Plain (Oscar II Land, NW Spitsbergen), is a polythermal valley glacier with an area of 4.1 km 2 , approximately 4 km in length and between 1 km (in its frontal zone) and 1.5 km (in the eastern section) in width (Sobota 2011(Sobota , 2014Sobota and Lankauf 2010;Sobota et al. 2016). The glacier's average retreat rate is 13 m per year (Sobota and Lankauf 2010). However, the rate of glacier recession is, apparently, not constant, and this process has clearly accelerated since the beginning of the twenty-first century (IPCC 2007(IPCC , 2014Sobota et al. 2016). The Irenebreen foreland is approximately 1700 m wide and 1200 m long (Fig. 2). The southern part of the foreland is relatively flat and covered by proglacial river sediments. In contrast, the northern part is characterised by hilly relief. The Irenebreen foreland consists of a series of Vendian rocks known as the Hecla Hoek formation, which also comprises the peak parts of the mountains surrounding the Kaffiøyra glaciers. The Hecla Hoek formation consists mainly of quartzite, phyllite, diamictite, marble, and green slate. Furthermore, the mountains east of the Kaffiøyra Plain consist of thick tilloids, phyllites with frequent intercalations of marble, green phyllites, and occasional metabasic rocks (metadiabase and metagabbro) (Sobota and Lankauf 2010). The immediate vicinity of the Irenebreen is composed of the Comfortlessbreen diamictite formation of Cryogenian and/or Ediacaran age.
The climate of the Kaffiøyra Plain is affected by the West Spitsbergen Current. Summer mean air temperature reaches 4.6°C; the mean annual air temperature is −5.8°C. January is the coldest month, while July is the warmest (mean monthly air temperature: −14.2°C and 4.9°C, respectively). Cloud cover, which averages more than 70% in summer and autumn, is lower in winter and spring, affecting the number of hours of sunshine, which in April and May reaches 250 h. In the Kaffiøyra region during the summer months, mean air humidity is 89% and mean precipitation amounts to 51 mm (Przybylak and Araźny 2006). The vegetation which surrounds the Irenebreen and its foreland belongs to the alliance Luzulion arcticae Gjaerev. 1950 and is composed mainly of the following species: Cerastium arcticum, Salix polaris, Saxifraga cernua, S. cespitosa, S. oppositifolia, Brachythecium glaciale, Polytrichastrum alpinum, Cetrariella delisei, Ochrolechia frigida, Cetraria islandica (Węgrzyn and Wietrzyk 2015). In the Irenebreen foreland, vegetation is mainly represented by initial stages of plant communities in which lichens and bryophytes are the dominant organisms.

Data sampling
Fieldwork was carried out in summer 2012. The study design included investigation along transects from the current glacier forehead to the terminal moraine. Preliminary fieldwork was carried out across the entire glacier foreland in order to determine the part with the most homogeneous site conditions. Three linear transects, each 1150 m long, were designated along the foreland. Sampling sites were designated roughly every 50 m along the transects and included various sites in terms of parent material, geomorphology, and wetness. At each sampling site, overall vegetation cover, including lichen, bryophytes, and vascular plants, was collectively estimated on a percentage scale within a 1 m 2 plot, and the percentage cover for each species was determined as well. In the case of taxonomically problematic specimens, samples were collected for precise identification in the laboratory. Simultaneously, soil samples from the surface to a depth of 5 cm were collected from each sampling site. Consequently, 72 sampling sites in total were taken into account, i.e. 24 per transect (Fig. 1). Moreover, three representative soil profiles were excavated at three different distances (i.e. 0 m, 550 m, and 1150 m) from the current glacier forehead in order to identify general changes in soil parameters with depth and soil samples from each horizon (AC, Cg1, Cg2) being taken. In addition, the depth of dead-ice occurrence was determined for each profile. The sampling depth differed from profile to profile, depending on the nature of the soil profile.
Glacier ranges in the past and estimated ages of foreland materials in the proglacial zone were assessed based on Sobota and Lankauf (2010). The time period between two different glacier ranges (isochrones) was divided into  equal parts appropriate for the number of years in each particular time period. Then, a specific age was assigned to given sampling site. This interpolation enabled estimation of the time elapsed since deglaciation for each sampling site.

Soil analyses
Soil samples were carefully cleaned of plant debris to avoid contamination by organic matter during analysis. The samples were air-dried, gently crushed, and sieved through a 2-mm sieve. All laboratory analyses were done on fine-earth material (fraction <2 mm). Concentrations of TOC, TN, and TS were determined for each soil sample. Previously, all calcareous soil samples had been treated with 10% HCl to remove carbonates. Three 1-g replicates were analysed by means of dry combustion technique using a LECO CNS 2000 automated analyser (LECO 1996), and the mean value was treated as one observation in further analysis. Soil pH was measured in distilled water in a 1:2.5 (w/v) ratio according to the potentiometric method (Tan 1995) using a pH meter (Elmetron CPI-505) with a glass electrode (Elmetron ERH-11).

Species identification
Standard light microscopy were used for taxonomic determination of specimens; the following publications were used: Atherton et al. (2010); Brodo et al. (2001); Osyczka and Węgrzyn (2008); Rønning (1996); Smith et al. (2009);and Thomson (1984;1997). Additionally, chemical analyses of lichen substances were conducted according to the technique of Orange et al. (2001). Herbal materials were deposited in the Herbarium of the Institute of Botany at Jagiellonian University in Kraków (KRA).

Statistical analyses
Two data matrices were created for statistical analyses. The first, which included data on soil chemical properties and vegetation cover in each analysed plot at a certain distance from the glacier forehead, will be referred to hereafter as the distance data matrix; the second, which included data on soil chemical properties and vegetation cover averaged for particular numbers of years elapsed after glacier retreat, will be referred to hereafter as the time data matrix. First, each sampling site was assigned to one of five distance intervals within each transect based on the geomorphological diversity of the studied moraine ( Fig. 1 and S1). The normality of the distribution of the variables TOC, TN, and TS contents and pH in the distinguished intervals was verified using the Kolmogorov-Smirnov test. Data which did not meet the assumptions of normality were Box-Cox transformed in order to find an optimal normalising transformation. Then, one-way analysis of variance (ANOVA) was performed, followed by Tukey's (HSD) test, to test differences in particular element contents and pH across the studied intervals. Prior to the analysis, Levene's test was performed to assess the equality of variances.
The relationships between element contents, soil pH, distance from the glacier forehead, time elapsed after glacier retreat, and vegetation cover were calculated separately, using Pearson correlation coefficients for distance matrix and time matrix separately. Next, scatter plots were created for an initial evaluation of the nature of the aforementioned relationships. Subsequently, various functions, i.e. linear, quadratic, power, exponential, logistic, Michaelis, and third-and fourth-degree polynomials, were fitted to the data using Levenberg-Marquardt optimisation. The Akaike Information Criterion (AIC) was used to select the best-fitting model. The analyses of relationships between soil parameters and time/distance were performed on the basis of values averaged for the same values of distance and time elapsed after glacier retreat defined for three transects; whereas the analyses of relationships between soil parameters and vegetation cover were performed on two separate data matrices separately.
Cluster analyses of all analysed soil parameters, distance from the glacier forehead, and time elapsed after glacier retreat were carried out using an unweighted pairgroup average (UPGMA) clustering algorithm and Euclidean distance separately for the distance and time matrices.
As regards epigeic vegetation data, canonical correspondence analysis (Legendre and Legendre 1998) was applied to relate the cover of a particular species in the studied plots to soil variables, distance from the glacier forehead and time elapsed after glacier retreat. Data from plots without epigeic vegetation were excluded from the analysis. A forward selection procedure of the explanatory variable from the set of all variables available was applied. A Monte Carlo permutation test (499 permutations) was used in order to assess the statistical significance of relationships between species, environmental factors, and axes (Ter Braak and Šmilauer, 2002). Furthermore, the effect of soil variables, distance from the glacier forehead and time elapsed after glacier retreat on the vegetation cover in the studied plots were evaluated by separate stepwise multiple linear regressions using forward variable selection with a threshold of p < 0.05 to entry. Strongly correlated variables (R > 0.95) were removed from the analysis. A detailed residual analysis was performed in order to obtain reliable regression coefficients and detect potential outliers (extreme cases). Finally, a seriation of all species (presence/absence values) recorded in the examined plots, combined for three transects, was performed using an unconstrained algorithm (Brower and Kile 1988).

Soil development in the glacier foreland
Along the studied transects, TOC in the surface soil layer increased from a very low content of 0.69% at the newly exposed sites, which were close to the current edge of the glacier, to 3.81% in the oldest foreland soils, which were approximately 100 years old; the content of TN increased from 0.04% to 0.15%. The content of TS was the most diversified along the transects, reaching the lowest values within the first 200 m (~0.02%) and the highest at 750 and 900 m (0.09%). Contrastingly, soil pH decreased along the transects, ranging from 8.1 in the vicinity of the glacier forehead to 6.3 at the furthest distance from the current glacier position. According to ANOVA, significant differences were found in soil parameters between particular distance intervals (TOC: F = 164.67; TN: F = 97.81; TS: F = 9.63; pH: F = 62.65; p < 0.05 for all cases). Moreover, the two intervals located furthest from the glacier forehead were characterised by the highest TOC and TN concentrations, in both cases significantly different from the others (Fig. 3a, b). Conversely, soil pH values were significantly lower in the last two intervals (Fig. 3d). Differences in TS concentration were not as distinct; however, the highest mean TS contents were determined in the last two intervals (Fig. 3c). Cluster analysis of soil parameters and distance from the glacier forehead showed that distance was most strongly associated with TOC content (Fig. S2a), whereas TOC and TN contents were highly correlated with time elapsed after glacier retreat (Fig. S2b).
We found that both TOC and TN contents positively depend on distance from the glacier forehead and time elapsed after glacier retreat (Table 1). Generally, both TOC and TN contents increased with time and distance from the glacier (Figs. 4, S2). Similarly, the content of TS was positively associated with both distance and time; however, the correlation was not as strong (Table 1). Contrastingly, soil pH decreased significantly with distance and time (Table 1). However, scatter plots provided a basis for believing that at least some of these relationships are non-linear in nature. The best-fitting functions based on AIC values for each relationship are provided in Table S1. Although different types of functions were defined for TOC-distance and TN-distance relationships, the shapes of the curves of these functions were generally similar (Fig. 4a, c). Fitting function graphs indicated that the increases in the contents of these elements with distance were initially slow, but gradually increased (Fig. 4a, c). A similar trend applied to TOC content and time (Fig. 4b). Contrastingly, TN content increased linearly with time elapsed after glacier retreat (Fig. 4d). Linear functions also proved to best fit relationships between TS content and pH with distance and time (Fig. 4e-h).
Generally, in the studied soil profiles, TOC and TN contents decreased with depth (Table S3). The highest concentrations of these elements were found in surface AC horizons. The greatest decline was recorded between the AC and Cg1 horizons, especially at the end of the studied transects (Profile 3). The content of TS was almost evenly distributed throughout all of the studied soil profiles (Table S3). The distinct increase in pH with depth was noticeable only at the site located furthest from the glacier forehead (Profile 3); however a slight increase in pH with depth in the soil profile excavated in the middle part of the transect (Profile 2) was also observed (Table S3). The dead-ice layer in the soil profile located at the beginning of the transects started at an approximately 20 cm depth; in the middle part of transects, dead-ice was recorded at an approximately 35 cm depth, and at the end of transects at an approximately 45 cm depth.

Vegetation development in the glacier foreland
In total, 52 species were recorded within the studied plots: 13 species of bryophytes, 30 species of lichens,  and 9 species of vascular plants (Fig. 7). In the analysis of the occurrence of particular species along the glacier foreland, groups of early and late colonisers may be distinguished: bryophytes, namely Bryum calophyllum and Encalypta alpina, were the first species to colonise the freshly deglaciated area, whereas lichens began to occur starting at 350 m from the glacier forehead and dominated towards the end of the foreland (Fig. 7). Vascular plants were present from a distance of 200 m from the glacier forehead through the rest of the foreland. The most abundant species of vascular plants were Saxifraga oppositifolia and Cerastium arcticum (Fig. 7). Relationships between species distribution and soil chemical properties, distance from the glacier forehead and time elapsed after glacier retreat were determined using CCA analysis (Fig. 9). The Monte Carlo permutation test showed that both the first axis and all canonical axes taken together were statistically significant (p = 0.02 and p = 0.04, respectively). Cladonia pyxidata, Lecanora epibryon, Lecidea ramulosa, Ochrolechia androgyna, Caloplaca nivalis and Polyblastia sendtneri were associated with high TOC, TN, and TS content. These species were also associated with lower soil pH and were present in the older part of Irenebreen foreland (Fig. 7). Bryophytes, such as Andreaea blyttii, Bryum pseudotriquetrum, B. turgidum, Ditrichum cylindricum, E. alpina, and Pohlia wahlenbergii exhibited opposite trend (Figs. 7 and 9). Nevertheless, the forward selection procedure of the variables showed that species distribution is associated only with distance from the glacier forehead and TOC content (F = 3.64; p = 0.001 and F = 2.94; p = 0.002, respectively). The obtained results indicated that distance from the glacier forehead and TOC content explained 7.75% of the overall variability of species occurrence. However, concerning overall vegetation cover, forward stepwise regression analysis with distance from the glacier forehead, time elapsed after glacier retreat and soil variables as predictor variables (F = 6.51; p < 0.05) showed that vegetation cover was mainly related to four of the analysed factors, i.e. soil pH (p = 0.027), distance from the glacier forehead (p = 0.030), content of TOC (p = 0.043) and TN (p = 0.024). This model explained 38.4% of the variability of vegetation cover.
Although there were some fluctuations in vegetation cover along the transects (Fig. 9), this parameter was evidently correlated with distance from the glacier forehead and time elapsed after glacier retreat (Table 1). Vegetation cover and soil parameters were also interrelated. Linear functions proved to be the best fitted for relationships between soil parameters and vegetation cover (Table S2). With increasing development of vegetation cover, contents of TOC, TN, and TS gradually increased, whereas pH decreased; however, the correlations between these parameters were stronger for the distance matrix than for the time matrix (Table 1, Fig. 5).

Discussion
In this paper, a modified chronosequence method is proposed. This novel approach provides nearcontinuous information on soil development, vegetation succession, and glacier retreat spanning one century. Despite the fact that the distance from the glacier forehead and time elapsed after glacier retreat are correlated, we examined their influence on the development of soil and succession of vegetation separately. Such an approach was applied because the rate of glacier retreat may be not constant over time. Because of that, the model of soil development and vegetation succession in relation to distance and time could be different. This enabled us to identify the nature of these relationships, which is particularly important in the context of the increasing rate of deglaciation process in recent years.

Soil development in the glacier foreland
Soil development is dependent upon various soilforming factors (parent material, topography, climate, biota, and time) and their fluctuations and interferences (Brady and Weil 2004). The obtained results confirmed our first hypothesis assuming that both increased distance from the current glacier forehead and time elapsed after deglaciation positively influence TOC, TN, and TS contents in soil and also they stand in line with our third hypothesis that changes in soil properties are mainly influenced by spontaneous development of vegetation.
The nature of changes in TOC and TN contents is similar, with slow increases at the outset followed by more rapid increases beginning at a distance of 800 m from the glacier. It was observed that soil development in the glacier forelands significantly accelerates under continuous vegetation cover (Alexander and Burt 1996;He and Tang 2008;. This may explain the nature of the obtained results, since vegetation cover did not exceed 30% up to a distance of 700 m from the glacier forehead, while at greater distances it reached significantly higher values (Figs. 7 and 9). As vegetation succession proceeded, an increase in TOC and TN contents in soil was recorded (Fig. 5). Apart from TOC and TN contents, TS content also follows a similar pattern. Sulphur is an essential nutrient for microorganisms, plants, and animals. In pristine regions with geological parent material categorised as S-poor due to little atmospheric S input and young soils, the bioavailability of S may limit ecosystem production and species succession (Schmalenberger and Noll 2010). We found that the TS content in the surface soil layer generally increased with time elapsed after glacier retreat and distance from the current glacier forehead; however, wide fluctuations were noticeable (Fig. 4). Moreover, in terms of distance intervals, the highest mean content of TS was recorded at 750-950 m from the current glacier forehead, compared to a lower mean content for the greatest distance interval, i.e. 1000-1150 m (Fig. 3). Similar observations were made by Prietzel et al. (2013) in foreland surface soil layers at the Damma Glacier in Switzerland, where the concentration of TS was at its lowest level in the fresh foreland material and its highest in mature soils, but generally was greater for 70-to 79-year-old than for 111-to 121year-old soils. In the foreland surface soil layer of the Damma Glacier, the content of TS was positively correlated with TOC content (Prietzel et al. 2013). In the present study we also found significant positive correlations between TS and TOC contents. However, older soils are not always characterised by higher TS contents than younger soils. Similar relationships were also observed in the soils of the forelands of Franz Josef Glacier and Damma Glacier (Allison et al. 2007;Prietzel et al. 2013).
The variable rate of glacier retreat over time is certainly one of the important factors influencing the observed pattern of dependence between distance and soil properties. In terms of distance, this rate was considerably higher over the past several years than previously (Fig. 6). The results indicate that the contents of TOC and TN in soils developed in the foreland of the Irenebreen are, in general, a function of time and distance from the glacier forehead. Nevertheless, the accumulation of these elements does not fit a simple linear model in accordance with distance from the glacier forehead. The patterns of changes with distance are very characteristic and the increase in the contents of these elements is expressed by non-linear functions. These Relationships between vegetation cover and TOC, TN, TS contents in surface soil layer and soil pH. Combined function graphs were prepared based on analyses performed on two separate data matrices for averaged values of soil parameters for the same values of the distance and time elapsed after glacier retreat defined for three transects. Fitted function formulas are presented in the graphs soil development measures are distance-related and fit power and exponential models for TOC and TN, respectively. Consequently, the accumulation of TOC and TN at distances such as the first 500 m from the current glacier forehead was lower than for the next 500 m; this was also connected with a shorter period of soil development. Nevertheless, the models created in this study stand in contradiction to relationships observed in the Werenskioldbreen foreland (south-western part of Spitsbergen), where TOC and TN accumulation fits a logarithmic model and, in the fourth/fifth decade after deglaciation, it has reached a quasi-steady state over a time span of 80 years (Kabala and Zapart 2012). These differences might be due to the different methodology used by these authors, namely, the classical chronosequence method, which also resulted in a smaller number of samples being collected in the foreland of the Werenskioldbreen glacier. In addition, Kabala and Zapart (2012) included in their study soils located only on typical moraine till, whereas soils formed from outwash sands and gravels, glaciofluvial and limnic deposits were excluded from their study. We collected soil samples from various sites in the foreland (not only at sites located on typical moraine till) without exclusion of any sites and this is most likely the second reason for a different pattern of TOC and TN content increase along the transects in our study in comparison with the pattern obtained by Kabala and Zapart (2012). Furthermore, the geomorphological character of each foreland, connected with thawing ice cores and thus surface stability, may influence the contents of TOC and TN along the foreland. The absence or deeper occurrence of thawing ice cores results in a more stable foreland surface, enabling faster vegetation colonisation and greater accumulation of elements such as TOC and TN. On the other hand, the general increasing trend in TN content found in the present study seems similar to that observed in the Midtre Lovénbreen foreland (Hodkinson et al. 2003), although the authors did not determine the best-fitting function for nitrogen content. The increase in the content of this element was relatively slow at first, but became more pronounced 60 years after glacier retreat. Different relationships were observed in the foreland of the Midtre Lovénbreen by Dong et al. (2016), where linear functions provided the best fit and TOC and TN contents increased steadily over a time span greater than 80 years after deglaciation. Consequently, our second hypothesis assuming that the relationship between time elapsed after glacier retreat and soil properties is close to being linear was confirmed only in the case of TN content, because the relationship between TOC content and time elapsed after glacier retreat was non-linear. Similarly, the above-mentioned hypothesis was not fully confirmed in the case of TS content, because the character of TS changes was linear with time elapsed after glacier retreat and also with the distance from the glacier forehead. The patterns and time trends of the surface soil layer TS status of young foreland-derived soils in cool, humid periglacial environments during initial stages of pedogenesis and Fig. 6 The relatioship between the time elapsed after glacier retreat and the distance from the glacier forehead for particular transects vegetation succession were caused by the combined effects of mineral weathering, accumulation and progressive microbial turnover of soil organic matter, and type of vegetation succession (Prietzel et al. 2013). Since the parent material of soils in the foreland of Irenebreen is S-poor (Bjørnerud 2010), we supposed that organic matter and vegetation cover were the major factors affecting the increase in TS content along with soil development in the studied foreland. The wide variations in TS content with time and distance could be somehow explained by parent material heterogeneity, influence of atmospheric aerosols containing sulphates (Udisti et al. 2016) or volcanic eruptions (e.g. Icelandic volcanos) during which large volume of ashes containing sulphur are produced, transported into the atmosphere and spread widely (Kekonen et al. 2005). Nevertheless, an increasing trend in the TS content is still observable due to the large contribution of the accumulation of TS in the surface soil layer as a consequence of the establishment of vegetation.
The expected negative correlation of soil pH with distance and time factors as well as with spontaneous development of vegetation was also confirmed by our results, which is in line with our first and third hypothesis. The pH of fresh foreland glacial till is quite high, fluctuating around 8. This is mainly due to the types of rocks representing the Hecla Hoek Formation and the mineralogical composition of glacial till, which includes a high proportion of calcium and magnesium carbonate (Chlebowski 1989). Since vegetation cover develops in forelands, organic acids from decaying plant debris cause a progressive decrease in soil pH (Fig. 5). In addition, the influence of summer precipitation and thawing of snow cover may also be responsible for leaching of carbonates from upper soil horizons leading to a decrease in soil pH. The increase in pH with depth in the studied soil profiles located in the middle and especially at the furthest sites from the current position of the glacier forehead indicates the leaching of carbonates (Table S3). However, it is not possible to determine which factor (i.e. development of vegetation cover or leaching) plays a more important role in decreasing soil pH in the studied foreland. We also found that the decrease in soil pH is a linear function of time in the foreland soils of Irenebreen (Fig. 4). Similar trends in pH were also found for the soils of the Midtre Lovénbreen foreland (Dong et al. 2016). On the other hand, Kabala and Zapart (2012) found that pH decreased with time and became finally stabilised approximately 80 years after glacier retreat. He and Tang (2008) also observed that pH in surface soil horizons tended to change asymptotically with time and thus progressive decreases in soil pH slowed with soil development and finally stabilised. The same pH changes were found in the recessional forelands of the Herbert and Mendenhall Glaciers in south-eastern Alaska where, after an initial rapid rate of decline in pH, changes become less rapid and finally stabilise (Crocker and Dickson 1957).
Although a linear function of pH decline proved to fit our data best, the increase in soil pH recorded at the last sampling sites may indicate that pH tends to stabilise in the Irenebreen foreland as well. Nevertheless, our second hypothesis assuming that the relationship between time elapsed after glacier retreat and soil pH is linear was confirmed; however, the results stand in contrast to the hypothesis that the relationship between distance and soil pH is non-linear in character.

Vegetation development in the glacier foreland
The increase in the number of species and their cover with distance from the glacier forehead indicates that vegetation gradually changes along the studied transect ( Figs. 7 and 9). The colonisation order is as follows: first appear bryophytes, then vascular plants and lichens. However, in the case of lichens, the first recorded species were epilithic taxa which overgrow rocks uncovered by the retreating ice, while epigeic species occur last (Fig. 7). According to phytosociological research on the climax Minuartia biflora community (Węgrzyn and Wietrzyk 2015), which surrounds the Irenebreen foreland, a gradual transition of species composition, from the pioneer stage to near-climax vegetation communities, can be observed in the glacier foreland. According to a study carried out in the Midtre Lovénbreen and Austre Lovénbreen (Moreau et al. 2008), the increase found in the frequency and cover of Stereocaulon alpinum is an indicator of climax vegetation. In the Irenebreen foreland, S. alpinum starts to occur approximately 500 m from the current glacier forehead. This may suggest that the presence of this lichen is not related strictly to climax communities. S. alpinum, such as another species of terricolous lichen, Rostania ceranisca (dominant taxon in the Irenebreen foreland), is a host for cyanobacterial symbionts with the ability to bind atmospheric nitrogen. This feature increases the colonisation potential of these lichens (Sancho et al. 2011).
In the majority of Earth's ecosystems, vegetation cover is considered to be a response to the gradient of soil parameters; however, in the Arctic region, in areas surrounded by glaciers, an important factor influencing vegetation cover is also the microclimatic impact of a glacier's presence. According to multiple regression analysis, except for soil properties such as pH, TOC and TN, vegetation cover depends on distance from the glacier forehead. This result partly confirmed our fourth hypothesis assuming that the development of vegetation cover is dependent mainly on the distance from the glacier forehead (Fig. 9). With regard to statistical significance, distance was the second significant factor (after soil pH) affecting the vegetation cover. Nevertheless, this model explained only 38.4% of the variability of vegetation cover. This indicates that further research focused on factors influencing the development of vegetation cover in a glacier foreland such as nutrients available to plants, soil humidity, soil texture, and presence of biological soil crusts (BSCs) is needed. Concerning the occurrence of particular species, the analysis of CCA along with forward selection procedure shows that distance from the glacier forehead and TOC content in the soil affect species distribution the most (Fig. 9). Nevertheless, regarding TOC content, it is mainly related to lichens, such as C. pyxidata, L. epibryon, L. ramulosa, and O. androgyna (Fig. 9). The occurence of these species is caused by the presence of plant debris and BSCs, which provide them with the substrate to grow on (Smith et al. 2009). Due to this fact, it seems that the positive relationship between the mentioned species and TOC content is indirect and connected with readily available plant debris and BSCs, which directly result in high TOC content in the soil. The majority of epigeic lichens dominate in the older part of the studied foreland (Fig. 7). This indicates the highest level of species diversity in the areas where ice disappeared early, the species had enough time to colonise the substrate, and where the influence of dead-ice thawing and glacier-affected microclimate is lower. A great abundance of lichens in the parts of the foreland exposed earliest was also observed in the Midtre Lovénbreen (Moreau et al. 2008). In studies conducted in the foreland of the Werenskioldbreen, Tishkov (1986) emphasised the importance of socalled micro-depressions, which accelerate the process of primary succession due to optimal snow accumulation conditions. Apart from that, various geomorphological processes transform glacial till in forelands. Among them, the most important geomorphological processes are related to the melting of ice-cores and activity of proglacial rivers. Many studies (e.g. Borysiak et al. 2015;Laffly and Mercier 2002;Mercier and Laffly Fig. 7 The seriation diagram using unconstrained algorithm of species presence/absence in the examined plots combined for three transects 2005; Moreau et al. 2008;Östrem 1959;Prach and Rachlewicz 2012;) indicate that the above-mentioned processes are the main and most influential in the formation of lowlands and mounds in glacier forelands, and thus appear to exert an important impact on vegetation development. In contrast to epigeic lichens, bryophytes dominate in the younger part of the studied foreland (Fig. 7). These include mainly minerotrophic species, which inhabit moist areas (Atherton et al. 2010;Walker et al. 2001) and seem to tolerate the erosional activity of proglacial rivers which occurs periodically in the areas located near the current glacier forehead (Fig. 7). Ditrichum cylindricum, in particular, prefers a distinctly wet habitat that undergoes periodic flooding (Walker et al. 2001). This is confirmed by CCA analysis, which shows that the majority of bryophytes are associated with substrates with lesser content of TOC, TN, and TS, but higher soil pH (Fig. 9). Bryophyte species recorded in the older part of foreland were associated with local micro-depressions, which provide higher moisture content in the substrate in comparison with surrounding mounds (Fig. S1). In the case of vascular plants, those whose presence is Fig. 8 The relationship between vegetation cover and the distance from the glacier forehead averaged for three transects Fig. 9 Canonical correspondence analysis ordination diagram showing relathionship between species (species abbreviations are explained in Fig. 7), soil variables (TOC, TN, TS, pH), distance from the current glacier forehead (distance) and time elapsed after glacier retreat (time). The red arrows and underlined variables mark the statistically significant factors. The variance explained by the first canonical axis reaches 45% whereas by all canonical axes 79% connected to the microclimatic influence of glacier and ice-cores, occur first in the colonisation process and are present only in the younger part of the foreland, e.g.: S. cesposita, S. cernua and P. alpina (Figs. 7 and 9). Conversely, S. oppositifolia and C. arcticum, which appear as the first vascular plant species, dominate towards the end of the foreland and show a greater tolerance for habitat conditions that change along the foreland (Fig. 7). The important role of S. oppositifolia and its domination in the succession process was also noted in the forelands of other glaciers in the northern (Hodkinson et al. 2003;Moreau et al. 2005Moreau et al. , 2009Tishkov 1986), southern (Pirożnikow and Górniak 1992;Wojtuń et al. 2008), north-eastern (Yoshitake et al. 2010), and central (Prach and Rachlewicz 2012) parts of the Svalbard Archipelago.

Conclusions
This paper provides a novel approach to studies concerning soil development and plant succession in a glacier foreland. The applied methodology provides nearcontinuous information on the soil development and vegetation succession spanning one century. Our results show that the contents of TOC and TN in soils developed in the Irenebreen foreland, rather than fitting a simple linear model in terms of their variation with distance from the glacier forehead, fit power and exponential models, respectively. The gradual colonisation of the foreland by cryptogams and vascular plants also affects soil pH, which decreases with time according to a linear function; however the additional impact of soil leaching related to summer precipitation and the thawing of snow cover cannot be entirely excluded. Similarly, the content of TS in soils appears to be affected by the formation of vegetation cover and accumulation of organic matter; however, its wide variations with time and distance from the current glacier forehead may also be explained by the heterogeneity of parent material and/or the influence of atmospheric aerosols and/or volcanic deposits. Bryophytes and vascular plants dominate in the younger part of the studied foreland, whereas epigeic lichens prevail in the older part. This indicates that epigeic lichens need more time than bryophytes and vascular plants to colonise a glacier foreland. The study showed that a particular species distribution is associated with the distance from the glacier forehead and TOC content in soil. Nevertheless, concerning the overall vegetation cover, apart from the above-mentioned factors, pH and TN seem to also impact the development of vegetation cover.