
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70282
ECE370282
ECE-2024-01-00162.R2
Biodiversity Ecology
Community Ecology
Restoration Ecology
Research Article
Research Article
High species richness of sheep‐grazed sand pastures is driven by disturbance‐tolerant and weedy short‐lived species
Kovacsics‐Vári et al.
Kovacsics‐Vári Gergely https://orcid.org/0000-0002-9995-5733
1
Sonkoly Judit https://orcid.org/0000-0002-4301-5240
1 2
Tóth Katalin 1 2
McIntosh‐Buday Andrea 1 2
Díaz Cando Patricia Elizabeth 1
Törő‐Szijgyártó Viktória 1 2
Balogh Nóra https://orcid.org/0000-0002-6156-0752
1
Guallichico Suntaxi Luis Roberto 2
Espinoza Ami Francis David 1
Matus Gábor 3
Tóthmérész Béla https://orcid.org/0000-0002-4766-7668
1 4
Török Péter https://orcid.org/0000-0002-4428-3327
1 2 5 molinia@gmail.com

1 Department of Ecology University of Debrecen Debrecen Hungary
2 HUN‐REN‐UD Functional and Restoration Ecology Research Group Debrecen Hungary
3 Department of Botany University of Debrecen Debrecen Hungary
4 HUN‐REN‐UD Biodiversity and Ecosystem Services Research Group Debrecen Hungary
5 Polish Academy of Sciences Botanical Garden‐Centre for Biological Diversity Conservation in Powsin Warszawa Poland
* Correspondence
Péter Török, Department of Ecology, University of Debrecen, 4032 Debrecen, Egyetem sqr. 1, Hungary.
Email: molinia@gmail.com

08 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7028222 8 2024
31 1 2024
27 8 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

We selected 15 sheep‐grazed sand pastures along a gradient of increasing grazing intensity to study the fine‐scale patterns of main biomass fractions (green biomass, litter) and that of plant species and functional groups (life forms and social behaviour types). We classified them into five grazing intensity levels based on stocking density, proximity to drinking and resting places and the number of faeces. We aimed to answer the following questions: (i) How does increasing intensity of sheep grazing affect the amount of green biomass, the species richness and their relationship in sand pastures? (ii) How does increasing intensity of sheep grazing affect the biomass of perennial and short‐lived graminoids and forbs? (iii) How does the disturbance value—expressed in the biomass ratio of disturbance‐tolerant and ruderal species—change along the gradient of grazing intensity? A unimodal relationship between green biomass and species richness was detected; however, the ordination (canonical correspondence analysis, CCA) showed no clustering of pastures subjected to the same levels of grazing intensity. Along the grazing intensity gradient we found an increasing trend in species richness and significant differences in green biomass (decreasing trend), litter (decreasing trend), graminoids (decreasing trend) and short‐lived forbs (increasing trend). We found an increasing amount of disturbance‐tolerant and ruderal species with increasing grazing intensity. We suggest that we might need to use multiple scales for sampling and a fine‐scale assessment of grazing intensity. Our findings might be instructive for pastures in densely populated regions, which are prone to the encroachment of disturbance‐tolerant and ruderal species.

We analysed the fine scale biomass patterns of sheep‐grazed sand pastures along increasing grazing intensity. Along the grazing intensity gradient, we found an increasing trend for species richness and the biomass of short‐lived forbs; while green biomass, litter and graminoid biomass showed a decreasing trend. We found an increasing amount of disturbance tolerant and ruderal species with increasing grazing intensity.

disturbance
European steppes
functional groups
humped‐back curve
plant biomass
sheep grazing
Ministry of Culture and Innovation, HungaryÚNKP‐23‐3‐I. New National Excellence Program Magyar Tudományos Akadémia 10.13039/501100003825 Bolyai János Scholarship (BO/00587/23/8) Nemzeti Kutatási, Fejlesztési és Innovaciós Alap 10.13039/501100012550 K137573 KKP144068 PD137747 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:08.09.2024
Kovacsics‐Vári, G. , Sonkoly, J. , Tóth, K. , McIntosh‐Buday, A. , Díaz Cando, P. E. , Törő‐Szijgyártó, V. , Balogh, N. , Guallichico Suntaxi, L. R. , Espinoza Ami, F. D. , Matus, G. , Tóthmérész, B. , & Török, P. (2024). High species richness of sheep‐grazed sand pastures is driven by disturbance‐tolerant and weedy short‐lived species. Ecology and Evolution, 14 , e70282. 10.1002/ece3.70282
==== Body
pmc1 INTRODUCTION

It is a truism in ecology that grazing has a crucial role in maintaining grassland biodiversity (Briske, 1996; Metera et al., 2010), although it is a disturbance that affects both the morphological characteristics and the functional trait composition of plant communities (Díaz et al., 2007; WallisDeVries et al., 2002). To avoid overgrazing, it is essential to study the effects of grazing on species composition, vertical and horizontal structure, regeneration capacity and functional composition of plant communities (Dong et al., 2012; Hao & He, 2019). Overgrazing is considered to be one of the most pressing problems for implementing sustainable management of pastures in many regions of the world (Gao & Li, 2016; Li, Zhang, et al., 2018; Török & Dengler, 2018). Former research found that plant community responses along a gradient of grazing intensity showed marked changes, which led to altered stability and ecosystem functioning (Li et al., 2018b; Su & Xu, 2021; Xu et al., 2022). When studying the effects of grazing, grazing intensity, livestock type and habitat type can be regarded as the main approaches, and various combinations of these approaches can be found in the scientific literature.

Sheep grazing has some specific characteristics; for example, (i) there is a higher preference for forbs compared to cattle grazing, (ii) sheep can consume plant parts closer to the ground and (iii) sheep prefer vegetative plant parts (Jerrentrup et al., 2015; Metera et al., 2010; Tóth et al., 2018). Sheep grazing supports seedling establishment on bare soil surfaces in different ways: (i) Seeds lying on the soil surface can be buried to an optimal depth for germination by sheep trampling (Eichberg et al., 2005). (ii) Flocks usually consist of up to a few hundred sheep and are frequently herded by shepherds over relatively long distances; consequently, they can contribute to a wide dispersal of certain seeds (Rosenthal et al., 2012). (iii) Sheep trampling opens dense vegetation cover and creates safe sites for seedling emergence and establishment (Faust et al., 2011; Freund et al., 2014).

Species richness often shows a humped‐back curve along a gradient of increasing disturbance, which can be explained by the intermediate disturbance hypothesis (Connell, 1978; Gao & Carmel, 2020). This finding was also supported by studies dealing with sheep grazing (del Pozo et al., 2006; Lázaro et al., 2016; Süss et al., 2007). However, the spatial (Süss et al., 2007) and temporal (del Pozo et al., 2006) scale of the study can influence the relationship, and monotonously decreasing species richness has also been detected, e.g. in sheep‐grazed desert steppes (Zhang et al., 2018).

