Skip to main content

Disentangling direct and indirect effects of experimental grassland management and plant functional-group manipulation on plant and leafhopper diversity



Plant biodiversity can affect trophic interactions in many ways, including direct bottom-up effects on insects, but is negatively affected by agricultural intensification. Grassland intensification promotes plant productivity, resulting in changes in plant community composition, and impacts on higher trophic levels. Here, we use a novel grassland management experiment combining manipulations of cutting and fertilization with experimental changes in plant functional group composition (independent of management effects) to disentangle the direct and indirect effects of agricultural management on insect herbivore diversity and abundance. We used leafhoppers as model organisms as they are a key insect taxon in grasslands and react rapidly to management changes. Leafhoppers were sampled between May and September 2010 using standardized sweep netting and pan traps.


Plant diversity, functional group composition and management regime in grasslands affected leafhopper species richness and abundance. Higher cutting frequencies directly led to decreasing leafhopper species richness, presumably due to the higher disturbance frequency and the reduction in food-resource heterogeneity. In contrast, fertilizer application had only a small indirect negative effect via enhanced aboveground plant biomass, reduced plant diversity and changes in functional group composition. The manipulated increase in grass cover had contrasting direct and indirect effects on leafhopper species richness: grass cover directly increased leafhopper species richness, but negatively affected plant diversity, which in turn was positively related to leafhopper species richness. In conclusion, insect diversity is driven in complex direct and indirect ways by grassland management, including changes in functional group composition.


The availability of preferred food sources and the frequency of disturbance are important direct and indirect drivers of leafhopper species richness, interacting in complex ways with plant diversity and food resource heterogeneity.


Grasslands, such as permanent meadows and pastures, cover about 37% of the agricultural area in Europe[1] and harbour much of Europe’s overall biodiversity[2]. Many plant and animal species are restricted to this habitat type[3]. However, agricultural intensification and land-use change have caused major losses in grassland biodiversity[4, 5]. Large amounts of fertilizer are applied in grassland to increase yield[6], allowing for earlier and more frequent cuttings in the growing season. Additionally, herbicides may be used to suppress unwanted plant species ([3]). These management practices greatly affect both plant biomass and species composition. Plant species that are adapted to low nutrient levels and low cutting frequencies are replaced by more competitive, faster-growing species[7, 8], thereby altering the invertebrate communities of the grassland as well[9]. Frequent cutting disturbs the vegetation structure, removes food resources and kills many animals[3, 1012].

Here, we present results from a novel grassland management and plant functional group manipulation experiment (GrassMan)[13], combining experimental variations in cutting frequency (two levels, one cut or three cuts) and fertilizer application (two levels, fertilized or unfertilized) with a plant functional group manipulation treatment (three levels), testing the enhancement of grasses or herbs independent of management changes. We focus on insect responses to plant species composition and management intensity. The resulting 12 treatment combinations were replicated six times, resulting in 72 plots laid out experimentally as a Latin rectangle[10].

Other studies on trophic interactions in biodiversity experiments have largely relied on artificially sown gradients in plant diversity[1416]. To achieve more realistic results, we performed this study in an old grassland and only changed the relative importance of grasses and forbs[17].

We chose leafhoppers, planthoppers and froghoppers (Auchenorrhyncha, hereafter referred to as leafhoppers), as model organisms as they are a highly diverse plant-sucking insect group that has been shown to be strongly influenced by management regime, productivity, vegetation structure and plant species composition[1820]. Leafhoppers play an important role both as herbivores and as prey for higher trophic levels. Their rapid reaction to changes in management regime makes them highly appropriate for ecological studies such as the one presented here[19]. Nevertheless, they have rarely been studied in this context (but see Hollier et al.[21]).

We hypothesize that

  1. (1)

    Leafhopper abundance and diversity decrease with cutting frequency, as cutting acts as a mechanical disturbance.

  2. (2)

    Leafhopper abundance and diversity increase with fertilizer application, as this enhances quantity and nutritional quality of available food resources.

  3. (3)

    Leafhopper abundance and diversity increase under experimental enhancement of graminoid cover, because many species feed preferably on graminoids.


Treatment effects on vegetation

We found 61 plant species, 22 graminoids and 39 forbs (including four legumes). Plant Shannon diversity was positively affected by higher cutting frequency (three times/year) (Figure 1a; Table 1). In addition, FG manipulation significantly affected plant Shannon diversity (Figure 1a; Table 1). Plant Shannon diversity (numbers equivalent, eH’) was highest in unfertilized, forb-rich plots with three cuts (9.61 ± 0.35) and lowest in fertilized, graminoid-rich plots with one cut (5.28 ± 0.45; Table 2), but fertilizer effects were generally non-significant (Table 1).

Figure 1
figure 1

Changes in plant diversity due to functional group manipulation, fertilization and mowing frequency. (a) Effects of FG manipulation and mowing frequency on the numbers equivalent (e) of plant Shannon diversity; only main effects were significant; (b) Effects of proportion of graminoids on e; Spearman’s rho = -0.70; p < 0.001; the regression line represents a linear model (F1,72 = 69.43; P < 0.001); (c) Effect of FG manipulation and fertilization on the proportion of graminoids. Bars represent the mean; error bars show ±1 SE.

Table 1 Sequential F tests of model terms for generalized linear models with plant and leafhopper diversity as response variables
Table 2 Effects of experimental treatments on vegetation parameters (mean ± 1 SE)

Plant Shannon diversity was strongly negatively correlated with the proportion of graminoids (Spearman’s rho = -0.70; p < 0.001, Figure 1b, Table 2).

