Tropical land use alters functional diversity of soil food webs and leads to monopolization of the detrital energy channel

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Agricultural expansion is among the main threats to biodiversity and functions of tropical ecosystems. It has been shown that conversion of rainforest into plantations erodes biodiversity, but further consequences for food-web structure and energetics of belowground communities remains little explored. We used a unique combination of stable isotope analysis and food-web energetics to analyze in a comprehensive way consequences of the conversion of rainforest into oil palm and rubber plantations on the structure of and channeling of energy through soil animal food webs in Sumatra, Indonesia. Across the animal groups studied, most of the taxa had lower litter-calibrated Δ¹³C values in plantations than in rainforests, suggesting that they switched to freshly-fixed plant carbon ('fast' energy channeling) in plantations from the detrital C pathway ('slow' energy channeling) in rainforests. These shifts led to changes in isotopic divergence, dispersion, evenness, and uniqueness. However, earthworms as major detritivores stayed unchanged in their trophic niche and monopolized the detrital pathway in plantations, resulting in similar energetic metrics across land-use systems. Functional diversity metrics of soil food webs were associated with reduced amount of litter, tree density, and species richness in plantations, providing guidelines on how to improve the complexity of the structure of and channeling of energy through soil food webs. Our results highlight the strong restructuring of soil food webs with the conversion of rainforest into plantations threatening soil functioning and ecosystem stability in the long term.

Editor's evaluation

Zhou et al., provide a robust study on isotopic and metabolic changes of a soil community across a gradient of different land-use types in Sumatra, Indonesia. By mixing community-based analyses of stable isotopes and size-based metabolic measures, they are able to elucidate, for the first time, important links among plants and the soil food web in tropical ecosystems. This study is of importance to tropical biologists, ecosystem ecologists and biodiversity conservationists aiming to understand the impacts of humans on tropical forests.

https://doi.org/10.7554/eLife.75428.sa0

Introduction

Worldwide, land use changes the structure of ecological communities and is associated with losses in multiple ecosystem functions, which is at the core of sustainable development goals (Bommarco et al., 2013; Matson et al., 1997; Newbold et al., 2015). Many tropical ecosystems are affected by land-use changes, losing their biodiversity and multifunctionality (Barnes et al., 2014; Laurance, 2007). It is projected that tropical ecosystems will face even greater pressures due to land-use change in the future (Dobrovolski et al., 2011). Decreases in biodiversity and changes in trophic interactions in animal communities (Newbold et al., 2015; Tsiafouli et al., 2015; Wilkinson et al., 2021) are associated with changes in nutrient dynamics and energy fluxes (de Vries et al., 2012; McGrath et al., 2001; Potapov et al., 2020), which ultimately influence ecosystem functioning and stability (Rooney et al., 2006; Rooney and McCann, 2012). However, interrelationships between the loss of diversity and changes in energy pathways in food webs are poorly studied and this applies in particular to tropical ecosystems.

Soils harbor a large portion of terrestrial biodiversity (Guerra et al., 2021), are intimately linked with aboveground biodiversity (Bardgett and Putten, 2014; Hooper et al., 2000; Yang et al., 2018), and deliver vital ecosystem services (Bardgett and Wardle, 2010; de Vries et al., 2013). Energetically, 80–90% of the carbon fixed by plants in terrestrial ecosystems enters the belowground system (Gessner et al., 2010) and is processed in soil food webs by microorganisms and invertebrate decomposers, the latter then become prey for predators (Bardgett and Wardle, 2010; Schmitz and Leroux, 2020). Shifts in resource use in the decomposer system results in asymmetries in energy fluxes through soil food-web channels, which modulate the resistance and resilience of terrestrial ecosystems to perturbations (de Vries et al., 2006; de Vries et al., 2012; Rooney et al., 2006; Rooney and McCann, 2012). Studies in temperate regions showed that more intensive land use reduces the diversity of soil organisms (Tsiafouli et al., 2015) and shifts soil food webs toward the ‘fast’ bacterial energy channel at the expense of the ‘slow’ fungal energy channel (de Vries et al., 2006), potentially undermining food-web stability. However, knowledge on how the rapid land-use change in tropical regions, such as the conversion of rainforest into plantations, affects soil food-web structure and energy channeling is scarce (Clough et al., 2016; Dobrovolski et al., 2011).

The present study took place in Jambi province, Sumatra, Indonesia, which is a global hotspot of biodiversity (Koh and Ghazoul, 2010; Miettinen et al., 2011), where over last 25–35 years rainforests and agroforests have been largely replaced by intensively managed plantations, mostly oil palm and rubber (Clough et al., 2016; Margono et al., 2012). Results of previous studies showed that land-use change in this region is associated with changes in soil chemistry (Ballauff et al., 2021), shifts in microbial and plant communities (Krashevska et al., 2015; Rembold et al., 2017a; Schulz et al., 2019), and reduced multitrophic biodiversity and functionality of soil animal communities (Barnes et al., 2014; Krause et al., 2021; Potapov et al., 2019a; Krashevska et al., 2019). These changes are expected to affect trophic niches of soil animals and alter both structure and energetics of soil food webs. In certain soil invertebrate groups (e.g., centipedes, springtails, and mites), conversion of rainforest into plantation systems has been associated with trophic shifts toward the plant energy channel in plantations (Klarner et al., 2017; Krause et al., 2021; Susanti et al., 2021), whereas the bacterial channel was reduced, as suggested by fatty acid analysis (Susanti et al., 2019). A more complete assessment of soil food webs showed that plantations are energetically dominated by large decomposers (i.e., earthworms), but have largely reduced energy fluxes to predators (Barnes et al., 2014; Potapov et al., 2019a), however, these studies ignored potential shifts in the trophic niches of individual soil taxa with land-use change.

Progress in understanding soil food-web responses to environmental changes is hampered by the chronic lack of empirical data for complex soil food webs (Brose and Scheu, 2014). Insights in their structure became possible with introduction of stable isotope, molecular, and biochemical methods, which showed inaccuracies in the traditional reconstructions (Bradford, 2016; Brose and Scheu, 2014; Geisen et al., 2019). Stable isotope analysis is now widely used as a first-line explorative tool in trophic ecology (Peterson and Fry, 1987; Parnell et al., 2010), allowing for in situ assessment of soil food-web structure (Potapov et al., 2019c). The method is especially promising to provide insight into the structure of soil food webs in the tropics, where the biology of species is poorly known. The ¹³C/¹²C and ¹⁵N/¹⁴N ratios in consumers depend on their food and can be used to explore the trophic niches of animal species and communities (Post, 2002a; Pollierer et al., 2009; Potapov et al., 2019c). The ¹⁵N/¹⁴N isotope ratio is used to indicate the trophic position of species since it is enriched by about 3–4‰ per trophic level (Post, 2002a; Pollierer et al., 2009; Potapov et al., 2019c); ¹³C typically is little affected by trophic transfer and thus reflects basal food resources of the trophic chain (Peterson and Fry, 1987; Potapov et al., 2019c). In soil communities, animals with high ¹³C concentration are considered to use ‘older’ carbon that has higher ¹³C values due to decomposition processes and preferential incorporation of labile plant compounds by microbes (Pollierer et al., 2009; Potapov et al., 2019c), and those with lower ¹³C concentration are considered to feed on freshly fixed plant material (Fujii et al., 2021; Potapov et al., 2019c).

To assess food-web structure using stable isotope analysis, Layman et al., 2012; Layman et al., 2007, suggested a number of ‘isotopic metrics’, which have been widely used in aquatic ecology. These metrics consider all species as having the same importance for food-web structure, which has a binary perspective (i.e., presence/absence of species), ignoring potential asymmetries in the magnitude of trophic interactions. However, in biological communities often only few species dominate, forming the energetic core of the food web, therefore, a non-binary perspective is important. Recently, Cucherousset and Villéger, 2015, joined Layman’s metrics and the functional diversity framework (Petchey and Gaston, 2006; Villéger et al., 2008; Mouillot et al., 2013) to calculate functional diversity indices for food webs, accounting for the dominance of species. These indices include isotopic ‘richness’ which represents the volume of the trophic niche across all species, isotopic ‘evenness’ which represents the regularity of the distribution of species’ trophic niches, isotopic ‘divergence’ which reflects the dominance of species with the most extreme trophic niches, and isotopic ‘dispersion’ which reflects the balance of the species distribution in the trophic space (Cucherousset and Villéger, 2015). These isotopic indices represent basic components of the functional diversity of food webs. Nevertheless, to our knowledge they have never been used to analyze soil food-web characteristics, either temperate or tropical, except for one case study on oribatid mites (Krause et al., 2021). Moreover, abundance and biomass each are biased in reflecting the functional role of consumers covering wide body size ranges. While abundance is biased toward the importance of small organisms, biomass is biased toward that of large ones. Considering these limitations, energetic demands of consumers (i.e., metabolic rates) may be used as less biased metric (Brown et al., 2004). In recent years, the energy flux approach was successfully used to represent functional changes in food webs, and therefore to link multitrophic biodiversity to ecosystem functioning (Barnes et al., 2018; Barnes et al., 2014; Jochum et al., 2021). To the best of our knowledge, however, the energy flux approach has never been used in conjunction with stable isotope analysis.

Here, for the first time we use stable isotope analysis to comprehensively investigate changes in tropical soil food webs associated with changes in land use. We apply a functional diversity framework to stable isotope data to assess which structural dimensions of soil food webs vary most across rainforests, agroforests, and intensively managed plantations of oil palm and rubber in Jambi province, Sumatra, Indonesia (Clough et al., 2016; Drescher et al., 2016). Using data on 23 high-rank taxonomic groups (orders, families), we focus on two perspectives of the functional diversity of soil food webs: a ‘community perspective’ in which we treat all groups as being equally important and an ‘energetic perspective’ in which we weight groups according to their shares in community metabolism. For both perspectives, we tested the following hypotheses: (1) shifts in trophic niches are uniform across all studied animal groups through land-use changes, with animals in plantations being less enriched in ¹³C than in rainforest due to stronger plant and weaker detrital energy channel; (2) functional diversity of soil food webs declines with land-use intensity in plantation systems reflected by reduced isotopic richness, redundancy, evenness, and divergence; (3) from an ‘energetic perspective’ soil food webs are less affected by changes in land use than from a ‘community perspective’ as total energy flux changes little with conversion of rainforest into plantations, whereas biodiversity declines strongly. Lastly, we aim at identifying the environmental factors driving changes in functional diversity of soil food webs under land-use change, both from a community and energetic perspective.

Results

Isotopic shifts in individual animal taxa

If averaged across rainforest sites, the mean Δ¹³C values of taxonomic groups covered a range of 3.0‰, from 3.4‰ (Coleoptera and Oribatida) to 6.4‰ (Orthoptera and Pauropoda). The respective Δ¹⁵N values covered a range of 14.2‰, from –5‰ (Pauropoda) to 9.2‰ (Diplura; Figure 1a). The Δ¹³C values of all groups varied across the four land-use systems, but Chilopoda, Diplura, and Annelida were typically most enriched in ¹³C and Coleoptera were consistently among the most depleted in ¹³C among all groups. Diplura, Pseudoscorpiones, Chilopoda, and Isopoda had the highest Δ¹⁵N values among all groups across the four land-use systems. The group with lowest Δ¹⁵N values was Pauropoda in rainforest and jungle rubber, and Protura in rubber and oil palm plantations. Overall, micropredators, that is, Diplura and Pseudoscorpiones, had 2–3‰ higher Δ¹⁵N values than macropredators, that is, Chilopoda, Formicidae, and Araneae (Figure 1). The share of Annelida (earthworms) in community metabolism was 15.4% in rainforest, but represented more than 75% in jungle rubber, rubber, and oil palm plantations (Figure 1—figure supplement 1).

