Variations in the natural 13C and 15N abundance of plants and soils under long-term N addition and precipitation reduction:interpretation of C and N dynamics

Guoyong Yan,Shijie Han,Mingxin Zhou,Wenjing Sun,Binbin Huang,Honglin Wang,Yajuan Xing1,*and Qinggui Wang*

Abstract

Keywords: δ13C,δ15N, N addition, Precipitation reduction, Nutrient acquisition strategies

Background

Anthropogenic activities, such as fossil fuel combustion and fertilizer application, have greatly increased the amount of reactive nitrogen (N) deposited onto terrestrial ecosystems (Galloway et al. 2008; Liu et al. 2013).Increases in soil N availability can significantly alter the carbon (C) and N cycles by influencing the biochemical processes of plants and soils (Chen and Ruan 2018;Schulte-Uebbing and de Vries 2018). Similarly, precipitation reduction, which is expected to become more frequent as climate changes (Ciais et al. 2013), can directly or indirectly affect C and N cycling by altering the physical and biochemical processes of soil (Borken and Matzner 2009; Díaz-Guerra et al. 2018). Although the effects of N deposition and changes of precipitation pattern on ecosystem C and N dynamics have been the hotspot of ecological research for decades (Waldrop et al. 2004;Xiang et al. 2008; Thomas et al. 2010; Díaz-Guerra et al.2018; Souza et al. 2019), how N deposition, precipitation reduction, and their interaction regulate the ecosystem C and N cycling is still unclear.

Based on previous studies, δ13C and δ15N provide a valuable approach exploring the C and N dynamics in different ecosystem components (Schlesinger 2009;Werth and Kuzyakov 2010). By analyzing δ15N values in soils and plants, the fates and retention of deposited N in forest ecosystems can be better understood (Schlesinger 2009).Increased N availability caused by N deposition may lead to an increase in net nitrification, NO3?leaching and gaseous N loss in soil ecosystem (Wallenstein et al. 2006; Gao et al. 2015). Atmospheric N deposition and fertilizer applied are generally isotopically depleted sources in comparison to δ15N values of soils and plants (Amundson et al. 2003). Processes associated with N loss under N deposition, e.g., mineralization, nitrification, denitrification, and N leaching discriminate against the heavier N isotope (15N) in favor of the lighter N isotope (14N) in soil ecosystems (Robinson 2001;Brookshire et al. 2012). A high loss of lighter14N causes overall enrichment of15N in soil ecosystems and indicates the openness of the N cycle (Watzka et al. 2006).Consequently, soil and plants have been expected to be more15N enriched because of the massive losses of the15N-depleted N compounds (Gurmesa et al. 2016). Furthermore, the plant absorbs N from soil mainly through fine roots and mycorrhizal fungi. The mycorrhizal fungi preferentially transfer14N and against15N to their host plants (Hobbie et al. 2000). Therefore, plant’s δ15N value under N deposition may reflect the forms and status of soil N utilized, mycorrhizal association type, and the subsequent fractionations in different plant organs (Conrad et al. 2018; Chen et al. 2019). Besides, the close coupling between plant leaf N concentration and C gain.The changes of N cycling under N deposition may alter the enrichment of plant-soil δ13C. Nevertheless, our knowledge of how the N deposition affects plant-soil δ13C and its change mechanisms is limited.

In addition to N deposition, change of precipitation pattern also affects δ13C and δ15N values in plants and soils. Previous studies have suggested that plant and soil δ15N values decreased with increasing precipitation intensity because the forms and the isotopic composition of the N left the ecosystems change (Schuur and Matson 2001). Studies have reported that plant-soil was much more enriched with δ15N in drought regions than that in wet regions, which, as a result of N cycling, was likely to become more open as the precipitation was reduced(Swap et al. 2004; Cheng et al. 2009). However, it remains unclear how precipitation reduction causes the loss of14N and enrichment of15N in plants and soils.For δ13C values, a strong correlation has been found between plant δ13C values and soil water relations, which thus can be a powerful tool to evaluate plants waterefficient (Klaus et al. 2016). With precipitation reduction, plants absorb heavier13C due to closing stomata leading to less negative δ13C values compared to plants with better water use efficiency (Farquhar et al. 1989;Mariotte et al. 2013). Therefore, the δ13C values have been successfully used as an indicator of drought stress of plants (Mariotte et al. 2013). The effects of precipitation changes on δ13C values have been mainly focused on sensitive grasslands in previous studies (Wrage et al.2009; Mariotte et al. 2013), and it is still unclear how forest δ13C values respond to precipitation reduction.

