Evaluating an Automated Approach for Monitoring Forest Disturbances in the Pacific Northwest from Logging, Fire and Insect Outbreaks with Landsat Time Series Data

Forests are the largest aboveground sink for atmospheric carbon (C), and understanding how they change through time is critical to reduce our C-cycle uncertainties. We investigated a strong decline in Normalized Difference Vegetation Index (NDVI) from 1982 to 1991 in Pacific Northwest forests, observed with the National Ocean and Atmospheric Administration’s (NOAA) series of Advanced Very High Resolution Radiometers (AVHRRs). To understand the causal factors of this decline, we evaluated an automated classification method developed for Landsat time series stacks (LTSS) to map forest change. This method included: (1) multiple disturbance index thresholds; and (2) a spectral trajectory-based image analysis with multiple confidence thresholds. We produced 48 maps and verified their accuracy with air photos, monitoring trends in burn severity data and insect aerial detection survey data. Area-based accuracy estimates for change in forest cover resulted in producer’s and user’s accuracies of 0.21 ± 0.06 to 0.38 ± 0.05 for insect disturbance, 0.23 ± 0.07 to 1 ± 0 for burned area and 0.74 ± 0.03 to 0.76 ± 0.03 for logging. We believe that accuracy was low for insect disturbance because air photo reference data were temporally sparse, hence missing some outbreaks, and the annual anniversary time step is not dense enough to track defoliation and progressive stand mortality. Producer’s and user’s accuracy for burned area was low due to the temporally abrupt nature of fire and harvest with a similar response of spectral indices between the disturbance index and normalized burn ratio. We conclude that the spectral trajectory approach also captures multi-year stress that could be caused by climate, acid deposition, pathogens, partial harvest, thinning, etc. Our study focused on understanding the transferability of previously successful methods to new ecosystems and found that this automated method does not perform with the same accuracy in Pacific Northwest forests. Using a robust accuracy assessment, we demonstrate the difficulty of transferring change attribution methods to other ecosystems, which has implications for the development of automated detection/attribution approaches. Widespread disturbance was found within AVHRR-negative anomalies, but identifying causal factors in LTSS with adequate mapping accuracy for fire and insects proved to be elusive. Our results provide a background framework for future studies to improve methods for the accuracy assessment of automated LTSS classifications.

[1]  A. Solomon,et al.  Forest dynamics in Oregon landscapes: evaluation and application of an individual-based model. , 2007, Ecological applications : a publication of the Ecological Society of America.

[2]  Jeffrey G. Masek,et al.  Large Area Scene Selection Interface (LASSI): Methodology of Selecting Landsat Imagery for the Global Land Survey 2005. , 2009 .

[3]  S W Pacala,et al.  Contributions of land-use history to carbon accumulation in U.S. forests. , 2000, Science.

[4]  Chengquan Huang,et al.  Global characterization and monitoring of forest cover using Landsat data: opportunities and challenges , 2012, Int. J. Digit. Earth.

[5]  Russell G. Congalton,et al.  A review of assessing the accuracy of classifications of remotely sensed data , 1991 .

[6]  Joanne C. White,et al.  A new data fusion model for high spatial- and temporal-resolution mapping of forest disturbance based on Landsat and MODIS , 2009 .

[7]  N. Coops,et al.  Estimating the reduction in gross primary production due to mountain pine beetle infestation using satellite observations , 2010 .

[8]  E. Davidson,et al.  Classifying successional forests using Landsat spectral properties and ecological characteristics in eastern Amazônia , 2003 .

[9]  David P. Turner,et al.  A Carbon Budget for Forests of the Conterminous United States , 1995 .

[10]  Michael A. Wulder,et al.  Surveying mountain pine beetle damage of forests: A review of remote sensing opportunities , 2006 .

[11]  R. Houghton Revised estimates of the annual net flux of carbon to the atmosphere from changes in land use and land management 1850 – 2000 , 2003 .

[12]  Nicholas C. Coops,et al.  Estimating forest productivity in the eastern Siskiyou Mountains of southwestern Oregon using a satellite driven process model, 3-PGS , 2001 .

[13]  S. Sader,et al.  Detection of forest harvest type using multiple dates of Landsat TM imagery , 2002 .

[14]  Jan Verbesselt,et al.  Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology , 2013, Remote. Sens..

[15]  C. Tucker,et al.  NASA’s Global Orthorectified Landsat Data Set , 2004 .

[16]  Sean P. Healey,et al.  Application of two regression-based methods to estimate the effects of partial harvest on forest structure using Landsat data , 2006 .

[17]  C. Woodcock,et al.  Making better use of accuracy data in land change studies: Estimating accuracy and area and quantifying uncertainty using stratified estimation , 2013 .

