Radial Variation of Wood Anatomical Properties Determines the Demarcation of Juvenile-Mature Wood in Schima superba

: Schima superba is a precious broadleaf tree species that produces excellent timber. Analyzing the radial variation of wood anatomical properties within stems and determining the boundary between the juvenile and mature wood can provide a theoretical basis for the breeding and cultivation of S . superba . Herein, 30 54-year-old trees were used to analyze the radial variation of wood anatomical properties qualitatively. The correlation among wood properties was analyzed. On this basis, four types of model were used to predict the wood properties quantitatively, and the optimal segmentation model was used to determine the boundary between juvenile and mature wood. The results showed that the variation of the early stage (before the 20th year) was larger than that of the later stage (after the 20th year); moreover, the variation of the vessel features (9.56%) was greater than that of the ﬁber features (7.42%), the vessel lumen diameter (11.94%) and the ﬁber length (8.00%) had the most variation. There were three radial variation patterns in the wood properties. In the early stage, there was a signiﬁcant positive correlation between the ring width and wood air-dry density (0.78) and wood basic density (0.89), and a signiﬁcant negative correlation between the ring width and the ﬁber characters ( − 0.79~ − 0.93) and vessel characters ( − 0.64~ − 0.90). The growth models of Logistic and Richards were more suitable ( R 2 = 84.2%~95.2%) for the radial variation of wood anatomical properties than the nonlinear least square model ( R 2 = 67.7%~90.9%). The transition wood between juvenile and mature wood of S. superba was from the 11th to the 16th year based on aggregative indicators from the pith to bark at breast height. Finally, when taken together, the wood anatomical properties of S. superba had regular radial variation. Breeding programs need to consider both growth and wood anatomical properties in the early stage, but they can be neglected in the mature stage. The formation of mature wood of S. superba occurred after the 16th year.


Instruction
One main function of a forest is to produce wood. The wood properties directly determine wood utilization. Wood density is an important property due to its correlation with other wood properties, such as stiffness and ability to resist deformation [1,2]. It is a composite trait affected by the annual ring width, latewood percentage, cell wall thickness and lumen diameter of the fibers and vessels, etc. [3]. Fibers and vessels are the main components of the wood anatomical structure in broadleaf trees [4], and their size, quantity and shape are important indicators affecting wood properties [5]. Wood properties are affected not only by genes [6] and the environment [7,8], but also by cambial age [9]. The radial variation of wood properties within plants caused by the xylem formed at different cambium ages is the main source of wood property variation. This variation shows a certain pattern with cambial age increasing in different tree species. Studying this variation pattern is beneficial for the wood breeding programs and utilization of this tree species.
The radial variation of wood anatomical properties can be not only qualitatively analyzed but also quantified. The theoretical growth equation is used as a model to describe

Study Site
The sampling site is located in Jian'ou, Fujian province, China (117 • 58 -118 • 57 E, 26 • 38 -27 • 20 N). The altitude is 140~180 m. The climate is characterized by distinct midsubtropical maritime monsoon influence. The mean annual temperature is 18.9 • C with average temperatures of 8 • C in January and of 28.6 • C in July. The annual precipitation is 1627.3~1662.9 mm, the relative humidity is 81% and the frost-free period is 286 days.

Plant Material
Five pure S. superba plantations at different ages were selected ( Table 1). The afforestation land was all low hills, and the soil was mountainous red soil with above medium fertility. The soil physical and chemical properties in each sample plot were tested (Table 2) so that the growth conditions of the study samples were as similar as possible.  25.80 In each plantation of a different age, survey plots of 20 × 20 m were set up. The diameter at breast height (DBH) of all trees was measured. Here, 30 sample trees were randomly selected. The sample trees were selected based on the mean DBH and the straightness of stems showing no apparent defects, and the height and DBH were measured (Table 1). At the same time, 150 intact wood cores from pith to bark for wood anatomical properties analyses were taken from the upper slope side of each trunk, and the wood cores were subsequently taken back to the laboratory for disposal. The wood cores were polished and processed into wood strips with a width of 5 mm and a thickness of 2.5 mm.