In general, plant-soil δ13C and δ15N values are related to the availability of soil N and water and are indicative of C and N cycling on different spatial scales (from soil to plant). Nitrogen rich soil system caused by N deposition tends to be15N enriched where the N cycle is more open to N loss. In such N rich soils, precipitation reduction may support low losses of N through the reducing leaching, indicating that precipitation patterns may potentially influence N deposition effects. Furthermore, increased soil N availability was found to potentially increase water use efficiency in plants, increasing their δ13C enrichments (Ripullone et al. 2004). Thus, under the simultaneous scenarios of N deposition and precipitation reduction, it is necessary to understand how their interactions alter the δ13C and δ15N values of plant-soil for accurately grasping ecosystem C and N dynamics.However, little evidence is collected in the forest to support the authenticity of this interaction.

We identified some of the remaining questions needed to be answered in order to advance our understanding of plant-soil δ13C and δ15N values responding to N deposition and precipitation reduction and consequently the ecosystem’s C and N cycle. Therefore, in this study,we examined the variation of δ13C and δ15N values in fresh plant foliage, litterfall, fine roots, coarse roots, fungal sporophores, and different soil layers in order to explore the effects of long-term stimulation of N deposition and precipitation reduction on the spatial variation of δ13C and δ15N values. We sought to determine whether the measurement of δ13C and δ15N values could be used as indicators of the relative sensitivity of soil C and N cycling across different tree species in response to N deposition and precipitation reduction. We did this by measuring the δ13C and δ15N values of three dominant species,including Fraxinus mandshurica,Tilia amurensis and Pinus koraiensis. We hypothesized that N deposition and precipitation reduction might lead to more open C and N cycling.

Methods

The study site was located at an old broad-leaved Korean pine mixed forest (over 300 years old) in the Changbai Mountains Natural Reserve, northeastern China (42°24′ N, 128°06′ E, 738 m above the sea level).The region is a typical temperate-continental climate,characterized by cold and long winters and warm and short summers. The mean annual precipitation is about 740 mm with more than 80% falling between May and October. The annual average temperature is approximately 3.6°C with an average growing season (from May to October) temperature of 15°C. The soil of the study site is classified as Eutric Cambisol (FAO classification)with 31.54% sand, 42.18% silt, 26.28% clay, and 25.42%organic matter in the topsoil in the 0-20 cm soil layer.Dominant tree species in this forest include Pinus koraiensis, Fraxinus mandshurica and Tilia amurensis. The tree density on the site was 432 trees·ha?1. The main shrub species include Euonymus alatus, Philadelphus schrenkii, Corylus mandshurica, Lonicera japonica, and the main herbaceous species include Anemone cathayensis, Anemone raddeana, Funaria officinalis, Cyperus microiria, and Filipendula palmate (Zhou et al. 2019).

Experimental design

In May 2009, six 50 m×50 m plots were randomly established in the forest, with at least 100 m buffer strips between each plot. Previous observation data showed that the precipitation was approximately 550 mm in drought years, which was almost 30% less than the longterm average annual precipitation of 740 mm in the last 30 years (Chinese Ecosystem Research Net data). To test the effect of precipitation reduction on ecosystem C and N cycling, three plots were treated with a 30% precipitation reduction and the other three were set as the control. The high light V-shaped panels, with high light transmittance, were installed to intercept through-fall.The V-shaped panels cover 30% of the plot area. The precipitation interception facility was approximately 1 m aboveground and was available to ensure normal air flow. The panels were removed in winter to avoid the disturbance of uneven snow reduction. To study the effect of N deposition on ecosystem C and N cycle, each plot was divided into two subplots (25 m×50 m): one subplot was used to simulate N deposition and the other subplot was used as a control. The N addition level was 50 kg N·ha?1·yr?1and approximately two times the real N deposition rate in the study area (Lü and Tian 2007).Ammonium nitrate (NH4NO3) was used as a nitrogen source in the study. The solution of NH4NO3was sprayed on the forest floor monthly using a backpack sprayer with six equal applications over the entire growing season (May to October). In every application, the NH4NO3was weighed and mixed with 40 L of water and each control plot received 40 L water without N fertilizer.