The proportion of graminoids was highest in fertilized, graminoid-rich plots with one cut (86.67 ± 2.35 t ha-1) and lowest in unfertilized, forb-rich plots with three cuts (46.33 ± 5.07 t ha-1). Generalized linear models with logit-transformed proportion of graminoids (Table 1) showed that the proportion of graminoids increased 0.42-fold with fertilizer application (estimate: 0.42 ± 0.08, t = 4.95, df = 68,P < 0.001). Forb-rich communities had a significantly lower proportion of grasses than control plots (estimate: -0.44 ± 0.1, t = -4.27, df = 68, P < 0.001), and graminoid-rich plots had significantly more grasses than forb-rich plots (estimate: 0.93 ± 0.10, t = 9.04, df = 68, P < 0.001). The positive effect of fertilization on the proportion of graminoids was even higher under FG reduction of forbs (see also Everwand et al.[10]).

Treatment effects on leafhoppers

In total, we caught 6497 adult leafhopper specimens from 36 species (Additional file1: Table S4). Twenty-eight species (86.5% of the total abundance) were graminoid-feeders, and eight species were forb-feeders[22]. The four most common species were Arthaldeus pascuellus (Fall.) with 67.9% of total abundance, Philaenus spumarius (L.) (10.7%), Streptanus sordidus (ZETT.) (7.6%) and Macrosteles viridigriseus (EDWARDS) (2.7%).

Leafhopper Shannon diversity (e) was highest in the most natural and least disturbed plots (one cut, no herbicide) with an average of 4.02 ± 0.38 species per plot (Figure2a, Additional file2: Table S2). A consistently negative effect of three cuts per year (Figure 2a, Table 1, Additional file2: Table S2) and of FG manipulation (Figure 2a, Table 1, Additional file2: Table S2) was observed. Fertilizer application had no direct effect on leafhopper Shannon diversity. Untransformed leafhopper species richness was positively affected by the proportion of graminoids (F1,72 = 6.56; P = 0.013; Figure 2b). For a numerical summary of effects on leafhoppers, see Additional file3: Table S5.

Figure 2
figure 2

Changes in leafhopper diversity due to functional group manipulation and mowing frequency. (a) Effects of FG manipulation and cutting frequency on the numbers equivalent of Shannon diversity of leafhoppers (e); (b) Positive relationship between leafhopper diversity [species per plot] and proportion of graminoids [%] (Spearman’s rho = 0.35; p = 0.003); the regression line represents a linear model (F1,72 = 6.56; P = 0.013). Bars represent the mean; error bars show ±1 SE.

Direct and indirect treatment effects on vegetation

Structural equation modeling (Figure 3) showed that an increasing proportion of graminoids, due to FG manipulation, reduced plant diversity (standardized path coefficient ß = -0.61). Plant diversity (latent variable) was mainly driven by forb diversity (ß = 0.71). Fertilizer application and cutting frequency jointly influenced the latent variable “plant productivity”: fertilization increased plant biomass and height, while cutting reduced both. Finally, higher plant productivity negatively influenced plant diversity (ß = -0.2). Notably, alternative pathways, e.g. arrows from plant diversity to plant productivity, were not supported by our data.

Figure 3
figure 3

Structural equation model showing indirect and direct treatment effects on leafhoppers. The graph shows the minimal adequate structural equation model with N = 72, X2 = 18.7, P = 0.664, 22 degrees of freedom and a root mean squared error of approximation of 0.00 (90% confidence interval = [0, 0.081]). Rectangles represent observed variables (organism abundance and diversity = species richness), ellipses represent latent variables. Solid (dashed) arrows indicate positive (negative) relationships among variables. Numbers next to arrows and boxes are standardized path coefficients. Design variables were specified as numeric variables. Fertilizer application: no fertilizer = 0; NPK-fertilizer = 1; cutting frequency; one cut per year = 0; three cuts per year = 1; and Graminoid enhancement: graminoid-reduced = -1; control = 0; forb-reduced = 1.

Direct and indirect effects of treatments and vegetation on leafhoppers

Our structural equation model (Figure 3) indicated that there were both direct and indirect effects of treatments on leafhopper abundance and species richness.Plant functional group manipulation was the most important direct driver of leafhopper species richness (ß = 0.53): In forb-rich communities, leafhopper species richness was lowest, while it was highest in grass-rich communities. In contrast, a higher percentage of grasses affected plant diversity negatively, while the number of plant species (latent variable “plant diversity”) had a positive effect on leafhopper species richness (ß = 0.62, see Figure 3). Cutting frequency had a negative effect on leafhopper species richness (ß = -0.42), whereas fertilizer application exhibited no effect on leafhopper species richness. Finally, leafhopper species richness was strongly related to leafhopper abundance (ß = 0.66).

In addition to these direct effects, the treatments had indirect effects on leafhopper abundance and species richness. This was mediated by changes in plant productivity and plant diversity. An inspection of the standardized total effects (Table 3) shows that FG manipulation was the most important driver of leafhopper species richness (total effect ϵ = 0.379). In addition, the latent variable “plant diversity” positively affected leafhopper species richness (ϵ = 0.249).

Table 3 Standardized total effects from structural equation model (ϵ), combining indirect and direct effects
Table 4 Effects of experimental treatments on leafhopper Shannon diversity, abundance and species richness (mean ±1SE)


The combination of management measures such as cutting and fertilizer application with the manipulation of functional group composition and their interactions led to contrasting sward types, ranging from nutrient poor, forb-dominated plots harbouring a greater diversity of plants to highly productive, graminoid-dominated plots with lower plant diversity. This had both direct and indirect effects on leafhopper species richness and abundance. Leafhopper species richness profited directly from both a higher cover of graminoids (due to herbicide-induced reduction of forbs) and from a lower cutting frequency. However, an indirect (negative) effect on leafhopper species richness was caused by the higher proportion of graminoids due to forb reduction and management intensification, which both had a negative impact on plant diversity.

Plant diversity had a negative effect on leafhopper abundance in our study, likely because the most abundant leafhopper species such as Arthaldeus pascuellus were grass feeders[22]. Plant diversity was driven by forbs, whereas productivity was driven by graminoids, which benefit from fertilizer application. Therefore, fertilizer application led to lower plant diversity but higher amounts of harvestable aboveground biomass.