Measurements of Wood Anatomical Properties
The observation of the annual ring width (RW) and latewood percentage (LP) was made under a microscope (Zeiss Imager. A2, Germany) ( Figure 1). Every 2 annual rings of the wood strips were divided into 1 pith-to-bark oriented sample. Wood air-dry density (WADD) was measured by a soft X-ray microdensitometer (type 2405). The wood basic density (WBD) was measured by the maximum moisture content method (wet weight and dry weight of the samples were measured, respectively) [23]. The fiber and vessel morphology were determined using the mixture segregation of 30% hydrogen peroxide and glacial acetic acid (1:1) at 90 • C for 4 h. The macerated tissue was washed in running water until the tissue was free of all acid traces, so that the lengths of the fiber and vessel elements could be determined. The macerated tissue was then crushed in water to make a temporary section by safranine (1%) staining, and 50 intact fibers and vessels were randomly measured, respectively, under a biological projection microscope (Zeiss Imager. A2, Germany). Software ImageJ (v1.52) was used to analyze all images.
where n is the number of the sample, and k is the category of sample. The smaller the value of L[ ( , )] is, the more reasonable the classification is. The descriptive statistical analysis of wood properties was calculated using the MEAS procedure in SAS (v9.4). The model parameters were estimated using PROC NLIN and PROC NLMIXED in SAS (v9.4).

Radial Variation among Different Tree Ages in Wood Anatomical Properties
The radial variation patterns of different tree ages in wood anatomical properties were similar (Figures 2-8). Only the values of wood air-dry density and wood basic density with different ages were quite different. The wood air-dry density and basic density of 54-year-old specimens were larger than those of other ages near the pith. The wood anatomical properties in 4-year-old, 6-year-old, 10-year-old and 20-year-old specimens had not yet reached a steady state. Therefore, 54-year-old S. superba specimens were selected for analysis.

Description of the Wood Anatomical Properties
The annual ring width of 54-year-old S. superba fluctuated greatly ( Table 3). The maximum ring (5.94 mm yr −1 ) width was 3.83 times higher than the minimum ring width (1.23 mm yr −1 ). The latewood (54.71% ± 6.71%) accounted for a larger proportion than earlywood. The maxima of wood air-dry density and wood basic density were 23.08% and 24.00% greater than their minima, respectively. The wood anatomical properties showed longer fibers (1838.58 ± 147.07 μm) and a smaller fiber width (32.63 ± 2.28 μm), compared to vessel length (1297.34 ± 107.16 μm) and vessel width (79.02 ± 8.38 μm). The fibers had a thicker cell wall thickness (16.33 ± 1.28 μm). This value was close to the fiber lumen diameter (16.31 ± 1.12 μm). Wood anatomical properties including fiber length (FL), fiber width (FW), fiber lumen diameter (FLD), fiber cell wall thickness (FCWT), vessel length (VL), vessel width (VW), vessel lumen diameter (VLD) and vessel cell wall thickness (VCWT) were measured. Fiber and vessel wall-lumen ratio was calculated using the double cell wall thickness divided by the lumen diameter. Fiber and vessel lumen-width ratio was calculated as lumen diameter divided by width. Fiber and vessel length-width ratio was calculated as length divided by width. In this study, wood properties included wood air-dry density, wood basic density and wood anatomical properties.

Data Analyses
Wood basic density was calculated according to the following formula: where M is the weight (g) of the sample at saturated moisture content, M 0 is the weight (g) of the oven-dry sample and D w is the specific gravity of the wood material forming the cell wall with an average value of 1.53. Coefficient of variation: where σ is the standard deviation of the wood property, and X is the mean of the wood property. Correlation coefficient: where Cov (x, y) is the covariance of the wood properties x and y, and σ 2 x and σ 2 y are the variance of the wood properties x and y, respectively.
The wood properties of 54-year-old plants were used to model according to the following equations: Binomial regression [24]: Logistic [10]: Richards [1]: where CA ij , FW ij and FWT ij are the cambial age, fiber wall thickness, and vessel lumen diameter of the jth annual ring from the ith tree, respectively. β 1 , β 2 , β 3 and β 4 are the model parameters.
The determination coefficient (R 2 ), root mean square error (RMSE), and mean absolute deviation (MAD) were used to evaluate and test the model [10]. The larger the coefficient of determination was, the smaller the root mean square error and mean absolute error was, which indicates that the accuracy of the model simulation is higher.
where y i is the observed value, andŷ i is the predicted value, and y i is the mean value of the sample. Optimal segmentation model: Assume that G includes sample X i , X i+1 , ······X j (j>i), where n is the number of the sample, and k is the category of sample. The smaller the value of L[b(n, k)] is, the more reasonable the classification is. The descriptive statistical analysis of wood properties was calculated using the MEAS procedure in SAS (v9.4). The model parameters were estimated using PROC NLIN and PROC NLMIXED in SAS (v9.4).

