Biochemical Society Focused Meeting

Advances in 15N-tracing experiments: new labelling and data analysis approaches

Tobias Rütting, Dries Huygens, Jeroen Staelens, Christoph Müller, Pascal Boeckx


To obtain an in-depth understanding of soil nitrogen dynamics, it is necessary to quantify a variety of simultaneously occurring gross nitrogen transformation processes. In order to do so, most studies apply 15N in a disturbed soil–microbial–root system and quantify gross rates based on the principles of 15N isotope dilution. However, this approach has several shortcomings. First, studying disturbed soil provides only limited information on in situ soil nitrogen dynamics. Secondly, the analytical data analysis allows the quantification of total production and consumption rates of the labelled pool, but does not provide information on process-specific transformation rates. Combining in situ 15N isotope labelling over 1–2 weeks with numerical data analysis allows determining process-specific gross nitrogen transformations in undisturbed soils under field conditions in the presence of live roots and their associated microbial communities. This has the potential to increase our understanding of nitrogen dynamics in the soil environment.

  • 15N labelling
  • 15N tracing
  • nitrogen transformation
  • soil nitrogen cycle

In situ15N labelling

Measurements of gross nitrogen dynamics are commonly made by 15N isotope dilution studies. The principle of the isotope dilution technique is based on the dilution of the product pool that has been labelled with 15N, such as ammonium (NH4+) or nitrate (NO3) [1]. Three experimental set-ups with various degrees of rhizosphere distortion are often used to assess gross nitrogen transformations: (i) laboratory incubations of fully disturbed soils, i.e. after mixing, drying, sieving, storing and/or root-picking of the soil (e.g. [2]), (ii) laboratory incubations of intact field-collected soil cores, eventually after a short cold storage period (e.g. [3]) and (iii) in situ incubations of introduced soil cores (e.g. [4]). All three set-ups may provide valuable information since they can be used to derive potential gross nitrogen transformations and to assess the influence of external drivers on soil nitrogen dynamics. Nevertheless, the study set-ups are known to alter factors that influence soil nitrogen transformations such as the chemistry of soil organic matter [5], the concentration and mobility of nitrogen pools [6], the microbial community structure and the root biomass [7,8]. Considering the differential degree of alteration of these manyfold factors, substantial discrepancies in gross nitrogen transformation fluxes have been observed as a function of the used set-up. For instance, Booth et al. [9] performed gross nitrogen transformation rates by in situ incubations using introduced soil cores in which the soil was either mixed or left intact before starting the 15N isotopic dilution experiment. It was indicated that soil mixing promoted gross nitrogen mineralization and NH4+ consumption [9]. Arnold et al. [10] showed that cold storage and laboratory incubation of intact soil cores resulted in a lower NH4+ and higher NO3 turnover compared to in situ incubated soil cores.

While concerns related to soil organic matter and nitrogen pools are largely surmounted when using in situ incubations of introduced soil cores [1,11], most current experimental set-ups fail to keep the soil–microbial–root system intact [12]. The general practice of inserting cylinders before adding 15N to the soil excises roots and disturbs the exchange of resources between roots and their associated microbial communities. Altering these interactions leads to unrepresentative estimates of nitrogen availability in plant communities and an incomplete understanding of the environmental factors that control plant-available nitrogen [12]. Plant roots secrete a variety of low- and high-molecular-mass compounds into the soil as exudates [13] that influence the turnover, composition and activity of the soil microbial community [1417]. Additionally, mycorrhizal fungi associated with plant roots alter soil nitrogen availability by competing for organic and inorganic nitrogen sources in the soil matrix [18]. Consequently, laboratory experiments that manipulated root densities and functioning have indicated the overarching influence of live roots on gross nitrogen transformations [17,1921]. Isotope dilution methods applied to the intact soil–microbial–root system have rarely been reported. One example is the study by Templer et al. [22] which investigated nitrogen dynamics using in-growth cores in a tropical montane forest soil. However, owing to the installation of the in-growth cores 11 months before the 15N isotope dilution study, fine root biomass in the cores was only 59% of the background level in undisturbed soil [22].

We developed a new in situ 15N soil-labelling method, called ‘virtual soil core’ injection, which is based on a technique to examine the dynamics of soil amino acids (e.g. [23,24]), but that to our knowledge has only recently been used to study gross nitrogen transformations (J. Staelens, T. Rütting, D. Huygens, A. De Schrijver, C. Müller, K. Verheyen and P. Boeckx, unpublished work, and T. Rütting, P. Boeckx, J. Staelens and L. Klemedtsson, unpublished work). In this approach, gross nitrogen dynamics are studied under field conditions in undisturbed soils with live roots and their associated microbial communities. Within an experimental site, several representative plots (approx. 1 m×1.5 m) are selected as replicates. For each 15N labelling treatment and time step, a 10 cm×10 cm location is prepared in each plot (Figure 1). The four corners of the location are then marked out by small sticks so that templates for nitrogen addition and soil sampling can be used later on. To study mineral forest soils, it is helpful to separate the litter and mineral layer by placing a nylon mesh on top of the mineral layer. As such, the litter layer can easily be removed during the following injection and sampling steps. To ensure that the microclimatic conditions of the labelled soils remain as intact as possible, any natural soil cover (e.g. straw in arable fields or litter layers in forests) should only be removed for 15N injection, but replaced during the incubation time until soil sampling.

