Soil Biodiversity-Ecosystem Functioning Relationships: Mechanisms, Methods, and Applications in a Changing World

Abstract

Understanding the relationships between soil biodiversity and ecosystem functioning (BEF) is critical for predicting ecosystem responses to global change and for informing sustainable land management. This article synthesizes current research on soil BEF relationships, exploring the foundational roles of diverse soil organisms—from microbes to macrofauna—in driving key processes like nutrient cycling, carbon sequestration, and plant productivity. We examine advanced methodological frameworks for quantifying these relationships, address challenges in troubleshooting context-dependent outcomes, and validate findings through cross-ecosystem comparisons. By integrating theoretical and applied perspectives, this review provides a comprehensive resource for researchers aiming to decipher the complex interplay between soil life and ecosystem multifunctionality, with implications for conservation, restoration, and climate change mitigation.

The Living Soil: Unraveling the Foundational Links Between Biodiversity and Ecosystem Processes

Soil biodiversity encompasses the variety of all living organisms within the soil, representing one of the most critical reservoirs of biological diversity on Earth [1]. This complex web of life includes a vast spectrum of organisms ranging from microscopic bacteria and archaea to more visible macrofauna such as earthworms and insects [2]. Soil is considered a major compartment of terrestrial ecosystems, hosting approximately a quarter of the planet's total taxonomic diversity, with species numbers several times higher than those observed aboveground [1]. This biodiversity is not merely a collection of species but constitutes a functional engine that drives essential ecosystem processes. The biological functioning of soils is intrinsically linked to microbial activity, which gives these microorganisms a major role in numerous soil functions [1]. Understanding this diversity—both structural and functional—is fundamental to testing biodiversity-ecosystem functioning (BEF) relationships in soils, a research area that has gained significant traction in recent years [3] [4]. The central premise of BEF research in soils investigates how the variety and abundance of soil organisms influence key ecosystem processes such as nutrient cycling, organic matter decomposition, and climate regulation [5].

The Structural Spectrum of Soil Biodiversity

Soil organisms are typically classified into functional groups based on body size, which influences their mobility through soil pores and their ecological roles [1] [3]. This classification system provides a practical framework for understanding the complex soil food web and its constituent components.

Table 1: Functional Classification of Soil Organisms

Organism Group	Size Range	Key Taxa Examples	Primary Functional Roles
Microorganisms (Microflora & Microfauna)	<0.2 mm	Bacteria, Archaea, Fungi, Protists, Nematodes	Organic matter decomposition, nutrient mineralization, pathogen suppression, biogeochemical cycling [1] [3] [6].
Mesofauna	0.2 - 2 mm	Springtails (Collembola), Mites (Acari), Enchytraeids	Litter fragmentation, microbial grazing, nutrient dispersal, soil microaggregate formation [2] [3].
Macrofauna	2 - 20 mm	Earthworms, Millipedes, Woodlice, Insect Larvae	Bioturbation, macro-pore creation, litter comminution, ecosystem engineering [3] [7].
Megafauna	>20 mm	Larger mammals (e.g., rabbits, bilbies, prairie dogs)	Soil aeration through burrowing, nutrient mixing, creating new niches for other organisms [6].

Microbial Diversity: The Biochemical Powerhouses

Microorganisms represent the most abundant and functionally diverse component of soil life. A single gram of soil can contain approximately one billion bacterial individuals from thousands of different species, alongside complex communities of fungi, archaea, and protists [1] [2]. Metagenomic studies have revealed that bacterial communities in diverse soils are often dominated by phyla such as Proteobacteria (∼40%), Acidobacteria (∼20%), and Actinobacteria (∼13%) [1]. Fungi are commonly classified by their functional modes: saprophytic fungi decompose dead organic matter; mycorrhizal fungi form symbiotic associations with plant roots (colonizing 80-90% of plant species); and pathogenic fungi can cause diseases in plants and animals [1]. These microbial communities are fundamental to regulating major biogeochemical cycles, including carbon, nitrogen, phosphorus, and sulfur, through processes like nitrogen fixation, nitrification, and organic matter mineralization [1] [5].

Soil Fauna: The Bioturbators and Engineers

Soil fauna, spanning from microfauna to macrofauna, play critical roles in modifying soil structure and regulating microbial communities. Earthworms, termites, and ants are recognized as key ecosystem engineers that significantly alter soil physical and chemical properties [6]. Their activities, including burrowing, nesting, and foraging, enhance soil aeration, water infiltration, and nutrient mixing. The functional diversity of soil macrofauna has been shown to help stabilize microbial diversity and community composition during periods of environmental stress, such as severe and prolonged drought [7]. Global estimates indicate that ants and termites alone can turn over between 0.001 to 10 tonnes of soil per hectare per year, substantially impacting soil formation processes, particularly in drylands [6].

Functional Groups and Ecosystem Processes

Beyond taxonomic classification, soil biodiversity can be understood through a functional lens—describing "who does what" in the ecosystem rather than merely "who is present" [1]. This functional perspective is crucial for understanding BEF relationships, as it focuses on the processes and services derived from soil organisms.

Table 2: Key Ecosystem Functions Driven by Soil Biodiversity

Ecosystem Function	Description	Primary Organisms Involved	Experimental Evidence
Organic Matter Decomposition	Breakdown of complex organic compounds into simpler molecules	Saprophytic fungi, bacteria, earthworms, springtails, millipedes	Litterbag studies show reduced decomposition rates in fauna-excluded treatments [1] [7].
Nutrient Cycling	Transformation and mobilization of essential nutrients (C, N, P)	Nitrifying & nitrogen-fixing bacteria, mycorrhizal fungi, nematodes, protozoa	Metagenomics and isotope tracing link microbial diversity to nutrient flux rates [1] [5].
Soil Structure Formation	Creation and stabilization of soil aggregates	Arbuscular mycorrhizal fungi, earthworms, termites, enchytraeids	Microcosm experiments demonstrate improved aggregation with diverse inocula [5] [6].
Climate Regulation	Sequestration of carbon in soil organic matter	Fungi, bacteria, root symbionts, bioturbators	Global surveys correlate microbial composition with soil C stocks [5] [6].
Pathogen Suppression	Control of plant and animal pathogens through competition and predation	Antagonistic bacteria, fungivorous nematodes, predatory mites	Bioassays show reduced disease incidence in soils with higher biodiversity [3] [6].

Multifunctionality and Functional Redundancy

A key concept in BEF research is ecosystem multifunctionality—the simultaneous performance of multiple ecosystem processes [5]. Recent research demonstrates that soil biodiversity is often positively associated with soil multifunctionality across natural, urban, and agricultural ecosystems [6]. While soil biodiversity was traditionally considered highly functionally redundant (where the loss of certain taxa could be compensated by others with similar functions), contemporary research suggests that losses in microbial diversity can result in proportional or exponential declines in specialized soil functions [6]. The relationship between biodiversity and function is context-dependent, with soil biodiversity playing a particularly critical role in supporting function in nutrient-poor soils and drylands, where the daily contribution of diverse soil microbiota is essential for decomposing organic matter and releasing nutrients [6].

Experimental Approaches for Testing BEF Relationships

Research on soil biodiversity and ecosystem functioning employs a range of experimental protocols, from highly controlled microcosm studies to large-scale field observations. These approaches aim to establish causal relationships between organismal diversity and process rates.

Biodiversity Manipulation Experiments

A common experimental design involves creating gradients of soil biodiversity and measuring subsequent effects on ecosystem functions.

Table 3: Key Experimental Methodologies in Soil BEF Research

Method Type	Protocol Description	Key Measurements	Applications & Limitations
Size-Fraction Filtering	Soil inocula are sequentially filtered through meshes of decreasing size (e.g., from 5000 μm to 6 μm) to selectively exclude organism groups based on body size [3].	Plant biomass, soil respiration, nutrient leaching, microbial community composition via DNA sequencing.	Creates a biodiversity gradient; may cause collateral damage to microbial communities.
Trophic Simplification	Removal or addition of specific functional groups (e.g., predators, decomposers) to assess their role in food webs.	Decomposition rates, nutrient mineralization, prey population dynamics, gas fluxes.	Isulates trophic interactions; difficult to maintain in open systems.
Dilution-to-Extinction	Serial dilution of soil suspensions to reduce microbial diversity while maintaining similar abiotic conditions.	Enzyme activities, substrate utilization profiles, process rates (e.g., nitrification).	Tests diversity-effects independent of density; may not reflect natural extinction sequences.
Mesocosm Systems	Controlled-environment experiments (e.g., Ecotron facilities) with manipulated biodiversity and environmental factors [7].	Multifunctionality indices, community resilience, functional gene expression, gas fluxes.	High realism and control; expensive and limited spatial scale.

A representative experiment by Radujković et al. [3] investigated the importance of soil biodiversity across five European grasslands with differing soil properties. The experimental protocol involved:

• Soil Collection: Soils were collected from five grassland sites across Europe with varying fertility levels.
• Sterilization and Re-inoculation: Soils were sterilized by gamma irradiation to eliminate native biota while preserving soil physicochemical properties.
• Biodiversity Gradient Establishment: Soil communities were filtered by size to create four biodiversity treatments: (1) Max - containing microorganisms, microfauna, and mesofauna; (2) Medium - containing microorganisms and microfauna; (3) Low - containing only microorganisms; and (4) Min - a minimal microbial community.
• Ecosystem Function Monitoring: Over the experimental period, researchers measured plant productivity (shoot and root biomass), soil microbial respiration, nutrient cycling indicators (N and P availability), and microbial community composition via molecular methods.

The results demonstrated that soil biodiversity decline often induced gradual, stepwise changes in ecosystem functions, suggesting that decreasing richness of different soil organism groups progressively influences how ecosystems operate [3].

Figure 1: Experimental Workflow for Soil BEF Studies. This diagram illustrates the common protocol for establishing biodiversity gradients in soil experiments.

Measuring Biodiversity and Function Responses

Advanced molecular techniques have revolutionized our ability to characterize soil biodiversity. DNA metabarcoding allows for the comprehensive identification of species present in soil samples, providing detailed information on community composition across all domains of life [2]. For functional assessments, researchers employ:

• Physiological Profiling: Microbial respiration, substrate-induced respiration, and enzyme activity assays.
• Stable Isotope Probing: Using 13C- or 15N-labeled substrates to trace nutrient flow through food webs.
• Metatranscriptomics and Metaproteomics: Assessing gene expression and protein production in microbial communities.
• Gas Flux Measurements: Quantifying greenhouse gas emissions (CO2, N2O, CH4) to assess climate-relevant processes.

The Scientist's Toolkit: Essential Research Reagents and Platforms

Table 4: Key Research Reagents and Platforms for Soil BEF Studies

Tool Category	Specific Examples	Function & Application
Molecular Analysis Kits	DNA/RNA extraction kits (e.g., MoBio PowerSoil), PCR reagents, sequencing library prep kits	Extraction and preparation of nucleic acids from complex soil matrices for community analysis [2].
Stable Isotopes	13C-labeled plant litter, 15N-ammonium nitrate, 18O-water	Tracing nutrient pathways through food webs; quantifying process rates in situ.
Environmental Data Platforms	Edaphobase data warehouse [8]	Collating and using soil biodiversity datasets with environmental metadata for large-scale analyses.
Controlled Environment Facilities	Ecotron mesocosms [7]	Precisely controlling environmental conditions (e.g., drought) while manipulating biodiversity.
Bioinformatics Tools	QIIME2, MOTHUR, PICRUSt2, FUNGuild	Processing sequencing data, assigning taxonomy, predicting functional potential.

Soil biodiversity, spanning from microbes to macrofauna, represents a complex hierarchical system that drives essential ecosystem functions through specialized functional groups. Experimental evidence consistently demonstrates that reductions in soil biodiversity can impair multiple ecosystem processes, though these relationships are influenced by environmental context and the degree of functional redundancy [3] [6]. The emerging consensus from BEF research is that different components of soil biodiversity contribute uniquely to ecosystem multifunctionality, with larger organisms often playing disproportionate roles as ecosystem engineers while microbial communities drive biochemical transformations [5] [6].

Future research directions include better understanding how environmental autocorrelation influences BEF relationships across spatial and temporal scales [4], developing management strategies that enhance microbial functions rather than simply abundance [9], and addressing trade-offs between biodiversity conservation and agricultural productivity [9]. As soil biodiversity faces increasing threats from global environmental change, elucidating the precise mechanisms linking soil organism diversity to ecosystem functioning becomes increasingly critical for conservation and sustainable management.

Soil biodiversity harbors a substantial fraction of the world's biodiversity and is a fundamental driver of many crucial ecosystem functions [10]. The organisms within soil—including microorganisms like bacteria and fungi, and fauna such as nematodes, earthworms, and collembolans—form complex biological networks whose activities are essential for terrestrial ecosystem processes [10] [5]. The functional traits of these organisms and the multitrophic interactions between them support key ecosystem processes including the decomposition of organic matter, nutrient cycling, soil carbon sequestration, and the maintenance of plant health [5]. Historically, the importance of soil biota was recognized by early scientists like Darwin, who identified soil fauna as an "engine" of ecosystem functioning [11]. However, this understanding was largely overshadowed during the 20th century by agricultural practices emphasizing external inputs like fertilizers and pesticides [11]. A resurgence of interest in soil ecology began in the 1990s, recognizing that the biodiversity within soils plays a disproportionately large role in governing ecosystem multifunctionality—the simultaneous provision of multiple ecosystem functions [5] [11]. This article examines the critical keystone functions performed by soil communities, exploring the mechanistic relationships between biodiversity and ecosystem processes that are essential for environmental sustainability and human wellbeing.

Experimental Evidence: Linking Keystone Taxa to Ecosystem Processes

Keystone Taxa Regulation of Plant Residue Decomposition and Carbon Chemistry

A long-term, large-scale field investigation examined how microbial keystone taxa drive the decomposition of plant residues and the subsequent transformation of carbon chemistry over a 9-year period [12]. The study tracked the decomposition processes of maize straw and wheat straw across three distinct climate zones, analyzing residue composition dynamics and associated microbiomes to understand the temporal patterns of this critical process.

Table 1: Experimental Protocol for Keystone Taxa in Decomposition Study

Experimental Element	Specification	Purpose
Experimental Sites	Three agroecological stations across climate gradients: Hailun (cold temperate), Fengqiu (warm temperate), Yingtan (mid-subtropical)	To assess impact of varying environmental conditions on decomposition processes
Experimental Design	Soil transplantation with three typical soils: Phaeozem, Cambisol, Acrisol; Litter bags containing maize and wheat straw buried at 12cm depth	To separate effects of soil type, climate, and residue quality on decomposition
Sampling Timeline	7 sampling points over 9 years: 1 month, 3 months, 0.5, 1, 2, 3, and 9 years after burial	To capture temporal dynamics of chemical and microbial succession
Analytical Methods	Solid-state 13C CP-MAS NMR; Py-GCMS; Extracellular enzyme assays; Amino sugar analysis; Co-occurrence network analysis	To characterize chemical composition, microbial activity, and community interactions

The research revealed that residue chemistry followed a divergent-convergent trajectory during decomposition [12]. During the initial 0.5-3 year period, residue composition diverged under the combined influence of straw type and climate conditions, followed by a convergence toward a common set of compounds during the 3-9 year period. This trajectory was primarily governed by shifts in microbial metabolic strategies: the initial divergent phase was driven by microbial catabolism (extracellular enzymatic degradation), while the later convergent phase was dominated by microbial anabolism (assimilation into microbial biomass) [12]. Keystone taxa belonging to Alphaproteobacteria (particularly Rhizobiales) guided bacterial networks that regulated extracellular enzyme activity during the early decomposition phase, whereas fungi, particularly Chaetomium, became the main contributors to microbial assimilation in the later stages [12].

Figure 1: Keystone Taxa Regulation of Plant Residue Decomposition Trajectory. The diagram illustrates the temporal succession of microbial processes governing the divergent-convergent pathway of residue chemistry transformation during long-term decomposition.

Mangrove Restoration Effects on Soil Multifunctionality Through Microbial Communities

A separate investigation examined how microbial diversity and keystone species drive soil nutrient cycling and multifunctionality following mangrove restoration [13]. This study evaluated a chronosequence of restored mangrove sites of different ages (2, 10, and 20 years since restoration) alongside unrestored mudflats to understand ecosystem recovery patterns.

Table 2: Experimental Protocol for Mangrove Restoration Study

Experimental Element	Specification	Purpose
Study Sites	Quanzhou Bay mangrove Reserve, Fujian Province, China; Sites with restoration histories of 2 years (M2), 10 years (M10), 20 years (M20), plus unrestored mudflats (M0)	To create space-for-time substitution chronosequence to study restoration trajectory
Field Measurements	Soil sampling and analysis of chemical properties; Assessment of extracellular enzyme activities; Plant community surveys	To quantify changes in soil properties and ecosystem functions through restoration time
Microbial Analysis	High-throughput sequencing of 16S rRNA and ITS genes; Co-occurrence network analysis; Identification of keystone taxa	To characterize microbial community composition, structure, and identify key species
Statistical Analysis	Calculation of soil multifunctionality index; Regression analysis; Structural equation modeling (SEM)	To establish links between restoration age, microbial communities, and ecosystem functioning

The findings demonstrated that mangrove restoration significantly enhanced soil multifunctionality, with older restoration sites exhibiting greater functional enhancement [13]. Specifically, soil total carbon (TC), total nitrogen (TN), total phosphorus (TP), and soil organic matter (SOM) concentrations showed substantial increases at M10 and M20 sites compared to unrestored mudflats. The research identified specific keystone microbial taxa (including Methylomirabilota, Anaerolineaceae, and Acidobacteriota) that played disproportionately important roles in regulating nutrient cycling processes [13]. Structural equation modeling revealed that mangrove restoration enhanced soil multifunctionality primarily through direct pathways (β = 0.41) and indirectly by increasing microbial diversity (β = 0.30) and the abundance of keystone taxa (β = 0.44) [13]. This demonstrates that the recovery of ecosystem functions following restoration is mediated by complex interactions between plant communities, soil properties, and soil microbial networks.

Research Reagent Solutions for Soil Biodiversity-Function Studies

Table 3: Essential Research Tools for Investigating Soil Biodiversity-Ecosystem Function Relationships

Research Tool Category	Specific Examples	Research Application	Key Function in Experiments
Molecular Community Profiling	16S rRNA sequencing (bacteria); ITS sequencing (fungi); Metagenomics	Characterizing taxonomic and functional diversity of soil microbial communities	Identification of keystone taxa and functional genes; Assessment of community structure changes
Chemical Composition Analysis	Solid-state 13C CP-MAS NMR; Pyrolysis-GCMS (Py-GCMS)	Tracking chemical transformations of organic matter during decomposition	Quantification of carbon chemistry changes; Assessment of decomposition pathways
Enzyme Activity Assays	β-1,4-glucosidase (βG); Cellobiohydrolase (CBH); Phenol oxidases (PhOX); Peroxidases (Perox)	Measuring microbial functional potential in nutrient cycling	Assessment of decomposition rates; Evaluation of nutrient mineralization capacities
Stable Isotope Techniques	13C and 15N labeling; Amino sugar analysis	Tracing nutrient flows through soil food webs	Quantification of carbon sequestration pathways; Assessment of microbial nutrient assimilation
Network Analysis Tools	Co-occurrence network analysis; Structural equation modeling (SEM)	Identifying microbial interactions and ecosystem relationships	Detection of keystone taxa; Elucidation of direct and indirect pathways in ecosystem functioning

Discussion: Implications for Conservation and Management

The experimental evidence demonstrates that soil biodiversity, particularly specific keystone taxa, plays a fundamental role in driving essential ecosystem functions including decomposition, nutrient cycling, and carbon sequestration [5] [12] [13]. These relationships between biodiversity and functioning have critical implications for ecosystem management and conservation strategies. Despite their importance, significant gaps remain in our understanding of global soil biodiversity patterns and their relationship to ecosystem functioning [10]. Current data show strong biogeographical biases in soil biodiversity research, with most sampling sites concentrated in temperate regions and important spatial, environmental, taxonomic, and functional gaps in coverage [10]. Alarmingly, only 0.3% of all sampling sites in macroecological studies have simultaneous information on both biodiversity and ecosystem function, severely limiting our ability to establish general patterns in soil biodiversity-ecosystem functioning relationships across different ecosystems [10].

The emerging understanding of soil keystone functions underscores the importance of integrating soil biodiversity conservation into land management practices and environmental policy [5] [11]. Soil organisms contribute significantly to climate regulation through their roles in carbon cycling, with microbial anabolism during decomposition leading to the formation of microbial necromass that contributes to stable soil organic matter pools [5] [12]. The manipulation of soil communities, such as through the application of synthetic microbial communities (SynComs), shows promise for enhancing ecosystem functions in agricultural and restoration contexts [14]. Future research priorities should focus on filling critical data gaps, particularly in underrepresented ecosystems and taxonomic groups, and on developing a more mechanistic understanding of how global environmental changes affect the relationships between soil biodiversity and ecosystem functioning [10] [5]. By placing soil biodiversity at the center of ecosystem management, we can better harness the fundamental keystone functions that support terrestrial ecosystem sustainability and human wellbeing.

The paradigm of soil biodiversity has undergone a fundamental shift in ecological science. Historically, soil communities were considered functionally redundant, with the loss of specific taxa presumed to have minimal ecosystem consequences due to compensatory mechanisms from remaining organisms [15] [6]. Contemporary research has fundamentally overturned this perspective, revealing that soil biodiversity supports ecosystem multifunctionality—the simultaneous performance of multiple ecosystem processes—through mechanisms more complex and essential than previously recognized [16] [15] [17]. This synthesis examines the experimental evidence establishing diverse soil communities as crucial drivers of simultaneous ecosystem services, addressing a core thesis in biodiversity-ecosystem functioning (BEF) relationships within soil ecosystems.

The functional importance of soil biodiversity extends beyond simple species-function relationships to encompass multitrophic interactions, community assembly processes, and the distinct roles of rare versus abundant taxa [16] [15]. Understanding these mechanisms provides critical insights for managing soil communities to support agricultural productivity, climate change mitigation, and ecosystem restoration amid global environmental change.

Experimental Comparisons: Diverse vs. Depauperate Soil Communities

Quantitative Evidence from Observational and Manipulative Studies

Table 1: Multifunctionality Responses to Soil Biodiversity Gradients Across Ecosystems

Ecosystem Context	Biodiversity Metric	Key Functions Measured	Multifunctionality Response	Citation
Agricultural fields (228 sites, Eastern China)	Rare vs. abundant soil taxa diversity	16 functions related to nutrient provisioning, element cycling, pathogen control, plant-microbe symbiosis	Positive relationship with rare taxa diversity (explained most variation); No significant relationship with abundant taxa diversity	[16]
Atlantic Forest landscapes (natural vs. degraded)	Multiple soil organism groups	Primary production, ecosystem stability, nutrient cycling	Significant positive correlation (p<0.01); Degraded landscapes showed reduced multifunctionality	[18]
Temperate grasslands (Switzerland)	Management intensity modulation	22 ecosystem service indicators across provisioning, regulating, cultural services	Extensive management (higher biodiversity) enhanced multifunctionality, particularly cultural services	[19]
Global Change Experiment (BioCON, Minnesota)	Plant species richness	Root biomass, soil respiration, microbial biomass, soil aggregation	Positive diversity effects on 3/4 functions and multifunctionality across ambient and future environments (eCO₂, eN)	[17]

The Critical Role of Rare Taxa

A pivotal advancement in BEF research has been the recognition that rare soil taxa disproportionately drive multifunctionality relationships. In extensive agricultural surveys, the diversity of rare species (relative abundance <0.05%) demonstrated strong positive relationships with multifunctionality, while abundant species diversity showed no significant relationship [16]. This paradigm-shifting finding challenges conventional emphasis on dominant taxa and reveals that the "rare biosphere" provides essential, specialized functional support.

Experimental evidence indicates distinct functional roles between rare and abundant soil taxa:

• Abundant taxa maintain broader functional performance, with individual phylotypes supporting multiple functions simultaneously, particularly in nutrient provisioning and element cycling [16]
• Rare taxa contribute more unique phylotypes supporting individual ecosystem functions, with particular importance for specialized processes like pathogen control and plant-microbe symbiosis [16]
• Community assembly processes significantly influence functional performance, with stochastic processes strengthening positive diversity-multifunctionality relationships in rare subcommunities [16]

Methodological Framework: Experimental Protocols for Soil BEF Research

Protocol 1: Assessing Multifunctionality Across Biodiversity Gradients

Field Sampling Design:

• Site selection: Implement stratified sampling across management intensity gradients (e.g., conventional to organic systems) or natural biodiversity gradients [16] [19]
• Soil sampling: Collect composite soil samples (0-15cm depth) during consistent seasonal periods to control for temporal variation
• Biodiversity quantification: Extract DNA for high-throughput sequencing of multiple organismal groups (bacteria, archaea, fungi, protists) using standardized primer sets [16] [18]

Multifunctionality Assessment:

• Function measurement: Quantify 16+ ecosystem functions spanning nutrient cycling (N, P, C transformations), pathogen suppression, soil structure maintenance, and productivity metrics [16]
• Standardization: Normalize individual function measurements to scale 0-1 relative to maximum observed values
• Integration: Calculate multifunctionality indices using averaging approaches (mean of standardized functions) and multiple threshold methods (number of functions exceeding critical thresholds) [17]

Protocol 2: Disentangling Rare vs. Abundant Taxon Contributions

Community Partitioning:

• Sequence processing: Cluster sequences into operational taxonomic units (OTUs) at 97% similarity threshold
• Abundance classification: Define "abundant" taxa as >0.5% relative abundance; "rare" taxa as <0.05% relative abundance; "intermediate" as 0.05-0.5% [16]
• Subcommunity analysis: Analyze diversity-function relationships separately for each abundance cohort

Statistical Modeling:

• Structural Equation Modeling (SEM): Quantify direct and indirect paths linking biodiversity components to multifunctionality while accounting for environmental covariates (climate, soil properties) [16]
• Random Forest analysis: Validate the predictive importance of rare versus abundant taxa diversity for multifunctionality [16]
• Network analysis: Construct correlation networks to identify keystone taxa and their associations with multiple ecosystem functions [16]

Conceptual Synthesis: Mechanisms Underlying Biodiversity-Multifunctionality Relationships

Soil Biodiversity-Multifunctionality Framework

Context Dependency in BEF Relationships

The strength of biodiversity-multifunctionality relationships demonstrates significant context dependency, moderated by several environmental and management factors:

Land Use Intensity: Agricultural intensification consistently weakens positive BEF relationships, with simplified landscapes (deforested areas, intensive pastures) exhibiting reduced multifunctionality compared to natural ecosystems [18]. Extensive management practices enhance multifunctionality, particularly for cultural and regulating services, though often at the expense of provisioning services [19].

Environmental Change Scenarios: Plant diversity effects on soil multifunctionality persist under future climate conditions, including elevated CO₂ and nitrogen deposition, though nitrogen enrichment can weaken these relationships at higher function thresholds [17]. This resilience demonstrates the robustness of biodiversity as a conservation strategy for maintaining ecosystem services under changing environmental conditions.

Soil Properties and Climate: Soil biodiversity contributes more significantly to ecosystem functioning in resource-poor environments, including drylands and poorly developed soils, where microbial communities are essential for nutrient depolymerization and daily ecosystem functioning [15] [6].

The Scientist's Toolkit: Essential Research Reagents and Methodologies

Table 2: Essential Research Solutions for Soil Biodiversity-Function Studies

