Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Nitrogen Addition and Warming Independently Influence the Belowground Micro-Food Web in a Temperate Steppe

  • Qi Li ,

    liq@iae.ac.cn (QL); liangw@iae.ac.cn (WL)

    Affiliations State Key Laboratory of Forest and Soil Ecology, Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang, China, Netherlands Institute of Ecology, Terrestrial Ecology Department, Wageningen, The Netherlands

  • Huahua Bai,

    Affiliation State Key Laboratory of Forest and Soil Ecology, Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang, China

  • Wenju Liang ,

    liq@iae.ac.cn (QL); liangw@iae.ac.cn (WL)

    Affiliation State Key Laboratory of Forest and Soil Ecology, Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang, China

  • Jianyang Xia,

    Affiliation State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing, China

  • Shiqiang Wan,

    Affiliations State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing, China, State Key Laboratory of Cotton Biology, Henan Key Laboratory of Stress Plant Biology, School of Life Sciences, Henan University, Henan, China

  • Wim H. van der Putten

    Affiliations Netherlands Institute of Ecology, Terrestrial Ecology Department, Wageningen, The Netherlands, Laboratory of Nematology, Wageningen University and Research Centre, Wageningen, The Netherlands

Abstract

Climate warming and atmospheric nitrogen (N) deposition are known to influence ecosystem structure and functioning. However, our understanding of the interactive effect of these global changes on ecosystem functioning is relatively limited, especially when it concerns the responses of soils and soil organisms. We conducted a field experiment to study the interactive effects of warming and N addition on soil food web. The experiment was established in 2006 in a temperate steppe in northern China. After three to four years (2009–2010), we found that N addition positively affected microbial biomass and negatively influenced trophic group and ecological indices of soil nematodes. However, the warming effects were less obvious, only fungal PLFA showed a decreasing trend under warming. Interestingly, the influence of N addition did not depend on warming.

Structural equation modeling analysis suggested that the direct pathway between N addition and soil food web components were more important than the indirect connections through alterations in soil abiotic characters or plant growth. Nitrogen enrichment also affected the soil nematode community indirectly through changes in soil pH and PLFA. We conclude that experimental warming influenced soil food web components of the temperate steppe less than N addition, and there was little influence of warming on N addition effects under these experimental conditions.

Introduction

Climate warming has been predicted to increase the global surface temperature by 1.8–4.0°C at the end of this century [1]. The rise in temperature could have profound effects on terrestrial ecosystems, such as changes in competition between species [2], altering plant productivity [3], [4], and in turn, influencing the supply of carbohydrates to belowground subsystems through root growth [5]. In addition to climate change drivers, terrestrial ecosystems are also affected by other global change phenomena, such as nitrogen deposition. It is predicted that global deposition of reactive N to the environment will increase from 100 Tg N yr−1 (in 1995) to 200 Tg N yr −1 by 2050 [6]. The intensive alteration of global nitrogen (N) cycles due to anthropogenic activities could change plant species composition and community structure [7], [8], [9], with consequent impacts on the structure and functions of soil ecosystems. Although individual effect of warming and N enrichment on ecosystem functioning has received wide attention [10], [11], [12], their combined effects are still unknown [13], especially on the responses of belowground organisms [13], [14], [15], [16], [17].

The few observations on the interactive effects of warming and N addition are relatively inconsistent. For example, there were additive effects of N addition and winter warming on plant productivity and soil N availability in temperate old fields [18], [19], and N addition increased the temperature sensitivity of the slowly cycling soil C pool in tropical forest [20]. In the Harvard Forest Long Term Ecological Research Site (LTER), the N availability diminished the warming effect on soil respiration in autumn [21], [22]. In a subarctic heath ecosystem, warming was found to negate the N addition effect on plant and microbial biomass after fifteen years of climate change manipulations [23]. On the other hand, N deposition and climate warming influenced litter decomposition and the associated microbial communities independently in low-alpine heath [24]. Given the important roles of soil biota in terrestrial ecosystems [15], understanding the direction and magnitude of interactive effects of N enrichment and warming on soil food webs and their components is crucial for predicting the changes in ecosystem structure and functioning under global climate change.

N addition and warming can each alter the activity of microbial decomposers, influencing the quantity of C lost from soils via respiration, and the transport of C from the surface into soils as dissolved organic C [25], [26]. N deposition can directly change soil C-cycling rates by inhibiting the microbial production of ligninolytic enzymes and enhancing cellulolytic enzyme activity [27]. Moreover, N deposition can also influence soil microorganisms and decomposition processes indirectly through altering plant composition and productivity by alleviating nitrogen limitation of plant growth [28]. Altogether, warming and N deposition can alter the rates of heterotrophic microbial metabolism in soil, and consequently the flow of C and N through soil food webs. While both warming and N deposition can impact soil biogeochemical processes, most of the research to date only quantified these effects independently, and the majority of these studies have focused on aboveground subsystems [10], [29]. Until now, there is relatively little knowledge on how different global change drivers interactively influence soil food web composition and functioning [15], [17], [30]. Since soil biota may influence biogeochemical cycling and physical conditions in terrestrial ecosystems, their responses to global changes are important at the ecosystem scale [15], [31]. Recent reviews of global change effects on soil biota also suggested that there is an urgent need for long-term studies to investigate interactions between different agents acting in concert [32], [33].

