Prediction of Competition Indices in a Norway Spruce and Silver Fir-Dominated Forest Using Lidar Data

: Competitive interactions are important predictors of tree growth. Spatial and temporal changes in resource availability, and variation in species and spatial patterning of trees alter competitive interactions, thus a ﬀ ecting tree growth and, hence, biomass. Competition indices are used to quantify the level of competition among trees. As these indices are normally computed only over small areas, where ﬁeld measurements are done, it would be useful to have a tool to predict them over large areas. On this regard, remote sensing, and in particular light detection and ranging (lidar) data, could be the perfect tool. The objective of this study was to use lidar metrics to predict competition (on the basis of distance-dependent competition indices) of individual trees and to relate them with tree aboveground biomass (AGB). The selected study area was a mountain forest area located in the Italian Alps. The analyses focused on the two dominant species of the area: Silver ﬁr ( Abies alba Mill.) and Norway spruce ( Picea abies (L.) H. Karst). The results showed that lidar metrics could be used to predict competition indices of individual trees ( R 2 above 0.66). Moreover, AGB decreased as competition increased, suggesting that variations in the availability of resources in the soil, and the ability of plants to withstand competition for light may inﬂuence the partitioning of biomass.


Introduction
Tree growth is influenced by several factors, including climate patterns, site conditions, and competition processes [1][2][3][4] Among them, tree competition measures are the main predictors of individual trees' growth [5]. Competition among trees is defined as the negative effects that neighboring trees have on a subject tree. These negative effects depend on the interactions between trees in acquiring limited resources, such as light, water, and nutrients [6,7]. Quantifying the competitive effect of neighboring trees is difficult due to the co-occurring effects of various environmental factors on trees.
Many indices have been developed in numerous studies in order to quantify the level of competition that individual trees experience, and to assess how competition affects growth rates (e.g., [8][9][10][11][12]). Competition models, based on competition indices of single trees, can be classified into two al. [44] predicted tree growth in terms of an increase in height crown area and crown volume using bi-temporal airborne lidar data and they related this to some competition indices.
The objective of this study was to use lidar metrics to predict competition indices and to show how they relate with tree aboveground biomass (AGB). In particular, we focused on two competition indices, one related to height and one to the diameter at breast height (DBH). To the best of our knowledge, no study has explored the possibility of predicting DBH and height competition indices for individual trees detected on lidar data, using lidar metrics extracted both at the plot and ITC level.

Materials and Methods
The workflow in Figure 1 shows the analyses carried out in this work. In the following subsections, we present the dataset analyzed in this study and each step of the workflow.

