Low-density Airborne Laser Scanning Research Articles

Benchmarking of Individual Tree Segmentation Methods in Mediterranean Forest Based on Point Clouds from Unmanned Aerial Vehicle Imagery and Low-Density Airborne Laser Scanning

Three raster-based (RB) and one point cloud-based (PCB) algorithms were tested to segment individual Aleppo pine trees and extract their tree height (H) and crown diameter (CD) using two types of point clouds generated from two different techniques: (1) Low-Density (≈1.5 points/m2) Airborne Laser Scanning (LD-ALS) and (2) photogrammetry based on high-resolution unmanned aerial vehicle (UAV) images. Through intensive experiments, it was concluded that the tested RB algorithms performed best in the case of UAV point clouds (F1-score > 80.57%, H Pearson’s r > 0.97, and CD Pearson´s r > 0.73), while the PCB algorithm yielded the best results when working with LD-ALS point clouds (F1-score = 89.51%, H Pearson´s r = 0.94, and CD Pearson´s r = 0.57). The best set of algorithm parameters was applied to all plots, i.e., it was not optimized for each plot, in order to develop an automatic pipeline for mapping large areas of Mediterranean forests. In this case, tree detection and height estimation showed good results for both UAV and LD-ALS (F1-score > 85% and >76%, and H Pearson´s r > 0.96 and >0.93, respectively). However, very poor results were found when estimating crown diameter (CD Pearson´s r around 0.20 for both approaches).

Quantifying changes in forest structural complexity using bi-temporal airborne laser scanning

ABSTRACT Structural complexity is an important forest characteristic for habitat assessment, forest management and planning. However, monitoring how forest structural complexity evolves over time in various forest types has not been widely explored. In this study we investigate the feasibility of bi-temporal low-density airborne laser scanning (ALS) for the assessment of changes in light availability conditions within forest canopy, considered to imply changes in forest structural complexity. We used ALS data acquired in 2012 and 2019 to generate canopy vertical profiles by slicing the point clouds into 4 × 4 ×1 m voxels which were then rasterized and reclassified into four light availability categories. To understand structural development over time in different forest types we used field-measured tree heights and tree species information to stratify sample plots in different stand complexity categories. Stands with higher structural complexity represented increased proportions of space occupied by vegetation as well as decreased proportions of empty space below the canopy. The experiments showed the ability of low-density ALS to characterize the dynamics in canopy layering structure, implying changes in forest structural complexity. The presented methodology could potentially be upscaled and applied in the landscape-level monitoring of the development of boreal forest structural characteristics.

Improvements of Fire Fuels Attributes Maps by Integrating Field Inventories, Low Density ALS, and Satellite Data in Complex Mediterranean Forests

One of the most determining factors in forest fire behaviour is to characterize forest fuel attributes. We investigated a complex Mediterranean forest type—mountainous Abies pinsapo–Pinus–Quercus–Juniperus with distinct structures, such as broadleaf and needleleaf forests—to integrate field data, low density Airborne Laser Scanning (ALS), and multispectral satellite data for estimating forest fuel attributes. The three-step procedure consisted of: (i) estimating three key forest fuel attributes (biomass, structural complexity and hygroscopicity), (ii) proposing a synthetic index that encompasses the three attributes to quantify the potential capacity for fire propagation, and (iii) generating a cartograph of potential propagation capacity. Our main findings showed that Biomass–ALS calibration models performed well for Abies pinsapo (R2 = 0.69), Juniperus spp. (R2 = 0.70), Pinus halepensis (R2 = 0.59), Pinus spp. mixed (R2 = 0.80), and Pinus spp.–Juniperus spp. (R2 = 0.59) forests. The highest values of biomass were obtained for Pinus halepensis forests (190.43 Mg ha−1). The structural complexity of forest fuels was assessed by calculating the LiDAR Height Diversity Index (LHDI) with regard to the distribution and vertical diversity of the vegetation with the highest values of LHDI, which corresponded to Pinus spp.–evergreen (2.56), Quercus suber (2.54), and Pinus mixed (2.49) forests, with the minimum being obtained for Juniperus (1.37) and shrubs (1.11). High values of the Fuel Desiccation Index (IDM) were obtained for those areas dominated by shrubs (−396.71). Potential Behaviour Biomass Index (ICB) values were high or very high for 11.86% of the area and low or very low for 77.07%. The Potential Behaviour Structural Complexity Index (ICE) was high or very high for 37.23% of the area, and low or very low for 46.35%, and the Potential Behaviour Fuel Desiccation Index (ICD) was opposite to the ICB and ICE, with high or very high values for areas with low biomass and low structural complexity. Potential Fire Behaviour Index (ICP) values were high or very high for 38.25% of the area, and low or very low values for 45.96%. High or very high values of ICP were related to Pinus halepensis and Pinus pinaster forests. Remote sensing has been applied to improve fuel attribute characterisation and cartography, highlighting the utility of integrating multispectral and ALS data to estimate those attributes that are more closely related to the spatial organisation of vegetation.

