LIU Yulin, LI Jiwei, HAI Xuying, WU Jianzhao, DONG Lingbo, PAN Yingjie,SHANGGUAN Zhouping,, WANG Kaibo, DENG Lei,*
1 State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation,Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, China;
2 University of Chinese Academy of Sciences, Beijing 100049, China;
3 Institute of Soil and Water Conservation, Northwest Agriculture and Forestry University, Yangling 712100, China;
4 State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi'an 710075, China
Abstract: Litter and root activities may alter the temperature sensitivity (Q10) of soil respiration. However,existing studies have not provided a comprehensive understanding of the effects of litter and root carbon inputs on the Q10 of soil respiration in different seasons. In this study, we used the trench method under in situ conditions to measure the total soil respiration (Rtotal), litter-removed soil respiration (Rno-litter),root-removed soil respiration (Rno-root), and the decomposition of soil organic matter (i.e., both litter and root removal; RSOM) in different seasons of pioneer (Populus davidiana Dode) and climax (Quercus liaotungensis Mary) forests on the Loess Plateau, China. Soil temperature, soil moisture, litter biomass, fine root biomass, litter carbon, and root carbon were analyzed to obtain the drive mechanism of the Q10 of soil respiration in the two forests. The results showed that the Q10 of soil respiration exhibited seasonality,and the Q10 of soil respiration was higher in summer. The litter enhanced the Q10 of soil respiration considerably more than the root did. Soil temperature, soil moisture, fine root biomass, and litter carbon were the main factors used to predict the Q10 of different soil respiration components. These findings indicated that factors affecting the Q10 of soil respiration highly depended on soil temperature and soil moisture as well as related litter and root traits in the two forests, which can improve our understanding of soil carbon–climate feedback in global warming. The results of this study can provide reference for exploring soil respiration under temperate forest restoration.
Keywords: litter biomass; root carbon; soil respiration; temperate forests; Loess Plateau
A forest ecosystem is the main carbon sequestration body (Kukum?gi et al., 2017) that accounts for 80% of global soil carbon pools (He et al., 2021). Soil CO2efflux presents the largest carbon flux between terrestrial ecosystems and the atmosphere, which is important for the carbon cycle in forest ecosystem (Bond-Lamberty and Thomson, 2010; Ballantyne et al., 2015). The annual emission of CO2from soils into the atmosphere accounts for 5%–20% of the greenhouse gas emitted (Wang et al., 2019), which is the key factor leading to global warming. The temperature sensitivity (Q10) of soil respiration is commonly used to estimate its feedback strength for the temperature increase (Bond-Lamberty et al., 2018; Yan et al., 2019; Ge et al., 2020). Therefore,determining the Q10of soil respiration is essential for predicting soil carbon pools and their potential effects on climate change (Todd-Brown et al., 2013; Bond-Lamberty et al., 2018).
Soil CO2emissions are derived from the metabolism of plant roots (autotrophic respiration, Ra)and microbial communities (heterotrophic respiration, Rh) (Kukum?gi et al., 2017; Zhang et al.,2021), and it is important to determine the contribution of each component to soil respiration and its dependence on multiple environmental factors (Wang et al., 2019; Ge et al., 2020). Vegetation controls soil respiration and its components by influencing litter and root biomass, nutrient availability, soil temperature, and soil moisture (Deng et al., 2014a, 2018; Li et al., 2020). On a large spatial scale, the length of the plant growth season is a chief driving factor for soil respiration (Davidson et al., 2006; Uri et al., 2012; Varik et al., 2013; Machmuller et al., 2018).Moreover, for various seasonal scales, the Q10of soil respiration shows different influencing factors (Davidson et al., 2006; Ferreira et al., 2016). For example, seasonal changes in the root activity trigger environmental factors that have a considerable effect on Ra(Li et al., 2020a); and substrates decompose through heterotrophic soil organisms, whose number and nutrient availability vary with substrates (such as litter) and seasonal variations, respectively (Luan et al.,2011; Ma et al., 2014).
The forest type is another important factor that leads to the change in the Q10values of Raand Rh(Wang et al., 2016; Yu et al., 2017). Studies have shown that various forest types have different effects on the Q10values of total soil respiration (Rtotal), Ra, and Rh(Wang et al., 2013) and that the Q10values of both Raand Rhare significantly positively influenced by forest types (Yu et al.,2017; Li et al., 2020a). However, other studies have reported that no significant difference exists between forest types and the Q10of rhizosphere respiration (Shi et al., 2015). Therefore, there is no consensus on how forest types affect the Q10values of Rtotal, Ra, and Rh. Separately studying the Q10values of Raand Rhchanges is crucial and could enhance our understanding of total soil respiration changes (Gaumant-Guay et al., 2009; Chang et al., 2016).
Populus davidianaDode (pioneer) forests andQuercus liaotungensisMary (climax) forests are typical natural secondary forests species available in the Ziwuling forests region situated on the Chinese Loess Plateau (Deng et al., 2014b). In the two forests, we set up four treatments, namely Rtotal, litter-removed soil respiration (Rno-litter), root-removed soil respiration (Rno-root), and the decomposition of soil organic matter (i.e., both litter and root removal; RSOM), to observe the changes in soil respiration fluxes on different season scales; and then, seasonal dynamic changes and factors influencing the Q10of soil respiration were obtained. The objectives of this study were to: (1) classify Rtotal, Rno-litter, Rno-root, and RSOMto independently quantify the Q10of soil respiration in different seasons in the two forests; (2) explore how the Q10of Rtoatl, Rno-litter,Rno-root, and RSOMvary with forest types; and (3) reveal the main factors affecting the changes in the Q10of soil respiration.
This study was conducted at the Lianjiabian Forest Farmland (36°00′43″–36°06′55″N,108°26′23″–108°41′10″E; 1371–1453 m a.s.l.) in the Heshui County of Gansu Province, China,located in the Ziwuling forests region on the Loess Plateau (Deng et al., 2014b). The study area has a mid-temperate continental monsoon climate. The annual average temperature is 10°C and average annual precipitation is 587 mm (Liu et al., 2020). The study area has a typical loess hilly topography. As the main soil type, loessial soil (Calcic Cambisols) is developed from the primary or secondary loess parent materials, and is distributed evenly at a depth of 50–130 m (Wang et al.,2020). The study area is located in a temperate zone, with a forest density of 80%–95%. This area is suitable for the development of the deciduous broadleaved and temperate coniferous forest species, such asP.davidiana,Betula platyphyllaSuk., andQ.liaotungensis. All these species grow on the top of the hillside, and the slope of the hillside does not exceed 5° (Liu et al., 2020).
According to local forest farm records, vegetation almost disappeared because of human-caused destruction. After the local residents moved out around 1860, the vegetation based on abandoned farmlands naturally recovered, and the vegetation landscapes dominated by the secondary deciduous broadleaved forests formed gradually. Moreover, a relatively complete succession chronosequence of secondary vegetation, including abandoned farmland, herbage,shrub, pioneer forest, climax community, formed in this area. We selectedP.davidianaforests andQ.liaotungensisforests in this study, which are approximately 110 and 160 years old,respectively (Deng et al., 2018; Wang et al., 2020). The main understory plants associated withP.davidianaincludedCarex lanceolataBoott,Artemisia campestrisL.,Ulmus macrocarpaHance,Acer ginnalaMaxim.,Armeniaca sibirica(L.) Lam., andB.platyphylla. The understory plants associated withQ.liaotungensiscomprisedC.lanceolata,A.campestris,P.davidiana, andB.platyphylla.
MatureP.davidianaforests andQ.liaotungensisforests were selected in 2016. The two forests are situated within a distance of less than 2 km from each other. In October 2016, we performed a series of treatments at three plots (each plot is 30 m×30 m) in each forest, in which Rtotal, Rno-litter,Rno-root, and RSOMwere studied with five replicates per treatment. Root exclusion was performed using the trench method (Rey et al., 2002). An area of 1 m×1 m was selected and dug to a depth of 100 cm to remove the root systems. A 100-mesh nylon net was laid around it, and then, in the dug area, the soil was filled again. The soil respiration rate was measured each month from March to November (spring–autumn) in 2018; spring was from March to May, summer was from June to August, and autumn was from September to November.
The soil respiration rate was measured using a soil respiration monitoring system (LI-8100A,Li-COR Inc., Nebraska, USA). From March to November in 2018, at the end or beginning of each month, the soil respiration rate was measured between 09:00 and 14:00 (LST) during the day; the sunny weather was the condition preferred for analyses each time, with 3-min measurements for each chamber. To eliminate the influence of soil disturbance, weeds were cut a day before each measurement. The LI-8100 temperature and humidity sensor was used to simultaneously measure the soil temperature and soil volumetric water content at the depth of 5 cm.
Fine root biomass (RB) was obtained after each soil respiration measurement. In each plot of the two forests, the soil surface mulch was cleaned, a root drill with 9-cm diameter was used to collect the RB at a soil depth of 0–10 cm, 5 drills for each plot were obtained in the diagonal direction, and the dry sieve method was employed to acquire the root system after mixing. At the same time, litter biomass (LB) was collected, and 10 litter collection frames with an area of 1 m×1 m were randomly set up in each sample site to determine LB. Leaves, fallen branches, fallen fruits, fallen flowers, miscellaneous branches, and miscellaneous fruits (worms and faeces, bird droppings, and plant debris) were separated from other parts and taken to the laboratory. The collected samples were dried in an oven at 80°C to obtain a constant weight, and then, were weighed. Litter carbon (LC) and root carbon (RC) were determined through wet digestion by using the K2Cr2O7oxidation method (Nelson and Sommers, 1982).
The Q10values of soil respiration are evaluated using the following equations (Li et al., 2020b):
where Rsis the measured soil respiration rate (μmol/(m2?s)); T is the measured soil temperature(℃); andaandbare the fitted parameters.
Two-way analysis of variance (ANOVA) was used to obtain significant differences (P<0.05)among various treatments and seasons, followed by the least-significant difference (LSD) test,which was performed for the post hoc comparisons of means. Statistical significance was defined asP<0.05. We used linear models to analyse the effects of soil temperature and soil moisture on the Q10values, and utilized the Pearson correlation analysis to determine the relationship of the Q10values with LB, RB, LC, and RC. Structural equation model was performed to identify the possible pathways that allow variables to affect the Q10of soil respiration in the two forests. The structural equation model was constructed using Amos 24.0. Other figures were drawn using SigmaPlot 12.5.
For all the treatments, summer was the season with the highest Q10values inP.davidianaforests andQ.liaotungensisforests (Fig. 1). Among different seasons, the Q10values in summer were considerably different from that in spring and autumn (P<0.05; Table 1). No significant difference was observed in the Q10values between spring and autumn (P>0.05); however, the Q10values in spring were higher than that in autumn. Overall, the Q10values inP.davidianaforests significantly varied among different seasons and treatments (P<0.01), and both the seasons and treatments affected the Q10values inP.davidianaforests (P<0.05; Table 1). The Q10values inP.davidianaforests andQ.liaotungensis forests were 0.455–3.869 and 0.430–3.861, respectively (Fig. 1).
Fig. 1 Seasonal differences among the temperature sensitivity (Q10) of total soil respiration (Rtotal), litter-removed soil respiration (Rno-litter), root-removed soil respiration (Rno-root), and the decomposition of soil organic matter (i.e.,both litter and root removal; RSOM) in Populus davidiana Dode forests (a) and Quercus liaotungensis Mary forests(b). Lowercase letters indicate the variations observed with changing seasons for the same treatment (P<0.05).Uppercase letters indicate the variations observed with varying treatments in the same season (P<0.05). Bars represent standard errors.
Table 1 Analysis of variance for the factors affecting the temperature sensitivity (Q10) of soil respiration with seasons and treatments in Populus davidiana Dode forests and Quercus liaotungensis Mary forests
In all the treatments, the soil respiration rate and soil temperature showed a trend of first increased and then decreased. The maximum values of soil respiration and soil temperature both obtained in July in summer; however, the soil moisture was relatively lower in July (Fig. 2). The Q10values were driven by soil moisture and soil temperature (Fig. 3; Table 2). The combined effect of soil temperature and soil moisture showed an extremely significant correlation with the Q10values (P<0.001; Fig. 3; Table 2). Regardless of the season, soil temperature showed a positive response to the Q10values inP.davidianaforests andQ.liaotungensis forests, with the coefficient of determination (R2) was above 0.80, and the promotion effect of soil temperature on the Q10values were more obvious in summer (Fig. 4a and b). On the whole, soil moisture showed a negative response to the Q10values (Fig. 4c and d).
Fig. 2 Seasonal changes in soil respiration rate, soil temperature, and soil moisture in P. davidiana forests (a, c,and e) and Q. liaotungensis forests (b, d, and f). Bars represent standard errors.
Fig. 3 Responses of the Q10 of soil respiration to the combination of soil temperature (T) and soil moisture (M)with a bivariate non-linear model (Q10=a×M×T+b×M+c×T+d, where a, b, and c represent the coefficients of soil temperature and soil moisture, and d represents the constant) in P. davidiana forests (a, c, and e) and Q.liaotungensis forests (b, d, and f). The specific fitting equation is shown in Table 2.
Table 2 Equations of the Q10 of soil respiration with soil temperature (T) and soil moisture (M) in P. davidiana forests and Q. liaotungensis forests
The LB ofQ.liaotungensisforests was significantly higher than that ofP.davidianaforests(P<0.05), and the LC ofP.davidianaforests was significantly higher than that ofQ.liaotungensisforests (P<0.05; Table 3). The RB and RC ofP.davidianaforests attained the peak in autumn, and the RB ofQ.liaotungensisforests was the highest in summer. No significant difference was observed in RB and RC betweenP.davidianaforests andQ.liaotungensisforests in the same season or different seasons (P>0.05; Table 3).Under different treatments, the Q10values inP.davidianaforests were positively correlated with RB (P<0.05; Fig. 5), and the Q10values inQ.liaotungensisforests was only negatively correlated with LC for Rtotaland Rno-littertreatments (P<0.05; Fig. 5b). The Q10values inP.davidianaforests was negatively correlated with LB in summer (P<0.05; Fig. 5c), and the Q10values inQ.liaotungensisforests was significantly correlated with RB and LC in spring (P<0.05;Fig. 5d).
Fig. 4 Relationships of the Q10 of soil respiration with soil temperature and soil moisture in P. davidiana forests(a and c) and Q. liaotungensis forests (b and d)
Table 3 Litter biomass, litter carbon, fine root biomass, and root carbon in P. davidiana forests and Q.liaotungensis forests
In structural equation model, the two forests showed strong relationships among soil temperature,soil moisture, and the Q10values (Fig. 6). Among the examined plant characteristics, RB and LC can be effectively used to predict the Q10values inP.davidianaforests (Fig. 6a and c). In addition, RB and LC had the greatest influence on the prediction of the Q10values inQ.liaotungensisforests (Fig. 6b and d).
Fig. 5 Pearson correlation of the Q10 of soil respiration with litter biomass (LB), fine root biomass (RB), litter carbon (LC), and root carbon (RC) under different treatments and seasons. (a and b), the correlation of different treatments with LB, RB, LC, and RC; (c and d), the correlation of different seasons with LB, RB, LC, and RC. *,P<0.05.
Fig. 6 Structural equation model for evaluating the direct and indirect effects of T, M, LB, RB, LC, and RC on the Q10 of soil respiration and standardised total effect in P. davidiana forests (a and c) and Q. liaotungensis forests (b and d). The black and grey lines indicate that the correlation is significant and not significant,respectively; the dashed and solid lines represent negative and positive correlations, respectively (P=0.05); and the line thickness represents the magnitude of the path coefficient. CFI, comparative fit index; RMSEA, root mean square error of approximation. *, P<0.05; **, P<0.01; ***, P<0.001.
The Q10of soil respiration showed strong seasonality (Fig. 1; Table 1). In both the forests, the highest Q10values were observed in summer, and the lowest Q10values were obtained in spring and autumn (Fig. 1). A higher Q10value was obtained at higher temperature range mainly because the Q10of unstable carbon mineralisation increased with temperature, and the metabolic efficiency of the microbial community was high (Wang et al., 2013). In spring and autumn, when the temperature was low, the activities of litter, root, and microbial respiration were weak, and the Q10of soil respiration was mainly restricted by biochemical reactions (Wang et al., 2018).Therefore, the soil respiration rate was more sensitive to temperature in summer than that in spring and autumn, which indicated that an increase in temperature leads to an increase in carbon emission.
The Q10of soil respiration between the two forests exhibited no significant difference (P>0.05;Fig. 1). The Q10values were not significantly different in the two forests due to similar soil temperature and soil respiration (Kim et al., 2010). The insignificant differences also resulted from subtle changes in plant variables (root biomass and carbon content) between the two forests(Yan et al., 2013). In this study, the differences in plant characteristics were small, and the most variations between different forests were not obvious (Table 3). Therefore, above factors may be responsible for the insignificant difference between the Q10values inP.davidianaforests and the Q10values inQ.liaotungensisforests.
The Q10of soil respiration enhanced when the litter was included in soil respiration (especially in summer), because the Q10values of Rtotalwere significantly higher than that of Rno-litter. In addition, the Q10values of Rno-rootwere higher than that of Rtotaland Rno-litter(Figs. 1 and 2), which further indicated that the litter enhanced the sensitivity of soil respiration to the temperature increase (Li et al., 2020, 2020c). Some studies have shown that the presence of litter leads to an increase in the substrate input, which increased microbial respiration and substances available for microorganisms (Li et al., 2016; Qu et al., 2018; Goncharova et al., 2020). Other studies also have shown different the Q10values for Raand Rh(Wang et al., 2017; Hu et al., 2018; Zeng et al.,2018). The interactions between litter decomposition and root activities led to an important impact on soil respiration (Subke et al., 2011; Wu et al., 2014), and the mechanism for the contact between the litter and root and soil respiration may be the microbial activities (Subke et al., 2011;Xiao et al., 2017). Due to the variations in availability of the two substrates of litter and root, the microbial decomposition activities led to the difference in soil respiration activities of various components of Raand Rh(D'Acunto et al., 2018; Yan et al., 2018). Therefore, different major drivers of Rno-litter, Rno-root, and RSOMas well as the asynchronous effects of these factors may lead to the variations in the Q10of soil respiration (Wang et al., 2018; Li et al., 2020a, 2020c).
Due to the influence of environmental factors, root and microbial activities were weak in spring and autumn (Wang et al., 2018), and the Q10values were low. In summer, compared with other treatments, Rno-litterexhibited a low Q10value (P<0.05; Fig. 1b), indicating that litter removal strongly affected the Q10values (Kukum?gi et al., 2017). The results from SEM supported this finding (Fig. 6). The results of the Rno-rootand RSOMtreatments are consistent with the temperature quality hypothesis reported in the literature (Suseela et al., 2013; Wu et al., 2014). Although missing roots exudate and the readily decomposable sources of carbon for microbes can influence soil respiration by affecting the decomposition of highly recalcitrant litter and soil organic matter(Wu et al., 2014; Kukum?gi et al., 2017), long trenching experiments can change the soil conditions and microbial biomass and activity in the trenched plots (Kukum?gi et al., 2017).Thus, the aforementioned factors can affect the Q10values among various respiratory components. In addition, studies have reported the similar variations in the Q10values (Jones et al., 2005; Wu et al., 2014), which are consistent with the Rno-litterresults. The opposite results obtained for Rno-rootand Rno-littermay be explained by labile carbon (Jones et al., 2005; Li et al.,2021). The Rno-littertreatment allowed the root rhizosphere activity, thus permitting the input of labile carbon into the soil through photosynthesis (Wu et al., 2014). The Rno-roottreatment allowed litter input; however, a small amount of labile carbon entered the soil because trenching prevented root decomposition (Wu et al., 2014; Li et al., 2020b). The significant differences between Rno-litterand Rtotalwere due to the strong response of carbon input, microbial activities, and soil respiration under the co-existence of litter and root (Li et al., 2021). Therefore, the difference between Rno-litterand Rtotalwas large, especially in the long trench experiments (Jones et al., 2005; Li et al.,2020a). The Q10values for Rno-littertreatment were lower than that for other treatments in summer(Fig. 1), which indicated that the root activity reduced the Q10values (Subke et al., 2011; Wu et al., 2014; Li et al., 2020b).
Soil moisture directly and indirectly affected both autotrophic and heterotrophic respiration through microbial as well as root physiological processes and oxygen and substrate diffusion,respectively (Aaltonen et al., 2017). In this study, soil moisture decreased during the rainy season(Fig. 2), which led to a significant increase in Rtotal, Rno-litter, Rno-root, and RSOM(Aaltonen et al.,2017; Chen et al., 2019). The decreased soil moisture was beneficial to improve the oxygen utilisation rate of the soil (Aaltonen et al., 2017), which led to better plant growth and provided more carbon substrates for root respiration (Chen et al., 2019). Soil moisture exhibited a similar effect on the Q10of soil respiration with higher substrate utilisation, which led to higher Q10values (Brunner et al., 2015; Ferreira et al., 2016). When soil moisture decreased in summer, the soil moisture level was close to the optimal range required for microbial and root activities, which could increase the Q10of soil respiration and components (Ferreira et al., 2016; Chen et al., 2019).
Moreover, soil temperature strongly affected the seasonal variations in the Q10of soil respiration, and the seasonal variation in the Q10values synchronises with soil temperature (Figs.2 and 4). The correlation between soil temperature and the Q10of soil respiration was stronger than that between soil moisture and the Q10of soil respiration (Figs. 4and 5). Soil temperature can be more important than soil moisture in influencing the Q10of soil respiration in temperate forests(Xu et al., 2015). In addition, the Q10of soil respiration was higher in summer than in spring and autumn, indicating that soil respiration was more sensitive to the temperature increase in warmer summers than in cooler spring and autumn (Fig. 2c and d). The Q10values increased with an increase in temperature, which can be confirmed by the strong positive correlation between the Q10of soil respiration and soil temperature (Figs. 4 and 5). The strong positive correlation between the Q10of soil respiration and soil temperature can be explained as an increase in the enzyme activity with an increase in soil temperature (Kunito et al., 2018). The seasonal adaptation of the soil enzyme activity resulted from changes in isoenzymes, which adapt to the corresponding environmental temperature conditions (Di Nardo et al., 2004; Tucker et al., 2017;Guan et al., 2021). Thus, more relevant biological and biochemical activities were observed in summer (Qin et al., 2013).
From the comparison of different litter and root treatments, the Q10of soil respiration can be enhanced without litter removal (Fig. 1). The LC exhibited a higher predictive capacity inP.davidianaforestsandQ.liaotungensisforests (Fig. 6a and b). Changes in LC affected the substrate availability of soil respiration, thereby influencing the Q10of soil respiration (Meyer et al., 2018a; Yan et al., 2019). Litter decomposition can provide carbon that can be easily decomposed by microorganisms (Xu et al., 2018). However, the quantity and quality of litter varied with the forest types, which further affected soil carbon distribution and substrate quality for microbial respiration; therefore, LC exhibited different effects on the Q10values in the two forests (Liu et al., 2017). Simultaneously, other plant characteristics played a role in the different Q10values of respiration components in different forests. Among them, RB exhibited a strong predictive capacity for the Q10values inP.davidianaforests andQ.liaotungensisforests (Fig. 6)because of the large seasonal duration (from spring to autumn) of root metabolism (Li et al.,2020b). Moreover, interactions between the root activity and soil organic matter may affect the variations in the Q10values (Li et al., 2020a). With an increase in root death due to the trenching treatment, the respiration of root fragment decomposition increased, which may last until the second year (October 2016 to November 2018) and lead to an increase in the root respiration proportion (Wu et al., 2014; Kukum?gi et al., 2017). Furthermore, the main predictors of the Q10of soil respiration varied with seasons in different forests because of variations in tree species,ecosystem types, and experimental operation modes (Hu et al., 2016; Meyer et al., 2018b; Chen et al., 2019). Therefore, variations in the Q10of soil respiration were affected by both the biological and abiotic factors (Hu et al., 2016).
The Q10of soil respiration shows considerable seasonality in temperate forests, and the highest Q10values are reached in summer. Among different soil respiration components, the litter enhances soil respiration. Compare with the litter effect, root have a lower CO2emission. In addition,P.davidianaforests andQ.liaotungensisforests have a similar Q10of soil respiration.The interactions and synergistic effects of soil temperature, soil moisture, litter biomass, and fine root biomass play an important role in influencing the Q10of soil respiration in the two forests.These results can explain how soil temperature, soil moisture, litter biomass, fine root biomass,litter carbon, and root carbon influence the Q10of soil respiration and help us better understand how temperate forest systems respond to long-term climate dynamics.
Acknowledgements
This study was sponsored by the National Natural Science Foundation of China (41877538, 41671511), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB40000000, XDA23070201), and the Funding of Special Support Plan of Young Talents Project of China and National Forestry and Grassland Administration in China (20201326015).