Radial Variation among Different Tree Ages in Wood Anatomical Properties
The radial variation patterns of different tree ages in wood anatomical properties were similar (Figures 2-8). Only the values of wood air-dry density and wood basic density with different ages were quite different. The wood air-dry density and basic density of 54-year-old specimens were larger than those of other ages near the pith. The wood anatomical properties in 4-year-old, 6-year-old, 10-year-old and 20-year-old specimens had not yet reached a steady state. Therefore, 54-year-old S. superba specimens were selected for analysis. vessel properties, the largest variation was in the vessel lumen diameter (11.94%), and the smallest variation was in the vessel cell wall thickness (7.43%).
Overall, the annual ring width rapidly decreased from the pith and slightly decreased close to the bark ( Figure 2). This reflected that S. superba grew faster and had a larger diameter increment in the early stage. The wood density decreased rapidly in the first 20 years with subsequent stabilization ( Figure 3A,B). In the first 20 years, the mean wood airdry density and basic density were 0.75 g·cm −3 and 0.56 g·cm −3 , respectively. After 20 years, they decreased to a mean of 0.70 g·cm −3 and 0.52 g·cm −3 , respectively.    The value of the fiber and vessel length increased gradually from the pith and then decreased after reaching the maximum value ( Figures 4A and 5A). The mean value of fiber length was 1769.55 ± 207.22 μm from the 2nd year to the 20th year, reaching the maximum at the 20th year (2034.11 ± 107.86μm) with a growth increase of 35.97%, and the  Vessel lumen diameter and vessel cell wall thickness increased from the pith and approached an asymptotic value near the outer wood, but there was a larger fluctuation in the later stage. The vessel lumen diameter reached the maximum (73.13 μm) in the 24th year, with growth increasing to 56.58%, and the variation was 14.47%. It then fluctuated above and below 69.05 μm after the 24th year, and the variation decreased to 3.98% (Figure 5C). The vessel cell wall thickness increased from 11.63 to 15.75 μm, with growth increasing to 35.48%, from the 2nd year to the 20th year, and the variation was 8.27%. It remained at a stable value of 14.65 μm, with variation decreasing to 5.95% ( Figure 5D).       There were inter-annual fluctuations in the latewood percentage, fiber and vessel wall-lumen ratio, lumen-width ratio and length-width ratio ( Figures 6-8), but their variation trends with cambial age were not significantly distinct.  There were inter-annual fluctuations in the latewood percentage, fiber and vessel wall-lumen ratio, lumen-width ratio and length-width ratio ( Figures 6-8), but their variation trends with cambial age were not significantly distinct.  5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 51

Correlation of Wood Anatomical Properties
According to the analysis of the radial variation trend of wood properties, the correlation analysis was divided into two sections (Table 4). In the first 20 years, there was an extremely significant positive correlation between cambial age and fiber (0.95~0.99) and vessel (0.90~0.98) morphology. The correlation between cambial age and fiber width was the greatest (R = 0.99). The ring width had a significant positive correlation with wood airdry density (0.78) and wood basic density (0.89). There was an extremely significant positive correlation between fiber and vessel. The fiber indicators were positively correlated with each other, as were the vessel indicators. There was an extremely significant negative correlation between cambial age and ring width (−0.88), wood air-dry density (−0.84) and

Radial Variation in Wood Anatomical Properties
According to the wood property variation trends, we first chose the 20th year as the boundary of early and later stages for analysis ( Figure 9). The results showed that the annual ring width had the greatest total variation (46.29%) than all wood properties (8.52%). The mean value of the radial variation of the wood properties in the early stage (10.18%) was greater than that in the later stage (4.19%). The total variation of the wood air-dry density (5.74%) was similar to that of the wood basic density (5.29%) with minor variation. The total variation of the fiber properties (7.42%) was smaller than that (9.56%) of the vessel properties. For the fiber properties, the largest variation was in the fiber length (8.00%), and the smallest variation was in the fiber lumen diameter (6.84%). For the vessel properties, the largest variation was in the vessel lumen diameter (11.94%), and the smallest variation was in the vessel cell wall thickness (7.43%).
(C) Figure 8. The radial variation of vessel wall-lumen ratio (A), vessel lumen-width ratio (B) and vessel length-width ratio (C) versus cambial age and the mean of differently aged S. superba.

