Mapping whitebark pine mortality caused by a mountain pine beetle outbreak with high spatial resolution satellite imagery

Insect outbreaks cause significant tree mortality across western North America, including in high-elevation whitebark pine forests. These forests are under several threats, which include attack by insects and white pine blister rust, as well as conversion to other tree species as a result of fire suppression. Mapping tree mortality is critical to determining the status of whitebark pine as a species. Satellite remote sensing builds upon existing aerial surveys by using objective, repeatable methods that can result in high spatial resolution monitoring. Past studies concentrated on level terrain and only forest vegetation type. The objective of this study was to develop a means of classifying whitebark pine mortality caused by a mountain pine beetle infestation in rugged, remote terrain using high spatial resolution satellite imagery. We overcame three challenges of mapping mortality in this mountainous region: (1) separating non-vegetated cover types, green and brown herbaceous cover, green (live) tree cover, and red-attack (dead) tree cover; (2) variations in illumination as a result of variations in slope and aspect related to the mountainous terrain of the study site; and (3) the difficulty of georegistering the imagery for use in comparing field measurements. Quickbird multi-spectral imagery (2.4 m spatial resolution) was used, together with a maximum likelihood classification method, to classify vegetation cover types over a 6400 ha area. To train the classifier, we selected pixels in each cover class from the imagery guided by our knowledge of the study site. Variables used in the maximum likelihood classifier included the ratio of red reflectance to green reflectance as well as green reflectance. These variables were stratified by solar incidence angle to account for illumination variability. We evaluated the results of the classified image using a reserved set of image-derived class members and field measurements of live and dead trees. Classification results yielded high overall accuracy (86% and 91% using image-derived class members and field measurements respectively) and kappa statistics (0.82 and 0.82) and low commission (0.9% and 1.5%) and omission (6.5% and 15.9%) errors for the red-attack tree class. Across the scene, 700 ha or 31% of the forest was identified as in the red-attack stage. Severity (percent mortality by canopy cover) varied from nearly 100% for some areas to regions with little mortality. These results suggest that high spatial resolution satellite imagery can provide valuable information for mapping and monitoring tree mortality even in rugged, mountainous terrain.     

Assessment of QuickBird high spatial resolution imagery to detect red attack damage due to mountain pine beetle infestation

High spatial resolution remotely sensed data has the potential to complement existing forest health programs for both strategic planning over large areas, as well as for detailed and precise identification of tree crowns subject to stress and infestation. The area impacted by the current mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has increased 40-fold over the previous 5 years, with approximately 8.5 million ha of forest infested in 2005. As a result of the spatial extent and intensity of the outbreak, new technologies are being assessed to help detect, map, and monitor the damage caused by the beetle, and to inform mitigation of future beetle outbreaks. In this paper, we evaluate the capacity of high spatial resolution QuickBird multi-spectral imagery to detect mountain pine beetle red attack damage. ANOVA testing of individual spectral bands, as well as the Normalized Difference Vegetation Index (NDVI) and a ratio of red to green reflectance (Red-Green Index or RGI), indicated that the RGI was the most successful (p < 0.001) at separating non-attack crowns from red attack crowns. Based on this result, the RGI was subsequently used to develop a binary classification of red attack and non-attack pixels. The total number of QuickBird pixels classified as having red attack damage within a 50 m buffer of a known forest health survey point were compared to the number of red attack trees recorded at the time of the forest health survey. The relationship between the number of red attack pixels and observed red attack crowns was assessed using independent validation data and was found to be significant (r² = 0.48, p < 0.001, standard error = 2.8 crowns). A comparison of the number of QuickBird pixels classified as red attack, and a broader scale index of mountain pine beetle red attack damage (Enhanced Wetness Difference Index, calculated from a time series of Landsat imagery), was significant (r² = 0.61, p < 0.001, standard error = 1.3 crowns). These results suggest that high spatial resolution imagery, in particular QuickBird satellite imagery, has a valuable role to play in identifying tree crowns with red attack damage. This information could subsequently be used to augment existing detailed forest health surveys, calibrate synoptic estimates of red attack damage generated from overview surveys and/or coarse scale remotely sensed data, and facilitate the generation of value-added information products, such as estimates of timber volume impacts at the forest stand level.     

