Loss of biodiversity impacts ecosystem functions, such as carbon (C)
cycling. Soils are the largest terrestrial C reservoir, containing more C
globally than the biotic and atmospheric pools together. As such, soil C
cycling, and the processes controlling it, has the potential to affect
atmospheric CO
Reductions in biodiversity have been linked with anthropogenic global change drivers such as climatic change, land cover change, reduction and fragmentation of natural areas, and human dependence on synthetic fertilizers (Vitousek and Mooney, 1997; Sanderson et al., 2002; Stevens, 2004; Phoenix et al., 2006; Fischer and Lindenmayer, 2007; Clark and Tilman, 2008). Over the past few hundred years, human activities have driven the species extinction rate to ca. 1000 times the typical background extinction rate (MEA, 2006). This global decline in biodiversity impacts important ecosystem functions, such as net primary production (NPP) and biogeochemical cycles of carbon (C) and nutrients, threatening the services that ecosystems provide to the human population (Wardle et al., 2011; Cardinale et al., 2012).
The C cycle is a particularly important ecosystem service because the
dynamic balance between C stored in ecosystems and in the atmosphere plays a
key regulatory role in the global climate. Although vegetation stores a
significant amount of C, soils are the largest terrestrial C reservoir,
containing more C globally than the biotic and atmospheric pools combined
(Lal, 2004). As such, soil C dynamics, and the processes that influence
them, have the potential to impact atmospheric carbon dioxide (CO
Loss of biodiversity has the potential to influence climate change if it
alters SOC pools by reducing ecosystem C uptake or by increasing CO
Soil communities typically have high levels of both taxonomic and functional
diversity (e.g., De Deyn and Van der Putten, 2005). High taxonomic
diversity, small sizes of organisms, and large population sizes make
characterization of soil communities much less straightforward than that of
plant communities. As such, characterization of soil organisms is often
based on body size (e.g., Bradford et al., 2002), grouping organisms into
macrofauna (
The positive impact of plant species diversity on soil C cycling processes can be mirrored in the soil community, with reported positive relationships between soil biodiversity and C cycling processes (e.g., Setälä, 2002; Heemsbergen et al., 2004). However, evidence suggests that this positive relationship is not universal, as other studies have found no significant impacts of soil biodiversity on C cycling (e.g., Griffiths et al., 2000). Understanding the relationship between soil biodiversity and C cycling is thus not so much a question of “does diversity matter?” but rather “under which circumstances does soil diversity affect C cycling?”. One possibility for addressing this question is to consider the role of functional similarity among taxa, because relatively small or no responses of ecosystem processes to loss or gain of soil biodiversity would be expected in the case of functional redundancy among soil organisms (Bengtsson, 1998; Andrén and Balandreau, 1999; Setälä et al., 2005).
To date, studies have assessed soil community diversity impacts on soil C cycling by manipulating diversity within or across multiple organismal groups (specifics of grouping criteria differ among studies, but are often taxonomic, functional, or body size groups). For example, studies have manipulated the diversity within groups of microorganisms (e.g., bacteria (Bonkowski and Roy, 2005; Griffiths et al., 2000) and mycorrhizal fungi (van der Heijden et al., 1998; Maherali and Klironomos, 2007)), soil mesofauna (e.g., microarthropods; Liiri et al., 2002), and macrofauna (Heemsbergen et al., 2004; Zimmer et al., 2005). Other studies have manipulated the diversity across groups of soil organisms that differ in body size (i.e., microbes, and micro-, meso-, or macrofauna) or, alternatively, have manipulated diversity of trophic or functional groups (Hedlund and Ohrn, 2000; Ladygina et al., 2010). Although different taxa within soil microbial (Cox et al., 2001; Hanson et al., 2008; Orwin et al., 2006) or faunal (Bignell and Eggleton, 2000; Milcu et al., 2008; Heemsbergen et al., 2004) groups can have unique impacts on the C cycle, functional redundancy among taxa would be expected to be reduced when a more complex food web of organisms is manipulated (e.g., across different size classes or feeding guilds) (Setälä, 2002). Thus, studies assessing biodiversity impacts on ecosystem processes across multiple groups of soil organisms may yield very different answers than studies that probe for biodiversity impacts within single groups.
Although our knowledge of relationships between soil biodiversity and soil C cycling processes has increased with expanding research emphasis in this area, a comprehensive understanding to date is hampered by a lack of quantitative synthesis of existing studies. Nielsen et al. (2011) performed the most extensive synthesis on this topic to date, with a qualitative analysis. They found that diversity is often positively correlated with ecosystem function (e.g., soil respiration), although they cautioned that negative relationships between soil biodiversity and C cycling may be related to experimental limitations. In particular, Nielsen et al. (2011) found that strong relationships between soil biodiversity and C cycling were most common under unrealistically low levels of diversity. Further, their synthesis showed that the soil community composition, rather than species richness per se, had significant impacts on C cycling processes. This indicates high levels of functional redundancy among soil organisms and suggests that a loss of soil biodiversity may not necessarily impact the C cycle.
We aimed to quantitatively assess how soil biodiversity affects soil C cycling pools and processes using meta-analysis. We tested the general hypothesis that soil biodiversity positively impacts the soil C cycle, where reductions in diversity decrease soil C pools and process rates. Further, we tested the hypothesis that biodiversity manipulations across groups of organisms with different body sizes more strongly affect C cycling processes than manipulations within groups, due to a higher degree of functional redundancy within than across groups (Andrén and Balandreau, 1999; Setälä, 2002). In addition, we tested whether diversity of soil microbes versus soil fauna (including micro-, meso-, and macrofauna) impacts C cycling differently. Finally, since “biodiversity” is a metric that differs greatly in absolute numbers for soil organisms that differ in body size, we evaluated how the relative loss of diversity (in percent) within body size groups (i.e., microbes, soil fauna) affects soil C cycling.
We compiled published studies that explicitly manipulated soil biodiversity and measured responses of soil C cycling pools and/or processes. We systematically searched ISI Web of Science using all possible combinations of one soil C search term (plant biomass, soil C, decomposition, respiration, or NPP), one soil organism search term (microbes, bacteria, fungi, microorganism*, soil fauna, soil biota, soil organism*, microfauna, mesofauna, macrofauna, nematode*, collembola, acari, termite*, earthworm*), and the term “*diversity”. We used “*” as a wildcard character such that papers using either singular or plural terms were returned. Additional relevant studies referenced in those returned by the search engines were also included in the literature compilation. While biodiversity sensu stricto includes both species richness and abundance (Magurran, 2005), we follow the recent nomenclature used in plant and soil studies and assume that the number of species present in a community represents the diversity of the community.
Each study included in our analysis presented data on one or more commonly measured biogeochemical C pool and/or process. Biogeochemical pool measurements were plant biomass and soil C pools (either total soil C, dissolved organic C (DOC), or microbial biomass). Measured biogeochemical processes were soil C respiration and plant tissue decomposition. The duration of manipulative experiments included in this analysis ranged from 14 days to 3 years. More studies were conducted under controlled laboratory and greenhouse conditions rather than under field conditions (37 and 8 studies, respectively). When extracting data from these studies, we took values directly from published tables or the text whenever possible. If necessary, we estimated values from graphical data with image analysis software (ImageJ, National Institutes of Health, Bethesda, MD, USA).
Overview of studies used in the analyses.
In total we analyzed 45 published studies, of which 8 examined the effects
of soil biodiversity on total plant biomass, 10 examined effects on
aboveground plant biomass, 9 examined effects on root biomass, 13 examined
effects on C respiration, 25 examined effects on decomposition, and 3 used
laboratory microcosms to examine effects on soil C pools (Table 1). For soil
C respiration, we included data that were obtained from either laboratory or
in situ incubation studies in which the substrate was soil only or soil with organic
amendments other than plant tissue (e.g., glucose). Laboratory studies
typically estimated potential C mineralization rates, using temperature and
moisture conditions assumed to be optimal for microbial activity. These
measurements were made in closed microcosms with flux rates estimated from
two or more repeated measurements of headspace gas concentrations. In situ studies
used static or flow-through chambers to measure CO
Soil biodiversity impacts on C respiration and decomposition were assessed by manipulating biodiversity either within a single body size group (i.e., microbes (including bacteria and fungi), micro-, or meso-, or macrofauna) or across multiple body size groups (e.g., micro-, meso-, macrofauna; e.g., Bradford et al., 2002). We treated the within-body size and across-body size groupings as two separate categories for the analysis. For plant biomass, however, there were not enough studies to run meta-analyses for individual categories. We also categorized the studies by soil microorganisms or soil fauna (micro-, meso-, and macro fauna grouped together due to inadequate numbers of studies to split these up). Categorizing studies in this manner allowed us to assess whether species diversity within or across body size groups affected C cycling differently, while also enabling us to compare the relative impacts of diversity within the soil microbial community versus soil biodiversity within the soil faunal community.
Many of the studies we compiled assessed soil biodiversity impacts on C
cycling by quantifying responses to a diversity gradient (e.g.,
To test how soil microbial and/or soil faunal diversity affects ecosystem C
pools (plant biomass) and processes (C respiration and decomposition), and
to test whether biodiversity manipulations across multiple body size groups
affected C cycling differently from manipulations within groups, we analyzed
the data set with meta-analysis (Curtis and Wang, 1998; Hungate et al.,
2009), using the statistical software MetaWin 2.0 (Rosenberg et al., 2000).
We were unable to use meta-analysis for soil C pools because the number of
studies available (3) was inadequate for a meaningful analysis. The response
ratio (
Conventional meta-analyses weight each individual observation by the
reciprocal of the mixed model variance (Curtis and Wang, 1998). However,
such an analysis requires that the standard deviations of individual studies
be known. These data were not available for a large proportion of the
studies used in our analysis. Thus, we weighted individual values included
in the analysis by experimental replication (Hedges and Olkin, 1985; Adams
et al., 1997), assuming that better-replicated experiments resulted in data
with lower variance. We choose this metric because well-replicated studies
provide more reliable estimates of the response of individual variables
(e.g., Hungate et al., 1996, 2009). We used bootstrapping to
calculate confidence intervals on mean effect size estimates for the whole
data set and for individual categories (Adams et al., 1997). We considered
diversity effects significant if the 95 % confidence intervals did not
overlap with zero. In addition, we considered diversity effects for
individual categories different from each other if they varied significantly
at the
Further, we tested how a loss of belowground species diversity is linked to
changes in C pools and processes by performing linear regressions with
percent change in species diversity and the effect size (ln
The percent response of total plant biomass, shoot biomass, and root biomass to a change in soil community diversity (i.e., “high” vs. “low” diversity). Studies included in the analysis manipulated diversity of the soil microbial community, the soil micro-, meso, or macrofaunal community or a combination of these trophic groups. Data represent means with 95 % confidence intervals; numbers in brackets represent the total number of data points included in the analysis.
Results from our meta-analyses indicate very different responses to changes in soil biodiversity among C pools (plant biomass; soil C pools not included because of inadequate number of studies) and processes (soil C respiration and plant tissue decomposition). Plant biomass did not respond to changes in diversity, either when analyzed as total biomass or partitioned into aboveground and belowground biomass (Fig. 1). In contrast to the lack of impact on plant C pools, decreased soil biodiversity (including studies that manipulated diversity within and across body size groups) corresponded to a mean 27.5 % reduction in soil C respiration (Fig. 2) and a mean 18 % reduction in decomposition (Fig. 3).
The percent response of soil C respiration to a change in soil community diversity (i.e., “high” vs. “low” diversity). Studies included in the analysis manipulated diversity of the soil microbial community; the soil micro-, meso, or macrofaunal community; or a combination of these body size groups (“all studies”). Further studies are categorized by studies that manipulated species diversity within or across body size groups, and by studies that manipulated the soil microbial community (including fungi and bacteria) or the soil faunal community (including micro-, meso-, and macrofauna). Data represent means with 95 % confidence intervals; numbers in brackets represent the total number of data points included in the analysis.
The percent response of decomposition to a change in soil community diversity (i.e., “high” vs. “low” diversity). Studies included in the analysis manipulated diversity of the soil microbial community; the soil micro-, meso, or macrofaunal community; or a combination of these trophic groups (“all studies”). Further studies are categorized by studies that manipulated species diversity within or across body size groups, and by studies that manipulated the soil microbial community (including fungi and bacteria) or the soil faunal community (including micro-, meso-, and macrofauna). Data represent means with 95 % confidence intervals; numbers in brackets represent the total number of data points included in the analysis.
When soil C respiration responses were partitioned into studies that manipulated diversity within body size groups versus those that manipulated diversity across body size groups, we found a significant effect only for within-group manipulations (Fig. 2). Due to a lack of studies that manipulated solely soil fauna (there was just one study: Scheu et al., 2002), we were unable to compare how a change in soil faunal biodiversity versus soil microbial biodiversity impacts soil C respiration. However, when we omitted studies manipulating soil faunal biodiversity from our analysis and assessed impacts of reductions in soil microbial biodiversity alone on soil C respiration, we found that soil C respiration was reduced by 41 % (Fig. 2).
Plant tissue decomposition generally responded negatively to reductions in soil biodiversity (Fig. 3). Studies that manipulated diversity within soil organismal groups and those that manipulated diversity across organismal groups both led to similar reductions in decomposition (means of 23 % and 15 %, respectively; Fig. 3). Further, whereas reduced soil microbial diversity did not significantly reduce decomposition rates, reduced soil faunal diversity led to a 37 % reduction in mean decomposition (Fig. 3).
Regression analyses revealed a negative linear relationship between soil
biodiversity and ln
Regressions between a percent change in the soil
microbial, soil macrofaunal, or soil microbial and soil faunal communities
(i.e., multiple organismal groups) and the natural log of the response ratio
of soil C respiration (ln
Regressions between a percent change in the soil
microbial, mesofaunal, macrofaunal, or soil microbial and soil faunal
communities (i.e., multiple organismal groups) and the natural log of the
response ratio of litter decomposition (ln
Regression analysis also revealed a significant response in decomposition to
altered biodiversity when all studies were included (Fig. 5a), but not when
only the highest and lowest diversity treatments were included (Fig. 5b).
Reductions in biodiversity did not significantly affect decomposition in
studies that measured litter mass loss. Conversely, when decomposition was
measured via CO
Changes in biodiversity have been linked with changes in ecosystem functioning, but so far studies have largely focused on plant diversity (e.g., Isbell et al., 2011). Whether or not similar patterns exist for soil biodiversity remains largely unknown. Here, we provide the first quantitative synthesis of studies testing effects of soil biodiversity on C cycling. Using meta-analysis and regression analysis, we showed that loss of soil biodiversity can have negative consequences for the soil C cycle but that relationships between C cycling processes with soil biodiversity vary across groups of soil organisms and are process-dependent. Below we explore how our findings contribute to our knowledge of how soil biodiversity drives ecosystem functions. We also discuss experimental shortcomings and methodological challenges and suggest directions for future research.
Few studies have assessed the relationship between soil biodiversity and soil C pools. We found just three studies in our literature search, and these studies all used different indices of soil C pools: ergosterol, which is a measure of fungal biomass (Liebich et al., 2007); dissolved organic C (Cragg and Bardgett, 2001); and soil organic C concentration (Zimmer et al., 2005). All three studies were of short duration (range: 42 to 70 days) and were conducted in microcosms in which diversity of the microbial community (Liebich et al., 2007), microfauna (collembola; Cragg and Bardgett, 2001), or macrofauna (earthworms and woodlice; Zimmer et al., 2005) was manipulated. Due to the small number of studies, we were unable to conduct a quantitative analysis. However, none of the individual studies showed an impact of soil biodiversity on soil C pools. It is probably not surprising that very few studies attempted to relate soil community diversity to soil C pools, due to the difficulty of maintaining soil biodiversity manipulations in microcosms for long time periods. Bulk soil organic C pools are typically stable on the order of years to decades due the large pool sizes and the relatively slow rates of biological processes (Conen et al., 2003; Smith, 2004), so short-term effect of soil biodiversity on this pool would be expected to be low and undetectable. We expect that some metrics of the relatively labile fraction in the soil C pool (e.g., particulate organic matter, DOC) will be more temporally dynamic than the bulk pool (Six et al., 2002) and would be better target response variables for assessing biodiversity impacts. Even more likely to provide information on soil biodiversity impacts on soil C cycling are studies assessing diversity effects on short-term C fluxes.
Plant biomass, assessed as a whole or partitioned into root and shoot biomass, was not significantly affected by soil biodiversity. Our analysis included studies that manipulated diversity of mycorrhizal fungi (Baxter and Dighton, 2001), microarthropods (Liiri et al., 2002), meso- and macrofaunal decomposers (Eisenhauer and Schädler, 2011), or more complex foodwebs consisting of multiple groups of different body sizes (Sulvaka et al., 2001; Laakso and Setälä, 1999; Ladygina et al., 2010; Eisenhauer et al., 2010; Eisenhauer and Schädler, 2011). The lack of a response of plant biomass production to soil biodiversity results from the contradicting results generated across a number of studies, and indicates that soil biodiversity does not unequivocally promote plant production (reviewed in van der Heijden et al., 2008). With the exception of rhizosphere organisms such as mycorrhizal fungi, rhizobia, and root pathogens or herbivores, linkages between soil organisms and plant biomass are indirect, i.e., decomposer organisms break down organic compounds and make nutrients available for plant uptake (Wardle et al., 2004). This indirect link between plant growth and soil organisms may result in a rather weak relationship between soil biodiversity and plant production (Balvanera et al., 2006). In line with this, we did find a trend of decreased plant production with loss of soil biodiversity. However, the limited number of studies reduced our statistical power, restricting our ability to quantify soil biodiversity impacts on plant biomass production. Another complication in assessing biodiversity impacts on plant production is that to date studies have exclusively been conducted in laboratory and greenhouse settings. While laboratory manipulations can provide useful information about potential controls over ecological processes, these manipulations are by necessity short-term (range: 35 days–52 weeks for the studies we compiled) and may include only a subset of the complex food webs and biogeochemical processes that occur in natural field settings (Hol et al., 2010). Furthermore, diversity effects may become more apparent over time as functional redundancy declines (Reich et al., 2012).
Although soil biodiversity did not conclusively impact soil C pools or affect plant biomass production, biodiversity as a whole appears to play an important role in maintaining and enhancing plant biomass production and soil C pools. For example, plant diversity can promote plant biomass production and soil C storage (Tilman et al., 2006; Cardinale et al., 2012; Fornara and Tilman, 2008; Steinbeiss et al., 2008), and these benefits of aboveground biodiversity on ecosystem functions are often attributed to increases in plant nutrient uptake resulting from belowground spatial and temporal differentiation in resource use (e.g., McKane et al., 1990; Tilman et al., 1996; Casper and Jackson, 1997; Schenk et al., 1999; van Ruijven and Berendse, 2005; van der Heijden et al., 2003). As such, a single limiting resource (e.g., soil N) may be spatially partitioned among co-occurring plant species, which reduces inter-specific competition and thereby facilitates coexistence (McKane et al., 1990). It is reasonable to expect that similar mechanisms occur for soil organisms (Prosser and Nicol, 2012; Sechi et al., 2015), and effects of plant and soil organism diversity on ecosystem functions may not be independent of each other, as increased plant diversity may be accompanied by increased soil biodiversity (Scherber et al., 2010; Eisenhauer et al., 2011a, b). If this is the case, soil biodiversity could explain, at least in part, the observed positive relationship between plant diversity and ecosystem C pools and processes. It is worthy of mention here that soil fungal pathogens have been found an important driver of observed positive relationships between plant diversity and productivity (Maron et al., 2011).
Results from our meta-analysis show that loss of soil biodiversity
significantly reduces soil C respiration (
Although many studies have assessed the impact of soil microbial diversity on soil C respiration, only one study evaluated effects of soil faunal diversity (earthworms) on C respiration (Scheu et al., 2002). This study indicated that soil faunal richness impacts soil microbial community functioning, which may in turn alter soil C respiration. Soil fauna can strongly affect microbial processes and community composition by enhancing the availability of resources to the microbial community (Edwards, 2000; Heal et al., 1997; Petersen and Luxton, 1982). For example, invasive earthworms in North America have been associated with changes in soil respiration, although these effects may be mediated through changes in the abundance of other organisms, such as microbes (Szlavecz et al., 2011) or availability of soil and litter resources (Huang et al., 2010). Thus, a change in the diversity of soil fauna is expected to alter soil C respiration, but more studies that manipulate soil fauna are required to confirm this hypothesis.
Results from the meta-analysis and the regression analysis show that loss of
soil biodiversity significantly reduces plant tissue decomposition
(
We hypothesized that diversity across multiple organismal groups composed of different body sizes would impact soil C cycling processes to a greater extent than diversity within groups of organisms with a single body size, with the assumption that within-group diversity would be accompanied by a greater degree of functional redundancy (Laakso and Setälä, 1999; Wardle, 1999; Cragg and Bardgett, 2001). For example, a higher diversity of feeding guilds has been linked to more effective substrate use (Setälä and Huhta, 1991). Our analysis, however, suggests that, for plant tissue decomposition, diversity across multiple body size groups has similar impacts on soil C cycling to diversity within body size groups. For soil C respiration, the effect of reduced diversity within groups was even stronger than that of smaller diversity across groups. This may result from the approaches taken to assess community impacts on soil C cycling. First, there is a lack of consistency in approaches taken to study effects of soil biodiversity on C cycling, both for the within- and the across-body size group approach. As such, the level of functional diversity between “high-” and “low-” diversity treatments may have varied across studies, and it is unclear whether shifts in functional diversity were greater for across-group manipulations than for within-group manipulations. Except for Heemsbergen et al. (2004), no studies have explicitly tested for the functional dissimilarity among the species manipulated. Second, populations of soil organisms at lower trophic levels may show compensatory growth responses to loss of biomass predation by organisms of higher trophic levels (e.g., Ingham et al., 1985), resulting in no net effect of manipulation of trophic diversity on the processes regulated by lower-trophic-level soil organisms. Third, effects of functional or trophic groups of organisms may have opposing effects on the C cycling pools and process rates. For example, Ladygina et al. (2010) showed that arbuscular mycorrhizal fungi and decomposer (enchytraeids and collembola) canceled each other out in affecting plant community biomass.
While manipulating diversity of any group of organisms is fraught with
challenges, manipulation of soil organism diversity is particularly
challenging. A more comprehensive assessment of soil diversity impacts on C
cycling will require that some of these challenges are addressed. Nielsen et al. (2011) found that the response to a reduction in diversity was greater
if diversity levels were low (i.e.,
If we are to fully understand how anthropogenically induced changes in biodiversity affect the terrestrial C cycle, we must dig deeper and embrace the challenges associated with studying the belowground world. Understanding the complex relationships between soil biodiversity and C cycling processes is currently limited by the sheer number of methodological concerns associated with these studies, which can greatly overestimate or underestimate the impact of soil biodiversity on soil C cycling, challenging extrapolation to natural field conditions. Nonetheless, our data point towards a definite importance of soil community diversity on key C cycling processes. Our quantitative analysis revealed significant negative effects of loss of soil biodiversity on rates of soil respiration and litter decomposition. If this is the case, declines in soil biodiversity could significantly affect the rates and dynamics of C cycling. However, biodiversity effects were not always consistent across groups of organisms. Differential responses of groups of organisms could be related to their functional role in the respective processes. It is, however, important to emphasize that we are still lacking a full understanding of the underlying changes in soil community functioning (or the lack of) with shifts in soil biodiversity (Nielsen et al., 2011; Briones, 2014). Future studies should therefore attempt to further elucidate the relative importance of taxonomic diversity (species numbers) versus functional diversity. Effects of loss of soil biodiversity on ecosystem C cycling should depend on the degree of functional dissimilarity of the organisms involved (Heemsbergen et al., 2004). Hence, unraveling the level of variation in functional traits among soil organisms, both within and across feeding groups, should be a future research priority. Such studies would improve predictions on the global C cycling in the face of future environmental changes. Given the importance of the soil community in regulating the direction and magnitude of C fluxes between the atmosphere and terrestrial ecosystems, advancing our understanding of soil biodiversity impacts on biogeochemical cycles may enhance the efficacy of climate change mitigation efforts.
All authors contributed to conception of the idea, compilation of the data, and preparation of the manuscript. Analyses were performed by M. A. de Graaff and J. Adkins.
This work was supported by Agriculture and Food Research Initiative competitive grant no. 2012-67010-20069 from the USDA National Institute of Food and Agriculture to M. A. de Graaff, and H. L. Throop's participation was supported by US National Science Foundation grants DEB-0953864 and DEB-1235828. Edited by: K. Denef