Search

Abstract Peatlands in the Canadian boreal forest are being negatively impacted by anthropogenic climate change, the effects of which are expected to worsen. Peatland types and sub‐classes vary in their ecohydrological characteristics and are expected to have different responses to climate change. Large‐scale modelling frameworks such as the Canadian Model for Peatlands, the Canadian Fire Behaviour Prediction System and the Canadian Land Data Assimilation System require peatland maps including information on sub‐types and vegetation as critical inputs. Additionally, peatland class and vegetation height are critical variables for wildlife habitat management and are related to the carbon cycle and wildfire fuel loading. This research aimed to create a map of peatland sub‐classes (bog, poor fen, rich fen permafrost peat complex) for the Canadian boreal forest and create an inventory of peatland vegetation height characteristics using ICESat‐2. A three‐stage hierarchical classification framework was developed to map peatland sub‐classes within the Canadian boreal forest circa 2020. Training and validation data consisted of peatland locations derived from various sources (field data, aerial photo interpretation, measurements documented in literature). A combination of multispectral data, L‐band SAR backscatter and C‐Band interferometric SAR coherence, forest structure and ancillary variables was used as model predictors. Ancillary data were used to mask agricultural areas and urban regions and account for regions that may exhibit permafrost. In the first stage of the classification, wetlands, uplands and water were classified with 86.5% accuracy. In the second stage, within the wetland areas only, peatland and mineral wetlands were differentiated with 93.3% accuracy. In the third stage, constrained to only the peatland areas, bogs, rich fens, poor fens and permafrost peat complexes were classified with 71.5% accuracy. Then, ICESat‐2 ATL08 spaceborne lidar data were used to describe regional variations in peatland vegetation height characteristics and regional and class‐wise variations based on a boreal forest wide sample. This research introduced a comprehensive large‐scale peatland sub‐class mapping framework for the Canadian boreal forest, presenting the first moderate resolution map of its kind.

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.