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.     

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.     

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.     

Detection of foliage conditions and disturbance from multi-angular high spectral resolution remote sensing

Disturbance of forest ecosystems, an important component of the terrestrial carbon cycle, has become a focus of research over recent years, as global warming is about to increase the frequency and severity of natural disturbance events. Remote sensing offers unique opportunities for detection of forest disturbance at multiple scales; however, spatially and temporally continuous mapping of non-stand replacing disturbance remains challenging. First, most high spatial resolution satellite sensors have relatively broad spectral ranges with bandwidths unsuitable for detection of subtle, stress induced, features in canopy reflectance. Second, directional and background reflectance effects, induced by the interactions between the sun-sensor geometry and the observed canopy surface, make up-scaling of empirically derived relationships between changes in spectral reflectance and vegetation conditions difficult. Using an automated tower based spectroradiometer, we analyse the interactions between canopy level reflectance and different stages of disturbance occurring in a mountain pine beetle infested lodgepole pine stand in northern interior British Columbia, Canada, during the 2007 growing season. Directional reflectance effects were modelled using a bidirectional reflectance distribution function (BRDF) acquired from high frequency multi-angular spectral observations. Key wavebands for observing changes in directionally corrected canopy spectra were identified using discriminant analysis and highly significant correlations between canopy reflectance and field measured disturbance levels were found for several broad and narrow waveband vegetation indices (for instance, r²_NDVI = 0.90; r²_CHL3 = 0.85; p < 0.05). Results indicate that multi-angular observations are useful for extraction of disturbance related changes in canopy reflectance, in particular the temporally and spectrally dense data detected changes in chlorophyll content well. This study will help guide and inform future efforts to map forest health conditions at landscape and over increasingly coarse scales.     

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.     

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).     

A Landsat time series approach to characterize bark beetle and defoliator impacts on tree mortality and surface fuels in conifer forests

Insects are important forest disturbance agents, and mapping their effects on tree mortality and surface fuels represents a critical research challenge. Although various remote sensing approaches have been developed to monitor insect impacts, most studies have focused on single insect agents or single locations and have not related observed changes to ground-based measurements. This study presents a remote sensing framework to (1) characterize spectral trajectories associated with insect activity of varying duration and severity and (2) relate those trajectories to ground-based measurements of tree mortality and surface fuels in the Cascade Range, Oregon, USA. We leverage a Landsat time series change detection algorithm (LandTrendr), annual forest health aerial detection surveys (ADS), and field measurements to investigate two study landscapes broadly applicable to conifer forests and dominant insect agents of western North America. We distributed 38 plots across multiple forest types (ranging from mesic mixed-conifer to xeric lodgepole pine) and insect agents (defoliator [western spruce budworm] and bark beetle [mountain pine beetle]). Insect effects were evident in the Landsat time series as combinations of both short- and long-duration changes in the Normalized Burn Ratio spectral index. Western spruce budworm trajectories appeared to show a consistent temporal evolution of long-duration spectral decline (loss of vegetation) followed by recovery, whereas mountain pine beetle plots exhibited both short- and long-duration spectral declines and variable recovery rates. Although temporally variable, insect-affected stands generally conformed to four spectral trajectories: short-duration decline then recovery, short- then long-duration decline, long-duration decline, long-duration decline then recovery. When comparing remote sensing data with field measurements of insect impacts, we found that spectral changes were related to cover-based estimates (tree basal area mortality [R²_adj = 0.40, F_1,34 = 24.76, P < 0.0001] and down coarse woody detritus [R²_adj = 0.29, F_1,32 = 14.72, P = 0.0006]). In contrast, ADS changes were related to count-based estimates (e.g., ADS mortality from mountain pine beetle positively correlated with ground-based counts [R²_adj = 0.37, F_1,22 = 14.71, P = 0.0009]). Fine woody detritus and forest floor depth were not well correlated with Landsat- or aerial survey-based change metrics. By characterizing several distinct temporal manifestations of insect activity in conifer forests, this study demonstrates the utility of insect mapping methods that capture a wide range of spectral trajectories. This study also confirms the key role that satellite imagery can play in understanding the interactions among insects, fuels, and wildfire.     

Detection of mountain pine beetle-killed ponderosa pine in a heterogeneous landscape using high-resolution aerial imagery