As part of the Eurasian grassland biome, the temperate steppe in northern China has suffered from over-grazing since 1980s, resulting in severe land degradation and soil N deficiency [34]. In addition, temperature in these areas increased substantially during past half century [35], and warming could stimulate evaportanspiration and reduced soil water availability, consequently aggravating soil water deficits in the semiarid regions [36]. Thus, it is predicted that the grassland ecosystems in this region are sensitive to climate change [37], [38]. However, there has been little information on how belowground biota respond to climatic change in this region, particularly to the combined driving factors, i.e. climate warming and N deposition. Therefore, a field experiment has been set up to examine the effects of climate warming, atmospheric N deposition, and their interactions. Our objectives were to explore the respective and interactive effects of N addition and warming on the micro-food web in a temperate steppe of Northern China. We focused on microbial biomass and nematode trophic groups as important components of soil micro-food web, because the micro-food web governs nutrient cycling and mineralization processes, and strongly determines the responses of belowground subsystem to climate change [15]. Previous studies at the field site we conducted our research have demonstrated that warming negatively affected root production [5], soil respiration, and microbial biomass [39]. In addition, N addition was found to stimulate gross ecosystem productivity and ecosystem C exchange by increasing C assimilation in the temperate steppe [26]. Accordingly, where they occur separately, warming and N addition appear to have an opposite effect on soil C dynamics [21]. Based on these previous studies and other reports [20], [33], [40], we hypothesized that N enrichment will positively affect micro-food web components, while warming will negatively influence them. Further, we hypothesized that warming and N addition will be counteracting each other, so that their interaction effect will be neutral. Structural equation modeling (SEM) was used to test whether N addition or warming influenced the components of soil food web directly or indirectly through changes in soil abiotic characters and/or plant growth.

Materials and Methods

Study site

This study was conducted at a field site of the Duolun Restoration Ecology Experimentation and Demonstration station of the Institute of Botany, the Chinese Academy of Sciences. The station is located at a temperate steppe in Duolun County, Inner Mongolia (42°02′N, 116°17′E, 1324ma.s.l), China. The soil in the study site is Chestnut soil (Chinese classification), or Haplic Calcisols according to the FAO classification, with sand, silt, and clay being 62.7%, 20.3%, and 17.0%, respectively. Mean bulk density is 1.31 g·cm−3 and pH is 7.7. The local climate in our study area has increased 0.45 °C over the past half century (1953–2005) [35]. N deposition in 2005–2006 was estimated at about 20 kg ha−1 yr−1 in this region [34]. The dominant plant species are Stipa krylovii Roshev., Artemisia frigida Willd., Potentilla acaulis L., Cleistogenes squarrosa (Trin.) Keng., Allium bidentatum Fisch. ex Prokh., and Agropyron cristatum (L.) Gaertn. Traditional land uses in this region include livestock grazing and farming. From the late 1950s to the 1970s, the government put strict policies into effect to ban grazing. Since the economic reforms and open-door policy in 1978, lands were thrown open to private use. After 2000, new bans on grazing have been carried out by the local government, but cattle and sheep grazing are still conducted in some areas of this region.

Experimental design

The experiment consisted of a complete random block design with six treatments, of which we included four treatments in the present study. There were six replicates of each treatment. Thirty-six plots of 3 m×4 m were arranged in a 6×6 matrix. The distance between two adjacent plots was 3 m. One of the six plots in each row was randomly assigned to one of the six treatments. In our study, we examined effects of control (CK), diurnal (24 h) warming (W), N addition (N), and diurnal warming plus N addition (WN) treatments. Effects of day and night warming were also included in the experimental design, but these treatments were not evaluated in the present study.

All the warmed plots were heated continuously by MSR-2420 infrared radiators (Kalglo Electronics Inc., Bethlehem, PA, USA) suspended 2.25 m above the soil surface. Across the four growing seasons from 2006 to 2009, warming treatment significantly increased the soil temperature at the depth of 10 cm by 1.79 °C. More detailed information of warming treatments on soil microclimates can be found in [41]. In each control or N addition plot, one ‘dummy’ heater with the same shape and size as the infrared heater was installed as in the heated plots to simulate the shading effects of the infrared radiator. All the heaters of the warming treatments were set at a radiation output of approximately 1600 W. The warming treatment started on 23 April 2006. N additions were spread by hand before the first rain event in the rainy season, and were applied once a year with NH4NO3 in the form of pellets. Given that N effects on species composition and ecosystem production saturate at N addition rates of about 10.5 g N m−2 yr−1 in this region [37], the level of N addition in this study was 10 g N m−2 yr−1. The experiment was started in 2006 and soil samples were collected in September 2009 and 2010. From each plot, five soil cores of 2.5 cm diameter were collected from a depth of 0–15 cm below soil surface and stored at 4°C until further analysis. In late August of each year, we clipped a 1×1 m2 quadrat in each subplot. Aboveground and belowground biomasses were measured. Root biomass at 0–15 cm depth was measured by soil auger (8 cm in diameter). The soil was carefully removed from the root system and the roots were thoroughly rinsed. The dry masses of aboveground and belowground biomass were determined by over-drying at 70 °C to constant weight.

