Organic vs. Conventional Grassland Management: Do 15N and 13C Isotopic Signatures of Hay and Soil Samples Differ?

Abstract

Distinguishing organic and conventional products is a major issue of food security and authenticity. Previous studies successfully used stable isotopes to separate organic and conventional products, but up to now, this approach was not tested for organic grassland hay and soil. Moreover, isotopic abundances could be a powerful tool to elucidate differences in ecosystem functioning and driving mechanisms of element cycling in organic and conventional management systems. Here, we studied the δ¹⁵N and δ¹³C isotopic composition of soil and hay samples of 21 organic and 34 conventional grasslands in two German regions. We also used Δδ¹⁵N (δ¹⁵N plant - δ¹⁵N soil) to characterize nitrogen dynamics. In order to detect temporal trends, isotopic abundances in organic grasslands were related to the time since certification. Furthermore, discriminant analysis was used to test whether the respective management type can be deduced from observed isotopic abundances.

Isotopic analyses revealed no significant differences in δ¹³C in hay and δ¹⁵N in both soil and hay between management types, but showed that δ¹³C abundances were significantly lower in soil of organic compared to conventional grasslands. Δδ¹⁵N values implied that management types did not substantially differ in nitrogen cycling. Only δ¹³C in soil and hay showed significant negative relationships with the time since certification. Thus, our result suggest that organic grasslands suffered less from drought stress compared to conventional grasslands most likely due to a benefit of higher plant species richness, as previously shown by manipulative biodiversity experiments. Finally, it was possible to correctly classify about two third of the samples according to their management using isotopic abundances in soil and hay. However, as more than half of the organic samples were incorrectly classified, we infer that more research is needed to improve this approach before it can be efficiently used in practice.

About the Authors:

Valentin H. Klaus

* E-mail: [email protected]

Affiliation: University of Münster, Institute of Landscape Ecology, Münster, Germany

Norbert Hölzel

Affiliation: University of Münster, Institute of Landscape Ecology, Münster, Germany

Daniel Prati

Affiliation: University of Bern, Institute of Plant Sciences, Bern, Switzerland

Barbara Schmitt

Affiliation: University of Bern, Institute of Plant Sciences, Bern, Switzerland

Ingo Schöning

Affiliation: Max-Planck-Institute for Biogeochemistry, Jena, Germany

Marion Schrumpf

Affiliation: Max-Planck-Institute for Biogeochemistry, Jena, Germany

Markus Fischer

Affiliation: University of Bern, Institute of Plant Sciences, Bern, Switzerland

Till Kleinebecker

Affiliation: University of Münster, Institute of Landscape Ecology, Münster, Germany

Introduction

Distinguishing organic and conventional products is a major issue of food security and authenticity and much research on method development has been conducted to tackle this issue [1]. Stable isotope analysis was proven to give important insight in ecosystem functioning and was successfully used to detect differences between organic and conventional agriculture [2], [3], [4]. Since Nakano et al. [5] proposed the use of natural abundances of stable isotopes to separate organic and conventional products, several studies tested this approach successfully for fruits, vegetables and other plant products [1], [6], [7], [8], [9], [10] as well as for beef [11], [12] and milk [13], but not for grassland hay or soil samples. Differences among organic and conventional plant products were mostly attributed to differences in δ¹⁵N isotopic signatures of applied fertilizers [14], because organic farming abandons the use of synthetic mineral fertilizers. While such conventional (synthetic) N sources exhibit δ¹⁵N values close to 0‰, organic N sources such as cattle dung or slurry are strongly enriched in δ¹⁵N [15]. Consequently, organic farming products are mostly enriched in δ¹⁵N compared to conventional ones due to the replacement of synthetic N sources by organic fertilizers. However, in nature δ¹⁵N abundances in plants are affected by a multiplicity of factors such as type and degree of mycorrhization, the chemical type of N-compounds taken up or further soil characteristics, which can be barely separated from each other [16], [17].

