Quantifying the E ﬀ ect of Crown Vertical Position on Individual Tree Competition: Total Overlap Index and Its Application in Sustainable Forest Management

: Competition is an essential driving factor that inﬂuences forest community sustainability, yet measuring it poses several challenges. To date, the Competition Index (CI) has generally been the tool of choice for quantifying actual competition. In this study, we proposed using the Total Overlap Index (TOI), a CI in which the Area Overlap (AO) index has been adapted and modiﬁed to consider the “shading” and “crowding” e ﬀ ects in the vertical dimension. Next, based on six mixed forest plots in Xiaolong Mountain, Gansu, China, we assessed the results to determine the TOI’s evaluation capability. Individual-tree simulation results showed that compared to the modiﬁed Area Overlap index (AO M ), the TOI has superior quantiﬁcation capability in the vertical direction. The results of the basal area increment (BAI) model showed that the TOI o ﬀ ers the best evaluation capability among the four considered CIs in mixed forest (with Akaike Information Criterion (AIC) of 1041.60 and log-likelihood (LL) of − 511.80 in the model ﬁtting test, mean relative error of − 28.67%, mean absolute percent error of 117.11%, and root mean square error of 0.7993 in cross-validation). Finally, the TOI was applied in the Kaplan–Meier survival analysis and Cox proportional-hazards analysis. The Kaplan–Meier survival analysis showed a signiﬁcant di ﬀ erence between the low- (consisting of trees with the TOI lower than 1) and high-competition (consisting of trees with the TOI higher than 1) groups’ survival and hazard curves. Moreover, the results of the Cox proportional-hazards analysis exhibited that the trees in the low-competition group only su ﬀ ered 34.29% of the hazard risk that trees in the high-competition group su ﬀ ered. Overall, the TOI expresses more dimensional information than other CIs and appears to be an e ﬀ ective competition index for evaluating individual tree competition. Thus, the competition status quantiﬁed using this method may provide new information to guide policy- and decision-makers in sustainable forest management planning projects.


Introduction
Sustainable forest management (SFM) is gradually replacing traditional forest management and becoming a new trend for forestry [1]. To ensure the rationality of SFM policy and strategy, decision-makers have to deeply understand processes and impact factors of forest dynamics. As one and its neighbors in the horizontal dimension. Some researchers promote the idea that tree competition is related to horizontal and vertical dimensions [5,33]. Thus, a proposed modified Area Overlap index (AO M ) concerns incorporating the tree height ratio parameter to show the relationship between the subject tree and its neighbors in the vertical dimension. However, the actual competition is related to the vertical position of the crowns, not the tree height, due to the light and the space competition being size-asymmetric in the vertical dimension [34,35]. Furthermore, different combinations of crown positions and crown lengths in the vertical dimension may cause the competition status to vary among the trees, even though they share the same overlapped area in the horizontal dimension. For example (Figure 1), the relative position of crowns in the vertical dimension may generate different types of competition (Figure 1a,b), while the subject tree's crown length (Figure 1a,d) and that of its neighbor (Figure 1a,c) also influence the CIs.
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 18 The Area Overlap (AO) index [30] has been widely used for predicting forest dynamics, growth, and yield [23,31]. It has also been adapted and modified by numerous researchers [12,32] for various forest types. The AO index is calculated using the area of the influence zone between the subject tree and its neighbors in the horizontal dimension. Some researchers promote the idea that tree competition is related to horizontal and vertical dimensions [5,33]. Thus, a proposed modified Area Overlap index (AOM) concerns incorporating the tree height ratio parameter to show the relationship between the subject tree and its neighbors in the vertical dimension. However, the actual competition is related to the vertical position of the crowns, not the tree height, due to the light and the space competition being size-asymmetric in the vertical dimension [34,35]. Furthermore, different combinations of crown positions and crown lengths in the vertical dimension may cause the competition status to vary among the trees, even though they share the same overlapped area in the horizontal dimension. For example (Figure 1), the relative position of crowns in the vertical dimension may generate different types of competition (Figure 1a,b), while the subject tree's crown length (Figure 1a,d) and that of its neighbor (Figure 1a,c) also influence the CIs. Trees with different vertical competition statuses that exhibit the same crown overlapping area. The green cylinder represents the neighbor tree, while blue represents the subject tree. Main view: (a) the subject tree lower than a neighbor tree whose crown length is relatively small; (b) the subject tree lower than a neighbor tree and share the space with the neighbor tree; (c) the subject tree lower than a neighbor tree whose crown length is relatively big; (d) the subject tree showed a relatively big crown length and lower than a neighbor tree. Vertical view: (a-d) showed same overlapped area between the subject tree and neighbor tree.
Tree mortality is an important process in forest dynamics and a fundamental factor in SFM. It is influenced by the environment, competition tree size, etc., and is modeled by different tree mortality models. The traditional tree mortality model lacks time-dependent data, methodology in hypothesis testing, and censoring of observations [36]. Some scholars proposed using survival analyses in modeling forest mortaility [37]. Survival analysis was developed in medical science and recognized as the basis of new analytical techniques for forest mortality [38]. Time was defined as a critical component of traditional survival analysis. Due to the inherent lack of death times or death ages of Figure 1. Trees with different vertical competition statuses that exhibit the same crown overlapping area. The green cylinder represents the neighbor tree, while blue represents the subject tree. Main view: (a) the subject tree lower than a neighbor tree whose crown length is relatively small; (b) the subject tree lower than a neighbor tree and share the space with the neighbor tree; (c) the subject tree lower than a neighbor tree whose crown length is relatively big; (d) the subject tree showed a relatively big crown length and lower than a neighbor tree. Vertical view: (a-d) showed same overlapped area between the subject tree and neighbor tree.
Tree mortality is an important process in forest dynamics and a fundamental factor in SFM. It is influenced by the environment, competition tree size, etc., and is modeled by different tree mortality models. The traditional tree mortality model lacks time-dependent data, methodology in hypothesis testing, and censoring of observations [36]. Some scholars proposed using survival analyses in modeling forest mortaility [37]. Survival analysis was developed in medical science and recognized as the basis of new analytical techniques for forest mortality [38]. Time was defined as a critical component of traditional survival analysis. Due to the inherent lack of death times or death ages of individual trees in forest surveys, researchers found that any variables that indicate status changes between the first survey and the re-survey could replace the variable of time or age in Sustainability 2020, 12, 7498 4 of 18 traditional survival analysis [36]. They succeed in using the growth of diameter at breast height (∆DBH) as a hypothetical surrogate for tree age [36,39,40] in refining estimates of tree death.
In this study, we proposed using the Total Overlap Index (TOI), a CI developed by modifying the AO index. To accurately calculate the individual trees' competition in both the horizontal and vertical dimensions, we conceptualized the crown as a three-dimensional entity and evaluated the tree competition based on shading and crowding [41]. Next, we verified the vertical dimension of the TOI's competitive representation and evaluation ability based on six natural forest plots. Finally, we used the TOI and growth in diameter at breast height (∆DBH) in applying survival analysis. In summary, the primary objectives of this study were to: (a) propose the use of the TOI, a new CI based on the AO index, to measure tree competition in the horizontal and vertical dimensions; (b) verify the vertical dimension's representational capability and the assessment ability of the new index in simulative and actual competition evaluation; (c) apply the TOI in Kaplan-Meier survival analysis and Cox proportional-hazards analysis, as well as declare the effect of different-level competition on tree mortality trends.

