Effect of thinning intensity on the stem CO2 efflux of Larix principis-rupprechtii Mayr

Kuangji Zhao,Timothy J.Fahey,Xiangzhen Wang,Jie Wang,Fang He,Chuan Fan,Zhongkui Jia*and Xianwei Li*

Abstract

Keywords: Stem CO2 efflux, Thinning intensity, Environmental factors, Nutrient content, Wood structure,Nonstructural carbon, Larix principis-rupprechtii

Introduction

Forest ecosystems play a pivotal role in the global carbon cycle, and the balance between carbon assimilation by forests and its release through respiratory processes is important in regulating the rate of CO2accumulation in the atmosphere (Dixon et al. 1994). Among the largest C flux pathways in forests is emission of CO2from stem surfaces (ES; Ryan et al. 2009); for example, Yang et al.(2016) estimated that annual global ESof forests was 6.7+/-1.1 Pg per year,comparable in magnitude to burning of fossil fuels. Despite considerable recent research, the mechanisms and controls of ESat the individual tree and forest stand level are not fully understood, constraining our ability to predict responses of this flux to changes in climate, forest composition and structure, and other potential drivers (Trumbore et al. 2013; Tu et al. 2017).

Although respiration by living cells in the stem (RS) is the dominant process contributing to ESin trees, interpretation of ESmeasurements is complicated by CO2transport processes either or both to and away from the chamber location,primarily in xylem sap(Ft;Saveyn et al.2008; Salomón et al. 2016b). The principal environmental factor regulating cellular respiration is temperature because the rate of most enzyme-catalyzed biochemical reactions increases exponentially with temperature.This effect is often represented using the Q10convention which specifies the proportional increase in reaction rate (or flux)for every 10°C increase in temperature, i.e. the temperature sensitivity.Both cell maintenance and growth respiration can contribute to the total RS, and maintenance respiration(Rm)is much more temperature sensitive than growth respiration (Rg; Adu-Bredu et al. 1997). Soil nutrient availability also might influence RSbecause cellular respiration often depends upon tissue nutrient (especially N) concentration (Stockfors and Linder 1998; Vose and Ryan 2002).

Regarding internal biotic factors influencing ES, the mass of living tissues directly determines RS(Gruber et al. 2009;Fan et al.2017).Three principal tissues are inner bark,cambium and parenchyma in the sapwood. Typically inner bark is thought to comprise the largest proportion of RS(Gruber et al.2009).Bark thickness can also influence ESthrough its resistance to CO2diffusion from stem to atmosphere(Steppe et al.2007).On a whole-tree basis both ESand RSincreased with gross production and stem growth rate (Yang et al.2016), and phloem girdling studies have also demonstrated that locally higher stem growth above the girdle results in higher ESand RS(Maier et al.2010;Yang et al.2019).However, this relationship is complicated by the local increase in substrate availability above the girdle which may stimulate RSindependently of growth; Maier et al. (2010) suggested that increased RSwith higher substrate supply could reflect upregulation as a means to avoid end-product suppression of photosynthesis.In sum,the environmental and biotic factors affecting ESare complex,resulting in challenges to predicting the response of this important C flux pathway to global change and forest management(Trumbore et al.2013).

Pre-commercial forest thinning is a standard and widely adopted silvicultural practice to enhance the growth in value of a forest stand by directing the wood production into the most valuable trees while maintaining favorable stem form (Kerr and Haufe 2011). Routine thinning is conducted in young, even-aged forest stands of both natural and plantation origin. After thinning the growth of the remaining trees increases and the crowns of individual trees gradually expand increasing gross photosynthesis(Binkley et al. 2004). Opening the canopy increases the penetration of solar radiation and alters the microclimate in the forest(Aussenac 2000).These changes in forest environment, structure, growth and physiological activity would be expected to cause profound responses in ES,but to our knowledge few detailed studies have been reported(Lavigne 1987). Because a large proportion of the world’s forests is managed by thinning, it is important to understand both the response of ESto thinning and the underlying mechanisms.

