Ecological Archives E095-232-D1
Laura Kirwan, John Connolly, Caroline Brophy, Ole Baadshaug, Gilles Belanger, Alistair Black, Tim Carnus, Rosemary Collins, Jure Čop, Ignacio Delgado, Alex De Vliegher, Anjo Elgersma, Bodil Frankow-Lindberg, Piotr Golinski, Philippe Grieu, Anne-Maj Gustavsson, Įslaug Helgadóttir, Mats Höglind, Olivier Huguenin-Elie, Marit Jųrgensen, ˇydre Kad˛iulienė, Tor Lunnan, Andreas Lüscher, Päivi Kurki, Claudio Porqueddu, M.-Teresa Sebastia, Ulrich Thumm, David Walmsley, and John Finn.. 2014. The Agrodiversity Experiment: three years of data from a multisite study in intensively managed grasslands. Ecology 95:2680. http://dx.doi.org/10.1890/14-0170.1
Metadata
Class I. Data set descriptors
A. Data set identity: The Agrodiversity Experiment: three years of data from a multisite plant diversity experiment in intensively managed grasslands.
B. Data set identification code: site_info.csv biomass.csv forage_quality.csv climate.csv soils.csv
C. Data set description
Abstract: Intensively managed grasslands are globally prominent ecosystems. We investigated whether experimental increases in plant diversity in intensively managed grassland communities can increase their resource use efficiency. This work consisted of a coordinated, continental-scale 33-site experiment. The core design was 30 plots, representing 15 grassland communities at two seeding densities. The 15 communities were comprised of four monocultures (two grasses and two legumes) and 11 four-species mixtures that varied in the relative abundance of the four species at sowing. There were 1028 plots in the core experiment, with another 572 plots sown for additional treatments. Sites agreed a protocol and employed the same experimental methods with certain plot management factors, such as seeding rates and number of cuts, determined by local practice. The four species used at a site depended on geographical location, but the species were chosen according to four functional traits: a fast-establishing grass, a slow-establishing persistent grass, a fast-establishing legume, and a slow-establishing persistent legume. As the objective was to maximize yield for intensive grassland production, the species chosen were all high-yielding agronomic species. The data set contains species-specific biomass measurements (yield per species and of weeds) for all harvests for up to four years at 33 sites. Samples of harvested vegetation were also analyzed for forage quality at 26 sites. Analyses showed that the yield of the mixtures exceeded that of the average monoculture in >97% of comparisons. Mixture biomass also exceeded that of the best monoculture (transgressive overyielding) at about 60% of sites. There was also a positive relationship between the diversity of the communities and aboveground biomass that was consistent across sites and persisted for three years. Weed invasion in mixtures was very much less than that in monocultures.
These data should be of interest to ecologists studying relationships between diversity and ecosystem function and to agronomists interested in sustainable intensification. The large spatial scale of the sites provides opportunity for analyses across spatial (and temporal) scales. The database can also complement existing databases and meta-analyses on biodiversity–ecosystem function relationships in natural communities by focusing on those same relationships within intensively managed agricultural grasslands.
D. Key words: agricultural grasslands; biodiversity; ecosystem function; forage quality; mixtures; monocultures; overyielding; plant community; species biomass; yield.
E. Authors:
Laura Kirwan,1,25 John Connolly,2 Caroline Brophy,3 Ole Baadshaug,4 Gilles Belanger,5 Alistair Black,6 Tim Carnus,2,7 Rosemary Collins,8 Jure Čop,9 Ignacio Delgado,10 Alex De Vliegher,11 Anjo Elgersma,12 Bodil Frankow-Lindberg,13 Piotr Golinski,14 Philippe Grieu,15 Anne-Maj Gustavsson,16 Áslaug Helgadóttir,17 Mats Höglind,18 Olivier Huguenin-Elie,19 Marit Jørgensen,18 Žydrė Kadžiulienė,20 Tor Lunnan,18 Andreas Lüscher,19 PÄivi Kurki,21 Claudio Porqueddu,22 M.-Teresa Sebastia,23 Ulrich Thumm,24 David Walmsley,1,7 and John Finn7
1 Waterford Institute of Technology, Cork Road, Waterford, Ireland
2 School of Mathematical Sciences, University College Dublin, Dublin 4, Ireland
3 Department of Mathematics and Statistics, National University of Ireland Maynooth, Maynooth, County Kildare, Ireland
4 Department of Plant Sciences, Norwegian University of Life Sciences, P.O. Box 5003, N-1432 Ås, Norway
5 Agriculture and Agri-Food Canada, 2560, Hochelaga Blvd, Québec, Québec G1V 2J3, Canada
6 Faculty of Agriculture and Life Sciences, PO Box 85084, Lincoln University, Lincoln 7647, Christchurch, New Zealand
7 Teagasc, Environment Research Centre, Johnstown Castle, Wexford, Ireland
8 Aberystwyth University, IBERS, Plas Gogerddan, Aberystwyth SY23 3EE, Wales, UK
9 Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000 Ljubljana, Slovenia
10 CITA-DGA, Av. Montañana 930, 50059 Zaragoza, Spain
11 ILVO, Plant Sciences Unit, Crop Husbandry and Environment, B. Van Gansberghelaan 109, 9820 Merelbeke, Belgium
12 Plant Sciences Group, Wageningen University, PO Box 16, 6700 AA Wageningen, The Netherlands
13 Swedish University of Agricultural Sciences, Department of Crop Production Ecology, Box 7043, SE-750 07 Uppsala, Sweden
14 Department of Grassland and Natural Landscape Sciences, Poznan University of Life Sciences, Wojska Polskiego 38/42, 60-627 Poznan, Poland
15 University of Toulouse; UMR AGIR, INP-ENSAT, F-31326 Castanet Tolosan, France
16 Swedish University of Agricultural Sciences, Department of Agricultural Research for Northern Sweden, SE-901 83 Umeå, Sweden
17 Agricultural University of Iceland, Árleyni 22, IS-112 Reykjavík, Iceland
18 Bioforsk, Norwegian Institute for Agricultural and Environmental Research, Frederik A. Dahls vei 20, N-1430 ÅS
19 Agroscope Institute for Sustainability Sciences ISS, Reckenholzstrasse 191, 8046 Zurich, Switzerland
20 Institute of Agriculture, Lithuanian Research Centre for Agriculture and Forestry, Akademija, LT-58344, Kedainiai, Lithuania
21 MTT Agrifood Research Finland, Plant Production Research Lönnrotinkatu 5, FI - 50100 Mikkeli, Finland
22 CNR-ISPAAM, Traversa la Crucca 3, località Baldinca, 07100 Sassari, Italy
23 Centre Tecnologic Forestal de Catalunya, Ctra Sant Llorenç km 2, 25280 Solsona, Spain
24 University of Hohenheim, Department of Crop Science, 70593 Stuttgart, Germany
25 Corresponding author. E-mail: [email protected]
Class II. Research origin descriptors
The Agrodiversity experiment was a coordinated continental-scale field experiment. The coordination of the network was supported by EU COST Action 852: Quality Legume-Based Forage Systems for Contrasting Environments. Each of the 33 participating sites established at least 30 plots using the same experiment design and managed the plots according to an agreed protocol. The design of the experiment consisted of a common set of 30 plots across all sites, with additional optional plots for applying treatments. Aboveground biomass was measured at the common set of 30 plots for all harvests at all sites. Forage quality measurements were undertaken at 26 sites. Additional treatments were voluntarily undertaken by various sites. The treatments applied varied across sites and included increased genetic diversity in species, different levels of nitrogen fertilizer and different levels of cutting intensity.
The network did not have a central funding source to cover the costs of running the experiment, thus participating sites each secured funding to cover the costs incurred at their own site; COST852 provided funding for regular meetings and for scientific networking between partners. The collation of this database was supported by the Irish Research Council for Science, Engineering and Technology through a research fellowship to L. Kirwan, with additional support from Science Foundation Ireland Research Frontiers Programme (09/RFP/EOB2546).
The email addresses of people that have contributed to the data are included in the file site_info.csv. Those in the best position to answer questions concerning the data are Laura Kirwan ([email protected]) and Caroline Brophy ([email protected]); those in the best position to answer questions about experimental protocol are Andreas Lüscher ([email protected]) and Maria-Teresa Sebastià ([email protected]); those in the best position to answer questions concerning the experimental design are John Connolly ([email protected]) and Laura Kirwan ([email protected]).
Introduction
Ecological research on the relationship between plant diversity and ecosystem function generally shows that reductions in plant diversity of randomly-assembled communities reduce the yield of aboveground biomass (Cardinale et al. 2007). Mechanisms that underpin these relationships are attributed to improved utilization of resources in total niche space (niche differentiation), positive interspecific interactions, and selection effects (Hooper et al. 2005). Most such experiments (as reviewed in Cardinale et al. 2007) study the effects of reductions in plant richness from relatively species-rich communities and in low-nutrient systems. In contrast, conventional agriculture in many areas maximises forage yield by planting monocultures of grasses and applying large quantities of nitrogen fertilizer, and has associated negative environmental impacts. We are interested in the role of multi-species agricultural mixtures in improving the resource use efficiency of agricultural grassland systems. If the same diversity-function mechanisms and relationships prevail in simple multi-species plant communities and under higher applications of nutrients, then increased plant diversity in species-poor agricultural systems would be expected to improve their resource-efficient provision of forage and other ecosystem services. For generality and relevance to agricultural practice, investigations of multi-species mixtures need to be conducted at multiple sites, and require comparison against the best-performing component monoculture species (or the prevailing conventional system).
The species composition of multi-species mixtures can be strategically designed to include traits that maximize complementarity and interspecific interactions to improve resource utilization and yield of aboveground biomass. A classic example is the combination of grass and legume species to exploit the ability of legumes to fix atmospheric nitrogen and thereby reduce reliance on application of chemical fertilizers (Nyfeler et al. 2011). However, other traits are also likely to be relevant, but are not usually explicitly tested in multi-species experiments conducted across multiple sites, years and mixture types. Here we test a temporal development trait. Fast-establishing species exhibit fast germination and growth, thereby providing adequate cover of soil and high biomass yields in the first and second year after sowing. These species often lack persistency. Slow-establishing, but persistent species exhibit slower germination and growth rate during establishment, but in the long run are highly competitive, therefore increasing in cover and biomass yields over initial years and constituting the majority of yield in the third and fourth production year and thereafter.
In this data paper, we present data from the Agrodiversity experiment, a co-ordinated continental-scale field experiment across 33 sites that was used to compare the biomass yield of monocultures and four-species mixtures (designed on the basis of specific traits) associated with intensively managed agricultural grassland systems. We addressed the following main questions (modified from Finn et al. 2013):
Recent analyses of some of these data quantified the relationships between community diversity and aboveground biomass (Kirwan et al. 2007, Finn et al. 2013). We showed that aboveground biomass of the mixtures exceeded that of the average monoculture in >97% of comparisons, and mixture biomass also exceeded that of the best monoculture (transgressive overyielding) at about 60% of sites (Finn et al. 2013). There was also a positive relationship between the evenness of the communities and aboveground biomass that was consistent across sites (Kirwan et al. 2007) and persisted for three years (Finn et al. 2013). In an analysis across six sites, positive yield effects were not accompanied by a reduction in herbage digestibility and crude protein concentration (Sturludóttir et al. 2013).
Ultimately, we hope that further analyses of this data set will promote understanding about (a) the design of multi-species mixtures for application in more resource-efficient agricultural systems and (b) how the diversity (composition, richness and relative abundance) of plant communities influences the delivery of ecosystem processes across spatial and temporal scales.
Fig. 1. Map of site locations in Europe. Geographical locations of the sites are indicated for the Mid-European (ME), North-European (NE), Dry-Mediterranean (DM), Moist-Mediterranean (MM) and Other species groups. Note that there is an additional site in Canada. Coordinates for each site are contained in site_info.csv
Methods
Description of the study area and experimental design:
This was a coordinated multi-site experiment in which sites a covered a broad geographical area (Fig. 1). Details of site coordinates, altitude, species and varieties used, plot management and any treatments applied are in site_info.csv.
Design of the core experiment:
The species compositions of experimental plots were selected using a simplex design for mixtures experiments (Cornell [2002] Table 1). At each site, four species were chosen according to four functional traits: a fast-establishing grass (G1), a slow-establishing persistent grass (G2), a fast-establishing legume (L1) and a slow-establishing persistent legume (L2). The species selected by individual sites are detailed in site_info.csv. The species chosen generally fell into four groups according to geographical location: Mid-European (18 sites), North-European (6 sites), Dry Mediterranean (2 sites) and Moist Mediterranean (1 site), with seven sites selecting site-specific species (Other). As the objective was to maximise yield of highly digestible forage from intensive grasslands, the chosen species were all high-yielding agronomic species.
The main treatment consisted of a core set of 30 plots, representing fifteen grassland communities at two levels of seed density. At each site, the four species were sown in monoculture, with sowing rates determined by each site according to local practice. The sowing rates of the four species were systematically varied to produce eleven mixture communities. These fifteen communities were repeated at low and high levels of seeding density ('high' was represented by the standard seeding rate of a monoculture species at a site, with 'low' being 60% of the high seeding rate). Relative proportions of species in a mixture were manipulated by varying seeding rates of the four species at sowing, and resulted in four planned levels of evenness in the design (Table 1). This design resulted in a core set of 30 experimental plots per site (arranged in a completely randomized design).
The species richness of the communities was either one (monocultures) or four (mixtures). However within the four-species mixtures, the level of evenness was manipulated. The levels of evenness (E) were calculated using the formula
,
where s is the maximum number of species in a community (s = 4 for this experiment), and Pi is the sown relative abundance of the ith species (see Kirwan et al. 2007). This lies between 0 for monocultures and 1 for a community in which all species are represented in equal proportions.
Two sites (45 and 46) were not part of the initial network and did not use this core design. They used a variation of the design. The data from these sites is included because there is high overlap in communities with those in the core design. The same plot number is given to design points (community type) at sites 45 and 46 that are the same as those in the core design. This facilitates similar analyses and inclusion in meta-analysis. Sites 45 and 46 included two-species mixtures, whereas the core design does not.
Note on replication: Within a site, the design is not replicated. The design of the experiment was optimized for good coverage of species composition in order to facilitate the use of a response surface regression approach. In regression analyses, one usually does not require replication and residual variation is estimated from the lack of fit of individual points to the regression model selected. However, the use of two seeding densities adds effective replication for mixture communities. The large number of individual experiments across different sites generally adds a very high level of statistical power to the overall experimental design.
Additional treatments:
Additional plots were sown at 22 sites to facilitate the assessment of an experimental treatment. The most commonly applied treatment was a wide genetic base (WGB) treatment, applied at ten sites. A nitrogen fertilizer application treatment was applied at eight sites.
A treatment to investigate genetic diversity was applied at ten sites. At these sites, the treatment plots were sown with legume species with increased intraspecific genetic diversity. The single varieties of white and red clover selected at each site were compared with a wide genetic base (WGB) treatment that consisted of composite populations of white and red clover that were each constructed from commercial varieties plus unselected material obtained from germplasm collections. Different WGB composite populations were constructed for white and red clover species for the ME and NE regions. The seed material for the WGB treatment was supplied to the participating sites and the populations supplied depended on whether the site had used the ME or NE species group for the core experiment. Further details on the composite populations and an analysis of the temporal change in genetic diversity at three sites was published in Collins et al. (2012).
At eight sites, an additional nitrogen fertilizer application level was tested. Details of the amounts of N applied at core and treatment plots are detailed in site_info.csv. The treatments applied at the other sites were cutting frequency, different legume species, and local varieties of the species.
Table 1. Description of experiment design. G1, G2, L1, and L2 represent the sowing proportions of the four species. E is the planned evenness of the community. Density is the sowing density (low is 60% of high) and is determined by each site according to local practice. Type indicates whether the design point is part of the core design, or an additional optional treatment. The core 30 plots make up the main experiment (shaded). The additional 18 treatment plots were established at 18 sites and the treatments applied are detailed in site_info.csv. Plots 49–68 are additional plots that were sown at sites 45 and 46. # Sites indicates the number of sites in which a particular community was sown.
PLOT |
G1 |
G2 |
L1 |
L2 |
E |
Density |
Type |
# Sites |
1 |
0.7 |
0.1 |
0.1 |
0.1 |
0.64 |
High |
Core |
32 |
2 |
0.1 |
0.7 |
0.1 |
0.1 |
0.64 |
High |
Core |
32 |
3 |
0.1 |
0.1 |
0.7 |
0.1 |
0.64 |
High |
Core |
32 |
4 |
0.1 |
0.1 |
0.1 |
0.7 |
0.64 |
High |
Core |
32 |
5 |
0.25 |
0.25 |
0.25 |
0.25 |
1 |
High |
Core |
33 |
6 |
0.4 |
0.4 |
0.1 |
0.1 |
0.88 |
High |
Core |
32 |
7 |
0.4 |
0.1 |
0.4 |
0.1 |
0.88 |
High |
Core |
32 |
8 |
0.4 |
0.1 |
0.1 |
0.4 |
0.88 |
High |
Core |
32 |
9 |
0.1 |
0.4 |
0.4 |
0.1 |
0.88 |
High |
Core |
32 |
10 |
0.1 |
0.4 |
0.1 |
0.4 |
0.88 |
High |
Core |
32 |
11 |
0.1 |
0.1 |
0.4 |
0.4 |
0.88 |
High |
Core |
32 |
12 |
1 |
0 |
0 |
0 |
0 |
High |
Core |
33 |
13 |
0 |
1 |
0 |
0 |
0 |
High |
Core |
33 |
14 |
0 |
0 |
1 |
0 |
0 |
High |
Core |
33 |
15 |
0 |
0 |
0 |
1 |
0 |
High |
Core |
33 |
16 |
0.7 |
0.1 |
0.1 |
0.1 |
0.64 |
Low |
Core |
32 |
17 |
0.1 |
0.7 |
0.1 |
0.1 |
0.64 |
Low |
Core |
32 |
18 |
0.1 |
0.1 |
0.7 |
0.1 |
0.64 |
Low |
Core |
32 |
19 |
0.1 |
0.1 |
0.1 |
0.7 |
0.64 |
Low |
Core |
32 |
20 |
0.25 |
0.25 |
0.25 |
0.25 |
1 |
Low |
Core |
32 |
21 |
0.4 |
0.4 |
0.1 |
0.1 |
0.88 |
Low |
Core |
32 |
22 |
0.4 |
0.1 |
0.4 |
0.1 |
0.88 |
Low |
Core |
32 |
23 |
0.4 |
0.1 |
0.1 |
0.4 |
0.88 |
Low |
Core |
32 |
24 |
0.1 |
0.4 |
0.4 |
0.1 |
0.88 |
Low |
Core |
32 |
25 |
0.1 |
0.4 |
0.1 |
0.4 |
0.88 |
Low |
Core |
32 |
26 |
0.1 |
0.1 |
0.4 |
0.4 |
0.88 |
Low |
Core |
32 |
27 |
1 |
0 |
0 |
0 |
0 |
Low |
Core |
32 |
28 |
0 |
1 |
0 |
0 |
0 |
Low |
Core |
32 |
29 |
0 |
0 |
1 |
0 |
0 |
Low |
Core |
32 |
30 |
0 |
0 |
0 |
1 |
0 |
Low |
Core |
32 |
31 |
0.7 |
0.1 |
0.1 |
0.1 |
0.64 |
Low |
Treatment |
21 |
32 |
0.1 |
0.7 |
0.1 |
0.1 |
0.64 |
Low |
Treatment |
21 |
33 |
0.1 |
0.1 |
0.7 |
0.1 |
0.64 |
Low |
Treatment |
21 |
34 |
0.1 |
0.1 |
0.1 |
0.7 |
0.64 |
Low |
Treatment |
21 |
35 |
0.25 |
0.25 |
0.25 |
0.25 |
1 |
Low |
Treatment |
21 |
36 |
1 |
0 |
0 |
0 |
0 |
Low |
Treatment |
16 |
37 |
0 |
1 |
0 |
0 |
0 |
Low |
Treatment |
16 |
38 |
0 |
0 |
1 |
0 |
0 |
Low |
Treatment |
21 |
39 |
0 |
0 |
0 |
1 |
0 |
Low |
Treatment |
21 |
40 |
0.7 |
0.1 |
0.1 |
0.1 |
0.64 |
Low |
Treatment |
21 |
41 |
0.1 |
0.7 |
0.1 |
0.1 |
0.64 |
Low |
Treatment |
21 |
42 |
0.1 |
0.1 |
0.7 |
0.1 |
0.64 |
High |
Treatment |
21 |
43 |
0.1 |
0.1 |
0.1 |
0.7 |
0.64 |
High |
Treatment |
21 |
44 |
0.25 |
0.25 |
0.25 |
0.25 |
1 |
High |
Treatment |
21 |
45 |
1 |
0 |
0 |
0 |
0 |
High |
Treatment |
16 |
46 |
0 |
1 |
0 |
0 |
0 |
High |
Treatment |
16 |
47 |
0 |
0 |
1 |
0 |
0 |
High |
Treatment |
21 |
48 |
0 |
0 |
0 |
1 |
0 |
High |
Treatment |
21 |
49 |
0 |
0.5 |
0 |
0.5 |
0.6667 |
High |
Additional |
3 |
50 |
0 |
0 |
0.5 |
0.5 |
0.6667 |
High |
Additional |
3 |
51 |
0.5 |
0.5 |
0 |
0 |
0.6667 |
High |
Additional |
3 |
52 |
0.5 |
0 |
0.5 |
0 |
0.6667 |
High |
Additional |
3 |
53 |
0.5 |
0 |
0 |
0.5 |
0.6667 |
High |
Additional |
3 |
54 |
0 |
0.5 |
0.5 |
0 |
0.6667 |
High |
Additional |
3 |
55 |
0 |
0.5 |
0 |
0.5 |
0.6667 |
Low |
Additional |
2 |
56 |
0 |
0 |
0.5 |
0.5 |
0.6667 |
Low |
Additional |
2 |
57 |
0.5 |
0.5 |
0 |
0 |
0.6667 |
Low |
Additional |
2 |
58 |
0.5 |
0 |
0.5 |
0 |
0.6667 |
Low |
Additional |
2 |
59 |
0.5 |
0 |
0 |
0.5 |
0.6667 |
Low |
Additional |
2 |
60 |
0 |
0.5 |
0.5 |
0 |
0.6667 |
Low |
Additional |
2 |
61 |
0.88 |
0.04 |
0.04 |
0.04 |
0.2944 |
High |
Additional |
2 |
62 |
0.04 |
0.88 |
0.04 |
0.04 |
0.2944 |
High |
Additional |
2 |
63 |
0.04 |
0.04 |
0.88 |
0.04 |
0.2944 |
High |
Additional |
2 |
64 |
0.04 |
0.04 |
0.04 |
0.88 |
0.2944 |
High |
Additional |
2 |
65 |
0.88 |
0.04 |
0.04 |
0.04 |
0.2944 |
Low |
Additional |
2 |
66 |
0.04 |
0.88 |
0.04 |
0.04 |
0.2944 |
Low |
Additional |
2 |
67 |
0.04 |
0.04 |
0.88 |
0.04 |
0.2944 |
Low |
Additional |
1 |
68 |
0.04 |
0.04 |
0.04 |
0.88 |
0.2944 |
Low |
Additional |
1 |
Management protocol:
Plots were not grazed. The number of cuts per annum and fertilizer application levels were determined by local practice at individual sites. See Site_info.csv for details on management practices employed at each site.
The plots were not weeded and there was generally no herbicide application. However, some targeted weeding was required in some sites during the establishment phase. See Site_info.csv for details on weeding in the establishment year. Margins between plots were sprayed to prevent stoloniferous ingression from neighboring plots.
Year 1 was defined as the first complete year after sowing. Some sites took cleaning cuts in the year of establishment (prior to year 1), but data is not recorded for these cuts. See Site_info.csv for details on cleaning cuts.
Response variables:
The response variables included in the database are total biomass, biomass of the five harvest fractions (G1, G2, L1, L2, weed) and measurements of forage quality. Analyses of total plot biomass and biomass of weed species have been published in Kirwan et al. (2007), Lüscher et al. (2008), Frankow-Lindberg et al. (2009), Nyfeler et al. (2009, 2011) and Finn et al. (2013). An analysis of forage quality at 4 sites was published in Sturludóttir et al. (2013).
Total Biomass per plot
For each plot, biomass of aboveground vegetation was measured at each harvest. This was done by cutting the whole plot and determining the fresh weight of the 'whole plant material'. A subsample of this material was taken, its fresh weight determined and the material dried at 65°C to constant weight. From the dry weight of the sample the percentage dry matter was calculated. From this, the total dry matter yield for the plot (DM/m²) was calculated from the fresh mass of the 'whole plant material'.
Biomass of the five harvest fractions (G1, G2, L1, L2, weed)
Biomass separation was done by one of the following two methods: (see Site_info.csv for method selected for each site). The harvests at which separate determination was carried out are given in biomass.csv.
Note on indigenous plants of experimental species (G1, G2, L1, L2): In the mixture plots indigenous (not sown) G1, G2, L1, and L2 plants cannot be separated from the sown plants of these species and so do not form part of the weed fraction. However, in each monoculture plot, indigenous plants of the other three species are included in weed.
Forage quality
The primary forage quality analysis was carried out at the Christian-Albrechts-Universität Kiel, Kiel, Germany. Additional analysis for four North European sites was carried out at Agriculture and Agrifood Canada at Levis, Canada. This analysis of samples from the North-European sites has been published in Sturludóttir et al. (2013). In addition, N concentration of bulk samples was locally analyzed by 21 sites. Forage quality data is given in forage_quality.csv. Measurements are coded _K, _C or _L to indicate whether samples were analysed at Kiel, Canada or local laboratories.
Kiel analysis: Bulk samples (not separated into component species) were analysed for 17 sites. At seven of the sites, all plots were analyzed. At 10 sites, samples were analyzed for 18 out of the 30 experimental plots (the 12 plots in the design co-dominated by two species were omitted). Separation of the forage into the five fractions was also considered for those particularly interested in the forage quality of the different functional groups. Separated samples were analyzed for eight sites. Sample material of at least 5 g per sample was prepared by drying to a constant weight at 65°C and then grinding to pass through a sieve of 1-mm mesh size. Near-infrared spectroscopy (NIRS) analysis was carried out to determine the dry matter percentage of nitrogen concentration (N), ADF (acid detergent fibre) NDF (neutral detergent fibre), ELOS (enzymatic soluble organic dry matter), and ash which represents the mineral content of the sample.
Canadian analysis: Bulk samples were analyzed for four sites. Samples were analysed for all plots. Samples were dried and ground and then analyzed using NIRS (FOSS NIRsystems 6500, Silver Spring, MD) to determine nitrogen concentration (N), acid detergent fiber (ADF), neutral detergent fiber (NDF), in vitro true digestibility (IVTD) and in vitro cell wall digestibility (IVCWD). The latter was calculated using the following equation: IVCWD = 1000 - [(1000 – IVTD)/(NDF/1000)].
Site description variables:
Climate data
Each site provided daily climate data for the duration of their experiment. Where possible, the climate data series begins on the sowing date of the experiment at each site. The variables recorded were precipitation (mm per day), minimum daily air temperature (°C), mean daily air temperature (°C) and maximum daily air temperature (°C). The climate data is contained in climate.csv.
Soil analysis
The soil analysis was carried out at Centre Tecnològic Forestal de Catalunya, Solsona, Spain. Composite soil samples were formed by combining samples from four plots. In each sampled plot a soil volume of 5 × 5 cm² per 15 cm depth was taken in a systematic manner, as follows. To avoid soil contamination by litter, the first 0.5 cm of litter-soil was removed. The composite sample was dried (temperature between 20 and 40ºC) over two or three days. The soil aggregates were gently ground and passed through a 2-mm sieve, discarding the fraction >2 mm. A 300 g subsample was sent to the laboratory for analysis. The percentage of sand, silt, fine silt, and clay were measured and the soil type classified. In addition, percent organic matter, soil carbonates, soil electrical conductivity and soil pH were measured, along with the concentrations of calcium, potassium, nitrate, magnesium, sodium, and phosphorus.
Data set description
Data from individual sites was recorded in standardized data recording spreadsheets and then collated into a central database. Data was checked using numerical and graphical summary methods. Data range was checked for each variable and pivot tables were used to check counts of measurements. The datasets were manipulated and merged using the SORT, MERGE, and DATA procedures in SAS 9.1 (SAS Institute Inc., Cary, NC, USA).
Descriptions and units of measurement for the columns of data are presented in the tables below. Site-level information is contained in the comma-separated-value data files named site_info.csv, climate.csv and soils.csv. Data relating to plot-level measurements are contained in the comma-separated-value data files named biomass.csv and forage_quality.csv. Note that in all files missing data values are represented by empty cells.
Column numbers, variable names and variable descriptions for file: site_info.csv
Column |
Variable Name |
Variable Description |
Unit |
1 |
SITE |
Site ID number |
|
2 |
Country |
Country |
|
3 |
Location |
Location of site within country |
|
4 |
Institute |
Institute responsible for site |
|
5 |
Contact Email |
Contact email of individual responsible for site |
|
6 |
Lat_Deg |
Location of site (Latitude degrees) |
|
7 |
Lat_Min |
Location of site (Latitude minutes) |
|
8 |
Lat_NS |
Location of site (Latitude North South) |
|
9 |
Long_Deg |
Location of site (Longitude degrees) |
|
10 |
Long_Min |
Location of site (Longitude minutes) |
|
11 |
Long_EW |
Location of site (Longitude East West) |
|
12 |
Elevation |
Elevation of site |
m above sea level |
13 |
Mixture Type |
Seed mixture used: ME=Mid European, NE=Northern European, |
|
14 |
G1 Species |
Fast establishing grass species |
|
15 |
G1 Variety |
Fast establishing grass variety |
|
16 |
G2 Species |
Persistent grass species |
|
17 |
G2 Variety |
Persistent grass variety |
|
18 |
L1 Species |
Fast establishing legume species |
|
19 |
L1 Variety |
Fast establishing legume variety |
|
20 |
L2 Species |
Persistent legume species |
|
21 |
L2 Variety |
Persistent legume variety |
|
22 |
Sowing Date |
Date the plots were established |
|
23 |
Sowing Method |
Method of sowing – drilled / hand sown |
|
24 |
P at sowing |
P fertilizer applied at establishment |
Kg/ha |
25 |
K at sowing |
K fertilizer applied at establishment |
Kg/ha |
26 |
N at sowing |
N fertilizer applied at establishment |
Kg/ha |
27 |
Annual P |
P fertilizer applied per annum |
Kg/ha |
28 |
Annual K |
K fertilizer applied per annum |
Kg/ha |
29 |
Annual N |
N fertilizer applied per annum |
Kg/ha |
30 |
Year 1 number of harvests |
Number of harvests taken in the year 1 of the experiment (year after sowing) |
|
31 |
Year 2 number of harvests |
Number of harvests taken in the year 2 of the experiment |
|
32 |
Year 3 number of harvests |
Number of harvests taken in the year 3 of the experiment |
|
33 |
Year 4 number of harvests |
Number of harvests taken in the year 4 of the experiment |
|
34 |
Harvest Height |
Cutting height when taking harvests |
cm |
35 |
Method of separation |
Method of selecting subsample of biomass for separation into species components |
|
36 |
Plot Size |
Size of plots |
m² |
37 |
Area Sampled for Total Yield |
Area within plot that was sampled for total aboveground biomass |
m² |
38 |
Area Sampled for Composition |
Area within plot that was sampled for separation into species components (if fixed quadrat method was used) |
m² |
39 |
Harvesting Method |
Method used to harvest biomass (manually or by machine) |
|
40 |
Treatment details |
Details of treatment(s) applied (where relevant) |
|
41 |
Cleaning cut date |
Date of cleaning cut (if any) |
|
42 |
Weeding details |
Details of any weeding undertaken during establishment |
|
Column numbers, variable names, and variable descriptions for file: climate.csv
Column |
Variable Name |
Variable Description |
Unit |
1 |
Site |
Site ID number |
|
2 |
Day |
Day |
|
3 |
Month |
Month |
|
4 |
Year |
Year |
|
5 |
Date |
Date |
|
6 |
Precip |
Daily precipitation |
mm/ay |
7 |
air_min |
Minimum daily air temperature |
°C |
8 |
air_mean |
Mean daily air temperature |
°C |
9 |
air_max |
Maximum daily air temperature |
°C |
Column numbers, variable names, and variable descriptions for file: soils.csv
Column |
Variable Name |
Variable Description |
Unit |
1 |
SITE |
Site ID number |
|
2 |
Carbonates |
Soil Carbonates |
% |
3 |
EC |
Soil electrical conductivity |
ds m-1 |
4 |
Silt |
Percent silt content in soil |
% |
5 |
Silt_fine |
Percent fine silt content in soil |
% |
6 |
Clay |
Percent clay content in soil |
% |
7 |
Sand |
Percent sand content in soil |
% |
8 |
OM |
Percent organic matter |
% |
9 |
Soil_Type |
Soil type |
|
10 |
Humidity |
Percent humidity |
% |
11 |
Ca |
Calcium concentration |
ppm |
12 |
K |
Potassium concentration |
ppm |
13 |
N-NO3 |
Nitrate concentration |
ppm |
14 |
Mg |
Magnesium concentration |
ppm |
15 |
Na |
Sodium concentration |
ppm |
16 |
P |
Phosphorus concentration |
ppm |
17 |
pH |
Soil pH |
|
Column numbers, variable names, and variable descriptions for file: biomass.csv
Column |
Variable Name |
Variable Description |
Unit |
1 |
SITE |
Site ID number |
|
2 |
COUNTRY |
Country |
|
3 |
YEAR |
Year |
|
4 |
YEARN |
Experimental year number |
|
5 |
NH |
Number of harvests – number of times the whole plot was cut in a year |
|
6 |
HARVEST |
Harvest number (within year) |
|
7 |
HARVEST_DATE |
Date of harvest |
|
8 |
PLOT |
Plot number as per design (1–30 = core design; 31–48 = treatment plots; 49–68 = additional plots at sites 45 and 46) |
|
9 |
TREAT |
Indicator variable: 1=basic 30 plots; 2 and 3=additional treatment plots (some sites implemented two levels of additional treatments) |
|
10 |
REP |
Replicate number (applies only to sites 15, 45, and 46) |
|
11 |
G1 |
Initial sown proportion of fast-establishing grass |
|
12 |
G2 |
Initial sown proportion of persistent grass |
|
13 |
L1 |
Initial sown proportion of fast-establishing legume |
|
14 |
L2 |
Initial sown proportion of persistent legume |
|
15 |
E |
Initial sown evenness |
|
16 |
DENS |
Indicator variable: high=high level of initial sown biomass, low = low level (60%of high) |
|
17 |
G1_Y |
Harvest Dry Matter Yield of fast-establishing grass |
t/ha |
18 |
G2_Y |
Harvest Dry Matter Yield of persistent grass |
t/ha |
19 |
L1_Y |
Harvest Dry Matter Yield of fast-establishing legume |
t/ha |
20 |
L2_Y |
Harvest Dry Matter Yield of persistent legume |
t/ha |
21 |
WEED_Y |
Harvest Dry Matter Yield of weed species |
t/ha |
22 |
HARV_YIELD |
Total Harvest Dry Matter Yield |
t/ha |
Column numbers, variable names, and variable descriptions for file: forage_quality.csv
Column |
Variable Name |
Variable Description |
Unit |
1 |
SITE |
Site ID number |
|
2 |
country |
Country |
|
3 |
YEAR |
Year |
|
4 |
YEARN |
Experimental year number |
|
5 |
NH |
Number of harvests |
|
6 |
Harvest |
Harvest number (within year) |
|
7 |
HARVEST_DATE |
Date of harvest |
|
8 |
PLOT |
Plot number as per design (1–30 = core design; 31–48 = treatment plots; 49–68 = additional plots at sites 45 and 46) |
|
9 |
TREAT |
Indicator variable: 1=basic 30 plots; 2 and 3=additional treatment plots (some sites implemented two levels of additional treatments) |
|
10 |
REP |
Replicate number (applies only to sites 15 and 45) |
|
11 |
G1 |
Initial sown proportion of fast-establishing grass |
|
12 |
G2 |
Initial sown proportion of persistent grass |
|
13 |
L1 |
Initial sown proportion of fast-establishing legume |
|
14 |
L2 |
Initial sown proportion of persistent legume |
|
15 |
E |
Initial sown evenness |
|
16 |
DENS |
Indicator variable: high=high level of initial sown biomass, low = low level (60% of high) |
|
17 |
Local_N |
Indicator variable (Local lab analysis present=1, absent=0) |
|
18 |
Kiel |
Indicator variable (Kiel bulk sample present=1 , absent=0) |
|
19 |
Canada |
Indicator variable ( Canada bulk sample present=1 , absent=0) |
|
20 |
Kiel_sep |
Indicator variable ( Kiel separated sample present=1 , absent=0) |
|
21 |
N_L |
Nitrogen percent in total harvest yield (analysis performed by local lab) |
% of dry matter |
22 |
N_K |
Nitrogen percent in total harvest yield (Kiel data) |
% of dry matter |
23 |
Ash_K |
Ash in total harvest yield (Kiel data) |
% of dry matter |
24 |
NDF_K |
Neutral Detergent Fibre in total harvest yield (Kiel data) |
% of dry matter |
25 |
ADF_K |
Acid Detergent Fibre in total harvest yield (Kiel data) |
% of dry matter |
26 |
CDOMD_K |
Cellulase Digestible of Organic Matter of Dry Matter in total harvest yield (Kiel data) |
% of dry matter |
27 |
ME_K |
Metabolizable Energy in total harvest yield (Kiel data) |
MJ ME per kg of DM |
28 |
N_C |
Nitrogen percent in total harvest yield (Canadian data) |
% of dry matter |
29 |
NDF_C |
Neutral Detergent Fibre in total harvest yield (Canadian data) |
% of dry matter |
30 |
ADF_C |
Acid Detergent Fibre in total harvest yield (Canadian data) |
% of dry matter |
31 |
IVTD_C |
In Vitro True Digestibility in total harvest yield (Canadian data) |
|
32 |
IVCWD_C |
In Vitro Cell Wall Digestibility in total harvest yield (Canadian data) |
|
33 |
N_G1_K |
Nitrogen percent in G1 harvest yield (Kiel data) |
% of dry matter |
34 |
Ash_G1_K |
Ash in total G1 harvest yield (Kiel data) |
% of dry matter |
35 |
NDF_G1_K |
Neutral Detergent Fibre in total G1 harvest yield (Kiel data) |
% of dry matter |
36 |
ADF_G1_K |
Acid Detergent Fibre in total G1 harvest yield (Kiel data) |
% of dry matter |
37 |
CDOMD_G1_K |
Cellulase Digestible of Organic Matter of Dry Matter in G1 harvest yield (Kiel data) |
% of dry matter |
38 |
ME_G1_K |
Metabolizable Energy in G1 harvest yield (Kiel data) |
MJ ME per kg of DM |
39 |
N_G2_K |
Nitrogen percent in G2 harvest yield (Kiel data) |
% of dry matter |
40 |
Ash_G2_K |
Ash in G2 harvest yield (Kiel data) |
% of dry matter |
41 |
NDF_G2_K |
Neutral Detergent Fibre in G2 harvest yield (Kiel data) |
% of dry matter |
42 |
ADF_G2_K |
Acid Detergent Fibre in G2 harvest yield (Kiel data) |
% of dry matter |
43 |
CDOMD_G2_K |
Cellulase Digestible of Organic Matter of Dry Matter in G2 harvest yield (Kiel data) |
% of dry matter |
44 |
ME_G2_K |
Metabolizable Energy in G2 harvest yield (Kiel data) |
MJ ME per kg of DM |
45 |
N_L1_K |
Nitrogen percent in L1 harvest yield (Kiel data) |
% of dry matter |
46 |
Ash_L1_K |
Ash in L1 harvest yield (Kiel data) |
% of dry matter |
47 |
NDF_L1_K |
Neutral Detergent Fibre in L1 harvest yield (Kiel data) |
% of dry matter |
48 |
ADF_L1_K |
Acid Detergent Fibre in L1 harvest yield (Kiel data) |
% of dry matter |
49 |
CDOMD_L1_K |
Cellulase Digestible of Organic Matter of Dry Matter in L1 harvest yield (Kiel data) |
% of dry matter |
50 |
ME_L1_K |
Metabolizable Energy in L1 harvest yield (Kiel data) |
MJ ME per kg of DM |
51 |
N_L2_K |
Nitrogen percent in L2 harvest yield (Kiel data) |
% of dry matter |
52 |
Ash_L2_K |
Ash in L2 harvest yield (Kiel data) |
% of dry matter |
53 |
NDF_L2_K |
Neutral Detergent Fibre in L2 harvest yield (Kiel data) |
% of dry matter |
54 |
ADF_L2_K |
Acid Detergent Fibre in L2 harvest yield (Kiel data) |
% of dry matter |
55 |
CDOMD_L2_K |
Cellulase Digestible of Organic Matter of Dry Matter in L2 harvest yield (Kiel data) |
% of dry matter |
56 |
ME_L2_K |
Metabolizable Energy in L2 harvest yield (Kiel data) |
MJ ME per kg of DM |
57 |
N_weed_K |
Nitrogen percent in weed species harvest yield (Kiel data) |
% of dry matter |
58 |
Ash_weed_K |
Ash in weed species harvest yield (Kiel data) |
% of dry matter |
59 |
NDF_weed_K |
Neutral Detergent Fibre in weed species harvest yield (Kiel data) |
% of dry matter |
60 |
ADF_weed_K |
Acid Detergent Fibre in weed species harvest yield (Kiel data) |
% of dry matter |
61 |
CDOMD_weed_K |
Cellulase Digestible of Organic Matter of Dry Matter in weed species harvest yield (Kiel data) |
% of dry matter |
62 |
ME_weed_K |
Metabolizable Energy in weed species harvest yield (Kiel data) |
MJ ME per kg of DM |
DATA-USE POLICY
The data presented here are publicly available. Those wishing to publish results from this data set should read this metadata document. The data set should be cited as: Kirwan et al. 2014. The Agrodiversity Experiment: three years of data from a multisite plant diversity experiment in intensively managed grasslands. Ecology xx:xxx.
Three papers are currently in preparation, focusing on (a) the effect of diversity on the contribution of weed species to biomass yield, (b) changes in the relative abundances of the mixtures, and (c) grassland biodiversity effects across an environmental gradient.
Literature cited
Cardinale, B. J., J. P. Wright, M. W. Cadotte, I. T. Carroll, A. Hector, D. S. Srivastava, M. Loreau, and J. J. Weis. 2007. Impacts of plant diversity on biomass production increase through time because of species complementarity. Proceedings of the National Academy of Sciences of the United States of America 104:18123–18128.
Collins, R. P., A. Helgadóttir, B. E. Frankow-Lindberg, L. Skøt, C. Jones, and K. P. Skøt. 2012. Temporal changes in population genetic diversity and structure in red and white clover grown in three contrasting environments in northern Europe. 2012. Annals of Botany 110(6):1341–1350.
Connolly, J., J. A. Finn, A. D. Black, L. Kirwan, C. Brophy, and A. Lüscher. 2009. Effects of multi-species swards on biomass production and weed invasion at three Irish sites. Irish Journal of Agricultural and Food Research. 48:243–260.
Cornell, J. A. 2002. Experiment with Mixtures: Designs, Models, and the Analysis of Mixture Data, Third Edition, John Wiley & Sons, Inc., USA.
Finn, J. A., L. Kirwan, J. Connolly, M. T. Sebastià, A. Helgadottir, O. H. Baadshaug, G. Bélanger, A. Black, C. Brophy, R. P. Collins, J. Čop, S. Dalmannsdóttir, I. Delgado,A. Elgersma, M. Fothergill, B. E. Frankow-Lindberg, A. Ghesquiere, B. Golinska, P. Golinski, P. Grieu, A. -M. Gustavsson, M. Höglind, O. Huguenin-Elie, M. Jørgensen, Z. Kadziuliene, P. Kurki, R. Llurba, T. Lunnan, C. Porqueddu, M. Suter, U. Thumm, and A. Lüscher. 2013. Ecosystem function enhanced by combining four functional types of plant species in intensively managed grassland mixtures: a 3-year continental-scale field experiment. Journal of Applied Ecology 50:365–375.
Frankow-Lindberg, B. E., C. Brophy, R. P. Collins, and J. Connolly. 2009. Biodiversity effects on yield and unsown species invasion in a temperate forage ecosystem. Annals of Botany 103(6):913–921.
Hooper, D. U., F. S. Chapin, J. J. Ewel, A. Hector, P. Inchausti, S. Lavorel, J. H. Lawton, D. M. Lodge, M. Loreau, S. Naeem, B. Schmid, H. Setala, A.J. Symstad, J. Vandermeer, and D. A. Wardle. 2005. Effects of biodiversity on ecosystem functioning: A consensus of current knowledge. Ecological Monographs 75:3–35.
Isbell, F., V. Calcagno, A. Hector, J. Connolly, W. S. Harpole, P. B. Reich, M. Scherer-Lorenzen, B. Schmid, D. Tilman, J. Van Ruijven, et al. 2011. High plant diversity is needed to maintain ecosystem services. Nature 477(7363):199–202.
Kirwan, L., A. Luescher, M. T. Sebastia, J. A. Finn, R. P. Collins,C. Porqueddu, A. Helgadottir, O. H. Baadshaug, C. Brophy, C. Coran, S. Dalmannsdottir, I. Delgado, A. Elgersma, M. Fothergill, B. E. Frankow-Lindberg, P. Golinski, P. Grieu, A. M. Gustavsson, M. Hoglind, O. Huguenin-Elie, C. Iliadis, M. Jorgensen, Z. Kadziuliene, T. Karyotis, T. Lunnan, M. Malengier, S. Maltoni, V. Meyer, D. Nyfeler, P. Nykanen-Kurki, J. Parente, H. J. Smit, U. Thumm, and J. Connolly. 2007. Evenness drives consistent diversity effects in intensive grassland systems across 28 European sites. Journal of Ecology. 95:530–539.
Lüscher, A., J. A. Finn, J. Connolly, M. T. Sebastià, R. Collins, M. Fothergill, C. Porqueddu, C. Brophy, O. Huguenin-Elie, L. Kirwan, D. Nyfeler, and A. Helgadóttir. 2008. Benefits of sward diversity for Agricultural grasslands. Biodiversity 9:29–32.
Nyfeler, D., O. Huguenin-Elie, M. Suter, E. Frossard,J. Connolly, and A. Lüscher, A. 2009. Strong mixture effects among four species in fertilized agricultural grassland led to persistent and consistent transgressive overyielding. Journal of Applied Ecology 46:683–691.
Nyfeler, D., O. Huguenin-Elie, M. Suter, E. Frossard, and A. Lüscher. 2011. Grass–legume mixtures can yield more nitrogen than legume pure stands due to mutual stimulation of nitrogen uptake from symbiotic and non-symbiotic sources. Agriculture, Ecosystems & Environment 140:155–163.
Sturludóttir, E., C. Brophy, G. Bélanger, A. -M. Gustavsson, M. Jørgensen, T. Lunnan, and Á. Helgadóttir. 2013. Benefits of mixing grasses and legumes for herbage yield and nutritive value in Northern Europe and Canada. Grass and Forage Science. doi: 10.1111/gfs.12037