Below-ground residue of agroforestry trees is an important N source for associated crops. Several studies have shown that its isotopic signature (δ15N) may change after tree pruning, which makes it difficult to study below-ground N inputs from pruned trees by isotopic techniques. We studied how temporal variation of legume root residue δ15N could be explained by considering differential decomposition kinetics and 15N content of residue fractions. A mathematical model on the isotopic patterns of soil and a N recipient plant during root decomposition was developed and applied for testing assumptions about residue characteristics against two experimental datasets. Observed 15N patterns of the recipient plants could be satisfactorily simulated only when the residue was assumed to consist of at least two fractions with distinct δ15N and decomposition rates depending on their C : N ratio. Assuming δ15N of residue constant over time resulted in substantial underestimates of N derived from low-quality residue (%Ndfr) by the recipient plant when compared with experimental data. Results of this study suggest that residue fractionation can help improve estimation of %Ndfr in isotopic studies, as an alternative or complementary method to assuming or aiming at homogenous isotopic composition of N sources.

1. Introduction

Pruning of trees is a common practice in agroforestry systems. In legume-based systems the main purpose of pruning is to provide nitrogen to soil and the associated crops from green manure and below-ground residue. Timing and intensity of pruning can be varied to adjust the N inputs with the requirements of the crops, and amount and optimal timing of these inputs have been the interest of numerous studies [1, 2]. Although most studies have concentrated on N release from above-ground biomass, management of N release from roots may be much more important for crop nutrition, since up to 50–60% of total plant N of frequently pruned agroforestry trees may occur in roots [3].

Isotopic techniques are commonly applied for studying the fate of N in agroforestry systems and the mechanisms involved in its cycling. They can be especially useful for studying below-ground N cycling processes that are difficult to trace otherwise. The techniques require measurement or a reasonable estimate of the N isotopic compositions of the N sources. However, recent research suggests that estimating N cycling in intercropping systems with isotopic techniques becomes difficult after management interventions which affect root turnover. After shoot harvest or pruning of the N donor plant, isotopic ratio of the N recipient plants peaks rapidly within a few days and then decreases slowly over time. Such patterns were observed in different experimental setups and N donor species: in studies applying both 15N natural abundance and 15N enrichment methods for the legume tree Gliricidia sepium under field conditions [4, 5] and in pot culture [6], and in a study applying the 15N natural abundance method for the herbaceous legume Canavalia ensiformis in pot culture [7]. The similarity of the results from different setups suggests that the unexpected isotopic patterns of the N recipient plants represent a general phenomenon independent of N donor species and the 15N tracing technique applied.

A key reason for difficulties in analysing the fate of N from organic sources with isotopic techniques certainly is the nonuniform isotopic composition of these sources, which as such is well established. Discrimination against 15N in biochemical processes results in variation of 15N natural abundance between plant organs and their nitrogenous compounds, including soluble proteins, amino acids, and nitrate [8, 9]. Numerous 15N enrichment studies have also demonstrated that 15N labelling techniques and 15N allocation within the plant may result in different plant parts becoming very differently enriched with 15N [10, 11].

Isotopic variation between nitrogenous compounds of a plant could be further enhanced by their different decomposition rates, which are well known from residue decomposition research [12, 13] yet not frequently considered in 15N tracer studies. Together these two factors may result in considerable temporal variation in the isotopic patterns of N released from organic inputs (e.g., roots decomposing after pruning), which is subsequently reflected in isotopic composition of the N recipient pools such as soil compartments, soil microbial biomass, or associated N recipient crops. Studies applying modelling approaches suggest that isotopic composition of N released from decomposing roots of a legume tree differs from that of living roots just before pruning, both with 15N natural abundance and 15N enrichment methods [6] and that N uptake from 15N-enriched, surface-applied residues is best explained when decomposition rate constants are estimated separately for labile and stable residue fractions [14]. However, because isotopic heterogeneity is considered difficult to evaluate over time, N uptake or N transfer studies often assume that isotopic composition of an organic N source remains constant over the length of an experiment. Although homogeneity of 15N labelling has improved through methodological development (reviewed in [15]), it remains difficult to achieve over time, especially at compound level (cf. e.g. [10, 16]). Moreover, inconsistent isotopic signatures can impair N cycling estimates also in 15N natural abundance studies where isotopic compositions of N sources cannot be controlled [4, 6]. Alternative or complementary approaches to more reliably quantify N cycles in agroforestry systems are, therefore, needed.

