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

Background: Stem CO2 efflux (ES) plays a critical role in the carbon budget of forest ecosystems. Thinning is a core practice for sustainable management of plantations. It is therefore necessary and urgent to study the effect and mechanism of thinning intensity (TI) on ES. Methods: In this study, five TIs were applied in Larix principis-rupprechtii Mayr 21-, 25-, and 41-year-old stands in North China in 2010. Portable infrared gas analyzer (Li-8100 A) was used to measure ES and its association with environmental factors at monthly intervals from May to October in 2013 to 2015. In addition, nutrients, wood structure and nonstructural carbon (NSC) data were measured in August 2016. Results: The results show that ES increased with increasing TI. The maximum ES values occurred at a TI of 35 % (3.29, 4.57 and 2.98 μmol∙m∙s) and were 1.54-, 1.94and 2.89-fold greater than the minimum ES value in the CK stands (2.14, 2.35 and 1.03 μmol∙m∙s) in July for the 21-, 25and 41-year-old forests, respectively. The ES of the trees in low-density stands was more sensitive to temperature than that of the trees in high-density stands. Soluble sugars (SS) and temperature are the main factors affecting ES. When the stand density is low enough as 41-year-old L. principis-rupprechtii forests with TI 35 %, bark thickness (BT) and humidity should be considered in addition to air temperature (Ta), wood temperature (Tw), sapwood width (SW), nitrogen concentration (N) and SS in the evaluation of ES. If a change in stand density is ignored, the CO2 released from individual 21-, 25and 41-year-old trees could be underestimated by 168.89 %, 101.94 % and 200.49 %, respectively. CO2 release was estimated based on the stem equation in combination with the factors influencing ES for reference. Conclusions: We suggest that it is not sufficient to conventional models which quantify ES only by temperature and that incorporating the associated drivers (e.g. density, SS, SW and N) based on stand density into conventional models can improve the accuracy of ES estimates.


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 CO 2 accumulation in the atmosphere (Dixon et al. 1994). Among the largest C flux pathways in forests is emission of CO 2 from stem surfaces (E S ; Ryan et al. 2009); for example, Yang et al. (2016) estimated that annual global E S of 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 E S at 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 (R S ) is the dominant process contributing to E S in trees, interpretation of E S measurements is complicated by CO 2 transport 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 Q 10 convention 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 R S , and maintenance respiration (Rm) is much more temperature sensitive than growth respiration (Rg; Adu-Bredu et al. 1997). Soil nutrient availability also might influence R S because cellular respiration often depends upon tissue nutrient (especially N) concentration (Stockfors and Linder 1998;Vose and Ryan 2002).
Regarding internal biotic factors influencing E S , the mass of living tissues directly determines R S (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 R S (Gruber et al. 2009). Bark thickness can also influence E S through its resistance to CO 2 diffusion from stem to atmosphere ). On a whole-tree basis both E S and R S increased 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 E S and R S (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 R S independently of growth; Maier et al. (2010) suggested that increased R S with higher substrate supply could reflect upregulation as a means to avoid end-product suppression of photosynthesis. In sum, the environmental and biotic factors affecting E S are 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 E S , 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 E S to 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 E S of 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 E S was measured monthly during the warm season in 2013-2015 using a Li-8100 A CO 2 flux 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 E S with thinning intensity (TI); (2) to explore the mechanism by which TI influences E S ; (3) to evaluate the effect of various influencing factors on E S ; (4) to determine the TI-related response and mechanism of temperature sensitivity in E S ; and (5) to evaluate the degree of error in the estimations of annual tree-scale and stand-scale CO 2 efflux if the effect of stand density on E S in the three forest ages is not considered and establish an estimation equation for E S that incorporates the factors influencing E S .

Site description
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 m 2 (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 E S was monitored in three model L. principisrupprechtii trees per standard plot.

Stem CO 2 efflux measurements
We used a portable infrared gas analyzer-Li-8100 A (LI-COR, Lincoln, Nebraska, USA) with 8100 − 103 chamber to monitor E S at 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 CO 2 concentration in real time (in a properly sealed chamber, the CO 2 concentration should increase by no more than 2 ppm) and by the leak test function of the Li8100A. After the measurement, the change in the CO 2 concentration 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 E S could be determined from measurements taken from 9:00 to 11:00 during the daily dynamic monitoring of E S over a one-year period (Zhao et al. 2018a). During the study period, nine days with roughly homogeneous weather conditions were selected from every month, and the E S at each measurement point was measured from 9:00 to 11:00. The air temperature (T a ), wood temperature (T w , 1 cm below bark) and relative humidity (RH, %) were simultaneously measured with a Li-8100 A instantaneous digital thermometer. The air moisture content (M a , mmol•mol -1 ) was measured using the Li-8100 A system. The temperature and moisture data were acquired simultaneously with the E S measurements. 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.

Analysis of wood structure, nutrient content and nonstructural carbon
In August 2016, before the leaves turned yellow, E S values 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 Lu-mina3300 (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, E S is 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 E S . Thus, the E S values in October were used to represent the E S values during dormancy. The period from May to October 2013-2015 was selected for monitoring E S , and the E S values during the period from November to April were estimated using the mature tissue method with T a and the change in E S resulting from a 10°C increase in temperature (Q 10 ) (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 E S .
Pearson's rank correlations and ordinary least square regression (OLS) were used to determine the relationships between the E S and 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 E S and 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 (Q 10 ) was computed by E S and 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 Q 10 were 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.

Thinning intensity variation in E S and influencing factors
As shown in Fig. 2, the E S increased with increasing TI for the 21-, 25-and 41-year-old forests. The E S of each stand was the highest in July of each year; the maximum E S values 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 E S values 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 E S variation 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 E S under 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 E S under 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 cm 2 •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 T a (P = 0.020), HS (P = 0.015) and SR (P = 0.060); the maximum T a values 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 T a (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 onE S in forests of different ages In the 21-year-old forest, E S was significantly correlated with T a (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 T a , density, BT, SW, HS, N and starch to explain the TI variation in E S in the 21-year-old forest stands (R 2 = 0.905) ( Table 2).
The T w (r = 0.66, P < 0.01), M a (r = 0.48, P = 0.01) and SS (r = 0.49, P = 0.01) were significantly correlated with E S in 25-year-old forests (Fig. 6b). T w was significantly correlated with T a (r = 0.43, P < 0.01), M a (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 T w , density, BT, HS, SR, N, K and SS to explain the TI variation in E S in the 25-year-old forest stands (R 2 = 0.863) ( Table 2).
In the 41-year-old forest, the M a (r = 0.53, P < 0.01), T a (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 E S (Fig. 6c). RH was significantly correlated with M a (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 T a , M a , density, BT, HS, N and SS to explain the TI variation in E S for the 41-year-old forests (R 2 = 0.907) ( Table 2).

Q 10 change in thinning intensity
Exponential fitting of the annual E S and T a at different TIs was conducted to evaluate the mechanism by which temperature affects E S (Table 3). For the three stand ages, the E S of 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 E S of TI 35 % was the most sensitive to temperature (Q 10 = 3.33 and 3.86). In the 21-year-old forest stands, the E S of the TI 30 % stand was the most sensitive to temperature (Q 10 = 3.26), and the Q 10 of the TI 35 % stand was 3.22, which was only 1.23 % lower than that of the TI 30 % stand. The Q 10 values of TI 35 % in the three forest ages were 1.25, 1.23 and 1.45 times higher than those in control, respectively.
Relative weight of impact factor to E S As shown in Fig. 7, fourteen factors affecting E S , including environmental factors, nutrient content, wood structure and NSC, could explain 90.9 %, 87.5 % and 89.3 % of the E S in the three forest ages, respectively. The five variables with the highest weight for E S in the 21-year-old stand were SS (14.51), starch (14.46), T a (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 E S were T a , SW, N and SS, respectively. In the 25-year-old forest, the five highest weighted factors for E S were T w (28.37), SS (15.21), M a (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 E S were T w , SW (5.03), N and SS, respectively. In the 41-year-old forest, the five highest weighted factors for E S were SS (17.16), BT (16.48), M a (16.07), RH (10.53) and T a (10.37); thus, the main environmental, nutrient content, wood structure and NSC factors that influenced the TI-related variations in E S were M a , BT, N (7.57) and SS, respectively.
The amount of CO 2 released from different density stands Table 4, the maximum stand-scale CO 2 efflux value occurred at TI 35 % (10.86 ± 0.37 t CO 2 •ha -1 •yr -1 ) and was 1.45-fold greater than the minimum E S value at control (7.47 ± 0.23 t CO 2 •ha -1 •yr -1 ) in the 21-yearold forest. In the 25-year-old forest, the maximum stand-scale CO 2 efflux occurred at control, with a value of 15.89 ± 0.54 t CO 2 •ha -1 •yr -1 , and the minimum occurred at TI 30 %, with a value of 10.08 ± 0.42 t CO 2 •ha -1 •yr -1 . In the 41-year-old forest, the maximum stand-scale CO 2 efflux occurred at TI 35 % (5.87 ± 0.31 t CO 2 •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 CO 2 •ha -1 •yr -1 . Under five TIs, the stand-scale CO 2 release of the 25-year-old stand was the highest among the three forest ages, and the CO 2 release of the 41-year-old stand was the lowest. There were significant differences in the stand-scale CO 2 release among the five TIs for all three forest ages (P < 0 01). The tree-scale CO 2 efflux increased with increasing TI for all three forest ages. The tree-scale CO 2 release of TI 35 % (4.84 ± 0.22, 10.40 ± 0.49 and 12.35 ± 0.39 kg CO 2 •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 CO 2 •plant -1 •yr -1 ) for the 21-, 25-and 41-year-old forests, respectively. Additionally, the tree-scale CO 2 emissions were significantly different among the five TIs for three forest ages (P < 0.01).

Seasonality and TI variation in E S
The current study found that E S increased 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 E S increase  There are several factors accounting for these complex relationships between E S and 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. E S and 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 E S . 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;  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 CO 2 in the catheter. The E S at a given location includes the CO 2 generated by trunk and root respiration and transported in the sap flow Bloemen et al. 2014) and the CO 2 exhaled by trunk cells at the measured location (Araki et al. 2015). Thus, an increase in CO 2 transport accelerates the outward discharge of CO 2 in trunk sap flow Bužková et al. 2015;Salomón et al. 2016bSalomón et al. , 2019, eventually leading to an increase in E S . Therefore, we speculate that changes in the microenvironment (temperature and humidity), transpiration tension and NSC content are key to the increase in E S caused 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 E S .
The present research shows that T a , T w , N and SW are the predominant environmental, nutrient content and wood structure factors that influenced on the TI variations in E S in 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 E S under different TIs in the 41-year-old forest showed that the main influencing environmental factor was M a , 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 41year-old forest, which is more water deficient than other stands. When the forest environment is arid, the xylem CO 2 concentration 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 . These studies demonstrate the dependency of the xylem CO 2 concentration on water availability; Fig. 6 Correlogram of the correlations among stem CO 2 efflux (E S ) and the variables in the environmental factors (air temperature (T a ), air moisture content (M a ), relative humidity (RH) and wood temperature (T w )), 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 change in environmental humidity in the forest can greatly increase the content of CO 2 in 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 . However, the change in BT also hinders the diffusion of xylem CO 2 to the atmosphere (Aubrey and Teskey 2009;Tu et al. 2017), which makes the BT a limiting factor for E S and makes E S more sensitive to changes in BT. Therefore, we can conclude that the main environmental factors affecting forest E S are T a and T w and that the main tree structure factor is SW. However, when the environmental humidity in a forest decreases to a certain threshold, the xylem CO 2 concentration of trees sharply decreases (Salomón et al. 2019), resulting in a decrease in forest E S . At this time, the main environmental factor affecting E S changes from temperature to M a , 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 E S . According to the weight ranking of the influencing factors of E S , SS is not only the main influencing factor in NSC but also plays an important role in all the influencing factors of E S in our study. The study of loblolly pine by Maier et al. (2010) showed that SS is an important factor controlling the rates of E S and explains much of the temporal and spatial variation in E S in 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 E S .
Response of Q 10 to TI Q 10 for the three forest ages ranged from 2.58 to 3.86, comparable with that observed in other species (Yang et al. 2012(Yang et al. , 2019Rodríguez-Calcerrada et al. 2014;Tarvainen et al. 2014;Darenova et al. 2018). Close inspection of Table 3 shows that Q 10 basically increases with increasing TI. In contrast to our study, a previous study suggested that thinning reduces Q 10 (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 E S (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 Q 10 (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 Table 3 Equations for fitting E S with air temperature and monthly mean temperature sensitivity (Q 10 ) at different thinning intensities for 21-, 25-and 41-year-old Larix principis-rupprechtii forests from May to October 2013-2015  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 ), carbohydrate availability (Yang et al. 2019) and soluble sugar content in stem partly explain the increase in Q 10 . 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 E S . Therefore, the thinninginduced increase in Q 10 might 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.

CO 2 effluxes for each TI
The range of our annual tree-scale estimates of E S for 21-to 41-year-old L. principis-rupprechtii trees (1.80 to 12.35 kg CO 2 •plant -1 •yr -1 ; Table 4) are similar to those for yellow poplar , red oak and red maple (Reinmann et al. 2016). Our study shows that if the change in stand density is ignored, the CO 2 releases 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 CO 2 via stem respiration but also increases the diameter of the trees and therefore the area of the trunk. Thus, the release of CO 2 from the trunk increases considerably. On the stand scale, the promoting effect of TI on the CO 2 release from the 21-and 41-year-old stands was weakened, and the CO 2 release at the stand scale in 25year-old stands decreased with the increase in TI. This trend is mainly because although thinning promotes E S and 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 CO 2 under different TIs varies significantly. Therefore, in the estimation of stem CO 2 release, 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 E S than other factors, but T a and T w can only explain 11.60 %-35.56 % of R 2 . Therefore, temperature is not sufficient to quantify E S accurately, and the effect of other factors, such as stand density, on E S must be integrated into models simulating E S in forest ecosystems. To estimate the carbon release from the trunk more accurately, we used the E S estimation equation combined with stand density and other factors for reference in Table 2.

Conclusions
In this study, to clarify the effects of TI on E S and Q 10 of Larix principis-rupprechtii and the related mechanisms and to evaluate the error in CO 2 emission estimates at the annual tree scale and stand scale without considering the effects of stand density on E S , five TIs were applied to Larix principis-rupprechtii stands with three forest ages of 21, 25 and 41 years, and the E S , environmental, nutritional, wood structure and NSC factors were monitored. The main findings are as follows: (1) Thinning can increase E S by improving the microenvironment (temperature and humidity), transpiration tension and NSC content.
(2) SS, like temperature, is a main factor affecting E S . When the stand density is low enough, BT and humidity should be considered in addition to T a , T w , SW, N and SS in the evaluation of E S . (3)

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