[18]  D. Schimel,et al.  CONTINENTAL SCALE VARIABILITY IN ECOSYSTEM PROCESSES: MODELS, DATA, AND THE ROLE OF DISTURBANCE , 1997 .

[19]  Zhiqiang Yang,et al.  Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr — Temporal segmentation algorithms , 2010 .

[20]  Alan R. Gillespie,et al.  Structural stage in Pacific Northwest forests estimated using simple mixing models of multispectral images , 2002 .

[21]  Jay D. Miller,et al.  Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR) , 2007 .

[22]  Arjan J. H. Meddens,et al.  Carbon stocks of trees killed by bark beetles and wildfire in the western United States , 2013 .

[23]  Chengquan Huang,et al.  Remote Sensing ofEnvironment Snow-covered Landsat time series stacks improve automated disturbance mapping accuracy in forested landscapes , 2022 .

[24]  Wei Gao,et al.  Modeling the land surface heat exchange process with the aid of moderate resolution imaging spectroradiomer images , 2009 .

[25]  L. Dale,et al.  FORESTS , 1994, Restoration & Management Notes.

[26]  Joseph W. Sherlock,et al.  Calibration and validation of the relative differenced Normalized Burn Ratio (RdNBR) to three measures of fire severity in the Sierra Nevada and Klamath Mountains, California, USA , 2009 .

[27]  Robert E. Wolfe,et al.  Automated registration and orthorectification package for Landsat and Landsat-like data processing , 2009 .

[28]  R. Pontius,et al.  Death to Kappa: birth of quantity disagreement and allocation disagreement for accuracy assessment , 2011 .

[29]  P. Ciais,et al.  Consistent Land- and Atmosphere-Based U.S. Carbon Sink Estimates , 2001, Science.

[30]  Christopher S. R. Neigh,et al.  An Automated Approach to Map the History of Forest Disturbance from Insect Mortality and Harvest with Landsat Time-Series Data , 2014, Remote. Sens..

[31]  K. Shadan,et al.  Available online: , 2012 .

[32]  Warren B. Cohen,et al.  Trajectory-based change detection for automated characterization of forest disturbance dynamics , 2007 .

[33]  James D. Wickham,et al.  Estimating accuracy of land-cover composition from two-stage cluster sampling , 2009 .

[34]  W. Cohen,et al.  North American forest disturbance mapped from a decadal Landsat record , 2008 .

[35]  Christopher Justice,et al.  The impact of misregistration on change detection , 1992, IEEE Trans. Geosci. Remote. Sens..

[36]  Warren B. Cohen,et al.  A Method to Efficiently Apply a Biogeochemical Model to a Landscape , 2006, Landscape Ecology.

[37]  S. Goward,et al.  An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks , 2010 .

[38]  Nicholas C. Coops,et al.  Curve fitting of time-series Landsat imagery for characterizing a mountain pine beetle infestation , 2010 .

[39]  W. Cohen,et al.  Estimating the age and structure of forests in a multi-ownership landscape of western Oregon, U.S.A. , 1995 .

[40]  B. Quayle,et al.  A Project for Monitoring Trends in Burn Severity , 2007 .

[41]  T. Wilbanks,et al.  The first state of the carbon cycle report (SOCCR): The North American carbon budget and implications for the global carbon cycle. , 2007 .

[42]  R. Congalton,et al.  Accuracy assessment: a user's perspective , 1986 .

[43]  K. Hibbard,et al.  Postfire carbon pools and fluxes in semiarid ponderosa pine in Central Oregon , 2007 .

[44]  T. Hofacker,et al.  Insect and disease conditions in the United States, 1979-83 : , 1985 .

[45]  Steven E. Franklin,et al.  Pixel sampling of remotely sensed digital imagery , 1991 .

[46]  Dirk Pflugmacher,et al.  The Relative Impact of Harvest and Fire upon Landscape-Level Dynamics of Older Forests: Lessons from the Northwest Forest Plan , 2008, Ecosystems.

[47]  D. Verbyla,et al.  Evaluation of remotely sensed indices for assessing burn severity in interior Alaska using Landsat TM and ETM , 2005 .

[48]  Peter E. Thornton,et al.  Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests , 2002 .

[49]  W. Cohen,et al.  Detection of deforestation and land conversion in Rondônia, Brazil using change detection techniques , 2004 .

[50]  Jesse A. Logan,et al.  Changing temperatures influence suitability for modeled mountain pine beetle (Dendroctonus ponderosae) outbreaks in the western United States , 2006 .

[51]  Russell G. Congalton,et al.  Assessing the accuracy of remotely sensed data : principles and practices , 1998 .

[52]  C. Potter,et al.  Satellite-derived estimates of potential carbon sequestration through afforestation of agricultural lands in the United States , 2007 .