The objective of this study was to quantify the effect of different levels of pre-commer cial thinning on individual tree and whole-stand ESof Larix principisrupprechtii, a widely-distributed tree species in both natural and plantation forests of northern China (Yao et al.2013). Three plantations were thinned at four different intensities (15%-35% removal) in 2010, and ESwas measured monthly during the warm season in 2013-2015 using a Li-8100 A CO2flux system. Due to the loose bark of Larix principis- rupprechtii, DBH growth rate can not accurately reflect the stem growth rate.Therefore, we use the sapwood growth rate to represent the stem growth rate. There were five primary aims of this study: (1) to clarify the variation in ESwith thinning intensity (TI); (2) to explore the mechanism by which TI influences ES; (3) to evaluate the effect of various influencing factors on ES; (4) to determine the TI-related response and mechanism of temperature sensitivity in ES;and (5) to evaluate the degree of error in the estimations of annual tree-scale and stand-scale CO2efflux if the effect of stand density on ESin the three forest ages is not considered and establish an estimation equation for ESthat incorporates the factors influencing ES.

Materials and methods

Site description

Fig.1 Monthly mean precipitation,air temperature(Ta)and soil temperature(Ts) in the Saihanba National Forest Park(SNFP)(1996-2013).Error bars indicate the standard error,n=9 (Zhao et al.2018b)

The study was conducted at the Saihanba National Forest Park (SNFP, 42°02′-42°36′ N; 116°51′-117°39′ E, altitude 1600-1800 m) of Weichang Manchu and Mongolian Autonomous county in Hebei Province,northern China.The climate is semiarid monsoon climate. The mean annual temperature was -1.5 °C; and average annual rainfall was 433 mm. The precipitation from June to August accounts for 55% of the total annual precipitation (Fig. 1).The soil is classified as gray forest soil (Xi 1994). The soil parent materials are eluvium, alluvium, and saprolite. The thickness of the organic surface horizon is approximately 3-8 cm in all stands(Zhao et al. 2018b).The stand is dominated by L.principis-rupprechtii(Yao et al.2013).

Experimental design

In our study, to be more convincing, the experimental site included 21- (young larch), 25- (immature larch),and 41-year-old (mature larch) artificial forests of L.principis-rupprechtii. Fifteen permanent plots of 400 m2(20 m × 20 m) with essentially the same site and stand conditions were randomly established in larch plantations for each forest age, a total of forty-five plots. According to stand volume calculation, the five TIs were expressed as the percentage of felled trees to the total volume of trees in each plot: 0% (unmanaged, control),15% (weak thinning, TI 15%), 20% (moderate thinning,TI 20%), 30% (strong thinning, TI 30%) and 35% (very strong thinning, TI 35%; Table 1; Shen 2001). In 2010,we randomly implemented these five TIs in fifteen plots for each forest age, with three replicates for each intensity. To minimize edge effects, a 6-m-wide buffer area was established around the inner edge of each plot. The rate of ESwas monitored in three model L. principisrupprechtii trees per standard plot.

Stem CO2 efflux measurements

We used a portable infrared gas analyzer-Li-8100 A (LICOR, Lincoln, Nebraska, USA) with 8100-103 chamber to monitor ESat breast height(approximately 1.3 m)of 45 sample trees.The 8100-103 chamber was sealed over the collar with gasket, and the chamber volume seals with its own gasket when it closes. Use gaskets to minimize leaks and wind effects. The installation steps were as follows:polyvinyl chloride (PVC) collars (with diameters of 7.5 cm) were fixed with environmentally friendly neutral silicone adhesive onto the tested parts of the stems. After drying, the residual silicone adhesive was removed, and the areas between the stem and the PVC collars were sealed. During the measurement, the Li-8100 A chamber was fixed onto the PVC collar with a string. Before the measurement, the seal was checked by exhaling on and around the chamber and the collar to ensure no rapid change in CO2concentration in real time (in a properly sealed chamber, the CO2concentration should increase by no more than 2 ppm) and by the leak test function of the Li8100A. After the measurement, the change in the CO2concentration and the coefficient of variation (CV)value recorded in the Li-8100 A were analysed to ensure the sealing of the Li-8100 A and collar and the accuracy of the data. More details of instrument tightness testing methods are described in Reference (Zhao et al. 2018b).Our previous study found that the daily mean EScould be determined from measurements taken from 9:00 to 11:00 during the daily dynamic monitoring of ESover a one-yearperiod (Zhao et al. 2018a). During the study period, nine days with roughly homogeneous weather conditions were selected from every month, and the ESat each measurement point was measured from 9:00 to 11:00. The air temperature (Ta), wood temperature (Tw, 1 cm below bark) and relative humidity (RH, %) were simultaneously measured with a Li-8100 A instantaneous digital thermometer. The air moisture content (Ma, mmol∙mol-1) was measured using the Li-8100 A system. The temperature and moisture data were acquired simultaneously with the ESmeasurements.The measured area of the stem was calculated by the ellipse area formula: fix a piece of white paper on the surface of the stem to be measured and press it with the treated PVC pipe. The indentation area on the white paper is the measured area of the stem.