Similarly, δ¹³C in plant and soil are also of broad ecological interest [18]. In C3 plants, which represent Central European grassland vegetation, δ¹³C abundances in biomass are first of all affected by water availability and drought stress, but show also significant interactions with nutrient availability and fertilization [2], [19]. Additionally, δ¹³C values are related to a different contribution of CO₂ from soil respiration to plant photosynthesis y [6], [20] and thus contain valuable ecological information related to agricultural management. Furthermore, was shown to be related to functional aspects of plant communities [18].

Although grasslands play a central role in the production of organic meat and dairy products [21], and proportions of organic grasslands have increased significantly during the last decade [22], stable isotope analysis was so far not used to distinguish between soils and yield (hay) of organically and conventionally managed grasslands. As organic fertilizers can even in grasslands lead to higher δ¹⁵N values in soil and vegetation [23], this might give the ability to classify organic and conventional plant products using isotopic abundances [7]. Moreover, isotopic abundances are related to important ecosystem processes affecting nutrient cycling and balances [24] and thus bear the potential to elucidate possible differences in ecosystem functioning of organic vs. conventional grasslands, which are otherwise difficult to detect.

Here, we studied 21 organic and 34 conventional grasslands in two German regions and analyzed δ¹³C and δ¹⁵N of soil and hay (plant biomass). Furthermore, Δδ¹⁵N values (δ¹⁵N plant - δ¹⁵N soil) were calculated to estimate differences in nitrogen dynamics. We also assessed the time since organic certification to test for temporal trends. In detail, we analyzed whether (a) differences in isotopic abundances among organic and conventional grasslands exist and whether there are (b) significant trends in isotopic composition with time since conversion to organic management. Additionally (c), we used discriminant analysis to deduce the respective management type from the isotopic composition of hay and/or soil samples.

Methods

Ethics statement

Field work permits were given by the responsible state environmental offices of Baden-Württemberg, Thüringen, and Brandenburg (according to § 72 BbgNatSchG).

Study design

We studied agriculturally used permanent grasslands in two regions in Germany which belong to the Biodiversity Exploratories project [25]: (I) Hainich-Dün in Thuringia in central Germany situated in and around the National Park Hainich and (II) the Biosphere Reserve Schwäbische Alb in Baden-Württemberg in south-west Germany. In grasslands of both regions Cambisols occur, while in the Schwäbsiche Alb Leptosols and in Hainich-Dün Stagnosols and Vertisols can also be found. Grassland types could be categorized as pastures, meadows and mown pastures [25]. To get information on land use for each grassland, farmers and land owners were annually questioned about the amount and type of fertilizer (kg N ha⁻¹) from 2006 to 2010 [26]. We chose organic and conventional grasslands from a randomly selected dataset of 50 plots in each region. Organic management of grasslands abandons pesticides and synthetic fertilizers, restricts livestock density and the use of organic fertilizers from animal husbandry to a maximum of 170 kg N*ha⁻¹*a⁻¹ (European Union, 2008). Accordingly, study plots can be distinguished in two sub-sets: organic plots which are managed according to an official organic farming certificate [27] and uncertified (conventional) plots, where management goes against certification criteria. Please note that all unfertilized but not certified grasslands were excluded from the analysis. Finally, we used 17 conventional plots per study region as well as 17 organic plots at Hanich-Dün and 4 organic plots at the Schwäbische Alb for comparison. Duration of organic management of the grasslands differs from 3 up to 20 years. Although we have no detailed data on the management prior conversion to organic farming, it seems to be likely that at least some of the organic grasslands were already previously managed at a low to medium intensity. The proportion of legumes varied widely among study plots (from 0.0 to 60.5%) but not among organic and conventional management (data not shown). C4 plants are generally no regular component of Central European grassland vegetation.

Field work and chemical analyses