Our objectives in this study were to (i) quantify how temporal variation in the isotopic composition of decomposing, 15N-enriched legume roots affects the estimates of the fate of recycled N and (ii) study how this variation could be accounted for by considering both the isotopic heterogeneity of residue compounds and their differential decomposition kinetics. To analyse the effects of heterogeneity of decomposing roots, we measured the isotopic patterns of a N recipient plant, which served as integrator of the biological processes involved in N cycling in the soil-plant system. A dynamic model on the isotopic composition of soil pools and the N recipient plant during residue decomposition was developed and then applied for testing the study assumptions against two experimental datasets. The objective of the simulations was not to find the actual parameter values for the specific residues and soil used in the experimental studies, but instead explore which overall patterns and general assumptions about residue characteristics were necessary to explain the observed temporal changes in isotopic composition of the N recipient plants. Options for improving experimental design and interpretation of results in 15N tracer studies and achieving more reliable estimates of the amount and fate of recycled N are then discussed.

2. Material and Methods

2.1. Experimental Design

Two separate datasets were used for evaluating the performance of the residue decomposition and fractionation model under different experimental designs. The first dataset was from a pot experiment where a fodder grass Dichanthium aristatum Poir C.E. Hubbard was grown alone in greenhouse, and 15N-enriched fine root and nodule residue of the legume tree Gliricidia sepium (Jacq.) Kunth ex Walp. were applied in the grass pots (residue application experiment, RA in the following). Evolution of the N isotopic signature (15N) of D. aristatum was then studied for 10 weeks. The experiment is described in the following. The second dataset was that described by Sierra et al. [6], where D. aristatum and G. sepium were grown together in pots in greenhouse (full interaction experiment, FI, in the following). Gliricidia sepium trees were labelled with 15N through foliar application and pruned 12 weeks thereafter to induce root turnover. Evolution of 15N of D. aristatum was then studied for 48 weeks. For this study we used the data for the first 24 weeks, after which the N source from the decomposing residue appeared exhausted [6]. Dichanthium aristatum shoots were harvested every 4–8 weeks during the FI experiment.

The RA experiment was conducted at the greenhouse facilities of the Antillean Research Centre of INRA, Guadeloupe, Lesser Antilles (16°12′N, 61°39′, 125 m a.s.l.). Soil and plant material originated from a cut-and-carry fodder production system of G. sepium and D. aristatum. The site was established in 1989 and managed thereafter by frequent tree pruning and grass cutting. The soil on the site is Vertisol with 80% of clay, pH of 7.8, organic carbon content 33.1 g kg−1, and organic N content 3.1 g kg−1 [6]. For detailed descriptions of the field site and the soil, see Daudin and Sierra [4] and Sierra et al. [17].

The experiment consisted of four potted G. sepium trees for providing the decomposing root and nodule residue, and eight pots of D. aristatum grass for studying the uptake of 15N mineralised from residue. Soil for the pot experiment was collected from the topsoil layer of the field site and sieved to <1 cm aggregates while removing plant roots. Mineral N content of the soil was analysed as described by Sierra et al. [17]. Gliricidia sepium trees were established from cuttings in nursery bags filled with the soil. One month after their propagation, the trees were transferred to pots of 14 L. Another two months later, D. aristatum swards were transplanted from the field site to similar pots as for the trees. Both series of pots were fertilised with 2 g of triple superphosphate and 2 g of K2SO4 at the time of planting and were irrigated daily throughout the experiment. Leaching was assumed negligible because soil inorganic N is mainly in the form of [18] and fixed on soil particles in the clayey soil [19]. Daily mean of air temperature in the greenhouse varied between 25.5 and 30.5°C, with a decreasing trend towards the end of the experiment.

