Herbivore pressure modulates soil multifunctionality in a Mediterranean landscape

Ongoing global environmental changes demand cost-effective strategies to restore ecosystem functioning. Within this context, the reintroduction of large mammalian herbivores has emerged as a promising approach to recover ecological processes and foster self-sustaining, biodiverse ecosystems. However, the impacts of such initiatives remain poorly understood, particularly regarding their effects on soil functions and ecological processes. We tested whether large herbivore densities (wild horses (Garrano breed) and cattle (Maronesa breed, an ancient lineage close to aurochs) modulates soil multifunctionality (the simultaneous provision of multiple ecosystem functions) across habitats (grasslands, shrublands, and forests) and seasons (fall and spring) in a Mediterranean landscape. Using composite multivariate indices describing soil functions related to biogeochemical cycles (i.e., enzymatic activities) and physicochemical properties (pH, conductivity, organic matter, and nutrient loads), we found that soil enzymatic activities varied seasonally, being higher in spring, and interacted with habitat and herbivore pressure gradients. Habitat and herbivore pressure explained the variation in soil multifunctionality during spring, while in fall it was mainly driven by herbivore pressure. Enzymatic stoichiometry, particularly C:P ratios, strongly predicted soil fertility and multifunctionality, showing positive relationships in forest habitats and under herbivore pressure. A mechanistic approach confirmed that herbivore impacts on soil functioning operated primarily through changes in soil properties and nutrient cycling rather than direct effects. Mediterranean landscapes are rapidly changing, and our results highlight herbivore management as a key tool to sustain soil ecological processes under the rewilding framework.