A Forest Attribute Mapping Framework: A Pilot Study in a Northern Boreal Forest, Northwest Territories, Canada

Improved k-NN Mapping of Forest Attributes in Northern Canada Using Spaceborne L-Band SAR, Multispectral and LiDAR Data

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 above-ground 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.

The Multisource Vegetation Inventory (MVI): A Satellite-Based Forest Inventory for the Northwest Territories Taiga Plains

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.

An Improved Generalized Hierarchical Estimation Framework with Geostatistics for Mapping Forest Parameters and Its Uncertainty: A Case Study of Forest Canopy Height

Forest canopy height is an essential parameter in estimating forest aboveground biomass (AGB), growing stock volume (GSV), and carbon storage, and it can provide necessary information in forest management activities. Light direction and ranging (LiDAR) is widely used for estimating canopy height. Considering the high cost of acquiring LiDAR data over large areas, we took a two-stage up-scaling approach in estimating forest canopy height and aimed to develop a method for quantifying the uncertainty of the estimation result. Based on the generalized hierarchical model-based (GHMB) estimation framework, a new estimation framework named RK-GHMB that makes use of a geostatistical method (regression kriging, RK) was developed. In this framework, the wall-to-wall forest canopy height and corresponding uncertainty in map unit scale are generated. This study was carried out by integrating plot data, sampled airborne LiDAR data, and wall-to-wall Ziyuan-3 satellite (ZY3) stereo images. The result shows that RK-GHMB can obtain a similar estimation accuracy (r = 0.92, MAE = 1.50 m) to GHMB (r = 0.92, MAE = 1.52 m) with plot-based reference data. For LiDAR-based reference data, the accuracy of RK-GHMB (r = 0.78, MAE = 1.75 m) is higher than that of GHMB (r = 0.75, MAE = 1.85 m). The uncertainties for all map units range from 1.54 to 3.60 m for the RK-GHMB results. The values change between 1.84 and 3.60 m for GHMB. This study demonstrates that this two-stage up-scaling approach can be used to monitor forest canopy height. The proposed RK-GHMB approach considers the spatial autocorrelation of neighboring data in the second modeling stage and can achieve a higher accuracy.

Shrub changes with proximity to anthropogenic disturbance in boreal wetlands determined using bi-temporal airborne lidar in the Oil Sands Region, Alberta Canada

In this study, we used bi-temporal airborne lidar data to compare changes in vegetation height proximal to anthropogenic disturbances in the Oil Sands Region of Alberta, Canada. We hypothesize that relatively low-impact disturbances such as seismic lines will increase the fragmentation of wetlands, resulting in shrub growth. Bi-temporal lidar data collected circa 2008 and 2018 were used to identify correspondence between the density of anthropogenic disturbances, wetland shape complexity and changes in vegetation height within >1800 wetlands near Fort McKay, Alberta, Canada. We found that up to 50% of wetlands were disturbed by anthropogenic disturbance in some parts of the region, with the highest proportional disturbance occurring within fens. Areas of dense anthropogenic disturbance in bogs resulted in increased growth and expansion of shrubs, while we found the opposite to occur in fens and swamps during the 10-year period. Up to 30% of bogs had increased shrubification, while shrub changes in fens and swamps varied depending on density of disturbance and did not necessarily correspond with shrub growth. As wetland shapes became increasingly elongated, the prevalence of shrubs declined between the two time periods, which may be associated with hydrological drivers (e.g. elongated may indicate surface and ground-water discharge influences). The results of this study indicate that linear disturbances such as seismic lines, considered to have relatively minimal impacts on ecosystems, can impact proximal wetland shape, fragmentation and vegetation community changes, especially in bogs.

Remote Sensing of Boreal Wetlands 2: Methods for Evaluating Boreal Wetland Ecosystem State and Drivers of Change

The following review is the second part of a two part series on the use of remotely sensed data for quantifying wetland extent and inferring or measuring condition for monitoring drivers of change on wetland environments. In the first part, we introduce policy makers and non-users of remotely sensed data with an effective feasibility guide on how data can be used. In the current review, we explore the more technical aspects of remotely sensed data processing and analysis using case studies within the literature. Here we describe: (a) current technologies used for wetland assessment and monitoring; (b) the latest algorithmic developments for wetland assessment; (c) new technologies; and (d) a framework for wetland sampling in support of remotely sensed data collection. Results illustrate that high or fine spatial resolution pixels (≤10 m) are critical for identifying wetland boundaries and extent, and wetland class, form and type, but are not required for all wetland sizes. Average accuracies can be up to 11% better (on average) than medium resolution (11–30 m) data pixels when compared with field validation. Wetland size is also a critical factor such that large wetlands may be almost as accurately classified using medium-resolution data (average = 76% accuracy, stdev = 21%). Decision-tree and machine learning algorithms provide the most accurate wetland classification methods currently available, however, these also require sampling of all permutations of variability. Hydroperiod accuracy, which is dependent on instantaneous water extent for single time period datasets does not vary greatly with pixel resolution when compared with field data (average = 87%, 86%) for high and medium resolution pixels, respectively. The results of this review provide users with a guideline for optimal use of remotely sensed data and suggested field methods for boreal and global wetland studies.

