El abandono del pastoreo afecta negativamente a la calidad del pasto en pastizales atlánticos ibéricos ; Grazing abandonment negatively affects forage quality in iberian atlantic grasslands

Forage quality in mountain grasslands is an important factor for maintaining or improving animal performance. Previous studies identified plant traits beneficial for animal nutrition; however, little is known regarding their maintenance in plant communities. Our main hypothesis was that grazing favors plants with high regrowth capacity and low resource utilization efficiency, promoting the production of high-quality forage. To test the hypothesis, we conducted a grazer exclusion experiment in an Iberian Atlantic grassland combined with a response-effect framework based on various plant traits and an appropriate statistical methodology to measure the relative contribution of species turnover and intraspecific trait variability (ITV) after grazing abandonment. Our results showed that forage quality declined after short-term grazing abandonment, via strong ITV effects and weak species turnover effects. Strong species turnover effects might appear after longer period of abandonment if, as expected, tall grasses with low tissue quality outcompete species with high tissue quality.


Introduction
Cattle and small ruminants, such as sheep and goats, consume 80% of the total animal feed produced worldwide (Herrero et al., 2013). Grasslands occupy approximately 40% of the global terrestrial surface and provide many ecosystem services, including most of the biomass used by livestock (Lee et al., 2017;Moreno García et al., 2014). Therefore, forage quality (i.e., digestibility and nutritional value) in grassland systems influences the amount of consumed feed, efficiency of rumination, rate of animal weight gain, quality and volume of produced milk, and animal reproductive success (Herrero et al., 2015). The quality of forage grasses is improved by reducing the proportion of fiber and increasing the proportion of protein (Lee et al., 2017;Waghorn & Clark, 2004). Although much research has been carried out focusing on the introduction of traits beneficial for animal nutrition into plants, little is known regarding the maintenance of such traits in plant communities (Parsons et al., 2011).
Grazing by domestic and wild ungulates affects ecosystem functions related to forage quality by promoting some plant traits over others (Deléglise et al., 2015;Peco et al., 2017). Thus, a response-effect framework based on plant traits (Lavorel & Garnier, 2002) may be a useful tool to understand the response of specific plant traits to grazing as well as their effect on the nutritional quality of produced forage. Plants use two main strategies against herbivores: avoidance (or defense) and tolerance (Díaz et al., 2007;Van Der Meijden et al., 1988). The latter, which is common in grasslands with high resource availability and long evolutionary history of grazing (Cingolani et al., 2005), favors plant species with fast regrowth rate after defoliation (Peco et al., 2017) and produces high-quality forage (Díaz et al., 2007). Additionally, herbivores enhance nitrogen (N) inputs to the soil through their urine and feces, increasing the abundance of grass species with low fiber content and high N content (Pontes et al., 2010(Pontes et al., , 2007. In this type of grassland, grazing abandonment promotes plant species that more efficiently compete for light and utilize resources (Díaz et al., 2001;Mc-Naughton, 1984;Milchunas et al., 1988), producing forage of low nutritional value.
Previous studies on Iberian Atlantic grasslands have shown that a reduction in forage quality occurred after cessation of grazing; specifically, lower protein concentration and digestibility, and higher fiber content were present (Odriozola et al., 2014). In this study, we propose a hypothesis for the mechanism behind these changes in forage quality that is based on the tolerance strategy of the plants: grazing favors plants with high regrowth capacity and low resource utilization efficiency, promoting the production of forage with high nutritional value; conversely, grazing abandonment favors plants with the ability to compete for light and that have high resource utilization efficiency, promoting the production of low-quality forage. To test this hypothesis, grazing abandonment by mixed livestock (sheep, cattle, and horses) was experimentally simulated (five years of exclusion) using grazerexclusion fences in productive grasslands of the Atlantic Iberian Peninsula. Then, in line with the mass-ratio hypothesis, plant traits of the most abundant (structural) species (covering at least 80% of the area) were measured to assess the effect of grazing abandonment on forage nutritional quality. The mass-ratio hypothesis states that the traits of dominant species largely determine the effect of plant communities on ecosystem functions (Grime, 1998). On the basis of this hypothesis, previous studies have successfully used the community weighted mean (CWM) of structural species to assess the effects of landuse changes on ecosystem processes (e.g., Garnier et al., 2004Garnier et al., , 2007Quétier et al., 2007). Usually, functional trait studies consider plant species as functionally homogeneous entities with no phenotypic plasticity; thus, the effect of intraspecific trait variability (ITV) on response capacity (plasticity) is disregarded. Nevertheless, most quantitative plant traits are highly variable within species (Albert et al., 2010;Cornelissen et al., 2003;Lepš et al., 2011;Westoby et al., 2002). Moreover, Lee et al. (2017) reported that intraspecific physiological responses were an im-RESUMEN: La calidad del pasto en los pastizales de montaña es un factor importante para mantener o mejorar el rendimiento animal. Estudios previos han identificado atributos vegetales beneficiosos para la nutrición animal; sin embargo, poco se sabe sobre su mantenimiento en comunidades vegetales. Nuestra principal hipótesis sostiene que el pastoreo favorece a las especies con una alta capacidad de rebrote y una baja eficiencia de utilización de recursos, promoviendo así la producción de pasto de alta calidad. Para testar la hipótesis, realizamos un experimento de exclusión del pastoreo en un pastizal de la montaña atlántica, combinado con un marco teórico de respuesta-efecto basado en varios atributos vegetales y una metodología estadística apropiada para medir la contribución relativa del recambio de especies y la variabilidad intraespecífica de los atributos (ITV) tras el abandono del pastoreo. Nuestros resultados mostraron que la calidad del pasto disminuyó tras el abandono del pastoreo a corto plazo, principalmente debido al efecto de la ITV. Los efectos por recambio de especies podrían aparecer tras un período más largo de abandono si, tal y como es esperable, las gramíneas altas con tejidos de baja calidad nutritiva sustituyen a las especies de mayor calidad.
PALABRAS CLAVE: Respuesta de la comunidad; atributos de efecto, calidad del forraje, exclusión de herbívoros, variabilidad intraespecífica, recambio de especies. portant driver of forage quality worldwide. In order to include plasticity effects on community response to grazing abandonment and overcome the limitations of previous studies, plant traits were measured in two grassland sites in this study (Alotza & Uzkuiti) under grazing and nongrazing conditions. Additionally, the relative contribution of species turnover and ITV after grazing abandonment was measured using the statistical methodology developed by Lepš et al. (2011).
The objectives of this study were to test 1) whether grazing abandonment promotes the production of forage with low nutritional value; and 2) whether species turnover and ITV after grazing reinforce each other through selection of similar plant traits. Although previous studies have attempted to associate the response of plant traits to grazing (or cutting and nutrient addition) with their effects on herbivore nutrition (Cingolani et al., 2005;Cruz et al., 2010;Deléglise et al., 2015;Moreno García et al., 2014;Pontes et al., 2010Pontes et al., , 2007, to our knowledge, this is the first study to measure the relative contribution of species turnover and ITV after grazing abandonment using a specific statistical methodology.

Study area
Experiments were conducted in two semi-natural grasslands, Alotza (43°0′10.6′′ N, 2°5′22′′ W; 1,223 m elevation) and Uzkuiti (43°0′50′′ N, 2°4′3′′ W; 1,300 m elevation), located in the Aralar Natural Park, which is an 11,000 ha protected area located in the Atlantic Basque Country, Spain. The climate is oceanic with a mean annual temperature of 7°C and a mean annual precipitation of 1,330 mm. Aralar grasslands are used by mixed livestock in a moderate/high stocking rate (3.2 livestock units ha -1 d -1 ; 13% beef cattle, 52% dairy sheep, and 35% horses; Odriozola et al., 2014). Livestock herds are managed based on a short-distance transhumance system using lowland farms in winter/spring (December to April) and, more extensively, the upland grasslands in summer/fall (May to November). The park comprises productive native grasslands (mean aboveground net primary production of 1.97 ± 0.39 t dry mass ha −1 yr −1 ; unpublished data) with a calcareous substrate (Gibbons & Moreno, 2002). The grasslands (Jasiono laevis-Danthonietum decumbentis Loidi;Loidi, 1982) have a high conservation value (code 6230, 92/43/EEC;European Commission, 2013) and are associated with the production of a local cheese, Idiazabal, from the milk of the Latxa sheep breed (Batalla, 2015).

Experimental design
In May 2005, we installed a permanent 50 m × 50 m fence plot at each site (Alotza, Uzkuiti) to establish grazer exclusion plots (E plots). Attached to each E plot, we delineated a 50 m × 50 m plot (G plots), in which sheep, cattle, and horses were allowed to graze continuously during the vegetative period (May-November). Both sites are located in a flat terrain and have similar environmental characteristics and grazing pressures (Table 1). Species composition was not significantly different between E and G plots prior to fence establishment in May 2005 (Odriozola et al., 2017). Soil from each plot was analyzed in July 2010 to measure pH, Kjeldahl total N (mg L −1 ), available phosphorus (P; mg L −1 ), available potassium (K; mg L −1 ), calcium (Ca; mg L −1 ), and magnesium (Mg; mg L −1 ). Soil water content was measured at 15 cm soil depth and at 2 h intervals from June 2010 to October 2011, using Em50 data loggers connected to an ECH2O sensor system (Decagon Devices Inc., Pullman, WA, USA).

Plant sampling and traits
In June 2010, the floristic composition in each plot was determined in 15 randomly placed 1 × 1 m 2 quadrats. Species abundance was measured as canopy cover using the following index: presence (<1%) and a scale from 1 (1-10%) to 10 (91-100%). Plant species nomenclature followed that of the standard floras (Table A1).
The following traits were used to assess the ecological response to grazing (i.e., response traits): plant height (H, cm) and specific leaf area (SLA, m 2 kg −1 ), since they each represent two key dimensions of plant function (Díaz et al., 2007(Díaz et al., , 2016. H is related to successional status and response to grazing (Díaz et al., 2007), and SLA is related to relative growth rate and resource utilization efficiency (Westoby et al., 2002). The traits used to estimate forage nutritional value (i.e., effect traits) were leaf N content (LNC, mg g −1 ) and carbon (C) to N ratio (CN). LNC is positively correlated to forage nutritional value and livestock productivity, and CN is related to fiber content and is negatively correlated with livestock productivity (Lee et al., 2017). SLA, LNC, and CN were measured in seven randomly selected, fully grown, undamaged individuals from each species, whereas H was measured in 25 individuals from each species using a metallic ruler (Cornelissen et al., 2003). SLA was defined as the one-sided area of a fresh leaf divided by its oven-dry mass. Hydrated and expanded leaves were scanned using an Epson Expression 10000 XL scanner (Epson, Suwa, Japan), and the area of each image (containing all the leaves of an individual) was measured using Win Folia 2008 (Regent Instruments, Ontario, Canada). The mass of the same leaves was measured using a precision balance after oven drying at 60°C for at least 72 h (Cornelissen et al., 2003). After removing the petiole or rachis, powdered leaf samples were prepared and used to determine the C and N contents using combustion elemental analysis (AOAC, 1990). CN was obtained by dividing the C content by the N content per leaf dry mass.

Data analysis
Distance-based redundancy analysis (dbRDA) was applied to test the null hypothesis of no effect of grazing abandonment on the abundances of the five structural species (Legendre & Legendre, 2012). The response matrix was created by Hellinger-transformation of the fivespecies composition matrix and computing Euclidean distances (Legendre & Gallagher, 2001). Analysis of variance (ANOVA) was fitted to each of the five structural species and four traits (H, SLA, LNC, and CN) to test the null hypothesis of no effect of grazer exclusion in the ITV of measured traits. Tukey's post-hoc test was used to separate the means. In all analyses, site (Alotza & Uzkuiti), treatment (G and E), and their interactions were used as explanatory variables.
The relative contribution of species composition and ITV to community response under grazer exclusion conditions was determined as described by Lepš et al. (2011). CWMs calculated using fixed trait averages (i.e., fixed mean trait value of each species for all individuals measured in the four plots) were only affected by changes in community floristic composition (species turnover), whereas CWMs calculated using specific trait averages for each field plot (i.e., specific mean trait value of each species in each plot) were affected by both changes in composition and ITV (i.e., ITV effect = specific average -fixed average). Linear models were used to test the null hypothesis of no effect of grazer exclusion on each community variable. Spearman correlation was performed to test the relationship between fixed and specific CWMs of response (H and SLA) and effect (LNC and CN) traits.
Variations in each trait were partitioned into fixed effects (species turnover), ITV effects, and their covariation by using the sum of squares (SS) decomposition of ANO-VAs or linear models. Total variation in specific averages is the total variation in community trait averages (i.e., variation in fixed averages and ITV effects). Thus, if SS specific = SS fixed + SS ITV , then changes in composition and ITV operate independently; if SS specific > SS fixed + SS ITV , both effects are positively correlated; and if SS specific < SSfixed + SS ITV , both effects are negatively correlated (Lepš et al., 2011). The R code and data are provided in Supplementary Data B and C, respectively.

Effect of grazing abandonment on structural species composition
The permutation test of dbRDA revealed significant effects of grazing abandonment on the abundance of the five structural species (Table A2). T. repens was clearly associated with G plots, whereas G. saxatile was associated with E plots (Figure 1). A. capillaris tended to be associated with E plots, whereas L. campestris and F. microphylla were not affected by the plot type ( Figure 1).

Effect of grazing abandonment on ITV
The five structural species had high ITV in all plots regarding the studied traits (Tables 2, A3). In E plots, H and SLA showed high ITV since the individuals of the same species were considerably taller and had lower SLA (Table 2). In G plots, species also tended to have higher LNC and lower CN; however, the results were not consistent. For instance, G. saxatile had low LNC in the G plot of Uzkuiti, but high LNC in the G plot of Alotza (Table 2). Moreover, the LNC and CN of F. microphylla and T. repens were not affected by the plot type (Table A3c, 3d); in all plots, F. microphylla had low LNC and high CN, whereas T. repens had the highest LNC and the lowest CN (Table 2).

Relative contribution of species turnover and ITV to the response to grazing abandonment
Weak and inconsistent effects of grazing abandonment were observed when CWMs based on fixed trait averages were considered (Figure 2). H and CN tended to be higher, whereas LNC tended to be lower in E plots (Figure 2A, C, and D). However, SLA was reduced in the E plot of Alotza but increased in the E plot of Uzkuiti ( Figure 2B). However, strong and consistent effects of grazing abandonment were observed when specific CWMs were considered (Figure 2). H and CN were markedly higher, whereas SLA and LNC were markedly lower in E plots ( Figure 2). Additionally, strong correlations were observed between response (H and SLA) and effect (LNC and CN) traits using both fixed and specific CWMs (Table 3). Plots with high H had low SLA and LNC but high CN; accordingly, plots with low H had
high SLA and LNC but low CN (Table 3). Species turnover and ITV selected similar H values after grazing abandonment, i.e., turnover and ITV effects were positively correlated; however, the contribution of turnover was almost negligible ( Figure 3A). The difference in SLA between grazing regimes was mainly due to intraspecific variability; however, the interaction fraction (Site × Treatment bar in Figure 3B) showed negative covariance between turnover and ITV. This reflects the fact that while individuals of the same species presented lower SLA in the E plot of Uzkuiti, species with higher fixed SLA were more abundant in the same plot ( Figure 2B). This inconsistency between turnover and ITV was not observed in Alotza ( Figure 2B). For LNC and CN, turnover and ITV effects were important and consistent, i.e., both effects were positively correlated, and selected similar trait values (i.e., low LNC and high CN) after grazing abandonment ( Figures 3C-3D).

Discussion
Grasslands provide multiple ecosystem services, including forage supply for livestock (Herrero et al., 2015). However, information on the maintenance of plant traits that support forage nutritional value in grasslands is limited (Parsons et al., 2011). Consequently, decision making is conducted based on empirical observations of the performance of the grass crop, but the mechanisms behind this performance are rarely Table 3: Spearman correlation between fixed and specific community weighted means (CWMs) of measured traits in five structural species under grazing and grazer-exclusion conditions in two Iberian Atlantic grasslands (Alotza and Uzkuiti). A) Spearman correlation coefficients for fixed CWMs; B) Spearman correlation coefficients for specific CWMs. Tabla 3: Correlación Spearman entre las medias ponderadas a nivel de comunidad (CWMs) fijas y específicas de los atributos medidos en las cinco especies estructurales bajo condiciones de pastoreo y exclusión en dos pastos ibéricos atlánticos (Alotza y Uzkuiti). A) Coeficiente de correlación Spearman para CWMs fijas; B) Coeficiente de correlación Spearman para CWMs específicas.  Figure 3. Decomposition of specific community weighted mean (CWM) variability into species turnover and intraspecific variability effects. Dark grey, changes in species composition; pale grey, changes in intraspecific variability. Black horizontal bars denote total variation (variation in specific averages). Fractions correspond to mean squares taken from ANOVA tables in Table A4. Space between the top of the column and the bar corresponds to covariation. Bar above the column indicates positive covariation/positive correlation between species turnover and intraspecific variability; bar across a column indicates negative covariation/negative correlation between species turnover and intraspecific variability; bar at the top of a column indicates independent effects. Values were standardized by total variation in specific averages. Figura 3. Descomposición de la variabilidad de la media ponderada a nivel de comunidad (CWM) específica en los efectos debidos al recambio de especies y a la variabilidad intraespecífica. Gris oscuro, cambios en la composición de especies; gris claro, cambios en la variabilidad intraespecífica. Las barras negras horizontales indican la variación total (variación en medias específicas). Las fracciones corresponden a las medias cuadradas tomadas de las tablas de ANOVA de la Tabla A4. El espacio entre la parte superior de la columna y la barra corresponde a la covariación. La barra por encima de la columna indica una covariación positiva/ correlación positiva entre el recambio de especies y la variabilidad intraespecífica; la barra que atraviesa la columna indica una covariación negativa/correlación negativa entre el recambio de especies y la variabilidad intraespecífica; La barra a la par de la parte superior de la columna indica efectos independientes. Los valores han sido estandarizados en base a la variación total de las medias específicas.
considered (Pontes et al., 2015). In the present study, we aimed to investigate the underlying mechanisms of grassland community response to grazing using a framework based on plant traits (Cingolani et al., 2005;Cruz et al., 2010;Deléglise et al., 2015;Moreno García et al., 2014;Pontes et al., 2010Pontes et al., , 2007. To our knowledge, this is the first study of forage quality that has measured the relative contribution of species turnover and ITV after grazing abandonment by using a specific statistical methodology.

Community response to grazing abandonment and effects on forage quality
The present study showed that the plant community had low aggregated values of H and high aggregated values of SLA under grazing conditions; however, this trend was reversed after grazing abandonment. These results were in agreement with those reported in previous studies (e.g., Cingolani et al., 2005;Cruz et al., 2010;Evju et al., 2009;Moreno García et al., 2014;Peco et al., 2017), suggesting that tolerance to grazing is the main strategy against defoliation by herbivores in productive grasslands with a long history of grazing (Cingolani et al., 2005;Díaz et al., 2007). Response traits H and SLA were strongly associated with the effect traits LNC and CN, revealing that the shift in trait composition following the change in grassland strategy after grazing abandonment (from tolerance to grazing, to aboveground competition) explains the decline in forage quality reported in these grassland sites (Odriozola et al., 2014). Theoretically, grazing tolerance requires high SLA, which increases shoot regrowth capacity (Westoby, 1999), and high LNC, which improves leaf quality and selectivity by herbivores (Cingolani et al., 2005). However, light competition requires increased H to capture the light and a conservative resource utilization (i.e., low SLA) for a higher investment in structural components (i.e., high CN) (Cingolani et al., 2005).
The production of forage of high quality under grazing conditions has been also reported in other sub-humid and humid grasslands where tolerance is the plant community strategy against herbivores (Cingolani et al., 2005;Cruz et al., 2010;Deléglise et al., 2015). This has also been observed in the semi-arid grasslands of South Africa (Moreno García et al., 2014), although intense grazing areas coincided with water points in this case. By contrast, in the mesic grasslands of France, where tolerance is combined with avoidance, some species showed reduced SLA and LNC with cutting (Pontes et al., 2010(Pontes et al., , 2007. Our results were in agreement with those of humid grasslands since the mean annual precipitation in Iberian Atlantic grasslands is 1,330 mm. Fertilization favors a resource acquisition strategy and may also improve forage quality (Pontes et al., 2010(Pontes et al., , 2007. Since defoliation by herbivores is always coupled with fertilization by urine and feces (Lezama & Paruelo, 2016), community response is probably a combination of tolerance to defolia-tion and response to fertilization, but these two effects cannot be disentangled in our study.

Species turnover and ITV effects on the response to grazing abandonment
In agreement with the results of previous studies, the traits measured in the present study showed high ITV (Albert et al., 2010;Cornelissen et al., 2003;Lepš et al., 2011;Westoby et al., 2002). Pontes et al. (2007) observed higher between-species variation than within-species variation. We also found that species could be distinguished based on the mean values of traits; however, the exclusion of ITV resulted in a clear underestimation of community response to grazing abandonment for all traits. Similarly, Lee et al. (2017) claimed that intraspecific physiological response is an important driver of forage quality worldwide.
We observed a mismatch between species turnover and ITV effects on the response to grazing abandonment for H and SLA; the effects of species turnover were very weak, whereas those of ITV were very strong. The lack of clear association of the tall grasses F. microphylla and A. capillaris with E plots along with the abundance of the short dicot G. saxatile might explain the weak species turnover effect. Five years after grazing abandonment, opposite effects of species turnover and ITV were observed for SLA in Uzkuiti, as indicated by the negative covariance between species turnover and ITV ( Figure  3B). This response was probably a result of the high abundance of G. saxatile (very high SLA) in the E plot. However, long-term changes in floristic composition were observed after grazing abandonment (nine years in both sites): F. microphylla and A. capillaris became the most abundant species in the E plots (Odriozola et al., 2017). Our results are in agreement with those reported by Lepš et al. (2011), who suggested that ITV determines the resistance and stability of grasslands to short-term changes, while species turnover reflects processes associated with the resilience of grasslands to long-term changes.
The contributions of species turnover and ITV were more equilibrated for LNC and CN. T. repens, which showed very high LNC and very low CN (Table 2), was the most sensitive species to grazing abandonment and was clearly more abundant on G plots. This species is highly digestible and nutrient-rich, and is widely cultivated for animal nutrition (Parsons et al., 2011). Our results demonstrated that T. repens is fundamental for maintaining the high nutritional value of Iberian Atlantic grasslands.

Management implications
Previous studies have reported reductions in forage quantity along with improvement in quality under grazing conditions (Cruz et al., 2010;Deléglise et al., 2015;Moreno García et al., 2014). Accordingly, an increase in the annual aboveground net primary productivity has been observed in Iberian Atlantic grasslands after grazing abandonment (unpublished data). Cruz et al. (2010) claimed that the reduction in forage availability under grazing conditions is incompatible with the high animal productivity, whereas Moreno  argued that the improved forage quality might compensate the reduction in quantity. It is widely accepted that both forage digestibility and nutritive value have direct effects on animal performance (Herrero et al., 2015;Lee et al., 2017). Additionally, although short-term grazing abandonment positively affects grassland productivity, longterm grazing abandonment has negative effects because litter accumulation inhibits plant growth (Patton et al., 2007).
Our study revealed that ITV responded quickly to grazing abandonment and determined the short-term changes in forage quality. Strong species turnover effects may appear after long-term grazing abandonment when tall grasses with tissue of low quality outcompete high forage quality species (Odriozola et al., 2017). The use of mountain grasslands by domestic livestock has declined slowly in most of the Pyrenean-Basque-Cantabrian grassland systems (Ruiz et al., 2009), and significant management changes, including shorter duration in upland grasslands, lower stocking rate, and less shepherd control, have been identified (Aldezabal et al., 2015). Our study suggested that an active management system of livestock herds is necessary to prevent damage to grasslands associated with long periods of grazing abandonment. Shortterm ITV-based changes in forage nutritional value are easily reversible; however, recovery times are much longer when changes are caused by species turnover effects, and high tissue quality species may be reduced in abundance or eliminated. This study also identified T. repens as a key indicator of forage quality and successional status in Iberian Atlantic grasslands.