The trees were labelled using foliar feeding of 99%  15N-enriched KNO3 four months after their propagation. The label was applied on tree leaves with a small paintbrush at three events in equal amounts with two day intervals, allowing time for absorption of the solution [6]. The total amount of 15N applied per tree was 30 mg. Foliar feeding was used because the study focused on N recycling from root and nodule turnover induced by pruning or shoot harvest of the N donor plant, common management methods in agroforestry and other intercropping systems. Labelling of above-ground plant parts with 15N is the most feasible 15N labelling methods to study N dynamics related to management practices in such systems. Foliar feeding also enabled comparison of the results with those of a larger study where various below-ground N transfer pathways between N donor and recipient plants were studied applying foliar feeding of 15N [10].

Application of fine root residue took place when two months had passed since grass transplanting in the pots, and three weeks since tree 15N-labelling. The grass shoots were first cut to approximately 2 cm height. The aim was to homogenise the initial situation and to minimise the dilution of assimilated 15N within grass biomass. Cutting was also consistent with grass management at the field site, where the grass is customarily cut every 40–50 days. After grass cutting, the trees were harvested for collecting fine roots (<2 mm of diameter) and nodules which were not detached from the roots. Coarse roots were not used because they decompose slowly, and roots recycled after tree pruning are mainly fine roots [20]. Twelve holes were then carefully drilled in the soil of the grass pots, and fresh fine roots and nodules mixed with a small amount of soil were applied in the holes. The residue was not mixed homogeneously within the soil in order to avoid destructing the grass. Each pot received 3.2 ± 0.6 g of residue, which corresponded to approximately half of the fine root and nodule mass of the trees. All remaining residue material was weighed, oven-dried at 70°C for 72 h, and ground to <0.2 mm for isotopic analysis and for determining the relation of residue fresh and dry weight. Coarse roots were analysed separately.

Grass shoots were sampled for 15N analysis immediately before residue application, and on days 28, 49, and 70 following it. Sampling was limited to these four events in order to avoid excessive disturbance of grass growth. At the end of the experiment on day 70, the grass was harvested and compartmented to shoot, stubble, and roots. All plant material was weighed, dried, and ground for isotopic analysis as described above. Characteristics of the labile and stable fractions were determined through model simulations as explained below. Biochemical properties of residue generally fail to predict N and C dynamics during residue decomposition [21], and the labile and stable fractions of residues are currently determined from laboratory incubations using a model to fit experimental data (e.g., [22]).

Sample N contents and isotopic ratios were determined at the Stable Isotope Facility of the University of California-Davis, US, using an element analyser interfaced to an isotope ratio mass spectrometer (Europa Integra CN; Sercon Ltd., Cheshire, UK).

The differences between the isotopic signatures (15N) of shoot subsample and shoot total biomass at the end of the experiment were tested with Student’s -test, in order to estimate the sampling error for the isotopic signatures during the experiment. Correlation between grass shoot 15N and N concentration during the experiment was calculated as described by Hamlett et al. [23], taking into account that values of the variables were obtained as repeated measures of individual plants and were linked over time. Correlations were expressed as Pearson’s correlation coefficients. Values of were interpreted as indicating statistically significant differences. Results were analysed with SAS statistical analysis software, version 9.1 (SAS Institute Inc., Cary, NC, USA).

2.2. Model Structure

A dynamic model was developed for simulating N mineralisation from decomposing 15N-labelled residue, its subsequent uptake by an adjacent plant, and the plant’s isotopic signature (Figure 1). Residue decomposition was modelled according to the STICS residue decomposition model [24]. The STICS model considers a single fraction for the residue, which decomposes with first order kinetics at a rate determined by its C : N ratio. It has reasonably well-explained overall C and N mineralisation from plant residues in tropical environments [25], including for roots of G. sepium [7]. It can also be easily combined with crop or ecosystem models, as it only requires residue C : N ratio as input [26].