Identification of Old-Growth Mediterranean Forests Using Airborne Laser Scanning and Geostatistical Analysis

The protection and conservation of old-growth forests (OGFs) are becoming a global concern due to their irreplaceability and high biodiversity. Nonetheless, there has been little research into the identification and characterization of OGFs of the oldest tree species in Mediterranean areas. We used forest inventory data, low-density airborne laser scanning (ALS) metrics, and geostatistical analysis to estimate old-growth indices (OGIs) as indicators of old-growth forest conditions. We selected a pilot area in European black pine (Pinus nigra subsp. salzmannii) ecosystems where the oldest known living trees in the Iberian Peninsula are found. A total of 756 inventory plots were established to characterize standard live tree and stand attributes. We estimated several structural attributes that discriminate old growth from younger age classes and calculated different types of OGI for each plot. The best OGI was based on mean tree diameter, standard deviation of tree diameter, and stand density of large trees (diameter > 50 cm). This index is useful for assessing old-growthness at different successional stages (young and OGFs) in Mediterranean black pine forests. Our results confirm that the estimation of OGIs based on a combination of forest inventory data, geostatistical analysis, and ALS is useful for identifying OGFs.

Estimating Stand and Fire-Related Surface and Canopy Fuel Variables in Pine Stands Using Low-Density Airborne and Single-Scan Terrestrial Laser Scanning Data

In this study, we used data from a thinning trial conducted on 34 different sites and 102 sample plots established in pure and even-aged Pinus radiata and Pinus pinaster stands, to test the potential use of low-density airborne laser scanning (ALS) metrics and terrestrial laser scanning (TLS) metrics to provide accurate estimates of variables related to surface and canopy fires. An exhaustive field inventory was carried out in each plot to estimate the main stand variables and the main variables related to fire hazard: surface fuel loads by layers, fuel strata gap, surface fuel height, stand mean height, canopy base height, canopy fuel load and canopy bulk density. In addition, the point clouds from low-density ALS and single-scan TLS of each sample plot were used to calculate metrics related to the vertical and horizontal distribution of forest fuels. The comparative performance of the following three non-parametric machine learning techniques used to estimate the main stand- and fire-related variables from those metrics was evaluated: (i) multivariate adaptive regression splines (MARS), (ii) support vector machine (SVM), and (iii) random forest (RF). The selection of the best modeling approach was based on a comparison of the root mean square error (RMSE), obtained by optimizing the parameters of each technique and performing cross-validation. Overall, the best results were obtained with the MARS techniques for data from both sensors. The TLS data provided the best results for variables associated with the internal characteristics of canopy structure and understory fuel but were less reliable for estimating variables associated with the upper canopy, due to occlusion by mid-canopy foliage. The combination of ALS and TLS metrics improved the accuracy of estimates for all variables analyzed, except the height and the biomass of the understory shrubs. The variability demonstrated by the combined use of both types of metrics ranged from 43.11% for the biomass of duff litter layers to 94.25% for dominant height. The results suggest that the combination of machine learning techniques and metrics derived from low-density ALS data, drawn from a single-scan TLS or a combination of both metrics, may represent a promising alternative to traditional field inventories for obtaining valuable information about surface and canopy fuel variables at large scales.

Modeling Diameter Distributions with Six Probability Density Functions in Pinus halepensis Mill. Plantations Using Low-Density Airborne Laser Scanning Data in Aragón (Northeast Spain)

