Responses of Fine Roots at Different Soil Depths to Different Thinning Intensities in a Secondary Forest in the Qinling Mountains, China

Simple Summary Fine roots of different sizes make critical contributions to carbon stocks and terrestrial productivity, but the changed characteristics of fine roots with different diameters at different soil depths under thinning disturbances remain poorly understood. In our study, we aimed to elucidate the response characteristics of fine roots with different diameters to thinning intensities at different soil depths, and to explore the driving mechanism of the change in the fine-root characteristics. We found that higher thinning intensities negatively affected 0.5–1 mm and 1–2 mm fine-root biomass, while the <0.5 mm fine-root characteristics fluctuated with increasing thinning intensities. Our results suggest that the thinning intensity exhibits varied influential mechanisms on the changed characteristics of <0.5 mm fine roots and thicker fine roots (0.5–2 mm). Collectively, our findings provide important insights into the effects of forest management on changes in fine-root characteristics, and supplement meaningful data on fine-root productivity to improve the parameterization of future ecological models. Abstract Fine roots make critical contributions to carbon stocks and terrestrial productivity, and fine roots with different diameters exhibit functional heterogeneity. However, the changed characteristics of fine roots with different diameters at different soil depths following thinning disturbances are poorly understood. We investigated the biomass, production, mortality and turnover rate of <0.5 mm, 0.5–1 mm, and 1–2 mm fine roots at 0–20 cm, 20–40 cm, and 40–60 cm soil depths under five thinning intensities (0%, 15%, 30%, 45%, and 60%) in a secondary forest in the Qinling Mountains. The biomass, production and turnover rate of <0.5 mm fine roots fluctuated with increasing thinning intensities, while the 0.5–1 mm and 1–2 mm fine-root biomass significantly decreased. The thinning intensities had no effects on the fine-root necromass or mortality. The change in the fine-root characteristics in deeper soils was more sensitive to the thinning intensities. The principal component analysis results showed that increased <0.5 mm fine-root biomass and production resulted from increased shrub and herb diversity and biomass and decreased soil nutrient availability, stand volume, and litter biomass, whereas the 0.5–1 mm and 1–2 mm fine-root biomass showed the opposite trends and changes. Our results suggest that different thinning intensities exhibit varied influential mechanisms on the changed characteristics of fine roots with different diameters.


Introduction
Thinning operations, a key forest management practice, are widely employed to maintain secondary forest ecosystem stability and promote stand productivity [1]. Thinning decreases stand density, regulates stand structure, reduces self-thinning mortality and wildfire risk, accelerates the rapid growth of remaining trees, increases resistance to natural disasters, stimulates ecosystem nutrient cycling processes, and improves forest health [2][3][4][5]. Additionally, reasonable thinning can also improve the stand environment and promote the regeneration of understory vegetation, and transform energy consumption caused by competition into effective productivity [6,7]. Nevertheless, little is known about the impact of thinning intensity on belowground productivity, although the belowground parts of vegetation make critical contributions to terrestrial productivity and carbon stock [3]. This is especially true for fine roots (Ø ≤ 2 mm). Fine roots account for less than 5% of all forest biomass but for 30-50% of total net primary production [8], and for approximately 50-70% of the soil C flux in forest ecosystems [9]. Unlike the sampling of other aboveground organic parts, the sampling of plant roots at the stand level is destructive, laborious and technically challenging, leading to incomplete estimations and understanding of fine-root processes [10].
Abiotic and biotic factors are the main drivers of change in fine-root characteristics (biomass, production, necromass, mortality and the turnover rate) [11,12]. Thinning treatments can directly modify the physical structure of trees and understory vegetation, and indirectly influence belowground environmental conditions (e.g., soil fertility and microclimate), potentially impacting changes in fine-root characteristics [13,14]. However, a few studies have reported contrasting results regarding thinning. For instance, some studies have shown that the biomass and production of fine roots generally increase following thinning due to improved soil nutrient ability and soil moisture levels [15,16]. In contrast, fine-root biomass and production responses to thinning have been reported in other studies, including negative and negligible responses [17,18]. These inconsistencies may result from the variations in sampling heterogeneity, plant species composition and climatic conditions. Undoubtedly, the exploration of the underlying mechanisms of change in fine-root characteristics' responses to thinning intensity will provide a theoretical basis for the further understanding of productivity and nutrient cycling in terrestrial ecosystems.
Fine roots are traditionally defined as roots with a diameter of <2 mm, though this definition has recently been debated [19]. Even in this small-diameter category (Ø ≤ 2 mm), fine roots of different sizes exhibit functional heterogeneity due to their different physiological activities [20,21]. Fine roots of <0.5 mm (very fine) are mainly responsible for the acquisition and absorption of resources in the soil, and 0.5-2 mm fine roots (thicker) are responsible for resource transfer and storage [22,23]. Very fine roots are typically more dynamic than thicker roots, and are extremely sensitive to environmental change, as the ratio of nonstructural to structural mass is much higher in smaller roots [15,24]. For instance, the responses of fine roots to soil nutrient and water changes are diameter dependent [25,26]. Furthermore, Ma et al. (2013) reported that very fine roots have a higher turnover rate and productivity levels than thicker roots following thinning.
Plants may also adjust their root systems to use resources at different soil depths [27]. When they are stimulated by more resource competition in the surface soil, plants can adjust their roots at deep soil depths to obtain more soil resources [10]. When competition for surface soil resources decreases, the exploration of deep-layer soil resources is also reduced, and the plant root system will again shift to the utilization of surface resources [28]. Thinning alters the original vegetation configuration, resulting in a change in resource competition among species that in turn leads to varied resource use strategies being employed by fine roots at different soil depths [29]. Light-and high-intensity thinning slightly or significantly reduce resource competition, leading deep-soil resource exploration to decrease [16]. Under a suitable thinning treatment, understory species diversity increases and resource competition intensifies, and fine roots gradually explore deep-soil resources, resulting in increased fine-root biomass and production at deeper soil depths [14]. It can be seen that there may be larger fluctuations in the characteristics of fine roots at the deep soil under disturbance. Thus, collecting samples at deeper soil depths is needed in order to fully quantify the change in fine-root characteristics and spatial distribution variability [30], and to clarify the resource acquisition strategies of fine roots under different thinning treatments.
In the present study, we determined the biomass, necromass, production, mortality, and turnover rate characteristics of <0.5 mm, 0.5-1 mm, and 1-2 mm fine roots at 0-20 cm, 20-40 cm and 40-60 cm soil depths under five thinning intensities (0%, 15%, 30%, 45%, and 60% of the stand volume removed) in secondary forests. Our hypotheses are as follows: (a) changes in very fine-root characteristics are more sensitive to thinning intensities than thicker fine roots due to their greater responsiveness to changing environments; (b) fineroot necromass and mortality increase when thinning intensities increase, as cutting trees would lead to the death of root systems [31,32]; and (c) changes in fine-root characteristics in deep soils are more sensitive to thinning intensities than those in shallow soils, and exhibit fluctuations because plants must flexibly adjust the root system to use resources at deep soil layers under the changing pressures of soil resource competition. Furthermore, in order to understand the driving effects of abiotic and biotic conditions on the biomass, production, necromass, mortality and turnover rate of fine roots, we explored the linkages between fine-root characteristics, soil properties and stand characteristics.