The model developed in this study applies a modified decomposition component with regard to STICS, in order to allow simulation of heterogeneity of residue and its isotopic composition, and the implications to soil and recipient plant in the system. Residue is divided into fractions, and 15N, N content, and C : N ratio can be specified individually for each of these. Uptake of N and its partitioning within the recipient plant was modelled according to a box model [18] which has explained well biomass accumulation and 15N in D. aristatum shoots [6, 18]. The final model consisted of seven types of N pools, namely, residue (NR), microbial biomass involved in decomposition (NB), humified organic matter (NH), soil inorganic N originating from residue (NSR), soil native inorganic N (NSN), and plant roots (NRO) and shoot (NS). All pools were divided to 14N and 15N according to their initial 15N values. Isotopic signatures of N flows were defined by 15N of the source pool.

All residue fractions and the associated microbial biomass pools were set to either deplete or contribute to the same soil inorganic N pool, depending on the C : N ratios of each residue fraction [24]. The model enables simulation of differential N uptake from the NSR and NSN pools by the recipient plant, by including a N source factor, Sf: where u, , and are total N uptake rate of soil N, and N uptake rates from NSR and NSN pools, respectively. Value of Sf in the model is fixed over time. If Sf > 0, then N is first absorbed from the NSR pool and only secondarily from NSN, if supply from NSR is insufficient. If Sf < 0, then N is first absorbed from the NSN pool.

Net mineralisation of soil native organic N () was simulated as a function of soil temperature according to the equation proposed for the same soil and site as in the experimental datasets by Sierra et al. [17]: Because NSN represents net N mineralisation from soil organic matter, it takes into account N immobilisation during mineralisation. However, further N immobilisation from NSN may occur when N demand for growth of the microbial biomass decomposing residues is greater than N supplied from NSR (Figure 1) [24]. Therefore, growth of microbial biomass decomposing residues is primarily supported by organic N in residues, then by NSR and then by NSN.

Ambient temperature and soil humidity affect mineralisation of both residue and microbial N pools [24] and soil native organic N. Mineralisation of humified organic matter originating from the residue was considered negligible as N source for the recipient plant, in comparison to the decomposing residue and soil inorganic N (cf. [27]). Approximately 95–98% of inorganic N in the soil used for the experimental studies is in the form of , possibly because some tropical grasses, including D. aristatum, release compounds which reduce the populations or activity of nitrifying bacteria [18, 28]. Negligible volatilisation of NH3 has been observed in a non-N-fertilised Vertisol similar to that used in the experimental studies [29], and frequent irrigation was assumed to further reduce concentration of NH4 in the soil. Therefore, denitrification and volatilisation were excluded from the model.

Daily uptake of soil inorganic N by the recipient plant was modelled as linear over the time span of the experiments and considering that grass shoots were regularly harvested in the FI experiment. Linear N uptake was based on the results of the FI experiment and previous observations for the N recipient grass [6, 30].

2.3. Calculations of N Uptake by the Recipient Plant

The enrichment of a sample with 15N in relation to the atmospheric standard is expressed as where the subscripts sa and atm refer to the sample and the atmospheric 15N atom-% (0.3663%), respectively. Proportion of N derived from residue (%Ndfr) per total N of the recipient plant was calculated as Subscripts P, SN, and SR refer to the recipient plant and to N derived from soil native N and residue, respectively, 0 to the initial values before 15N enrichment, and t to the point in time. Amount of N derived from residue () in the recipient plant was expressed as where denotes plant N content. Proportion of residue N uptaken by the recipient plant (%) was obtained from where NR is the total N content of the residue.

2.4. Model Parameterisation and Simulations

We assumed that the fine root residue consisted of two fractions, labile and stable. The following constraints were imposed for the labile fraction: (i)  the minimum C : N ratio was set to 3.0, and the maximum N content to 30% of total residue N. These are in the range of values measured for water-soluble N in plant roots [22, 31, 32]. (ii)  The maximum 15N was set to 600 for the RA experiment, and to 495 for the FI experiment. With the maximum N content and respective 15N of the labile residue fraction, 15N of the stable fraction would equal to 1 in each dataset.

