Validation of the Two Standard MODIS Satellite Burned-Area Products and an Empirically-Derived Merged Product in South Africa

The 500-m Moderate Resolution Imaging Spectroradiometer (MODIS) burned area products, MCD45A1, MCD64A1, and a merged product were validated across six study sites in South Africa using independently-derived Landsat burned-area reference data during the fire season of 2007. The objectives of this study were to: (i) investigate the likelihood of the improved detection of small burns through an empirically-derived merged product; (ii) quantify the probability of detection by each product using sub-pixel burned area measures; and, (iii) compare the mean percent concurrence of burned pixels between the standard products over a ten-year time series in each site. Results show that MCD45A1 presented higher detection probabilities (i.e., 3.0%–37.9%) for small fractions ≤50%, whereas MCD64A1 appeared more reliable (i.e., 12.0%–89.2%) in detecting large fractions >50% of a burned MODIS pixel, respectively. Overall, the merged product demonstrated improved detection of the burned area in all fractions. This paper also demonstrates that, on average, >50% of MODIS burned pixels temporally concur between OPEN ACCESS Remote Sens. 2014, 6 1276 the MCD45A1 and MCD64A1 products in each site. These findings have significant implications for fire monitoring in southern Africa and contribute toward the understanding of the range and of the sources of errors present in the MODIS burned area products.