Figure 1 with 1 supplement see all

Download asset Open asset

Mean litter-calibrated Δ¹³C and Δ¹⁵N values of soil animal taxa in rainforest (a), jungle rubber (b), rubber, (c) and oil palm plantations (d).

Error bars represent standard errors across sampling plots (n = 1–8 per land-use system). Size of the points is scaled to the total share of the taxonomic group in the community metabolism in the corresponding land-use system (metabolism was log₁₀-transformed to show trends in rarer groups).

Figure 1—source data 1 Metabolism data of each group in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig1-data1-v3.xlsx
Download elife-75428-fig1-data1-v3.xlsx
Figure 1—source data 2 Stable isotope data of groups in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig1-data2-v3.xlsx
Download elife-75428-fig1-data2-v3.xlsx

The Δ¹³C values were significantly higher in rainforest than in the other land-use systems in Coleoptera, Diplopoda, Hemiptera, Orthoptera, Pauropoda, Protura, Pseudoscorpiones, and Thysanoptera (Figure 2, Figure 2—figure supplement 1). Chilopoda, Diplura, Formicidae, Isopoda, Mesostigmata, and Symphyla were significantly more enriched in ¹³C in rainforest than in oil palm, but not significantly different from those in jungle rubber and rubber plantations. In general, most groups in rainforest were higher in Δ¹³C by 1–3‰ than in the other land-use systems, but this shift was only significant for two out of six macrodecomposer groups. Annelida, which accounted for much of the community metabolism in each of the land-use systems, had similar Δ¹³C values across land-use systems.

Figure 2 with 2 supplements see all

Download asset Open asset

Average Δ¹³C and Δ¹⁵N values of taxonomic groups in rainforest (F), jungle rubber (J), rubber, (R) and oil palm plantations (O).

Numbers show means, asterisks indicate significant differences between the mean value in the corresponding land-use system and in rainforest (Student’s t-test *p < 0.05, **p < 0.01). Color represents the direction (red – increase, blue – decrease) and magnitude (darker color indicate stronger change) of the difference between rainforest and other land-use systems.

Figure 2—source data 1 Metabolism data of each group in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig2-data1-v3.xlsx
Download elife-75428-fig2-data1-v3.xlsx

The Δ¹⁵N values were by 1.5–2.5‰ lower in rainforest than in the other land-use systems in Oribatida and Blattodea (except for similar Δ¹⁵N values in oil palm and rainforest for Blattodea). By contrast, Δ¹⁵N values of Hemiptera, Orthoptera, Isoptera, and Pseudoscorpiones were lower in oil palm than in rainforest, whereas in jungle rubber this was only true for Orthoptera (Figure 2, Figure 2—figure supplement 2).

One-dimensional isotopic metrics

One-dimensional isotopic metrics described the overall range and average Δ¹³C and Δ¹⁵N values of each community. The maxima of Δ¹³C values were by 1–2‰ higher in forest than in jungle rubber and oil palm plantations, but minima and the overall range did not differ significantly (Figure 3). The unweighted average Δ¹³C values of communities were by 1–2‰ higher in rainforest than in the other land-use systems and were also higher in rubber than in oil palm plantations. However, the energetic average positions did not differ significantly due to similar Δ¹³C values of Annelida (dominant invertebrate group) across land-use systems (Figure 2e, Figure 1—figure supplement 1).

Figure 3

Download asset Open asset

One-dimensional metrics for Δ¹³C (upper panel) and Δ¹⁵N values (lower panel) of communities in rainforest (F, green), jungle rubber (J, blue), rubber (R, red), and oil palm plantations (O, yellow).

Each point represents one community (n = 8 per land-use system). For the calculation of the weighted average values, species were weighted according to their contribution to the total community metabolism per plot. Means sharing the same letter within each pane are not significantly different (Tukey’s HSD test following ANOVA, p > 0.05).

Figure 3—source data 1 Community metrics of soil food webs in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig3-data1-v3.csv
Download elife-75428-fig3-data1-v3.csv
Figure 3—source data 2 Energetic metrics of soil food webs in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig3-data2-v3.csv
Download elife-75428-fig3-data2-v3.csv

Extreme values of Δ¹⁵N were most pronounced in jungle rubber, both maximum and minimum, resulting in the largest range in Δ¹⁵N values (16.5‰) among the four land-use systems. In the other land-use systems, maxima, minima, and ranges of Δ¹⁵N values were similar. The unweighted average Δ¹⁵N values of communities were lowest in oil palm, being significantly lower than in rubber (Figure 3i). By contrast, the energetic average Δ¹⁵N values did not differ significantly.

Multidimensional isotopic metrics

Among the unweighted metrics (i.e., community metrics), isotopic dispersion was significantly higher in oil palm than in each of the other land-use systems; isotopic divergence and uniqueness were significantly higher in oil palm than in jungle rubber; isotopic evenness was significantly lower in jungle rubber than in each of the other land-use systems; only isotopic richness showed no significant differences between land-use systems, but in trend the two monoculture systems had lower values than in rainforest and jungle rubber. For detailed information on the plot-level metrics values, see Appendix (Figure 4—figure supplements 1 and 2).

By contrast, the weighted multidimensional metrics (i.e., energetic metrics) did not differ among land-use systems for isotopic dispersion, isotopic evenness, isotopic richness, and isotopic uniqueness (Figure 4). Only isotopic divergence was significantly lower in rainforest than in the other land-use systems, showing an opposite trend to isotopic dispersion. For detailed information on plot level metrics values see Appendix (Figure 4—figure supplements 3–6).

Figure 4 with 64 supplements see all

Download asset Open asset

Multidimensional isotopic metrics of soil animal communities in rainforest (F, green), jungle rubber (J, blue), rubber (R, red), and oil palm plantations (O, yellow).

Community (upper panel) and energetic metrics (lower panel) are shown. Each point represents one community (n = 8 per land-use system). Means sharing the same letter within each pane are not significantly different (Tukey’s HSD test following ANOVA, p < 0.05).

Figure 4—source data 1 Community metrics of soil food webs in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig4-data1-v3.csv
Download elife-75428-fig4-data1-v3.csv
Figure 4—source data 2 Energetic metrics of soil food webs in each plot.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig4-data2-v3.csv
Download elife-75428-fig4-data2-v3.csv

Environmental effects on functional diversity of soil food webs

As indicated by multivariate analysis, community metrics differed strongly between the four land-use systems (anosim R = 0.404, p < 0.001), whereas differences for the energetic metrics were less pronounced (anosim R = 0.138, p = 0.014). Among all tested environmental factors, soil pH, tree species richness, litter amount, and understory density had the strongest correlations with both community metrics (p = 0.003, p = 0.037, p < 0.001, p < 0.001, respectively) and energetic metrics (p = 0.004, p = 0.002, p = 0.003, p = 0.002, respectively). These variables were subsequently selected for the structural equation model (SEM) analysis (see Materials and methods). SEM indicated that the changes in the community metrics (PC1_unweighted) were induced directly by tree properties and litter amount (tree density: p < 0.05, effect size = 0.72; tree species richness: p < 0.001, effect size = –0.83; litter amount: p < 0.001, effect size = 0.71), while changes in the energetic metrics (PC1_weighted) were indirectly driven by soil pH via increased metabolism of earthworms (p < 0.05, effect size = –0.36; Figure 5; Figure 6).

Figure 5

Download asset Open asset

Environmental drivers of community and energetic soil food-web metrics.

Community and energetic food-web metrics were related to environmental factors using multivariate analysis of variance (MANOVA); the thickness of connection lines shows statistical significance, dashed line for p > 0.05. Pairwise Spearman’s correlations among environmental factors are shown with a tile chart (blue – negative, red – positive). The vegetation parameters included tree species richness (TreeRich), tree density (TreeDen), understory species richness (UnderRich), understory density (UnderDen), and average understory height (UnderHeight). Parameters of litter and soil include soil pH, litter amount, soil carbon concentration (Csoil), carbon-to-nitrogen ratio of litter (CtoNlitter), soil microbial biomass C (Soil Cmic), and soil humidity (Water soil) (Krashevska et al., 2015; Rembold et al., 2017a).

Figure 5—source data 1 Data of environmental factors.: https://cdn.elifesciences.org/articles/75428/elife-75428-fig5-data1-v3.csv
Download elife-75428-fig5-data1-v3.csv

Figure 6

Download asset Open asset

Structural equation model on the effects of environmental change on food-web metrics.

Numbers adjacent to arrows are standardized path coefficients that show effect sizes and directions (blue – positive, red – negative) of the relationship, arrow width is proportional to the strength of path coefficients. Gray arrows represent paths that were not significant; *p < 0.05, **p < 0.01, and ***p < 0.001. Numbers above every response variable in the model denotes the proportion of variance explained. For abbreviations, see Figure 5.

Figure 6—source data 1 Data for building structural equation model (SEM).: https://cdn.elifesciences.org/articles/75428/elife-75428-fig6-data1-v3.xlsx
Download elife-75428-fig6-data1-v3.xlsx

Discussion

We used stable isotope data of 23 high-rank animal taxa to comprehensively assess changes in functional diversity of soil food webs under tropical land-use change. We found shifts in basal resource use for most taxonomic groups in plantations compared to rainforest, and responses of food-web diversity metrics to land-use change were more pronounced for community than for energetic metrics. In agreement to our first hypothesis, ¹³C values of animal taxa and communities were more enriched in rainforest than in plantations, but this shift vanished if the average Δ¹³C values were weighted by metabolism. Soil animals in jungle rubber had the largest range of Δ¹⁵N values among all land-use systems, which suggests the longest food chain in this system. Refuting our second hypothesis, when considering all taxonomic groups being equally important (‘community perspective’), soil food webs in oil palm had a significantly higher community dispersion than in the other land-use systems, and in trend also had a higher community divergence and uniqueness (i.e., a lower redundancy). By contrast, most energetic isotopic metrics (‘energetic perspective’) varied less between the four land-use systems. Conform to our third hypothesis, community isotopic metrics were more sensitive to changes in land use than energetic isotopic metrics, suggesting soil food webs are more sensitive to land-use change from a community than from an energetic perspective. Further, community metrics of soil food webs were influenced directly by land use (tree properties and litter amount), whereas energetic metrics were influenced indirectly via pH-induced changes in earthworm abundance.

The structure of tropical soil food webs