Evaluating the potential of multispectral imagery to map multiple stages of tree mortality

Insect outbreaks are major forest disturbances, causing tree mortality across millions of ha in North America. Resultant spatial and temporal patterns of tree mortality can profoundly affect ecosystem structure and function. In this study, we evaluated the classification accuracy of multispectral imagery at different spatial resolutions. We used four-band digital aerial imagery (30-cm spatial resolution and aggregated to coarser resolutions) acquired over lodgepole pine-dominated stands in central Colorado recently attacked by mountain pine beetle. Classes of interest included green trees and multiple stages of post-insect attack tree mortality, including dead trees with red needles (“red-attack”), dead trees without needles (“gray-attack”), and non-forest. The 30-cm resolution image facilitated delineation of trees located in the field, which were used in image classification. We employed a maximum likelihood classifier using the green band, Red-Green Index (RGI), and Normalized Difference Vegetation Index (NDVI). Pixel-level classification accuracies using this imagery were good (overall accuracy of 87%, kappa = 0.84), although misclassification occurred between a) sunlit crowns of live (green) trees and herbaceous vegetation, and b) sunlit crowns of gray- and red-attack trees and bare soil. We explored the capability of coarser resolution imagery, aggregated from the 30-cm resolution to 1.2, 2.4, and 4.2 m, to improve classification accuracy. We found the highest accuracy at the 2.4-m resolution, where reduction in omission and commission errors and increases in overall accuracy (90%) and kappa (0.88) were achieved, and visual inspection indicated improved mapping. Pixels at this resolution included more shadow in forested regions than pixels in finer resolution imagery, thereby reducing forest canopy reflectance and allowing improved separation between forest and non-forest classes, yet were fine enough to resolve individual tree crowns better than the 4.2-m imagery. Our results illustrate that a classification of an image with a spatial resolution similar to the area of a tree crown outperforms that of finer and coarser resolution imagery for mapping tree mortality and non-forest classes. We also demonstrate that multispectral imagery can be used to separate multiple postoutbreak attack stages (i.e., red-attack and gray-attack) from other classes in the image.     

Continuous wavelet analysis for the detection of green attack damage due to mountain pine beetle infestation

Mountain pine beetle (Dendroctonus ponderosae Hopkins) is the most destructive insect infesting mature pine forests in North America and has devastated millions of hectares of forest in western Canada. Past studies have demonstrated the use of multispectral imagery for remote identification and mapping of visible or red attack damage in forests. This study aims to detect pre-visual or green attack damage in lodgepole pine needles by means of hyperspectral measurements, particularly via continuous wavelet analysis. Field measurements of lodgepole pine stands were conducted at two sites located northwest of Edmonton, Alberta, Canada. In June and August of 2007, reflectance spectra (350-2500 nm) were collected for 16 pairs of trees. Each of the 16 tree pairs included one control tree (healthy), and one stressed tree (girdled to simulate the effects of beetle infestation). In addition, during the period of June through October 2008, spectra were collected from 15 pairs of control- and beetle-infested trees. Spectra derived from these 31 tree pairs were subjected to a continuous wavelet transform, generating a scalogram that compiles the wavelet power as a function of wavelength location and scale of decomposition. Linear relationships were then explored between the wavelet scalograms and chemical properties or class labels (control and non-control) of the sample populations in order to isolate the most useful distinguishing spectral features that related to infested or girdled trees vs. control trees.A deficit in water content is observed in infested trees while an additional deficit in chlorophyll content is seen for girdled trees. The measurable water deficit of infested and girdled tree samples was detectable from the wavelet analysis of the reflectance spectra providing a novel method for the detection of green attack. The spectral features distinguishing control and infested trees are predominantly located between 950 and 1390 nm from scales 1 to 8. Of those, five features between 1318 to 1322 nm at scale 7 are consistently found in the July and August 2008 datasets. These features are located at longer wavelengths than those investigated in previous studies (below 1100 nm) and provide new insights into the potential remote detection of green attack. Spectral features that distinguish control and girdled trees were mostly observed between 1550 and 2370 nm from scales 1 to 5. The differing response of girdled and infested trees appears to indicate that the girdling process does not provide a perfect simulation of the effects caused by beetle infestation.It remains to be determined if the location of the 1318-1322 nm features, near the edge of a strong atmospheric water absorption band, will be sufficiently separable for use in airborne detection of green attack. A plot comparing needle water content and wavelet power at 1320 nm reveals considerable overlap between data derived from both infested and control samples, though the groups are statistically separable. This obstacle may preclude a high accuracy separation of healthy and infected single individuals, but establishing threshold identification levels may provide an economical, efficient and expeditious method for discriminating between healthy and infested tree populations.     