Study Area and Data Collection
The study area ( Figure 2) was established in Xiaolong Mountain (104 • 23 -106 • 43 E, 33 • 31 -34 • 41 N), which is part of the western Qinling Mountain range located in Gansu Province, China. The area has an average altitude of 1600 m and is characterized by a continental monsoon climate. Furthermore, it lies in a transitional region between the temperate and warm temperate zones, with a mean annual temperature of 7-12 • C, mean annual precipitation of 750 mm, and mean annual sunshine duration of 1520-2313 h, as well as an annual average frost-free period of 175 days. The vertical distribution of forest soil showed a significant difference, which showed a mean soil thickness of 30-60 cm and a mean soil pH of 6.5-7.5. The forest communities mainly consist of Quercus aliena var. acuteserrata, Quercus liaotungensis, Lauraceae. Obtusiloba BI., Pinus armandii Franch. and Tilia tuan Szyszyl. Quercus aliena var. acuteserrata dominates the study area's vegetation and represents >50% of the forest composition. individual trees in forest surveys, researchers found that any variables that indicate status changes between the first survey and the re-survey could replace the variable of time or age in traditional survival analysis [36]. They succeed in using the growth of diameter at breast height (∆DBH) as a hypothetical surrogate for tree age [36,39,40] in refining estimates of tree death. In this study, we proposed using the Total Overlap Index (TOI), a CI developed by modifying the AO index. To accurately calculate the individual trees' competition in both the horizontal and vertical dimensions, we conceptualized the crown as a three-dimensional entity and evaluated the tree competition based on shading and crowding [41]. Next, we verified the vertical dimension of the TOI's competitive representation and evaluation ability based on six natural forest plots. Finally, we used the TOI and growth in diameter at breast height (∆DBH) in applying survival analysis. In summary, the primary objectives of this study were to: (a) propose the use of the TOI, a new CI based on the AO index, to measure tree competition in the horizontal and vertical dimensions; (b) verify the vertical dimension's representational capability and the assessment ability of the new index in simulative and actual competition evaluation; (c) apply the TOI in Kaplan--Meier survival analysis and Cox proportional-hazards analysis, as well as declare the effect of different-level competition on tree mortality trends.