Table 1 Stand characteristics for 21-, 25-and 41-year-old Larix principis-rupprechtii plantations

Analysis of wood structure,nutrient content and nonstructural carbon

In August 2016, before the leaves turned yellow, ESvalues were measured, and nine increment cores were taken at a height of 1.3 m from the stems of forty-five trees in each stand age using a wood core sampler.Three cores measured the wood structure by LINTAB 5 tree-ring measurement station (RINNTECH, Heidelberg,Germany) after polished smooth. The growth rate of sapwood is equal to the width of sapwood after thinning divided by the number of annual rings. The other three cores were dried at 70 °C for 48 h, milled and analysed of the nutrient content. The N and P were analyzed by means of an AA3 HR AutoAnalyzer (Seal Analytical Ltd., Southampton, UK); and K were analyzed by Lumina3300 (Aurora Biomed Inc., Vancouver, Canada).The remaining three cores were immediately placed in a cooler and transported to our laboratory. SS and starch were analyzed through an enzymatic assay (McCleary et al. 1997; Yang et al. 2019).

Data analysis

During dormancy, ESis dominated by maintenance respiration, and growth respiration is negligible(Adu-Bredu et al. 1997). In this study, during the cold season (November to April), the respiration rates were too low and the variation too high to meet the standard criterion,which states that the CV of the measurements should be below 2, so we were unable to obtain reliable measurements of ES. Thus, the ESvalues in October were used to represent the ESvalues during dormancy. The period from May to October 2013-2015 was selected for monitoring ES, and the ESvalues during the period from November to April were estimated using the mature tissue method with Taand the change in ESresulting from a 10 °C increase in temperature (Q10) (Amthor 1989).

The assumptions of normality and homogeneity of the variance in the data for each variable were verified using Shapiro-Wilk test and the Levene test. The variables that did not conform to these assumptions were mathematically transformed using logarithms or reciprocal functions. We performed statistical analysis separately for each stand age to statistically analyze the effects of TI and each factor (environmental factors, nutrient content,wood structure and NSC) on ES.

Pearson’s rank correlations and ordinary least square regression (OLS) were used to determine the relationships between the ESand each factor (including the stand density caused by the five TIs). The relative weight method was used to compare the weights of each prediction variable and closely approximated the average increase in R-square obtained by adding a predictor variable across all possible submodels (Johnson 2004;Johnson and Lebreton 2004; Lebreton and Tonidandel 2008). The quantitative relationships between ESand the factors were determined and modeled by all-subsets regression analysis, during which the variance inflation factor (VIF) was used to test the multicollinearity of predicted variables. Based on a comparison of the weights of predictive variables with multicollinearity problems,the predictive variables with smaller weights were excluded, and a new equation was established. Then, the goodness of fit between the old and new models was compared by analysis of variance (ANOVA) and the Akaike information criterion (AIC).The temperature coefficient (Q10) was computed by ESand temperature; the surface area was estimated by the circular truncated cone surface area equation. Estimation of yearly standscale and tree-scale Es and calculation of Q10were done according to Zhao et al. (2018b).

Statistical analysis was performed with the R v.3.5.0 statistical software (R Foundation for Statistical Computing, Vienna, Austria) and at a P ＜0.05 significance level.All figures were created using Origin 8 (OriginLab,Northampton, MA, USA) and the R v.3.5.0 statistical software.

Results

Thinning intensity variation in ES and influencing factors

As shown in Fig. 2, the ESincreased with increasing TI for the 21-, 25- and 41-year-old forests. The ESof each stand was the highest in July of each year; the maximum ESvalues occurred at TI 35% (3.29, 4.57 and 2.98 μmol∙m-2∙s-1) and were 1.54-, 1.94- and 2.89-fold greater than the minimum ESvalues at control (2.14, 2.35 and 1.03 μmol∙m-2∙s-1) in July for the 21-, 25- and 41-yearold stands, respectively. The seasonal ESvariation range in the 21-, 25- and 41-year-old forests increased with increasing TI. The amplitudes of variation at TI 35% were the highest (2.26, 3.57, and 2.55 μmol∙m-2∙s-1), and the amplitudes of variation at control were the lowest (1.74,1.69, and 1.12 μmol∙m-2∙s-1). In the 21- and 25-year-old stands, there were significant differences in ESunder different TIs (P ＜0.05) in July, August and September. In the 41-year-old stands, except for in May and October,there were significant differences in ESunder different TIs (P ＜0.05).