Input parameters for total N content, C : N ratio, and 15N of the residue in the RA experiment were the values measured for fine roots and nodules. In the FI experiment the input parameters corresponded to the amount of N recycled after pruning, and total C : N ratio and 15N of roots (including nodules; Table 1). Rate constant equations were common for the two datasets (Table 2). Temperature in the simulations was daily mean temperature in the greenhouse during the experimental studies, and soil humidity was set equal to field capacity as the grass pots were irrigated daily.

The model was run for a period of 70 days for the RA experiment and 168 days for the FI experiment. Three simulation steps were performed for each dataset to evaluate the assumptions about residue characteristics. The steps considered (i) one single residue fraction, that is homogenous residue, (ii) two residue fractions with equal isotopic composition, and (iii) two residue fractions with different isotopic composition. The aim was to test the hypotheses that simulating the observed trends of 15N of the N recipient plants requires more than one residue fraction and that the fractions differ in terms of 15N. Each simulation step was run with two options for the N source factor: N uptake proportionally equal to the size of each inorganic N pool (Sf = 0, (1)), and optimising the factor in the simulation. In simulation step 3 for the RA experiment, N content of the labile fraction was fixed to 29.0, the value obtained in all previous simulations, in order to limit the number of optimised parameters.

Agreement with simulations and experimental observations was evaluated using the coefficient of variation of the root mean square error [34]: and values <0.05 were defined to represent satisfactory model performance. Sensitivity of the model outputs was then studied by varying the values of input parameters and analysing the subsequent changes in the model outputs.

After finding the parameter values which allowed satisfactory simulation of the experimental observations, a fourth simulation step was performed to study the impacts of isotopic heterogeneity of decomposing residue on the N uptake estimates by the N recipient plant. Nitrogen uptake from the residue was estimated according to (5), using two options for 15N of N derived from the residue: (i) total 15N of the residue at the time of residue application, which corresponds to the assumption of a homogenous residue over time (uptake denoted %), and (ii) simulated 15N of N derived from residue (%NdfrS), which can vary over time if residue fractions are differently enriched with 15N and decompose at different rates.

The model was built using the Simile software, version 4.9 (Simulistics Ltd., Edinburgh, UK). Optimal fit was searched by using the PEST software (Model-Independent Parameter Estimation and Uncertainty Analysis, version 11.8; Watermark Numerical Computing, Australia). For detailed description of the software and the optimisation method, see the PEST manual [35].

3. Results

3.1. Experimental Data of the RA Experiment

Fine root and nodule residue of G. sepium applied in the grass pots had C : N ratio of 14.5 ± 0.6, N content of 97.1 ± 2.0 mg pot−1, N concentration of 3.1 ± 0.1%, and 15N of 180 ± 15. Shoot 15N of the N recipient grass increased rapidly at the beginning of the experiment from 4.3 to 33.8 ± 2.1 at the end of the first study period and decreased thereafter during the last two periods (Figure 2). The coefficient of variation for shoot 15N varied between 12 and 22% during the experiment. Grass shoot N concentration peaked simultaneously with shoot 15N (Figure 2) and correlated statistically significantly with it (, ; ). The values of 15N of shoot subsamples did not differ statistically significantly from 15N of total shoot biomass at the end of the experiment (23.5 ± 1.8 versus 21.7 ± 1.9, resp.). Value of 15N of coarse roots of G. sepium at the time of root harvest and residue application was 499 ± 38.

3.2. Simulations

Simulations with a single residue fraction resulted in low 15N values for the N recipient plant in both experimental datasets (Sf > 0, Figures 3(a), and 3(b)). Neither the trends nor the ranges of plant shoot 15N corresponded to the observations. When the residues were divided into two fractions with different C : N but equal 15N, a decreasing trend of shoot 15N after an initial peak was obtained for both datasets, as also observed in the experimental studies. However, simulations deviated too much from the observations for satisfactory results (Figures 3(c), and 3(d); Table 3). Best agreement with observations for both datasets was reached when 15N of the two residue fractions were allowed to differ (Sf > 0, Figures 3(e), and 3(f); Table 3). Characteristics of the residue fractions in those simulations are given in Table 3.

