Detecting Growth Phase Shifts Based on Leaf Trait Variation of a Canopy Dipterocarp Tree Species (Parashorea chinensis)

Canopy species need to shift their adaptive strategy to acclimate to very different light environments as they grow from seedlings in the understory to adult trees in the canopy. However, research on how to quantitively detect ecological strategy shifts in plant ontogeny is scarce. In this study, we hypothesize that changes in light and tree height levels induce transitions in ecological strategies, and growth phases representing different adaptive strategies can be classified by leaf trait variation. We examined variations in leaf morphological and physiological traits across a vertical ambient light (represented by the transmittance of diffuse light, %TRANS) and tree height gradient in Parashorea chinensis, a large canopy tree species in tropical seasonal rainforest in Southwestern China. Multivariate regression trees (MRTs) were used to detect the split points in light and height gradients and classify ontogenetic phases. Linear piecewise regression and quadratic regression were used to detect the transition point in leaf trait responses to environmental variation and explain the shifts in growth phases and adaptive strategies. Five growth phases of P. chinensis were identified based on MRT results: (i) the vulnerable phase, with tree height at less than 8.3 m; (ii) the suppressed phase, with tree height between 8.3 and 14.9 m; (iii) the growth release phase, with tree height between 14.9 and 24.3 m; (iv) the canopy phase, with tree height between 24.3 and 60.9 m; and (v) the emergent phase, with tree height above 60.9 m. The suppressed phase and canopy phase represent “stress-tolerant” and “competitive” strategies, respectively. Light conditions drive the shift from the “stress-tolerant” to the “competitive” strategy. These findings help us to better understand the regeneration mechanisms of canopy species in forests.


Introduction
Ecological strategies represent adaptations to the intensity of environmental filters involving competition, stress and disturbance [1]. When considered as an operating process for resource allocation, plants' ecological strategies can be expected to shift from seedling to mature stages [1]. Differences in P. chinensis is a typical canopy species with shade-tolerant seedlings; previous studies suggested that irradiance is not a limiting factor for seed germination of P. chinensis, and partial shading favored seedling growth [24]. He et al. [25] reported that leaf morphological traits such as LMA have a positive relation with tree height of P. chinensis in the canopy layer, but we still need more ontogenetic evidence to explain the variations in strategy from understory to canopy. In this study, a canopy crane was used to measure morphology and photosynthesis-related leaf traits of sample trees of P. chinensis in the forest at different life stages. The following hypotheses were tested: (1) insufficient ambient light in the understory leads to a stress-tolerant strategy in P. chinensis, and with the growth in tree height, the increased light conditions induce a shift from stress-tolerant strategy to release growth; (2) with the continuous height growth above the main canopy, hydraulic limitation limits plant growth and changes the adaptive strategy to relatively competitive; (3) all these shifts in adaptive strategy can be detected by variation in leaf traits, and based on the threshold level of light or height, the ontogenetic growth of P. chinensis can be classified into several phases. To test these hypotheses, variations in both morphological traits (LMA, leaf area and leaf thickness) and physiological traits (leaf dark respiration rate (R d ), A max , LCP, LSP, leaf nitrogen content and carbon-nitrogen ratio (C:N ratio)) across a vertical ambient light and tree height (sapling to adult) gradient in P. chinensis were examined, and different growth phases of P. chinensis were also detected by quantitative evaluation of leaf traits in response to height and light variance.

Study Area
This study was performed in a tropical seasonal rainforest in Xishuangbanna, Southwestern China (101 • 34 E, 21 • 36 N). This region is predominated by a typical monsoon climate that alternates between a dry season (November-April) and a rainy season (May-October). The mean air temperature is 21.8°C, and the mean annual precipitation is 1493 mm [20].
The tropical seasonal rainforest in this area is distributed in lowlands, valleys and hills with sufficient water supply [19]. The emergent layer of the forest in the present study was dominated by P. chinensis; its canopy can reach an average height > 45 m and approaches approximately 60 m for individuals in the emergent layer. P. chinensis also showed the largest basal area among all the tree species in the 20 ha Xishuangbanna tropical seasonal rainforest dynamics plot (XTRDP) [20]. The main canopy layer (approximately 30 m high) is dominated by Sloanea tomentosa, Pometia pinnata (formerly Pometia tomentosa) and Semecarpus reticulatus [20].
The canopy crane in the same tropical seasonal rainforest (TCT7015-10E, Zoomlion Heavy Industry, Changsha, China) was built in December 2014, approximately 300 m away from the XTRDP. The maximum operational height of this crane is 80 m (the maximum canopy height around the crane is approximately 60 m), and the jib length is 60 m (Figure 1). Centered at the base of the crane, a 120 m × 120 m ancillary plot was established for plant species inventory in 2014. A total of 6928 individuals of 217 woody tree species with a DBH ≥ 1 cm were recorded in the plot of 1.44 ha, and P. chinensis had the highest relative importance values (21.73%) ( Table S1). The forest structure and species composition in the crane plot are similar to those in the XTRDP [20].
A total of 21 P. chinensis individuals from 0.6 to 66.5 m tall (with diffuse transmittance from 0.30% to 100%) in the plot were sampled and measured from October to November in 2016 and September to October in 2017, at the end of the rainy season, which is also the main growing season for most local plants. The weather had enough sunny days so we could operate the crane and conduct measurements in the field. In total, 34 workdays were actually taken for field measurements. The tree height of each selected individual was measured with a tape measure from the crane. A total of 21 P. chinensis individuals from 0.6 to 66.5 m tall (with diffuse transmittance from 0.30% to 100%) in the plot were sampled and measured from October to November in 2016 and September to October in 2017, at the end of the rainy season, which is also the main growing season for most local plants. The weather had enough sunny days so we could operate the crane and conduct measurements in the field. In total, 34 workdays were actually taken for field measurements. The tree height of each selected individual was measured with a tape measure from the crane.

