Burned Area Mapping of an Escaped Fire into Tropical Dry Forest in Western Madagascar Using Multi-Season Landsat OLI Data

A human-induced fire cleared a large area of tropical dry forest near the Ankoatsifaka Research Station at Kirindy Mitea National Park in western Madagascar over several weeks in 2013. Fire is a major factor in the disturbance and loss of global tropical dry forests, yet remotely sensed mapping studies of fire-impacted tropical dry forests lag behind fire research of other forest types. Methods used to map burns in temperature forests may not perform as well in tropical dry forests where it can be difficult to distinguish between multiple-age burn scars and between bare soil and burns. In this study, the extent of forest lost to stand-replacing fire in Kirindy Mitea National Park was quantified using both spectral and textural information derived from multi-date satellite imagery. The total area of the burn was 18,034 ha. It is estimated that 6% (4761 ha) of the Park’s primary tropical dry forest burned over the period 23 June to 27 September. Half of the forest burned (2333 ha) was lost in the large conflagration adjacent to the Research Station. The best model for burn scar mapping in this highly-seasonal tropical forest and pastoral landscape included the differenced Normalized Burn Ratio (dNBR) and both uniand multi-temporal measures of greenness. Lessons for burn mapping of tropical dry forest are much the same as for tropical dry forest mapping—highly seasonal vegetation combined with a mixture of background spectral information make for a complicated analysis and may require multi-temporal imagery and site specific techniques.


