Herbivore species and grazing intensity regulate community composition and an encroaching woody plant in semi-arid rangeland

Grazing by livestock can influence ecosystems in various ways, including altering plant communities, influencing woody plant encroachment, and determining livestock productivity. Evaluating long term effects of grazing on plant composition is valuable not only to understand herbivory on rangelands but to be able to address the primary factors that can threaten long term livestock productivity. We examined plant species composition and woody plant encroachment 45 years after the initiation of differing grazing treatments within a semiarid savanna of the southern Great Plains, USA. Grazing treatments varied in herbivore type (domestic cattle, sheep, and goats vs. goats only) and grazing intensity (heavy, moderate, and no-herbivory). All individual trees of Juniperus ashei Buchholz, the encroaching woody plant of the area, were removed prior to treatment initiation. Moderate and heavy grazing by a combination of species resulted in similar plant communities, while a history of heavy browsing by goats only and no-herbivory resulted in more distinct communities. Cover of J. ashei did not differ between mixed grazing and no-herbivory treatments, indicating that grazing was not responsible for woody plant encroachment. J. ashei cover within the browsed treatment was a third less compared to other treatments; compositional differences within this treatment are possibly due to reduced cover of woody vegetation. Declines in livestock productivity of the area are likely related to compositional changes resulting from increased woody plants. Livestock production within this semi-arid rangeland is likely unsustainable without management of woody plant encroachment, as communities tend to a closed canopy woodland.     

Response of plant species and life form diversity to variable fire histories and biomass in the jarrah forest of south‐west Australia

Frequent fires reduce the abundance of woody plant species and favour herbaceous species. Plant species richness also tends to increase with decreasing vegetation biomass and cover due to reduced competition for light. We assessed the influence of variable fire histories and site biomass on the following diversity measures: woody and herbaceous species richness, overall species richness and evenness, and life form evenness (i.e. the relative abundance or dominance among six herbaceous and six woody plant life forms), across 16 mixed jarrah (Eucalyptus marginata) and marri (Corymbia calophylla) forest stands in south‐west Australia. Fire frequency was defined as the total number of fires over a 30‐year period. Overall species richness and species evenness did not vary with fire frequency or biomass. However, there were more herbaceous species (particularly rushes, geophytes and herbs) where there were fewer shrubs and low biomass, suggesting that more herbaceous species coexist where dominance by shrubs is low. Frequently burnt plots also had lower number and abundance of shrub species. Life form evenness was also higher at both high fire frequency and low biomass sites. These results suggest that the impact of fire frequency and biomass on vegetation composition is mediated by local interactions among different life forms rather than among individual species. Our results demonstrate that measuring the variation in the relative diversity of different woody and herbaceous life forms is crucial to understanding the compositional response of forests and other structurally complex vegetation communities to changes in disturbance regime such as increased fire frequency.     

Climate and vegetation structure determine plant diversity in Quercus ilex woodlands along an aridity and human-use gradient in Northern Algeria

We studied the influence of environmental factors relating to climate, soil and vegetation cover on total species richness, species richness of different life-forms and species composition of plant communities occurring in Quercus ilex woodlands, across a 450-km long transect in Northern Algeria constituting a gradient of aridity and human use. We sampled vegetation and collected environmental data in 81 10 m × 10 m plots in five zones representing the largest Q. ilex woodlands throughout the study area, analysing them within an a priori hypothesis framework with the use of Path Analysis. Changes in plant diversity were mainly influenced by environmental factors related to precipitation and temperature regimes, as well as by total plant cover. In particular, changes in species composition were determined by factors associated with the temperature regime through their influence on both woody and annual herbaceous plant richness, and by factors related to the precipitation regime through their influence on perennial herbaceous plant richness, likely due to the differential tolerances of these functional groups to cold and water stress. Our results emphasize the importance of differences in environmental adaptability of the most important life-forms with regard to explaining compositional change (beta diversity) along aridity gradients, and the mediator role of total plant cover in relation to the effects of soil conditions on plant diversity.     