Soil analyses

Soil temperature at the depth of 10 cm was recorded automatically with a Datalogger (STM-01 Soil Temperature Measurement System; Henan Electronic Institute, Zhengzhou, China). Six soil temperature measurements in each plot were collected with 10-min intervals and averages of the six measurements were stored as the hourly averages. The average soil temperature of each month is presented in (Fig. 1). The total soil carbon (SOC) and nitrogen (TN) were determined by a TruSpec CN Elemental Analyzer (Leco Corporation, USA). Soil pH was determined with a glass electrode in 1∶2.5 soil:water solution (w/v). Soil moisture (SM) was measured by weight loss after drying at 105°C for 48 h.

thumbnail
Figure 1. Monthly mean soil temperature in 2009 and 2010 as affected by nitrogen addition and warming in a temperate steppe.

https://doi.org/10.1371/journal.pone.0060441.g001

PLFA analysis

The soil microbial community was characterized using phospholipid fatty acid (PLFA) analysis as described by Bossio & Scow (1998) with slight modifications [42]. PLFA was extracted from 8 g freeze-dried soil, and the extracted fatty acid methyl esters were identified with a standard qualitative mix ranging from C9 to C30 and a MIDI peak identification system (Microbial ID. Inc., Newark, DE).

The sum of the following PLFAs was used as the measure of bacterial biomass: i15∶0, a15∶0, 15∶0, i16∶0, 16∶1ω7c, i17∶0, a17∶0, 17∶0, cy17∶0, and cy19∶0. The fatty acid 18∶2ω6 and 18∶1ω9c were used as indicator of fungal biomass [43] and 16∶1ω5c for AM fungi (AMF) [44], [45]. Taken together, all of the PLFAs indicated above were considered to be representative of the total PLFAs of soil microbial community.

Nematode community analysis

Nematodes were extracted from 100 g of soil (fresh weight) by a modified cotton-wool filter method [46], [47]. Nematode populations were expressed as number of nematodes per 100 g dry soil, and at least 150 nematodes from each sample were identified to genus level using an inverted compound microscope. The nematodes were assigned to the following trophic groups characterized by feeding habits: (1) bacterivores (Ba); (2) fungivores (Fu); (3) omnivore-predators (Om-Ca) and (4) plant parasites (PP) following Yeates et al. (1993) [48].

The following nematode community indices were calculated:

  1. Trophic diversity TD = 1/Σpi2; where pi is the proportion of trophic group i;
  2. Generic richness SR = (S−1)/ln(N), where S is the number of taxa and N is the number of total nematodes [49];
  3. Maturity index MI = ∑v(i)×f(i), where v(i) is the c-p value of taxon i according to their r and K characteristics following Bongers (1990) [44], f(i) is the frequency of taxon i in a sample;
  4. Plant parasite index PPI = ∑v(i)×f(i), which was determined in a similar manner with maturity index for plant parasitic genera [50];
  5. Structure index SI = 100×(∑ksns/(∑ksns+∑kbnb));
  6. Enrichment index EI = 100×(∑kene/(∑kene+∑kbnb)),

where kb is the weight assigned to guilds Ba2 and Fu2 and nb is the abundance of nematodes in guilds Ba2 and Fu2, which indicate basal characteristics of the food web; ks the weight assigned to guilds Ba3–Ba5, Fu3–Fu5, Om4–Om5 and Ca2–Ca5, ns is the abundance of nematodes in these guilds, which represent the structure condition of the food web; ke the weight assigned to guilds Ba1 and Fu2, and ne is the abundance of nematodes in these guilds, which represent an enriched condition of the food web [51]. Bax, Fux, Cax, Omx, (where x = 1–5) represent the functional guilds of nematodes that are bacterivores, fungivores, predators and omnivores where the guilds have the characters indicated by x on the colonizer-persister (cp) scale (1–5) following Bongers & Bongers (1998) [52].

Statistical analysis

Nematode abundances and the PLFA-derived biomass were ln(x+1) transformed prior to statistical analysis. General linear model analysis of variance was used to test the main effects and interactions of date, N addition and warming on soil properties and biota. The least significant difference (LSD) test was used when the treatment effects were significant. All statistical analyses were performed by SPSS statistical software (SPSS Inc., Chicago, IL). Difference at P<0.05 level was considered to be statistically significant. For the PLFA and nematode data, no interactive effects with sampling dates were observed in the ANOVA analysis, so the data were combined and presented for the average of two years.

