Abstract Introduction Research methods and design Results Discussion Conclusions Acknowledgements References About the Author(s) Katlego K. Mashiane Department of Geography, Faculty of Natural and Agricultural Sciences, University of the Free State, Qwaqwa, South Africa Abel Ramoelo Department of Geography, Geoinformatics and Meteorology, Faculty of Natural and Agricultural Sciences, University of Pretoria, Hatfield, South Africa Samuel Adelebu Department of Geography, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein, South Africa Ernest Daemane Department of Geography, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein, South Africa Citation Mashiane, K.K., Ramoelo, A., Adelebu, S. & Daemane, E., 2023, ‘Assessing environmental factors contributing to plant species richness in mountainous mesic grasslands’, Koedoe 65(1), a1732. https://doi.org/10.4102/koedoe.v65i1.1732 Original Research Assessing environmental factors contributing to plant species richness in mountainous mesic grasslands Katlego K. Mashiane, Abel Ramoelo, Samuel Adelebu, Ernest Daemane Received: 27 May 2022; Accepted: 22 June 2023; Published: 04 Aug. 2023 Copyright: © 2023. The Author(s). Licensee: AOSIS. This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Southern African mountain ranges are characterised by rich and diverse plant species thriving in different habitats with pronounced endemic species. However, globally, biodiversity at the species level is deteriorating rapidly because of environmental change leading to habitat degradation and fragmentation. Mountainous grassland communities are particularly vulnerable to rapid ecosystem changes because of their specialised niches and sensitivity to global warming. Understanding the determinants of vegetation is necessary for effective and efficient management. This study aimed to determine significant environmental drivers influencing plant species richness in mountainous mesic grasslands. Topographical variables (slope and elevation) were derived by using a 30 m resolution Digital Elevation Model. Soil variables such as bulk density, silt fragments, pH, coarse fragments, soil organic carbon, sand and nitrogen were acquired from the International Soil Conference and Information Centre (ISCI), and species richness and diversity were derived from vascular plant species inventory data compiled using a field survey. Species richness was influenced by soil bulk density, and the interaction between elevation and soil bulk density; higher species richness was associated with lower bulk density and higher elevations. Similarly, species composition changed as edaphic factors and elevation changed. Conservation implications: Species richness increased with soil bulk density, which decreased with elevation. Fire severity had little effect on species richness and diversity, which may indicate that management actions do not affect species. However, the impact of grazers still needs to be better understood at this time. Moreover, the strong positive relationship between species richness and diversity in this study shows that species richness can be used as a surrogate for diversity and conservation monitoring, especially in mesic highland grasslands. Keywords: mountains; subalpine; edaphic; topography; protected grasslands. Introduction Plant diversity controls processes are pivotal to the functioning of ecosystems (Cardinale et al. 2011). However, species-level biodiversity is deteriorating rapidly because of global environmental change drivers such as climate change, pollution, invasive species and land-use changes (Hooper et al. 2012). It has been stressed that reduced species richness and concomitant changes in species composition could alter ecosystem services (Tilman, Reich & Isbell 2012). In addition, studies showed that reducing plant diversity in grasslands may lead to lower production of above-ground plant biomass (Tilman, Wedin & Knops 1996). Modelling the distribution of mountainous species diversity is essential for identifying priority areas for conservation (Barros et al. 2015), mainly because of the sensitivity of these species to climate change (Picó & Van Groenendael 2007; White et al. 2012). For example, subtropical highland grasslands in Brazil are home to several endemic groups of plants (Iganci et al. 2011), and some taxa are represented by an island-like distribution (Lorenz-Lemke et al. 2010). Furthermore, in mountainous areas of Europe, plant diversity is increased by disturbance regimes such as grazing and soil erosion; as a result, the subalpine grasslands are mosaics comprising varying herbaceous communities (Armas-Herrera et al. 2020). This enhances species composition for grazing while restricting the expansion of unwanted species, that is, shrub encroachment. Understanding of biodiversity indicators and factors driving heterogeneity is crucial for conservation interventions and modelling future distributions in response to global change (Barros et al. 2015; Kuzemko et al. 2016). In addition, the composition of grassland species is controlled by interacting environmental factors such as moisture regime, soil characteristics and other disturbances, including fire. Despite interventions such as policy reform and increasing protected areas, most pressures on biodiversity are intensifying (Butchart et al. 2010; Hooper et al. 2012). As a result, biodiversity loss is accelerating, and thus, this loss of species could be detrimental to the ecosystem and its functions (Butchart et al. 2010; Meyer et al. 2018; Oliver et al. 2015). In addition, richness in grassland species is controlled by parameters that increase heterogeneity such as disturbances, management practices and topography (Zulka et al. 2014). Spatio-temporal studies of grassland vegetation are required to examine species’ response to heterogeneity and ecosystem processes (Lundholm 2009). Besides evidence supporting a positive relationship between environmental heterogeneity and species richness, some studies report negative and monotonous relationships (Lundholm 2009; Yang et al. 2015). Establishing heterogeneity-diversity relationships (HDRs) is fundamental for developing management strategies and abating biodiversity loss (Filibeck et al. 2019). For example, forage nutrition and spatiotemporal vegetation variations because of geographic rainfall gradients influence the occurrence of large mammals in Serengeti National Park, Tanzania (McNaughton 1985, 1988). Soil catena characteristics control species composition in sub-alpine grassland communities of Monte Perdido National Park in Spain (Armas-Herrera et al. 2020). Although potentially susceptible to global environmental change, the vegetation of temperate mountains in Southern Africa is influenced by ecological, physical and anthropogenic processes (Brown & Du Preez 2020). Information and studies on the mountains in southern Africa still need to be documented. Therefore, by intensifying research and monitoring in this region, information on the vegetation and resilience of this ecosystem would become available (Brown & Du Preez 2020). Studies that combine multiple datasets based on ecological traits are encouraged to better understand what determines spatial patterns of biodiversity at explicit scales (Barros et al. 2015). Grassland ecosystems rank among the most threatened ecosystems in South Africa. Thus, conservation plans should be put in place to abate the loss of rare and endangered species. In mesic Afromontane grasslands of southern Africa, vegetation is subjected to intense land-use activities such as grazing by livestock and wildlife (O’Connor 2008). The Golden Gate Highlands National Park (GGHNP) is a protected conservation area experiencing intense grazing by livestock and wild animals, with the potential to change the composition of grass species. The GGHNP is highly variable in topography and climate, which results in a complex mosaic of plant communities (Kay, Bredenkamp & Theron 1993). The existing monitoring programmes for grassland vegetation and data ought to be evaluated and analysed for effective conservation policies. The objective of the study was to determine significant environmental drivers influencing plant species richness in the mountainous mesic grasslands of GGHNP. Research methods and design Study area The study was conducted in GGHNP, situated in the north-eastern part of the Free State province, South Africa (Figure 1). The park covers 32758.35 ha and lies between 28°27ʹ S – 28°37ʹ S and 28°33ʹ E – 28°42ʹ E. The GGHNP is a mountainous grassland located at the foothills of the Drakensberg and forms part of the mesic highveld grassland with marked variation in geology and topography. The following soil types were identified in the study area: shallow rocky soils (Glenrosa and Mispah), deep soil along drainage lines (Oakleaf), well-developed sandy soils (Hutton and Clovelly), as well as Clayey structured soils (Milkwood and Tambakulu) (SANParks 2020). The Park has a rich geology divided into stratigraphic units; four sedimentary units (Tarkastad subgroup, Molteno, Elliot and Clarens formations) and an igneous Drakensberg group that dates 195–210 million years ago. Golden Gate Highlands National Park is characterised by summer rainfall, temperate summers and cold winters. The rainfall season stretches from September to April with a mean annual ranging from 1800 mm to 2000 mm (Kay et al. 1993). The GGHNP lies between 1654 m and 2829 m above sea level (Figure 1) and comprises the following grasslands units: Eastern Free State sandy grasslands (Gm 4), Basotho montane shrubland (Gm 5), Lesotho highveld basalt grassland (Gd 8) and Northern Drakensberg highveld (Gd5) (Mucina & Rutherford 2006). FIGURE 1: A map of the Golden Gate Highlands National Park study area (the park) with names of sampling sites, plots and elevation (metres above sea level [masl]). Sampling design and field data collection The researchers stratified GGHNP according to old crop farming fields that are distributed evenly, these sites have quadrats, were sampled in the study area. Data were collected in March 2019, April 2019 and May 2019 during the peak growing season. At each quadrat, all herbaceous plant species were identified to species level where possible. The 1 m × 1 m was gridded into 100 cells of 1 cm2 each (representing 1% cover) and summed different land types and elevation (Figure 1). Golden Gate Highlands National Park was proclaimed in 1962 with additional proclamation of eight farms on the border of QwaQwa in 1989 (SANParks 2020). Vegetation sampling was undertaken in 30 m × 30 m plots placed purposively within homogeneous patches of grass communities. Each sampling plot was divided into four transects of 30 m. Within each 30 m transect, a nested 1 m × 1 m quadrat was sampled at 10 m intervals. In total 142 plots each with 16 quadrats, species cover was averaged per species from each plot. Ancillary data acquisition and preparation Chemical and physical properties of soil such as soil bulk density (BD), silt fragments (SF), pH, coarse fragments (CFs), soil organic carbon (SOC), sand (SD) and nitrogen of topsoil (15 cm) were downloaded from the International Soil Conference and Information Centre (https://soilgrids.org); the soil data were validated using ground data from across the globe (Hengl et al. 2017). Elevation data and Shuttle Radar Topography Mission (SRTM), Digital Elevation Model (DEM) at 30 m resolution were obtained from US Geological Survey’s EROS data centre (https://earthexplorer.usgs.gov). The slope was derived from a DEM using the Spatial Analyst Tool in ArcGIS 10. 4. Fire severity data were acquired from a study conducted in the study area (Adagbasa, Adelabu & Okello 2018). The study estimated fire severity using the Normalized Burn Ratio (NBR) index by analysing pre and post fire season (April 2017–September 2017) from remote sensing images with burnt and unburnt pixels. To prepare the data for analysis, the raster images were clipped to the extent of the study area and imported into R to extract pixel values representing environmental variables in each plot using the ‘extract’ function. The taxonomic composition of the vegetation and its cover were used to derive species richness and diversity for each plot. The former was determined by averaging the 1 m2 quadrats to obtain mean plot species richness. Total species diversity was calculated using the Shannon–Wiener Index: where pi is the proportion of the species within the sampling unit, these units were averaged to obtain plot level value of species diversity (alpha diversity), which was used as a predictor of species richness. The species richness and diversity were computed using ‘vegan’ and ‘plyr’ found in R statistical software packages (Oksanen 2017) by implementing the diversity and apply function for species diversity and richness, respectively. Data analysis The function ‘stat. desc’ (in pastecs package) in R was used to determine descriptive statistics of all the variables in this study. Exploratory analysis, that is, a correlogram, was used to test for intercorrelations and the relationships between biotic (species richness and diversity), abiotic (BD, SF, sand, soil pH, SOC, CF), fire severity, topographic data (elevation and slope) and an interaction of soil and topo-edaphic variables. Simple multiple linear regression was used to determine the influence of each environmental variable on species richness. To ensure that multicollinearity was low inflation factors (VIFs) were calculated and the variable with the highest VIF was removed. This was repeated until the VIF of all variables in the final model was at acceptable levels (VIF < 5). Hierarchical partitioning was used to estimate the contribution of each explanatory variable in the multivariate linear models because the goal was to identify variables with strong direct influence. Hierarchical partitioning calculates all possible regression subsets and averages model improvements across all hierarchies, thus permitting a ranking of variables by their independent effects. This analysis was calculated using the R package ‘hier.part’. Permutational multivariate analysis of variance (PERMANOVA) was used to determine the influence of environmental predictor variables on species composition among different sites in Golden Gate Highlands National Park (GGHNP). Subsequently, environmental variables were fitted onto an ordination of species composition to determine variables that significantly affected species composition across 9 of 13 sites based on the ample species composition data > 5 species per quadrat. These relationships were displayed using non-metric multidimensional scaling (NMDS). Ethical considerations This article followed all ethical standards for research without direct contact with human or animal subjects. Results Descriptive statistics The maximum number of species recorded in a single plot was 17 and the lowest was 3. The highest record of species diversity was 2.45 while the lowest was 0.64 in GGHNP (Table 1). Out of the ecological variables, species richness had the highest coefficient of variation (38.9%) compared with fire severity (25.9%) and species diversity (24.7%) (Table 1). Among all edaphic variables, SOC showed the highest coefficient of variation (31.6%). Slope (80.4%) had the highest coefficient of variation out of the topographic variables when compared with elevation (5.72%) (Table 1). TABLE 1: Descriptive statistics including species richness, species diversity, soil variables (list them), topographic variables (list them) and fire severity. Relationship between species richness and environmental variables There was a significantly negative relationship between grass species richness and SF, soil organic content and CF; however, there was a positive relationship between species richness and sand. Although significant, the other variables had a weak relationship with species richness. Species diversity was not significantly related to topographical and soil variables (see Figure 2). FIGURE 2: Correlogram for 142 herbaceous samples in the Golden Gate Highlands National Park. The multiple linear regression showed that species richness was positively and significantly associated with species diversity within the biotic-based model (Table 2). However, after removing species diversity, the model with abiotic variables showed soil BD as the only significant variable associated with species richness, albeit negatively. With the topo-edaphic-based model, interaction variables of elevation with BD, CF and sand fragments contributed significantly to variation. Only soil BD was positively associated with species richness. TABLE 2: Herbaceous species richness versus biotic, abiotic and topo-edaphic variables. Hierarchical partitioning showed species diversity as the variable having the highest independent effects followed by edaphic variables (Figure 3). The abiotic-based model exhibited silt fragments as the variable with the highest explanation of variation in species richness. The interaction term of elevation and BD was more pronounced as the variable with the strongest direct influence on species richness. FIGURE 3: Hierarchical partitioning regression analysis for herbaceous species richness of (a) biotic-based (b) abiotic-based and (c) topo-edaphic-based variables. Grass species composition was significantly related to the environmental predictors, that is, BD, SF, sand fragments, pH and elevation. Annashope communities were associated with low sand and pH values. Maluti and Heuwelkop (Figure 4) communities separated from the other sites and thrived on soils with high BD. FIGURE 4: Non-metric multidimensional scaling ordination of Golden Gate Highlands National Park, grass species composition at different sites. Discussion In this study, species richness was not significantly related to any of the topographic variables. Although a positive HDR was expected, these are common but not universal. A metanalysis conducted by Lundholm (2009) found positive HDRs to be mostly lacking. The lack of a relationship may have been because of the low levels of topographic variation across the sites. For instance, the variation in elevation range was very low, this may be a result of our sampling technique of ‘stratified random’, which used old farming lands as a stratum and not elevation. The full mountain range in GGHNP has highly variable and unique topography, which results in diverse habitats (Kay et al. 1993). Therefore, stratification based on elevation maybe important for testing positive HDR of mountainous areas. The variations in mountain topography create micro-environmental conditions that influence habitat heterogeneity and fragmentation (Brown & Du Preez 2020; Gaury & Devi 2017). Consequently, these micro-habitats harbour and shelter a diversity of species serving as essential biodiversity hotspots, providing ecosystem services such as pollination and controlling biological invasion (Zulka et al. 2014). Studies investigating the influence of elevational gradients on species composition would be a boon for conservation management. A model incorporating species diversity and a set of environmental variables performed well in predicting species richness. In addition, species diversity was strongly related to species richness in Spearman’s correlation analysis. The results highlight the importance of incorporating biotic variables in establishing HDRs, which require an understanding of individual species’ response to gradients of spatial environmental heterogeneity factors (Lundholm 2009). For example, local patterns in soil moisture seem to be affected by topography through more mechanisms than other topographically controlled factors and with a strong influence on plant diversity. Lundholm (2009) suggests that species diversity allows for compartmentalisation and complementary resource use because no single species can fully exploit the conditions of the entire ecosystem, a concept known as niche differentiation (Tilman 1999). Therefore, species richness is merely the number of species in an ecosystem, while diversity considers the evenness; thus, considering species diversity and richness as independent and dependent variables helps to understand their influence on the GGHNP ecosystems’ resilience. Soil BD was the only variable negatively affecting species richness within the abiotic-based model and the results of the topo-edaphic-based model showed that elevation and BD jointly affected species richness. This is not surprising because the mountain soil is affected by topography and thus influences local terrestrial plant diversity patterns (Moeslund et al. 2013). Edaphic features such as soil BD can impair or amplify water loss and spatial accumulation patterns (Moeslund et al. 2013). As such, the extent to which portions of a landscape capture or retain water and nutrient resources is an essential determinant of vegetation patterns (Moeslund et al. 2013). This emphasises the role of topographically controlled edaphic features on landscapes in GGHNP. This could indicate that soil with a lower BD can be vital for plant biodiversity in grassland terrestrial plant species. The results suggest that stable landscapes at high elevations harbour niche species diversity. Mid and foot-slope catena are often the recipients of all upslope water and nutrient resources. This causes enrichment that either limits or enhances species richness and diversity, depending on soil quality (McNaughton 1985). Steep slopes are prone to soil erosion, a disturbance could either promote the growth of early-succession species or preclude disturbance-intolerant species in mountainous grassland ecosystems (Armas-Herrera et al. 2020; Filibeck et al. 2019). Species composition was significantly related to the environmental predictors of soil variables and elevation. The ordination showed that species composition turned over along sand, pH and BD gradients. This is consistent with the results of Armas-Herrera et al. (2020), who suggested that grassland vegetation mosaics are maintained by soil–plant feedbacks, especially in subalpine grassland communities where soil disturbances and inevitable succession drive plant co-existence. This finding emphasises the importance of soil drainage in maintaining species diversity in African Mountain vegetation (Brown & Du Preez 2020). A species composition shift will likely occur associated with increasing soil resources such as water and nutrients. As defined by particle size, soil texture, is important for controlling soil–plant relationships, especially seedling germination and establishment. Soil texture especially, silt and clay content, strongly influences the establishment of many plant species (Fair, Peters & Lauenroth 2001). The risk of species extinction because of global environmental change is imminent, propelling the modelling and identification of conservation hotspots (Barros et al. 2015). Subtropical highland grasslands have co-evolved with climate and specific microhabitat conditions (Iganci et al. 2011). Changes in climatic conditions may lead to environmental changes in microhabitats with dire consequences for species diversity, especially in mountainous grasslands, which are sensitive and thus susceptible to environmental changes. The conservation of topographically controlled soils that support and enhance plant diversity can be crucial for maintaining grassland plant communities. In fact, McNaughton (1985) showed that the resilience for grass species to environmental changes will likely depend on their physiological traits in rangelands. Monitoring aspects of biodiversity, that is, richness and diversity as well as species composition are critical for managing ecosystem functions and services for rangelands. These parameters of diversity are crucial for determining carrying capacity and stocking rates to ensure sustainable management of natural resources. Conclusions The research corroborates the fact that drivers of species diversity in grassland communities are multidimensional by nature (Shoemaker et al. 2020). This is evident by the variations in importance and direction of multiple environmental drivers and their influence on species richness. To understand community dynamics, the incorporation of biotic and abiotic factors is beneficial (Shoemaker et al. 2020). The notion that diverse habitats, which are influenced by unique and heterogeneous topography in mountainous biomes, enhance species composition should be investigated. These habitats could be meta-populations fragmented into islands of similar species composition. The research shows that slope and soil interactions are important influencers of species richness in this montane grassland, particularly slope and BD. Soil conservation practices such as revegetation and erosion control blankets on bare soil should therefore be a priority, especially in mesic mountainous grasslands where the soil is prone to erosion and nutrient leaching. Acknowledgements The authors would like to thank the South African National Parks for affording the opportunity to conduct this research. They also acknowledge scientific services intern Mpho Mbedzi, and the staff at the Golden Gate Highlands National Park, especially the rangers who aided immensely in the field data collection. Competing interests The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article. Authors’ contributions All authors conceived and designed the research. K.K.M. conducted the fieldwork, analysed the data and wrote the article. All authors reviewed and edited the manuscript. Funding information Funding for this research was provided by the University of the Free State and SANParks Scientific Service Division. Data availability Data used in the analysis for the results are available upon request. Disclaimer The views and opinions expressed in this article are those of the authors and do not necessarily reflect the official policy or position of any affiliated agency of the authors. References Adagbasa, G.E., Adelabu, S.A. & Okello, T.W., 2018, ‘Spatio-Temporal Assessment of Fire Severity in a Protected and Mountainous Ecosystem’, IGARSS 2018 - 2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 2018, pp. 6572–6575. https://doi.org/10.1109/IGARSS.2018.8518268 Armas-Herrera, C.M., Badía-Villas, D., Mora, J.L. & Gómez, D., 2020, ‘Plant-topsoil relationships underlying subalpine grassland patchiness’, Science of the Total Environment 712, 134483. https://doi.org/10.1016/j.scitotenv.2019.134483 Barros, M.J.F., Silva-Arias, G.A., Fregonezi, J.N., Turchetto-Zolet, A.C., Iganci, J.R.V., Diniz-Filho, J.A.F. et al., 2015, ‘Environmental drivers of diversity in Subtropical Highland Grasslands’, Perspectives in Plant Ecology, Evolution and Systematics 17(5), 360–368. https://doi.org/10.1016/j.ppees.2015.08.001 Brown, L.R. & Du Preez, J., 2020, ‘Alpine vegetation of temperate mountains of Southern Africa’, in M.I. Goldstein & D.A. Dellasala (eds.), Encyclopedia of the world’s biomes, pp. 395–404, Elsevier, Amsterdam. Butchart, S.H., Walpole, M., Collen, B., Van Strien, A., Scharlemann, J.P., Almond, R.E. et al., 2010, ‘Global biodiversity: Indicators of recent declines’, Science 328(May), 1164–1169. Cardinale, B.J., Matulich, K.L., Hooper, D.U., Byrnes, J.E., Duffy, E., Gamfeldt, L. et al., 2011, ‘The functional role of producer diversity in ecosystems’, American Journal of Botany 98(3), 572–592. https://doi.org/10.3732/ajb.1000364 Fair, J.L., Peters, D.P.C. & Lauenroth, W.K., 2001, ‘Response of individual Bouteloua gracilis (Gramineae) plants and tillers to small disturbances’, American Midland Naturalist 145(1), 147–158. https://doi.org/10.1674/0003-0031(2001)145[0147:ROIBGG]2.0.CO;2 Filibeck, G., Sperandii, M.G., Bazzichetto, M., Mancini, L.D., Rossini, F. & Cancellieri, L., 2019, ‘Exploring the drivers of vascular plant richness at very fine spatial scale in sub-Mediterranean limestone grasslands (Central Apennines, Italy)’, Biodiversity and Conservation 28(10), 2701–2725. https://doi.org/10.1007/s10531-019-01788-7 Gaury, P.K. & Devi, R., 2017, ‘Plant species composition and diversity at the Aravalli Mountain Range in Haryana, India’, Journal of Biodiversity 8(1), 34–43. https://doi.org/10.1080/09766901.2017.1336306 Hengl, T., Mendes de Jesus, J., Heuvelink, G.B., Ruiperez Gonzalez, M., Kilibarda, M., Blagotić, A. et al., 2017, ‘SoilGrids250m: Global gridded soil information based on machine learning’, PLoS One 12(2), e0169748. https://doi.org/10.1371/journal.pone.0169748 Hooper, D.U., Adair, E.C., Cardinale, B.J., Byrnes, J.E., Hungate, B.A., Matulich, K.L. et al., 2012, ‘A global synthesis reveals biodiversity loss as a major driver of ecosystem change’, Nature 486(7401), 105–108. https://doi.org/10.1038/nature11118 Iganci, J.R.V., Heiden, G., Miotto, S.T.S. & Pennington, R.T., 2011, ‘Campos de Cima da Serra: The Brazilian Subtropical Highland Grasslands show an unexpected level of plant endemism’, Botanical Journal of the Linnean Society 167(4), 378–393. https://doi.org/10.1111/j.1095-8339.2011.01182.x Kay, C., Bredenkamp, G.J. & Theron, G.K., 1993, ‘The plant communities of the Golden Gate Highlands National Park in the north-eastern Orange Free State’, South African Journal of Botany 59(4), 442–449. https://doi.org/10.1016/S0254-6299(16)30717-7 Kuzemko, A.A., Steinbauer, M.J., Becker, T., Didukh, Y.P., Dolnik, C., Jeschke, M. et al., 2016, ‘Patterns and drivers of phytodiversity in steppe grasslands of Central Podolia (Ukraine)’, Biodiversity and Conservation 25(12), 2233–2250. https://doi.org/10.1007/s10531-016-1060-7 Lorenz-Lemke, A.P., Togni, P.D., Mäder, G., Kriedt, R.A., Stehmann, J.R., Salzano, F.M. et al., 2010, ‘Diversification of plant species in a subtropical region of eastern South American highlands: A phylogeographic perspective on native Petunia (Solanaceae)’, Molecular Ecology 19(23), 5240–5251. https://doi.org/10.1111/j.1365-294X.2010.04871.x Lundholm, J.T., 2009, ‘Plant species diversity and environmental heterogeneity: Spatial scale and competing hypotheses’, Journal of Vegetation Science 20(3), 377–391. https://doi.org/10.1111/j.1654-1103.2009.05577.x McNaughton, S.J., 1985, ‘Ecology of a grazing ecosystem: The Serengeti’, Ecological Monographs 55(3), 259–294. https://doi.org/10.2307/1942578 McNaughton, S.J., 1988, ‘Mineral nutrition and spatial concentrations of African ungulates’, Nature 334(6180), 343–345. https://doi.org/10.1038/334343a0 Mashiane, K., Ramoelo, A. & Adelabu, S., 2023, ‘Diversifying modeling techniques to disentangle the complex patterns of species richness and diversity in the protected Afromontane grasslands’, Biodiversity and Conservation 32, 1423–1436. https://doi.org/10.1007/s10531-023-02560-8 Meyer, S.T., Ptacnik, R., Hillebrand, H., Bessler, H., Buchmann, N., Ebeling, A. et al., 2018, ‘Biodiversity-multifunctionality relationships depend on identity and number of measured functions’, Nature Ecology and Evolution 2(1), 44–49. https://doi.org/10.1038/s41559-017-0391-4 Moeslund, J.E., Arge, L., Bøcher, P.K., Dalgaard, T. & Svenning, J.-C., 2013, ‘Topography as a driver of local terrestrial vascular plant diversity patterns’, Nordic Journal of Botany 31(2), 129–144. https://doi.org/10.1111/j.1756-1051.2013.00082.x Mucina, L. & Rutherford, M.C., 2006, ‘The vegetation of South Africa, Lesotho and Swaziland’, Strelitzia 19(December), 1–30, viewed 23 February 2022, from http://ebooks.cambridge.org/ref/id/CBO9781107415324A009. O’Connor, T., 2008, ‘Influence of land use on phytomass accumulation in Highland Sourveld grassland in the southern Drakensberg, South Africa’, African Journal of Range & Forage Science 25(1), 17–27. https://doi.org/10.2989/AJRFS.2008.25.1.3.381 Oksanen, J., 2017, Vegan: Ecological diversity, R Package Version 2.4-4, p. 11, viewed 10 January 2022, from https://cran.r-project.org/package=vegan. Oliver, T.H., Heard, M.S., Isaac, N.J.B., Roy, D.B., Procter, D., Eigenbrod, F. et al., 2015, ‘Biodiversity and resilience of ecosystem functions’, Trends in Ecology & Evolution 30(11), 673–684. https://doi.org/10.1016/j.tree.2015.08.009 Picó, F.X. & Van Groenendael, J., 2007, ‘Large-scale plant conservation in European semi-natural grasslands: A population genetic perspective’, Diversity and Distributions 13(6), 920–926. https://doi.org/10.1111/j.1472-4642.2007.00349.x SANParks, 2020, Golden Gate Highlands National Park Park Management Plan, South African National Parks. Shoemaker, L.G., Sullivan, L.L., Donohue, I., Cabral, J.S., Williams, R.J., Mayfield, M.M. et al., 2020, ‘Integrating the underlying structure of stochasticity into community ecology’, Ecology 101(2), 1–17. https://doi.org/10.1002/ecy.2922 Tilman, D., Wedin, D. & Knops, J., 1996, ‘Productivity and sustainability in grassland ecosystems’, Nature, 379(3), 718–720. https://doi.org/10.1038/379718a0 Tilman, D., 1999, ‘The ecological consequences of changes in biodiversity: A search for general principles’, Ecology 80(5), 1455–1474. https://doi.org/10.1890/0012-9658(1999)080[1455:TECOCI]2.0.CO;2 Tilman, D., Reich, P.B. & Isbell, F., 2012, ‘Biodiversity impacts ecosystem productivity as much as resources, disturbance, or herbivory’, Proceedings of the National Academy of Sciences 109(26), 10394–10397. https://doi.org/10.1073/pnas.1208240109 White, S.R., Carlyle, C.N., Fraser, L.H. & Cahill, Jr. J.F., 2012, ‘Climate change experiments in temperate grasslands: Synthesis and future directions’, Biology Letters 8(4), 484–487. https://doi.org/10.1098/rsbl.2011.0956 Yang, Z., Liu, X., Zhou, M., Ai, D., Wang, G., Wang, Y. et al., 2015, ‘The effect of environmental heterogeneity on species richness depends on community position along the environmental gradient’, Scientific Reports 5(1), 15723. https://doi.org/10.1038/srep15723 Zulka, K.P., Abensperg-Traun, M., Milasowszky, N., Bieringer, G., Gereben-Krenn, B.A., Holzinger, W. et al., 2014, ‘Species richness in dry grassland patches of eastern Austria: A multi-taxon study on the role of local, landscape and habitat quality variables’, Agriculture, Ecosystems and Environment 182, 25–36. https://doi.org/10.1016/j.agee.2013.11.016