The diameter distributions of trees in 50 temporary sample plots (TSPs) established in Pinus halepensis Mill. stands were recovered from LiDAR metrics by using six probability density functions (PDFs): the Weibull (2P and 3P), Johnson’s SB, beta, generalized beta and gamma-2P functions. The parameters were recovered from the first and the second moments of the distributions (mean and variance, respectively) by using parameter recovery models (PRM). Linear models were used to predict both moments from LiDAR data. In recovering the functions, the location parameters of the distributions were predetermined as the minimum diameter inventoried, and scale parameters were established as the maximum diameters predicted from LiDAR metrics. The Kolmogorov–Smirnov (KS) statistic (Dn), number of acceptances by the KS test, the Cramér von Misses (W2) statistic, bias and mean square error (MSE) were used to evaluate the goodness of fits. The fits for the six recovered functions were compared with the fits to all measured data from 58 TSPs (LiDAR metrics could only be extracted from 50 of the plots). In the fitting phase, the location parameters were fixed at a suitable value determined according to the forestry literature (0.75·dmin). The linear models used to recover the two moments of the distributions and the maximum diameters determined from LiDAR data were accurate, with R2 values of 0.750, 0.724 and 0.873 for dg, dmed and dmax. Reasonable results were obtained with all six recovered functions. The goodness-of-fit statistics indicated that the beta function was the most accurate, followed by the generalized beta function. The Weibull-3P function provided the poorest fits and the Weibull-2P and Johnson’s SB also yielded poor fits to the data.

Assessing the Potential of the DART Model to Discrete Return LiDAR Simulation—Application to Fuel Type Mapping

Fuel type is one of the key factors for analyzing the potential of fire ignition and propagation in agricultural and forest environments. The increase of three-dimensional datasets provided by active sensors, such as LiDAR (Light Detection and Ranging), has improved the classification of fuel types through empirical modelling. Empirical methods are site and sensor specific while Radiative Transfer Models (RTM) approaches provide broader universality. The aim of this work is to analyze the suitability of Discrete Anisotropic Radiative Transfer (DART) model to replicate low density small-footprint Airborne Laser Scanning (ALS) measurements and subsequent fuel type classification. Field data measured in 104 plots are used as ground truth to simulate LiDAR response based on the sensor and flight characteristics of low-density ALS data captured by the Spanish National Plan for Aerial Orthophotography (PNOA) in two different dates (2011 and 2016). The accuracy assessment of the DART simulations is performed using Spearman rank correlation coefficients between the simulated metrics and the ALS-PNOA ones. The results show that 32% of the computed metrics overpassed a correlation value of 0.80 between simulated and ALS-PNOA metrics in 2011 and 28% in 2016. The highest correlations were related to high height percentiles, canopy variability metrics as for example standard deviation and Rumple diversity index, reaching correlation values over 0.94. Two metric selection approaches and Support Vector Machine classification method with variants were compared to classify fuel types. The best-fitted classification model, trained with the DART simulated sample and validated with ALS-PNOA data, was obtained using Support Vector Machine method with radial kernel. The overall accuracy of the classification after validation was 88% and 91% for the 2011 and 2016 years, respectively. The use of DART demonstrates its value for simulating generalizable 3D data for fuel type classification providing relevant information for forest managers in fire prevention and extinction.

Developing a site index model for P. Pinaster stands in NW Spain by combining bi-temporal ALS data and environmental data

Site index (SI) is a common measure of forest site productivity, serving as a valuable baseline for forest management. The main objective of this study was to develop a SI model for Pinus pinaster Ait. in north-west Spain by combining bi–temporal, low–density airborne laser scanning (ALS) data (acquired in the periods 2009–2011 and 2015–2017) with climatic, edaphic and physiographical data. Site productivity, assessed by site quality curves, was modelled using an age-independent difference equation method based on ALS metrics and environmental variables. For the model development process, we used data from 156 sample plots in pure and even-aged P. pinaster stands distributed throughout Galicia (NW Spain) and measured in the Spanish National Forest Inventory (SNFI). The generalized algebraic difference approach (GADA) formulation was tested by using two different base equations for modelling the dominant height growth (ΔH) from ALS variables. The GADA formulation derived from the Bertalanffy’s base model produced the best estimates of dominant height (H) for P. pinaster stands in Galicia. Use of the proposed model to estimate ΔH for a new pine stand requires two ALS data sets for estimating site-specific (local) parameters. To enable use of the model when such information is not available, the relationship between the values of the site-specific parameter and environmental variables was described using Multivariate Adaptive Regression Splines (MARS). Use of the MARS equation enabled us to develop spatially-explicit predictive maps of the site-specific parameter values, which can be used together with the GADA model to derive ΔH curves and SI estimates for P.pinaster stands in the whole study region.