The grassland in the experimental field site has been used for cattle grazing and hay making for at least a century. We therefore also expect the leafhopper community to be adapted to the long-term managed, grass-dominated vegetation, which may have led to the selection of the dominant species in the pool associated with grasses[22].

The increase in plant Shannon diversity (e) with graminoid reduction, the decrease in plant Shannon diversity with forb reduction, as well as the strong decrease in plant Shannon diversity with increasing proportion of graminoids highlights the impact of the functional group manipulation on plant diversity and composition and the strong contribution of forbs to plant species richness. Some graminoid species, such as Agrostis capillaris, Festuca rubra and Dactylis glomerata, play a dominant role in productivity in our study, resulting in a negative species richness–biomass relationship[23].

Plant diversity was particularly high under a higher frequency of three cuts per year, confirming that cutting can increase plant species diversity due to removal of nutrients from the soil[24]. However, a moderate frequency of two cuts per year may improve plant species richness in our experimental field site even more, because the disturbance rate is lower but removal of nutrients is still high[25]. This nutrient limitation only occurs when no fertilizer is applied and the harvested AGB is removed. Recent studies have shown that such grasslands may harbour a higher species number and proportion of forbs[26]. This is in line with our finding of a higher number of plant species (mainly driven by forbs) under the regime of three cuts per year.

When fertilizer is applied, many herbs cannot efficiently use higher nitrogen inputs[25] and are out-competed by more competitive species[27], which are mainly graminoids[28]. However, disturbance events such as mowing create niches for weaker competitors[29].

Similar to Morris[30], we observed a negative direct effect of higher cutting frequency on leafhopper species richness, but a positive effect on leafhopper abundance. However, we did not observe any significant direct responses of leafhopper species richness, abundance or Shannon diversity (e) to fertilizer application, which are often described in the literature (e.g.[31]).

As indicated by our structural equation model, fertilizer application had a weak indirect negative effect on leafhoppers, mediated by the strong increase in aboveground primary production and its negative association with plant species diversity.

NPK fertilizer application has been shown to result in a strong increase in aboveground biomass production[32] along with an increase in the proportion of graminoids[29] and a decrease in plant species diversity[33]. In our study, this may have been due to the higher tolerance of the dominant graminoid species (Dactylis glomerata, Festuca rubra) to high cutting frequencies and their faster re-growth capacity after cutting, especially under fertilizer application[34].

Many generalist leafhoppers can benefit from the improved performance of some graminoids under more frequent cutting, as the opposing effect of cutting frequency on leafhopper abundance (positive) and leafhopper species richness (negative) in the path diagram shows. Some very abundant species (e.g. Arthaldeus pascuellus, Streptanus sordidus, Deltocephalus pulicaris), which are generalists on graminoids[22], show a clear preference for plots cut three times per year.

The negative direct effect of higher cutting frequency on Shannon diversity (e) and species richness of leafhoppers indicates that the majority of leafhopper species found within our study site showed a clear preference for plots cut only once a year, which is in line with findings of other studies[31, 35]. This preference of leafhoppers for plots cut only once a year[31, 35] is also supported by our finding of a negative impact of higher cutting frequency on vegetation height and biomass. This shows that reduced disturbance rate results in larger amounts of food resources and shelter due to higher vegetation and aboveground biomass. Higher vegetation and aboveground biomass were shown by Kőrösi et al.[36] to have a positive effect on leafhopper abundance and species richness. We expect the same to apply here, since lower vegetation height and biomass on plots cut three times per year also implies a lower amount of available food resources and therefore directly results in lower leafhopper species richness.

The strong positive direct effect of the herbicide-induced increase of graminoids on leafhopper species richness can be explained by the preference for graminoids displayed by most leafhopper species found in our study (see Additional file4: Figure S3 and Additional file1: Table S4) which is in accordance with other studies[37].

Higher plant diversity has also been shown to lead to more diverse leafhopper communities[19], which, in our study, is owed to less common but more specialized species (e.g. Acanthodelphax spinosa, Cicadula persimilis, Conomelus anceps).

The herbicide treatment was applied only once in 2008. Nevertheless, this manipulation of plant functional group composition was efficient, since plant diversity and functional group composition were clearly affected even three years after the treatment.

Although most leafhopper species are mobile[18, 38] and our plots were small, the differences in management and plant functional group composition caused species sorting even on a small scale, which has also been found to be the case for highly mobile bees and wasps[39].

Our experimental manipulation of established grassland shows very strong design effects on vegetation, therefore we expect the direct and indirect effects of management regime and vegetation parameters on leafhoppers to be even stronger in larger and unconnected areas[40]. Our findings also indicate that leafhoppers prefer certain microhabitats within a defined community.

To gain a deeper understanding of the interactions behind the design effects (cutting frequency, fertilizer application and FG manipulation), we included vegetation parameters, such as vegetation height, harvested peak biomass and plant species richness (for forbs and graminoids separately) in our structural equation model.

The strong direct negative effect of higher cutting frequency and the positive direct effect of higher proportion of graminoids (due to FG manipulation) on leafhopper species richness (and also Shannon diversity, e) may be explained by several causes: (i) Plant diversity was greater with herbicide-induced reduction of graminoids, higher cutting frequency and without application of NPK-fertilizer, since it was mainly driven by forbs. But species richness of graminoids, which are the preferred food source for the majority of the leafhopper species caught[22] did not increase. Therefore, a higher proportion of graminoids indirectly reduced leafhopper species richness via its negative effect on plant diversity, but also had a direct positive effect due to greater availability of the preferred food source.

(ii) Cutting is a disturbance event and reduces food resource heterogeneity. The higher its frequency, the higher the direct impact on leafhoppers, since many leafhoppers are killed and removed during the process of harvesting [11]. This is comparable to the negative effects of cutting on slug abundance within the same experiment [10].

