Study area
The N-addition experiment was conducted in Mt. Gongga (29°20′−30°20′ N, 101°30′–102°15′ E), which is on the southeastern fringe of the Tibetan Plateau. Mt. Gongga has concentrated N deposition, an ideal area for studying the changes in forest soil microbial CUE at different concentrations of N deposition. The N-addition experiment was established in an evergreen broad-leaf forest (the soil type is yellow brown soil, 2257 m) and a coniferous forest (the soil type is dark brown soil, 2839 m), these two vegetation types are typical and widely distributed in Mt. Gongga. In the evergreen broad-leaf forest, the predominant species is Lithocarpus cleistocarpus and is about 25 ~ 30m high, and the stand density is about 525 trees per hectare. In the coniferous forest, the predominant species is Abies fabri and is about 20 ~ 30m high, and the stand density is about 800 trees per hectare. East Asian monsoon is the dominant climate in Gongga Mountain, and the mean annul precipitation was 1403 mm at 2200 m and 1938 mm at 2900 m. The mean annual temperature was 13–14 ℃ at 2200 m and 3.5-5.0 ℃ at 2900 m. The soil is mostly acid soil with a pH value ranging from 4.5 to 6.
Nitrogen addition experiment
The N addition experiment was established in 2015, and we designed three gradients of N-addition rate (applicated 10 kg N ha− 1 yr− 1, 20 kg N ha− 1 yr− 1and 40 kg N ha− 1 yr− 1, respectively) and a control (0 kg N ha− 1 yr− 1), marked as N10, N20, N40 and N0, respectively. There are 3 uniform plots (1 m×2 m) per treatment for each forest type (24 plots in total). We applied the N fertilizer (urea) once a month from May to October on yearly basis (6 times a year in total), by mixing it in 500 ml deionized water, and then spreading it on the soil surface.
Samples collection
Soil samples were collected in August 2020, and the litters and humus were removed before sampling at each site. The samples were taken with cutting rings randomly (diameter = 5 cm, volume = 100 cm3) ranging from 0–10 cm in each of the replicate for measuring several soil physical indicators (e.g., water content, bulk density, and porosity, n = 3). Given the acreage of the plots, we collected four samples in S-shape from 0–10 cm in each of the replicate and sieved the (2 mm) immediately. Ice bags were used to keep the samples at a lower temperature during the transportation. Each sample was divided into three parts: i) For incubation, the subsample was stored in the frozen state (-20°C) in Ziploc bags until further use; ii) For microbial properties testing, it was stored under 4°C until further use; iii) For soil physical and chemical properties testing, the subsample was air-dried at room temperature.
Analysis of soil properties
Soil water holding capacity was determined by gravimetric method. In brief, weighed the samples after soaking in water for 12 hours (soil samples were in the cutting rings), and then calculate the water-holding capacity. Take a part of each soil sample, weighed it and dry it to a consistent weight at 105℃, then calculate the soil density based on its proportion to the total soil sample and the volume of the cutting ring. A pH meter (Precision and Scientific Corp, China) was use to measuring soil pH (soil: water = 1:2.5). An elemental analyzer (Elementar Vario MACRO cube, Germany) was used to determine soil organic carbon (SOC) and total nitrogen (TN). Soil total phosphorus (TP) was extracted into sulfuric acid (98%) and then measured by an Auto Discrete Analyzer (Smartchem 200, AMS, Italy).
Soil microbial biomass carbon (MBC), phosphorus (MBP) and nitrogen (MBN) were determined by fumigation extraction method. K2SO4 solution (0.5 M) was used to extract dissolved carbon (DOC) and dissolved nitrogen (DON) in fumigated (fumigated with chloroform for 24 h) and non-fumigated soil (Brookes et al. 1985; Vance et al. 1987). The conversion factors of calculating the MBC and MBN were 0.45 and 0.54, respectively (Wu et al. 1990). Dissolved phosphorus (DOP) was extracted into NaHCO3 solution (0.5 M) from fumigated (fumigated with chloroform for 24 h) and non-fumigated soil, and the conversion factor was 0.40 for calculating MBP (Brookes et al. 1985).
Determination of soil microbial CUE
For the understanding of the metabolic process of soil microorganisms, we applied the 13C isotope tracing method to determine soil microbial 13C-respiration and 13C-growth rate following Davey L. Jones et al. (2018). Before incubating, samples were kept for 3 days at 15℃ (average soil temperature at the sites during the growing season) to reactive the microorganisms and the water content of soils was kept at 60% of soil’s water holding capacity to maintain an optimal conditions for microbial activity. At the end of the pre-incubation stage, 5g of the pre-incubated soil was placed into a polypropylene tube (50cm3) and then received a quantified glucose solution (0.2 mg glucose-13C g− 1 soil) that included other reagents (0.1% MgCl2, 0.2% KH2PO4, and 0.1% K2SO4) and sodium nitrate to make it C:N ratio = 40 to mimic the circumstances of natural environment (Wadsö 2009). 2 ml sodium hydroxide solution (NaOH, 1M) was injected into a polypropylene scintillation vial, and placed the scintillation vial upright on the soil surface to catch the respired 13CO2, sealed the tubes and put them into a thermostat with a temperature of 15°C. After 72 h, the polypropylene tubes and scintillation vials were harvested and the remaining 13C-glucose was extracted in ice-cold 1M NaCl (25 ml), the cultivation time was selected based on Glanville et al. (2016) that most of the glucose have been consumed by the microorganisms. The content of 13C in NaCl and NaOH was determined by Delta V Advantage (Thermo, America) and ISOPRIME100 (Elementar, German), respectively.
Microbial 13C-uptake (13Cu) can be estimated as follows:
$${}_{0}{}^{13}{C}_{u}={}_{0}{}^{13}{C}_{t}-{}_{0}{}^{13}{C}_{NaCl}$$
where 13Ct is the total amount of 13C-glucose added to each of the sample, 13CNaCl is the amount of 13C extracted in the 1M NaCl. Microbial CUE for 13C-glucose was estimated by:
$$CUE=\left({}_{0}{}^{13}{C}_{u}-{}_{0}{}^{13}{C}_{r}\right)/{}_{0}{}^{13}{C}_{u}$$
where 13Cr is the amount of 13C trapped by the 1M NaOH which is represents the 13C-respiration. The threshold (TER) of C:N ratio was estimated by (Sterner and Elser 2002; Soong et al. 2020):
$$TER\approx {C:N}_{MB}\times \frac{NUE}{CUE}$$
where C:NMB is microbial biomass C:N ratio, NUE is the microbial N efficiency and it can be estimated by (Zhong et al. 2015):
$$\frac{CUE}{NUE}=\frac{{C:N}_{MB}}{{C:N}_{soil}}$$
Where C:Nsoil is C:N ratio of the soil.
Statistical analysis
The data were checked for normality and homogeneity of variance, and transformed in the case that they were not normally distributed. After that, a one-way ANOVA was performed to test the differences of the properties of soil (soil bulk density, water-holding capacity, pH, SOC, TN and TP) and the microorganisms (MBC, MBN and MBP) in the treatments, which was followed by a Tukey post-hoc test for multiple comparisons. The significant differences were calculated at the confidence interval of 95%. Similarly, the differences of the soil microbial metabolic parameters (13C-uptake, 13C-growth, 13C-respiration and CUE) at varying N-addition rates and soil types was analyzed by using a two-way ANOVA. Furthermore, a liner regression model was established to describe the relationships between soil C:P (or N:P) ratio and microbial CUE. Besides, a redundancy analysis (RDA) was performed to reveal the influences of soil properties on microbial biomass and parameters in metabolic process. The one-way and two-way ANOVA analyses were carried out by using SPSS 25.0 software, and the RDA analysis was conducted using with the vegan package (version 2.5.7) of R software (version 4.4.1).