Reagent/Technology	Primary Application	Key Function in BEF Research	Experimental Considerations
High-throughput DNA sequencing kits (16S, 18S, ITS, COI markers)	Characterization of microbial and faunal diversity across taxa	Provides comprehensive biodiversity inventories for correlation with function measures	Multi-marker approach essential for capturing full community diversity; requires standardization across studies
Standardized DNA extraction kits (e.g., MoBio PowerSoil)	Nucleic acid extraction from diverse soil types	Ensures comparable molecular analysis across environmental gradients	Critical for removing PCR inhibitors; extraction efficiency varies with soil properties
Microbial metabolic profiling plates (e.g., EcoPlates, Biolog)	Functional potential assessment of soil communities	Measures catabolic diversity and substrate utilization patterns	Provides functional complement to taxonomic data; laboratory conditions may not reflect field reality
Stable isotope tracers (¹³C, ¹⁵N)	Nutrient flow and process rate quantification	Tracks element cycling pathways through soil food webs	Enables direct measurement of process rates rather than potential functions
SynComs (Synthetic Communities)	Manipulative experiments testing BEF mechanisms	Isolates effects of specific taxa or functional groups on ecosystem processes	Allows controlled testing of biodiversity hypotheses; complexity reduction may limit real-world applicability
Soil enzyme assay kits	Biochemical process measurement	Quantifies potential decomposition and nutrient cycling rates	Standardized colorimetric methods enable cross-study comparisons; represents potential rather than in situ rates

The accumulating evidence unequivocally demonstrates that diverse soil communities are fundamental to ecosystem multifunctionality across terrestrial ecosystems. The functional importance of soil biodiversity extends beyond simple redundancy to encompass specialized contributions from rare taxa, multitrophic interactions, and context-dependent responses to global change drivers [16] [15] [17]. These insights provide a mechanistic basis for managing soil communities as essential components of sustainable agricultural systems, climate-resilient ecosystems, and restoration initiatives.

Future research directions should prioritize (1) understanding the genomic and metabolic basis of multifunctionality, (2) elucidating the dynamics of soil biodiversity in response to simultaneous global change drivers, and (3) developing management practices that specifically enhance the diversity-function relationships identified in experimental systems. Integrating these insights into land management and policy frameworks offers a pathway to enhance multiple ecosystem services simultaneously through targeted conservation of soil biodiversity.

Understanding the mechanisms through which soil biodiversity influences ecosystem processes is a central goal in soil ecology. Three foundational theoretical frameworks have emerged to explain this relationship: the complementarity effect, selection effect, and functional redundancy. These frameworks offer distinct yet interconnected explanations for how the immense diversity of soil organisms—from bacteria and fungi to nematodes and earthworms—collectively drives essential functions like nutrient cycling, carbon sequestration, and organic matter decomposition [15] [6]. Research demonstrates that soil biodiversity supports multiple ecosystem functions simultaneously, a concept known as ecosystem multifunctionality [20] [18]. As global change stressors threaten soil ecosystems, deciphering the relative contributions of these mechanisms becomes crucial for predicting ecosystem responses and informing conservation strategies [21]. This guide provides a comparative analysis of these frameworks, synthesizing experimental evidence and methodologies to equip researchers with tools for testing biodiversity-ecosystem functioning relationships in soil systems.

Theoretical Framework Comparison

The table below defines, describes the underlying mechanisms, and outlines the research support for each of the three core theoretical frameworks.

Framework	Definition & Core Principle	Underlying Mechanism	Key Research Support
Complementarity Effect	Positive biodiversity effect arising from niche differentiation and facilitative interactions among species [21].	- Resource partitioning: Different species utilize distinct resources or the same resource in different ways [22].- Abiotic facilitation: One species modifies the environment to benefit another [21].- Biotic feedbacks: Interactions between species enhance overall community performance.	Long-term grassland experiments show complementarity effects strengthen over time and drive positive biodiversity-productivity relationships [22] [21].
Selection Effect	Biodiversity effect occurring when particular high-performing species become dominant in diverse communities [21].	- Probability effect: Diverse communities have a higher chance of containing a highly productive species.- Dominance effect: This high-performing species outcompetes others and disproportionately contributes to ecosystem function.	Observed in plant communities where a dominant plant species drives community biomass production; its influence can be altered by global change factors like N and CO₂ enrichment [21].
Functional Redundancy	The concept that multiple species perform similar ecological roles, and the loss of one can be compensated for by others [22].	- Functional similarity: Different species possess overlapping functional traits and can maintain specific processes if one is lost [22].- Ecosystem resilience: Redundancy provides insurance against species loss, stabilizing ecosystem functions.	Traditionally applied to soil microbes for broad functions like decomposition; however, the term is now debated and may be context-dependent [22].

Experimental Evidence and Data

Empirical studies across different ecosystems have quantified the roles of complementarity, selection, and redundancy. The following table summarizes key experimental findings.

Experimental Context	Key Findings on Theoretical Frameworks	Implications
Global Change Experiment (BioCON)24-year grassland study [21]	- Complementarity was the primary driver of the positive biodiversity-ecosystem functioning relationship.- Selection effects were not significantly altered by biodiversity alone.- N addition and CO₂ enrichment interacted to diminish the positive effects of biodiversity on complementarity.	Long-term global change can disrupt the key mechanisms through which biodiversity supports ecosystem productivity.
Soil Community ManipulationModel grassland microcosms [20]	- Soil biodiversity loss impaired multiple ecosystem functions (e.g., plant diversity, nutrient retention, decomposition).- Ecosystem multifunctionality exhibited a strong positive linear relationship with a composite soil biodiversity index.	Supports the complementarity effect, showing that diverse soil communities are needed to maintain multiple functions simultaneously.
Microbial Diversity Loss SynthesisLiterature review [23]	- Loss of microbial diversity decreased both specific and general functions.- Challenged the assumption of widespread functional redundancy, even for general processes like C and N mineralization.	Functional redundancy in soil microbial communities may be overestimated; diversity is critical for a wide range of functions.
Atlantic Forest Landscape StudyField assessment [18]	- Soil biodiversity was positively correlated with primary production, ecosystem stability, and nutrient cycling.- Land-use simplification reduced biodiversity and these related functions.	Provides field evidence for the role of soil biodiversity, via complementarity, in supporting multifunctionality in natural ecosystems.

Detailed Experimental Protocols

Soil Community Fractionation and Microcosm Setup

This protocol, adapted from a foundational study [20], tests the effects of soil biodiversity loss on ecosystem multifunctionality by creating a gradient of soil community complexity.

• Soil Inoculum Collection: Collect bulk soil from a representative grassland site. Remove large rocks and visible plant debris.
• Community Fractionation: Process the soil inoculum through a series of sieves with decreasing mesh sizes (e.g., 5000 μm, 250 μm, 50 μm, 20 μm) to create different community treatments.
  • Rationale: Larger mesh sizes retain a more complete soil community (including macrofauna, mesofauna, and microbes). Smaller mesh sizes sequentially remove specific groups (e.g., nematodes, mycorrhizal fungi, and larger microorganisms), simplifying the community composition and reducing overall diversity [20].
• Microcosm Establishment: Fill sterile microcosms (e.g., self-contained pots) with a sterilized background soil matrix. Inoculate each microcosm with one of the fractionated soil communities.
• Plant Community Setup: Sow each microcosm with a standardized model grassland plant community.
• Maintenance and Monitoring: Maintain the microcosms under controlled environmental conditions. Prevent cross-contamination between treatments.
• Ecosystem Function Measurement: After a designated growth period, measure a suite of ecosystem functions, including:
  • Plant Productivity: Above-ground and below-ground biomass.
  • Nutrient Cycling: Litter decomposition rates, nitrogen leaching, and nitrous oxide (N₂O) emissions.
  • Soil Carbon: Carbon sequestration in soil.
  • Plant Diversity: Number and abundance of plant species.

Field-Based Assessment of Biodiversity-Multifunctionality Relationships

This protocol outlines a correlative approach to link soil biodiversity with ecosystem multifunctionality across natural and managed landscapes [18].

• Site Selection: Select paired field sites representing different land-use types (e.g., natural forest, pasture, deforested area) to capture a gradient of disturbance and biodiversity.
• Soil Sampling: Conduct seasonal soil sampling (e.g., dry and rainy seasons) over multiple years to account for temporal variation. Collect multiple soil cores per site and combine them into a composite sample for analysis.
• Biodiversity Quantification:
  • Macrofauna: Extract and identify insects, arachnids, and myriapods from soil samples.
  • Microfauna: Extract nematodes using centrifugal-flotation techniques.
  • Microbes: Quantify bacterial and fungal diversity and abundance via molecular techniques (e.g., DNA sequencing and qPCR).
  • Mycorrhizal Fungi: Quantify spore density and root colonization rates.
• Ecosystem Function Measurement: Assess key ecosystem functions from the same sites:
  • Primary Production (PP): Estimate via above-ground plant biomass.
  • Nutrient Cycling (NC): Measure soil nitrogen and phosphorus availability.
  • Ecosystem Stability (ES): Calculate the temporal stability of biomass production.
• Data Analysis: Use Pearson's correlation analysis and linear mixed-effects models (LMMs) to relate soil biodiversity indices (e.g., overall diversity or diversity of specific groups) to ecosystem multifunctionality, while accounting for random effects like site and season [18].

Conceptual Diagrams of Frameworks and Mechanisms

Biodiversity-Ecosystem Functioning Theoretical Frameworks

Global Change Impacts on BEF Relationships

The Scientist's Toolkit: Key Research Reagents & Materials

The following table details essential reagents, materials, and tools used in experimental soil biodiversity and ecosystem functioning research.

Category / Item	Specific Examples	Function / Application in Research
Soil Fractionation	Sieve sets (e.g., 5000μm, 250μm, 50μm) [20]	Creates a biodiversity gradient by mechanically separating soil organisms based on body size.
DNA/RNA Analysis	DNA extraction kits, PCR reagents, primers (e.g., 16S/18S/ITS), sequencing platforms [18]	Quantifies microbial and fungal diversity, community composition, and functional gene abundance.
Microcosm Systems	Sterile pots, growth chambers, controlled-environment rooms [20]	Provides a self-contained, replicable experimental system to test cause-effect relationships.
Fauna Extraction	Berlese-Tullgren funnels, centrifugal-flotation systems [18]	Extracts active microarthropods and nematodes from soil samples for identification and counting.
Mycorrhizal Analysis	Staining solutions (e.g., trypan blue), microscopes [18]	Quantifies mycorrhizal root colonization and spore density, key to plant nutrient uptake.
Gas Flux Measurement	Gas chromatographs, infrared gas analyzers (IRGA) [20]	Measures ecosystem functions like soil respiration (CO₂ flux) and denitrification (N₂O flux).
Nutrient Analysis	Colorimetric assays, ion chromatographs, elemental analyzers [20] [18]	Quantifies nutrient cycling (e.g., N, P) and carbon sequestration in soil and plant tissues.

The theoretical frameworks of complementarity, selection effects, and functional redundancy provide powerful, complementary lenses through which to understand the soil biodiversity-ecosystem functioning relationship. Empirical evidence strongly indicates that complementarity effects are a primary mechanism underpinning ecosystem multifunctionality, while the concept of strict functional redundancy may be an overgeneralization [22]. Critically, long-term studies reveal that global change factors can interact to disrupt these fundamental mechanisms [21]. Future research must address critical gaps, including the functioning of bare soils, intra-specific trait variability, and multi-trophic interactions [24]. Integrating these frameworks and employing robust, standardized experimental protocols will be essential for advancing predictive soil ecology and developing effective strategies for soil conservation and restoration in an era of global change.

Soil biodiversity, encompassing organisms from bacteria and fungi to earthworms and nematodes, forms the foundation of virtually all terrestrial ecosystems. The relationship between this biodiversity and ecosystem functioning (BEF)—termed the "soil biodiversity-ecosystem functioning" relationship—has emerged as a critical research frontier in ecology [3]. Soils represent the most biodiverse systems on Earth, harboring complex communities of microorganisms (archaea, bacteria, fungi, protists), microfauna (e.g., nematodes), mesofauna (e.g., microarthropods), and macrofauna (e.g., earthworms, beetles) [3]. These organisms collectively drive essential ecosystem processes including nutrient cycling, organic matter decomposition, carbon sequestration, and plant productivity [6]. Understanding how global change factors—particularly land-use shifts and climate change—alter these relationships is paramount for predicting ecosystem responses and developing conservation strategies.

In the face of accelerated global change, including climate change, elevated carbon dioxide, nitrogen deposition, and land-use intensification, the stability of ecosystem functions has become a major concern [25]. The 21st century has witnessed growing recognition that soil biodiversity supports almost every ecosystem function we know, with profound implications for human well-being through its contributions to food security, climate change mitigation, and the One Health concept linking human, animal, and environmental health [6]. Simultaneously, a third of all soils worldwide are already significantly degraded, with less than 40% of European soils in healthy condition—a concerning statistic given that soil is essentially a non-renewable resource requiring up to thousands of years to form a single centimeter [6]. This review synthesizes current understanding of how land-use and climate changes are reshaping soil BEF relationships, with particular emphasis on experimental evidence, mechanistic insights, and conservation strategies.

Theoretical Foundation: Multidimensional Biodiversity and Ecosystem Functioning

The Multidimensional Nature of Biodiversity-Ecosystem Functioning Relationships

The relationship between biodiversity and ecosystem functioning extends beyond simple species counts to encompass multiple dimensions of biological organization. Contemporary BEF research recognizes that assessing biodiversity effects on ecosystem functioning requires consideration of multiple dimensions of diversity: diversity across multiple trophic levels (plants, animals, and microbes), multiple facets (taxonomic, functional traits, and phylogeny), and multiple spatial scales (local, regional, and landscape) [25]. This multidimensional perspective reveals that high levels of biodiversity can enhance ecosystem stability and mitigate the negative impacts of global change on ecosystem functioning [25].

Soil biodiversity was traditionally considered highly functionally redundant, meaning the loss of certain soil taxa was thought to have little impact on soil functions because their roles would be replaced by others [6]. However, emerging research challenges this view, demonstrating that losses in microbial diversity can result in proportional or exponential losses of soil function, particularly for specialized processes such as denitrification [6]. This is especially evident in extremely complex processes like organic matter decomposition, which require the cooperation of multiple organisms bringing together diverse metabolic tools to tackle the numerous metabolic routes involved in this process [6].

Table 1: Key Concepts in Multidimensional Soil Biodiversity-Ecosystem Functioning Relationships

Dimension	Components	Ecosystem Significance
Trophic Levels	Plants, animals, microbes	Determines energy flow, nutrient cycling efficiency, and trophic interactions
Facets	Taxonomic, functional, phylogenetic	Enhances functional redundancy, response diversity, and evolutionary potential
Spatial Scales	Local, regional, landscape	Influences meta-community dynamics, dispersal, and cross-system subsidies
Temporal Scales	Seasonal, interannual, decadal	Affects ecosystem stability, resilience, and adaptive capacity

The Critical Role of Soil Organisms in Ecosystem Processes

Soil organisms operate as a complex, interconnected network that sustains critical ecosystem functions through multiple pathways. Primary decomposers, principally bacteria and fungi, transform soil organic matter into mineral nutrients necessary for plant growth [3]. Many soil microorganisms form symbioses with plant roots—such as nitrogen-fixing bacteria and mycorrhizal fungi—directly influencing plant growth and performance [3]. Larger organisms, including protists, nematodes, and microarthropods, feed on plant roots, organic matter, or smaller organisms, thereby affecting the abundance of decomposer and mutualist communities and releasing mineral nutrients locked in microbial biomass [3]. This intricate web of interactions creates a self-sustaining cycle that maintains ecosystem productivity.

Beyond decomposition and nutrient cycling, soil biodiversity contributes to numerous other ecosystem services. Soils perform five essential functions: regulating water, sustaining plant and animal life, filtering and buffering potential pollutants, cycling nutrients, and providing physical stability and support [26]. The diversity and productivity of living things depend fundamentally on soil, while the minerals and microbes in soil are responsible for filtering, buffering, degrading, immobilizing, and detoxifying organic and inorganic materials [26]. These functions highlight why soil biodiversity conservation is essential for the sustainability of life on Earth.

Global Change Drivers: Land-Use and Climate Change Impacts

Land-Use Change Patterns and Impacts on Soil Biodiversity

Global land use has undergone dramatic transformations in recent decades, with profound implications for soil biodiversity and functioning. Recent analyses reveal that land use change has affected almost a third (32%) of the global land area in just six decades (1960-2019)—approximately four times greater in extent than previously estimated from long-term land change assessments [27]. This conversion has followed divergent patterns across the globe: afforestation and cropland abandonment have characterized the Global North, while deforestation and agricultural expansion have dominated the Global South [27]. These geographically diverging land use change processes create distinctly different pressures on soil biodiversity in different regions.

The temporal dynamics of land use change reveal an acceleration phase with increasing rates of change from 1960 to 2004, followed by a deceleration phase from 2005 to 2019 [27]. This transition correlates with shifting global food regimes and the increasing influence of tele-connected markets, particularly the offshoring of land use change from the Global North to the South [27]. The growing proportion of cropland in Global South countries used for export and consumption outside their territories represents a significant driver of soil biodiversity loss, as natural ecosystems are converted to intensive agriculture [27]. These land use changes disproportionately affect some soil organism groups over others, potentially leading to the impoverishment or complete exclusion of entire functional groups [3].

Visual Summary 1: Conceptual framework showing how global change drivers, including land-use change and climate change, affect soil communities and ecosystem functioning through multiple pathways.

Climate Change Effects on Environmental Suitability

Climate change is altering the environmental suitability of land use and land cover classes globally, with particularly pronounced impacts in semiarid regions. Research from the southern Brazilian semiarid region projects significant shifts in land use suitability under future climate scenarios, with croplands expected to lose approximately 8% of their suitable area under the worst-case scenario (RCP 8.5), while pastures may expand by up to 30% [28]. Areas suitable for savannas are expected to increase under both RCP scenarios, potentially expanding into lands historically occupied by forests, grasslands, and eucalyptus plantations [28]. These projected changes will likely lead to biodiversity loss and socioeconomic disruptions while fundamentally altering soil BEF relationships.

The mechanisms through which climate change affects soil BEF relationships involve both direct and indirect pathways. Direct effects include physiological stresses on soil organisms from increased temperatures and altered moisture regimes, while indirect effects operate through climate-induced vegetation changes that modify resource inputs and habitat structure [28]. In semiarid regions characterized by low water availability, projected substantial increases in air temperature and decreases in precipitation will increase aridity and alter the environmental suitability dynamics of different land use and land cover classes [28]. These changes may disproportionately impact specialized soil taxa with limited dispersal capabilities, potentially disrupting critical ecosystem functions.

Table 2: Projected Changes in Land Use/Land Cover Suitability under Climate Change Scenarios in the Brazilian Semiarid Region

Land Use/Land Cover Class	RCP 2.6 (Optimistic)	RCP 8.5 (Pessimistic)	Primary Drivers
Croplands	Moderate decrease	8% loss of suitable area	Reduced precipitation, increased aridity
Pastures	Moderate expansion	Up to 30% expansion	Replacement of croplands, vegetation changes
Savannas	Increase	Significant increase	Expansion into forest/grassland areas
Forests	Decrease	Decrease	Climate-mediated habitat suitability loss
Eucalyptus Plantations	Decrease	Decrease	Water stress, economic viability

Experimental Evidence: Soil Biodiversity-Functioning Relationships Across Ecosystems

Methodology of a Multi-Grassland Soil Biodiversity Experiment

A sophisticated experimental approach examining soil BEF relationships across five European grasslands with differing soil properties provides compelling evidence for context-dependent biodiversity effects. The study employed a filtered biodiversity gradient design, where soil communities from each grassland were filtered by size to create four biodiversity treatments: (1) High (complete soil community), (2) Medium (excluding macrofauna), (3) Low (excluding macrofauna and mesofauna), and (4) Min (bacteria and fungi only) [3]. This design allowed researchers to test how progressive loss of soil biodiversity groups affects ecosystem functioning across grassland types classified as relatively fertile (rich in organic matter) or relatively poor [3].

The experimental protocol involved collecting soils from five European grasslands (Belgium, two in Hungary, Germany, and Italy) encompassing a range of soil carbon and nitrogen concentrations, pH levels, and soil textures [3]. Sterilized soil substrates from each site were inoculated with size-filtered soil organisms extracted from respective field soils, and also received the plant seed community specific to each grassland [3]. Researchers then measured multiple ecosystem functions, including plant shoot and root biomass, plant diversity, microbial respiration, microbial carbon and nitrogen, and the abundance of functional groups such as nitrifiers, arbuscular mycorrhizal fungi, and plant pathogens [3]. This comprehensive approach enabled direct testing of how soil biodiversity effects vary across environmental contexts.

Key Findings on Biodiversity-Functioning Relationships

The experimental results revealed complex, context-dependent relationships between soil biodiversity and ecosystem functioning. Contrary to the initial hypothesis that soil biodiversity decline would universally reduce ecosystem functions, the study found that increased soil biodiversity may promote some but suppress other specific functions related to plant productivity and nutrient cycling across different grassland soils [3]. This suggests that intricate direct and indirect interactions among soil organisms can result in compensatory increases in certain ecosystem processes at the expense of others, creating trade-offs in functioning.

The majority of soil biodiversity-function relationships did not vary systematically with variation in initial soil properties across the five grasslands [3]. However, the study identified that the presence of particular functional groups, rather than overall diversity, drove specific processes. For instance, the presence of arbuscular mycorrhizal fungi was particularly important for plant phosphorus uptake, especially in nutrient-poor conditions [3]. Additionally, the reduction in predator pressure (microfauna) in the lower diversity treatments resulted in the accumulation of fungal plant pathogens and bacterial parasites, demonstrating how trophic interactions influence ecosystem health [3].

Visual Summary 2: Experimental workflow of a multi-grassland soil biodiversity experiment testing how filtered soil communities affect ecosystem functions across different environmental contexts.

Context Dependency: Soil Properties and Environmental Conditions Modulate BEF Relationships

The Role of Soil Fertility in Mediating BEF Relationships

Soil fertility emerges as a critical factor modulating the relationship between soil biodiversity and ecosystem functioning. Experimental evidence indicates that fertile soils with higher organic matter and nutrient content can support more diverse soil communities, harboring a variety of different taxa capable of performing the same function [3]. This creates higher levels of functional redundancy, particularly for general functions such as nutrient mineralisation, making fertile soils less sensitive to biodiversity decline compared to poorer soils [3]. This finding challenges universal predictions about biodiversity effects and highlights the need to consider environmental context.

The dependence of ecosystem functioning on specific soil functional groups varies with soil fertility. In nutrient-poor soils, functioning may be largely dependent on the presence and richness of particular soil functional groups, making these systems more vulnerable to biodiversity loss [3]. This aligns with broader patterns showing that soil biodiversity is particularly critical for supporting function in drylands and poorly developed soils, which have small stocks of nutrients and organic matter [6]. In these systems, nutrient inputs depend heavily on the daily contribution of diverse soil microbiota to decompose and depolymerize litter and organic matter [6].

The Importance of Multidimensional Biodiversity

The conservation of multidimensional biodiversity—encompassing multiple trophic levels, taxonomic groups, and spatial scales—proves essential for maintaining ecosystem functioning under global change. Research demonstrates that multidimensional biodiversity regulates the response of ecosystem functioning to global change factors, with high levels of multidimensional biodiversity mitigating the negative impacts of global change on ecosystem functioning [25]. This buffering capacity arises from greater functional redundancy, response diversity, and complementary resource use across diverse biological communities.

The interaction of multiple global change factors may lead to greater reductions in biodiversity and ecosystem functioning than single global change factors alone [25]. This synergistic effect underscores the importance of comprehensive approaches to biodiversity conservation that address the full suite of anthropogenic pressures. Moreover, different dimensions of biodiversity may exhibit varying sensitivities to global change drivers, creating complex response patterns that necessitate multidimensional assessment frameworks [25]. Prioritizing conservation efforts to maintain and enhance all dimensions of biodiversity is therefore essential to address the challenges of future global change.

Management Implications and Conservation Strategies

Soil Health Management Principles

Effective soil management can enhance biodiversity and buffer ecosystems against global change impacts. The Natural Resources Conservation Service outlines four key principles for soil health management: (1) maximize presence of living roots, (2) minimize disturbance, (3) maximize soil cover, and (4) maximize biodiversity [26]. These principles work synergistically to create conditions that support diverse soil communities and enhance ecosystem functioning. For instance, living roots maintain a rhizosphere—an area of concentrated microbial activity—where peak nutrient and water cycling occurs [26].

Disturbance minimization, particularly through reduced tillage practices, preserves soil organic matter and structure while protecting habitat for soil organisms [26]. Tillage can destroy soil organic matter and structure, reduce water infiltration, increase runoff, and make soil more susceptible to erosion [26]. By contrast, no-till systems maintain more organic matter and moisture for plant use, support better nutrient cycling and root growth, and improve carbon sequestration efficiency [26]. These management approaches create more favorable conditions for diverse soil communities to persist and maintain ecosystem functions.

Biodiversity Conservation in Agricultural Landscapes

Maintaining biodiversity in agricultural landscapes requires specific strategies that support soil organisms while meeting production needs. Diversified crop rotations and cover crop mixtures enhance biodiversity above and below ground, helping to prevent disease and pest problems associated with monocultures [26]. Cover crops provide multiple benefits, including increased organic matter, improved water infiltration, natural nutrient cycling, and protection against erosion [26]. These practices collectively support more diverse soil food webs that provide for nutrient, energy, and water cycling, allowing soils to express their full potential.

Integrating conservation strategies with agricultural production represents a promising pathway for sustaining soil BEF relationships in human-dominated landscapes. Lack of biodiversity severely limits the potential of any cropping system and increases disease and pest problems [26]. Ultimately, biodiversity is the key to the success of any agricultural system, creating resilient, productive soils capable of withstanding environmental fluctuations [26]. By implementing soil health management systems, farmers can increase organic matter, improve microbial activity, sequester more carbon, enhance water infiltration, and improve wildlife habitat—all while maintaining productivity and profitability [26].

Table 3: Research Reagent Solutions for Soil Biodiversity-Ecosystem Functioning Studies

Research Tool Category	Specific Methods/Approaches	Application in Soil BEF Research
Molecular Analysis	DNA/RNA extraction, amplicon sequencing, metagenomics	Characterizing microbial, fungal, protist diversity and functional potential
Stable Isotope Probing	13C, 15N, 18O labeling	Tracing nutrient flows, identifying active microbial taxa
Microcosm/Mesocosm Systems	Size-filtered inocula, sterilization techniques	Experimental manipulation of soil community complexity
Enzyme Assays	Fluorometric substrates, colorimetric detection	Measuring functional potential for nutrient cycling
Respiratory Measurements	Microrespiration systems, gas chromatography	Quantifying microbial metabolic activity
Microscopy Techniques	Fluorescence in situ hybridization (FISH)	Visualizing spatial organization of microbial communities
Physical Fractionation	Density separation, size filtering	Isolating soil fractions and organism size classes

Knowledge Gaps and Future Research Directions

Despite significant advances in understanding soil BEF relationships under global change, critical knowledge gaps remain. Essential unresolved questions include how different dimensions of biodiversity (taxonomic, functional, phylogenetic) contribute to ecosystem stability across environmental gradients, and how global change factors interact to affect soil biodiversity and functioning [6]. Understanding these interactions is crucial for predicting ecosystem responses to simultaneous stressors and developing effective conservation strategies.

Future research should prioritize several key areas: (1) assessing the effects of multiple interacting global change factors on soil BEF relationships; (2) understanding the role of underrepresented organism groups (e.g., viruses) in ecosystem processes; (3) elucidating the mechanisms underlying context-dependent biodiversity effects; and (4) developing manipulative experiments that simulate future global change scenarios [6] [25]. Addressing these knowledge gaps will enhance our ability to conserve soil biodiversity and maintain ecosystem functions in the face of ongoing global change.

Long-term experimental platforms, such as the Jena Experiment—one of the longest-running biodiversity experiments in the world (running since 2002)—provide valuable insights into the mechanisms that determine BEF relationships over extended timeframes [29]. These long-term studies reveal that the relationship between plant diversity and several ecosystem functions strengthens over time, suggesting that both selection effects and complementary resource use contribute to biodiversity effects [29]. Extending such long-term research to include more comprehensive soil biodiversity assessments and global change manipulations will further enhance our understanding of soil BEF dynamics.

From Theory to Practice: Methodological Approaches for Quantifying Soil BEF Relationships

Understanding the intricate relationships between biodiversity and ecosystem functioning (BEF) represents one of the most pressing challenges in soil ecology. Researchers employ sophisticated experimental designs to unravel how the composition and diversity of soil organisms drive essential processes like nutrient cycling, carbon sequestration, and primary production. This guide compares two pivotal approaches in modern soil research: size-fractionation inoculation experiments, which manipulate soil community complexity, and long-term field experiments, which observe ecosystem development over extended periods. Each methodology offers distinct advantages and limitations, providing complementary insights into the mechanisms sustaining terrestrial ecosystems.