Previous studies have used remote-sensing images to map tree mortality caused by mountain pine beetle (Dendroctonus ponderosae; MPB) in relatively homogeneous lodgepole pine (Pinus contorta) forests; however, classification methods have not been tested for the patchy landscape of ponderosa pine-dominated (Pinus ponderosae) montane forests characterized by highly variable tree density. This study explores two supervised classification methods to identify MPB-caused mortality (red attack) in heterogeneous montane forests of the Colorado Front Range using 1 m-resolution 2011 imagery of the National Agriculture Imagery Program (NAIP): maximum likelihood using the red, green, and blue bands, and the red-green index (RGI), and a thresholding technique using the RGI. Two variations of the RGI threshold method were also explored: the addition of a green-band threshold and the incorporation of a focal analysis. Evaluation pixels were used to assess the accuracy of the classification methods. The maximum likelihood (97 Percentage Correctly Classified (PCC); 11% error of commission for red attack) and RGI threshold (85 PCC; 46% error of commission for red attack) classification methods overestimated the red attack. The RGI and green band threshold classification reduced the error of commission (5%) and had high overall accuracy (97 PCC). In a comparison of classification methods across tree-density sites, we found the maximum likelihood classification had a very high accuracy in the high-density site (95 PCC), but substantially lower accuracy in the low-density site (85 PCC) due to the presence of more visible cover types. The RGI threshold classification with the green band constraint produced more consistent PCCs across tree densities: high (93.7 PCC), moderate (95.2 PCC), and low (92.0 PCC). Our results indicate forest structure may affect the classification accuracy and should be considered when selecting a classification method for a landscape.     

Modeling mountain pine beetle infestation with an agent-based approach at two spatial scales

Extensive outbreaks of tree-killing insects have been occurring in many parts of North America, including the province of British Columbia, raising concerns about the health of pine forest ecosystems. The dynamic phenomenon of mountain pine beetle (MPB), Dendroctonus ponderosae Hopkins, infestation outbreaks is an inherent spatial and temporal complex process. Agent-based modeling (ABM) facilitates simulating spatial interactions that describe the ecological context in which insect populations spread. The main objective of this study was to develop a model of the MPB forest infestation dynamics. This spatially explicit model integrates geographic information systems (GISs) and ABM to simulate MPB outbreaks at the tree and landscape scales, providing spatiotemporal information of annual distribution and patterns of MPB outbreaks. This prototype was implemented with geographic data generated from aerial overview surveys carried out by the B.C. Ministry of Forests and Range, for the study site in Kamloops, Canada. Results show the direct influence that vigorous forest stands and trees have on higher breeding rates, and therefore in the MPB population increment at a tree scale, in a period of 5 years. The simulation results at the landscape level help to determine the most probable locations of future MPB infestations in a time frame of 10 years.     

A comparison of multi-spectral, multi-angular, and multi-temporal remote sensing datasets for fractional shrub canopy mapping in Arctic Alaska

Shrub cover appears to be increasing across many areas of the Arctic tundra biome, and increasing shrub cover in the Arctic has the potential to significantly impact global carbon budgets and the global climate system. For most of the Arctic, however, there is no existing baseline inventory of shrub canopy cover, as existing maps of Arctic vegetation provide little information about the density of shrub cover at a moderate spatial resolution across the region. Remotely-sensed fractional shrub canopy maps can provide this necessary baseline inventory of shrub cover. In this study, we compare the accuracy of fractional shrub canopy (> 0.5 m tall) maps derived from multi-spectral, multi-angular, and multi-temporal datasets from Landsat imagery at 30 m spatial resolution, Moderate Resolution Imaging SpectroRadiometer (MODIS) imagery at 250 m and 500 m spatial resolution, and MultiAngle Imaging Spectroradiometer (MISR) imagery at 275 m spatial resolution for a 1067 km² study area in Arctic Alaska. The study area is centered at 69 °N, ranges in elevation from 130 to 770 m, is composed primarily of rolling topography with gentle slopes less than 10°, and is free of glaciers and perennial snow cover. Shrubs > 0.5 m in height cover 2.9% of the study area and are primarily confined to patches associated with specific landscape features. Reference fractional shrub canopy is determined from in situ shrub canopy measurements and a high spatial resolution IKONOS image swath. Regression tree models are constructed to estimate fractional canopy cover at 250 m using different combinations of input data from Landsat, MODIS, and MISR. Results indicate that multi-spectral data provide substantially more accurate estimates of fractional shrub canopy cover than multi-angular or multi-temporal data. Higher spatial resolution datasets also provide more accurate estimates of fractional shrub canopy cover (aggregated to moderate spatial resolutions) than lower spatial resolution datasets, an expected result for a study area where most shrub cover is concentrated in narrow patches associated with rivers, drainages, and slopes. Including the middle infrared bands available from Landsat and MODIS in the regression tree models (in addition to the four standard visible and near-infrared spectral bands) typically results in a slight boost in accuracy. Including the multi-angular red band data available from MISR in the regression tree models, however, typically boosts accuracy more substantially, resulting in moderate resolution fractional shrub canopy estimates approaching the accuracy of estimates derived from the much higher spatial resolution Landsat sensor. Given the poor availability of snow and cloud-free Landsat scenes in many areas of the Arctic and the promising results demonstrated here by the MISR sensor, MISR may be the best choice for large area fractional shrub canopy mapping in the Alaskan Arctic for the period 2000-2009.     

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.     

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.     