Study Area and Data Collection
The study area ( Figure 2) was established in Xiaolong Mountain (104°23′-106°43′ E, 33°31′-34°41′ N), which is part of the western Qinling Mountain range located in Gansu Province, China. The area has an average altitude of 1600 m and is characterized by a continental monsoon climate. Furthermore, it lies in a transitional region between the temperate and warm temperate zones, with a mean annual temperature of 7-12 °C, mean annual precipitation of 750 mm, and mean annual sunshine duration of 1520-2313 h, as well as an annual average frost-free period of 175 days. The vertical distribution of forest soil showed a significant difference, which showed a mean soil thickness of 30-60 cm and a mean soil pH of 6.5-7.5. The forest communities mainly consist of Quercus aliena var. acuteserrata, Quercus liaotungensis, Lauraceae. Obtusiloba BI., Pinus armandii Franch. and Tilia tuan Szyszyl. Quercus aliena var. acuteserrata dominates the study area's vegetation and represents >50% of the forest composition.   Compared to the pure forest, the mixed forest showed more complexity in the crowns' vertical structure. Herein, to verify the TOI's vertical evaluation capability, we selected six natural mixed forest plots (Table 1). All six plots were established in May 2008 at an altitude of 1700-1800 m and were resurveyed in 2013. Based on the method proposed by Husch et al. [42], each square plot had a diagonal distance of 40 m and an overall area of 800 m 2 . To avoid the edge effect, a 2 m buffer zone was implemented, which divided the square plot into the buffer area and central plot ( Figure 3). In addition, we measured and recorded the DBH (diameter tape, accurate to 0.1 cm), tree height (Blume-Leiss height indicator, accurate to 0.5 m), clear bole height (sliding staff, accurate to 0.1 m), crown diameter (tape measure, accurate to 0.1 m), tree age (increment borer), and the relative coordinates of trees in each plot (total station) of all trees with a DBH ≥ 5 cm. The altitude, aspect, slope, soil layer thickness, and coordinates of each plot were also recorded. The basal area increment (BAI) was calculated based on initial DBH and resurveyed DBH [42]. Compared to the pure forest, the mixed forest showed more complexity in the crowns' vertical structure. Herein, to verify the TOI's vertical evaluation capability, we selected six natural mixed forest plots (Table 1). All six plots were established in May 2008 at an altitude of 1700-1800 m and were resurveyed in 2013. Based on the method proposed by Husch et al. [42], each square plot had a diagonal distance of 40 m and an overall area of 800 m 2 . To avoid the edge effect, a 2 m buffer zone was implemented, which divided the square plot into the buffer area and central plot ( Figure 3). In addition, we measured and recorded the DBH (diameter tape, accurate to 0.1 cm), tree height (Blume-Leiss height indicator, accurate to 0.5 m), clear bole height (sliding staff, accurate to 0.1 m), crown diameter (tape measure, accurate to 0.1 m), tree age (increment borer), and the relative coordinates of trees in each plot (total station) of all trees with a DBH ≥ 5 cm. The altitude, aspect, slope, soil layer thickness, and coordinates of each plot were also recorded. The basal area increment (BAI) was calculated based on initial DBH and resurveyed DBH [42]. The mean dominant height was used to indicate the average tree height in each plot.

Selection of the Subject Trees and Neighbor Trees
To calculate the CIs of each individual tree within each plot, all trees in the central plot (as declared in Section 2.1 and shown in Figure 3) were recognized as subject trees. Based on the method proposed by Bella [23] and Opie [30], we defined the neighbor trees as those whose crown projection

Selection of the Subject Trees and Neighbor Trees
To calculate the CIs of each individual tree within each plot, all trees in the central plot (as declared in Section 2.1 and shown in Figure 3) were recognized as subject trees. Based on the method proposed by Bella [23] and Opie [30], we defined the neighbor trees as those whose crown projection connected to the subject tree's crown projection. The canopy projection relationship between trees was assessed using their relative coordinates and crown diameter, as follows:

TOI Model Algorithm
The Area Overlap index (AO) is calculated using a horizontal area ratio of the influence zone of trees and represents the resource competition among trees in the horizontal direction [12,23]. The equation of the AO index was defined as: where a ij is the overlapped area between crowns of subject tree i and its neighbor tree j, CA i is the crown area of subject tree i, and AO i is the Area Overlap index of tree i. However, the competitive pressure on trees is also present in the vertical dimension. Due to individual differences in crown length and vertical position, subject trees compete for vertical dimensional light and space in different competition levels. Thus, as proposed by Canham et al. [43], we divided the competitive states of individual trees based on the concepts of "shading" and "crowding", and modified the AO index using two crown-relationship-based parameters: shading and crowding modification coefficients. The shading modification coefficient was defined as a representation of the effective photosynthetic thickness and shaded thickness of a crown, while the crowding modification coefficient was defined as a representation of space packing and physical crowding. We defined the TOI's formula as follows: where a ij is the overlapped area between the crowns of subject tree i and its neighbor tree j; CA i is the crown area of subject tree i; θ is the modification coefficient of shading; ϕ is the modification coefficient of crowding; shading ij and crowding ij indicate the shading and crowding models that subject tree i suffers from neighbor tree j, respectively. Both shading ij and crowding ij are the volume ratios in three-dimensional space rather than area ratios in the AO index. Because separating the light radiation transfer is hard [44] and generating a three-dimensional model is challenging, the shading and crowding models are based on two hypotheses: (a) The crown shapes of all the trees were recognized as symmetrical cylinders [45], and (b) the cylinders were considered homogeneous, suggesting that light radiation losses were the same when light passed through crowns of the same thickness.

