Abstract

Quantifying the effects of nutrient additions on soil microbial respiration () and its contribution to soil respiration () are of great importance for accurate assessment ecosystem carbon (C) flux. Nitrogen (N) addition either alone (coded as LN and HN) or in combination with phosphorus (P) (coded as LN + P and HN + P) were manipulated in a semiarid alpine meadow on the Tibetan Plateau since 2008. Either LN or HN did not affect , while LN + P enhanced during peak growing periods, but HN + P did not affect . Nutrient addition also significantly affected , and the correlations of with climatic factors varied with years. Soil water content (Sw) was the main factor controlling the variations of . During the years with large rainfall variations, was negatively correlated with Sw, while, in years with even rainfall, was positively correlated with Sw. Meanwhile, in N + P treatments the controlling effects of climatic factors on were more significant than those in CK. Our results indicate that the sensitivity of soil microbes to climatic factors is regulated by nutrient enrichment. The divergent effects of Sw on suggest that precipitation distribution patterns are key factors controlling soil microbial activities and ecosystem C fluxes in semiarid alpine meadow ecosystems.

1. Introduction

Soil nutrient, especially nitrogen (N) availability, is an important factor determining primary productivity in many terrestrial ecosystems [13]. N enrichment by either N deposition or anthropogenic fertilization could enhance plant growth and promote net ecosystem primary productivity [4, 5]. Obviously, N enrichment can increase ecosystem carbon (C) sequestration by photosynthetic fixation on atmospheric carbon. However, the effect of N enrichment on CO2-C emissions from ecosystems, especially from the belowground portion, is inconsistent, being found to be positive [69], negative [1012], or neutral [1315]. Extensive research has focused on the effects of nutrient enrichment on soil respiration (), while the effects of nutrient enrichment on soil respiration () and the contribution of to have received less attention. As a key component of , investigations on the effects of nutrient enrichment on soil microbial activity improve our understanding of potential effects of global change.

Principally, consists of autotrophic respiration (primarily from roots) and heterotrophic respiration (mainly from soil microbes) [16]. Each of these components generally accounts for approximately 50% of the soil CO2 efflux [1618]. However, this proportion varies dramatically in different ecosystems because root and soil microbial respiration differ in their sensitivity to changes of nutrient conditions [19, 20]. N enrichment probably either enhances root respiration by increasing fine root production [21, 22] or reduces root respiration by stimulating fine root turnover [23] while N amendment can either stimulate by the increasing soil available N or reduce due to the suppression on soil microbes after input of the high N amount [8, 24]. Therefore, accurately quantifying can help us identify the source of variation in [25, 26]. Meanwhile partitioning of is essential to further detect the sensitivity of and to the increasing nutrient availability. Knowledge about the sensitivity is critical for accurately assessing the responses of belowground carbon fluxes to nutrient enrichment in the face of increasing of atmospheric N deposition, agricultural fertilization, and soil nutrient availability due to soil warming.

In the short term is strongly and positively related to temperature [27, 28]. Exponential function has been widely used to describe the sensitivity of to soil temperature. However, it has been proved that this relationship is regulated by soil water availability, especially in long-term and large scale [29]. In semiarid and arid region, precipitation is the main factor regulating plant growth and soil microbial activity [30, 31] due to the great contribution of precipitation to soil water availability [32]. Nutrient additions to soil can affect microbial activity [24, 33], yet it is not clear if these impacts can interact with the regulation effects of soil water availability.

The Tibetan Plateau covers about 2.5 million km2 with an average altitude of more than 4,000 m a.s.l. and 35% of that area is occupied by alpine meadows [34]. As the hinterland of the Tibetan Plateau, most of the meadows are located in the semiarid area. Most N is in organic form, since low temperature restricts decomposition of soil organic matter [35]. Therefore, plant growth and soil microbial activity are limited by the low soil N availability [7]. Moreover, phosphorus (P) is another essential element for plant growth and previous studies indicated that addition of N could induce P deficiency in grasslands [36]. N and P fertilization has been widely used as an efficient grassland management technique to recover degraded alpine meadow and to increase primary productivity to meet the needs of livestock [37, 38]. The N deposition rate in this area is about 7 kg N hm−2 yr−1 in the 2000s [39, 40] but is projected to increase to 40 kg N hm−2 yr−1 by 2050 [41]. In addition, the Tibetan Plateau is experiencing climatic warming [42, 43] and is predicted to experience “much greater than average” increases in surface temperature in the future [44] which mean that soil available nutrients may increase as the soil warming. Several lines of evidence have shown that the soil nutrient condition is improving in this meadow [29, 45]. Here we conducted a continuation of N and P addition study from 2008 and measured and during the growing season from 2010 to 2012. Our objectives were to (1) examine the effects of N and the combination of N and P on and the contribution of to during growing seasons and (2) detect how nutrient enrichment regulates the correlations of and in years with different variations in precipitation patterns. We hypothesize that (1) nutrient enrichment may decrease the contribution of to as belowground biomass was promoted by N + P treatments [46] and (2) the regulation effects of climatic factors on may intensify as soil microbial biomass was high in N + P treatments [47].