Methodological Comparison: Core Experimental Approaches

The investigation of soil biodiversity-ecosystem functioning relationships relies on carefully controlled manipulations and observational frameworks. The table below summarizes the fundamental characteristics of two primary experimental designs.

Table 1: Comparison of Core Experimental Designs in Soil BEF Research

Experimental Design	Key Manipulation	Primary Objective	Typical Scale	Temporal Scope
Size-Fractionation Inoculation	Physical sieving of soil communities to create biodiversity gradients [20]	Isolate the effect of soil biodiversity from abiotic factors; test causality [20]	Microcosms/Laboratory [20]	Short to Medium-term (weeks to months) [20]
Long-Term Field Experiments	Establishment of plots with varying plant diversity and management [30] [29]	Understand temporal development of BEF relationships under realistic conditions [30]	Field Plots (square meters to hectares) [31] [29]	Long-term (years to decades) [30] [29]

Size-Fractionation Inoculation Protocol

The size-fractionation inoculation method mechanically reduces soil biodiversity through sequential filtration, creating a gradient of soil community complexity for establishing causal relationships.

• Step 1: Soil Collection: Collect bulk soil from a representative field site, ensuring appropriate depth and habitat representation.
• Step 2: Fractionation Procedure: Process soil through a series of sieves with decreasing mesh sizes [20]:
  • >5,000 μm: Retains large fauna, organic debris, and most soil organisms
  • 250 μm: Removes larger mesofauna but retains microbes, microfauna, and some smaller mesofauna
  • 50 μm: Eliminates most nematodes and microarthropods, retaining primarily microbial communities
  • <20 μm: Greatly reduces microbial diversity, creating highly simplified communities [20]
• Step 3: Inoculation: Use the size-fractionated soil as inoculum for sterile soil in experimental microcosms [20].
• Step 4: Ecosystem Monitoring: Measure response variables including plant diversity, nutrient leaching, decomposition rates, and ecosystem multifunctionality [20].

This approach successfully creates broad soil biodiversity gradients, with some groups (e.g., nematodes, mycorrhizal fungi) entirely eliminated in finer fractions, while fungal and bacterial communities show reduced abundance and richness [20].

Long-Term Field Experiment Establishment

Long-term field experiments investigate BEF relationships under realistic environmental conditions across extended timeframes, revealing temporal dynamics inaccessible through short-term studies.

• Step 1: Site Selection: Choose representative ecosystem with appropriate size (e.g., 2.8 hectares for mining reclamation studies [31]).
• Step 2: Experimental Design:
  • Implement random block designs with multiple diversity treatments [30]
  • Include gradient of species richness (e.g., monocultures to 16-species mixtures [30])
  • Incorporate various planting configurations (monocultures, specific species combinations [31])
• Step 3: Ongoing Maintenance and Monitoring:
  • Conduct regular censuses of species composition, survival rates, and growth parameters [31]
  • Measure ecosystem functions (e.g., aboveground net primary productivity, nutrient cycling, soil carbon sequestration) at regular intervals [30]
  • Document environmental conditions and management practices consistently

The Jena Experiment, established in 2002, exemplifies this approach, examining how plant diversity affects ecosystem processes in Central European grassland [30] [29].

Quantitative Outcomes and Temporal Dynamics

Each experimental approach generates distinct quantitative insights into BEF relationships, particularly regarding the strength and mechanisms of biodiversity effects.

Table 2: Quantitative Outcomes from Different Experimental Approaches

Measured Parameter	Size-Fractionation Findings	Long-Term Field Experiment Findings	Key Implications
Ecosystem Multifunctionality	Strong positive relationship with soil biodiversity index (explaining major variation) [20]	Strengthening positive relationship with species richness over 17 years [30]	Soil biodiversity drives multiple functions; effects strengthen over time
Plant Diversity	Declined strongly with soil biodiversity reduction [20]	Positive relationship with species richness maintained long-term [30]	Aboveground-belowground diversity linkages are crucial
Nutrient Retention	Phosphorus leaching increased threefold in most simplified communities [20]	Complementarity effects increased over time, enhancing nutrient use [30]	Diverse communities improve nutrient conservation
Temporal Stability	Not directly measured in short-term studies	Community stability increased with diversity; effect strengthened over time [30]	Biodiversity provides insurance against environmental fluctuations

Visualizing Experimental Workflows

The diagrams below illustrate the key procedural pathways for both experimental approaches.

Size-Fractionation Inoculation Workflow: This reductionist approach creates biodiversity gradients through sequential filtration.

Long-Term Field Experiment Workflow: This approach monitors ecosystem development across extended timescales.

The Scientist's Toolkit: Essential Research Reagents and Materials

Successful implementation of soil BEF research requires specialized materials and reagents for manipulating and measuring soil communities.

Table 3: Essential Research Reagents and Materials for Soil BEF Experiments

Category	Specific Items	Research Function	Application Context
Soil Fractionation	Sieve series (5000μm, 250μm, 50μm, 20μm), Low-energy sonication equipment, Centrifuge [32] [20]	Create soil biodiversity gradients via physical separation	Size-fractionation inoculation experiments
Molecular Analysis	DNA extraction kits, PCR reagents, 16S/ITS primers, Sequencing platforms [2] [33]	Characterize microbial community composition and diversity	Both approaches; essential for soil microbiome analysis
Field Equipment	Soil corers, Root ingrowth bags, Litter traps, Gas flux chambers [30]	Measure ecosystem processes under field conditions	Long-term field experiments
Microscopy	Fluorescence stains (e.g., WGA-FITC for fungi), Epifluorescence microscope [33]	Quantify fungal colonization and microbial biomass	Verification of inoculation success
Chemical Analysis	PLFA/NLFA reagents, Elemental analyzers, Nutrient extraction solutions [2]	Assess microbial biomass and nutrient cycling processes	Both experimental approaches

Key Research Insights and Applications

Complementary Evidence for Biodiversity Effects

Both experimental approaches provide strong, complementary evidence that soil biodiversity drives ecosystem functioning, though they reveal different aspects of this relationship:

• Size-fractionation studies demonstrate that simplified soil communities show impaired ecosystem functioning, with reduced plant diversity, increased nutrient leaching, and decreased decomposition rates [20]. The relationship between soil biodiversity and ecosystem multifunctionality is strikingly positive and linear [20].

• Long-term experiments reveal that biodiversity-ecosystem functioning relationships strengthen over time [30]. In the Jena Experiment, the positive effect of species richness on productivity significantly increased over 17 years, with diverse communities maintaining more stable productivity [30].

Divergent Insights into Mechanisms

The different approaches illuminate complementary mechanisms underlying BEF relationships:

• Size-fractionation experiments highlight the importance of specific functional groups. The loss of key groups like mycorrhizal fungi and nematodes corresponds with abrupt shifts in ecosystem functioning [20].

• Long-term field experiments demonstrate how temporal mechanisms develop over time. Species asynchrony and complementarity take more than a decade to develop strong stabilizing effects on ecosystem functioning [30].

Predictive Applications in Restoration

Mining restoration research demonstrates how BEF principles guide effective intervention strategies:

• Species compatibility critically determines restoration success. In Pingshuo open-pit mine reclamation, oil pine and locust exhibited mutually beneficial interactions, while combinations involving sea buckthorn and caragana microphylla showed competitive suppression [31].

• Functional trait complementarity enhances ecosystem development. The locust-sea buckthorn configuration emerged as a strategic model for promoting multi-species growth in degraded landscapes [31].

Size-fractionation inoculation and long-term field experiments represent complementary methodological paradigms in soil BEF research. The reductionist approach of size-fractionation enables precise identification of causal mechanisms by isolating biodiversity effects from confounding environmental variables. Conversely, long-term field experiments capture the temporal dynamics and real-world complexity of ecosystem development, revealing patterns inaccessible in short-term studies.

Together, these approaches demonstrate that soil biodiversity is not merely a response variable but a fundamental driver of ecosystem processes including nutrient cycling, productivity, and stability. Their integration provides a robust evidence base predicting that ongoing losses of soil biodiversity will impair multiple ecosystem functions, with consequences for food security, climate regulation, and ecosystem sustainability [15] [20] [34]. Future research combining mechanistic isolation with ecological realism will further illuminate the processes maintaining terrestrial ecosystems and guide effective restoration of degraded soils worldwide.

Understanding the complex relationships between biodiversity and ecosystem functioning is a central goal in soil ecological research. The development of omics technologies has provided scientists with powerful tools to probe the vast diversity of soil microbial and faunal communities, which were previously difficult to study using traditional culture-based methods. Among these tools, metagenomics and metabarcoding have emerged as the two primary high-throughput sequencing approaches for profiling soil communities. While both methods extract and sequence DNA directly from environmental samples, they differ fundamentally in their scope, biases, and applications. This guide provides an objective comparison of these technologies, supported by experimental data, to help researchers select the appropriate method for testing specific biodiversity-ecosystem functioning relationships in soils.

Metabarcoding and metagenomics represent distinct approaches for assessing biodiversity in soil communities, each with characteristic workflows and analytical outputs.

Metabarcoding employs polymerase chain reaction (PCR) to amplify specific marker genes from environmental DNA that serve as "barcodes" for taxonomic identification. For bacteria and archaea, the 16S ribosomal RNA (rRNA) gene is predominantly targeted, with universal primers typically amplifying hypervariable regions (e.g., V3-V4) [35]. For fungi, the internal transcribed spacer (ITS) region is the standard marker, while the cytochrome c oxidase subunit I (COI) gene is commonly used for soil fauna such as mites and springtails [36] [37]. This targeted approach allows for cost-effective sequencing and provides data primarily suited for taxonomic profiling and estimating relative abundances of different organismal groups.

In contrast, metagenomics involves the random sequencing of all DNA fragments present in an environmental sample without a prior PCR amplification step [38]. This method captures genomic information from all domains of life simultaneously—bacteria, archaea, fungi, protists, and microfauna—and enables functional profiling by identifying genes involved in biogeochemical cycling and other metabolic processes [39] [40]. While computationally more intensive and requiring greater sequencing depth, metagenomics provides a more comprehensive view of both the taxonomic composition and functional potential of soil communities.

The following diagram illustrates the key procedural differences between these two approaches:

Performance Comparison: Experimental Data and Technical Considerations

Multiple studies have directly compared the performance of metagenomics and metabarcoding for assessing soil biodiversity. The table below summarizes key comparative findings from empirical studies:

Table 1: Comparative performance of metagenomics and metabarcoding for soil community profiling

Performance Metric	Metabarcoding	Metagenomics	Experimental Evidence
Taxonomic Coverage	Limited to taxa amplified by primers	Recoveres additional taxa missed by primers	Metagenomics identified large number of additional OTUs absent in metabarcoding [38]
Primer Bias	High: preferential amplification of certain taxa	None: no PCR step	Metabarcoding can dramatically change community composition even at lower taxonomic levels [38]
Detection of Rare Taxa	Lower sensitivity due to amplification bias	Higher sensitivity with sufficient sequencing depth	Metagenomics allows detection of low-abundance community members [38]
Taxonomic Resolution	Varies with marker gene; species-level often challenging	Potentially higher with sufficient genomic coverage	PacBio long-read metabarcoding provided more reliable community profiles than Illumina [36]
Functional Profiling	Limited to inference from taxonomy	Direct assessment of functional genes	Metagenomics enables analysis of genes encoding enzymes in biogeochemical cycling [39] [40]
Sequencing Depth	Lower requirement	Higher requirement for comparable coverage	Classification rates: Metabarcoding (77-97%) vs Metagenomics (17-32%) with same database [38]

Beyond these technical comparisons, the choice between methods also depends on the specific research questions. A novel Two-Step Metabarcoding (TSM) approach has been developed to address primer bias issues in standard metabarcoding. This method combines initial sequencing with universal 16S rDNA primers to identify dominant taxa, followed by a second sequencing step using taxa-specific primers designed for the most abundant phyla [35]. This approach provides more detailed resolution within dominant taxonomic groups while remaining more cost-effective than metagenomics.

For ecosystem functioning studies, metagenomics offers distinct advantages by enabling direct assessment of functional genes. Research on successional gradients has demonstrated the capability of metagenomics to track changes in functional genes encoding enzymes involved in carbon (C), nitrogen (N), and phosphorus (P) cycling [39]. This approach revealed that grassland afforestation increases functional diversity while decreasing taxonomic diversity, highlighting a decoupling between taxonomic and functional measures of biodiversity [39].

Experimental Protocols for Soil Community Profiling

Metabarcoding Protocol for Soil Microbes

The following protocol outlines a standardized approach for metabarcoding analysis of soil bacterial and fungal communities:

• Soil Sampling and DNA Extraction: Collect soil cores from the study site using sterile corers. For homogeneous representation, combine multiple subsamples from each sampling location. Store samples immediately at -80°C until DNA extraction. Extract DNA using commercial soil DNA extraction kits (e.g., FastDNA SPIN Kit or DNeasy PowerSoil Pro Kit), following manufacturer's instructions with bead-beating step for efficient cell lysis [35] [36]. Quantify DNA concentration and quality using fluorometry and gel electrophoresis.

• PCR Amplification and Library Preparation: For bacteria, amplify the V3-V4 region of the 16S rRNA gene using primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) [41]. For fungi, target the ITS1 region using primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′) [41]. Perform PCR in triplicate reactions to minimize stochastic amplification bias, then pool products. Use a high-fidelity polymerase to reduce amplification errors. Index PCR adds dual indices and sequencing adapters for multiplexing.

• Sequencing and Bioinformatic Analysis: Sequence amplified libraries on Illumina MiSeq or HiSeq platforms with 2×250 bp or 2×300 bp paired-end reads. Process raw sequences using QIIME 2 or DADA2 pipelines to demultiplex, quality filter, denoise, and remove chimeras [41]. Cluster sequences into amplicon sequence variants (ASVs) or operational taxonomic units (OTUs) at 100% similarity. Assign taxonomy using reference databases (Silva for bacteria, UNITE for fungi) [41].

Metagenomic Protocol for Soil Communities

• DNA Extraction and Quality Control: Extract high-quality, high-molecular-weight DNA from 2-10g of soil using kits designed for complex environmental samples (e.g., PowerMax Soil Kit) [36]. Include additional purification steps to remove humic acids and other PCR inhibitors. Verify DNA integrity by pulse-field gel electrophoresis and quantify using fluorometric methods suitable for complex mixtures.

• Library Preparation and Sequencing: Fragment DNA to desired size (typically 350-500 bp) using acoustic shearing or enzymatic fragmentation. Prepare sequencing libraries using Illumina TruSeq, NovaSeq, or PacBio SMRTbell kits following manufacturer's protocols. For comprehensive community profiling, aim for minimum sequencing depth of 20-50 million reads per sample for complex soil communities [38]. Consider long-read sequencing technologies (PacBio, Oxford Nanopore) for improved assembly and taxonomic classification.

• Bioinformatic Processing and Analysis: Quality filter raw reads using Trimmomatic or similar tools. For taxonomic profiling, use Kraken 2 or MetaPhlAn with customized databases [38]. For functional analysis, assemble reads into contigs using metaSPAdes or MEGAHIT, then predict open reading frames and annotate against functional databases (Subsystems, KEGG, COG, CAZy) [40]. Calculate functional diversity metrics and map genes to specific biogeochemical pathways.

Applications in Biodiversity-Ecosystem Functioning Research

The relationship between soil biodiversity and ecosystem functioning represents a critical research area where omics technologies provide unique insights. Experimental evidence demonstrates that soil biodiversity enhances the temporal stability of multiple ecosystem functions, including plant productivity, litter decomposition, and soil carbon assimilation [42]. This stabilizing effect occurs through asynchrony among microbial taxa, where different fungi and bacteria promote different functions at different times [42].

Metagenomic approaches have been particularly valuable for documenting functional redundancy—the phenomenon where multiple taxa perform similar ecosystem functions—in soil microbial communities. Research shows that despite significant taxonomic diversity loss (up to 72% species reduction), the relative abundance of most functional categories remains remarkably stable [40]. This functional resilience demonstrates a decoupling between taxonomy and function in soil ecosystems, with implications for ecosystem management and conservation.

Metabarcoding studies have revealed how land-use changes affect soil communities and their functions. For example, research on grassland afforestation successional gradients showed threshold dynamics in soil microbial communities, with abrupt decreases in bacterial diversity and marked shifts in community composition [39]. These taxonomic changes were linked to alterations in functional genes involved in nutrient cycling, highlighting the complex relationships between biodiversity and ecosystem processes.

For soil fauna, metabarcoding approaches have enabled broader assessment of microarthropod communities (mites and springtails) that function as decomposers, bacterivores, fungivores, and carnivores in soil food webs [37]. Comparisons between morphological identification and metabarcoding show that while metabarcoding may detect a different subset of taxa, it reveals comparable ecological patterns in response to environmental changes such as irrigation treatments [37].

Essential Research Reagent Solutions

The table below outlines key reagents and kits used in metagenomics and metabarcoding workflows for soil community profiling:

Table 2: Essential research reagents and materials for soil omics studies

Reagent/Kits	Application	Function	Examples/Alternatives
Soil DNA Extraction Kits	Both methods	Isolation of high-quality DNA from complex soil matrices	DNeasy PowerSoil Pro Kit (0.2g samples) [36], PowerMax Soil Kit (2g samples) [36], FastDNA SPIN Kit [35]
PCR Master Mix	Metabarcoding	Amplification of target marker genes	HOT FIREPol Blend Master Mix [36]
16S rRNA Primers	Metabarcoding	Amplification of bacterial/archaeal marker genes	515F/806R for V4 region [41], 341F/785R for V3-V4 region [35]
ITS Primers	Metabarcoding	Amplification of fungal marker genes	ITS1F/ITS2 [41]
COI Primers	Metabarcoding	Amplification of animal/fauna marker genes	LCO1490/HCO2198 [36], mlCOIintF/jgHCO2198 [36]
Library Prep Kits	Both methods	Preparation of sequencing libraries	Illumina Nextera XT, TruSeq DNA PCR-Free, PacBio SMRTbell
Sequence Databases	Both methods	Taxonomic and functional annotation	SILVA (rRNA genes) [41], UNITE (fungal ITS) [41], Greengenes [35], RefSeq [40], Subsystems [40]

Metagenomics and metabarcoding offer complementary approaches for profiling soil communities in biodiversity-ecosystem functioning research. Metabarcoding provides a cost-effective method for taxonomic profiling of specific organismal groups, while metagenomics delivers comprehensive insights into both taxonomic composition and functional potential without amplification biases. The choice between these methods should be guided by research questions, resources, and specific aspects of biodiversity-ecosystem functioning relationships under investigation. As sequencing technologies continue to advance and decrease in cost, metagenomics is likely to become more accessible for routine analysis of soil communities, particularly for studies requiring functional insights. However, metabarcoding will remain valuable for large-scale surveys targeting specific taxonomic groups or when resources are limited. Methodological innovations such as the two-step metabarcoding approach and long-read sequencing continue to enhance the resolution and accuracy of both techniques, promising new discoveries in the complex relationships between soil biodiversity and ecosystem functioning.

Understanding and quantifying ecosystem functions is a paramount challenge in soil ecology, particularly for testing biodiversity-ecosystem functioning (BEF) relationships in complex natural systems [43]. Soils harbor a substantial fraction of the world's biodiversity, contributing to many crucial ecosystem functions and services, including climate regulation, nutrient cycling, and food production [10]. However, major gaps persist in our understanding of how to accurately measure and link the diversity of soil organisms to the multiple functions they support [10] [44].

The field currently faces a significant limitation: the lack of a common measurement framework that enables systematic comparison across different soils and ecosystem types [44]. This comparison is hindered by the existence of numerous measurement methods and indicators across research groups. For instance, a review of 65 soil health monitoring schemes revealed that only three indicator variables (soil organic matter, acidity, and available phosphorus) were commonly measured, while no biological measure was implemented in more than 30% of the schemes [44]. This status quo prevents robust comparisons and undermines efforts to establish general principles in BEF relationships.

This guide objectively compares current methodologies for quantifying key ecosystem functions, including biomass, respiration, and nutrient fluxes, while framing them within the broader context of assessing soil multifunctionality. By providing standardized protocols and comparative data, we aim to support researchers in selecting appropriate methodologies for their specific BEF research questions.

Theoretical Frameworks for Ecosystem Function Assessment

The Three Major Axes of Terrestrial Ecosystem Function

Recent research has revealed that most variability in ecosystem functions (71.8%) can be captured by three key axes derived from surface gas exchange measurements across major terrestrial biomes [45]. Understanding these axes provides a crucial framework for designing targeted measurement approaches in BEF research.

• Axis 1: Maximum Ecosystem Productivity - This primary axis explains 39.3% of variance and is dominated by maximum gross CO2 uptake (GPPsat) and maximum net ecosystem productivity (NEPmax). It is largely explained by vegetation structure (leaf area index, above-ground biomass, canopy height) and atmospheric aridity (vapor pressure deficit) [45]. This axis runs from cold, arid shrublands with low productivity to forests in continental, tropical, and temperate climates with high photosynthesis rates.

• Axis 2: Ecosystem Water-Use Strategies - Explaining 21.4% of variance, this axis reflects water-use efficiency metrics and evaporative fraction. It is jointly explained by variation in vegetation height and climate, running from temperate forests and dry subtropical sites with higher water-use efficiency to cold or tropical climates and wetlands with high evaporative fraction [45].

• Axis 3: Ecosystem Carbon-Use Efficiency - This axis explains 11.1% of variance and is dominated by apparent carbon-use efficiency (aCUE) and basal ecosystem respiration. It features a gradient related to aridity and temperature, running from Arctic and boreal sites with high carbon-use efficiency to hot, dry climates with lower efficiency [45].

Energy Flux as an Integrative Framework

An emerging framework proposes energy flux dynamics in food webs as a universal tool for understanding BEF relationships in multitrophic systems [43]. This approach integrates concepts from network theory, metabolic ecology, BEF theory, and ecological stoichiometry to provide a common currency for ecosystem functioning across trophic levels and ecosystem types.

Energy flux characterizes the rate of energy flow among nodes in food webs, expressing energy consumption by different trophic groups and describing the energetic structure of communities [43]. This approach is particularly valuable because fluxes reflect most of the functions commonly measured in biodiversity experiments and could provide a standardized index of ecosystem multifunctionality.

Table 1: Comparative Advantages of Theoretical Frameworks for BEF Research

Framework	Key Variables Measured	Spatial Scale Applicability	Integration of Multitrophic Diversity	Implementation Complexity
Three-Axis Framework [45]	Gas exchange (CO2, H2O); Vegetation structure; Climate variables	Ecosystem to landscape	Limited	Moderate
Energy Flux Framework [43]	Trophic interactions; Metabolic demands; Biomass stocks	Local to regional	Excellent	High
Latent-Variable Modeling [44]	Multiple soil biological, chemical, and physical indicators	Plot to ecosystem	Moderate	High

Methodological Approaches for Measuring Key Ecosystem Functions

Quantifying Soil Carbon Stocks and Greenhouse Gas Fluxes

A standardized five-stage protocol has been developed for systematically measuring soil carbon stocks and greenhouse gas fluxes in complex environments with heterogeneous land uses, soil types, and management practices [46]. This approach is particularly valuable for BEF research as it enables comparable measurements across different biodiversity treatments or sampling sites.

The protocol employs the SCORPAN framework (Soil, Climate, Organisms, Relief, Parent material, Age, and Spatial position) to account for environmental covariates that influence soil formation and properties [46]. The five stages include: (1) planning and developing a soil sampling design, (2) conducting a time-zero baseline assessment, (3) strategic selection of sampling sites, (4) soil sampling and GHG flux measurements, and (5) data analysis and monitoring design.

For GHG flux measurements, the flux gradient (FG) method enables near-continuous measurements of trace gas exchanges at multiple plots with a single laser spectrometer, facilitating treatment comparison and replication in experimental designs [47]. This method has been successfully deployed to address "pollution swapping" concerns in agricultural management, such as assessing whether practices that reduce nitrate leaching might inadvertently increase emissions of potent greenhouse gases like N2O and CH4 [47].

Partitioning Soil Respiration Components

Understanding the different components of soil respiration is crucial for BEF studies, as autotrophic (root-driven) and heterotrophic (microbe-driven) respiration may respond differently to environmental changes and biodiversity manipulations [48].

The root exclusion method (using bare plots or trenching) is commonly employed to partition soil respiration (RS) into autotrophic (RA) and heterotrophic (RH) components [48]. However, this approach carries an important methodological assumption that must be validated: microbial respiration in root-free plots should be similar to microbial respiration in planted plots. Recent research challenges this assumption, showing that microbial respiration in planted plots can be nearly four times higher than in root-free plots due to rhizosphere stimulation effects [48].

Temperature sensitivity (Q10) differs between respiration components, with RH generally exhibiting greater sensitivity (Q10 = 2.67) than RS (Q10 = 2.29) [48]. This has important implications for BEF research under climate change scenarios, as carbon cycling may be more sensitive to temperature changes than previously assumed based on bulk respiration measurements alone.

Table 2: Methodological Comparison for Soil Respiration Partitioning

Method	Key Principle	Critical Assumptions	Advantages	Limitations
Root Exclusion (Trenching) [48]	Physical separation of roots from soil	Microbial respiration similar in rooted and root-free soils; No root regrowth	Direct field measurement; Long-term monitoring possible	Disturbs soil structure; Rhizosphere priming effects ignored
Bare Soil Plots [48]	Elimination of living roots through vegetation removal	Same as trenching method; Minimal edge effects	Less destructive than trenching; Suitable for annual crops	Limited to systems where vegetation can be permanently excluded
Isotopic Approaches	Use of 13C or 14C labels to trace root-derived C	Isotopic label does not alter microbial processes; Homogeneous label distribution	Non-destructive; Can distinguish recent vs. old C sources	Technically complex; Expensive; Requires specialized equipment

Advanced Sensing Technologies for Nutrient Fluxes

Smart sensing technologies are revolutionizing our ability to monitor nutrient fluxes in both open-field and controlled cultivation systems [49]. These technologies enable real-time, high-resolution measurements that capture the spatial and temporal variability crucial for understanding BEF relationships.

Electrochemical sensors measure specific ions (N, P, K, S) by detecting electrical conductivity or potential changes, while optical sensors using near-infrared (NIR) or mid-infrared (MIR) spectroscopy provide non-invasive analysis of soil properties including organic matter, moisture, and nutrient content [49]. These sensors are increasingly deployed on integrated platforms including handheld devices for spot-checking, drones for large-area mapping, and tractors or autonomous vehicles for continuous, real-time data collection.

For hydroponic systems or soil solution monitoring, sensors such as pH meters, electrical conductivity (EC) sensors, and ion-selective electrodes (ISE) enable continuous monitoring of nutrient levels [49]. When integrated with automated dosing systems, these sensors can adjust nutrient concentrations in real-time, providing precise experimental control for BEF research investigating nutrient cycling processes.

Assessing Ecosystem Multifunctionality: Integrated Approaches

The Challenge of Multifunctionality Quantification

Ecosystem multifunctionality reflects the simultaneous provisioning of multiple ecosystem functions [43]. Measuring this integrated concept presents significant statistical and methodological challenges in BEF research, particularly in soil ecosystems where multiple functions operate across different spatial and temporal scales.

A critical limitation in current soil macroecological studies is the extreme scarcity of sites with co-located biodiversity and ecosystem function data—only 0.3% of all sampling sites have both information about biodiversity and function [10]. This data gap severely constrains our ability to explore general patterns in soil BEF relationships at large spatial scales.

Novel Approaches: Latent-Variable Modeling

To address the challenge of multifunctionality assessment, a novel approach proposes using latent-variable modeling to develop a common "IQ test for soils" [44]. This method, based on factor analysis with a history in social and economic sciences, treats soil functions as latent variables—complex processes that cannot be measured directly but can be inferred from measurable indicators.

The framework formally separates the causes, consequences, predictors, and indicators of soil functioning, linking them to underlying processes and environmental context [44]. This approach represents a significant advancement over current practices where drivers of soil functions (e.g., nutrient content as a proxy for fertility) are often used implicitly as incomplete proxies for the functions themselves.

