Estimation of forest structural and compositional variables using ALS data and multi-seasonal satellite imagery

Advanced forest resource inventory (FRI) information is of critical importance for sustainable forest management. FRIs are dependent on remote sensing data and processing methods, along with field calibration/validation to generate cost-effective options for modelling forest inventory and biophysical variables over large areas. The objective of this study was to examine the impact of combining multi-seasonal multispectral satellite imagery with airborne laser scanning (ALS) data for estimating basal area, species mixture and stem density for an uneven-aged tolerant hardwood forest in Ontario, Canada. Using random forest (RF) regression as a non-parametric diagnostic technique, three multispectral optical sensors (i.e., Landsat-5 TM, Sentinel-2 A and WorldView-2) were compared to examine the most cost-effective sensor configuration for modelling FRI variables. The contribution of spectral predictors derived from these optical sensors as well as ALS height and intensity metrics were evaluated using RF variable importance. As part of our variable selection framework, all predictor variables were grouped into relatively independent clusters using a hierarchical variable clustering technique, which revealed the distinctiveness between information contained in spectral predictors, height- and intensity-based metrics. This indicates that ALS intensity data carry unique information complementary to passive near-infrared data for forest characterization. ALS data alone did not result in accurate models for basal area and species mixture, but predictive accuracies were improved significantly with the addition of spectral predictors. Compared to single-date images, multi-seasonal imagery proved to be more accurate for modelling FRI variables, especially when combined with ALS data. Despite its limited spatial resolution, Sentinel-2 A was found to be the most cost-effective image source for enhancing ALS-based FRI models. Using variables identified by the variable selection procedure, best subsets regression outperformed the RF models developed for diagnostic analysis, resulting in a suite of accurate and parsimonious predictive models, with coefficients of determination of 0.73, 0.90 and 0.67, for basal area, species mixture, and stem density, respectively.Financial support for this research was provided by Natural Resources Canada (NRCan) through an ecoENERGY Innovative Initiative Grant [BIO 063] to the authors and a Natural Sciences and Engineering Research Council Discovery Grant to Treitz [RGPIN/03822]. The authors are also grateful to the Haliburton Forest, Tembec Inc., MITACS and Ontario Power Generation for their financial and logistical support of various components of this research