Search

Differential interferometry of synthetic aperture radar (DInSAR) can be used to generate high-precision surface displacement maps in continuous permafrost environments, capturing isotropic surface subsidence and uplift associated with the seasonal freeze and thaw cycle. We generated seasonal displacement maps using DInSAR with ultrafine-beam Radarsat-2 data for the summers of 2013, 2015, and 2019 at Cape Bounty, Melville Island, and examined them in combination with a land-cover classification, meteorological data, topographic data, optical satellite imagery, and in situ measures of soil moisture, soil temperature, and depth to the frost table. Over the three years studied, displacement magnitudes (estimated uncertainty ± 1 cm) of up to 10 cm per 48-day DInSAR stack were detected. However, generally, the displacement was far smaller (up to 4 cm). Surface displacement was found to be most extensive and of the greatest magnitude in low-lying, wet, and steeply sloping areas. The few areas where large vertical displacements (>2.5 cm) were detected in multiple years were clustered in wet, low lying areas, on steep slopes or ridges, or close to the coast. DInSAR also captured the expansion of two medium-sized retrogressive thaw slumps (RTS), exhibiting widespread negative surface change in the slump floor.

Airborne laser scanning (ALS) systems tuned to the near-infrared (NIR; 1064 nm) wavelength have become the best available data source for characterizing vegetation structure. Proliferation of multi-spectral ALS (M-ALS) data with lasers tuned at two additional wavelengths (commonly 532 nm; green, and 1550 nm; short-wave infrared (SWIR)) has promoted interest in the benefit of additional wavelengths for forest inventory modelling. In this study, structural and intensity based M-ALS metrics were derived from wavelengths independently and combined to assess their value for modelling forest inventory attributes (Lorey’s height (HL), gross volume (V), and basal area (BA)) and overstorey species diversity (Shannon index (H), Simpson index (D), and species richness (R)) in a diverse mixed-wood forest in Ontario, Canada. The area-based approach (ABA) to forest attribute modelling was used, where structural- and intensity-based metrics were calculated and used as inputs for random forest models. Structural metrics from the SWIR channel (SWIRstruc) were found to be the most accurate for H and R (%RMSE = 14.3 and 14.9), and NIRstruc were most accurate for V (%RMSE = 20.4). The addition of intensity metrics marginally increased the accuracy of HL models for SWIR and combined channels (%RMSE = 7.5). Additionally, a multi-resolution (0.5, 1, 2 m) voxel analysis was performed, where intensity data were used to calculate a suite of spectral indices. Plot-level summaries of spectral indices from each voxel resolution alone, as well as combined with structural metrics from the NIR wavelength, were used as random forest predictors. The addition of structural metrics from the NIR band reduced %RMSE for all models with HL, BA, and V realizing the largest improvements. Intensity metrics were found to be important variables in the 1 m and 2 m voxel models for D and H. Overall, results indicated that structural metrics were the most appropriate. However, the inclusion of intensity metrics, and continued testing of their potential for modelling diversity indices is warranted, given minor improvements when included. Continued analyses using M-ALS intensity metrics and voxel-based indices would help to better understand the value of these data, and their future role in forest management.

Climate warming has led to an urgent need for improved estimates of carbon accumulation in uneven-aged, mixed temperate forests, where high uncertainty remains. We investigated the feasibility of using LiDAR-derived forest attributes to initialize a growth and yield (G&Y) model in complex stands at the Petawawa Research Forest (PRF) in eastern Ontario, Canada; i.e., can G&Y models based on LiDAR provide accurate predictions of aboveground carbon accumulation in complex forests compared to traditional inventory-based estimates? Applying a local G&Y model, we forecasted aboveground carbon stock (tons/ha) and accumulation (tons/ha/yr) using recurring plot measurements from 2012−2016, FVS1. We applied statistical predictors derived from LiDAR to predict stem density (SD), stem diameter distribution (SDD), and basal area distribution (BA_dist). These data, along with measured species abundance, were used to initialize a second model (FVS2). A third model was tested using LiDAR-initialized tree lists and photo-interpreted estimates of species abundance (i.e., FVS3). The carbon stock projections for 2016 from the inventory-based G&Y model) were equivalent to validation carbon stocks measured in 2016 at all size-class levels (p < 0.05), while LiDAR-based G&Y models were not. None of the models were equivalent to validation data for accumulation (p > 0.05). At the plot level, LiDAR-based predictions of carbon accumulation over a nine-year period did not differ when using either inventory or photo-interpreted species (p < 0.05). Using a constant mortality rate, we also found statistical equivalency of inventory and photo-interpreted accumulation models for all size classes ≥17 cm. These results suggest that more precise information is needed on tree characteristics than we could derive from LiDAR, but that plot-level species information is not as critical for predictions of carbon accumulation in mixed-species forests. Further work is needed on the use of LiDAR to quantify stand properties before this technique can be used to replace recurring plot measurements to quantify carbon accumulation.