2. Materials and Methods

2.1. Study Site

The study site is located in the midsouth portion of the Tibetan Plateau in the grassland station of Damxung County (91°05′E, 30°51′N, 4333 m a.s.l.). This site is characterized as semiarid continental climate influenced by monsoons from the Pacific Ocean. Mean annual temperature is 1.3°C with a minimum of −10.4°C in January and a maximum of 10.7°C in July. Difference between diurnal temperatures is 18.0°C. Annual precipitation is 477 mm, 85% of which is concentrated from June to August. The data given here were mean values based on climate data collected from 1962 to 2006 [48]. This experiment was conducted in an alpine meadow dominated by Kobresia pygmaea, Stipa capillacea and Carex montis-everestii, and about 15 species present in per square meter. The dominant species contribute over 40% of the aboveground biomass in this meadow. Vegetation cover ranges from 30% to 50% depending on yearly precipitation [48]. The soil is classified as Mat-Gryic Cambisol, corresponding to Gelic Cambisol, with a depth of about 0.3–0.5 m. Soil particle composition is 67.02% of sand, 18.24% of silt, and 14.74% of clay [49]. Detailed soil characteristics can be found in literature [50].

2.2. Experimental Design

An area of 40 m × 40 m alpine meadow with uniform vegetation cover was selected as the field fertilization experiment site. Twenty-five 5 m × 5 m split plots were laid out in a complete randomized block design with 5 replicates for each of the 5 treatments which included CK (i.e., no nutrient addition), 2 levels of N enrichment, and 2 levels of combinations of N and P enrichment. Plots were separated by 2-meter aisles as buffering zones. A study carried out in temperate grassland indicated that N addition rate more than 10.5 g N m−2 yr−1 did not affect plant production either in mature or in degraded grasslands [38]. Thus, two levels of N, 5 and 10 g N m−2 yr−1 (hereafter coded as LN for low N and HN for high N, resp.), were manipulated in the alpine meadow. Considering that alpine communities in this semiarid region are probably colimited by N and P availability [51, 52], we chose to combine constant 5 g P m−2 yr−1 with the LN and HN treatments (hereafter coded as LN + P and HN + P, resp.). Granular CO(NH2)2 and (NH4)2HPO4 fertilizers were directly applied before plant seedling establishment in each year since 2008 (June 15 in 2010, June 8 in 2011, and June 15 in 2012). Fertilizers were applied in the evening in fine weather to reduce leaching or volatilization. Plots were located in winter rangelands which meant no grazing from May to September and grazed in other months. For each treatment, four replicate plots were randomly chosen for measurements.

2.3. Field Sampling and Measurements

Root exclusion was used to estimate soil microbial respiration by comparing CO2 efflux rates from soil surfaces with and without living roots [16, 53]. A previous study in our meadow site showed that roots in the top 0–0.15 m of soil accounted for more than 95% of the total root biomass (0–0.5 m) [54]. Therefore, we created a root-free soil quadrat in each plot, with an area of 0.5 m × 0.5 m and depth of 0.15 m, by removing both aboveground and belowground plant materials in every 5-cm layer and backfilling the soil according to the original order [55, 56], in May before fertilization. Then one polyvinyl chloride (PVC) collar (20 cm in diameter and 15 cm in height) was placed on the plant-free quadrat in each plot. The PVC collars were inserted into the soil to a depth of 13 cm which prevented root growth in them but allowed the movement of water and soil microbes [57]. The quadrats with collars were kept free of seedlings and plant growth by frequent manual removal during the growing season. Thus, it was assumed that CO2 efflux measured within the collars was derived only from soil microbes and this respiration was coded as [55]. Similar quadrats without root exclusion were also created in each plot and collars 5 cm in height were inserted into soil to a depth of about 3 cm. These collars retained intact plant root, and CO2 efflux measured in these collars was treated as . To minimize disturbance, the deep collars were installed one month before the first measurement and the aboveground plant material within the shallow collars was clipped to ground level and litter was removed 24 hours prior to each measurement [57]. A previous study in our meadow site indicated that daily mean values of and were quite close to the values measured between 09:00 and 11:00 a.m. [58]. Therefore, and were directly measured by a portable soil CO2 flux system (LI-8100, LI-COR Biosciences, Lincoln, NE, USA) between 09:00 and 11:00 a.m. on each measuring date. The interval between measurements was 10–15 days during the growing seasons from late June to late September in 2010, 2011, and 2012. We tried to select measuring dates at least two days after a rainfall event to avoid any pulse effect of precipitation on and .