Study Area
The experiment was implemented at the Qinling National Forest Ecosystem Research Station (33 • 18 -33 • 28 N, 108 • 21 -108 • 39 E) located on the Huoditang Experimental Forest Farm of Northwest A&F University, in Ningshaan County, Shaanxi Province, China. This study area covers an area of 22.25 square kilometers, and has a subtropical humid montane climate with altitudes ranging from 1900 to 2200 m and a mean slope of approximately 35 • . Across the site, Cambisols, Umbrisols and Podzols (FAO) are the dominant soil types, reaching a mean soil depth of 50 cm [33]. The annual average temperature and humidity are approximately 10.5 • C and 77%, respectively. The annual mean precipitation level is 1000 mm, with 50% occurring from July to September. Local plants have a 177-day growth period, while the average frost-free period of the study area is approximately 199 days [34]. During the 1960s and 1970s, forests around the Qinling National Forest Ecosystem Research Station underwent extensive logging. After decades of regeneration, the forest coverage rate reached 93.8%, and secondary growth dominated the area. The dominant tree species are Quercus aliena var. acutiserrata, Quercus variabilis, Pinus armandii, Betula albosinensis, Picea asperata, and Populus davidiana. Shrubs (i.e., Lonicera tragophylla, Cerasus stipulacea, and Symplocos paniculata) and herbs (i.e., Lysimachia christinae, Rubus parvifolius, Saussurea mutabilis, and Rubia cordifolia) occupy the understory space [35][36][37].

Experimental Design and Treatments
For the purpose of promoting high-quality trees, forest farm staff cut down disturbing trees in the secondary forests according to the target tree operation system. Whole-tree harvesting was used during the felling, which removed most of the aboveground tree parts (stems, needles, branches, and twigs) to a greater extent than conventional stemonly harvesting [38,39]. These thinning treatments in secondary forest stands (dominated by Pinus armandii, Betula albosinensis, and Picea asperata) were carried out from July to September 2013. All of the selected plots were of the same stand age (35 years), occupied a similar topography, and had no history of fertilization. A randomized complete block design was used in the study. Five 20 × 20 m plots were randomly installed within each secondary forest block. Based on the target tree operation system, the thinning intensities included: (1) no thinning (CK), (2) the 15% removal of the stand volume (T1), (3) the 30% removal of the stand volume (T2), (4) the 45% removal of the stand volume (T3), and (5) the 60% removal of the stand volume (T4). Given that the study area had undergone large-scale logging, the forest community is now at a young stage. Therefore, combined with the Chinese forestry industry standard, a buffer zone of 5 m was set to avoid potential edge effects. All of the harvesting debris was removed from the plots. Each of the five thinning treatments was replicated into four blocks, totaling 20 sampling plots (five thinning intensities × four blocks). The layout of the experimental design is shown in Figure 1 (including the blocks and plots).
the Chinese forestry industry standard, a buffer zone of 5 m was set to avoid potential edge effects. All of the harvesting debris was removed from the plots. Each of the five thinning treatments was replicated into four blocks, totaling 20 sampling plots (five thinning intensities × four blocks). The layout of the experimental design is shown in Figure 1 (including the blocks and plots).