Shading Modification Cofficient
Light interception can be an effective competition predictor [46]. Neighbor and subject crowns' projected areas, effective photosynthetic thickness, and shaded thickness control the total amount of light radiation that passes through the subject crown. Thus, we improved the AO index by adding the vertical dimensional shading status of subject trees. Moreover, we divided the "shading" competition status into three different groups ( Figure 4): competition status (a)-the neighbor tree's clear bole height is not lower than the subject tree's height; competition status (b)-the neighbor tree's height is Sustainability 2020, 12, 7498 7 of 18 higher, while its clear bole height is lower than the subject tree's height; and competition status (c)-the tree height of the subject tree is not lower than the tree height of its neighbor. A shading coefficient (θ), which demonstrates how these three groups of competition statuses were used to modify the AO index, is expressed as follows: where h 1 and h 3 are the neighbor trees' crown lengths under different competition statuses, and h 2 and h 4 are the lengths from the subject tree's treetop to the neighbor tree's treetop in different competition statuses. θ is a modification coefficient and was calculated based on parameters that could influence the amount of light radiation. For example, in competition status (a), the calculation of θ considered the crown thickness and the thickness of gaps between crowns. In competition status (b), only the neighbor crown's effective thickness was considered. In competition status (c), the subject tree suffered no shading competition pressure from the neighbor tree.
where h1 and h3 are the neighbor trees' crown lengths under different competition statuses, and h2 and h4 are the lengths from the subject tree's treetop to the neighbor tree's treetop in different competition statuses. is a modification coefficient and was calculated based on parameters that could influence the amount of light radiation. For example, in competition status (a), the calculation of considered the crown thickness and the thickness of gaps between crowns. In competition status (b), only the neighbor crown's effective thickness was considered. In competition status (c), the subject tree suffered no shading competition pressure from the neighbor tree.

Crowding Modification Cofficient
The subject tree's growth is also highly influenced by competitive pressure from space packing and physical crowding [41]. These competitive pressures originate from branch interlacing, and may reduce the subject tree's vigor by facilitating foliage and branch loss. To precisely describe these competitive pressures, we used a three-dimensional cylinder based on crown abstraction in the "crowding" model, as opposed to using projected circles in the AO index. Using the crown length and vertical position of each tree, we modified the AO index by adding the trees' vertical crowding status and dividing the "crowding" competition status into two distinct groups ( Figure 4): competition status (d)-where the subject tree's crown shared space with its neighbor in the vertical dimension; and competition status (e)-where the subject tree's crown shared no space with its neighbor in the vertical dimension. To quantitatively describe these different crowding statuses, we

Crowding Modification Cofficient
The subject tree's growth is also highly influenced by competitive pressure from space packing and physical crowding [41]. These competitive pressures originate from branch interlacing, and may reduce the subject tree's vigor by facilitating foliage and branch loss. To precisely describe these competitive pressures, we used a three-dimensional cylinder based on crown abstraction in the "crowding" model, as opposed to using projected circles in the AO index. Using the crown length and vertical position of each tree, we modified the AO index by adding the trees' vertical crowding status and dividing the "crowding" competition status into two distinct groups ( Figure 4): competition status (d)-where the subject tree's crown shared space with its neighbor in the vertical dimension; and competition status (e)-where the subject tree's crown shared no space with its neighbor in the vertical dimension. To quantitatively describe these different crowding statuses, we proposed a crowding coefficient (ϕ), as follows: where h 5 is the length of the area where the crowns overlap, and h 6 is the crown length of the subject tree. We defined ϕ as another modification coefficient that indicates the proportion of shared volume.

Individual-Tree Simulation
Trees suffer competition pressure in the vertical direction. Although the subject tree and its neighbor share the same horizontal overlapping area, the subject tree's competition status may change in response to the crown relationship being altered in the vertical direction. To test the possibility of the TOI expressing changes in the vertical competition status, we proposed individual-tree simulations. In these simulations, a subject tree and a neighbor tree were generated based on vertical crown relationships. The subject tree exhibited a random tree height (from 5.1 to 20 m), as well as a random clear bole height (from 5 to 19.9 m). The neighbor tree had a fixed tree height and clear bole height of 15 m and 10 m, respectively. In addition, we assumed that the horizontal overlapping area between the subject tree and its neighbor was constant to ensure the same horizontal competition status. To show the different effects of the vertical dimension's algorithm on the competition index, we calculated the TOI and a modified area overlap index (AO M ), proposed by Ek and Monserud [24], for the subject tree in each simulation. The formula of AO M was defined as: where AO ij is the area overlap index of subject tree i from neighbor tree j; H i and H j are the tree height of subject tree i and neighbor tree j, respectively; and AO Mi is the modified area overlap index of subject tree i.