(iii) The lower vegetation height and biomass on plots cut three times per year also implies a lower amount of available food resources and therefore results in lower leafhopper species richness.

The proportion of graminoids increased with fertilizer application and also with application of forb-specific herbicides. This increased proportion of graminoids may explain the negative effect of vegetation biomass and height on plant diversity, as indicated by our structural equation model, demonstrating how a continuously managed agricultural system results in a high-yield but low-diversity system. This can be associated with less resilience towards sudden perturbation, microclimatic change and invasions[29].


We have shown that contrasting sward types may result from direct and indirect interactions among management regimes (cutting frequency and fertilizer application) and the manipulation of functional group composition.

Leafhopper species richness increased indirectly in nutrient-poor, forb-dominated plots, since these harboured a greater diversity of plants and therefore greater food-resource heterogeneity. On the other hand, leafhopper species richness increased directly with higher graminoid cover in highly productive plots with lower plant diversity, likely due to a greater availability of graminoids as a preferred food source for many leafhopper species.

With our novel approach of combining variation in cutting and fertilization with a manipulation of plant functional group composition (independent of management effects), we were able to disentangle the complexity of direct and indirect effects within a 100-year old, moderately species-rich and continuously managed grassland, based on the example of the easily accessible group of leafhoppers.

This study allows insights into the effects of grassland management on plant diversity, productivity and functional group composition and the trophic links and feedbacks on leafhopper species richness and abundance. Furthermore, we show that - apart from food resource heterogeneity represented by plant diversity - the availability of preferred food sources and the frequency of disturbance are important drivers of leafhopper species richness, but interactions between biodiversity and management are highly complex.


Description of study site

This study was performed in 2010 as part of the “GrassMan”-Experiment[13] near Neuhaus (Solling) in the Solling Mountains in Northern Germany (51°440 N, 9°320 E, 490 m a.s.l.).

Prior to the start of the experiment the study site was a nutrient poor, moderately wet Lolio-Cynosuretum grassland with high abundances of Agrostis capillaris (L.), Festuca rubra (L.), Rumex acetosa (L.), Veronica chamaedrys (L.) and Ranunculus repens (L.)[13].

Mean annual precipitation is 1028 mm and mean annual temperature is 6.9°C (Deutscher Wetterdienst, 1961–1990, station Holzminden-Silberborn, 440 m a.s.l.). In 2010, the year of the study, mean annual temperature was 8.0°C and annual precipitation was 1110 mm. The dominant soil type in the experimental area is a shallow (40–60 cm), stony Haplic Cambisol, developed on sediments of loess on the Middle Bunter (Triassic sandstone) formation with a loamy silt texture[41].

Study design

The experiment was established in 2008 in a permanent, formerly extensively used, cattle-grazed grassland. It was laid out as a three-factorial Latin rectangle[42] with the following factors (Additional file5: Figure S1): (i) plant functional group manipulation (three levels) using herbicides, (ii) fertilizer application (two levels) and (iii) cutting frequency (two levels), resulting in twelve treatment combinations.

To manipulate plant functional group presence, we applied (i) a combination of the forb-specific herbicides Fluroxypyr (Starane; Dow AgroSciences, Munich, Germany; 3 L ha-1) and Mecoprop-P (Duplosan; KV, Du Pont de Nemours, Neu-Isenburg, Germany; 3 L ha-1) or (ii) the graminoid-specific herbicide Select 240 EC (Stähler Int., Stade, Germany; 0.5 L ha-1), resulting in three levels of plant diversity: (i) forb-reduced (=graminoid-rich), (ii) graminoid-reduced (=forb-rich) and (iii) control. Herbicides were applied once in June 2008 (a “pulse” experiment sensu Bender et al.[43]).

In 2009 and 2010, plots were fertilized with N (Calcium ammonium nitrate N27: 13.5% NH4-N, 13.5% NO3-N, 4% MgO, 6% Ca) at two equal doses (2 × 90 kg ha-1) in April and May/June; in addition, fertilized plots received 30 kg P ha-1 and 100 kg K ha-1 in early June (Thomaskali®, 8% P2O5, 15% K2O, 20% CaO).

Control plots were not fertilized. Plots were cut either once (in July) or three times a year (May, July, September) using a Haldrup® forage combine harvester (INOTEC Engineering GmbH, Ilshofen, Germany) at a cutting height of 7 cm. The resulting twelve treatment combinations (equalling one block of the Latin rectangle; see Everwand et al.[10], Figure 4 and Additional file5: Figure S1) were arranged randomly and replicated six times, resulting in 72 plots. Each plot was a 15 × 15 m square surrounded by at least 3 m of frequently cut grass between plots, and 5 m between blocks.

Figure 4
figure 4

Aerial photograph of the GrassMan experiment. Photograph taken on 13th June 2009 by Laura Rose. The horizon was levelled using the ruler tool in Adobe Photoshop CS6.

Plant functional groups were not entirely removed, but target plant species were strongly reduced in abundance. Plant functional groups slowly recovered following herbicide application, but all FG manipulation treatments significantly affected vegetation parameters, such as compressed vegetation height, harvested biomass, functional group composition and plant species richness. More details on the experimental design, setup and treatment effects on vegetation can be found in Petersen et al.[13, 44] and Rose et al.[23, 45, 46].

Leafhopper sampling

Leafhoppers were sampled using two methods: i) by sweep netting (Heavy Duty Sweep Net, 7215HS, BioQuip, diameter: 38 cm), while walking a circular transect with a diameter of 8 m around the centre of each plot (30 sweeps each), in dry weather on two occasions (at the beginning of July and at the end of August 2010). Transects length was approximately 20 m, and there was a distance of at least 4 m to the edge of each plot. In addition, ii) we sampled leafhoppers by placing two transparent pan traps, containing an ethylene glycol/water mixture (1:3), 1 m apart, near to the centre of each plot. Pan traps were about 5 cm above vegetation height and were active for one week in five time intervals in 2010 (end of June, mid-July, early August, mid-August, end of August).