Soil, Litter and Vegetation Survey
The soil sampling was conducted at three soil depths (0-20 cm, 20-40 cm and 40-60 cm) using a soil auger (40-mm diameter) in August 2018. We collected nine replicate soil samples following an "S"-shaped pattern at three depths in each plot ( Figure 1). Then, the collected soil samples were fully homogenized from the same depth to form a composite soil sample. In total, 60 composite soil samples (5 treatments × 4 blocks × 3 depths) were collected. Plant and fauna residues were manually removed, and the soil was then passed through a 2-mm screen. The soil samples were then divided into two portions: the first part was air-dried to measure the soil organic carbon, pH, available nitrogen, available phosphorus and available potassium. The second portion was used to measure the water content after oven-drying at 105 °C for 48 h. Soil samples of three duplicates were collected at depths of 0-20 cm, 20-40 cm and 40-60 cm by volumetric rings (100 cm 3 ) after continuously sunny conditions in order to measure the soil's bulk density [34,40].
The tree height (H) and diameter at breast height (DBH ≥ 5 cm, 1.3 m) in each plot were measured. The understory species diversity was investigated in five shrub subplots (2 × 2 m) and five herb subplots (1 × 1 m) established along the diagonals in each plot. Whole-plant sampling techniques were used to determine the shrub and herb biomass [33]. For the litter sampling, all of the organic material (the undecomposed and decomposed parts on the ground) in five 1 × 1 m subplots were collected. Herb and litter subplots were located on larger shrub subplots, and all of the vegetation surveys were carried out in August at peak vegetation coverage, as previously described [41].

Fine-root Sampling
The sequential soil coring method was used to collect fine-root biomass, production, mortality and turnover rate data using a previously described method [42]. Because the change in fine-root characteristics exhibits strong seasonal variations, fine roots were sampled throughout the year [43]. Furthermore, we expanded the soil depth interval to 0-60 cm based on the average soil layer thickness in the study area.
In each sampling plot, we randomly collected eight soil cores (90-mm inner diameter) over the first three days of September (autumn) and November (winter) 2018, and of April Figure 1. Location of the study area and the layout of the experimental plots. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively.

Soil, Litter and Vegetation Survey
The soil sampling was conducted at three soil depths (0-20 cm, 20-40 cm and 40-60 cm) using a soil auger (40-mm diameter) in August 2018. We collected nine replicate soil samples following an "S"-shaped pattern at three depths in each plot ( Figure 1). Then, the collected soil samples were fully homogenized from the same depth to form a composite soil sample. In total, 60 composite soil samples (5 treatments × 4 blocks × 3 depths) were collected. Plant and fauna residues were manually removed, and the soil was then passed through a 2-mm screen. The soil samples were then divided into two portions: the first part was air-dried to measure the soil organic carbon, pH, available nitrogen, available phosphorus and available potassium. The second portion was used to measure the water content after oven-drying at 105 • C for 48 h. Soil samples of three duplicates were collected at depths of 0-20 cm, 20-40 cm and 40-60 cm by volumetric rings (100 cm 3 ) after continuously sunny conditions in order to measure the soil's bulk density [34,40].
The tree height (H) and diameter at breast height (DBH ≥ 5 cm, 1.3 m) in each plot were measured. The understory species diversity was investigated in five shrub subplots (2 × 2 m) and five herb subplots (1 × 1 m) established along the diagonals in each plot. Whole-plant sampling techniques were used to determine the shrub and herb biomass [33]. For the litter sampling, all of the organic material (the undecomposed and decomposed parts on the ground) in five 1 × 1 m subplots were collected. Herb and litter subplots were located on larger shrub subplots, and all of the vegetation surveys were carried out in August at peak vegetation coverage, as previously described [41].