Assessing canopy mortality during a mountain pine beetle outbreak using GeoEye-1 high spatial resolution satellite data

Continuing, severe outbreaks of mountain pine beetle (Dendroctonus ponderosae) across western North America have resulted in widespread mortality of lodgepole pine (Pinus contorta). Multiple studies have used high spatial resolution satellite data to map areas of beetle kill; these studies have largely focused on mapping red canopy cover associated with recent tree mortality and have not examined mapping gray canopy cover that occurs after red needles have dropped. The work presented here examines the use of newly available GeoEye-1 data for mapping both red and gray canopy area in southeastern Wyoming lodgepole pine forest. A 0.5 m spatial resolution, pan-sharpened GeoEye-1 image was used to classify areas of green, red, and gray canopy cover. Reference data were collected at twelve 500 m² field plots. Shadow-normalized green, red, and gray canopy area from classified GeoEye-1 data closely agreed with field-estimated green, red, and gray canopy area. Mean absolute error in canopy cover for the twelve sample plots was 8.3% for the green class, 5.4% for the red class, and 7.2% for the gray class. When all twelve plots were aggregated, remotely sensed estimates of green, red, and gray cover were within 1.7% of the field-estimated cover. Our results demonstrate that high spatial resolution spaceborne multispectral data are a promising tool for mapping canopy mortality caused by mountain pine beetle outbreaks.     

Integrating remotely sensed and ancillary data sources to characterize a mountain pine beetle infestation

The ongoing mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has reached epidemic proportions, with the beetle expanding into geographic areas outside its known biological range. In this study, estimates of red attack damage were derived from a logistic regression model using multi-date Landsat imagery, and ancillary information including terrain attributes and solar radiation. The model estimates were found to be approximately 70% accurate using an independent set of beetle survey data as validation. This probability surface of red attack damage, along with forest inventory and terrain attributes, were used as inputs to decision tree analyses, in order to identify which forest attributes were associated with stands that had a greater likelihood of mountain pine beetle red attack damage. Three distinct decision tree models were developed, with each having a different set of input variables. The results of the analyses indicated that site index (an indicator of the quality of a forest site) and slope were the principal discriminators of the current mountain pine beetle attack, followed by basal area of pine dominated stands, and to a lesser extent, crown closure and stem density. The results suggest that indicators of site quality, particularly site index, could be a complementary addition to existing stand susceptibility rating models.     

Multi-temporal analysis of high spatial resolution imagery for disturbance monitoring