When N uptake was considered proportional to the size of each soil inorganic N pool (i.e., Sf = 0), plant shoot 15N remained low in all simulations for both experimental datasets, regardless of the number of residue fractions and their 15N values. When the N source factor was optimised in the simulations, best fit for both datasets was obtained when proportionally more N was absorbed from the inorganic N pool originating from residue, in comparison to soil native inorganic N (Sf > 0; Table 3).

According to the simulations, 10.3% of N in the N recipient plant originated from the residue at the end of the RA experiment (Figure 4(a)). This corresponded to 25 mg of N and 26% of initial residue N. Nitrogen uptake estimates calculated assuming homogenous residue (%NdfrR) differed only slightly from the simulated N uptake (%NdfrS), except for the first week after residue application (Figure 4(a)). In the FI experiment N originating from decomposing residue constituted 8.1% of N of the recipient plant 25 days after tree pruning (Figure 4(b)), corresponding to 12 mg and 3% of initial residue N. Nitrogen originating from residue rapidly diluted in the plant biomass thereafter as a result of frequent shoot harvests. Nitrogen uptake estimates calculated assuming homogenous residue were 40–51% lower than simulated N uptake during the first 10 days, and 36–39% lower during the rest of the experiment, with the difference slowly decreasing with time (Figure 4(b)).

The model showed a short period of net N immobilisation between days 6 and 16 in the RA experiment. At the end of the simulations, all residue N would have decomposed, and 40% and 34% of it remained in soil microbial biomass and humified organic matter, respectively (Figure 5(a)). In the FI experiment the simulations indicated net N immobilisation from day 0 to day 47. Mineralisation of SOM exceeded immobilisation from day 12 onwards. At the end of the experiment all initial residue N would have decomposed, and 17% of it remained in soil microbial biomass and 76% in humified organic matter (Figure 5(b)). Values presented in Figure 5 correspond to total N content of each pool and include N immobilisation from NSN in soil microbial biomass and humified organic matter.

The sensitivity analyses for the RA experiment indicated that model output (15N of the N recipient plant) was sensitive to C : N ratio of the decomposing microbial biomass (RB), especially during the initial phase of residue decomposition, and to N uptake rate of the recipient plant (u), especially towards the end of the experiment (Figures 6(b), and 6(d)). The model output was negligibly affected by changes in the mineralisation rate constant of soil organic matter (; data not shown), decomposition rate constant of residue (k), and N source factor (Sf; Figures 6(a), and 6(c)) within the tested range of values. Relative values of k for the two residue fractions were retained the same in the simulations, that is, 15N of N released from the residue was not altered. The model was very sensitive to changes in Sf values between 0.0 and approximately 0.3 (data not shown) but was little affected by changes at higher Sf values because of the rapid decomposition and exhaustion of N from the labile residue fraction.

4. Discussion

4.1. Residue Characteristics

We observed a rapid initial peak and a slow subsequent decrease in the isotopic composition of the N recipient plant after incorporation of organic residues in the RA experiment, similar to that obtained in other studies applying both 15N natural abundance and 15N enrichment methods [47, 14]. The simulation results show that the isotopic patterns of the N recipient plants may be explained by the combined effects of a rapid initial release of N from a labile residue fraction, and its differential isotopic composition with regard to total residue. It is noteworthy that the temporal pattern of the 15N content of the recipient plant was not affected by the fact that in the RA experiment the tree roots were harvested for residue application 3 weeks after 15N labelling, whereas in the FI experiment tree pruning which induced root decomposition occurred 12 weeks after labelling. Moreover, the fact that similar trends in the isotopic composition of the N recipient plant were observed in studies applying 15N natural abundance [6, 7] indicates that the observed trends cannot be explained by the method of 15N labelling alone.