Growth Model
Compared with other competition indexes, the TOI can express both the horizontal and the vertical competition status of the trees. To verify the plausibility of adding vertical competition status in CIs to evaluate the actual competition of the trees, we compared the TOI, AO M index, Simple Competition index (Hegyi index) [9], and SCI [11] based on growth models. The formulas of the Hegyi index [9] and SCI [11] are defined as: where dbh i and dbh j are the diameters at breast height (DBH) of subject tree i and neighbor tree j, respectively; U i and C i are the dominance and crowding indexes of subject tree i, respectively; λ W i and λ M i are the weighting factors for the neighbor tree's spatial distribution and species identity, respectively. The basal area increment (BAI) model is a widely used tree growth model, which is highly influenced by competition among trees. The BAI model's model-fitting goodness can be improved by incorporating CIs [47]. Thus, we verified the CIs' evaluation capability by comparing the BAI models after incorporating the AO M , Hegyi index, SCI, and TOI. The BAI models were then applied in our Sustainability 2020, 12, 7498 9 of 18 plots. Based on the Wykoff model [48], the BAI was represented as a function of the independent variables of tree size, site condition, and competition. It is mathematically expressed as follows: The tree size data used in modeling were collected and re-collected at different periods and were recognized as mutually dependent observations that showed a strong autocorrelation. Thus, the mixed linear model was selected for our study to avoid the error caused by data aggregation. Herein, based on the correlation relationship and the biological process associated with tree growth, three fixed variables were selected to represent tree size, including DBH, tree height, and crown-diameter-DBH ratio. In addition, the soil thickness ratio was selected as a representation of the site condition of each plot. The structure-based competition index (SCI) was selected as the benchmark CI in the BAI models to estimate the goodness of different CIs. Furthermore, the tree species was treated as the random variable in the model. The initial BAI model was presented as: where i represents tree i in BAI model; BAI is the five-year basal area increment; DHB is the initial diameter at breast height; H is the tree height; DC is the crown-diameter-DBH ratio; ST OA is the total thickness from the soil's O-layer to the A-layer; ST OB is the total thickness from the soil's O-layer to the B-layer; SCI is the structure-based competition index; ε is the random error; µ sp is the random variable of tree species; and b 0 ∼b 5 are the parameters to be estimated. The within-subject covariance structure effectively eliminated the autocorrelation of mutually dependent observations and improved the accuracy of mixed models [49]. In our study, we compared the initial BAI model, which incorporated no covariance structure, and the initial model, which incorporated the autoregressive covariance structure (AR(1)) based on the fitting goodness parameters-including the Akaike Information Criterion (AIC), Likelihood Ratio Test (LRT), and log-likelihood (LL). The model with better fit goodness was selected as the benchmark BAI model. The comparison among the CIs was based on the benchmark BAI models. By replacing the benchmark CI, a group of BAI models incorporating the abovementioned CIs was established. These models were validated via 600 replications using the 0.632 bootstrap method. This technique, which was proposed by Efron [50], has been recognized as a very effective mixed model cross-validation method [51] and uses a validation dataset that was established based on data resampling. In contrast to the conventional bootstrap method, the 0.632 bootstrap method prediction errors were calculated by a subset of original data not included in the bootstrap sample to eliminate its downward bias [51]. The evaluation was based on mean relative error (MRE), mean absolute percent error (MAPE), and root mean square error (RMSE). The MRE, MAPE and RMSE were calculated as follows: where y i is the observed BAI, y i is the predicted BAI, and n is the number of observations in the bootstrapping validation dataset.

Survival and Hazard Analysis
In our study, we applied the TOI in survival analysis to declare the effects of different levels of competition on tree mortality trends. We firstly divided dead trees in our plots into two different groups. The threshold value used in grouping was calculated based on the TOI value of all trees by K-means clustering. Then, a Kaplan-Meier analysis was performed to explore if mortality curves are significantly different between the two groups. The Kaplan-Meier functions were used to calculate the probability distribution of mortality, which comprises a survival function (Equation (14)) and hazard function (Equation (15)), and is expressed as [36,52]: where S(t) is the probability that death occurs at time T. The time T is at least as great as time t. h(t) is the probability that death occurrs exactly at time t. Finally, a Cox proportional-hazards analysis was performed between two competition groups, which showed the different effects of competition status on tree mortality. The Cox function (Equation (16)) is expressed by the hazard function and can be written as a linear regression: where h(t) is the hazard function, h 0 (t) is the baseline hazard, which is equal to the value of the hazard when covariate x is equal to zero, and b is the coefficient of the covariate variable. The time in our research was replaced by ∆DBH.