[53]  Giles M. Foody,et al.  Good practices for estimating area and assessing accuracy of land change , 2014 .

[54]  N. Coops,et al.  Estimation of insect infestation dynamics using a temporal sequence of Landsat data , 2008 .

[55]  Dongmei Chen,et al.  The effect of spatial autocorrelation and class proportion on the accuracy measures from different sampling designs , 2009 .

[56]  S. Gower Patterns and Mechanisms of the Forest Carbon Cycle1 , 2003 .

[57]  Steven W. Running,et al.  Testing a MODIS Global Disturbance Index across North America , 2009 .

[58]  W. Kurz,et al.  Mountain pine beetle and forest carbon feedback to climate change , 2008, Nature.

[59]  J. Grace,et al.  Delayed conifer mortality after fuel reduction treatments: interactive effects of fuel, fire intensity, and bark beetles. , 2009, Ecological applications : a publication of the Ecological Society of America.

[60]  C. Tucker,et al.  New vegetation index data set available to monitor global change , 2004 .

[61]  W. Cohen,et al.  Characterizing 23 Years (1972–95) of Stand Replacement Disturbance in Western Oregon Forests with Landsat Imagery , 2002, Ecosystems.

[62]  John R. G. Townshend,et al.  North American vegetation dynamics observed with multi-resolution satellite data , 2008 .

[63]  J. Hicke,et al.  Evaluating methods to detect bark beetle-caused tree mortality using single-date and multi-date Landsat imagery , 2013 .

[64]  Philip A. Townsend,et al.  Estimating the effect of gypsy moth defoliation using MODIS , 2008 .

[65]  S. VAN TUY,et al.  Disturbance and climate effects on carbon stocks and fluxes across Western Oregon USA , 2004 .

[66]  Steven E. Franklin,et al.  Forest Structure Classification in the North Columbia Mountains Using the Landsat TM Tasseled Cap Wetness Component , 2001 .

[67]  C. H. R I S T O P H E R P O T T E R,et al.  Major Disturbance Events in Terrestrial Ecosystems Detected Using Global Satellite Data Sets , 2003 .

[68]  Qihao Weng,et al.  A survey of image classification methods and techniques for improving classification performance , 2007 .

[69]  R. B. Jackson,et al.  CO 2 emissions from forest loss , 2009 .

[70]  S. Stehman,et al.  Statistical sampling to characterize recent United States land-cover change , 2003 .

[71]  Houghton,et al.  The U.S. Carbon budget: contributions from land-Use change , 1999, Science.

[72]  Stephen V. Stehman,et al.  Sampling designs for accuracy assessment of land cover , 2009 .

[73]  Molly E. Brown,et al.  EMD CORRECTION OF ORBITAL DRIFT ARTIFACTS IN SATELLITE DATA STREAM , 2010 .

[74]  Valerie A. Thomas,et al.  On-the-Fly Massively Multitemporal Change Detection Using Statistical Quality Control Charts and Landsat Data , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[75]  Arjan J. H. Meddens,et al.  Spatial and temporal patterns of Landsat-based detection of tree mortality caused by a mountain pine beetle outbreak in Colorado, USA , 2014 .

[76]  S. Goward,et al.  Characterization of the Landsat-7 ETM Automated Cloud-Cover Assessment (ACCA) Algorithm , 2006 .

[77]  H. Maffei,et al.  Forest health decline in a central Oregon mixed-conifer forest revisited after wildfire: a 25-year case study , 2007 .

[78]  R. B. Jackson,et al.  A Large and Persistent Carbon Sink in the World’s Forests , 2011, Science.

[79]  Jeffrey G. Masek,et al.  Carbon consequences of forest disturbance and recovery across the conterminous United States , 2012 .

[80]  G. Hurtt,et al.  NASA's Earth Science Use of Commercially Available Remote Sensing Datasets , 2008 .

[81]  Jakub Langhammer,et al.  Comparison of two types of forest disturbance using multitemporal Landsat TM/ETM+ imagery and field vegetation data. , 2009 .

[82]  E. Crist,et al.  Application of the Tasseled Cap concept to simulated thematic mapper data , 1984 .

[83]  C. Tucker,et al.  Northern hemisphere photosynthetic trends 1982–99 , 2003 .

[84]  Warren B. Cohen,et al.  A Landsat time series approach to characterize bark beetle and defoliator impacts on tree mortality and surface fuels in conifer forests , 2011 .

[85]  Edwin W. Pak,et al.  An extended AVHRR 8‐km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data , 2005 .

[86]  Martin Hais,et al.  Surface temperature change of spruce forest as a result of bark beetle attack: remote sensing and GIS approach , 2008, European Journal of Forest Research.