Soil sampling was conducted in early May 2011. On each plot mixed samples of 14 soil cores from 0 to 10 cm depth were collected using a split tube sampler with a diameter of 5 cm. Cores were taken along two 20 m transects at each plot. Roots were removed from the samples in the field and soil samples were air-dried, sieved to <2 mm and ground. For δ13C analyses, soil samples were weighted into tin capsules and treated with sulphurous acid inside the capsules to remove carbonates. Samples were dried again at 70°C before combustion in an oxygen stream using an elemental analyser (NA 1110, CE Instruments, Milan, Italy). Evolved CO2 was analysed using an isotope ratio mass spectrometer (IRMS; Delta C or DELTA+XL, Thermo Finnigan MAT, Bremen, Germany). For biomass sampling, we harvested aboveground community biomass in four quadrates of 0.25 m² from mid-May to mid-June 2011 in both regions simultaneously. Temporary fences ensured that no mowing or grazing took place before yield was sampled. Plant material was dried for 48 h at 80°C and ground to fine powder for lab analyses. Both isotopic abundances of biomass samples and δ¹⁵N of soil samples were determined by mass spectrometry (Finnigan MAT DeltaPlus with Carlo Erba Elementar Analysator with ConFlo II Interface). Isotope ratios are given in per mille (‰), whereby δ¹³C is relative to the international reference standard v-PDB using NBS19 and δ¹⁵N relative to AIR-N₂ [28].

The nearby climate station in Hainich-Dün, located in Schönstedt (193 m A.S.L.), revealed that the precipitation prior to sampling (April and May 2011) was 70% lower than the mean of the same month from 2003–2010 (33.8 mm instead of 116.2 mm) [29]. At the Schwäbische Alb, the next climate station is situated in Münsingen-Rietheim (732 m A.S.L). Compared to the same month in the period 2005–2010, precipitation in April and May 2011 was 20% lower (153.8 mm instead of 191.5 mm) [30].

Data analysis

Multiple analysis of variance (ANOVA) models were used to examine differences among management types (organic vs. conventional) while also accounting for farm or land owner (n = 7), soil type (n = 4), study region (n = 2) and grasslands type (meadow, pasture, mown pasture) including all two-way interactions with management type. Therefore, the lm() and step() functions in R for stepwise reduction of explanatory variables were used. Obtaining model results by using the anova() output assured that the order of variables entering the model had no effect on significance gained for the respective variable, because all other variables were taken into account previously as co-variables. Model assumptions were checked using diagnostic plots function. Similarly, also applying the lm() function analysis of co-variance (ANCOVA) models were calculated to relate the isotopic composition to the time since certification in all organic plots including farm, soil type and study region. Soil texture (proportions of clay, silt and sand) was also incorporated in lm() models but removed for final analysis as they did not add to explained variance. Additionally, we employed the qda() function from the MASS package [31] to perform Quadratic Discriminant Analysis for the classification of organic and conventional grasslands based on ¹⁵N and ¹³C isotopic values in soil and biomass samples. To ensure normal distribution of variables log10 transformation was applied prior to all analyses where necessary. All statistical tests were performed with R [32].

Results

Farmer's questionnaires revealed that organic grasslands received average fertilizer applications of 21.6 (±37.0) kg N a⁻¹, while conventional grasslands received more than the threefold amount: 72.3 (±41.5) kg N a⁻¹. However, both management systems show wide variation in fertilization intensity at the plot level. Organic plots received only organic fertilizers, as prescribed by organic management guidelines, whereas conventional plots received mostly mineral fertilizers (75% of total fertilizer N).

Grasslands of both management types overlap widely in ¹³C and ¹⁵N isotopic abundances in soil and hay (Figure 1). Nevertheless, organic grasslands were characterized by significantly lower δ¹³C abundances in soil (Figure 2a). Meanwhile, δ¹³C in hay, δ¹⁵N in hay and soil and Δδ¹⁵N did not differ between organic and conventional grasslands but only among farms, study regions and/or grassland types (Table 1). In none of the models the soil type was a significant predictor of ¹³C and ¹⁵N isotopic signals. Although this does not mean that isotopic signals are independent of further soil characteristics, it, nevertheless, underlines the comparability of the selected plots. Generally, analyses explained only 6 to 34% of the variance (Table 1).