Effects of plant species richness and evenness on soil microbial community diversity and function

Understanding the links between plant diversity and soil communities is critical to disentangling the mechanisms by which plant communities modulate ecosystem function. Experimental plant communities varying in species richness, evenness, and density were established using a response surface design and soil community properties including bacterial and archaeal abundance, richness, and evenness were measured. The potential to perform a representative soil ecosystem function, oxidation of ammonium to nitrite, was measured via archaeal and bacterial amoA genes. Structural equation modeling was used to explore the direct and indirect effects of the plant community on soil diversity and potential function. Plant communities influenced archaea and bacteria via different pathways. Species richness and evenness had significant direct effects on soil microbial community structure, but the mechanisms driving these effects did not include either root biomass or the pools of carbon and nitrogen available to the soil microbial community. Species richness had direct positive effects on archaeal amoA prevalence, but only indirect impacts on bacterial communities through modulation of plant evenness. Increased plant evenness increased bacterial abundance which in turn increased bacterial amoA abundance. These results suggest that plant community evenness may have a strong impact on some aspects of soil ecosystem function. We show that a more even plant community increased bacterial abundance, which then increased the potential for bacterial nitrification. A more even plant community also increased total dissolved nitrogen in the soil, which decreased the potential for archaeal nitrification. The role of plant evenness in structuring the soil community suggests mechanisms including complementarity in root exudate profiles or root foraging patterns.     

Complex facilitation and competition in a temperate grassland: loss of plant diversity and elevated CO2 have divergent and opposite effects on oak establishment

Encroachment of woody vegetation into grasslands is a widespread phenomenon that alters plant community composition and ecosystem function. Woody encroachment is often the result of fire suppression, but it may also be related to changes in resource availability associated with global environmental change. We tested the relative strength of three important global change factors (CO₂ enrichment, nitrogen deposition, and loss of herbaceous plant diversity) on the first 3 years of bur oak (Quercus macrocarpa) seedling performance in a field experiment in central Minnesota, USA. We found that loss of plant diversity decreased initial oak survival but increased overall oak growth. Conversely, elevated CO₂ increased initial oak seedling survival and reduced overall growth, especially at low levels of diversity. Nitrogen deposition surprisingly had no net effect on survival or growth. The magnitude of these effects indicates that long-term woody encroachment trends may be most strongly associated with those few individuals that survive, but grow much larger in lower diversity patches. Further, while the CO₂ results and the species richness results appear to describe opposing trends, this is due only to the fact that the natural drivers are moving in opposite directions (decreasing species richness and increasing CO₂). Interestingly, the mechanisms that underlie both patterns are very similar, increased CO₂ and increased species richness both increase herbaceous biomass which (1) increases belowground competition for resources and (2) increases facilitation of early plant survival under a more diverse plant canopy; in other words, both competition and facilitation help determine community composition in these grasslands.     

Biodiversity of semi‐arid Mediterranean grasslands: Impact of grazing and afforestation

Question: What is the impact of grazing and/or afforestation on grassland diversity, species composition and cover parameters? Location: Semi‐arid Mediterranean grasslands of Jordan. Methods: Vegetation, litter, bare soil and rock cover were compared among four management types – free grazing and protected from grazing with three levels of tree cover. Species composition, plant cover, species richness and evenness were used to evaluate differences in vegetation among management types. Species composition differences among management types were also investigated. Results: Semi‐arid Mediterranean grasslands harbour appreciable levels of plant biodiversity. Grazing did not affect plant diversity, indicating the high resilience against and adaptation to grazing; however,grazing affected species composition and cover parameters. Afforestation seems to protect soil through higher litter cover but its impact on plant biodiversity was negative and markedly affected species composition. Conclusions: Neither protection from grazing or massive afforestation alone are sufficient for conserving biodiversity in this system. A management model is suggested where the landscape should be maintained as a mosaic of four management types: complete protection from grazing, grazing rotation, planting sparse trees in eroded areas and revegetating degraded areas using native, herbaceous and grazing tolerant species.     