Remote Sensing of Boreal Wetlands 1: Data Use for Policy and Management

Wetlands have and continue to undergo rapid environmental and anthropogenic modification and change to their extent, condition, and therefore, ecosystem services. In this first part of a two-part review, we provide decision-makers with an overview on the use of remote sensing technologies for the ‘wise use of wetlands’, following Ramsar Convention protocols. The objectives of this review are to provide: (1) a synthesis of the history of remote sensing of wetlands, (2) a feasibility study to quantify the accuracy of remotely sensed data products when compared with field data based on 286 comparisons found in the literature from 209 articles, (3) recommendations for best approaches based on case studies, and (4) a decision tree to assist users and policymakers at numerous governmental levels and industrial agencies to identify optimal remote sensing approaches based on needs, feasibility, and cost. We argue that in order for remote sensing approaches to be adopted by wetland scientists, land-use managers, and policymakers, there is a need for greater understanding of the use of remote sensing for wetland inventory, condition, and underlying processes at scales relevant for management and policy decisions. The literature review focuses on boreal wetlands primarily from a Canadian perspective, but the results are broadly applicable to policymakers and wetland scientists globally, providing knowledge on how to best incorporate remotely sensed data into their monitoring and measurement procedures. This is the first review quantifying the accuracy and feasibility of remotely sensed data and data combinations needed for monitoring and assessment. These include, baseline classification for wetland inventory, monitoring through time, and prediction of ecosystem processes from individual wetlands to a national scale.

Extending ALS-Based Mapping of Forest Attributes with Medium Resolution Satellite and Environmental Data

Airborne laser scanner (ALS) data are used to map a range of forest inventory attributes at operational scales. However, when wall-to-wall ALS coverage is cost prohibitive or logistically challenging, alternative approaches are needed for forest mapping. We evaluated an indirect approach for extending ALS-based maps of forest attributes using medium resolution satellite and environmental data. First, we developed ALS-based models and predicted a suite of forest attributes for a 950 km2 study area covered by wall-to-wall ALS data. Then, we used samples extracted from the ALS-based predictions to model and map these attributes with satellite and environmental data for an extended 5600 km2 area with similar forest and ecological conditions. All attributes were predicted well with the ALS data (R2 ≥ 0.83; RMSD% < 26). The satellite and environmental models developed using the ALS-based predictions resulted in increased correspondence between observed and predicted values by 13–49% and decreased prediction errors by 8–28% compared with models developed directly with the ground plots. Improvements were observed for both multiple regression and random forest models, and for the suite of forest attributes assessed. We concluded that the use of ALS-based predictions in this study improved the estimation of forest attributes beyond an approach linking ground plots directly to the satellite and environmental data.

SAR and Lidar Temporal Data Fusion Approaches to Boreal Wetland Ecosystem Monitoring

The objective of this study was to develop a decision-based methodology, focused on data fusion for wetland classification based on surface water hydroperiod and associated riparian (transitional area between aquatic and upland zones) vegetation community attributes. Multi-temporal, multi-mode data were examined from airborne Lidar (Teledyne Optech, Inc., Toronto, ON, Canada, Titan), synthetic aperture radar (Radarsat-2, single and quad polarization), and optical (SPOT) sensors with near-coincident acquisition dates. Results were compared with 31 field measurement points for six wetlands at riparian transition zones and surface water extents in the Utikuma Regional Study Area (URSA). The methodology was repeated in the Peace-Athabasca Delta (PAD) to determine the transferability of the methods to other boreal environments. Water mask frequency analysis showed accuracies of 93% to 97%, and kappa values of 0.8–0.9 when compared to optical data. Concordance results comparing the semi-permanent/permanent hydroperiod between 2015 and 2016 were found to be 98% similar, suggesting little change in wetland surface water extent between these two years. The results illustrate that the decision-based methodology and data fusion could be applied to a wide range of boreal wetland types and, so far, is not geographically limited. This provides a platform for land use permitting, reclamation monitoring, and wetland regulation in a region of rapid development and uncertainty due to climate change. The methodology offers an innovative time series-based boreal wetland classification approach using data fusion of multiple remote sensing data sources.