Soil temperature (°C) and volumetric soil moisture (m3 m−3) at 5 cm depth were collected from the nearby eddy covariance system (100 m away from our experimental site) automatically recording every 30 minutes. The representative range of this observation system is about 200 m [58].

2.4. Statistical Analysis

Repeated measures ANOVAs were used to analyze effects of nutrient fertilization and sampling year on and , with measuring date as the repeated variables. Since the effects of year and date and fertilization and the interactions between year and date on and were significant (), and the interactions between year and fertilization significantly affected (, Table 1), we used a one-way ANOVA to test the differences in and among fertilization treatments followed by Tukey’s test for multiple comparisons in each year. Regression analyses were also used to test the correlations of with Sw and both in CK and the combination of N and P treatments (LN + P, HN + P) in the three consecutive years, respectively. In these regression analyses, only N + P (combined treatments of LN + P and HN + P) and CK were concerned due to significant effects of nutrient addition on occurring only in N + P treatments. All the analyses were performed in SPSS 16.0 (SPSS for Windows, version 16.0, Chicago, USA).

3. Results

3.1. Effects of Nutrient Enrichment on Seasonal Variations of

showed significant seasonal variations (Figures 1(A), 1(B), and 1(C); Table 1, P < 0.001). LN or HN addition did not affect relative to CK in each year of the experiment with the exception of higher under LN treatment in September 2010 (Figure 1(a)). However, LN + P and HN + P enhanced relative to CK in September 2010 (Figure 1(a)) and in August 2011 (Figure 1(b)). Similarly, LN + P increased relative to CK in August and September of 2012 (Figure 1(c)).

3.2. Effect of Nutrient Enrichment on Contribution of to

The contribution of to () varied with measured dates and years (Table 1, P < 0.001). In 2010, decreased sharply in September, while in 2011 and 2012 was relatively steady throughout the growing seasons (Figure 2). Nutrient enrichment also significantly affected but varied with years (Figure 2; Table 1, P < 0.001). In 2010, LN significantly increased in September, and HN increased it in July and August in 2011 (Figures 2(a) and 2(b)) while in 2011 HN + P decreased it compared with CK (Figures 2(a) and 2(b)), and this pattern also occurred in August in 2012 (Figure 2(c), P = 0.001). The average during the entire growing season under N + P and CK was 45% and 46.4% in 2010, 53% and 57.2% in 2011, and 73.5% and 80.6% in 2012.

3.3. Patterns of Precipitation, Soil Water Content, and Soil Temperature

Total precipitation from June to September was 360.6, 397.1, and 299.6 mm during 2010, 2011, and 2012, respectively, and the distribution patterns differed among these three years (Figure 3). In 2010 high precipitation events mainly concentrated at the end of the growing season (from mid-August to early September) (Figure 3(A), (a)) whereas, in 2011 and 2012, high precipitation events synchronized with the peak plant growth (from July to early August, Figures 3(B) and 3(C)). Specifically, precipitation in July was 63.6, 183.9, and 159.9 mm in 2010, 2011, and 2012, respectively, which amounted to 17.5%, 46.3%, and 53.4% of total precipitation from June to September (Figures 3(a), 3(b), and 3(c)). Moreover, precipitation in August was 169.8, 61.5, and 54.2 mm in 2010, 2011, and 2012, respectively, amounting to 47.1%, 15.4%, and 18.1% of total precipitation from June to September (Figures 3(a), 3(b), and 3(c)). In 2010 Sw remained with very low values till August and then increased sharply as rainfall rising, but in 2011 and 2012 they were very high from July to August and then decreased in late growing season. All the variations in Sw corresponded to the patterns of precipitation.