Positive versus negative environmental impacts of tree encroachment in South Africa

Woody plant encroachment in grasslands is a worldwide phenomenon. Despite many studies, the consequences of woody plant encroachment on sub-canopy vegetation and soil properties are still unclear. To better understand the impacts of trees on grassland properties we examined the following questions using a mountainous sub-tropical grassland of South Africa encroached by an indigenous tree, Acacia sieberiana as a case study: (1) Do trees increase sub-canopy herbaceous diversity, quality and biomass and soil nitrogen content? (2) Do large trees have a stronger effect than medium-sized trees on grass and soil properties? (3) Does the impact of trees change with the presence of livestock and position of trees in a catena? We studied grass and non-graminoid species diversity and biomass, grass quality and soil properties during the wet season of 2009. Nitrogen in grass leaves, soil cation exchange capacity and calcium and magnesium ion concentrations in the soil increased under tall Acacia versus open areas. Medium-sized Acacia decreased the gross energy content, digestibility and neutral detergent fibre of grasses but increased the species richness of non-graminoids. Tall and medium Acacia trees were associated with the presence of Senecio inaequidens, an indigenous species that is toxic to horses and cattle. The presence of livestock resulted in a decrease in herbaceous root biomass and an increase in soil carbon and leaf biomass of grass under Acacia. Tree position in the catena did not modify the impact of trees on the herbaceous layer and soil properties. For management of livestock we recommend retaining tall Acacia trees and partially removing medium-sized Acacia trees because the latter had negative effects on grass quality.     

Encroachment of woody plants and its impact on pastoral livestock production in the Borana lowlands, southern Oromia, Ethiopia

Encroachment of woody plants has been among the major threats to the livelihoods of Borana pastoralists and their ecosystem. An approach that integrated vegetation survey and pastoralists’ perception was followed to study the impacts of encroachment of woody plants in the Borana lowlands, Ethiopia. Density of woody species was determined in 192 plots of 500 m². Canopy cover of woody plants was estimated in 123 quadrates of 400 m². Pastoralists’ perception was assessed through group discussions and a semi‐structured questionnaire. Results showed that plant density was 3014 woody plants ha^?1. Cover of woody plants was 52%, indicating an increasing trend from ≤40% cover reported in the early 1990s. It was concluded that the increase of woody plants density and cover has crossed the critical threshold and has entered into the encroached condition. Principal components analysis (PCA) and redundancy analysis (RDA) also showed that woody plants were negatively correlated with herbaceous biomass. Commiphora africana, Acacia melliphera, A. drepanolobium, A. brevispica and Lannea rivae were among the dominant encroachers. RDA revealed that soil nutrients were positively correlated with woody plants density and cover. The pastoralists perceived that encroachment of woody plants had decreased the production of their grazingland. A ban on fire was perceived as the major factor that caused encroachment of woody plants. Re‐utilization of fire and strengthening of traditional rangeland management strategies are recommended.     

Trade‐offs between savanna woody plant diversity and carbon storage in the Brazilian Cerrado

Incentivizing carbon storage can be a win‐win pathway to conserving biodiversity and mitigating climate change. In savannas, however, the situation is more complex. Promoting carbon storage through woody encroachment may reduce plant diversity of savanna endemics, even as the diversity of encroaching forest species increases. This trade‐off has important implications for the management of biodiversity and carbon in savanna habitats, but has rarely been evaluated empirically. We quantified the nature of carbon‐diversity relationships in the Brazilian Cerrado by analyzing how woody plant species richness changed with carbon storage in 206 sites across the 2.2 million km² region at two spatial scales. We show that total woody plant species diversity increases with carbon storage, as expected, but that the richness of endemic savanna woody plant species declines with carbon storage both at the local scale, as woody biomass accumulates within plots, and at the landscape scale, as forest replaces savanna. The sharpest trade‐offs between carbon storage and savanna diversity occurred at the early stages of carbon accumulation at the local scale but the final stages of forest encroachment at the landscape scale. Furthermore, the loss of savanna species quickens in the final stages of forest encroachment, and beyond a point, savanna species losses outpace forest species gains with increasing carbon accumulation. Our results suggest that although woody encroachment in savanna ecosystems may provide substantial carbon benefits, it comes at the rapidly accruing cost of woody plant species adapted to the open savanna environment. Moreover, the dependence of carbon‐diversity trade‐offs on the amount of savanna area remaining requires land managers to carefully consider local conditions. Widespread woody encroachment in both Australian and African savannas and grasslands may present similar threats to biodiversity.     

