Next Article in Journal
Effects of Water Content and Mesh Size on Tea Bag Decomposition
Next Article in Special Issue
Genus-Physiognomy-Ecosystem (GPE) System for Satellite-Based Classification of Plant Communities
Previous Article in Journal / Special Issue
Response of the Zygopteran Community (Odonata: Insecta) to Change in Environmental Integrity Driven by Urbanization in Eastern Amazonian Streams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plant Diversity in Sardinian Mountain Rangelands: Analysis of Its Relationships with Grazing, Land Management, and Pastoral Value

1
Department of Agricultural Sciences, University of Sassari, Viale Italia 39, 07100 Sassari, Italy
2
CNR (Italian National Research Centre), ISPAAM, Institute for Animal Production System in Mediterranean Environment, Via La Crucca 3, Località Baldinca, 07100 Sassari, Italy
*
Author to whom correspondence should be addressed.
Submission received: 11 January 2021 / Revised: 20 February 2021 / Accepted: 22 February 2021 / Published: 4 March 2021
(This article belongs to the Special Issue Feature Papers of Ecologies 2021)

Abstract

:
In this study, we analyzed the effects of grazing on native and endemic plant diversity, as well as its relationship with pastoral value along a gradient of abiotic and biotic factors and types of land management in a mountainous area of central-eastern Sardinia, Italy. Plant diversity was estimated by conducting a floristic survey within plots. In total, 231 plant species were recorded in 63 plots distributed within the study area, and this total number included 20 endemic species. Species richness was mainly affected by the type of management, soil attributes, altitude, and bioclimate. Pastoral value was strongly affected by nutrient availability and bioclimate. Our results suggest that the cover of endemic species increases with altitude. Finally, in Sardinian rangelands, a negative effect of grazing pressure on endemic species was observed.

1. Introduction

Rangelands cover approximately one-third of the Earth’s land area and at least one billion people depend on these lands to live, both directly through livestock production and indirectly from other resources and ecosystem services [1,2]. They are traditionally characterized by highly heterogeneous natural vegetation communities, which often have a high conservation value, frequently grow in harsh environments [3], are associated with wildlife and/or domestic grazing, and are managed by ecological/traditional rather than agronomic methods [4]. One of the most common forms of use of rangelands with domestic livestock is continuous grazing [5]. Nevertheless, rotational grazing (i.e., moving livestock among paddocks) has become common and has provoked discussion among rangeland ecologists and managers regarding the best management practices [6,7]. Consequently, the management and preservation of rangelands, which are important ecosystems for the maintenance of biodiversity [8], must be supported by a deep knowledge of natural resources and by an accurate level of grazing [9].
Rangelands in the Mediterranean area and on Sardinia are predominantly man-made grasslands, oak-wood grasslands, and oak-dominated woodlands [10], of which the stability and productivity largely depend on the amount of species diversity [11]. Importantly, there is increasing evidence that maintaining low or sustainable levels of grazing livestock is essential to preserve plant species and habitat types linked to pastures [11]. Furthermore, particularly in Sardinia, the presence of wild animals (e.g., mouflons and wild boars), overgrazing, soil-ploughing in natural or semi-natural land, mechanical bush removal, and high frequency of fires can lead to degradation of natural vegetation, particularly where soils are shallow, topography is mountainous, and the climate is thermo-Mediterranean [12]. In addition, the intensive exploitation of rangelands under high stocking rates promotes an increase in unpalatable species (e.g., poisonous, toxic, and spiny species), and a decrease in forage biomass and value, which are common indicators of land degradation [12,13].
The Mediterranean Basin is a renowned global biodiversity area with high plant diversity, sheltering 20% of the world’s total floristic richness [14]. Thus, understanding the drivers that modify plant diversity, particularly in rangelands, and its relationships with grazing and land management is a fundamental question in applied ecology and an important factor for conservation purposes, management practices, and the adoption of suitable strategies for the preservation of plant diversity. The aim of the present study was to evaluate changes in plant diversity (native and endemic species richness) and in pastoral value (based on species pabularity) due to grazing pressure, abiotic and biotic factors, and land management types in a mountain area on a Mediterranean island. More precisely, the present research addresses the following questions: (1) How does plant diversity, plant community composition, and pastoral value respond to the effect of grazing? (2) Which abiotic, biotic, and management variables drive the highest differences in plant community composition?

2. Experimental Section

2.1. Study Area

The research was carried out on Sardinia, Italy, the second largest island in the Mediterranean basin (24,100 km2), in the mountainous area of the Ogliastra region, which is located in central-eastern Sardinia (Figure 1), at an altitude ranging between 290 and 1350 m a.s.l. The climate is typically Mediterranean, and is characterized by two main seasons, a hot-dry and a cold-humid season [15]. The vegetation is composed of meso-Mediterranean silicicolous shrublands, Sardinian holm-oak forests and supra-Mediterranean holm-oak forests, Mediterranean riparian elm forests, Southern and Sicilian Italian Quercus pubescens woods, Mediterranean xeric grasslands, evergreen oak matorral and Juniper matorral, and conifer plantations, according to the Map of the Habitats of Sardinia [16].
The Ogliastra region was selected as a suitable study area due to its rich flora and due to its being predominantly characterized by extensively managed rangelands, which contain agro-forestry land used for grazing various animals (e.g., cattle, donkeys, goats, horses, pigs, and sheep) throughout the year. The investigated rangelands are partly included in territories managed by the Sardinian Forest Agency (FoReSTAS). These rangelands are also partly included within Special Areas of Conservation (SAC) and Special Protection Areas (SPA) according to the Council Directive no. 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora.

2.2. Vegetation Surveys