Besides patterns of biodiversity, a further important aspect is how changes in species richness are reflected in the abundance of plant life and growth forms. Díaz et al. (2007) analysed plant trait responses to grazing in a meta‐analysis. They found that increasing grazing intensity favoured stoloniferous plants, rosette‐formation likeliness and short height and increased the abundance of short‐lived and fast‐growing species when climatic conditions and grazing history were both taken into consideration. Some findings of the above‐mentioned meta‐analysis were partly confirmed for sheep‐grazed pastures by Pettit et al. (1995), Yang et al. (2022) and Farmilo et al. (2023). However, none of the above papers studied the biomass of life and growth form groups at multiple levels of grazing intensity in sheep‐grazed pastures. Compared to classical analyses of species composition and richness, functional trait‐ and functional group‐based analyses provide a more direct view for understanding the link between community assembly and functioning and how it is affected by livestock grazing (Carboni et al., 2023; Carmona et al., 2015).

Avoiding the encroachment of weedy and/or unwanted species in grazed grasslands is among the most important challenges of sustainable management (Bretas et al., 2023). While the formation of vegetation gaps by grazing is essential for the establishment of short‐lived pioneer species, it also provides colonisation windows for weeds and invasive species. To evaluate the effect of grazing intensity on the encroachment of weedy species is especially crucial for sand pastures as these types of grasslands are highly exposed to invasion (Botta‐Dukát, 2008). In some ecosystems, sheep grazing is an established method for supressing weedy forbs, mostly because sheep are considered to be highly selective grazers, which prefer forbs over grasses and shrubs (Jerrentrup et al., 2015; Olson & Lacey, 1994). However, it has also been stressed that the effect of livestock grazing is strongly context‐dependent and varying in different habitat types (Liu et al., 2015; Mládek et al., 2013; Török et al., 2024). Many sand grasslands in Central Europe are classified in the EU Habitats Directive as Pannonian and Pontic sandy steppes (E1.1a). These grasslands are located in Central and Southeast Europe and critically endangered according to the European red list of habitats (EC Directorate‐General for Environment et al., 2017), which underlines their conservation importance. For their sustainable management, it is crucial to understand how their vegetation composition and functional diversity respond to different grazing regimes.

With the study of sand pastures subjected to increasing intensity of sheep grazing, we aimed to address the following questions: (i) How does increasing intensity of sheep grazing affect the amount of green biomass, the species richness and their relationship in sand pastures? (ii) How does increasing intensity of sheep grazing affect the biomass of perennial and short‐lived graminoids and forbs? (iii) How does the disturbance value (expressed in the biomass ratio of disturbance‐tolerant and ruderal species) change along the gradient of grazing intensity?

2 MATERIALS AND METHODS

2.1 Study area

For our study, we selected 15 sand pastures (Table 1) in the Nyírség region, East Hungary, where there is a high proportion of man‐made habitats such as croplands and tree plantations (Botta‐Dukát, 2008). The Nyírség region is characterised by an annual rainfall ranging between 530 and 680 mm and an average annual temperature between 9.4 and 9.8°C (Dövényi, 2010). In some years, the annual rainfall is even less than 400 mm so that serious drought events occur (Négyesi, 2018). The soil types are Dystric to Brunic Arenosols, which are characterised by a rather low humus content (0.6%–2.6% by mass) and a pH ranging typically from 4.45 to 5.71. Site codes are given in Table 1, further soil characteristics in Appendix S1. All investigated pastures are leased by the shepherds from the Hortobágy National Park Directorate except pasture 13, which is leased from the village authorities. The vegetation of the sand pastures is characterised by a high cover of tussock‐forming (Festuca pseudovina, F. vaginata and Corynephorus canescens) and stoloniferous (Poa angustifolia, Cynodon dactylon, Carex praecox and C. stenophylla) graminoids and a relatively high cover of short‐lived and perennial sand grassland forbs (Scleranthus annuus, Vicia lathyroides, Potentilla arenaria, Rumex acetosella, Eryngium campestre and Chondrilla juncea) (Figure 1).

TABLE 1 Main characteristics of the studied sites. The site codes were used in the canonical correspondence analysis (CCA) (Figure 5). The levels of grazing intensity (increasing from ‘1’ to ‘5’) are explained in Table 2.

Site code	Settlement name	GPS coordinates	Elevation (m)	Area of pasture (ha)	Level of grazing intensity	
1	Létavértes	N47.41797	E21.89950	116	200	5	
2	Létavértes	N47.42269	E21.91110	117	200	4	
3	Létavértes	N47.44133	E21.92817	120	130	2	
4	Hajdúbagos	N47.41580	E21.68023	107	140	2	
5	Hajdúbagos	N47.41183	E21.68336	109	140	2	
6	Létavértes	N47.42478	E21.86181	116	200	5	
7	Létavértes	N47.42928	E21.86382	114	200	5	
8	Hajdúsámson (Martinka)	N47.57292	E21.77960	129	126	3	
9	Hajdúsámson (Martinka)	N47.57391	E21.78137	130	126	3	
10	Hajdúsámson (Martinka)	N47.57486	E21.79247	130	76	3	
11	Monostorpályi	N47.40945	E21.77691	112	45	3	
12	Monostorpályi	N47.41554	E21.78521	112	45	3	
13	Hajdúsámson (Martinka)	N47.58068	E21.77093	132	57	4	
14	Hajdúsámson (Martinka)	N47.57541	E21.79281	131	Grazing exclusion	1	
15	Vámospércs	N47.53241	E21.95022	132	Grazing exclusion	1	

FIGURE 1 Typical view of sampled sand pastures in the surroundings of Hajdúbagos (a), Létavértes (b) and Hajdúsámson (Martinka district) (c, d), East Hungary. Photo credit: Judit Sonkoly.

2.2 Sampling

We sampled the biomass of 15 pastures from late May to early June 2021. The pastures were managed by seasonal sheep herding (Merino breed, typically from early April to the end of October). According to the national park rangers, the grazing intensity in any given site was similar between 2017 and 2021. Two sites had been fenced for 13 years (since the summer of 2008) to exclude livestock grazing (Aszalósné Balogh et al., 2023). The pastures were selected to cover a broad range of grazing intensity (Table 1) considering the stocking density, the proximity of (frequently used) drinking and/or resting places, the number of faeces (dropping events) and the physiognomy of the vegetation (for details see Tables 1 and 2). We also used satellite imaging to assess the distance to the resting/drinking places. We avoided sites where the prevailing disturbance was not grazing and sites with markedly different habitat types (e.g. forest plantations, dirt roads and drained low‐lying areas). Information on approximate grazing intensities provided by national park rangers were refined during the field sampling by interviewing the shepherds. National park rangers also helped in site selection and confirmed that the stocking density alone is not sufficient to assess the grazing intensity, an assumption also supported by some former studies (e.g. Tonn et al., 2019). Taking all this information into account, we classified our sampling sites into five grazing intensity categories (Table 2).