Intraspecific variation overrides origin effects in impacts of litter-derived secondary compounds on larval amphibians

Vertebrate herbivores can be key determinants of grassland plant species richness, although the magnitude of their effects can largely depend on ecosystem and herbivore characteristics. It has been demonstrated that the combined effect of primary productivity and body size is critical when assessing the impact of herbivores on plant richness of perennial-dominated grasslands; however, the interaction of site productivity and herbivore size as determinants of plant richness in annual-dominated pastures remains unknown. We experimentally partitioned primary productivity and herbivore body size (sheep and wild rabbits) to study the effect of herbivores on the plant species richness of a Mediterranean semiarid annual plant community in central Spain over six years. We also analyzed the effect of grazing and productivity on the evenness and species composition of the plant community, and green cover, litter, and plant height. We found that plant richness was higher where the large herbivore was present at high-productivity sites but barely changed at low productivity. The small herbivore did not affect species richness at either productivity site despite its large effects on species composition. We propose that adaptations to resource scarcity and herbivory prevented plant richness changes at low-productivity sites, whereas litter accumulation in the absence of herbivores decreased plant richness at high productivity. Our results are consistent with predictions arising from a long history of grazing and highlight the importance of both large and small herbivores to the maintenance of plant diversity of Mediterranean annual-dominated pastures.     

Assessment of woody species encroachment in the grasslands of Nechisar National Park,Ethiopia

Ecological survey was executed to assess woody species encroachment into the grassland plain of Nechisar National Park (NNP). Forty‐one woody species were recorded. Dichrostachys cinerea Wight & Arn., Acacia mellifera (Vahl) Benth., Acacia nilotica (L) Willd., Acacia senegal (L.) Willd., Acacia seyal Del. and Acacia tortilis (Forssk.) Hayne were among the major encroaching woody species. The majority of the woody species were found to be highly aggregated in their pattern of distribution, while only few species showed some degree of randomness. The mean woody species density was ca. 1995 woody plants ha^?1. Mean cover of woody, grass, unpalatable forbs and total herbaceous species were 31%, 58%, 68% and 121%, respectively. The woody species density and cover, unpalatable forbs and bare land cover were significantly higher in the highly grazed and fire‐suppressed part of the grassland plain. Pearson correlation coefficient matrix indicated that woody species cover and density were negatively correlated with total herbaceous and grass cover. The high woody, unpalatable forbs and bare land cover indicated the progressively increasing perennial grass species diversity deterioration in the grass plain of the Park. Decline in the grassland condition, unless reversed, will jeopardize the biological diversity as well as the aesthetic value of the NNP.     

The Importance of Phytogenic Mounds (Nebkhas) for Restoration of Arid Degraded Rangelands in Northern Sinai