In the 21-year-old forest stands, the variables with significant differences among the five TIs were SR(P=0.012)and ST (P=0.002); the maximum SR values occurred at TI 35% (13.73 cm2∙yr-1), which was significantly higher than other plots(Fig.3);the maximum ST values occurred at TI 15%(19.25 mg∙g-1)(Fig.4).In the 25-year-old forest stands,the variables with significant differences under different TIs were Ta(P=0.020),HS(P=0.015)and SR(P=0.060); the maximum Tavalues occurred at TI 35%(14.44°C) (Fig. 5); HS showed a single peak curve with thinning intensity,and the maximum value appeared at TI 20% (0.64) (Fig. 5). In the 41-year-old forest stands, the variables with significant differences were Ta(P=0.029),BT (P=0.001), SW (P=0.023), HS (P=0.000) and SR(P=0.054); BT showed “V” shape with the increase of thinning intensity, the minimum value is 0.40 cm, which appears in TI 20%(Fig.3).

Effects of environmental factors, nutrient content, wood structure and nonstructural carbon onESin forests of different ages

In the 21-year-old forest, ESwas significantly correlated with Ta(r=0.57, P ＜0.01), SW (r=0.39, P=0.03), N (r=0.51, P ＜0.01) and starch (r=0.37, P=0.04) when the data from all TIs over the 21-year-old forests were included, as shown in Fig. 6a. The results of the model selections based on the all-subsets regression analysis and VIF testing (VIF of SW=12.82,VIF of SR=7.28, SW and SR exhibit multicollinearity)demonstrated that the best model used Ta, density,BT, SW, HS, N and starch to explain the TI variation in ESin the 21-year-old forest stands (R2=0.905)(Table 2).

The Tw(r=0.66, P ＜0.01), Ma(r=0.48, P=0.01) and SS (r=0.49, P=0.01) were significantly correlated with ESin 25-year-old forests (Fig. 6b). Twwas significantly correlated with Ta(r=0.43, P ＜0.01), Ma(r=0.86, P=0.02) and RH (r=0.76, P ＜0.01). The results of the model selections based on the all-subsets regression analysis and VIF testing demonstrated that the best model used Tw, density, BT, HS, SR, N, K and SS to explain the TI variation in ESin the 25-year-old forest stands (R2=0.863) (Table 2).

In the 41-year-old forest, the Ma(r=0.53, P ＜0.01), Ta(r=0.37, P=0.04), density (r = -0.65, P=0.01), N (r=0.46, P=0.01) and SS (r=0.54, P ＜0.01) were significantly correlated with ES(Fig. 6c). RH was significantly correlated with Ma(r=0.87, P ＜0.01). The results of the model selections based on the all-subsets regression analysis and VIF testing demonstrated that the best model used Ta, Ma, density, BT, HS, N and SS to explain the TI variation in ESfor the 41-year-old forests (R2=0.907)(Table 2).

The above models were compared with the equation fitting ESwith all 14 predictive variables by ANOVA and AIC for 21- (AIC: 28.03 VS 39.60, P=0.60), 25- (AIC: -6.94 VS -2.13, P=0.78) and 41-year-old (AIC: -70.24 VS-54.70, P=0.75) forests.

Fig.2 Thinning intensity(TI)variations in the monthly stem CO2 efflux(ES) of 21-(a),25-(b) and 41-year-old(c) Larix principis-rupprechtii forests from May to October 2013-2015.Error bars indicate the standard error,n=3

Q10 change in thinning intensity

Exponential fitting of the annual ESand Taat different TIs was conducted to evaluate the mechanism by which temperature affects ES(Table 3). For the three stand ages, the ESof the trees with lower density was more sensitive to temperature than that of the trees with higher density. In the 25- and 41-year-old forest stands,the ESof TI 35% was the most sensitive to temperature(Q10=3.33 and 3.86). In the 21-year-old forest stands,the ESof the TI 30% stand was the most sensitive to temperature (Q10=3.26), and the Q10of the TI 35%stand was 3.22, which was only 1.23% lower than that of the TI 30% stand. The Q10values of TI 35% in the three forest ages were 1.25, 1.23 and 1.45 times higher than those in control, respectively.

