New GOES imager algorithms for cloud and active fire detection and fire radiative power assessment across North, South and Central America

Vegetation fires are a key global terrestrial disturbance factor and a major source of atmospheric trace gases and aerosols. Therefore, many earth-system science and operational monitoring applications require access to repetitive, frequent and well-characterized information on fire emissions source strengths. Geostationary imagers offer important temporal advantages when studying rapidly changing phenomena such as vegetation fires. Here we present a new algorithm for detecting and characterising active fires burning within the imager footprints of the Geostationary Operational Environmental Satellites (GOES), including consideration of cloud-cover and calculation of fire radiative power (FRP), a metric shown to be strongly related to fuel consumption and smoke emission rates. The approach is based on a set of algorithms now delivering near real time (NRT) operational FRP products from the Meteosat Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) imager (available from http://landsaf.meteo.pt/), and the GOES processing chain presented here is designed to deliver a compatible fire product to complete geostationary coverage of the Western hemisphere. Results from the two GOES imagers are intercompared, and are independently verified against the well regarded MODIS cloud mask and active fire products. We find that the detection of cloud and active fires from GOES matches that of MODIS very well for fire pixels having FRP > 30 MW, when the GOES omission error falls to less than 10%. The FRP of fire clusters detected near simultaneously by both GOES and MODIS have a bias of only 22 MW, and a similar bias is found when comparing near-simultaneous GOES East and GOES West FRP observations. However, many fire pixels having FRP < 30 MW remain undetected by GOES, probably unavoidably since it has a much coarser spatial resolution than MODIS. Adjustment using data from the less frequent but more accurate views obtained from high spatial resolution polar orbiting imagers could be used to bias correct regional FRP totals. Temporal integration of the GOES FRP record indicates that during the summer months, biomass burning combusts thousands of millions of tonnes of fuel daily across the Americas. Comparison of these results to those of the Global Fire Emissions Database (GFEDv2) indicate strong linear relationships (r² > 0.9), suggesting that the timely FRP data available from a GOES real-time data feed is likely to be a suitable fire emissions source strength term for inclusion in schemes aiming to forecast the concentrations of atmospheric constituents affected by biomass burning.     

Development of a virtual active fire product for Africa through a synthesis of geostationary and polar orbiting satellite data

We explore the ability to enhance landscape fire detection and characterization by constructing a ‘virtual’ fire product from a synthesis of geostationary and polar orbiting satellite data. Active fire pixels detected by the Spinning Enhanced Visible and Infrared Imager (SEVIRI) and the Moderate Resolution Imaging Spectroradiometer (MODIS) were spatially and temporally collated across Africa between February 2004 and January 2005. Coincident fire pixels detected by SEVIRI and MODIS were used to populate an empirical database of frequency density (f-D) distributions of fire radiative power (FRP). Frequency density distributions of FRP measured by SEVIRI at 5.0° grid cell resolution and 15-minute temporal resolution were then cross referenced in the database to a set of counterpart f-D distributions of FRP measured by MODIS. This procedure resulted in the first generation of a ‘virtual’ fire product that exhibits the full continental coverage and high temporal resolution of SEVIRI whilst quantifying fire pixel counts and FRP with accuracies approaching those of MODIS. Diurnal cycles extracted from the virtual fire product indicate that SEVIRI measures a greater proportion of the active fire pixels and FRP potentially detectable by MODIS during the day due to the increased prevalence and stronger radiant contribution of highly energetic fire pixels. On a daily basis (sample size n = 365) the peak magnitude in the diurnal cycle of the virtual FRP occurred within the same 15-minute timeslot as in the native SEVIRI fire product. Continental-scale ignition and extinction events, however, were detected on average 44 min earlier (standard deviation s.d. = 40 min) and 137 min later (s.d. = 92 min), respectively. It is anticipated that the methodology developed here can be used to cross-calibrate active fire products between a variety of different satellite platforms.     

Active fire detection and characterization with the advanced spaceborne thermal emission and reflection radiometer (ASTER)

We present an automated fire detection algorithm for the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor capable of mapping actively burning fires at 30-m spatial resolution. For daytime scenes, our approach uses near infrared and short-wave infrared reflectance imagery. For nighttime scenes a simple short wave infrared radiance threshold is applied. Based on a statistical analysis of 100 ASTER scenes, we established omission and commission error rates for nine different regions. In most regions the probability of detection was between 0.8 and 0.9. Probabilities of false alarm varied between 9 × 10^− 8 (India) and 2 × 10^− 5 (USA/Canada). In most cases, the majority of false fire pixels were linked to clusters of true fire pixels, suggesting that most false fire pixels occur along ambiguous fire boundaries. We next consider fire characterization, and formulate an empirical method for estimating fire radiative power (FRP), a measure of fire intensity, using three ASTER thermal infrared channels. We performed a preliminary evaluation of our retrieval approach using four prescribed fires which were active at the time of the Terra overpass for which limited ground-truth data were collected. Retrieved FRP was accurate to within 20%, with the exception of one fire partially obscured by heavy soot.     

Fire radiative energy for quantitative study of biomass burning: derivation from the BIRD experimental satellite and comparison to MODIS fire products

A major focus in global change research is to quantify the amount of gaseous and particulate pollutants emitted from terrestrial vegetation fires. Determination of the emitted radiant energy released during biomass combustion episodes (the so-called fire radiative energy or FRE) has been suggested as a new tool for determining variations in biomass combustion rates and the rate of production of atmospheric pollutants. We review the physical principals behind the remote determination of FRE and present an alternative method for its derivation via analysis of ‘fire pixel’ radiances in the middle infrared spectral region. We compare our method to the existing FRE retrieval approach used in the EOS Moderate Resolution Imaging Spectro-radiometer (MODIS) fire products, and to retrievals of FRE based on derived fire temperature and area made via the so-called Bi-spectral method. We test each FRE retrieval method using both simulated data and imagery from a new experimental space mission, the Bi-spectral InfraRed Detection (BIRD) small satellite, which has sensors specifically designed for the study of active fires. We analyse near simultaneous MODIS and BIRD data of the fires that burned around Sydney, Australia in January 2002. Despite the markedly different pixel size and spectral coverage of these sensors, where the spatial extent of the fire pixel groups detected by MODIS and BIRD are similar, the derived values of FRE for these fires agree to within ±15 %. However, in certain fires, the lower spatial resolution of MODIS appears to prevent many of the less intensely radiating fire pixels being detected as such, meaning MODIS underestimates FRE for these fires by up to 46% in comparison to BIRD. Though the FRE release of each of these low intensity fire pixels is relatively low, their comparatively large number makes their overall FRE significant. Thus, total FRE release of the Sydney fires on 5 January 2002 is estimated to be 6.5×10⁹ J s⁻¹ via BIRD but 4.0×10⁹ J s⁻¹ via MODIS. The ability of BIRD to resolve individual fire fronts further allows the first accurate calculation of ‘radiative’ fireline intensity from spaceborne measurements, providing values of 15-75 kJ s⁻¹ m⁻¹ for fire fronts that are up to 9 km in length. Finally, we analyse the effectiveness of the satellite-based FRE retrieval methods in estimating the FRE from the active flaming and smouldering components only (FRE_Active, believed to be proportional to the rate of biomass combustion), despite the sensor receiving additional radiance from the ‘cooling ground’. The MIR radiance method appears particularly strong in this regard, allowing FRE_Active to be estimated to within ±30% in the range 100-100,000 J s⁻¹ m⁻². These results provide further confidence in the ability of spaceborne missions to derive physically meaningful values of FRE that could be used to support biomass burning emissions inventories. Future comparisons between FRE derived via MODIS and those from higher spatial resolution BIRD or airborne imagery may allow the MODIS-derived FRE values to be ‘calibrated’ for any systematic underestimation. We therefore expect FRE to become an important tool for enhancing global studies of terrestrial vegetation fires with infrared remote sensing, particularly as the majority of large fires are now imaged four times per day via the MODIS instruments on the Terra and Aqua spacecraft.     

Monitoring forest fires over the Indian region using Defense Meteorological Satellite Program-Operational Linescan System nighttime satellite data

Monitoring and management of forest fires is very important in countries like India where 55% of the total forest cover is prone to fires annually. The present study aims at effective monitoring of forest fires over the Indian region using Defense Meteorological Satellite Program-Operational Linescan System (DMSP-OLS) nighttime satellite data and to evaluate the active fire detection capabilities of the sensor. Nightly DMSP-OLS fire products were generated from February to May 2005 (peak fire season) and analyzed to study the occurrence and behavior of fires over different forest physiognomies in Indian region. Fire products generated from DMSP-OLS were validated with ground observations of fire records from state forest departments to evaluate the accuracy of fire products. Further, inter-comparison of the DMSP-OLS derived fire products with contemporary fire products from Moderate resolution Imaging Spectroradiometer (MODIS) (both daytime and nighttime products) in addition to fires and burnt areas derived from Indian Remote sensing Satellite (IRS-P6) Advanced Wide Field Sensor (AWiFS) data has been done to analyze spatial agreement of fire locations given by the above sensors.Results from the DMSP-OLS fire products (derived from February to May 2005) over Indian region showed high forest fires in southern dry deciduous forests during February-March; central Indian dry and mixed deciduous forests during March-April; northeastern tropical forests during February-April and northern pine forests during May. Spatial pattern in fires showed a typical seasonal shift in fire activity from the southern dry deciduous forests to the northern pine forests and temperate forests as the fire season progressed. Statistical evaluation of DMSP-OLS fire products with ground observations showed an over all accuracy of 98%. Comparison of DMSP-OLS derived fires with consecutive MODIS and AWiFS derived fires for individual days indicated that 69% of the fires continued from current day (DMSP-OLS pass around ∼ 7 pm to ∼ 10 pm local time) to the next day (MODIS and AWiFS pass ∼ 10:30 am local time). Comparison of DMSP-OLS derived fires with burnt areas estimated from AWiFS showed that 98% of DMSP-OLS derived fires on the current day fell within the burnt area of AWiFS on subsequent day. Since the worst forest fires are those that extend from the current to the consecutive days, DMSP-OLS derived fires provide a valuable augmentation to the fires derived from other sensors operating in daytime.     

A MODIS direct broadcast algorithm for mapping wildfire burned area in the western United States

Improved wildland fire emission inventory methods are needed to support air quality forecasting and guide the development of air shed management strategies. Air quality forecasting requires dynamic fire emission estimates that are generated in a timely manner to support real-time operations. In the regulatory and planning realm, emission inventories are essential for quantitatively assessing the contribution of wildfire to air pollution. The development of wildland fire emission inventories depends on burned area as a critical input. This study presents a Moderate Resolution Imaging Spectroradiometer (MODIS) - direct broadcast (DB) burned area mapping algorithm designed to support air quality forecasting and emission inventory development. The algorithm combines active fire locations and single satellite scene burn scar detections to provide a rapid yet robust mapping of burned area. Using the U.S. Forest Service Fire Sciences Laboratory (FiSL) MODIS-DB receiving station in Missoula, Montana, the algorithm provided daily measurements of burned area for wildfire events in the western U.S. in 2006 and 2007. We evaluated the algorithm's fire detection rate and burned area mapping using fire perimeter data and burn scar information derived from high resolution satellite imagery. The FiSL MODIS-DB system detected 87% of all reference fires > 4 km², and 93% of all reference fires > 10 km². The burned area was highly correlated (R² = 0.93) with a high resolution imagery reference burn scar dataset, but exhibited a large over estimation of burned area (56%). The reference burn scar dataset was used to calibrate the algorithm response and quantify the uncertainty in the burned area measurement at the fire incident level. An objective, empirical error based approach was employed to quantify the uncertainty of our burned area measurement and provide a metric that is meaningful in context of remotely sensed burned area and emission inventories. The algorithm uncertainty is ± 36% for fires 50 km² in size, improving to ± 31% at a fire size of 100 km². Fires in this size range account for a substantial portion of burned area in the western U.S. (77% of burned area is due to fires > 50 km², and 66% results from fires > 100 km²). The dominance of these large wildfires in burned area, duration, and emissions makes these events a significant concern of air quality forecasters and regulators. With daily coverage at 1-km² spatial resolution, and a quantified measurement uncertainty, the burned area mapping algorithm presented in this paper is well suited for the development of wildfire emission inventories. Furthermore, the algorithm's DB implementation enables time sensitive burned area mapping to support operational air quality forecasting.     

Addressing the spatiotemporal sampling design of MODIS to provide estimates of the fire radiative energy emitted from Africa

Satellite-based estimates of the fire radiative power (FRP) and energy (FRE) emitted from open biomass burning are affected by the spatiotemporal resolution of polar-orbiting and geostationary sensors. Here the impacts of the MODIS sampling design on estimates of FRE are characterized by superimposing the timing and extents of the Terra and Aqua granules onto the SEVIRI active fire product. Results for different land-cover types across Africa indicate that the FRE measured by SEVIRI during eight days is linearly related to the sum of FRP measured by SEVIRI within the MODIS granules. These relationships are least variable during the height of the fire season when diurnal cycles of FRP measured by SEVIRI are most consistent. Relationships between FRE and the sum of FRP developed using the SEVIRI active fire product are directly applied to the sum of FRP retrieved from the MODIS Terra and Aqua Climate Modeling Grid (CMG) fire products. Estimates of FRE from MODIS herein agree within 5% of those obtained from previously published methods, but remain a factor of 0.72 times those obtained by adjusting SEVIRI measurements of FRE to account for low spatial resolution detection biases. An examination of the MODIS scan geometry suggests that the latter underestimation is attributed to the coupling between a MODIS imaging artefact referred to as the “bow-tie” effect and the typical calculation used to retrieve the sum of FRP from the MODIS CMG fire products. Depending on the availability of MODIS scan angle information, we offer rigorous and simplified calculations to account for the bow-tie effect. Applying the simplified adjustment to the MODIS CMG fire products yields national estimates of monthly FRE that are 1.44 times greater than originally predicted.     

Mapping canopy damage from understory fires in Amazon forests using annual time series of Landsat and MODIS data

Understory fires in Amazon forests alter forest structure, species composition, and the likelihood of future disturbance. The annual extent of fire-damaged forest in Amazonia remains uncertain due to difficulties in separating burning from other types of forest damage in satellite data. We developed a new approach, the Burn Damage and Recovery (BDR) algorithm, to identify fire-related canopy damages using spatial and spectral information from multi-year time series of satellite data. The BDR approach identifies understory fires in intact and logged Amazon forests based on the reduction and recovery of live canopy cover in the years following fire damages and the size and shape of individual understory burn scars. The BDR algorithm was applied to time series of Landsat (1997-2004) and MODIS (2000-2005) data covering one Landsat scene (path/row 226/068) in southern Amazonia and the results were compared to field observations, image-derived burn scars, and independent data on selective logging and deforestation. Landsat resolution was essential for detection of burn scars < 50 ha, yet these small burns contributed only 12% of all burned forest detected during 1997-2002. MODIS data were suitable for mapping medium (50-500 ha) and large (> 500 ha) burn scars that accounted for the majority of all fire-damaged forests in this study. Therefore, moderate resolution satellite data may be suitable to provide estimates of the extent of fire-damaged Amazon forest at a regional scale. In the study region, Landsat-based understory fire damages in 1999 (1508 km²) were an order of magnitude higher than during the 1997-1998 El Niño event (124 km² and 39 km², respectively), suggesting a different link between climate and understory fires than previously reported for other Amazon regions. The results in this study illustrate the potential to address critical questions concerning climate and fire risk in Amazon forests by applying the BDR algorithm over larger areas and longer image time series.     

Field work and statistical analyses for enhanced interpretation of satellite fire data

Because their broad spatial and temporal coverage, satellites provide the main source of fire data for Amazonia. A key to the application of these tools for environmental studies is the appropriate interpretation of the data they provide. To enhance the interpretation of satellite fire data for this region, we collected ground-based data on fires in 2001 and 2002 using a simple and passive method, and statistically related these data to corresponding estimates from AVHRR and MODIS fire products using error matrices. Multiple methods of analyses from simple to complex produced qualitatively similar results. Total accuracies for both fire products were very high (> 99%) and dominated by accurate (> 99%) non-fire detection. Kappa statistics and fire-detection accuracies were substantially lower, with omission errors higher than commission errors. Results calculated using several different sets of spatial-matching parameters of analysis showed that Kappa was 1-10.6% for AVHRR, and 0-1.4% for MODIS. User's accuracy for fires was 0-40% for AVHRR and 3-100% for MODIS. Producer's accuracy for fires was 0-8% for AVHRR and 0-1% for MODIS. Statistical evaluations of potential explanatory factors showed that fire size and sampling time were dominant factors for low accuracies. Results from this study indicate that current satellite fire products are providing a limited sample of the fire activity in the region, and that ground-based analyses can substantially contribute to the interpretation of these products.     

Estimation of instantaneous net surface longwave radiation from MODIS cloud-free data

This paper develops a statistical regression method to estimate the instantaneous Downwelling Surface Longwave Radiation (DSLR) for cloud-free skies using only the satellite-based radiances measured at the Top Of the Atmosphere (TOA), and subsequently combines the DSLR with the MODIS land surface temperature/emissivity products (MOD11_L2) to estimate the instantaneous Net Surface Longwave Radiation (NSLR). The proposed method relates the DSLR directly to the TOA radiances in the MODIS Thermal InfraRed (TIR) channels provided that the terrain altitude and the satellite Viewing Zenith Angle (VZA) are known. The simulation analysis shows that the instantaneous DSLR could be estimated by the proposed method with the Root Mean Square Error (RMSE) of 12.4 W/m² for VZA = 0 and terrain altitude z = 0 km. Similar results are obtained for the other VZAs and altitudes. Considering the MODIS instrumental errors of 0.25 K for the TOA brightness temperatures in channels 28, 33 and 34, and of 0.05 K for channels 29 and 31, and of 0.35 K for channel 36, the overall retrieval accuracy in terms of the RMSE is decreased to 13.1 W/m² for the instantaneous DSLR. Moreover, a comparison of MODIS derived DSLR and NSLR are done with the field measurements made at six sites of the Surface Radiation Budget Network (SURFRAD) in the United States for days with cloud-free conditions at the moment of MODIS overpass in 2006. The results show that the bias, RMSE and the square of the correlation coefficient (R²) between the MODIS derived DSLR with the proposed method and the field measured DSLR are 20.3 W/m², 30.1 W/m² and 0.91 respectively, and bias = 11.7 W/m², RMSE = 26.1 W/m² and R² = 0.94 for NSLR. In addition, the scheme proposed by Bisht et al. [Bisht, G., Venturini, V., Islam, S., & Jiang, L. (2005). Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear-sky days. Remote Sensing of Environment, 97, 52-67], which requires the MODIS atmospheric profile product (MOD07) and also the MODIS land surface temperature/emissivity products (MOD11_L2) as inputs, is used to estimate the instantaneous DSLR and NSLR for comparison with the field measurements as well as the MODIS derived DSLR and NSLR using our proposed method. The results of the comparisons show that, at least for our cases, our proposed method for estimating DSLR from the MODIS radiances at the TOA and the resultant NSLR gives results comparable to those estimated with Bisht et al.'s scheme [Bisht, G., Venturini, V., Islam, S., & Jiang, L. (2005). Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear-sky days. Remote Sensing of Environment, 97, 52-67].     

Multi-resolution spectral analysis of wildfire potassium emission signatures using laboratory, airborne and spaceborne remote sensing

Thermal remote sensing studies of actively burning wildfires are usually based on the detection of Planckian energy emissions in the MIR (3-5 μm), LWIR (8-14 μm) and/or SWIR (1.0-2.5 μm) spectral regions. However, vegetation also contains a series of trace elements which present unique narrowband spectral emission lines in the visible and near infrared wavelength range when the biomass is heated to high temperatures during the process of flaming combustion. These spectral lines can be discriminated by detector systems that are less costly than the longer wavelength, actively cooled instruments more typically used in EO-based active fire studies. The main trace element resulting in the appearance of spectral emission lines appears to be potassium (K), with features at 766.5 nm and 769.9 nm. Here we study K-emission line spectral signature in laboratory scale fires using a field spectrometer, at a series of moderately-sized woodland and shrubland fires using airborne imagery from a new compact hyperspectral imager (HYPER-SIM.GA) operating at a relatively fine spectral sampling interval (1.2 nm), and at large open wildfires using the EO-1 satellite's Hyperion sensor. We derive a metric based on band differencing of the spectral signal both close to and outside of the K-line region in order to quantify the magnitude of the K-emission signature, and find that variations in this metric appear to track quite well with the commonly used measures of fire radiometric temperature and fire radiative power (FRP). We find that substantial flaming activity is required to generate a potassium emission signature, but that once present this can be detected using airborne remote sensing even through a substantial smoke layer that apparently obscures fire across the remainder of the VIS spectral range. Being specific to flaming combustion, detection of the K-emission line signature could prove useful in refining estimates of the gases released in open wildfires, since trace gas emission factors can vary substantially between flaming and smouldering stages. Finally, we demonstrate the first identification of the K-emission line signature from space using the EO-1 Hyperion instrument, but find it detectable only in certain instances. We conclude that a finer spectral and spatial resolution than that offered by Hyperion is required for improved detection performance. Nevertheless, our results point to the potential effectiveness of airborne and spaceborne K-emission signature detection as a complement to the more common thermal remote sensing approaches to wildfire detection and analysis. Sensors targeting this application should consider careful placement of the measurement wavelengths around the location of the K-line wavelengths, in part to minimise influences from the nearby oxygen A-band features.     

Using multiple endmember spectral mixture analysis to retrieve subpixel fire properties from MODIS

The Moderate-Resolution Imaging Spectroradiometer (MODIS) sensors on NASA's Terra and Aqua satellites image most of the Earth multiple times each day, providing useful data on fires that cannot be practically acquired using other means. Unfortunately, current fire products from MODIS and other sensors leave large uncertainties in measurements of fire sizes and temperatures, which strongly influence how fires spread, the amount and chemistry of their gas and aerosol emissions, and their impacts on ecosystems. In this study, we use multiple endmember spectral mixture analysis (MESMA) to retrieve subpixel fire sizes and temperatures from MODIS, possibly overcoming some limitations of existing methods for characterizing fire intensities such as estimating the fire radiative power (FRP). MESMA is evaluated using data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) to assess the performance of FRP and MESMA retrievals of fire properties from a simultaneously acquired MODIS image, for a complex of fires in Ukraine from August 21, 2002. The MESMA retrievals of fire size described in this paper show a slightly stronger correlation than FRP does to fire pixel counts from the coincident ASTER image. Prior to this work, few studies, if any, had used MESMA for retrieving fire properties from a broad-band sensor like MODIS, or compared MESMA to higher-resolution fire data or other measures of fire properties like FRP. In the future, MESMA retrievals could be useful for fire spread modeling and forecasting, reducing hazards that fires pose to property and health, and enhancing scientific understanding of fires and their effects on ecosystems and atmospheric composition.     

Evaluation of Collections 4 and 5 of the MODIS Gross Primary Productivity product and algorithm improvement at a tropical savanna site in northern Australia

In this study, we assessed the accuracy of the MODIS (Moderate Resolution Imaging Spectroradiometer) GPP (gross primary productivity) Collections 4.5, 4.8 and 5 along with Leaf Area Index (LAI), fraction of absorbed Photosynthetically Active Radiation (fPAR), light use efficiency (LUE) and meteorological variables that are used to estimate GPP for a northern Australian savanna site. Results of this study indicated that the MODIS products captured the seasonal variation in GPP, LAI and fPAR well. Using the index of agreement (IOA), it was found that Collections 4.5 and 4.8 (IOA 0.89 respectively) agreed reasonably well with flux tower measurements between 2001 and 2006. It was also found that MODIS Collection 4.5 predicted the dry season GPP well (Relative Predictive Error (RPE) 4.17%, IOA 0.72 and Root Mean Square Error (RMSE) of 1.05 g C m^− 2 day^− 1), whilst Collection 4.8 performed better in capturing wet season dynamics (RPE 1.11%, IOA 0.80 and RMSE of 0.91 g C m^− 2 day^− 1). Although the wet season magnitude of GPP was predicted well by Collection 4.8, an examination of the inputs to the GPP algorithm revealed that MODIS fPAR was too high, but this was compensated by PAR and LUE that was too low. Although LAI and fPAR estimated by Collection 5 were more accurate, GPP for this Collection resulted in a much lower value (RPE 25%) due to errors in other factors. Recalculation of MODIS GPP using site specific input parameters indicated that MODIS fPAR was the main reason for the differences between MODIS and tower derived GPP followed by LUE and meteorological inputs. GPP calculated using all site specific values agreed very well with tower data on an annual basis (IOA 0.94, RPE 6.06% and RMSE 0.83 g C m^− 2 day^− 1) but the early initiation of the growing season calculated by the MODIS algorithm was improved when the vapor pressure deficit (VPD) function was replaced with a soil water deficit function. The results of this study however, reinforce previous findings in water limited regions, like Australia, and incorporation of soil moisture in a LUE model is needed to accurately estimate the productivity.     

An intercomparison of three remote sensing-based energy balance models using Large Aperture Scintillometer measurements over a wheat-corn production region

This paper compares three remote sensing-based models for estimating evapotranspiration (ET), namely the Surface Energy Balance System (SEBS), the Two-Source Energy Balance (TSEB) model, and the surface temperature-vegetation index Triangle (TVT). The models used as input MODIS/TERRA products and ground measurements collected during the wheat and corn growth period in a subhumid climate at a measurement station in Yucheng, China. MODIS land surface temperature (LST) and leaf area index (LAI) products, corrected using ground-truth observations, were used in the three models. The TSEB model output of sensible (H) and latent (LE) heat fluxes were in good agreement with Large Aperture Scintillometer (LAS)-measured H and LE derived by residual (RMSD < 45 W/m²). Reasonable agreement was also obtained with the SEBS model output yielding RMSD for H of ~ 40 W/m² and LE ~ 55 W/m². However, the TVT model output resulted in poor agreement with the LAS-estimated H and LE with RMSD-values > 110 W/m². Using the uncorrected MODIS LST and LAI products resulted in a deterioration of the agreement in H and LE with LAS-estimated values for both the TSEB and SEBS models, whereas TVT performance improved marginally. These results indicate that the TSEB model yielded the closest agreement with the LAS-estimated fluxes using either the corrected or uncorrected MODIS inputs (LST and LAI). The SEBS model also computed reasonable H and LE values but was significantly more sensitive to errors in MODIS LST and LAI inputs than the TSEB model. In the TVT model, output of H and LE was unacceptable in either scenario of MODIS input which was attributable to errors in selection of the dry edge. With the TVT method, accurate determination of the dry edge end member is critical in regional ET estimation, but for humid and subhumid regions this end member may often be quite difficult to identify or encompass within a satellite scene.     

Monitoring turbidity in Tampa Bay using MODIS/Aqua 250-m imagery

We developed an approach to map turbidity in estuaries using a time series (May 2003 to April 2006) of 250-m resolution images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA's Aqua satellite, using Tampa Bay as a case study. Cross-calibration of the MODIS 250-m data (originally designed for land use) with the well-calibrated MODIS 1-km ocean data showed that the pre-launch radiometric calibration of the 250-m bands was adequate. A simple single scattering atmospheric correction provided reliable retrievals of remote sensing reflectance at 645 nm (0.002 < R_rs(645) < 0.015 sr^− 1, median bias = − 7%, slope = 0.95, intercept = 0.00, r² = 0.97, n = 15). A more rigorous approach, using a multiple scattering atmospheric correction of the cross-calibrated at-sensor radiances, retrieved similar R_rs(645). R_rs(645) estimates, after stringent data quality control, showed a close correlation with in situ turbidity (turbidity = 1203.9 × R_rs(645)^1.087, 0.9 < turbidity < 8.0 NTU, r² = 0.73, n = 43). MODIS turbidity imagery derived using the developed approach showed that turbidity in Hillsborough Bay (HB) was consistently higher than that in other sub-regions except in August and September, when higher concentrations of colored dissolved organic matter seem to have caused underestimates of turbidity. In comparison, turbidity in Middle Tampa Bay (MTB) was generally lowest among the Bay throughout the year. Both Old Tampa Bay (OTB) and Low Tampa Bay (LTB) showed marked seasonal variations with higher turbidity in LTB during the dry season and in OTB during the wet season, respectively. This seasonality is linked to wind-driven bottom resuspension events in lower portion of the Bay and river inputs of sediments in the upper portion of the Bay. The Bay also experiences significant interannual variation in turbidity, which was attributed primarily to changes in wind forcing. Compared with the once-per-month, non-synoptic in situ surveys, synoptic and frequent sampling facilitated by satellite remote sensing provides improved assessments of turbidity patterns and thus a valuable tool for operational monitoring of water quality of estuarine and coastal waters such as in Tampa Bay.     

Satellite mapping of CO emission from forest fires in Northwest America using MOPITT measurements

We present a study on MOPITT (Measurements Of Pollution In The Troposphere) detection of CO emission from large forest fires in the year 2000 in the northwest United States. Fire data used are from the space-borne Advanced Very High Resolution Radiometer (AVHRR) at 1-km resolution. The study shows that MOPITT can reliably detect CO plumes from forest fires whenever there are >30 AVHRR hotspots in a 0.25° × 0.25° grid, which is comparable to the pixel area of MOPITT in the region. The spatial CO pattern during the fire events is found to be consistent with the location and density of AVHRR hotspots and wind direction. While the increase of CO abundance inside the study area is closely correlated to the AVHRR-derived hotspot number in general (R > 0.75), the non-linearity of fire emission with fuel consumption is also observed. MOPITT can also capture the temporal variation in CO emission from forest fires through 3-day composites so it may offer an opportunity to enhance our knowledge of temporal fire emission over large areas. The CO emission is quantitatively estimated with a one-box model. The result is compared with a bottom-up approach using surface data including burnt area, biomass density, and fire emission factors. If mean emission factors for the region are used, the bottom-up approach results in total emission estimates which are 10%-50% lower than the MOPITT-based estimate. In spite of the limitations and uncertainties addressed in this study, MOPITT data may provide a useful constraint on uncertain ground-based fire emission estimates.     

Detection and monitoring of African vegetation fires using MSG-SEVIRI imagery

An operational procedure is presented that allows detecting active fires based on information from Meteosat-8/SEVIRI over Africa. The procedure takes advantage of the temporal resolution of SEVIRI (one image every 15 min), and relies on information from SEVIRI channels (namely 0.6, 0.8, 3.9, 10.8 and 12.0 µm) together with information on illumination angles. The method is based on heritage from contextual algorithms designed for polar, sun-synchronous instruments, namely NOAA/AVHRR and MODIS/TERRA-AQUA. A potential fire pixel is compared with the neighboring ones and the decision is made based on relative thresholds as derived from the pixels in the neighborhood.An overview is provided of results obtained for January and July 2007, respectively over Northern Africa (NAfr) and Southern Africa (SAfr), paying special attention to the spatial and temporal distribution of active fires. In both NAfr and SAfr, two types of vegetation clearly dominate in terms of fire activity, namely tree-covered areas, containing 40% of total fires observed, and shrub-covered areas, with 25% (19%) of total fires in NAfr (SAfr). However, marked differences were also to be found between the two regions; more than two-thirds (70%) of fires in SAfr were observed in land cover classes dominated by trees but the proportion is much lower (40%) in the case of NAfr. The duration of active fires in both regions tends to follow two-parameter generalized Pareto distributions, with both the scale and the shape parameters presenting very similar values for NAfr and SAfr.An assessment of the robustness of the algorithm, consistency of results and added value of the product was made by studying the daily cycle of fire activity over two regions located in northern and southern hemisphere Africa and by means of systematic comparisons against fire incidence reported in previous works and against hot spots extracted from the global daily active fire product developed by the MODIS Fire Team. The observed fire incidence by land cover class compares well with the results reported in previous works and it is shown that there is an overall coherence between results obtained from SEVIRI and MODIS when adequate spatial and temporal scales are chosen when performing the comparison. Data from MODIS and SEVIRI may be viewed as complementary, the latter having the added value of providing a much finer temporal resolution that allows uncovering certain aspects of fire behavior, namely the characterization of daily fire cycles.     

Development and analysis of a 12-year daily 1-km forest fire dataset across North America from NOAA/AVHRR data

Fires in boreal and temperate forests play a significant role in the global carbon cycle. While forest fires in North America (NA) have been surveyed extensively by U.S. and Canadian forest services, most fire records are limited to seasonal statistics without information on temporal evolution and spatial expansion. Such dynamic information is crucial for modeling fire emissions. Using the daily Advanced Very High Resolution Radiometer (AVHRR) data archived from 1989 to 2000, an extensive and consistent fire product was developed across the entire NA forest regions on a daily basis at 1-km resolution. The product was generated following data calibration, geo-referencing, and the application of an active fire detection algorithm and a burned area mapping algorithm. The spatial-temporal variation of forest fire in NA is analyzed in terms of (1) annual and monthly patterns of fire occurrences in different eco-domains, (2) the influence of topographic factors (elevation zones, aspect classes, and slope classes), and (3) major forest types and eco-regions in NA. It was found that 1) among the 12 years analyzed, 1989 and 1995 were the most severe fire years in NA; 2) the majority of burning occurred during June-July and in low elevation zones (< 500 m) with gentle slopes (< 10°), except in the dry eco-domain where more fires occurred in higher elevation zones (> 2000 m); 3) most fires occurred in the polar eco-domain, sub-arctic eco-division, and in the taiga ( boreal forests), forest-tundras and open woodlands eco-provinces in the boreal forests of Canada. The tendency for multiple burns to occur increases with elevation and slope until about 2500 m elevation and 24° slope, and decreases therefore. In comparison with ground observations, the omission and commission errors are on the order of 20%.     