3.4. Correlations of with and Sw

Correlations between and climatic factors showed that the variations of were mainly controlled by Sw, but these correlations differed with years (Figure 4). was negatively correlated with Sw in 2010 with larger variations of rainfall than the other two years (Figure 4(d), , P < 0.001 in CK, , P < 0.001 in N + P treatments) while in 2011 and 2012 was positively correlated with Sw, and these correlations were significant only in N + P treatments (Figures 4(e) and 4(f)). Soil temperature had little effects on , and the correlation was only significant in N + P treatments in 2012 (Figure 4(c)).

4. Discussion

Our results showed that N addition alone at a rate greater than 5 g m−2 yr−1 did not affect during the growing season from the third to the fifth experimental year (2010 to 2012) of nutrient enrichment. However, both LN + P and HN + P treatments increased in most measuring dates, especially during the peak growth periods. Contrary to our assumption, LN + P did not decrease the contribution of to compared with CK. While consistent with the hypothesis, in N + P treatments the relationship between and Sw was more significant than in CK, and this relationship depended on precipitation distribution patterns. Our results indicate that the responses of to climatic factors were regulated by exogenous nutrient enrichment in this semiarid alpine meadow.

4.1. Effects of Nutrient Enrichment on and the Contribution of to

Positive effects of N fertilization on and have been found in some terrestrial ecosystems [10, 59]. As N availability is often limited in most terrestrial ecosystems, exogenous N enrichment can stimulate soil microbial activity and thus enhance CO2 flux from soil. Another N fertilization experiment conducted in the same meadow in 2010 found N addition at rates of 1, 2, and 4 g N m−2 yr−1 significantly increased plant aboveground biomass, , and [60]. However, in our N fertilization experiment adding N at the rate 10 g m−2 yr−1 did not significantly affect and (Figures 1 and 5). Additionally, N addition at the rate of 5 g m−2 yr−1 in this meadow did not significantly affect plant aboveground biomass (Table 2). Therefore, we presume that 5 g N m−2 yr−1 could be the saturation threshold for this alpine meadow ecosystem. Previous studies also found that ceased to continue to increase after years of fertilization [61, 62]. The neutral or suppression of and under long-term and high N addition treatments could be due to labile C depletion [63, 64], reduction of microbial biomass [64, 65], inhibition of microbial activity [66], and/or reduction of the belowground allocation [59]. Nitrogen additions significantly accelerated decomposition of light soil carbon fractions, which caused the depletion of labile C [63]. High rates of N addition could also lead to toxicity and reduction to plants and soil microbes by soil acidification, which resulted from the depletion of base cations and the release of ammonium from the soil [67, 68]. In addition, increased N availability may reduce the investment of photosynthetic labile carbon to root systems and subsequently decreased the supply of labile carbon to soil microbes [24]. In our study, although N addition at 10 g m−2 yr−1 did not significantly decrease and , the soil microbial biomass carbon indeed was much lower in HN than in CK (data were not presented). We presume that HN would decrease if we continue HN addition in the coming years.

In our study LN + P enhanced at the end of 2010 and in the middle of 2011 and 2012. Previous nutrient amendment studies showed that N availability was the main limitation of primary production in dry alpine meadows, while N and P availability colimited production in a wet alpine meadow in Colorado, USA [51, 69]. Another N and P fertilization experiment conducted in New Zealand showed that N addition changed the composition and activity of soil microbes, but the addition of P alone did not [70]. Additionally, it was found that P addition typically increased microbial biomass in the short term but generally decreased biomass over longer terms [71]. P addition alone was not concerned in our experiment and we failed to know if the positive effect derived from N + P or P addition alone. Therefore, the critical next step is to carry out P addition treatment in order to examine whether P limits plant growth and soil microbe activity in this meadow and to help in our understanding of the effects of N + P.