Light Conditions
The spatial and temporal distribution of light in the understory is quite heterogeneous, and the measurement of direct light becomes difficult for long time spans or with too many samples [26]. Indirect measurements, such as the transmittance of diffuse light (%TRANS, sometimes presented as the diffuse non-interceptance, DIFN), are also acceptable to estimate light conditions. A methodological comparison previously showed that %TRANS measured with LAI-2000 has a good fit with directly measured long-term PPFD in the understory of a tropical forest [26]. In the present study, we measured %TRANS above each sampled tree with LAI-2200 (LI-COR, Inc., Lincoln, NE, USA), an upgraded vision of LAI-2000 (LI-COR, Inc., Lincoln, NE, USA). Open sky measurements were collected on top of the crane, and each individual was measured at the crown top with five repeats. A 270° view cap was used to prevent the crane jib from obstructing the view. The repeated measurements of each individual were used to adjust the open sky measurements with FV2200 software (LI-COR, Inc., Lincoln, NE, USA) and to estimate diffuse transmittance.

Leaf Dark Respiration Rate (Rd)
Five leaves on sun-exposed terminal shoots were selected from the upper-and outermost crown of each individual, and all measurements (including Rd, light-response curve, leaf area, thickness, LMA, leaf nitrogen and carbon) were performed in same leaf. Prior to sunset (approximately 6.00 p.m.) on the day before the respiration measurements, the leaves were covered with thin aluminum foil so that they were not exposed to sunlight until the measurements. The abaxial side was not completely covered to allow free gas exchange. Because the leaves were covered overnight, the effects on respiration of light-enhanced dark respiration [27] and light-induced metabolites and respiratory gene expression [28] were avoided in the measurements.
The Rd was measured using a portable infrared gas analyzer (LI-6400, LI-COR Inc., Lincoln, NE, USA) at ambient humidity of 60%-90%, while the CO2 concentration was maintained at 400 ppm with the built-in CO2 mixer. The block temperature of the LI-6400 instrument was set to the ambient temperature in this study, approximately 26.9 ± 3.5 ℃ (mean value ± standard deviation), and the vapor pressure deficit (VPD) in the cuvette was maintained at 1.1 ± 0.3 kPa. The system stability was assessed by "stableF" in LI-6400 output, which indicated the comprehensive stability of CO2, H2O

Light Conditions
The spatial and temporal distribution of light in the understory is quite heterogeneous, and the measurement of direct light becomes difficult for long time spans or with too many samples [26]. Indirect measurements, such as the transmittance of diffuse light (%TRANS, sometimes presented as the diffuse non-interceptance, DIFN), are also acceptable to estimate light conditions. A methodological comparison previously showed that %TRANS measured with LAI-2000 has a good fit with directly measured long-term PPFD in the understory of a tropical forest [26]. In the present study, we measured %TRANS above each sampled tree with LAI-2200 (LI-COR, Inc., Lincoln, NE, USA), an upgraded vision of LAI-2000 (LI-COR, Inc., Lincoln, NE, USA). Open sky measurements were collected on top of the crane, and each individual was measured at the crown top with five repeats. A 270 • view cap was used to prevent the crane jib from obstructing the view. The repeated measurements of each individual were used to adjust the open sky measurements with FV2200 software (LI-COR, Inc., Lincoln, NE, USA) and to estimate diffuse transmittance.

Leaf Dark Respiration Rate (R d )
Five leaves on sun-exposed terminal shoots were selected from the upper-and outermost crown of each individual, and all measurements (including R d , light-response curve, leaf area, thickness, LMA, leaf nitrogen and carbon) were performed in same leaf. Prior to sunset (approximately 6.00 p.m.) on the day before the respiration measurements, the leaves were covered with thin aluminum foil so that they were not exposed to sunlight until the measurements. The abaxial side was not completely covered to allow free gas exchange. Because the leaves were covered overnight, the effects on respiration of light-enhanced dark respiration [27] and light-induced metabolites and respiratory gene expression [28] were avoided in the measurements.
The R d was measured using a portable infrared gas analyzer (LI-6400, LI-COR Inc., Lincoln, NE, USA) at ambient humidity of 60%-90%, while the CO 2 concentration was maintained at 400 ppm with the built-in CO 2 mixer. The block temperature of the LI-6400 instrument was set to the ambient temperature in this study, approximately 26.9 ± 3.5°C (mean value ± standard deviation), and the vapor pressure deficit (VPD) in the cuvette was maintained at 1.1 ± 0.3 kPa. The system stability was assessed by "stableF" in LI-6400 output, which indicated the comprehensive stability of CO 2 , H 2 O and flow in the cuvette, and the value varied from 0 to 1 (1 = stable, 0 = not) (LI-COR 2011). In this measurement, the stability is 0.9 ± 0.1.
Each leaf was measured for 1 min (recorded every 3 s) the following morning (between 8:00 a.m. and 11.00 a.m.) at ambient temperature, and a mean value of 1 min was used as the R d for each leaf.
After the R d measurement, leaves with their shoots were cut, and the bases were immersed in deaerated water and transported to the laboratory for additional trait measurements.

