Underestimation of LiDAR heights is widely known but has never been evaluated for several sensors and for diverse types of ecological conditions. This underestimation is mainly linked to the probability of the pulse to reach the ground and the top of vegetation. Main causes of this underestimation are pulse density, pattern of scan (sensors), scan angles, specific contract parameters (flying altitude, pulse repetition frequency) and characteristics of the territory (slopes, stand density and species composition). This study, carried out at a resolution of 1 × 1 m, first assessed the possibility of making an adjustment model to correct the bias of the digital terrain model (DTM), and then proposed a global adjustment model to correct the bias on the canopy height model (CHM). For this study, the bias of both DTM and CHM were calculated by subtracting two LiDAR datasets: high-density pixels with 21 pulses/m² (first return) and more (DTM or CHM reference value pixels) and low-density pixels (DTM or CHM value to correct). After preliminary analyses, it was concluded that the DTM did not need specific adjustment. In contrast, the CHM needed adjustments. Among the variables studied, three were selected for the final CHM adjustment model: the maximum height of the pixel (H2Corr); the density of first returns by m² (D_first); and the standard deviation of nine maximum heights of the neighborhood cells (H_STD9). The modeling occurred in three steps. The first two steps enabled the determination of significant variables and the shape of the equation to be defined (linear mixed model and non-linear model). The third step made it possible to propose an empirical equation using a non-linear mixed model that can be applied to a 1 × 1 m CHM. The CHM underestimation correction could be used for a preliminary step to several uses of the CHM such as volume calculation, forest growth models or multi-temporal analysis.

So far, large-scale projects aiming to map forest attributes using aerial LiDAR data have been developed using ground sample plots acquired synchronously with LiDAR. No large projects have been developed using aerial LiDAR acquired independent of ground sample plot datasets. The goal of this study was to develop and validate large-scale parametric merchantable wood volume estimation models using existing historical ground sample plots. The models can be applied to large LiDAR datasets to map merchantable wood volume as a 10 × 10 m raster. The study demonstrated that a relative density index (RDI) based on a self-thinning equation and dominant height were suitable variables that can be calculated both for ground sample plots and LiDAR datasets. The resulting volume raster showed sound accuracy rates when compared to validation zones: R², 82.25%; RMSE, 13.7 m³/ha; and bias, −4.09 m³/ha. The results show that ground sample plot datasets acquired synchronously with LiDAR can be used to calculate the RDI and dominant height. These variables can consequently be used to map forest attributes over a large area with a high level of accuracy, thus not requiring the implementation of new costly sample plots.

A shadow fraction method was developed previously for mapping forest attributes of northern black spruce ( Picea mariana (Mill.) Britton, Sterns & Poggenb.) forests. This paper evaluates application of the method for (i) balsam fir stands ( Abies balsamea (L.) Mill.), (ii) stands with higher volume and biomass than those of previous studies, and (iii) stands with a higher composition of deciduous trees and steeper slopes. Models developed for new test sites in (i) central Labrador and (ii) western Newfoundland were not statistically different from previous models for biomass, volume, and basal area. Relative root mean square errors (RMSE_r) for central Labrador were slightly lower than those found in other test sites (RMSE_r: 24%–29%) but higher for western Newfoundland (RMSE_r = 37%–43%), attributed to the higher upper limit of measured attributes and increased presence of deciduous trees. Results suggest that reasonable estimates can be generated for conifer forests of northeastern Canada; however, an alternative solution is needed where mixed and deciduous stands are prevalent. Measurement of ground plots over a wider range of species composition and forest structure is recommended for broader application to northern boreal forests and to further assess the potential role of the shadow fraction method in national-scale inventory programs.

Monitoring crown closure evolution using multi-temporal Light Detection and Ranging (LiDAR) surveys is a method that we expect to be increasingly adopted given the availability of LiDAR sensors and the accumulating survey archives. However, little attention was devoted to comparing crown closure estimates from independent surveys. Although survey parameters cannot be modified after the data collection, we speculate that the error associated to crown closure estimates comparison can be reduced by selecting optimal post-survey parameters. In this study, we compared crown closure estimates of three airborne LiDAR surveys from 2018 (40 pt/m²) used as a reference, and two lower-density surveys from 2016 (4.5 pt/m²) and 2018 (2 pt/m²). We studied the effect of the height threshold used to separate canopy points and the grid resolution, using skewness and variance of lagged difference of crown closure. Crown closure estimates using low height thresholds were more different across surveys, resulting in higher root mean squared error (RMSE), bias and more different variograms. Results show that optimal height threshold was 3 m and grid resolution was 25 m, although there was room for decision (RMSE of 7% and 5%, and bias of 4% and 0% for 2016 and 2018 low-density surveys).