The specimens caught with both methods were transferred into ethanol (70% vol.) separately and identified to species level in the laboratory using Biedermann & Niedringhaus[47] and Kunz et al.[48]. One species with woody host plants was excluded, as we assumed that it had been swept off its host tree by wind and was not a true member of the grassland fauna. Species whose larvae used herbs or grasses as host plants and whose imagines fed on trees were also included in the analysis.

For female specimens of several genera, identification to species level is not possible (e.g. Psammotettix)[47, 48]. Thus, if male specimens were present, female specimens were assumed to belong to the same species. If not, they were only identified to genus level. If males of more than one species of a genus were present, the proportion of females was assumed to mirror that of males.

We found no interaction effects of the two sampling methods with the management variables (cutting frequency, fertilizer application) on leafhopper species richness (see e.g. Figure 4, Additional file6: Table S1 and Additional file7: Figure S2). In addition, vegetation height (a proxy for vegetation density) did not affect the number of leafhoppers caught by sweep netting.

We therefore pooled the data of both methods, which allowed us to cover the growing season of 2010 from early May until late September. For all diversity assessments, we used species richness, Shannon’s diversity index (H’) or its numbers equivalent exp(H’)[49].

Assessment of vegetation parameters

Because our treatments are likely to have affected plant productivity and vegetation structure, possibly indirectly affecting leafhopper species richness, we additionally measured a series of vegetation parameters:

(i) We conducted vegetation surveys on two quadrates, each of 1 m2 size, twice (in May before the first harvest and also again in August) on each plot. We recorded the percentage of cover, proportional yield of each species [50], plant species richness, functional group composition and presence-absence data of the functional groups (graminoids and forbs).

(ii) Plant aboveground biomass (AGB) was estimated as follows: First, during harvest, fresh weight of two 1.50 × 15 m strips per plot was measured using the harvester’s built-in scale.

To determine the water content of this sample, we took four subsamples that were homogenized, weighed, dried for 48 h at 65°C and subsequently weighed again. We then multiplied fresh weight by water content to obtain the total aboveground dry biomass (t ha-1) for every plot.

(iii) Proportions of graminoids and forbs (%) were determined from the vegetation surveys (derived as described above). Harvest was performed on all plots once a year in the end of June and additionally in mid-May and mid-September for the 3-cut treatment [13].

(iv)Compressed sward height (cm) was measured using a rising plate meter according to Castle [51] and the average value of 25 measures per plot was calculated. This was performed every 2-3 weeks, resulting in eleven time points throughout the growing season of 2010.

Statistical analysis

Data were analysed using the statistical software package R (version 2.15.2)[52]. In addition, we performed structural equation modelling using AMOS 20.0 (SPSS, Inc.). Treatment effects on vegetation and leafhoppers were assessed using generalized linear models (GLMs;[53]).

Models contained row- and column effects (fitted as factors, column was nested within block), sward composition (factor with three levels), cutting frequency and nutrient input (two levels each) with up to two-way-interactions. For abundance data (e.g. Table 4) we used quasipoisson GLMs, for proportion data LMs with a logit link[54, 55] and for exp(H’) we used GLMs with Gamma errors and an inverse link. Corresponding alternative models (e.g. quasipoisson or Gamma with log link) had higher residual deviance and were therefore not considered.

Continuous response variables (e.g. biomass or vegetation height) were log-transformed and analyzed using GLMs with an identity link. For each response variable in turn, maximal models containing all possible terms were manually simplified into models containing fewer explanatory variables. We compared the resulting nested models using F-tests (and Chi2 for quasipoisson models), until a minimal adequate model that only contained significant effects was obtained.

Significance of terms was assessed in two ways: (i) each parameter estimate from linear models was compared to zero using marginal t-tests; and (ii) terms in the models were additionally tested by sequential addition to a null model (sequential analysis of deviance tables; Additional file8: Table S3).

In addition to traditional GLM-based analyses, we employed structural equation models (SEMs), allowing us to test more complex hypotheses on indirect effects of treatments, plant productivity and plant diversity on leafhoppers[5658]. SEMs are particularly well suited for experimental contexts, i.e. where some variables are deliberately manipulated experimentally[59]. Furthermore, SEMs “can be used to develop accurate and meaningful final multiple regression models when collinearities among explanatory variables are thought to be present”[60], which was clearly the case for the vegetation properties measured here.

SEMs contained all three treatment variables, as well as latent variables[56] for plant productivity and plant diversity. For the SEMs we specified our design variables as numeric variables as follows:

Fertilizer treatment: no fertilizer = 0; NPK-fertilizer application = 1

Cutting frequency: one cut/year = 0; three cuts/year = 1

FG manipulation: FG graminoid-reduced = -1; FG control = 0; FG forb-reduced = 1