TABLE 2 Levels of grazing intensity depending on stocking density, proximity to resting and drinking places and units of faeces found per sampling site. One sheep is equivalent to 0.2 livestock units per hectare (LU/ha), and one unit of faeces is equal to one dropping event. Grazing intensity increases from ‘1’ to ‘5’, where ‘1’ means grazing exclosure and ‘5’ is the highest level of grazing intensity.

Level of grazing intensity	Stocking density (LU/ha)	Proximity to resting and drinking places (m)	Units of faeces (10 m × 10 m area)	
1	Grazing exclusion	
2	0.5–0.8	>150	0–20	
3	0.5–0.8	<150	>20	
4	1.1–4	>150	0–20	
5	1.1–4	<150	>20	

In each pasture, we designated a 10 m × 10 m sampling site to ensure uniform biomass sample heterogeneity. Before the biomass sampling, we recorded the complete list of vascular plant species in each site to facilitate the biomass sorting in the laboratory. In each sampling site, we harvested the total aboveground biomass of ten 20 cm × 20 cm sampling plots (altogether 150 samples) using secateurs. Standing litter and the litter layer were also included in the samples. The samples were dried using a drying oven (65°C for 48 h). After drying the biomass was sorted into the main fractions moss, lichen, litter (including both the litter layer and standing litter) and green biomass. Green biomass was further sorted to vascular plant species level, whereas the moss and lichen fractions were not further differentiated. The sorted biomass fractions were weighed using an electronic balance (accuracy: ±0.01 g).

During the biomass sampling, we also collected soil samples from the upper 5 cm soil layer of each biomass sampling plot to characterise the average site properties. The ten samples per site were pooled to eventually obtain at least 500 g air‐dried soil for each sampling site. Soil samples were analysed by an accredited laboratory for physical soil type, pH, humus, NO2 − and NO3 −, K2O, P4 O10, CaCO3 and total of water‐soluble salts (Appendix S1).

2.3 Data processing and analyses

For the functional analysis of the sorted biomass of vascular plant species, we obtained regional plant trait data from the Pannonian Database of Plant Traits (PADAPT, Sonkoly et al., 2023). We classified the species into simplified morpho‐functional groups of short‐lived forbs, short‐lived graminoids, perennial forbs and perennial graminoids according to PADAPT and Király (2009) and into social behaviour types (SBTs) following Borhidi (1995). SBT is a refined Grime's CSR classification scheme adapted to the Hungarian flora, which has a robust and expanded subclassification for the ruderal strategy type. Using the SBT classification system we grouped the species into three categories along a gradient of increasing disturbance tolerance: (1) sand grassland species including the categories competitors (C), specialists (S), generalists (G) and natural pioneers (NP) of sand grasslands, (2) natural disturbance‐tolerant species (DT) and (3) ruderal weedy species including the categories ruderal competitors (RC), 0000 (AC) and weeds (W). For each plot we calculated community‐weighted means (CWMs) of this ordinal variable (Groups 1, 2 and 3, respectively) weighted by biomass and used it as an ecological indicator of disturbance (disturbance index) in the analyses. We used generalised linear mixed‐effect models (GLMMs) to assess the impact of sheep grazing on dependent variables (intensity level included as fixed factor, site identity as random factor and dependent variables listed in Table 3; SPSS 26.0 program package, IBM Corp, 2019). We plotted species richness against green biomass and analysed their relationship using a second‐order polynomial fit. The biomass composition of sites and grazing intensities were explored by canonical correspondence analysis (CCA) using the CANOCO 5.0 programme package (Šmilauer & Lepš, 2014). We included seven variables (soil phosphorous content, pH, soil compactness, soil nitrogen content, soil potassium content, disturbance index and soil humus content) in the secondary explanatory matrix of the CCA and selected significant predictors by a Monte‐Carlo permutation test.

TABLE 3 Effects of the intensity of sheep grazing on species richness, disturbance and biomass fractions.

Variable	Grazing intensity	
F 4,145	p	
Species richness	5.66	<.001	
Disturbance	6.24	<.001	
Main biomass fraction	
Green biomass	2.49	.046	
Litter	10.97	<.001	
Mosses + lichens	0.953	.435	
Specific biomass fraction	
Perennial forbs	25.88	<.001	
Perennial graminoids	3.21	.015	
Short‐lived forbs	11.20	<.001	
Short‐lived graminoids	0.27	.897	
Note: Significant values (generalised linear mixed‐effect model (GLMM), p < .05) are highlighted in boldface.

3 RESULTS

We detected 84 species in the samples consisting of 24 graminoids (8 short‐lived and 16 perennial) and 60 forbs (36 short‐lived and 24 perennial). Grazing intensity had a significant effect on green biomass, litter, the biomass of perennial forbs, perennial graminoids and short‐lived forbs and on species richness and disturbance index (Table 3). The species richness at the fourth and fifth level of grazing intensity was significantly higher than that at the second level, and we detected the highest species richness at the fourth level of grazing intensity (Figure 2a). The disturbance index was significantly higher at the fifth level compared to all other levels (Figure 2b). We detected a humped‐back relationship between the amount of green biomass and species richness (Figure 3).

FIGURE 2 Species richness (a) and disturbance index (b) along a grazing intensity gradient obtained by a generalised linear mixed‐effect model (GLMM). Circles denote estimated means, whiskers standard errors. Different letters denote significant differences. The grazing intensity levels 1 to 5 are explained in Table 2.

FIGURE 3 Relationship between the green biomass and species richness obtained by a second‐order polynomial fit (p < .001, R 2 = 0.143). The ‘+’ signs denote the centroids of the samples originating from the same site and grazing intensity, and all other symbols denote the estimated means of individual samples collected on site. The green shade indicates the 95% confidence interval. Different symbols and their colours denote the level of grazing intensity as follows: 1: Green circles, 2: Blue triangles, 3: Orange diamonds, 4: Red rectangles and 5: Violet stars. The levels of grazing intensity are explained in Table 2.

We also detected significant differences for green biomass and litter. Both were highest at the first and lowest at the fifth level of grazing intensity (Figure 4a,b). Moss and lichen biomass did not show significant differences along the grazing intensity gradient. Perennial forb biomass was highest at the fourth level of grazing intensity (Figure 5a), whereas perennial graminoid biomass was significantly higher at the second than at the fourth and fifth level (Figure 5b). We found no significant differences in the biomass of short‐lived forbs between the first and fourth level of grazing intensity, whereas that at the fifth level was significantly higher compared to the other levels (Figure 5c). The biomass of short‐lived graminoids did not show significant differences (Table 3; Figure 5d).