Fig.3 Thinning intensity variations (TI)in the wood structures(bark thickness(BT),sapwood width(SW),width ratio of heartwood to sapwood(HS)and sapwood growth rate (SR))of 21-(a), 25-(b)and 41-year-old (c) Larix principis-rupprechtii forests.Error bars indicate the standard error,n= 9.

Relative weight of impact factor to ES

As shown in Fig. 7, fourteen factors affecting ES, including environmental factors, nutrient content, wood structure and NSC, could explain 90.9%, 87.5% and 89.3% of the ESin the three forest ages,respectively.The five variables with the highest weight for ESin the 21-year-old stand were SS (14.51), starch (14.46), Ta(10.72), SW(9.78) and N (9.37); thus, the main environmental, nutrient content, wood structure and NSC factors that influenced the TI-related variations in ESwere Ta, SW, N and SS, respectively. In the 25-year-old forest, the five highest weighted factors for ESwere Tw(28.37), SS(15.21), Ma(13.81), RH (9.73) and N (8.06); thus, the main environmental, nutrient content, wood structure and NSC factors that influenced the TI-related variations in ESwere Tw, SW (5.03), N and SS, respectively. In the 41-year-old forest, the five highest weighted factors for ESwere SS (17.16), BT (16.48), Ma(16.07), RH (10.53)and Ta(10.37); thus, the main environmental, nutrient content, wood structure and NSC factors that influenced the TI-related variations in ESwere Ma, BT, N (7.57)and SS, respectively.

The amount of CO2 released from different density stands

Fig.4 Thinning intensity(TI)variations in the stem nutrient concentrations(nitrogen concentration(N),phosphorus concentration(P)and potassium concentration(K))and nonstructural carbon(soluble sugar content(SS)and starch content(ST))of 21-(a),25-(b)and 41-year-old(c)Larix principisrupprechtii forests.Error bars indicate the standard error,n=9.

As shown in Table 4,the maximum stand-scale CO2efflux value occurred at TI 35% (10.86±0.37 t CO2∙ha-1∙yr-1)and was 1.45-fold greater than the minimum ESvalue at control (7.47±0.23 t CO2∙ha-1∙yr-1) in the 21-yearold forest. In the 25-year-old forest, the maximum stand-scale CO2efflux occurred at control, with a value of 15.89±0.54 t CO2∙ha-1∙yr-1, and the minimum occurred at TI 30%, with a value of 10.08±0.42 t CO2∙ha-1∙yr-1. In the 41-year-old forest, the maximum stand-scale CO2efflux occurred at TI 35% (5.87±0.31 t CO2∙ha-1∙yr-1), similar to the 21-year-old forest, and the minimum occurred at control,with a value of 3.80±0.10 t CO2∙ha-1∙yr-1. Under five TIs, the stand-scale CO2release of the 25-year-old stand was the highest among the three forest ages, and the CO2release of the 41-year-old stand was the lowest. There were significant differences in the stand-scale CO2release among the five TIs for all three forest ages (P ＜0 01). The tree-scale CO2efflux increased with increasing TI for all three forest ages. The tree-scale CO2release of TI 35% (4.84±0.22, 10.40±0.49 and 12.35±0.39 kg CO2∙plant-1∙yr-1) increased by 168.89%,101.94% and 200.49% compared with that of control(1.80±0.21,5.15±0.27 and 4.11±0.13 kg CO2∙plant-1∙yr-1)for the 21-, 25- and 41-year-old forests, respectively.Additionally, the tree-scale CO2emissions were significantly different among the five TIs for three forest ages(P ＜0.01).

Discussion

Seasonality and TI variation in ES

Fig.5 Thinning intensity(TI)variations in the environmental factors(air temperature (Ta), air moisture content(Ma),relative humidity(RH)and wood temperature (Tw)) of 21-(a),25-(b)and 41-year-old(c) Larix principis-rupprechtii forests.Error bars indicate the standard error,n= 9.