The sorting of FG manipulation was according to its effect on plant diversity and proportion of graminoids (see Figure 1). Plant productivity had two indicator variables: harvested aboveground biomass in July (AGB, t ha-1), and average compressed sward height. Plant diversity had the indicator variables “forbs” and “graminoids”; since only four legumes species (Lotus corniculatus, numeric variables as followsL. pedunculatus , Trifolium repens, Lathyrus pratensis) were present in a very low cover on 61 plots only, and none of the leafhopper species found had been categorized as preferentially feeding on legumes, we did not take legumes into account separately for the SEMs. Leafhopper abundance and species richness were taken separately (instead of (eH`) Shannon diversity) for the SEM to identify effects of design variables and vegetation parameters on leafhoppers.


  1. Food and agricultural organization statistics database. []

  2. WallisDeVries MF, Poschlod P, Willems JH: Challenges for the conservation of calcareous grasslands in northwestern Europe: integrating the requirements of flora and fauna. Biol Conserv. 2002, 104 (3): 265-273. 10.1016/S0006-3207(01)00191-4.

    Article  Google Scholar 

  3. Isselstein J, Jeangros B, Pavlu V: Agronomic aspects of biodiversity targeted management of temperate grasslands in Europe–a review. Agron Res. 2005, 3 (2): 139-151.

    Google Scholar 

  4. Sala OE, Chapin FS, Armesto JJ, Berlow E, Bloomfield J, Dirzo R, Huber-Sanwald E, Huenneke LF, Jackson RB, Kinzig A, et al: Biodiversity - global biodiversity scenarios for the year 2100. Science. 2000, 287 (5459): 1770-1774. 10.1126/science.287.5459.1770.

    Article  CAS  PubMed  Google Scholar 

  5. Allan E, Bossdorf O, Dormann CF, Prati D, Gossner MM, Tscharntke T, Blüthgen N, Bellach M, Birkhofer K, Boch S, et al: Interannual variation in land-use intensity enhances grassland multidiversity. Proc Natl Acad Sci. 2014, 111 (1): 308-313. 10.1073/pnas.1312213111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Galloway JN, Dentener FJ, Capone DG, Boyer EW, Howarth RW, Seitzinger SP, Asner GP, Cleveland CC, Green PA, Holland EA, et al: Nitrogen cycles: past, present, and future. Biogeochemistry. 2004, 70 (2): 153-226. 10.1007/s10533-004-0370-0.

    Article  CAS  Google Scholar 

  7. Smart SM, Thompson K, Marrs RH, Le Duc MG, Maskell LC, Firbank LG: Biotic homogenization and changes in species diversity across human-modified ecosystems. Proc R Soc B Biol Sci. 2006, 273 (1601): 2659-2665. 10.1098/rspb.2006.3630.

    Article  Google Scholar 

  8. Wesche K, Krause B, Culmsee H, Leuschner C: Fifty years of change in central European grassland vegetation: large losses in species richness and animal-pollinated plants. Biol Conserv. 2012, 150 (1): 76-85. 10.1016/j.biocon.2012.02.015.

    Article  Google Scholar 

  9. Attwood S, Maron M, House A, Zammit C: Do arthropod assemblages display globally consistent responses to intensified agricultural land use and management?. Global Ecol Biogeogr. 2008, 17 (5): 585-599. 10.1111/j.1466-8238.2008.00399.x.

    Article  Google Scholar 

  10. Everwand G, Scherber C, Tscharntke T: Slug responses to grassland cutting and fertilizer application under plant functional group removal. Acta Oecol. 2013, 48: 62-68.

    Article  Google Scholar 

  11. Humbert JY, Ghazoul J, Sauter GJ, Walter T: Impact of different meadow mowing techniques on field invertebrates. J Appl Entomol. 2010, 134 (7): 592-599.

    Google Scholar 

  12. Humbert J-Y, Ghazoul J, Walter T: Meadow harvesting techniques and their impacts on field fauna. Agric Ecosyst Environ. 2009, 130 (1–2): 1-8.

    Article  Google Scholar 

  13. Petersen U, Wrage N, Köhler L, Leuschner C, Isselstein J: Manipulating the species composition of permanent grasslands—a new approach to biodiversity experiments. Basic Appl Ecol. 2012, 13 (1): 1-9. 10.1016/j.baae.2011.10.003.

    Article  Google Scholar 

  14. Scherber C, Mwangi PN, Temperton VM, Roscher C, Schumacher J, Schmid B, Weisser WW: Effects of plant diversity on invertebrate herbivory in experimental grassland. Oecologia. 2006, 147 (3): 489-500. 10.1007/s00442-005-0281-3.

    Article  PubMed  Google Scholar 

  15. Scherber C, Eisenhauer N, Weisser WW, Schmid B, Voigt W, Fischer M, Schulze E-D, Roscher C, Weigelt A, Allan E, et al: Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. Nature. 2010, 468 (7323): 553-556. 10.1038/nature09492.

    Article  CAS  PubMed  Google Scholar 

  16. Hudewenz A, Klein A-M, Scherber C, Stanke L, Tscharntke T, Vogel A, Weigelt A, Weisser WW, Ebeling A: Herbivore and pollinator responses to grassland management intensity along experimental changes in plant species richness. Biol Conserv. 2012, 150 (1): 42-52. 10.1016/j.biocon.2012.02.024.

    Article  Google Scholar 

  17. Diaz S, Symstad AJ, Chapin FS, Wardle DA, Huenneke LF: Functional diversity revealed by removal experiments. Trends Ecol Evol. 2003, 18 (3): 140-146. 10.1016/S0169-5347(03)00007-7.

    Article  Google Scholar 

  18. Nickel H: The leafhoppers and planthoppers of Germany (Hemiptera, Auchenorrhyncha): patterns and strategies in a highly diverse group of phytophagous insects. 2003, Keltern: Goecke & Evers

    Google Scholar 

  19. Biedermann R, Achtziger R, Nickel H, Stewart A: Conservation of grassland leafhoppers: a brief review. J Insect Conserv. 2005, 9 (4): 229-243. 10.1007/s10841-005-0531-z.

    Article  Google Scholar 

  20. Rösch V, Tscharntke T, Scherber C, Batáry P, Osborne J: Landscape composition, connectivity and fragment size drive effects of grassland fragmentation on insect communities. J Appl Ecol. 2013, 50 (2): 387-394. 10.1111/1365-2664.12056.

    Article  Google Scholar 

  21. Hollier J, Maczey N, Masters G, Mortimer S: Grassland leafhoppers (Hemiptera: Auchenorrhyncha) as indicators of habitat condition–a comparison of between-site and between-year differences in assemblage composition. J Insect Conserv. 2005, 9 (4): 299-307. 10.1007/s10841-005-8821-z.

    Article  Google Scholar 

  22. Nickel H, Remane R: Artenliste der Zikaden Deutschlands, mit Angabe von Nährpflanzen, Nahrungsbreite, Lebenszyklus, Areal und Gefährdung (Hemiptera, Fulgoromorpha et Cicadomorpha). Beiträge zur Zikadenkunde. 2002, 5: 27-64.

    Google Scholar 

  23. Rose L, Leuschner C: The diversity–productivity relationship in a permanent temperate grassland: negative diversity effect, dominant influence of management regime. Plant Ecol Diver. 2012, 5 (3): 265-274. 10.1080/17550874.2012.723763.

    Article  Google Scholar 

  24. Oomes M, Mooi H: The effect of cutting and fertilizing on the floristic composition and production of an Arrhenatherion elatioris grassland. Vegetatio. 1981, 46 (1): 233-239.

    Article  Google Scholar 

  25. Zechmeister HG, Schmitzberger I, Steurer B, Peterseil J, Wrbka T: The influence of land-use practices and economics on plant species richness in meadows. Biol Conserv. 2003, 114 (2): 165-177. 10.1016/S0006-3207(03)00020-X.

    Article  Google Scholar 

  26. Fischer M, Wipf S: Effect of low-intensity grazing on the species-rich vegetation of traditionally mown subalpine meadows. Biol Conserv. 2002, 104 (1): 1-11. 10.1016/S0006-3207(01)00149-5.

    Article  Google Scholar 

  27. Chapin FS: The mineral-nutrition of wild plants. Annu Rev Ecol Syst. 1980, 11: 233-260. 10.1146/

    Article  CAS  Google Scholar 

  28. Pearce IS, van der Wal R: Effects of nitrogen deposition on growth and survival of montane Racomitrium lanuginosum heath. Biol Conserv. 2002, 104 (1): 83-89. 10.1016/S0006-3207(01)00156-2.

    Article  Google Scholar 

  29. MacDougall AS, McCann KS, Gellner G, Turkington R: Diversity loss with persistent human disturbance increases vulnerability to ecosystem collapse. Nature. 2013, 494 (7435): 86-89. 10.1038/nature11869.

    Article  CAS  PubMed  Google Scholar 

  30. Morris MG: Responses of Grassland Invertebrates to Management by Cutting. IV. Positive Responses of Auchenorhyncha. J Appl Ecol. 1981, 18 (3): 763-771. 10.2307/2402368.

    Article  Google Scholar 

  31. Nickel H, Achtziger R: Wiesen bewohnende Zikaden (Auchenorrhynca) im Gradienten von Nutzungsintensität und Feuchte. Beiträge zur Zikadenkunde. 1999, 3: 65-80.

    Google Scholar 

  32. Harpole WS, Tilman D: Grassland species loss resulting from reduced niche dimension. Nature. 2007, 446 (7137): 791-793. 10.1038/nature05684.

    Article  CAS  PubMed  Google Scholar 

  33. Pan Q, Bai Y, Wu J, Han X: Hierarchical plant responses and diversity loss after nitrogen addition: testing three functionally-based hypotheses in the Inner Mongolia grassland. Plos One. 2011, 6 (5): e20078-10.1371/journal.pone.0020078.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Leuschner C, Gebel S, Rose L: Root trait responses of six temperate grassland species to intensive mowing and NPK fertilisation: a field study in a temperate grassland. Plant Soil. 2013, 373: 687-698. 10.1007/s11104-013-1836-4.

    Article  CAS  Google Scholar 

  35. Morris MG, Lakhani KH: Responses of grassland invertebrates to management by cutting.1. Species-Diversity of Hemiptera. J Appl Ecol. 1979, 16 (1): 77-98. 10.2307/2402730.

    Article  Google Scholar 

  36. Koroesi A, Batary P, Orosz A, Redei D, Baldi A: Effects of grazing, vegetation structure and landscape complexity on grassland leafhoppers (Hemiptera: Auchenorrhyncha) and true bugs (Hemiptera: Heteroptera) in Hungary. Insect Conserv Diver. 2011, 5 (1): 57-66.

    Article  Google Scholar 

  37. Nickel H, Hildebrandt J: Auchenorrhyncha communities as indicators of disturbance in grasslands (Insecta, Hemiptera)— a case study from the Elbe flood plains (Northern Germany). Environment. 2003, 98: 183-199.

    Google Scholar 

  38. Waloff N: Dispersal by flight of leafhoppers (Auchenorrhyncha: Homoptera). J Appl Ecol. 1973, 10: 705-730. 10.2307/2401864.

    Article  Google Scholar 

  39. Ebeling A, Klein A-M, Weisser WW, Tscharntke T: Multitrophic effects of experimental changes in plant diversity on cavity-nesting bees, wasps, and their parasitoids. Oecologia. 2012, 169 (2): 453-465. 10.1007/s00442-011-2205-8.

    Article  PubMed  Google Scholar 

  40. Rosch V, Tscharntke T, Scherber C, Batary P: Landscape composition, connectivity and fragment size drive effects of grassland fragmentation on insect communities. J Appl Ecol. 2013, 50 (2): 387-394. 10.1111/1365-2664.12056.

    Article  Google Scholar 

  41. Keuter A, Hoeft I, Veldkamp E, Corre M: Nitrogen response efficiency of a managed and phytodiverse temperate grassland. Plant Soil. 2013, 364: 193-206. 10.1007/s11104-012-1344-y.

    Article  CAS  Google Scholar 

  42. Clewer AG, Scarisbrick DH: Practical statistics and experimental design for plant and crop science. 2001, Chichester: John Wiley & Sons Ltd

    Google Scholar 

  43. Bender EA, Case TJ, Gilpin ME: Perturbation experiments in community ecology: theory and practice. Ecology. 1984, 65 (1): 1-13. 10.2307/1939452.

    Article  Google Scholar 

  44. Petersen U, Wrage-Mönnig N, Isselstein J: Effects of herbicide application to control sward composition in different management variants. International Journal of Biodiversity Science, Ecosystem Services & Management. 2013, 9 (2): 1-11.

    Article  Google Scholar 

  45. Rose L, Coners H, Leuschner C: Effects of fertilization and cutting frequency on the water balance of a temperate grassland. Ecohydrology. 2012, 5 (1): 64-72. 10.1002/eco.201.

    Article  Google Scholar 

  46. Rose L, Rubarth MC, Hertel D, Leuschner C: Management alters interspecific leaf trait relationships and trait-based species rankings in permanent meadows. J Veg Sci. 2012, doi:10.1111/j.1654-1103.2012.01455.x

    Google Scholar 

  47. Biedermann R, Niedringhaus R: Die Zikaden Deutschlands: Bestimmungstafeln für alle Arten. 2004, Scheeßel: Wissenschaftlicher Akademischer Buchvertrieb Fründ

    Google Scholar 

  48. Kunz G, Nickel H, Niedringhaus R: Fotoatlas der Zikaden Deutschlands – photographic atlas of the planthoppers and leafhoppers of Germany. 2011, Scheeßel: Wissenschaftlicher Akademischer Buchvertrieb Fründ

    Google Scholar 

  49. Jost L: Entropy and diversity. Oikos. 2006, 113 (2): 363-375. 10.1111/j.2006.0030-1299.14714.x.

    Article  Google Scholar 

  50. Klapp E, Stählin A: Standorte, Pflanzengesellschaften und Leistung des Grünlandes. 1936, Stuttgart: Ulmer

    Google Scholar 

  51. Castle ME: Simple disk instrument for estimating herbage yield. Journal of the British Grassland Society. 1976, 31 (1): 37-40. 10.1111/j.1365-2494.1976.tb01113.x.

    Article  Google Scholar 

  52. R Core Team: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, []

    Google Scholar 

  53. Crawley MJ: The R book. 2007, Southern Gate, Chichester, West Sussex: John Wiley & Sons Ltd.

    Book  Google Scholar 

  54. Warton DI, Hui FKC: The arcsine is asinine: the analysis of proportions in ecology. Ecology. 2010, 92 (1): 3-10.

    Article  Google Scholar 

  55. Scherber C, Gladbach DJ, Stevnbak K, Karsten RJ, Schmidt IK, Michelsen A, Albert KR, Larsen KS, Mikkelsen TN, Beier C, et al: Multi-factor climate change effects on insect herbivore performance. Ecology and Evolution. 2013, 3 (6): 1449-1460. 10.1002/ece3.564.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Bollen KA: Structural equations with latent variables. 1989

    Book  Google Scholar 

  57. Shipley B: Cause and correlation in biology. 2000, Cambridge University Press

    Book  Google Scholar 

  58. Scherber C, Mwangi PN, Schmitz M, Scherer-Lorenzen M, Bessler H, Engels C, Eisenhauer N, Migunova VD, Scheu S, Weisser WW, et al: Biodiversity and belowground interactions mediate community invasion resistance against a tall herb invader. J Plant Ecol-Uk. 2010, 3 (2): 99-108. 10.1093/jpe/rtq003.

    Article  Google Scholar 

  59. Grace JB: Structural equation modeling and natural systems. 2006, Cambridge University Press

    Book  Google Scholar 

  60. Graham MH: Confronting multicollinearity in ecological multiple regression. Ecology. 2003, 84 (11): 2809-2815. 10.1890/02-3114.

    Article  Google Scholar 

Download references


We thank Tatiana From and Andreas Keuter for providing data on vegetation, biomass and soil properties, Lars Köhler and Johannes Isselstein for coordination of the GrassMan project as well as David Perovic for proof reading. Furthermore we thank the Ministry of Science and Culture of Lower Saxony and the “Niedersächsisches Vorab” for financial support.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christoph Scherber.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GE, TT and CS designed the study. GE sampled leafhoppers in the field. GE and VR determined leafhopper species. GE and CS did all statistical analyses. All authors contributed to manuscript drafting, writing and revisions. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S4: Species list of the Auchenorrhyncha within the GrassMan experiment with food preferences according to Nickel & Remane (2002) and their abundance. (PDF 83 KB)

Additional file 2: Table S2: ANOVA-table: Leafhopper species richness vs. Design. (PDF 29 KB)


Additional file 3: Table S5: Mean and standard error of leafhopper responses to experimental design (treatments). (PDF 12 KB)


Additional file 4: Figure S3: Box plot showing the effects of herbicide application and cutting frequency on Shannon diversity of grass specialist leafhopper. (PDF 29 KB)


Additional file 5: Figure S1: Experimental design of the Grassman Experiment, showing the Latin rectangle of 12 treatments in 6 replications. Gra- = graminoid reduced plots (=forb enhanced); Forb- =forb reduced (=graminoid enhanced); Con = Control (no herbicide application). The grey area around and between the plots is mown monthly. Plot size 15 m x 15 m, space between plots 3 m, between blocks 5 m (PDF 125 KB)


Additional file 6: Table S1: F and p-values of glm’s testing for significant effects of sampling method in combination with design treatments. (PDF 12 KB)


Additional file 7: Figure S2: Comparison of the two different sampling methods in combination with functional group manipulation and cutting frequency. (a) leafhopper species richness; (b) leafhopper Shannon diversity (eH´); (c) leafhopper abundance. (PDF 80 KB)

Additional file 8: Table S3: ANOVA-table: Leafhopper abundance vs. Design. (PDF 29 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Cite this article

Everwand, G., Rösch, V., Tscharntke, T. et al. Disentangling direct and indirect effects of experimental grassland management and plant functional-group manipulation on plant and leafhopper diversity. BMC Ecol 14, 1 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Auchenorrhyncha
  • Management intensity
  • Mowing
  • Plant species composition
  • Forb
  • Graminoid
  • Biodiversity experiment
  • Removal experiment