Vegetation surveys were carried out in 63 plots, located according to a random sampling design, with field surveys from May 2013 to July 2015, when the floristic richness is higher and it is possible to observe higher plant diversity. The surveys were scheduled taking into account the altitudinal gradient, starting from the plots at the lower altitudes. The rangelands in the study area are grazed throughout the year, but more frequently from the beginning of spring until the end of autumn (i.e., for about 200–250 days of grazing per year). Each plot, sized 1 × 1 m, was georeferenced with a portable GPS device. Plant species presence and abundance were visually estimated. Plant abundance was recorded according to the Braun-Blanquet scale (r, +, 1–5).
Plant specimens were collected only if was not possible to identify the species in the field, then pressed, dried, and identified in the laboratory according to Arrigoni’s Flora [15]. Endemic species were identified following Arrigoni’s Flora [15]. This was also the reference used for the diverse biogeographical categories of the endemic species observed.

2.3. Estimation of Response Variables

Plant Diversity and Pastoral Value

Species richness (native and endemic), and Pielou’s evenness (J) were used to express plant diversity. Species richness was the number of plant species observed in each vegetation plot. The function ‘diversity’ finds the most used diversity indices and the R package ’vegan’ was used to calculate the diversity indices, estimated in each plot [17].
The pastoral value (PV) for each plot was estimated using the methodology described by Vacca et al. [18]. The abundance value assigned to each plant species according to the Braun–Blanquet scale was converted into an average cover value per class, as proposed by [19], r = 0.1; + = 2.5; 1 = 7.5; 2 = 17.5; 3 = 37.5; 4 = 62.5; 5 = 87.5, to calculate the specific coverage coefficient (CSR). Therefore, PV was calculated as: PV = 0.2 × (ƩCSPi × SI), where CSPi is the presence specific contribution of a single species and SI is the specific index (from 0—species of no forage interest, to 5—species excellent for quality, palatability, and productivity) [20,21,22]. CSPi = (CRSI/CRStot) × 100, where CSRi is the specific coverage coefficient of eaten species and CSRtot is the specific coverage coefficient of the community. PV ranges from 0 to 100 and estimates the forage potentiality of a pasture area [21].
Finally, we considered for each plot the following functional plant traits: life form, life span, presence of thorns, and pollination type.

2.4. Estimation of Explanatory Variables

Abiotic and Biotic Predictors