[Figure omitted. See PDF.]

Figure 1. δ¹⁵N and δ¹³C composition of a) hay and b) soil samples of organic and conventional grasslands.

https://doi.org/10.1371/journal.pone.0078134.g001

[Figure omitted. See PDF.]

Figure 2. δ¹³C abundances in soil of a) organic vs. conventional grasslands and b) organic grasslands in relation to the time since certification (r_s = −0.70; p<0.001).

Letters indicate significant group differences according to ANOVA analyses (for details see Table 1 & 2).

https://doi.org/10.1371/journal.pone.0078134.g002

[Figure omitted. See PDF.]

Table 1. Summary of multiple ANOVA models of isotopic abundances (no interaction with management was significant) (n = 55). “Grassland type”  =  pasture, meadow or mown pasture.^a

https://doi.org/10.1371/journal.pone.0078134.t001

When related to the time since certification, ¹³C in hay and soil showed significant negative relationships, while the explained variance of the models increased considerably compared to previous models (Table 2). This was most significant for δ¹³C in soil (Figure 2b). Neither δ¹⁵N in soil or hay nor Δδ¹⁵N showed significant relationships with the time since certification (Table 2).

[Figure omitted. See PDF.]

Table 2. Summary of ANCOVA models of isotopic abundances testing for relationships with the time since certification (only organic plots, n = 21).

https://doi.org/10.1371/journal.pone.0078134.t002

Due to significant effects of study region on some of the isotopic abundances (Table 3), quadratic discriminant analysis (QDA) was carried out using data which was centered (standardized) according to the regional mean. Using QDA, 60% of the samples could be correctly classified as organic or conventional using δ¹³C and δ¹⁵N isotopic information of soils, while 70% could be correctly classified using the isotopic information of the hay, assumingly due to higher variation in isotopic values of hay compared to soil samples. Combining both, increased correct classifications only slightly up to 73%. In all cases analyses mismatched organic samples to a higher degree than conventional samples and classified significantly more organic samples wrongly as conventional ones (Table 3).

[Figure omitted. See PDF.]

Table 3. Results of quadratic discriminant analysis (QDA) of management types (organic vs. conventional) deduced from regionally standardized δ¹⁵N and δ¹³C isotopic abundances of soil and/or hay samples of grasslands.

https://doi.org/10.1371/journal.pone.0078134.t003

Discussion

Isotopic abundances of ¹⁵N and ¹³C can give insight in ecosystem functioning and were often shown to be a useful tool to separate between organic and conventional products [1], [6]. In the case of organic management in grasslands, this is only partly true for δ¹⁵N and δ ¹³C in soil and hay.

¹⁵N in organic grasslands

Organic grasslands received only organic fertilizers and were characterized by significantly lower fertilization intensity, in line with Klaus et al. [33]. Nevertheless, this difference in fertilization regime did not imprint in the δ¹⁵N signal of soil and hay samples. While δ¹⁵N abundances of conventional synthetic fertilizers vary between −2 and 2‰, organic fertilizers such as cattle dung or slurry are strongly enriched in δ¹⁵N (5 to 35‰) [15]. While experimental studies revealed an effect of organic compared to synthetic fertilization on the δ¹⁵N isotopic composition of grassland soils and vegetation [23], it was likewise shown that both organic but also mineral fertilizers can lead to increasing δ¹⁵N vales in soil and vegetation probably due to increased microbial activity and subsequent losses of ¹⁵N depleted N from the system [17]. Furthermore, conventional management includes both the use of organic and mineral fertilizers and farmers may considerably change respective proportions among years [26]. In our study, the naturally diverse abiotic conditions such as strong variation in soil properties and land-use history among plots might have additionally impeded clear patterns of δ¹⁵N in organic and conventional grasslands. This is supported by findings from Wrage et al. [34] giving poor to missing associations among N balances and δ¹⁵N in soil and vegetation due to spatially heterogeneous conditions and short-term changes in stocking densities in pastures. Moreover, as only a certain proportion of the N stored in biomass and soil originates from fertilization, biologically fixed N (e.g. through legumes) have additionally diluted isotopic signals of fertilizer applications [35].

