Integrating Landsat pixel composites and change metrics with lidar plots to predictively map forest structure and aboveground biomass in Saskatchewan, Canada

Forest inventory and monitoring programs are needed to provide timely, spatially complete (i.e. mapped), and verifiable information to support forest management, policy formulation, and reporting obligations. Satellite images, in particular data from the Landsat Thematic Mapper and Enhanced Thematic Mapper (TM/ETM +) sensors, are often integrated with field plots from forest inventory programs, leveraging the complete spatial coverage of imagery with detailed ecological information from a sample of plots to spatially model forest conditions and resources. However, in remote and unmanaged areas such as Canada's northern forests, financial and logistic constraints can severely limit the availability of inventory plot data. Additionally, Landsat spectral information has known limitations for characterizing vertical vegetation structure and biomass; while clouds, snow, and short growing seasons can limit development of large area image mosaics that are spectrally and phenologically consistent across space and time. In this study we predict and map forest structure and aboveground biomass over 37 million ha of forestland in Saskatchewan, Canada. We utilize lidar plots—observations of forest structure collected from airborne discrete-return lidar transects acquired in 2010—as a surrogate for traditional field and photo plots. Mapped explanatory data included Tasseled Cap indices and multi-temporal change metrics derived from Landsat TM/ETM + pixel-based image composites. Maps of forest structure and total aboveground biomass were created using a Random Forest (RF) implementation of Nearest Neighbor (NN) imputation. The imputation model had moderate to high plot-level accuracy across all forest attributes (R2 values of 0.42–0.69), as well as reasonable attribute predictions and error estimates (for example, canopy cover above 2 m on validation plots averaged 35.77%, with an RMSE of 13.45%, while unsystematic and systematic agreement coefficients (ACuns and ACsys) had values of 0.63 and 0.97 respectively). Additionally, forest attributes displayed consistent trends in relation to the time since and magnitude of wildfires, indicating model predictions captured the dominant ecological patterns and processes in these forests. Acknowledging methodological and conceptual challenges based upon the use of lidar plots from transects, this study demonstrates that using lidar plots and pixel compositing in imputation mapping can provide forest inventory and monitoring information for regions lacking ongoing or up-to-date field data collection programs

Ecosystem classifications based on summer and winter conditions

Ecosystem classifications map an area into relatively homogenous units for environmental research, monitoring, and management. However, their effectiveness is rarely tested. Here, three classifications are (1) defined and characterized for Canada along summertime productivity (moderate-resolution imaging spectrometer fraction of absorbed photosynthetically active radiation) and wintertime snow conditions (special sensor microwave/imager snow water equivalent), independently and in combination, and (2) comparatively evaluated to determine the ability of each classification to represent the spatial and environmental patterns of alternative schemes, including the Canadian ecozone framework. All classifications depicted similar patterns across Canada, but detailed class distributions differed. Class spatial characteristics varied with environmental conditions within classifications, but were comparable between classifications. There was moderate correspondence between classifications. The strongest association was between productivity classes and ecozones. The classification along both productivity and snow balanced these two sets of variables, yielding intermediate levels of association in all pairwise comparisons. Despite relatively low spatial agreement between classifications, they successfully captured patterns of the environmental conditions underlying alternate schemes (e.g.; snow classes explained variation in productivity and vice versa). The performance of ecosystem classifications and the relevance of their input variables depend on the environmental patterns and processes used for applications and evaluation. Productivity or snow regimes, as constructed here, may be desirable when summarizing patterns controlled by summer- or wintertime conditions, respectively, or of climate change responses. General purpose ecosystem classifications should include both sets of drivers. Classifications should be carefully, quantitatively, and comparatively evaluated relative to a particular application prior to their implementation as monitoring and assessment frameworks

Forest inventory stand height estimates from very high spatial resolution satellite imagery calibrated with lidar plots

Many areas of forest across northern Canada are challenging to monitor on a regular basis as a result of their large extent and remoteness. Although no forest inventory data typically exist for these northern areas, detailed and timely forest information for these areas is required to support national and international reporting obligations. We developed and tested a sample-based approach that could be used to estimate forest stand height in these remote forests using panchromatic Very High Spatial Resolution (VHSR, <1 m) optical imagery and light detection and ranging (lidar) data. Using a study area in central British Columbia, Canada, to test our approach, we compared four different methods for estimating stand height using stand-level and crown-level metrics generated from the VHSR imagery. 'Lidar plots' (voxel-based samples of lidar data) are used for calibration and validation of the VHSR-based stand height estimates, similar to the way that field plots are used to calibrate photogrammetric estimates of stand height in a conventional forest inventory or to make empirical attribute estimates from multispectral digital remotely sensed data. A k-nearest neighbours (k-NN) method provided the best estimate of mean stand height (R-2 = 0.69; RMSE = 2.3 m, RMSE normalized by the mean value of the estimates (RMSE-%) = 21) compared with linear regression, random forests, and regression tree methods. The approach presented herein demonstrates the potential of VHSR panchromatic imagery and lidar to provide robust and representative estimates of stand height in remote forest areas where conventional forest inventory approaches are either too costly or are not logistically feasible. While further evaluation of the methods is required to generalize these results over Canada to provide robust and representative estimation, VHSR and lidar data provide an opportunity for monitoring in areas for which there is no detailed forest inventory information available