Plant, soil and fungal sporophore sampling

The fresh foliage, litterfall, fine roots (<2 mm), and coarse roots (>2 mm) of Pinus koraiensis (PKS), Tilia amurensis (TAR) and Fraxinus mandshurica (FMR)were collected. In July 2017, all the plant samples from the three separate plants per species were collected in each subplot at the seasonal peak of aboveground biomass. Meanwhile, soil samples were collected from beneath each tree in four directions using a stainless soil corer (5-cm inner diameter), then divided into three soil layers (0-10, 10-20 and 20-30 cm soil layer). The six soil cores from each soil layer were mixed to form a composite sample in each subplot. All of the soil samples were put through a 2-mm soil sieve while the roots,stones,and other debris were removed.Six fungal sporophores were collected in the forest floor in each subplot.In the laboratory, all samples for C concentration,13C,N concentration, and15N analysis were dried at 65°C to a constant weight and then ground into fine powders using a Tecator sample mill (Subang, Shanghai, China).Subsequently, C%, δ13C, N% and δ15N of all samples were measured on an isotope ratio mass spectrometer(Isoprime 100, Isoprime Ltd.) coupled to an automatic,online elemental analyzer (Vario ISOTOPE cube). When13C abundance of soil was measured, the soil needs acidification with 1 N HCL for 24 h at room temperature to remove any inorganic carbon. The soil is acid, pH range is 4.77 to 5.44 (Zhou et al. 2019). The reproducibility of the isotope determination was analyzed in multiple runs at several laboratories using an internal standard soil and plant sample from temperate forests. Natural abundance δ15N and δ13C in samples were reported in the conventional delta (δ) notation, with units of per mil(‰). The natural abundance of the N fertilizer(NH4NO3) in this study was 0.21%±0.43‰.

Statistical analysis

All data were tested with the Shapiro-Wilk test and Levene’s tests for normality and homogeneity of variances respectively. Data were analyzed using the general linear mixed model (GLMM) for a split-plot design with Tukey’s multiple comparisons. Nitrogen addition, precipitation reduction, and their interaction were considered fixed factors. Blocks, as random effects, were used to determine the effects of N addition, precipitation reduction, and their interaction on C concentration, δ13C,δ15N, and N concentration values of plant and soil samples. When a significant effect was detected, Tukey’s honestly significant difference test was applied to elucidate treatment differences. Pearson correlation was used to test for correlation among C%, δ13C, N%, and δ15N values of fresh plant foliage, fine roots, and surface soil layer (0-10 cm). All statistical analyses were performed using the SPSS software (version 17.0).

Results

Simulated nitrogen deposition effects

There was a statistical significant effect of different treatments (N addition, precipitation reduction, and combined effect) (p<0.05), tree species (p<0.05),sample types (P<0.05) and their interaction on C concentration, δ13C, N concentration, and δ15N of plant and soil (Table 1). The mean δ15N values of fresh foliage and litterfall were significantly increased by N addition (Table 2, and Tables S1 and S2). When different plant species were considered, under N addition, δ15N values of fresh foliage and litterfall were only significantly higher in Tilia amurensis and Fraxinus mandshurica (p<0.05, Fig. 1). Contrary to foliage and litterfall, the δ15N value of fine roots was significantly lower in N addition compared to that in control (p=0.006, Table 2 and Table S3). There was also a trend for a lower δ15N value of coarse roots in N addition compared to that in control, but this pattern was not statistically significant (p=0.963) (Table 2 and Table S3). Similarly, roots δ15N value of Tilia amurensis and Fraxinus mandshurica, including fine roots and coarse roots, was significantly decreased by N addition (p<0.05, Fig. 1). The δ15N value in the surface soil layer was increased with N addition among all plant species but did not increase in deeper soil layers (Fig. 1). For N concentration, N addition significantly increased the N concentration of litterfall and significantly decreased the N concentration of fine roots and coarse roots (Table 1, Fig. 1, and Tables S1-S7).

