Potential of VIIRS Time Series Data for Aiding the USDA Forest Service Early Warning System for Forest Health Threats: A Gypsy Moth Defoliation Case Study

This report details one of three experiments performed during FY 2007 for the NASA RPC (Rapid Prototyping Capability) at Stennis Space Center. This RPC experiment assesses the potential of VIIRS (Visible/Infrared Imager/Radiometer Suite) and MODIS (Moderate Resolution Imaging Spectroradiometer) data for detecting and monitoring forest defoliation from the non-native Eurasian gypsy moth (Lymantria dispar). The intent of the RPC experiment was to assess the degree to which VIIRS data can provide forest disturbance monitoring information as an input to a forest threat EWS (Early Warning System) as compared to the level of information that can be obtained from MODIS data. The USDA Forest Service (USFS) plans to use MODIS products for generating broad-scaled, regional monitoring products as input to an EWS for forest health threat assessment. NASA SSC is helping the USFS to evaluate and integrate currently available satellite remote sensing technologies and data products for the EWS, including the use of MODIS products for regional monitoring of forest disturbance. Gypsy moth defoliation of the mid-Appalachian highland region was selected as a case study. Gypsy moth is one of eight major forest insect threats listed in the Healthy Forest Restoration Act (HFRA) of 2003; the gypsy moth threatens eastern U.S. hardwood forests, which are also a concern highlighted in the HFRA of 2003. This region was selected for the project because extensive gypsy moth defoliation occurred there over multiple years during the MODIS operational period. This RPC experiment is relevant to several nationally important mapping applications, including agricultural efficiency, coastal management, ecological forecasting, disaster management, and carbon management. In this experiment, MODIS data and VIIRS data simulated from MODIS were assessed for their ability to contribute broad, regional geospatial information on gypsy moth defoliation. Landsat and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data were used to assess the quality of gypsy moth defoliation mapping products derived from MODIS data and from simulated VIIRS data. The project focused on use of data from MODIS Terra as opposed to MODIS Aqua mainly because only MODIS Terra data was collected during 2000 and 2001-years with comparatively high amounts of gypsy moth defoliation within the study area. The project assessed the quality of VIIRS data simulation products. Hyperion data was employed to assess the quality of MODIS-based VIIRS simulation datasets using image correlation analysis techniques. The ART (Application Research Toolbox) software was used for data simulation. Correlation analysis between MODIS-simulated VIIRS data and Hyperion-simulated VIIRS data for red, NIR (near-infrared), and NDVI (Normalized Difference Vegetation Index) image data products collectively indicate that useful, effective VIIRS simulations can be produced using Hyperion and MODIS data sources. The r(exp 2) for red, NIR, and NDVI products were 0.56, 0.63, and 0.62, respectively, indicating a moderately high correlation between the 2 data sources. Temporal decorrelation from different data acquisition times and image misregistration may have lowered correlation results. The RPC experiment also generated MODIS-based time series data products using the TSPT (Time Series Product Tool) software. Time series of simulated VIIRS NDVI products were produced at approximately 400-meter resolution GSD (Ground Sampling Distance) at nadir for comparison to MODIS NDVI products at either 250- or 500-meter GSD. The project also computed MODIS (MOD02) NDMI (Normalized Difference Moisture Index) products at 500-meter GSD for comparison to NDVI-based products. For each year during 2000-2006, MODIS and VIIRS (simulated from MOD02) time series were computed during the peak gypsy moth defoliation time frame in the study area (approximately June 10 through July 27). Gypsy moth defoliation mapping products from simated VIIRS and MOD02 time series were produced using multiple methods, including image classification and change detection via image differencing. The latter enabled an automated defoliation detection product computed using percent change in maximum NDVI for a peak defoliation period during 2001 compared to maximum NDVI across the entire 2000-2006 time frame. Final gypsy moth defoliation mapping products were assessed for accuracy using randomly sampled locations found on available geospatial reference data (Landsat and ASTER data in conjunction with defoliation map data from the USFS). Extensive gypsy moth defoliation patches were evident on screen displays of multitemporal color composites derived from MODIS data and from simulated VIIRS vegetation index data. Such defoliation was particularly evident for 2001, although widespread denuded forests were also seen for 2000 and 2003. These visualizations were validated using aforementioned reference data. Defoliation patches were visible on displays of MODIS-based NDVI and NDMI data. The viewing of apparent defoliation patches on all of these products necessitated adoption of a specialized temporal data processing method (e.g., maximum NDVI during the peak defoliation time frame). The frequency of cloud cover necessitated this approach. Multitemporal simulated VIIRS and MODIS Terra data both produced effective general classifications of defoliated forest versus other land cover. For 2001, the MOD02-simulated VIIRS 400-meter NDVI classification produced a similar yet slightly lower overall accuracy (87.28 percent with 0.72 Kappa) than the MOD02 250-meter NDVI classification (88.44 percent with 0.75 Kappa). The MOD13 250-meter NDVI classification had a lower overall accuracy (79.13 percent) and a much lower Kappa (0.46). The report discusses accuracy assessment results in much more detail, comparing overall classification and individual class accuracy statistics for simulated VIIRS 400-meter NDVI, MOD02 250-meter NDVI, MOD02-500 meter NDVI, MOD13 250-meter NDVI, and MOD02 500-meter NDMI classifications. Automated defoliation detection products from simulated VIIRS and MOD02 data for 2001 also yielded similar, relatively high overall classification accuracy (85.55 percent for the VIIRS 400-meter NDVI versus 87.28 percent for the MOD02 250-meter NDVI). In contrast, the USFS aerial sketch map of gypsy moth defoliation showed a lower overall classification accuracy at 73.64 percent. The overall classification Kappa values were also similar for the VIIRS (approximately 0.67 Kappa) versus the MOD02 (approximately 0.72 Kappa) automated defoliation detection product, which were much higher than the values exhibited by the USFS sketch map product (overall Kappa of approximately 0.47). The report provides additional details on the accuracy of automated gypsy moth defoliation detection products compared with USFS sketch maps. The results suggest that VIIRS data can be effectively simulated from MODIS data and that VIIRS data will produce gypsy moth defoliation mapping products that are similar to MODIS-based products. The results of the RPC experiment indicate that VIIRS and MODIS data products have good potential for integration into the forest threat EWS. The accuracy assessment was performed only for 2001 because of time constraints and a relative scarcity of cloud-free Landsat and ASTER data for the peak defoliation period of the other years in the 2000-2006 time series. Additional work should be performed to assess the accuracy of gypsy moth defoliation detection products for additional years.The study area (mid-Appalachian highlands) and application (gypsy moth forest defoliation) are not necessarily representative of all forested regions and of all forest threat disturbance agents. Additional work should be performed on other inland and coastal regions as well as for other major forest threats