Structural equation modeling (SEM) was used to investigate how N addition and warming affected soil microbe and nematode trophic groups in each year. We started the SEM procedure with an initial SEM model, based on our predictions [53]. In this model we hypothesized that N addition and climate warming may alter soil abiotic properties and plant biomass, which in turn may affect microbial community and soil nematode trophic groups. Soil abiotic characteristics, microbial community and nematode community were treated as latent variables [54]. We used soil moisture,soil organic C, total nitrogen and soil pH as the indicators of soil abiotic characteristics; bacterial, fungal and AMF PLFA as the indicators of microbial community; and the abundance of nematode trophic groups as the indicators of nematode community. In the initial model, the relationships between plant biomass and micro-food web were not significant and contributed little for the fit of the model, so plant biomass was removed from the final model for simplification. In the SEM analysis, we compared the model-implied variance-covariance matrix against the observed variance-covariance matrix with maximum likelihood estimation. We used the χ2 and its associated P-value to judge the model fit to the data. In the first phase, we constructed models containing all soil abiotic variables proposed to have influences on soil micro-food web. Based on the results of goodness-of-fit tests, we excluded less predictive measures (soil moisture, soil organic carbon and total nitrogen) and retained soil pH as the most informative abiotic variables in the final model. By stepwise removal of non-significant paths from the model we selected the model that fit our data best. SEM analyses were performed using AMOS 5.0 (Amos Development, Spring House, Pennsylvania, USA).

Results

Soil charateristics and plant biomass

Although soil characteristics (SM and pH) and plant biomass varied between 2009 and 2010 (Table 1), no warming effects were observed on the tested parameters. Soil pH was decreased by 8% and 1% in 2009 and 2010, respectively, following N addition (P<0.01).

thumbnail
Table 1. Soil characteristics and above and belowground plant biomass as affected by nitrogen addition (N) and warming (W) in a temperate steppe (Means ± SD).

https://doi.org/10.1371/journal.pone.0060441.t001

Soil microbial community

No sampling date effect was observed in the total and subgroup of PLFA biomass (Table 2). N addition effects on the PLFA biomass were more obvious than the warming effects. The average biomasses of total PLFA and the bacterial, fungal and AMF PLFA were enhanced 47.4%, 48.1%, 74.5% and 34.3%, respectively, by N addition (Fig. 2 a, b, c, d; P<0.05). Only fungal PLFA responded to warming, which resulted in lower fungal biomass under warming (Fig. 2c; P<0.05).

thumbnail
Figure 2. The PLFA biomasses for total (a) and indicator subgroups (b, c, d) and the abundance of nematode trophic groups (e–h) as affected by nitrogen addition and warming in a temperate steppe soil (Means

± SE). Bars indicate standard errors. Horizontal lines indicate the N treatment effects derived from General linear model analysis of variance. * and ** indicate N treatment effects significant at P<0.05 and P<0.01, respectively. Capital letters indicate significant differences among different treatments derived from LSD.

https://doi.org/10.1371/journal.pone.0060441.g002

thumbnail
Table 2. Results (P-values) of ANOVA on the effects of Date, warming (W), nitrogen addition (N) and their interactions on the PLFA biomass, nematode trophic groups and ecological indices.

https://doi.org/10.1371/journal.pone.0060441.t002

Soil nematode community

Warming and N addition didn't influence the total nematode abundance (data not shown). No warming effects were observed on nematode trophic groups and ecological indices. Among nematode trophic groups, only bacterivores varied between different sampling years (P<0.05), but no treatment effects were observed on the abundance of bacterivores and plant parasites (Fig. 2 e, f). Surprisingly, although the fungal PLFA was increased by N addition, the abundances of fungivores and omnivore-predators showed opposite patterns, as N addition reduced their abundances with 33.7% and 40.3%, respectively (Fig. 2 g, h; P<0.05). Among nematode ecological indices, trophic diversity, generic richness, maturity and structural indices were reduced by N addition (P<0.05), whereas plant-parasite index and enrichment index were not affected by the experimental treatments (Fig. 3).

thumbnail
Figure 3. Nematode ecological indices as affected by nitrogen addition and warming in a temperate steppe soil (Means

± SE). Bars indicate standard errors. * and ** indicate N treatment effects significant at P<0.05 and P<0.01, respectively. Capital letters indicate significant differences among different treatments derived from LSD.

https://doi.org/10.1371/journal.pone.0060441.g003

SEM results

Since no significant path way was observed in the final model of 2010, we only presented the SEM results of 2009. SEM analysis of 2009 data suggested that N addition altered soil nematode community composition indirectly through changes in soil microbial biomass (PLFA) and soil pH. However, warming effects were less obvious than those of N addition. The final models adequately fit the data on soil food web components in our study (standardized path coefficients are given in Fig. 4). The model of 2009 explained 35% of the variance in soil microbial biomass and 51% of variance in nematode abundance. The final SEM models indicated that soil food web components were stronger influenced by N addition directly rather than indirectly through changes in soil abiotic parameters (soil pH).

thumbnail
Figure 4. Structural equation models of N addition and warming effects on soil micro-food web components in a temperate steppe in 2009 (χ2 = 27.511; df = 24, P = 0.281; CFI = 0.975; RMSEA = 0.080).

Numbers on arrows are standardized path coefficients. Width of the arrows indicates the strength of the causal influence (non-significant pathways are dashed). N represents N addition effects; W, warming effects; Soil Char, soil characteristics; Microbe, soil microbial PLFA biomass; and Nema, soil nematode trophic groups.

https://doi.org/10.1371/journal.pone.0060441.g004

Discussion

N addition effects on soil micro-food web