FIGURE 4 Green biomass (a) and litter (b) along a grazing intensity gradient obtained by a generalised linear mixed‐effect model (GLMM). Circles denote estimated means, whiskers standard errors. Different letters denote significant differences. The grazing intensity levels 1 to 5 are explained in Table 2.

FIGURE 5 Biomass of perennial forbs (a), perennial graminoids (b), short‐lived forbs (c) and short‐lived graminoids (d) along a grazing intensity gradient obtained by a generalised linear mixed‐effect model (GLMM). Circles denote estimated means, whiskers standard errors. Different letters denote significant differences. The grazing intensity levels 1 to 5 are explained in Table 2.

The CCA identified only two significant predictors for the species compositional patterns, soil humus content and disturbance index, but these two predictors were weakly correlated with each other (Figure 6). We found no clear separation of pastures grazed at different levels of intensity based on the biomass composition, and sites with increasing levels of grazing intensity were also not separated along the predictors. Similarly, the samples of the different sites did not show any clear separation along either the axis of green biomass quantity or the axis of species richness (Figure 3).

FIGURE 6 Relationship between species composition, humus content and level of disturbance (community‐weighted means (CWM) of social behaviour types (SBT), see text). For the canonical correspondence analysis (CCA), biomass per species was used. Eigenvalues were 0.617 and 0.504 for the first and second axes, respectively. The cumulative percentage variance of the species‐environment relation was 82.2 for the first four axes. Numbers denote the sampling sites (see Table 1) and their colours the levels of grazing intensity (1: Green, 2: Blue, 3: Orange, 4: Red and 5: Violet). Species names are abbreviated to the first four letters of the genus and the first four letters of the specific epithet. Full names of species are given in Appendix S3. Only significant explanatory variables are shown in the figure (499 permutations, p < .008).

4 DISCUSSION

4.1 Species richness and green biomass

We found a humped‐back relationship between species richness and green biomass. Moderate levels of disturbance provide favourable conditions for a wider range of species (Metera et al., 2010). The highest species richness was detected at the fourth and fifth grazing intensity levels. In accordance with our results several studies detected that species richness is higher under moderate grazing pressure than, for example, in ungrazed pastures (Deng et al., 2013; Fensham et al., 2011; Sasaki et al., 2009). Deng et al. (2013) found the density, the height and the cover of vegetation to be highest in the ungrazed plot, where also the densities of dominant, good competitor species were highest. When grazing intensity is low, more competitive plants can grow taller as they acquire more resources and can thus overcome resource limitation more effectively and grow efficiently in dense, ungrazed vegetation (He et al., 2021; Westoby, 1999). Livestock grazing decreases green biomass and litter (Magnano et al., 2019), and this effect is obviously more pronounced under more intensive grazing. This is also reflected in the high proportion of biomass samples belonging to the first and second level of grazing intensity at higher biomass values (Figure 3). Although we detected a humped‐back relationship between species richness and green biomass in the studied pastures, the possibility of a decreasing trend should not be excluded in unusually arid years (Gao & Carmel, 2020; Milchunas et al., 1988).

One would expect that the samples originating from the same grazing intensity levels would cluster in the ordination, but this was not validated by our results (Figure 6). The reason for this might be a fine‐scale heterogeneity of vegetation. Studying the vegetation composition of alkaline and sand pastures at multiple‐scales, Godó et al. (2017) observed that there is a significant relationship between plot size and grazing effects: With increasing plot size they found decreasing levels of differences in species composition, ie the small‐sized plots showed higher beta diversity than larger ones. We assume this to be the reason why our fine‐scale samples (harvested in 20 cm × 20 cm plots) were not clearly separated along the gradient of green biomass and species richness.

4.2 Main biomass fractions and the biomass of life forms

We found significant differences along the grazing intensity gradient for the main biomass fractions and the biomass of functional groups. Both the green biomass and the litter fraction were significantly lower at higher grazing intensities. Obviously, the consumption of plants by livestock contributes to a decrease in green biomass and litter (Magnano et al., 2019). According to Kemp et al. (2000) perennial graminoids are most sensitive to grazing, and with increasing grazing disturbance subordinated species are able to spread (Grime & Mackey, 2002). This might explain why we detected significantly lower amounts of perennial graminoids at higher grazing intensities. Green biomass showed significant differences, and though the biomass of short‐lived forbs was significantly higher at the highest grazing intensity, total green biomass decreased with increasing grazing intensity (as one would expect). There are several possible explanations for this pattern. First, sheep grazing may be less selective for forbs at higher grazing intensities (Golodets et al., 2009; Tóth et al., 2018), and the feeding strategy of sheep is flexible depending on the available biomass and/or forb richness of the pasture, displaying a mass foraging strategy in communities characterised by graminoids (Liu et al., 2015; Mládek et al., 2013). This causes a net loss of total green biomass, but at the same time favours short‐lived forbs due to their fast regrowth rate and colonisation ability in gaps (Hofmann & Isselstein, 2004; Westoby, 1999). Second, the mean height of species is typically lower at higher grazing intensities (Deng et al., 2013; Török et al., 2016) so that they may be represented by less biomass. Third, annuals have a higher specific leaf area (SLA) and a lower leaf dry matter content (LDMC) and thus a lower dry weight than other plants (with a lower SLA and a higher LDMC) (E‐Vojtkó et al., 2020). In their meta‐analysis of plant responses to grazing, Díaz et al. (2007) found that the abundance of perennial plant species decreased with increasing grazing intensity. This was also confirmed by our results to some extent, but a striking leap can be observed at the fourth level of grazing intensity for perennial forbs (Figure 5a). Among the detected perennial forbs Thymus glabrescens had remarkably high values at this grazing intensity level. Without this species there would be a continuous decline of perennial forb biomass along the gradient of grazing intensity. Their stout woody stems increase the biomass of perennial forbs, and their large amount of biomass shapes the differences of the total green biomass scores. Thymus glabrescens produces monoterpenes, which are known to act as feeding deterrents against sheep (Linhart & Thompson, 1999). It is therefore probably one of the less favoured food plants on the pastures we investigated. The meta‐analysis mentioned above (Díaz et al., 2007) details the response of annuals to grazing. It found them to increase together with grazing intensity, which is confirmed by our results in case of forbs. Another explanation for the spreading success of short‐lived plants is their generally higher SLA compared to perennials. A high SLA is linked to fast re‐growth ability (Helm et al., 2019), and higher SLA scores at higher grazing intensities were confirmed by a former study addressing livestock grazing in sand grasslands (Kovacsics‐Vári et al., 2023). We found no significant differences for short‐lived graminoids, which is presumably due to their low species number and the even distribution of the species along the gradient of grazing intensity. For example, Apera spica‐venti had higher biomass at lower levels of grazing intensity, whereas Bromus hordeaceus was stronger at higher grazing intensities.