A set of predictor variables in GIS format was retrieved from the Sardinian regional geo-database (Sardegna Geoportale, http://www.sardegnageoportale.it, accessed on April 2020). Topographic variables—namely, altitude, aspect, and slope—were determined for each plot from a digital terrain model (DTM, raster 10 × 10 m resolution), using QGIS 3.16.3 software (QGIS Development team, 2018, Open Source Geospatial Foundation). The type of geological substrate was selectedfrom the following: intrusive rocks (granites, granodiorites, and leucogranites), metamorphic rocks (schists, arenaceous, and shales), effusive rocks (basalts), and acid-effusive rocks (andesites and rhyolites) [23]. The bioclimate map of Sardinia was used to identify the main bioclimatic categories found in the 63 plots of the study area (Supplementary Materials, Table S1) [24]. The Special Areas of Conservation (SAC) vector map (https://www.regione.sardegna.it, accessed on April 2020) was used to identify the plots that were inside SACs. The visual assessment of the presence and proximity of trees and shrubs for each plot was carried out and classified as presence and absence. The effects of past fires on the plant community were visually evaluated for each plot and classified as absent or very low, low, medium, or high.
Topsoil samples were collected at a depth of approximately 30 cm for each plot after removing litter, air-dried, and the analyses were conducted on the <2 mm soil fraction after sieving. Soil chemical characteristics were measured in the laboratory for pH, organic carbon (OC) using the Walkley-Black method, nitrogen (N) using the Kjeldahl procedures, and phosphorous (P) as P2O5 using the Olsen method [25].
Grazing pressure was visually assessed (VLU) for each plot and classified as absent or very low, low, medium, or high. In addition to this visual assessment, the spatial variability of livestock units (LU) was evaluated through 49 pasture visual random observation sites (GPS located) within the study area. In each of the 49 sites, all grazing animals within a radius of about 1 km were recorded. Stocking rate by pasture was expressed in livestock units (LU), using the conversion values reported in Regione Autonoma della Sardegna-Programma di Sviluppo Rurale 2007–2013 for the different types of livestock (e.g., cattle value as 1 LU, pigs 0.3 LU, sheep and goats 0.15 LU). Spatialized livestock unit pressure (SLU), used hereafter as a proxy for grazing pressure, was performed between interpolated maps by livestock unit (LU) and ordinary kriging (contains the predicted values with the coordinate covariates) with the R packages ‘gstat’ [26] and ‘kriging’ [27]. Spatial interpolation used the data from the 49 sampling (observation) points to predict values of the observation variable at the 63 plots where data was not collected. We checked the relationship between the visual assessment and the spatial interpolation (ordinary kriging) and retrieved a non-significant association (p-value = 0.0688). Therefore, we decided to consider only the spatial value (SLU), which is a quantitative variable, in all of the following statistical analyses involving grazing pressure.

2.5. Statistical Analysis

Plant community composition at each plot was analyzed using canonical correspondence analysis (CCA) of the vegetation matrix (231 species × 63 plots). The species scores were weighted average plot scores. CCA is a gradient analysis that shows the relationship amongst vegetation patterns under the influence of the abiotic, biotic, and management predictors described above (i.e., analyzing the ecological gradient of vegetation).
Analysis of variance (ANOVA) was performed to test the relationship between the species richness, evenness diversity index, pastoral value, and the number of endemic species, calculated in the 63 sampled plots of the study area.
We used generalized linear models (GLMs) to assess the effects of abiotic, biotic, and management predictors (explanatory variables) on species richness, pastoral value, and the number of endemic species (response variables). Response and explanatory variables are listed in Table 1 and explained in more detail in Table S2 (see Supplementary Materials). A stepwise variable selection procedure was adopted. We used an information-theoretic model selection approach based on Akaike’s information criterion (AIC) to identify the final model among all of the subsets of reduced models. Models considered best had a delta AIC higher than two compared to the next-ranked model.
In addition, we also used GLM to evaluate the effect of the SLU (i.e., a proxy for grazing pressure, considered as an explanatory variable) on functional traits (e.g., life forms, life span, and the presence of thorns) (all considered as response variables) within the study area.
All statistical analyses were performed using R software, version 3.5.2. (R Core Team, 2018, Vienna, Austria). GLMs were performed with the package ‘lme4’ [28]. CCA was run with the ‘cca’ and ‘envfit’ functions of the ‘vegan’ package [17]. Akaike’s information criterion (AIC) was performed with the package ‘My.stepwise’ and function ‘My.stepwise.glm’ [29].

3. Results

3.1. Floristic Richness, Native and Endemic Species

A total of 231 plant species were recorded in the 63 plots, including 20 endemic species (Supplementary Materials, Table S3). These endemics included Sardinian endemics, i.e., endemic species exclusive to the island of Sardinia (8), Sardinian–Corsican endemics (7), and Sardinian–Mediterranean islands endemics (5). Endemic species were found in 45 of the 63 plots, and their number ranged from one to four per plot. Importantly, we observed with significant high frequencies the following endemic species: Armeria sardoa Spreng., Astragalus genargenteus Moris, Genista morisii Colla, Pancratium illyricum L., Plantago sarda C.Presl, and Ptilostemon casabonae (L.) Greuter.
The ANOVA results showed that the number of endemic species per plot was positively correlated with species richness. In addition, we found a negative correlation between endemic species and evenness index (Supplementary Materials, Table S4). Additionally, the GLM procedure highlighted that the percentage cover of endemic species was significantly positively affected by altitude and negatively by SLU and the number of endemic species was positively influenced by both altitude and SLU (p-values > 0.05).
Through the GLM, we found that species richness was positively affected by the abiotic predictors of altitude, soil attributes (OC and N), and iso-bioclimatic type 20. Conversely, soil P2O5 and iso-bioclimatic type 31 resulted in a negatively correlated species richness. Interestingly, species richness showed a highly significant positive correlation with SLU.
In general, recorded species were mainly therophytes (111), hemicryptophytes (62), chamaephytes (18), and geophytes (16), belonging to 47 families, mainly Asteraceae (37), Poaceae (31), and Fabaceae (28). The most common lifespans were annual (111) and perennial (110). Entomophily was the most common pollination type between species (180), followed by anemophily (51) (Supplementary Materials, Figure S1, Table S3).

3.2. Ordination Analysis

The overall results of the CCA ordination plot are shown in Figure 2. A correlation among plant species and the explanatory variables was observed, with 48.73% of the variance explained by the first two canonical axes of CCA. A plot of the first two axes shows two major trends. Axis 1 (eigenvalue = 0.262) was mainly dominated, with a positive correlation, by abiotic factors such as N and the iso-bioclimatic type 37, and negatively correlated with iso-bioclimatic type 28. Soil N is the variable with a long arrow, indicating that it was more strongly correlated with the ordination axis 1 and more closely related to the pattern of community variation than those with short arrows, whereas axis 2 (eigenvalue = 0.224) was distinguished, with a negative correlation, by abiotic factors such as P2O5 and biotic factors such as the proximity of plots to trees. In addition, the plot shows that the plant communities’ composition had a significant negative relationship with grazing pressure (Table 2, Figure 2).

3.3. Pastoral Value

The pastoral value evaluated in the 63 plots of the study area ranged between 18.05 and 46.27 (mean = 29.43 ± 6.23). The results of the GLM are presented in Table 3. The most significant abiotic predictors that influenced PV were soil N and iso-bioclimate “lower mesomediterranean–lower subhumid–euoceanic weak” (isobioclim type 20, Table S1). Pastoral value was not correlated with the altitudinal gradient. The ANOVA results showed that the pastoral value was not significantly correlated with the number of endemic species per plot (p-value = 0.83) (Supplementary Materials, Table S4), or with the cover of endemic species per plot (p-value = 0.902). In fact, the endemic species recorded in the plots had a lower specific index (used to assess the PV) compared to non-endemic species, but often a lower percentage of cover (Supplementary Materials, Figure S2, Table S4).

3.4. Management Type

Of the total number of 63 plots, 38 were under FoReSTAS management and 25 plots were on common land, i.e., they were grazed without any management plan. The GLM procedure highlighted that the land managed by FoReSTAS had a significant negative influence on species richness (Table 3), whereas it did not have a significant effect on the cover of endemic species and pastoral value, or on plant community composition (Figure 2).

3.5. Effect of Grazing on Plant Traits

The grazing values (SLU) ranged between 0.11 and 115.9 (mean = 35.32 ± 23.63). We found, through the GLM, a significant positive effect of SLU on the percentage cover of hemicryptophytes and biannual species per plot (Table 4).

4. Discussion

The range of abiotic, biotic, and management conditions explored in this study allowed us to evaluate the relative importance of these factors in affecting plant diversity (native and endemic species richness) and pastoral value in a Mediterranean mountain rangeland, and reinforced previous studies that showed the effect of grazing pressure on plant community diversity in other parts of Sardinia.

4.1. Plant Diversity, Endemic Species, and Pastoral Value

In the present study, higher species richness was evident in plots at higher altitudes, with a high content of organic carbon and high nutrient availability (soil N). Soil attributes can have a positive relationship with plant diversity due to the creation of optimal growth conditions, and species richness is prone to change due to organic matter, providing a suitable place for the activity of soil microorganisms [30].
Mountains areas in the Mediterranean are characterized by the presence of endemic species [31,32,33], including narrow endemics. As already described for Sardinia [34,35] and other Mediterranean islands [36], and according to our results, altitude was the factor that mainly affected the cover of endemic species per plot, suggesting that the increase in endemic species is highly dependent on the increase in altitude.
We found a clear predominance of therophytes in all of the investigated plots, which is in accordance with previous studies in the Mediterranean region [37]. On the contrary, although the authors of a study [38] in a Mediterranean dwarf shrubland found a positive correlation between pH and species richness concerning annuals, geophytes, and chamaephytes, we found no significant correlation among species richness and pH.
Ordination analysis (CCA) showed that the plant community composition was mainly positively or negatively affected by abiotic (P2O5, N, and iso-bioclimatic type) and biotic factors (grazing pressure and proximity of plots to trees). Some of these variables have a prevailing effect on plant community structure on the island of Sardinia, according to Drissen et al. [39]. Nevertheless, the results of the CCA indicated that the abiotic and biotic factors considered in the study do not explain all of the floristic differences in the study area. Furthermore, our results highlighted that the pastoral value was influenced by different abiotic factors, e.g., nutrient availability (soil N) and the lower mesomediterranean–lower subhumid–euoceanic weak iso-bioclimatic type.
Management type—in our case land managed by FoReSTAS, and therefore subject to some level of regulation of the grazing activities—had a significant importance in conditioning species richness.

4.2. Effect of Grazing

Grazing is known to be a key factor for ecosystem dynamics in Mediterranean island landscapes [40]. Accordingly, in our study and that reported by Noy-Meir and Oron [41], grazing pressure promotes the abundance of geophytes and hemicryptophytes. The abundance of these two lifeforms is related to the production of bulbs, rhizomes, and buds protected in different ways at the ground level, thus acting as an effective barrier against grazing [42]. In addition, recurrent fires can promote these plant traits. However, we did not find any significant relationship in the present study, which can be related to the fact that land managed by FoReSTAS is also actively protected by the risk of wildfires.
In general, moderate grazing pressure optimizes pasture productivity, increases plant diversity [43], and might lead some species to become more dominant and others to become more relegated. On the contrary, overgrazing threatens plant diversity, especially on islands [44]. On the island of Sardinia, a negative effect of high grazing pressure on endemic species was demonstrated by Pisanu et al. [45], which is consistent with our results. However, in a study carried out by Camarda et al. [42], the authors found that 137 (45.4%) endemic species are not affected by grazing, whereas 116 (38.4%) endemic species are promoted by grazing. Grazing reduces the abundance of herbaceous palatable species, whereas it can promote the presence and abundance of species with defensive traits (e.g., spiny), and this may also promote a number of the endemic species, such as Astragalus genargenteus and Genista morisii. Therefore, the effects of grazing pressure can be twofold, both promoting the local richness of endemic species and reducing their cover.

5. Conclusions

Our results highlighted that altitude, soil attributes, and grazing are the main drivers of the plant diversity in the investigated Sardinian mountain rangelands, although the type of management had a notable influence, and pastoral value was highly affected by soil N. Hemicryptophytes were significantly affected by grazing, which must be considered an important component in the maintenance of the landscapes. The effects of altitude and grazing on endemic species should also be taken into account in rangeland management plans in Mediterranean mountain areas.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2673-4133/2/1/9/s1, Figure S1: The pie charts show the number of species distributed according to the lifeforms, lifespan, presence of thorns, pollination type, and the number of endemic species in the 63 sampled plots of the study area. Lifeforms—Ch: Chamaephytes, G: geophytes, H: hemicryptophytes, He: helophytes, NP: nano-phanerophytes, P: phanerophytes, T: therophytes. Lifespan—A: annual, B: biennial, P: perennial. Presence of thorns—Y: yes, N: no. Pollination—A: anemophily, E: entomophily. Endemic species—Y: yes, N: no. Figure S2: Correlation between the presence of endemic species and the specific index for the evaluation of forage pastoral value (from 0—species of no forage interest, to 5—species excellent for quality, palatability, and productivity, according to Delpech, 1960). Table S1: The main iso-bioclimatic types, found in the 63 plots of the study area, extracted from the bioclimate map of Sardinia. The code defines the categories in the original vector map of the bioclimates in Sardinia (Italy). Table S2: Predictors used in the generalized linear models (GLMs), with their detailed descriptions. Table S3: Floristic list and plant traits. Table S4: Analysis of variance (ANOVA) was performed to test the effect of the number of endemic species (explanatory variable) on species richness, evenness diversity index, and pastoral value, calculated in the 63 sampled plots of the study area.

Author Contributions

Conceptualization, G.B., F.S. and G.A.R.; methodology, I.C., G.B., V.L., G.V. and G.A.R.; writing—original draft preparation, V.L. and G.V.; writing—review and editing, V.L., G.B., F.S. and G.A.R.; supervision, I.C., G.B., G.A.R. and S.C.; project administration, G.A.R. and F.S.; funding acquisition, G.A.R. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Autonomous Region of Sardinia (Italy), in the framework of the PO Sardegna FSE 2007-2013, L.R. n. 7/2007 “Promotion of Scientific Research and Innovative Technologies in Sardinia”.

Acknowledgments

The authors gratefully acknowledge Salvatore Nieddu, Daniele Dettori, Piero Saba, and Maddalena Sassu for their continuous and useful technical assistance in field surveys, handling of samples, and laboratory analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ragab, R.; Prudhomme, C. SW—Soil and water. Biosyst. Eng. 2002, 81, 3–34. [Google Scholar] [CrossRef]
  2. Follett, R.F.; Reed, D.A. Soil carbon sequestration in grazing lands: Societal benefits and policy implications. Rangel. Ecol. Manag. 2010, 63, 4–15. [Google Scholar] [CrossRef] [Green Version]
  3. Jouven, M.; Lapeyronie, P.; Moulin, C.H.; Bocquier, F. Rangeland utilization in Mediterranean farming systems. Animal 2010, 4, 1746–1757. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wang, H.-H.; Grant, W.E.; Teague, R. Modeling rangelands as spatially-explicit complex adaptive systems. J. Environ. Manag. 2020, 269, 110762. [Google Scholar] [CrossRef]
  5. Jacobo, E.J.; Rodriguez, A.M.; Bartoloni, N.; Deregibus, V.A. Rotational grazing effects on rangeland vegetation at a farm scale. Rangelands 2006, 59, 249–257. [Google Scholar] [CrossRef] [Green Version]
  6. Briske, D.; Derner, J.; Brown, J.; Fuhlendorf, S.; Teague, W.; Havstad, K.; Gillen, R.; Ash, A.; Willms, W. Rotational grazing on rangelands: Reconciliation of perception and experimental evidence. Rangel. Ecol. Manag. 2008, 61, 3–17. [Google Scholar] [CrossRef] [Green Version]
  7. Teague, R.; Provenza, F.; Kreuter, U.; Steffens, T.; Barnes, M. Multi-paddock grazing on rangelands: Why the perceptual dichotomy between research results and rancher experience? J. Environ. Manag. 2013, 128, 699–717. [Google Scholar] [CrossRef]
  8. Sabatini, S.; Argenti, G.; Staglianò, N.; Bianchetto, E. Il Monitoraggio delle Risorse Prative e Pascolive per la Definizione di Idonee Linee di Gestione Pastorale Sostenibile; Comunicazioni di Ricerca 2001/2; Istituto Sperimentale per l’Assestamento Forestale e l’Alpicoltura (ISAFA): Trento, Italy, 2001; pp. 93–99. [Google Scholar]
  9. Papanastasis, V.P.; Bautista, S.; Chouvardas, D.; Mantzanas, K.; Papadimitriou, M.; Mayor, A.G.; Koukioumi, P.; Papaioannou, A.; Vallejo, R.V. Comparative assessment of goods and services provided by grazing regulation and reforestation in degraded Mediterranean rangelands. Land Degrad. Dev. 2017, 28, 1178–1187. [Google Scholar] [CrossRef] [Green Version]
  10. Bagella, S.; Caria, M.C.; Seddaiu, G.; Leites, L.; Roggero, P.P. Patchy landscapes support more plant diversity and ecosystem services than wood grasslands in Mediterranean silvopastoral agroforestry systems. Agric. Syst. 2020, 185, 102945. [Google Scholar] [CrossRef]
  11. Farris, E.; Secchi, Z.; Rosati, L.; Filigheddu, R. Are all pastures eligible for conservation? A phytosociological survey of the Sardinian–Corsican Province as a basic tool for the Habitats Directive. Plant Biosyst. 2013, 147, 931–946. [Google Scholar] [CrossRef]
  12. Camarda, I. Indicators of Degradation in Agro-ecosystems. In Indicators for Assessing Desertification in the Mediterranean; Os-servatorio Nazionale sulla Desertificazione: Porto Torres, Italy, 1998; pp. 159–163. [Google Scholar]
  13. Merdas, S.; Menad, A.; Mostephaoui, T.; Sakaa, B. Plant community structure and diversity under grazing gradient in arid Mediterranean steppe of Algeria. J. Mater. Environ. Sci. 2017, 8, 4329–4338. [Google Scholar] [CrossRef]
  14. Medail, F.; Quezel, P. Biodiversity hotspots in the Mediterranean Basin: Setting global conservation priorities. Conserv. Biol. 1999, 13, 1510–1513. [Google Scholar] [CrossRef]
  15. Arrigoni, P.V. Flora Dell’Isola di Sardegna; Carlo Delfino Editore: Sassari, Italy, 2006; Volumes 1–6. [Google Scholar]
  16. Camarda, I.; Carta, L.; Brunu, A.; Angelini, P.; Laureti, L. Carta degli Habitat della Regione Sardegna per il Sistema Informativo di Carta della Natura alla Scala 1:50.000; ISPRA-Regione Sardegna Università degli Studi di Sassari: Sassari, Italy, 2010. [Google Scholar]
  17. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R package version 2.5-7. 2019. Available online: https://CRAN.R-project.org/package=vegan (accessed on 15 December 2020).
  18. Vacca, G.; Camarda, I.; Brundu, G.; Sanna, F.; Caredda, S.; Re, G. Comparison of quadrate and transect survey methods to evaluate Pastoral Value (PV) in SE Sardinian rangelands. Options Méditerranéennes Série A Séminaires Méditerranéens 2016, 114, 135–138. [Google Scholar]
  19. Tomaselli, R.I. Introduzione allo Studio della Fitosociologia; Industria Poligrafica Lombarda: Milano, Italy, 1956; 319p. [Google Scholar]
  20. Roggero, P.P.; Bagella, S.; Farina, R. Un archivio dati di Indici Specifici per la valutazione integrata del Valore Pastorale. Rivista di Agronomia 2002, 36, 149–156. [Google Scholar]
  21. Cavallero, A.; Aceto, P.; Gorlier, A.; Lombardi, G.; Lonati, M.; Martinasso, B.; Tagliatori, C. I Tipi Pastorali delle Alpi Piemontesi; Alberto Perdisa Editore: Bologna, Italy, 2007. [Google Scholar]
  22. Bagella, S. Valore Pastorale delle associazioni vegetali: Un esempio di applicazione nell’Appennino Umbro-Marchigiano, Italia. Fitosociologia 2001, 38, 153–165. [Google Scholar]
  23. Aru, A.; Baldaccini, P.; Delogu, G.; Dessena, M.A.; Madrau, S.; Melis, S.T.; Vacca, S.; Vacca, A. Carta dei Suoli della Sardegna; Regione Autonoma della Sardegna: Cagliari, Italy, 1989; Available online: http://www.sardegnaportalesuolo.it/ (accessed on 15 December 2020).
  24. Canu, S.; Rosati, L.; Fiori, M.; Motroni, A.; Filigheddu, R.; Farris, E. Bioclimate map of Sardinia, Italy. J. Maps 2015, 11, 711–718. [Google Scholar] [CrossRef] [Green Version]
  25. Ministero delle Politiche Agricole e Forestali. Decreto Ministeriale 13 Settembre 1999. Approvazione dei Metodi ufficiali di analisi chimica del suolo. (GU Serie Generale n. 248 del 21-10-1999-Suppl. Ordinario n. 185). Rome (Italy). Available online: https://www.gazzettaufficiale.it/eli/id/1999/10/21/099A8497/sg (accessed on 15 December 2020).
  26. Pebesma, E.J. Multivariable geostatistics in S: The gstat package. Comput. Geosci. 2004, 30, 683–691. [Google Scholar] [CrossRef]
  27. Olmedo, O.E. Kriging: Ordinary Kriging. R package Version 1.1. 2014. Available online: https://CRAN.R-project.org/package=kriging (accessed on 15 December 2020).
  28. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Usinglme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  29. International-Harvard Statistical Consulting Company. My Stepwise: Stepwise Variable Selection Procedures for Regression Analysis; R package Version 0.1.0; International-Harvard Statistical Consulting Company: Taipei, Taiwan, 2017. [Google Scholar]
  30. Abbasi-Kesbi, M.; Tataian, M.R.; Tamartash, R.; Fattahi, B. Relationships between soil properties and plant diversity indices (Case Study: Lashgardar Protected Rangeland, Malyer, Iran). J. Rangel. Sci. 2016, 7, 79–89. [Google Scholar]
  31. Hobohm, C.; Janišová, M.; Jansen, J.; Bruchmann, I.; Deppe, U. Biogeography of endemic vascular plants–overview. In Endemism in Vascular Plants; Springer: Dordrecht, The Netherlands, 2014; pp. 85–163. [Google Scholar]
  32. Thompson, J.D.; Lavergne, S.; Affre, L.; Gaudeul, M.; Debussche, M. Ecological differentiation of Mediterranean endemic plants. Taxon 2005, 54, 967–976. [Google Scholar] [CrossRef] [Green Version]
  33. Brundu, G.; Peruzzi, L.; Domina, G.; Bartolucci, F.; Galasso, G.; Peccenini, S.; Raimondo, F.M.; Albano, A.; Alessandrini, A.; Banfi, E.; et al. At the intersection of cultural and natural heritage: Distribution and conservation of the type localities of Italian endemic vascular plants. Biol. Conserv. 2017, 214, 109–118. [Google Scholar] [CrossRef]
  34. Cañadas, E.M.; Fenu, G.; Peñas, J.; Lorite, J.; Mattana, E.; Bacchetta, G. Hotspots within hotspots: Endemic plant richness, environmental drivers, and implications for conservation. Biol. Conserv. 2014, 170, 282–291. [Google Scholar] [CrossRef]
  35. Fois, M.; Fenu, G.; Cañadas, E.M.; Bacchetta, G. Disentangling the influence of environmental and anthropogenic factors on the distribution of endemic vascular plants in Sardinia. PLoS ONE 2017, 12, e0182539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Trigas, P.; Panitsa, M.; Tsiftsis, S. Elevational gradient of vascular plant species richness and endemism in Crete–The effect of post-isolation mountain uplift on a continental island system. PLoS ONE 2013, 8, e59425. [Google Scholar] [CrossRef]
  37. Aronson, J.; Shmida, A. Plant species diversity along a Mediterranean-desert gradient and its correlation with interannual rainfall fluctuations. J. Arid. Environ. 1992, 23, 235–247. [Google Scholar] [CrossRef]
  38. Drissen, T.; Faust, C.; Stadtmann, R.; Treitler, J.T.; Zerbe, S.; Mantilla-Contreras, J. Plant composition and diversity in a semi-natural Mediterranean island landscape: The importance of environmental factors. Plant Biosyst. 2019, 153, 756–766. [Google Scholar] [CrossRef]
  39. Chytrý, M.; Danihelka, J.; Axmanová, I.; Božková, J.; Hettenbergerova, E.; Li, C.; Rozbrojova, Z.; Sekulová, L.; Tichý, L.; Vymazalova, M. Floristic diversity of an Eastern Mediterranean dwarf shrubland: The importance of soil pH. J. Veg. Sci. 2010, 21, 1125–1137. [Google Scholar] [CrossRef]
  40. Médail, F. The specific vulnerability of plant biodiversity and vegetation on Mediterranean islands in the face of global change. Reg. Environ. Chang. 2017, 17, 1775–1790. [Google Scholar] [CrossRef] [Green Version]
  41. Noy-Meir, I.; Oron, T. Effects of grazing on geophytes in Mediterranean vegetation. J. Veg. Sci. 2001, 12, 749–760. [Google Scholar] [CrossRef]
  42. Camarda, I.; Brundu, G.; Carta, L.; Vacca, G.; Brunu, A. Plants and grazing: An evaluation of the effects on Sardinian endemic plants conservation. Flora Mediterranea 2015, 25, 263–275. [Google Scholar]
  43. Naveh, Z.; Whittaker, R.H. Structural and floristic diversity of shrublands and woodlands in Northern Israel and other Mediterranean areas. Vegetatio 1980, 41, 171–190. [Google Scholar] [CrossRef]
  44. Ahlborn, J.; Von Wehrden, H.; Lang, B.; Römermann, C.; Oyunbileg, M.; Oyuntsetseg, B.; Wesche, K. Climate–grazing interactions in Mongolian rangelands: Effects of grazing change along a large-scale environmental gradient. J. Arid. Environ. 2020, 173, 104043. [Google Scholar] [CrossRef]
  45. Pisanu, S.; Farris, E.; Filigheddu, R.; García, M.B. Demographic effects of large, introduced herbivores on a long-lived endemic Plant. Plant Ecol. 2012, 213, 1543–1553. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of the study area. Distribution of the 63 vegetation plots (red dots) in central-eastern Sardinia (plotted on a digital terrain model).
Figure 1. Map of the study area. Distribution of the 63 vegetation plots (red dots) in central-eastern Sardinia (plotted on a digital terrain model).
Ecologies 02 00009 g001
Figure 2. Canonical correspondence analysis (CCA) triplot. The ordination triplot shows the relationship between abiotic, biotic, and management predictors (black arrows), and the plant species (grey crosses) within the different plots (red dots) on the vegetation matrix (231 species × 63 plots). Species codes: All_sub—Allium subhirsutum; Arm_gen—Armeria sardoa subsp. genargentea Bel_bel—Bellium bellidioides; Car_hir—Cardamine hirsuta; Dit_vis—Dittrichia viscosa; Hell_arg—Helleborus argutifolius; Jun_acu—Juncus acutus; Lag_ova—Lagurus ovatus; Pot_rep—Potentilla reptans; Pyr_spi—Pyrus spinosa; Ros_ser—Rosa serafinii; Rub_per—Rubia peregrina; Rub_ulm—Rubus ulmifolius; Umb_rup—Umbilicus rupestris; Vic_lat—Vicia lathyroides; Vio_lim—Viola limbarae.
Figure 2. Canonical correspondence analysis (CCA) triplot. The ordination triplot shows the relationship between abiotic, biotic, and management predictors (black arrows), and the plant species (grey crosses) within the different plots (red dots) on the vegetation matrix (231 species × 63 plots). Species codes: All_sub—Allium subhirsutum; Arm_gen—Armeria sardoa subsp. genargentea Bel_bel—Bellium bellidioides; Car_hir—Cardamine hirsuta; Dit_vis—Dittrichia viscosa; Hell_arg—Helleborus argutifolius; Jun_acu—Juncus acutus; Lag_ova—Lagurus ovatus; Pot_rep—Potentilla reptans; Pyr_spi—Pyrus spinosa; Ros_ser—Rosa serafinii; Rub_per—Rubia peregrina; Rub_ulm—Rubus ulmifolius; Umb_rup—Umbilicus rupestris; Vic_lat—Vicia lathyroides; Vio_lim—Viola limbarae.
Ecologies 02 00009 g002
Table 1. The response and explanatory variables used in the present study with mean, standard deviation (Sd), minimum and maximum values for quantitative variables, and category used for qualitative variables.
Table 1. The response and explanatory variables used in the present study with mean, standard deviation (Sd), minimum and maximum values for quantitative variables, and category used for qualitative variables.
Type of VariableVariableUnitMeanSdMinMax
Response variables
Pastoral value (PV)No.29.436.2318.0546.27
Species richnessNo.42.357.9127.0067.00
Evenness (J)No.0.810.040.710.92
Endemic species richnessNo.2.391.051.004.00
Endemic species—plot cover%5.883.791.0039.5
Life forms—plot cover%Seven types found in the study area
Life span—plot cover%Three types considered
Pollination type—plot cover%Two types considered
Presence of thorns—plot cover%7.168.051.0041.00
Explanatory variables
AbioticpHNo.5.910.544.807.75
OC%4.901.461.327.80
P2O5ppm8.386.151.8635.25
Ng·kg−14.522.060.4110.54
Altitudem a.s.l.10111952801362
Slope°11.638.25041.62
Geological substrateFive types found in the study area
FireClassified into four categories
Iso-bioclimatic typeEight types found in the study area
BioticSpatialized livestock unit pressure (SLU)LU·km−235.3223.630.11115.90
Trees01
Shrubs01
ManagementLand managementCommon land vs. land managed by FoReSTAS
Special Areas of Conservation (SAC)01
N: nitrogen; OC: organic carbon; P2O5: phosphoric anhydride.
Table 2. Canonical correspondence analysis (CCA) model. The table shows the fitting values of abiotic, biotic, and management predictors.
Table 2. Canonical correspondence analysis (CCA) model. The table shows the fitting values of abiotic, biotic, and management predictors.
Predictorsr2Pr (>r)
Trees0.070.01 *
P2O5 (ppm)0.140.01 *
N (g·kg−1)0.020.02 *
SLU (kg·km−2)0.120.02 *
Iso-bioclimatic type0.180.09
pH0.060.14
Altitude (m a.s.l.)0.060.15
Slope (°)0.040.26
Land management0.020.26
Geological substrate0.100.28
SAC0.020.31
OC (%)0.020.60
Fire0.020.88
Shrubs0.000.96
* N: nitrogen; OC: organic carbon; P2O5: phosphoric anhydride; SAC: special areas of conservation; SLU: spatialized livestock unit pressure.
Table 3. Stepwise generalized linear model (GLM). The table shows the GLM analysis of pastoral value, species richness, and the cover of endemic species per plot (response variables) using abiotic, biotic, and management predictors within the 63 sampled plots of the study area. Significance codes: *** 0.001; ** 0.01; * 0.05. Iso-bioclimatic type 20: lower mesomediterranean, lower subhumid, euoceanic weak; Iso-bioclimatic 31: upper mesomediterranean, lower humid, semicontinental weak; N: nitrogen; OC: organic carbon; P2O5: phosphoric anhydride; SLU: spatialized livestock unit pressure.
Table 3. Stepwise generalized linear model (GLM). The table shows the GLM analysis of pastoral value, species richness, and the cover of endemic species per plot (response variables) using abiotic, biotic, and management predictors within the 63 sampled plots of the study area. Significance codes: *** 0.001; ** 0.01; * 0.05. Iso-bioclimatic type 20: lower mesomediterranean, lower subhumid, euoceanic weak; Iso-bioclimatic 31: upper mesomediterranean, lower humid, semicontinental weak; N: nitrogen; OC: organic carbon; P2O5: phosphoric anhydride; SLU: spatialized livestock unit pressure.
Response VariablesPredictorsEstimateStd. Errorz Valuep-Value
Pastoral valueN (g·kg−1)1.030.382.690.01 **
SLU (kg·km−2)0.080.032.270.03 *
Iso-bioclimatic type 200.550.252.230.03 *
Species richnessSLU (kg·km−2)0.010.0013.81<2 × 10−16 ***
Management: FoReSTAS−0.190.02−7.864.00 × 10−15 ***
N (g·kg−1)0.060.017.111.17 × 10−12 ***
Altitude (m a.s.l.)0.000.005.942.93 × 10−9 ***
OC (%)13.143.393.881.00 × 10−4 ***
P2O5 (ppm)−0.010.02−3.301.00 × 10−3 ***
Iso-bioclimatic type 31−0.460.15−3.003.00 × 10−3 **
Iso-bioclimatic type 200.300.132.280.02 *
Endemic speciesAltitude (m a.s.l.)0.000.002.913.00 × 10−3 **
SLU (kg·km−2)−0.010.01−2.350.02 *
Table 4. Generalized linear model (GLM) analysis. The table shows the GLM analysis of the effects of spatialized livestock unit pressure (SLU) on plant life forms (Ch: chamaephytes, G: geophytes, H: hemicryptophytes, He: helophytes, NP: nano-phanerophytes, P: phanerophytes, T: therophytes), life span (A: annual, B: biennial, P: perennial), presence of thorns, and pollination type (A: anemophily, E: entomophily) (all considered response variables). Significance codes: ** 0.01; * 0.05.
Table 4. Generalized linear model (GLM) analysis. The table shows the GLM analysis of the effects of spatialized livestock unit pressure (SLU) on plant life forms (Ch: chamaephytes, G: geophytes, H: hemicryptophytes, He: helophytes, NP: nano-phanerophytes, P: phanerophytes, T: therophytes), life span (A: annual, B: biennial, P: perennial), presence of thorns, and pollination type (A: anemophily, E: entomophily) (all considered response variables). Significance codes: ** 0.01; * 0.05.
Response VariablesEstimateStd. Errorz Valuep-Value
Chamaephytes0.020.100.170.87
Geophytes0.010.040.180.86
Hemicryptophytes0.260.132.010.05 *
Helophytes−0.020.01−1.590.12
Nanophanerophytes−0.000.05−0.020.98
Phanerophytes−0.020.07−0.240.81
Therophytes0.170.200.860.39
Annual0.170.200.860.39
Biennial0.080.032.680.01 **
Perennial0.180.190.950.35
Thorns0.060.090.640.53
Anemophily0.230.161.490.14
Entomophily0.190.210.910.37
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lozano, V.; Vacca, G.; Camarda, I.; Brundu, G.; Sanna, F.; Caredda, S.; Re, G.A. Plant Diversity in Sardinian Mountain Rangelands: Analysis of Its Relationships with Grazing, Land Management, and Pastoral Value. Ecologies 2021, 2, 164-174. https://0-doi-org.brum.beds.ac.uk/10.3390/ecologies2010009

AMA Style

Lozano V, Vacca G, Camarda I, Brundu G, Sanna F, Caredda S, Re GA. Plant Diversity in Sardinian Mountain Rangelands: Analysis of Its Relationships with Grazing, Land Management, and Pastoral Value. Ecologies. 2021; 2(1):164-174. https://0-doi-org.brum.beds.ac.uk/10.3390/ecologies2010009

Chicago/Turabian Style

Lozano, Vanessa, Gabriella Vacca, Ignazio Camarda, Giuseppe Brundu, Federico Sanna, Salvatore Caredda, and Giovanni Antonio Re. 2021. "Plant Diversity in Sardinian Mountain Rangelands: Analysis of Its Relationships with Grazing, Land Management, and Pastoral Value" Ecologies 2, no. 1: 164-174. https://0-doi-org.brum.beds.ac.uk/10.3390/ecologies2010009

Article Metrics

Back to TopTop