Mountain pine beetle red attack damage has been successfully detected and mapped using single-date high spatial resolution (< 4 m) satellite multi-spectral data. Forest managers; however, need to monitor locations for changes in beetle populations over time. Specifically, counts of individual trees attacked in successive years provide an indication of beetle population growth and dynamics. Surveys are typically used to estimate the ratio of green (current) attack trees to red (previous) attack trees or G:R. In this study, we estimate average stand-level G:R using a time-series of QuickBird multi-spectral and panchromatic satellite data, combined with field data for three forested stands near Merritt, British Columbia, Canada. Using a ratio of QuickBird red to green wavelengths (Red-Green Index or RGI), the change in RGI (ΔRGI) in successive image pairs is used to estimate red attack damage in 2004, 2005, and 2006, with true positive accuracies ranging from 89 to 93%. To overcome issues associated with differing viewing geometry and illumination angles that impair tracking of individual trees through time, segments are generated from the QuickBird multi-spectral data to identify small groups of trees. These segments then serve as the vehicle for monitoring changes in red attack damage over time. A local maxima filter is applied to the panchromatic data to estimate stem counts, thereby allowing an indication of the total stand population at risk of attack. By combining the red attack damage estimates with the local maxima stem counts, predictions are made of the number of attacked trees in a given year. Backcasting the current year's red attack damaged trees as the previous year's green attack facilitates the estimation of an average stand G:R. In this study area, these retrospective G:R values closely match those generated from field surveys. The results of this study indicate that a monitoring program using a time series of high spatial resolution remotely sensed data (multi-spectral and panchromatic) over select sample locations, could be used to estimate G:R over large areas, facilitating landscape level management strategies and/or providing a mechanism for assessing the efficacy of previously implemented strategies.     

Using hyperspectral remote sensing to detect and quantify southeastern pine senescence effects in red-cockaded woodpecker (Picoides borealis) habitat

Conservation of threatened and endangered species requires maintenance of critical habitat. The red-cockaded woodpecker Picoides borealis (RCW) is a threatened bird species, endemic to the mixed conifer forests of the southeastern United States. RCW nests and forages preferentially in mature longleaf pine Pinus palustris, but also uses mature loblolly pine Pinus taeda and shortleaf pine Pinus echinata forests. In the last century, the extent of mature pine forests has been greatly reduced by logging. The RCW, in contrast to other woodpeckers, excavates nest cavities in living trees and senescence symptoms (year round leaf chlorosis and leaf mortality) have been observed in mature pine stands across the southeast. Widespread mortality of the mature pine forests would threaten the long-term survival of the RCW. We used airborne hyperspectral data across a portion of Ft. Benning Military Installation, Georgia, U.S.A., to determine if senescent trees can be identified and mapped and assess the likely persistence of mature pines in the RCW habitat. Univariate analysis of variance showed good separation between asymptomatic, senesced and dead physiological conditions with asymptomatic trees having significantly higher reflectance for all bands in the wavelength range between 0.719 and 1.1676 µm, senescent trees having significantly lower reflectance for bands in the range between 1.1927 and 1.3122 µm, and dead trees having significantly higher reflectance for bands in the range between 1.8151 and 1.9471 µm. Classification and Regression Tree (CART) models achieved correct classification rates and kappa statistics above 70%. CART models selected information from wavelength regions similar to those identified from the ANOVA, which likely explains their performance. Our aggregated CART model of tree senescence estimated that 141.4 ha (3%) of the study area is affected. RCW nests occurred in areas with significantly higher tree cover, and trees within foraging and home ranges did not have significantly more senescence than areas without RCW. These results indicate that although tree senescence is widespread, mortality is yet to significantly affect RCW habitat. Results and analysis of critical habitat similar to those exemplified in this study can extend our knowledge about the stressors of these important and imperiled components of biodiversity.     

Research on Crop Classification Method based on Sentinel-2 Time Series Combined Vegetation Index