Our study is among the first comprehensive assessment of tropical soil food webs based on stable isotope analysis. Collembola, Symphyla, and Isopoda showed a much higher ¹⁵N enrichment than, for example, Oribatida, but all three groups occupy similar trophic positions in temperate forests and predominantly function as decomposers (Potapov et al., 2019c). This difference may be caused by low litter quality in tropical forests forcing decomposers to switch to more microbial or even animal diet (Illig et al., 2005). Protura in temperate forests are enriched in ¹⁵N and feed on ectomycorrhizal fungi (Bluhm et al., 2019), whereas in the studied tropical forests, Protura had the lowest Δ¹⁵N values among all groups, suggesting that they feed on saprotrophic rather than mycorrhizal fungi. The low Δ¹⁵N and high Δ¹³C values of Pauropoda, reported for the first time for this group, indicate that they function as decomposers by feeding on saprotrophic microorganisms (Tiunov et al., 2015), confirming earlier suggestions (Starling, 1944). Low Δ¹⁵N values in Protura, Diplopoda, Isoptera, Psocoptera, and Blattodea may be associated with feeding on algae (Potapov et al., 2018), shown to be important food for mesofauna in tropical soil food webs (Susanti et al., 2019; Semenina et al., 2020). Unexpectedly, micropredators (e.g., Diplura and Pseudoscorpiones) had higher trophic positions (Δ¹⁵N values) than macropredators (e.g., Araneae and Formicidae) across all land-use systems, and Diplura had the highest Δ¹⁵N values among all taxa studied. Diplura were represented mostly by predatory Japygidae, which may hunt springtails, mites, and other small invertebrates (Sendra et al., 2021). The higher trophic position of small-sized predators suggests that they form part of a different energy channel than macropredators. In fact, the micro-food web in soil has been shown to be based mainly on microbial resources channeled to higher trophic levels by microarthropod predators, whereas the macro-food web is based more on litter and detritus consumed by macrofauna taxa with the energy channeled to higher trophic levels by macroarthropod predators (Potapov et al., 2021a). This implies more trophic transactions in the micro-food web (Pollierer et al., 2009; Steffan et al., 2015) explaining the higher trophic position of micro- than macroarthropod predators. Among other groups with a wide food spectrum, Diptera had higher isotopic values (both Δ¹⁵N and Δ¹³C) than Coleoptera, indicating that flies are more closely linked to detrital and microbial food chains in tropical soil food webs than beetles. Isopoda had strikingly high Δ¹⁵N values for macrodecomposers, possibly due to intense coprophagy (Potapov et al., 2022). Overall, despite general similarities, we also found consistent differences between tropical and temperate soil food webs. Further studies comparing differences in soil food-web structure and associated soil functions across temperate and tropical ecosystems are needed to prove the generality of these differences and their consequences for biodiversity – ecosystem functioning relationships.

Shifts in trophic positions of soil invertebrates and energy channels in soil food webs

In agreement with our first hypothesis, Δ¹³C values of most of the studied soil animal taxa were higher in rainforest than in plantations. The ¹³C concentration in dead plant material is increasing during decomposition compared to fresh leaf litter (Ågren et al., 1996; Boström et al., 2007; Potapov et al., 2019c), and high Δ¹³C values in soil fauna in forest likely indicate feeding on saprotrophic fungi and bacteria that assimilate predominantly labile ¹³C-enriched plant compounds (Pollierer et al., 2009; Potapov et al., 2013; Hyodo, 2015). Vascular and non-vascular plants have generally lower Δ¹³C values than saprotrophic microorganisms and animals (Hyodo et al., 2010; Potapov et al., 2019c), therefore, the high Δ¹³C values in soil invertebrates in rainforest point to a more pronounced detritus-based ‘brown’ food web relying heavily on saprotrophic fungi and bacteria based on litter material. Among the plantations, the unweighted average Δ¹³C values were lowest in oil palm suggesting a shift toward a more plant-based ‘green’ food web relying more heavily on the consumption of living plant tissue (Fujii et al., 2021), which has been previously shown for Chilopoda, Oribatida, Collembola, and Pseudoscorpiones at the same study sites (Klarner et al., 2017; Krause et al., 2019; Liebke et al., 2021; Susanti et al., 2021). Results of the study of Susanti et al., 2019, using fatty acids as trophic biomarkers at our study sites further support the conclusion of a more pronounced plant- and reduced detritus-based energy channel in soil food webs of plantations compared to rainforest. Compared to rainforest the herb layer is much more developed in plantations due to more open canopy and coverage by weeds (Rembold et al., 2017a), presumably providing high-quality resources for plant and litter feeding soil animals. By contrast, litter in rainforest is high in lignin and therefore of low food quality (Krashevska et al., 2018), increasing the use of saprotrophic microorganisms rather than litter by detritivores (Illig et al., 2005).

In soil food-web models the bacterial energy channel typically is considered to be ‘fast’, while the fungal energy channel is considered to be ‘slow’ due to differences in the turnover rates of bacteria and their consumers, and fungi and their consumers, respectively (Coleman et al., 1983; Moore et al., 2005). Generally, the fast energy channel is defined by fast-growing populations with short turnover rates (Rooney et al., 2006). Based on this concept, animals feeding on living plants or fresh litter contribute to fast energy channeling, whereas animals feeding on recalcitrant detritus (being slowly decomposed belowground) are associated with slow energy channeling. This perspective complements the green (plant) and brown (detrital) energy channels in soil food webs. Since ¹³C accumulates in plant-derived materials with microbial decomposition (Pollierer et al., 2009; Potapov et al., 2013; Fujii et al., 2021), Δ¹³C values in soil animals may indicate slow versus fast carbon cycling in soil. Therefore, the shift toward plant-based energy channeling in food webs of plantations may reflect accelerated carbon cycling at the ecosystem level, contributing to carbon losses of plantations (Guillaume et al., 2018) and thereby potentially compromising ecosystem long-term stability (McCann et al., 1998; Rooney and McCann, 2012).

Contrasting the community perspective, energetic Δ¹³C values of communities did not vary significantly among land-use systems, which was due to similar Δ¹³C values of Annelida (earthworms) and some other macro-decomposers across land-use systems. Earthworms had the highest share in community metabolism among detritivores in plantations suggesting that they predominate animal-mediated decomposition and organic matter transformation processes (Potapov et al., 2019a). Notably, earthworms were among the most ¹³C-enriched animal groups in jungle rubber, rubber and oil palm plantations, but their Δ¹³C values were similar to other soil animal groups in rainforest. Earthworms feed on detritus and microorganisms and are able to efficiently use fresh litter carbon, but also ‘old’ microbially processed carbon (e.g., soil organic matter; Scheu and Falca, 2000; Hyodo et al., 2012; Blouin et al., 2013) reflected in high Δ¹³C values (Pollierer et al., 2009). Land-use change in tropical lowland landscapes is associated with the loss of biodiversity and reduced biomass of litter arthropods (Barnes et al., 2014), but the negative effect of biodiversity loss on soil functions may be at least in part counteracted by earthworms that monopolize the detrital channel in plantations. Thereby, earthworms also may counteract the destabilization of the system through sequestration of carbon in their large body and thus strengthening ‘slow’ energy channeling (Rooney and McCann, 2012; Schwarzmüller et al., 2015). Earthworms contribute only about 15.4% to community metabolism in rainforest, leaving vacated trophic space for other groups that vanish or are reduced in plantations. Combined, high Δ¹³C values and shift in dominance of detritivore taxa suggest that the detrital energy channel in rainforest is diversified and comprises a wider range of consumer groups than in plantations, whereas in plantations it comprises almost exclusively Annelida. The similar weighted average Δ¹³C values in plantations and rainforest suggest that from an energetic perspective, soil food webs in plantations are as efficient in processing old organic carbon as in rainforest, despite having a very different structure. However, the decomposer system in plantations may be more vulnerable against future changes since their functioning relies on a single detritivore group comprising predominantly a single invasive species (Potapov et al., 2021b).

Unlike Δ¹³C, changes in Δ¹⁵N with changes in land use were less consistent across animal groups. Lower Δ¹⁵N values in Oribatida and Blattodea in rainforest compared to plantations may be linked to their trophic plasticity, allowing them to shift their trophic position from primary decomposer in rainforest to secondary decomposer in litter-poor plantations (Krause et al., 2019). By contrast, many taxa had lower trophic positions in oil palm compared to rainforest, including Hemiptera, Orthoptera, Isoptera, Pseudoscorpiones, Diplopoda, and Psocoptera. The low Δ¹⁵N values in these taxa at least in part may be due to feeding on resources depleted in Δ¹⁵N relative to litter, such as algae and lichens (Chahartaghi et al., 2005; Schneider et al., 2004), suggesting that non-vascular plants may play a more important role for groups such as Hemiptera, Orthoptera, and Isoptera in plantations than in rainforest, potentially associated with the more open canopies in plantations.

The range of Δ¹⁵N values in soil communities reflects the length of food chains (Cabana and Rasmussen, 1994; Scheu and Falca, 2000), and was the largest in jungle rubber. This was caused by the very low Δ¹⁵N values of Pauropoda (–5.1‰) and Orthoptera (–2.1‰) and high Δ¹⁵N values of Diplura (10.0‰). Jungle rubber is a system that is highly heterogeneous in management practices and plant richness (Gouyon et al., 1993; Rembold et al., 2017b), with species richness in some arthropod predators even exceeding that in rainforest at our study sites (Junggebauer et al., 2021). Anthropogenic disturbances in jungle rubber are moderate compared to monoculture plantation systems (Barnes et al., 2014) and food chains have been found to be longest at intermediate levels of disturbance (Menge and Sutherland, 1987; Polis and Winemiller, 2018; Post, 2002b), which may explain the largest range of Δ¹⁵N values in jungle rubber. As a note of caution, however, the δ¹⁵N values of primary producers (vascular plants, algae, and lichens) may vary among our study systems, which may have affected the Δ¹⁵N values of consumers, but unlikely our overall conclusions.

Changes in functional diversity of soil food webs from community and energetic perspectives

Refuting our second hypothesis, neither isotopic diversity nor isotopic redundancy were higher in rainforest than in plantations. However, isotopic richness was slightly higher in the two more natural systems (i.e., rainforest and jungle rubber) than in rubber and oil palm plantations. Oil palm showed significantly higher community dispersion values than the other land-use systems and in trend had the highest community divergence (unweighted values) reflecting the proportion of groups with the most extreme trophic (isotopic) niches within the community (Cucherousset and Villéger, 2015; Mason et al., 2005; Villéger et al., 2008). At least in part this likely was due to feeding on non-vascular plants, such as algae and lichens, characterized by very different stable isotope values than C3 plants, that is, the dominant vegetation at our study sites (Potapov et al., 2019c). As discussed above, the more open canopy in plantations favors algae and lichens (Drescher et al., 2016; Schulz et al., 2019), together with the monopolization of detrital channel by earthworms (representing another ‘extreme’ isotopic niche), the use of non-vascular plants explains the high dispersion and divergence of energy channeling in oil palm plantations.

From the ‘energetic perspective’, soil food-web divergence in plantations was significantly higher than in rainforest. This contrasts previous evidence that functional divergence decreases with disturbance (Gerisch et al., 2012; Mouillot et al., 2013). However, contrary to divergence, energetic dispersion was in trend higher in rainforest than in the other land-use systems. Similar to the community metrics, the energetic metrics indicated that food-web characteristics in plantations deviate from those in rainforest (high divergence), with food webs being less balanced (low dispersion) with most of the energy being channeled and locked into earthworms.

