- Research article
- Open Access
Zn isotope fractionation in a pristine larch forest on permafrost-dominated soils in Central Siberia
© Viers et al.; licensee Springer. 2015
- Received: 17 September 2014
- Accepted: 1 April 2015
- Published: 16 April 2015
Stable Zn isotopes fractionation was studied in main biogeochemical compartments of a pristine larch forest of Central Siberia developed over continuous permafrost basalt rocks. Two north- and south-oriented watershed slopes having distinctly different vegetation biomass and active layer depth were used as natural proxy for predicting possible future climate changes occurring in this region. In addition, peat bog zone exhibiting totally different vegetation, hydrology and soil temperature regime has been studied.
The isotopic composition of soil profile from Central Siberia is rather constant with a δ66Zn value around 0.2‰ close to the value of various basalts. Zn isotopic composition in mosses (Sphagnum fuscum and Pleurozium schreberi) exhibits differences between surface layers presenting values from 0.14 to 0.2‰ and bottom layers presenting significantly higher values (0.5 – 0.7‰) than the underlain mineral surface. The humification of both dead moss and larch needles leads to retain the fraction where Zn bound most strongly thus releasing the lighter isotopes in solution and preserving the heavy isotopes in the humification products, in general accord with previous experimental and modeling works [GCA 75:7632–7643, 2011].
The larch (Larix gmelinii) from North and South-facing slopes is enriched in heavy isotopes compared to soil reservoir while larch from Sphagnum peatbog is enriched in light isotopes. This difference may result from stronger complexation of Zn by organic ligands and humification products in the peat bog compared to mineral surfaces in North- and South-facing slope.
During the course of the growing period, Zn followed the behavior of macronutrients with a decrease of concentration from June to September. During this period, an enrichment of larch needles by heavier Zn isotopes is observed in the various habitats. We suggest that the increase of the depth of rooting zone, and the decrease of DOC and Zn concentration in soil solution from the root uptake zone with progressively thawing soil could provoke heavy isotopes to become more available for the larch roots at the end of the vegetative season compared to the beginning of the season, because the decrease of DOC will facilitate the uptake of heavy isotope as it will be less retained in strong organic complexes.
In the context of climate warming boreal forested regions with continuous permafrost attract special attention because they are likely to encounter the most important changes. The latest technical study from the IPCC (www.ipcc.ch and references therein) on climate change and water resources reports abnormal permafrost thawing [1,2], a longer growing season [3-5], an increase of subsurface water drainage  and foresees significant changes in the distribution of plants and their productivity for arctic and sub-arctic regions . All these changes are likely to modify the fluxes of chemical elements, including carbon and metal micronutrients between the main reservoirs (soil-rock system, plants, aquatic domains, and atmosphere) [7,8]. Siberian forests are especially important in this regard because they store approximately 140 gigatonnes of carbon in above-ground biomass and at least three times more as soil organic matter [9-11]. However, the limiting factors in the development of these forests remain poorly understood. In contrast to the large number of articles devoted to macronutrients limitation in coniferous forests over permafrost, studies of micronutrients, especially the divalent metals, are notably scarce. Among the different metals, zinc is especially interesting because i) it is essential for plant growth similar to manganese, molybdenum and nickel , and ii) after iron, Zn is the most abundant metal in living organisms . In addition to the availability of nutrients in soils [14-18] internal mechanisms such as nutrients resorption or leaf-span longevity sustaining plant growth by increasing the mean residence time of nutrients within the plants are also important .
Until now, measuring Zn concentrations in the whole plant biomass, sap solution and individual organs remained the basic technique used by scientists for studying Zn biogeochemistry . Over the past decade, with the appearance of a new generation of MC-ICP-MS, stable isotope methods have been widely employed to characterize Zn transport mechanisms and sources within continental or marine environments [20-41] comprising field studies, in vitro experiments on plants, experimental studies on Zn isotope fractionation during interaction with mineral and organic compounds and paleo proxies. The present work extends our knowledge of Zn fractionation in soil-plant system to a previously unknown boreal larch forest that has developed on permafrost soils. Our general objectives were to constrain the use of Zn as a proxy for other metals in biogeochemical processes in Siberian larch forest and to test the possibility that Zn is a limiting micronutrient using a stable isotope approach . The specific questions to be addressed in this study were as follows: 1) To what degree are Zn isotopes fractionated between soil and tree in the permafrost environment? 2) Does this fractionation occur during soil – plant transfer, or within the plant itself? 3) How do the properties of contrasting habitats (in terms of nutrients availability, soil moisture and plant productivity) affect Zn isotopes fractionation? Finally, 4) Is the pool of Zn mobilized by the tree in these contrasted environments changing during the growing season (due to thawing)? We conducted our study within the Tura pilot-site (Central Siberia) to address these questions. This site is ideally suited for isotopic research in terrestrial biogeochemistry of permafrost forested environments under climate warming scenario [16,42-48]. As a first working hypothesis, and following results of previous studies, we assume that Zn uptake by plants growing in the permafrost environment is essentially controlled by Zn availability in the frozen soil. At the beginning of the active season, when the soil is frozen, Zn concentration and isotopic signature in the tree will be determined by available stocks of Zn within the plant. Upon progressive thawing of soils, deeper horizons may provide Zn to the roots. The second working hypothesis is that one can use the difference in Zn geochemistry between the slopes of southern and northern exposition, exhibiting dramatically different temperature and soil depth conditions, as a proxy for Zn isotope evolution in ecosystem of Central Siberia under on-going climate change.
Plot, stand (2006) and soil characteristics of the Tura site (64°19' N, 100°15'E) in Central Siberia, modified after Prokushkin et al. 
Stand density (tree/ha)
Active layer thickness** (cm)
Organic layer*** (cm)
Max rooting depth in mineral layer (cm)
South-facing slope (SF)
18 ± 8
86 ± 15
112 ± 26
123 ± 18
7 ± 2
North-facing slope (NF)
1 ± 1
25 ± 12
36 ± 11
42 ± 10
19 ± 5
Sphagnum bog (B)
1 ± 2
15 ± 8
20 ± 5
38 ± 8
Elemental concentrations (Zn, C, N) and Zn isotopic compositions for the different samples from the Kulingdakan catchment (Russia)
[Zn] 1 μg/g
[C] 2 %
[N] 2 %
δ 66 Zn(‰) 3
δ 66 Zn(‰) 4
δ 67 Zn(‰) 3
δ 68 Zn(‰) 3
0.27 ± 0.03
0.44 ± 0.04
0.55 ± 0.06
0.18 ± 0.06
0.25 ± 0.07
0.34 ± 0.12
0.24 ± 0.02
0.37 ± 0.03
0.47 ± 0.04
0.20 ± 0.05
0.30 ± 0.09
0.41 ± 0.08
0.21 ± 0.02
0.32 ± 0.06
0.41 ± 0.04
Sphagnum fuscum mosses
T1 (0–5 cm)
T1 (10–20 cm)
T1 (20–27 cm)
T1 (28–32 cm)
T2 (0–5 cm)
T2 (5–25 cm)
T2 (25–30 cm)
T2 (30–32 cm)
Mosses (Pleurozium schreberi), litter and organic layer
Live (green) mosses (0–3 cm)
0.12 ± 0.06
0.36 ± 0.09
0.35 ± 0.12
Dead (brown) mosses (3–6 cm)
Litter (7–10 cm)
0.41 ± 0.04
0.72 ± 0.06
0.87 ± 0.07
Organic layer (11–13 cm)
Elemental and isotopical composition of larch organs
Environment and date
[Zn] 1 μg/g
[C] 2 %
[N] 2 %
δ 66 Zn(‰) 3
δ 66 Zn(‰) 4
δ 67 Zn(‰) 3
δ 68 Zn(‰) 3
NF, June 7
38.4 ± 1.9
48.2 ± 0.8
2.66 ± 0.22
0.06 ± 0.05
0.17 ± 0.09
0.17 ± 0.08
NF, July 18
21.1 ± 6.9
48.6 ± 0.3
1.13 ± 0.09
NF, August 28
16.1 ± 3.6
48.7 ± 0.2
1.01 ± 0.12
0.22 ± 0.01
0.39 ± 0.01
0.48 ± 0.09
NF, Sept 10
11.8 ± 2.0
48.7 ± 1.2
0.23 ± 0.01
0.27 ± 0.08
0.37 ± 0.14
0.46 ± 0.18
SF, June 7
36.7 ± 13.6
47.8 ± 0.4
3.06 ± 0.40
0.20 ± 0.01
0.37 ± 0.04
0.44 ± 0.04
SF, July 18
12.6 ± 5.4
48.6 ± 0.5
1.47 ± 0.10
SF, August 28
14.5 ± 4.8
49.1 ± 0.4
1.21 ± 0.17
0.37 ± 0.03
0.59 ± 0.06
0.75 ± 0.07
SF, Sept 10
13.3 ± 3.7
48.7 ± 0.8
0.31 ± 0.11
0.28 ± 0.06
0.43 ± 0.10
0.56 ± 0.12
PB, June 7
36.1 ± 5.7
48.4 ± 0.5
2.23 ± 0.44
−0.21 ± 0.10
−0.22 ± 0.16
−0.35 ± 0.20
PB, July 18
20.9 ± 2.4
48.2 ± 0.8
1.27 ± 0.10
PB, August 28
13.3 ± 1.3
49.8 ± 0.3
1.04 ± 0.05
−0.15 ± 0.05
−0.16 ± 0.07
−0.24 ± 0.09
PB, Sept 10
18.2 ± 7.0
48.8 ± 0.7
0.26 ± 0.02
0.13 ± 0.09
0.21 ± 0.12
0.26 ± 0.18
SF, Sept 10
7.0 ± 2.3
48.4 ± 1.9
0.08 ± 0.01
Wood (bulk) 5
SF, Sept 10
4.8 ± 3.2
47.1 ± 3.3
0.10 ± 0.05
SF, Sept 10
PB, Sept 10
7.7 ± 2.4
43.9 ± 4.4
PB, Sept 10
15.8 ± 16.6
44.3 ± 3.4
0.07 ± 0.02
Wood (bulk) 5
PB, Sept 10
10.0 ± 7.1
47.1 ± 3.3
0.09 ± 0.02
Fine roots (<2 mm) n = 5
SF, Sept 10
riparian zone, Sept 10
In September, all 3 trees from each plot were harvested to determine their above ground biomass. Stem discs (3 cm thick) were taken from every tree at breast height (1.3 m) to analyze tree ages and then different portions (parts) of the disc (bark, sapwood, heartwood and integral wood samples) originating from the South-facing slope and Sphagnum peatbog environments were taken for Zn concentration and isotopes analysis. Again Zn concentrations were obtained on the analysis of three samples while Zn isotopic composition was obtained on the composited sample.
In addition, in the slopes we have collected 5 replicate columns (0–13 cm depth) of surficial organic material divided into moss, Pleurozium schreberii, (separately live and dead portions), litter and soil organic layer. Replicates were composited to one sample. The Zn concentrations and isotopic compositions were measured on the composited samples. Within the Sphagnum peatbog site we sampled 2 peat columns (0–32 cm depths), which were divided into live Sphagnum fuscum moss layer (0–5 cm) and three layers of peat at different stages of decomposition (between 5 and 32 cm, see Tables 2 and 3). The Zn concentrations and isotopic compositions were measured on these samples.
Soil samples were collected just beneath the sampled organic columns in the South-facing slope, North-facing slope and Sphagnum peat bog in mid-August of 2006, at the period of maximum active layer thickness. Mineral soil sampling was performed by 100 cm3 cylinder for the entire active layer to the permafrost (0 to 120 cm). Further in the text, soil samples are indicated by the locality identifier followed by the sampling depth in cm (e.g., North0-10). In this example, "North" is the aspect of the soil sampling site and 0–10 indicates that the sample was collected from the 0 to 10 cm soil layer. Soil samples were passed through a 2 mm sieve and air-dried on site. Back at the laboratory, samples were dried at 80°C for 24 hours and finely ground using an agate mixer mill (Retsch, Germany). Given the extremely high homogeneity of the soil profile and the similarity between the unaltered basalt and the deepest (C) soil horizon [42,45,47,49,52], the mineral soil horizon sampled below the rooting depth, 60–100 cm at the South-facing slope, can be considered as unaltered mineral substrate. The bog site has no clear plant-related mineral substrate as the peat is permanently frozen below the rooting depth.
Interstitial soil solutions were collected both from N- and S-facing slope at different depths using suction cups and plate lysimeters, respectively, as described elsewhere .
Finally, the water samples from the large river (Kochechum) that drains this region [46,49] was collected in August 2006, filtered on site using trace clean techniques (a pre-washed Nalgene disposable filter unit with 0.22 μm acetate cellulose filter) and acidified with doubly distilled HNO3. An untreated MilliQ water sample was filtered to produce field blanks for analysis.
Sample preparation and Zn concentration measurement
To measure the elemental concentrations, all samples were digested in Teflon vials within individual polycarbonate compartments (A100) containing Teflon hot plates in a clean room (ISO3). Between 100 and 200 mg of plant and soil was first reacted with hydrogen peroxide (H2O2) for 24 hours at ambient temperature and further digested in HNO3 + HF, for 36 hours at 80°C, then in HCl for 36 hours at 80°C, and finally treated with HCl-HNO3 for 36 hours at 80°C.
Zn concentrations were measured by ICP-MS (Agilent 7500ce). Indium and rhenium were used as internal standards to correct for instrumental drift and eventual matrix effects. The international geostandards of basaltic rock BE-N (from CRPG, France) and lichens BCR-CRM 482 (from BCR, Belgium) were used to check the accuracy of the analytical methodology. The relative standard deviation between the certified values and our measurements were expressed as ([Zn]certified - [Zn]measured)/(([Zn]certified + [Zn]measured)/2)*100, and they were lower than 10% (1 s.d., n = 10). The agreement between our Zn analyses in the SLRS-4 standard of natural water and the certified value was better than 15% (2 s.d., n = 20). For isotopic analysis of the river water sample, large volume of filtered solution (~700 mL) was evaporated to dryness in the clean room under the hood box (A100) and the residue was subjected to acid leaching as for larch needles.
Zn separation and isotopic measurements
Zn was separated from the matrix using the AG-MP1 anion exchange procedure following Maréchal et al.  inside a laminar hood box ISO3 placed in the clean room. Two successive separations were performed for each sample. An aliquot of sample containing 800 ng was loaded on the column. Column loading was limited to <15% of the total anion exchange capacity of the resin. After the first separation procedure Zn fraction was loaded on a new column to perform a second separation (purification). Column yields was checked for each sample by analyzing Zn concentrations by quadrupole ICP-MS. Yields were found to be 95 ± 11% for Zn. Two samples presenting bad yields (<70%) were discarded and not considered in this study. Zn was removed using 0.05 M HNO3, the matrix used during MC-ICP-MS measurements. Note that this matrix is different from that used by Maréchal et al.  Total Zn blank due to the whole procedure separation are <8 ng which is insignificant relative to the total amounts loaded in the columns.
Zn isotopic composition was measured on a Thermo Finnigan NEPTUNE MC-ICP-MS hosted in GET (Toulouse). The Neptune introduction system used consists of a tandem quartz glass spray chamber (cyclone + standard Scott double pass) coupled with a low flow PFA nebulizer (50 to μL/min). Ar gas flows of ~15 L/min, ~1.2 L/min, and ~0.6 L/min were used for coolant, nebulizer and auxiliary, respectively. For ICP operation, an RF generator power of 1300 W was used. The beams were collected on a multi collector module with eight movable and one fixed Faraday cup displaying a maximum relative mass range of 17%. This cup configuration (used for Zn isotopic measurements) allowed simultaneous collection of the four Zn masses (64, 66, 67, and 68) 63Cu and 65Cu for mass bias correction and mass 62 to check for eventual Ni interferences . After a minimum of 3 hours of instrument warm up and setting, analytical sequences were run automatically using a CETAC ASX-100 autosampler.
Between each sample and standard the machine was rinsed with 0.05 N HNO3 from two different vials for 1 minute each. Blank measurements consist of 1 block of 10 cycles (8 s) and sample and standard measurements consist of 2 blocks of 20 cycles of 8 seconds each. For sample and standard measurements the internal precision was found to be in the range 5–10 ppm (2σ err.) for both 65Cu/63Cu and 66Zn/64Zn. Instrumental mass fractionation was corrected by using two different methods: bracketing and internal standard [53,54]. A Cu standard (NIST 976) was added to the purified Zn fractions and a Cu (NIST 976) + Zn (JMC 3-0749 L) standard mixture was run as a brack standard, thus allowing to perform different corrections. The 62Ni signal was measured to evaluate the isobaric interferences of 64Ni on the 64Zn and no significant interferences were found in the samples. During one analytical session we checked that Cu and Zn fractionation factors remained constant. The Zn isotopic results in this study are given in the recommended delta notation for the 66Zn/64Zn ratio: δ66Zn = [((66Zn/64Zn)s / (66Zn/64Zn)JMC)-1]*1000 (in ‰) where s stands for sample and JMC represents the Zn isotopic standard solution (JMC 3–0749 L) . Note that this solution is not a referenced material but an elemental standard solution used by several laboratories . The δ66Zn has been also expressed relative to the new certified reference material (IRMM-3702) for Zn isotopes measurements (see Tables 2 and 3) .
Mass dependent isotopic fractionation, analytical control and error evaluation
During the samples preparation, two reference materials (BCR-1 and CRM 482) were involved to check the external reproducibility. Our data (δ66Zn = 0.20 ± 0.07 for BCR-1 and δ66Zn = 0.10 ± 0.08) are in agreement with published values [21,22,57]. In the tables of results an error (standard deviation) is given for some of the samples when several replicates have been performed. Further details on the analytical precision of δ66Zn on standards during the analytical sessions, the analytical repeatability of a Zn monoelemental standard and procedural repeatability of the CRM standard as well as the long term measurement of Zn standards in our laboratory are presented elsewhere .
Zn concentration and isotopic composition in soils and Zn concentration in soil solutions
The Zn concentrations and isotopic compositions of 5 selected soil samples are given in Table 2. We observed that both Zn concentrations (108 to 121 μg/g) and δ66Zn (0.18 to 0.27‰) are similar among different environments (south and north-facing slopes, and the bog) over the full depth of the soil column, except the upper mineral horizon of the north-facing slope exhibiting ca. < 0.1‰ heavier isotopic composition. The river water had a Zn isotopic composition (δ66Zn = 0.16 ± 0.05‰) similar to that of the measured soils.
Zn concentration and isotopic composition in mosses (Sphagnum fuscum and Pleurozium schreberi)
Samples from two microenvironments (trough and mound) of the Sphagnum fuscum peat bog (T1 and T2) were taken at different depths, ranging from the surface to 32 cm (Table 2). Zn concentrations ranged between 14 and 21 μg/g without any visible depth-related trend. In the T1 (trough) location, the peat Zn concentration was close to 20 μg/g in the most superficial sample and approximately 14 to 15 μg/g in the samples collected between 5 and 27 cm depth. In the T2 peat mound location, the Zn concentration was approximately 15 μg/g in the most superficial sample and between 14 and 21 μg/g in the other samples collected at greater depths. Again, we did not observe any relationship between Zn isotopic composition and Zn concentration in peat bog samples. However, δ66Zn increases steadily from 0.15-0.2‰ in the most superficial samples to around 0.6‰ at 32 cm depth.
The green (alive) and brown (dead) parts of Pleurozium schreberi moss exhibited a δ66Zn value close to 0‰. In comparison, litter and organic layer horizons developed beneath the moss exhibited distinctly heavier isotopic compositions from around 0.4‰, similar to the deeper layers of the sphagnum peat column.
Zn concentration and isotopic composition in larch needles and wood
Comparing the isotopic signature of larch needles between the three different environments revealed a distinct difference between the bog area and valley slopes (Table 3 and Figure 5). Larch needles collected in August at the North and South-facing slopes exhibited δ66Zn values (0.22 to 0.37‰) close to those of the soils (0.18 to 0.27‰). In contrast, larch needles from the peat bog area exhibited much lower isotopic signatures (−0.21 to 0.13‰), which were especially pronounced in August.
Zn isotopic composition in soils
The isotopic composition of soils from Central Siberia is notably constant with a δ66Zn value approximately 0.2‰. Although the isotopic composition of the basaltic rock constituting the geological substratum in Central Siberia is not available, the values of various basalts are highly homogeneous (δ66Zn = 0.20 ± 0.09‰ for BCR-1 standard ; δ66Zn = 0.24 ± 0.01‰ for 5 basaltic rock samples ; δ66Zn = 0.23‰ for BCR-2 standard ) and similar to those of the Tura bottom soil horizon (0.20 ± 0.05‰).
Results from the Tura site in Central Siberia indicate that there is no measurable Zn isotopes fractionation between parent rocks and soils assuming that i) the parent bedrock is similar to the C horizon of the soil, as it follows from other isotopical [45,47] and chemical  data, and ii) the Siberian basalts exhibit the same Zn isotopic signature as other basalts worldwide.
A homogeneous Zn isotopic composition over the full depth of the soil column corroborates previous studies on basalt chemical weathering in this region based on chemical and mineralogical analyses of rocks and soils and fluid solution chemistry [49,52]. These authors demonstrated that the weathering of basaltic rocks does not lead to chemical and mineralogical differentiation along the soil profiles under the current ecosystem state. As such, the fractionation of Zn isotopes in the soil profile is also minimal. The most likely mechanism for such a homogeneous composition is seasonal freezing front migration, downward in the beginning of the active season and both downward and upward in the autumn. This intensive mechanical mixing of the frozen/unfrozen soil horizon is further enhanced by sliding of the unfrozen soil layers along the ice boundary, especially visible at the N-facing slope. The observation that the whole soil profile is involved in seasonal cycles of freezing-thawing is also supported by high variability of Zn concentration in the interstitial soil solution in the course of the active season, sampled at different depths of the mineral layer, from 5 to 40 cm and from 5 to 60 cm on the NF and SF sites, respectively (Figure 3). Noteworthy a decrease of Zn concentration at the end of the active season, in September, observed on both environments. The averaged over all seasons Zn concentration in soil solution demonstrates the enrichment in surface horizons (5 cm) on the S-facing slope but virtually identical concentration over the full soil column in the N-facing slope, consistent with much higher water content of the NF habitats  and possible impact of the plant litter in productive SF habitats on Zn enrichment in the most surface soil horizon.
Zn isotopic composition in mosses
Surface samples of live (green) mosses from Sphagnum fuscum peat bogs exhibited much lighter isotopic compositions compared to deep (brown) samples representing decomposed biomass. For the Pleurozium schreberi moss column, heavy isotope enrichment was also recorded in the litter and organic layer relative to the live and dead biomass of surface moss samples. Because there is no local contamination from an industrial complex, road or urban areas in the vicinity of our site, the δ66Zn values of the moss surface layers (from 0.14 to 0.20‰) may represent the background atmospheric deposition for this region . Indeed, various studies have shown that mosses can be considered as a time-integrated measure of both wet and dry metal deposition from the atmosphere [60,61]. The δ66Zn values obtained from the bottom layers of the peat core (0.5 to 0.7‰) are significantly higher than the underlying mineral substrate. This heavy isotope enrichment may be related to the decay of organic matter induced by bacteria and fungi as may follow from Weiss et al. , although these authors did not evidenced this directly. A similar fractionation tendency towards heavy isotope enrichment in bottom horizons has been reported for carbon, nitrogen and sulfur in sphagnum peat bogs from a pristine area of the Czech Republic . According to the latter authors, various other processes, such as nutrient recycling, ion exchange and absorption, may be involved in changing the initial isotope ratio.
Even if the main processes controlling isotope fractionation could be different between major nutrients and Zn, the retention of heavy isotopes during the humification process appears to be consistent with recent studies on Zn complexation, adsorption, and incorporation within biological and organic compartments [39,64-66]. Jouvin et al.  showed that the complexation of Zn by high molecular weight organic compounds (i.e., humic and fulvic acids) favors the heavy isotopes. Similarly, adsorption onto biological surfaces leads to cell enrichment with heavy isotopes as a consequence of the change in Zn coordination numbers and bond distances . During humification process the labile products such as carbohydrates and amino acids are destroyed while aromatic structures appear [67,68]. Within that framework, and considering that chemical bonds Zn-O-C are stronger within the aromatic structures, in general accord with previous experimental and modeling works  we suggest that humification of both dead moss and larch needles leads to the enrichment of residual fraction by heavier isotope while releasing the lighter isotope in soil solution although direct field verification of this hypothesis is still lacking.
Zn isotopic fractionation between soil and plants in different habitats
A larch root sample exhibited 0.2-0.4‰ heavier isotopic composition relative to the needles. Despite limited sampling, this result is in agreement with previous studies showing that roots are preferentially enriched in heavier isotopes relative to plant tissue following the conceptual model for Cu and Zn isotope fractionation during uptake at the root surface (Jouvin et al.  and references therein). According to these studies it is mostly agreed that i) within the root epidermal cell, Zn2+ ions and its organic and inorganic complexes can move via non-specific transport pathways favoring light isotopes; and ii) transport via specific Zinc/Iron regulated Proteins (ZIP) or phytosiderophores (PS) favors heavy isotopes . A possible reason for the 66Zn enrichment in the roots is Zn-phosphate binding as supported by ab-initio calculations .
If we consider the root-shoot translocation, our results corroborate previous works showing heavy isotope depletion from root to shoot . The translocation of metal by the xylem from roots to shoots occurs in the form of free Zn2+ and its complexes with organic acids (amino acids and peptides) when produced by the plant itself. The diffusion of these Zn species induced by transpiration and by exchange with cell wall binding sites has been found to promote light isotopes migration along the transport path [12,71].
Because the stem accounts for between ~80 and ~92% (depending on the habitat) of the above-ground larch biomass, the isotopic signature of the bulk wood approximates the isotopic composition of the whole larch above-ground biomass. Our results (section 4.3) demonstrate that the larch bulk wood collected from the South-facing slope is ~ 0.2‰ heavier than the nutritive topsoil and the bulk wood of the larch from a sphagnum peat bog is 0.2‰ lighter than the Sphagnum peat bog organic substrate of the rooting zone in this habitat. Therefore, there is a slight but measurable fractionation between the larch tree and the nutrient reservoir (mineral or peat soil). Various studies have shown that the whole vegetation could be enriched in light or heavy isotopes relative to the soil substrate depending on the element (e.g., Ca, Mg and Zn), the physiological status of the plant and the environment [20,33,72-78]. The trends observed for the larch in the South-facing slope in this study are consistent with those reported for tropical watershed , namely the plant enrichment in heavier Zn isotopes by 0.2-0.4‰ relative to the soil reservoir.
Interestingly, the single river water sample exhibited an isotopic composition similar to that of the larch needles, upper organic layers of the peat bog and the upper mineral soil horizons (0.2 ± 0.05‰). The isotopic measurements thus did not allow for discrimination between mineral and organic sources of Zn in surface waters within a 0.1‰ resolution. Curiously, this result corroborates recent conclusions achieved for Ca47, Si45 and Mg48 in this region; the stable isotope approach does not allow for a straightforward quantification of the relative contribution of plant litter, soil organic and mineral horizons and the bedrocks to the riverine fluxes of elements.
Evolution of δ66Zn in the course of vegetative season
A recent study of nutrients and non-essential elements in larch needles during the vegetative season of the same boreal forest plot in Central Siberia  revealed that Zn micronutrient parallels the behavior of such macronutrients as nitrogen (N), phosphorus (P) or potassium (K) with decreasing concentration from June to September (see Figure 5). High concentrations of all these elements at the beginning of the growing period (June) are indicators of intense metabolic processes such as photosynthetic activity and protein synthesis when needles are first produced. The observed decrease of macro and micronutrient concentrations such as Zn between June and July can be explained by a dilution caused by C rather than Zn accumulation in the leafs whereas the decrease of Zn concentration between August and September results from a resorption . The latter resorption process is a key mechanism for deciduous plants to avoid essential nutrient losses due to the litterfall. In fact, these elements can be mobilized towards the phloem in September at the end of the active period and before senescence and redistributed the next year to the younger tissues by the xylem.
Note that the Zn content in June should originate from the pool of metal accumulated at the end of the previous seasons as soil is still frozen in June. In the end of summer, the bulk isotopic signature can be considered as a mixture between a Zn pool stored in shoots from the previous year and a pool of Zn taken from the progressively thawing soils.
A factor of 2 higher Zn concentrations in solution of SF soil compared to that of the NF slope (Figure 4) suggests more intense cycling of Zn in the warmer habitats. In deeper layers (BC/C, up to 60 cm depth which is close to permafrost) Zn concentrations in soil solution decrease to ca. 4 μg/L at both slope soils. Considering that all DOC of soil solution corresponds to fulvic acid we calculated the isotopic composition of the various Zn species when the DOC concentration increases from 10 to 50 mg/L with pH ranging from 3 to 6, encompassing the majority of soil solution composition. This calculation was performed for a constant Zn concentration of 10 μg/L considering that at the equilibrium the organic matter - complexed Zn will be 0.2‰ heavier than the free Zn2+ (see Jouvin et al. ). This first-order calculation shows that in acidic environment, the isotopic composition of the free Zn available for plants will be up to 0.1‰ lighter for a DOC concentration of 50 mg/L compared to 10 mg/L. This result is consistent with the measured difference of isotopic composition between the needles of the North-facing slope and the bog area.
The possible role of resorption processes on Zn isotope fractionation between June and September can be assessed from temporal evolution of Zn concentration in larch needles from various habitats. Zn concentration in the needles remained constant at the South-facing slope and bog and decreased between August and September at the North-facing slope. The carbon concentration in the needles decreased during the same period at the South-facing slope and bog or remains constant at the North-facing slope (Table 1) . It follows that Zn resorption process involved between 10 to 25% of total Zn stock in the plant . These values imply that for the resorption processes to influence the isotopic signatures at the level of 0.1 to 0.2‰, the δ66Zn isotopic shift corresponding to 10 to 25% of resorption should be on the order of 1.0 to 1.5‰, which is half a permil higher than the maximum Zn isotopic shifts reported during single-reaction adsorption or complexation processes [64,70]. As such, Rayleigh distillation processes should be invoked to explain isotopic fractionation by Zn resorption; however this could not be tested in this study.
Therefore, the most likely mechanism of Zn isotopic composition difference in the larch needles at the beginning and at the end of the vegetative season is the change of rooting depth in soil horizons upon soil thaw during active season, with rather uniform isotopic signature among contrasting habitats. At the same time, the obtained data did not allow to confirm the possibility that Zn is a limiting micronutrient for the larch. On slopes of contrasting nutritional status, the re-translocation affects Zn isotopic composition in needles in a variable degree, often comparable with the resolution of isotopic measurements.
Implying the previously developed approach of substituting space for time and assuming that the contemporary difference in local scales within the climate gradient created by slope orientation can serve as a proxy for future changes of a given system,42,45,82 results of the present study do not suggest any significant change in Zn isotopic pool of terrestrial biome in the Central Siberia at the most extreme scenario of climate warming in this region, i.e., in the case of complete replacement of N-facing cold and wet slope conditions by S-facing warm and dry habitats.
This study revealed that Zn isotopes are not significantly fractionated within the permafrost-dominated soils of this site in Central Siberia. In contrast, measurable isotope fractionations were observed between the above-ground larch biomass and the soil reservoir. The isotopic composition (δ66Zn) of the larch needles depends on the habitat (North and South-facing slope or bog area) with 0.2-0.4‰ heavier values in the North- and South-facing slope relative to the peat bog. This trend is not related to a different physiological response of the tree but rather to the isotopic composition of the solution available to the rooting zone. It was hypothesized that stronger Zn complexation by organic ligands and humification products in bog areas leaves a soil solution which is isotopically lighter than the bulk mineral horizon.
Zn isotopes help to trace which free Zn pool is mobilized by the tree via following the rooting depth and soil solution chemistry in contrasted local environments and how it changes during the soil thawing over growing season. Specifically, during the course of the vegetation we observed enrichment in Zn heavy isotopes of larch needles. We suggest that the increased rooting depth and the decreased DOC concentration in the root uptake zone resulting from progressively thawing soil lead to heavy isotopes becoming more and more available for the larch roots in the course of the vegetative season. This is consistent with physico-chemical calculations showing that the decrease of DOC will facilitate the uptake of heavy isotope as it will be less retained in strong organic complexes. With the extreme climate evolution in Central Siberia, the North-facing slopes may experience climatic variations that are typical of present-day-South-facing slopes but it is unlikely that any significant change of Zn isotopic composition occurs in the soil-tree-river system.
This work was supported by the GDRI Car-Wet-Sib, the ANR CESA “Artic Metals” (France), RSF 14-24-0013, and the grant “BIO-GEO-CLIM” of the Russian Ministry of Education and Research and Tomsk State University (N° 14.B25.31.0001).
- Schuur EAG, Bockheim J, Canadell JG, Euskirchen E, Field CB, Goryachkin SV, et al. Vulnerability of permafrost carbon to climate change: implications for the global carbon cycle. BioScience. 2008;58(8):701–14.View ArticleGoogle Scholar
- Schuur EAG, McGuire AD, Schädel C, Grosse G, Harden JW, Hayes DJ, et al. Climate change and the permafrost carbon feedback. Nature. 2015;520:171–9.View ArticleGoogle Scholar
- Linderholm HW. Growing season changes in the last century. Agricultural Forest Meteorol. 2006;137:1–14.View ArticleGoogle Scholar
- McGuire AD, Anderson LG, Christensen TR, Dallimore S, Guo L, Hayes DJ, et al. Sensitivity of the carbon cycle in the Arctic to climate change. Ecological Monographs. 2009;79(4):523–55.View ArticleGoogle Scholar
- Kim Y, Kimball JS, Zhang K, McDonald KC. Satellite detection of increasing Northern Hemisphere non-frozen seasons from 1979 to 2008: Implications for regional vegetation growth. Remote Sens Environ. 2012;121:472–87.View ArticleGoogle Scholar
- Serreze MC, Walsh JE, Chapin III FS, Osterkamp T, Dyurgerov M, Romanovsky V, et al. Observational evidence of recent change in the northern high-latitude environment. Climatic Change. 2000;46:159–207.View ArticleGoogle Scholar
- Holmes RM, McClelland JW, Peterson BJ, Tank SE. Seasonal and annual fluxes of nutrients and organic matter from large rivers to the Arctic Ocean and surrounding seas. Estuaries and Coasts. 2012;35:369–82.View ArticleGoogle Scholar
- Pokrovsky OS, Viers J, Dupré B, Chabaux F, Gaillardet J, Audry S, et al. Biogeochemistry of carbon, major and trace elements in watersheds of northern Eurasia drained to the Artic Ocean: the change of fluxes, sources and mechanisms under the climate warming prospective. Compte Rendus Geoscience. 2012;344:663–77.View ArticleGoogle Scholar
- Houghton RA, Skole DL. Carbon. In: Turner BL, Clark WC, Kates RW, Richards JF, Mathews JT, Meyer WB, editors. The Earth as Transformed by Human Action. Cambridge: Cambridge University Press; 1990. p. 393–408.Google Scholar
- Botch MS, Kobak KI, Vinson TS, Kolchugina TP. Carbon pools and accumulation in peatlands of the former Soviet Union. Global Biogeochem Cycles. 1995;9:37–46.View ArticleGoogle Scholar
- Jarvis PG, Saugier B, Detlef-Schulze E. Productivity of boreal forests. In: Ray J, Saugier B, Mooney HA, editors. Terrestrial Global Productivity. Amsterdam: Academic Press, Elsevier; 2001. p. 211–44.View ArticleGoogle Scholar
- Palmer CM, Guerinot ML. Facing the challenges of Cu, Fe and Zn homeostasis in plants. Nature Chem Biol. 2009;5:333–40.View ArticleGoogle Scholar
- Broadley MR, White PJ, Hammond JP, Zelko I, Lux A. Zinc in plants. New Phytol. 2007;173(4):677–702.View ArticleGoogle Scholar
- Aerts R, Chapin III FS. The mineral nutrition of wild plants revisited: a re-evaluation of processes and patterns. Adv Ecol Res. 2000;30:1–67.View ArticleGoogle Scholar
- Hobbie SE, Vitousek PM. Nutrient limitation of decomposition in Hawaiian forests. Ecology. 2000;81:1867–77.View ArticleGoogle Scholar
- Prokushkin AS, Knorre AA, Kirdyanov AV, Schulze ED. Productivity of mosses and organic matter accumulation in the litter of sphagnum larch forest in the permafrost zone. Russ J Ecol. 2006;37(4):225–32.View ArticleGoogle Scholar
- Wang YP, Law RM, Pak B. A global model of carbon, nitrogen and phosphorus cycles for the terrestrial biosphere. Biogeosci Disc. 2009;6:9891–944.View ArticleGoogle Scholar
- Tokuchi N, Hirobe MK, Arai H, Hobara S, Fukushima K, Matsuura Y. Soil Nitrogen Dynamics in Larch Ecosystem. In: Osawa A, Zyryanova OA, Matsuura Y, Kajimoto T, Wein RW, editors. Permafrost Ecosystems: Siberian Larch Forests, vol. 209. New York: Springer Ecological series; 2010. p. 229–43.View ArticleGoogle Scholar
- Eckstein RL, Karlsson PS, Weih M. Leaf life span and nutrient resorption as determinants of plant nutrient conservation in temperate-arctic regions. New Phytol. 1999;143:177–89.View ArticleGoogle Scholar
- Weiss DJ, Mason TFD, Zhao FJ, Kirk GJD, Coles BJ, Horstwood MSA. Isotopic discrimination of zinc in higher plants. New Phytol. 2005;165:703–10.View ArticleGoogle Scholar
- Cloquet C, Carignan J, Libourel G. Isotopic composition of Zn and Pb atmospheric depositions in an urban/periurban area of Northeastern France. Environ Sci Technol. 2006;40:6594–600.View ArticleGoogle Scholar
- Viers J, Oliva P, Sonke J, Nonell A, Freydier R, Gainville R, et al. Evidence of Zn isotopic fractionation in a soil-plant-system of a pristine tropical watershed (Nsimi, south Cameroon). Chem Geol. 2007;239:124–37.View ArticleGoogle Scholar
- Sivry Y, Riotte J, Sonke J, Audry S, Schäfer J, Viers J, et al. Zn isotopes as tracers on anthopogenic pollution from Zn-ore smelters. The case of the Riou Mort-Lot River system. Chem Geol. 2008;255:295–304.View ArticleGoogle Scholar
- Sonke JE, Sivry Y, Viers J, Freydier R, DeJonghe L, Andre L, et al. Historical variations in zinc stable isotope composition of smelter polluted sediments. Chem Geol. 2008;252:145–57.View ArticleGoogle Scholar
- Borrok DM, Nimick DA, Wanty RB, Ridley WI. Isotopic variations of dissolved copper and zinc in stream waters affected by historical mining. Geochim Cosmochim Acta. 2008;72:329–44.View ArticleGoogle Scholar
- Borrok DM, Wanty RB, Ridley WI, Lamothe PJ, Kimball BA, Verplanck PL, et al. Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed. Appl Geochem. 2009;24:1270–7.View ArticleGoogle Scholar
- Fernandez A, Borrok DM. Fractionation of Cu, Fe and Zn isotopes during the oxidative weathering of sulfide-rich rocks. Chem Geol. 2009;264:1–12.View ArticleGoogle Scholar
- Moynier F, Pichat S, Pons ML, Fike D, Balter V, Albarède F. Isotopic fractionation and transport mechanisms of Zn in plants. Chem Geol. 2009;267:125–30.View ArticleGoogle Scholar
- Mattielli N, Petit JCJ, Deboudt K, Flament F, Perdrix E, Taillez A, et al. Zn isotope study of atmospheric emissions and dry depositions within a 5 km radius of a Pb–Zn refinery. Atmos Environ. 2009;43(6):1265–72.View ArticleGoogle Scholar
- Arnold T, Kirk GJD, Wissuwa M, Frei M, Zhao FJ, Mason TFD, et al. Evidence for the mechanisms of zinc uptake by rice using isotopic fractionation. Plant Cell Environ. 2010;33:370–81.View ArticleGoogle Scholar
- Bigalke M, Weyer S, Kobza J, Wilcke W. Stable Cu and Zn isotope ratios as tracers of sources and transport of Cu and Zn in contaminated soil. Geochim Cosmochim Acta. 2010;74(23):6801–13.View ArticleGoogle Scholar
- Caldelas C, Dong S, Araus JL, Weiss DJ. Zinc isotopic fractionation in Phragmites australis in response to toxic level of zinc. J Exper Bot. 2011;62(6):2169–78.View ArticleGoogle Scholar
- Aucour AM, Pichat S, Macnair MR, Oger P. Fractionation of stable zinc isotopes in the zinc hyperaccumulator Arabidopsis halleri and nonaccumulator Arabidopsis petraea. Environ Sci Technol. 2011;45:9212–7.View ArticleGoogle Scholar
- Tang YT, Cloquet C, Sterckeman T, Echevarria G, Carignan J, Qiu RL, et al. Fractionation of stable zinc isotopes in the field-grown zinc hyperaccumulator Noccaea caerulescens and the zinc-tolerant plant Silene vulgaris. Environ Sci Technol. 2012;46(18):9972–9.Google Scholar
- Jouvin D, Weiss DJ, Mason TFM, Bravin MN, Louvat P, Zhao F, et al. Stable isotopes of Cu and Zn in higher plants: evidence for Cu reduction at the root surface and two conceptual models for isotopic fractionation processes. Environ Sci Technol. 2012;46:2652–60.View ArticleGoogle Scholar
- Aranda S, Borrok DM, Wanty RB, Balistrieri LS. Zinc isotope investigation of surface and pore waters in a mountain watershed impacted by acid rock drainage. Sci Total Environ. 2012;410:202–13.View ArticleGoogle Scholar
- Dong S, Weiss DJ, Strekopytov S, Kreissig K, Sun Y, Baker AR, et al. Stable isotope ratio measurements of Cu and Zn in mineral dust (bulk and size fractions) from the Taklimakan Desert and the Sahel and in aerosols from the eastern tropical North Atlantic Ocean. Talanta. 2013;114:103–9.View ArticleGoogle Scholar
- Zhao Y, Vance D, Abouchami W, de Baar HJW. Biogeochemical cycling of zinc and its isotopes in the Southern Ocean. Geochim Cosmochim Acta. 2014;125:653–72.View ArticleGoogle Scholar
- Wanty RB, Podda F, De Giudici G, Cidu R, Lattanzi P. Zinc isotope and transition-element dynamics accompanying hydrozincite biomineralization in the Rio Naracauli; Sardinia, Italy. Chem Geol. 2013;337/338:1–10.View ArticleGoogle Scholar
- Weiss DJ, Boye K, Caldelas C, Fendorf S. Isotope fractionation during early dissolution of biotite granite. Soil Sci Soc Amer J. 2014;78:171–9.View ArticleGoogle Scholar
- Jaouen K, Pons ML, Balter V. Iron, copper and zinc isotopic fractionation up mammal trophic chains. Earth Planet Sci Lett. 2013;274:164–72.View ArticleGoogle Scholar
- Viers J, Prokushkin AS, Pokrovsky OS, Kirdyanov AV, Beaulieu E, Zouiten C, et al. Seasonal and spatial variability of elemental concentrations in boreal forest larch foliage of Central Siberia on continuous permafrost. Biogeochemistry. 2013;113:435–49.View ArticleGoogle Scholar
- Prokushkin AS, Gleixner G, McDowell WH, Ruchlow S, Schulze ED. Source- and substrate-specific export of dissolved organic matter from permafrost-dominated forested watershed in central Siberia. Global Biogeochem Cycles. 2007;21:GB4003.View ArticleGoogle Scholar
- Prokushkin AS, Pokrovsky OS, Shirokova LS, Korets MA, Viers J, McDowell WH. Export of dissolved carbon from watersheds of the central Siberian plateau. Doklady Earth Sci. 2011;441:1568–71.View ArticleGoogle Scholar
- Pokrovsky OS, Reynolds BC, Prokushkin AS, Schott J, Viers J. Silicon isotope variations in Central Siberia rivers during basalt weathering in permafrost-dominated larch forests. Chem Geol. 2013;355:103–16.View ArticleGoogle Scholar
- Bagard ML, Chabaux F, Pokrovsky O, Viers J, Prokushkin A, Stille P, et al. Annual variations of the geochemical fluxes in rivers from high latitude permafrost areas (Kochechumo River and Nizhniya Tunguska River, Central Siberia). Geochim Cosmochim Acta. 2011;75:3335–57.View ArticleGoogle Scholar
- Bagard ML, Schmit AD, Chabaux F, Pokrovsky OS, Viers J, Stille P, et al. Biogeochemistry of stable Ca and radiogenic Sr isotopes in a larch-covered permafrost dominated watershed of Central Siberia. Geochim Cosmochim Acta. 2013;114:169–87.View ArticleGoogle Scholar
- Mavromatis VM, Prokushkin AS, Pokrovsky OS, Viers J, Koretz M. Mg isotopes in permafrost-dominated Central Siberian larch forest watersheds. Geochim Cosmochim Acta. 2014;147:76–89.View ArticleGoogle Scholar
- Pokrovsky OS, Schott J, Kudryavtzev DI, Dupré B. Basalt weathering in Central Siberia under permafrost conditions. Geochim Cosmochim Acta. 2005;69:5659–80.View ArticleGoogle Scholar
- Prokushkin SG, Bugaenko TN, Prokushkin AS, Shikunov VG. Succession-driven transformation of plant and soil cover on solifluction sites in the permafrost zone of Central Evenkia. Biol Bulletin. 2010;1:80–8.View ArticleGoogle Scholar
- Markert B. Sample preparation (cleaning, drying, homogenization) for trace element analysis in plant matrices. Sci Tot Environ. 1995;176:45–61.View ArticleGoogle Scholar
- Pokrovsky OS, Schott J, Dupré B. Trace elements fractionation and transport in boreal streams and soil solutions of basaltic terrain. Central Siberia Geochim Cosmochim Acta. 2006;70:3239–60.View ArticleGoogle Scholar
- Maréchal C, Télouk P, Albarède F. Precise analysis of copper and zinc isotopes. Geochim Cosmochim Acta. 1999;156:251–73.Google Scholar
- Maréchal CN, Nicolas E, Douchet C, Albarede F. Abundance of zinc isotopes as a marine biogeochemical tracer. Geochem Geophys Geosyst. 2000;1:1015. doi:10.1029/1999GC000029.Google Scholar
- Johnson CM, Beard BL, Albarède F. Geochemistry of nontraditional stable isotopes. Rev Mineral Geochem, Mineral Soc Amer. 2004;55:1–22. Washington, DC, USA.View ArticleGoogle Scholar
- Moeller K, Schoenberg R, Pedersen RB, Weiss D, Dong S. Calibration of the new certified reference material ERM-AE633 and ERM-AE647 for copper and IRMM-3702 for zinc isotope amount ratio determinations. Geostandards Geoanalytical Res. 2012;36/2:177–99.View ArticleGoogle Scholar
- Archer C, Vance D. Mass discrimination correction in multiple-collector plasma source mass spectrometry: an example using Cu and Zn isotopes. J Anal Atom Spec. 2004;19:656–65.View ArticleGoogle Scholar
- Toutain JP, Sonke J, Munoz M, Nonell A, Polvé M. Evidence for Zn isotopic variations at Merapi active volcanic system. Chem Geol. 2008;253:74–82.View ArticleGoogle Scholar
- Harmens H, Norris DA, Koerber GR, Buse A, Steinnes E, Rühling A. Temporal trends in the concentration of arsenic, chromium, copper, iron, nickel, vanadium and zinc in mosses across Europe between 1990 and 2000. Atmos Environ. 2007;41:6673–87.View ArticleGoogle Scholar
- Weiss D, Shotyk W, Kramers JD, Gloor M. Sphagnum mosses as archives of recent and past atmospheric lead deposition in Switzerland. Atmos Environ. 1999;33:3751–63.View ArticleGoogle Scholar
- Haack U, Kienholz B, Reimann C, Schneider J, Stumpfl EF. Isotopic composition of lead in moss and soil of the European Arctic. Geochim Cosmochim Acta. 2004;68(12):2613–22.View ArticleGoogle Scholar
- Weiss DJ, Rausch N, Mason TFD, Coles BJ, Wilkinson JJ, Ukonmaanaho L, et al. Atmospheric deposition and isotope biogeochemistry of zinc in ombrotrophic peat. Geochim Cosmochim Acta. 2007;71:3498–517.View ArticleGoogle Scholar
- Novák M, Buzek F, Adamova M. Vertical trends in δ13C, δ15N, and δ34S ratios in bulk Sphagnum peat. Soil Biol Biochem. 1999;31:1343–6.View ArticleGoogle Scholar
- Gélabert A, Pokrovsky OS, Viers J, Schott J, Boudou A, Feurtet-Mazel A. Interaction between zinc and freshwater and marine diatom species: Surface complexation and Zn isotope fractionation. Geochim Cosmochim Acta. 2006;70:839–57.View ArticleGoogle Scholar
- Jouvin D, Louvat P, Juillot F, Maréchal C, Benedetti MF. Zinc isotopic fractionation : why organic matters ? Environ Sci Technol. 2009;43:5747–54.View ArticleGoogle Scholar
- Kafantaris FC, Borrok DM. Zinc isotope fractionation during surface adsorption and intracelular incorporation by bacteria. Chem Geol. 2014;336:42–51.View ArticleGoogle Scholar
- Berg B, Wessen B. Changes in organic-chemical components and ingrowth of fungal mucelium in decomposing birch leaf litter as compared to pine needles. Pedobiologia. 1984;26:285–98.Google Scholar
- Klavins M, Purmalis O. Properties and structure of raised bog peat humic acids. J Mol Struc. 2013;1050:103–13.View ArticleGoogle Scholar
- Colangelo EP, Guerinot MM. Put the metal to the petal: metal uptake and transport throughout the plant. Current Opinion Plant Biol. 2006;9:322–30.View ArticleGoogle Scholar
- Fujii T, Albarede F. Ab initio calculation of the Zn isotope effect in phosphates, citrates, and malates and applications to plants and soil. PLoS ONE. 2012;7:e30726.View ArticleGoogle Scholar
- Kramer U, Clemens S. Functions and homeostasis of zinc, copper, and nickel in plants. In: Tamas MJ, Martinoia E, editors. Molecular Biology of Metal and Homeostasis and Detoxification. Berlin Heidelberg: Springer Verlag; 2005. p. 215–72.View ArticleGoogle Scholar
- Opfergel S, Cardinal D, Henriet C, André L, Delvaux B. Silicon isotope fractionation between plant parts in banana: in situ vis in vitro. J Geochem Expl. 2006;88:224–7.View ArticleGoogle Scholar
- Cenki-Tok B, Chabaux F, Lemarchand D, Schmidt AD, Pierret MC, Viville D, et al. The impact of water-rock interaction and vegetation on calcium isotope fractionation in soil- and stream waters of a small, forested catchment (the Strengbach case). Geochim Cosmochim Acta. 2006;73:2215–28.View ArticleGoogle Scholar
- Engström E, Rodushkin I, Öhlander B, Ingri J, Baxter DC. Silicon isotopic composition of boreal forest vegetation in Northern Sweden. Chem Geology. 2008;257:247–56.View ArticleGoogle Scholar
- Ding TP, Zhou JX, Wan DF, Chen ZY, Wang CY, Zhang F. Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon. Geochim Cosmochim Acta. 2009;72:1381–95.View ArticleGoogle Scholar
- Bolou-Bi EB, Poszwa A, Leyval C, Vigier N. Experimental determination of magnesium isotope fractionation during higher plant growth. Geochim Cosmochim Acta. 2010;74:2523–37.View ArticleGoogle Scholar
- Guelke M, von Blanckenburg F, Schoenberg R, Staubwasser M, Stuetzel H. Determining the stable Fe isotope signature of plant-available iron in soils. Chem Geol. 2010;277:269–80.View ArticleGoogle Scholar
- Weinstein C, Moynier F, Wang K, Paniello R, Foriel J, Catalano J, et al. Isotopic fractionation of Cu in plants. Chem Geol. 2011;286:266–71.Google Scholar
- Kajimoto T, Matsuura Y, Osawa A, Prokushkin AS, Sofronov MA, Abaimov AP. Root system development of Larix gmelinii trees affected by micro-scale conditions of permafrost soils in central Siberia. Plant Soil. 2003;255:281–92.View ArticleGoogle Scholar
- Kajimoto T, Osawa A, Matsuura Y, Abaimov AP, Zyryanova OA, Kondo K, et al. Individual-based measurement and analysis of root system development: case studies for Larix gmelinii trees growing on the permafrost region in Siberia. J Forestry Res. 2007;12:103–12.View ArticleGoogle Scholar
- Brown J, Ferrians OJ, Heginbottom JA, Melnikov ES. Circum-arctic map of permafrost and ground ice conditions. Boulder, CO, USA: National Snow and Ice Data Center/World Data Center for Glaciology. Digital media; 2002.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.