Information on foliar macronutrients is required in order to understand plant physiological and ecosystem processes such as photosynthesis, nutrient cycling, respiration and cell wall formation. The ability to measure, model and map foliar macronutrients (nitrogen (N), phosphorus (P), potassium (K), calcium (Ca) and magnesium (Mg)) at the forest canopy level provides information on the spatial patterns of ecosystem processes (e.g., carbon exchange) and provides insight on forest condition and stress. Imaging spectroscopy (IS) has been used particularly for modeling N, using airborne and satellite imagery mostly in temperate and tropical forests. However, there has been very little research conducted at these scales to model P, K, Ca, and Mg and few studies have focused on boreal forests. We report results of a study of macronutrient modeling using spaceborne IS and airborne light detection and ranging (LiDAR) data for a mixedwood boreal forest canopy in northern Ontario, Canada. Models incorporating Hyperion data explained approximately 90% of the variation in canopy concentrations of N, P, and Mg; whereas the inclusion of LiDAR data significantly improved the prediction of canopy concentration of Ca (R2 = 0.80). The combined used of IS and LiDAR data significantly improved the prediction accuracy of canopy Ca and K concentration but decreased the prediction accuracy of canopy P concentration. The results indicate that the variability of macronutrient concentration due to interspecific and functional type differences at the site provides the basis for the relationship observed between the remote sensing measurements (i.e., IS and LiDAR) and macronutrient concentration. Crown closure and canopy height are the structural metrics that establish the connection between macronutrient concentration and IS and LiDAR data, respectively. The spatial distribution of macronutrient concentration at the canopy scale mimics functional type distribution at the site. The ability to predict canopy N, P, K, Ca and Mg in this study using only IS, only LiDAR or their combination demonstrates the excellent potential for mapping these macronutrients at canopy scales across larger geographic areas into the next decade with the launch of new IS satellite missions and by using spaceborne LiDAR data.

Vegetation in the Arctic is often sparse, spatially heterogeneous, and difficult to model. Synthetic Aperture Radar (SAR) has shown some promise in above-ground phytomass estimation at sub-arctic latitudes, but the utility of this type of data is not known in the context of the unique environments of the Canadian High Arctic. In this paper, Artificial Neural Networks (ANNs) were created to model the relationship between variables derived from high resolution multi-incidence angle RADARSAT-2 SAR data and optically-derived (GeoEye-1) Soil Adjusted Vegetation Index (SAVI) values. The modeled SAVI values (i.e., from SAR variables) were then used to create maps of above-ground phytomass across the study area. SAVI model results for individual ecological classes of polar semi-desert, mesic heath, wet sedge, and felsenmeer were reasonable, with r2 values of 0.43, 0.43, 0.30, and 0.59, respectively. When the outputs of these models were combined to analyze the relationship between the model output and SAVI as a group, the r2 value was 0.60, with an 8% normalized root mean square error (% of the total range of phytomass values), a positive indicator of a relationship. The above-ground phytomass model also resulted in a very strong relationship (r2 = 0.87) between SAR-modeled and field-measured phytomass. A positive relationship was also found between optically derived SAVI values and field measured phytomass (r2 = 0.79). These relationships demonstrate the utility of SAR data, compared to using optical data alone, for modeling above-ground phytomass in a high arctic environment possessing relatively low levels of vegetation.