Forest structural diversity characterization in Mediterranean landscapes affected by fires using Airborne Laser Scanning data

ABSTRACT Forest fires can change forest structure and composition, and low-density Airborne Laser Scanning (ALS) can be a valuable tool for evaluating post-fire vegetation response. The aim of this study is to analyze the structural diversity differences in Mediterranean Pinus halepensis Mill. forests affected by wildfires on different dates from 1986 to 2009. Several types of ALS metrics, such as the Light Detection and Ranging (LiDAR) Height Diversity Index (LHDI), the LiDAR Height Evenness Index (LHEI), and vertical and horizontal continuity of vegetation, as well as topographic metrics, were obtained in raster format from low point density data. In order to map burned and unburned areas, differentiate fire occurrence dates, and distinguish between old and more recent fires, a sample of pixels was previously selected to assess the existence of differences in forest structure using the Kruskal–Wallis test. Then, k-nearest neighbors algorithm (k-NN), support vector machine (SVM) and random forest (RF) classifiers were compared to select the most accurate technique. The results showed that, in more recent fires, around 70% of the laser returns came from grass and shrub layers, yielding low LHDI and LHEI values (0.37–0.65 and 0.28–0.46, respectively). In contrast, the areas burned more than 20 years ago had higher LHDI and LHEI values due to the growth of the shrub and tree strata. The classification of burned and unburned areas yielded an overall accuracy of 89.64% using the RF method. SVM was the best classifier for identifying the structural differences between fires occurring on different dates, with an overall accuracy of 68.79%. Furthermore, SVM yielded an overall accuracy of 75.49% for the classification between old and more recent fires.

Assessment of the Carbon Stock in Pine Plantations in Southern Spain through ALS Data and K-Nearest Neighbor Algorithm Based Models

Accurate estimation of forest biomass to enable the mapping of forest C stocks over large areas is of considerable interest nowadays. Airborne laser scanning (ALS) systems bring a new perspective to forest inventories and subsequent biomass estimation. The objective of this research was to combine growth models used to update old inventory data to a reference year, low-density ALS data, and k-nearest neighbor (kNN) algorithm Random Forest to conduct biomass inventories aimed at estimating the C sequestration capacity in large Pinus plantations. We obtained a C stock in biomass (Wt-S) of 12.57 Mg·ha−1, ranging significantly from 19.93 Mg·ha−1 for P. halepensis to 49.05 Mg·ha−1 for P. nigra, and a soil organic C stock of the composite soil samples (0–40 cm) ranging from 20.41 Mg·ha−1 in P. sylvestris to 37.32 Mg·ha−1 in P. halepensis. When generalizing these data to the whole area, we obtained an overall C-stock value of 48.01 MgC·ha−1, ranging from 23.96 MgC·ha−1 for P. halepensis to 58.09 MgC·ha−1 for P. nigra. Considering the mean value of the on-site C stock, the study area sustains 1,289,604 Mg per hectare (corresponding to 4,732,869 Mg CO2), with a net increase of 4.79 Mg·ha−1·year−1. Such C cartography can help forest managers to improve forest silviculture with regard to C sequestration and, thus, climate change mitigation.

Challenges of Multi-Temporal and Multi-Sensor Forest Growth Analyses in a Highly Disturbed Boreal Mixedwood Forests