The current study found that ESincreased with TI. Previous studies have presented that thinning could significantly promote the growth of tree DBH (Zhao et al.2016) and that DBH growth could promote ESincrease(Robertson et al. 2010; Yang et al. 2012; Araki et al.2015; Tu et al. 2017), supporting our research results.There are several factors accounting for these complex relationships between ESand TI. Thinning reduces tree density and opens the canopy, thus changing the microclimatic conditions by allowing greater penetration of solar radiation and wind into forest stands (Brèda et al.1995; Aussenac 2000). In conifer forests, low-density stands with more open canopies are usually sunnier,hotter and drier than high-density stands (Meyer et al.2001; Rambo and North 2009; Ma et al. 2010; Martín-Benito et al. 2010;Moreno-Gutierrez et al.2011).Consequently, there are significant differences in temperature under different TIs, thereby supporting the above point of view. ESand temperature are strongly positively correlated (Bown and Watt 2016; Darenova et al. 2018; Salomón et al. 2019). Therefore, temperature is one of the mechanisms by which thinning affects ES. Additionally,thinning may increase photosynthetic production because of the lower competition for light, water and mineral nutrients (López-Serrano et al. 2005; Han et al.2006; Barbour 2007; Forrester et al. 2012; Forrester and Baker 2012). These changes could affect the carbohydrate levels of the remaining trees, and the accumulation of SS could increase the stem respiration rate (Yang et al. 2019). As shown in Fig. 6, thinning increased the SS and promoted the carbon emission rate of stem.Moreover, a decrease in stand density could lead to an enhancement of transpiration (Forrester 2015) and an acceleration of the transport of CO2in the catheter. The ESat a given location includes the CO2generated by trunk and root respiration and transported in the sap flow (Teskey et al. 2008; Bloemen et al. 2014) and the CO2exhaled by trunk cells at the measured location(Araki et al. 2015). Thus, an increase in CO2transport accelerates the outward discharge of CO2in trunk sap flow (Saveyn et al. 2008; Bužková et al. 2015; Salomón et al. 2016b, 2019), eventually leading to an increase in ES. Therefore, we speculate that changes in the microenvironment (temperature and humidity), transpiration tension and NSC content are key to the increase in EScaused by thinning. Future research exploring the relationship between growth respiration and maintenance respiration under different TIs and the above influencing factors could be of great significance for explaining the relationship between thinning and ES.

Fig.6 Correlogram of the correlations among stem CO2 efflux(ES)and the variables in the environmental factors(air temperature(Ta),air moisture content(Ma),relative humidity(RH)and wood temperature(Tw)),nutrient content(nitrogen concentration(N),phosphorus concentration(P)and potassium concentration(K)),wood structure(bark thickness(BT),sapwood width(SW),width ratio of heartwood to sapwood(HS)and sapwood growth rate(SR))and nonstructural carbon(soluble sugar content(SS)and starch content(ST))for 21-(a),25-(b)and 41-year-old(c)Larix principisrupprechtii.Rows and columns have been reordered using the angular order of eigenvectors.The blue color and hashing that goes from lower left to upper right represents a positive correlation between the two variables that meet at that cell.Conversely,a red color and hashing that goes from the upper left to the lower right represents a negative correlation.The darker and more saturated the color,the greater the magnitude of the correlation.The cross sign indicates that the P value of the correlation coefficient of two variables was less than 0.05

The present research shows that Ta, Tw, N and SW are the predominant environmental, nutrient content and wood structure factors that influenced on the TI variations in ESin the 21- and 25-year-old forests. These findings also agree with our earlier observations, which were calculated by a different method (Zhao et al.2018b). Contrary to expectations, the study of ESunder different TIs in the 41-year-old forest showed that the main influencing environmental factor was Ma, while the main influencing wood structural factor was BT. These inconsistencies might be explained as follows: the density of the 41-year-old forest stands was only approximately 25%-33% that of the 21-year-old and 25-yearold stands (Table 1). This lower stand density can result in drought conditions in the forest environment (Moreno-Gutierrez et al. 2011), especially in the thinned 41-year-old forest, which is more water deficient than other stands. When the forest environment is arid, the xylem CO2concentration becomes more sensitive to changes in air humidity, as has been reported for Eucalyptus tereticornis (Salomón et al. 2019), Quercus pyrenaica(Salomón et al. 2016a), and Populus deltoides (Saveyn et al. 2008). These studies demonstrate the dependency of the xylem CO2concentration on water availability;the change in environmental humidity in the forest can greatly increase the content of CO2in the xylem and promote the metabolism of the trunk. Furthermore,thinning can promote an increase in tree diameter (Zhao et al. 2016). Trunk BT increased by 1-2 mm per 10 cm increase in trunk diameter (Paine et al. 2010). The thickening of the bark helps improve trunk water retention(Teskey et al. 2008). However, the change in BT also hinders the diffusion of xylem CO2to the atmosphere(Aubrey and Teskey 2009; Tu et al. 2017), which makes the BT a limiting factor for ESand makes ESmore sensitive to changes in BT. Therefore, we can conclude that the main environmental factors affecting forest ESare Taand Twand that the main tree structure factor is SW.However, when the environmental humidity in a forest decreases to a certain threshold, the xylem CO2concentration of trees sharply decreases (Salomón et al. 2019),resulting in a decrease in forest ES. At this time, the main environmental factor affecting ESchanges from temperature to Ma, and the tree structure factor changes from SW to BT. Therefore, the effects of humidity and BT should be considered when a low-density forest is used to study the changes of ES.