Missing differences in Δδ¹⁵N between organic and conventional plots suggest that the studied grasslands do not substantially differ in nitrogen cycling [24]. Instead, the respective farm (or land owner) was responsible for most of the variation explained in Δδ¹⁵N indicating that nitrogen cycling is strongly driven by individual decisions and practice at the farm level. Because plant species richness was shown decease Δδ¹⁵N in plant and soil assumedly via nitrogen partitioning [35], missing or poor differences in plant diversity between organic and conventional grasslands [e.g. 33,36,37] could be another reason why Δδ¹⁵N did not respond to organic grassland management. Contrary, in our study organic grasslands had higher mean plant species richness compared to conventional ones (t-test; p<0.01; data not shown), apparently not leading to significant differences in Δδ¹⁵N.

Haas et al. [21] suggest that lower fertilization intensity under organic management is associated with fewer losses of N to neighbouring habitats or to the ground water. However in our case, even after up to 20 years of organic farming ¹⁵N isotopic abundances revealed no significant indication for fewer losses of N from the studied organic grasslands.

¹³C in organic grasslands

We found significant differences in δ¹³C abundances between organic and conventional soils and trends in decreasing δ¹³C abundances with time since certification in both organic soil and hay samples. Missing statistical significances of soil types in the analyses suggested that this is not only caused by differences in water availability among study plots. However, clear differences in δ¹³C in hay might only get apparent in dry years, because δ¹³C isotopic patterns can differ considerably among seasons [38]. Thus, δ¹³C abundances in soils are probably a more reliable indicator for water supply for plants, because they integrate δ¹³C signals over a longer period of time, while δ¹³C in hay strongly depends on the photosynthetic activity and stomatal conductance influenced by current weather conditions [2], [38]. However, it is questionable whether differences in δ¹³C in hay can get more pronounced than in the studied year as prior to sampling the weather was particularly dry.

It was shown that vegetation δ¹³C decrease with plant diversity in grassland experiments due to facilitation and/or complementarity among plant species [39], [40]. Assumedly, we also detected such a pattern of a positive feedback of higher plant diversity under organic management causing less water stress in organic grasslands. However, as mentioned above, organic and conventional grasslands do not always differ in plant diversity and more research is needed to give full evidence for reduced water stress under higher plant diversity in established permanent grasslands. Furthermore, changes in δ¹³C with time since certification cannot be traced back to changes in plant diversity, because diversity was stable over time (data not shown).

There are several other mechanisms, which could have also occurred, but cannot be easily separated from the whole complex of factors influencing isotopic abundances. Georgi et al. [6] found organic biomass (vegetables) to be depleted in δ¹³C. They suggest higher soil respiration rates in organic soils due to higher microbial activity might cause a lowering of the δ¹³C value of the soil CO₂ pool available to plants. As soil density fractions can have different ¹³C signals [41], differences in the distribution of soil carbon among density fractions could also have also led to isotopic patterns.

Other reasons for lower δ¹³C values in organic grasslands could be direct or indirect effects of (former) fertilization. As organic farms use less often corn (a C4 plant) to feed their livestock, conventional slurry is thus enriched in ¹³C and may have led to higher δ¹³C values in soil organic matter in conventional grasslands [42]. As the “conventional” organic matter deceases after time since certification, this effect is likely to have caused decreasing δ¹³C in soil with time since certification. Similarly, it was shown that high N availability can lead to higher δ¹³C in plant biomass mostly due to structural changes in plant tissue further tightening drought conditions [43]. However, differences in δ¹³C values of vegetables observed by Georgi et al. [6] were independent of optimal or reduced N supply. Finally, we can’t rule out that further factors such as the local rainfall distribution also influenced observed δ¹³C abundances.

