A Comparison of Three Airborne Laser Scanner Types for Species Identification of Individual Trees.

Species identification is a critical factor for obtaining accurate forest inventories. This paper compares the same method of tree species identification (at the individual crown level) across three different types of airborne laser scanning systems (ALS): two linear lidar systems (monospectral and multispectral) and one single-photon lidar (SPL) system to ascertain whether current individual tree crown (ITC) species classification methods are applicable across all sensors. SPL is a new type of sensor that promises comparable point densities from higher flight altitudes, thereby increasing lidar coverage. Initial results indicate that the methods are indeed applicable across all of the three sensor types with broadly similar overall accuracies (Hardwood/Softwood, 83-90%; 12 species, 46-54%; 4 species, 68-79%), with SPL being slightly lower in all cases. The additional intensity features that are provided by multispectral ALS appear to be more beneficial to overall accuracy than the higher point density of SPL. We also demonstrate the potential contribution of lidar time-series data in improving classification accuracy (Hardwood/Softwood, 91%; 12 species, 58%; 4 species, 84%). Possible causes for lower SPL accuracy are (a) differences in the nature of the intensity features and (b) differences in first and second return distributions between the two linear systems and SPL. We also show that segmentation (and field-identified training crowns deriving from segmentation) that is performed on an initial dataset can be used on subsequent datasets with similar overall accuracy. To our knowledge, this is the first study to compare these three types of ALS systems for species identification at the individual tree level.

Spatial modeling of daily concentrations of ground-level ozone in Montreal, Canada: A comparison of geostatistical approaches.

Ground-level ozone (O3) is a powerful oxidizing agent and a harmful pollutant affecting human health, forests and crops. Estimating O3 exposure is a challenge because it exhibits complex spatiotemporal patterns. The aim in this study was to provide high-resolution maps (100 × 100 m) of O3 for the metropolitan area of Montreal, Canada. We assessed the kriging with external drift (KED) model to estimate O3 concentration by synoptic weather classes for 2010. We compared these results with ordinary kriging (OK), and a simple average of 12 monitoring stations. We also compared the estimates obtained for the 2010 summer with those from a Bayesian maximum entropy (BME) model reported in the literature (Adam-Poupart et al., 2014). The KED model with road and vegetation density as covariates showed good performance for all six synoptic classes (daily R2 estimates ranging from 0.77 to 0.92 and RMSE from 2.79 to 3.37 ppb). For the summer of 2010, the model using KED demonstrated the best results (R2 = 0.92; RMSE = 3.14 ppb), followed by the OK model (R2 = 0.85, RMSE = 4 ppb). Our results showed that errors appear to be substantially reduced with the KED model. This may increase our capacity of linking O3 levels to health problems by means of improved assessments of ambient exposures. However, future work integrating the temporal dependency in the data is needed to not overstate the performance of the KED model.

Response of a boreal forest to canopy opening: assessing vertical and lateral tree growth with multi-temporal lidar data.

Fine-scale height-growth response of boreal trees to canopy openings is difficult to measure from the ground, and there are important limitations in using stereophotogrammetry in defining gaps and determining individual crowns and height. However, precise knowledge on height growth response to different openings is critical for refining partial harvesting techniques. In this study, we question whether conifers and hardwoods respond equally in terms of sapling growth or lateral growth to openings. We also ask to what distance gaps affect tree growth into the forest. We use multi-temporal lidar to characterize tree/sapling height and lateral growth responses over five years to canopy openings and high resolution images to identify conifers and hardwoods. Species-class-wise height-growth patterns of trees/saplings in various neighborhood contexts were determined across a 6-km matrix of Canadian boreal mixed deciduous coniferous forests. We then use statistical techniques to probe how these growth responses vary by spatial location with respect to the gap edge. Results confirm that both mechanisms of gap closure contribute to the closing of canopies at a rate of 1.2% per annum. Evidence also shows that both hardwood and conifer gap edge trees have a similar lateral growth (average of 22 cm/yr) and similar rates of height growth irrespective of their location and initial height. Height growth of all saplings, however, was strongly dependent on their position within the gap and the size of the gap. Results suggest that hardwood and softwood saplings in gaps have greatest growth rates at distances of 0.5-2 m and 1.5-4 m from the gap edge and in openings smaller than 800 m2 and 250 m2, respectively. Gap effects on the height growth of trees in the intact forest were evident up to 30 m and 20 m from gap edges for hardwood and softwood overstory trees, respectively. Our results thus suggest that foresters should consider silvicultural techniques that create many small openings in mixed coniferous deciduous boreal forests to maximize the growth response of both residual and regenerating trees.

Implications of differing input data sources and approaches upon forest carbon stock estimation.

Site index is an important forest inventory attribute that relates productivity and growth expectation of forests over time. In forest inventory programs, site index is used in conjunction with other forest inventory attributes (i.e., height, age) for the estimation of stand volume. In turn, stand volumes are used to estimate biomass (and biomass components) and enable conversion to carbon. In this research, we explore the implications and consequences of different estimates of site index on carbon stock characterization for a 2,500-ha Douglas-fir-dominated landscape located on Eastern Vancouver Island, British Columbia, Canada. We compared site index estimates from an existing forest inventory to estimates generated from a combination of forest inventory and light detection and ranging (LIDAR)-derived attributes and then examined the resultant differences in biomass estimates generated from a carbon budget model (Carbon Budget Model of the Canadian Forest Sector (CBM-CFS3)). Significant differences were found between the original and LIDAR-derived site indices for all species types and for the resulting 5-m site classes (p < 0.001). The LIDAR-derived site class was greater than the original site class for 42% of stands; however, 77% of stands were within +/-1 site class of the original class. Differences in biomass estimates between the model scenarios were significant for both total stand biomass and biomass per hectare (p < 0.001); differences for Douglas-fir-dominated stands (representing 85% of all stands) were not significant (p = 0.288). Overall, the relationship between the two biomass estimates was strong (R(2) = 0.92, p < 0.001), suggesting that in certain circumstances, LIDAR may have a role to play in site index estimation and biomass mapping.

Spatio-temporal models to estimate daily concentrations of fine particulate matter in Montreal: Kriging with external drift and inverse distance-weighted approaches.

Air pollution is a major environmental and health problem, especially in urban agglomerations. Estimating personal exposure to fine particulate matter (PM2.5) remains a great challenge because it requires numerous point measurements to explain the daily spatial variation in pollutant levels. Furthermore, meteorological variables have considerable effects on the dispersion and distribution of pollutants, which also depends on spatio-temporal emission patterns. In this study we developed a hybrid interpolation technique that combined the inverse distance-weighted (IDW) method with Kriging with external drift (KED), and applied it to daily PM2.5 levels observed at 10 monitoring stations. This provided us with downscaled high-resolution maps of PM2.5 for the Island of Montreal. For the KED interpolation, we used spatio-temporal daily meteorological estimates and spatial covariates as land use and vegetation density. Different KED and IDW daily estimation models for the year 2010 were developed for each of the six synoptic weather classes. These clusters were developed using principal component analysis and unsupervised hierarchical classification. The results of the interpolation models were assessed with a leave-one-station-out cross-validation. The performance of the hybrid model was better than that of the KED or the IDW alone for all six synoptic weather classes (the daily estimate for R(2) was 0.66-0.93 and for root mean square error (RMSE) 2.54-1.89 µg/m(3)).