Community isotopic uniqueness, defined as the inverse of the average isotopic redundancy, and community evenness (Cucherousset and Villéger, 2015) were low in jungle rubber. Most of the soil animal groups clustered in a small region in stable isotope space in jungle rubber, resulting in low community uniqueness, while Pauropoda and Orthoptera were far from this cluster, resulting in low community evenness. High functional redundancy may buffer against land-use impacts and promote stable food webs with long trophic chains (Brodie et al., 2014; Chua et al., 2021; Sanders et al., 2018), which is supported by the results of our study. The increase in isotopic uniqueness (both energetic and unweighted) in oil palm plantations may reflect eroded resilience of this system against future changes.

Overall, community food-web metrics were more sensitive to changes in land use than energetic food-web metrics, suggesting that compositional changes in soil food webs with land use are stronger than changes in energy channeling. This echoes earlier findings that land-use effects on soil animal biodiversity exceed those on functional diversity (Potapov et al., 2020). Results of our SEM indicated significant direct effects of land use-induced environmental changes (i.e., litter amount, tree density, and tree species richness) on community food-web metrics. The amount and quality of leaf litter are important drivers of soil fauna composition and soil food-web structure, being both the food and the habitat for soil animals (Fujii et al., 2020; Sayer et al., 2006). Apart from litter-mediated effects, changes in tree density and species richness are associated with changes in root-derived resources (Ballauff et al., 2021), which also fuel belowground food webs (Pollierer et al., 2007; Bradford, 2016), and this may explain the direct effects of tree communities on the food-web metrics. By contrast, energetic food-web metrics were not directly affected by changes in tree communities and pH, but were linked to the changes in earthworm abundance. High soil pH favors colonization of plantations by earthworms and this is common in the tropics (Marichal et al., 2010; Potapov et al., 2021b). The close association between energetic food-web metrics and the fraction earthworms contribute to community metabolism stems to a large extent from the mathematical dependence between these two variables. However, we intentionally wanted to illustrate that those strong shifts in the functional diversity of food webs may result from a single group benefiting from certain environmental changes.

In conclusion, our study is among the first comprehensive assessment of tropical soil food webs and their variation due to land-use changes. Low Δ¹³C values in most soil animal groups in plantations in comparison to rainforest indicate a shift toward using plant carbon and ‘fast’ energy channeling based on high-quality understory plants (weeds) as well as algae. On the other hand, the trophic niche of earthworms as major macrofauna detritivores stayed unchanged and they monopolized the ‘slow’ detrital channel in plantations. This resulted in systems with strong divergence and imbalance in energetic pathways potentially compromising functional stability of plantation systems. Other studies at our sites showed that these changes in soil food-web characteristics with transformation of rainforest into plantations are associated with reduced soil functioning (Grass et al., 2020) and litter invertebrate biodiversity (Barnes et al., 2014). Our analyses allowed to uncover the mechanisms responsible for these changes and demonstrated that land-use effects on soil biodiversity from a ‘community perspective’ are in part buffered from the perspective of energy channeling (‘energetic perspective’), but resistance of plantations against future changes in climate and land use may be compromised.

Materials and methods

Sampling sites

Request a detailed protocol

The study was conducted in the framework of the collaborative research project CRC990/EFForTS investigating ecological and socio‐economic changes associated with the transformation of lowland rainforest into agricultural systems (Drescher et al., 2016). Four land-use systems, rainforest, jungle rubber, rubber plantations, and oil palm plantations were investigated in two regions, that is, Harapan and Bukit Duabelas (Drescher et al., 2016). Jungle rubber sites were established by planting rubber trees (Hevea brasiliensis) into selectively logged rainforest and contain rainforest tree species. Jungle rubber sites represent low intensive land-use systems, lacking fertilizer input as well as herbicide application; the age of rubber trees varied between 15 and 40 years (Kotowska et al., 2015). Rubber and oil palm (Elaeis guineensis) monocultures represent high land-use intensity plantation systems managed by the addition of fertilizers as well as herbicides (Drescher et al., 2016). Each land-use system was replicated four times in each landscape, resulting in a total of 32 sites; for more details, see Drescher et al., 2016.

Sampling, extraction, and classification of soil fauna

Request a detailed protocol

Soil animals were sampled at each of 32 study sites during October and November 2013. Soil samples measuring 16 cm × 16 cm and including the litter layer and 0–5 cm of the mineral soil were taken in three 5 m × 5 m subplots within each of 50 m × 50 m plots established at each study site, resulting in a total of 96 samples. The samples were transported to the laboratory and animals were extracted by heat (Kempson et al., 1963) until the substrate was completely dry (6–8 days). Until further analysis, species were stored in 70% ethanol. For calibration of the animal stable isotope values, we used mixed litter samples that were taken from each site and analyzed in a previous study (Klarner et al., 2017).

Animals were classified into 23 high-rank taxonomic groups (Oribatida, Collembola, Symphyla, Protura, Annelida, Blattodea, Diplopoda, Isopoda, Isoptera, Psocoptera, Psocoptera, Lepidoptera, Orthoptera, Thysanoptera, Diptera, Coleoptera, Pseudoscorpiones, Mesostigmata, Diplura, Formicidae, Chilopoda, Araneae, Pauropoda). For stable isotope analysis, we adopted a group-level analysis representing the stable isotope niche at the level of taxonomic groups. Although this approach may miss the variability in stable isotope niches of species within high-rank taxonomic groups, it has the advantage that it integrates across species allowing generalizations on the trophic structure and energy flux of whole communities. The approach has been recently advocated for analyzing the channeling of energy through food webs using lipid profiling (Kühn et al., 2018), but has not been adopted yet in stable isotope analysis although it has been shown that at least in soil high-rank animal taxa typically represent the trophic niches of species in most taxa (Potapov et al., 2019b). To ensure that our samples reliably represent the trophic niche of the studied taxa, we included (if ever possible) 15 individuals per taxon per study site. Doing that we considered the turnover of species among sites and focused on dominant species representing the majority of biomass, which we considered most important for our approach. We further classified taxonomic groups into five major functional groups according to their trophic guild and body size class (Potapov et al., 2019b; Potapov et al., 2021a): herbivores including, for example, Hemiptera and Orthoptera, microdecomposers including, for example, Oribatida and Collembola, macrodecomposers including, for example, Annelida and Diplopoda, micropredators including, for example, Diplura and Mesotigmata, macropredators including, for example, Araneae and Chilopoda, and groups with mixed feeding habits including, for example, Diptera and Coleoptera.

Stable isotope analysis

Request a detailed protocol

To cover the entire community, for each sampling site we analyzed all taxa for which we were able to collect enough biomass for stable isotope analysis and which were represented by more than two individuals. We analyzed a minimum of 3 and a maximum of 15 individuals for each taxonomic group for each site as a single mixed sample to cover the species- and individual-level isotopic variation. We mixed individuals from different subplots whenever possible to cover spatial variation in stable isotope values. Animals from the litter and soil layer were analyzed separately, but were merged for data analysis since stable isotope values did not differ significantly between layers. Animal samples were dried at 60°C for 24 hr before stable isotope analysis, weighed and wrapped into tin capsules; sample weight varied between 0.01 and 1.00 mg. For small-sized animal groups we used bulk individuals, for large-sized animal groups we used body parts dominated by muscle tissue (e.g., legs) from different individuals and pooled them (Tsurikov et al., 2015). In total, 626 samples of 23 taxonomic groups were analyzed across 32 sites. For Collembola, Oribatida, and Chilopoda, we additionally used stable isotope data collected at species level (Klarner et al., 2017; Krause et al., 2019; Susanti et al., 2021) to calculate a single average value for each group at each site, which were collected at the same sampling campaign. The number of analyzed taxonomic groups varied between 6 and 17 per site (i.e., per one soil food web) and was on average 12.3.

Animal samples were analyzed using a coupled system of an elemental analyzer (NA 1500, Carlo Erba, Milan, Italy) and a mass spectrometer (MAT 251, Finnigan, Bremen, Germany) adopted for the analysis of small sample sizes (Langel and Dyckmans, 2014). Ratios of the heavy isotope to the light isotope (¹³C/¹²C,¹⁵N/¹⁴N, denoted as R) were expressed in parts per thousand relative to the standard using the delta notation with δ¹³C or δ¹⁵N = (R_sample/R_standard − 1) × 1000 (‰). Vienna PD Belemnite and atmospheric nitrogen were used as standard for ¹³C and ¹⁵N, respectively. Acetanilid was used for internal calibration.

Environmental parameters of the study sites were used as given in Potapov et al., 2020, Krashevska et al., 2015, and Rembold et al., 2017a, which included tree species richness, tree density, understory species richness, understory density, and average understory height, soil pH, litter amount, soil carbon concentration, carbon-to-nitrogen ratio of litter, soil microbial biomass C, and soil humidity.

Statistical analysis

Request a detailed protocol

The stable isotope compositions of animals were calibrated to that of the local leaf litter. Calibrated δ¹³C and δ¹⁵N values were calculated as the difference between the plot-specific litter δ¹³C and δ¹⁵N values and the δ¹³C and δ¹⁵N values of each group, and given as Δ¹³C and Δ¹⁵N values, respectively. Statistical analyses were done in R v 4.0 (R Development Core Team, 2020) with R studio interface (RStudio Team, 2020).

To characterize the trophic structure of soil animal communities, we calculated isotopic metrics as given in Cucherousset and Villéger, 2015. One-dimensional metrics describe the isotopic parameters of the communities based on Δ¹⁵N or Δ¹³C values. Multidimensional metrics combine both Δ¹³C and Δ¹⁵N values, and join the ones from Layman et al., 2007, with functional diversity framework (Villéger et al., 2008; Laliberté and Legendre, 2010). The Δ¹³C and Δ¹⁵N values were scaled between 0 and 1 based on maximum and minimum across all communities to ensure equal contribution of two isotopes prior to calculation of multidimensional metrics. Multidimensional metrics were calculated from two perspectives: (1) a ‘community perspective’, assuming all taxonomic groups being equally important, that is, unweighted metrics, and (2) an ‘energetic perspective’, assuming that groups that have higher contribution to total community metabolism are also more functionally important, that is, metrics were weighted by community metabolism. We used metabolism instead of biomass because it better reflects the contribution of organisms to energy processing and thus their importance in the food web (Brown et al., 2004; Barnes et al., 2018). Community metabolism for each group at each plot was taken from Potapov et al., 2019a; it was based on length and width measurements of all individuals and using body size to body mass ratios and group-specific allometric regressions to calculate metabolic rates (Ehnes et al., 2011). Individual metabolic rates were then summed up for groups to estimate contribution of each taxonomic group to the total community metabolism per plot (Supplementary file 1, Figure 1—figure supplement 1).

