Comparative analysis of taper models for Pinus nigra Arn. using terrestrial laser scanner acquired data
iForest - Biogeosciences and Forestry, Volume 17, Issue 4, Pages 203-212 (2024)
doi: https://doi.org/10.3832/ifor4525-017
Published: Jul 22, 2024 - Copyright © 2024 SISEF
Research Articles
Abstract
Taper equations are indispensable tools for characterizing the stem profile of trees, providing valuable insights for forest management, timber inventory, and optimal assortments allocation. The recent progress in Terrestrial Laser Scanning (TLS) has revolutionized forest inventory practices by enabling non-destructive data collection. In this study, four taper models from three different model categories were established based on point cloud data of 219 Pinus nigra trees. The taper equations fitted with TLS data were used to predict the diameter at specific stem heights and the total stem volume. The results show that among fitted models, the Max and Burkhart segmented model calibrated by the means of a mixed-effects approach provided the best estimate of the diameter at different heights and the total stem volume evaluated for different diameter at breast height (DBH) classes. In numerical terms, this model estimated the diameter and the volume with a respective overall error of 0.781 cm and 0.021 m3. The predicted profile also shows that above a relative height of 0.7, the diameter error tends to increase due to the low reliability of data collected beyond the base of the crown primarily caused by interference from branches and leaves. Nevertheless, this study shows that TLS technology presents a compelling opportunity and a promising non-destructive alternative for generating taper profiles and estimating tree volume.
Keywords
Taper Equations, Volume Equations, Forest Mensuration, Forest Assessment, Environmental Management, Max and Burkhart, B-Splines, Random Forest, TLS
Introduction
Precision forestry is becoming increasingly important in the face of the urgent challenges posed by climate change and the rising call for sustainable forest management (SFM) practices ([22]). This emerging trend depends on leveraging precise data obtained using advanced technologies to guide well-informed decision-making ([15]). In the realm of SFM, a noteworthy potential exists for value creation through implementing improved practices. Besides the ecological advantages related to enhanced productivity, which contribute to carbon sequestration and alleviate pressure on forests, the adoption of this paradigm shift is fundamentally linked with substantial economic and social value ([10]).
One concrete and increasingly prevalent instance of the potential being realized is precision harvesting (PH), which aims to maximize the efficiency of the harvesting operations while minimizing the impact on the environment and overcoming several limitations inherent to conventional logging practices (e.g., lack of precision and inefficient resource utilization due to lack of tree-level data - [37]). Within the scope of PH, the optimal allocation of assortments is a crucial component of the wood products supply chain and the carbon stock projections since different wood products hold different economical values and carbon storage potential ([5], [32]). Stem taper equations, which predict the change in stem form from ground to tip can accurately disaggregate trees into specific products based on certain specifications like log length and diameter and can help reach a better optimization of wood products ([6], [7], [43]).
In the academic literature, a panoply of taper model formulations and methods of parameter estimation is currently in use and the selection of a proper model form is more important than the actual fitting method ([49]). Following the classification of McTague & Weiskittel ([34]), taper models could be grouped into three main categories: (i) parametric models (e.g., simple-taper, variable-exponent, and segmented equations) which include most taper equations, are calibrated by the means of parametric approaches (e.g., nonlinear least squares - NLS - and nonlinear mixed effects - NLME - methods) and can ensure biologically consistent behavior ([35]). The Max and Burkhart model stands as an instance of this class, particularly within the segmented equations category, and has been demonstrated to exhibit high accuracy in predicting stem taper ([17], [23], [33]). (ii) Semiparametric models (e.g., B-splines and P-splines) which offer greater flexibility in the fit without requiring an extensive addition of parameters, may provide a better representation for complex form species ([26]). The B-spline regression model is a particular instance of semi-parametric model that was, first introduced by Kublin et al. ([24]), and has been successfully used to predict stem taper ([14], [19]). (iii) Nonparametric models (e.g., Random Forest - RF - and Artificial Neural Network -ANN) can offer a strong predictive performance without requiring the testing of statistical hypotheses ([1]). In particular, RF, an ensemble learning method in machine learning, has been gaining popularity in modeling taper by providing precise predictions and demonstrating generalizable performance ([36], [46]).
Developing taper equations requires multiple observations on upper-stem diameters and other covariates collected on standing or felled trees ([18]). Precision forestry in general offers substantial opportunities that may be beneficial for the forest sector including the development of decision tools such as taper equations. The advancements in Terrestrial Laser Scanning (TLS) technology present an unprecedented opportunity for acquiring essential data with high precision and non-destructive methods, thereby significantly reducing the time required compared to conventional approaches ([6], [42]). Data acquisition for TLS-based forest missions can be performed following single or multi-scan approaches. The single-scan approach has the simplest data acquisition setting and the fastest speed, however, a major problem with that scan method is that only parts of the trees are covered due to occlusion effects by other objects. On the other hand, the multi-scan approach appears the most accurate technique and has the potential to fully cover trees within the sampling area, since scans are performed from multiple directions to overcome the occlusion issue ([29]). Many studies using the multi-scan mode reported a stem detection rate between 62.1% and 100% depending on the forest structure and the scanning setup ([28], [31]).
Pinus nigra Arn. is a fast-growing conifer with a wide though fragmented distribution across Europe and Asia Minor, predominantly in mountainous areas. Due to its ecological flexibility, it is one of the most widely used tree species for reforestation worldwide, and it is considered a potential substitute for native coniferous species in Central Europe under future climate change ([48]). Economically, it is one of the most important conifers in Southern Europe, as its wood is highly suitable for general construction, indoor flooring, furniture industry, fuelwood, and paper pulp ([13]). In Italy, Pinus nigra Arn. stands cover an area of 444.785 hectares, representing 4.25% of the total national forest area, with age ranging from 50 to 95 years ([32]). The previously developed yield equations for Pinus nigra Arn. were based on data collected from young trees and have been shown to underestimate the volume of mature trees in the current conditions ([2]). Therefore, the primary aim of this study was to compare the performance of four taper models from tree class categories in predicting diameter over bark (DOB) and total stem volume of Pinus nigra Arn. specifically using data collected from multi-scan approach with TLS. The developed tool will provide forest managers with more accurate estimates of volume and enable better optimization of assortments for industrial purposes.
Materials and methods
Study area and sampling protocol
The study was conducted in the Vallombrosa forest (43° 44′ N, 11° 34′ E), a biogenetic reserve located about 50 km east-southeast of Florence, Tuscany, Italy. The native vegetation of this forest is mainly represented by beech (Fagus sylvatica) at higher elevations, oak-hornbeam stands (Quercus spp. mixed with Carpinus betulus L. and Ostrya carpinifolia Scop.), and chestnut (Castanea sativa Mill.) at lower altitudes ([12]). As part of the main reforestation program that occurred after the World War II, the European black pine (Pinus nigra Arn.) was introduced in the extreme western side of this forest ([3]). Covering 11% of the total forest area, Pinus nigra Arn. stands were specifically studied here.
The Pinus nigra Arn. area was stratified as part of the sampling scheme based on two criteria: stand density and total tree height. Each criterion was divided into three levels to ensure optimal representation of forest stand variability, resulting in a total of nine strata. An airborne lidar flight was utilized to achieve this stratification, providing high-resolution data on forest structure and height. The ALS survey was carried out under leaf-on condition using a Eurocopter AS350 B3 equipped with a LiDAR RIEGL LMS-Q680i sensor. The flight height was 1000 m above ground level. Full-waveform LiDAR data were registered and discretized to a point density of 10 points m-2. Standard procedures for pre-processing ALS data (e.g., outliers and noise removal, classification of ground/non-ground, and computation of height on the ground) were carried out with the LAStools software ([20]) to obtain ALS normalized point clouds. Subsequently, one scan center was randomly placed in each stratum; then, a multi-scan TLS mission was performed on each scan center to collect detailed information on forest characteristics (Fig. 1).
Fig. 1 - Map of the geographical situation of the Vallombrosa forest and the location of Pinus nigra Arn. stands in this forest (Orange color). Stratification of those stands by density (D1-D3) and canopy height (H1-H3; down left sub-figure).
TLS data collection and processing
TLS data were acquired using a FARO Focus 3D × 130 (FARO Technologies Inc., Lake Mary, FL, USA). The instrument uses a phase-shift-based technology with a maximum range of 120 m and acquires data with an azimuth scan angle of 360°. It collects the x, y, and z coordinates, and the intensity of the laser returns with a scan-ranging noise of ± 1 mm.
Using the TLS device, the data acquisition was set to a reasonable scan configuration providing a good trade-off between a sufficient density and the required time for a single scan (resolution of 7.6 mm at 10 m and 1/5 4× overall quality), for a total of about 28 million pulses per scan. With the aim of reducing occlusion due to other obstacles or vegetation and ensuring sufficient coverage, scanning was performed in multiple positions, at least 8 scans per plot subjectively distributed based on the density and structure. Each of the trees in the stand was at least scanned from 3 positions and up to 12 white polystyrene registration spheres (14 cm diameter) were placed in each plot to aid in the digital registration of individual scans.
Fig. 2 - Schematic illustration of the point cloud data exploration. (a): The resulting point cloud after the pre-processing steps from a small region within one stratification unit; (b): identification of a representative tree within the point cloud; (c): zoomed-in view of the stem of the identified tree; (d): removal of branches for stem cleaning; (e): fitting a cylinder to a specific section of the tree; (f): cross-sectional view of the point cloud for the section depicted in panel (e).
Individual scans were merged on the plot level using the automatic registration algorithm implemented in Trimble Real Works® (TRW) software. The software automatically joins overlapping redundant points to create one seamless 3D point cloud suited for the analysis. Details of the operation (settings, criteria, and thresholds) performed by the software are not declared nor accessible. The process attained a very low plot-level registration error, and TRW achieved a high precision scan placement (mean = 2.33 mm, sd = 3.18 mm, max = 8.21 mm). Using the TRW tool, each pine tree stem was manually separated from the rest (i.e., from the ground, the crown, and other non-pine points). Diameter measurements were acquired at approximately 1-meter intervals through the fitting of cylinders to point cloud data (Fig. 2). To determine the volume, Huber’s formula was employed, which entails using the diameter measured at the midpoint of a log segment (eqn. 1):
where V is the log volume, H is the log length, and s0.5 is the mid-section sectional area.
Detailed descriptive statistics are provided for the trees inventoried in this study (Tab. 1). Statistics are derived from stem attributes estimated from point cloud data collected by the TLS device. The diameters and total height of the stems were estimated using the method described earlier, while the stem volume was determined by applying the Huber formula to each of the 1m-length cylinders. In total, 219 trees were extracted from the inventoried area and were systematically allocated to 5 folds. Four folds were exclusively designated for model training, while the remaining fold was reserved for rigorous testing during each iteration of the cross-validation procedure.
Tab. 1 - Mean, standard deviation (SD), and range for tree characteristics from all the samples. (*): Disk refers to a given cross-section of the stem, at a given height level, for which the diameter is measured.
Parameter | Mean | SD | min | max |
---|---|---|---|---|
DBH (cm) | 45.05 | 11.22 | 21.0 | 82.0 |
Total height (m) | 30.7 | 4.0 | 20.5 | 40.5 |
Disk* dob (cm) | 29.3 | 12.5 | 0 | 82.5 |
Disk height (m) | 14.5 | 8.7 | 0.9 | 40.5 |
Crown base height (m) | 21.0 | 4.0 | 11.7 | 31.3 |
Volume (m3) | 2.27 | 1.32 | 0.38 | 8.04 |
Used models
The data acquired by TLS was used to fit 4 taper models from 3 different model categories, namely parametric, semi-parametric, and non-parametric. Thereafter, the fitted models were evaluated in terms of their ability to predict diameters at different heights and the total stem volume as well. In the following, the different used taper models, as well as the evaluation metrics, are explicitly introduced.
Max and Burkhart (1976)
The segmented stem profile model developed by Max & Burkhart ([33]) was used to represent the parametric model’s category. Numerous studies have shown that this model is accurate in estimating upper stem diameters of many species ([21], [38], [39]). The model consists of three quadratic functions representing the stem’s lower, middle, and upper parts grafted together at two points (a1 and a2) with continuous polynomials and continuous first derivatives at each point. In the present study the taper equation was constrained to pass through diameter at breast height (dbh) while it is by construction constrained for the total height (ht).
Max & Burkhart ([33]) taper equation can be expressed as follows (eqn. 2):
where yij = (dij / dbhi)2, dbhi is the dbh of tree i, with i = [1, 2, …, N], N is the number of trees in the sample, dij is the i-th upper-stem diameter over bark (dob) at height hij on tree i, with j = [1, 2, …, ni], ni is the number of diameter measurements for tree i, zij = 1 - (hij/ht) is the complement of the relative height, a1 and a2 are the join points to be estimated from the data, Ik = 1 if z > ai, and 0 otherwise with k = [1, 2], bp’s are regression coefficients with p = [1, 2, 3, 4] and εij is the residual error term associated to the prediction of yij.
To constrain the taper equation to pass through dbh at the breast height, only one regression coefficient is needed to be changed. In this case, the b1 parameter was modified as suggested in Cao ([8]) and replaced by the modified b1* which can be expressed as follows (eqn. 3):
where zbh = 1 - (1.3/ht). To estimate the parameters of eqn. 2, we used the Levenberg-Marquardt algorithm, which is implemented within the “nls” function in the R stats package ([44]).
Mixed-effects model
The stem taper models may show a high degree of multicollinearity, autocorrelation, and heteroscedasticity, all of which can be addressed with appropriate statistical methods. Many studies reported a reduction in the within-tree autocorrelation when including a tree-level random effect ([16], [9]). Using a mixed model for example, the autocorrelation within the stem can be accounted for at each height measurement level.
In the mixed effect formulation, all parameters of eqn. 2 are expressed as fixed effect (population level coefficients) while some parameters may contain additional random effect (for each tree within the population). In the case of the Max & Burkhart ([33]) taper equation, different combinations of the parameters (b1, b2, b3, b4, a1, a2) could be a potential candidate for random parameters.
Here is an expression for the stem profile including the fixed and random effects (eqn. 4):
where a1i and a2i are join points to be estimated from the data for each tree i, Ik = 1 if z > aki, and 0 otherwise with k = [1, 2], and bpi’s are the regression coefficients to be estimated for each tree i, with p = [1, 2, 3, 4].
To define the statistical properties of the terms in eqn. 4, its vectorial formulation presented in the following offers a more simplified notation and will be considered (eqn. 5):
where Yi is a vector of predictions, Xi is a matrix of predictors, Ai is a matrix for the fixed effects, Bi is a matrix for random effects, B is a vector of fixed parameters, Ui is the vector for random effects and εi is a vector of random errors, with the assumptions that εi = N(0,R) and Ui = N(0,D) where R is the diagonal variance-covariance matrix of εi and D is a diagonal variance-covariance matrix of Ui.
Using the “saemix” package in R ([11]), the stochastic approximation expectation maximization algorithm was used to estimate the parameters of eqn. 4.
The volume equation derived from the integration of the Max and Burkhart taper equation can be expressed as follows (eqn. 6):
where zy = 1 - (hy/ht) and zx = 1 - (hx/ht), hx and hy being respectively the lower height of interest (m) and the upper height of interest, K is a conversion factor (= π/40.000), Ik = 1 if zy > ai, and 0 otherwise, Jk = 1 if zx > ai, and 0 otherwise, with k = [1, 2].
B-splines
Semiparametric models are a class of statistical models that provide a flexible compromise between the rigidity of fully parametric models and the complexity of fully non-parametric models. By avoiding rigid assumptions about functional form, semi-parametric models offer the potential to capture the underlying relationship between variables more accurately, particularly for unknown or complex functional forms, while parametric approaches can be used for estimating model coefficients. Basis-splines (B-splines) are a type of semiparametric model that possess favorable numerical properties. To construct a B-spline model, a fixed number of splines, also called polynomials, are connected at specific points known as knots while ensuring the condition of continuity of their second derivatives across the knots.
To formalize, here we consider that the population mean diameter for a given relative height (hr) may be approximated by a B-spline function of degree p and can be expressed as follows (eqn. 7):
where f(hr) is the response of the model which renders the diameter value di at a given hr, βl’s are the parameters of the model, d1 is the number of parameters in the B-spline function with the condition that d1 ≤ k1 + 1 + p; k1 corresponds to the number of internal knots, Bl,p(1) (hr) is the B-spline basis function computed based on the “de Boor recurrence relation”.
Due to the nonlinear nature of stem tapering and the hierarchy of data acquired from sampled trees, semiparametric models can be calibrated more efficiently using mixed approaches, which can also consider the random deviation of an individual tree from the population average. Linear mixed effect models can be used for the calibration if knot values are fixed a priori, which has the advantage of being less computationally expensive and more numerically stable than nonlinear mixed models.
The B-spline based taper equation considering both fixed and random effects can be written as (eqn. 8):
where f(hrij) corresponds to the population mean response, g(hrij) represents the tree-specific deviation from f(hrij), εij is the tree-specific residual error, dij corresponds to the j-th diameter from the i-th sample tree, hrij corresponds to the j-th relative height from the i-th sample tree, βl’s are the fixed regression coefficient of the model, θi,l’s are the random effect parameters to be estimated for each tree, d2 is the number of random parameters in the g function with the condition that d2 ≤ k2 + 1 + p; k2 corresponds to the number of internal knots in the random effect part of the model, Bl,p(2)(hrij) is the B-spline basis function for the random effect part of the model.
The vectorial formulation of eqn. 8 can be expressed as follows (eqn. 9):
where Yi is the vector of diameter predictions at the level of the i-th tree, Xi is the matrix of B-spine basis values for the fixed effect part of the model, Zi is the matrix of B-spine basis values for the random effect part of the model, β is a vector of fixed parameters, θi is the vector of random effect parameters, εi is the vector of residual error with the assumption that εi = N(0, R) and θi = N(0, G) where R is the diagonal variance-covariance matrix of εi and G is a positive non-diagonal matrix of θi. For more details on this method, see Kublin et al. ([24]).
According to the suggestions by Kublin et al. ([24]). The parameters of eqn. 8 were estimated considering knots [0.0, 0.1, 0.75, 1.0] corresponding to the B-spline basis function B(1)(hr) for the population mean function and [0.0, 0.1, 1.0] for B(2)(hr) representing the deviation from the population average. Hence, the regression model depends only on relative height. Measurements of dbh and further diameter measurements of a specific tree are merely used to calibrate the tree specific deviation from the average taper curve by estimating the corresponding random effects. To ensure diameter estimates of zero at tree top, the last spline in both B(1)(hr) and B(2)(hr) are omitted (see Fig. S1 in Supplementary material for a graphical visualization and [24]).
The R package “TapeR” was used to calibrate the model by the Restricted Maximum Likelihood (REML) method (default in the internally used R package “nlme” - [41]) providing the dbh values for each sampling tree of training data. To estimate stem volumes, the implemented numerical integration method included in the “TapeR” package was used ([25]).
Random Forest
Nonparametric methods have been recently used for taper predictions providing reasonable performances ([36], [50]). The present study used Random Forest to estimate the diameter at different heights of the stem. This algorithm works by training many decision trees on random subsets of the features, then averaging out their predictions. Regarding the structure of the model, the diameter was set as the dependent variable, which were predicted by four independent predictors namely dbh, h, ht and h/ht (abbreviations as above). For fine-tuning the model, the grid search method was used to evaluate all possible combinations of hyperparameters (viz., the number of decision trees; the number of features to consider when looking for the best split; the minimum number of samples required to split an internal node) within the search space as well as the cross-validation method to prevent overfitting. The Python library “scikit-learn” (ver. 0.24.2) which includes an implementation of the Random Forest algorithm was used here ([40]).
Fig. 3 - Schematic tree form representation for the three applied models. (a) the Max and Burkhart parametric model; (b) B-splines semi-parametric model; (c) Random Forest non-parametric model.
Since nonparametric methods are not algebraically integrable, including the random forest algorithm, the Monte Carlo integration, which is a numerical integration approach, was self-implemented and used to predict the stem volume for each sample tree (Fig. 3). This technique relies on a random choice of points at which the integrand is evaluated. Monte Carlo integration could be expressed as follows (eqn. 10):
where ha and hb being respectively the lower and the upper heights of interest (m), K is the conversion factor (= π/40.000), hi corresponds to the i-th height on the tree generated randomly from a uniform distribution U(ha, hb), i = [1, 2, …, N], and N is the number of samples generated with the distribution U.
Evaluation of the methods
Taper models should yield unbiased estimates with minimal variance in both diameter outside bark (dob) and stem volume. A single overall measure of bias or residual error of estimation calculated for dob and stem volume does not provide an exclusive indicator of the goodness in evaluating several taper equations for a set of trees of a given species. This holds true when calculating bias metrics at different parts of stems and for different trees can result in a zero or close-to-zero average given that bias can be positive or negative. Alternatively, the standard error which can report the variability of the biases is considered of better reliability. Moreover, when dealing with small sample size, the residual variance can be influenced by the degrees of freedom, making the correlation index, also known as the fitness index, a useful metric to address this concern. All those metrics were used to rank the fitted models.
The most reliable measure of mean bias (B), overall residual variance (SSE) and the fitness index (FI) can be expressed as follows (eqn. 11 - eqn. 13):
where Yi is the actual observation, Ŷi is the predicted value of the actual observation, bar{Y}i is the mean of the actual observations, n is the number of observation and k corresponds to the number of estimated parameters.
As part of our evaluation process, we employed a grouping methodology to assess the performance of the taper models in conjunction with the overall assessment presented earlier. For dob, biases and standard errors were calculated for each 10% step of the height along the stem, while for the total stem volume (from ground to top), the biases and standard errors were calculated by DBH classes. This allowed us to gain a comprehensive understanding of the accuracy of the taper equations across various tree dimensions as well as across different parts within the trees.
Results
Taper prediction
Tab. 2 - Estimates of parameters (± standard errors) for fixed-effects and mixed-effects taper models, based on the validation data set.
Parameters | Fixed-effect model | Mixed-effects model |
---|---|---|
a 1 | 0.6384 ± 0.0242 | 0.71 ± 0.007 |
a 2 | 0.8756 ± 0.0043 | 0.89 ± 0.001 |
b 1 | 0.9577 ± 0.0085 | 0.96 ± 0.014 |
b 2 | -0.2546 ± 0.0166 | -0.25 ± 0.017 |
b 3 | 1.1693 ± 0.2123 | 2.06 ± 0.083 |
b 4 | 25.5313 ± 1.9310 | 38.67 ± 1.428 |
var(a1) | - | 0.006 ± 0.00074 |
var(a2) | - | 0.0001 ± 0.00002 |
var(b1) | - | 0.042 ± 0.0042 |
var(b2) | - | 0.063 ± 0.0065 |
var(b4) | - | 52.250 ± 20.10 |
The model by Max & Burkhart ([33]) was fit using both fixed and mixed-effects approaches. Different mixed-effects combinations of parameters were tested (Tab. S1 in Supplementary material), and the combination involving a1, a2, b1, b2 and b4 produced the lowest values of Akaike’s information criterion (AIC) and Bayesian information criterion (BIC), resulting in the following parameter values (Tab. 2) and model form (eqn. 14):
The mixed B-splines model was fit using the linear mixed-effects approach implemented in the TapeR package and reached convergence. In order to constrain the model predictions to equal zero at the tree top, d1 which is the number of fixed parameters was set to 5, while d2 which corresponds to the number of random parameters was set to 4, omitting the topmost spline, respectively (see Fig. S1 in Supplementary material). The results of B-splines model fitting are represented in Tab. 3.
Tab. 3 - Estimates of parameters (± standard errors) for B-splines model, based on the validation data set.
Parameters | B-splines model |
---|---|
β 1 | 59.10 ± 0.89 |
β 2 | 41.58 ± 0.52 |
β 3 | 31.76 ± 0.22 |
β 4 | 29.27 ± 0.44 |
β 5 | 13.88 ± 0.05 |
var(θ1) | 172.95 ± 13.15 |
var(θ2) | 114.94 ± 10.72 |
var(θ3) | 38.71 ± 6.22 |
var(θ4) | 136.16 ± 11.66 |
Regarding the RF model, after rigorous fine-tuning, the model reached a state of convergence as the number of decision trees stabilized at an optimal value of 267. In addition, the hyperparameter tuning process led to the following optimized configuration: number of features to best split = 3; minimum samples to split internal nodes = 4.
For each taper equation, overall statistics of fit (Bias, SEE and FI) for the entire stem were calculated and are represented in Tab. 4 for dob by model. Among the four models under investigation, discernible trends in bias emerged. The Max and Burkhart models exhibited a noteworthy positive bias, suggestive of their propensity to overestimate diameters. In contrast, the B-spline and Random Forest models displayed a discernible negative bias, indicating a tendency to underestimate diameters. Notably, the Max and Burkhart fixed effect model demonstrated a particularly pronounced inclination towards overestimation. The results suggest that mixed-effect models outperformed models considering only fixed effects. In particular, the Max and Burkhart and B-splines mixed-effects models showed relatively comparable SEE values of 0.78 cm and 0.96 cm, respectively. In contrast, the Max and Burkhart fixed-effect and random forest models had comparable SEE values of 1.90 cm and 2.26 cm, respectively. This comparison holds true for FI as well, although to a lesser extent than for SEE.
Tab. 4 - Total stem fit statistics for the four taper models from cross-validation results. The underlined values refer to the best statistics. M&B is used here as the acronym of Max and Burkhart.
Model | Bias ± SD (cm) |
SEE ± SD (cm) |
FI ± SD |
---|---|---|---|
M&B fixed effect | 0.274 ± 0.292 | 1.90 ± 0.213 | 0.97 ± 0.030 |
M&B mixed effects | 0.027 ± 0.073 | 0.78 ± 0.094 | 0.99 ± 0.021 |
B-splines | -0.015 ± 0.082 | 0.96 ± 0.102 | 0.99 ± 0.006 |
Random forest | -0.010 ± 0.357 | 2.26 ± 0.286 | 0.96 ± 0.046 |
To assess the performance of taper models at different positions throughout the stem, the statistics of fit, derived from the 5-fold cross-validation results, were analyzed by relative height classes. These statistics are documented in Tab. S2 (Supplementary material) and graphically depicted in Fig. 1. To comprehensively assess the taper of all models, the average bias, SEE, and FI were calculated for each model by relative height class along the stem. The results show that models including mixed effects outperformed those considering only fixed effect roughly along all the section heights. Noticeably, a significant increase of SEE was observed for all the models from a relative height class of (70-80%) and was shown to coincide with the crown base height (CBH) for most of the trees (Fig. 4).
Fig. 4 - Distribution of the relative crown base height (in gray color). The distribution fit curve (in red color), The boxplot of the relative crown base height (in blue color).
Bias measurements indicate that the Max and Burkhart fixed-effect model tends to consistently overestimate upper stem diameters throughout the entire stem. On the other hand, the random forest model tends to underestimate upper stem diameters. B-spline-based predictions tend to overestimate upper stem diameters in the lower and upper portions while underestimating in the middle portion. In contrast, predictions from the Max and Burkhart mixed-effects models exhibit an opposing bias pattern in the lower and middle sections but demonstrate similar behavior in the upper portion of the stem.
Volume prediction
The predictions of volume (over bark) by four taper models were compared with each other. Overall statistics of fit (Bias, SEE, and FI) for the total stem volume are represented in Tab. 5. The results show that mixed-effects models outperformed models considering only fixed effect. Specifically, Max and Burkhart and B-splines models had a comparable SEE values of respectively 0.021 and 0.034 m3, while the Max and Burkhart fixed effect and the Random Forest models had comparable SEE values of respectively 0.221 and 0.273 m3. In addition, the overall results show that all the four models tend to generally overestimate the volume as the bias resulted in positive values for all the candidate models. Notably, the Max and Burkhart models (both fixed and mixed effects) exhibited the smallest overprediction, with an absolute bias value of less than 0.01 m3.
Tab. 5 - Total stem volume fit statistics for the four taper models from cross-validation results. The underlined values refer to the best statistics. M&B is used here as the acronym of Max and Burkhart.
Model | Bias ± SD (m3) | SEE ± SD (m3) | FI ± SD |
---|---|---|---|
M&B fixed effect | 0.0085 ± 0.065 | 0.221 ± 0.044 | 0.97 ± 0.207 |
M&B mixed effects | 0.0080 ± 0.001 | 0.021 ± 0.002 | 0.99 ± 0.001 |
B-splines | 0.0240 ± 0.006 | 0.034 ± 0.005 | 0.99 ± 0.005 |
Random forest | 0.0560 ± 0.066 | 0.273 ± 0.066 | 0.91 ± 0.169 |
To assess the performance of stem volume prediction for different tree dimensions, the statistics of fit (SEE, bias, and FI) were analyzed by DBH class categories. These statistics are documented in Tab. S3 (Supplementary material) and graphically depicted in Fig. 6. In general, results show that prediction errors are larger in large diameter trees for both Max and Burkhart fixed effect and Random Forest models, while the prediction error tend to be constant over all the DBH classes for Max and Burkhart mixed effects and B-splines models. The results also demonstrate a strong correlation between volume predictions and volume observations (FI > ~ 0.5) for all four models across all DBH classes, except for the DBH class category (55-65 cm) and (> 65 cm), where both the M&B fixed effect and Random Forest models exhibit limitations (FI ~ 0.42 and 0.21, SEE ~ 0.39 m3 and 0.71 m3, respectively).
Model ranking
A taper-estimating system should provide an unbiased estimation with minimum variance of both diameter outside bark (dob) and total stem volume outside bark. Using the results from Tab. 4 and Tab. 5, it is possible to rank the four estimating systems. The rank sums were created by ranking the performance by dob and stem volume which were equally weighted here in this study. The rank sums reported in Tab. 6 were generated in such a way that each estimating system was assigned a rank separately for every relative height for dob and for every DBH class for the volume. These ranks were also summed for standard errors and biases. As a result, Max and Burkhart mixed-effects system ranked first in most of all the ranking categories. B-splines based system ranked second in estimating diameter and volume by standard error of estimates while random forest system ranked third and Max and Burkhart fixed effect ranked fourth with regards to the same previously discussed estimates.
Tab. 6 - Rank sum analysis for the four estimating systems.
Model Parameters | Max and Burkhart | B-splines | Random forest | |
---|---|---|---|---|
Fixed effect | Mixed effects | |||
Overall SEE for diameter and volume | 6 (3) | 2 (1) | 4 (2) | 8 (4) |
Overall bias for diameter and volume | 6 (3) | 4 (1) | 5 (2) | 5 (2) |
SEE for diameter by relative height | 31 (3) | 17 (1) | 16 (2) | 36 (4) |
Bias for diameter by relative height | 29 (3) | 21 (2) | 30 (4) | 20 (1) |
SEE for volume by DBH classes | 15 (3) | 6 (1) | 8 (2) | 18 (4) |
Bias for volume by DBH classes | 16 (4) | 9 (1) | 13 (3) | 12 (2) |
Total | 103 (4) | 59 (1) | 76 (2) | 99 (3) |
Discussion
Taper equations are an indispensable tool in the field of forest management and industry. However, the development of an accurate taper model present inherent challenges, mainly due to the need for comprehensive customization to account for the characteristics exhibited by each species and site-specific environmental conditions, and the availability of high-quality data required for the model development. In the present study, four taper models from three different class categories were compared in terms of their capacity to predict diameter over bark (dob) and total stem volume of Pinus nigra Arn. Across models, the Max and Burkhart and B-splines mixed-effects models demonstrated superior accuracy and reliability in predicting both dob and total stem volume, with a small advance for the Max and Burkhart model.
The Terrestrial Laser Scanner (TLS) has previously been employed as an effective tool for measuring trees in forest environments and generating taper equations ([45], [47]). The case study conducted on black pine forests, which is presented here, demonstrated the feasibility of using TLS for surveys and highlighted the high quality of the acquired data. In the Italian context, where tree measurements for quantitative purposes have significantly decreased in recent decades, technological tools like TLS represent a pivotal advancement for rejuvenating our understanding of the Italian forest sector.
In the present study, biases and standard errors of estimate were evaluated for diameter outside bark (dob) and total stem volume. These biases and standard errors were also evaluated for different heights within the trees and across different tree sizes. The overall goodness-of-fit statistics (Bias, SEE, and FI) were calculated and reported in Tab. 4, Tab. 5 and Fig. 5, Fig. 6 for the four tested models. The results indicate that the Max and Burkhart and B-splines mixed-effects models explained more than 99% of the total variation in predicting the upper-stem diameter and total stem volume. The Max and Burkhart fixed-effect model explained more than 99% of this total variation, while the Random Forest model explained more than 96% for the diameter variation and 96% for the total volume variation.
Fig. 5 - Bias, SEE and FI of the estimate of diameter over bark (dob) by relative class for the four models.
The superiority of models that include random effects over those that only consider fixed effects has been demonstrated in various previous studies, and our findings are in agreement with this. For example, a study conducted by Leites & Robinson ([27]) showed that the inclusion of random effects considerably improved the fit of the Max and Burkhart taper equation for Pinus taeda. In a different study by Bronisz & Zasada ([4]), the Kozak’s taper equation fitted using a mixed-effects approach provided better results for the upper diameter and total stem volume for Pinus sylvestris. In our study, the Max & Burkhart ([33]) model was fitted using both fixed and mixed effects approaches. Notably, the error was reduced by 52% for the upper-stem diameter and by 89% for the total stem volume when random effects were additionally considered.
It is also noteworthy that the evaluation of the proposed models exhibited a higher Standard Error of Estimation (SEE) in predicting the upper-stem diameter at a relative height of [0.7, 0.8] (Tab. S2 - Supplementary material). The analysis of the distribution of the tree’s relative crown base height reveals that 50% of the observations are located between a value of 0.7 and 0.8, and the maximum relative crown base height was found to be approximately equal to 0.8 (Fig. 4). This implies that all observations above a relative height of 0.8 were necessarily measured inside the crown. This last result may partly explain the notable increase in error in predicting the upper-stem diameter at a relative height of 0.7. In fact, many studies have reported the inability of Terrestrial Laser Scanners to collect reliable data beyond the base of the crown due to interference from branches and leaves as well as the mutual occlusion between canopies. Practically, when approaching higher parts of the stem, occlusion starts to deteriorate the quality of cylinder fitting ([30], [42]). Additionally, as the distance from the scanner increases, spatial resolution decreases, which is a well-known effect in scanning that particularly impact the upper parts of trees ([29]). The interference from the crown, compounded with the increased distance, may explain the models’ inability to detect a well-defined taper pattern above the crown base.
In this study, Monte Carlo numerical integration was employed for estimating the volume of tree stems using the non-differentiable and complex Random Forest model. This decision was justified by the model’s inherent characteristics mentioned earlier, which render traditional integration methods less suitable. It is imperative to underscore the significance of comparing our chosen method with alternative numerical integration techniques to gain insights into result robustness and accuracy, potentially informing future research. Notably, the RF model, unlike the parametric and semi-parametric models used in this study, typically demands a substantial amount of training data. In our analysis, the RF model obtained the lowest score in terms of SEE (dob and volume). This highlights the potential impact of our relatively small dataset (219 trees) on RF performance. In addition, the use of various integration methods (Huber, Monte Carlo, and analytical integration) in our study may introduce minor biases affecting our results. Further research is needed to assess and potentially mitigate these biases, as their combined impact on the findings remains uncertain.
Considering the results of this study, the Max & Burkhart ([33]) taper model, calibrated by means of a mixed-effects approach, and its compatible volume integration could be considered for operational use in estimating diameter outside bark at different heights, as well as total stem volume, respectively, with an overall precision of 0.781 cm and 0.021 m³ in Pinus nigra Arn.
Conclusion
This study explored the application of precise and non-destructive techniques, such as Terrestrial Laser Scanning, for collecting data to develop taper equations specific to Pinus nigra Arn., in line with the principles of precision forestry. The collected data were used to fit four taper models from three different class categories, evaluated for their ability to predict diameter at various heights and total stem volume. The Max & Burkhart ([33]) taper model, calibrated with a mixed-effects approach, proved superior in estimating diameter and volume compared to other models. However, a common limitation among all models was their inability to accurately predict diameter beyond a relative height of 0.7 due to crown interference. The selected model, within the scope of this study, can be operationally utilized by forest managers to enhance predictions of diameter and stem volume, aiding in optimizing wood allocation to different harvested products. Future research perspectives may focus on refining algorithms and methods for extracting diameter at the upper parts of the stem.
Data, code availability and setup
The whole dataset is made available at ⇒ https://doi.org/10.5281/zenodo.8414408. The study was performed on a laptop with Win10 22H2 operating system, an Intel Core i5 7200U CPU with maximal capacity of 2.70 GHz, and 8 Gb of RAM.
Authors contribution
NP and CV contributed equally to the manuscript.
References
Gscholar
CrossRef | Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Authors’ Info
Authors’ Affiliation
Riccardo Valentini 0000-0002-6756-5634
Department for Innovation in Biological, Agri-Food and Forest Systems (DIBAF), University of Tuscia, 01100 Viterbo (Italy)
Monia Santini 0000-0002-8041-8241
Riccardo Valentini 0000-0002-6756-5634
CMCC Foundation - Euro-Mediterranean Center on Climate Change, Italy, 01100 Viterbo (Italy)
Matteo Guasti 0000-0002-0352-0760
CREA, Research Centre for Forestry and Wood, v.le Santa Margherita 80, IT-52100 Arezzo (Italy)
Forest Research Institute Baden-Württemberg, 79100, Freiburg (Germany)
Chair of Forest Growth and Dendroecology, University of Freiburg, Tennenbacherstaße 4, 79106 Freiburg (Germany)
Department of Forest Development, National School of Forest Engineers, 11000 Salé (Morocco)
Corresponding author
Paper Info
Citation
Boukhris I, Puletti N, Vonderach C, Guasti M, Lahssini S, Santini M, Valentini R (2024). Comparative analysis of taper models for Pinus nigra Arn. using terrestrial laser scanner acquired data. iForest 17: 203-212. - doi: 10.3832/ifor4525-017
Academic Editor
Angelo Nolè
Paper history
Received: Nov 17, 2023
Accepted: Jun 25, 2024
First online: Jul 22, 2024
Publication Date: Aug 31, 2024
Publication Time: 0.90 months
Copyright Information
© SISEF - The Italian Society of Silviculture and Forest Ecology 2024
Open Access
This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Web Metrics
Breakdown by View Type
Article Usage
Total Article Views: 178
(from publication date up to now)
Breakdown by View Type
HTML Page Views: 0
Abstract Page Views: 0
PDF Downloads: 74
Citation/Reference Downloads: 0
XML Downloads: 104
Web Metrics
Days since publication: 141
Overall contacts: 178
Avg. contacts per week: 8.84
Article Citations
Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Nov 2020)
(No citations were found up to date. Please come back later)
Publication Metrics
by Dimensions ©
Articles citing this article
List of the papers citing this article based on CrossRef Cited-by.
Related Contents
iForest Similar Articles
Research Articles
Modeling compatible taper and stem volume of pure Scots pine stands in Northeastern Turkey
vol. 16, pp. 38-46 (online: 22 January 2023)
Research Articles
Modelling taper and stem volume considering stand density in Eucalyptus grandis and Eucalyptus dunnii
vol. 14, pp. 127-136 (online: 16 March 2021)
Research Articles
Compatible taper-volume models of Quercus variabilis Blume forests in north China
vol. 10, pp. 567-575 (online: 08 May 2017)
Research Articles
Nonlinear mixed model approaches to estimating merchantable bole volume for Pinus occidentalis
vol. 5, pp. 247-254 (online: 24 October 2012)
Technical Reports
GlobAllomeTree: international platform for tree allometric equations to support volume, biomass and carbon assessment
vol. 6, pp. 326-330 (online: 18 July 2013)
Research Articles
Allometric equations to estimate above-ground biomass of small-diameter mixed tree species in secondary tropical forests
vol. 13, pp. 165-174 (online: 02 May 2020)
Research Articles
Allometric equations to assess biomass, carbon and nitrogen content of black pine and red pine trees in southern Korea
vol. 10, pp. 483-490 (online: 12 April 2017)
Research Articles
Equations for estimating belowground biomass of Silver Birch, Oak and Scots Pine in Germany
vol. 12, pp. 166-172 (online: 15 March 2019)
Technical Reports
Biomass equations for European beech growing on dry sites
vol. 9, pp. 751-757 (online: 17 June 2016)
Research Articles
Alternative methods of scaling Eucalyptus urophylla trees in forest stands: compatibility and accuracy of volume equations
vol. 11, pp. 275-283 (online: 29 March 2018)
iForest Database Search
Google Scholar Search
Citing Articles
Search By Author
Search By Keywords
PubMed Search
Search By Author
Search By Keyword