Leaf Area Index (LAI) is an important input variable for forest ecosystem modeling as it is a factor in predicting productivity and biomass, two key aspects of forest health. Current in situ methods of determining LAI are sometimes destructive and generally very time consuming. Other LAI derivation methods, mainly satellite-based in nature, do not provide sufficient spatial resolution or the precision required by forest managers for tactical planning. This paper focuses on estimating LAI from: (i) height and density metrics derived from Light Detection and Ranging (LiDAR); (ii) spectral vegetation indices (SVIs), in particular the Normalized Difference Vegetation Index (NDVI); and (iii) a combination of these methods. For the Hearst Forest of Northern Ontario, in situ measurements of LAI were derived from digital hemispherical photographs (DHPs) while remote sensing variables were derived from low density LiDAR (i.e., 1 m−2) and high spatial resolution WorldView-2 data (2 m). Multiple Linear Regression (MLR) models were generated using these variables. Results from these analyses demonstrate: (i) moderate explanatory power (i.e., R2 = 0.53) for LiDAR height and density metrics that have proven to be related to canopy structure; (ii) no relationship when using SVIs; and (iii) no significant improvement of LiDAR models when combining them with SVI variables. The results suggest that LiDAR models in boreal forest environments provide satisfactory estimations of LAI, even with narrow ranges of LAI for model calibration. Models derived from low point density LiDAR in a mixedwood boreal environment seem to offer a reliable method of estimating LAI at high spatial resolution for decision makers in the forestry community. This method can be easily incorporated into simultaneous modeling efforts for forest inventory variables using LiDAR.

As a result of the warming observed at high latitudes, there is significant potential for the balance of ecosystem processes to change, i.e., the balance between carbon sequestration and respiration may be altered, giving rise to the release of soil carbon through elevated ecosystem respiration. Gross ecosystem productivity and ecosystem respiration vary in relation to the pattern of vegetation community type and associated biophysical traits (e.g., percent cover, biomass, chlorophyll concentration, etc.). In an arctic environment where vegetation is highly variable across the landscape, the use of high spatial resolution imagery can assist in discerning complex patterns of vegetation and biophysical variables. The research presented here examines the relationship between ecological and spectral variables in order to generate an ecologically meaningful vegetation classification from high spatial resolution remote sensing data. Our methodology integrates ordination and image classifications techniques for two non-overlapping Arctic sites across a 5° latitudinal gradient (approximately 70° to 75°N). Ordination techniques were applied to determine the arrangement of sample sites, in relation to environmental variables, followed by cluster analysis to create ecological classes. The derived classes were then used to classify high spatial resolution IKONOS multispectral data. The results demonstrate moderate levels of success. Classifications had overall accuracies between 69%–79% and Kappa values of 0.54–0.69. Vegetation classes were generally distinct at each site with the exception of sedge wetlands. Based on the results presented here, the combination of ecological and remote sensing techniques can produce classifications that have ecological meaning and are spectrally separable in an arctic environment. These classification schemes are critical for modeling ecosystem processes.

Over the past two decades there has been an abundance of research demonstrating the utility of airborne light detection and ranging (LiDAR) for predicting forest biophysical/inventory variables at the plot and stand levels. However, to date there has been little effort to develop a set of protocols for data acquisition and processing that would move governments or the forest industry towards cost-effective implementation of this technology for strategic and tactical (i.e., operational) forest resource inventories. The goal of this paper is to initiate this process by examining the significance of LiDAR data acquisition (i.e., point density) for modeling forest inventory variables for the range of species and stand conditions representing much of Ontario, Canada. Field data for approximately 200 plots, sampling a broad range of forest types and conditions across Ontario, were collected for three study sites. Airborne LiDAR data, characterized by a mean density of 3.2 pulses m−2 were systematically decimated to produce additional datasets with densities of approximately 1.6 and 0.5 pulses m−2. Stepwise regression models, incorporating LiDAR height and density metrics, were developed for each of the three LiDAR datasets across a range of forest types to estimate the following forest inventory variables: (1) average height (R2(adj) = 0.75–0.95); (2) top height (R2(adj) = 0.74–0.98); (3) quadratic mean diameter (R2(adj) = 0.55–0.85); (4) basal area (R2(adj) = 0.22–0.93); (5) gross total volume (R2(adj) = 0.42–0.94); (6) gross merchantable volume (R2(adj) = 0.35–0.93); (7) total aboveground biomass (R2(adj) = 0.23–0.93); and (8) stem density (R2(adj) = 0.17–0.86). Aside from a few cases (i.e., average height and density for some stand types), no decimation effect was observed with respect to the precision of the prediction of the majority of forest variables, which suggests that a mean density of 0.5 pulses m−2 is sufficient for plot and stand level modeling under these diverse forest conditions across Ontario.