Light-Response Curve
The swinging of the crane basket made it difficult to measure a light-response curve in situ, because this measurement requires clamping a leaf for more than one hour and the photosynthesis of P. chinensis decreased at midday [29]; the use of detached shoots could eliminate these immediate effects of path length and gravity on gas exchange [30]. In this study, we used detached shoots to isolate the influence of any height-related trends in foliar structure on the various gas exchange parameters, a method applied in previous studies on P. chinensis [31]. Although excised branches lead to a reduction in photosynthetic rate and respiration rate [32], the purpose of this study was to detect the tendency of variation in leaf traits between ontogenetic phases rather than measuring the absolute value of traits, so the systematic bias from excising branches was acceptable. The R d_area in situ and the respiration rate under 0 µmol·m −2 ·s −1 PPFD in the light-response curve (R d_area ') were consistent, with little bias due to excision ( Figure S1).
The light-response curve (A-Q curve) was measured with a LI-6400 instrument equipped with a 6400-02B LED source and CO 2 injector. The cuvette CO 2 was set to 400 ppm, the relative humidity was 40%-70%, the air flow was 500 µmol·s −1 , the leaf temperature was approximately 30°C, the VPD was 1.3 ± 0.4 kPa, and the system stability was 0.9 ± 0.1. After light adaptation under 1000 µmol·m −2 s −1 (the adaptation time depended on the stomatal situation and was often more than 30 min), the PPFD gradient was set as follows: 1800, 1500, 1200, 900, 600, 300, 150, 70, 30, 15, and 0 µmol·m −2 s −1 . The leaves were allowed to equilibrate for at least 3 min under each PPFD gradient based on the stability of stomatal conductance (g s ) (the fluctuation of g s in same leaf in same gradient during 1 min was 1.9 ± 1.7 mmol m −2 ·s −1 ). The instrument was matched before the measurement in each PPFD gradient. The branch samples were measured within 8 h of excision.
The following model was used to analyze the A-Q curve [33]: where A max is the area-based maximal net assimilation rate, α is the efficiency of the quantum yield in weak light, and C 0 is a constant. The LSP and LCP calculations were performed using the following equations [33]:

Leaf Mass per Area (LMA), Leaf Nitrogen, and Carbon
After the gas exchange measurements, the leaf lamina thickness of fresh leaves was determined using a digital caliper (accuracy: 0.01 mm, No. 111-101-40, Guanglu Co., Guilin, China) based on the mean value of three repetitive measurements across the lamina while avoiding major veins. The areas of the leaves were measured with a leaf area meter (LI-3000C, LI-COR Inc., Lincoln, NE, USA). After then, all leaves were dried at 60 • C in an oven for 3 days and the leaf dry mass was measured. The petiole of P. chinensis is short (1-3 cm), and we included it in the measurement of leaf area and leaf mass.
The dried leaves were then pulverized, and the leaf nitrogen and carbon contents were determined with an elemental analyzer (Vario ISOTOPE Cube, Elementar Analysensysteme GmbH, Langenselbold, Germany) at the Central Laboratory of Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences. Photosynthetic nitrogen use efficiency (PNUE) was calculated by the ratio of the area-based maximal net assimilation rate to leaf nitrogen content per area (N area ) [34].

Statistical Analysis
The value of light intensity was log-transformed in this study, for light often shows a curved, asymptotic relationship with height [13] and LMA [5,35].
We used multivariate regression trees (MRTs) to cluster growth phases. This method is applied to analyze complex ecological data that may include imbalanced or missing values, nonlinear relationships between variables, and high-order interactions [36]. Height and ln(%TRANS) were explanatory variables, and leaf traits including LCP, LSP, PNUE, C:N ratio, leaf thickness, LMA, R d_area , A max_area , N area , R d_mass , A max_mass , N mass , and leaf area were response variables in MRTs. All variables were standardized in this analysis.
Regression trees are summarized by their size (the number of leaves) and overall fit. Fit is defined by relative error (RE): the total impurity of the leaves divided by the impurity of the root node (the undivided data). RE gives an over-optimistic estimate of how accurately an MRT will predict new data, and predictive accuracy is better estimated from the cross-validated relative error (CVRE). CVRE varies from zero for a perfect predictor to close to one for a poor predictor [36].
The size of the MRT was selected by cross-validation that generates a series of MRTs with only marginally worse predictive ability than the best predictive MRT. From this series, the smallest MRT within one standard error of the best is often selected (the 1-SE rule; [37]).
Because of the discontinuous variation of light availability with increasing tree height, traits such as LMA and leaf area may be highly correlated with light in the lower crown, but the relationships become weaker in the upper crown [5,38]. Hence, linear approaches such as multiple regression could not be used to distinguish the effects of height versus light on these variables [38]. Therefore, linear piecewise regression was used to model traits vs. height (or ln(%TRANS)) to determine the stage at which the trait value may not respond to increasing height (or ln(%TRANS)) levels. Model parameters are presented in Tables S2 and S3. To better predict the "hump-shaped" variation in traits such as area-based maximal net assimilation rate (A max_area ), leaf area, and the C:N ratio with increasing tree height, quadratic regressions (traits = β 0 + β 1 × height + β 2 × height 2 ) were also used, and the model parameters are presented in Table S4. R program v3.6.0 [39] was used for this data analysis, and MRTs were constructed using the "mvpart" package in R.

Tree Size and Light Availability
Quadratic regressions described the positive relationship between DBH and tree height of P. chinensis very well (height = 2.8598 + 0.8863 × DBH -0.0030 × DBH 2 , r 2 = 0.973, F = 362.1, P < 0.01; Figure 2). Considering that the tree height could better reflect the effect of hydrostatic limitation, we selected tree height instead of DBH to represent tree size in this study.

Statistical Analysis
The value of light intensity was log-transformed in this study, for light often shows a curved, asymptotic relationship with height [13] and LMA [5,35].
We used multivariate regression trees (MRTs) to cluster growth phases. This method is applied to analyze complex ecological data that may include imbalanced or missing values, nonlinear relationships between variables, and high-order interactions [36]. Height and ln(%TRANS) were explanatory variables, and leaf traits including LCP, LSP, PNUE, C:N ratio, leaf thickness, LMA, Rd_area, Amax_area, Narea, Rd_mass, Amax_mass, Nmass, and leaf area were response variables in MRTs. All variables were standardized in this analysis.
Regression trees are summarized by their size (the number of leaves) and overall fit. Fit is defined by relative error (RE): the total impurity of the leaves divided by the impurity of the root node (the undivided data). RE gives an over-optimistic estimate of how accurately an MRT will predict new data, and predictive accuracy is better estimated from the cross-validated relative error (CVRE). CVRE varies from zero for a perfect predictor to close to one for a poor predictor [36].
The size of the MRT was selected by cross-validation that generates a series of MRTs with only marginally worse predictive ability than the best predictive MRT. From this series, the smallest MRT within one standard error of the best is often selected (the 1-SE rule; [37]).
Because of the discontinuous variation of light availability with increasing tree height, traits such as LMA and leaf area may be highly correlated with light in the lower crown, but the relationships become weaker in the upper crown [5,38]. Hence, linear approaches such as multiple regression could not be used to distinguish the effects of height versus light on these variables [38]. Therefore, linear piecewise regression was used to model traits vs. height (or ln(%TRANS)) to determine the stage at which the trait value may not respond to increasing height (or ln(%TRANS)) levels. Model parameters are presented in Tables S2 and S3. To better predict the "hump-shaped" variation in traits such as area-based maximal net assimilation rate (Amax_area), leaf area, and the C:N ratio with increasing tree height, quadratic regressions (traits = β0 + β1 × height + β2 × height 2 ) were also used, and the model parameters are presented in Table S4. R program v3.6.0 [39] was used for this data analysis, and MRTs were constructed using the "mvpart" package in R.

Tree Size and Light Availability
Quadratic regressions described the positive relationship between DBH and tree height of P. chinensis very well (height = 2.8598 + 0.8863 × DBH -0.0030 × DBH 2 , r 2 = 0.973, F = 362.1, P < 0.01; Figure 2). Considering that the tree height could better reflect the effect of hydrostatic limitation, we selected tree height instead of DBH to represent tree size in this study.   Light availability was promoted as a result of increasing tree height up to 27.4 m (r 2 = 0.804, F = 54.14, P < 0.01; Figure 3; Table S2), and the transmittance of diffuse light varied from 0.3% at 3.7 m to 100% at ≥ 21.2 m. When a tree grew higher than 27.4 m, the linear model became nonsignificant (r 2 = 0.262, F = 3.12, P = 0.13).
Forests 2020, 11, x FOR PEER REVIEW 7 of 20 Light availability was promoted as a result of increasing tree height up to 27.4 m (r 2 = 0.804, F = 54.14, P < 0.01; Figure 3; Table S2), and the transmittance of diffuse light varied from 0.3% at 3.7 m to 100% at ≥ 21.2 m. When a tree grew higher than 27.4 m, the linear model became nonsignificant (r 2 = 0.262, F = 3.12, P = 0.13).

Figure 3.
Relationship between height (x) and ln(%TRANS) (y). Lines indicate linear piecewise regression results; solid lines denote significant relationships (P < 0.05), and broken lines denote nonsignificant relationships. Phase partitions were determined based on the results of multivariate regression tree (MRT). Note that a natural log scale is used on the Y-axis.

Growth Phases
Based on the 1-SE rule, the CVRE result generated an MRT size of five leaves (when size of tree = 5, CVRE ± SE = 0.544 ± 0.086, RE = 0.204; Figure 4). Four splits were presented in MRT (     Light availability was promoted as a result of increasing tree height up to 27.4 m (r 2 = 0.804, F = 54.14, P < 0.01; Figure 3; Table S2), and the transmittance of diffuse light varied from 0.3% at 3.7 m to 100% at ≥ 21.2 m. When a tree grew higher than 27.4 m, the linear model became nonsignificant (r 2 = 0.262, F = 3.12, P = 0.13).

Figure 3.
Relationship between height (x) and ln(%TRANS) (y). Lines indicate linear piecewise regression results; solid lines denote significant relationships (P < 0.05), and broken lines denote nonsignificant relationships. Phase partitions were determined based on the results of multivariate regression tree (MRT). Note that a natural log scale is used on the Y-axis.
Although the negative relationship between leaf area and tree height above the main canopy layer appeared to be nonsignificant in the linear piecewise regression (with height < 9.1 m, r 2 = 0.864, F = 45.27, P < 0.01; with height ≥ 9.1 m, P > 0.05; Figure 6c), quadratic regression of leaf area vs. tree height showed a significant effect (r 2 = 0.671, F = 21.35, P < 0.001) and indicated that the highest leaf area appeared at a height of 37.1 m (Table S4).
Leaf area had a significant linear relationship with increasing diffuse light in the understory (with ln(%TRANS) < −0.132 or diffuse transmittance < 87.63%, slope = 10.766, r 2 = 0.601, F = 22.06, P < 0.001). After the tree crown was exposed to high light conditions with diffuse transmittance ≥ 87.63%, the leaf area decreased with increasing ln(%TRANS) (Figure 6d). The quadratic regression detected that leaf area vs. ln(%TRANS) had a lower r 2 (r 2 = 0.449, F = 9.13, P < 0.01) than leaf area vs. tree height (Table S4). Although the negative relationship between leaf area and tree height above the main canopy layer appeared to be nonsignificant in the linear piecewise regression (with height < 9.1 m, r 2 = 0.864, F = 45.27, P < 0.01; with height ≥ 9.1 m, P > 0.05; Figure 6c), quadratic regression of leaf area vs. tree height showed a significant effect (r 2 = 0.671, F = 21.35, P < 0.001) and indicated that the highest leaf area appeared at a height of 37.1 m (Table S4).
Leaf thickness could increase linearly with increasing tree height (with height < 59.6 m, r 2 = 0.595, F = 24.54, P < 0.001 (Figure 6e)), but the linear relationship between leaf thickness and diffuse transmittance was not significant (Figure 6f). Quadratic regression also indicated that leaf thickness and tree height had positive relationship during the whole life history (Figure 6e), and the lowest thickness appeared at 2.50% diffuse transmittance (Figure 6f; Table S4) at the height of around 6.3 m (converted by ln(%TRANS)-height model in Table S2).

Physiological Trait Transition Point of P. chinensis
Rd_area showed the lowest value at 57.3 m in the quadratic model (Figure 7a), and the heightrelated piecewise linear model of Rd_area was divided into two parts at 17.5 m (Figure 7a). Rd_area had Leaf thickness could increase linearly with increasing tree height (with height < 59.6 m, r 2 = 0.595, F = 24.54, P < 0.001 (Figure 6e)), but the linear relationship between leaf thickness and diffuse transmittance was not significant (Figure 6f). Quadratic regression also indicated that leaf thickness and tree height had positive relationship during the whole life history (Figure 6e), and the lowest thickness appeared at 2.50% diffuse transmittance (Figure 6f; Table S4) at the height of around 6.3 m (converted by ln(%TRANS)-height model in Table S2).   (Table S4).
N mass had the highest value at 39.9 m in the quadratic model (r 2 = 0.494, F = 10.77, P < 0.001; Table S4). Piecewise linear regression showed a significant linear increase in N mass with height growth up to 27.4 m (slope = 0.148, r 2 = 0.478, F = 12.9, P < 0.01). Although the quadratic term in the quadratic model of N mass vs. light was nonsignificant (P > 0.05, Table S4), N mass only had a significant linear increase with ln(%TRANS) up to −1.808 (representing 16.4% transmittance; slope = 0.7143, r 2 = 0.337, F = 6.074, P < 0.05), and the relationship between N mass and higher height or light gradient was not significant in piecewise regression (Figure 8e,f).
The tendency of PNUE in terms of both height and light was quite similar to that of A max_mass (Figure 8g,h), and the transition point of the linear piecewise regression in PNUE was 4.6 m in the height model and −3.65 ln(%TRANS) (representing 2.60% transmittance) in the light model. In quadratic regression, the transition point in PNUE was 25.3 m in the height model and 7.23% transmittance in the light model (Table S4).
LCP showed a positive relationship with height ( Figure 9a) and light (Figure 9b) in either the linear or quadratic model. In piecewise regression, LCP was stable before the tree height reached 13.5 m (Figure 9a). Subsequently, it increased linearly with tree height growth (r 2 = 0.706, F = 24.95, P < 0.001). Although LCP increased with improved light conditions and a transition point at 16.40% transmittance of diffuse light (ln(%TRANS) = −1.808) was observed in piecewise regression, no significant piecewise linear relationship was found between LCP and ln(%TRANS) (Figure 9b). increase with ln(%TRANS) up to −1.808 (representing 16.4% transmittance; slope = 0.7143, r 2 = 0.337, F = 6.074, P < 0.05), and the relationship between Nmass and higher height or light gradient was not significant in piecewise regression (Figure 8e,f). The tendency of PNUE in terms of both height and light was quite similar to that of Amax_mass (Figure 8g,h), and the transition point of the linear piecewise regression in PNUE was 4.6 m in the height model and −3.65 ln(%TRANS) (representing 2.60% transmittance) in the light model. In quadratic regression, the transition point in PNUE was 25.3 m in the height model and 7.23% transmittance in the light model (Table S4).
LCP showed a positive relationship with height ( Figure 9a) and light (Figure 9b) in either the linear or quadratic model. In piecewise regression, LCP was stable before the tree height reached 13.5 m (Figure 9a). Subsequently, it increased linearly with tree height growth (r 2 = 0.706, F = 24.95, P < 0.001). Although LCP increased with improved light conditions and a transition point at 16.40% transmittance of diffuse light (ln(%TRANS) = −1.808) was observed in piecewise regression, no significant piecewise linear relationship was found between LCP and ln(%TRANS) (Figure 9b).  (Table S4). High light could significantly promote LSP in quadratic regression (r 2 = 0.572, F = 14.36, P < 0.001), but the piecewise linear regression was nonsignificant (Figure 9d).  (Table S4). High light could significantly promote LSP in quadratic regression (r 2 = 0.572, F = 14.36, P < 0.001), but the piecewise linear regression was nonsignificant (Figure 9d).
Linear piecewise regression showed a C:N ratio transition point at 27.4 m in the height-related model or 16.40% transmittance (or −1.808 in ln(%TRANS)) in the light-related model, but all the linear models were nonsignificant (P > 0.05; Figure 9e,f). By contrast, the relationship of the C:N ratio vs. height could be predicted by a significant quadratic model (r 2 = 0.333, F = 5.99, P < 0.05), and the lowest predicted C:N ratio (18.30) appeared at a height of 37.8 m. In fact, the quadratic model of the C:N ratio vs. height had a higher r 2 than the quadratic model of the C:N ratio vs. ln(%TRANS) (r 2 = 0.198, F = 3.47, P > 0.05).

Relationships between Variations in Leaf Traits
Photosynthetic traits such as respiration rate and LCP play an important role in determining a plant's tolerance to low light levels [40]. P. chinensis has low LCP and R d_area in the understory, which allows its seedlings and treelets to conserve carbon better than trees in the canopy stage and to gain carbon better, having net photosynthesis at lower light levels [41]. The trade-off between low respiration rates and high maximum photosynthetic rates of P. chinensis in habitat shift from understory to canopy reflected its change in adaptive strategy: conservation of carbon in shaded habitats, which indicates a "stress-tolerant" strategy; increased photosynthesis in high light environments, representing a competitive strategy. All the strategies described here were inferred only relative to ontogenetic adaptation in P. chinensis, which only reflected the similarity between the leaf traits in the different phases with those of the Grimes' strategies.
Investing more nitrogen in proteins to limit the light saturation rate of photosynthesis can improve photosynthesis under high light conditions [42]; therefore, P. chinensis has higher N mass and LSP under high light transmittance (Figure 5f). These increased photosynthetic activities could promote leaf respiration rate in the dark [43], leading to a close negative relationship between R d_mass and the light gradient (Figure 8b).
A max_area and LSP increased with tree height in understory but decreased in very tall trees, while A max_mass and PNUE peaked at low %TRANS (3.50% for A max_mass and 2.60% for PNUE) and then monotonically decreased with height (Figure 8c,g). This difference between area-and mass-based characteristics was similar to the findings in a previous report on other tropical canopy trees [44] and can probably be attributed to the finding that area-based characteristics mainly depend on LMA, while mass-based characteristics mainly depend on chlorophyll content per unit dry mass [42].
Although He et al. [25] reported that both leaf thickness and LMA have positive correlations with tree height of P. chinensis in the canopy layer, the increasing tendency with increased tree height differs between leaf thickness and LMA. Increasing height seems to be the main factor affecting leaf area throughout ontogeny (Figure 6e) and is correlated with a linear increase in leaf thickness of P. chinensis from seedlings to very tall trees. This pattern corresponds to an increase in palisade and mesophyll cells [42] and results in greater chlorophyll per unit dry mass with tree height growth. However, the leaf area of P. chinensis displayed a quadratic relation with height growth (Figure 6c), and the piecewise linear regression of LMA vs. tree height showed different slopes beginning at a height of 27.4 m (Figure 6a), which demonstrates that the effect of light conditions on LMA must be taken into account in the understory. Although hydraulic limitations with the height gradient can potentially limit leaf development as indicated by increasing LMA [14,35], high light availability could also lead to thicker leaves as a result of thicker palisade mesophyll cell layers, which maximize light capture [45]. The light gradient in the understory could enhance the sensitivity of LMA to the change in tree height, but after the tree crown overtopped the main canopy layer, the saturated light conditions could not further increase LMA. Both light and height conditions contributed to the ontogenetic variability of LMA.

Growth Phases and Adaptive Strategies of P. chinensis
MRT analysis showed that the leaf trait combinations throughout the ontogeny of P. chinensis could be divided into five phases. The first split was divided by the diffuse transmittance at 18.33%, which indicated the importance of ambient light conditions for leaf trait acclimation in the ontogeny of P. chinensis.
The other three splits were mainly driven by height, at 8.3, 24.3, and 60.9 m. The effects of tree size at 8.3 m may be combined with light acclimation because the height and light are always positively correlated in the understory and within the canopy [5]. The other two splits, at 24.3 and 60.9 m, occurred in trees exposed to a light-saturated environment; thus, other age-specific environmental filters, such as hydraulic limit and reproductive onset, are expected to drive selection towards specific resource-investment strategies for these specific height gradients.
With further combination with the transition points of linear piecewise regression (Tables S2 and S3) (or peak value of the quadratic regression, Table S4) in leaf traits vs. environmental factors, five phases in MRT were redefined as the vulnerable phase, suppressed phase, growth release phase, canopy phase, and emergent phase to better describe the growth strategies in different phases of P. chinensis. Split levels of height and corresponding diffuse transmittance are listed in Table 1. Table 1. Growth phases and changes in the adaptive strategy of leaf traits of P. chinensis. The letters in the strategy column indicate the following: C-relatively "competitive" strategy; S-relatively "stress-tolerant" strategy.  Table S4). Below these values, P. chinensis had a lower PNUE (Figure 8d,h) and A max_mass (Figure 8b,f), probably due to the extremely low light intensity in the understory, which led to a larger fraction of investment of leaf nitrogen in light harvesting [42]. PNUE is quite closely related to the core leaf economic traits [34]. Low PNUE in the vulnerable phase might indicate that P. chinensis seedlings invested much less photosynthetic product into growth. Some studies have reported that minimum light requirements for species in dipterocarp forests is 1.7-5.2 mol m −2 d −1 or 6%-19% for sun [46]. We also observed that the height growth of P. chinensis seedlings after germination remained unchanged under low transmittance in a greenhouse [24], although high mortality caused by drought, herbivory, and missing occurred within 10 months after germination in the field [47], indicating that P. chinensis seedlings were quite vulnerable in this phase.

Growth
The highest PNUE and A max_mass occurred in the transitional zone between the vulnerable phase and the suppressed phase and subsequently decreased with increasing tree height growth and light. A high PNUE indicates that treelet leaves partition more nitrogen to photosynthetic enzymes (such as Rubisco) [48] and/or have lower mesophyll resistance because of large air spaces and thin cell walls within the lamina [49]. Treelet leaves in the suppressed phase had the lowest A max_area (Figure 7c,d) and the lowest N area (Figure 7e,f), similar to saplings in the vulnerable phase, suggesting that the A max_area in treelet leaves in the suppressed phase is severely limited by low nitrogen content [50]. In the suppressed phase, PNUE and A max_mass exhibited a continuous and significant linear relationship with light (or height) variation until the canopy phase. That is, P. chinensis was physiologically flexible in its adaptation to light environmental change in the understory. Previous studies have reported that treelets of P. chinensis can survive in the understory for several years [51], which would be a consequence of this strategy, just as many other tropical canopy species can survive in the understory as seedlings or treelets for decades [52][53][54]. Tall seedlings of Shorea leprosula, S. parvifolia, and S. pauciflora in a Sabah rainforest also showed distinctly lower mortality than seedlings shorter than 50 cm [54]. Su [55] observed that seedlings of Pometia pinnata (formerly, P. tomentosa), another canopy tree species in the Xishuangbanna rainforest, had 99.6% mortality before reaching 1.5 m in height; beyond this height, mortality rapidly decreased and remained very low until physiological death, which Silvertown [56] described as Oscar Syndrome. In conclusion, in the suppressed phase, P. chinensis has conservative leaf economic trait values such as a smaller leaf area, lower A max_area , and lower area-based leaf nitrogen concentration during ontogenesis, and these characters may indicate that the suppressed phase is an S-selected ("stress-tolerant" strategy) phase [10].
P. chinensis had a lower but relatively stable LCP below 16.40% transmittance or 13.5 m in height (piecewise linear regression; Figure 9a,d), which is close to the split point between the suppressed phase and the growth release phase (18.33% transmittance or 14.9 m in height; Table 1). The lower LCP allowed the treelets and saplings to achieve positive carbon gains under lower light intensities [57]. Based on these findings, diffuse transmittance at 18.33% could be suggested as a threshold value dividing the suppressed phase from the growth release phase.
Most leaf traits, including LMA, N area , N mass , LSP, leaf area, R d_area , and R d_mass , showed significant variations with increasing tree height and light intensity in the growth release phase, in accordance with the hypothesis that tree juveniles exhibit greater intraspecific variability in ecological strategies than adults [1]. The rapid change in leaf traits caused by the improvement of the light environment is the most important feature for the growth release phase of P. chinensis and indicates the conversion of the adaptive strategy from the understory to the canopy.
When the tree crown completely entered the main canopy layer at 24.3 m, the growth release phase ended and the canopy phase began. Ambient light in the canopy phase was close to saturation, and hydraulic limitation was achieved because of increasing trunk height, leaf tissue density, and cell wall thickness and decreased mesophyll air space [14]. The decrease in R d_mass (Figure 8a) and increase in N mass (Figure 8e) with height growth were terminated in this phase (transition points: 17.5 m in R d_mass and 27.4 m in N mass , piecewise linear regression; Table 1). Furthermore, leaf area-related traits such as LMA (Figure 6a Table 1). The canopy phase had the largest leaf area and the highest A max_area and N mass throughout the life history of P. chinensis, implying the highest photosynthetic efficiency in this phase. The leaves of P. chinensis in the canopy layer displayed higher A max_mass than other local Dipterocarp species [29]; this species does not suffer from irreversible photoinhibition in the uppermost canopy leaves [31]. Although hydraulic limitation becomes a major determinant of photosynthesis in the canopy layer in other tree species [58], Kenzo et al. [59,60] observed that high leaf nitrogen concentration and a well-developed mesophyll structure help to maintain high A max in the upper-canopy leaves of dipterocarp canopy trees, which explains why the peak value of A max_area occurred in the canopy phase. In short, the canopy phase was more similar to a C-selected ("competitive" strategy) phase, and higher LMA and A max_area in this phase implied conservative growth and maximized resource acquisition to improve competitiveness in stable and productive habitats [6,10].
MRT suggests a final phase (emergent phase) at tree height over 60.9 m, which is correlated with the decrease in some leaf traits associated with carbon-gain, including A max_area (Figure 7c), leaf area (Figure 6c), and leaf C:N ratios (Figure 9e); these traits showed "hump-shaped" patterns with increasing tree height. Quadratic regression predicted that these traits had peak values at specific heights (A max_area at 37.7 m, leaf area at 37.1 m, and leaf C:N ratios at 37.8 m; Table S4) that were already emergent over the main canopy layer (approximately 30 m) and then decreased with the further growth in tree height. This strategy is similar to other mid-to late-successional species [15,61]. Age-and size-dependent declines in leaf nitrogen concentration could lead to a reduction in A max in very tall trees [16,17]; this could be a reason for the decreasing A max_area in the emergent phase. Increased reproductive allocation to nonstructural carbohydrates among very large trees [61] could explain the reduction in nitrogen and A max , and the increase in leaf C:N ratios (Figure 9e) provides potential support for this hypothesis. Former research in P. chinensis reported that reproduction of this species occurs only after its crown emerges from the main canopy, with DBH approximately 30-40 cm [22], and dipterocarps in Malaysia have a reproductive size threshold of approximately 60 cm in primary forest [62]. This empirical evidence indicated that P. chinensis had a reproductive phase only when the individual became huge and tall enough. In addition, increased mechanical perturbation due to herbivory [63] and wind exposure [64] among large canopy trees is also a potential cause for the declines in leaf nitrogen content and photosynthetic capacity with size. Due to the limit of the crane jib (60 m long), we could sample very few individuals of emergent trees in this phase. More experiments are still required to better support our explanations of the emergent phase in the future.
P. chinensis has larger leaves in the canopy phase than in the suppressed phase and vulnerable phase, but in the emergent phase, leaf area again decreases (Figure 6b). Similar phenomena have been reported for canopy tree species in other tropical forests [18]. Energy-demanding and resource-limited old trees could undergo more morphological than biochemical acclimation [65], and a limited leaf area could better control the leaf temperature, rates of carbon dioxide uptake, and water loss for adaptive acclimation to ambient light promotion in the canopy [66]. The decrease in leaf size in the emergent phase may also be a result of the physiological effects of carbon allocation to reproductive structures [18,61].

Conclusions
As a canopy tree species, P. chinensis adapts to light environments both in the understory and canopy throughout its life history. Based on the different strategies of leaf traits during ontogenetic development in response to light and height, five phases of P. chinensis could be recognized: vulnerable phase, suppressed phase, growth release phase, canopy phase, and emergent phase. According to CSR theory [6] and ontogenetic shifts in leaf traits [1,8,10], the suppressed phase and canopy phase could be classified as "stress-tolerant" strategy and "competitive" strategy, respectively, although these strategies were inferred only relative to ontogenetic adaptation in P. chinensis. The vulnerable phase and growth release phase may also not fall into any specific strategies because the trait values in these two phases were quite dynamic: the former was the key phase during which the plant achieved the transformation from seeds to seedlings, and the latter served as a spanning period during which light-related traits overcame suppression by low light in the understory and responded rapidly to change within a narrow height range (14.9-24.3 m) with a large change in ambient light (18.33%-100% transmittance) in the canopy. P. chinensis had lower carbon assimilation in the emergent phase than in the canopy phase, which may be attributed to reproductive allocation strategy, and this deserves more experiments in the future. Our results suggest that this tropical forest canopy species adjusts its adaptive strategy in leaf functional traits to adapt to variation in the ambient light environment and tree size in different life-history stages. The ecological implications of this pattern of change throughout the lifetime of a tree species for the forest growth cycle deserve further investigation in terms of maintaining forest biodiversity.