Introduction
Fire is a major cause of land surface changes with impacts on the climate system, vegetation composition and the chemical composition of the atmosphere [1][2][3].The frequency of the occurrence, distribution and severity of biomass burning are fundamental drivers of global climate, ecosystem dynamics and land cover change.Global satellite observations of biomass burning reveal that Africa has by far the highest occurrence of vegetation fires [4,5], and this is believed to constitute approximately 40% of annual carbon emissions globally [6,7].Most of the global satellite observations are converted to fire detection proxies by use of specific algorithms.Two distinct satellite-based detection methods have been developed, namely, active fire [8] and burned area mapping [3,9].Only a fraction of the actively burning fires during satellite overpass is detected [10], and thus, this methodology is limited in its ability to estimate the actual extent of burning.On the other hand, burned area mapping provides a more direct estimate of the impacts of fires, required to address various environmental questions and management decisions.
In the last few years, several global burned-area thematic products derived from low-to moderate-resolution satellite data based on various algorithms have become publicly available [1,5,[11][12][13][14].Due to the differences in the data and algorithms, the derived products exhibit inherent systematic biases, which would therefore require their validation.Numerous validation studies of the burned area products have been reported in the literature [2,[15][16][17][18][19][20][21].Overall, the validation studies focused on three major inter-related aspects: (a) quantitative analyses of the relationship between the spectral reflectance of the burned areas and product mapping accuracy; (b) quantification of burned area detection by the product; and (c) analysis of the influence of the sub-pixel burned area on the detection accuracy of the product.This paper addresses the last two aspects.A key feature of accuracy assessment was in the use of a limited selection of reference data and time periods that would not sufficiently characterize the errors accompanying a particular product [13], and this called for the need for rigorous regional validation, particularly those defined using dissimilar algorithms [19].In this regard, product accuracy information could be more useful if product discrepancies are described in association to biome type and fire season [22].This would aid in the selection and interpretation of the most suitable product for a particular region.
Studies validating the Moderate Resolution Imaging Spectroradiometer (MODIS) burned area product (MCD45A1) and the MODIS direct broadcast (DB) burned area product (MCD64A1) in South Africa were limited to independently-derived reference datasets located in the savanna and fynbos biomes [9,19,21].The present study expands the validation of both standard MODIS products (i.e., MCD45A1 and MCD64A1) across six sites, encompassing not only savanna and fynbos, but also grassland and exotic pine forest plantations in South Africa during the fire season of 2007; by applying the traditional thematic accuracy reporting method, i.e., error matrix [9,19] and MODIS sub-pixel burned area analysis.Furthermore, this study introduces and validates a 500-m merged MODIS burned area product, which was recently derived by temporally aggregating the existing MCD45A1 and MCD64A1 products.The objectives of this study were to: (i) investigate the likelihood of improved detection (i.e., lower omission error) of a small burned area through an empirically-derived merged product; (ii) quantify the probability of detection by each of the three products using sub-pixel burned area measures; and (iii) compare the mean percent concurrence of burned pixels between the standard products over a ten-year time series in each site.These findings have significant implications for fire monitoring in southern Africa and contribute toward the understanding of a range and sources of errors present in the MODIS burned area products.

Description of Study Area
The January-April and June-September fire seasons of 2007 were considered in the analysis for various study sites (Figure 1), while the active fire data [23] and expert knowledge of the geographic regions susceptible to biomass burning were applied in the selection of six validation sites (Figure 1) in the savanna, grassland, fynbos biomes [24], as well as the exotic pine forest plantations [25,26].Mean annual rainfall [27,28] in the study areas ranges from approximately 260-975 mm per annum (Table 1).
These sites encompass the main fire-prone biomes that typically burn almost every year, partly as a result of their dominant grassy layer, climate and topography [29,30].

Ancillary Datasets Used
In the current study, the national vegetation map [24] and the national land cover (NLC 2000) product [25,26] were used as ancillary datasets to provide the spatial detail of the region's biome types.Furthermore, the 30-m Landsat 5 Thematic Mapper (TM) images (Table 1) were provided by the South African National Space Agency (SANSA) in tiff format processed to Level 1T (terrain and precision corrected).The temporal differences of the before and after fire images had different time intervals ranging from 16 to 64 days.At least 16 ground control points (GCPs) were selected uniformly throughout both TM images of each scene (with a total root mean square (RMS) error threshold of <0.5 pixel) and used to re-register the TM pairs precisely.Consequently, burned area reference datasets were independently derived from multi-date TM pairs (Figure 2) following the validation protocol described in [31] to validate the MODIS burned area products.

MODIS Burned Area Products
The MCD45A1 product is generated from 500-m MODIS cloud-free surface reflectance time series data using a predictive bidirectional reflectance modeling approach [3].On the other hand, the MCD64A1 product uses the -active-fire based burned area mapping algorithm‖ developed by [9] applied on the 500-m MODIS cloud-free surface reflectance imagery [32].It is important to note that active fire (training) data used in the latter algorithm may not fully represent regions dominated by small and fragmented burned areas, and therefore, such regions would generally result in large omission errors according to MCD64A1 [9].Unlike MCD64A1, the approach used to generate the MCD45A1 product makes no use of training samples, but rather, relies heavily on change in surface reflectance and a static wavelength threshold [3,5,13].Both MCD45A1 and MCD64A1 are released by the University of Maryland as monthly products that depict the location and approximate ordinal day (1-366) of burned areas at 500 m × 500 m spatial resolution [9,13].For the purpose of comparison, these products were re-projected to correspond with the Landsat Universal Transverse Mercator (UTM) projection of each scene while retaining their output pixel-dimension at 500-m spatial resolution [19].Unlike the MCD45A1 and MCD64A1 products, the 500-m empirically-derived MODIS merged burned-area product was produced in this study by temporally aggregating the monthly MCD45A1 and MCD64A1 products.Only the burn dates falling within each month were selected, and the dates outside of that particular month were discarded, as this data overlapped between consecutive months [3].Further, the individual pixel burn dates of both the MCD45A1 and MCD64A1 products were analyzed for temporal proximity (Table 2).Table 2 reveals that both products contained a significant number of pixels with burn dates that were too close together (i.e., eight or less days apart) to be two independent fire events.A temporal separation of zero means that both products had detected burned area for the same pixel on the same day.Nearly all the other burn dates with n > 0 occurred towards the end of one month and the beginning of the next month (Table 2).Therefore, if two burn dates for a particular pixel were within eight days of each other and spanned a month boundary, then they were treated as a single fire event, and the earlier of the two dates was chosen as definitive for the event.

Assessing the Error of MCD45A1 and MCD64A1 Products
A classification error assessment was conducted by using an error matrix [33] to compute the commission and omission errors (hereafter, Ce and Oe, respectively) present in the MODIS satellite-derived burned area products.In particular, Ce and Oe were derived for each study site by comparing the 500-m MODIS burned area products with the 30-m Landsat reference data using Equation ( 1), (1) where represents the intersection of the MODIS and Landsat burned area at 30-m resolution, are the pixels at 30-m resolution of the resampled MODIS burned area and are the pixels at 30-m resolution of the Landsat burned area.
An example of classification accuracy assessment using Equation (1) based on the MCD45A1 product is illustrated in Figure 3.It is evident that according to the error matrix, when the burned area product does not map the same extent of area burned as the reference data (Landsat interpreted burned area), a Ce and/or Oe is introduced in the product being evaluated [34].As a result, this method does not take into account the fractional burned area and is therefore susceptible to inherent classification biases.

Classification Error Assessment of the MCD45A1, MCD64A1 and Merged Products in Relation to the Fractional Sub-Pixel Burned Area
The classification error assessment [19] (Figure 3) considered the entire 500-m MODIS pixel as either burned or unburned and assessed the 500-m pixel in relation to all the 30-m Landsat reference pixels it contained.As a result, the fractional sub-pixel burned area within a 500-m pixel was not considered, and an overestimation of error during the assessment of the moderate-resolution burned-area products could be inherent [34].In the polygon-based analysis procedure considered here, a fixed 500-m reference grid was generated from the MODIS imagery and overlaid onto the Landsat interpreted burned areas (Figure 4).These grids were then utilized to compute proportions (fractional sub-pixel) of the burned area within each 500-m MODIS grid cell according to the Landsat data.
The computed proportions were subsequently categorized into four classes of ≤25% (6.25 ha), 25%-50% (6.25-12.5 ha), 50%-75% (12.5-18.75ha) and ≥75% (18.75 ha) of a MODIS pixel burned according to the Landsat TM data.The burned area products were analyzed and compared across each class in order to determine their probability of detection, derived as 100%; Oe for the different sizes of burned area fractions across the sites.Further, in order to account for the inherent error in burned areas due to the difference in resolution (30-m vs. 500-m pixels), a classification accuracy assessment of the MCD45A1, MCD64A1 and merged product based on 500-m MODIS pixels that are ≥50% burned according to the Landsat reference data was carried out based on indices derived from Equations ( 2) and (3) as follows: (2) Commission Error Correctly classified pixels Landsat burned area MCD45A1 Omission Error (3) where ba denotes the 500-m MODIS burned pixel and ref denotes the 500-m MODIS pixel burned to a specified fraction according to the Landsat data.These two equations were also applied to 500-m MODIS pixels that are ≥75% burned according to the Landsat reference data.The analysis method presented in Equations ( 2) and ( 3) addresses the concerns raised [9] about the reliability of detecting burned areas that correspond to at least 50% (~13 ha) of a 500-m MODIS pixel.

Analysis of Mean Percent Concurrence of Burned Pixels between MCD45A1 and MCD64A1 Products
The concurrence between the two MODIS burned area products was analyzed per pixel through time (January 2002 to December 2011) in each of the sites.For each pixel, the burned area dates of the two products were compared as follows: If the pixel was detected as burned by both products within n ≤ 8 days of each other, it was treated as a concurrence (i.e., an agreement).The concurrence could be expressed as the binary, AND, of the two products, while the merged product could be expressed as the binary, OR, of the two products (Table 3).The merged product was used to determine the total number of burns per pixel.Table 3.An example of the agreement and the result of a merge between the MCD45A1 and MCD64A1 products for a particular pixel.

MCD45A1 MCD64A1 Concurrence (Binary AND)
Merge (Binary OR) Further, the percentage concurrence per pixel was therefore computed as the total count of concurrences in proportion to the total number of times a pixel burned.Hereafter, the mean percent concurrence of burned area pixels between the MCD45A1 and MCD64A1 products was derived in each of the sites.

Assessing the Error of MCD45A1 and MCD64A1 Products
Table 4 presents the accuracy assessment of the MCD45A1 and MCD64A1 products per 30-m pixel through the derived Oe and Ce, respectively.In addition, the percentage of individual burns in different size classes (i.e., <50 ha, 50-100 ha, >100-500 ha and >500 ha, according to the Landsat reference data) intersecting with 500-m MODIS pixels of each burned area product is also presented.This percentage was computed as the count of individual burns per class intersecting (fully or partially) with the MODIS pixels of each product, relative to the total count of individual burns in that particular class.

Classification Error Assessment of the MCD45A1, MCD64A1 and Merged Products in Relation to the Fractional Sub-Pixel Burned Area
Table 5 presents the Oe and Ce of the MCD45A1, MCD64A1 and merged burned area products for 500-m MODIS pixels that were ≥50% and ≥75% burned according to the Landsat reference data.On average, when considering only the MODIS pixels ≥50% (13 ha) burned across all sites (Table 5), the Ce reduces by 7.82% and 15.24% (in absolute terms) for the MCD64A1 and MCD45A1 products, respectively, in comparison to the results in Table 4, where all fractional burned pixels were considered.However, the Oe was less affected than Ce (Table 5), and with increasing fractional cover (i.e., MODIS pixels ≥75% burned) the Oe reduces, while Ce increases [34].Overall, the reductions in Oe and Ce suggest improvement in the detection of burned areas corresponding to ≥50% of a MODIS pixel [9,34,35].
The Oe of the merged product (Table 5) based on MODIS pixels ≥50% burned had an improved detection of burned area across all study sites with an average reduction in Oe of 34.17% in contrast with 48.67% (MCD45A1) and 42.49% (MCD64A1).This suggests that the merged product appears to capture more burns, but at the cost of higher Ce (18.91%) than the individual MCD45A1 (Ce = 11.35%) and MCD64A1 (Ce = 14.83%) averaged across all sites.Overall, increasing the fraction burned in the MODIS pixels to ≥50% and ≥75% burned causes a drop in Oe coupled with an increase in Ce.

Analysis of Mean Percent Concurrence of Burned Pixels between MCD45A1 and MCD64A1 Products
Figure 6 illustrates the mean percent concurrence of burned pixels between the MCD45A1 and MCD64A1 products across the individual study sites during the period 2002-2011.Overall, more than 50% of the burned pixels concur between the MCD45A1 and MCD64A1 products in all the sites.In decreasing order, the products revealed the mean percentage concurrences of 71.13% (Sabie), 70.65% (Thabazimbi) 63.99% (Southern Kruger NP) 63.17% (Middelburg) 62.89% (Western Cape) and 54.92% (Free Sate).Further investigation is needed on the non-overlapping burned pixels.

Discussion
In the forestry plantations of Sabie (Figure 2), approximately 85% of the burned areas are small (<50 ha), and only a few spatially extensive burned areas (120-1,000 ha; >1,000 ha) were present.The MCD45A1 product did not show reliable detection of the burned areas over the pine forest plantations (Oe of 80.50%; Table 4).The variable post-fire reflectance changes on account of the (partially burned and/or unburned) high tree canopy may have compromised the detection accuracy of the MCD45A1 product in the forest plantations.This finding agrees with results of the comparisons with 1-km activefire data across low and high forest cover in South America [3].Their study reported higher burned area detection accuracies by the MCD45A1 product over low forest cover, as opposed to the high forest cover results, where most of the area burned was detected by the active-fire data and not MCD45A1.Furthermore, part of the algorithm used to define the MCD45A1 product assumes that post-fire reflectance of MODIS band 6 (1.640 µm) will drop following a forest fire [3], which may not be the case across all forest regions [3].Conversely, the MCD64A1 product detected most of the burned forest plantations in the Sabie site with a lower Oe of 38.50% (see Table 4), and that could be due to the exploitation of active-fire data in its approach to calibrate the search for burned areas [9].The spatial complexity of the burned area shapes in Sabie (Figure 2) resulted in unburned fractions within the MODIS pixels labeled as burned by the MCD64A1 product and, consequently, could have inflated Ce in this product.Additionally, the Oe of MCD64A1 reduces further to 18.02% (Table 5) over MODIS pixels with fractional burns ≥50% and, thus, showed an adequate probability of burned area detection (Figure 5).Despite the difference in the burned area performance of MCD45A1 and MCD64A1 in Sabie, the overlapping MODIS burned pixels of the two products were high over a ten year period in the Sabie site (Figure 6).
In the Western Cape site (Figure 2), the majority of burned area sizes were in the range of 0.36-1,000 ha, as well as a few exceptional burns greater than 1,000 ha, which were characterized by simple homogeneous shapes.This site encompass the fynbos biome that includes a mixture of fine-leaved low shrubs coupled with a sparse grass layer and is marked by winter rainfall and a summer season of burning, generally between December and April [30].The derived Oe for both products suggests that a significant portion of the burned area in the fynbos biome remained undetected (Table 4).This could be because spectral changes pre-to post-fire are not very large in this vegetation.For example, Figure 7 depicts a quantitative measure of the brightness level for the reference burned pixels over four selected sites, where core burned pixels have high values of the brightness level (close to one).The Western Cape site, above all, suggests a rather variable contrast from burns showing uni-modal distribution of values centered on 0.7 (Figure 7).In particular, the validation results (Table 5) based on MODIS pixels with fractional burns ≥50% show that the MCD45A1 and MCD64A1 product detected a burned area with low Ce of 4.08% and 3.12%, respectively.This outcome indicated a minor difference of about 4% and 12% compared to the derived Ce [21], corresponding to the MCD45A1 and MCD64A1 products during the 2007 fire season.Similarly, our validation results (Table 5) indicated that MCD45A1 detected more (i.e., lower Oe) of the burned fynbos than MCD64A1, which is also depicted in the detection rates across the individual burned areas below 500 ha (Table 4) and the probability of detection for MODIS pixel burned area proportions of 25%, 50% and 75% (Figure 5).Furthermore, similar detection rates (Table 4) are observed for both products over burned area sizes greater than 500 ha (and since only seven out of 668 individual burns in the Western Cape site were greater than 1,000 ha, with 4,105 ha being the largest burn), and this implies that we could expect the Oe for the MCD64A1 product to improve over larger burned area sizes in the fynbos biome [21].The burned areas in the grassland biome appeared (Figure 2) to be widely spaced with large burns of sizes 120-1,000 ha and >1,000 ha, constituting most of the total burned surface.It is evident that both MODIS products detected all burned areas larger than 500 ha over the grassland biome with variable detection rates of 43%-89% for burned area sizes between 100 and 500 ha (Table 4).The MCD45A1 product detected the area burned with fewer omissions (Tables 4 and 5) over the grassland biome than the MCD64A1 product; though a significant fraction dominated by small burned area sizes below 100 ha (i.e., the detection rates in Table 4) remained undetected by both products.In addition, the MCD45A1 product demonstrated the highest detection probabilities across all MODIS pixel burned area proportions (Figure 5) and, thus, provides adequate burned area mapping for both small and large burns in the respective grassland sites.In the Middelburg and Free State sites, the MCD64A1 product had a lower Ce of 25.14% and 20.85% (Table 4), which decreases further to 3.53% and 14%, respectively, over MODIS pixels ≥50% burned (Table 5) compared to MCD45A1.Given the burned area spatial characteristics (Figure 2) and sufficient brightness level (e.g., more than 50% available core burned pixels (Figure 7)) in the grassland sites, there were 45.08% and 36.83%(Figure 6) non-concurring burned pixels between the products in the Free State and Middelburg sites, respectively.As a result, the two MODIS products perform very differently in the grassland sites, regardless of the presence of numerous large burned areas.
The MODIS burned area products indicated variable detection rates (Table 4) for burned areas below 100 ha in the Southern Kruger NP and Thabazimbi sites.Generally, both products clearly missed all the burned areas between 50 and 100 ha in Thabazimbi compared to Southern Kruger NP (i.e., 66.67% detection rate).This observation, especially in Southern Kruger, demonstrates that depending on the degree of brightness, even small burned areas could reasonably be detected at lower resolution [36] (e.g., 500-m MODIS pixel).Nonetheless, in the savanna biome, the detection rate of the MCD45A1 and MCD64A1 products gradually increases with a burn size between 100 and 500 ha; and further, reaches a 100% detection rate for all burned areas greater than 500 ha, with the exception of only one burned area, which was missed by the MCD64A1 product in Thabazimbi (Table 4).The Southern Kruger NP site had few (i.e., 25) individual burned areas (in the reference data), which revealed complex shapes [18] and, thus, contributed toward the Oe (Table 4) of the products.Furthermore, the high Oe (~89%) for both products in Thabazimbi is mainly the consequence of about 868 individual burned areas that were not detected.In the aforementioned site, when considering 50% of fractional burned areas (Table 5), the Oe of the MCD45A1 and MCD64A1 products is significantly decreased by 34% and 35%, respectively.In addition, the Ce decreased by 12% for both preceding products.A similar trend is observed in the Southern Kruger NP (Table 5).These findings therefore suggest that MCD64A1 provides superior detection of area burned, especially with increasing fractional cover ≥50% (Figure 5) than the MCD45A1 product in the savanna biome.
A quantitative analysis of the influence of MODIS sub-pixel burned area on the probability of detection by the MCD45A1, MCD64A1 and merged products was also performed.It is evident from Figure 5 that both burned area products display a positive linear trend between the proportion of a burned MODIS pixel and the probability of being detected by MODIS across all validation sites.There are also large differences in the size of fires that are detected between sites (Table 4).These results suggest that a burned area has to be at least 13 ha in size (≥50% of a 500-m MODIS pixel) for the products to detect it with a probability of approximately 31%, 36% and 46% (for MCD45A1, MCD64A1 and merged, respectively) averaged over all six validation sites.These findings corroborate previous conclusions that a burned area proportion of ~50% can be considered an appropriate threshold for medium resolution burned area detection [9,18,34,37,38].The adopted minimum mapping unit of four pixels at Landsat 30-m resolution corresponding to 240 m or less during the production of burned area reference maps, reduced the inflation of Oe to some extent (Table 4).The differences in Oe and Ce (Tables 4 and 5) show that the sensitivity of the MODIS products becomes much higher over burned areas corresponding to ≥75% of a MODIS burned pixel.

Conclusions
This paper presented validation results for three MODIS satellite burned area products (including the merged product) across the four fire-prone biomes (encompassing savanna, grassland, fynbos and exotic pine forest plantations) in South Africa.The accuracy assessment methods were based on the error matrix and MODIS sub-pixel analysis, to compare the 500-m burned area products with the independently-derived reference data from multi-temporal Landsat TM imagery distributed over six validation sites.In particular, the classification procedures of the MCD45A1 and MCD64A1 products have important differences, which might contribute to differences in their performances across the sites.These differences could be caused by the following: (i) the MCD45A1 product approach does not exploit any active fire information, whereas MCD64A1 makes use of active fire data to enhance the search for probable burned pixels in its approach; and (ii) the MCD45A1 product primarily relies on a bi-directional reflectance change detection approach, whereas the MCD64A1 product approach hinges on a burn-sensitive vegetation index (VI) as a primary indicator to discover the burned areas [3,5,9,39].
Overall, this study has demonstrated that merging the two approaches [3,9] into a single product holds the potential for the improved detection of burned area, particularly for small burns.Efforts are currently underway to include a merged product into the official MODIS products in collection 6 [40], and the method for merging the two presented in this paper could be adopted.Small errors in the MODIS burned area product collection 5 data used here will also be fixed in collection 6.Studies that employed locally-adapted MODIS burned area algorithms instead of the global algorithms tested here achieved a very high accuracy level [41][42][43].This suggest that there is merit in designing algorithms that may allow local experts to set specific parameters in the algorithm to better account for local vegetation properties and the spatial and spectral properties of local burned area.Efforts are also underway to develop automated burned area algorithms that operate on high-resolution imagery, such as Satellite Pour l'Observation de la Terre (SPOT)-4 and Landsat-7 Enhanced Thematic Mapper plus (ETM+), have been successfully demonstrated on a regional scale [44][45][46][47].Developing operational automated burned area products based on Sentinel-2 and Landsat-8 could alleviate issues related to the detection of small and spatially fragmented burns with coarser resolution satellite data.

Figure 1 .
Figure 1.The biome map showing the location of the validation sites and multi-temporal Landsat 5 Thematic Mapper (TM) acquisitions distributed across the four fire-prone biomes in South Africa.

Figure 2 .
Figure 2. Spatial analysis of the burned areas interpreted and produced from multi-temporal Landsat TM data over the six validation sites.The burned area patterns show a diversity of burned shapes, from narrow and/or fragmented (0.36, 12.5, 50 ha) to large and compact (>1,000 ha).

Figure 3 .
Figure 3.Comparison of Landsat (interpreted) burned-area reference data and the MCD45A1 product to derive the Ce and Oe by use of a confusion matrix.

Figure 4 .
Figure 4.The 500-m MODIS reference grids used to compute the fractional sub-pixel burned areas corresponding to the Landsat TM interpreted burns.

Figure 5 .
Figure 5. Histograms of MODIS sub-pixel burned area analysis used to depict the relationship between burned area proportions according to Landsat TM and the probability of detection by the three MODIS products over the six validation sites.

Figure 6 .
Figure 6.Mean percent concurrence of burned pixels between the MCD45A1 and MCD64A1 products based on a ten-year time series, January 2002 to December 2011, across the six study sites.

Figure 7 .
Figure 7. Histograms of the brightness levels for Landsat reference pixels interpreted as burned from the TM temporal-difference between vegetation index (VI) date2 and VI date1 across the selected sites.

Table 1 .
Landsat 5 TM scene pairs used to produce high-resolution burned-area reference datasets and the corresponding biome characteristics per site.

Table 2 .
Analysis of the temporal proximity of the burned pixel dates of the MCD45A1 and MCD64A1 products within monthly eight-day intervals for the period April 2000-April 2012, over southern Africa encompassing four Moderate Resolution Imaging Spectroradiometer (MODIS) tiles.

Table 4 .
Accuracy assessment of the official MODIS burned area product (MCD45A1) and the MODIS direct broadcast (DB) burned area product (MCD64A1) against the Landsat TM burned area reference data at 30-m resolution for individual validation sites.The percentage of individual burns in different size classes detected by the MODIS burned area product.

Table 5 .
Accuracy assessment of the official MODIS burned area product (MCD45A1), MODIS direct broadcast (DB) burned area product (MCD64A1) and merged MODIS burned area product based on the 500-m MODIS pixels that are ≥50% and ≥75% burned according to Landsat.