The earth observing one (EO-1) Hyperion and Advanced land imager sensors for use in tundra classification studies within the Upper Kuparuk river basin, Alaska

The heterogeneity of Arctic vegetation can make land cover classification very difficult when using medium to small resolution imagery (Schneider et al., 2009; Muller et al., 1999). Using high radiometric and spatial resolution imagery, such as the SPOT 5 and IKONOS satellites, have helped arctic land cover classification accuracies rise into the 80 and 90 percentiles (Allard, 2003; Stine et al., 2010; Muller et al., 1999). However, those increases usually come at a high price. High resolution imagery is very expensive and can often add tens of thousands of dollars onto the cost of the research. The EO-1 satellite launched in 2002 carries two sensors that have high spectral and/or high spatial resolutions and can be an acceptable compromise between the resolution versus cost issues. The Hyperion is a hyperspectral sensor with the capability of collecting 242 spectral bands of information. The Advanced Land Imager (ALI) is an advanced multispectral sensor whose spatial resolution can be sharpened to 10 meters. This dissertation compares the accuracies of arctic land cover classifications produced by the Hyperion and ALI sensors to the classification accuracies produced by the Systeme Pour l' Observation de le Terre (SPOT), the Landsat Thematic Mapper (TM) and the Landsat Enhanced Thematic Mapper Plus (ETM+) sensors. Hyperion and ALI images from August 2004 were collected over the Upper Kuparuk River Basin, Alaska. Image processing included the stepwise discriminant analysis of pixels that were positively classified from coinciding ground control points, geometric and radiometric correction, and principle component analysis. Finally, stratified random sampling was used to perform accuracy assessments on satellite derived land cover classifications. Accuracy was estimated from an error matrix (confusion matrix) that provided the overall, producer's and user's accuracies. This research found that while the Hyperion sensor produced classification accuracies that were equivalent to the TM and ETM+ sensor (approximately 78%), the Hyperion could not obtain the accuracy of the SPOT 5 HRV sensor. However, the land cover classifications derived from the ALI sensor exceeded most classification accuracies derived from the TM and ETM+ sensors and were even comparable to most SPOT 5 HRV classifications (87%). With the deactivation of the Landsat series satellites, the monitoring of remote locations such as in the Arctic on an uninterrupted basis throughout the world is in jeopardy. The utilization of the Hyperion and ALI sensors are a way to keep that endeavor operational. By keeping the ALI sensor active at all times, uninterrupted observation of the entire Earth can be accomplished. Keeping the Hyperion sensor as a "tasked" sensor can provide scientists with additional imagery and options for their studies without overburdening storage issues