Fine-Root Sampling
The sequential soil coring method was used to collect fine-root biomass, production, mortality and turnover rate data using a previously described method [42]. Because the change in fine-root characteristics exhibits strong seasonal variations, fine roots were sampled throughout the year [43]. Furthermore, we expanded the soil depth interval to 0-60 cm based on the average soil layer thickness in the study area.
In each sampling plot, we randomly collected eight soil cores (90-mm inner diameter) over the first three days of September (autumn) and November (winter) 2018, and of April (spring) and June (summer) 2019 (no samples were collected from December to March because low temperatures had caused the soil layer to freeze). We collected samples at depths of 0-20, 20-40 and 40-60 cm using a soil auger, producing 480 (5 treatments × 4 blocks × 8 cores × 3 depths) samples seasonally and 1920 (480 samples × 4 seasons) samples for the four seasons of sampling used for the laboratory analysis. The root samples collected at each soil depth were thoroughly mixed to form one fine-root sample. In total, 240 root samples (5 treatments × 4 blocks × 3 depths × 4 seasons) were collected over the four seasons.
The composite fine-root samples were transported to the laboratory in an icebox. In order to separate the roots from the soil, we first soaked the fine-root samples in water. Then, three diameter classes of fine roots (<0.5, 0.5-1, and 1-2 mm, determined using electronic calipers) were carefully washed and sorted into living and dead groups according to their status using the method described by Brassard et al. (2013) [44]. Live roots were classified as having a pale exterior, as elastic and flexible, and as being free of decay, with a whitish cortex, while dead roots were brown or black in color, and were inflexible. Finally, all of the live and dead fine roots of the three diameters were oven-dried at 65 • C to a constant mass.

Chemical and Biochemical Analyses
All of the soil chemical indicators were determined following a previously described method [40]. The soil's organic carbon content of the soils was measured using the K 2 Cr 2 O 7 oxidation method. The soil's available nitrogen was identified by alkaline hydrolysis diffusion, and the available phosphorus was measured by colorimetry after extraction with NaHCO 3 . The soil's available potassium was extracted in ammonium acetate (pH 7.0) and identified on a flame photometer. The soil pH was determined in a 1:2.5 soil:water suspension. The soil's bulk density was obtained by calculating the ratio of the soil's mass to the total volume (g·cm −3 ) after oven drying at 105 • C to a constant weight [45].

Data Calculation and Analysis
The fine-root biomass (g m −2 ) and necromass (g m −2 ) were calculated for each sampling season in each plot by summing the dry weight of the live and dead fine roots in each soil core. The fine-root production (g m −2 year −1 ) and mortality (g m −2 year −1 ) were determined using a simplified decision matrix method (Table 1) [14]. The fine-root turnover rate (year −1 ) was defined as the ratio of annual fine-root production (g m −2 year −1 ) to the mean biomass (g m −2 ) of fine roots over a year [46]. Table 1. Simplified decision matrix for the calculation of the production and mortality of the fine roots. In the present study, the biomass characteristics of fine roots were repeatedly measured across sampling seasons according to the soil depth within each plot. Therefore, we performed a linear mixed-model analysis with three fixed effects (thinning intensity (T), sampling season (S), and soil depth (D)) and random effects (plot and block), as described by where Y ijkl is the fine-root biomass, necromass, or total mass (g m −2 ); T i (i = 0, 15, 30, 45, 60) is the thinning intensity; S j(l) is the sample season (i.e., autumn, winter, spring, or summer); D k(l) is the soil depth (0-20 cm, 20-40 cm, and 40-60 cm); and π l is a random plot or block effect (l = 1, 2, . . . , 20). Under different thinning intensities, the changes in the fine-root characteristics were highly variable at different soil layers. Thus, the fine-root characteristic percentage change (compared to CK) of the average for all of the thinning intensities at a certain soil depth could be seen as an indicator for the evaluation of the fine-root characteristic responses at different soil depths to thinning intervention. A higher percentage change of one soil depth indicates a more sensitive response to thinning intensities, which was calculated as follows: where C p denotes the fine-root biomass, necromass, total mass, production, mortality or turnover rate percentage change value; T n (n = 1, 2, 3, 4) is the thinning intensity; C Tn and C CK are the fine-root characteristic mean values for the thinning and control conditions, respectively; and |C Tn − C CK | represents the absolute value of C Tn − C CK . The effects of thinning intensities on the soil properties (water content, bulk density, soil organic carbon, available nitrogen, available phosphorus, available potassium and pH) and stand characteristics (tree density, height, DBH and volume, understory vegetation biomass and species diversity index, and litter biomass) were also tested using a linear mixed-effects model ANOVA. For all of the models, the significance of fixed effects was assessed using Satterthwaite approximations for degrees of freedom. When the fixed effects or interactions were significant, the least square means differences test was performed for multiple comparisons (main effect or simple effect analysis). The statistical value F was used to evaluate the sensitivity differences of three diameter classes of fine roots to thinning intensities. The linear mixed-effects model was obtained with the 'lmerTest' and 'lme4' packages [47,48]. Multiple comparisons were drawn using the 'emmeans' package. Principal component analysis (PCA) was performed to determine the relationships between the fine-root characteristics, stand characteristics and soil properties, using the 'FactoMineR' package [49]. All of the analyses were implemented using R for Windows statistical software, version 4.1.1 [50].