Energy Flux as a Multifunctionality Metric

The quantification of energy fluxes in food webs provides an alternative approach to assessing ecosystem multifunctionality [43]. This method calculates energy flux between network nodes using measured biomass stocks, energetic expenditure, and ecological efficiencies to balance energetic demands of biomass stocks with energy outflow.

The energy flux approach can be implemented using the equation:

F = (1/e) · (X + L)

Where F is energy flux to each consumer node, e is diet-specific assimilation efficiency, X is the summed metabolic demands of individuals in a consumer node, and L is loss of energy to higher trophic levels [43]. This approach can incorporate individual-level metabolic demands into network nodes, accounting for community composition, body size structure, trophic topology, and temperature effects on metabolism and assimilation efficiency.

The Scientist's Toolkit: Essential Reagents and Materials

Table 3: Essential Research Reagent Solutions for Ecosystem Function Quantification

Reagent/Material	Primary Function	Application Context	Key Considerations
Laser Spectrometers	Quantification of N2O and CH4 concentrations	Flux gradient method for GHG measurements [47]	Requires calibration with standard gases; Sensitivity to environmental conditions
Ion-Selective Electrodes	Detection of specific nutrient ions (NO3-, NH4+, K+)	Soil solution nutrient monitoring [49]	Ion interference can affect accuracy; Regular calibration required
Infrared Gas Analyzers	Measurement of CO2 and H2O fluxes	Soil respiration chambers; Ecosystem-scale flux towers [45] [48]	Requires temperature and pressure corrections; Sensitive to flow rates
Soil Cores	Collection of intact soil samples	Soil carbon stock assessment [46]	Diameter and depth must be standardized; Preservation methods affect microbial activity
Root Exclusion Materials	Physical barriers to root growth (mesh, trenches)	Separation of autotrophic and heterotrophic respiration [48]	Mesh size critical for excluding roots but allowing mycorrhizal hyphae; Trenching disturbs soil structure
Chemical Extractants (K2SO4, KCl, etc.)	Solubilization of soil nutrients and organic carbon	Microbial biomass; Available nutrient pools [48]	Extraction efficiency varies by soil type; Filter pore size affects results
Isotopic Tracers (13C, 15N, 18O)	Tracing element pathways through ecosystems	Partitioning respiration sources; Nutrient uptake studies	Requires mass spectrometry; Homogeneous labeling challenging in field studies

Quantifying ecosystem functions—including biomass, respiration, nutrient fluxes, and multifunctionality—requires careful methodological selection tailored to specific research questions in biodiversity-ecosystem functioning relationships. The comparative analysis presented in this guide reveals that while no single approach perfectly captures the complexity of soil ecological processes, emerging frameworks like the three-axis model [45], energy flux quantification [43], and latent-variable modeling [44] offer promising avenues for standardization.

Critical methodological considerations include: (1) the spatial and temporal mismatch between biodiversity and ecosystem function measurements [10], (2) unverified assumptions in common methods like root exclusion [48], and (3) the trade-off between methodological simplicity and comprehensive function assessment. Future methodological development should focus on harmonizing protocols across studies to enable meaningful cross-site comparisons and more robust tests of BEF theory in soil ecosystems.

As the field advances, integrating multiple measurement approaches within coherent theoretical frameworks will be essential for unraveling the complex relationships between soil biodiversity and ecosystem functioning—a crucial step for predicting ecosystem responses to global environmental change and informing sustainable soil management practices.

Biodiversity-ecosystem functioning (BEF) research has emerged as a critical field in ecology, investigating how the diversity of life influences ecosystem processes and services. Statistical and modeling tools provide the foundation for quantifying these complex relationships, testing hypotheses, and predicting outcomes under different environmental scenarios. In soil ecosystems, where immense biodiversity drives essential functions like nutrient cycling, carbon sequestration, and climate regulation, advanced analytical approaches are particularly vital for uncovering patterns and mechanisms that operate across multiple trophic levels [50] [51].

The historical context of BEF research reveals an ongoing evolution of analytical approaches. Beginning with early observations of natural communities, the field has progressed through theoretical models to sophisticated experiments that manipulate biodiversity as an independent variable [51]. Over the past 25 years, BEF research has demonstrated that the identity and combinations of species are powerful drivers of ecosystem processes, with positive biodiversity effects observed across spatial and temporal scales [51]. This research has moved conservation arguments beyond ethical motivations to include utilitarian perspectives grounded in empirical evidence of biodiversity's functional importance [51].

Foundational Statistical Frameworks for BEF Analysis

Core Statistical Approaches

BEF researchers employ a suite of statistical tools to analyze relationships between biodiversity measures and ecosystem functioning parameters. These tools range from basic descriptive statistics to advanced multivariate techniques, each serving specific analytical purposes in soil research.

Measures of central tendency (mean, median, mode) and dispersion (standard deviation, variance, range) provide fundamental characterization of datasets, forming the basis for more complex analyses [52]. In soil BEF studies, these basic statistics help researchers understand the distribution of soil organic carbon levels, microbial biomass measurements, and nutrient cycling rates across different management practices [53]. The t-test serves as a fundamental statistical tool for comparing means between two groups, such as assessing differences in soil carbon stocks between conventional and conservation agricultural practices [52].

Analysis of Variance (ANOVA) extends beyond t-tests by comparing means across multiple groups simultaneously, helping researchers avoid the increased error risk associated with multiple t-tests [52]. In soil BEF research, one-way ANOVA might examine the effect of a single factor (e.g., tillage practice) on a dependent variable (e.g., soil respiration), while two-way ANOVA investigates the effects of two independent variables and their interaction (e.g., tillage practice and cover cropping on microbial diversity) [52].

Regression analysis explores relationships between variables, predicting outcomes based on one or more predictors [52]. Linear regression models relationships between continuous variables, while multiple regression incorporates several independent variables to explain variations in the dependent variable. For example, researchers might use multiple regression to understand how soil texture, pH, and organic matter content collectively influence microbial community composition and function [52].

Table 1: Foundational Statistical Methods for BEF Research

Statistical Method	Primary Function	Soil BEF Application Example
T-test	Compare means between two groups	Compare mycorrhizal diversity between organic and conventional farms
ANOVA	Compare means across multiple groups	Test differences in soil carbon across cropping systems (monoculture, 2-species, 4-species mixtures)
Linear Regression	Model relationship between continuous variables	Predict soil organic carbon based on plant diversity metrics
Multiple Regression	Model multi-factor influences on an outcome	Predict nutrient cycling rates using microbial diversity, soil moisture, and temperature
Chi-Square Test	Analyze categorical data relationships	Examine association between soil type and presence/absence of key ecosystem engineers

Advanced Multivariate Techniques

As BEF research addresses increasingly complex questions, multivariate statistical techniques have become essential for analyzing the interconnected nature of soil ecological systems.

Factor analysis reduces complex datasets by identifying underlying variables that explain observed patterns [52]. This statistical tool proves valuable in soil BEF research for identifying latent variables that represent composite measures of soil health or ecosystem functioning. Principal component analysis (PCA) and related ordination techniques help visualize and interpret high-dimensional ecological data, revealing patterns that might not be apparent through univariate approaches [54].

The Generalized Linear Model (GLM) framework represents a more complex extension of traditional linear models, accommodating various data distributions and response types [55]. In soil BEF studies, GLMs can analyze count data (e.g., numbers of soil organisms), binary outcomes (e.g., presence/absence of nutrient limitation), and non-normal distributions common in ecological measurements [55].

Predictive Modeling Approaches for BEF Relationships

Classification Models

Classification models categorize data based on learned patterns from historical data, making them ideal for answering categorical questions in soil BEF research [55]. These models can address questions such as "Will this agricultural management practice improve soil health classification?" or "Does this soil community composition indicate degraded or healthy functioning?" The Random Forest algorithm stands as perhaps the most popular classification approach, capable of both classification and regression tasks with high accuracy on large volumes of data [55]. As an ensemble method, it combines multiple decision trees to create a robust predictive model that reduces the variance and bias associated with individual trees [55].

Clustering Models

Clustering models sort data into nested groups based on similar attributes, enabling pattern discovery without pre-defined categories [55]. In soil BEF research, clustering can identify distinct soil microbial communities across land-use gradients or group similar ecosystem functioning profiles across landscape positions. These approaches help researchers identify natural groupings in complex soil ecological data that might reflect underlying processes or states [55].

Forecast and Time Series Models

Forecast models estimate numeric values for new data based on historical patterns, applying to any research context with historical numerical data [55]. In soil BEF studies, forecast models might predict future soil carbon stocks under different climate scenarios or project nutrient cycling rates following biodiversity changes.

Time series models analyze data sequences using time as the input parameter, making them particularly valuable for understanding ecological dynamics [55]. These approaches can model seasonal patterns in soil processes, legacy effects of land-use history, and trajectory of ecosystem recovery following restoration interventions. Unlike simple averages, time series models better capture non-linear dynamics and exponential trends common in ecological systems [55].

Outlier Detection Models

Outlier models identify anomalous data entries within datasets, either individually or in conjunction with other values [55]. In soil BEF research, these approaches can detect unusual measurements that might indicate experimental error, unique ecological phenomena, or critical transition points in ecosystem state. For example, outlier detection might identify sites with unexpectedly high multifunctionality despite low biodiversity, suggesting the presence of unmeasured explanatory factors [55].

Table 2: Predictive Modeling Approaches for Soil BEF Research

Model Type	Primary Function	Strengths for Soil BEF Research
Classification	Categorize data into discrete classes	Identifying soil health status; Predicting management outcomes
Clustering	Group similar observations without pre-defined categories	Discovering natural patterns in microbial communities; Identifying functional groups
Forecast	Predict numerical values based on historical data	Projecting soil carbon sequestration; Estimating future nutrient availability
Time Series	Analyze temporal patterns and trends	Modeling seasonal dynamics; Detecting long-term trends in soil indicators
Outlier Detection	Identify anomalous observations	Flagging unusual measurements; Detecting ecosystem state transitions

Machine Learning Algorithms for BEF Prediction

Algorithm Selection and Comparison

Machine learning algorithms provide powerful approaches for analyzing complex BEF relationships in soil systems. Random Forest operates by creating an ensemble of decision trees, each depending on the values of a random vector sampled independently [55]. This algorithm offers numerous advantages for soil BEF research, including resistance to overfitting, ability to handle thousands of input variables without variable deletion, and capacity to estimate variable importance in classification [55].

The Gradient Boosted Model (GBM) also produces a prediction model composed of an ensemble of decision trees, but uses a different approach to combine these "weak learners" into a strong predictive model [55]. While both Random Forest and GBM employ multiple trees, they differ in their construction approach—Random Forest uses "bagging" (bootstrap aggregating) while GBM uses "boosting" that sequentially improves model performance [55].

Research comparing machine learning algorithms for ecological prediction demonstrates that optimal algorithm selection depends on the specific classification problem. A study comparing eight different machine learning algorithms for classifying beef quality attributes found that the best-performing algorithm varied depending on the specific classification task, suggesting that a "one size fits all" approach is not appropriate [54]. The highest-performing models for different classifications achieved prediction accuracies between 81.5–99%, indicating the potential of tailored machine learning approaches for complex biological classification problems [54].

Deep Learning Approaches

Deep learning represents a specialized subset of machine learning that is particularly valuable for analyzing unstructured data like images, audio, video, and text [55]. In soil BEF research, deep learning approaches can extract features from soil imagery, classify soil organisms using microscopic images, or analyze sensor data from environmental monitoring systems. AI-enhanced microscopy is emerging as a particularly promising approach for soil biodiversity assessment, offering scalable pathways to connect soil condition data with decision-making [56]. By making soil life visible and measurable, these technologies bridge the gap between expert knowledge and public understanding, enabling translation of scientific insights into actionable conservation strategies [56].

Experimental Design and Protocols for BEF Studies

Establishing BEF Experiments

Robust BEF research requires carefully designed experiments that can isolate biodiversity effects from confounding factors. The Jena Experiment, a prominent BEF study, exemplifies the long-term, multi-factorial approach needed to understand complex relationships in grassland ecosystems [51]. Such experiments typically manipulate species richness as an independent variable while controlling for environmental conditions, then measure multiple ecosystem functions simultaneously to assess multifunctionality.

Methodologies for global BEF syntheses involve systematic literature reviews following established protocols like the Population-Intervention-Control-Outcome (PICO) framework [53]. For example, a global dataset compiling results from 232 articles that experimentally compare effects of agricultural management practices on soil organic carbon accrual employed rigorous eligibility criteria and search terms to ensure data quality and comparability [53]. The resulting dataset included 570 experimental effects from 254 experiments across 38 countries, providing a robust foundation for analyzing patterns and context dependencies in BEF relationships [53].

Measuring Soil Biodiversity and Function

Comprehensive soil BEF assessment requires measuring biodiversity across multiple trophic levels and linking it to ecosystem processes. Advanced analytical techniques include rapid evaporative ionization mass spectrometry (REIMS), which enables rapid chemical fingerprinting of samples without extensive preparation [54]. This ambient mass spectrometry approach has applications in both human health and food science, with potential for adaptation to soil ecological research.

High-throughput sequencing of marker genes (e.g., 16S rRNA for bacteria, ITS for fungi) enables characterization of microbial communities, while metagenomics and metatranscriptomics provide insights into functional potential and activity. For soil fauna, extraction methods followed by morphological or molecular identification generate data on nematode, microarthropod, and enchytraeid communities. Integrating these diverse measurements requires sophisticated statistical approaches that can handle different data types and scales.

Table 3: Essential Methodologies for Soil BEF Experiments

Method Category	Specific Techniques	Measured Parameters
Biodiversity Assessment	High-throughput sequencing, AI-enhanced microscopy, Morphological identification	Species richness, Community composition, Phylogenetic diversity, Food web structure
Ecosystem Function Measurement	Soil respiration chambers, Nutrient leaching collection, Enzyme assays, Stable isotope tracing	Carbon cycling, Nutrient mineralization, Decomposition rates, Multifunctionality
Abiotic Factor Quantification	Soil coring and analysis, Climate station monitoring, Moisture sensors	Soil texture, pH, Organic matter, Temperature, Moisture, Bulk density
Statistical Analysis	Multivariate ordination, Structural equation modeling, Network analysis, Machine learning	Biodiversity effects, Context dependencies, Interaction pathways, Predictive relationships

Conceptual Framework and Analytical Workflow

The complexity of soil BEF relationships necessitates conceptual frameworks that accommodate multitrophic interactions, environmental context dependencies, and spatial-temporal dynamics. A multitrophic perspective recognizes that BEF relationships emerge from interactions across trophic levels, with soil network complexity enhancing links between biodiversity and ecosystem multifunctionality [50] [51].

The following diagram illustrates a generalized analytical workflow for soil BEF studies, from experimental design through interpretation:

Research Reagent Solutions and Essential Materials

Conducting robust soil BEF research requires specialized reagents and materials for accurate measurement of biodiversity and ecosystem functions. The following table details key solutions and their applications in soil BEF studies:

Table 4: Essential Research Reagents and Materials for Soil BEF Studies

Reagent/Material	Primary Function	Application in Soil BEF Research
DNA/RNA Extraction Kits	Nucleic acid isolation from soil matrices	Molecular characterization of microbial, fungal, and microfaunal communities
Stable Isotope Tracers (e.g., ¹⁵N, ¹³C)	Tracking element movement through food webs	Quantifying nutrient cycling rates, energy flows, and trophic interactions
Enzyme Assay Substrates	Measuring extracellular enzyme activities	Assessing microbial functional potential for carbon, nitrogen, and phosphorus cycling
PCR Reagents and Primers	Target gene amplification for diversity assessment	Taxonomic and functional gene characterization through amplicon sequencing
Soil Respiration Chambers	Contained measurement of CO₂ flux	Quantifying soil heterotrophic respiration and carbon mineralization rates
Lysimeter Systems	Collection of soil pore water	Measuring nutrient leaching losses and water quality impacts
Mass Spectrometry Reagents	Chemical fingerprinting of samples	Rapid characterization of soil metabolic profiles using REIMS approaches
AI-Enhanced Microscopy Platforms	Automated image analysis of soil organisms	High-throughput characterization of soil fauna and root-microbe interactions

Context Dependencies and Scaling in BEF Relationships

A critical insight from BEF research is that relationships between biodiversity and ecosystem functioning are strongly context-dependent, varying with environmental conditions, disturbance regimes, and spatial-temporal scales [51]. Soil BEF relationships have been shown to depend on climatic conditions, local site characteristics, and management practices, which interact with biodiversity in complex ways [51].

The multitrophic perspective has been particularly valuable for understanding context dependencies in soil systems, where interactions across trophic levels strongly influence ecosystem processes [51]. Research has demonstrated that soil phylotypes with larger sizes or at higher trophic levels (e.g., invertebrates or protist predators) may exhibit weaker or no BEF relationships compared to those with smaller sizes or at lower trophic levels (e.g., archaea, bacteria, fungi, and protist phototrophs) [50].

The role of soil network complexity, reflected by co-occurrence patterns among multitrophic-level organisms, appears crucial in enhancing the link between soil biodiversity and ecosystem functions [50]. This insight represents a significant advance in forecasting the impacts of belowground multitrophic organisms on ecosystem functions in agricultural systems, suggesting that soil multitrophic network complexity should be considered a key factor in enhancing ecosystem productivity and sustainability under land-use intensification [50].

The following conceptual diagram illustrates the complex, context-dependent nature of soil BEF relationships:

Future Directions in BEF Statistical Modeling

Soil BEF research continues to evolve with emerging analytical approaches and priority research areas. Six key frontiers have been identified for future BEF research: (1) non-random biodiversity change across trophic levels; (2) predicting the strength of BEF relationships across environmental contexts; (3) spatial scaling of BEF relationships; (4) eco-evolutionary implications of multitrophic BEF; (5) FAIR data (findable, accessible, interoperable, reusable) and reproducible processing; and (6) operationalizing BEF insights for ecosystem management, society, and decision making [51].

AI-enhanced microscopy represents a promising technological development for soil biodiversity assessment, providing scalable pathways to connect soil condition data with decision-makers and inform evidence-based governance [56]. By making soil life visible and measurable, these technologies bridge the gap between expert knowledge and public understanding, enabling translation of scientific insights into actionable conservation and management strategies [56].

The integration of FAIR data principles and reproducible processing will be key to advancing soil BEF research by making it more integrative and collaborative [51]. As datasets grow in size and complexity, standardized data management and analytical workflows will ensure that research findings are transparent, reproducible, and accessible to diverse stakeholders including researchers, policymakers, and land managers.

Understanding the relationship between plant diversity and soil fungal communities is fundamental to testing biodiversity-ecosystem functioning (BEF) relationships in soils research. Grassland ecosystems, which provide critical services including carbon sequestration, nutrient cycling, and soil fertility maintenance, serve as ideal models for investigating these interactions [57]. Soil fungi, encompassing mutualistic mycorrhizal fungi, decomposers, and pathogens, constitute key biological interfaces driving ecosystem multifunctionality—the simultaneous performance of multiple ecosystem processes [58] [59]. This case study synthesizes experimental data from recent research to compare how different experimental approaches elucidate the mechanisms underlying plant diversity-soil fungal interactions and their consequences for ecosystem functioning.

Comparative Experimental Data on Plant-Fungal Interactions

Table 1: Key quantitative findings from plant diversity-fungal community studies

Study System & Approach	Plant Diversity Effect	Fungal Response	Ecosystem Multifunctionality Impact
Experimental Plant Communities (Manipulation) [58]	Low (1-2 species) vs. High (up to 8 species)	21% reduction in fungal diversity with fungicide decreased multifunctionality only in low-diversity communities	Fungal diversity buffers ecosystem multifunctionality against plant diversity loss
Natural Grassland Gradient (Observational) [57]	High-diversity plots (>20 species/0.1m²) vs. Low-diversity (<10 species/0.1m²)	Higher fungal diversity, different community composition, and increased AMF biomass in high-diversity plots	Soil chemistry drove bacterial communities; plant community type primarily shaped fungal communities
Root Economics Space (Trait-Based) [60]	Variation along collaboration ("outsourcing" vs "DIY") and conservation ("fast-slow") gradients	"Slow" root traits favored saprotrophs; "fast" and "DIY" traits favored pathogenic fungi	Root economics space and plant richness jointly determined fungal guild composition
Vegetation Restoration Chronosequence (Observational) [61]	Land use types (Artificial forest, Natural shrub, Grassland, Cropland)	Fungal diversity followed: Natural shrub > Artificial forest > Natural grassland > Artificial grassland > Cropland	Soil organic carbon and plant richness were primary drivers of fungal diversity

Table 2: Fungal guild-specific responses to plant diversity

Fungal Guild	Response to Plant Diversity	Driving Factors	Functional Consequences
Arbuscular Mycorrhizal Fungi (AMF)	Positive correlation with plant diversity [57] [60]	Root collaboration gradient ("outsourcing" strategy with thicker roots) [60]	Enhanced plant nutrition, pathogen protection [60]
Saprotrophic Fungi	Increased diversity in species-rich communities with "slow" root traits [60]	Litter quality and quantity [60]	Organic matter decomposition, carbon cycling [62]
Plant Pathogenic Fungi	Highest in communities with "fast" and "DIY" root traits [60]	Reduced host density in diverse communities [60]	Plant fitness regulation, diversity maintenance [62]
Ectomycorrhizal Fungi	Affected by host identity and abundance [62]	Specific tree hosts [62]	Host plant nutrient acquisition [62]

Methodological Approaches: Experimental Protocols

Plant Diversity Manipulation Experiment

The most direct approach for establishing causality involves experimental manipulation of plant diversity combined with fungal suppression:

Plant Community Assembly: Researchers established 190 experimental plant communities with controlled species richness gradients (1, 2, 4, 6, and 8 species) in greenhouse conditions [58].

Fungal Diversity Manipulation: Soil fungal diversity was experimentally reduced through fungicide addition (applied prior to planting), which successfully reduced fungal diversity by 21% without affecting fungal abundance, creating a crossed design of plant diversity × fungal diversity [58].

Ecosystem Function Assessment: Ten ecosystem functions were quantified, including primary productivity, floral abundance, crown interception of light/water/wind, water conservation, microbial biomass, litter decomposition, soil enzyme activity, and carbon, nitrogen, and phosphorus stocks [58].

Statistical Analysis: Ecosystem multifunctionality was calculated using both average value and multiple thresholds methods. Complementarity and selection effects were quantified using plant aboveground biomass data [58].

Natural Gradient Observational Approach

For established ecosystems, researchers employ comparative sampling along existing plant diversity gradients:

Paired-Plot Design: In Central European grasslands, researchers identified and sampled 12 plots (0.1m²) distributed in six pairs, each containing one low-diversity plot (<10 plant species) near forest edges or solitary trees, and one high-diversity plot (>20 species) in treeless patches [57].

Comprehensive Sampling: Each plot underwent complete botanical survey, shoot biomass collection, soil core extraction (10cm depth), and root system retrieval. Large roots (>2mm diameter) and debris were removed before soil sieving (5mm mesh) and homogenization [57].

Multi-Omic Analysis: Soil samples were processed for DNA extraction, followed by high-throughput sequencing of fungal ITS regions for community composition, phospholipid fatty acid (PLFA) analysis for microbial biomass, and standard soil chemistry analysis (pH, organic carbon, nitrogen, etc.) [57].

Data Integration: Multivariate statistics (PERMANOVA, RDA) linked microbial community data to environmental variables, while structural equation modeling tested hypothesized pathways between plant diversity, soil properties, and microbial communities [61] [57].

Root Trait-Based Approach

Linking plant functional traits to fungal communities provides mechanistic understanding:

Root Economics Space Characterization: Researchers collected and analyzed root traits along two orthogonal axes: the conservation gradient (root tissue density, nitrogen content - "fast-slow" spectrum) and the collaboration gradient (specific root length, root diameter - "outsourcing" vs "DIY" spectrum) [60].

Fungal Guild Partitioning: Soil samples were subjected to DNA extraction and high-throughput sequencing with subsequent functional classification of fungi into guilds (saprotrophic, pathogenic, mycorrhizal) using curated databases [60].

Community-Level Trait Aggregation: Plant community-weighted mean traits were calculated based on species abundances and their trait values [60].

Conceptual Workflow and Signaling Pathways

The following diagram illustrates the conceptual framework and experimental workflow for investigating plant diversity-soil fungal interactions in grassland ecosystems:

Conceptual Workflow for Studying Plant-Fungal Interactions This diagram outlines the major methodological pathways for investigating biodiversity-ecosystem functioning relationships between plants and soil fungi, highlighting complementary experimental and observational approaches.

The following diagram illustrates the mechanistic pathways through which plant diversity influences ecosystem multifunctionality via soil fungal communities:

Mechanistic Pathways from Plant Diversity to Ecosystem Multifunctionality This diagram visualizes the key mechanisms and pathways through which plant diversity influences ecosystem multifunctionality via alterations in soil fungal communities, highlighting both direct and indirect effects.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential research reagents and methodologies for plant-fungal interaction studies

Reagent/Method	Application	Specific Function	Example Implementation
ITS2 rRNA Sequencing	Fungal community profiling	Amplification and sequencing of ITS2 region for fungal identification and diversity assessment [61] [63]	Illumina HiSeq platform with ITS2-specific primers [63]
Fungicide Treatments	Experimental reduction of fungal diversity	Selective suppression of fungal communities to test causality in plant-fungal relationships [58]	Soil application prior to planting; reduces diversity by 21% without affecting abundance [58]
PLFA/NLFA Analysis	Microbial biomass quantification	Measurement of phospholipid and neutral lipid fatty acids to determine fungal biomass and activity [57]	Soil sample extraction and GC-MS analysis for AMF-specific lipid biomarkers [57]
Root Economics Space Traits	Plant functional strategy characterization	Quantification of root traits along collaboration and conservation gradients [60]	Measurement of specific root length, root diameter, tissue density, nitrogen content [60]
Structural Equation Modeling	Causal pathway analysis	Statistical modeling to test direct and indirect effects in complex plant-soil-fungal systems [61] [63]	Analysis of climate → soil properties → fungal diversity → multifunctionality pathways [63]
Co-occurrence Network Analysis	Microbial interaction mapping	Construction of correlation networks to visualize fungal community structure and identify keystone taxa [63] [64]	Spearman correlation networks of abundant taxa; module identification [64]

This comparison of experimental approaches demonstrates that plant diversity consistently shapes soil fungal communities across grassland ecosystems, with significant consequences for ecosystem multifunctionality. Manipulative experiments provide the strongest evidence for causality, revealing that fungal diversity buffers ecosystem functioning against plant diversity loss [58]. Natural gradient studies offer ecological realism and identify soil properties and plant community composition as major drivers of fungal communities [61] [57]. The emerging root economics framework provides mechanistic understanding by linking specific root trait strategies to fungal guild composition [60]. Together, these approaches advance our understanding of BEF relationships in soils and highlight the critical importance of conserving both plant and fungal diversity for maintaining ecosystem services. Future research should integrate these approaches across wider environmental gradients and temporal scales to predict ecosystem responses to global change.

Navigating Complexity: Context-Dependency and Troubleshooting in Soil BEF Research

The relationship between biodiversity and ecosystem functioning (BEF) is a cornerstone of ecological research. While a positive correlation is well-established, this relationship is not universal. In soil ecosystems, it is profoundly mediated by abiotic factors, creating a significant "context-dependency challenge." Understanding how soil type, fertility, and climate modulate BEF relationships is crucial for predicting ecosystem responses to global change and for developing effective soil management and conservation strategies. This guide synthesizes current experimental data to objectively compare how these contextual factors shape BEF relationships in soil environments.

How Climate Sets the Boundaries for BEF Relationships

Climate, particularly temperature and precipitation, operates as a master regulator that sets the broad thermal and moisture boundaries within which soil organisms function. Recent global analyses reveal that its influence is not linear but can involve critical thresholds.

Temperature Thresholds and Global Projections

A 2025 global analysis of belowground ecosystem multifunctionality (BEMF) identified a mean annual temperature (MAT) threshold of 16.4°C [65]. This threshold separates two distinct regimes:

• Below 16.4°C (High BEMF regime): BEMF is generally higher but decreases rapidly with rising temperature. Soil nutrient pools are the primary drivers of multifunctionality in these cooler regions [65].
• Above 16.4°C (Low BEMF regime): BEMF is lower and less sensitive to temperature changes. In these warmer regions, precipitation and plant species richness become the dominant positive drivers of multifunctionality [65].

This threshold effect has significant consequences. Climate models project that ongoing climate change could lead to a 20.8% loss of global BEMF by 2100 under a high-emission scenario (SSP585), with temperate and continental biomes expected to suffer the most significant losses [65].

Precipitation and Aridity

The influence of precipitation is often intertwined with temperature. In drylands, which cover over 41% of the Earth's surface, low water availability strongly limits ecosystem functioning and soil fertility [24]. Research in West African semi-arid savannas demonstrates that the positive relationship between tree species richness and aboveground carbon (AGC) is stronger in areas with higher mean annual precipitation (MAP) [66]. This suggests that water availability can modulate the strength of biodiversity effects on key ecosystem functions like carbon storage.

Table 1: Climate-Mediated Shifts in Belowground Ecosystem Multifunctionality (BEMF)

Climatic Factor	Observed Effect on BEF Relationship	Key Data/Threshold	Underlying Mechanism
Mean Annual Temperature (MAT)	Abrupt shift in global BEMF patterns [65]	Threshold at 16.4°C [65]	Shift from temperature/pH-dominated controls (cooler regions) to precipitation/biodiversity-dominated controls (warmer regions) [65]
Mean Annual Precipitation (MAP)	Modulates strength of diversity-carbon relationship [66]	Positive relationship between species richness and aboveground carbon strengthens with higher MAP [66]	Enhanced water availability supports greater stem density and growth, enabling fuller expression of biodiversity effects [66]
Aridity	Limits soil fertility and functioning [24]	Low water availability restricts nutrient cycling and organism activity [24]	Scarcity of water directly stresses biota and reduces the input and decomposition of organic matter [24]

The Mediating Role of Soil Type and Fertility

At a more local scale, soil properties such as texture, pH, and nutrient content act as filters on soil communities and modify the relationship between biodiversity and function.

Soil Fertility and Functional Redundancy

Soil fertility can determine the sensitivity of an ecosystem to biodiversity loss. A 2025 study across five European grasslands found that the initial soil organic matter content (a key indicator of fertility) influenced how ecosystem functions responded to changes in soil biodiversity [3].

• Fertile Soils: Soils rich in organic matter and nutrients often support more diverse communities with higher functional redundancy. This means multiple taxa can perform the same function, making the ecosystem more resilient to biodiversity decline, particularly for general processes like nutrient mineralization [3].
• Poor Soils: In contrast, nutrient-poor soils are more dependent on specific soil functional groups. Their functioning is more sensitive to biodiversity loss, as the loss of key taxa is less likely to be compensated for by others [3]. For instance, the functioning of poorly developed soils in drylands relies heavily on the daily contribution of diverse microbiota to decompose organic matter and introduce nutrients [15] [6].

Soil Texture and Chemistry

The physical and chemical composition of soil directly influences BEF relationships:

• Texture: The "inverse-texture hypothesis" suggests that in arid and semi-arid regions, coarser-textured (sandy) soils may support higher primary productivity than fine-textured (clay) soils, a relationship that reverses in humid regions [66]. This occurs because sandy soils in dry areas allow for better water infiltration and reduced evaporation, making limited water more available to plants and soil organisms.
• pH and Nutrient Pools: In cool regions (MAT ≤ 16.4°C), soil pH generates strong negative effects on BEMF alongside temperature [65]. Furthermore, the vast nutrient pools stored in cold, often acidic, polar and continental soils are primary drivers of their high multifunctionality, despite having lower productivity and decomposition rates than tropical biomes [65].

Table 2: Experimental Evidence for Soil Type and Fertility Mediating BEF Relationships

Experimental Context	Key Finding on BEF Relationship	Methodology Summary	Citation
Five European Grasslands (varying soil fertility)	Soil organic matter content influenced functional redundancy; effects of biodiversity decline were context-dependent, not systematically stronger in poorer soils [3].	Soil Biodiversity Filtering: Sterilized soils from five sites were inoculated with soil organisms extracted from donor grasslands and filtered by size to create a biodiversity gradient. Plant communities were established, and functions like productivity and nutrient cycling were measured [3].	[3]
Global Drylands & Poorly Developed Soils	Soil biodiversity is particularly critical for supporting function in drylands and nutrient-poor soils [15] [6].	Field Sampling & Correlation Analysis: Global-scale field surveys measuring soil biodiversity (e.g., bacteria, fungi, protists, invertebrates) and multiple ecosystem functions (e.g., nutrient cycling, carbon storage, decomposition) [15] [6].	[15] [6]
Atlantic Forest Landscape (natural vs. degraded)	Land-use simplification for agriculture (e.g., pastures) reduces soil biodiversity, which directly and indirectly diminishes ecosystem multifunctionality [18].	Land-Use Comparison: Linear mixed-effects models applied to field data from natural forests, pastures, and deforested areas, linking soil biodiversity to ecosystem functions across seasons [18].	[18]

Conceptual Framework and Pathways

The interaction between climate, soil, and biodiversity in driving ecosystem function involves complex, interconnected pathways. The following diagram synthesizes these relationships into a conceptual model.

Diagram 1: Context-Dependency of Soil BEF Relationships. This model shows how climate and soil properties directly influence ecosystem function, but also operate indirectly by mediating soil biodiversity and stand structure. Threshold effects (e.g., temperature) and soil properties like fertility and texture are key sources of context-dependency.

The Scientist's Toolkit: Key Research Reagents and Methods

To investigate these complex BEF relationships, researchers employ a suite of field, laboratory, and molecular tools. The table below details essential reagents, materials, and methods used in the cited studies.

Table 3: Research Reagent Solutions for Soil BEF Studies

Reagent / Material / Method	Primary Function in BEF Research	Experimental Application Example
Size-Based Fractionation	To experimentally create gradients of soil biodiversity by filtering out organism groups based on body size [3].	Separating soil organisms into size classes (e.g., >40μm, >20μm, >5μm) to prepare inocula for microcosm experiments, allowing testing of the role of different functional groups (microfauna, mesofauna, microbes) [3].
"-Omics" Techniques (e.g., metagenomics)	To characterize the taxonomic and functional diversity of soil microbial communities (bacteria, archaea, fungi, protists, viruses) [15] [6].	Profiling soil microbiomes from global drylands to link microbial diversity with ecosystem multifunctionality [15] [6].
Soil Sterilization	To create a sterile background substrate for inoculation experiments, isolating the effect of the added soil community.	Autoclaving or gamma-irradiating soil to eliminate native biota before inoculation with filtered biodiversity treatments in microcosms [3].
Standardized Soil Assays	To quantify key ecosystem functions related to nutrient cycling and microbial activity.	Measuring microbial biomass carbon/nitrogen (chloroform fumigation), gross nitrogen mineralization rates, soil respiration, and available phosphorus [65] [3].
Structural Equation Modeling (SEM)	A statistical method to test and quantify complex networks of direct and indirect effects between climate, soil, biodiversity, and ecosystem functions [65] [66].	Used in global and landscape studies to disentangle how climate and soil properties directly influence BEMF versus indirectly influencing it via shifts in biodiversity [65] [66].

The relationship between soil biodiversity and ecosystem functioning is not a simple one-size-fits-all paradigm. It is a context-dependent relationship powerfully mediated by climate—which can create abrupt thresholds like the 16.4°C MAT tipping point—and by local soil properties like fertility and texture, which determine functional redundancy and resource availability. Overcoming the context-dependency challenge requires research that explicitly integrates these factors across scales, from global models that predict broad thresholds to targeted experiments that unravel the mechanisms in specific soils. This integrated understanding is paramount for forecasting the impacts of global change and for managing soil ecosystems to maintain their vital functions for future generations.

Soil biodiversity represents one of the largest reservoirs of biological diversity on Earth, hosting organisms ranging from microorganisms (e.g., bacteria, fungi, protists) to microfauna (e.g., nematodes), mesofauna (e.g., microarthropods), and macrofauna (e.g., earthworms, beetles) [3]. This biodiversity plays fundamental roles in driving essential ecosystem functions and services, including climate regulation, nutrient cycling, food production, and even human health through the One Health framework [6]. However, our understanding of soil biodiversity-ecosystem functioning (BEF) relationships remains critically limited by systematic biases in research focus and geographical coverage. Recent syntheses reveal that soil macroecological studies suffer from significant blind spots that constrain our ability to draw general conclusions about BEF relationships across different environmental conditions and taxonomic groups [10]. These biases impede the development of effective conservation policies and sustainable land management practices, particularly in the context of global environmental change. This review synthesizes evidence of these biases, quantifies their extent, and provides methodological frameworks for overcoming these limitations in future soil health research.

Quantifying the Geographic and Taxonomic Biases in Current Research

Documented Geographic Biases in Soil Research

Table 1: Documented Geographic Biases in Soil Biodiversity and Function Research

Region/System Type	Representation in Studies	Specific Gaps Identified	Consequences
Tropical Zones	Highly unrepresentative field observations [67]	Skewed toward high-precipitation regions with allophanic clay mineralogy [67]	Inaccurate extrapolation of land-cover change effects on soil C stocks [67]
Temperate Systems	Overrepresented in global datasets [10]	-	Limited understanding of global patterns; temperate-centric models [10]
Global Soil Sampling	Only 17,186 sampling sites documented [10]	Far fewer than aboveground databases (e.g., PREDICTS: ~29,000 sites) [10]	Reduced power for global analyses and macroecological relationships [10]
Seven High-Income Countries	Prevail in evidence base [68]	USA, China, and Brazil particularly dominant [68]	Lack of knowledge for regions with different management practices and environmental conditions [68]

The geographic distribution of soil biodiversity research exhibits pronounced systematic biases. A meta-analysis of studies quantifying changes in soil carbon stocks with land use in the tropics found that field observations are "highly unrepresentative of most tropical landscapes" [67]. This bias is particularly concerning given that tropical latitudes account for the bulk of current CO₂ emissions from land-cover change. The distribution of field observations was found to be highly skewed toward high-precipitation regions with allophanic clay mineralogy, while historically, land-conversion activities in the tropics have focused on high-activity clay soils in lower precipitation regions [67]. This mismatch between research focus and actual land-use patterns strongly cautions against extrapolating average values of land-cover change effects on soil carbon stocks to regions that differ in biophysical conditions.

At a global scale, an analysis of 17,186 sampling sites from macroecological studies on soil biodiversity and ecosystem functions revealed that data are not evenly distributed worldwide [10]. There is a clear concentration of studies in temperate systems, while tropical, polar, and arid ecosystems remain significantly underrepresented. This bias is further exacerbated in secondary research, with a recent systematic map of 200 meta-analyses (gathering over 9,000 primary studies) confirming that evidence predominantly stems from a limited set of high-income countries, notably the United States, China, and Brazil [68]. This unequal representation creates significant blind spots in our understanding of soil biodiversity and functioning across diverse global ecosystems.

Taxonomic and Methodological Biases in Soil Studies

Table 2: Taxonomic and Methodological Biases in Soil Biodiversity Research

Category	Current Focus	Underrepresented Elements	Impact on BEF Understanding
Taxonomic Groups	Bacteria, fungi, and Formicoidea (48.8% of all records) [10]	Rotifera, Collembola, Acari, annelids, vertebrates, plants [68] [10]	Incomplete picture of soil food webs and functional diversity
Ecosystem Functions	Soil respiration (78.8% of all function records) [10]	Nutrient cycling, secondary productivity, multiple simultaneous functions [10]	Limited understanding of multifunctionality and interacting processes
Biodiversity Metrics	Abundance data [68]	Functional and phylogenetic diversity metrics [68]	Oversimplification of biodiversity-ecosystem function relationships
Management Focus	Individual practices (fertilizer, phytosanitary, diversification) [68]	Farm and landscape levels, practice combinations [68]	Limited relevance to real-world farming contexts with multiple interventions
Temporal Data	Single sampling events dominate [10]	Repeated measurements over multiple years [10]	Inability to assess trends and vulnerability to global change

Taxonomic biases in soil biodiversity research are equally pronounced. Analysis of global soil macroecological studies reveals that bacteria, fungi, and Formicoidea (ants) together constitute 48.8% of all soil biodiversity records, while other ecologically important groups such as Rotifera, Collembola, and Acari have substantially lower representation [10]. This bias extends to secondary research, where arthropods and microorganisms are most frequently studied, while annelids, vertebrates, and plants are less represented [68]. The taxonomic bias is particularly problematic given that different organism groups contribute uniquely to ecosystem functioning, with larger organisms (e.g., earthworms, nematodes) being crucial for maintaining soil health through organic matter decomposition, nutrient cycling, and soil aeration, while smaller organisms (e.g., bacteria, fungi) support a large number of functions at different activity levels [6].

Methodologically, soil biodiversity research shows a strong preference for simple abundance data, with substantial gaps in studies on functional and phylogenetic diversity [68]. This limits our understanding of how biodiversity influences ecosystem functioning beyond simple species numbers. Furthermore, research predominantly focuses on individual agricultural practices rather than their combinations, overshadowing investigations at the farm and landscape levels that would better reflect real-world farming contexts [68]. Perhaps most critically, there is an "almost complete absence of temporally explicit data" in soil macroecological studies [10], with most studies based on single sampling events rather than repeated measurements over time, severely limiting our ability to assess trends and vulnerability to global change.

Critical Consequences: How Biases Impair BEF Relationship Understanding

The geographic and taxonomic biases in soil health research have profound implications for our understanding of biodiversity-ecosystem functioning relationships. The most significant consequence is the limited overlap between biodiversity and ecosystem function data—only 0.3% of all sampling sites have both information about biodiversity and function, and even these few sites have non-systematic coverage of taxa and functions [10]. This data gap severely constrains macroecological analyses of soil BEF relationships at global scales.

The interaction between biophysical factors and land-use impacts further complicates extrapolation from biased datasets. Research has demonstrated that the direction and magnitude of changes in soil carbon pools with land-use conversion depend strongly on biophysical factors such as mean annual precipitation and dominant soil clay mineralogy [67]. For example, conversion of forest to pasture can increase soil carbon stocks on soils with low-activity clay in moderate precipitation regions (1,501-2,500 mm annually) but decrease stocks on allophanic soils in high precipitation regions (>3,501 mm annually) [67]. When research focus is biased toward specific environmental conditions, predictions for underrepresented regions become highly uncertain.

The focus on individual practices rather than their combinations limits practical applications, as real-world agricultural management typically involves implementing multiple practices simultaneously [68]. This oversight reduces the relevance of research findings for farmers and policymakers attempting to implement sustainable land management strategies. Furthermore, the concentration of studies in temperate systems and a few high-income countries means that soil biodiversity conservation and management strategies may be ill-suited for regions with different environmental conditions and socioeconomic contexts.

Methodological Frameworks for Mitigating Research Biases

Sampling and Analysis Approaches

Representativeness Heuristic for Spatial Bias Mitigation Digital Soil Mapping (DSM) approaches offer promising frameworks for addressing spatial bias in existing soil samples. One proposed method involves quantifying sample representativeness as the "goodness-of-coverage" of sample locations over the environmental covariate space, then improving this representativeness through weighting samples toward maximal goodness-of-coverage [69]. This approach mitigates spatial bias by assigning optimal weights to existing samples based on their coverage of environmental gradients, rather than discarding valuable data. Case studies have demonstrated that this weighting approach can improve DSM accuracy, particularly for mapping soil organic matter [69].

Standardized Soil Health Assessment Frameworks Comparative analyses of soil health assessment frameworks provide guidance for standardizing methodologies across studies. Research comparing the Soil Management Assessment Framework (SMAF) with the Haney Soil Health Test (HSHT) found that SMAF provides a more comprehensive assessment through evaluation of physical, chemical, and biological indicators, while HSHT primarily reflects biological health through the Solvita CO₂-C burst (r = 0.82) [70]. The SMAF methodology includes ten indicators: bulk density (Bd), water-stable aggregates (WSAs), β-glucosidase activity (BG), microbial biomass carbon (MBC), potential mineralizable nitrogen (PMN), soil organic carbon (SOC), pH, electrical conductivity (EC), and extractable P and K [70]. This comprehensive approach facilitates more comparable results across different regions and soil types.

Experimental Protocols for BEF Relationship Elucidation Well-designed experiments can specifically address biodiversity-ecosystem function relationships while controlling for biases. One experimental approach involves establishing gradients of soil biodiversity through size-based filtering of soil organisms to create different biodiversity treatments (e.g., High, Medium, Low, and Minimum diversity levels) [3]. These treatments are then used to inoculate sterilized soil substrates from different grassland systems, allowing researchers to test BEF relationships across varying environmental contexts. Measurements typically include plant biomass production, soil microbial respiration, microbial biomass carbon and nitrogen, nutrient cycling rates, and molecular analyses of microbial communities [3].

Diagram Title: Research Blind Spot Mitigation Workflow

The Scientist's Toolkit: Essential Methodologies and Reagents

Table 3: Research Reagent Solutions for Comprehensive Soil Health Assessment

Category	Tool/Method	Specific Function	Considerations for Bias Reduction
Molecular Analysis	High-throughput sequencing	Characterizes microbial communities (bacteria, fungi, protists)	Enables inclusion of underrepresented microbial taxa [3]
Physical Assessment	Bulk density measurement	Quantifies soil compaction and pore space	Standardized in SMAF framework for cross-study comparison [70]
Chemical Analysis	Water-stable aggregates (WSA)	Measures soil structural stability	Physical indicator in comprehensive frameworks [70]
Biological Activity	Solvita CO₂-C burst test	Assesses microbial respiration	Primary driver of HSHT score; biological focus [70]
Enzyme Assays	β-glucosidase (BG) activity	Measures carbon cycling potential	Included in SMAF as key biological indicator [70]
Nutrient Cycling	Potential mineralizable N (PMN)	Estimates nitrogen availability	Biological indicator in SMAF [70]
Microbial Biomass	Microbial biomass C (MBC)	Quantifies living microbial component	Biological indicator in SMAF [70]
Experimental Systems	Size-based filtering	Creates biodiversity gradients	Allows testing BEF relationships across systems [3]

The identified geographic and taxonomic biases in soil health research represent significant impediments to understanding biodiversity-ecosystem functioning relationships and developing effective conservation strategies. The concentration of studies in temperate systems, focus on limited taxonomic groups, and absence of temporal data collectively constrain our ability to predict how soil ecosystems will respond to global environmental changes. Addressing these blind spots requires concerted efforts to expand research into underrepresented regions and ecosystems, incorporate a wider range of taxonomic groups, develop standardized assessment methodologies that facilitate cross-study comparisons, and implement long-term monitoring programs. Only through such comprehensive approaches can we develop a truly global understanding of soil biodiversity-ecosystem functioning relationships and formulate effective strategies for conserving these vital ecosystems under changing environmental conditions.

Understanding the relationship between biodiversity and ecosystem functioning represents a central goal in ecology. In soil ecosystems, this relationship is primarily mediated through complex trophic interactions and food webs. While it is established that soil biodiversity contributes significantly to ecosystem functions such as decomposition, nutrient cycling, and carbon sequestration, the specific mechanisms through which direct and indirect effects operate within soil food webs remain poorly characterized. This review synthesizes current experimental approaches and findings regarding the disentanglement of these effects, focusing specifically on methodological frameworks and their applications in soil ecosystems. We examine how trophic interactions modulate biodiversity-ecosystem functioning relationships and provide a comparative analysis of key experimental protocols used to quantify these relationships.

Theoretical Framework: Trophic Interaction Modifications (TIMs) in Soil Ecosystems

Defining Direct and Indirect Trophic Effects

In soil food webs, species are connected through consumer-resource relationships that often involve third-party species that modify the strength and outcome of these interactions. These trophic interaction modifications (TIMs) represent a fundamental class of indirect effects where one species alters how another species interacts with its resources or predators [71]. The TIM framework provides a powerful approach for understanding soil ecosystem dynamics because it explicitly recognizes the multi-species nature of these interactions rather than attempting to reduce them to pairwise relationships [71].

Soil ecosystems present particular challenges for studying TIMs due to the immense diversity of organisms, the physical complexity of the soil matrix, and the difficulty in observing interactions directly. A single gram of soil can contain millions of bacterial individuals representing thousands of species, alongside diverse populations of fungi, protozoa, nematodes, and microarthropods [6]. Within this complex environment, TIMs emerge through various mechanisms including behavioral changes in response to predators (e.g., nematodes altering feeding behavior in presence of microarthropod predators), ecosystem engineering (e.g., earthworms creating burrows that affect microbial communities), and chemical mediation (e.g., microbial production of compounds that influence predator-prey interactions) [5] [71].

Conceptual Model of Interaction Modifications

The following diagram illustrates how trophic interaction modifications operate within a three-component soil food web, where one species modifies the interaction between two others:

Diagram 1: Trophic Interaction Modification Framework. This illustrates how a modifier species alters the base trophic interaction between consumer and resource species, resulting in a modified interaction that affects feeding rates or behavior.

In soil systems, these modification networks become exponentially more complex. For instance, earthworms (modifier) can alter the interaction between bacteria (resource) and bacterivorous nematodes (consumer) through their burrowing activities that redistribute bacterial communities and change their accessibility to nematodes [5]. Simultaneously, these nematodes may themselves modify the interaction between fungi and their fungal-feeding predators through chemical signaling [71]. This multi-layered network of direct and indirect effects creates the complex dynamics observed in soil biodiversity-ecosystem functioning relationships.

Quantitative Approaches for Measuring Interaction Strengths

Metrics for Quantifying Trophic Interaction Modifications

Researchers have developed multiple metrics to quantify the strength of TIMs in ecological systems, each incorporating different levels of contextual information about the system. These metrics can be applied to soil food webs to measure how species losses, environmental changes, or management practices alter interaction strengths and subsequent ecosystem functions [71].

Table 1: Metrics for Quantifying Trophic Interaction Modification Strength

Metric Type	Definition	Application in Soil Systems	Key References
Relative Modification Strength	Proportional change in consumption rate	Measuring how predator presence alters nematode grazing on bacteria	[71]
Density-Mediated TIM	Effect driven by changes in modifier density	Earthworm density effects on microbe-nematode interactions	[71]
Trait-Mediated TIM	Effect driven by changes in functional traits	Microbial trait shifts under predation pressure	[71]
Integrated TIM Strength	Combined effect size across multiple functions	Multifunctionality response to cross-trophic interactions	[5] [72]
Network Fragmentation Index	Structural changes following species loss	Food web robustness to habitat-specific extinctions	[73]

The application of these metrics in soil systems requires careful experimental design. For example, relative modification strength can be quantified by comparing organic matter decomposition rates in presence and absence of specific predator groups, while network fragmentation indices have been used to model how the loss of wetland-associated species disproportionately disrupts regional food web connectivity [73].

Experimental Data on Soil Biodiversity-Function Relationships

Recent empirical studies have generated quantitative data on the relationships between soil biodiversity and ecosystem functions, highlighting how trophic interactions moderate these relationships across different ecosystems.

Table 2: Experimental Evidence for Biodiversity-Function Relationships Moderated by Trophic Interactions

Study System	Biodiversity Measure	Ecosystem Function(s)	Key Finding	Trophic Interaction Role
Global drylands [10]	Bacterial & fungal diversity	Multifunctionality (C, N, P cycling)	Positive correlation (R²=0.62)	Microbial trophic complementarity
Urban greenspaces [72]	Cross-domain soil biodiversity	18 ecosystem functions	Soil biodiversity > plant diversity as predictor	Multi-trophic interactions drive multifunctionality
Agricultural soils [5]	Microbial & faunal diversity	Organic matter decomposition	2.4x faster with diverse decomposers	Functional redundancy across trophic levels
Swiss metaweb [73]	Network connectance	Food web robustness	Targeted species loss increased fragmentation	Habitat-specific trophic dependencies
Forest soils [74]	Microbial community composition	Nutrient cycling rates	30% variation explained by community structure	Trophic cascade effects on process rates

The data from global studies consistently show that soil biodiversity accounts for significant variation in ecosystem functioning, with coefficients of determination (R²) ranging from 0.42 to 0.76 for different functions [10] [72]. Importantly, these relationships are often non-linear, with threshold effects observed where biodiversity loss beyond certain points leads to disproportionate functional declines [72]. The strength of these relationships is consistently moderated by the complexity of trophic interactions, with more diverse food webs generally exhibiting greater functional stability and resilience to perturbations [75] [73].

Methodological Approaches: Experimental Protocols

Mesocosm Experiments with Trophic Manipulations

Controlled mesocosm experiments represent a powerful approach for disentangling direct and indirect trophic effects in soil systems. These experiments typically involve manipulating the presence/absence or diversity of specific trophic groups while measuring resulting ecosystem processes.

Protocol 1: Trophic Exclusion Experiment

• Objective: Quantify the relative contributions of different trophic groups to ecosystem functions
• Procedure:
  • Establish soil microcosms with standardized soil properties
  • Apply selective biocides to eliminate specific trophic groups (e.g., antibiotics for bacteria, fungicides for fungi, nematicides for nematodes)
  • Inoculate with characterized microbial and faunal communities
  • Measure process rates over time (decomposition, nutrient mineralization, respiration)
  • Use statistical partitioning to attribute functions to different trophic groups

Applications: This approach has revealed that microbial communities (bacteria and fungi) contribute approximately 60-75% of decomposition and nutrient cycling functions, while soil fauna contribute 25-40% through regulation and dispersal activities [5]. The protocol allows researchers to quantify both the direct effects of excluded groups and the indirect effects mediated through trophic cascades.

Molecular Characterization of Trophic Interactions

Advanced molecular techniques now enable direct tracking of energy flows and trophic relationships in soil food webs, moving beyond inferences based on coexistence.

Protocol 2: Stable Isotope Probing (SIP) with High-Throughput Sequencing

• Objective: Directly link taxonomic identity to trophic function and resource use
• Procedure:
  • Label substrates with stable isotopes (¹³C, ¹⁵N)
  • Introduce labeled substrates to soil mesocosms
  • Harvest samples at multiple time points
  • Separate density-fractionated nucleic acids (DNA-SIP) or analyze bulk tissue
  • Sequence marker genes (16S rRNA for bacteria, 18S rRNA for eukaryotes, ITS for fungi)
  • Quantify isotope incorporation to identify active consumers

Applications: SIP approaches have revealed that approximately 20-30% of bacterial and fungal taxa account for 70-80% of organic matter processing in soil, challenging assumptions of high functional redundancy [6] [5]. This method has been particularly valuable for identifying keystone taxa that disproportionately influence ecosystem processes through their trophic activities.

Network Reconstruction from Metawebs

Regional-scale metawebs provide a framework for predicting how species losses affect food web structure and function across habitat types.

Protocol 3: Metaweb Inference and Robustness Testing

• Objective: Predict food web responses to non-random species loss scenarios
• Procedure:
  • Compile comprehensive database of potential trophic interactions (metaweb)
  • Document species distributions across habitat types
  • Infer realized food webs for specific regions/habitats
  • Simulate extinction scenarios based on species habitat associations and abundances
  • Calculate robustness metrics (secondary extinctions, network fragmentation)

Applications: Application of this protocol to the Swiss trophiCH metaweb (7808 species, 281,023 interactions) demonstrated that targeted loss of wetland-associated species caused 45% greater network fragmentation than random species loss, highlighting the disproportionate importance of specific habitats for maintaining food web integrity [73].

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Key Research Reagents for Disentangling Trophic Interactions in Soils

Reagent/Material	Primary Function	Application Example	Technical Considerations
Selective biocides	Trophic group exclusion	Isolating microbial vs. faunal contributions	Non-target effects must be quantified
Stable isotope labels (¹³C, ¹⁵N)	Tracking energy flow	Identifying active decomposers via SIP	Resolution depends on incorporation time
DNA/RNA extraction kits	Community characterization	Molecular analysis of food web structure	Yield varies with soil type and organism size
Species-specific primers	Taxonomic identification	Quantifying predator-prey relationships	Specificity must be empirically validated
Microcosm/mesocosm setups	Controlled experimentation	Manipulating trophic complexity	Scaling from microcosms to field conditions
Metaweb databases	Food web reconstruction	Predicting responses to species loss	Distinguishing potential vs. realized interactions