Statistical Analysis
The establishment of the TOI, the individual-tree simulation, Spearman correlation analysis, the regression of the growth model, and the stand simulation were all performed in R software (ver-3.6.1) using the basic module. The model fitting and cross-validation in the BAI model were calculated using the "nlme" and "sjstats" packages in R software, while the Kaplan-Meier survival analysis and Cox proportional-hazards analysis were based on the "survival" and "survminer" packages. Th relevant figures in this study were created in the "corrplot" package, "ggplot2" package, and "vioplot" packages in R software, ArcGIS/ArcGIS version 10.2 (ESRI, Redlands, CA, USA) and Origin/OriginPro 2018C (OriginLab Corporation, Northampton, MA, USA).

Vertical Dimension Representational Ability
In the individual-tree simulation, we tested the efficacy of the AO M index and TOI to express vertical competition. In general, the TOI and AO M index value distribution exhibited significant discrepancies throughout the simulations ( Figure 5). The AO M index showed a negative correlation with the subject tree's height and was not influenced by the clear bole height (Figure 5a). The maximum and minimum of the AO M index values appeared when the tree height of the subject tree was 5 and 20 m, respectively. Overall, the AO M value distribution was unable to express vertical changes in the relationship between crowns. The subject tree's TOI value is highly influenced by both tree height and clear bole height (Figure 5b). The minimum TOI value appeared when the subject tree's clear bole height and tree height were (15,20]. Under this condition, the lowest part of the subject crown is higher than the highest part of its neighbor's crown; thus, the subject tree did not suffer competitive pressure from the neighbor's crown. The maximum part of the TOI value appeared when the subject tree's height was [10,15], and the clear bole height was around 10 m. Under this condition, the subject tree suffered both relatively high "shading" and "crowding" pressure from the neighbor tree. Overall, the TOI values' distribution trend indicates that the TOI can express changes in the relationship between crowns.

Growth Model
In selecting the benchmark BAI model's covariance structure, we tested model fitting parameters of the initial BAI models with no covariance structure and with the AR(1) covariance structure, respectively ( Table 2). The initial model with the AR(1) covariance structure showed an AIC value of 1042.83 and an LL value of −512.42. These results indicate that this model performed better than the initial model with no covariance structure, which depicted an AIC value of 1050.27 and an LL value of −517.13. Furthermore, the likelihood ratio test between initial models with the AR(1) covariance structure and no covariance structure exhibited an LRT value of 9.437 (p < 0.01), which indicates a significant improvement after the AR(1) covariance structure was incorporated into the initial model. We also compared models that incorporated the AOM, TOI, SCI, and Hegyi index based on the benchmark BAI model, respectively. Table 3 shows the model fitting and the model validation statistics for the benchmark BAI model and models that incorporated the other three CIs. According to model fitting statistics, the model incorporating the TOI showed the best performance, with an AIC value of 1041. 60

Growth Model
In selecting the benchmark BAI model's covariance structure, we tested model fitting parameters of the initial BAI models with no covariance structure and with the AR(1) covariance structure, respectively ( Table 2). The initial model with the AR(1) covariance structure showed an AIC value of 1042.83 and an LL value of −512.42. These results indicate that this model performed better than the initial model with no covariance structure, which depicted an AIC value of 1050.27 and an LL value of −517.13. Furthermore, the likelihood ratio test between initial models with the AR(1) covariance structure and no covariance structure exhibited an LRT value of 9.437 (p < 0.01), which indicates a significant improvement after the AR(1) covariance structure was incorporated into the initial model. We also compared models that incorporated the AO M , TOI, SCI, and Hegyi index based on the benchmark BAI model, respectively. Table 3 shows the model fitting and the model validation statistics for the benchmark BAI model and models that incorporated the other three CIs. According to model fitting statistics, the model incorporating the TOI showed the best performance, with an AIC value of 1041.

Kaplan-Meier Survival Analysis and Cox Proportional-Hazards Analysis
Based on the K-means cluster results of the TOI values, we divided the dead trees into two groups. The dead trees whose TOI values were lower than 1 belonged to the low-competition group; otherwise, they belonged to the high-competition group. The results of the Kaplan-Meier analysis (Figure 6a) showed that the survival probability of the low-competition group is higher than that of the high-competition group, which indicated that higher competition might induce the increase of trees' death rate. The 50% survival probability showed when ∆DBH reached 0.4 in the high-competition group and 1.1 in the low-competition group (Figure 6a). For both the low-and high-competition groups, the curve of cumulative hazard increases while ∆DBH increases (Figure 6b), which indicated that trees with a faster growth rate faced a higher hazard of death. The coefficient of the Cox analysis (Table 4) was −1.0703 and passed three different tests, indicating that trees in the high-competition group suffer a higher risk of death than trees in the low-competition group. The hazard ratio in the Cox analysis was 0.3429, which indicated that trees in the low-competition group had a reduction of 65.71% (or suffered only 34.29%) of hazard risk compared to the high-competition group.