Time series is a widely used phenological research method. A new time series vegetation indices which takes full advantage of the red edge information of Sentinel 2 data were used for crop classification to improve the classification accuracy. The NDVI, EVI, and red edge NDVI were combined to construct a time series vegetation index image. Then, four different algorithms (support vector machine, random forest, CART decision tree and maximum likelihood) were used to classify four crops, three forest grasses, bare land, and water bodies. Among the original classification results, the random forest with the highest overall accuracy is 87.92%, and the maximum likelihood with the lowest overall accuracy is 80.07%. In the classification details, the boundaries of random forest and support vector machine are the clearest. Among the four classification results, the classification accuracy of crops is higher than other land types, just smaller than water body. The error mainly comes from the mixture of three forests. It indicates that the time series combined vegetation index is feasible and accurate for crop classification.     

基于哨兵2时间序列组合植被指数的作物分类研究

时间序列是一种常用的物候研究方法。为充分利用哨兵2数据在红边波段的丰富信息,本文利用多种植被指数组合成时间序列进行作物分类。将NDVI、EVI、红边NDVI三种植被指数进行组合,构建时序植被指数图像,然后使用支持向量机、随机森林、CART决策树和最大似然4种不同的算法对四种作物、三种林草、裸露地表、水体进行分类。原始分类结果中,总体精度最高的随机森林为87.92%,最低的最大似然为80.07%,在分类细节上,随机森林和支持向量机的边界最清晰,4种分类结果中,农作物的分类精度均高于其他地类,仅次于水体的精度,误差主要来自三种林草的混分,表明时间序列组合植被指数用于农作物分类是可行的。     

A hyperspectral GA-PLSR model for prediction of pine wilt disease

Pine wilt disease caused by a forest-invasive alien species, the pine wood nematode (Bursaphelenchus xylophilus) is considered as one of the most destructive pest problems. In recent years, spectroscopic technologies have shown great potentials for the assessment of forest damage due to their nondestructive, noninvasive, cost-effective, and rapidly responsive nature. This paper first identified the hyperspectral characteristics of pine wilt disease by measuring and analyzing the changes in spectral reflectance of healthy and infected Pinus massoniana trees. Then 16 spectral features were extracted from the spectral bands covering the green region (510~580 nm), the red region (620~680 nm), the red edge (680~760 nm), the near-infrared region (780~1100 nm), and coded as genes composing the chromosome of a genetic algorithm (GA). Based on the optimal spectral features with suitable fitness from the GA, a partial least squares regression (PLSR) prediction model was built with highest determination coefficient R²_c?=?0.91, R²_v?=?0.82, relative prediction deviation RPD?=?3.3 and lowest root mean square error RMSE_c?=?0.23, RMSE_v?=?0.33 on the calibration and validation datasets. Compared with other PLSR models, our proposed GA-based approach significantly improves the prediction accuracy with few input spectral features.

     

Mapping health levels of Robinia pseudoacacia forests in the Yellow River delta,China, using IKONOS and Landsat 8 OLI imagery

The largest artificial Robinia pseudoacacia forests in the Yellow River delta of China have been infected by dieback diseases. Over the past several decades, this has caused a large amount of mortality of Robinia pseudoacacia forests in this area. Timely and accurate information on the health levels of the forests is crucial to improving local ecological and economic conditions. Remote sensing has been demonstrated to be a useful tool to map forest diseases over a large area. In this study, IKONOS and Landsat 8 Operational Land Imager (OLI) sensor data were collected for comparing their capability of accurately mapping health levels of the artificial forests. There were three health levels (i.e. healthy, medium dieback, and severe dieback) based on explicit tree crown symptoms. After the IKONOS and OLI images were preprocessed, both spatial and spectral features were extracted from the IKONOS and OLI imagery, and a maximum likelihood classification method was used to identify and map health levels of Robinia pseudoacacia forests. The experimental results indicate that the IKONOS sensor has greater potential for identifying and mapping forest health levels. Furthermore, texture features, especially texture variance, derived from the IKONOS panchromatic band, contributed greatly to the accuracy of classification results, achieving an overall accuracy (OA) of 96% for the IKONOS sensor and an OA of 88% for the OLI 2, which used both OLI spectral and IKONOS spatial features, compared with an OA of 74% for the OLI sensor alone. Our results indicate that the texture features extracted from high resolution imagery can improve the classification accuracy of health levels of planted forests with a regular spatial pattern. Our experimental results also demonstrate that classification of an image with a spatial resolution similar to, or finer than, tree crown diameter outperforms that of relatively coarse resolution imagery for differentiating living tree crowns and understorey dense green grass.     