Figure 1 Example layout of 15N labelling plot for two 15N treatments and six time steps

The square plots indicate the locations for injection of 15N-enriched solution at t0 containing either 15NH4NO3 or NH415NO3, and the circular areas within the squares indicate the ‘virtual soil cores’ that will be sampled at different times (t=time of soil sampling). Example distances between plots (within and between rows) and times for extraction are given, but they will depend on the actual field conditions.

For injection of 15N labelling solutions, where one of the moieties is enriched in 15N (e.g. NH4+, NO2 and NO3), a 10 cm×10 cm template is slid over the sticks. The template contains a sufficient number of holes in order to inject the 15N enriched solutions in the studied soil layer. Hart et al. [1] recommended at least six injection locations within a soil core of 4 cm diameter. However, to minimize edge effects, it is advisable to inject the 15N solution into a larger area than the final sample core size (T. Rütting, P. Boeckx, J. Staelens and L. Klemedtsson, unpublished work). Similarly to isotope dilution studies, the 15N solutions should be added as uniformly as possible both horizontally and vertically [27]. Soil sampling is carried out at defined time steps after 15N injection (e.g. 15 min, 1, 2, 5, 9 and 12 days) by pushing a PVC sampling tube into the soil layer through the fitting central hole of a second template. The templates with the marking sticks are used to ensure the sampling of the 15N labelled soil. The sampled soil cores are then processed as quickly as possible (weighing, sieving and extraction) and the size and 15N content of nitrogen pools are determined (e.g. NH4+, NO3 and organic nitrogen). To quantify potential nitrogen losses during the field experiment in the open ‘virtual’ soil cores, it is recommended to measure the size and 15N content of total soil nitrogen at the different time points.

Numerical data analysis

In 15N labelling experiment, gross soil nitrogen dynamics are typically determined by analytically solved equations with data from 15N dilution experiments. The 15N isotope dilution model of Kirkham and Bartholomew [28] is applied to determine the gross rates between two nitrogen pools. For this analysis, it is assumed that (i) no isotopic fractionation occurs during microbial transformations of soil nitrogen, (ii) no remineralization of added 15N takes place and (iii) nitrogen transformation rates are usually assumed to be constant [1], although first-order models have also be developed (e.g. [29]). Thus analytical equations can most realistically only be applied for short-term incubations when remineralization is negligible [3032]. Moreover, analytical solutions do actually quantify only the total gross production and consumption of the labelled pool and do not provide information on gross rates for specific processes [33].

To overcome the limitations of analytical solutions, it was recommended to use a numerical data analysis via so-called 15N-tracing models [32,34,35]. Barraclough and Puri [36] used such an approach to separate NO3 production, quantified by an analytical model, into heterotrophic and autotrophic nitrification pathways. Generally, numerical approaches combine the principles of the 15N dilution, where a product pool is labelled, with 15N-tracing techniques, where a substrate pool is labelled and the movement of the isotope is followed in various product pools over a longer time period [1]. These 15N-tracing models consist of at least three nitrogen pools (NH4+, NO3 and organic nitrogen) and various nitrogen transformations that simultaneously transfer nitrogen between the pools. Process-specific gross rates for the nitrogen transformations can be quantified via non-linear optimization routines [35,37]. However, 15N-tracing models often face the problem of over-parameterization [35], for which reason only a low number of nitrogen processes and simple kinetic settings (e.g. no Michaelis–Menten kinetics) can be considered [35,37,38]. To apply more complex, and hence probably more realistic, 15N-tracing models, a robust optimization technique, not limited by the degree of freedom, was developed based on Monte Carlo sampling techniques (Figure 2, [37]). It was shown that this 15N analysis tool was able to reliably quantify a range of gross nitrogen transformation rates, with the additional benefit that a suitable kinetic setting could be implemented (i.e. consideration of Michaelis–Menten kinetics if appropriate).

Figure 2 Example of a conceptual 15N-tracing model to analyse gross nitrogen transformation rates

The considered pools in this example are ammonium (NH4+), nitrate (NO3) as well as two organic nitrogen pools (labile Nlab and more recalcitrant Nrec). Eight nitrogen transformations are considered: mineralization of labile and recalcitrant organic nitrogen (MNlab and MNrec), immobilization of NH4+ (INH4→Nlab and INH4→Nrec) and of NO3 (INO3), oxidation of NH4+ (ONH4) and organic nitrogen (ONrec) to NO3 as well as dissimilatory NO3 reduction to NH4+ (DNO3). Modified from Soil Biology and Biochemistry, vol. 39, C. Müller, T. Rütting, J. Kattge, R.J., Laughlin and R.J. Stevens, Estimation of parameters in complex 15N tracing models via Monte Carlo sampling, pp. 715–726, © 2007, with permission from Elsevier.