Overall, 13 isotopic metrics were calculated for each of 32 communities (i.e., sampling plots). One-dimensional metrics included average position, range, minimum, and maximum. The unweighted and metabolism-weighted average position of communities (mean isotopic value across groups) represent mean community-level isotopic trait values. The isotopic range represents the difference between minimum and maximum values of both Δ¹³C and Δ¹⁵N. Range, minimum, and maximum could not be weighted and are given unweighted. Multidimensional metrics included isotopic divergence, isotopic dispersion, isotopic evenness, isotopic uniqueness, and isotopic richness, which were calculated as both unweighted and metabolism-weighed. Isotopic divergence represents the distance between all species and the center of the convex hull area. Isotopic divergence values close to 0 indicate that groups with extreme stable isotope values are rare (community divergence) or contribute little to the community metabolism (energetic divergence), whereas isotopic divergence values close to 1 indicate that there are many groups with extreme stable isotope values (community divergence) or they contribute considerably to the community metabolism (energetic divergence). Isotopic dispersion combines convex hull area with isotopic divergence values and can be interpreted as scaled multidimensional variance. Isotopic dispersion approaches 1 when species with contrasting stable isotope values have similar abundance, which is a more functionally diverse and balanced system, whereas it approaches 0 when most groups (community dispersion) or community metabolism (energetic dispersion) are concentrated near the ‘center of gravity’ of the community in stable isotope space. Isotopic evenness quantifies the distribution of groups or metabolism in stable isotope space. Isotopic evenness values close to 1 indicate that the isotope values of the groups/metabolism are evenly distributed, while values close to 0 indicate that the groups/metabolism cluster together. Isotopic uniqueness reflects the closeness of stable isotope values of the studied groups/metabolism within the community, which is defined as the inverse of the average isotopic redundancy. Finally, isotopic richness is the volume occupied by all groups in isotopic space (convex hull area in two-dimensional isotopic space) and reflects functional richness of the food web; it is the only multidimensional metric that cannot be weighted since it considers the total isotopic space (Mason et al., 2005; Villéger et al., 2008).

To assess differences in food-web structure among land-use systems, we used a set of analyses of variance (aov function) with the Δ¹³C and Δ¹⁵N values of each taxonomic group, one-dimensional isotopic metrics, and multidimensional community and energetic isotopic metrics as response variables, and land-use system (rainforest, jungle rubber, rubber, oil palm) and landscape (Harapan or Bukit Duabelas) as factors (total n = 32, 8 plots as replicates per land-use system). Pairwise comparisons of means among land-use systems were done using post hoc HSD.test function from the package agricolae Margur, 2020 following analyses of variance. Differences in Δ¹³C and Δ¹⁵N values between rainforest and other land-use systems for each taxonomic group were analyzed with Student’s t.test function in R. Results were visualized using the ggplot2 package (Wickham, 2016).

To assess effect size of land use on all food-web metrics combined, we used analysis of similarities based on community and energetic metrics with land use as the grouping variable (anosim in package vegan). Besides, we used multivariate analyses of variance (MANOVAs) to inspect the effects of environmental factors on community and energetic metrics, and additionally explored pairwise correlations between environmental factors and food-web metrics using Spearman’s correlation from the package agricolae (Margur, 2020).

Finally, a SEM based on generalized least squares was constructed to provide insight into how land use affected soil food webs from both community and energetic perspectives. The analysis was performed with the lavaan package in R (Rosseel, 2012). The model included tree, understory, and soil properties selected according to permutation tests based on R² which were used to quantify the land-use effects (ordiR2stepin package vegan), and before permutation tests, the environmental factors were filtered based on the MANOVAs and Spearman’s correlation. The final model included soil pH, tree density, and tree richness as the three most important variables that represented direct land-use effects (i.e., logging and liming; Drescher et al., 2016). Furthermore, we included litter amount, understory density, and earthworm metabolism as the three mediators that are affected by changes in tree density, tree richness, and pH, and have strong impacts on tropical soil invertebrate communities (Darras et al., 2019; Potapov et al., 2019a). Food-web metrics (i.e., isotopic divergence, dispersion, uniqueness, evenness, and average position) were first combined using principal component analysis (prcomp function) and the PC1 was used as the response variable in SEM; PC1 explained 66% and 42% of the variance for community and energetic metrics, respectively. To determine the goodness of fit of the model, we used χ²-test associated p-value ≥ 0.05, the comparative fit index (CFI) > 0.95, the root-mean-square error of approximation (RMSEA) and the standardized root-mean-square residual (SRMR) with values ≤ 0.05 (Schermelleh-Engel et al., 2003). Our SEM adequately described the data (χ² = 6.23, p = 0.40, df = 5, CFI = 0.98, RMSEA = 0.04, SRMR = 0.05).

Data availability

All data generated or analysed during this study are included in the manuscript and supporting file.

References

(1996) Isotope Discrimination during Decomposition of Organic Matter: A Theoretical Analysis
Soil Science Society of America Journal 60:1121–1126.

https://doi.org/10.2136/sssaj1996.03615995006000040023x
- Google Scholar
1. Ballauff J
2. Schneider D
3. Edy N
4. Irawan B
5. Daniel R
6. Polle A
(2021) Shifts in root and soil chemistry drive the assembly of belowground fungal communities in tropical land-use systems
Soil Biology and Biochemistry 154:108140.

https://doi.org/10.1016/j.soilbio.2021.108140
- Google Scholar
Book
1. Bardgett RD
2. Wardle DA
(2010)
Aboveground-Belowground Linkages: Biotic Interactions, Ecosystem Processes, and Global Change. Oxford Series in Ecology and Evolution

Oxford: Oxford University Press.
- Google Scholar
1. Bardgett RD
2. Putten WH
(2014) Belowground biodiversity and ecosystem functioning
Nature 515:505–511.

https://doi.org/10.1038/nature13855
- Google Scholar
1. Barnes AD
2. Jochum M
3. Mumme S
4. Haneda NF
5. Farajallah A
6. Widarto TH
7. Brose U
(2014) Consequences of tropical land use for multitrophic biodiversity and ecosystem functioning
Nature Communications 5:5351.

https://doi.org/10.1038/ncomms6351
- PubMed
- Google Scholar
1. Barnes AD
2. Jochum M
3. Lefcheck JS
4. Eisenhauer N
5. Scherber C
6. O’Connor MI
7. Ruiter P
8. Brose U
(2018) Energy Flux: The Link between Multitrophic Biodiversity and Ecosystem Functioning
Trends Ecol Evol 33:186–197.

https://doi.org/10.1016/j.tree.2017.12.007
- Google Scholar
1. Blouin M
2. Hodson ME
3. Delgado EA
4. Baker G
5. Brussaard L
6. Butt KR
7. Dai J
8. Dendooven L
9. Peres G
10. Tondoh JE
11. Cluzeau D
12. Brun JJ
(2013) A review of earthworm impact on soil function and ecosystem services: Earthworm impact on ecosystem services
European Journal of Soil Science 64:161–182.

https://doi.org/10.1111/ejss.12025
- Google Scholar
(2019) Protura are unique: first evidence of specialized feeding on ectomycorrhizal fungi in soil invertebrates
BMC Ecology 19:10.

https://doi.org/10.1186/s12898-019-0227-y
- PubMed
- Google Scholar
(2013) Ecological intensification: harnessing ecosystem services for food security
Trends in Ecology & Evolution 28:230–238.

https://doi.org/10.1016/j.tree.2012.10.012
- PubMed
- Google Scholar
(2007) Isotope fractionation and 13C enrichment in soil profiles during the decomposition of soil organic matter
Oecologia 153:89–98.

https://doi.org/10.1007/s00442-007-0700-8
- PubMed
- Google Scholar
1. Bradford MA
(2016) Re-visioning soil food webs
Soil Biology and Biochemistry 102:1–3.

https://doi.org/10.1016/j.soilbio.2016.08.010
- Google Scholar
1. Brodie JF
2. Aslan CE
3. Rogers HS
4. Redford KH
5. Maron JL
6. Bronstein JL
7. Groves CR
(2014) Secondary extinctions of biodiversity
Trends in Ecology & Evolution 29:664–672.

https://doi.org/10.1016/j.tree.2014.09.012
- PubMed
- Google Scholar
1. Brose U
2. Scheu S
(2014) Into darkness: unravelling the structure of soil food webs
Oikos 123:1153–1156.

https://doi.org/10.1111/oik.01768
- Google Scholar
1. Brown JH
2. Gillooly JF
3. Allen AP
4. Savage VM
5. West GB
(2004) TOWARD A METABOLIC THEORY OF ECOLOGY
Ecology 85:1771–1789.

https://doi.org/10.1890/03-9000
- Google Scholar
1. Cabana G
2. Rasmussen JB
(1994) Modelling food chain structure and contaminant bioaccumulation using stable nitrogen isotopes
Nature 372:255–257.

https://doi.org/10.1038/372255a0
- Google Scholar
(2005) Feeding guilds in Collembola based on nitrogen stable isotope ratios
Soil Biology and Biochemistry 37:1718–1725.

https://doi.org/10.1016/j.soilbio.2005.02.006
- Google Scholar
1. Chua KWJ
2. Liew JH
3. Wilkinson CL
4. Ahmad AB
5. Tan HH
6. Yeo DCJ
(2021) Land‐use change erodes trophic redundancy in tropical forest streams: Evidence from amino acid stable isotope analysis
The Journal of Animal Ecology 1:1365–2656.

https://doi.org/10.1111/1365-2656.13462
- Google Scholar
1. Clough Y
2. Krishna VV
3. Corre MD
4. Darras K
5. Denmead LH
6. Meijide A
7. Moser S
8. Musshoff O
9. Steinebach S
10. Veldkamp E
11. Allen K
12. Barnes AD
13. Breidenbach N
14. Brose U
15. Buchori D
16. Daniel R
17. Finkeldey R
18. Harahap I
19. Hertel D
20. Holtkamp AM
21. Hörandl E
22. Irawan B
23. Jaya INS
24. Jochum M
25. Klarner B
26. Knohl A
27. Kotowska MM
28. Krashevska V
29. Kreft H
30. Kurniawan S
31. Leuschner C
32. Maraun M
33. Melati DN
34. Opfermann N
35. Pérez-Cruzado C
36. Prabowo WE
37. Rembold K
38. Rizali A
39. Rubiana R
40. Schneider D
41. Tjitrosoedirdjo SS
42. Tjoa A
43. Tscharntke T
44. Scheu S
(2016) Land-use choices follow profitability at the expense of ecological functions in Indonesian smallholder landscapes
Nature Communications 7:13137.

https://doi.org/10.1038/ncomms13137
- PubMed
- Google Scholar
Book
(1983) Biological Strategies of Nutrient Cycling in Soil Systems
In: Coleman DC, editors. Advances in Ecological Research. Amsterdam: Elsevier. pp. 1–55.

https://doi.org/10.1016/S0065-2504(08)60107-5
- Google Scholar
1. Cucherousset J
2. Villéger S
(2015) Quantifying the multiple facets of isotopic diversity: New metrics for stable isotope ecology
Ecological Indicators 56:152–160.

https://doi.org/10.1016/j.ecolind.2015.03.032
- Google Scholar
1. Darras KFA
2. Corre MD
3. Formaglio G
4. Tjoa A
5. Potapov A
6. Brambach F
7. Sibhatu KT
8. Grass I
9. Rubiano AA
10. Buchori D
11. Drescher J
12. Fardiansah R
13. Hölscher D
14. Irawan B
15. Kneib T
16. Krashevska V
17. Krause A
18. Kreft H
19. Li K
20. Maraun M
21. Polle A
22. Ryadin AR
23. Rembold K
24. Stiegler C
25. Scheu S
26. Tarigan S
27. Valdés-Uribe A
28. Yadi S
29. Tscharntke T
30. Veldkamp E
(2019) Reducing Fertilizer and Avoiding Herbicides in Oil Palm Plantations—Ecological and Economic Valuations
Frontiers in Forests and Global Change 2:65.

https://doi.org/10.3389/ffgc.2019.00065
- Google Scholar
(2006) Fungal/bacterial ratios in grasslands with contrasting nitrogen management
Soil Biology and Biochemistry 38:2092–2103.