Abstract Natural accumulation of wind‐borne sediments within or around the canopies of plants plays an important role in the ecological and evolutionary dynamics of many coastal and desert ecosystems. The formation of such phytogenic mounds (nebkhas) creates patches that can strongly influence the spatial distribution of plant and soil resources. In land restoration of arid and semiarid environments it is important to study the potential role of such biological patchiness that may provide sites for coexistence of species with different life and growth forms. Our main objective was to test whether the nebkhas of a leguminous shrub, Retama raetam (white broom), promote restoration of herbaceous vegetation and soil in the degraded rangelands of northern Sinai. Vegetation and microclimatic and edaphic characteristics within the nebkhas, as well as within internebkha spaces, were compared for ungrazed and grazed sites. Abundance and richness of herbaceous plants were positively related to nebkha area, which explained more of the variance of abundance and richness in the grazed site than in the ungrazed one. Protection from grazing, especially on nebkhas, was associated with an increase in abundance and richness of herbaceous plants, improved soil microclimate, and increased soil fine particles and nutrient concentrations. The results suggest that management (in casu protection from grazing) of nebkhas of woody perennial shrubs changes rangeland conditions and improves the resource regulatory processes. Furthermore, nebkhas of unpalatable plants have the potential to preserve plant diversity in overgrazed plant communities, because they are effective in capturing and retaining water, soil materials, and propagules within and from nearby areas, resources that would otherwise be lost.     

Effects on soil,microclimate and vegetation of the native-invasive <Emphasis Type="Italic">Retama monosperma</Emphasis> (L.) in coastal dunes

The aim of this study was to analyse and quantify the effects of the canopy of the native-invasive N-fixer woody shrub Retama monosperma in the dune ecosystem, affecting the structure and function of the dune environment as well as plant community, in the context of the facilitation mechanism. Air temperature and relative humidity; soil pH, electric conductivity, organic matter (OM) and nutrient content; above and below-ground vegetation biomass, litter mass, species richness and Shannon diversity were determined and compared from sampling plots below the R. monosperma canopy and in canopy gaps within a coastal dune system in SW Spain. The relationships between soil OM and nutrient contents and above and below-ground vegetation biomass, litter mass, species richness and Shannon diversity were also assessed. A predominance of positive interactions was confirmed. The canopy of R. monosperma ameliorated temperature extremes beneath, and soil OM and nutrient concentrations were increased by 188–466%, compared to those found in gaps. Plant biomass increased by 442% beneath the canopy and was composed almost exclusively of herbaceous annuals. Plant diversity was not affected. Plant communities were clearly structured as fertility islands, distributed in an environmentally stressful dune matrix characterized by scarce vegetation cover and low biomass.     

Large-scale manipulation of plant litter and fertilizer in a managed successional temperate grassland

Plant litter may play an important role in herbaceous plant communities by limiting primary production and influencing plant species richness. However, it is not known how the effect of litter interacts with fertilization. We tested for the role of litter and fertilization in a large-scale experiment to investigate effects on diversity and biomass of plant species, growth forms, native vs. non-native groups, and abiotic ecosystem components (e.g., soil moisture, PAR). We manipulated plant litter (removed vs. left in situ) and nutrient availability (NPK-fertilized vs. unfertilized) for 4 years in 314-m² plots, replicated six times, in an old-field grassland. While many of our species-level results supported previously published studies and theory, our plant group results generally did not. Specifically, grass species richness and forb biomass was not affected by either fertilization or plant litter. Moreover, plant litter removal significantly increased non-native plant species richness. Relative to native plant species, all of our experimental manipulations significantly increased both the biomass and the species richness of non-native plant species. Thus, this grassland system was sensitive to management treatments through the facilitation of non-native plant species. We coupled biotic and abiotic components within a nonmetric multidimensional scaling (NMS) analysis to investigate treatment effects, which revealed that specific treatments altered ecosystem development. These results suggest that fertilization and plant litter may have larger impacts on plant communities and on ecosystem properties than previously understood, underscoring the need for larger-scale and longer-term experiments.     

The role of area enclosures and fallow age in the restoration of plant diversity in northern Ethiopia