Prediction and assessment of bark beetle-induced mortality of lodgepole pine using estimates of stand vigor derived from remotely sensed data

The current outbreak of mountain pine beetle (Dendroctonus ponderosae Hopkins) in British Columbia (BC), Canada, has led forest managers to consider thinning as a means of decreasing residual tree susceptibility to attack and subsequent mortality. Previous research indicates that susceptibility to mountain pine beetle is a function of a tree's physiological vigor and the intensity of attack. Trees able to produce ≥ 80 g (g) of wood per m² of projected leaf area annually are highly resistant, because they are able to shift resource allocation locally from wood to resin production to isolate blue-stain fungi introduced by attacking beetles. Typically, the leaf area of susceptible stands must be reduced by two-thirds to permit most residual trees to increase their vigor to a safe level. We evaluate whether Landsat Thematic Mapper (TM) imagery (30 × 30 m) provides a means to assess the maximum leaf area index (LAI) of unthinned stands and the extent that thinning reduces LAI. The extent that residual trees in thinned stands may have increased their resistance to attack from mountain pine beetle is predicted from a non-linear relationship between % maximum LAI and mean tree vigor.We investigated the merits of this approach in the vicinity of Parson, British Columbia using four stands of lodgepole pine (Pinus contorta Dougl.), two of which were heavily thinned (stands were spaced to 4 and 5 m, approximately 70% reduction in stand density). An analysis of archived Landsat TM imagery indicated that prior to thinning in 1993, all four stands had full canopy, which, for mature stands, would translate to mean tree vigor between 40 and 70 g of annual wood production per m² of foliage. By 1995, based on estimated changes in LAI derived from a second data of Landsat TM imagery, stand vigor in the unthinned stands had not changed; however, in the thinned stands, a nearly two third reduction in LAI resulted in a predicted increase in vigor to between 100 and 160 g wood m^− 2 of leaf area. A subsequent assessment in 2001 indicated that stand vigor remained higher in the thinned stands relative to the control stands. Following an infestation of mountain pine beetle in the study area in 2002, mortality data indicated that the thinned stands experienced no mortality relative to the unthinned stands which experienced 5.5% mortality in the initial years of the attack. In the larger area surrounding the study site, a general relationship was found between predicted stand vigor and mountain pine beetle-induced mortality as estimated from aerial overview survey data (r² = 0.43, p < 0.01).     

The effects of bark beetle stress on the foliar spectral reflectance of lodgepole pine

Abstract