Testing the potential of multi-spectral remote sensing for retrospectively estimating fire severity in African Savannahs

The remote sensing of fire severity is a noted goal in studies of forest and grassland wildfires. Experiments were conducted to discover and evaluate potential relationships between the characteristics of African savannah fires and post-fire surface spectral reflectance in the visible to shortwave infrared spectral region. Nine instrumented experimental fires were conducted in semi-arid woodland savannah of Chobe National Park (Botswana), where fire temperature (T_max) and duration (dt) were recorded using thermocouples positioned at different heights and locations. These variables, along with measures of fireline intensity (FLI), integrated temperature with time (T_sum) and biomass (and carbon/nitrogen) volatilised were compared to post-fire surface spectral reflectance. Statistically significant relationships were observed between (i) the fireline intensity and total nitrogen volatilised (r² = 0.54, n = 36, p < 0.001), (ii) integrated temperature (T_sum−μ) and total biomass combusted (r² = 0.72, n = 32, p < 0.001), and (iii) fire duration as measured at the top-of-grass sward thermocouple (dt_T) and total biomass combusted (r² = 0.74, n = 34, p < 0.001) and total nitrogen volatilised (r² = 0.73, n = 34, p < 0.001). The post-fire surface spectral reflectance was found to be related to dt and T_sum via a quadratic relationship that varied with wavelength. The use of visible and shortwave infrared band ratios produced statistically significant linear relationships with fire duration as measured by the top thermocouple (dt_T) (r² = 0.76, n = 34, p < 0.001) and the mean of T_sum (r² = 0.82, n = 34, p < 0.001). The results identify fire duration as a versatile measure that relates directly to the fire severity, and also illustrate the potential of spectrally-based fire severity measures. However, the results also point to difficulties when applying such spectrally-based techniques to Earth Observation satellite imagery, due to the small-scale variability noted on the ground. Results also indicate the potential for surface spectral reflectance to increase following higher severity fires, due to the laying down of high albedo white mineral ash. Most current techniques for mapping burned area rely on the general assumption that surface albedo decreases following a fire, and so if the image spatial resolution was high enough such methods may fail. Determination of the effect of spatial resolution on a sensor's ability to detect white ash was investigated using a validated optical mixture modelling approach. The most appropriate mixing model to use (linear or non-linear) was assessed using laboratory experiments. A linear mixing model was shown most appropriate, with results suggesting that sensors having spatial resolutions significantly higher than those of Landsat ETM+ will be required if patches of white ash are to be used to provide EO-derived information on the spatial variation of fire severity.     