4.3 Impact of grazing disturbance on species composition

Significantly higher biomass scores of disturbance‐tolerant and ruderal species were detected at the highest level of grazing intensity, where their species numbers were high, but they suppressed the biomass of characteristic species of sand grasslands. The process leading to overgrazing has been succinctly summarised by Schulze (2019): First, the vegetation composition changes; second, vegetation cover decreases; third, bare soil surfaces are formed and fourth, soil erosion becomes more severe. Our results might partly support this description, but we presume that these stages are not reached one after another, but occur in parallel. For example, when bare soil surfaces are formed due to the suppression of characteristic perennial grass and forb species, disturbance‐tolerant and ruderal species will take the opportunity and establish quickly, which changes the vegetation composition. Midolo et al. (2023) assessed plant species by disturbance categories and found that annuals are favoured by disturbance as their ability to grow fits to circumstances that do not provide stable biotic and abiotic features, which would otherwise favour better competitors using resources efficiently on the long‐run (Salguero‐Gómez, 2017; Schulze, 2019). According to Botta‐Dukát (2008), open sand grasslands are among the habitats most exposed to disturbance in Hungary. The Nyírség region has been densely populated for centuries, and despite the low productivity of sand grasslands, large areas were cultivated, resulting in high covers of disturbance‐tolerant and ruderal species. Pastures with open surfaces can quickly become colonised by these species, which typically have a good dispersal ability (Schulze, 2019), and according to our findings they even increased in species richness. This draws the attention to landscapes where pastures have to be sustainably managed in a densely populated region.

Our results showed that open surfaces and colonisation by disturbance‐tolerant and ruderal species are to be expected in intensively grazed pastures. The fourth and fifth level of grazing intensity included sampling sites with relatively high stocking densities (1.1 to 4 LU/ha), and the highest scores of species richness were observed in these pastures due to the higher number of short‐lived – mainly disturbance‐tolerant and ruderal – species. Increasing species richness in more intensively grazed sites was also confirmed by Kiss et al. (2006). In their assessment of grazing intensity,0 they considered the proximity of study sites to stables, and similar to us they detected a high proportion of disturbance‐tolerant and ruderal species contributing to the higher species richness. We detected significantly higher species richness at the fifth and fourth level of intensity compared to the second level. At the fourth level, the presence of grassland species was still substantial, but disturbance‐tolerant and ruderal species already occurred more frequently. One can find significant differences in most of the studied characteristics if the second and fourth levels of grazing intensity are compared to the third and fifth (Figures 1 and 3), which implies that the sites representing the third and fifth level were grazed and trampled more frequently. These significant differences suggest that besides the stocking density the proximity to resting and drinking places has a noticeable effect on the vegetation as well. We assume that a closer proximity to resting and watering points results in a higher grazing frequency and intensity. In their study on the effect of proximity to resting and drinking places on plant characteristics, Kovacsics‐Vári et al. (2023) found similar results for flowering period and life forms. In this context, it should be noted that the findings by Kiss et al. (2006) and Tonn et al. (2019) suggest that the effects of livestock grazing should be studied on finer scales than the scale of a pasture since the stocking density, as an important metric in grazing regimes, is an average figure for the entire pasture area and does not take local variations into account.

5 CONCLUSIONS

We found that disturbance‐tolerant and ruderal species contributed to the increase of species richness at higher levels of grazing intensity, but also suppressed some of the characteristic species, which might lead to decreasing stability and ecosystem functions. To interpret this finding data on proximity to frequently used places seemed useful. We hope that our results will help in the arrangement of pastures and be a decision‐making aid with regard to a site‐adapted stocking density in situations where conservation grazing is necessary to maintain characteristic and reduce the encroachment of non‐desired species. One implication of our findings is that resting and drinking places should not be placed close to the field margins and/or to roads or ruderal habitats, which could act as sources of weed propagules. We detected a humped‐back relationship between green biomass and species richness in the studied sand pastures, but the plots characterised by different grazing intensities were not clearly separated from each other along the biomass gradient. This clearly indicates that (i) levels of grazing intensity based on the stocking density on pasture scale are too robust to assess effects of grazing intensity on the vegetation and (ii) even sheep grazing created local variations in the vegetation with respect to species richness and biomass at the scale of the biomass sampling. These results suggest that we might need to use multiple scales for sampling and a fine‐scale assessment of grazing intensity. We formulated assumptions regarding the impact of proximity on grazing frequency, but to find clearer effects, a broader knowledge about each pasture would be needed.

AUTHOR CONTRIBUTIONS

Gergely Kovacsics‐Vári: Conceptualization (supporting); data curation (equal); investigation (equal); methodology (supporting); project administration (supporting); validation (supporting); visualization (supporting); writing – original draft (supporting); writing – review and editing (supporting). Judit Sonkoly: Conceptualization (supporting); data curation (supporting); funding acquisition (supporting); investigation (supporting); validation (supporting); writing – review and editing (supporting). Katalin Tóth: Validation (lead); writing – review and editing (supporting). Andrea McIntosh‐Buday: Data curation (supporting); investigation (supporting); validation (supporting); writing – review and editing (supporting). Patricia Elizabeth Díaz Cando: Investigation (supporting); writing – review and editing (supporting). Viktória Törő‐Szijgyártó: Data curation (supporting); investigation (supporting); writing – review and editing (supporting). Nóra Balogh: Investigation (supporting); writing – review and editing (supporting). Luis Roberto Guallichico Suntaxi: Investigation (supporting); writing – review and editing (supporting). Francis David Espinoza Ami: Investigation (supporting); writing – review and editing (supporting). Gábor Matus: Resources (supporting); writing – review and editing (supporting). Béla Tóthmérész: Conceptualization (supporting); writing – review and editing (supporting). Péter Török: Conceptualization (lead); data curation (equal); formal analysis (equal); funding acquisition (lead); investigation (equal); methodology (lead); project administration (lead); validation (supporting); visualization (lead); writing – original draft (lead); writing – review and editing (equal).

CONFLICT OF INTEREST STATEMENT

None declared.

Supporting information

Appendix S1. Soil properties of the studied sand pastures.

Appendix S2. Raw biomass data of the grazed pastures.

Appendix S3. Species composition data of the grazed pastures.

ACKNOWLEDGEMENTS

