A rule-based semi-automatic method to map burned areas: exploring the USGS historical Landsat archives to reconstruct recent fire history

This study presents a new semi-automatic method to map burned areas by using multi-temporal Land Remote Sensing Satellite Program (Landsat) Thematic Mapper (TM) and Enhanced TM Plus (ETM+) images. The method consists of a set of rules that are valid especially when the post-fire satellite image has been captured shortly after the fire event. The overall accuracy of the method when applied to two case studies in Mt Parnitha and Samos Island in Greece were 95.69% and 93.98%, respectively. The commission and omission errors for Mt Parnitha were 6.92% and 10.24%, while those for Samos Island were 3.97% and 8.80%, respectively. Between the two types of error, it is preferred to minimize omission errors, since commission errors can be easily identified as part of product quality assessment and algorithm tuning procedures. The rule-based approach minimizes human interventions and makes it possible to run the mapping algorithm for a series of images that would otherwise need extensive time investment. In case of failure to capture burned areas correctly, it is possible either to make some adjustments by modifying the thresholding coefficients of the rules, or to discard some of the rules, since some editing is usually required to correct errors following the automated extraction procedures. When this method was applied to a series of US Geological Survey (USGS) Landsat TM and ETM+ archived satellite images covering the periods 1984–1991 and 1999–2009, a total of 1773 fires were identified and mapped from six different scenes that covered Attica and the Peloponnese in Greece. The majority of uncaptured burned areas corresponded to fires with size classes of 0–1 ha and 1–5 ha, where the loss in capturing fire scars is generally significant. This was expected since it is possible that small fires, identified and recorded by forest authorities, may not have been captured by satellite data due to limitations arising either from the spatial resolution of the sensor or imposed by the temporal series, which do not systematically cover the full period.