Estimating the area of stubble burning from the number of active fires detected by satellite

To estimate the area affected by stubble burning in southern Australia, use was made of observations from the MODIS (MODerate resolution Imaging Spectroradiometer) on the Terra and Aqua satellites. The burnt area (BA) was calculated from the number of active fires, known as fire hot spots (FHS) using parameters estimated from a survey of farms in the agricultural area of south-western Australia. The study also served as a ground validation of the capability and limitations of the MODIS sensor and associated algorithm for detection of small agricultural fires.During the period from 1 March to 23 May 2005, 3240 unique FHS associated with stubble burning were detected. The majority of these FHS occurred in the afternoon in the last 3 weeks of April. To estimate the total area of stubble burning associated with these FHS, a survey was sent to 2066 farmers. This survey determined for each farm, the number of fields in which stubble was burnt, average size of field burnt (A), crop types burnt, dates and time of day of burning. Responses were received from 273 farms, 38% of whom reported over 500 stubble burns. The 3240 FHS were intersected with the polygons of the farm boundaries to determine the proportion (P) of stubble burns detected using MODIS. Only 13% (± 3%) of the stubble burns recorded in the farm survey were detected. Average field size burnt was 75 ha (± 6 ha). Total BA was calculated as: BA = A ^? 3240/P, which gave an estimated area of 1.87 million ha. This area was 27% of the total cropped area in south-western Australia. This level of stubble burning was similar to that determined by a 1990s survey in the State of Victoria.Neither cloud cover nor field size was significantly correlated with the low number of stubble fires that were detected. Therefore it was concluded that many stubble burns went undetected because of the lack of coincidence between the time of the MODIS overpass and when stubble burns were initiated. Also the use of additional sensors such as the Advanced Very High Resolution Radiometer (AVHRR) on NOAA satellites with afternoon overpasses would improve the fraction of stubble burns detected.Across the whole of Australia where winter cropping occurs, there was a high correlation (r² = 0.96) between FHS and total cropped area in each State. This provided the basis for extrapolating the field results from south-western Australia, to estimate the total area of stubble burning in southern Australia for 2005.