Stand and Soil Properties
The stand and soil properties changed in different ways following the thinning treatments ( Table 2). The stand volume and litter biomass significantly decreased with increasing thinning intensities (p < 0.05), whereas the biomass and Shannon-Wiener index of the shrubs and herbs showed the opposite trend (p < 0.05).
The soil's organic carbon, available nitrogen, available phosphorus, available potassium, water content and pH values decreased with increasing soil depths, while the soil's bulk density increased. T2, T3 and T4 had significantly lower soil organic carbon, available nitrogen, available phosphorus and available potassium values than their CK values, decreasing by 47.1-72.3%, 23.2-77.6%, 18.7-61.8 and 6.3-53.6% (p < 0.05), respectively. Thinning had no significant effects on the soil's water content, bulk density or pH (p > 0.05) (except for the pH at 40-60 cm) ( Table 2).

Fine-Root Biomass, Necromass, and Total Mass of Different Diameter Classes
The whole fine-root system (Ø ≤ 2 mm) biomass (F = 2.11, p > 0.05), necromass and total mass (F = 2.33, p > 0.05) values of the total soil depth (0-60 cm), on average, did not differ among the treatments, except for the significant increase in necromass occurring in T4 (F = 25.31, p < 0.01) (Figure 2d,h,l and Table 3). The very fine-root (<0.5 mm) biomass and total mass fluctuated with the increasing thinning intensity. The very-fine root biomass (F = 10.91, p < 0.01) and total mass (F = 10.62, p < 0.01) in T1 and T3 were significantly lower than those in CK, while those in T2 and T4 exceeded those in CK (although not significantly) (Figure 2d,l and Table 3). The 0.5-1 mm (F = 4.42, p < 0.05) and 1-2 mm (F = 9.01, p < 0.01) fine-root biomass, and the 1-2 mm (F = 4.15, p < 0.05) total mass decreased with the thinning intensity, showing a significant difference at T3 and T4 (Figure 2d,l and Table 3). The 0.5-1 mm and 1-2 mm fine-root necromass mirrored the overall fine-root Biology 2022, 11, 351 7 of 18 system necromass characteristics, whereas very fine-root necromass was rarely observed (Figure 2h). The fine-root biomass, necromass, and total mass decreased with the soil depth, and the values in the topsoil (0-20 cm) accounted for more than 70% of the total values for all of the soil depths ( Figure 2). The biomass and the total mass of the very fine roots were mainly distributed in the topsoil, and accounted for approximately 50% of the values in the topsoil and 40% of the total values. The biomass and total mass of 0.5-1 mm and 1-2 mm roots dominated the deep soil depths, while the necromass of 0.5-1 mm and 1-2 mm roots was dominant at all of the soil depths ( Figure 2). The very fine-root biomass and total mass were more sensitive to the thinning intensities than those of other size classes across all of the soil depths (Table 3 and Figure 2). The fine-root biomass, necromass, and total mass at the middle and deep soil depths exhibited greater percentage changes than those in the topsoil following thinning, and these characteristics in the deep soil increased again in T3 and T2 (Figure 2 and Table 4).

Figure 2.
The effects of thinning intensities on the fine-root biomass, necromass, and total mass (average over four seasons) among different diameter classes at different soil depths (a-c,e-g,i-k), and the entire soil profile (d,h,l). The values are the mean of four replicates ± SE. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively. Different lowercase letters (a, b, c, d) indicate significant differences for the same diameter class fine roots among the thinning intensities (p < 0.05). Table 3. The effects of the thinning intensity (T) on the mean annual biomass, necromass, and total biomass of fine roots over the entire soil profile, and at specific depths (D).  . The values are the mean of four replicates ± SE. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively. Different lowercase letters (a, b, c, d) indicate significant differences for the same diameter class fine roots among the thinning intensities (p < 0.05). Table 3. The effects of the thinning intensity (T) on the mean annual biomass, necromass, and total biomass of fine roots over the entire soil profile, and at specific depths (D).  Note-NA: no data were detected in the control group.