Kaplan-Meier Survival Analysis and Cox Proportional-Hazards Analysis
Based on the K-means cluster results of the TOI values, we divided the dead trees into two groups. The dead trees whose TOI values were lower than 1 belonged to the low-competition group; otherwise, they belonged to the high-competition group. The results of the Kaplan-Meier analysis (Figure 6a) showed that the survival probability of the low-competition group is higher than that of the high-competition group, which indicated that higher competition might induce the increase of trees' death rate. The 50% survival probability showed when ∆DBH reached 0.4 in the highcompetition group and 1.1 in the low-competition group (Figure 6a). For both the low-and highcompetition groups, the curve of cumulative hazard increases while ∆DBH increases (Figure 6b), which indicated that trees with a faster growth rate faced a higher hazard of death. The coefficient of the Cox analysis (Table 4) was −1.0703 and passed three different tests, indicating that trees in the high-competition group suffer a higher risk of death than trees in the low-competition group. The hazard ratio in the Cox analysis was 0.3429, which indicated that trees in the low-competition group had a reduction of 65.71% (or suffered only 34.29%) of hazard risk compared to the high-competition group.   The coefficient is estimated by Cox regression. A positive coefficient indicates that the high-competition group suffered a higher hazard risk than the low-competition group, while a negative coefficient indicates that the low-competition group suffered a higher hazard risk than the high-competition group; hazard ratio is the ratio of the low-competition group versus the high-competition group; * indicates p < 0.05.