In the present study, N addition enhanced microbial biomass, but negatively affected soil nematode trophic groups and community composition. The latter effect did not support our hypothesis. Low-levels of N addition have been reported to stimulate microbial growth by ameliorating N and C limitation through improving soil N availability and stimulating plant growth and litter decomposition [55]. However, the increase in microbial biomass did not lead to an increase in the abundance of nematode trophic groups. The differential responses of microbial biomass and nematode community composition might be driven by the complex interactions between bottom-up (resource quantity and quality) and top-down forces (regulation by predation) in the soil food web [56]. In the present study, microbial biomass increased while the abundances of fungivores and omnivore-predators decreased with N addition. One possibility was that microbial biomass in the temperate steppe was controlled by bottom-up forces, such as the quality of organic material and plant exudates [57], which increased under N additions. Whereas the fungivores and bacterivores were not controlled by bottom-up forces, because fungivores declined even though their resources (total PLFA) increased. Interactions between soil organisms are inherently complex and these bottom-up and/or top-down regulations in the soil food web may also depend on climate and physical conditions.

On the other hand, Bai et al. (2009) found that there is a threshold in the effects of nitrogen addition on aboveground biomass and plant functional composition in the Inner Mongolia grasslands [37]. Based on their findings, we inferred that there may be also a threshold for the N addition effects on the soil micro-food webs. Since high N additions might have a direct toxic effect on some saprophytic fungi by inhibiting their enzymes [58], which might also affect soil nematodes indirectly. SEM analysis also revealed that N addition could affect nematode community through changes in soil pH. Previous reports have suggested that higher N levels are detrimental to omnivores–predators [59], thus precluding any effects of temperature change. Soil acidification or ammonium suppression following N addition has been reported as important factors inhibiting soil nematodes [60].

Corresponding to the decrease in the abundances of fungivores and omnivore-predators, nematode ecological indices such as trophic diversity and generic richness were decreased by N additions, indicating a soil food web with lower diversity in both generic level and trophic community composition. The relatively lower structural and maturity indices under N additions also revealed a degraded food web relative to the treatments without N additions. Due to the trophic complexity associated with soil food webs, our knowledge is still limited about the processes controlling resource acquisition and use in belowground subsystems.

SEM results in this study revealed that the direct effects of N addition on soil food web components were more important than the indirect effects through alterations in soil abiotic characteristics or plant growth. Our observations were not in line with previous findings that changes in nematode community composition operated largely through effects on aboveground vegetation [16], [61]. Since we only tested above- and belowground biomass in the present study, and no significant treatment effects were found on plant biomass, the difference between our study and results by Kardol et al. (2010) and Veen et al. (2010) may be due to the species-specific relationships [16], [61]. De Deyn et al. (2007) found that nematode community composition responded differently to specific plants [62], and previous studies also suggest the important effects of plant species identity on multiple trophic levels in the soil food web [63]. In addition, it is also possible that soil biota may have been affected directly by living plant roots or indirectly through root exudates [64], which were not reflected in the plant biomass during the experimental period (3–4 years). In a long-term biodiversity experiment, Scherber et al. (2010) also found that plant species richness effects generally mediated the belowground responses through changes in root production or root exudates, but not through the vegetation biomass or the amount of litter input [65]. Our results underline suggestions made in other studies that responses of soil biological processes to N addition mainly resulted from direct effects other than indirect effects through changes in soil abiotic characters or plant growth [66], [67].

Warming effects on soil micro-food web

In our experiment, fungal PLFA biomass decreased following warming. In the same experimental field, the increased temperature associated with climate warming was found to stimulate evapotranspiration, leading to decreases in soil water content and supply of carbohydrates to belowground ecosystems [2], [5]. These changes may help to explain the decrease of fungal PLFA biomass. Similarly, in four heathland ecosystems along a climatic gradient in Europe, experimental warming was found to increase soil microbial biomass at the coldest and wettest site, and decrease soil microbial biomass at the warmest and driest site [40]. While the root-associated fungal community at tundra sites were quite resilient to warming effects [68]. The meta-analysis of 75 manipulative experiments on the responses of soil biota abundance to global change also revealed that effects of warming did not depend on taxon or body size of soil biota, and negative effects of warming were more likely to occur at the colder and drier sites [33].

Warming did not influence other micro-food web components tested in the present study, such as bacterial PLFA biomass or nematode trophic groups. This may be because bacteria and organisms feeding on them are less sensitive to climatic change than fungi, since bacteria mainly live inside soil aggregates and experience less extreme fluctuations in microclimate [69], [70]. In addition, other studies also found that soil moisture mediated the warming effects, thus null and negative effects of experimental warming have been attributed to the reductions in soil moisture [71], [72].

Interactive effects of warming and N addition on soil micro-food web

One of our objectives in this experiment was to investigate potential interactions between N deposition and climate change. In our study, no interactive effects of warming and N addition were observed on the micro-food web components, which was in support to our hypothesis. Only fungal biomass decreased with warming, so that the absence of interactive effects might be due to the minimal effects of warming on the micro-food web components. Similarly, in a low-alpine heathland, Papanikolaou et al. (2010) also found that there was lack of interactions between warming and N addition [24]. They concluded that N addition and warming acted independently on decomposing and the associated microbial community. Based on our observations in soil micro-food web, the effects of N addition did not depend on warming. Therefore, we conclude that in the temperate steppe ecosystem of our study climate warming and N addition influence soil micro-food web independently.

