1. Measuring ecosystem multifunctionality (EMF) could be a valuable method for assessing whether nature-based solutions, such as introducing non-native species to benefit ecosystem functioning, yield net positive outcomes. While EMF metrics may provide a coarse overview of the functioning of a system, such metrics can obscure the interactions among multiple simultaneously occurring ecosystem functions. This could limit the predictive value and overall efficacy of EMF metrics for the management of environmental change impacts on ecosystem functioning.
2. We carried out a mesocosm experiment with dung beetle temporal species richness (diversity) treatments in a pasture ecosystem in Aotearoa New Zealand and quantified cascading biotic and abiotic effects of introduced dung beetles (Onthophagus binodis and Onitis alexis) on multiple ecosystem functions. We measured nine individual ecosystem functions affected by dung beetles, including dung removal, total soil carbon, nitrogen, phosphorus, soil microbial respiration, nitrate N, ammonium N, anaerobically mineralisable N and pasture biomass production.
3. We first tested the effects of temporal dung beetle diversity on the nine individual ecosystem functions. We then tested how introduced dung beetle temporal diversity affects ecosystem multifunctionality by calculating EMF metrics using the averaging, threshold and the Hill numbers approaches. Finally, to investigate the underlying interactions among individual functions, we tested the direct and indirect effects of introduced dung beetles on the nine ecosystem functions using structural equation modelling (SEM).
4. Only one out of the nine individual ecosystem functions was significantly affected by dung beetle introductions, whereas five of the six calculated EMF metrics increased significantly with temporal diversity of introduced dung beetles. However, when modelled within an SEM framework, dung beetle temporal.