EO-1 Data Quality and Sensor Stability with Changing Orbital Precession at the End of a 16 Year Mission

The Earth Observing One (EO-1) satellite has completed 16 years of Earth observations in early 2017. What started as a technology mission to test various new advancements turned into a science and application mission that extended many years beyond the satellites planned life expectancy. EO-1s primary instruments are spectral imagers: Hyperion, the only civilian full spectrum spectrometer (430-2400 nm) in orbit; and the Advanced Land Imager (ALI), the prototype for Landsat-8s pushbroom imaging technology. Both Hyperion and ALI instruments have continued to perform well, but in February 2011 the satellite ran out of the fuel necessary to maintain orbit, which initiated a change in precession rate that led to increasingly earlier equatorial crossing times during its last five years. The change from EO-1s original orbit, when it was formation flying with Landsat-7 at a 10:01am equatorial overpass time, to earlier overpass times results in image acquisitions with increasing solar zenith angles (SZAs). In this study, we take several approaches to characterize data quality as SZAs increased. Our results show that for both EO-1 sensors, atmospherically corrected reflectance products are within 5 to 10 of mean pre-drift products. No marked trend in decreasing quality in ALI or Hyperion is apparent through 2016, and these data remain a high quality resource through the end of the mission

Estimation of leaf area index and its sunlit portion from DSCOVR EPIC data: theoretical basis

This paper presents the theoretical basis of the algorithm designed for the generation of leaf area index and diurnal course of its sunlit portion from NASA's Earth Polychromatic Imaging Camera (EPIC) onboard NOAA's Deep Space Climate Observatory (DSCOVR). The Look-up-Table (LUT) approach implemented in the MODIS operational LAI/FPAR algorithm is adopted. The LUT, which is the heart of the approach, has been significantly modified. First, its parameterization incorporates the canopy hot spot phenomenon and recent advances in the theory of canopy spectral invariants. This allows more accurate decoupling of the structural and radiometric components of the measured Bidirectional Reflectance Factor (BRF), improves scaling properties of the LUT and consequently simplifies adjustments of the algorithm for data spatial resolution and spectral band compositions. Second, the stochastic radiative transfer equations are used to generate the LUT for all biome types. The equations naturally account for radiative effects of the three-dimensional canopy structure on the BRF and allow for an accurate discrimination between sunlit and shaded leaf areas. Third, the LUT entries are measurable, i.e., they can be independently derived from both below canopy measurements of the transmitted and above canopy measurements of reflected radiation fields. This feature makes possible direct validation of the LUT, facilitates identification of its deficiencies and development of refinements. Analyses of field data on canopy structure and leaf optics collected at 18 sites in the Hyytiälä forest in southern boreal zone in Finland and hyperspectral images acquired by the EO-1 Hyperion sensor support the theoretical basis.Shared Services Center NAS