Nevertheless, N addition only significantly increased mean δ13C values of fine root and significantly decreased C concentration of fine roots (Table 3). However, δ13C values and C concentration of different plant species and organs showed different responses to N addition (Fig. 2). For Pinus koraiensis, N addition had a significantly positive effect on δ13C values of fresh foliage and had a significant negative effect on δ13C values of litterfall and subsurface soil layer (10-20 cm) (Fig. 2). Nitrogen addition significantly increased δ13C values of litterfall and roots for Tilia amurensis (p<0.05, Fig. 2). Moreover, a significant increase in C concentration of litterfall and δ13C values of surface soil layer (0-10 cm) as well as asignificant decrease in δ13C values of fresh foliage and C concentration of surface soil layer (0-10 cm) for Fraxinus mandshurica were found (p<0.05, Fig. 2).Nitrogen addition also significantly decreased the δ13C value, C and N concentration of fungal sporophore (p<0.05), and mean N return proportion compared with control (Figs. 3 and 4 and Table S8).

Table 1 Effect of different treatments,plant species, sample types and their interaction on C concentration,N concentration,δ13C and δ15N plant and soil

Table 2 The δ15N values (‰) and N concentration (%)of trees and soils as affected by N addition, precipitation reduction and their combination (mean±S.E)

Precipitation reduction effects

Our analysis shows that mean δ13C and δ15N values of foliage significantly increased with precipitation reduction (Tables 2 and 3 and Table S1). Precipitation reduction also significantly increased mean C concentration of foliage, the δ15N value of fine roots, the δ15N value of litterfall, and the C and N concentration of mineral soil layer (20-30 cm), but significantly decreased N concentration of foliage and litterfall,and δ13C value of mineral soil layer (Tables 2 and 3 and Tables S1-S7).

Different plant species and different plant organs showed different responses to precipitation reduction(Figs. 1 and 2). The δ15N value of litterfall, fine roots,coarse roots, and surface soil layer was significantly increased by precipitation reduction in Pinus koraiensis. Precipitation reduction significantly increased the δ15N value of litterfall and surface soil layer but significantly decreased the δ15N value of fine roots and coarse roots and N concentration of foliage and litterfall in Tilia amurensis. There was a significant increase in the δ15N value of foliage and surface soil layer and N concentration of litterfall with precipitation reduction in Fraxinus mandshurica. Besides, precipitation reduction significantly increased the δ13C value of foliage among all plant species. The C concentration of foliage and litterfall was significantly decreased by precipitation reduction in Pinus koraiensis,while the C concentration of foliage and roots was significantly increased by precipitation reduction in

Tilia amurensis. For Fraxinus mandshurica, precipitation reduction had a positive effect on C concentration of foliage and a negative effect on C concentration of fine roots.

Interaction of simulated nitrogen deposition and precipitation reduction

Mean δ15N value of foliage, litterfall, fine roots, and subsurface soil layer (10-20 cm) were significantly increased under the interaction of N addition and precipitation reduction (Tables 1 and 2 and Tables S1-S6). Moreover,the interaction of N addition and precipitation reduction significantly increased the δ13C value of foliage but significantly decreased the δ13C value of litterfall. Interaction of N addition and precipitation reduction had a significantly positive effect on the δ15N value of foliage and litterfall among all plant species. Interaction of N addition and precipitation reduction had a positive effect on the δ15N value of fine root and coarse roots in Pinus koraiensis and Fraxinus mandshurica, but a negative effect in Tilia amurensis. Furthermore, the interaction of N addition and precipitation reduction significantly increased the δ13C value of foliage in Pinus koraiensis, the δ13C value of coarse roots in Tilia amurensis, and the δ13C value of fine roots in Fraxinus mandshurica, but significantly decreased the δ13C value of foliage in Pinus koraiensis.

Correlation analysis

To elucidate the relationship between input and output of C and N values in the plant-soil system, Pearson correlation coefficients were calculated for C%, δ13C,δ15N, and N% in plants foliage, fine roots and surface soil layer (0-10 cm) (Table 4). This showed a strongpositive correlation between the δ15N value of foliage and δ15N value of the surface soil layer (p=0.001) and δ13C value of foliage (p=0.026). The δ15N value of foliage was also negatively correlated with the N concentration of the surface soil layer. A significant negative correlation was found between N concentration of foliage and δ15N value and N concentration of fine roots(p=0.018 and p<0.001, respectively). The δ13C value of foliage was positively correlated with the δ13C value of fine roots and negatively correlated with the δ15N value and N concentration of fine roots. Significant positive correlations were found between C concentration of foliage and δ15N value and N concentration of fine roots.The δ15N value of fine roots was negatively correlated with the δ15N value of the surface soil layer, while the N concentration of fine roots was negatively correlated with C and N concentration of the surface soil layer.