Correlation of Wood Anatomical Properties
According to the analysis of the radial variation trend of wood properties, the correlation analysis was divided into two sections (Table 4). In the first 20 years, there was an extremely significant positive correlation between cambial age and fiber (0.95~0.99) and vessel (0.90~0.98) morphology. The correlation between cambial age and fiber width was the greatest (R = 0.99). The ring width had a significant positive correlation with wood airdry density (0.78) and wood basic density (0.89). There was an extremely significant positive correlation between fiber and vessel. The fiber indicators were positively correlated with each other, as were the vessel indicators. There was an extremely significant negative correlation between cambial age and ring width (−0.88), wood air-dry density (−0.84) and  Overall, the annual ring width rapidly decreased from the pith and slightly decreased close to the bark (Figure 2). This reflected that S. superba grew faster and had a larger diameter increment in the early stage. The wood density decreased rapidly in the first 20 years with subsequent stabilization ( Figure 3A,B). In the first 20 years, the mean wood air-dry density and basic density were 0.75 g·cm −3 and 0.56 g·cm −3 , respectively. After 20 years, they decreased to a mean of 0.70 g·cm −3 and 0.52 g·cm −3 , respectively.
The value of the fiber and vessel length increased gradually from the pith and then decreased after reaching the maximum value ( Figures 4A and 5A). The mean value of fiber length was 1769.55 ± 207.22 µm from the 2nd year to the 20th year, reaching the maximum at the 20th year (2034.11 ± 107.86µm) with a growth increase of 35.97%, and the variation was 11.71%. After the 20th year, the value of the fiber length fluctuated with a downward trend with cambial age, and the variation decreased to 3.80%. The vessel length increased rapidly from the 2nd year to the 14th year with growth increasing to 31.35%. The variation was 10.64% and remained stable from the 14th year to the 26th year with variation decreasing to 1.69%, and vessel length then decreased slightly with variation increasing to 3.33%.
The radial variation patterns of the fiber and vessel width, fiber lumen diameter and fiber cell wall thickness were slightly different from those of the fiber and vessel length. They increased rapidly with cambial age from the pith and remained stable with minor fluctuation. The fiber width varied between 27.11 and 34.28 µm, with growth increasing to 26.45% in the first 20 years and varying around a mean value of 33.87 µm in the later stage ( Figure 4B). The value of vessel width increased from 58.72 to 87.66 µm gradually from the 2nd year to 24th year, with growth increasing to 49.28% and then approached an asymptotic value (83.68 µm) after the 24th year ( Figure 5B). The fiber lumen diameter increased to 29.18%, and its variation was 8.29% from the 2nd year to the 20th year and then remained at a stable value (16.85 µm), with variation decreasing to 3.05% ( Figure 4C). The fiber cell wall thickness reached the maximum (17.42 µm) in the 22nd year, with growth increasing to 25.26% and then fluctuated above and below 16.99 µm after the 22nd year ( Figure 4D).
Vessel lumen diameter and vessel cell wall thickness increased from the pith and approached an asymptotic value near the outer wood, but there was a larger fluctuation in the later stage. The vessel lumen diameter reached the maximum (73.13 µm) in the 24th year, with growth increasing to 56.58%, and the variation was 14.47%. It then fluctuated above and below 69.05 µm after the 24th year, and the variation decreased to 3.98% ( Figure 5C). The vessel cell wall thickness increased from 11.63 to 15.75 µm, with growth increasing to 35.48%, from the 2nd year to the 20th year, and the variation was 8.27%. It remained at a stable value of 14.65 µm, with variation decreasing to 5.95% ( Figure 5D).
There were inter-annual fluctuations in the latewood percentage, fiber and vessel walllumen ratio, lumen-width ratio and length-width ratio (Figures 6-8), but their variation trends with cambial age were not significantly distinct.