Study Area
This study was conducted in the Autonomous Province of Trento (Italy), in the municipality of Lavarone (45°57′30.09"N, 11°16′25.17′'E). The study area of approximately 4 km 2 ( Figure 2) presents an altitude between 1200 and 1600 m above sea level, and it is composed of an uneven-aged forest with patches of mixed-species and pure-forest stands. The average number of trees per hectare is 839.4. In particular, dominant tree species are Norway spruce (Picea abies (L.) H. Karst.), about 47% of the total stem volume, and silver fir (Abies alba Mill.), about 36% of the total stem volume. Other tree species are present in the study area, although with a low percentage of the total stem volume: European beech (Fagus sylvatica L.) with about 13%, and European larch (Larix decidua Mill.) and Scots pine (Pinus sylvestris L.) with an overall percentage of about 4%.

Study Area
This study was conducted in the Autonomous Province of Trento (Italy), in the municipality of Lavarone (45 • 57 30.09"N, 11 • 16 25.17 'E). The study area of approximately 4 km 2 ( Figure 2) presents an altitude between 1200 and 1600 m above sea level, and it is composed of an uneven-aged forest with patches of mixed-species and pure-forest stands. The average number of trees per hectare is 839.4. In particular, dominant tree species are Norway spruce (Picea abies (L.) H. Karst.), about 47% of the total stem volume, and silver fir (Abies alba Mill.), about 36% of the total stem volume. Other tree species are present in the study area, although with a low percentage of the total stem volume: European beech (Fagus sylvatica L.) with about 13%, and European larch (Larix decidua Mill.) and Scots pine (Pinus sylvestris L.) with an overall percentage of about 4%.

Lidar Data
Lidar data were acquired in 2015 by an Optech ALTM 3100EA sensor with a maximum scan angle of 21 degrees. The mean point density was 21.5 points per square meter for the first return. Up to four returns per pulse were measured. A digital terrain model (DTM) was generated from the lidar points by the vendor using the TerraScan software with a grid size of 0.5 m. The lidar point cloud was normalized to create a canopy height model (CHM) by subtracting the DTM from the z values of the lidar pulses. This operation was carried out using the software lasground of the package LAStools (https://rapidlasso.com/). The intensity value of each lidar point was range calibrated using the following equation: where I C is the calibrated intensity, I is the raw intensity, R is the sensor-to-target range, and Rs is the reference range or average flying height. An exponential factor, α, of 2.5 was considered since the environmental factors can be considered stable and the same acquisition parameters and instruments were maintained during the survey [45].

Field Data
Inside the study area, 49 circular plots of a 15-m radius were placed (see Figure 2). In summer 2016, within each plot, diameter at breast height (DBH), species, and positions were measured for all the trees with DBH ≥ 7 cm. Height was measured only for certain trees while for the others, it was predicted using height-DBH models defined on the basis of the measured trees. Above-ground biomass was obtained for all trees using local stem volume equations [46,47] multiplied by a conversion factor [48]. In Table 1, a summary of the collected field data and plot characteristics is presented.

Extraction of Competition Indices
Competition dynamics were analyzed using distance-dependent individual competition indices (CIs). These indices provide spatial information of the competitive status of an individual tree [49]. In this study, we used the indices proposed by Hegyi [50] in order to represent different combinations of tree characteristics (i.e., DBH and height).
The first step to calculate individual-based competition indices for a tree (called the subject tree) is the identification of the trees actively competing with it. A search radius around the subject tree needs to be defined: The radius influences the number of competitors to consider and thus the indices. A too small or too large radius may lead to underestimates or overestimates of the competitive effects between the subject tree and its real competitors [49]. In the literature, different methods to define the neighborhood area are presented. In our case, to evaluate the level of competition for each tree, we considered only the competitors located within a certain radius from the subject tree. Previous studies showed that a search radius of 10 m is sufficiently wide to capture all the competitive effects of neighboring trees [51,52]. In particular Szwagrzyk et al. [53] used a radius of 10 m in an area with similar structural parameters to the area in this study. Thus, we decided to use such a value. Moreover, to avoid edge effects, the competition indices were calculated only for trees positioned less than 10 m from any of the plot borders.
Two competition indices were considered in this study, one related to stem DBH and one to tree height. The two indices were calculated using the following formulas [50]: where CI i DBH and CI i H are the DBH and height competition indices for the subject tree, i; DBH i is the diameter at breast height of the subject tree, i; DBH j is the diameter at breast height of the competitor tree, j; H i is the height of the subject tree, i; H j is the height of the competitor tree, j; dist ij is the distance between the subject tree, i, and the competitor tree, j; and n is the number of competitors in the neighborhood zone.
Because the most abundant species in our study area are silver fir and Norway spruce, we considered only these species as subject trees, while all trees were competitors.

ITCs Delineation
The delineation of the ITCs was carried out on the lidar data using the delineation algorithm of the R package itcSegment [54]. In particular, we used the function itcLiDAR. A detailed description of the method adopted can be found in [55]. This approach takes as input the canopy height model (CHM) on which local maxima (i.e., treetops) are located, and around them the crowns of the trees are delineated. The approach included three phases: (1) Smoothing of the canopy height model for which a Gaussian low-pass filter is applied to the rasterized CHM to smooth the surface and to reduce the number of potential local maxima; (2) local maxima extraction: A circular moving window of variable size is applied to the smoothed CHM to find a set of potential treetops (local maxima). A pixel of the CHM is identified as a local maximum when its value is greater than the other values contained in the moving window. The window size is defined according to the height of its central pixel and it spans in a range of odd values defined by the user (e.g., 3,5,7,9). If the height of the central pixel is low, a small value of the window size is used and vice versa; and 3) crown region growing: The crown of a tree is identified by the algorithm through the proximity of the pixels to that particular local maximum. A pixel is considered to belong to a specific region when its vertical distance is less than a percentage of a default difference given by the height of the local maximum. This process is repeated until no pixel is added to a region. Once the region is fully grown, a 2D convex hull is applied, resulting in polygons that represent individual trees (ITCs). The algorithm assigns to each delineated ITC a value of height (i.e., the value of the maximum elevation value of the lidar points inside the ITC) and a value of the crown area derived from the convex hull. The input parameters of the function itcLiDAR used in this study were: Resolution 0.5, MinSearchFilSize 3, MaxSearchFilSize 9, TRESHSeed 0.55, TRESHCrown 0.6, minDIST 5, maxDIST 40, HeightThreshold 2, and cw 1. For each ITC, DBH and AGB were predicted using the equations of Jucker et al. [56] implemented in the R package itcSegment considering the temperate coniferous forest as the biome.
To generate the ITCs dataset to use in the modelling part, a matching process between delineated ITCs and reference ground observations was done. The matching procedure followed two steps: (1) Candidate search: all ground reference trees falling inside an ITC were considered as matching candidates; (2) candidate vote: selected candidates were ranked by their difference in height with the delineated ITCs and their Euclidean distance to the treetop. A distance metric, D, was estimated by considering both parameters to select the best candidate as follows: where x and y denote the locations and h the heights of the field-measured trees and the delineated ITCs, respectively; w is a user-defined weight (set to 0.5 in this study) [57].

Lidar Metrics Extraction
Lidar metrics were extracted for each delineated ITC. As explained before, the competition indices computed on the field tree data refer to the competing trees in a radius of 10 m from the subject tree. Thus, for each ITC, two sets of metrics were defined: (i) plot metrics: 93 elevation and intensity metrics (see Figure 3) extracted from the first and last return of the lidar point cloud data. For each ITC, we considered lidar points located in an area of radius of 10 m from the ITC location. These metrics were extracted from the entire set of points in the 10-m radius from the ITC location, and they did not depend on the ITCs characteristics. Each metric was extracted from the first return points only (_F), and last return points only (_L); and (ii) ITC metrics: 23 metrics computed on the basis of the ITCs located in a radius of 10 m around the considered ITC (see Figure 4). These metrics were computed combining the characteristics (height, crown area, location, DBH) of the ITCs located in a radius of 10 m from the reference ITC.

Prediction
To evaluate the efficiency of lidar metrics in modelling the competition indices, ordinary least square (OLS) models were built, in which the dependent variables were the competition indices estimated in the field and the independent variables were the lidar metrics. We developed three models for each competition index: One model using all the ITCs matched with the field data, one using only the ones of silver fir, and one using only the ones of Norway spruce. Before building the regression models, the presence of multicollinearity between the independent variables (the lidar metrics) was evaluated with the function findCorrelation of the R package caret [58]. In order to reduce the number of lidar metrics and to remove the ones most correlated among each other, we set the correlation threshold to 0.9. After this, the OLS models were built using the stepAIC function of the R package MASS. A natural logarithmic transformation of the original independent variables (the competition indices) was performed in order to avoid non-normality. The stepAIC function was restricted in order to reduce/avoid overfitting of the models. In particular, the selection was restricted in order to have at least 10 field samples for each metric selected and to have a value of the accuracy indices, R2R (Equation (13)) and SSR (Equation (14)), close to one.
The models were validated using a leave-one-out cross-validation and the accuracy statistics used in [59]. In particular, we considered: (1) The mean difference (MD) between the predicted and the observed values: where n is the total number of samples, pre CV i is the predicted value of the sample, i, obtained by cross-validation, and obs i is the observed value of the sample, i.
(2) The mean absolute difference (MAD): (3) The root mean squared differences (RMSDs) of the predicted values: where the SS CV is the sum of the squared differences between the observed values and the predicted values obtained by cross-validation: (4) The coefficient of determination obtained from the models' residuals: where the SS f it is the sum of squares of the model residuals: and SS tot is the sum of squared differences of each observation from the overall mean: (5) The coefficient of determination obtained from the cross-validation: (6) The R 2 ratio: and (7) the sum of squares ratio: Each one of these statistics measures a different aspect of the prediction accuracy: MD measures the prediction bias, MAD and RMSD the prediction precision, R 2 CV the agreement, and R2R and SSR the overfitting. Regarding these last statistics, a desirable value for R2R and SSR in order to avoid overfitting is below 1.1 [59].

Relationship between AGB and Competition Indices
The relationship between competition and AGB was evaluated by developing two linear models between the AGB of the subject trees and the corresponding competition indices. In particular, we built: (i) an OLS model linking the logarithm of the individual trees AGB estimated in the field with the competition indices estimated in the field, and (ii) an OLS model linking the logarithm of the individual trees AGB predicted using lidar with the competition indices predicted using lidar. The logarithm was used in order to avoid non-linearity. The coefficient of determination was used to evaluate the relationships.

ITC Crown Delineation
Out of 464 trees used as subject trees to calculate the competition indices, only 115 matched with a delineated ITC. The detection rate was not very high at only 24.7%, but it is worth noting that, as we wanted to use these data to build up a reliable model, we excluded all the matched trees for which the field height and lidar height differed by more than 2 m. Among the 115 matched trees, 100 belonged to Norway spruce (34 ITCs) and silver fir (66 ITCs), and in the following analyses, only these ones were considered. In Figure 5, a scatterplot of the field-measured/estimated versus lidar-predicted values of DBH and AGB is reported along with the R 2 .
between the AGB of the subject trees and the corresponding competition indices. In particular, we built: (i) an OLS model linking the logarithm of the individual trees AGB estimated in the field with the competition indices estimated in the field, and (ii) an OLS model linking the logarithm of the individual trees AGB predicted using lidar with the competition indices predicted using lidar. The logarithm was used in order to avoid non-linearity. The coefficient of determination was used to evaluate the relationships.

ITC Crown Delineation
Out of 464 trees used as subject trees to calculate the competition indices, only 115 matched with a delineated ITC. The detection rate was not very high at only 24.7%, but it is worth noting that, as we wanted to use these data to build up a reliable model, we excluded all the matched trees for which the field height and lidar height differed by more than 2 m. Among the 115 matched trees, 100 belonged to Norway spruce (34 ITCs) and silver fir (66 ITCs), and in the following analyses, only these ones were considered. In Figure 5, a scatterplot of the field-measured/estimated versus lidarpredicted values of DBH and AGB is reported along with the R 2 . Figure 5. Field-measured DBH and field-estimated AGB versus lidar-predicted DBH and AGB for the 100 detected trees used to build the competition models.

Prediction of Competition Indices
Regarding the models comprising both species, the model ( Table 2) was made up of two plot metrics and four ITC metrics. Among the plot metrics, one was an elevation metric (Zpcum1_F) and one was an intensity metric (Iskew_L). The model (Table 3) was made up of four plot metrics and one ITC metric. Among the plot metrics, two were elevation metrics (Z) and two were intensity

Prediction of Competition Indices
Regarding the models comprising both species, the CI DBH model ( Table 2) was made up of two plot metrics and four ITC metrics. Among the plot metrics, one was an elevation metric (Zpcum1_F) and one was an intensity metric (Iskew_L). The CI H model (Table 3) was made up of four plot metrics and one ITC metric. Among the plot metrics, two were elevation metrics (Z) and two were intensity metrics (I). Regarding the species specific models, for the silver fir (Table 4), the CI DBH model was made up of two elevation plot metrics and one ITC metric, while the CI H model was made up of one elevation plot metric and two ITC metrics. Contrastingly, the models for Norway spruce (Table 5) were made up by plot metrics only: three intensity metrics for the CI DBH model and one elevation metric and two intensity metrics for the CI H model.
Among the ITC metrics, five of them were used in the models: CI_H_ITC, CI_DBH_ITC, DBHsumITC, CAmeanITC, and DsdITC. CI_H_ITC is the height competition index (Equation (3)) computed using only the detected ITCs in the 10-m radius from the subject tree and the lidar-predicted heights; CI_DBH_ITC is the DBH competition index (Equation (2)) computed using only the detected ITCs in the 10-m radius from the subject tree and the lidar-predicted DBH; DBHsumITC is the sum of the predicted DBH values of the delineated ITCs in the 10-m radius around the subject tree; CAmeanITC is the mean value of the crown areas of the delineated ITCs in the 10-m radius around the subject tree; and DsdITC is the standard deviation of the distances among the delineated ITCs in the 10-m radius around the subject tree. Among the elevation plot metrics, two percentile metrics were used (Zq20_F, and Zq95_L) and three cumulative percentage of points metrics (Zpcum1_F, Zpcum2_F, and Zpcum2_L). Six intensity plot metrics were used: Iskew_L, Iskew_F, Isd_L, Imean_F, Ipcumzq50_L, and Ipcumzq90_F.    The accuracy statistics of the models obtained with the leave-one-out cross-validation are shown in Table 6. As can be seen, the behaviors are quite different for each model. Regarding the generic models, the model for the prediction of CI H obtained slightly better results compared to the one for the prediction of CI DBH : All the statistics except for the R 2 f it and R 2 CV obtained better values. Among the species-specific models, the ones for the silver fir competition indices had better performances compared to the ones of Norway spruce. All models experienced quite high values of MAD% (over 25%) and RMSD% (over 36%), and all had a negative bias (negative value of MD). Regarding the overfitting statistics (R2R and SSR), it can be seen that only the generic models have values below 1.1 while for the other models, these values can be higher even if (except for the one for CI DBH of Norway spruce) the values are just slightly above 1.1. Figures 6-8 show the scatterplots between the competition indices predicted by the lidar models and those calculated from the data in the field. Table 6. Summary of models' accuracy statistics. MD: mean difference (Equation (5)). MD: relative MD. MAD: mean absolute difference (Equation (6)). MAD%: relative MAD. RMSD: root mean squared differences (Equation (7)). RMSD%: relative RMSD. R 2 f it : coefficient of determination obtained from the models' residuals (Equation (9)). R 2 CV : coefficient of determination obtained from the cross validation (Equation (12)). R2R: the R 2 ratio (Equation (13)). SSR: sum of squares ratio (Equation (14)).

All Trees
Silver

Relationship between Competition Indices and AGB
In Figure 9 scatterplots of the field-estimated AGB and field-estimated competition indices are shown. This figure shows that high levels of competition resulted in low biomass values while when the competition was below a certain value, it did not influence the biomass. In the figures, the linear models relating the AGB and competition indices are also shown. The model relating log(AGB) to showed an R 2 of 0.65 while the one relating log(AGB) to was 0.44. The same trend was found using the lidar-predicted AGB and competition indices (Figure 10), but with lower values of correlation (R 2 of 0.43 and 0.16, respectively).

Relationship between Competition Indices and AGB
In Figure 9 scatterplots of the field-estimated AGB and field-estimated competition indices are shown. This figure shows that high levels of competition resulted in low biomass values while when the competition was below a certain value, it did not influence the biomass. In the figures, the linear models relating the AGB and competition indices are also shown. The model relating log(AGB) to CI DBH showed an R 2 of 0.65 while the one relating log(AGB) to CI H was 0.44. The same trend was found using the lidar-predicted AGB and competition indices ( shown. This figure shows that high levels of competition resulted in low biomass values while when the competition was below a certain value, it did not influence the biomass. In the figures, the linear models relating the AGB and competition indices are also shown. The model relating log(AGB) to showed an R 2 of 0.65 while the one relating log(AGB) to was 0.44. The same trend was found using the lidar-predicted AGB and competition indices (Figure 10), but with lower values of correlation (R 2 of 0.43 and 0.16, respectively).

Discussion
In this study, we demonstrated that it is possible to predict DBH and height competition indices using lidar metrics. We also showed how competition affects the AGB of individual trees. The results showed that no real improvement is gained in using a species-specific model with respect to a general model. It is worth noting that for the species-specific models, we had quite a low number of samples (especially for Norway spruce) and this could have influenced the results. In terms of overfitting, all models showed reasonable values of R2R and SSR; in particular, only the model for of Norway spruce showed values much above 1.1, which was suggested by Lipovetsky [60] as a desirable limit in order to not have overfitting.
Analysis of the five ITC metrics selected indicated that they are all related to the competition indices. Some of them are clearly related, like the and indices computed using the ITCs (CI_DBH_ITC, and CI_H_ITC), while the others are representative of a part of the competition index equations (DBHsumITC, and CAmeanITC) and of the density of the forest (DsdITC), which is related to competition. Regarding the plot metrics, it was slightly harder to find a direct relation to the competition indices. Metrics based on the distribution of Z are likely related to competition even if

Discussion
In this study, we demonstrated that it is possible to predict DBH and height competition indices using lidar metrics. We also showed how competition affects the AGB of individual trees. The results showed that no real improvement is gained in using a species-specific model with respect to a general model. It is worth noting that for the species-specific models, we had quite a low number of samples (especially for Norway spruce) and this could have influenced the results. In terms of overfitting, all models showed reasonable values of R2R and SSR; in particular, only the model for CI DBH of Norway spruce showed values much above 1.1, which was suggested by Lipovetsky [60] as a desirable limit in order to not have overfitting.
Analysis of the five ITC metrics selected indicated that they are all related to the competition indices. Some of them are clearly related, like the CI DBH and CI H indices computed using the ITCs (CI_DBH_ITC, and CI_H_ITC), while the others are representative of a part of the competition index equations (DBHsumITC, and CAmeanITC) and of the density of the forest (DsdITC), which is related to competition. Regarding the plot metrics, it was slightly harder to find a direct relation to the competition indices. Metrics based on the distribution of Z are likely related to competition even if not directly. The Zq20_F, Zq95_L, Zpcum1_F, Zpcum2_F, and Zpcum2_L metrics describe the vertical distribution of the lidar points, and the vertical distribution of the points is related to the forest structure and density, which are related to competition. In contrast, the intensity metrics could be related to the species. As an example, Imean_F has a quite different distribution of values for the two species: It has a mean value of 23.47 (standard deviation of 3.63) for silver fir compared to 20.15 (standard deviation of 3.82) for Norway spruce.
The effectiveness of lidar metrics in predicting both AGB and competition indices was also found in the study conducted by Lin et al. [43]. In particular, Lin et al. [43] showed that the height competition index estimated by lidar, especially when combined with other lidar metrics (crown radius and height) of the trees, is capable of effectively estimating above-ground carbon (AGC) at both the stand and tree level. In our case, the competition indices were used to assess the influence that high or low competition values have on biomass. The results showed that high competition values led to a decrease in biomass. Therefore, the competitive pressure of neighboring trees is probably an important factor influencing tree growth and biomass partitioning, especially for small trees. Indeed, according to Litton et al. [61] and Poorter et al. [62], biomass partitioning may vary with soil resource availability and with the ability of plants to withstand competition for light. Furthermore, according to the theory of biomass allocation, high competition may increase or reduce biomass allocation in plants [63,64]. Zhou et al. [33] found that the biomass ratio of roots and stems decreased with increasing intensity of competition from neighboring trees, while biomass at the level of branches and leaves increased.
Few previous studies were found in the literature that combined lidar and competition indices. Among the ones present in the literature, the ones of Lo et al. [42], Lin et al. [43], and Ma et al. [44] are the only ones slightly related to this work. In all these studies, competition indices were computed using ITCs automatically delineated on lidar data, in a similar way to our computation of the lidar metrics CI_DBH_ITC and CI_H_ITC. None of these studies analyzed the accuracy of the prediction of competition indices using lidar metrics or validated the predictions using field data. In contrast to the present study, these studies used lidar-predicted competition indices as metrics to predict trees attributes, such as DBH, volume, and carbon density.
Several studies have shown the relationship between radial growth and height growth of trees [65,66] and that the crown:height ratio may quantify competition among trees [67,68]. Therefore, the diameter and height of a tree are not only closely linked to light capture but also to the effects of water, nutrients, and soil conditions [33]. Moreover, according to the results of Zhou et al. [33], during growth, plants change how they are affected by the competition of neighboring plants, and their competitive effect on other plants. This suggests a close relationship between individual competition and tree growth.
In the workflow proposed in this study, some parameters were fixed in a way that could have influenced the final results. The main one was the search radius used to compute the competition indices. Many studies in the literature focusing just on the computation of indices using field data used a different radius for each area analyzed. This is possible if all the tree crowns are measured on the ground; however, that was not our case. Moreover, as we wanted to relate the indices to the lidar data in order to have the possibility of also predicting competition indices in areas not covered by field data, we needed to have a fixed value of the search radius. We chose 10 m as it was used before in other studies investigating forests with similar characteristics [53]. It is worth noting that we also carried out the same analyses using other values, computed in other ways, but the used radius was quite close to 10 m and the final results were very similar or the same.
The ITCs delineation could also have had an effect on the final results. Indeed, the higher the accuracy of the delineation, the more valuable ITC metrics are, and the more trees can be considered as subject trees in the area. The algorithm selected is a simple method when compared with the many algorithms in the literature [69], and has been used successfully in many other studies on forests with similar characteristics (e.g., [47,55]). It was effective for different forest scenarios in a previous study conducted on various forest sites in the Alps [70].
Despite the high potential of lidar technology for the estimation of vegetation parameters, it must be considered that lidar also has limitations. According to Rosette et al. [71], the ability to estimate vegetation parameters (tree height and DTM) decreases in the presence of high terrain slopes and high canopy coverage. Moreover, in very dense forests, it is only able to identify dominant trees, as in our study. Values of competition indices calculated using metrics derived from lidar can be biased due to the fact that small trees are not detected.

Conclusions
Our results showed that lidar metrics have a good capacity to predict competition indices. We developed a system that, after detecting individual tree crowns (ITCs) in the forest, on the basis of lidar metrics extracted in the neighborhood of the detected ITC, predicts two competition indices related to height and diameter at breast height (DBH). From the analyses, it emerged that the use of lidar metrics based on ITCs is important in prediction models. Moreover, we showed that all the information that can be extracted from lidar data should be used, as both plots and ITCs metrics were used in our models. Regarding the relationship between AGB and competition indices, it was observed that the AGB value decreases, increasing competition at both the DBH and height level. These results are probably related to the variation in the availability of soil-level resources and the ability of plants to withstand competition for light. The possibility of predicting competition indices in large areas opens interesting perspectives for forest management practices aimed at regulating species mixture, in particular for forests managed extensively with selective logging, typical of many mountain areas. Funding: This study is linked to activities conducted within the COST (European Cooperation in Science and Technology) Action CLIMO (Climate-Smart Forestry in Mountain Regions-CA15226) financially supported by the EU Framework Programme for Research and Innovation HORIZON 2020.

Conflicts of Interest:
The authors declare no conflict of interest.