Abstract. This paper reports the results of a study of the foliar reflectance properties of lodgepole pine trees under mountain pine beetle attack which aimed to determine the sequence of changes which occur. The motivation for this study was to identify spectral regions showing the earliest signs of attack in order to choose spectral bands for airborne linear array sensors. Normal illumination, diffuse hemispherical reflectance spectra were obtained of needles from trees showing varying degrees of stress from beetle attack and for needles from unattacked trees for comparison. The more pronounced changes in the spectra were interpreted visually and compared to changes reported by other authors. The more subtle changes were studied by analysing variance methods. Three spectral bands (the green peak, red edge and near-infrared (NIR) shoulder regions) have been identified as most promising for detecting early effects of bark beetle attack. Three additional bands (the blue, red and NIR plateau regions) are identified as reference bands for calibration and comparison purposes. The difference between foliar reflectance of attacked and unattacked trees is found to be most significant in the visible and red edge regions for current foliage and in the NIR for previous foliage. A red edge red shift was observed in the spectra of current foliage from attacked trees, in contrast to red edge blue shifts associated with stress in other studies. The observed sequence of subtle to more pronounced changes cannot be explained qualitatively using current knowledge of the plant pigment and anatomical changes which occur at the cellular and needle levels in stressed conifers. Further study of the detailed changes in pigments anel cellular and foliar anatomy is recommended, both to elucidate on the cause-and-effect relationships which occur, and to indicate the extent to which this paper's findings can be generalized. The findings suggest that multispectral linear array airborne scanners may be able to detect stressed conifers long before the red attack stage, but further investigation is required to determine whether the differences at foliar level between attacked and unattacked trees can also be detected at the whole-tree level, and whether confusion with other ground cover types can be avoided.     

Annual seasonality in Sentinel-1 signal for forest mapping and forest type classification

ABSTRACT

The Sentinel-1 satellites provide the formerly unprecedented combination of high spatial and temporal resolution of dual polarization synthetic aperture radar data. The availability of dense time series enables the derivation and analysis of temporally filtered annual backscatter signals. The study concentrates on the use of Sentinel-1 seasonal backscatter signatures for forest area estimation and forest type classification. A classification method based on time series similarity measures is introduced and tested in three test areas covered by various forest types including broadleaf temperate, boreal and montane forests. The results are compared with two European-wide Copernicus high resolution layers, namely forest type and tree cover density (TCD). The correspondence of forest/non-forest maps and TCD is high in all test areas, with overall accuracies for forest/non-forest classification between 86% and 91% and Pearson correlation coefficients for TCD between 0.68 and 0.74. The forest type classification (non-forest, coniferous and broadleaf forest classes) provides best results in temperate forests with an overall accuracy of 85%; in boreal forest, the accuracy decreases to only 65%. Generally, the method provides reliable results for forest area estimation, including regions where methods based on static parameters are often problematic (mountainous areas), and it enables forest type classification in temperate forests.     

Estimation of insect infestation dynamics using a temporal sequence of Landsat data

The current outbreak of mountain pine beetle (Dendroctonus ponderosae) in western Canada has been increasing over the past decade and is currently estimated to be impacting 9.2 million hectares, with varying levels of severity. Large area insect monitoring is typically undertaken using manual aerial overview sketch mapping, whereby an interpreter depicts areas of homogenous insect attack conditions onto 1:250,000 or 1:100,000 scale maps. These surveys provide valuable strategic data for management at the provincial scale. The coarse spatial and attribute resolution of these data however, make them inappropriate for fine-scale monitoring and operational planning. For instance, it is not possible to estimate the initial timing of attack and year of stand death. In this study, we utilise eight Landsat scenes collected over a 14 year period in north-central British Columbia, Canada, where the infestation has gradually developed both spatially and temporally. After pre-processing and normalising the eight scenes using a relative normalisation procedure, decision tree analysis was applied to classify spectral trajectories of the Normalised Difference Moisture Index (NDMI). From the classified temporal sequence of images, key parameters were extracted including the presence of beetle disturbance and timing of stand decline. The accuracy of discriminating beetle attack from healthy forest stands was assessed both spatially and temporally using three years of aerial survey data (1996, 2003, and 2004) with results indicating overall classification accuracies varying between 71 and 86%. As expected, the earliest and least severe attack year (1996), recorded the lowest overall accuracy. The relationship between the timing of stand attack (i.e. moderate to severe beetle infestation) and NDMI (initial year of detected disturbance) was also explored. The results suggest that there is potential for deriving regional estimates of the year of stand death using Landsat data and decision tree analysis however, a higher temporal frequency of images is required to quantify the timing of mountain pine beetle attack.     