Correlation of Wood Anatomical Properties
According to the analysis of the radial variation trend of wood properties, the correlation analysis was divided into two sections (Table 4). In the first 20 years, there was an extremely significant positive correlation between cambial age and fiber (0.95~0.99) and vessel (0.90~0.98) morphology. The correlation between cambial age and fiber width was the greatest (R = 0.99). The ring width had a significant positive correlation with wood air-dry density (0.78) and wood basic density (0.89). There was an extremely significant positive correlation between fiber and vessel. The fiber indicators were positively correlated with each other, as were the vessel indicators. There was an extremely significant negative correlation between cambial age and ring width (−0.88), wood air-dry density (−0.84) and wood basic density (−0.98). The wood air-dry density and basic density had a significant negative correlation with fiber and vessel, especially the correlations between wood air-dry density, wood basic density, fiber cell wall thickness and vessel lumen diameter were the greatest. There was an extremely significant negative correlation between ring width and fiber (−0.79~−0.93) and vessel (−0.64~−0.90).
After 20 years, only the ring width (−0.86), fiber length (−0.85) and vessel length (−0.70) were significantly negatively correlated with cambial age (Table 4). This showed that ring width, fiber length and vessel length continued to decrease with age, while other indicators such as wood air-dry density, wood basic density and width, lumen diameter, cell wall thickness of fiber and vessel were in a stable state.

Models for Predicting Wood Anatomical Properties
The wood air-dry density and basic density are complex wood characteristics determining wood quality. Based on the above correlation results, the fiber width, fiber cell wall thickness and vessel lumen diameter had the greatest correlation with cambial age, wood air-dry density, wood basic density, respectively. Therefore, wood air-dry density, wood basic density, fiber width, fiber cell wall thickness and vessel lumen diameter were selected for model fitting, using fiber width, fiber cell wall thickness and vessel lumen diameter as variables and using multiple linear regression for wood basic density and wood air-dry density (Table 5). Using the Richards model to simulate negative growth was difficult, so this model of wood basic density and air-dry density is not represented in the table. The fitting accuracy of growth model was higher than that of the nonlinear least square model. The wood basic density, wood air-dry density and fiber cell wall thickness had a higher fitting accuracy in the Logistic model, which improved the R 2 to 82.4%, 83.2% and 90.2% and reduced the RESE to 0.012, 0.017 and 0.400, respectively. The Richards model in the fiber width and vessel lumen diameter had a higher fitting accuracy, which accounted for 95.2% and 91.2% of the variation. This shows that the Logistic and Richards model had good predictability and practicality for wood basic density, wood air-dry density, fiber cell wall thickness, fiber width and vessel lumen diameter. In the multiple linear regression of wood basic density and wood air-dry density, the parameter accounting for a larger proportion was fiber cell wall thickness (40.76%) and vessel lumen diameter (62.46%).

Demarcation of Juvenile and Mature Wood
According to the variation trend of each wood property, since the latewood percentage had no obvious variation pattern with cambial age, it was not used to analyze for demarcation. The variation was divided into three sections using the optimal segmentation model, namely juvenile wood, transition wood and mature wood ( Figure 10). The fiber width (R 2 = 0.92), ring width (R 2 = 0.91) and fiber cell wall thickness (R 2 = 0.90) had larger accuracy, while vessel cell wall thickness (R 2 = 0.50) had smaller accuracy. The model accuracy of the aggregative indicator was 0.76 ( Figure 11). The juvenile period of each wood property was different and the calculation of aggregative indicators showed that the juvenile period was before the 10th year. Among all wood properties, the juvenile period of the fiber lumen diameter and vessel cell wall thickness was shorter (before the 4th year) and the juvenile period of wood air-dry density was longest (before the 14th year). Overall, the juvenile period of fiber properties was relatively shorter (before the 8th year), while the vessel properties gradually completed the juvenile period before the 12th year and entered the transition period. Most properties took 4-14 years to complete the transition period, except for the annual ring width (14 years), vessel length (18 years) and fiber length (28 years). This showed that the continuous change in growth rate, vessel length and fiber length took longer and, as a result, reached maturity relatively late. Overall, the transitional period was from the 11th to the 16th year, and the mature period was after the 16th year through calculation of the aggregative indicators.
There were significant differences among juvenile wood, transitional wood and mature wood for each wood property, showing that the above division was effective. The 12 wood properties showed three main patterns ( Table 6). The first is where the values of juvenile wood were relatively larger, and those of the transitional wood and mature wood gradually became smaller, e.g., annual ring width, wood air-dry density and wood basic density. The second was where the values of juvenile wood were smaller, and the transition wood gradually increased towards mature wood, e.g., the width, lumen diameter and cell wall thickness of the fiber and vessel. The third was where the juvenile wood was smaller, the transition wood increased and the mature wood decreased, e.g., the length of fiber and vessel.     There were significant differences among juvenile wood, transitional wood and mature wood for each wood property, showing that the above division was effective. The 12 wood properties showed three main patterns ( Table 6). The first is where the values of juvenile wood were relatively larger, and those of the transitional wood and mature wood gradually became smaller, e.g., annual ring width, wood air-dry density and wood basic density. The second was where the values of juvenile wood were smaller, and the transition wood gradually increased towards mature wood, e.g., the width, lumen diameter and cell wall thickness of the fiber and vessel. The third was where the juvenile wood was  S. superba is a precious and high-quality broadleaf tree species. Understanding its wood anatomical properties is essential for wood utilization and genetic improvement. The wood air-dry density and basic density were, respectively, 0.73 g·cm −3 and 0.54 g·cm −3 , which is moderately dense, similar to loblolly pine [10,11]. There was little variation in wood air-dry density and wood basic density among rings, indicating that the wood properties of S. superba were uniform and are convenient for processing and utilization. According to the regulations of the International Wood Anatomy Society, a wood fiber larger than 1600 µm is long, so the fiber of S. superba is long [25].