Classification of organic and conventional samples

Although Rapisarda et al. [6] showed that approximately 90% of fruit samples could be correctly categorized as organic or conventional products using isotopic abundances, classification of grassland samples revealed rather weak results. Only one third of organic and conventional samples were correctly classified and a high proportion especially of organic samples were mismatched. This can be largely explained by missing differences in especially δ¹⁵N abundances. These were significantly higher in organic fruits used in the study by Rapisarda et al. [6]. Thus, we have to conclude that separating organic and conventional soil and hay samples is barely possible using only δ¹⁵N and δ¹³C abundances. This is especially true, if we would also include uncertified but unfertilized grasslands, which occur frequently [33], [44]. However, including further stable isotopes of S, O and H or other indicative substances might offer additional possibilities to improve the classification of organic vs. conventional samples from permanent grasslands [11], [12].

Acknowledgments

We thank Harald Strauβ and Artur Fugman for friendly isotopic analyses and Svenja Agethen and Max Moenikes for help during laboratory work. Willi A. Brand, Heike Geilmann, Jessica Schäfer and Theresa Klötzing are acknowledged for ¹³C analyses of soil samples. Furthermore, our thank goes to the managers of the exploratories, Swen Renner, Sonja Gockel annd Kerstin Wiesner for their work in maintaining the plot and project infrastructure; Simone Pfeiffer and Christiane Fischer giving support through the central office, Michael Owonibi for managing the central data base, and Markus Fischer, Eduard Linsenmair, Dominik Hessenmöller, Jens Nieschulze, Daniel Prati, Ingo Schöning, François Buscot, Ernst-Detlef Schulze, Wolfgang W. Weisser and the late Elisabeth Kalko for their role in setting up the Biodiversity Exploratories project.

Author Contributions

Conceived and designed the experiments: NH DP TK IS MS MF. Performed the experiments: VHK BS IS MS. Analyzed the data: VHK. Contributed reagents/materials/analysis tools: IS TK. Wrote the paper: VHK TK NH DP MF BS MS IS.

Citation: Klaus VH, Hölzel N, Prati D, Schmitt B, Schöning I, Schrumpf M, et al. (2013) Organic vs. Conventional Grassland Management: Do 15N and 13C Isotopic Signatures of Hay and Soil Samples Differ? PLoS ONE 8(10): e78134. https://doi.org/10.1371/journal.pone.0078134

References

1. Camin F, Perini M, Bontempo L, Fabroni S, Faedi W, et al. (2011) Potential isotopic and chemical markers for characterizing organic fruits. Food Chemistry 125: 1072–1082.

2. Adams MA, Grierson PF (2001) Stable isotopes at natural abundance in terrestrial plant ecology and ecophysiology: An update. Plant Biology 3: 299–310.

3. Franke BM, Gremaud G, Hadorn R, Kreuzer M (2005) Geographic origin of meat - elements of an analytical approachto its authentication. European Food Research and Technology 221: 493–503.

4. Schwertl M, Auerswald K, Schäufele R, Schnyder H (2005) Carbon and nitrogen stable isotope composition of cattle hair: ecological fingerprints of production systems? Agriculture, Ecosystems and Environment 109: 153–165.

5. Nakano A, Uehara Y, Yamauchi A (2003) Effect of organic and inorganic fertigation on yields, δ15N values, and δ13C values of tomato (Lycopersicon esculentum Mill. cv. Saturn). Plant and Soil 255: 343–349.

6. Georgi M, Voerkelius S, Rossmann A, Grassmann J, Schnitzler WH (2005) Multielement isotope ratios of vegetables from integrated and organic production. Plant and Soil 275: 93–100.