Fine-root Production, Mortality and Turnover rate in Different Diameter Classes
The whole fine-root system (Ø ≤ 2 mm) production and mortality of the entire soil profile (0-60 cm) and the averaged turnover rate did not differ among the different thinning intensities (all p > 0.05) (Figure 4d,h,l and Table 5). The very fine-root production also fluctuated following thinning, and significantly exceeded that of the CK in T2 and T4 (F = Figure 3. The effects of the season on the fine-root biomass (live) among different diameter classes in different soil depths, and the entire soil profile. The values are the mean of four replicates ± SE. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively.

Fine-Root Production, Mortality and Turnover Rate in Different Diameter Classes
The whole fine-root system (Ø ≤ 2 mm) production and mortality of the entire soil profile (0-60 cm) and the averaged turnover rate did not differ among the different thinning intensities (all p > 0.05) (Figure 4d,h,l and Table 5). The very fine-root production also fluctuated following thinning, and significantly exceeded that of the CK in T2 and T4 (F = 5.7, p < 0.01) (Figure 4d and Table 5). The very fine-root turnover rate in T3 was significantly lower than those at the other thinning intensities (F = 2.71, p < 0.05) (Figure 4l and Table 5). The thinning intensities did not affect the production and turnover rates of 0.5-1 mm and 1-2 mm fine roots (all p > 0.05) (Figure 4d,l and Table 5). The 0.5-1 mm fine-root mortality value in T1 was higher than that under the other treatments (F = 5.27, p < 0.01). The 0.5-1 mm and 1-2 mm fine-root mortality levels mirrored the total mortality levels ( Figure 4h and Table 5).  . The effects of the thinning intensities on the fine-root production, mortality and turnover rate among different diameter classes at different soil depths (a-c,e-g,i-k), and the entire soil profile (d,h,l). The values are the mean of four replicates ± SE. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively. Different lowercase letters (a, b, c, d) indicate significant differences for the same diameter class of fine roots among the thinning intensities (p < 0.05). Table 5. The effects of the thinning intensity (T) on the annual production, mortality, and turnover rate of fine roots over the entire soil profile, and at specific depths (D).  Figure 4. The effects of the thinning intensities on the fine-root production, mortality and turnover rate among different diameter classes at different soil depths (a-c,e-g,i-k), and the entire soil profile (d,h,l). The values are the mean of four replicates ± SE. CK, T1, T2, T3, and T4 represent 0%, 15%, 30%, 45% and 60% thinning intensities, respectively. Different lowercase letters (a, b, c, d) indicate significant differences for the same diameter class of fine roots among the thinning intensities (p < 0.05). The production and mortality of fine roots decreased with the soil depth, and the values in the topsoil accounted for approximately 69% of the total production and 71% of the total mortality (Figure 4). The production of very fine roots also occupied the topsoil, and accounted for approximately 51% of the topsoil and 35% of the total value. The mortality of 1-2 mm roots was greater at all of the soil depths (Figure 4e-g). The fine-root turnover rate generally increased with deepening soil depths (Figure 4i-k). The very fine-root production and turnover rate were more sensitive to thinning intensities than those of other size classes among the soil depths (Table 5 and Figure 4). The percentage changes in the fine-root production, mortality and turnover rates in the deep soil layers were also generally higher than those in the topsoil following thinning, and these characteristics in the deep soil were enhanced again in T2 (Figure 4 and Table 4).

The Linkages between Fine-Root Characteristics and Stand and Soil Attributes
The PCA of the entire soil profile according to the standardized data shows that the first two trait axes accounted for 32.3% and 14.3% of the total variation, respectively ( Figure 5). We found that the 0.5-1 mm and 1-2 mm fine-root biomasses were highly positively related to the decreased soil properties (soil organic carbon, available nitrogen, available phosphorus and available potassium) and the decreased stand characteristics (volume and litter biomass) but were negatively correlated with the increased Shannon-Wiener index and biomass of shrubs and herbs ( Figure 5 and Table S4). Conversely, the biomass, production, and turnover rate of very fine roots and the necromass of 0.5-1 mm and 1-2 mm fine roots were negatively correlated with the decreased soil nutrient availability, stand volume and litter biomass, whereas they exhibited strong positive correlations with the increased Shannon-Wiener index and the biomass of shrubs and herbs ( Figure 5 and Table S4). The PCA of individual soil depths shows that the association between the understory vegetation characteristics and the very fine-root portion was stronger in the topsoil, whereas this correlation disappeared with the thicker fine-root portion appearing at deeper soil depths ( Figure S3).  Tables S2 and S3). AN, AP, AK, SOC, WC and BD: the available nitrogen, available phosphorus, available potassium, soil organic carbon, water content, and soil bulk density; DBH, H, SD, SS, HS, SB, HB, LB: diameter at breast height, tree height, stem density, shrub Shannon-Wiener index, herb Shannon-Wiener index, shrub biomass, herb biomass, and litter biomass; P<0.5, P0.5-1, and P1-2: <0.5, 0.5-1, and 1-2 mm fine-root production; M<0.5, M0.5-1, and M1-2: <0.5, 0.5-1, and 1-2 mm fine-root mortality; T<0.5, T0.5-1, and T1-2: <0.5, 0.5-1, and 1-2 mm fine-root turnover rate; B<0.5, B0.5-1, and B1-2: <0.5, 0.5-1, and 1-2 mm fine-root biomass; N<0.5, N0.5-1, and N1-2: <0.5, 0.5-1, and 1-2 mm fine-root necromass.