Radial Variation of Wood Anatomical Properties
Variation is the basis of genetic improvement. Effective variation determines the potential for improvement [6]. The annual ring width is the primary indicator representing the radial growth rate of trees and one of the main factors causing radial variation. The results showed that the annual ring width had a larger variation than wood properties. This meant that the growth had greater potential for improvement than wood properties [26]. It may be that ring width is easily affected by environmental factors, and the general wood anatomical features are mainly controlled by genetics. In general, the variation of the vessels was greater than that of the fibers, and vessel lumen diameter and fiber length had the largest variation, which had greater potential for improvement by selective breeding. The radial variation of wood properties in the early stage was greater than that in the later stage, which revealed that cultivating a short juvenile period for wood properties is a strategy for wood property breeding. This may be due to the vigorous division of cambium cells in the juvenile stage and the subsequent decline in the mature stage.
Researching the radial variation pattern of wood anatomical properties is the basic way to understand wood formation and utilization [5]. Previous research had found that wood density increased in the early stage and then remained stable [8,27]. Other research found that wood density was largest in the pith and then decreased and stabilized to a constant value [28][29][30]. In our study, the radial variation of wood density of S. superba belonged to the latter pattern; a decrease followed by stabilization. Wang et al. [31] also found that 20-year-old and 42-year-old S. superba had the same variation pattern within plants. It is worth noting that the wood air-dry density and basic density of 54-year-old specimens were larger than at other ages near the pith. This may be because more extractives were deposited in the wood cell walls of heartwood. Cell death in the heartwood increases with age, and the extractives gradually increase [32].
The fibers and vessels were the main components of broadleaf wood, and their changes are affected by many factors, especially the age of the tree. They showed a certain radial variation pattern with cambial age. A previous study found that the fiber length, fiber cell wall thickness and fiber cell wall to lumen ratio of earlywood and latewood of 13-year-old Catalpa bungei increased significantly with cambial age [5]. Similarly, fiber length and fiber width of 12-year-old poplar were found to gradually increase from the pith to the bark [27]. In another study, the fiber length of 40-year-old Betula pendula increased from the pith to the bark in all discs [33]. Similarly, according to our results, there were two radial variation patterns of fiber and vessel. One pattern showed an increase and a decrease, e.g., fiber and vessel length. These wood properties therefore needed continuous observation until they stabilized. The other pattern showed an increase and then stabilization, e.g., the width, lumen diameter and the cell wall thickness of the fiber and vessel. The variation pattern of the fiber and vessel among plants of different ages was similar to that within plants. This verifies the radial variation of the fiber and vessel. Furthermore, the fiber and vessel varied regularly, and this meant that it was feasible to distinguish the demarcation of juvenile and mature wood according to their radial variation pattern.