Sensitivity of the thematic mapper enhanced wetness difference index to detect mountain pine beetle red-attack damage

Red-attack damage caused by mountain pine beetle (Dentroctonus ponderosa Hopkins) infestation in stands of lodgepole pine (Pinus contorta) in the Prince George Forest Region of British Columbia was examined using multitemporal Landsat-7 ETM+ imagery acquired in 1999, 2000, and 2001. The image data were geometrically and atmospherically corrected, and processed using the Tasseled Cap Transformation (TCT) to obtain wetness indices. The final steps included pixel subtraction, enhancement, and thresholding of the wetness index differences. The resulting enhanced wetness difference index (EWDI) was used to interpret spectral patterns in stands with confirmed (through aerial survey) red-attack damage in 2001, and these EWDI patterns were compared to the patterns of reflectance in normal-colour composites. We stratified the aerial survey dataset into two levels and used the EWDI to discriminate classes of 10-29 red-attack trees and 30-50 red-attack trees, and a sample of healthy forest collected from inventory data. Classification accuracy of red-attack damage based on the EWDI ranged from 67% to 78% correct.     

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.     

Assessing the potential of VEGETATION sensor data for mapping snow and ice cover: A Normalized Difference Snow and Ice Index

The VEGETATION (VGT) sensor in SPOT 4 has four spectral bands that are equivalent to Landsat Thematic Mapper (TM) bands (blue, red, near-infrared and mid-infrared spectral bands) and provides daily images of the global land surface at a 1-km spatial resolution. We propose a new index for identifying and mapping of snow/ice cover, namely the Normalized Difference Snow/Ice Index (NDSII), which uses reflectance values of red and mid-infrared spectral bands of Landsat TM and VGT. For Landsat TM data, NDSII is calculated as NDSII_TM=(TM3-TM5)/(TM3+TM5); for VGT data, NDSII is calculated as NDSII_VGT=(B2-MIR)/(B2+MIR). As a case study we used a Landsat TM image that covers the eastern part of the Qilian mountain range in the Qinghai-Xizang (Tibetan) plateau of China. NDSIITM gave similar estimates of the area and spatial distribution of snow/ice cover to the Normalized Difference Snow Index (NDSI=(TM2-TM5)/(TM2+TM5)) which has been proposed by Hall et al. The results indicated that the VGT sensor might have the potential for operational monitoring and mapping of snow/ice cover from regional to global scales, when using NDSII_VGT.     

Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model

A new method is described for the retrieval of fractional cover of large woody plants (shrubs) at the landscape scale using moderate resolution multi-angle remote sensing data from the Multiangle Imaging SpectroRadiometer (MISR) and a hybrid geometric-optical (GO) canopy reflectance model. Remote sensing from space is the only feasible method for regularly mapping woody shrub cover over large areas, an important application because extensive woody shrub encroachment into former grasslands has been seen in arid and semi-arid grasslands around the world during the last 150 years. The major difficulty in applying GO models in desert grasslands is the spatially dynamic nature of the combined soil and understory background reflectance: the background is important and cannot be modeled as either a Lambertian scatterer or by using a fixed bidirectional reflectance distribution function (BRDF). Candidate predictors of the background BRDF at the Sun-target-MISR angular sampling configurations included the volume scattering kernel weight from a Li-Ross BRDF model; diffuse brightness (ρ0) from the Modified Rahman-Pinty-Verstraete (MRPV) BRDF model; other Li-Ross kernel weights (isotropic, geometric); and MISR near-nadir bidirectional reflectance factors (BRFs) in the blue, green, and near infra-red bands. The best method was multiple regression on the weights of a kernel-driven model and MISR nadir camera blue, green, and near infra-red bidirectional reflectance factors. The results of forward modeling BRFs for a 5.25 km² area in the USDA, ARS Jornada Experimental Range using the Simple Geometric Model (SGM) with this background showed good agreement with the MISR data in both shape and magnitude, with only minor spatial discrepancies. The simulations were shown to be accurate in terms of both absolute value and reflectance anisotropy over all 9 MISR views and for a wide range of canopy configurations (r² = 0.78, RMSE = 0.013, N = 3969). Inversion of the SGM allowed estimation of fractional shrub cover with a root mean square error (RMSE) of 0.03 but a relatively weak correlation (r² = 0.19) with the reference data (shrub cover estimated from high resolution IKONOS panchromatic imagery). The map of retrieved fractional shrub cover was an approximate spatial match to the reference map. Deviations reflect the first-order approximation of the understory BRDF in the MISR viewing plane; errors in the shrub statistics; and the 12 month lag between the two data sets.     