Table 3 The δ13C values (‰) and C concentration(%) of plants and soils as affected by N addition,precipitation reduction and their combination (mean±S.E)

Discussion

Effects of simulated nitrogen deposition on ecosystem C and N cycling

In this study, the N addition (0.21%±0.43‰) with a δ15N abundance close to that of atmospheric N clearly affected the δ15N values of both plants and soils(Table 2 and Fig. 1). We observed a significant increase in surface soil layer δ15N value with N addition, which is consistent with a previous study (Liu et al. 2017). Increasing soil δ15N value may be explained by (1) input of litterfall δ15N value and (2) processes of N mineralization,immobilization, nitrification, denitrification leading to losses of15N-depleted mineral N with N addition(Watzka et al. 2006; Kleinebecker et al. 2014). In addition, plant N uptake could also cause isotope fractionating processes leading to N losses depleted in15N, enriching soil N with the heavier isotope (Evans 2001). The previous study in the experiment site showed that N addition significantly increased N2O emission (Geng et al. 2017), which could partially explain the increase of soil δ15N value.

Although N addition changed the δ15N value of foliage and roots, it did not significantly change their total N concentration among all tree species, which indicated that N in plants had internal homeostasis. Nitrogen addition also significantly increased the mean N concentration of litterfall and decreased the mean N reabsorption of plant foliage, indicating that N addition increased available N in soil and altered plant N uptake strategies.For the δ13C value of plant and soil, we found a positive correlation between foliar δ13C values and foliar δ15N values and N concentrations (Table 4), which indicated that N addition could alter C dynamic by changing N dynamic. However, compared with the change of N and δ15N value,N addition had no or little effect on C and δ13C value of plant and soil. The δ13C value of the plant is relatively stable,and change in plant input would alter the δ13C value of soil. However, the addition of N does not change the total N status in plant and soil,although it could change the form of N uptake by plants, and therefore has no less effect on the C status of plants.Moreover,different species often respond differently to increases in N inputs(Templer et al.2005).In this study,Pinus koraiensis was less sensitive to N addition than others, suggesting C and N cycling in conifers may be less susceptible than broadleaved plants to N deposition in the temperate forest.

Effects of precipitation reduction on ecosystem C and N cycling

In addition to N deposition, change of precipitation also exerts strong controls on15N natural abundance on an ecosystem level (Cheng et al. 2009). Previous studies have shown that plants adapted to precipitation reduction environments are expected to have higher δ15N value (Amundson et al. 2003; Swap et al. 2004; Cheng et al. 2009), which is consistent with our result that δ15N value of foliage and fine root significantly increased with precipitation reduction. The soil is the main N source for plants; thus soil processes play a key role in plant δ15N enrichment under precipitation reduction (Swap et al. 2004). In this study, the δ15N value of the surface soil layer was also increased by precipitation reduction,which may lead to a positive effect on plant δ15N values.Nevertheless, several possible mechanisms may explain this effect of precipitation reduction on the increase of plant-soil δ15N value. Firstly, precipitation reduction may control15N abundance by influencing NH3volatilization (Cheng et al. 2009). Precipitation reduction limited leaching of base cations, and thus might raise surface soil pH value (Cheng et al. 2009). The higher pH values under precipitation reduction permit larger NH3volatilization than in control, thus increasing soil δ15N value. Secondly, gross N mineralization reduced with precipitation reduction (Wu et al. 2011), which may have improved δ15N values of plant-soil in this study. N mineralization discriminates against the heavier N isotope (15N) (Brookshire et al. 2012), which may cause overall enrichment of15N in the soil ecosystem. Consequently, plants have been expected to be more15N enriched. The previous studies at the study site may indirectly confirm that precipitation reduction restrains mineralization by depressing litter decomposition(Zheng et al. 2017). Thus, the reduced N mineralization may be an essential process responsible for the enhancement of plant-soil δ15N values in this study.Thirdly, under precipitation reduction, the increase of plant δ15N values may be interpreted as the reduction in a shift in the N source for plants (Amundson et al.2003). Precipitation reduction may decrease NO3?leaching (Nielsen and Ball 2015), and the enhancement of plant δ15N values could be caused by increasing uptake and assimilation of NO3?in plants. However, the enhancement of plant-soil δ15N values may be caused by an N cycle that becomes more open as precipitation reduction.