The results are in line with the observations of Sierra et al. [6] who showed that 15N of N released from decomposing roots may not correspond to that measured from the living roots just before pruning. By fractioning 15N-enriched residue into labile and stable components, Hadas et al. [14] managed to simulate fairly well the isotopic composition of the N recipient plant, except for its simulated initial peak soon after residue application. They measured isotopic composition of the recipient plant first time 47 days after residue application. This may have concealed the role of isotopic heterogeneity of the residue fractions, which became evident in the RA experiment involving earlier measurements. Recently, it was shown that nitrate content and isotopic fractionation during metabolic processes explain variation in 15N of nitrogenous compounds in rapeseed leaves (Brassica napus L.), where nitrate and amino acids were considerably more enriched than soluble proteins [9]. While fractionation is not of relevance in 15N labelling studies, the results demonstrate that N allocation and plant metabolic processes are capable of considerably altering 15N of nitrogenous compounds and that of residue fractions as a consequence.

The fact that satisfactory simulation of the observations for both experimental datasets in this study was reached with low C : N values of the labile fraction of 3.3 suggests that the fraction consisted of water-soluble components, for example amino acids and possibly inorganic N [31, 32]. Simulated proportion of labile N was much smaller in the FI experiment than in the RA experiment because residue characteristics in the FI experiment corresponded to those of total roots instead of fine roots only. Total N concentration of the root residue in the FI experiment also was approximately a third lower than in the RA experiment (data not shown).

Labile residue fraction appeared more enriched than the stable fraction in the RA experiment, but the opposite was true for the FI experiment in which 15N labelling of the trees took place much earlier. The 15N label in plants may initially occur as metabolically active N or be stored predominantly in coarse roots. In contrast, it can be assumed that N bound in structural components during their formation is generally not replaced over their life time. Therefore, labelling would not affect the isotopic composition of these components, if they were formed before 15N labelling, as presumably was the case for majority of the fine roots of G. sepium in the RA experiment. Relative 15N enrichment of root fractions in the RA experiment (coarse roots > labile fraction of fine roots > stable fraction of fine roots) is in line with these assumptions about the fate of applied 15N within the plant. Relative 15N enrichment of the fractions can be assumed to change over time, as metabolically active and stored 15N are converted to structural components during biomass growth. Differential 15N enrichment in the FI experiment may also partly be explained by the fact that stable fraction corresponded to both the structural components of the fine roots, and coarse roots as a whole, including stored 15N. Method of 15N labelling may importantly affect the relative 15N enrichment of residue fractions [15]. An early start of 15N labelling would help to enrich structural components of the plant more evenly, but may in turn disturb the system more through repeated interventions and need of higher 15N inputs.

4.2. Nitrogen Availability from the Residue to the Associated Plant

The observed correlation between shoot N concentration and 15N of the N recipient grass after the incorporation of 15N-labelled residue in the RA experiment indicates that grass obtained proportionally more N from decomposing residue than from soil native N pools. This was supported by the results of the modelling study, where proportionally higher N uptake of N mineralised from residue compared to soil native N was necessary for obtaining a satisfactory agreement with the experimental observations (Figure 3). Residue application appeared, therefore, to enhance grass N nutrition instead of merely substituting soil native N as N source (cf. [36]).

Plant N uptake is affected by the availability of N in the proximity of roots [37]. The physicochemical properties of the very clayey soil in the experimental studies restrict solute flow [38, 39]. Therefore, it can be assumed that spots of high N concentration emerged in soil as a result of root decay, and that the N recipient plant absorbed proportionally more N from these spots. Such spots are also more porous as the surrounding soil, which may have facilitated colonisation by grass roots and subsequent uptake of N originating from residue [40, 41]. Soil in the holes where the residue in the RA experiment was applied was more loosely packed than the elsewhere in the pots, which may have simulated the effect of increased soil porosity. Small diameter of the holes, and the fact that grass roots had already effectively colonised the pots, likely contributed to the rapid capture of N from the residue. In the FI experiment roots of G. sepium and D. aristatum were mixed together in the pots, and N release from decomposing root residue could be assumed to have been spatially more homogenous than in the RA experiment. Intraspecific competition for soil N could explain why according to the simulations grass would have absorbed proportionally more N originating from residue than in the RA experiment. Grasses are generally more effective competitors for soil N than legumes [42]. Although heterogenous distribution of residue influences N uptake from decomposing roots also in natural systems [41], it was probably more pronounced in the experimental datasets because of experimental factors. The differences in N uptake between the two sources do not, however, affect the results about residue fractionation and distinctive 15N enrichment of labile and stable residue fractions discussed above (cf. Figures 3(e), and 3(f)).