We are thankful to the rangers of the Hortobágy National Park Directorate (Krisztián Pompola, László Demeter, László Széll, Norbert Patalenszki, Zoltán Kovács) for the information they provided on the studied pastures. We are also thankful to Luca Di Vita (BSc student, University of Palermo) for the help field works and data collection. Aiko Huckauf (Kiel) kindly polished the paper linguistically. The authors were supported by the ÚNKP‐23‐3‐I. New National Excellence Program of the Ministry for Culture and Innovation from the source of the National Research, Development and Innovation Found (GK‐V, AM‐B, VT‐SZ). The authors were supported by the National Research, Development and Innovation Office [PT: KKP 144068, K 137573; JS: PD 137747] during manuscript preparation. The work of JS was supported by the Bolyai János Scholarship of the Hungarian Academy of Sciences [BO/00587/23/8]. This project has received funding from the HUN‐REN Hungarian Research Network.

DATA AVAILABILITY STATEMENT

Underlying data are provided in Appendices [Link], [Link].
==== Refs
REFERENCES

Aszalósné Balogh, R. , Farkas, E. , Tüdősné Budai, J. , Lőkös, L. , & Matus, G. (2023). Cryptogamic biomass in Pannonic acidic sand steppes subject to changing land‐use. Plants, 12 , 2972. 10.3390/plants12162972 37631183
Borhidi, A. (1995). Social behaviour types, the naturalness and relative ecological indicator values of the higher plants in the Hungarian Flora. Acta Botanica Hungarica, 39 , 97–181.
Botta‐Dukát, Z. (2008). Invasion of alien species to Hungarian (semi‐)natural habitats. Acta Botanica Hungarica, 50 (Suppl ), 219–227. 10.1556/abot.50.2008.suppl.11
Bretas, I. L. , Dubeux, J. C. B. , Cruz, P. J. R. , Queiroz, L. M. D. , Ruiz‐Moreno, M. , Knight, C. , Flynn, S. , Ingram, S. , Pereira Neto, J. D. , Oduor, K. T. , Loures, D. R. S. , Novo, S. F. , Trumpp, K. R. , Acuña, J. P. , & Bernardini, M. A. (2023). Monitoring the effect of weed encroachment on cattle behavior in grazing systems using GPS tracking collars. Animals, 13 (21 ), 3353. 10.3390/ani13213353 37958108
Briske, D. (1996). Strategies of plant survival in grazed systems: A functional interpretation. In J. Hodgson & A. Illius (Eds.), The ecology and Management of Grazed Systems (pp. 37–67). CAB International.
Carboni, L. J. , Yahdjian, L. , & Oñatibia, G. R. (2023). Effects of livestock grazing intensification on plant communities of Patagonian drylands increase with increasing aridity. Applied Vegetation Science, 26 , e12754. 10.1111/avsc.12754
Carmona, C. P. , Mason, N. W. H. , Azcárate, F. M. , & Peco, B. (2015). Inter‐annual fluctuations in rainfall shift the functional structure of Mediterranean grasslands across gradients of productivity and disturbance. Journal of Vegetation Science, 26 , 538–551. 10.1111/jvs.12260
Connell, J. H. (1978). Diversity in tropical rain forests and coral reefs. Science, 199 , 1302–1310. 10.1126/science.199.4335.1302 17840770
del Pozo, A. , Ovalle, C. , Casado, M. A. , Acosta, B. , & de Miguel, J. M. (2006). Effects of grazing intensity in grasslands of the Espinal of central Chile. Journal of Vegetation Science, 17 , 791–798. 10.1111/j.1654-1103.2006.tb02502.x
Deng, L. , Sweeney, S. , & Shangguan, Z.‐P. (2013). Grassland responses to grazing disturbance: Plant diversity changes with grazing intensity in a desert steppe. Grass and Forage Science, 69 , 524–533. 10.1111/gfs.12065
Díaz, S. , Lavorel, S. , McIntyre, S. , Falczuk, V. , Casanoves, F. , Milchunas, D. G. , Skarpe, C. , Rusch, G. , Sternberg, M. , Noy‐Meir, I. , Landsberg, J. , Zhang, W. , Clark, H. , & Cambell, B. D. (2007). Plant trait responses to grazing – A global synthesis. Global Change Biology, 13 , 313–341. 10.1111/j.1365-2486.2006.01288.x
Dong, S. K. , Wen, L. , Li, Y. Y. , Wang, X. X. , Zhu, L. , & Li, X. Y. (2012). Soil‐quality effects of grassland degradation and restoration on the Qinghai‐Tibetan plateau. Soil Science Society of America Journal, 76 , 2256–2265. 10.2136/sssaj2012.0092
Dövényi, Z. (Ed.). (2010). Cadastral inventory of landscapes in Hungary (Hungarian) (p. 876). Budapest.
EC Directorate‐General for Environment , Tsiripidis, I. , Piernik, A. , Janssen, J. A. M. , Tahvanainen, T. , Molina, J. A. , Giusso del Galdo, G. , Gardfjell, H. , Dimopoulos, P. , Šumberová, K. , Acosta, A. , Biurrun, I. , Poulin, B. , Hájek, M. , Bioret, F. , Essl, F. , Rodwell, J. R. , García Criado, M. , Schaminée, J. H. J. , … Gigante, D. (2017). European red list of habitats. Part 2, Terrestrial and freshwater habitats. European Union. 10.2779/091372
Eichberg, C. , Storm, C. , & Schwabe, A. (2005). Epizoochorous and post‐dispersal processes in a rare plant species: Jurinea cyanoides (L.) Rchb. (Asteraceae). Flora, 200 , 477–489. 10.1016/j.flora.2005.02.005
E‐Vojtkó, A. , Balogh, N. , Deák, B. , Kelemen, A. , Kis, S. Z. , Kiss, R. , Lovas‐Kiss, Á. , Löki, V. , Lukács, K. , Molnár, V. A. , Nagy, T. , Sonkoly, J. , Süveges, K. , Takács, A. , Tóth, E. , Tóth, K. , Tóthmérész, B. , Török, P. , Valkó, O. , … Lukács, B. A. (2020). Leaf trait records of vascular plant species in the Pannonian flora with special focus on endemics and rarities. Folia Geobotanica, 55 , 73–79. 10.1007/s12224-020-09363-7
Farmilo, B. , Duncan, D. , Moxham, C. , Buchan, A. , Turner, V. , Ballentine, M. , White, M. , & Schultz, N. (2023). Transient shifts in composition of degraded temperate native grassland following grazing exclusion. Applied Vegetation Science, 26 , e12731. 10.1111/avsc.12731
Faust, C. , Eichberg, C. , Storm, C. , & Schwabe, A. (2011). Post‐dispersal impact on seed fate by livestock trampling – A gap of knowledge. Basic and Applied Ecology, 12 , 215–226. 10.1016/j.baae.2011.02.009
Fensham, R. J. , Silcock, J. L. , & Dwyer, J. M. (2011). Plant species richness responses to grazing protection and degradation history in a low productivity landscape. Journal of Vegetation Science, 22 , 997–1008. 10.1111/j.1654-1103.2011.01305.x
Freund, L. , Eichberg, C. , Retta, I. , & Schwabe, A. (2014). Seed addition via epizoochorous dispersal in restoration: An experimental approach mimicking the colonization of bare soil patches. Applied Vegetation Science, 17 , 74–85. 10.1111/avsc.12052
Gao, J. , & Carmel, Y. (2020). Can the intermediate disturbance hypothesis explain grazing–diversity relations at a global scale? Oikos, 129 , 493–502. 10.1111/oik.06338
Gao, J. , & Li, X. (2016). Degradation of frigid swampy meadows on the Qinghai‐Tibet plateau: Current status and future directions of research. Progress in Physical Geography, 40 , 794–810. 10.1177/0309133316659283
Godó, L. , Valkó, O. , Tóthmérész, B. , Török, P. , Kelemen, A. , & Deák, B. (2017). Scale‐dependent effects of grazing on the species richness of alkaline and sand grasslands. Tuexenia, 37 , 229–246. 10.14471/2017.37.016
Golodets, C. , Sternberg, M. , & Kigel, J. (2009). A community‐level test of the leaf‐height‐seed ecology strategy scheme in relation to grazing conditions. Journal of Vegetation Science, 20 , 392–402. 10.1111/j.1654-1103.2009.01071.x
Grime, J. P. , & Mackey, J. M. L. (2002). The role of plasticity in resource capture by plants. Evolutionary Ecology, 16 , 299–307. 10.1023/A:1019640813676
Hao, Y. , & He, Z. (2019). Effects of grazing patterns on grassland biomass and soil environments in China: A meta‐analysis. PLoS One, 14 , e0215223. 10.1371/journal.pone.0215223 31009490
He, Q. , Jiang, K. , Hou, W. , Zhao, Y. , Sun, X. , Wang, L. , Zou, Y. , Zhu, Z. , & Zhang, H. (2021). Grazing alters species relative abundance by affecting plant functional traits in a Tibetan subalpine meadow. Ecology and Evolution, 11 , 11028–11037. 10.1002/ece3.7891 34429900
Helm, J. , Dutoit, T. , Saatkamp, A. , Bucher, S. F. , Leiterer, M. , & Römermann, C. (2019). Recovery of Mediterranean steppe vegetation after cultivation: Legacy effects on plant composition, soil properties and functional traits. Applied Vegetation Science, 22 , 71–84. 10.1111/avsc.12415
Hofmann, M. , & Isselstein, J. (2004). Seedling recruitment on agriculturally improved mesic grassland: The influence of disturbance and management schemes. Applied Vegetation Science, 7 , 193–200. 10.1111/j.1654-109X.2004.tb00610.x
IBM Corp . (2019). IBM SPSS Statistics for Windows (Version 26.0) [Computer software]. IBM Corp.
Jerrentrup, J. S. , Seither, M. , Petersen, U. , & Isselstein, J. (2015). Little grazer species effect on the vegetation in a rotational grazing system. Agriculture, Ecosystems & Environment, 202 , 243–250. 10.1016/j.agee.2015.01.007
Kemp, D. R. , Michalk, D. L. , & Virgona, J. M. (2000). Towards more sustainable pastures: lessons learnt. Australian Journal of Experimental Agriculture, 40 , 343–356. 10.1071/EA99001
Király, G. (Ed) (2009). New Hungarian Herbal. The Vascular Plants of Hungary (Hungarian). Identification key. Jósvafő, Aggtelek National Park Directorate, 616 p.
Kiss, T. , Malatinszky, Á. , & Penksza, K. (2006). Comparative coenological examinations on pastures of the great Hungarian plain I. (horse and cattle pasture near Hódmezővásárhely). Tájökológiai Lapok, 4 , 339–346. 10.56617/tl.4471
Kovacsics‐Vári, G. , Sonkoly, J. , Tóth, K. , McIntosh‐Buday, A. , Díaz Cando, P. , Törő‐Szijgyártó, V. , Balogh, N. , Guallichico Suntaxi, L. R. , Espinoza Ami, F. D. , Demeter, L. , Tóthmérész, B. , & Török, P. (2023). Intensity‐dependent effects of cattle and sheep grazing in sand grasslands ‐ does livestock type really matter? Applied Vegetation Science, 26 , e12727. 10.1111/avsc.12727
Lázaro, A. , Tscheulin, T. , Devalez, J. , Nakas, G. , Stefanaki, A. , Hanlidou, E. , & Petanidou, T. (2016). Moderation is best: Effects of grazing intensity on plant‐flower visitor networks in Mediterranean communities. Ecological Applications, 26 , 796–807. 10.1890/15-0202 27411251
Li, L. , Zhang, Y. , Liu, L. , Wu, J. , Li, S. , Zhang, H. , Zhang, B. , Ding, M. , Wang, Z. , & Paudel, B. (2018). Current challenges in distinguishing climatic and anthropogenic contributions to alpine grassland variation on the Tibetan plateau. Ecology and Evolution, 8 , 5949–5963. 10.1002/ece3.4099 29938105
Li, W. , Li, X. , Zhao, Y. , Zheng, S. , & Bai, Y. (2018). Ecosystem structure, functioning and stability under climate change and grazing in grasslands: Current status and future prospects. Current Opinion in Environmental Sustainability, 33 , 124–135. 10.1016/j.cosust.2018.05.008
Linhart, Y. B. , & Thompson, J. D. (1999). Thyme is of the essence: Biochemical polymorphism and multi‐species deterrence. Evolutionary Ecology Research, 1 , 151–171.
Liu, J. , Feng, C. , Wang, D. , Wang, L. , Wilsey, B. J. , & Zhong, Z. (2015). Impacts of grazing by different large herbivores in grassland depend on plant species diversity. Journal of Applied Ecology, 52 , 1053–1062. 10.1111/1365-2664.12456
Magnano, A. L. , Krug, P. , Casa, V. , & Quintana, R. D. (2019). Changes in vegetation composition and structure following livestock exclusion in a temperate fluvial wetland. Applied Vegetation Science, 22 , 484–493. 10.1111/avsc.12453
Metera, E. , Sakowski, T. , Słoniewski, K. , & Romanowicz, B. (2010). Grazing as a tool to maintain biodiversity of grassland – A review. Animal Science Papers and Reports, 28 , 315–334.
Midolo, G. , Herben, T. , Axmanová, I. , Marcenò, C. , Pätsch, R. , Bruelheide, H. , Karger, D. N. , Aćić, S. , Bergamini, A. , Bergmeier, E. , Biurrun, I. , Bonari, G. , Čarni, A. , Chiarucci, A. , De Sanctis, M. , Demina, O. , Dengler, J. , Dziuba, T. , Fanelli, G. , … Chytrý, M. (2023). Disturbance indicator values for European plants. Global Ecology and Biogeography, 32 , 24–34. 10.1111/geb.13603
Milchunas, D. G. , Sala, O. E. , & Lauenroth, W. K. (1988). A generalized model of the effects of grazing by large herbivores on grassland community structure. The American Naturalist, 132 , 87–106. 10.1086/284839
Mládek, J. , Mládková, P. , Hejcmanová, P. , Dvorský, M. , Pavlu, V. , De Bello, F. , Duchoslav, M. , Hejcman, M. , & Pakeman, R. J. (2013). Plant trait assembly affects superiority of Grazer's foraging strategies in species‐rich grasslands. PLoS One, 8 (7 ), e69800. 10.1371/journal.pone.0069800 23922805
Négyesi, G. (2018). Investigation of spatial and temporal variation of field hedgerows in the shear zone ‐ from the perspective of wind erosion (Hungarian). Tájökológiai Lapok, 16 , 113–128.
Olson, B. E. , & Lacey, J. R. (1994). Sheep: A method for controlling rangeland weeds. Sheep Research Journal, Special Issue, 1994 , 105–112.
Pettit, N. E. , Froend, R. H. , & Ladd, P. G. (1995). Grazing in remnant woodland vegetation: Changes in species composition and life form groups. Journal of Vegetation Science, 6 , 121–130. 10.2307/3236263
Rosenthal, G. , Schrautzer, J. , & Eichberg, C. (2012). Low‐intensity grazing with domestic herbivores: A tool for maintaining and restoring plant diversity in temperate Europe. Tuexenia, 32 , 167–205.
Salguero‐Gómez, R. (2017). Applications of the fast‐slow continuum and reproductive strategy framework of plant life histories. New Phytologist, 213 , 1618–1624. 10.1111/nph.14289 27864957
Sasaki, T. , Okubo, S. , Okayasu, T. , Jamsran, U. , Ohkuro, T. , & Takeuchi, K. (2009). Management applicability of the intermediate disturbance hypothesis across Mongolian rangeland ecosystems. Ecological Applications, 19 , 423–432. 10.1890/08-0144.1 19323200
Schulze, E.‐D. (Ed.). (2019). Plant ecology (second edition). In Springer‐Verlag GmbH, Germany, part of springer nature (p. 926). Berlin. 10.1007/978-3-662-56233-8
Šmilauer, P. , & Lepš, J. (2014). Multivariate analysis of ecological data using Canoco 5 (p. 362). Cambridge University Press. 10.1017/CBO9781139627061
Sonkoly, J. , Tóth, E. , Balogh, N. , Balogh, L. , Bartha, D. , Csendesné Bata, K. , Bátori, Z. , Békefi, N. , Botta‐Dukát, Z. , Bölöni, J. , Csecserits, A. , Csiky, J. , Csontos, P. , Dancza, I. , Deák, B. , Dobolyi, Z. K. , E‐Vojtkó, A. , Gyulai, F. , Hábenczyus, A. A. , … Török, P. (2023). PADAPT 1.0 – The Pannonian dataset of plant traits. Scientific Data, 10 , 742. 10.1038/s41597-023-02619-9 37880224
Su, J. , & Xu, F. (2021). Root, not aboveground litter, controls soil carbon storage under grazing exclusion across grasslands worldwide. Land Degradation & Development, 32 (11 ), 3326–3337. 10.1002/ldr.4008
Süss, K. , Storm, C. , Zimmermann, K. , & Schwabe, A. (2007). The interrelationship between productivity, plant species richness and livestock diet: A question of scale? Applied Vegetation Science, 10 , 169–182. 10.1111/j.1654-109X.2007.tb00515.x
Tonn, B. , Densing, E. M. , Gabler, J. , & Isselstein, J. (2019). Grazing‐induced patchiness, not grazing intensity, drives plant diversity in European low‐input pastures. Journal of Applied Ecology, 56 , 1624–1636. 10.1111/1365-2664.13416
Török, P. , & Dengler, J. (2018). Palearctic grasslands in transition: Overarching patterns and future prospects. In V. R. Squires , J. Dengler , H. Feng , & L. Hua (Eds.), Grasslands of the world: Diversity, management, and conservation (pp. 15–26). CRC Press.
Török, P. , Lindborg, R. , Eldridge, D. , & Pakeman, R. (2024). Grazing effects on vegetation: Biodiversity, management, and restoration. Applied Vegetation Science, 27 , e12794. 10.1111/avsc.12794
Török, P. , Valkó, O. , Deák, B. , Kelemen, A. , Tóth, E. & Tóthmérész, B. (2016). Managing for species composition or diversity? Pastoral and free grazing systems in alkali steppes. Agriculture, Ecosystems & Environment, 234 , 23–30. http://dx.doi.org/10.1016/j.agee.2016.01.010
Tóth, E. , Deák, B. , Valkó, O. , Kelemen, A. , Miglécz, T. , Tóthmérész, B. , & Török, P. (2018). Livestock type is more crucial than grazing intensity: Traditional cattle and sheep grazing in short‐grass steppes. Land Degradation & Development, 29 , 231–239. 10.1002/ldr.2514
WallisDeVries, M. F. , Poschlod, P. , & Willems, J. H. (2002). Challenges for the conservation of calcareous grasslands in northwestern Europe: Integrating the requirements of flora and fauna. Biological Conservation, 104 , 265–273. 10.1016/S0006-3207(01)00191-4
Westoby, M. (1999). The LHS strategy scheme in relation to grazing and fire. In D. Eldridge & D. Freudenberger (Eds.), People and rangelands building the future, vols 1 and 2 (pp. 893–896). VI INT RANGELAND CONGRESS INC.
Xu, F. , Li, J. , Wu, L. , Su, J. , Wang, Y. , Chen, D. , & Bai, Y. (2022). Linking leaf traits to the temporal stability of above‐ and belowground productivity under global change and land use scenarios in a semi‐arid grassland of Inner Mongolia. Science of the Total Environment, 818 , 151858. 10.1016/j.scitotenv.2021.151858 34822882
Yang, H. , Auerswald, K. , Gong, X. , Schnyder, H. , & Bai, Y. (2022). Climate and anthropogenic drivers of changes in abundance of C4 annuals and perennials in grasslands on the Mongolian plateau. Grassland Research, 1 , 131–141.
Zhang, R. , Wang, Z. , Han, G. , Schellenberg, M. P. , Wu, Q. , & Gu, C. (2018). Grazing induced changes in plant diversity is a critical factor controlling grassland productivity in the desert steppe, northern China. Agriculture, Ecosystems & Environment, 265 , 73–83. 10.1016/j.agee.2018.05.014