Effects of Thinning on the Fine-Root Biomass, Production and Turnover Rate
As is consistent with our hypothesis, we found that the very fine-root biomass, production and turnover rates were more sensitive to thinning intensities than thicker fine roots. The results align with a recent observation that forest cutting has more significant effects on fine roots than on thicker roots [18]. These results may be attributed to the fact that fine roots of different diameter classes exhibit heterogeneous physiological functions and structural compositions, leading to discrepant responses following thinning intensities [21]. Moreover, a similar study found that very fine roots with a higher ratio of nonstructural to structural mass are more sensitive to changes in abiotic or biotic factors caused by thinning treatments [24].  Tables S2 and S3). AN, AP, AK, SOC, WC and BD: the available nitrogen, available phosphorus, available potassium, soil organic carbon, water content, and soil bulk density; DBH, H, SD, SS, HS, SB, HB, LB: diameter at breast height, tree height, stem density, shrub Shannon-Wiener index, herb Shannon-Wiener index, shrub biomass, herb biomass, and litter biomass; P <0.5 , P 0.5-1 , and P 1-2 : <0.5, 0.5-1, and 1-2 mm fine-root production; M <0.5 , M 0.5-1 , and M 1-2 : <0.5, 0.5-1, and 1-2 mm fine-root mortality; T <0.5 , T 0.5-1 , and T 1-2 : <0.5, 0.5-1, and 1-2 mm fine-root turnover rate; B <0.5 , B 0.5-1 , and B 1-2 : <0.5, 0.5-1, and 1-2 mm fine-root biomass; N <0.5 , N 0.5-1 , and N 1-2 : <0.5, 0.5-1, and 1-2 mm fine-root necromass.

Effects of Thinning on the Fine-Root Biomass, Production and Turnover Rate
As is consistent with our hypothesis, we found that the very fine-root biomass, production and turnover rates were more sensitive to thinning intensities than thicker fine roots. The results align with a recent observation that forest cutting has more significant effects on fine roots than on thicker roots [18]. These results may be attributed to the fact that fine roots of different diameter classes exhibit heterogeneous physiological functions and structural compositions, leading to discrepant responses following thinning intensities [21]. Moreover, a similar study found that very fine roots with a higher ratio of nonstructural to structural mass are more sensitive to changes in abiotic or biotic factors caused by thinning treatments [24].
The very fine-root biomass, production and turnover rate levels fluctuated (positive or negative effects) with increasing thinning intensities, which is inconsistent with the positive or negative results obtained following thinning in previous studies [3,32]. Regarding positive effects, our PCA results suggest that the increase in the biomass and productivity of very fine roots resulted from increases in herb and shrub layer species diversity and biomass (especially in the topsoil) ( Figure 5 and Figure S3a), which compensated for decreased fine-root biomass and production resulting from the cutting of canopy trees, as is consistent with early studies [14,31]. In addition, previous research has demonstrated that low-nutrient conditions can stimulate the growth of fine-root biomass and productivity [51]. In the present study, our whole-tree harvesting measures (reduced stand volume and liter biomass) increased the export of nutrients and reduced soil nutrient availability (Table 2), as is consistent with an early study [52]. Furthermore, the PCA results show that very fine-root biomass and production were negatively correlated with soil nutrient availability, stand volume and litter biomass, supporting this view. Regarding negative effects, a possible explanation is that thinning reduces resource competition pressure, and low fine-root biomass, production and turnover rate values may satisfy vegetation resource absorption and utilization requirements [53].
In contrast, we found that higher thinning intensities significantly reduced the biomass of 0.5-1 mm and 1-2 mm fine roots, supporting the negative effect of thinning found in previous studies [54]. On the one hand, our PCA results indicate that thinning reduced the stand volume and litter biomass, which may decrease the photosynthate partitioning to the root system, and may ultimately lead to the reduced biomass and productivity of 0.5-1 mm and 1-2 mm fine roots [18]. Early studies reported that the nutrient acquisition strategy for thicker fine roots may be achieved by increasing their lifespan and extending the period of nutrient absorption, which is strongly dependent on nutrient availability [20,55,56]. Our PCA results show that 0.5-1 mm and 1-2 mm fine-root biomass was positively correlated with soil nutrient availability, and negatively correlated with herb and shrub species diversity and biomass. Therefore, the decreased stand volume, litter biomass and nutrient availability following thinning together reduced the 0.5-1 mm and 1-2 mm fine-root biomass. Overall, our findings provide a reasonable explanation for the inconsistent impacts of thinning on fine-root biomass and production observed in previous studies [3,32].
In terms of fine-root biomass and productivity, our results reveal that a 30% approach (T2: moderately reduced aboveground biomass) may be a more suitable thinning strategy to promote productivity and increase stand heterogeneity in forest ecosystems. Previous studies have also reported that root gaps in forests recover faster and are more ephemeral than canopy gaps after thinning [57,58], demonstrating that fine-root characteristics serve as very valuable data for the selection of a suitable thinning intensity. Moreover, our findings align with traditional studies showing that a moderate thinning intensity could better improve stand productivity and increase the diversity of understory plants compared to other thinning intensities [59,60].