Integrated Conceptual Framework of Soil Food Web functioning

The complex relationships between soil biodiversity, trophic interactions, and ecosystem functions can be integrated into a conceptual framework that highlights both direct and indirect pathways of influence:

Diagram 2: Integrated Framework of Soil Food Web Functioning. This conceptual model illustrates how soil biodiversity directly and indirectly (via trophic interactions) influences ecosystem multifunctionality, with modifying factors affecting these relationships.

Knowledge Gaps and Future Research Directions

Despite significant advances, critical knowledge gaps remain in our understanding of trophic interactions in soil ecosystems. Current global soil biodiversity and ecosystem function research exhibits important spatial, environmental, taxonomic, and functional biases, with only 0.3% of sampling sites having concurrent biodiversity and function measurements [10]. These limitations constrain our ability to develop general models of soil food web functioning.

Priority research directions include:

• Temporally explicit studies: Current data overwhelmingly represent single time points, limiting understanding of feedback dynamics [10]
• Integration across scales: Linking micro-scale interaction data with ecosystem-level processes [75]
• Standardized methodologies: Developing comparable protocols for quantifying interaction strengths across systems [71]
• Network stability assessment: Moving beyond structural analysis to dynamic robustness testing [73]

Addressing these gaps will require coordinated international efforts, such as the global soil biodiversity monitoring networks, and the application of emerging technologies including environmental DNA meta-barcoding, sensor networks, and machine learning approaches to predict interaction outcomes under changing environmental conditions.

The long-held paradigm in ecology that "more biodiversity is always better" is being fundamentally challenged by cutting-edge research. Modern soil and ecosystem science reveals a far more complex reality: the relationships between biodiversity and ecosystem functioning are characterized by significant trade-offs, non-linear responses, and critical thresholds. Understanding these complex dynamics is crucial for predicting ecosystem responses to global change and for developing effective conservation and restoration strategies.

This paradigm shift recognizes that while biodiversity generally enhances ecosystem multifunctionality—the simultaneous performance of multiple ecosystem functions—these relationships are rarely straightforward or linear [76]. Instead, threshold effects and context dependencies dictate when and how biodiversity supports functions ranging from carbon sequestration to nutrient cycling [77] [39]. This article synthesizes recent experimental evidence to guide researchers in testing these complex biodiversity-ecosystem functioning relationships in soils, with particular emphasis on methodological approaches and analytical frameworks.

Experimental Evidence: Key Studies and Findings

Quantitative Evidence of Non-Linear Relationships

Table 1: Documented Non-Linear Responses and Threshold Effects in Ecosystem Studies

Study System/Location	Key Driver(s)	Ecosystem Function(s)	Threshold/Non-linear Pattern	Citation
City Belt along Yellow River, Ningxia	Landscape fragmentation (CONTAG, SHDI, FRAC_MN)	Crop production, carbon sequestration, nutrient retention, sand fixation	3 km identified as optimal scale for analysis; trade-offs altered at high/low fragmentation	[77]
Atlantic Forest, Brazil	Soil biodiversity (multi-trophic)	Primary production, ecosystem stability, nutrient cycling	Land-use simplification reduces multifunctionality via biodiversity loss	[18]
Hubei Province, China	Urbanization intensity	Water yield, carbon storage, soil conservation, food supply	Spatial heterogeneity in trade-offs: CS-SC-NPP synergies vs. CS/SC-NPP-FS trade-offs	[78]
Nationwide grassland-afforestation gradient	Soil pH, C:N ratio, leaf dry matter content	Microbial taxonomic/functional diversity	Threshold transitions between microbial generalists/specialists	[39]
BioCON experiment, Minnesota	Plant diversity under elevated CO₂ and N	Root biomass, soil respiration, microbial biomass, soil aggregation	Weakened biodiversity-multifunctionality relationships under N enrichment	[76]

Documented Trade-offs and Synergies Among Ecosystem Functions

Table 2: Experimentally Documented Trade-offs and Synergies Among Ecosystem Functions

Ecosystem Functions Compared	Relationship Type	Context Dependencies	Citation
Agricultural production vs. Carbon sequestration	Synergy	Becomes trade-off in mountain regions with higher/lower fragmentation	[77]
Sand fixation vs. Regulating services	Conditional synergy	Altered under higher/lower diversity conditions	[77]
Carbon storage vs. Food supply	Trade-off	Consistent across Hubei Province landscapes	[78]
Functional diversity vs. Functional redundancy	Trade-off	Microbial genetic specialization during ecosystem development	[39]
Soil respiration vs. Microbial biomass	Positive correlation	Strong association (r² = 0.45) across environmental contexts	[76]

Methodological Approaches for Detecting Non-Linearities

Experimental Designs and Scaling Considerations

Detecting non-linear responses requires carefully designed experiments and appropriate spatial scaling. The optimal spatial scale for analyzing landscape pattern impacts on ecosystem services has been identified at 3 km across a range of 1.5 km to 30 km [77]. This scale captures the relevant ecological processes without introducing excessive noise or oversimplification.

Paired-site approaches that establish gradients of management intensity (e.g., managed grasslands to abandoned and afforested sites) effectively reveal threshold dynamics [39]. These designs should incorporate multiple successional stages to capture transitions rather than simple endpoint comparisons. For plant diversity experiments, orthogonal manipulations of global change factors (e.g., elevated CO₂, nitrogen enrichment) enable testing of biodiversity-ecosystem functioning relationships across future environmental scenarios [76].

Analytical Frameworks for Non-Linear Relationships

• Random Forest Analysis: Effectively quantifies relative contributions of multiple drivers and identifies dominant variables through built-in importance metrics [77] [79].

• Partial Dependence Analysis (PDA): Characterizes non-linear response curves and detects impact thresholds for dominant drivers, revealing "single threshold," "monotonic impact," and "complex curve" effects (including S-shape, inverted U-shape, and inverted S-shape) [79].

• Multivariable Regression Trees: Identifies critical thresholds in multivariate data where relationships between drivers and ecosystem functions shift abruptly [77].

• Structural Equation Modeling: Uncovers direct and indirect pathways through which biodiversity influences multifunctionality, separating effects of species richness, functional diversity, and evenness [76].

• Multiple Threshold Approach: Assesses ecosystem multifunctionality across different performance thresholds (e.g., 20%, 40%, 60%, 80%) to determine how many functions are maintained above critical levels [76].

Figure 1: Experimental workflow for detecting non-linear responses in ecosystem functions, showing the progression from study design through analysis to pattern identification.

The Scientist's Toolkit: Essential Research Solutions

Table 3: Key Research Reagents and Methodological Solutions for Ecosystem Function Studies

Tool/Category	Specific Examples	Primary Application	Experimental Consideration
Ecological Modeling	InVEST model water yield module	Quantifying water yield services based on water balance principle	Requires precipitation, evapotranspiration, and soil data inputs	[78]
Landscape Metrics	Fragstats 4.3	Calculating landscape-level metrics (CONTAG, SHDI, FRAC_MN)	3 km scale recommended for ecosystem service analyses	[77]
Molecular Tools	16S/ITS amplicon sequencing; Metagenomics	Assessing microbial taxonomic/functional diversity	Enables tracking of genetic redundancy vs. specialization	[39]
Soil Bioassays	Soil respiration measurements; Microbial biomass C; Decomposition rates	Quantifying multiple soil functions simultaneously	Critical for multifunctionality assessments	[76]
Statistical Packages	R packages for Random Forest, PDA, SEM	Analyzing non-linear relationships and threshold effects	Partial Dependence Analysis essential for threshold detection	[79]
Field Instruments	Soil core samplers; Aggregate stability kits; Root scanners	Standardized measurement of root biomass and soil structure	Enables cross-study comparisons	[76]

Mechanisms Underlying Non-linear Responses

Soil Microbial Transitions During Ecosystem Development

Recent research on land abandonment reveals threshold dynamics in soil microbial communities during the grassland-to-forest transition. These transitions involve:

• Functional-Redundancy Trade-offs: Afforestation leads to increasing functional diversity but decreasing taxonomic diversity, creating a putative trade-off between two desirable ecosystem properties [39]. This specialization decreases genetic redundancy as microbial communities adapt to more complex, recalcitrant plant litter.

• Carbon Cycling Specialization: Fungal functional diversity underpins higher microbial carbon-cycling capacity, with threshold effects notable in the response of fungal C-cycling genes, which increase in forests compared to grasslands [39].

• Abiotic Trigger Points: Thresholds in diversity changes coincide with a sharp decline in soil pH, increasing soil C:N ratio, and higher levels of leaf dry matter content, with soil pH identified as the primary determinant of bacterial diversity decline during afforestation [39].

Plant-Soil Feedback Mechanisms

Plant-soil feedback represents a critical mechanism generating non-linear responses in ecosystem functions:

• Soil Legacy Effects: Plants condition soil abiotic and biotic properties via litter inputs and root activity, creating legacies that influence subsequent plant performance [80]. These feedbacks can be positive or negative and strongly influence plant community dynamics.

• Environmental Context Dependency: The strength and direction of plant-soil feedback depends on environmental factors such as nitrogen and water availability [80]. Understanding this context-dependency is essential for predicting ecosystem responses to global change.

• Trait-Based Frameworks: Incorporating root trait gradients of the root economics spectrum enhances predictability of plant-soil feedback effects [80]. This functional trait approach provides mechanistic understanding beyond correlational patterns.

Figure 2: Mechanism of threshold effects in ecosystem multifunctionality, showing how environmental drivers trigger soil property shifts that alter microbial communities, leading to non-linear responses in individual functions and overall multifunctionality.

Understanding the non-linear responses and trade-offs in ecosystem functions provides a more realistic framework for predicting how ecosystems will respond to global environmental change. The experimental evidence synthesized here demonstrates that biodiversity-ecosystem functioning relationships are context-dependent and often exhibit critical thresholds beyond which system behavior changes dramatically.

For researchers designing experiments to test these relationships, key recommendations emerge: (1) incorporate environmental gradients rather than simple comparisons; (2) measure multiple functions simultaneously to detect trade-offs; (3) employ analytical methods capable of detecting non-linearities and thresholds; and (4) consider both taxonomic and functional dimensions of biodiversity. As we face unprecedented environmental change, recognizing that "more isn't always better" provides the nuanced understanding necessary to conserve and manage ecosystems effectively.

Robust assessments of soil biodiversity and its relationship to ecosystem functioning (BEF) hinge on optimized sampling strategies that account for depth, scale, and temporal dynamics. Soils represent one of Earth's most complex biological habitats, harboring immense biodiversity that drives crucial ecosystem functions including nutrient cycling, carbon sequestration, and climate regulation [81]. However, significant gaps persist in global soil biodiversity data due to methodological inconsistencies and biogeographical biases across studies [81]. Research demonstrates that BEF relationships generally strengthen over time in both grassland and forest ecosystems, with soil characteristics significantly influencing these temporal trajectories [82]. This guide systematically compares sampling approaches to address these challenges, providing researchers with evidence-based protocols for generating comparable, high-quality data across studies and ecosystems. By optimizing sampling strategies across spatial and temporal dimensions, scientists can more accurately quantify how soil biodiversity sustains ecosystem functioning amid global environmental change.

Comparing Core Sampling Strategies: Methodologies and Applications

Depth-Specific Sampling Approaches

Table 1: Depth-Specific Sampling Recommendations for Different Ecosystems

Ecosystem/Crop Type	Sampling Depth	Rationale	Key Functions Assessed
Shallow-rooted crops (rice, millet, groundnut)	0-15 cm (6 inches)	Captures majority of root activity and associated rhizosphere processes [83]	Nutrient uptake, microbial activity
Deep-rooted crops (cotton, sugarcane, orchard crops)	Multiple depths: 0-30 cm, 30-60 cm, 60-90 cm	Profiles nutrient distribution across root zones [83]	Water retention, nutrient leaching, carbon storage
Grasses and grasslands	0-5 cm (2 inches)	Focuses on organic-rich surface layer [83]	Surface stability, organic matter decomposition
General agricultural assessments	0-6 inches (standard)	Balances practicality with representation of plow layer [84]	Fertility assessment, pH, macronutrients
Construction/contamination studies	12-24 inches	Reaches below surface to assess subsurface risks [84]	Contaminant profiling, structural stability

Depth-specific sampling involves collecting soil at predetermined intervals to profile how properties and processes vary vertically through the soil column. This approach reveals that nutrient availability, microbial communities, and soil constraints differ significantly across depths, requiring tailored sampling protocols for different research objectives [84] [83]. For instance, standard agricultural assessments typically target the 0-6 inch depth where most nutrient uptake occurs, while construction and contamination studies require sampling down to 12-24 inches to assess subsurface risks [84]. Research on perennial crops and orchards benefits from multi-depth sampling at 30, 60, and 90 cm to capture root activity across the soil profile [83].

Spatial Sampling Designs

Table 2: Comparison of Spatial Sampling Strategies

Method	Protocol	Optimal Use Cases	Efficiency Metrics	Limitations
Grid Sampling	Divides field into equal cells (e.g., 1-acre); samples at intersections [84] [83]	High-precision agriculture; heterogeneous fields	Captures 80% of spatial variability [84]	Cost-intensive: $500-$1,000 per 100 acres [84]
Zone Sampling	Delineates areas with similar properties (soil type, yield history); samples each zone separately [83]	Sites with known management history; distinct soil zones	Improves accuracy by 20% in heterogeneous fields [84]	Requires prior knowledge; GIS data dependency
Composite Sampling	Combines 10-15 subsamples per area into one representative sample [84] [83]	Large-area assessments; initial screening	Reduces analysis costs by 30% while maintaining 90% accuracy [84]	Masks micro-scale variability
Transect Sampling	Linear sampling along environmental gradients [84]	Gradient studies (e.g., pollution plumes, moisture gradients)	Detects 15% more variation than random points [84]	Limited areal coverage
Stratified Random Sampling	Splits area into strata based on soil type/land use; samples each stratum [84]	Region-scale studies; heterogeneous landscapes	Cost-effective for large areas [84]	Complex statistical design required

Spatial sampling strategies range from highly systematic grid-based approaches to more targeted zone sampling, each with distinct advantages depending on research objectives, spatial heterogeneity, and available resources [84] [83]. Grid sampling provides high-resolution data but at greater cost, while composite sampling offers a practical balance between cost-efficiency and representativeness for larger areas [84]. The selection of an appropriate spatial design fundamentally influences the detection of biodiversity-ecosystem functioning relationships, as soil properties can vary by 15-20% within individual fields [84].

Temporal Sampling Considerations

Table 3: Temporal Dynamics in Biodiversity-Ecosystem Functioning Relationships

Ecosystem Type	Temporal Pattern	Key Mechanisms	Sampling Implications
Grassland Systems	Strengthening positive BEF relationships over 17+ years [30]	Increasing complementarity effects; species asynchrony [30]	Multi-year studies essential (minimum 3 years) [82]
Forest Ecosystems	Consistent positive diversity effects on biomass accumulation [82]	Early complementarity effects; niche partitioning [82]	Longer duration required (minimum 5 years) [82]
Soil Microbial Communities	Diversity stabilizes ecosystem functioning via species asynchrony [42]	Different taxa support different functions at different times [42]	Seasonal sampling across growing cycles; multiple time points
Agricultural Systems	Post-harvest/pre-planting optimal [83]	Avoids management practice distortions [83]	Consistent annual timing for comparability

Temporal dynamics introduce critical considerations for sampling regime design. Long-term studies reveal that biodiversity-ecosystem functioning relationships typically strengthen over time, with grassland experiments demonstrating increasingly positive diversity-productivity relationships over 17-year periods [30]. These temporal patterns emerge through different mechanisms across ecosystems: in grasslands, diverse communities show greater stability while low-diversity communities exhibit functional declines, whereas forests demonstrate consistent positive diversity effects on biomass accumulation [82]. Soil microbial communities stabilize ecosystem functioning through asynchrony, where different taxa support different functions at different times [42]. These findings underscore the necessity of multi-year studies with consistent seasonal timing to accurately capture BEF relationships.

Integrated Methodological Framework

Experimental Protocols for Robust BEF Assessment

Protocol 1: Depth-Explicit Soil Biodiversity Sampling

• Site Characterization: Document land use history, vegetation cover, and soil type using standardized forms [84]
• Sample Collection: Using sterile soil probes or augers, collect minimum 10-15 cores per defined area [84] [83]
• Depth Stratification: Process samples by predetermined depth intervals (e.g., 0-10 cm, 10-20 cm, 20-30 cm) using the V-shaped method for precise depth control [83]
• Composite Strategy: For homogeneous areas, combine depth-specific samples from multiple cores; maintain separate samples for heterogeneous zones [84]
• Preservation: For molecular analyses, immediately freeze samples at -20°C; for chemical analyses, air-dry in shade within 24 hours [84] [83]

Protocol 2: Temporal Sampling for BEF Dynamics

• Baseline Establishment: Conduct comprehensive initial sampling characterizing all biodiversity compartments and ecosystem functions [30]
• Seasonal Captures: Sample at key phenological stages (e.g., early growth, peak biomass, senescence) to capture intra-annual variation [42]
• Inter-annual Replication: Maintain consistent protocols across years with sampling at the same seasonal points [83] [30]
• Event-Responsive Sampling: Implement additional sampling following extreme events (drought, flooding) to assess resilience [30]

Protocol 3: Spatial Design Implementation

• Pre-Sampling Stratification: Using historical data, aerial imagery, or geophysical surveys, delineate sampling zones [84] [83]
• Grid Establishment: For precision studies, establish permanent grid points with GPS coordinates (accuracy <3m) [84]
• Sample Intensity Determination: Allocate 10-20 samples per 20 acres to capture 85% of variability in heterogeneous systems [84]
• Reference Samples: Include adjacent natural ecosystems as reference points where possible [85]

Workflow Visualization: Integrated Sampling Design

Research Reagent Solutions and Essential Materials

Table 4: Essential Research Reagents and Equipment for Soil BEF Studies

Category	Specific Products/Tools	Research Function	Application Notes
Sampling Equipment	Soil probes; augers; stainless steel corers; GPS devices	Standardized collection; precise spatial referencing	Use sterile tools to prevent cross-contamination (reduces errors by 10%) [84]
Preservation Solutions	DNA/RNA stabilization buffers; cryogenic containers; silica gel	Maintains sample integrity for molecular and chemical analyses	Store at 4°C to preserve microbial community structure [84]
Molecular Analysis Kits	DNA extraction kits (e.g., PowerSoil Kit); PCR reagents; sequencing libraries	Biodiversity assessment via metabarcoding/eDNA	Standardized across samples for cross-study comparisons [86]
Chemical Analysis Reagents	Spectrometry standards; pH buffers; nutrient extraction solutions	Quantification of soil properties and ecosystem functions	70% of labs use spectrometry for nutrient analysis (95% accuracy) [84]
Stable Isotope Tracers	¹⁵N, ¹³C labeled compounds; isotope ratio mass spectrometry	Tracking nutrient pathways and biogeochemical processes	Critical for quantifying process rates in BEF studies [87]

Discussion: Integration and Future Directions

Optimizing sampling strategies requires careful balancing of practical constraints with scientific objectives. While comprehensive sampling across depth, space, and time provides the most robust BEF assessments, researchers must often make strategic tradeoffs. Composite sampling reduces costs by 30% while maintaining 90% accuracy, offering a practical approach for large-scale studies [84]. Emerging technologies like eDNA metabarcoding are revolutionizing soil biodiversity assessment, with ongoing research optimizing core sampling strategies for capturing invertebrate diversity through composite sampling approaches [86].

Critical gaps remain in global soil biodiversity data, particularly in tropical regions and for specific taxonomic groups [81]. Future methodological development should focus on standardizing protocols across studies to enable meaningful cross-site comparisons. Integrating remote sensing with targeted ground sampling represents a promising approach for scaling BEF relationships across landscapes. Furthermore, research must address the temporal mismatch between rapid microbial community dynamics and slower ecosystem processes by developing multi-scale temporal sampling frameworks [42].

By implementing these optimized sampling strategies—carefully considering depth stratification, spatial design, and temporal dynamics—researchers can generate the high-quality, comparable data needed to advance our understanding of critical biodiversity-ecosystem functioning relationships in soil environments.

Validation and Synthesis: Cross-System Comparisons and Emerging Universal Principles

The relationship between biodiversity and ecosystem functioning (BEF) is a cornerstone of modern ecology, with critical implications for conservation, restoration, and climate change mitigation. While BEF relationships have been demonstrated across diverse ecosystems, understanding how these dynamics differ between major biome types remains essential for developing targeted management strategies. This review synthesizes current evidence to compare BEF dynamics in grassland versus forest ecosystems, with particular focus on patterns across successional gradients and under different land-use types. We examine how fundamental differences in physiology, life history strategies, and temporal dynamics between these ecosystems shape their functional response to biodiversity change, providing a framework for predicting ecosystem behavior in an era of global change.

Theoretical Foundations of BEF Relationships

BEF research investigates how the variety and abundance of species influence ecosystem processes such as productivity, nutrient cycling, and stability. Two primary mechanisms underlie positive biodiversity effects: complementarity effects, where diverse species partitions resources more efficiently through niche differentiation, and selection effects, where dominant species with particular traits drive ecosystem processes [88] [89]. The relative importance of these mechanisms varies between ecosystems and across temporal scales.

In grasslands, complementarity often arises from differences in rooting depth, photosynthetic pathways, or phenological patterns that allow species to utilize resources at different times or spaces [88]. In forests, niche differentiation frequently occurs through canopy stratification, variation in shade tolerance, or differences in leaf phenology [89]. The temporal development of these mechanisms differs substantially between these ecosystems due to their distinct life history strategies and developmental trajectories.

Comparative Analysis of BEF Dynamics

Temporal Dynamics and Successional Patterns

Table 1: Comparative Temporal Dynamics of BEF Relationships in Grasslands and Forests

Aspect	Grassland Ecosystems	Forest Ecosystems
Temporal strengthening	Effects strengthen over 1-10 years [88]	Effects strengthen over decades [88] [89]
Early succession drivers	Acquisitive species with fast resource capture [88]	Light-capturing species with rapid height growth [88]
Late succession drivers	Conservative species with resource efficiency [88]	Shade-tolerant species with high resource acquisition capacity [88]
Key mechanisms	Rapid changes in population densities [88]	Slow processes of individual growth and canopy space filling [88]
Successional study approach	Direct observation of succession [90]	Chronosequence studies and modeling [90] [89]

The temporal dynamics of BEF relationships differ markedly between grasslands and forests. In grassland ecosystems, biodiversity effects on productivity typically strengthen over relatively short timeframes (1-10 years) [88]. This rapid response is facilitated by the fast life cycles of herbaceous species and their ability to quickly adjust population densities in response to changing conditions [88]. During early succession, fast-growing acquisitive species with traits like high specific leaf area (SLA) and leaf nitrogen content (LNC) dominate biodiversity effects, while in later successional stages, resource-conservative species with greater efficiency in resource use become increasingly important contributors to overyielding [88].

In contrast, forest ecosystems exhibit much slower development of BEF relationships, with effects strengthening over decades rather than years [88] [89]. The slow growth and longevity of trees means that early successional dynamics are dominated by light-capturing species with rapid height growth, while later successional stages see increasing contributions from shade-tolerant species with high soil resource acquisition capacity [88]. The slower dynamics in forests necessitate alternative research approaches, including chronosequence studies that substitute space for time [90] and modeling techniques that project short-term observations across full rotation periods [89].

Trait-Based Mechanisms and Functional Diversity

Plant functional traits provide a mechanistic link between biodiversity and ecosystem functioning, with different traits being important in grassland versus forest ecosystems.

Table 2: Key Functional Traits Influencing BEF Relationships in Grasslands and Forests

Trait Category	Grassland Ecosystems	Forest Ecosystems
Aboveground traits	Specific Leaf Area (SLA), Leaf Dry Matter Content (LDMC), Leaf Nitrogen Content (LNC) [88]	Canopy structure, leaf phenology, shade tolerance [89]
Belowground traits	Specific Root Length (SRL), Root Nitrogen Content (RNC) [88]	Rooting depth, mycorrhizal associations [88]
Resource strategy	Acquisition-conservation tradeoff along economics spectrum [88]	Light capture vs. soil resource acquisition tradeoffs [88]
Response to mixing	Plasticity in rooting patterns and nutrient uptake [88]	Crown plasticity, vertical stratification [89]

In grasslands, the resource economics spectrum—ranging from acquisitive to conservative strategies—strongly influences species contributions to BEF relationships [88]. Acquisitive species with high SLA, LNC, SRL, and RNC typically dominate early biodiversity effects, while conservative species with opposite trait values become increasingly important over time [88]. This temporal shift in species importance reflects changing resource availability and competitive interactions as communities develop.

Forest BEF relationships are strongly influenced by traits related to light capture (e.g., canopy structure, phenology) and belowground resource acquisition (e.g., rooting depth, mycorrhizal associations) [88] [89]. The overyielding of species with high light capture capacity tends to decrease over time, while overyielding of species with high soil resource acquisition capacity increases [88]. This reflects the changing light environment as canopies close and increased importance of belowground resources as stands mature.

Methodological Approaches in BEF Research

Table 3: Research Methods for Studying BEF Relationships in Grasslands and Forests

Method	Grassland Applications	Forest Applications
Experimental designs	Manipulative diversity experiments [88]	Tree diversity experiments [89]
Observational approaches	Permanent plots across successional gradients [90]	Triplet studies (paired monoculture-mixture stands) [89]
Modeling techniques	Not specified in results	Forest gap models [89]
Temporal scope	Direct observation over years to decades [88]	Modeling projections over decades to centuries [89]
Scale considerations	Small plots (m² to hm²) [31] [3]	Large plots (hm²) [31]

Grassland BEF research has heavily relied on manipulative experiments where species richness is directly controlled, allowing for strong causal inference [88]. These experiments benefit from the relatively small size and fast growth of grassland species, enabling researchers to test multiple diversity levels with adequate replication over meaningful timeframes [3]. The ORPHEE experiment, for instance, demonstrates how combining experimental data with modeling can explore long-term effects of diversity [89].

Forest BEF research employs more diverse methodologies due to the practical challenges of manipulating tree diversity. While tree diversity experiments have been established [89], they are typically younger and smaller than their grassland counterparts. Observational approaches like triplet studies (comparing monocultures and mixtures) [89] and chronosequence studies [90] are common. Modeling approaches, particularly forest gap models, have proven valuable for projecting diversity effects across full rotation periods [89]. These models simulate community dynamics by modeling establishment, growth, and mortality of individual trees in small patches, allowing virtual experiments over century timescales [89].

Context Dependency and Global Change Implications

BEF relationships in both grasslands and forests are strongly context-dependent, influenced by abiotic factors, land use history, and global change drivers. In drylands, which include both grassland and forest ecosystems, BEF relationships are particularly sensitive to climate change and human activities [24]. Soil biodiversity plays a crucial role in mediating BEF relationships across ecosystems, with soil organisms supporting multiple ecosystem functions including carbon sequestration, nutrient cycling, and plant productivity [15] [3] [6].

Climate change alters BEF dynamics by shifting temperature and precipitation regimes, which subsequently affect species interactions and ecosystem processes [91] [92]. In dryland shrublands and woodlands of the western US, climate warming has been associated with vegetation changes that decrease ecological resilience and increase fire risk [92]. Such changes may fundamentally alter BEF relationships by favoring species with different functional strategies and changing the nature of species interactions.

Land use intensification represents another major driver of BEF relationships, often reducing biodiversity and simplifying ecosystem structure in both grasslands and forests [91]. Understanding how management practices can optimize BEF relationships while meeting human needs represents a critical research frontier.

Research Gaps and Future Directions

Despite significant advances in BEF research, important knowledge gaps remain in comparing grassland and forest ecosystems. Dryland ecosystems, which encompass both grasslands and forests, are particularly understudied regarding BEF relationships in bare soils devoid of perennial vegetation [24]. The role of intra-specific trait variability in modulating BEF relationships represents another understudied area across ecosystems [24]. Similarly, understanding how biotic interactions beyond competition—including plant-animal interactions and microbial relationships—influence BEF relationships requires additional research [24]. Finally, integrating temporal variability and human activities into BEF research remains challenging but essential for predicting real-world ecosystem dynamics [24].