Table 4 Pearson correlation coefficients between δ15N, N, δ13C and C in trees foliage,fine roots and surface soil layer(0-10 cm)

According to our hypothesis, we also expected significant changes in δ13C values of plant and soil due to the precipitation reduction. In this study, precipitation reduction significantly increased δ13C values of foliage and fine root but did not significantly change the δ13C values of the soil, which partly supported our hypothesis. Foliage δ13C values vary across changes of precipitation,with the general pattern being that plants in precipitation reduction tend to have more enriched δ13C in their foliage than those in moist environments (Werth and Kuzyakov 2010). Many previous studies reported that plants adapted to precipitation reduction environments are expected to have higher δ13C values (Swap et al.2004; Du et al. 2014; Blessing et al. 2016), which is consistent with our study. This result is mainly due to plant13C discrimination. According to a general notion, under precipitation reduction, plants tend to have low stomatal conductance, high water use efficiency, and low intercellular CO2, which results in decreasing discrimination against13CO2during photosynthesis leading to the increase of foliage δ13C (Blessing et al. 2016). Compared to the increase of foliage δ13C value, the increase in root δ13C value was probably caused by processes occurring during the transport of photoassimilates from foliage to roots (Du et al. 2014). A recent study has also proved that plants maintain the transport of fixed photoassimilates into the roots under moderate drought stress(Hommel et al. 2016), which could lead to an increase in root δ13C value. However, the change of soil δ13C mainly derives from the plant litterfall. The foliage with high concentrations of sugars, starch, protein, organic acids,and cellulose displays high δ13C values (Bowling et al.2008). The increase of foliage δ13C under precipitation reduction may be mainly stored in these organic matters.These organic matters with high δ13C values may decrease with foliage senescence and consequently lead to depleted in litterfall δ13C, especially under the pressure of reduced precipitation. Thus, precipitation reduction did not significantly change in soil δ13C. Moreover, precipitation reduction may reduce the input of plant C into the soil, resulting in a slight decrease in soil C concentration, suggesting that reduced precipitation may lead to a more closed C cycle.

Effect of the interaction between simulated nitrogen deposition and precipitation reduction on ecosystem C and N cycling

In this study, our results directly showed that the interaction of N deposition and precipitation reduction might significantly alter the C and N cycle of plant-soil in a temperate forest. Firstly, in the study, δ15N value of foliage, litterfall, fine root, coarse root and subsurface soil layer (10-20 cm) was higher in the interaction between N addition and precipitation reduction than that in others treatments, providing evidence that the interaction might increase the openness of N cycle compared with N addition alone. Precipitation reduction increased N addition effects mostly indirectly via changing the N status of soil and the physiology of the plants. Mechanistically it is evident that the reduction of precipitation inhibits mineralization and alters soil N status, causing plants to absorb more inorganic N from N addition. The N addition was enriched in15N relative to the control,and decreasing mineralization could also increase the15N enrichment of soil in the study, which may lead to increased plant δ15N value. Furthermore, plants increase water uptake by regulating their osmotic pressure under precipitation reduction conditions, and plants may increase N uptake to regulate osmotic pressure under N addition, which could also increase plant δ15N value.However, increasing plants δ15N values under N addition are always associated with higher levels of water consumption causing plant water stress, being particularly pronounced in precipitation reduction. Previous studies found that N addition could induce an increase in water use efficiency, resulting in a decrease in plant δ13C values (Gong et al. 2011). The result is consistent with the result that foliage δ13C values in the interaction of N addition and precipitation reduction was a slightly lower than that in precipitation reduction alone but still significantly higher than that in the control, which indicated that N addition might increase water use efficiency but precipitation reduction had a stronger effect on δ13C values than N addition, although we have found that interaction of N addition and precipitation reduction may affect ecosystem C and N dynamics.

Conclusions