Correlation of Wood Anatomical Properties
Wood traits are controlled by multiple factors, causing one effect or one causing multiple effects [34]. Therefore, the correlation among wood properties was analyzed to provide a theoretical basis for comprehensive improvement of wood properties. Concurrently, rapid growth and excellent wood properties are an important goal for genetic improvement in wood [35]. According to our results, there was a significant negative correlation between annual ring width and fiber vessel in the early stage. This indicated that excessive growth would lose the qualities of the fiber and vessel in the early stage. Therefore, both growth and wood anatomical properties should be taken into account in the early improvement, but these correlations can be neglected in the mature stage. This result was also confirmed by Pacheco et al. [14], who found that, in Pinus halepensis, the lower the radial growth rate was, the thicker the cell wall became. In addition, Lasserre et al. [36] also found that a narrower ring width had a longer fiber length and a thicker cell wall in Pinus radiata. Furthermore, cambial age had a persistent and extremely significant negative affect on fiber length and vessel length. This also indicated that the fiber length and vessel length of S. superba were greatly affected by cambial age for a long time. Therefore, continuous observation of fiber length and vessel length was needed to determine their stable period.

Models for Predicting Wood Anatomical Properties
The tree theoretical growth equation can reflect the pattern of tree growth. Logistic and Richards are the most widely used models and have a good biological theoretical basis. Dahlen et al. [10,11] used a four-parameter logistic function to model the variation of specific gravity, tracheid length and width for Pinus taeda, and found that their initial model of cambial age explained 45%, 55% and 46% of the variation, respectively. The effect of disk height and ring width was then incorporated into the fixed model, improving the model fit to 56%, 71% and 57%. Auty et al. [37] used an exponential function of cambial age, with the additional inclusion of ring width and sample height as explanatory variables, to model the wood's basic specific gravity of Pinus sylvestris, and found that the fixed model could explain 57% of the variation, with an RESE of 0.043. There was also a non-linear least square method to simulate wood density. It also achieved the effect of quantifying the extent and pattern of radial variation in wood density [1,13]. Our results showed that the growth model was more suitable for the variation of wood anatomical properties than the nonlinear least square model, which laterally indicated that the wood anatomical properties accord with the "S" curve of the growth pattern. The simulation accuracy of the growth model reached more than 82.4%, showing that the model had relatively high precision.

Demarcation of Juvenile and Mature Wood
The demarcation of juvenile wood and mature wood is of great significance for understanding wood properties and for evaluating wood quality and use value. Different tree species have different boundaries between juvenile and mature wood. The fiber length of Betula pendula of 40-year-old specimens tended to be stable at the 15th growth ring from the pith to the bark, so the 15th ring was assumed to be the demarcation point of juvenile-mature wood [33]. The demarcation of juvenile and mature wood of loblolly pine is based on ring specific gravity, for a juvenile period averaging 4.3 years, followed by an 8.6-year transition period, with mature wood produced by the 13th year [8]. The juvenile wood of Hevea brasiliensis occurred between 40 and 55 mm from the pith [38]. From previous research, it was feasible to distinguish the demarcation using different wood properties from the pith to the bark. In our study, ring width, each wood property and aggregative indicators were used to distinguish between juvenile wood and mature wood. According to the different indicators, the results were different; overall, the juvenile wood was before the 10th year, transition age was between the 11th and the 16th years and the mature age was after the 16th year. Therefore, the wood of S. superba had better workability after the 16th year. Jordan et al. [8] found that the demarcation between juvenile and mature wood of Pinus taeda based on specific gravity was not significantly different among different regions. This indicates that the demarcation between juvenile and mature wood in S. superba is representative.

Conclusions
The wood properties had smaller variation than the growth rate. The vessel lumen diameter and fiber length had greater potential for improvement. The juvenile stage was the best period for improvement. There were three radial variation patterns in the wood properties. One pattern showed a decrease and then a stabilization, e.g., wood air-dry density and wood basic density. Another pattern showed an increase and then a decrease, e.g., fiber length and vessel length. The third pattern showed an increase and then a stabilization, e.g., the width, lumen diameter and the cell wall thickness of the fiber and vessel. Both growth and wood anatomical properties should be taken into account in early improvement (before the 20th year), and they could be improved independently in the mature stage. The growth models of Logistic and Richards were more suitable for the radial variation of wood anatomical properties than the nonlinear least square model. Vessel lumen diameter and fiber cell wall thickness were the key factors affecting wood air-dry density and wood basic density. The transition wood between juvenile wood and mature wood was from the 11th to the 16th year based on aggregative indicators from the pith to the bark at breast height. Therefore, the wood of S. superba had better workability after the 16th year.