Applying this novel 15N methodology helped us to deepen our understanding of specific pathways within the soil nitrogen cycle. Experiments in grassland soils have revealed more insight into gross nitrogen dynamics in conjunction with nitrite (NO2) dynamics, a precursor of nitrous oxide, indicating the importance of oxidation of organic nitrogen for the soil NO2 pool [39]. For a pristine forest soil, Rütting et al. [40] showed that dissimilatory NO3 reduction to NH4+ (DNRA) was a more important NO3 consumption process than denitrification. Applying the virtual soil core approach indicated the importance of heterotrophic nitrification (oxidation of organic nitrogen to NO3) in forest soils (J. Staelens, T. Rütting, D. Huygens, A. De Schrijver, C. Müller, K. Verheyen and P. Boeckx, unpublished work). The use of the virtual soil core set-up might have been crucial in this study, since soil disturbing processes such as sieving or delimiting soil cores have been observed to disrupt the functioning of heterotrophic nitrifiers such as soil fungi [41,42].

Advantages and drawbacks

The first advantage of the described 15N application and numerical analysis methods is that the obtained gross nitrogen dynamics are representative of the actual field conditions because the soil–microbial–root system is only marginally disturbed. Hence, soil properties are not changed, and heat, nutrients and gas can exchange according to actual field conditions. Nevertheless, short-term 15N labelling studies face the problem that nitrogen transformations can be altered as a result of the 15N addition itself, as the mineral soil nitrogen and the water contents are increased. It is therefore recommended to add only small amounts of highly enriched 15N solutions [1]. The potential stimulation of the production and consumption rates by adding inorganic nitrogen will be greatest at the start of an experiment. As the rates are averaged over time, this drawback is largely compensated by the fact that the 15N-tracing approach with numerical data analysis allows longer study periods (1–2 weeks) than in common isotope dilution experiments with analytical data analysis (usually 1–2 days). The second advantage is that a numerical 15N-tracing model provides information on many simultaneously occurring gross nitrogen transformation rates. This clearly increases our understanding of soil nitrogen dynamics because the contribution of specific nitrogen transformation pathways to the total production and consumption of organic and inorganic nitrogen pools can be quantified.

The numerical 15N-tracing analysis requires, however, more time steps and thus more chemical analyses of soil samples in comparison with the 15N isotope dilutions. This has earlier been discussed as an disadvantage of the tracing techniques [34]. Additionally, the inherent soil heterogeneity may result in a high variation in gross nitrogen cycling fluxes between replicate samples (e.g. [43]), for which a high number of spatial replicates in the field should be taken.

The nitrogen transformation rates obtained depend, among other things, on the actual soil temperature, moisture and microclimatic conditions during the study. Precipitation during the experimental period should be controlled, since changes in soil moisture will affect microbial activity and nitrogen process rates. The dependency on actual soil conditions implies that the comparison of nitrogen dynamics between sites or treatments should occur as synchronized as possible. This requires appropriate coordination during field campaigns.

Similarly to classical 15N experiments, there could be an incomplete recovery of 15N in the sampled soil due to various processes, including gaseous nitrogen loss through ammonia volatilization or denitrification, uptake and transport of inorganic nitrogen by roots and mycorrhizae, leaching of dissolved nitrogen below the sampled soil, and diffusion of injected 15N to non-labelled soil around the sampled virtual cores. In order to quantify nitrogen loss pathways in the tracing model, the 15N content of the bulk soil should be analysed (J. Staelens, T. Rütting, D. Huygens, A. De Schrijver, C. Müller, K. Verheyen and P. Boeckx, unpublished work). The inflow of unlabelled mineral nitrogen via passive or active transport processes may lead to an overestimation of the gross nitrogen production fluxes. This diffusion can be accounted for by injecting 15N into an edge zone around the sampled soil (T. Rütting, P. Boeckx, J. Staelens and L. Klemedtsson, unpublished work).


It is concluded that the presented methodology has the potential to increase our understanding in gross nitrogen cycling studies that focus on the soil–microbial–root system. The incorporation of these methodologies into ecosystem biogeochemistry will allow discerning the importance of rhizosphere nitrogen transformation pathways other than the classical mineralization–immobilization turnover.


T.R. was supported by the NitroEurope IP under the European Commission Sixth Framework Programme [contract number 017841]. D.H. and J.S. were funded as postdoctoral fellows by the Research Foundation-Flanders (FWO). A Short Visit Grant of the European Science Foundation Nitrogen in Europe (NinE) research network was granted to J.S.


  • Enzymology and Ecology of the Nitrogen Cycle: A Biochemical Society Focused Meeting held at University of Birmingham, U.K., 15–17 September 2010. Organized and Edited by Jeff Cole (University of Birmingham, U.K.), Rosa María Martínez-Espinosa (University of Alicante, Spain), David Richardson (University of East Anglia, Norwich, U.K.) and Nick Watmough (University of East Anglia, Norwich, U.K.).


View Abstract