The Scientist's Toolkit

Table 4: Essential Research Reagents and Methodologies for BEF Studies

Tool Category	Specific Methods/Technologies	Application in BEF Research
Field Survey	CTFS (Center for Tropical Forest Science) protocols [31]	Standardized measurement of tree growth and survival
Trait Measurements	Specific Leaf Area (SLA), Leaf Dry Matter Content (LDMC) [88]	Quantifying plant functional strategies
Soil Analysis	Chloroform fumigation extraction [3]	Measuring microbial biomass nitrogen
Molecular Techniques	"Omics" approaches (e.g., metagenomics) [15] [6]	Characterizing soil microbial diversity
Remote Sensing	Satellite imagery analysis [92]	Tracking vegetation change over large scales
Modeling Tools	Forest gap models (e.g., ForCEEPS) [89]	Projecting long-term forest dynamics

Conceptual Diagrams

Graph 1: Comparative BEF dynamics show grassland ecosystems characterized by rapid, population-level adjustments, while forests exhibit slower dynamics driven by individual-level processes.

Graph 2: Successional shifts show divergent patterns, with grasslands and forests displaying different trajectories in species contributions to ecosystem functioning over time.

Grassland and forest ecosystems demonstrate fundamentally different BEF dynamics shaped by their distinct biological attributes, temporal scales, and methodological constraints. Grasslands exhibit rapid BEF relationships driven by population-level adjustments and clear shifts from acquisitive to conservative species, while forests display slower dynamics governed by individual-level processes and transitions from light-capturing to soil resource-acquiring species. These differences necessitate distinct research approaches, with grasslands favoring manipulative experiments and forests relying more on modeling and observational studies. Understanding these divergent patterns is crucial for developing ecosystem-specific management strategies that maintain biodiversity and ecosystem functioning under global change. Future research should address critical knowledge gaps regarding intra-specific trait variation, multi-trophic interactions, and the impacts of global change drivers on BEF relationships across ecosystems.

The relationship between soil biodiversity and ecosystem functioning (BEF) is a cornerstone of soil ecological research. A critical test of this relationship lies in comparing resource-rich (fertile) and resource-poor (nutrient-limited) environments. In fertile soils, high nutrient availability can allow a few dominant, fast-growing organisms to drive key processes, suggesting a degree of functional redundancy. In contrast, theory and growing evidence indicate that in nutrient-poor soils, the diversity of soil life becomes paramount; a wider range of species with complementary functions is required to capture and recycle scarce resources efficiently [6] [93]. This article compares soil biodiversity-ecosystem function relationships under contrasting fertility conditions, validating the critical role of biological diversity in low-resource environments through experimental data and mechanistic insights.

---

Comparative Analysis: Biodiversity and Ecosystem Functioning

The following tables synthesize quantitative data from key studies, comparing biodiversity metrics and ecosystem functions between fertile and nutrient-poor soils under various experimental conditions.

Table 1: Soil Biota Diversity and Ecosystem Multifunctionality Under Nutrient Enrichment

This table summarizes findings from a long-term gradient nutrient addition experiment, demonstrating how nutrient-induced changes impact biodiversity and function [94].

Metric / Parameter	Unfertilized Control (Low Fertility)	Medium NP Addition (NP90)	High NP Addition (NP120) (High Fertility)
Soil Chemical Properties
Soil pH	7.20	Significantly decreased	6.54
Soil Mineral N (mg N·kg⁻¹)	3.67	Significantly increased	7.78
Soil Available P (mg P·kg⁻¹)	8.12	Significantly increased	173.23
Soil Biota Diversity (Shannon Index)
Bacteria Diversity	Highest	Significantly decreased	Significantly decreased
Fungi Diversity	Highest	Significantly decreased	Significantly decreased
Total Nematode Diversity	Highest	Significantly decreased	Significantly decreased
Soil Biota Multidiversity	Highest	Significantly decreased	Significantly decreased
Ecosystem Functions & Multifunctionality
C & N Cycling Functions (e.g., Microbial Biomass C & N)	Highest	Significantly decreased	Significantly decreased
Ecosystem Multifunctionality (EMF) Index	Highest (0% reduction)	Decreased by 28%	Decreased by 36%

Table 2: Soil Community and Functional Responses to Legume Intercropping in a Nutrient-Poor System

This table details how introducing legume plants into a degraded, nutrient-poor ecosystem alters the soil micro-food web and related functions, highlighting a phosphorus-limited context [93].

Parameter / Relationship	Camellia Oleifera Monoculture (Nutrient-Poor)	With Legume Intercropping (e.g., Peanut)
Soil Nutrient Context	Strong soil P limitation	P limitation aggravated by legume input
Soil Biota Response
Bacterial Diversity	Lower	Increased
Protist Consumer Diversity	Lower	Increased
Omnivore-Predator Nematode Abundance	Lower	Increased
Trophic Cascade & Function
Bottom-Up Effect	Weaker	Stronger: P limitation → ↑ Bacterial diversity → ↑ Protists & Nematodes
Top-Down Effect	Weaker	Stronger: Omnivore-predator nematodes & protists ↓ soil P metabolism via bacteria
Overall P-related Function	Baseline	Down-regulated

---

Detailed Experimental Protocols

To validate the critical role of biodiversity in low-resource environments, researchers employ carefully designed experiments. Below are the methodologies for two key studies cited in this article.

Protocol 1: Long-Term Nutrient Enrichment Experiment

This protocol is designed to test the mechanisms by which nutrient enrichment affects soil biodiversity and ecosystem multifunctionality [94].

• 1. Site Description: The experiment was established in a Tibetan alpine meadow.
• 2. Experimental Design:
  • A long-term (13-year) field experiment with a gradient of nutrient additions was used.
  • Treatments included a control (no fertilizer) and various levels of NP fertilizer ((NH₄)₂HPO₄), such as NP30, NP90, and NP120 (i.e., 30, 90, and 120 g m⁻²).
• 3. Data Collection and Sampling:
  • Soil Physicochemical Properties: Soil samples were analyzed for labile C, mineral N, available P, and pH.
  • Soil Biodiversity Assessment:
    • DNA Extraction and Sequencing: Soil microbial DNA was extracted. High-throughput sequencing of 16S rRNA genes (for bacteria) and ITS regions (for fungi) was performed to characterize community composition and calculate Shannon diversity indices.
    • Nematode Extraction and Identification: Nematodes were extracted from soil samples via centrifugation-flotation. They were identified to genus level and classified into trophic groups (bacterivores, fungivores, plant parasites, omnivores, predators) under a microscope. Diversity indices were calculated.
  • Ecosystem Function Measurement: Fourteen ecosystem functions related to C, N, and P cycling were measured. These included:
    • C & Nutrient Stocks: Total soil C and N, microbial biomass C, N, and P.
    • C & Nutrient Turnover: Soil basal respiration (microbial activity), degradation rates of various carbon compounds (sugar, chitin, lignin, polymer), and P mineralization.
    • Ecosystem Stability: Aggregate stability and resistance to plant-parasitic nematodes.
  • Ecosystem Multifunctionality (EMF) Calculation: EMF was quantified using the average approach (averaging the standardized values of all measured functions) and the multi-threshold approach (counting the number of functions passing certain percentage thresholds of their maximum observed values).
• 4. Statistical Analysis: Linear models and correlation analyses were used to relate changes in soil properties (especially pH) to changes in biodiversity and EMF.

Protocol 2: Legume Intercropping Experiment in a Degraded Ecosystem

This protocol tests how alleviating nutrient limitation through legume intercropping triggers trophic cascades in the soil micro-food web [93].

• 1. Site and Soil: A field experiment was conducted in a degraded agroecosystem with nutrient-poor purple soil.
• 2. Experimental Treatments:
  • Camellia oleifera monoculture (CK, control).
  • C. oleifera intercropped with peanut (Arachis hypogaea, CP).
  • C. oleifera intercropped with Senna tora (CS).
• 3. Data Collection and Analysis:
  • Soil Nutrient Limitation Assessment: Soil nutrient limitation, particularly for phosphorus (P), was assessed as a foundational premise of the experiment.
  • Soil Biota Community Analysis:
    • Microbial Community: Soil DNA was sequenced to analyze the diversity and composition of bacterial communities.
    • Protists and Nematodes: Protist consumer diversity and nematode communities were characterized, with nematodes identified to trophic groups.
  • Metabolomic Analysis: Untargeted soil metabolomic analysis was conducted to profile the metabolic pathways and confirm nutrient limitation states.
  • Statistical Modeling: Partial Least Squares Path Modeling (PLS-PM) was applied to test the causal pathways (trophic cascades) linking legume intercropping, soil P limitation, bacterial and protist diversity, nematode abundance, and soil P metabolism functions.
• 4. Key Measured Variables: Biodiversity of bacteria, protists, and nematodes; abundance of omnivore-predator nematodes; soil P-related metabolic functions.

---

Visualizing Mechanisms and Workflows

The following diagrams illustrate the key mechanistic pathways and experimental workflows discussed in this article.

Diagram 1: Mechanisms Linking Soil Fertility and the Biodiversity-Ecosystem Function (BEF) Relationship. In nutrient-poor soils, high complementarity and strong trophic cascades make biodiversity critical for function. In fertile soils, species dominance and acidification can weaken the BEF relationship.

Diagram 2: Trophic Cascade Triggered by Legume Intercropping in a Nutrient-Poor Soil. This pathway shows how a management practice creates a bottom-up cascade that increases diversity at multiple trophic levels, which in turn top-down regulates soil function.

---

The Scientist's Toolkit: Key Research Reagents & Materials

This table details essential reagents, materials, and tools for conducting research on soil biodiversity and ecosystem function.

Item / Category	Primary Function in Research	Application Example / Notes
Molecular Biology Reagents
DNA Extraction Kits (e.g., MoBio PowerSoil)	Extract high-quality genomic DNA from complex soil matrices.	Fundamental for downstream sequencing of bacterial, fungal, and other microbial communities [94] [95].
16S rRNA & ITS Primers	Amplify specific gene regions for identifying bacteria/archaea (16S) and fungi (ITS) via PCR.	Enables taxonomic classification and diversity estimation using high-throughput sequencing.
Field & Lab Equipment
Soil Corers	Collect standardized, minimally disturbed soil samples at specific depths.	Essential for consistent spatial and temporal sampling [94] [93].
Centrifuges	Separate soil particles, microorganisms, or nematodes from suspensions.	Used in nematode extraction protocols and sample processing [94].
High-Throughput Sequencer (e.g., Illumina)	Determine the sequence of amplified DNA fragments from many samples in parallel.	Generates the raw data for characterizing soil microbial and fungal diversity [94] [95].
Analytical Tools & Software
QIIME 2, MOTHUR	Process and analyze raw sequencing data: quality filtering, clustering sequences into OTUs/ASVs, taxonomic assignment.	Standard bioinformatics pipelines for microbial ecology [94].
Partial Least Squares Path Modeling (PLS-PM)	Test complex causal hypotheses and path relationships (e.g., trophic cascades) from observational data.	Used to model the effect of legume intercropping on the soil micro-food web [93].
Chemical Assays
Microplate Enzymatic Assays	Measure the potential activity of extracellular enzymes (e.g., for C, N, P cycling) in soil samples.	Key for linking microbial community to ecosystem function [94].
Chloroform Fumigation	Determine soil microbial biomass carbon, nitrogen, and phosphorus by measuring lysed cell content.	Standard method for estimating the size of the soil microbial pool.

---

The comparative evidence robustly validates that soil biodiversity is not merely a passive consequence of fertility but an active driver of ecosystem function, whose role becomes critically indispensable in low-resource environments. In nutrient-poor soils, the loss of species can directly and disproportionately impair essential processes like nutrient cycling and carbon sequestration due to low functional redundancy and strong trophic interactions [6] [93]. Conversely, the common anthropogenic practice of nutrient enrichment, while seemingly boosting short-term fertility, can erode this biodiversity through mechanisms like soil acidification, thereby weakening the natural resilience of the soil ecosystem [94].

For researchers, this underscores the necessity of incorporating biodiversity metrics into soil health assessments and the development of sustainable management strategies. In agricultural and restoration contexts, practices like legume intercropping that manipulate soil nutrient stoichiometry and foster diverse soil food webs offer promising pathways to enhance ecosystem functioning without compromising soil biological integrity [93]. Future research must continue to unravel the complex feedbacks between specific nutrient limitations, distinct trophic levels, and ecosystem multifunctionality to better predict and manage soil health under global change.

Biodiversity-ecosystem functioning (BEF) research seeks to unravel the complex relationships between biological diversity and the stability and efficiency of ecological processes. In soil ecosystems, this relationship is paramount, underpinning critical functions from nutrient cycling to carbon sequestration. A central thesis in this field posits that specific taxa exert a disproportionate influence on these functions, their effects shaped by unique physiological traits and interspecific interactions. This guide provides a comparative examination of three such taxa—mycorrhizal fungi, nematodes, and earthworms—evaluating their validated roles through the lens of experimental data and mechanistic studies. By synthesizing quantitative evidence and standard methodologies, this analysis aims to inform research strategies and reagent selection for testing BEF hypotheses in subterranean environments.

Comparative Impact of Key Soil Taxa

The disproportionate impact of mycorrhizal fungi, nematodes, and earthworms is demonstrated through their contributions to various soil processes. The table below provides a comparative summary of their validated roles and effect sizes based on experimental data.

Table 1: Comparative Quantitative Impacts of Mycorrhizal Fungi, Nematodes, and Earthworms on Soil Processes

Taxon	Key Ecosystem Functions	Quantified Impact (Effect Size/Direction)	Key Experimental Context
Mycorrhizal Fungi	Soil aggregation, macroaggregate formation	+24% overall soil aggregation [96]; fungi contribute more strongly to macroaggregates than bacteria or animals [96]	Global meta-analysis of 183 studies [96]
	Nutrient uptake (Phosphorus, Nitrogen)	Regulates ~80% of plant nitrogen and 50% of soil carbon [97]; solubilizes bound phosphorus [97]	Synthesis of field and pot studies [97]
	Plant growth and stress tolerance	Improves host plant water use efficiency and detoxification under heavy metal stress [97]	Inoculation experiments with tree seedlings [97]
Earthworms	Soil aggregation (varies by functional group)	Endogeic species show significant positive effect; epigeic and anecic effects are neutral or less consistent [96]	Global meta-analysis [96]
	Nutrient cycling (Nitrogen)	Enhance denitrification through increased organic compounds [98]	Literature review and experimental studies [98]
	Plant growth promotion	Effects are species-specific; can be positive, negative, or neutral [99] [98]	Factorial pot experiments (e.g., with turf grasses) [99]
Nematodes	Nutrient cycling, nitrogen mineralization	Free-living nematodes contribute significantly to nitrogen mineralisation [18]	Field studies in Atlantic Forest landscapes [18]

Detailed Experimental Protocols and Methodologies

Validating the disproportionate impact of soil taxa requires robust, replicable experimental designs. Below are detailed methodologies for key experiments cited in this guide.

Table 2: Key Experimental Protocols for Investigating Soil Taxa Functionality

Experiment Focus	Core Protocol Summary	Key Reagents & Measurements	Supporting Citations
Soil Aggregation Meta-Analysis	Global compilation of 183 studies (345 trials) measuring soil biota effects on aggregate stability. Data grouped by taxon (animals, bacteria, fungi) and aggregate size (macro vs. micro).	Measurement: Mean weight diameter (MWD) or water-stable aggregates (WSA).Analysis: Random effects meta-analysis to calculate overall effect sizes and confidence intervals [96].	[96]
AMF & Earthworm Interactions	Factorial pot experiment: Five turfgrass species inoculated with AMF (Glomus mosseae), earthworms (Pheretima tschiliensis), both, or none.	Inocula: AMF spores, hyphae, colonized root fragments; endogeic earthworms.Plant Metrics: Height, biomass.Soil Metrics: Physical/chemical properties. Data analyzed via PCA and weighted subordinate function for comprehensive evaluation [99].	[99]
Soil Biodiversity-Multifunctionality	Field sampling in natural and degraded Atlantic Forest. Soil biodiversity (microbes, fungi, nematodes, invertebrates) and ecosystem functions (PP, ES, NC) measured across seasons.	Organism Extraction: Wet sieving and decanting for mycorrhizae [18]; centrifugal-flotation for nematodes [18].Statistical Analysis: Pearson’s correlation; Linear Mixed-Effects Models (LMMs) to disentangle direct/indirect effects [18].	[18]
Earthworm-Microbiota Interactions	Review synthesizing studies on earthworm effects on microbial community structure (via PLFA, sequencing) and functional genes (e.g., for N cycling).	Methods: Phospholipid fatty acid (PLFA) analysis, DNA sequencing, qPCR for functional genes.Focus: Compare effects across drilosphere components (gut, casts, burrows) and earthworm ecological categories [98].	[98]

Research Reagent Solutions for Key Taxa Studies

Selecting appropriate reagents and materials is fundamental for the experimental validation of soil taxa functions.

Table 3: Essential Research Reagents and Materials for Soil Taxa Functional Studies

Reagent/Material	Function in Experiment	Specific Application Example
AMF Inoculum (e.g., Glomus mosseae)	Establish symbiotic relationship with plant roots to test effects on nutrient uptake, growth, and stress tolerance.	Inoculum consists of spores (~28 per gram), hyphae, and colonized root fragments [99].
Earthworm Species (e.g., Pheretima tschiliensis)	Test the role of soil engineers from specific functional groups (endogeic) on soil structure and nutrient dynamics.	Endogeic earthworms are used to study their role as major agents of soil aggregation and organic matter stabilization [99].
Extraction Solutions (e.g., Sterilized Water)	Safely handle and prepare biological specimens without introducing contaminants or causing harm.	Used for washing earthworms before their addition to experimental pots [99].
Molecular Reagents (for PLFA, DNA Sequencing)	Quantify and characterize soil microbial community structure, abundance, and functional potential.	Used to analyze how earthworms modify soil microbiota in their gut, casts, and the surrounding drilosphere [98].
Enzyme Assay Kits	Measure the activity of microbial enzymes (e.g., proteinases, phosphatases) critical to nutrient cycling.	To quantify the extracellular enzymes produced by EMF that break down organic nitrogen or solubilize phosphorus [97].

Mechanisms and Pathways of Influence

The disproportionate impact of these taxa is mediated through distinct but sometimes interconnected mechanistic pathways.

Diagram 1: Key pathways showing the direct (&) and indirect (§) mechanisms through which mycorrhizal fungi, earthworms, and nematodes influence the soil system and plant health. The 'Sleeping Beauty Paradox' describes how earthworm mucus awakens dormant microbes [98].

Integrated Experimental Workflow

Testing BEF relationships for specific taxa requires an integrated approach that spans from field sampling to molecular analysis.

Diagram 2: A generalized workflow for designing and executing experiments to test biodiversity-ecosystem functioning relationships, integrating field and lab methods [18] [99].

The experimental data and mechanistic pathways synthesized in this guide validate the thesis that mycorrhizal fungi, nematodes, and earthworms are not merely residents of the soil matrix but are disproportionate drivers of ecosystem functioning. The evidence confirms that their impacts are quantifiable, taxon-specific, and often mediated through complex interactions with the soil microbiota and plant hosts. Mycorrhizal fungi, particularly ectomycorrhizal types, stand out for their dual role in soil aggregation and nutrient cycling in forest systems [96] [97]. Earthworms, especially endogeic species, function as physical ecosystem engineers but their effects on plant growth are highly context-dependent [96] [99]. Nematodes contribute significantly to nutrient mineralization, linking microbial activity to plant-available nutrients [18].

Future research should leverage emerging technologies, such as AI-enhanced microscopy and high-throughput molecular tools, to scale up these mechanistic understandings from controlled experiments to landscape-level monitoring [56]. Furthermore, adopting a holistic "plant-soil system" evaluation framework, as demonstrated in recent multi-factorial experiments [99], is crucial for accurately predicting the outcomes of biodiversity changes on the myriad functions upon which terrestrial life depends.

Soil biodiversity represents one of the most complex and abundant assemblages of life on Earth, encompassing microorganisms (bacteria, fungi, archaea, protists) and micro-, meso-, and macrofauna (nematodes, collembola, earthworms, and larger organisms) [6] [15]. The relationship between this biodiversity and the functions it performs—the soil Biodiversity-Ecosystem Functioning (BEF) relationship—has emerged as a critical frontier in ecological research. Historically, soil systems were thought to exhibit high functional redundancy, where the loss of specific taxa would have minimal impact on ecosystem processes because other organisms could perform the same roles [6] [15]. However, a paradigm shift is underway, driven by growing experimental and observational evidence. Recent research suggests that losses in soil microbial diversity can result in proportional or even exponential declines in certain soil functions, particularly for specialized processes like denitrification or complex ones like organic matter decomposition that require coordinated action from multiple organisms [6] [15].

This review synthesizes evidence from global studies to examine the consistency and variability of soil BEF relationships. We compare findings across biomes, taxonomic groups, and spatial scales, providing a comprehensive analysis of the experimental data, methodological approaches, and contextual factors that determine how soil biodiversity sustains ecosystem functions crucial for human wellbeing and planetary health.

Global Patterns and Key Research Gaps in Soil BEF Research

Documented Consistency: Positive Biodiversity-Multifunctionality Relationships

Strong evidence from multiple studies confirms that soil biodiversity is significantly and positively associated with multiple ecosystem functions (multifunctionality) across diverse biomes. A landmark global study analyzing soils from across the world demonstrated that soil biodiversity (including bacteria, fungi, protists, and invertebrates) is positively correlated with nutrient cycling, decomposition, plant production, and pathogen control [100]. This relationship held across broad environmental gradients, suggesting a general pattern where diverse soil communities support a wider array of functions simultaneously. Similarly, research in grassland ecosystems has confirmed that the richness of multiple soil organism groups promotes ecosystem multifunctionality, though the strength of these relationships varies with environmental context [3] [100].

The mechanisms underlying these consistent positive relationships primarily involve functional complementarity and asynchronous responses among soil taxa. Different soil organisms perform distinct functions at different times or under varying conditions, creating stability in ecosystem processes. Experimental evidence shows that diverse soil microbial communities exhibit asynchrony, where different fungi and bacteria promote different ecosystem functions at different times, thereby stabilizing functioning over time [101]. This asynchrony provides ecological insurance, ensuring that some taxa capable of supporting critical functions are present despite environmental fluctuations [101].

Documented Variability: Context-Dependent Nature of Soil BEF

Despite these consistent patterns, soil BEF relationships exhibit substantial variability influenced by environmental factors, spatial scales, and ecosystem characteristics. The strength of biodiversity-multifunctionality relationships displays marked context dependency [100]. For instance, soil biodiversity has proven particularly critical for supporting function in drylands and poorly developed soils with limited nutrient stocks [6] [15]. In these systems, the daily contribution of diverse soil microbiota to decompose and depolymerize litter represents a essential nutrient input, making these ecosystems more sensitive to biodiversity loss [15].

Spatial scale further influences observed BEF relationships. Experimental work with bacterial metacommunities found positive diversity-productivity relationships at regional scales driven by niche differentiation, but no correlation at local scales where diversity was maintained primarily through immigration [102]. This highlights how mechanisms maintaining diversity (niche partitioning vs. mass effects) operate differently across scales, generating variable BEF patterns [102].

Table 1: Key Drivers of Variability in Soil BEF Relationships

Factor	Effect on BEF Relationship	Supporting Evidence
Environmental Context	Stronger in nutrient-poor soils; weaker in fertile soils	[3] [6] [15]
Spatial Scale	Positive at regional scales; weak or absent at local scales	[102]
Soil Community History	Long-term assembly enhances functioning; recent disturbance weakens it	[103]
Trophic Complexity	Multitrophic diversity stronger predictor than single-group diversity	[100]
Temporal Scale	Relationships strengthen over time as communities assemble	[101] [103]

Comparative Analysis of Key Experimental Findings

Global Observational Studies vs. Controlled Experiments

Research approaches to soil BEF relationships broadly fall into two categories: large-scale observational studies and controlled experiments, each yielding complementary insights. Global field surveys provide unparalleled evidence for naturally occurring patterns but face challenges in establishing causality and controlling confounding variables [10] [100]. These studies have revealed that the most dominant bacterial phylotypes in soils worldwide are consistently linked to multiple ecosystem functions, suggesting certain keystone taxa play disproportionate roles in maintaining functioning [100].

Controlled experiments, in contrast, enable researchers to isolate biodiversity effects from environmental covariates and establish mechanistic links. A pioneering microcosm experiment that manipulated soil microbial diversity through size-based filtering demonstrated that microbial diversity enhanced the temporal stability of multiple ecosystem functions, including plant biomass production, plant diversity, litter decomposition, and soil carbon assimilation [101]. The stabilizing effect was particularly strong in treatments with reduced microbial richness where over 50% of microbial taxa were lost, highlighting the vulnerability of simplified systems [101].

The Plant Diversity-Soil BEF Controversy

The relationship between plant diversity and soil functioning represents an area of ongoing debate and variable findings. Some studies report indirect effects of plant diversity on soil multifunctionality mediated through soil biodiversity [100], while others find minimal plant diversity effects compared to other factors. A recent Ecotron experiment found that soil community history rather than plant species richness (1-6 species) enhanced belowground multitrophic functioning [103]. After four months, soil fauna-driven ecosystem functioning was promoted in communities with shared soil history regardless of plant diversity level, suggesting that soil community assembly and biomass accumulation may outweigh short-term plant diversity effects [103].

Table 2: Comparison of Major Soil BEF Study Approaches and Findings

Study Type	Key Findings	Limitations	Representative Evidence
Global Field Survey	Positive soil biodiversity-multifunctionality relationship across biomes; context-dependent effects	Correlation not causation; environmental co-variation	[100]
Microcosm Experiment	Diversity stabilizes ecosystem functioning; asynchrony mechanisms	Simplified communities; limited environmental realism	[101]
Grassland Biodiversity Experiment	Soil biodiversity decline reduces some functions; fertility mediates effects	Specific to grassland systems	[3]
Ecotron Experiment	Soil history outweighs plant diversity; body mass distributions matter	Short duration; controlled conditions	[103]

Methodological Approaches and Experimental Protocols

Standardized Experimental Designs for Soil BEF Research

Experimental research on soil BEF relationships has converged on several methodological approaches that enable comparability across studies while allowing investigation of specific mechanisms.

Soil Community Dilution or Filtering Approach: A widely employed method involves creating a gradient of soil biodiversity through serial dilution or size-based filtering of soil organisms [3] [101]. In a typical protocol, soil samples are suspended in sterile solution and passed through filters of decreasing pore sizes (e.g., 5000 μm, 100 μm, 25 μm, and finally sterile through 0 μm) to systematically exclude larger organisms while allowing smaller ones to pass through [101]. This approach successfully creates a monotonically declining gradient in fungal and bacterial richness, with the sterile treatment typically resulting in approximately 60% loss of bacterial and 55% loss of fungal richness compared to the highest biodiversity treatment [101].

Common Garden and Mesocosm Experiments: These controlled environment studies allow researchers to test the effects of manipulated soil communities on ecosystem functions while holding constant other environmental variables. Standard protocols involve: (1) collecting soils from different field sites or under different management histories; (2) optionally, sterilizing soils (via gamma irradiation or autoclaving) to create a microbial-free baseline; (3) inoculating sterilized soils with different microbial communities; and (4) measuring ecosystem functions over time [3] [103]. These experiments typically measure functions including plant productivity, litter decomposition, soil respiration, nutrient cycling rates, and pathogen suppression [3] [101] [100].

DNA-Based Molecular Techniques: Standardized molecular approaches now enable comprehensive characterization of soil biodiversity. Typical workflows include: (1) DNA extraction from soil samples using commercial kits; (2) PCR amplification of taxonomic marker genes (e.g., 16S rRNA for bacteria, ITS for fungi, 18S rRNA for protists and microfauna); (3) high-throughput sequencing; and (4) bioinformatic processing using pipelines like UPARSE or QIIME to assign sequences to operational taxonomic units (OTUs) or amplicon sequence variants (ASVs) [100]. These methods allow researchers to characterize microbial and microfaunal communities with unprecedented depth, enabling tests of how specific taxonomic groups correlate with ecosystem functioning.

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Key Research Reagents and Materials for Soil BEF Studies