In summary, the results from our experiment indicated that N addition has led to significant changes in the soil micro-food web of temperate steppe, resulting in a lower diversity belowground micro-food web. However, effects of warming were less obvious than N addition and there were no interactive effects of N addition and warming on soil food web components. Direct effects of N addition on soil food web components explained the observed responses more than indirect effects through changes of soil abiotic conditions.

Acknowledgments

We are grateful to G.F. (Ciska) Veen and Shixiu Zhang for their help with structural equation modeling. We thank the staff of Duolun Restoration Ecology Experimentation and Demonstration Station. Comments and suggestions from Dr Amy Treonis and two anonymous reviewers greatly improved the quality of this manuscript.

Author Contributions

Conceived and designed the experiments: SW WL. Performed the experiments: QL HB JX. Analyzed the data: QL WvdP. Contributed reagents/materials/analysis tools: HB JX QL. Wrote the paper: QL WvdP WL.

References

  1. 1. Intergovernmental Panel on Climate Change (IPCC) (2007) Climate change 2007: the physical science basis. In: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (eds Solomon S, Qin D, Manning M, et al.), pp. 235–335; 747–846. Cambridge University Press, Cambridge, UK and New York .
  2. 2. Niu SL, Li ZX, Xia JY, Han Y, Wu MY, et al. (2008) Climatic warming changes plant photosynthesis and its temperature dependence in a temperate steppe of northern China. Environ Exp Bot 63: 91–101.
  3. 3. Weltzin J, Bridgham SD, Pastor J, Chen JQ, Harth C (2003) Potential effects of warming and drying on peatland plant community composition. Global Change Biol 9: 141–151.
  4. 4. Schuur EAG, Crummer KG, Vogel JG, Mack MC (2007) Plant species composition and productivity following permafrost thaw and thermokarst in Alaskan tundra. Ecosystems 10: 280–292.
  5. 5. Bai WM, Wan SQ, Niu SL, Liu WX, Chen QS, et al. (2009) Increased temperature and precipitation interact to affect root production, mortality and turnover in a temperate steppe: implications for ecosystem C cyling. Global Change Biol 16: 1306–1316.
  6. 6. Galloway JN, Townsend AR, Erisman JW, Bekunda M, Cai ZC, et al. (2008) Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science 320: 889–892.
  7. 7. Gruber N, Galloway JN (2008) An earth-system perspective of the global nitrogen cycle. Nature 451: 293–296.
  8. 8. Stevens CJ, Dise NB, Mountford JO, Gowing DJ (2004) Impact of nitrogen deposition on the species richness of grasslands. Science 303: 1876–1879.
  9. 9. Suding KN, Collins SL, Gough L, Clark C, Cleland EE, et al. (2005) Functional- and abundance-based mechanisms explain diversity loss due to N fertilization. Proc Natl Acad Sci USA 102: 4387–4392.
  10. 10. Xia JY, Wan SQ (2008) Global response patterns of terrestrial plant species to nitrogen addition. New Phytol 179: 428–439.
  11. 11. Bleeker A, Hicks WK, Dentener F, Galloway J, Erisman JW (2011) N deposition as a threat to the word's protected areas under the convention on biological diversity. Environ Pollut 159: 2280–2288.
  12. 12. Zhang NL, Xia JY, Yu XJ, Ma KP, Wan SQ (2011) Soil microbial community changes and their linkages with ecosystem carbon exchange under asymmetrically diurnal warming. Soil Biol Biochem 43: 2053–2059.
  13. 13. Clark CM, Tilman D (2008) Loss of plant species after chronic low-level nitrogen deposition to prairie grasslands. Nature 451: 712–715.
  14. 14. Antoninka A, Wolf JE, Bowker M, Classen AT, Johnson NC (2009) Linking above-and belowground responses to global change at community and ecosystem scales. Global Change Biol 15: 914–929.
  15. 15. Bardgett RD, Wardle DA (2010) Aboveground–Belowground Linkages, Biotic Interactions, Ecosystem Processes, and Global Change. Oxford Series in Ecology and Evolution, Oxford University Press, New York.
  16. 16. Kardol P, Cregger MA, Campany CE, Classen AT (2010) Soil ecosystem functioning under climate change: plant species and community effects. Ecology 91: 767–781.
  17. 17. Eisenhauer N, Cesarz S, Koller R, Worm K, Reich PB (2012) Global change belowground: impacts of elevated CO2, nitrogen, and summer drought on soil food webs and biodiversity. Global Change Biol 18: 435–447.
  18. 18. Hutchison JS, Henry HAL (2010) Additive effects of warming and increased nitrogen depostion in a temperate old field: plant productivity and the importance of winter. Ecosystems 13: 661–672.
  19. 19. Turner MM, Henry HAL (2010) Net nitrogen mineralization and leaching in response to warming and nitrogen deposition in a temperate old field: the importance of winter temperature. Oecologia 162: 227–236.
  20. 20. Cusack DF, Torn MS, McDowell WH, Silver WL (2010) The responses of heterotrophic activity and carbon cycling to nitrogen additions and warming in two tropical soils. Global Change Biol 16: 2555–2572.
  21. 21. Contosta AR, Frey SD, Cooper AB (2011) Seasonal dynamics of soil respiration and N mineralization in chronically warmed and fertilized soils. Ecosphere 2: 1–21.
  22. 22. Contosta AR, Frey SD, Ollinger SV, Cooper AB (2012) Soil respiration does not acclimatize to warmer temperatures when modeled over seasonal timescales. Biogeochemistry doi:10.1007/s10533-012-9749-6
  23. 23. Rinnan R, Michelsen A, Bååth E, Jonasson S (2007) Fifteen years of climate change manipulations alter soil microbial communities in a subarctic heath ecosystem. Global Change Biol 13: 28–39.
  24. 24. Papanikolaou N, Britton AJ, Helliwell RC, Johnson D (2010) Nitrogen deposition, vegetation burning and climate warming act independently on microbial community structure and enzyme activity associated with decomposing litter in low-alpine heath. Global Change Biol 16: 3120–3132.
  25. 25. Fenner N, Freeman C, Lock MA, Harmens H, Reynolds B, et al. (2007) Interactions between elevated CO2 and warming could amplify DOC exports from peatland catchments. Environ Sci Technol 41: 3146–3152.
  26. 26. Xia JY, Niu SL, Wan SQ (2009) Response of ecosystem carbon exchange to warming and nitrogen addition during two hydrologically contrasting growing seasons in a temperate steppe. Global Change Biol 15: 1544–1556.
  27. 27. Carreiro MM, Sinsabaugh Rl, Repert DA, Parkhurst DF (2000) Microbial enzyme shifts explain litter decay responses to simulated nitrogen deposition. Ecology 81: 2359–2365.
  28. 28. Bardgett RD, Freeman C, Ostle NJ (2008) Microbial contributions to climate change through carbon cycle feedbacks. ISME J 2: 805–814.
  29. 29. Wu ZT, Koch GW, Dijkstra P, Bowker MA, Hungate BA (2011) Responses of ecosystem carbon cycling to climate change treatments along an elevation gradient. Ecosystems 14: 1066–1080.
  30. 30. Fierer N, Strickland MS, Liptzin D, Bradford MA, Cleveland CC (2009) Global patterns in belowground communities. Ecol Lett 12: 1238–1249.
  31. 31. Bradford MA, Tordoff GM, Eggers T, Jones TH, Newington JE (2002) Microbiota, fauna, and mesh size interactions in litter decomposition. Oikos 99: 317–323.
  32. 32. Maraldo K, Holmstrup M (2010) Enchytraeids in a changing climate: a mini-review. Pedobiologia 53: 161–167.
  33. 33. Blankinship JC, Niklaus PA, Hungate BA (2011) A meta-analysis of responses of soil biota to global change. Oecologia 165: 553–565.
  34. 34. Zhang Y, Zheng LX, Liu XJ, Jickells T, Cape JN, et al. (2008) Evidence for organic N deposition and its anthropogenic sources in China. Atmos Environ 42: 1035–1041.
  35. 35. Wan SQ, Xia JY, Liu WX, Niu SL (2009) Photosynthetic overcompensation under nocturnal warming enhances grassland carbon sequestration. Ecology 99: 2700–2710.
  36. 36. Wan SQ, Luo Y, Wallace LL (2002) Changes in microclimate induced by experimental warming and clipping in tallgrass prairie. Global Change Biol 8: 754–768.
  37. 37. Bai Y, Wu J, Clark CM, Naeem S, Pan Q, et al. (2010) Tradeoffs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning: evidence from inner Mongolia Grasslands. Global Change Biol 16: 358–372.
  38. 38. Christensen L, Coughenour MB, Ellis JE, Chen ZZ (2004) Vulnerability of the Asian typical steppe to grazing and climate change. Climatic Change 63: 351–368.
  39. 39. Liu WX, Zhang Z, Wan SQ (2009) Predominant role of water in regulating soil and microbial respiration and their responses to climate change in a semiarid grassland. Global Change Biol 15: 184–195.
  40. 40. Sowerby A, Emmett B, Beier C, Tietema A, Penuelas J, et al. (2005) Microbial community changes in heathland soil communities along a geographical gradient: interaction with climate change manipulations. Soil Biol Biochem 37: 1805–1813.
  41. 41. Xia J, Chen S, Wan S (2010) Impacts of day versus night warming on soil microclimate: Results from a semiarid temperate steppe. Sci Total Environ 408: 2807–2816.
  42. 42. Bossio DA, Scow KM (1998) Impacts of carbon and flooding on soil microbial communities: phospholipids fatty acid profiles and substrate utilization patterns. Microb Ecol 35: 265–278.
  43. 43. Bååth E, Anderson TH (2003) Comparison of soil fungal/bacterial ratios in a pH gradient using physiological and PLFA-based techniques. Soil Biol Biochem 35: 955–963.
  44. 44. Bardgett RD, Hobbs PJ, Frostegård Å (1996) Changes in soil fungal:bacterial biomass ratios following reductions in the intensity of management of an upland grassland. Biol Fertil Soils 22: 261–264.
  45. 45. Frostegård Å, Bååth E (1996) The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil. Biol Fertil Soils 22: 59–65.
  46. 46. Oostenbrink M (1960) Estimating nematode populations by some selected methods. In: Sasser, J N, Jenkins, W R. (Eds) Nematology. University of North Carolina Press, Chapel Hill.
  47. 47. Townshend JL (1963) A modification and evaluation of the apparatus for the Oostenbrink direct cottonwool filter extraction method. Nematologica 9: 106–110.
  48. 48. Yeates GW, Bongers T, de Goede RGM, Freckman DW, Georgieva SS (1993) Feeding habits in soil nematode families and genera - an outline for soil ecologists. J Nematol 25: 315–331.
  49. 49. Yeates GW, Bongers T (1999) Nematode diversity in agroecosystems. Agric Ecosyst Environ 74: 113–135.
  50. 50. Bongers T (1990) The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 83: 14–19.
  51. 51. Ferris H, Bongers T, de Goede RGM (2001) A framework for soil food web diagnostics: extension of the nematode faunal analysis concept. Appl Soil Ecol 18: 13–29.
  52. 52. Bongers T, Bongers M (1998) Functional diversity of nematodes. Appl Soil Ecol 10: 239–251.
  53. 53. Grace JB (2006) Structural Equation Modeling and Natural Systems. Cambridge University Press, Cambridge.
  54. 54. Shipley B (2002) Cause and correlation in biology. A user's guide to path analysis, structural equa-tions and causal inference. Cambridge University PressCambridge.
  55. 55. Zhang NL, Wan SQ, Li LH, Bi J, Zhao MM, et al. (2008) Impacts of urea N addition on soil microbial community in a semi-arid temperate steppe in northern China. Plant Soil 311: 19–28.
  56. 56. Hoeksema JD, Lussenhop J, Teeri JA (2000) Soil nematodes indicate food web responses to elevated atmospheric CO2. Pedobiologia 44: 725–735.
  57. 57. Bardgett RD, Chan KF (1999) Experimental evidence that soil fauna enhance nutrient mineralization and plant uptake in montane grassland ecosystems. Soil Biol Biochem 31: 1007–1014.
  58. 58. Fog K (1988) The effect of added nitrogen on the rate of decomposition of organic matter.Biol Rev Camb Phil Soc. 63: 433–462.
  59. 59. Liang WJ, Lou YL, Li Q, Zhong S, Zhang XK, et al. (2009) Nematode faunal response to long-term application of nitrogen fertilizer and organic manure in Northeast China. Soil Biol Biochem 41: , 883–890.
  60. 60. Wei CZ, Zheng HF, Li Q, Lü XT, et al. (2012) Nitrogen addition regulates soil nematode community compostion through ammonium suppression. PLoS One 7: e43384.
  61. 61. Veen GF, Olff H, Duyts H, Van der Putten WH (2010) Vertebrate herbivores influence soil nematodes by modifying plant communities. Ecology 91: 828–835.
  62. 62. De Deyn GB, van Ruijven J, Raaijmakers CE, de Ruiter PC, van der Putten WH (2007) Above- and belowground insect herbivores differentially affect soil nematode communities in species-rich plant communities. Oikos 116: 923–930.
  63. 63. Wardle DA, Yeates GW, Williamson W, Bonner K (2003) The response of a three trophic level soil food web to the identity and diversity of plant species and functional groups. Oikos 102: 45–56.
  64. 64. Wardle DA, Bardgett RD, Klironomos JN, Setälä H, Van der Putten WH, et al. (2004) Ecological linkages between aboveground and belowground biota. Science 304: 1629–1633.
  65. 65. Scherber C, Eisenhauer N, Weisser WW, Schmid B, Voigt W, et al. (2010) Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. Nature 468: 553–556.
  66. 66. Manning P, Newington JE, Robson HR, Saunders M, Eggers T, et al. (2006) Decoupling the direct and indirect effects of nitrogen deposition on ecosystem function. Ecol Lett 9: 1015–1024.
  67. 67. Manning P, Morrison SA, Bonkowski M, Bardgett RD (2008) Nitrogen enrichment modifies plant community structure via changes to plant-soil feedback. Oecologia 157: 661–673.
  68. 68. Fujimura KE, Egger KN, Henry GHR (2008) The effect of experimental warming on the root-associated fungal community of Salix arctica. ISME J 2: 105–114.
  69. 69. Hattori T (1988) Soil aggregates in microhabitats of microorganisms. Rep Inst Agric Res, Tohoku Univ 37: 23–36.
  70. 70. Denef K, Six J, Bossuyt H, Frey SD, Elliott ET, et al. (2001) Influence of dry-wet cycles on the interrelationship between aggregate, particulate organic matter, and microbial community dynamics. Soil Biol Biochem 33: 1599–1611.
  71. 71. De Boeck HJ, Lemmens CMHM, Zavalloni C, Gielen B, Malchair S, et al. (2008) Biomass production in experimental grasslands of different species richness during three years of climate warming. Biogeosciences 5: 585–94.
  72. 72. Verburg PSJ, Johnson DW, Schorran DE, Wallace LL, Luo Y, et al. (2009) Impacts of an anomalously warm year on soil nitrogen availability in experimentally manipulated intact tallgrass prairie ecosystems. Global Change Biol 15: 888–900.