Spectral contribution of understory to forest reflectance in a boreal site : an analysis of EO-1 Hyperion data

Boreal forests exhibit strong seasonal dynamics in their reflectance spectra during the short, snow-free growing period. This short communication paper reports an analysis of the seasonality of boreal forest spectra from the end of snowmelt until the time of maximal leaf area. We apply a forest reflectance model (FRT) to estimate the seasonal contribution of understow vegetation to forest reflectance from a time series of three Earth Observing 1 (EO-1) Hyperion images acquired in May, June and July. The reflectance simulations are based on detailed seasonal series of leaf area index and understory spectra measurements carried out in ten stands at the Hyytiala Forestry Field Station in Finland. Our results show that the contribution of understory to boreal forest reflectance is high in the visible domain, but it drops at the red edge and stays relatively low and constant in near infrared (NIR). Throughout the growing season, the contribution of the understory remains approximately the same in the NIR domain, whereas larger changes can be observed in the visible domain. (C) 2015 The Authors. Published by Elsevier Inc.Peer reviewe

Two-level model for land degradation mapping on multispectral satellite imagery

Two-level model for land degradation mapping on multispectral satellite imagery of low and medium spatial resolution is shown. First level model applies several different thematic classifications of source multispectral images, e.g. vegetation change, soil erosion etc. Second level gives data fusion of specific thematic classifications of the first level into final thematic map to improve accuracy and reliability of wing to joint interpretation. Proposed model will be useful for information support systems to provide land management

Comparative analysis of Landsat TM, ETM+, OLI and EO-1 ALI satellite images at the Tisza-tó area, Hungary

L

Detecting the effects of hydrocarbon pollution in the Amazon forest using hyperspectral satellite images

The global demand for fossil energy is triggering oil exploration and production projects in remote areas of the world. During the last few decades hydrocarbon production has caused pollution in the Amazon forest inflicting considerable environmental impact. Until now it is not clear how hydrocarbon pollution affects the health of the tropical forest flora. During a field campaign in polluted and pristine forest, more than 1100 leaf samples were collected and analysed for biophysical and biochemical parameters. The results revealed that tropical forests exposed to hydrocarbon pollution show reduced levels of chlorophyll content, higher levels of foliar water content and leaf structural changes. In order to map this impact over wider geographical areas, vegetation indices were applied to hyperspectral Hyperion satellite imagery. Three vegetation indices (SR, NDVI and NDVI705) were found to be the most appropriate indices to detect the effects of petroleum pollution in the Amazon forest

Applications of Remote Sensing to Alien Invasive Plant Studies

Biological invasions can affect ecosystems across a wide spectrum of bioclimatic conditions. Therefore, it is often important to systematically monitor the spread of species over a broad region. Remote sensing has been an important tool for large-scale ecological studies in the past three decades, but it was not commonly used to study alien invasive plants until the mid 1990s. We synthesize previous research efforts on remote sensing of invasive plants from spatial, temporal and spectral perspectives. We also highlight a recently developed state-of-the-art image fusion technique that integrates passive and active energies concurrently collected by an imaging spectrometer and a scanning-waveform light detection and ranging (LiDAR) system, respectively. This approach provides a means to detect the structure and functional properties of invasive plants of different canopy levels. Finally, we summarize regional studies of biological invasions using remote sensing, discuss the limitations of remote sensing approaches, and highlight current research needs and future directions