https://doi.org/10.1016/j.soilbio.2006.01.008
- Google Scholar
(2012) Land use alters the resistance and resilience of soil food webs to drought
Nature Climate Change 2:276–280.

https://doi.org/10.1038/nclimate1368
- Google Scholar
(2013) Soil food web properties explain ecosystem services across European land use systems
PNAS 110:14296–14301.

https://doi.org/10.1073/pnas.1305198110
- PubMed
- Google Scholar
(2011) Agricultural expansion and the fate of global conservation priorities
Biodiversity and Conservation 20:2445–2459.

https://doi.org/10.1007/s10531-011-9997-z
- Google Scholar
1. Drescher J
2. Rembold K
3. Allen K
4. Beckschäfer P
5. Buchori D
6. Clough Y
7. Faust H
8. Fauzi AM
9. Gunawan D
10. Hertel D
11. Irawan B
12. Jaya INS
13. Klarner B
14. Kleinn C
15. Knohl A
16. Kotowska MM
17. Krashevska V
18. Krishna V
19. Leuschner C
20. Lorenz W
21. Meijide A
22. Melati D
23. Nomura M
24. Pérez-Cruzado C
25. Qaim M
26. Siregar IZ
27. Steinebach S
28. Tjoa A
29. Tscharntke T
30. Wick B
31. Wiegand K
32. Kreft H
33. Scheu S
(2016) Ecological and socio-economic functions across tropical land use systems after rainforest conversion
Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 371:20150275.

https://doi.org/10.1098/rstb.2015.0275
- PubMed
- Google Scholar
1. Ehnes RB
2. Rall BC
3. Brose U
(2011) Phylogenetic grouping, curvature and metabolic scaling in terrestrial invertebrates
Ecology Letters 14:993–1000.

https://doi.org/10.1111/j.1461-0248.2011.01660.x
- PubMed
- Google Scholar
(2020) Living Litter: Dynamic Trait Spectra Predict Fauna Composition
Trends in Ecology & Evolution 35:886–896.

https://doi.org/10.1016/j.tree.2020.05.007
- PubMed
- Google Scholar
(2021) Radiocarbon signature reveals that most springtails depend on carbon from living plants
Biology Letters 17:20210353.

https://doi.org/10.1098/rsbl.2021.0353
- PubMed
- Google Scholar
1. Geisen S
2. Briones MJI
3. Gan H
4. Behan-Pelletier VM
5. Friman VP
6. de Groot GA
7. Hannula SE
8. Lindo Z
9. Philippot L
10. Tiunov AV
11. Wall DH
(2019) A methodological framework to embrace soil biodiversity
Soil Biology and Biochemistry 136:107536.

https://doi.org/10.1016/j.soilbio.2019.107536
- Google Scholar
(2012) More species, but all do the same: contrasting effects of flood disturbance on ground beetle functional and species diversity
Oikos 121:508–515.

https://doi.org/10.1111/j.1600-0706.2011.19749.x
- Google Scholar
1. Gessner MO
2. Swan CM
3. Dang CK
4. McKie BG
5. Bardgett RD
6. Wall DH
7. Hättenschwiler S
(2010) Diversity meets decomposition
Trends in Ecology & Evolution 25:372–380.

https://doi.org/10.1016/j.tree.2010.01.010
- PubMed
- Google Scholar
(1993) Does ‘jungle rubber’ deserve its name? An analysis of rubber agroforestry systems in southeast Sumatra
Agroforestry Systems 22:181–206.

https://doi.org/10.1007/BF00705233
- Google Scholar
1. Grass I
2. Kubitza C
3. Krishna VV
4. Corre MD
5. Mußhoff O
6. Pütz P
7. Drescher J
8. Rembold K
9. Ariyanti ES
10. Barnes AD
11. Brinkmann N
12. Brose U
13. Brümmer B
14. Buchori D
15. Daniel R
16. Darras KFA
17. Faust H
18. Fehrmann L
19. Hein J
20. Hennings N
21. Hidayat P
22. Hölscher D
23. Jochum M
24. Knohl A
25. Kotowska MM
26. Krashevska V
27. Kreft H
28. Leuschner C
29. Lobite NJS
30. Panjaitan R
31. Polle A
32. Potapov AM
33. Purnama E
34. Qaim M
35. Röll A
36. Scheu S
37. Schneider D
38. Tjoa A
39. Tscharntke T
40. Veldkamp E
41. Wollni M
(2020) Trade-offs between multifunctionality and profit in tropical smallholder landscapes
Nature Communications 11:1186.

https://doi.org/10.1038/s41467-020-15013-5
- PubMed
- Google Scholar
1. Guerra CA
2. Bardgett RD
3. Caon L
4. Crowther TW
5. Delgado-Baquerizo M
6. Montanarella L
7. Navarro LM
8. Orgiazzi A
9. Singh BK
10. Tedersoo L
11. Vargas-Rojas R
12. Briones MJI
13. Buscot F
14. Cameron EK
15. Cesarz S
16. Chatzinotas A
17. Cowan DA
18. Djukic I
19. van den Hoogen J
20. Lehmann A
21. Maestre FT
22. Marín C
23. Reitz T
24. Rillig MC
25. Smith LC
26. de Vries FT
27. Weigelt A
28. Wall DH
29. Eisenhauer N
(2021) Tracking, targeting, and conserving soil biodiversity
Science (New York, N.Y.) 371:239–241.

https://doi.org/10.1126/science.abd7926
- PubMed
- Google Scholar
(2018) Carbon costs and benefits of Indonesian rainforest conversion to plantations
Nature Communications 9:2388.

https://doi.org/10.1038/s41467-018-04755-y
- PubMed
- Google Scholar
1. Hooper DU
2. Bignell DE
3. Brown VK
4. Brussard L
5. Mark Dangerfield J
6. Wall DH
7. Wardle DA
8. Coleman DC
9. Giller KE
10. Lavelle P
11. Van Der Putten WH
12. De Ruiter PC
13. Rusek J
14. Silver WL
15. Tiedje JM
16. Wolters V
(2000) Interactions between Aboveground and Belowground Biodiversity in Terrestrial Ecosystems: Patterns, Mechanisms, and Feedbacks
Bioscience 50:1049.

https://doi.org/10.1641/0006-3568(2000)050[1049:IBAABB]2.0.CO;2
- Google Scholar
1. Hyodo F
2. Matsumoto T
3. Takematsu Y
4. Kamoi T
5. Fukuda D
6. Nakagawa M
7. Itioka T
(2010) The structure of a food web in a tropical rain forest in Malaysia based on carbon and nitrogen stable isotope ratios
Journal of Tropical Ecology 26:205–214.

https://doi.org/10.1017/S0266467409990502
- Google Scholar
1. Hyodo F
2. Uchida T
3. Kaneko N
4. Tayasu I
(2012) Use of radiocarbon to estimate diet ages of earthworms across different climate regions
Applied Soil Ecology 62:178–183.

https://doi.org/10.1016/j.apsoil.2012.09.014
- Google Scholar
1. Hyodo F
(2015) Use of stable carbon and nitrogen isotopes in insect trophic ecology: Use of isotope in insect trophic ecology
Entomol Sci 18:295–312.

https://doi.org/10.1111/ens.12128
- Google Scholar
1. Illig J
2. Langel R
3. Norton RA
4. Scheu S
5. Maraun M
(2005) Where are the decomposers? Uncovering the soil food web of a tropical montane rain forest in southern Ecuador using stable isotopes (¹⁵ N)
Journal of Tropical Ecology 21:589–593.

https://doi.org/10.1017/S0266467405002646
- Google Scholar
1. Jochum M
2. Barnes AD
3. Brose U
4. Gauzens B
5. Sünnemann M
6. Amyntas A
7. Eisenhauer N
(2021) For flux’s sake: General considerations for energy-flux calculations in ecological communities
Ecology and Evolution 11:12948–12969.

https://doi.org/10.1002/ece3.8060
- PubMed
- Google Scholar
1. Junggebauer A
2. Hartke TR
3. Ramos D
4. Schaefer I
5. Buchori D
6. Hidayat P
7. Scheu S
8. Drescher J
(2021) Changes in diversity and community assembly of jumping spiders (Araneae: Salticidae) after rainforest conversion to rubber and oil palm plantations
PeerJ 9:e11012.

https://doi.org/10.7717/peerj.11012
- PubMed
- Google Scholar
(1963)
A new extractor for woodland litter