7. Rapisarda P, Calabretta ML, Romano G, Intrigliolo F (2005) Nitrogen metabolism components as a tool to discriminate between organic and conventional citrus fruits. Journal of Agricultural and Food Chemistry 53: 2664–2669.

8. Bateman AS, Kelly SD, Woolfe M (2007) Nitrogen Isotope Composition of Organically and Conventionally Grown Crops. Journal of Agricultural and Food Chemistry 55: 2664–2670.

9. del Amor FM, Navarro N, Aparicio PM (2008) Isotopic Discrimination as a Tool for Organic Farming Certification in Sweet Pepper. Journal of Envionmental Quality 37: 182–185.

10. Camin F, Moschella A, Miselli F, Parisi B, Versini G, et al. (2007) Evaluation of markers for the traceability of potato tubers grown in an organic versus conventional regime. Journal of the Science of Food and Agriculture 87: 1330–1336.

11. Boner M, Förstel H (2004) Stable isotope variation as a tool to trace the authenticity of beef. Analytical & Bioanalytical Chemistry 378: 301–310.

12. Schmidt O, Quilter JM, Bahar B, Moloney AP, Scrimgeour CM, et al. (2005) Inferring the origin and dietary history of beef from C, N and S stable isotope ratio analysis. Food Chemistry 91: 545–549.

13. Molkentin J (2009) Authentication of Organic Milk Using δ13C and the r-Linolenic Acid Content of Milk Fat. Journal of Agricultural and Food Chemistry 2009 57: 785–790.

14. Bateman AS, Kelly SD, Jickells TD (2005) Nitrogen Isotope Relationships between Crops and Fertilizer:Implications for Using Nitrogen Isotope Analysis as an Indicator of Agricultural Regime. Journal of Agricultural and Food Chemistry 53: 5760–5765.

15. Bateman AS, Kelly SD (2007) Fertilizer nitrogen isotope signatures. Isotopes in Environmental and Health Studies 43: 237–247.

16. Hobbie EA, Högberg P (2012) Nitrogen isotopes link mycorrhizal fungi and plants to nitrogen dynamics. New Phytologist 196: 367–382.

17. Kleinebecker T, et al. (submitted) 15N natural abundances reveal plant diversity effects on nitrogen dynamics in permanent grassland ecosystems. Submitted

18. de Bello F, Buchmann N, Casals P, Leps J, Sebastia MT (2009) Relating plant species and functional diversity to community δ13C in NE Spain pastures. Agriculture, Ecosystems and Environment 131: 303–307.

19. Högberg P, Johannisson C, Hallgren J-E (1993) Studies of 13Cin the foliage reveal interactions between nutrients and water in forest fertilization experiments. Plant and Soil 152: 207–214.

20. Šantrůčková H, Bird MI, Lloyd J (2000) Microbial processes and carbon-isotope fractionation in tropical and temperate grassland soils. Functional Ecology 14: 108–114.

21. Haas G, Wetterich F, Köpke U (2001) Comparing intensive, extensified and organic grassland farming in southern Germany by process life cycle assessment. Agriculture, Ecosystems and Environment 83: 43–53.

22. Schaack D, Iller S, Würtenberger E (2010) AMI-Marktbilanz Öko-Landbau 2010, Bonn: Agrarmarkt Informationsgesellschaft mbH.

23. Watzka M, Buchgraber K, Wanek W (2006) Natural N-15 abundance of plants and soils under different management practices in a montane grassland. Soil Biology & Biochemistry 38: 1564–1576.

24. Kahmen A, Wanek W, Buchmann N (2008) Foliar delta 15N values characterize soil N cycling and reflect nitrate or ammonium preference of plants along a temperate grassland gradient. Oecologia 156: 861–870.