This study investigated the roles played by area enclosures and fallow age in the restoration of plant species richness and soil seed bank species richness in degraded mountain rangelands in northern Ethiopia. Management types (enclosures versus grazed) influenced woody and herbaceous species richness, while fallow age showed no effect on the woody species. Management, age and the doubling of fallow age influenced the herbaceous species richness and species diversity. Management showed no effect on soil seed bank species richness. Fallow age and the doubling of fallow time also showed no influence on the soil seed bank of grass species, but they were influential on the forbs species soil seed bank. The trends for restoration of plant species richness and diversity and grass seed bank in response to fallow age were positive‐linear, but they declined when the fallow ages were doubled. The exception was the forbs seed bank showed linear trends when age of restoration was doubled. The data suggest that the restoration of degraded rangelands in the high mountain zones of northern Ethiopia was still in the weedy succession stages. Long‐term monitoring will be required to gain an informed understanding of the roles played by area enclosures and fallow age in the restoration of plant biodiversity.     

Effects of termite mounds on composition,functional types and traits of plant communities in Pendjari Biosphere Reserve (Benin,West Africa)

Understanding the role of termite mounds in biodiversity and ecosystem functioning is a priority for the management of tropical terrestrial protected areas dominated by savannahs. This study aimed to assess the effects of termite mounds on the diversity of plant functional types (PFTs) and herbaceous’ net aboveground primary productivity (NAPP) in plant communities (PCs) of the Pendjari Biosphere Reserve. PCs were identified through canonical correspondence analysis performed on 96 phytosociological ‘relevés’ realized in plots of 900 m². PFTs’ diversity was compared between savannahs and mounds’ plots using generalized linear models. In each plot, 7 m² subplots were harvested and NAPP was determined. Linear mixed models were performed to assess change in herbaceous NAPP regarding species richness, graminoids’ richness, specific leaf area and termite mounds. There is no specific plant community related to mounds. However, the occurrence of termite mounds induced an increase of woody and forbs diversity while the diversity of legumes and graminoids decreased. These diversity patterns led to decreasing of PCs’ NAPP. This study confirms that termite‐induced resource heterogeneity supports niche differentiation theory and increased savannah encroachment by woody species.     

Understory response to varying fire frequencies after 20 years of prescribed burning in an upland oak forest

Ecosystems in the eastern United States that were shaped by fire over thousands of years of anthropogenic burning recently have been subjected to fire suppression resulting in significant changes in vegetation composition and structure and encroachment by invasive species. Renewed interest in use of fire to manage such ecosystems will require knowledge of effects of fire regime on vegetation. We studied the effects of one aspect of the fire regime, fire frequency, on biomass, cover and diversity of understory vegetation in upland oak forests prescribe-burned for 20 years at different frequencies ranging from zero to five fires per decade. Overstory canopy closure ranged from 88 to 96% and was not affected by fire frequency indicating high tolerance of large trees for even the most frequent burning. Understory species richness and cover was dominated by woody reproduction followed in descending order by forbs, C3 graminoids, C4 grasses, and legumes. Woody plant understory cover did not change with fire frequency and increased 30% from one to three years after a burn. Both forbs and C3 graminoids showed a linear increase in species richness and cover as fire frequency increased. In contrast, C4 grasses and legumes did not show a response to fire frequency. The reduction of litter by fire may have encouraged regeneration of herbaceous plants and helped explain the positive response of forbs and C3 graminoids to increasing fire frequency. Our results showed that herbaceous biomass, cover, and diversity can be managed with long-term prescribed fire under the closed canopy of upland oak forests.     

Woody overstorey effects on soil carbon and nitrogen pools in South African savanna

Abstract Woody plant encroachment in savannas may alter carbon (C) and nitrogen (N) pools over the long‐term, which could have regional or global biogeochemical implications given the widespread encroachment observed in the vast savanna biome. Soil and litter %C and %N were surveyed across four soil types in two encroached, semi‐arid savanna landscapes in northern South Africa. Litter at sampling points with a woody component had a higher %N and lower C : N ratio than litter at solely herbaceous points. Severely encroached areas had lower C : N ratios throughout the soil profile than less encroached areas. Soil %C and %N were highly influenced by soil texture but were also influenced by the presence of a woody overstorey, which increased surface soil %C on three soil types but decreased it on the most heavily encroached soil type. Soil C sequestration may initially increase with bush encroachment but then decline if bush densities become so high as to inhibit understorey grass growth.     