Table 2 Equations for fitting ES with influence factors in 21-, 25- and 41-year-old Larix principis-rupprechtii plantations in August 2016

Table 3 Equations for fitting ES with air temperature and monthly mean temperature sensitivity (Q10) at different thinning intensities for 21-, 25- and 41-year-old Larix principis-rupprechtii forests from May to October 2013–2015

According to the weight ranking of the influencing factors of ES, SS is not only the main influencing factor in NSC but also plays an important role in all the influencing factors of ESin our study. The study of loblolly pine by Maier et al. (2010) showed that SS is an important factor controlling the rates of ESand explains much of the temporal and spatial variation in ESin forest ecosystems, which is consistent with our results. Thinning enhances photosynthesis, and photosynthates are retained in the form of NSC (primarily SS and starch)(Hoch and Körner 2003). Although the photosynthesis and NSC concentration of the needles increased after thinning (Li et al. 2018), the responses of the NSC concentrations of the inner bark and xylem to thinning were not statistically significant. According to theory, the plant NSC pool increases when the assimilated C exceeds the demand for C but decreases when the assimilated C fails to meet demand (Sala et al. 2012). Hence,when the tree growth is slow and the SS is in surplus,the NSC could be stored in the form of starch, whereas when the growth is rapid or conditions are dry, part of the starch could be converted to SS to promote respiration (Rodríguez-Calcerrada et al. 2014; Sevanto et al.2014; Cui et al. 2017), and the allocation of a plant’s internal carbohydrates could present a sink-driven pattern(Lambers et al. 2008). Therefore, we speculate that in a stand with a high TI, trees give priority to assigning more NSC to growth and therefore increase the availability of trunk NSC to occupy more growth space. Although the trunks of trees in stands with a high TI exhibit more NSC, the respiratory consumption of NSC also increases, resulting in lower proportions but constant concentrations of NSC in the trunk. An increase in TI can promote an increase in the DBH of trees, which also confirms this explanation. Therefore, although thinning does not cause a difference in SS, SS, like temperature, is an important factor affecting ES.

Response of Q10 to TI

Q10for the three forest ages ranged from 2.58 to 3.86, comparable with that observed in other species(Yang et al. 2012, 2019; Rodríguez-Calcerrada et al.2014; Tarvainen et al. 2014; Darenova et al. 2018).Close inspection of Table 3 shows that Q10basically increases with increasing TI. In contrast to our study,a previous study suggested that thinning reduces Q10(Lavigne 1987). The differences in the trends probably result from the limited measuring periods (July and September) and TI levels in Lavigne (1987) and from differences among tree species. Thinning can improve the light conditions of a stand, and enhancement of photosynthesis can increase the temperature sensitivity of ES(De Roo et al. 2019). Thinning can increase the transpiration of trees by reducing stand density and then accelerate the flow velocity of trunk sap flow; faster sap flow velocity can increase Q10(Maier et al. 2010; Zhu et al. 2012; Yang et al. 2019). In addition, thinning can enhance the metabolic activity of trunk cells and carbohydrate availability (Lambers et al. 2008; Sala et al. 2012) and increase the SS content in the trunk (Fig. 3), thereby promote the rapid growth of tree diameter, so the improvements in cell metabolic activity (Saveyn et al. 2007), carbohydrate availability (Yang et al. 2019) and soluble sugar content in stem partly explain the increase in Q10. Significantly, thinning not only improved the growth respiration of trees, but also greatly improved the maintenance respiration, and the temperature sensitivity of maintenance respiration was far greater than that of growth respiration (Adu-Bredu et al.1997). We speculate that the increase of maintenance respiration after thinning led to the increase of temperature sensitivity of ES. Therefore, the thinninginduced increase in Q10might be partly due to the potential effects of photosynthesis, sap flow, soluble sugar content, metabolic activity of living cells and carbohydrate availability and the increase of Rm.

Fig.7 Relative weight of predictive variables of 21-(a), 25-(b) and 41-year-old (c) Larix principis-rupprechtii plantations for stem CO2 efflux(ES)

