Improving estimates of forest disturbance by combining observations from Landsat time series with U.S. Forest Service Forest Inventory and Analysis data

To the best of our knowledge, one or more authors of this paper were federal employees when contributing to this work. This is the publisher’s final pdf. The published article is copyrighted by Elsevier and can be found at: http://www.journals.elsevier.com/remote-sensing-of-environment/.With earth's surface temperature and human population both on the rise a new emphasis has been placed on\ud monitoring changes to forested ecosystems the world over. In the United States the U.S. Forest Service Forest\ud Inventory and Analysis (FIA) program monitors the forested land base with field data collected over a permanent\ud network of sample plots. Although these plots are visited repeatedly through time there are large temporal gaps\ud (e.g. 5–10 years) between remeasurements such that many forest canopy disturbances go undetected. In this\ud paper we demonstrate how Landsat time series (LTS) can help improve FIA's capacity to estimate disturbance\ud by 1.) incorporating a new, downward looking response variable which is more sensitive to picking up change\ud and 2.) providing historical disturbance maps which can reduce the variance of design-based estimates via\ud post-stratification. To develop the LTS response variable a trained analyst was used to manually interpret 449\ud forested FIA plots located in the Uinta Mountains of northern Utah, USA. This involved recording cause and timing\ud of disturbances based on evidence gathered from a 26-year annual stack of Landsat images and an 18-year,\ud periodically spaced set of high resolution (~1 m) aerial photographs (e.g. National Aerial Image Program, NAIP\ud and Google Earth). In general, the Landsat data captured major disturbances (e.g. harvests, fires) while the air\ud photos allowed more detailed estimates of the number of trees impacted by recent insect outbreaks. Comparing\ud the LTS and FIA field observations, we found that overall agreement was 73%, although when only disturbed\ud plots were considered agreement dropped to 40%. Using the non-parametric Mann–Whitney test, we compared\ud distributions of live and disturbed tree size (height and DBH) and found that when LTS and FIA both found non-stand\ud clearing disturbance the median disturbed tree size was significantly larger than undisturbed trees,\ud whereas no significant difference was found on plots where only FIA detected disturbance. This suggests\ud that LTS interpretation and FIA field crews both detect upper canopy disturbances while FIA crews alone add\ud disturbances occurring at or below canopy level. The analysis also showed that plots with only LTS disturbance\ud had a significantly greater median number of years since last FIA measurement (6 years) than plots with both\ud FIA and LTS disturbances (2.5 years), indicating that LTS improved detection on plots which had not been field\ud sampled for several years. Next, to gauge the impact of incorporating LTS disturbances into the FIA estimation\ud process we calculated design-based estimates of disturbance (for the period 1995–2011) using three response\ud populations 1.) LTS observations, 2.) FIA field observations, and 3.) Combination of FIA and LTS observations.\ud The results showed that combining the FIA and LTS observations led to the largest and most precise (i.e. smallest\ud percent standard error) estimates of disturbance. In fact, the estimate based on the combined observations\ud (486,458 ha, +/−47,101) was approximately 65% more than the estimate derived solely with FIA data\ud (294,295 ha, +/−44,242). Lastly, a Landsat forest disturbance map was developed and tested for its ability to\ud post-stratify the design-based estimates. Based on relative efficiency (RE), we found that stratification mostly improved\ud the estimates derived with the LTS response data. Aside from insects (RE = 1.26), the estimates of area affected\ud by individual agents saw minimal gain, whereas the LTS and combined FIA + LTS estimates of total disturbance saw modest improvement, with REs of 1.43 and 1.50 respectively. Overall, our results successfully demonstrate\ud two ways LTS can improve the completeness and precision of disturbance estimates derived from FIA inventory\ud data