Airborne laser scanning (ALS) and digital aerial photogrammetry (DAP) have both been demonstrated as reliable sources of information on forest stand inventory attributes. The increasing availability of both datasets provides a means for improving stand dynamics information over time; however, the cost of multi-temporal ALS can be prohibitive in some circumstances. As a result, a combination of ALS at an initial time step and subsequent updates using DAP has been proposed as a cost-effective alternative for maintaining forest inventories. In this study we used low density ALS and DAP point clouds acquired in 2007 and 2015, respectively, to quantify changes in forest structure, in a highly disturbed boreal mixedwood forest in Alberta, Canada. We examined the capacity of the two technologies to model changes in top height (H), volume (V), and basal area (BA) using both direct and indirect approaches for estimation. Results indicate that the proportion of explained variance (adjusted R2) for the models derived from the ALS (Time 1; T1) and DAP (Time 2; T2) data were highest for models predicting H at T1, and lowest for BA at T1 and T2 (R2 was 0.66–0.70). The indirect estimates of change in H, BA, and V were calculated by subtracting the T1 and T2 predictions. For the direct approach, separate regression models were developed that used the differences in point cloud metrics between T1 and T2 as predictors. Results indicated that the accuracy of the estimates generated using the indirect approach were markedly lower than the estimates generated using the direct approach, with especially poor results for ∆BA and ∆V. Best results were achieved for ∆H using the direct approach with an R2 coefficient of 0.65 and an root mean square error (RMSE)% of 190.06%. We found that the error associated with change estimates of H, BA and V increased with the increase or decrease in mortality. We conclude that forest managers should act carefully when applying the multi-temporal and multi-sensor analysis of forest growth if forest growth is slow and the level of mortality is high.

Temporal Transferability of Pine Forest Attributes Modeling Using Low-Density Airborne Laser Scanning Data

This study assesses model temporal transferability using airborne laser scanning (ALS) data acquired over two different dates. Seven forest attributes (i.e. stand density, basal area, squared mean diameter, dominant diameter, tree dominant height, timber volume, and total tree biomass) were estimated using an area-based approach in Mediterranean Aleppo pine forests. Low-density ALS data were acquired in 2011 and 2016 while 147 forest inventory plots were measured in 2013, 2014, and 2016. Single-tree growth models were used to generate concomitant field data for 2011 and 2016. A comparison of five selection techniques and five regression methods were performed to regress field observations against ALS metrics. The selection of the best regression models fitted for each stand attribute, and separately for both 2011 and 2016, was performed following an indirect approach. Model performance and temporal transferability were analyzed by extrapolating the best fitted models from 2011 to 2016 and inversely from 2016 to 2011 using the direct approach. Non-parametric support vector machine with radial kernel was the best regression method with average relative % root mean square error differences of 2.13% for 2011 models and 1.58% for 2016 ones.

Airborne Laser Scanning Cartography of On-Site Carbon Stocks as a Basis for the Silviculture of Pinus Halepensis Plantations

Forest managers are interested in forest-monitoring strategies using low density Airborne Laser Scanning (ALS). However, little research has used ALS to estimate soil organic carbon (SOC) as a criterion for operational thinning. Our objective was to compare three different thinning intensities in terms of the on-site C stock after 13 years (2004–2017) and to develop models of biomass (Wt, Mg ha−1) and SOC (Mg ha−1) in Pinus halepensis forest, based on low density ALS in southern Spain. ALS was performed for the area and stand metrics were measured within 83 plots. Non-parametric kNN models were developed to estimate Wt and SOC. The overall C stock was significantly higher in plots subjected to heavy or moderate thinning (101.17 Mg ha−1 and 100.94 Mg ha−1, respectively) than in the control plots (91.83 Mg ha−1). The best Wt and SOC models provided R2 values of 0.82 (Wt, MSNPP) and 0.82 (SOC-S10, RAW). The study area will be able to stock 134,850 Mg of C under a non-intervention scenario and 157,958 Mg of C under the heavy thinning scenario. High-resolution cartography of the predicted C stock is useful for silvicultural planning and may be used for proper management to increase C sequestration in dry P. halepensis forests.

Estimation of Total Biomass in Aleppo Pine Forest Stands Applying Parametric and Nonparametric Methods to Low-Density Airborne Laser Scanning Data