Introduction
Globally, tropical forests account for the majority of the earth's remaining forest resources [1]. Thirty years ago, tropical dry forests (TDF) outnumbered tropical humid forests nearly two-to-one [2], but today, land occupied by tropical humid forests surpasses that of TDF [1]. TDFs support a higher human population density than do tropical humid forests [3] and they are at high risk of disturbance from a variety of human activities such as agricultural conversion, hunting and gathering, livestock grazing, fuelwood collection, and selective logging [1,2,[4][5][6]. Natural disturbances such as fire, drought, insect outbreaks, and cyclones also take a toll on TDF leading to detrimental positive feedback loops between natural and human-induced disturbances. TDF are also often overlooked in forest protection schemes [7].
Remote sensing is a cost effective means to map and monitor changes in forest landscapes, but research in mapping TDF lags behind other forest types [7] due to (1) inconsistency in defining what is TDF [7][8][9]; (2) phenological challenges of this seasonal forest type [8,10]; and (3) TDF's numerous vegetation types having similar spectral properties (e.g., spiny forest, deciduous forest, broadleaf forest) [7,10,11]. The challenges are compounded in disturbed TDF landscapes where there may be even more vegetation types. In particular, it can be difficult differentiating between recent and past burns scars, as well as between burn scars and agricultural lands. even more vegetation types. In particular, it can be difficult differentiating between recent and past burns scars, as well as between burn scars and agricultural lands.
Despite the "widespread and potential significant threat" of fire to global TDF [4], there is a paucity of remote sensing research concerning fire in TDF, and fire mapping research of TDF lags behind fire mapping of other forest types. Tropical forest fire studies, to date, have focused largely on escaped fires in tropical evergreen rain forests [12][13][14][15].
Madagascar's dry forests experienced one of the largest decreases of any forest type worldwide between 1980 and 2000 [16]. TDF in Madagascar are undergoing rapid loss and fragmentation [17] largely due to land clearing for agriculture and pasture. During the period 1990-2000, rates of deforestation in dry and spiny forest were nearly as high or higher (7.5% and 8.6% respectively) than that of Malagasy humid forest (7.8%) [18]. Smith [19] estimates that western dry forest extent decreased from 12.5% in 1950 to 2.8% in 1990. In the western region of Central Menabe, annual forest loss is estimated to be 1% between 1999 and 2005 [20]. Annual rates of deforestation in the Kirindy-Ambadira Forest Complex in Central Menabe (an area greater than 100,000 ha) vary by location but have steadily increased, from 0.42% in the period 1973-1992, to 0.8-1.68% between 2000 and 2008, and 2.55% from 2008 to 2010; the most extreme forest losses occurred as incursions from existing agricultural clearings [21].
Kirindy Mitea National Park (KMNP) was established to protect animals and plants of dry deciduous forests in the transition zone between the sub-arid and dry ecoregions of western Madagascar in Southern Menabe. The Menabe region of western Madagascar ( Figure 1) is comprised of TDF, shrubland, pasture, savanna, and agricultural lands. Forests in Menabe have a long history of human modification through agriculture, logging, and grazing [20,[22][23][24]. Fire is commonly used in Menabe to clear lands for agriculture and to elicit fresh growth of pasture grasses for cattle [17]; both savannas and abandoned agricultural clearings are maintained as zebu pasture [25,26]. A pasture-renewal fire escaped into the forest of KMNP and burned for several weeks in August-September 2013 coming to within 600 m of the Ankoatsifaka Research Station, home of Lewis' long-term study of the endangered Propithecus verreauxi (sifaka) [27]. While the fire did not enter the 1 km 2 research trail system, all intact forests in the Park support high levels of plant and animal diversity and any forest loss will have ensuing impacts on the biodiversity of the park [28,29]. A pasture-renewal fire escaped into the forest of KMNP and burned for several weeks in August-September 2013 coming to within 600 m of the Ankoatsifaka Research Station, home of Lewis' long-term study of the endangered Propithecus verreauxi (sifaka) [27]. While the fire did not enter the 1 km 2 research trail system, all intact forests in the Park support high levels of plant and Remote Sens. 2018, 10, 371 3 of 17 animal diversity and any forest loss will have ensuing impacts on the biodiversity of the park [28,29]. The purpose of this study was to map the location of a large stand-clearing fire and calculate the amount of forest lost in the major KMNP fire of 2013.
There are multiple methods of satellite-based wildfire mapping, but most rely on spectral indices calculated from remotely-sensed satellite images. The first commonly adopted spectral index designed especially for fire studies is the Normalized Burn Ratio (NBR) [30] which has been used for burn scar mapping [31][32][33] and fire severity mapping [34,35]. Classification of burn areas using NBR can be challenging when there are multi-aged burn scars in the study area [31,36]. A multi-temporal differenced image of NBR (dNBR), which represents fire occurrence within a time frame defined by two image dates can aid in identification of burned areas from specific fire events. Differenced NBR (dNBR and RdNBR) has been used exclusively in some studies to map burn areas [35,37,38].
Automatic image segmentation of burned areas using NBR/dNBR alone can be especially challenging when there are a variety of non-vegetated land classes such as fallow agriculture fields, bare soil, and multi-aged burned areas of varying degrees of reflectance (black to grey charcoal) [31,36], as was the case in this study site. To reduce the spectral confusion between multi-color soil patches, additional spectral information sensitive to other parameters such as wetness and greenness can be included [36]. Patterson and Yool [39] used components (brightness, greenness, and wetness) from a single date Kauth-Thomas Tasselled Cap Transformation (KT) to map burned areas in Arizona, and Henry [33] used the same to map fire scars in pine and scrub vegetation in Florida. Barret et al. [40] used KT in combination with other information to map burn severity in Alaska.
Healey [41] combined variants of the Tasseled Cap Transform into a suite of indices designed to detect stand-replacing disturbance. Their Disturbance Index (DI) is a combination of the three KT components arranged so as to exploit the observed differences in brightness value compared to greenness and wetness in pixels of cleared forest. Compared to simple vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the Second Modified Soil Adjusted Vegetation Index (MSAVI2), DI incorporates information from shortwave infrared spectrum (i.e., wetness) which is especially useful for assessing change in forest structure. Also, DI is resilient to variability in vegetation phenology between image dates [42].
The highly seasonal nature of this study site posed an additional challenge to mapping burn scars from the large fire event. Background phenological changes in vegetation and soil moisture are registered along with changes due to fire. Inclusion of multi-temporal spectral bands can aid in the differentiation between phenological changes and land cover changes by capturing a baseline change due only to seasonal phenology. In addition to uni-temporal spectral information described above, land cover change can be captured also by using multi-temporal data in the form of a linear transformation of spectral bands.
This study used both spectral and textural information from multi-date satellite imagery to determine the extent of forest loss during the 2013 fire. Burn mapping methods typically used for temperate forests were incapable of pinpointing burn scars from a single fire event. The best model for burn scar mapping in this highly-seasonal tropical forest and pastoral landscape included dNBR and both uni-and multi-temporal measures of greenness.

Study Area
Kirindy Mitea National Park (KMNP) is located in the southern Menabe region of western Madagascar along the Mozambique Channel, about 70 km south of Morondava ( Figure 1). Ankoatsifaka Research Station (20 • 47 11 S, 44 • 09 53 E) is located approximately 1.5 km south of the 1 ha forest study area which has trails cut at 25 m intervals. Forest distribution in both Central and Southern Menabe is closely tied to soils and geological substrate [43] with dense, dry deciduous forests located largely in the west nearer the coast, and savannas occupying large swaths in the east along the central highlands [43]. Forest likely has never occupied the entire bioclimatic region, as there is evidence here for early grasslands [44]. Large areas of forest have been converted to agriculture and pasture.
KMNP, located approximately 100 km south of the similarly-named, Kirindy Forest/Centre National de Formation, d'Etudes et de Recherches en Environnement et Foresterie (CNFEREF) (also known as the Swiss Forest), was established to protect animals and plants of dry deciduous forests in the transition zone between the sub-arid and dry ecoregions of western Madagascar in Southern Menabe. KMNP is characterized by high biotic diversity and high local endemism [45].
The park encompasses a variety of ecological features including mangroves, coral reefs, beaches, spiny forest, baobab groves, and dry deciduous forest. Mean annual temperature is 25 • C (range: 9-40 • C) [46], and the site is highly seasonal with a long dry season followed by a shorter rainy season. Average rainfall is 700 mm/year, with the majority falling between December and February [47].
The dominant forest vegetation in the park is tropical dry forest (Commiphora sp., Dalbergia sp., Tamarindus indica, Euphorbia sp., Grewia sp.) including bottle-trunked trees (pachypodiums, Delonix adansonoides, Euphoribia antso, and Givotia madagascariensisi) and baobabs (Adansonia fony and Adansonia grandidieri). There are also small bands of gallery forest along rivers, as well as patches of spiny thicket of the endemic family Didieraceae near the coast. The eastern edge of the park is characterized by wooded savanna that grades to grassland.
Forest loss within KMNP is especially prominent at the eastern edge of the park and along roads with incursions due to agriculture and grazing. Forest lost to agricultural activities is also concentrated along riparian areas. A fairly extensive road network within park boundaries, pre-dating park inauguration, also facilitates wood and animal extraction. In 2009, Cyclone Fanele caused substantial damage to the forests of KMNP resulting in the downing of large trees and the creation of forest gaps [46].
For each of the scenes, digital numbers (DN) were first converted to top-of-atmosphere reflectance using the reflectance rescaling factors provided in Landsat OLI metadata file [49]. Adjacent north and south scenes were then mosaicked in ERDAS Imagine. Finally, the mosaicked image was atmospherically corrected using dark object subtraction.

Burn Mapping
Fire maps may represent burn scars, fire perimeters, or fire severity. For this study, burned area was defined as fire that results in significant change in forest biomass and structure that can be detected by space-borne optical instruments. That is, burn scars resulting from a fire so severe that is a stand-clearing event [50]. Burned areas within remaining standing forest where understory was cleared and surviving trees were scorched up to heights of 1-2 m were not identified as burned area for this study, as these areas could not be reliably detected by 30 m remotely sensed imagery. Therefore, the total burn area estimated in this study is an underestimate of the actual area impacted by this fire event. Field data on burn area was first collected on the author's behalf in April 2014 by KMNP park staff present during the fire event, and then again more extensively by the author in June 2014 with the assistance of the same park staff members.
One difficulty with mapping burns in this study area using NBR alone was discriminating between areas that had burned just prior to June 2013 (and still appeared burned in the post-fire image), and those areas that burned after June 2013 during the large fire of interest. A differenced NBR (dNBR) was used to help distinguish between fires that occurred between 23 June and 27 September 2013 and those that occurred prior to 23 June 2013.
For calculation of multi-temporal indices, both a pre-fire image (23 June 2013) and post-fire image (27 September 2013) from the dry season were selected. The post-fire image was the first image date that pixels in the study area were no longer obscured by smoke given that the fire burned for weeks, and then smoldered even longer.
To help resolve confusion between charcoal and exposed soil, the KT [51] and the MSAVI2 (in lieu of NDVI because of its better performance in regions with high degree of exposed soil [52]) were calculated on the post-fire image [31,36].

Multi-Temporal Kauth-Thomas Transform
Common multi-temporal linear transformation techniques for fire mapping include Principal Components Analysis [53,54] and Multitemporal Kauth-Thomas (MKT) Transform [39,51,55,56]. MKT components represent inter-date change in brightness, greenness, and wetness components. Rogan and Yool [36] combined MKT brightness, greenness, and wetness change components in a supervised classification of fire in semi-arid lands in the American southwest. They achieved higher accuracies using a combination of components than by classifying single components alone. Potere et al. [57] used a combination of a single date KT with MKT components to identify forest clearings along the Appalachian Trail in eastern United States.
Coefficients for the KT Transform were initially developed for Landsat MSS [58] and later updated for Landsat TM by Crist and Cicone [59]. Output bands (brightness, greenness, and wetness) are well-described in Patterson & Yool [39]. To extend the single-date KT Transform to a multi-temporal linear transform approach, bands 2-7 of both dry season Landsat OLI images were layer-stacked into a 12 layer file. A 12 × 12 reflectance factor coefficient matrix was created using published Landsat OLI coefficients [60] and normalized as in Collins and Woodcock [51] ( Table 1). The MKT transformation matrix is an extension of Fung's MKT Gramm-Schmidt orthonormalization technique [55]. The matrix coefficients were used to transform the two images into both stable and change components of brightness, greenness, and wetness.

Burn Classification
Threshold partitioning (the use of hard threshold values to define categories) is generally used to identify burnt pixels or burn severity with NBR. The U.S. Forest Service adopted a thresholding methodology [30] for their Fire Effects Monitoring and Inventory system (FIREMON) [61]; however, there are limitations to this method. First, defined "universal" threshold values for burn severity are specific to geographic areas [31], and there has been little research using dNBR in tropical dry forests. Second, thresholds derived from statistical regression models require field severity measurements which may be difficult to collect, especially directly after the fire [62]. Additionally, threshold partitioning can result in misclassification of pixels when classifying burn severity in sites having both pre-fire low and high vegetated pixels [63].
An alternative to thresholding is classification. Given the relatively small extent of area sampled on the ground for this project, an unsupervised classification of burn metrics was deemed a superior method to threshold partitioning. In this way, field data could inform the identification of classes provided by the classification.
An initial classification using only dNBR confirmed that additional information would be necessary to reliably identify burned areas. The next step was to identify additional spectral bands or indices that might aid in the discrimination of burn patches when combined with dNBR. This list included: dMSAVI2 (June to September); differenced DI (June to September); post-fire Landsat OLI bands 2-7; Greenness, Brightness, and Wetness bands of Tasseled Cap Transform of post-fire image; change components of Greenness, Brightness, and Wetness of MKT using both dry season images (June and September). A series of unsupervised classifications were run using a combination of spectral indices in an additive process starting with dNBR. Three indices from the original list were selected for the final classification: dNBR (June to September), uni-temporal Tasseled Cap Greenness (September), and MKT change-in-Greenness (June and September). These were layer-stacked and then classified using k-means clustering with 30 classes and an Iterative Self-Organizing Data Analysis Technique (ISODATA) with 15 maximum iterations and a convergence threshold of 96%.
Multiple fire classes were output corresponding to decreasing fire severity; however, in this analysis, only those pixels that resulted in bare surfaces were classified as fire, as this was the most conservative estimate of fire on the landscape, and field data were not extensive enough to reliably classify under-canopy fire. A majority smoothing filter (3 × 3) was passed over the classification and the image was clipped to the land portion of KMNP to exclude the Park's marine islands. A total of 11 classes represented burned area where fire burned to bare soil. The image was reclassified to a binary burned/not burned image. Total area of burned and non-burned pixels was calculated in ArcGIS 10.1.
The accuracy of the burn classification was evaluated using 244 validation data points from three sources: locations geolocated on the ground in 2014; locations extracted from post-fire Google Earth Very High Resolution Imagery (VHRI) dated October and November 2013; and a random set of points placed in locations not covered by ground data or VHRI extracted from post-fire Landsat OLI raw spectral bands dated 27 September 2013. Google Earth VHRI are widely used for validation of map classifications [64][65][66][67].
Post-fire VHRI was not available across the entire study area, so those areas of KMNP with VHRI imagery dated October and November 2013 were identified in Google Earth. Polygons were then drawn in those regions of VHRI coverage in Google Earth Pro. These polygons were imported into the University of New Hampshire Cooperative Extension KML Tools Project (http://extension.unh.edu/ kmltools/) and using the Cruise Tool, a random set of points was created within the polygons.

Forest Classification
Dry season images provide the best discrimination between vegetation classes having different canopy structure [68,69], and difference images between wet season and dry season TDF imagery can greatly improve classification accuracy [70,71]. To control for vegetation differences due to seasonal changes, two pre-fire Landsat OLI images (Bands 2-7) were utilized in the pre-fire forest classification. The wet season pre-fire image was dated 20 April 2013 and the dry season pre-fire image was dated 23 June 2013. Pre-processing of the images included conversion of digital numbers (DN) to Remote Sens. 2018, 10, 371 7 of 17 top-of-atmosphere reflectance, mosaicking of the two Landsat scenes, and radiometric normalization by dark object subtraction.
Texture is a measure of heterogeneity in the tonal values of pixels, and texture measures have been used to distinguish between land cover types having different vegetation structure [73]. Textural measures including mean Euclidian distance and skewness (asymmetry of the distribution around its mean [74]) were also calculated using a 3 × 3 moving window. The two texture measures were calculated on bands 5-6 of Landsat OLI image dated 23 June 2013. Textural mathematical operators implemented in ERDAS Imagine followed those in Irons and Peterson [75].
The spectral indices and textural measures selected were layer-stacked with wet and dry season Landsat 8 OLI spectral bands 2-7 ( Table 2) and then classified using k-means clustering with 22 classes and an Iterative Self-Organizing Data Analysis Technique (ISODATA) with 15 maximum iterations and a convergence threshold of 96%. Of the 22 classes, six were identified as forest and 16 as non-forest. The image was reclassified to forest/non-forest raster, a majority smoothing filter (3 × 3) was passed over the classification, and the image was clipped to the land portion of KMNP. Total area of forest and non-forest pixels was calculated in ArcGIS 10.1.
The accuracy of the forest classification was evaluated using 150 validation data points from three sources: locations geolocated on the ground in 2014; locations extracted from post-fire Google Earth Very High Resolution Imagery (VHRI) dated June 2012; and a random set of points placed in locations not covered by ground data or VHRI extracted from post-fire Landsat OLI raw spectral bands dated 20 April 2013.
Areas of KMNP with VHRI imagery dated June 2012 were identified in Google Earth and polygons were then drawn in those regions. These polygons were imported into the University of New Hampshire Cooperative Extension KML Tools Project (http://extension.unh.edu/kmlTools/) and using the Cruise Tool, a random set of points was created within the polygons.

Burn Classification Map
The total area of stand-replacing burn within KMNP during the large August-September fire was 18,034 ha (Figure 2A). Because so much of the forest on the eastern edge had been previously cleared, just 4761 ha (26%) of the 18,034 ha burned area was identified as forest prior to the fire ( Figure 2B). In total, it is estimated that 6% of primary TDF burned in KMNP over the period 23 June to 27 September.
A confusion matrix was calculated and both producer's and user's accuracy of the burn classification map was well above 80% with an overall Kappa statistic of 0.9031 ( Table 2). The Kappa index of agreement has been a staple of remote sensing accuracy assessment for decades [76][77][78][79]; however, it is criticized as being flawed and difficult to interpret [80,81]. Pontius and Millones [80] Remote Sens. 2018, 10, 371 8 of 17 recommend using quantity disagreement and allocation disagreement, measures that can also be summarized from the cross-tabulation matrix. These measures take into the account both the relative proportions of categories in the accuracy assessment (quantity) and the spatial allocation of the categories (allocation).
The total disagreement of the burn classification was 2% representing equal amounts of quantity and allocation disagreement at 1% each. Low total disagreement suggests the classification properly allocated cover classes to the correct spatial location and in the correct proportions with respect to reference data. The vast majority of the KMNP pixels were unburned from this fire event (although, many were burned from previous events) (Figure 3). There was strong agreement between the map classification and test data. Burned area commission disagreement is much higher than unburned area suggesting there may be an overestimate of burned lands; however, the overall percent disagreement is low. The largest burn scar, located in close proximity to Ankoatsifaka Research Station, can definitely be attributed to the escaped fire and accounts for half of the forest lost in KMNP over the period (2333 ha). It is unknown the true extent of the burn attributable to the original escaped fire, although the large conflagration extending from center to the east is known to have originated from that source. Forest burns just to the north and south of the large burn area could be due to the same escaped fire (particularly the one extending from center to western edge); it is difficult to determine just how far that burn extends. Burns along the park boundary in the far north and far south cannot be attributed to the same source but they are included in this estimate as they did burn during the same period of study. Weather conditions during that season suggest that fires set in other locations may have also spread beyond target areas as well. The largest burn scar, located in close proximity to Ankoatsifaka Research Station, can definitely be attributed to the escaped fire and accounts for half of the forest lost in KMNP over the period (2333 ha). It is unknown the true extent of the burn attributable to the original escaped fire, although the large conflagration extending from center to the east is known to have originated from that source. Forest burns just to the north and south of the large burn area could be due to the same escaped fire (particularly the one extending from center to western edge); it is difficult to determine just how far that burn extends. Burns along the park boundary in the far north and far south cannot be attributed to the same source but they are included in this estimate as they did burn during the same period of study. Weather conditions during that season suggest that fires set in other locations may have also spread beyond target areas as well.
(particularly the one extending from center to western edge); it is difficult to determine just how far that burn extends. Burns along the park boundary in the far north and far south cannot be attributed to the same source but they are included in this estimate as they did burn during the same period of study. Weather conditions during that season suggest that fires set in other locations may have also spread beyond target areas as well.  It is important to note that the extent of fire during this period was greater than reported here, as some areas of forest were burned in the understory and in the canopy without stand removal. Without extensive ground-truthing of these understory burns, it is very difficult to accurately model the full extent of this type of fire. Furthermore, at this time, neither spectral nor thermal remote sensing methods can accurately detect changes in lower forest strata [15,82]. However, changes in the understory are arguably an especially important component of this forest's ecology.
Cyclones that frequent this region remove large trees and create forest gaps which then "release" suppressed subcanopy vegetation resulting in a dramatic increase in understory growth [83,84]. A forest that has lost understory seedlings and saplings to a fire may not be able to recover from a windstorm that removes overstory trees; rather these gaps will likely be invaded by grasses instead.
Even undisturbed forests are at risk from grass fires originating from pasture or agricultural lands which burn up to and just into forest edges. Forest edges recolonized by grass fuel further grass encroachment by increasing fire susceptibility, fuel loading, and fire intensity in a feedback process that can quickly transform large areas of tropical dry forest into savanna, shrubland, or grassland [14,85].
Windfall gaps colonized by grasses within TDF may represent an even higher fire risk to forests as these gap grasses are susceptible to fire ignition by wind-borne fire embers, and they may be subject to intentional clearing by fire for agricultural activities. Consequently, the interplay of fire and cyclone disturbance may be fueling grass and forest fuel regimes (typical of tropical forest edges) from within large blocks of forest.
While dNBR was instrumental in identifying burned areas, it underestimated the amount of the burn, particularly in areas where old clearings in June had a spectral signature similar to burned areas that occurred after June. The dNBR also falsely identified some unburned areas undergoing phenological change as having burned. A difference image of MSAVI2 of the two dry season images did not significantly differ from dNBR and provided no additional information.
A major challenge was separating fresh burn scars appearing on 23 June with those that appeared after. Given the frequency of fire on the landscape, a classification method was sought that would reliably differentiate between burns created at such close time intervals.
The addition of KT greenness to dNBR helped differentiate between fresh burn patches created just before 23 June and those attributed to the large fire after that date. Areas of forest that burned after 23 June were covered in ash leaving a unique gray spectral signature, while grassy areas that burned left a different signature. Bare patches that burned between April and June 2013 and then burned again in July/August 2013 left yet another spectral signature. The ultimate goal was to quantify the amount of forest burned, so it was not critical to identify re-burnt patches (as they would not be classified as forest in the pre-fire forest classification). Therefore, making this distinction did not alter the calculation of the extent of forest lost; however, if one were interested in modeling the trajectory of the fire or quantify the full extent of burn, regardless of land cover type, the addition of the Greenness band would facilitate such an analysis.
Whereas KT greenness helped resolve differences in the timing of burns in open clearings, MKT ∆G (change in greenness) helped resolve burns that left scars with a spectral signature not significantly different from background color which would have been otherwise overlooked by other indices. So while KT helped to avoid overestimation of burn scars, MKT aided in avoiding underestimating them.
The MKT ∆W (change in wetness) band failed to differentiate between burns just prior to 23 June and those just after, yet it did successfully differentiate between burns after 23 June and those that occurred at least a season prior to that date. The MKT ∆B (change in brightness) easily differentiated between dark burn scars just prior to 23 June and those afterwards, but it falsely labeled some old burn scars (before April 2013) located within recent burned areas as recent burns, and it failed to discriminate between white ash spots from 23 June from white ash spots after 23 June.
Together, dNBR and DI performed well in distinguishing forest burned since 23 June nearly obviating the need to perform a separate classification of pre-fire forest. Although the DI/dNBR combination showed strong promise for rapid classification, DI failed to distinguish between vegetation regrowth in old previously cleared patches (perhaps agriculture or early shrub growth) and tropical dry forest; thus, while DI did appropriately identify this as burn since 23 June, it did not differentiate between burned tropical dry forest and burned secondary vegetation. Furthermore, DI misidentified some clearings burned before 23 June as having been burned after that date. Conversely, both multi-temporal and uni-temporal measures of Greenness easily distinguished regions of regrowth in old cleared patches from tropical dry forest. In study areas without cleared patches in various stages of regrowth, DI may identify forest burned in a stand-replacing fire in a single step (without having to first classify pre-fire forest cover).

Forest Classification Map
The total land area of KMNP (excluding the small islands off the coast) is 128,666 ha, and the total area of TDF within this area, as of 23 June 2013, was estimated to be 80,339 ha.
A confusion matrix was calculated and both producer's and user's accuracy of the forest classification map was well above 80% with an overall Kappa statistic of 0.8541 (Table 3). The total disagreement of the forest classification was 6.8% which was comprised of quantity disagreement (1.6%) and allocation disagreement (5.2%). The difference in agreement measures suggest that while the classification performed very well in detecting the total number of forest pixels (quantity), the classification was not as adept at allocating them all to the correct locations on the map (allocation). Non-forest commission disagreement is slightly higher than forest suggesting there may be an overestimate of non-forested lands. Likewise, omission disagreement is slightly higher for forest suggesting an underestimate of forested lands (Figure 4).
Non-forest commission disagreement is slightly higher than forest suggesting there may be an overestimate of non-forested lands. Likewise, omission disagreement is slightly higher for forest suggesting an underestimate of forested lands (Figure 4). A major challenge to mapping forest burns in TDF landscapes is the initial classification of the TDF. Mapping of TDF requires considerable attention to both phenological changes and the presence of regenerating secondary vegetation. Classification of tropical dry forest is especially challenging because vegetation undergoes substantial spectral changes over the course of the year [7]. Dry season A major challenge to mapping forest burns in TDF landscapes is the initial classification of the TDF. Mapping of TDF requires considerable attention to both phenological changes and the presence of regenerating secondary vegetation. Classification of tropical dry forest is especially challenging because vegetation undergoes substantial spectral changes over the course of the year [7]. Dry season imagery is deemed more useful for TDF mapping because differences in forest structural properties are more distinctive during at this time. This is especially important when trying to differentiate between evergreen gallery forests and seasonally dry deciduous forest [68,86]. Inclusion of multi-season imagery can further improve classification by helping to distinguish between tropical dry forests having different degrees of deciduousness (e.g., xerophytic vs. seasonally dry deciduous forest, forest vs. pastureland) [87].
The commonly used vegetation index NDVI, has been shown to be significantly related to TDF canopy closure, plant diversity, plant richness and plant species composition [88,89]; however, NDVI saturates at high vegetation densities and researchers have thus developed "enhanced" and "soil-adjusted" options [90,91]. MSAVI2 is a modified version of the original soil-adjusted vegetation index which outperforms NDVI in regions with a large amount of exposed soil [52], and MSAVI2 performed substantially better than NDVI in this landscape. On the other hand, indices representing wetness such as NDWI and TCAP wetness did not improve the classification.
Haboudane et al. [72] modified existing indices to create the Modified Triangular Vegetation Index, an index that is more responsive to variations in green leaf area index (LAI), less resistant to the effects of chlorophyll absorption on green reflectance, and more resistant to the effects of soil properties on vegetation reflectance. In this landscape, MTVI2 was particularly helpful in helping to distinguish between secondary growth in old burn patches and TDF.
Texture measures were included in this classification to combat the difficulties of classifying land cover in the highly heterogeneous semi-arid environment; texture has been shown to improve upon spectral classification alone in these environments [92]. Textural measures have improved land cover classifications of various landscape types [93], as they are especially effective in distinguishing among "spectrally indistinct" cover classes [94]. The addition of texture measures to the forest classification was especially helpful in teasing apart those areas having similar spectral signal but different structures.
Only texture measures of near-infrared and shortwave infrared bands (Landsat OLI Bands 5 and 6) improved the classification. Texture measures have been shown to correlate with stem density [69], aboveground biomass [67], and canopy structure [70,71]. In this case, high mean Euclidian distance corresponded to high vegetation heterogeneity typical of disturbed lands; tropical dry forest was homogeneous compared to other vegetation types and registered very low measures of texture. In a study of textural measures (standard deviation), Nagendra et al. [89] also found a negative relationship between texture and plant diversity (e.g., undisturbed tropical forest) in a human-disturbed landscape.
While there are different types of TDF located within KMNP, the primary goal here was to map primary forest cover of any type, so it was therefore not important to distinguish between spiny thicket, gallery forest, baobab groves, and seasonal dry deciduous forest. Use of multi-season imagery and vegetation indices did improve this forest classification by allowing all forest types (merged into a single class) to be readily distinguished from pasture and other non-forest classes.

Conclusions
Despite the fact that the majority of land within the large burn scar adjacent to Ankoatsifaka Research Station (8441 ha) was not identified as forest, it is important to recognize that most of this burned area was recently forested, as evidenced by satellite imagery. Forest clearings materialize patch-by-patch over seasons and years. As such, it is easy to overlook the collective loss of forest.
This research is one of few burn mapping studies of tropical dry forests. Fire is a major factor in the disturbance and loss of global tropical dry forests, yet this avenue of research has received little attention. Remote sensing mapping studies of fire-impacted TDF lag behind studies of other forest types, including boreal forest. Furthermore, as of yet, there are no established standardized techniques for tropical dry forest burn mapping. Active fire map products derived from instruments such as NASA's MODIS (Moderate Resolution Imaging Spectroradiometer) and VIIRS (Visible Infrared Imaging Radiometer Suite) and ESA's ATSR and AATSR (no longer operational) identify active fires through detection of thermal anomalies at large spatial resolution (~550-1000 m). These instruments provide some insight into the number, and location of fires, but they do not provide any information on the extent of fire or forest lost.
The research presented here addresses a critical gap in understanding the challenges of analyzing fire through remotely sensed satellite imagery in TDF systems. Spectral information based primarily on change in greenness metrics resulted in highest accuracy suggesting that accounting for phenological changes in TDF (particularly senescence) significantly aids in burn mapping. Standards established to calculate burn severity in other areas do not necessarily hold true in tropical dry forest. With fire representing a major risk to the biome, there is an urgent need for additional studies on tropical dry forest fires. Land stewards, conservation planners, and wildlife managers all require accurate figures of forest lost to fire in order to manage resources sensibly.
While TDF may share common spectral characteristics across the globe, they are frequently embedded in a host of different-and largely anthropogenic-landscapes rendering them each fairly spectrally distinctive at the landscape scale. This further complicates efforts to establish standardized techniques for tropical dry forest burn mapping. Lessons for burn mapping of TDF are much the same as for TDF forest mapping-highly seasonal vegetation combined with a mixture of background spectral information complicate analysis and may require site specific techniques.