25. Fischer M, Bossdorf O, Gockel S, Hansel F, Hemp A, et al. (2010) Implementing large-scale and long-term functional biodiversity research: The Biodiversity Exploratories. Basic and Applied Ecology 11: 473–485.

26. Bluethgen N, Dormann CF, Prati D, Klaus VH, Kleinebecker T, et al. (2012) A quantitative index of land-use intensity in grasslands: Integrating mowing, grazing and fertilization. Basic and Applied Ecology 13: 207–220.

27. European Union (2008) Commission Regulations (EC) No 889/2008, Brussels.

28. Werner RA, Brand WA (2001) Referencing strategies and techniques in stable isotope ratio analysis. Rapid Communications in Mass Spectrometry 15: 501–519.

29. LUFTGEIST website. Available: http://www.luftgeist.de. Accessed 2013 Jun 3.

30. Wetterstation Riethein Lichse website. Available: http://www.wetterstation-rietheim-lichse.de/html/archiv.html. Accessed 2013 Jun 3.

31. Venables WN, Ripley BD (2002) Modern Applied Statistics with S. Fourth Edition. New York: Springer.

32. R Development Core Team (2011) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna.

33. Klaus VH, Kleinebecker T, Prati D, Fischer M, Alt F, et al. (2013) Does organic grassland farming benefit plant and arthropod diversity at the expense of yield and soil fertility? Agriculture, Ecosystems and Environment 177: 1–9.

34. Wrage N, Küchenmeister F, Isselstein J (2011) Isotopic composition of soil, plant or cattle hair no suitable indicator of nitrogen balances in permanent pasture. Nutrient Cycling in Agroecosystems 90: 189–199.

35. Gubsch M, Roscher C, Gleixner G, Habekost M, Lipowsky A, et al. (2011) Foliar and soil delta 15N values reveal increased nitrogen partitioning among species in diverse grassland communities. Plant Cell and Environment 34: 895–908.

36. Hole DG, Perkins AJ, Wilson JD, Alexander IH, Grice PV, et al. (2005) Does organic farming benefit biodiversity? Biological Conservation 122: 113–130.

37. Batáry P, Báldi A, Sárospataki M, Kohler F, Verhulst J, et al. (2010) Effect of conservation management on bees and insect-pollinated grassland plant communities in three European countries. Agriculture, Ecosystems and Environment 136: 35–39.

38. Neilson R, Hamilton D, Wishart J, Marriott CA, Boag B, et al. (1998) Stable isotope natural abundances of soil, plants and soil invertebrates in an upland pasture. Soil Biology and Biochemistry 30: 1773–1782.

39. Caldeira MC, Ryel RJ, Lawton JH, Pereira JS (2001) Mechanisms of positive biodiversity-production relationships: insights provided by δ13C analysis in experimental Mediterranean grassland plots. Ecology Letters 4: 439–443.

40. Jumpponen A, Mulder CPH, Huss-Danell K, Högberg P (2005) Winners and losers in herbaceous plant communities: insights from foliar carbon isotope composition in monocultures and mixtures. Journal of Ecology 93: 1136–1147.

41. Baisden WT, Amundson R, Cook AC, Brenner DL (2002) Turnover and storage of C and N in five density fractions from California annual grassland surface soils. Global Biogeochemical Cycles 16: 64–71.

42. Bol R, Moering J, Preedy N, Glaser B (2004) Short-term sequestration of slurry-derived carbon into particle size fractions of a temperate grassland soil. Isotopes in Environmental and Health Studies 40: 81–87.

43. Shangguan ZP, Shao MA, Dyckmans J (2000) Nitrogen nutrition and water stress effects on leaf photosynthetic gas exchange and water use efficiency in winter wheat. Environmental and Experimental Botany 44: 141–149.

44. Socher S, Prati D, Müller J, Klaus VH, Hölzel N, et al. (2012) Direct and productivity-mediated indirect effects of fertilization, mowing and grazing intensities on grassland plant species richness. Journal of Ecology 100: 1391–1399.