Mapping insect defoliation in Scots pine with MODIS time-series data

Insect damage is a general problem that disturbs the growth of forests, causing economic losses and affecting carbon sequestration. Coarse-resolution data from satellites are potentially useful for national and regional mapping of forest damage, but the accuracy of these methods has not been fully examined. In this study, a method was tested for the mapping of defoliation in Scots pine [Pinus silvestris] forests in southeast Norway caused by the pine sawfly [Neodiprion sertifer], with the use of multi-temporal MODIS 16-day composite vegetation index data and the TIMESAT processing method. The damage mapping method used differences in summer mean values and angles of the seasonal profiles, indicating decreasing foliage density, to identify pixels that represent areas containing forest damage. In addition to 16-day NDVI the Wide Dynamic Range Vegetation Index (WDRVI) was tested. Damage areas were identified by classifying data into pixels representing damaged versus undamaged forest areas using a boolean combination of thresholded parameters. Classification results were evaluated against the change in LAI estimated from airplane LIDAR measurements, as an indicator of defoliation. The damage classifications detected 71% to 82% of the pixels with damage, and had kappa coefficients varying between 0.48 and 0.63, indicating some overestimation. This was due e.g. to failure to include clear-cut areas in the evaluation data. Damage classification with WDRVI only resulted in slight improvement compared to the NDVI. Only weak relationships were found between the LIDAR-estimated defoliation and the change parameters obtained from MODIS. Consequently, mapping of the degree of defoliation from MODIS was abandoned. In conclusion, the damage detection method based on MODIS data was found to be useful for locating insect damage, but not for estimating its intensity. Control of the detected damage areas using high-resolution remote sensing data, aerial survey, or fieldwork is recommended for accurate delineation in operational applications.     

Examining pine spectral separability using hyperspectral data from an airborne sensor: An extension of field‐based results

Three southern USA forestry species, loblolly pine (Pinus taeda), Virginia pine (Pinus virginiana), and shortleaf pine (Pinus echinata), were previously shown to be spectrally separable (83% accuracy) using data from a full‐range spectroradiometer (400–2500 nm) acquired above tree canopies. This study focused on whether these same species are also separable using hyperspectral data acquired using the airborne visible/infrared imaging spectrometer (AVIRIS). Stepwise discriminant techniques were used to reduce data dimensionality to a maximum of 10 spectral bands, followed by discriminant techniques to measure separability. Discriminatory variables were largely located in the visible and near‐infrared regions of the spectrum. Cross‐validation accuracies ranged from 65% (1 pixel radiance data) to as high as 85% (3×3 pixel radiance data), indicating that these species have strong potential to be classified accurately using hyperspectral data from air‐ or space‐borne sensors.     

Linear mixture model applied to Amazonian vegetation classification

Many research projects require accurate delineation of different secondary succession (SS) stages over large regions/subregions of the Amazon basin. However, the complexity of vegetation stand structure, abundant vegetation species, and the smooth transition between different SS stages make vegetation classification difficult when using traditional approaches such as the maximum likelihood classifier (MLC). Most of the time, classification distinguishes only between forest/non-forest. It has been difficult to accurately distinguish stages of SS. In this paper, a linear mixture model (LMM) approach is applied to classify successional and mature forests using Thematic Mapper (TM) imagery in the Rondônia region of the Brazilian Amazon. Three endmembers (i.e., shade, soil, and green vegetation or GV) were identified based on the image itself and a constrained least-squares solution was used to unmix the image. This study indicates that the LMM approach is a promising method for distinguishing successional and mature forests in the Amazon basin using TM data. It improved vegetation classification accuracy over that of the MLC. Initial, intermediate, and advanced successional and mature forests were classified with overall accuracy of 78.2% using a threshold method on the ratio of shade to GV fractions, a 7.4% increase over the MLC. The GV and shade fractions are sensitive to the change of vegetation stand structures and better capture biophysical structure information.