Our study showed was 46.4%, 57.2%, and 80.6% in CK during the 2010, 2011, and 2012 growing seasons, respectively, which are right within the range of 10%–90% across different ecosystems [16] and within the range of 35%–90% in a tallgrass prairie ecosystem over many years [72]. Generally, exhibits interannual variations with changes in climate and relevant physiological and ecological processes year by year [7274]. Contrary to our assumption, LN + P treatments did not decrease the contribution of to , while, in some sampling months, the effects of HN + P partly conformed to the assumption. Our previous study demonstrated that LN + P significantly enhanced belowground biomass in September [46]. It is well known that belowground biomass is the direct source of root respiration and this increase could enhance the contribution of root respiration to . However, the supply of liable carbon to soil microbes is mainly from root rhizodeposits and exogenous nutrient enrichment could stimulate soil microbial activities [24, 33]. Therefore, the contribution of to lies in the effects of nutrient enrichment on belowground biomass and soil microbial activities. In this study, LN + P treatment did not decrease the contribution of to , indicating that nutrient enrichment has a greater impact on soil microbial activities than belowground biomass. However, in some sampling months HN + P decreased relative to CK and we infer that high N addition could inhibit soil microbial activity.

4.2. Effects of Soil Temperature and Moisture on

CO2 emission from soil is regulated by several abiotic and biotic factors [53], such as temperature, rainfall events, soil moisture, soil physiochemical properties, and plant and soil microbial activities. In our study, the significant correlation between and Sw during growing season suggests soil water availability is an important factor limiting ecosystem C fluxes in this semiarid alpine meadow. However, differential responses of to soil water availability were also observed in these three years. The relationship between and Sw was negative in 2010 but positive in 2011 and 2012. This divergent correlation among years was attributed to the different precipitation patterns. Temperature varies usually in synchrony with precipitation during the growing season in this region; that is, high temperature and high precipitation corresponded to occurrence at the peak plant growth period, such as in 2011 and 2012, while high precipitation events mainly concentrated at the end of the growing season in 2010 owing to the late monsoon from the Pacific Ocean. For instance, 17.5% of total precipitation in July 2010 is far less than 46.3% and 53.4% of total precipitation in July 2011 and 2012, respectively. The continual and sharp increasing precipitation led to increase of Sw but decrease of in August 2010, which led to low soil microbial activities. Correspondingly, maintained high levels from July to mid-August and decreased gradually from then on (Figure 2). Therefore, the asynchronous variations of and precipitation accounted for the negative correlation between and Sw in 2010. However, in 2011 and 2012 with even distribution of precipitation during growing seasons, was positively correlated with Sw. The divergent effects of Sw on the contribution of to in years with different precipitation patterns suggest that precipitation distribution patterns are the key factors controlling soil microbial activities and ecosystem C fluxes in semiarid alpine meadow ecosystems.

As N availability is limited in alpine meadow ecosystems [35, 51], exogenous nutrient enrichment can stimulate soil microbial activities, and this has been proved in many previous studies [46, 75]. Although we did not measure soil microbial activities in this experiment, soil microbial biomass carbon was higher in N + P treatments than that in CK, especially in LN + P treatments [47], indicating that soil microbial activities were stimulated by exogenous nutrient enrichment. In 2011 and 2012 with even distribution of precipitation during growing seasons, the variations of did not depend on soil water availability in CK, while, in the N + P treatments with stimulated soil microbial activities, the variations of depended on soil water availability. This result indicated that nutrient enrichment mediates the relationships between soil microbial respiration and climatic factors through the stimulation of soil microbial activities.

5. Conclusions

N addition at a rate greater than 5 g N m−2 yr−1 did not significantly affect and plant aboveground biomass; we presume that 5 g N m−2 yr−1 could be the saturation threshold for this alpine meadow ecosystem. During the years with large variations of rainfall, was negatively correlated with Sw, while, in years with even distribution of rainfall, was positively correlated with Sw. The divergent effects of Sw on the contribution of to in years with different precipitation patterns suggest that precipitation distribution patterns are the key factors controlling soil microbial activities and ecosystem C fluxes in semiarid alpine meadow ecosystems. In the future climate change scenarios, spatial and temporal changes in precipitation patterns may have great impacts on semiarid alpine meadow ecosystems. Meanwhile, our results also indicate that the increase of exogenous N deposition in the future climate change scenarios may mediate the controlling effects of climatic factors on .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Elaine Stebler from Department of Natural Resource Ecology and Management at Oklahoma State University for language improvement. This research was supported by National Basic Research Program (973 Program) of China (no. 2010CB833502), the Strategic Priority Research Program of the Chinese Academy of Sciences “Climate Change: Carbon Budget and Relevant Issues” (Grant no. XDA05060700), the Excellent Young Scientists grant from the Institute of Geographic Sciences and Natural Resources Research, CAS (2011RC101), the Innovation Project grant from the Institute of Geographic Sciences and Natural Resources Research, CAS (201003012), and the National Basic Research Program of China (Grant no. 2010CB951704).