Evaluating MODIS data for mapping wildlife habitat distribution

Habitat distribution models have a long history in ecological research. With the development of geospatial information technology, including remote sensing, these models are now applied to an ever-increasing number of species, particularly those located in areas in which it is logistically difficult to collect habitat data in the field. Many habitat studies have used data acquired by multi-spectral sensor systems such as the Landsat Thematic Mapper (TM), due mostly to their availability and relatively high spatial resolution (30 m/pixel). The use of data collected by other sensor systems with lower spatial resolutions but high frequency of acquisitions has largely been neglected, due to the perception that such low spatial resolution data are too coarse for habitat mapping. In this study we compare two models using data from different satellite sensor systems for mapping the spatial distribution of giant panda habitat in Wolong Nature Reserve, China. The first one is a four-category scheme model based on combining forest cover (derived from a digital land cover classification of Landsat TM imagery acquired in June, 2001) with information on elevation and slope (derived from a digital elevation model obtained from topographic maps of the study area). The second model is based on the Ecological Niche Factor Analysis (ENFA) of a time series of weekly composites of WDRVI (Wide Dynamic Range Vegetation Index) images derived from MODIS (Moderate Resolution Imaging Spectroradiometer – 250 m/pixel) for 2001. A series of field plots was established in the reserve during the summer–autumn months of 2001–2003. The locations of the plots with panda feces were used to calibrate the ENFA model and to validate the results of both models. Results showed that the model using the seasonal variability of MODIS-WDRVI had a similar prediction success to that using Landsat TM and digital elevation model data, albeit having a coarser spatial resolution. This suggests that the phenological characterization of the land surface provides an appropriate environmental predictor for giant panda habitat mapping. Therefore, the information contained in remotely sensed data acquired with low spatial resolution but high frequency of acquisitions has considerable potential for mapping the habitat distribution of wildlife species.     

Super-resolution land cover mapping with indicator geostatistics

Many satellite images have a coarser spatial resolution than the extent of land cover patterns on the ground, leading to mixed pixels whose composite spectral response consists of responses from multiple land cover classes. Spectral unmixing procedures only determine the fractions of such classes within a coarse pixel without locating them in space. Super-resolution or sub-pixel mapping aims at providing a fine resolution map of class labels, one that displays realistic spatial structure (without artifact discontinuities) and reproduces the coarse resolution fractions. In this paper, existing approaches for super-resolution mapping are placed within an inverse problem framework, and a geostatistical method is proposed for generating alternative synthetic land cover maps at the fine (target) spatial resolution; these super-resolution realizations are consistent with all the information available.More precisely, indicator coKriging is used to approximate the probability that a pixel at the fine spatial resolution belongs to a particular class, given the coarse resolution fractions and (if available) a sparse set of class labels at some informed fine pixels. Such Kriging-derived probabilities are used in sequential indicator simulation to generate synthetic maps of class labels at the fine resolution pixels. This non-iterative and fast simulation procedure yields alternative super-resolution land cover maps that reproduce: (i) the observed coarse fractions, (ii) the fine resolution class labels that might be available, and (iii) the prior structural information encapsulated in a set of indicator variogram models at the fine resolution. A case study is provided to illustrate the proposed methodology using Landsat TM data from SE China.     

Forest mapping by geoinformatics for landscape fire behaviour modelling in coastal forests,Greece

This study aims at quantifying and mapping fire-related characteristics of forest structure through field inventories, statistics, remote sensing, and geographical information systems in the island of Lesvos, northeast Aegean Sea, Greece. Simulation of fire behaviour requires forest biomass inputs that describe surface fuel types/models along with canopy fuel properties, such as canopy cover, stand height, crown base height, and crown bulk density, to accurately predict surface and crown fire spread and spotting potential. Forest canopy characteristics and other vegetation attributes were sampled and derived in over 100 field plots, the majority of which were located in coastal pine forest stands. Regression models involving four dependent forest stand variables (stand height, canopy cover, crown base height, and crown bulk density) were developed using generalized additive models. The values of adjusted R ² were 0.72 for stand height, 0.68 for canopy cover, 0.51 for crown base height, and 0.33 for crown bulk density. These regression models were used to create forest fuel characteristics layers, which can be used as inputs to fire management applications and state-of-the-art landscape-scale fire behaviour models.