Reagent/Material	Function in Soil BEF Research	Example Application
DNA Extraction Kits	Isolation of genomic DNA from complex soil matrices	Characterizing bacterial, fungal, and protist communities [100]
PCR Primers	Amplification of taxonomic marker genes	16S rRNA genes (bacteria), ITS regions (fungi), 18S rRNA (protists) [100]
Sterile Soil Substrates	Create microbial-free baseline for inoculation experiments	Assessing function of specific microbial inocula [3] [101]
Biolog Plates	Assess community metabolic potential	Measuring carbon substrate utilization profiles [102]
Tetrazolium Dyes	Indicator of metabolic activity	Colorimetric measure of substrate metabolism [102]
Size-Specific Filters	Create biodiversity gradients via size exclusion	Separating soil organisms by body size [101]
Enzyme Assay Reagents	Measure extracellular enzyme activities	Indicators of nutrient cycling processes [100]

Conceptual Framework and Visual Synthesis

The relationship between soil biodiversity and ecosystem functioning operates through multiple interconnected pathways that operate across temporal and spatial scales. The following conceptual diagram synthesizes these key relationships and mechanisms based on current evidence:

Conceptual Framework of Soil Biodiversity-Ecosystem Functioning Relationships

The diagram illustrates how soil biodiversity influences ecosystem functioning through multiple mechanisms (complementarity, asynchrony, trophic interactions, and ecosystem engineering), with these relationships moderated by contextual factors including spatial scale, soil properties, community history, and climate. These pathways collectively enhance multiple ecosystem functions and their stability over time.

The synthesis of current evidence reveals both consistent patterns and important variability in soil BEF relationships. Consistently, soil biodiversity enhances ecosystem multifunctionality and stability across diverse systems, primarily through mechanisms of functional complementarity and asynchronous responses among taxa. However, the strength of these relationships varies substantially with environmental context, spatial scale, soil properties, and community history.

Critical knowledge gaps remain in our understanding of soil BEF relationships. The spatial and environmental representativeness of current studies is limited, with significant blind spots in tropical regions, drylands, and specific soil types [10]. Furthermore, a striking disparity exists in data availability across taxonomic groups, with bacteria and fungi being considerably better studied than microfauna and mesofauna [10]. Perhaps most importantly, temporally explicit data is nearly absent from soil macroecological studies, limiting our understanding of how BEF relationships evolve as communities assemble and respond to global changes [10] [103].

Future research should prioritize: (1) expanding global coverage of soil biodiversity monitoring to underrepresented regions and ecosystems; (2) integrating multiple trophic levels to understand how food web structure influences BEF relationships; (3) establishing long-term observational and experimental networks to capture temporal dynamics; and (4) translating mechanistic understanding into management practices that conserve and enhance soil biodiversity for ecosystem service provision. As soils face increasing pressures from global change, understanding the conditions under which soil biodiversity sustains ecosystem functioning becomes not merely an academic pursuit but an imperative for ecosystem conservation and human wellbeing.

Understanding the relationships between soil biodiversity and ecosystem functioning (BEF) is a central goal in soil ecology, with profound implications for agriculture, climate change mitigation, and ecosystem restoration [15]. This field faces a unique methodological challenge: bridging the gap between controlled, small-scale experimental validation and the complex reality of landscape-level ecosystem dynamics. Research conducted in microcosms allows for precise manipulation and causal inference, while landscape-level studies capture emergent properties and real-world contextual factors that determine how biodiversity influences ecosystem multifunctionality—the simultaneous performance of multiple ecosystem functions [18] [5]. This review systematically compares these complementary research approaches, their experimental protocols, and their respective contributions to advancing BEF theory and application. By synthesizing findings from controlled manipulations to large-scale observational studies, we provide researchers with a comprehensive framework for selecting appropriate methodologies across the spectrum of scientific inquiry in soil ecology.

Comparative Analysis of Experimental Approaches

Soil biodiversity research employs a hierarchy of experimental approaches, each with distinct advantages, limitations, and appropriate applications. The following table summarizes the key characteristics of the primary methodologies used across the scale spectrum.

Table 1: Comparison of Experimental Approaches in Soil Biodiversity Research

Experimental Approach	Spatial Scale	Key Control Level	Primary Data Types	Ecological Context	Major Advantages	Principal Limitations
Microcosm Experiments	Laboratory (mg to kg of soil)	High	Taxon abundance, process rates	Highly simplified	Established causation, mechanistic insights, high replication	Limited realism, artificial conditions
Mesocosm Studies	Intermediate (field enclosures, lysimeters)	Moderate	Community composition, multiple functions	Semi-natural	Balance of control and realism, multi-trophic studies	Constrained spatial/temporal scales
Landscape Observational Studies	Field (plots to regions)	Low (natural gradients)	Biodiversity indices, multifunctionality metrics, environmental covariates	Natural ecosystems	Real-world relevance, contextual factors, spatial patterns	Correlation rather than causation, confounding factors
Land-Use Change Studies	Landscape	Natural experiments	Biodiversity, nutrient cycling, primary production, stability [18]	Anthropogenically modified	Direct policy relevance, restoration insights	Historical contingencies, multiple stressors

Detailed Methodological Protocols

Microcosm and Mesocosm Experimental Designs

Controlled experiments using microcosms and mesocosms form the foundation for establishing causal relationships in soil BEF research. These approaches typically involve manipulating biodiversity levels while holding environmental conditions constant.

Protocol for Biodiversity Manipulation Studies:

• Soil Inoculum Preparation: Collect field soil and subject it to gradient sterilization or sieving to create biodiversity gradients, or use synthetic communities (SynComs) assembled from cultured isolates [15].
• Experimental Units: Use sterile containers with controlled drainage, typically containing 100g to 1kg of sterilized background soil.
• Biodiversity Treatments: Inoculate with soil suspensions representing different biodiversity levels (e.g., from 1 to 10+ taxonomic groups).
• Incubation Conditions: Maintain under controlled temperature and moisture conditions relevant to the ecosystem of interest.
• Function Measurements: Quantify process rates including organic matter decomposition (mass loss of standardized litter), nutrient mineralization (soil extractions), and microbial respiration (gas chromatography) [5].

Landscape-Level Observational Approaches

Landscape-scale studies examine BEF relationships across natural environmental gradients or land-use types, capturing the complexity of real-world ecosystems.

Protocol for Atlantic Forest Landscape Study [18]:

• Site Selection: Establish study sites across natural forests, pastures, and deforested areas to represent a land-use intensity gradient.
• Soil Sampling: Collect soil cores during both rainy (May-August) and dry seasons (November-February) on a monthly schedule across multiple years to capture temporal variation.
• Biodiversity Quantification:
  • Extract and identify soil fauna (insects, arachnids, myriapods) using standardized methods [18].
  • Quantify microbial communities through DNA sequencing targeting 16S rRNA (bacteria) and ITS (fungi) regions.
  • Assess mycorrhizal fungi using wet sieving and decanting techniques [18].
• Ecosystem Function Assessment: Measure multiple simultaneous functions:
  • Primary production (plant biomass)
  • Nutrient cycling (N, P availability)
  • Ecosystem stability (resistance to disturbance)
• Statistical Analysis: Apply Pearson's correlation analysis to identify relationships between biodiversity and functions, followed by linear mixed-effects models (LMMs) to account for spatial and temporal autocorrelation [18].

Protocol for Mine Reclamation Chronosequence Study [104]:

• Space-for-Time Substitution: Establish study plots across a 30-year restoration gradient in an opencast coal mine, including naturally restored grasslands and forests, artificially reclaimed vegetation, and unreclaimed bare land.
• Multidimensional Biodiversity Assessment:
  • Characterize α-diversity (within-sample), β-diversity (between-sample), and network complexity for bacteria, archaea, fungi, and eukaryotic communities.
  • Construct co-occurrence networks to quantify topological features.
• Multifunctionality Quantification: Calculate soil multifunctionality (SMF) index as the average Z-score of 19 standardized variables representing soil nutrients, enzyme activities, and microbial biomass.
• Statistical Modeling: Apply Random Forest and multiple regression analyses to identify dominant biodiversity drivers of SMF recovery [104].

Key Research Findings Across Scales

Consistent Patterns Emerging from Multiple Approaches

Despite methodological differences, research across scales reveals several consistent patterns regarding soil BEF relationships:

• Positive Biodiversity-Multifunctionality Correlation: Studies from microcosms to landscapes consistently report positive relationships between soil biodiversity and ecosystem multifunctionality [18] [15] [104]. In the Atlantic Forest landscape study, soil biodiversity (including insects, arachnids, myriapods, fungi, nematodes, and bacteria) showed highly significant correlations (p < 0.01) with primary production, ecosystem stability, and nutrient cycling [18].

• Land-Use Impacts: Both observational studies and land-use experiments demonstrate that agricultural intensification and deforestation reduce soil biodiversity and ecosystem functioning. Synthetic N fertilization negatively impacts arbuscular mycorrhizal fungal (AMF) and faunal diversity, while tillage reduces soil faunal and bacterial diversity [105].

• Functional Redundancy Reassessment: Contemporary research challenges the historical view of high functional redundancy in soil communities. Diversity losses can result in proportional or exponential reductions in specialized processes like denitrification and organic matter decomposition [15].

Scale-Dependent Insights

Different methodological approaches reveal unique insights that are often scale-dependent:

• Microcosm Studies demonstrate mechanistic causality, showing that microbial diversity losses directly reduce specific process rates like litter decomposition and nitrogen mineralization [15].

• Landscape Observations reveal contextual importance, finding that soil biodiversity is particularly critical for supporting function in nutrient-poor soils where daily microbial contributions to nutrient cycling are essential [15].

• Restoration Studies show that different reclamation approaches yield varying outcomes, with mixed coniferous and broad-leaved forests exhibiting greater soil multifunctionality than pure forests or grasslands in mine reclamation contexts [104].

Table 2: Quantitative Findings from Major Soil BEF Studies

Study Type	Location/System	Key Quantitative Findings	Statistical Approaches
Landscape observational [18]	Atlantic Forest, Brazil	Natural ecosystems showed significantly higher soil biodiversity and function than pastures and deforested areas; Soil biodiversity mediated nutrient cycling under seasonal variability	Pearson correlation, Linear mixed-effects models
Mine reclamation chronosequence [104]	Antaibao opencast coal mine, China	Mixed forests had highest SMF; Bacterial and fungal β-diversity were dominant factors affecting SMF recovery; Network complexity of bacteria, archaea, and eukaryotes significantly influenced SMF	Random Forest, Multiple regression analysis
Agricultural intensification meta-analysis [105]	Global agricultural systems	Low synthetic N inputs (<150 kg N/ha/year) increased bacterial diversity; Organic inputs and >5 year application duration enhanced bacterial biodiversity; Tillage negatively impacted faunal and bacterial diversity	Meta-analysis of 85 studies

Conceptual Framework and Signaling Pathways

The relationship between methodological approaches in soil BEF research can be visualized as a hierarchical framework where findings at smaller scales inform understanding at larger scales, and landscape observations generate hypotheses for mechanistic testing.

Diagram 1: Multiscale Research Integration

The complex interactions between soil organisms and ecosystem functions can be conceptualized as a network of signaling pathways and feedback mechanisms that operate across spatial and temporal scales.

Diagram 2: Soil Biodiversity-Function Signaling Pathways

The Scientist's Toolkit: Essential Research Reagent Solutions

Contemporary soil BEF research requires specialized reagents, analytical tools, and computational approaches to effectively characterize biodiversity and ecosystem processes across scales.

Table 3: Essential Research Solutions for Soil BEF Studies

Tool Category	Specific Solutions	Research Application	Technical Considerations
Molecular Analysis	16S rRNA gene sequencing (V3-V4 region); ITS sequencing for fungi	Taxonomic characterization of bacterial and fungal communities	Primer selection critical for taxonomic resolution; standardization enables cross-study comparisons
Bioinformatics	QIIME2, mothur for amplicon analysis; Network analysis pipelines (e.g., SPIEC-EASI)	Processing sequencing data; constructing co-occurrence networks	Parameter selection significantly impacts network topology; null model correction essential
Soil Fauna Extraction	Berlese-Tullgren funnels; Wet sieving and decanting [18]	Standardized extraction of microarthropods; mycorrhizal fungi quantification	Extraction efficiency varies with soil type; morphological identification requires specialist expertise
Function Assays	Microplate enzyme assays; Substrate-induced respiration; Litter decomposition bags	Quantification of process rates; functional potential assessments	Standardized substrates (e.g., MUB-linked substrates for enzymes) enable cross-study comparison
Isotopic Tracers	13C-labeled substrates; 15N tracing experiments	Tracking carbon and nitrogen fluxes through food webs	Requires specialized mass spectrometry facilities; powerful for process rate quantification
Statistical Modeling	Linear mixed-effects models (LMMs); Random Forest; Structural equation modeling (SEM)	Accounting for spatial/temporal autocorrelation; identifying key predictors	Model selection critical for robust inference; cross-validation essential for machine learning approaches
Geospatial Analysis	GIS integration; Remote sensing (NDVI, EVI); Digital soil mapping [106]	Landscape-level pattern analysis; spatial prediction of soil properties	Resolution limitations for microbial patterns; useful for scaling plot-level measurements

The integration of microcosm validations with landscape-level observations represents the most promising path forward for soil BEF research. While microcosm studies provide essential mechanistic understanding and causal evidence, landscape approaches capture the contextual factors and emergent properties that determine real-world ecosystem functioning [18] [104]. Future research should prioritize designs that explicitly bridge these scales, such as nested studies that combine broad-scale surveys with intensive mechanistic investigations at subset locations. The increasing sophistication of digital soil mapping [106] [107], molecular tools, and computational approaches [108] offers unprecedented opportunities to overcome traditional limitations of scale and complexity. By strategically combining these approaches, researchers can generate robust, predictive understanding of how soil biodiversity supports ecosystem functions across changing landscapes—knowledge that is essential for addressing pressing global challenges including climate change, food security, and ecosystem restoration.

Conclusion

The body of evidence unequivocally confirms that soil biodiversity is a fundamental driver of ecosystem functioning, but this relationship is profoundly complex and context-dependent. Key takeaways reveal that while higher soil biodiversity generally promotes multifunctionality, its effects are modulated by environmental factors like soil fertility and land-use history. Methodological advances now allow us to move beyond simple correlations to decipher mechanistic links. However, significant challenges remain, including geographic research blind spots and the need to better understand the trade-offs between functional diversity and redundancy. Future research must prioritize interdisciplinary approaches that integrate soil ecology with global change science, empower underrepresented regions, and translate findings into practical strategies for soil conservation and sustainable management. This knowledge is not just academically vital but is essential for informing policies that ensure the health of our planet's soils and the countless services they provide.

References

• 1. Soil biodiversity - Encyclopedia of the Environment [https://www.encyclopedie-environnement.org/en/soil/soil-biodiversity/]
• 2. Soil Biodiversity: State‐of‐the‐Art and Possible ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8246784/]
• 3. Soil biodiversity and ecosystem functions in grasslands [https://www.sciencedirect.com/science/article/abs/pii/S0038071725002822]
• 4. Scaling up biodiversity–ecosystem functioning relationships [https://pmc.ncbi.nlm.nih.gov/articles/PMC7944106/]
• 5. Soil biodiversity promotes key ecosystem functions by its ... [https://soil-organisms.org/index.php/SO/article/view/450]
• 6. Soil biodiversity and function under global change - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC12135966/]
• 7. Functional diversity of soil macrofauna may contribute to ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12202550/]
• 8. Edaphobase 2.0: Advanced international data warehouse ... [https://www.sciencedirect.com/science/article/pii/S0929139324004414]
• 9. Soil biodiversity management [https://www.biodiversa.eu/2025/09/25/soil-biodiversity-management/]
• 10. Blind spots in global soil biodiversity and ecosystem ... [https://www.nature.com/articles/s41467-020-17688-2]
• 11. Biodiversity and ecosystem functioning in soil - Ambio [https://link.springer.com/article/10.1007/s13280-021-01507-z]
• 12. Microbial keystone taxa drive succession of plant residue ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10119086/]
• 13. Microbial diversity and keystone species drive soil nutrient ... [https://www.sciencedirect.com/science/article/abs/pii/S0013935124006194]
• 14. Global Soil Biodiversity Initiative [https://www.globalsoilbiodiversity.org/]
• 15. Soil biodiversity and function under global change [https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3003093]
• 16. Rare Species-Driven Diversity–Ecosystem Multifunctionality ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9239226/]
• 17. Plant diversity maintains multiple soil functions in future ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6296783/]
• 18. Soil biodiversity drives ecosystem multifunctionality across ... [https://www.sciencedirect.com/science/article/abs/pii/S0031405625000605]
• 19. Effects of management practices on the ecosystem-service ... [https://www.nature.com/articles/s41467-024-48049-y]
• 20. Soil biodiversity and soil community composition determine ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3986181/]
• 21. Nitrogen and CO2 enrichment interact to decrease ... [https://www.nature.com/articles/s41467-025-62691-0]
• 22. Reconsidering functional redundancy in biodiversity research [https://pmc.ncbi.nlm.nih.gov/articles/PMC11332098/]
• 23. Response of soil general and specific functions following ... [https://www.sciencedirect.com/science/article/pii/S266700622400025X]
• 24. Research needs on the biodiversity–ecosystem functioning ... [https://www.nature.com/articles/s44185-024-00046-6]
• 25. The critical role of multidimensional biodiversity in... | CoLab [https://colab.ws/articles/10.1016%2Fj.geosus.2023.05.002]
• 26. Soil Health | Natural Resources Conservation Service [https://www.nrcs.usda.gov/conservation-basics/natural-resource-concerns/soil/soil-health]
• 27. Global land use changes are four times greater than ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8113269/]
• 28. Climate Change Effects on Land Use and Land Cover ... [https://www.mdpi.com/2073-445X/13/12/2008]
• 29. The Jena Experiment [https://the-jena-experiment.de/]
• 30. Biodiversity–stability relationships strengthen over time in a ... [https://www.nature.com/articles/s41467-022-35189-2]
• 31. Population dynamics of an experimental plot of biodiversity ... [https://www.frontiersin.org/journals/ecology-and-evolution/articles/10.3389/fevo.2025.1581588/full]
• 32. Microbial population structures in soil particle size fractions ... [https://pubmed.ncbi.nlm.nih.gov/11526026/]
• 33. Soil microbiome indicators can predict crop growth ... [https://www.nature.com/articles/s41564-023-01520-w]
• 34. Restoring soil biodiversity [https://www.sciencedirect.com/science/article/pii/S0960982224001775]
• 35. A novel two-step metabarcoding approach improves soil ... [https://www.nature.com/articles/s41598-025-18936-5]
• 36. DNA metabarcoding of mites from small soil samples [https://pmc.ncbi.nlm.nih.gov/articles/PMC12548636/]
• 37. Comparing soil microarthropod communities derived ... [https://www.sciencedirect.com/science/article/pii/S0929139325001805]
• 38. Metagenomics Provides a Deeper Assessment of the ... [https://www.mdpi.com/2073-4425/14/4/812]
• 39. Functional diversity of soil microbial communities increases ... [https://www.nature.com/articles/s41467-025-66544-8]
• 40. Stable Soil Microbial Functional Structure Responding to ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2021.716764/full]
• 41. 长期覆膜条件下农田土壤微生物群落的响应特征 [http://www.hjkx.ac.cn/hjkx/ch/html/20221045.htm]
• 42. Diversity and asynchrony in soil microbial ... [https://elifesciences.org/articles/62813]
• 43. Energy flux: The link between multitrophic biodiversity and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6181201/]
• 44. SOIL - Benchmarking soil multifunctionality - Copernicus.org [https://soil.copernicus.org/articles/11/609/2025/]
• 45. The three major axes of terrestrial ecosystem function [https://www.nature.com/articles/s41586-021-03939-9]
• 46. A five-stage protocol for systematic measuring soil carbon ... [https://link.springer.com/article/10.1007/s11104-025-07703-0]
• 47. Quantification of soil N 2 O and CH 4 fluxes using the ... [https://www.sciencedirect.com/science/article/abs/pii/S0168192323004525]
• 48. Soil heterotrophic and autotrophic respiration respond ... [https://www.nature.com/articles/s41598-025-06082-x]
• 49. Trends of Soil and Solution Nutrient Sensing for Open Field ... [https://www.mdpi.com/1424-8220/25/2/453]
• 50. Soil multitrophic network complexity enhances the link between biodiversity and multifunctionality in agricultural systems - PubMed [https://pubmed.ncbi.nlm.nih.gov/34610173/]
• 51. A multitrophic perspective on biodiversity–ecosystem functioning research - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC6944504/]
• 52. Best Statistical Tools for Data Analysis and Visualization [https://www.6sigma.us/six-sigma-in-focus/statistical-tools/]
• 53. A global dataset of experimental agricultural management ... [https://www.nature.com/articles/s41597-025-05238-8]
• 54. Comparison of Machine Learning Algorithms for Predictive ... [https://www.nature.com/articles/s41598-019-40927-6]
• 55. Top Predictive Analytics Models and Algorithms to Know [https://insightsoftware.com/blog/top-5-predictive-analytics-models-and-algorithms/]
• 56. The role of AI-enhanced microscopy in soil biodiversity ... [https://www.sciencedirect.com/science/article/pii/S2667006225000280]
• 57. Effect of plant communities on bacterial and fungal ... [https://environmentalmicrobiome.biomedcentral.com/articles/10.1186/s40793-024-00583-4]
• 58. Soil fungi influence the relationship between plant diversity ... [https://www.nature.com/articles/s41467-025-60661-0]
• 59. Fungal diversity drives soil multifunctionality after new tree ... [https://www.sciencedirect.com/science/article/abs/pii/S0929139325005037]
• 60. Plant Species Richness and the Root Economics Space ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11687415/]
• 61. Links between Soil Fungal Diversity and Plant ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5682006/]
• 62. Experimental Evidence of Functional Group-Dependent ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2018.02312/full]
• 63. Soil Fungal Diversity, Community Structure, and Network ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12113623/]
• 64. Biodiversity of network modules drives ecosystem functioning ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10847562/]
• 65. Climate change is predicted to reduce global belowground ... [https://www.nature.com/articles/s41467-025-64453-4]
• 66. Climate and soil effects on tree species diversity and ... [https://www.nature.com/articles/s41598-023-38225-3]
• 67. Geographic bias of field observations of soil carbon stocks ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3076837/]
• 68. Geographic, Taxonomic and Metric Gaps in Biodiversity ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12513691/]
• 69. A representativeness heuristic for mitigating spatial bias in ... [https://www.sciencedirect.com/science/article/abs/pii/S0016706118319268]
• 70. Comparing the Soil Management Assessment Framework ... [https://www.mdpi.com/2073-4395/15/3/643]
• 71. Trophic interaction modifications: an empirical and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6849598/]
• 72. 褚海燕组-土壤生物多样性与城市绿地生态系统功能 [https://agri.ucas.ac.cn/index.php/zh/xueshukeyan/keyandongtai/1041-2023-01-11-02-08-29]
• 73. Species loss in key habitats accelerates regional food web ... [https://www.nature.com/articles/s42003-025-08396-y]
• 74. Linking species richness, biodiversity and ecosystem function ... [https://research.fs.usda.gov/treesearch/25438]
• 75. 土壤生态系统微生物多样性-稳定性关系的思考 [https://www.biodiversity-science.net/CN/Y2013/V21/I4/411]
• 76. Plant diversity maintains multiple soil functions in future ... [https://elifesciences.org/articles/41228]
• 77. Non-linearity impacts of landscape pattern on ecosystem ... [https://www.sciencedirect.com/science/article/pii/S1470160X22000796]
• 78. Spatiotemporal evolution and trade-offs/synergies of ... [https://www.nature.com/articles/s41598-025-19570-x]
• 79. Non-Linear Impacts of Social and Ecological Drivers on ... [https://www.mdpi.com/2413-8851/9/10/390]
• 80. Plant-soil feedback: the next generation [https://link.springer.com/article/10.1007/s11104-023-06000-y]
• 81. Blind spots in global soil biodiversity and ecosystem function ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7400591/]
• 82. Diversity-dependent temporal divergence of ecosystem ... - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC5659383/]
• 83. Top 6 methods for soil sampling efficiency [https://www.cropway.com/soil-sampling/]
• 84. Soil Testing and Collection Tips and Techniques for 2025 [https://alluvialsoillab.com/blogs/news/soil-testing-and-collection-tips-and-techniques-for-2025?srsltid=AfmBOoqUwBnBU_XqLfAxSCTeDVMxf8_mm3cCVw3Qc53Ii3qayMvdE7pN]
• 85. Positive plant diversity-soil stability relationships are mediated ... [https://research.fs.usda.gov/treesearch/45572]
• 86. Optimizing Soil Core Sampling Strategies for Enhanced ... [https://pure.au.dk/portal/en/projects/optimizing-soil-core-sampling-strategies-for-enhanced-soil-invert-2/]
• 87. Crop management strategies shape the shared temporal ... [https://www.sciencedirect.com/science/article/pii/S0167880924001762]
• 88. Effects of plant diversity on productivity strengthen over ... [https://www.nature.com/articles/s41467-024-46355-z]
• 89. Using forest gap models and experimental data to explore ... [https://annforsci.biomedcentral.com/articles/10.1007/s13595-020-00954-0]
• 90. Field observed relationships between biodiversity and ... [https://www.sciencedirect.com/science/article/abs/pii/S1146609X13001525]
• 91. Biodiversity and Ecosystem Function Under Global Change [https://www.mdpi.com/2079-7737/14/11/1503]
• 92. Major changes in climate, vegetation, and ... - Fire Ecology [https://fireecology.springeropen.com/articles/10.1186/s42408-025-00388-x]
• 93. Soil nutrient limitation controls trophic cascade effects of ... [https://www.sciencedirect.com/science/article/pii/S2090123225000372]
• 94. Nutrient-induced acidification modulates soil biodiversity- ... [https://www.nature.com/articles/s41467-024-47323-3]
• 95. 草原生态系统地上下生物多样性关系研究取得进展 [https://www.cas.cn/syky/202203/t20220303_4826932.shtml]
• 96. Soil biota contributions to soil aggregation - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC5701735/]
• 97. PP387/PP387: Biology, Ecology, and Benefits of Ectomycorrhizal Fungi in Natural and Managed Forest Ecosystems [https://edis.ifas.ufl.edu/publication/PP387]
• 98. Earthworms Building Up Soil Microbiota, a Review [https://www.frontiersin.org/journals/environmental-science/articles/10.3389/fenvs.2019.00081/full]
• 99. Effects of arbuscular mycorrhizal fungi and earthworms on ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12619942/]
• 100. Multiple elements of soil biodiversity drive ecosystem ... functions [https://www.nature.com/articles/s41559-019-1084-y?error=cookies_not_supported&code=0171d70c-6444-48fe-bfc6-8a91e000b654]
• 101. and asynchrony in Diversity microbial communities stabilizes... soil [https://pmc.ncbi.nlm.nih.gov/articles/PMC7987343/]
• 102. Dispersal scales up the biodiversity–productivity relationship ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC2894895/]
• 103. Soil community history strengthens belowground ... [https://www.nature.com/articles/s41467-024-54401-z]
• 104. Soil biodiversity and network complexity jointly drive ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2025.1668494/full]
• 105. Effects of agricultural intensification on soil biodiversity and ... [https://asu.elsevierpure.com/en/publications/effects-of-agricultural-intensification-on-soil-biodiversity-and-/]
• 106. Digital Soil Mapping Platforms Market 2025 Growth Rate [https://farmonaut.com/precision-farming/digital-soil-mapping-platforms-market-2025-growth-rate]
• 107. Soil Organic Carbon Projections and Climate Adaptation ... [https://egusphere.copernicus.org/preprints/2025/egusphere-2025-4258/]
• 108. Predictive Modeling of Soil Carbon Sequestration with AI [https://prism.sustainability-directory.com/scenario/predictive-modeling-of-soil-carbon-sequestration-with-ai/]