The account of total biomass can assist with the evaluation of climate regulation policies from local to global scales. This study estimates total biomass (TB), including tree and shrub biomass fractions, in Pinus halepensis Miller forest stands located in the Aragon Region (Spain) using Airborne Laser Scanning (ALS) data and fieldwork. A comparison of five selection methods and five regression models was performed to relate the TB, estimated in 83 field plots through allometric equations, to several independent variables extracted from ALS point cloud. A height threshold was used to include returns above 0.2 m when calculating ALS variables. The sample was divided into training and test sets composed of 62 and 21 plots, respectively. The model with the lower root mean square error (15.14 tons/ha) after validation was the multiple linear regression model including three ALS variables: the 25th percentile of the return heights, the variance, and the percentage of first returns above the mean. This study confirms the usefulness of low-density ALS data to accurately estimate total biomass, and thus better assess the availability of biomass and carbon content in a Mediterranean Aleppo pine forest.

Modeling diameter distributions in radiata pine plantations in Spain with existing countrywide LiDAR data

We evaluated the use of low-density airborne laser scanning data to estimate diameter distributions in radiata pine plantations. The moment-based parameter recovery method was used to estimate the diameter distributions, considering LiDAR metrics as explanatory variables. The fitted models explained more than 77% of the observed variability. This approach can be replicated every 6 years (temporal cover planned for countrywide LiDAR flights in Spain). The estimation of stand diameter distribution is informative for forest managers in terms of stand structure, forest growth model inputs, and economic timber value. In this sense, airborne LiDAR may represent an adequate source of information. The objective was to evaluate the use of low-density Spanish countrywide LiDAR data for estimating diameter distributions in Pinus radiata D. Don stands in NW Spain. The empirical distributions were obtained from 25 sample plots. We applied the moment-based parameter recovery method in combination with the Weibull function to estimate the diameter distributions, considering LiDAR metrics as explanatory variables. We evaluated the results by using the Kolmogorov–Smirnov (KS) test and a classification tree and random forest (RF) to relate the KS test result for each plot to stand-level variables. The models used to estimate average (dm) and quadratic (dg) mean diameters from LiDAR metrics, required for recovery of the Weibull parameters, explained a high percentage of the observed variance (77 and 80%, respectively), with RMSE values of 3.626 and 3.422 cm for the same variables. However, the proportion of plots accepted by the KS was low. This poor performance may be due to the strictness of the KS test and/or by the characteristics of the LiDAR flight. The results justify the assessment of this approach over different species and forest types in regional or countrywide surveys.

A NEW APPROACH FOR PROGRESSIVE DENSE RECONSTRUCTION FROM CONSECUTIVE IMAGES BASED ON PRIOR LOW-DENSITY 3D POINT CLOUDS

Abstract. In recent years, the increasing incidence of climate-related disasters has tremendously affected our environment. In order to effectively manage and reduce dramatic impacts of such events, the development of timely disaster management plans is essential. Since these disasters are spatial phenomena, timely provision of geospatial information is crucial for effective development of response and management plans. Due to inaccessibility of the affected areas and limited budget of first-responders, timely acquisition of the required geospatial data for these applications is usually possible only using low-cost imaging and georefencing sensors mounted on unmanned platforms. Despite rapid collection of the required data using these systems, available processing techniques are not yet capable of delivering geospatial information to responders and decision makers in a timely manner. To address this issue, this paper introduces a new technique for dense 3D reconstruction of the affected scenes which can deliver and improve the needed geospatial information incrementally. This approach is implemented based on prior 3D knowledge of the scene and employs computationally-efficient 2D triangulation, feature descriptor, feature matching and point verification techniques to optimize and speed up 3D dense scene reconstruction procedure. To verify the feasibility and computational efficiency of the proposed approach, an experiment using a set of consecutive images collected onboard a UAV platform and prior low-density airborne laser scanning over the same area is conducted and step by step results are provided. A comparative analysis of the proposed approach and an available image-based dense reconstruction technique is also conducted to prove the computational efficiency and competency of this technique for delivering geospatial information with pre-specified accuracy.

Modelling the vertical distribution of canopy fuel load using national forest inventory and low-density airbone laser scanning data.

