Search

Sustainable forest management requires information on the spatial distribution, composition, and structure of forests. However, jurisdictions with large tracts of noncommercial forest, such as the Northwest Territories (NWT) of Canada, often lack detailed forest information across their land base. The goal of the Multisource Vegetation Inventory (MVI) project was to create a large area forest inventory (FI) map that could support strategic forest management in the NWT using optical, radar, and light detection and ranging (LiDAR) satellite remote sensing anchored on limited field plots and airborne LiDAR data. A new landcover map based on Landsat imagery was the first step to stratify forestland into broad forest types. A modelling chain linking FI plots to airborne and spaceborne LiDAR was then developed to circumvent the scarcity of field data in the region. The developed models allowed the estimation of forest attributes in thousands of surrogate FI plots corresponding to spaceborne LiDAR footprints distributed across the project area. The surrogate plots were used as a reference dataset for estimating each forest attribute in each 30 m forest cell within the project area. The estimation was based on the k-nearest neighbour (k-NN) algorithm, where the selection of the four most similar surrogate FI plots to each cell was based on satellite, topographic, and climatic data. Wall-to-wall 30 m raster maps of broad forest type, stand height, crown closure, stand volume, total volume, aboveground biomass, and stand age were created for a ~400,000 km2 area, validated with independent data, and generalized into a polygon GIS layer resembling a traditional FI map. The MVI project showed that a reasonably accurate FI map for large, remote, predominantly non-inventoried boreal regions can be obtained at a low cost by combining limited field data with remote sensing data from multiple sources.

Satellite forest inventories are the only feasible way to map Canada’s vast, remote forest regions, such as those in the Northwest Territories (NWT). A method used to create such inventories is the k-nearest neighbour (k-NN) algorithm, which spatially extends information from forest inventory (FI) plots to the entire forest land base using wall-to-wall features typically derived from Landsat data. However, the benefits of integrating L-band synthetic aperture radar (SAR) data, strongly correlated to forest biomass, have not been assessed for Canadian northern boreal forests. Here we describe an optimized multivariate k-NN implementation of a 151,700 km2 area in southern NWT that included ca. 2007 Landsat and dual-polarized Phased Array type L-band SAR (PALSAR) data on board the Advanced Land Observing Satellite (ALOS). Five forest attributes were mapped at 30 m cells: stand height, crown closure, stand/total volume and aboveground biomass (AGB). We assessed accuracy gains compared to Landsat-based maps. To circumvent the scarcity of FI plots, we used 3600 footprints from the Geoscience Laser Altimeter System (GLAS) as surrogate FI plots, where forest attributes were estimated using Light Detection and Ranging (LiDAR) metrics as predictors. After optimization, k-NN predicted forest attribute values for each pixel as the average of the 4 nearest (k = 4) surrogate FI plots within the Euclidian space of 9 best features (selected among 6 PALSAR, 10 Landsat, and 6 environmental features). Accuracy comparisons were based on 31 National Forest Inventory ground plots and over 1 million airborne LiDAR plots. Maps that included PALSAR HV backscatter resulted in forest attribute predictions with higher goodness of fit (adj. R2), lower percent mean error (ME%), and percent root mean square error (RMSE%), and lower underestimation for larger attribute values. Predictions were most accurate for conifer stand height (RMSE% = 32.1%, adj. R2 = 0.58) and AGB (RMSE% = 47.8%, adj. R2 = 0.74), which is much more abundant in the area than mixedwood or broadleaf. Our study demonstrates that optimizing k-NN parameters and feature space, including PALSAR, Landsat, and environmental variables, is a viable approach for inventory mapping of the northern boreal forest regions of Canada.

A study was undertaken to assess Envisat Advanced Synthetic Aperture Radar (ASAR) ScanSAR data for quantifying forest growing stock volume (GSV) across three boreal regions with varying forest types, composition, and structure (Sweden, Central Siberia, and Québec). Estimates of GSV were obtained using hyper-temporal observations of the radar backscatter acquired by Envisat ASAR with the BIOMASAR algorithm. In total, 5.3×106 km2 were mapped with a 0.01° pixel size to obtain estimates representative for the year of 2005. Comparing the SAR-based estimates to spatially explicit datasets of GSV, generated from forest field inventory and/or Earth Observation data, revealed similar spatial distributions of GSV. Nonetheless, the weak sensitivity of C-band backscatter to forest structural parameters introduced significant uncertainty to the estimated GSV at full resolution. Further discrepancies were observed in the case of different scales of the ASAR and the reference GSV and in areas of fragmented landscapes. Aggregation to 0.1° and 0.5° was then undertaken to generate coarse scale estimates of GSV. The agreement between ASAR and the reference GSV datasets improved; the relative difference at 0.5° was consistently within a magnitude of 20–30%. The results indicate an improvement of the characterization of forest GSV in the boreal zone with respect to currently available information.

A methods framework is presented that utilizes field plots, airborne light detection and ranging (LiDAR), and spaceborne Geoscience Laser Altimeter System (GLAS) data to estimate forest attributes over a 20 Mha area in Northern Canada. The framework was implemented to scale up forest attribute models from field data to intersecting airborne LiDAR data, and then to GLAS footprints. GLAS data were sequentially filtered and submitted to the k-nearest neighbour (k-NN) imputation algorithm to yield regional estimates of stand height and crown closure at a 30 m resolution. Resulting outputs were assessed against independent airborne LiDAR data to evaluate regional estimates of stand height (mean difference = −1 m, RMSE = 5 m) and crown closure (mean difference = −5%, RMSE = 9%). Additional assessments were performed as a function of dominant vegetation type and ecoregion to further evaluate regional products. These attributes form the primary descriptive structure attributes that are typical of forest inventory mapping programs, and provide insight into how they can be derived in northern boreal regions where field information and physical access is often limited.