Analysis of the endophytic lifestyle and plant growth promotion of Burkholderia terricola ZR2-12

Nitrogen (N) enrichment has been pinpointed as a main driver for biodiversity change. Most of our knowledge of effects of increased N availability on ecosystems comes from northern Europe and America. Most other ecosystem types have been neglected. In contribution to filling this gap, our study examined the short-term effects of N enrichment in a N-manipulation (doses and forms) field study of a severely nutrient-limited Mediterranean ecosystem located in a Natura 2000 site in Portugal. Our aims were to (a) understand the effects of N enrichment on plant diversity, and to (b) link N-driven plant community changes with changes in soil inorganic N availability. In general, the standing plant community responded to short-term N enrichment with increased richness and evenness. Changes in the plant community occurred through changes in species composition and cover, and were correlated with soil N, and N and phosphorus availability. Fertilization with 80 kg NH₄NO₃ ha^?1 y^?1 was the treatment which changed plant composition the most, while geophytes, hemicryptophytes and therophytes were the biological types more responsive to N enrichment. Dittrichia viscosa was the only species that responded significantly to increased N, i.e., its cover decreased in control plots, but increased in fertilized plots, suggesting that it could be used as an indicator of N enrichment in Mediterranean maquis. Changes in plant richness and evenness were correlated with the mean and/or the variation (standard deviation) of soil inorganic N parameters (e.g. nitrate concentration in the soil solution and the soil??s ratio of bioavailable N and phosphorus) measured along the time between the two plant community assessments. However, short- and long-term effects can be quite distinct, thus highlighting the need for further studies.     

On the Importance of Shrub Encroachment by Sprouters, Climate, Species Richness and Anthropic Factors for Ecosystem Multifunctionality in Semi-arid Mediterranean Ecosystems

One of the most important changes taking place in drylands worldwide is the increase of the cover and dominance of shrubs in areas formerly devoid of them (shrub encroachment). A large body of research has evaluated the causes and consequences of shrub encroachment for both ecosystem structure and functioning. However, there are virtually no studies evaluating how shrub encroachment affects the ability of ecosystems to maintain multiple functions and services simultaneously (multifunctionality). We aimed to do so by gathering data from ten ecosystem functions linked to the maintenance of primary production and nutrient cycling and storage (organic C, activity of β-glucosidase, pentoses, hexoses, total N, total available N, amino acids, proteins, available inorganic P, and phosphatase activity), and summarizing them in a multifunctionality index (M). We assessed how climate, species richness, anthropic factors (distance to the nearest town, sandy and asphalted road, and human population in the nearest town at several historical periods) and encroachment by sprouting shrubs impacted both the functions in isolation and M along a regional (ca. 350 km) gradient in Mediterranean grasslands and shrublands dominated by a non-sprouting shrub. Values of M were higher in those grasslands and shrublands containing sprouting shrubs (43 and 62%, respectively). A similar response was found when analyzing the different functions in isolation, as encroachment by sprouting shrubs increased functions by 2–80% compared to unencroached areas. Encroachment was the main driver of changes in M along the regional gradient evaluated, followed by anthropic factors and species richness. Climate had little effects on M in comparison to the other factors studied. Similar responses were observed when evaluating the functions in isolation. Overall, our results showed that M was higher at sites with higher sprouting shrub cover, longer distance to roads and higher perennial plant species richness. Our study is the first documenting that ecosystem multifunctionality in shrublands is enhanced by encroaching shrubs differing in size and leaf attributes. Our findings reinforce the idea that encroachment effects on ecosystem functioning cannot be generalized, and that are largely dependent on the traits of the encroaching shrub relative to those of the species being replaced.