Table 4 Annual stem CO2 efflux for 21-,25- and 41-year-old Larix principis-rupprechtii forests.The means±the standard deviations of the between tree variations are shown(n=3)

CO2 effluxes for each TI

The range of our annual tree-scale estimates of ESfor 21- to 41-year-old L. principis-rupprechtii trees (1.80 to 12.35 kg CO2∙plant-1∙yr-1; Table 4) are similar to those for yellow poplar (Fan et al. 2017), red oak and red maple (Reinmann et al. 2016). Our study shows that if the change in stand density is ignored, the CO2releases of individual trees in the three forest age groups could be underestimated by 168.89%, 101.94% and 200.49%,respectively. Thinning not only promotes the release of CO2via stem respiration but also increases the diameter of the trees and therefore the area of the trunk. Thus,the release of CO2from the trunk increases considerably. On the stand scale, the promoting effect of TI on the CO2release from the 21- and 41-year-old stands was weakened, and the CO2release at the stand scale in 25-year-old stands decreased with the increase in TI. This trend is mainly because although thinning promotes ESand increases the trunk area of individual trees, it reduces the number of trees per unit area. Nevertheless,we know from the results of this study that, regardless of tree scale or stand scale, the release of CO2under different TIs varies significantly. Therefore, in the estimation of stem CO2release, stand density is a factor that cannot be ignored,especially in the calculation of individual tree trunk respiration carbon release. Temperature has a greater effect on ESthan other factors, but Taand Twcan only explain 11.60%-35.56% of R2. Therefore,temperature is not sufficient to quantify ESaccurately,and the effect of other factors, such as stand density, on ESmust be integrated into models simulating ESin forest ecosystems. To estimate the carbon release from the trunk more accurately, we used the ESestimation equation combined with stand density and other factors for reference in Table 2.

Conclusions

In this study, to clarify the effects of TI on ESand Q10of Larix principis-rupprechtii and the related mechanisms and to evaluate the error in CO2emission estimates at the annual tree scale and stand scale without considering the effects of stand density on ES, five TIs were applied to Larix principis-rupprechtii stands with three forest ages of 21, 25 and 41 years, and the ES, environmental,nutritional, wood structure and NSC factors were monitored. The main findings are as follows: (1) Thinning can increase ESby improving the microenvironment(temperature and humidity), transpiration tension and NSC content. (2) SS, like temperature, is a main factor affecting ES. When the stand density is low enough, BT and humidity should be considered in addition to Ta,Tw, SW, N and SS in the evaluation of ES. (3) The temperature sensitivity of ESincreases with increasing TI, and photosynthesis, sap flow, metabolic activity of living cells and carbohydrate availability may be potential influencing factors. (4) If the change in stand density is ignored, the CO2release of individual 21-, 25- and 41-year-old trees could be underestimated by 168.89%,101.94% and 200.49%, respectively. (5) The CO2estimation based on the stem equation was simulated for reference.

Abbreviations

TI: Thinning intensity; ES: Stem CO2efflux; Ta: Air temperature; Tw: Wood temperature; Ma: Air moisture content; RH: Relative humidity; BT:Bark thickness; SW: Sapwood width; HS: Width ratio of heartwood to sapwood;SR: Sapwood growth rate; N: Nitrogen concentration; P: Phosphorus concentration;K: Potassium concentration;SS: Soluble sugar content;ST: Starch content

Authors' contributions

Kuangji Zhao: Design of the study, data collection, data analysis, and manuscript writing. Timothy J. Fahey: Guide the revision of manuscript.Xiangzhen Wang and Jie Wang: Data collection;Fang He and Chuan Fan:Data analysis and polish the manuscript.Zhongkui Jia and Xianwei Li:Supervision of the scientific experiments. The author(s)read and approved the final manuscript.

Funding

This project received funding from National Natural Science Foundation of China (No.31870387),China Scholarship Council.

Declarations

Competing interests

The authors declare that they have no conflict of interest to disclose.

Author details

1Sichuan Province Key Laboratory of Ecological Forestry Engineering on the Upper Reaches of the Yangtze River, College of Forestry, Sichuan Agricultural University,, 611130 Chengdu, China.2Key Laboratory for Silviculture and Conservation of the Ministry of Education, College of Forestry, Beijing Forestry University, 100083 Beijing, China.3Department of Natural Resources,Cornell University, 14853 New York, USA.

Received: 28 December 2020 Accepted: 11 May 2021