Our study reveals that N addition, reduced precipitation,and their interaction altered C and N cycles in temperate forest ecosystems. Nitrogen addition and precipitation reduction increased plant-soil δ15N values,suggesting that N cycling becomes more open as precipitation decreases and N deposition increases. Furthermore, precipitation reduction caused an increase in the δ13C values of plants and soils as a result of improved water use efficiency, indicating that C-cycling could become tighter with decreased precipitation. Interaction of N addition and precipitation reduction showed that precipitation reduction might aggravate the effect of N addition on ecosystem N dynamics by regulating N status of soil and the physiology of the plants, while N addition may alleviate the effect of precipitation reduction on ecosystem C cycling by improving water use efficiency. These data give evidence that N deposition and reduced precipitation change the C and N cycles in forest ecosystems; however, the mechanisms that were responsible for these effects, and their intensity were not quantitatively assessed in this study.

Supplementary information

Supplementary informationaccompanies this paper at https://doi.org/10.1186/s40663-020-00257-w.

Additional file 1:Table S1.Effects of N addition, precipitation reduction and their combination on δ15N values (‰) and N concentration (%) of plants and soils among different plants. (HS): Pinus koraiensis; (ZD): Tilia amurensis;and (SQL):Fraxinus mandshurica.Table S2.Effects of N addition, precipitation reduction and their combination on δ13C values (‰)and C concentration (%) of plants and soils among different plants.(HS): Pinus koraiensis;(ZD): Tilia amurensis; and (SQL):Fraxinus mandshurica. Table S3. Effects of N addition, precipitation reduction and their combination on mean δ15N (‰), N.concentration(%), δ13C (‰)and C concentration (%) of foliage. Table S4. Effects of N addition, precipitation reduction and their combination on mean δ15N(‰),N concentration (%), δ13C (‰) and C concentration (%) of litterfall.Table S5. Effects of N addition, precipitation reduction and their combination on mean δ15N (‰),N concentration (%), δ13C (‰)and C concentration (%) of fine roots. Table S6. Effects of N addition,precipitation reduction and their combination on mean δ15N (‰),N concentration (%), δ13C (‰)and C concentration (%) of coarse roots.Table S7. Effects of N addition, precipitation reduction and their combination on mean δ15N (‰),N concentration (%), δ13C (‰)and C concentration (%) of surface soil layer(0-10 cm).Table S8.Effects of N addition, precipitation reduction and their combination on mean δ15N(‰),N concentration (%), δ13C (‰) and C concentration (%) of subsurface soil layer(10-20 cm). Table S9.Effects of N addition,precipitation reduction and their combination on mean δ15N (‰),N concentration (%), δ13C (‰)and C concentration (%) of mineral soil layer(20-30 cm). Table S10.Effects of N addition, precipitation reduction and their combination on mean N return proportion (%).

Acknowledgements

We gratefully acknowledge Professor Xingliang Xu from the Institute of

Geographic Sciences and Natural Resource Research, Chinese Academy of Science for his advice about field experiment design and Professor Yunting Fang and Dr. Feifei Zhu from the Institute of Applied Ecology, Chinese

Academy of Science for their comments and suggestions on an earlier draft of this manuscript.

Authors’contributions

QW and YX designed the study, were awarded funding, supervised data collection and contributed to and edited manuscripts. QW, GY,YX, SH and MZ contributed the whole manuscript preparation and design and wrote the main manuscript text. QW,GY, YX,SH, MZ and WS prepared all figures,YX, HW, HS, BH, WS and QW prepared field experiments, tables and collected literatures. All authors read and approved the final manuscript.

Funding

This research was supported by grants from National Natural Science Foundation of China (Grant Nos: 41773075, 41575137, 31370494, 31170421).

Availability of data and materials

Data are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1School of Life Sciences, Qufu Normal University, Qufu 273165,Shandong,China.2College of Agricultural Resource and Environment, Heilongjiang University, 74 Xuefu Road, Harbin 150080, Heilongjiang,China.3School of Forestry, Northeast Forestry University, 26 Hexing Road, Harbin 150040,Heilongjiang, China.4School of Life Sciences, Henan University, Kaifeng 475004, Henan,China.5Heilongjiang Institute of Construction Technology,999 Xueyuan Road, Harbin 150025, Heilongjiang,China.

Received: 2 November 2019 Accepted: 26 June 2020