Necromass and Mortality Changes following Thinning
In this study, the thinning intensity had no influence on the fine-root necromass of the entire soil profile (except for T4). Our observed characteristics of the fine-root necromass are inconsistent with the results of Wang et al. (2019) and Ma et al. (2013) [14], who found thinning intensities to decrease or increase necromass. On the one hand, previous studies have demonstrated that the biomass and necromass of very fine roots account for a greater proportion of the fine-root system than other root size classes [24], and our very fine-root biomass results support this conclusion, which may theoretically determine the necromass. However, very fine roots also possess higher nonstructural carbohydrate concentrations and decompose more easily; thus, these roots are rarely observed [44,61,62]. On the other hand, higher thinning intensities may improve the soil conditions (e.g., soil temperature, Figure S4), and may also further increase the fine-root decomposition rate [32]. Thus, the fine-root necromass levels did not differ between the treatments, potentially due to a large amount of fine-root necromass decomposing and disappearing [63,64]. We note, however, that T4 exhibited significantly increased necromass levels. This finding may be attributable to the fact that high-intensity thinning and species replacement produced a large amount of fine-root necromass [28,30], especially for 0.5-1 mm and 1-2 mm fine roots, which had not completely decomposed, resulting in a significant increase in necromass in T4.

Response of Deeper-Soil Fine-Root Characteristics to Thinning
As hypothesized, we found the thinning effects on the change in the fine-root characteristics to be stronger at deep soil depths, and the fine-root characteristics exhibited fluctuating patterns. These findings are consistent with the results of a recent study showing dramatically altered fine-root biomass and necromass levels at deeper soil depths following thinning [18]. This response may reflect the resource acquisition strategy of the root system following thinning. In the topsoil layer, the fine-root systems of shallow root understory plants and the remaining trees quickly colonize and recover in the area liberated by the disturbance; thus, these roots are less affected by thinning intensities [17,65]. It is well known that trees and understories jointly determine the characteristics of roots in the surface soil layer due to belowground niche partitioning, while trees determine such characteristics at deep soil depths ( Figure S3) [66,67]. In our study, all of the four thinning intensities reduced the fine-root biomass and productivity at deep soil depths (Figures 2 and 4), which may be because these intensities reduced the fine-root densities of trees, and consequently alleviated resource competition pressure at this depth [68]. Compared to low-and high-intensity thinning at deep soil depths, higher fine-root biomass and productivity levels were observed in the T2 or T3 thinning treatments (Figures 2 and 4). This change could be attributed not only to thinning interventions but also to the regeneration of understory plants [10]. In the T2 or T3 thinning treatments, the faster regeneration of understory plants excessively consumes resources, which results in greater resource competition pressure in the surface soil, and requires trees to adjust their rooting depth and increase their fine-root growth at deeper soil depths with less root competition [14]. Thus, we observed the fluctuating phenomena of the fine-root characteristics in the deep soil layers. However, the fine-root characteristics at deep soil depths could not recover in a short period of time after thinning (Figures 2 and 4), reflecting their more sensitive response to thinning intensities. Moreover, 0.5-1 and 1-2 mm fine roots dominated the fine-root characteristics at the deep soil depths, and exhibited an increased turnover rate, indicating that thicker fine roots could better mirror potential carbon pools at deeper soil depths in forest ecosystems.
Some studies have reported that fine roots can be categorized according to their functions (e.g., the first order of the root system plays a role in absorption) [44,69]. Furthermore, fine-root characteristics may vary among different functional plant groups [14]. Although the current study uses a more nuanced classification method that diverges from the traditional definition (≤2 mm) to study the stand-level characteristics of fine roots, future studies could build on root functional approaches and plant functional group distinctions in order to better understand how root function specificity and species diversity impact belowground processes at the ecosystem level.

Conclusions
Our study suggests that different thinning intensities have substantial effects on the changed characteristics of fine roots with different diameters at different soil depths. The positive effect of thinning on very fine roots and the negative effect on thicker fine roots and all of the diameter classes of fine roots provide reasonable explanations for the inconsistent effects of thinning intensities. Here, a 30% (T2) thinning intensity moderately reduced the aboveground biomass and yielded increased biomass and productivity among very fine roots compared to the other treatments, suggesting that a 30% approach is a more suitable thinning strategy to promote productivity and increase stand heterogeneity in forest ecosystems. The fine-root characteristics at deeper soil depths are more sensitive to thinning intensities. The sampled 0.5-1 and 1-2 mm fine roots dominated the fine-root characteristics at deep soil layers and exhibited a higher turnover rate, indicating that thicker fine roots could better mirror the potential carbon pools of deeper soils in forest ecosystems. Collectively, our findings provide important insights into the effects of forest management on change in fine-root characteristics, and supplement meaningful data on fine-root productivity to improve the parameterization of future ecological models.