The fuel complex variables canopy bulk density and canopy base height are often used to predict crown fire initiation and spread. Direct measurement of these variables is impractical, and they are usually estimated indirectly by modelling. Recent advances in predicting crown fire behaviour require accurate estimates of the complete vertical distribution of canopy fuels. The objectives of the present study were to model the vertical profile of available canopy fuel in pine stands by using data from the Spanish national forest inventory plus low-density airborne laser scanning (ALS) metrics. In a first step, the vertical distribution of the canopy fuel load was modelled using the Weibull probability density function. In a second step, two different systems of models were fitted to estimate the canopy variables defining the vertical distributions; the first system related these variables to stand variables obtained in a field inventory, and the second system related the canopy variables to airborne laser scanning metrics. The models of each system were fitted simultaneously to compensate the effects of the inherent cross-model correlation between the canopy variables. Heteroscedasticity was also analyzed, but no correction in the fitting process was necessary. The estimated canopy fuel load profiles from field variables explained 84% and 86% of the variation in canopy fuel load for maritime pine and radiata pine respectively; whereas the estimated canopy fuel load profiles from ALS metrics explained 52% and 49% of the variation for the same species. The proposed models can be used to assess the effectiveness of different forest management alternatives for reducing crown fire hazard.

Comparison of regression models to estimate biomass losses and CO2 emissions using low-density airborne laser scanning data in a burnt Aleppo pine forest

ABSTRACTThe knowledge of the forest biomass reduction produced by a wildfire can assist in the estimation of greenhouse gases to the atmosphere. This study focuses on the estimation of biomass losses and CO2 emissions by combustion of Aleppo pine forest in a wildfire occurred in the municipality of Luna (Spain). The availability of low point density airborne laser scanning (ALS) data allowed the estimation of pre-fire aboveground forest biomass. A comparison of nine regression models was performed in order to relate the biomass, estimated in 46 field plots, to several independent variables extracted from the ALS data. The multivariate linear regression selected model, including the percentage of first returns above 2 m and 40th percentile of the return heights, was validated using a leave-one-out cross-validation technique (6.1 ton/ha root mean square error). Biomass losses were estimated in a three-phase approach: (i) wildfire severity was obtained using the difference normalized burn ratio , (ii) Aleppo pine forest was delimited using the National Forest Map and ALS data and (iii) burning efficiency factors were applied considering severity levels. Post-fire biomass was then transformed into CO2 emissions (426,754.8 ton). This study evidences the usefulness of low-density ALS data to accurately estimate pre-fire biomass, in order to assess CO2 emissions in a Mediterranean Aleppo pine forest.

Estimating Tree Height Distribution Using Low-Density ALS Data With and Without Training Data

This study applies an approach based on stochastic geometry for retrieval of forest characteristics from airborne laser scanning (ALS) in two situations: 1) without ground-measured training data and 2) with training data. The applied model treats the ALS echo heights as an outcome of a random process, expressing the observed heights of canopy envelope as a function of stand density, the parameters of the tree height distribution, and the shape of the individual tree crown. The model was applied to a eucalyptus plantation dataset with known spacing, where the main interest was to estimate the plot-specific tree height distribution. Estimation without training data resulted in RMSEs of 2.9 and 0.9 m for mean and dominant heights, respectively. Estimation using training data resulted in RMSE's of 1.4 and 0.8 m, respectively. In both cases, the estimates of dominant height were more accurate than with the reference method, but the estimates of mean height were less accurate (area-based approach; RMSEs 1.1 and 0.9 m, respectively). The model-based method was robust to substantial decrease in echo density from 1.4 echoes/m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> to 0.14 echoes/m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> .

Retrieving Suppressed Trees from Model-Based Height Distribution by Combining High- and Low-Density Airborne Laser Scanning Data

The individual tree detection (ITD) using airborne laser scanning (ALS) data suffers from the problem of underestimating forest total attributes, caused by deficient detectability of suppressed trees. A new method is proposed in this study to retrieve height information for suppressed trees in order to improve the estimation of forest total stem volume at plot level. Area-based approach (ABA) and ITD were, respectively, applied to predict plot-level height distributions. Automatic detection of the cut points that distinguish dominant and suppressed trees was developed. Suppressed trees in the ITD-derived height distribution were retrieved from the corresponding height classes of the ABA-derived height distribution. The 2 basic and 5 calibrated height distributions were assessed using the Reynolds error index and stem volume estimates. ITD, although exposed to systematic underestimation, obtained better accuracy for the stem volume than ABA did. After calibration, the relative root mean square error (RMSE) for the estimated stem volume decreased from 18.66% to 16.70%, and the bias decreased from 10.39% to −0.41%, at best.