According to the simulations, the N recipient plant had obtained 26% of initial residue N at the end of the RA experiment, and 3% of initial residue N during the first 25 days of the FI experiment. These were on the lower side of the values measured in previous studies, which range from 10 to 100% of total residue N taken up by the associated plants [4346]. Variation in the observations may result for example from residue type and method of application, biological or environmental conditions facilitating or restricting decomposition, capacity of the soil to sequester organic inputs, synchrony with the nutrient needs of the recipient plant, methods for estimating N uptake, and the length of the observation period. Nitrogen from below-ground residues appears to decompose initially faster than that of above-ground residues, and more of it is also taken up by subsequent crops [12, 43, 47]. Very little research, however, has focused on recycling of plant below-ground N, especially in tropical agroecosystems.

4.3. Implications of the Findings to N Cycling Estimates

Quantification of N cycling in agroforestry systems using isotope techniques crucially depends on the determination of the isotopic signature of the N sources. We used simulation methods to evaluate how the common assumption of homogenous isotopic composition and decomposition rates of organic N sources affects estimates of N recycling from decomposing residue, when the residue fractions actually differ in both parameters. Accuracy of the estimates depends on the quality of the residue, relative sizes, and 15N enrichment of its fractions. Nitrogen recycling estimates may not be largely affected by isotopic heterogeneity when only relatively high-quality residues (e.g. fine roots of legume trees) are concerned, which decompose entirely over a short time and N originating from the residue fractions become effectively mixed in the system. In contrast, substantial errors in N recycling estimates can be obtained for lower-quality residues with large, differentially 15N-enriched stable fractions, if isotopic composition of released N is assumed constant. Role of the residues in crop N nutrition in such cases may be underestimated as was the case for the FI experiment where the labile residue fraction was less enriched with 15N than the stable fraction. Depending on the timing or method of 15N labelling, N recycling might also be overestimated if labile residue fraction is more enriched with 15N than the stable fraction (cf. the RA experiment).

While homogenous 15N labelling of the N source is commonly assumed in 15N enrichment studies, it is generally difficult to achieve over time. Assimilation of N from soil and its allocation within the plant during growth, and in legumes also N2-fixation, result in dilution of the label at different rates for different organs [16, 48]. Moreover, in 15N natural abundance studies the isotopic composition of the N sources cannot be controlled [9]. Analysis of temporal isotopic variation of the N sources and reasons affecting it could provide an alternative or complementary method aiming at, or simply assuming, homogenous isotopic composition in 15N tracer studies. Measuring 15N and N content of labile and stable residue fractions separately can be a useful first step in evaluating whether isotopic heterogeneity is likely to result in too small or large estimates for N recycling from organic sources in the concerned system. Water-soluble N may be a good approximation for the labile fraction for this purpose (cf. [2, 31]).

If no information of 15N of residue fractions in 15N enrichment studies is available, caution should be exerted in estimating N recycling and uptake from the 15N values of the recipient pools in the short term. This is especially the case after discrete events in time such as pruning or green manure application, or under pronounced seasonality, if there are large inputs of residue from senescing biomass over a short-time period. Previous studies indicate that these considerations are important also for studies using the 15N natural abundance method [4, 6, 7, 9], and the results explained here reveal that the role of residue fractionation in explaining isotopic patterns in 15N natural abundance studies merits research.


The authors thank Dr. Pekka Nygren for his valuable comments on the paper, and Saint-Ange Sophie for his skillful technical assistance during the experiment. They also thank two anonymous reviewers for their valuable comments. Contribution of R. Jalonen was financed by the Academy of Finland (Grants 111796 and 129166). The authors declare no conflict of interests.