Discussion
As stated above, the competition among trees can be characterized as aboveground and belowground. The aboveground competition primarily focuses on the competition for light and growing space [10,23,24]. Thus, concerning the horizontal dimension, many CIs have been developed to focus on the crown's effective projection area [53] to evaluate the actual light and space competition. However, the light and space competition among trees in the vertical dimension is more complicated and assumed to be size-asymmetric [10]. As such, the vertical dimension has not been legitimately evaluated with previous CIs [10]. For example, a shorter tree undergoes disproportionate pressure from light competition due to the shading effect of taller trees [34]. In this scenario, the taller and shorter trees were exposed to the same magnitude of competitive pressure from each other, which is recognized as symmetrical. However, this competition pressure may require higher resistance from the shorter trees, which is manifested in longitudinal growth, and, subsequently, an asymmetric effect is produced. Thus, incorporating only the size ratio (e.g., DBH ratio according to Bella [23] and volume ratio according to Ek and Monserud [24]), expressing the symmetric competition in a CI may not express the trees' actual vertical competition status.
As demonstrated in individual tree simulation, the AO M index, which incorporated the tree height size ratio, did not reflect the competition status changes when the crown position and crown length changed in the vertical dimension. This is because the function and structure of the crown and bole remain highly divergent. In the horizontal dimension, tree competition could be represented by crown diameter ratio or bole diameter ratio (always represented as diameter at breast height) due their strong correlation relationship. However, in the vertical dimension, trees shared different bole lengths and crown lengths, even though they showed similar tree heights (e.g., neighbor trees in conditions 1 and 2 of Figure 7). Considering the competition mode difference in the horizontal and vertical dimensions, the light and space competition of trees in the vertical direction may not be simply represented by using tree-size-based parameters, but rather by adding crown-relationship-based parameters. For example, with the conditions shown in Figure 7, the subject trees in conditions 1 and 2 suffered different competition pressures in the vertical dimension. This difference could be accurately explained by the TOI (including the crown-relationship-based parameter). However, the CIs that used only the tree-size-based parameter showed no difference between these conditions. Thus, we modified the AO index by incorporating two modification coefficients that reduce the excessive effect of smaller neighbor trees [10] to improve its vertical-dimensional competition evaluation ability. In addition, considering the differences between light and growing space competition, two modification coefficients were based on different calculation methods of the crown relationship. The individual-tree simulation results showed that adding modification coefficients significantly contributed to the TOI's vertical-dimensional competition evaluation ability compared to AO M , and enabled it to express the subject tree's competition status changes when the relationship among the crowns changed in the vertical dimension. Sustainability 2020, 12, x FOR PEER REVIEW 14 of 18 The BAI models that incorporate individual trees' competition status can be more accurate than models that concentrate only on tree size [54]. Therefore, researchers always evaluated CIs by using a BAI model that incorporates CIs [16,55,56]. In this work, the BAI model that incorporated the AOM index showed better performance than the benchmark BAI model. This result is in good agreement with Pretzsch [10], who promoted the idea that the AO-based index always outperforms most other CIs and can perform satisfactorily in most forest types. In addition, the BAI model that incorporated the TOI showed the best model fitting performance and model validation performance. This indicates that the TOI, as a modification of the AO index, was more efficient in quantifying individual trees' actual competition status in mixed forests. This is because adding the TOI, with its crownrelationship-based parameters, can effectively improve the accuracy of the AO index.
The competition also influenced tree mortality in both natural and planted forests [57]. Researchers hold different standpoints for the effect of competition on tree mortality and estimated mortality by different methods [58,59]. In forestry, competition is widely used as the representation of the social positions of trees in mortality models [10]-for example, basal area in lager trees (BAL) in Pukkala's model [60], or crown closure in higher trees (CCH) in Hann and Wang's model [61]. These models could represent the contributions of different variables to tree death, and qualitatively described the effects of different levels of specific variables on tree mortality. For example, Pretzsch [10] reported that trees in lower social positions might have a higher mortality rate than trees in higher social positions, which agreed with our Kaplan-Meier analysis result in mixed forests. On the one hand, trees in the low-competition group live in environments with relatively adequate resources, and have a low probability of dying [62]. On the other hand, trees in the high-competition group showed a lower growth rate and a bigger tree height-diameter ratio, reperesenting a high risk of being attacked by insects or disease and wind, respectively [10]. In addition, even though these traditional tree mortality models could describe the relationship between different levels of specific variables and tree death, the quantitative relationship between different levels of competition groups and tree mortality is still unknown. In this study, we used Cox proportional-hazards analysis based on the value of the TOI to explore the effect of different levels of competition on tree mortality. This method helped in quantifying the effect of a specific variable's different level on tree mortality. As our results showed in mixed forests, trees in the low-competition group only suffered 34.29% of hazard risk of that which trees in the high-competition group suffered. This result indicated that keeping the competition pressure of trees in the lower level (TOI value less than 1) could effectively reduce the death rate of trees in mixed forests.
Overall, due to the incorporation of crown-relationship-based parameters, the TOI showed some improvements over the AOM index in actual competition evaluation capability in mixed forests, and The BAI models that incorporate individual trees' competition status can be more accurate than models that concentrate only on tree size [54]. Therefore, researchers always evaluated CIs by using a BAI model that incorporates CIs [16,55,56]. In this work, the BAI model that incorporated the AO M index showed better performance than the benchmark BAI model. This result is in good agreement with Pretzsch [10], who promoted the idea that the AO-based index always outperforms most other CIs and can perform satisfactorily in most forest types. In addition, the BAI model that incorporated the TOI showed the best model fitting performance and model validation performance. This indicates that the TOI, as a modification of the AO index, was more efficient in quantifying individual trees' actual competition status in mixed forests. This is because adding the TOI, with its crown-relationship-based parameters, can effectively improve the accuracy of the AO index.
The competition also influenced tree mortality in both natural and planted forests [57]. Researchers hold different standpoints for the effect of competition on tree mortality and estimated mortality by different methods [58,59]. In forestry, competition is widely used as the representation of the social positions of trees in mortality models [10]-for example, basal area in lager trees (BAL) in Pukkala's model [60], or crown closure in higher trees (CCH) in Hann and Wang's model [61]. These models could represent the contributions of different variables to tree death, and qualitatively described the effects of different levels of specific variables on tree mortality. For example, Pretzsch [10] reported that trees in lower social positions might have a higher mortality rate than trees in higher social positions, which agreed with our Kaplan-Meier analysis result in mixed forests. On the one hand, trees in the low-competition group live in environments with relatively adequate resources, and have a low probability of dying [62]. On the other hand, trees in the high-competition group showed a lower growth rate and a bigger tree height-diameter ratio, reperesenting a high risk of being attacked by insects or disease and wind, respectively [10]. In addition, even though these traditional tree mortality models could describe the relationship between different levels of specific variables and tree death, the quantitative relationship between different levels of competition groups and tree mortality is still unknown. In this study, we used Cox proportional-hazards analysis based on the value of the TOI to explore the effect of different levels of competition on tree mortality. This method helped in quantifying the effect of a specific variable's different level on tree mortality. As our results showed in mixed forests, trees in the low-competition group only suffered 34.29% of hazard risk of that which trees in the high-competition group suffered. This result indicated that keeping the competition pressure of trees in the lower level (TOI value less than 1) could effectively reduce the death rate of trees in mixed forests.
Overall, due to the incorporation of crown-relationship-based parameters, the TOI showed some improvements over the AO M index in actual competition evaluation capability in mixed forests, and could be effectively used in survival analysis. However, the total numbers of research plots and forest types in our study are limited. Ideally, more forest types and research plots should be included in future three-dimensional competition research.

Conclusions
Competition index research is crucial for understanding actual tree competition status and could provide new evidence in sustainable forest management. In this work, we proposed and verified the TOI, a new CI based on the modification of AO, and applied it in Kaplan-Meier survival and Cox proportional-hazards analyses. The results showed that: (1) compared to the AO M , the TOI has improved the quantization ability in the vertical direction; (2) in mixed forests, the TOI showed a better evaluation capability than the AO M , SCI, and Hegyi index in the fitting goodness test and cross-validation of growth models; (3) the low-competition group (consisting of trees with TOI values of lower than 1) only suffered 34.29% of the hazard risk that the high-competition group (consisting of trees with TOI values higher than 1) suffered in mixed forests. These results substantiate that the TOI showed excellent performance in actual competition evaluation mixed forests and, thus, could offer a more accurate understanding of tree competition, which helps sustainable forest management planners make more informed decisions.