Pedobiologia 3:1–21.
- Google Scholar
(2017) Trophic niches, diversity and community composition of invertebrate top predators (Chilopoda) as affected by conversion of tropical lowland rainforest in Sumatra (Indonesia
PLOS ONE 12:e0180915.

https://doi.org/10.1371/journal.pone.0180915
- PubMed
- Google Scholar
1. Koh LP
2. Ghazoul J
(2010) Spatially explicit scenario analysis for reconciling agricultural expansion, forest protection, and carbon conservation in Indonesia
PNAS 107:11140–11144.

https://doi.org/10.1073/pnas.1000530107
- PubMed
- Google Scholar
(2015) Quantifying above- and belowground biomass carbon loss with forest conversion in tropical lowlands of Sumatra (Indonesia
Global Change Biology 21:3620–3634.

https://doi.org/10.1111/gcb.12979
- PubMed
- Google Scholar
(2015) Impact of tropical lowland rainforest conversion into rubber and oil palm plantations on soil microbial communities
Biology and Fertility of Soils 51:697–705.

https://doi.org/10.1007/s00374-015-1021-4
- Google Scholar
1. Krashevska V
2. Malysheva E
3. Klarner B
4. Mazei Y
5. Maraun M
6. Widyastuti R
7. Scheu S
(2018) Micro-decomposer communities and decomposition processes in tropical lowlands as affected by land use and litter type
Oecologia 187:255–266.

https://doi.org/10.1007/s00442-018-4103-9
- PubMed
- Google Scholar
(2019) Changes in Nematode Communities and Functional Diversity With the Conversion of Rainforest Into Rubber and Oil Palm Plantations
Frontiers in Ecology and Evolution 7:487.

https://doi.org/10.3389/fevo.2019.00487
- Google Scholar
1. Krause A
2. Sandmann D
3. Bluhm SL
4. Ermilov S
5. Widyastuti R
6. Haneda NF
7. Scheu S
8. Maraun M
(2019) Shift in trophic niches of soil microarthropods with conversion of tropical rainforest into plantations as indicated by stable isotopes (15N, 13C
PLOS ONE 14:e0224520.

https://doi.org/10.1371/journal.pone.0224520
- PubMed
- Google Scholar
1. Krause A
2. Sandmann D
3. Potapov A
4. Ermilov S
5. Widyastuti R
6. Haneda NF
7. Scheu S
8. Maraun M
(2021) Variation in Community-Level Trophic Niches of Soil Microarthropods With Conversion of Tropical Rainforest Into Plantation Systems as Indicated by Stable Isotopes (15N, 13C
Frontiers in Ecology and Evolution 9:e92149.

https://doi.org/10.3389/fevo.2021.592149
- Google Scholar
1. Kühn J
2. Richter A
3. Kahl T
4. Bauhus J
5. Schöning I
6. Ruess L
7. Davey M
(2018) Community level lipid profiling of consumers as a tool for soil food web diagnostics
Methods in Ecology and Evolution 9:1265–1275.

https://doi.org/10.1111/2041-210X.12966
- Google Scholar
1. Laliberté E
2. Legendre P
(2010) A distance-based framework for measuring functional diversity from multiple traits
Ecology 91:299–305.

https://doi.org/10.1890/08-2244.1
- PubMed
- Google Scholar
1. Langel R
2. Dyckmans J
(2014) Combined 13C and 15N isotope analysis on small samples using a near-conventional elemental analyzer/isotope ratio mass spectrometer setup
Rapid Communications in Mass Spectrometry 28:1019–1022.

https://doi.org/10.1002/rcm.6878
- PubMed
- Google Scholar
1. Laurance WF
(2007) Have we overstated the tropical biodiversity crisis?
Trends in Ecology & Evolution 22:65–70.

https://doi.org/10.1016/j.tree.2006.09.014
- PubMed
- Google Scholar
(2007) Can stable isotope ratios provide for community-wide measures of trophic structure?
Ecology 88:42–48.

https://doi.org/10.1890/0012-9658(2007)88[42:csirpf]2.0.co;2
- PubMed
- Google Scholar
1. Layman CA
2. Araujo MS
3. Boucek R
4. Hammerschlag-Peyer CM
5. Harrison E
6. Jud ZR
7. Matich P
8. Rosenblatt AE
9. Vaudo JJ
10. Yeager LA
11. Post DM
12. Bearhop S
(2012) Applying stable isotopes to examine food-web structure: an overview of analytical tools
Biological Reviews of the Cambridge Philosophical Society 87:545–562.

https://doi.org/10.1111/j.1469-185X.2011.00208.x
- PubMed
- Google Scholar
1. Liebke DF
2. Harms D
3. Widyastuti R
4. Scheu S
5. Potapov AM
(2021) Impact of rainforest conversion into monoculture plantation systems on pseudoscorpion density, diversity and trophic niches
Soil Organisms 93:83–95.

https://doi.org/10.25674/SO93ISS2ID147
- Google Scholar
1. Margono BA
2. Turubanova S
3. Zhuravleva I
4. Potapov P
5. Tyukavina A
6. Baccini A
7. Goetz S
8. Hansen MC
(2012) Mapping and monitoring deforestation and forest degradation in Sumatra (Indonesia) using Landsat time series data sets from 1990 to 2010
Environmental Research Letters 7:034010.

https://doi.org/10.1088/1748-9326/7/3/034010
- Google Scholar
1. Margur F
(2020)
Statistical Procedures for Agricultural Research

Agricolae 1:1–704.
- Google Scholar
1. Marichal R
2. Martinez AF
3. Praxedes C
4. Ruiz D
5. Carvajal AF
6. Oszwald J
7. del Pilar Hurtado M
8. Brown GG
9. Grimaldi M
10. Desjardins T
11. Sarrazin M
12. Decaëns T
13. Velasquez E
14. Lavelle P
(2010) Invasion of Pontoscolex corethrurus (Glossoscolecidae, Oligochaeta) in landscapes of the Amazonian deforestation arc
Applied Soil Ecology 46:443–449.

https://doi.org/10.1016/j.apsoil.2010.09.001
- Google Scholar
1. Mason NWH
2. Mouillot D
3. Lee WG
4. Wilson JB
(2005) Functional richness, functional evenness and functional divergence: the primary components of functional diversity
Oikos 111:112–118.

https://doi.org/10.1111/j.0030-1299.2005.13886.x
- Google Scholar
1. Matson PA
2. Parton WJ
3. Power AG
4. Swift MJ
(1997) Agricultural intensification and ecosystem properties
Science (New York, N.Y.) 277:504–509.

https://doi.org/10.1126/science.277.5325.504
- PubMed
- Google Scholar
(1998) Weak trophic interactions and the balance of nature
Nature 395:794–798.

https://doi.org/10.1038/27427
- Google Scholar
(2001) Effects of Land-Use Change on Soil Nutrient Dynamics in Amazônia
Ecosystems 4:625–645.

https://doi.org/10.1007/s10021-001-0033-0
- Google Scholar
1. Menge BA
2. Sutherland JP
(1987) Community Regulation: Variation in Disturbance, Competition, and Predation in Relation to Environmental Stress and Recruitment
The American Naturalist 130:730–757.

https://doi.org/10.1086/284741
- Google Scholar
(2011) DEFORESTATION rates IN INSULAR SOUTHEAST ASIA between 2000 and 2010: DEFORESTATION IN INSULAR SOUTHEAST ASIA 2000-2010
Glob Change Biol 17:2261–2270.

https://doi.org/10.1111/j.1365-2486.2011.02398.x
- Google Scholar
(2005) Modeling trophic pathways, nutrient cycling, and dynamic stability in soils
Pedobiologia 49:499–510.

https://doi.org/10.1016/j.pedobi.2005.05.008
- Google Scholar
(2013) A functional approach reveals community responses to disturbances
Trends in Ecology & Evolution 28:167–177.

https://doi.org/10.1016/j.tree.2012.10.004
- PubMed
- Google Scholar
1. Newbold T
2. Hudson LN
3. Hill SLL
4. Contu S
5. Lysenko I
6. Senior RA
7. Börger L
8. Bennett DJ
9. Choimes A
10. Collen B
11. Day J
12. De Palma A
13. Díaz S
14. Echeverria-Londoño S
15. Edgar MJ
16. Feldman A
17. Garon M
18. Harrison MLK
19. Alhusseini T
20. Ingram DJ
21. Itescu Y
22. Kattge J
23. Kemp V
24. Kirkpatrick L
25. Kleyer M
26. Correia DLP
27. Martin CD
28. Meiri S
29. Novosolov M
30. Pan Y
31. Phillips HRP
32. Purves DW
33. Robinson A
34. Simpson J
35. Tuck SL
36. Weiher E
37. White HJ
38. Ewers RM
39. Mace GM
40. Scharlemann JPW
41. Purvis A
(2015) Global effects of land use on local terrestrial biodiversity
Nature 520:45–50.

https://doi.org/10.1038/nature14324
- PubMed
- Google Scholar
(2010) Source partitioning using stable isotopes: coping with too much variation
PLOS ONE 5:e9672.

https://doi.org/10.1371/journal.pone.0009672
- PubMed
- Google Scholar
1. Petchey OL
2. Gaston KJ
(2006) Functional diversity: back to basics and looking forward
Ecology Letters 9:741–758.

https://doi.org/10.1111/j.1461-0248.2006.00924.x
- PubMed
- Google Scholar
1. Peterson BJ
2. Fry B
(1987) STABLE ISOTOPES IN ECOSYSTEM STUDIES
Annual Review of Ecology and Systematics 18:293–320.

https://doi.org/10.1146/annurev.es.18.110187.001453
- Google Scholar
Book
1. Polis GA
2. Winemiller KO
(editors) (2018)
Food Webs: Integration of Patterns & Dynamics, Softcover Reprint of the Harcover (1st edition 1996. ed)

Dordrecht: Springer-Science+Business.
- Google Scholar
1. Pollierer MM
2. Langel R
3. Körner C
4. Maraun M
5. Scheu S
(2007) The underestimated importance of belowground carbon input for forest soil animal food webs
Ecology Letters 10:729–736.

https://doi.org/10.1111/j.1461-0248.2007.01064.x
- PubMed
- Google Scholar
(2009) Compartmentalization of the soil animal food web as indicated by dual analysis of stable isotope ratios (15N/14N and 13C/12C
Soil Biology and Biochemistry 41:1221–1226.

https://doi.org/10.1016/j.soilbio.2009.03.002
- Google Scholar
1. Post DM
(2002a) The long and short of food-chain length
Trends in Ecology & Evolution 17:269–277.

https://doi.org/10.1016/S0169-5347(02)02455-2
- Google Scholar
1. Post DM
(2002b) USING STABLE ISOTOPES TO ESTIMATE TROPHIC POSITION: MODELS, METHODS, AND ASSUMPTIONS
Ecology 83:703–718.

https://doi.org/10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2
- Google Scholar
(2013) Large 13C/12C and small 15N/14N isotope fractionation in an experimental detrital foodweb (litter–fungi–collembolans
Ecological Research 28:1069–1079.

https://doi.org/10.1007/s11284-013-1088-z
- Google Scholar
(2018) Non-vascular plants as a food source for litter-dwelling Collembola: Field evidence
Pedobiologia 66:11–17.

https://doi.org/10.1016/j.pedobi.2017.12.005
- Google Scholar
(2019a) Linking size spectrum, energy flux and trophic multifunctionality in soil food webs of tropical land-use systems
The Journal of Animal Ecology 88:1845–1859.

https://doi.org/10.1111/1365-2656.13027
- PubMed
- Google Scholar
(2019b) Trophic consistency of supraspecific taxa in below-ground invertebrate communities: Comparison across lineages and taxonomic ranks
Functional Ecology 33:1172–1183.

https://doi.org/10.1111/1365-2435.13309
- Google Scholar
(2019c) Uncovering trophic positions and food resources of soil animals using bulk natural stable isotope composition
Biological Reviews of the Cambridge Philosophical Society 94:37–59.

https://doi.org/10.1111/brv.12434
- PubMed
- Google Scholar
1. Potapov AM
2. Dupérré N
3. Jochum M
4. Dreczko K
5. Klarner B
6. Barnes AD
7. Krashevska V
8. Rembold K
9. Kreft H
10. Brose U
11. Widyastuti R
12. Harms D
13. Scheu S
(2020) Functional losses in ground spider communities due to habitat structure degradation under tropical land-use change
Ecology 101:e02957.

https://doi.org/10.1002/ecy.2957
- PubMed
- Google Scholar
1. Potapov AM
2. Rozanova OL
3. Semenina EE
4. Leonov VD
5. Belyakova OI
6. Bogatyreva VY
7. Degtyarev MI
8. Esaulov AS
9. Korotkevich AY
10. Kudrin AA
11. Malysheva EA
12. Mazei YA
13. Tsurikov SM
14. Zuev AG
15. Tiunov AV
(2021a) Size compartmentalization of energy channeling in terrestrial belowground food webs
Ecology 102:e03421.

https://doi.org/10.1002/ecy.3421
- PubMed
- Google Scholar
(2021b) Oil palm and rubber expansion facilitates earthworm invasion in Indonesia
Biological Invasions 23:2783–2795.

https://doi.org/10.1007/s10530-021-02539-y
- Google Scholar
1. Potapov AM
2. Beaulieu F
3. Birkhofer K
4. Bluhm SL
5. Degtyarev MI
6. Devetter M
7. Goncharov AA
8. Gongalsky KB
9. Klarner B
10. Korobushkin DI
11. Liebke DF
12. Maraun M
13. Mc Donnell RJ
14. Pollierer MM
15. Schaefer I
16. Shrubovych J
17. Semenyuk II
18. Sendra A
19. Tuma J
20. Tůmová M
21. Vassilieva AB
22. Chen TW
23. Geisen S
24. Schmidt O
25. Tiunov AV
26. Scheu S
(2022) Feeding habits and multifunctional classification of soil-associated consumers from protists to vertebrates
Biological Reviews of the Cambridge Philosophical Society 1:brv.12832.

https://doi.org/10.1111/brv.12832
- PubMed
- Google Scholar
Software
1. R Development Core Team
(2020) R: A Language and Environment for Statistical Computing
R Foundation for Statistical Computing, Vienna, Austria.

http://www.r-project.org
(2017a) Plant diversity, forest dependency, and alien plant invasions in tropical agricultural landscapes
Biological Conservation 213:234–242.

https://doi.org/10.1016/j.biocon.2017.07.020
- Google Scholar
(2017b) Biodiversity, Macroecology and Biogeography, Faculty of Forest Sciences and Forest Ecology of the University of Goettingen
Common Wayside Plants of Jambi Province 1:3979.

https://doi.org/10.3249/WEBDOC-3979
- Google Scholar
1. Rooney N
2. McCann K
3. Gellner G
4. Moore JC
(2006) Structural asymmetry and the stability of diverse food webs
Nature 442:265–269.

https://doi.org/10.1038/nature04887
- PubMed
- Google Scholar
1. Rooney N
2. McCann KS
(2012) Integrating food web diversity, structure and stability
Trends in Ecology & Evolution 27:40–46.

https://doi.org/10.1016/j.tree.2011.09.001
- PubMed
- Google Scholar
1. Rosseel Y
(2012)
lavaan: An R Package for Structural Equation Modeling

Journal of Statistical Software 48:1–36.
- Google Scholar
Book
1. RStudio Team
(2020)
RStudio: Integrated Development Environment for R

Boston, MA: RStudio, PBC.
- Google Scholar
(2018) Trophic redundancy reduces vulnerability to extinction cascades
PNAS 115:2419–2424.

https://doi.org/10.1073/pnas.1716825115
- PubMed
- Google Scholar
(2006) Effects of litter manipulation on early-stage decomposition and meso-arthropod abundance in a tropical moist forest
Forest Ecology and Management 229:285–293.

https://doi.org/10.1016/j.foreco.2006.04.007
- Google Scholar
(2003)
Evaluating the Fit of Structural Equation Models: Tests of Significance and Descriptive Goodness-of-Fit Measures

MPR-Online 8:e53.
- Google Scholar
1. Scheu S
2. Falca M
(2000) The soil food web of two beech forests (Fagus sylvatica) of contrasting humus type: stable isotope analysis of a macro- and a mesofauna-dominated community
Oecologia 123:285–296.

https://doi.org/10.1007/s004420051015
- PubMed
- Google Scholar
1. Schmitz OJ
2. Leroux SJ
(2020) Food Webs and Ecosystems: Linking Species Interactions to the Carbon Cycle
Annual Review of Ecology, Evolution, and Systematics 51:271–295.

https://doi.org/10.1146/annurev-ecolsys-011720-104730
- Google Scholar
1. Schneider K
2. Migge S
3. Norton RA
4. Scheu S
5. Langel R
6. Reineking A
7. Maraun M
(2004) Trophic niche differentiation in soil microarthropods (Oribatida, Acari): evidence from stable isotope ratios (15N/14N
Soil Biology and Biochemistry 36:1769–1774.

https://doi.org/10.1016/j.soilbio.2004.04.033
- Google Scholar
1. Schulz G
2. Schneider D
3. Brinkmann N
4. Edy N
5. Daniel R
6. Polle A
7. Scheu S
8. Krashevska V
(2019) Changes in Trophic Groups of Protists With Conversion of Rainforest Into Rubber and Oil Palm Plantations
Frontiers in Microbiology 10:240.

https://doi.org/10.3389/fmicb.2019.00240
- PubMed
- Google Scholar
(2015) “Trophic whales” as biotic buffers: weak interactions stabilize ecosystems against nutrient enrichment
The Journal of Animal Ecology 84:680–691.

https://doi.org/10.1111/1365-2656.12324
- PubMed
- Google Scholar
(2020) Trophic Structure of Small Invertebrates Inhabiting Litter of a Monsoon Tropical Forest
Russian Journal of Ecology 51:492–496.

https://doi.org/10.1134/S1067413620050112
- Google Scholar
(2021) Diversity, ecology, distribution and biogeography of Diplura
Insect Conservation and Diversity 14:415–425.

https://doi.org/10.1111/icad.12480
- Google Scholar
1. Starling JH
(1944) Ecological Studies of the Pauropoda of the Duke Forest
Ecological Monographs 14:291–310.

https://doi.org/10.2307/1948445
- Google Scholar
1. Steffan SA
2. Chikaraishi Y
3. Currie CR
4. Horn H
5. Gaines-Day HR
6. Pauli JN
7. Zalapa JE
8. Ohkouchi N
(2015) Microbes are trophic analogs of animals
PNAS 112:15119–15124.

https://doi.org/10.1073/pnas.1508782112
- PubMed
- Google Scholar
(2019) Conversion of rainforest to oil palm and rubber plantations alters energy channels in soil food webs
Ecology and Evolution 9:9027–9039.

https://doi.org/10.1002/ece3.5449
- PubMed
- Google Scholar
(2021) Trophic niche differentiation and utilisation of food resources in Collembola is altered by rainforest conversion to plantation systems
PeerJ 9:e10971.

https://doi.org/10.7717/peerj.10971
- PubMed
- Google Scholar
(2015) Stable isotope composition (δ(13)C and δ(15)N values) of slime molds: placing bacterivorous soil protozoans in the food web context
Rapid Communications in Mass Spectrometry 29:1465–1472.

https://doi.org/10.1002/rcm.7238
- PubMed
- Google Scholar
1. Tsiafouli MA
2. Thébault E
3. Sgardelis SP
4. de Ruiter PC
5. van der Putten WH
6. Birkhofer K
7. Hemerik L
8. de Vries FT
9. Bardgett RD
10. Brady MV
11. Bjornlund L
12. Jørgensen HB
13. Christensen S
14. Hertefeldt TD
15. Hotes S
16. Gera Hol WH
17. Frouz J
18. Liiri M
19. Mortimer SR
20. Setälä H
21. Tzanopoulos J
22. Uteseny K
23. Pižl V
24. Stary J
25. Wolters V
26. Hedlund K
(2015) Intensive agriculture reduces soil biodiversity across Europe
Global Change Biology 21:973–985.

https://doi.org/10.1111/gcb.12752
- PubMed
- Google Scholar
(2015) Intra-body variation and ontogenetic changes in the isotopic composition (13C/12C and 15N/14N) of beetles (Coleoptera
Entomological Review 95:326–333.

https://doi.org/10.1134/S0013873815030057
- Google Scholar
(2008) New multidimensional functional diversity indices for a multifaceted framework in functional ecology
Ecology 89:2290–2301.

https://doi.org/10.1890/07-1206.1
- PubMed
- Google Scholar
Book
1. Wickham H
(2016) Ggplot2: Elegant Graphics for Data Analysis
Cham: Springer.

https://doi.org/10.1007/978-3-319-24277-4
- Google Scholar
1. Wilkinson CL
2. Chua KWJ
3. Fiala R
4. Liew JH
5. Kemp V
6. Hadi Fikri A
7. Ewers RM
8. Kratina P
9. Yeo DCJ
(2021) Forest conversion to oil palm compresses food chain length in tropical streams
Ecology 102:e03199.

https://doi.org/10.1002/ecy.3199
- PubMed
- Google Scholar
1. Yang G
2. Wagg C
3. Veresoglou SD
4. Hempel S
5. Rillig MC
(2018) How Soil Biota Drive Ecosystem Stability
Trends in Plant Science 23:1057–1067.

https://doi.org/10.1016/j.tplants.2018.09.007
- PubMed
- Google Scholar

Article and author information

Author details

Zheng Zhou

JF Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany

Contribution
Data curation, Formal analysis, Visualization, Writing - original draft, Writing – review and editing

For correspondence
zzhou@gwdg.de

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8078-6378
Valentyna Krashevska

JF Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany

Contribution
Data curation, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9765-5833
Rahayu Widyastuti

Department of Soil Sciences and Land Resources, Institut Pertanian Bogor, Bogor, Indonesia

Contribution
Investigation

Competing interests
No competing interests declared
Stefan Scheu
1. JF Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany
2. Centre of Biodiversity and Sustainable Land Use, Göttingen, Germany
Contribution
Conceptualization, Funding acquisition, Supervision, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4350-9520
Anton Potapov
1. JF Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany
2. Russian Academy of Sciences, A.N. Severtsov Institute of Ecology and Evolution, Moscow, Russian Federation
Contribution
Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4456-1710

Funding

Deutsche Forschungsgemeinschaft (192626868-SFB 990)

Valentyna Krashevska
Rahayu Widyastuti
Stefan Scheu
Anton Potapov

China Scholarship Council (202004910314)

Zheng Zhou

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This study was funded by the Deutsche Forschungsgemeinschaft (DFG), project number 192626868-SFB 990 in the framework of the collaborative German-Indonesian research project CRC990. ZZ are supported by China Scholarship Council (CSC) (202004910314). We thank Dr Katja Rembold and Prof Holger Kreft for providing vegetation parameters; we also thank Zhijing Xie and Haifeng Yin for discussion. Special gratitude goes to Svenja Meyer for the animal silhouettes. We acknowledge support by the Open Access Publication Funds of the University of Göttingen.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,783

views
307

downloads
19

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Zheng Zhou
Valentyna Krashevska
Rahayu Widyastuti
Stefan Scheu
Anton Potapov

(2022)

Tropical land use alters functional diversity of soil food webs and leads to monopolization of the detrital energy channel

eLife 11:e75428.

https://doi.org/10.7554/eLife.75428

Categories and tags

Research organism

Arthropods

Share this article

Cite this article

Mean litter-calibrated Δ13C and Δ15N values of soil animal taxa in rainforest (a), jungle rubber (b), rubber, (c) and oil palm plantations (d).

Figure 1—source data 1

Figure 1—source data 2

Average Δ13C and Δ15N values of taxonomic groups in rainforest (F), jungle rubber (J), rubber, (R) and oil palm plantations (O).

Figure 2—source data 1

One-dimensional metrics for Δ13C (upper panel) and Δ15N values (lower panel) of communities in rainforest (F, green), jungle rubber (J, blue), rubber (R, red), and oil palm plantations (O, yellow).

Figure 3—source data 1

Figure 3—source data 2

Multidimensional isotopic metrics of soil animal communities in rainforest (F, green), jungle rubber (J, blue), rubber (R, red), and oil palm plantations (O, yellow).

Figure 4—source data 1

Figure 4—source data 2

Environmental drivers of community and energetic soil food-web metrics.

Figure 5—source data 1

Structural equation model on the effects of environmental change on food-web metrics.

Figure 6—source data 1

Author details

Zheng Zhou

Contribution

For correspondence

Competing interests

Valentyna Krashevska

Contribution

Competing interests

Rahayu Widyastuti

Contribution

Competing interests

Stefan Scheu

Contribution

Competing interests

Anton Potapov

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading

Mean litter-calibrated Δ¹³C and Δ¹⁵N values of soil animal taxa in rainforest (a), jungle rubber (b), rubber, (c) and oil palm plantations (d).

Average Δ¹³C and Δ¹⁵N values of taxonomic groups in rainforest (F), jungle rubber (J), rubber, (R) and oil palm plantations (O).

One-dimensional metrics for Δ¹³C (upper panel) and Δ¹⁵N values (lower panel) of communities in rainforest (F, green), jungle rubber (J, blue), rubber (R, red), and oil palm plantations (O, yellow).