Seasonal and interannual variations in CO and ... - Wiley Online Library

2 downloads 173 Views 510KB Size Report
[5] In southern Africa, Swap et al. [2003] showed ..... [22] The NDVI and cloud mask data set were derived from the ... and desert masks determined by NDVI to remove pixels ..... Goode, and W. M. Hao (2003), Trace gas and particle emissions from .... Hoelzemann, J. J., M. G. Schultz, G. P. Brasseur, C. Granier, and M. Simon.
GLOBAL BIOGEOCHEMICAL CYCLES, VOL. 21, GB2011, doi:10.1029/2006GB002848, 2007

Seasonal and interannual variations in CO and BC emissions from open biomass burning in Southern Africa during 1998–2005 Akinori Ito,1 Akihiko Ito,1 and Hajime Akimoto1 Received 26 September 2006; revised 2 February 2007; accepted 16 February 2007; published 9 May 2007.

[1] We estimate the emissions of carbon monoxide (CO) and black carbon (BC) from

open vegetation fires in the Southern Hemisphere Africa from 1998 to 2005 using satellite information in conjunction with a biogeochemical model. Monthly burned areas at a 0.5degree resolution are estimated from the Visible InfraRed Scanner (VIRS) fire count product and the MODerate resolution Imaging Spectroradiometer (MODIS) burned area data set associated with the MODIS tree cover imagery in grasslands and woodlands. The monthly fuel load distributions are derived from a 0.5-degree terrestrial carbon cycle model in conjunction with satellite data. The monthly maps of combustion factors and emission factors are estimated using empirical models that predict the effects of fuel conditions on these factors in grasslands and woodlands. Our annually averaged effective CO and BC emissions per area burned are 27 g CO m2 and 0.17 g BC m2 which are consistent with the products of fuel consumption and emission factors typically measured in southern Africa. The CO and BC emissions from open vegetation burning in southern Africa range from 45 Tg CO yr1 and 0.26 Tg BC yr1 for 2002 to 75 Tg CO yr1 and 0.42 Tg BC yr1 for 1998. The monthly averaged burned areas from VIRS fire counts peak earlier than modeled CO emissions. This characteristic delay between burned areas and emissions is mainly explained by significant changes in combustion factors for woodlands in our model. Consequently, the peaks in CO and BC emissions from our bottom-up approach are identical to those from previous top-down estimates using the Measurement Of the Pollution In The Troposphere (MOPITT) and the Total Ozone Mapping Spectrometer (TOMS) Aerosol Index (AI) data. Citation: Ito, A., A. Ito, and H. Akimoto (2007), Seasonal and interannual variations in CO and BC emissions from open biomass burning in Southern Africa during 1998 – 2005, Global Biogeochem. Cycles, 21, GB2011, doi:10.1029/2006GB002848.

1. Introduction [2] Biomass burning releases many trace gases and aerosols, such as carbon monoxide, oxygenated volatile organic compounds and carbonaceous aerosols, into the atmosphere, which have a significant influence on ambient air quality and climate [Crutzen et al., 1979; Andreae and Crutzen, 1997]. Most of the biomass burning takes place in the tropical belt [Crutzen and Andreae, 1990]. However, fire statistics are rather insufficient in many areas of the tropical regions, even though the savanna and grassland fires are mainly associated with human practices. Thus many researchers have utilized satellite information to quantify the emissions from the biomass burning [Chuvieco and Kasischke, 2007]. [3] Global satellite information of burned areas [Simon et al., 2004; Tansey et al., 2004; Giglio et al., 2006a] have been applied to estimates of open biomass burning emissions [Hoelzemann et al., 2004; Ito and Penner, 2004; van 1 Frontier Research Center for Global Change, JAMSTEC, Yokohama, Japan.

Copyright 2007 by the American Geophysical Union. 0886-6236/07/2006GB002848

der Werf et al., 2006]. While a number of satellite-based burned area products are currently under development [e.g., Carmona-Moreno et al., 2005; Roy et al., 2005b; Plummer et al., 2006], multiyear data of burned area at a coarse resolution on a continental scale provided useful information for analyses on the seasonal and interannual variations to emphasize important scientific issues that remained important subjects for discussion [Barbosa et al., 1999]. The active fire observations have also been used for multiyear estimates of burned areas [van der Werf et al., 2003; Giglio et al., 2006a]. However, many studies have shown that the uncertainty of the burned area estimates is significantly large without a quantitative validation and calibration using high-resolution data [e.g., Kasischke et al., 2003; Boschetti et al., 2004]. [4] The MODerate resolution Imaging Spectroradiometer (MODIS) burned area data product is available at a 500 m resolution only for 2 months for the year of 2000 [Roy et al., 2002]. Using the MODIS burned area data, Ito and Penner [2005b] showed that the estimated ranges from regional modeling approaches for carbon monoxide (CO) emissions from open biomass burning were within the range of the estimates constrained by chemical transport models and measurements in 2000 [Arellano et al., 2004; Pe´tron et al.,

GB2011

1 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

2004]. Thus the MODIS burned area product may serve as a useful data to calibrate the active fire counts to the burned areas in the Southern Hemisphere Africa until long-term burned area data sets become available. [5] In southern Africa, Swap et al. [2003] showed the apparent inconsistency of the results between the commonly held understanding maximum burning seasons of August and September and the satellite measurements of burned areas, which indicated that the largest burned areas were observed during June and July [Silva et al., 2003]. The updated version of this burned area data set [Tansey et al., 2004] was used to estimate the seasonal variations of CO emissions, which also showed the peak in July [Ito and Penner, 2004; Sinha et al., 2004; Korontzi, 2005; Jain et al., 2006]. Here the peak in monthly data is defined as the calendar month during which the maximum value is estimated for a year. The 1- to 2-month offset was also reported between the peak in the MODIS fire counts [Justice et al., 2002] and the peaks in the Measurement Of the Pollution In The Troposphere (MOPITT) [Edwards et al., 2006] and the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) [Gloudemans et al., 2006] CO retrievals in southern Africa. [6] In addition to the fire affected areas, seasonal variations in the CO emissions are also influenced by the available biomass for combustion (i.e., fuel load) and burning conditions (i.e., flaming and smoldering fires). Litter fall occurs as the dry season progress, so that the amount of leaf litter increases seasonally [Hoffa et al., 1999]. van der Werf et al. [2003] used a biogeochemical model to represent the seasonality in the delivery of leaves to the litter pools. [7] The combustion factor (CF) and the emission factor (EF) are mainly characterized by composition, availability, and moisture content of biomass. To exploit the spectral signature of green vegetation for a measure of the fuel moisture content, the strong absorption of chlorophyll in the red and the high reflectance of vegetation in the nearinfrared portion provides a useful information [Tucker, 1979]. The Normalized Difference Vegetation Index (NDVI) is calculated as the difference in reflectance between the near-infrared and visible red wavelengths divided by their sum. Korontzi [2005] used the monthly NDVI data composites at a spatial resolution of approximately 64 km2 to represent the temporal and spatial variations of the regional fuel moisture condition, and related them to CF and EF. Verbesselt et al. [2006] used the measurements of time series greenness from VEGETATION instrument on board the Syste`m Pour l’Observation de la Terre (SPOT) at finer spatial resolution of about 1 km2 to examine them for a fire risk assessment. [8] The purpose of this paper is to improve the mismatch of 1 – 2 months, by incorporating the seasonal variations of burned areas with fuel loads, combustion factors, and emission factors at a fine resolution. For this purpose, satellite information and a biogeochemical model are used to estimate CO and black carbon (BC) emissions from open biomass burning in southern Africa during the 8 year period from 1998 to 2005. Section 2 describes data sets used for estimating the CO and BC emissions from open vegetation

GB2011

fires using satellite data and a biogeochemical model. Section 3 shows multiyear estimates of CO and BC emissions and comparisons with previous estimates. We also present an analysis of the sensitivity of the CO emissions to the fire products and effective CO emissions per area burned to explore the relative importance of temporal and spatial variations in different treatments of burned areas and seasonality in fuel conditions. Section 4 presents a summary of our findings.

2. Materials and Methods [9] In a modeling approach, the emissions of gases and particles per month ([Q]lit g species month1) within a land cover type (l) at location (i) for month (t) from open vegetation fires can be described by the following equation [Seiler and Crutzen, 1980]: ½QðXÞlit ¼ ½ BAlit  ½FLlit  ½CF lit  ½ EF lit ;

ð1Þ

where X is chemical species, BA is the burned area per month (m2 month1), FL is the fuel load (kg m2) expressed on a dry weight (DM) basis within each grid, CF is the combustion factor, and EF is the emission factor in grams of species per kilogram of dry matter burned. 2.1. Area Burned [10] Monthly burned areas at a 0.5-degree resolution are estimated from the Visible InfraRed Scanner (VIRS) fire count product [Giglio et al., 2003] and the MODIS burned area data set [Roy et al., 2002] associated with the MODIS tree cover (TC) imagery [Hansen et al., 2003] in grasslands and woodlands. Monthly 0.5-degree VIRS fire product is available for the Tropical Rainfall Measuring Mission (TRMM) latitudes from 38°N to 38°S, which covers the entire African continent [Giglio et al., 2003]. The MODIS burned area data set was obtained on the approximate day of burning at a 0.5 km resolution for 2 months in 2000 [Roy et al., 2002]. We used the pixels with burns occurring in July and September and on the last day of the previous months [Roy, 2003]. The MODIS tree and herbaceous cover data sets were collected at a resolution of 0.5 km  0.5 km from November 2000 to December 2001 [Hansen et al., 2002]. Comparing the validation database in Zambia with a subset of the global MODIS TC map in a validation exercise showed that the general tree cover distribution was well represented [Hansen et al., 2006]. [11] In southern Africa, two classifications of land cover types (i.e., grasslands and woodlands) have been used to estimate the emissions from open vegetation fires [Korontzi et al., 2003; Sinha et al., 2004; Ito and Penner, 2005b]. Following these studies, land cover types were classified into two classes, based on the MODIS percentage of tree cover. Grasslands were defined as those areas with less than or equal to 15% tree cover, and woodlands were defined as those regions with greater than 15% tree cover. [12] van der Werf et al. [2003] related burned area at a spatial resolution of 1° to the VIRS active fire counts and fractional tree cover using a linear regression model. Giglio et al. [2006a] calibrated the MODIS active fire observations

2 of 13

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

GB2011

Figure 1. Relationships between MODIS tree cover percentages and the ratios of VIRS fire counts to MODIS burned areas in (a) grasslands and (b) woodlands. The solid lines represent the regression lines. to the burned areas at 1° spatial resolution as a function of tree and herbaceous vegetation cover, and the mean size of monthly cumulative fire (i.e., pixel clusters). On the basis of these studies, we developed equations between the ratio of fire counts to burned areas and tree cover in grasslands and woodlands, separately. The 0.5-km resolution MODIS burned areas in grasslands and woodlands were aggregated to 0.5° resolution. The ratio between burned areas and fire counts was binned in each one percentage of tree cover, and then the average of the ratio was related to the percentage tree cover, which was the average value for the corresponding 0.5° grid cells. [13] For grasslands (Figure 1a), we find the following linearly decreasing function of the mean tree cover (%) for the ratio of fire count (FC) to burned area (km2), ½ BAgit =½FCgit ¼ 0:43  ½TC gi þ 57:74

 TC 15; r2 ¼ 0:79 : ð2Þ

The negative slope in the equation (2) is explained by one possible factor that increasing tree cover reduces the fire spread rate and increases the detection probability of hot spots and thus decreases the effective burned area per fire count [van der Werf et al., 2003; Giglio et al., 2006a]. [14] For woodlands (Figure 1b), we find the following two linear functions of tree cover, ½ BAwit =½FCwit ¼ 2:77  ½TC wi  40:89  15 < TC 28; r2 ¼ 0:96

ð3Þ

½ BAwit =½FCwit ¼ 0:80  ½TC wi þ 54:05  28 < TC 61; r2 ¼ 0:65 :

ð4Þ

In woodlands, intensity of fires is likely increased by the amount of fine fuel load [Hoffa et al., 1999], which may increase as tree cover increases because more fine litters fallen from trees are available for open fires. Thus the positive relationship in the equation (3) and negative in (4) may be the combined results of the two opposite factors of

fuel injections for fire intensity and of obstacles against fire spread. We note that the uncertainties associated with any calibration methods are likely to be comparatively large. The MODIS burned area of Roy et al. [2002] (0.81  106 km2) for the sum of July and September was smaller (30%) than our estimate (1.16  106 km2). This amount of difference was also found in the study between the different calibration techniques of Giglio et al. [2006a] (0.83  106 km2 yr1) and van der Werf et al. [2003] (1.16  106 km2 yr1). [15] In addition to the burned area observations, the errors in the estimated burned areas also depend on the accuracy of the fire count data. Thus we examine the sensitivity study to investigate the effect of different seasonal variability of fire counts on modeled CO emissions. In southern Africa, VIRS fire counts [Giglio et al., 2003] peak about two months earlier than the Along Track Scanning Radiometer (ATSR) [Arino and Plummer, 2001] and the MODIS [Justice et al., 2002]. We used the Collection 4, version 3 Terra MODIS monthly Climate Modeling Grid (CMG) fire product at a 0.5° spatial resolution from January 2001 through December 2005 [Giglio et al., 2006b]. The method used to produce the seasonal variations is based on the overlap period between the VIRS and the MODIS fire data sets. The monthly 0.5-degree burned areas and the MODIS fire counts were summed to annual data in southern Africa, and then the monthly variations of burned areas were scaled and distributed by the MODIS fire counts at 0.5-degree resolution. 2.2. Fuel Consumption and Emission Factors [16] Our fuel consumption model considers herbaceous vegetation in grasslands and that plus the fine litter of leaves fallen from trees in woodlands. The primary biomass density was derived from a 0.5-degree terrestrial carbon cycle model from 1998 to 2003. The ecosystem model used in this study is a modified version of Simulation model of Carbon cYCle in Land Ecosystems (Sim-CYCLE) [Ito and Oikawa, 2002; Ito et al., 2006], which is composed of 18 carbon pools representing the canopy tree layer, floor C3 grass layer, floor C4 grass layer, soil surface litter, and mineral soil. Carbon flows such as photosynthesis, respira-

3 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

tion, and decomposition are calculated using physiological relationships with moderate simplicity for the global-scale applications. The model parameters were calibrated using field data from a number of sites, so that the model validation shows a good agreement with observations of carbon storage and productivity [Ito and Oikawa, 2002]. The carbon pools for herbaceous living biomass and tree fine litter are taken from the calculations of the processbased model [Ito et al., 2006]. The majority of fires in southern Africa are surface fires and the woody debris is generally collected for use as fuels [Shea et al., 1996; Desanker et al., 1997; Hoffa et al., 1999]. Thus the coarse fuels (i.e., stems and coarse woody debris) are not treated in this work and implementation of this assumption is discussed in section 3. [17] Monthly meteorological data for the carbon cycle model were obtained from the United States National Centers for Environmental Prediction and the United States National Center for Atmospheric Research (NCEP/NCAR) [Kistler et al., 2001] for the period from 1948 to 2003. Since the meteorological data were available to 2003, we used the monthly averages of the output at each grid for 1998 – 2003 for the years of 2004 and 2005. Trenberth and Guillemot [1998] evaluated the NCEP/NCAR reanalysis precipitation data set. They found a reasonable agreement in annual mean precipitation over southern Africa, but excessive precipitation, for example, in June 1995. In future work, correction using observation-based precipitation data may help to improve the carbon cycle model [Sheffield et al., 2006], which is currently postprocessed by incorporating the satellite data of tree cover [Hansen et al., 2003] and the Global Precipitation Climatology Project (GPCP) product [Adler et al., 2003], as described below. [18] The herbaceous biomass density from the prognostic model generally showed good agreement during the wet season but overestimated during the dry season in southern Africa, mainly because the herbs continued to grow even during the dry season owing to the excessive precipitation from the reanalysis data. Thus we used the living biomass density at the end of growing season from the model for the fuel load of the herbs during the burning season. [19] When the herbaceous biomass density from the carbon cycle model was calculated to be zero owing to the light limitation for the under-story vegetation but the MODIS herbaceous cover was plus, the biomass density based on the relationship between the annual rainfall (mm) and in situ biomass density (kg DM m2) measurements at about 150 sites in West African savannas from Menaut et al. [1991] (slope: 4.9  104, intercept: 0.058, correlation coefficient: 0.7) was used to fill the gap in each grid. Thus we applied this relationship to the annual rainfall data from 1998 to 2005 at a 2.5-degree resolution using monthly rainfall data from Adler et al. [2003] [Ito and Penner, 2004]. [20] Fractional tree cover from the process-based model, which simulates potential vegetation, is generally higher than that derived from the satellite observation. Thus we scaled the biomass density for tree fine litter from the model by the ratio of tree cover between the model and the MODIS. We also scaled the biomass density for herbaceous

GB2011

vegetation from the model by the ratio of herbaceous cover between the model and the MODIS. [21] The combustion factors and emission factors are modeled as functions of grass fuel moisture and composition in grasslands and woodlands [Ward et al., 1996; Hoffa et al., 1999; Korontzi, 2005]. The fuel moisture is derived from 1-km satellite vegetation index time series. We used the relationship between the percentage of green grass out of the total grass (PGREEN) and NDVI [Korontzi, 2005; Verbesselt et al., 2006]. [22] The NDVI and cloud mask data set were derived from the sensor VEGETATION on board the SPOT satellite platforms by the Flemish Institute for Technological Research (VITO). The images had a spatial resolution of approximately 1 km2 after geometric, radiometric, and atmospheric corrections [Maisongrande et al., 2004]. Clouds obscure the land surface, reducing NDVI considerably. The cloud screening for VEGETATION products is based mainly on a threshold approach using the contrast between clouds and Earth surface in the blue band [Kempeneers et al., 2000]. We used the cloud mask data to identify cloud-free pixels. We also used evergreen forests and desert masks determined by NDVI to remove pixels where the NDVI-based method to assess the PGREEN was unreliable [Korontzi, 2005]. The cloud-free NDVI data in grasslands and woodlands are calculated from the 1998– 2005 time series of NDVI and cloud mask data. Then the relative greenness indices, PGREEN, are calculated using the following equation [Kogan, 1997; Burgan et al., 1998],  ½PGREENlit ¼ ½NDVIlit  ½NDVIlimin =  ½NDVIlimax  ½NDVIlimin ;

ð5Þ

where [NDVI]limin is the minimum of all NDVI values at the pixel over the preceding growing season, [NDVI]limax is the maximum of all NDVI values at the pixel over the preceding growing season. If fine-scale data are averaged over coarse-grid data from different land cover types, the results are often biased. Thus we calculated the averages from the cloud-free NDVI data over a 0.5-degree resolution in grasslands and woodlands, separately. Since the NDVI data were available from April 1998, we used the monthly averages of the PGREEN from 1999 to 2005 for the period from January to March in 1998. [23] For the fine fuels in grasslands and woodlands, we adopted combustion factors that depend on the fuel moisture. Thus we used the relationship between the combustion factors and the PGREEN [Hoffa et al., 1999]. Ito and Penner [2005b] recalculated the relationships between CF and PGREEN in grasslands (slope: 2.0, intercept: 1.4, correlation coefficient: 0.7) and woodlands (slope: 2.1, intercept: 0.87, correlation coefficient: 0.6), based on the measurements at 8 sites presented by Shea et al. [1996], Ward et al. [1996], Hoffa et al. [1999], and Korontzi et al. [2003]. The relationships include measurements obtained in both the early and late burning season. [24] Following Hoffa et al. [1999], Ito and Penner [2005b] related the modified combustion efficiency (MCE), which was defined as the molar ratio of carbon emitted as CO2 to the sum of carbon emitted as CO2 and CO

4 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

Figure 2. Comparison of PGREEN from open vegetation fires in southern Africa in (a) grasslands and (b) woodlands. The dashed line with circle represents the model estimates, while the solid line with error bar denotes the medians and ranges of the observations. (DCO2/(DCO2 + DCO)) [Ward and Radke, 1993], to PGREEN in grasslands (slope: 0.21, intercept: 1.0, correlation coefficient: 0.7). This formula was developed based on the measurements at 8 sites taken from Shea et al. [1996], Ward et al. [1996], Hoffa et al. [1999], and Korontzi et al. [2003]. [25] Fire extent and intensity in southern Africa are strongly influenced by biomass availability for fires, and hence rainfall [Lindesay et al., 1996; Swap et al., 2003]. Fires in woodlands tend to be more intense in areas of low tree cover and high mean annual rainfall, where grass production is high but where grass quality and therefore grazing pressure is low [Desanker et al., 1997]. The grasses dry faster and are consumed during the flaming combustion which increases the MCE, while the litter fuels dry slower and tend to involve more smoldering combustion, which can decrease the MCE [Ward et al., 1996; Bertschi et al., 2003]. The dependence of MCE on fuel types has been determined for southern Africa in the late dry season by Ward et al. [1996] and in both the early and late dry season by Korontzi [2005]. A reanalysis of the measurements in woodlands presented by Shea et al. [1996], Ward et al. [1996], Hoffa et al. [1999], and Korontzi et al. [2003] provide the following relationship for MCE as a function of fuel types,     ½MCEwit ¼ 1:25  Pgg wit þ 1:00  Pdg wit þ 0:89  ½Pfl wit  R2 ¼ 0:83 ; ð6Þ

where P denotes the proportion of each fuel type of green grass (gg), dead grass (dg) and fine litter (fl) in the fuel mixture and the three proportions sum to 1. These measurements took place at 9 sites in Zambia from June to September. The fuel loads of green and dead grasses were calculated from the NDVI-based PGREEN values and the fuel-load modeled estimates of grass [Korontzi, 2005]. [26] We used the regression model between MCE and EF for CO (g-CO kg-DM1, slope: 1134.6, intercept: 1135.4) from Ito and Penner [2004], based on the approach developed by Ward et al. [1996]. The correlation coefficient

between MCE and EF for CO is unity, because EF for CO is not an independent variable against MCE (R. Yokelson, personal communication, 2004). These measurements of MCE and EF for CO were taken from a wide variety of independent 39 fires varying from tropical to boreal zone fuels [Hurst et al., 1994a, 1994b; Ward et al., 1996; Goode et al., 2000; Korontzi et al., 2003; Shirai et al., 2003; Yokelson et al., 2003]. [27] We recalculated the regression equation between MCE and EF for BC (g-BC kg-DM1) from 10 fires measured in southern Africa by Sinha et al. [2003, 2004]. ½ EF ðBCÞit ¼ 9:66  ½MCEit þ 9:54

 R2 ¼ 0:65 :

ð7Þ

[28] To quantify the relative importance of different aspects of our model in explaining the seasonal delay of the CO emissions peak relative to burned area, we set two different scenarios. For case 1, the monthly averaged fuel loads, combustion factors, and emission factors were used to estimate the monthly CO emissions. The annually averaged effective CO emissions per area burned were calculated in grasslands (18 g CO m2) and woodlands (30 g CO m2) in southern Africa from the calculations of case 1. For case 2, these averages for each land cover type are multiplied by the monthly burned areas to estimate the monthly CO emissions.

3. Results and Discussion 3.1. Comparison of Relative Greenness Indices, Fuel Consumption, and Emission Factors With Observations [29] We compare our PGREEN, FC and EF model results with the field measurements as a function of time and land cover type. The comparisons of monthly variations of PGREEN in the area affected by open vegetation fires are shown in Figure 2. Monthly averaged observations were taken from Shea et al. [1996] and Hoffa et al. [1999] in Zambia. The median and range of observations are shown because of the limited number of the observations for different fires. Overall, our estimates of PGREEN were in good agreement with observations. The monthly averages

5 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

Figure 3. Comparison of fuel consumption (kg DM m2) from open vegetation fires in southern Africa in (a) grasslands and (b) woodlands. The dashed line with circle represents the model estimates, while the solid line with error bar denotes the medians and ranges of the observations. of PGREEN in grasslands (Figure 2a) decreased from June to July. From July to September, the dead grasses in grasslands were already dry and the monthly averages of PGREEN were almost constant in time. The monthly averages of PGREEN in woodlands (Figure 2b) decreased from the early to late burning season. [30] The comparisons of monthly variations of FC from open vegetation fires in southern Africa are shown in Figure 3. Monthly averaged observations were taken from Shea et al. [1996], Hoffa et al. [1999], and He´ly et al. [2003]. The ecology of miombo ecosystems is closely related to that of the savanna-like cerrado in South America [Desanker et al., 1997]. Thus the two measurements of fuel consumption in the cerrado [Ward et al., 1992] were also used for the averaged FC in woodlands. Our fuel consumption in grasslands (Figure 3a) was generally higher than the measurement. These overestimates in grasslands do not cause a significant error in the total emissions of CO because woodland fires (80%) are a much larger source of CO than are grassland fires (20%) in southern Africa. Our fuel consumption in woodlands (Figure 3b) increased from 0.24 (kg DM m2) in June to 0.68 in September, while the

measurement increased from 0.01– 0.25 to 0.45 – 0.84. The seasonal increase in FC is mainly caused by the variation in fuel moisture during the burning season [Hoffa et al., 1999]. At the beginning of the burning season, the fuel is wetter and burns less completely. As it dries out, it burns more completely. [31] The comparisons of monthly variations of EF for CO from open vegetation fires in southern Africa are shown in Figure 4. Monthly averaged observations were taken from Ward et al. [1992, 1996], Korontzi et al. [2003], and Yokelson et al. [2003]. Our emission factor in grasslands (Figure 4a) decreased from 74 (g-CO kg-DM1) in June to 48 in July, while that of measurements decreased from 52– 101 to 33– 54. The seasonal decrease in EF for grasslands is associated with the burning of moister vegetation at the beginning of the dry season [Korontzi et al., 2003]. The burning of moist vegetation produces a higher fraction of the products of incomplete combustion in the early season, which increases the production of CO emissions. Our emission factor in woodlands (Figure 4b) did not significantly change with month, which was consistent with the measurements in southern Africa. Our averaged EF in

Figure 4. Comparison of emission factors (g-CO kg-DM1) from open vegetation fires in southern Africa in (a) grasslands and (b) woodlands. The dashed line with circle represents the model estimates, while the solid line with error bar denotes the medians and ranges of the observations. 6 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

Figure 5. Seasonal and interannual variations of burned areas and CO emissions in southern Africa. The dashed line with circle represents the burned areas, while the solid line with square denotes the CO emissions.

GB2011

one pixel in late burning season are not differentiated owing to the coarser VIRS spatial resolution (4.8 km2 before 7 August 2001 and 5.8 km2 after 24 August 2001) [Kummerow et al., 1998] than MODIS (1 km2). [35] The fire season in southern Africa typically starts in May in the northwest, moves southeast, and ends in October near the east coast of South Africa [Cahoon et al., 1992]. The peaks in burned areas, combustion factors and CO emissions for 2000 are shown in Figure 7. The peaks in burned area occurred on the west side of the continent around between 5°S and 10°S in June and July, and on the east side in July – September. In contrast, the peaks in combustion factors and CO emissions were found on the west side of the continent in July – September, and the east side below 10°S in September. This characteristic delay between burned areas and CO emissions is mainly caused by significant changes in CF. In woodlands, the averaged combustion factor increased from 0.25 in June to 0.71 in September by a factor of 2.80, while the burned area

woodlands (69) was within the range of the observations for grassland fires in June and significantly smaller than the average for tropical forest fires (104 ± 20) [Andreae and Merlet, 2001]. This indicates that the woodland fires involve smoldering smokes as much as in early grassland fires but more flaming than the forest fires. [32] Our averaged BC emission factors in grasslands and woodlands were 0.31 and 0.45 (g-BC kg-DM1) and in good agreement with those of measurements (0.24 and 0.47) in southern Africa [Sinha et al., 2003, 2004] and with the global average (0.48 ± 0.18) [Andreae and Merlet, 2001]. We note that the measurement data is not available for comparisons of monthly variations of EF for BC. 3.2. Seasonal Variability [33] The monthly variations of burned areas and CO emissions are shown in Figure 5. The peaks in BA often occurred in early burning season of June and July, while those in CO emissions were found in late burning season of August and September. [34] The monthly averages of burned areas and CO emissions in different scenarios are shown in Figure 6. Two different fire counts of (Figure 6a) VIRS and (Figure 6b) MODIS were used to estimate the burned areas and CO emissions. The monthly averaged burned area from VIRS fire counts peaked earlier than MODIS, as well as monthly mean CO emissions when the CO emission were calculated from annually averaged effective CO emissions per area burned (case 2). In contrast, the peaks in monthly averaged CO emissions were identical between two different fire counts of VIRS and MODIS when the CO emissions were derived from time-dependent CO emissions per area burned (case 1). Even though more fire affected areas were observed in June and July by VIRS, the effective CO emissions per burned area were substantially low in the early burning season when the fires might cease owing to the wet fuels. These results suggest that the VIRS is sensitive to subpixel small fires as well as large ones, but large-scale fires within

Figure 6. Seasonal variations of burned areas and CO emissions in southern Africa. The dashed line with circle represents the burned areas, while the solid line with square denotes the CO emissions derived from time-dependent CO emissions per area burned (case 1) and the solid line with triangle shows the CO emissions from annually averaged CO emissions per area burned (case 2). Two different fire counts of (a) VIRS and (b) MODIS were used to estimate the burned areas.

7 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

ranged from 1.8  106 km2 for 2002 to 3.0  106 km2 for 1998, while those of CO and BC emissions ranged from 45 Tg CO and 0.26 Tg BC for 2002 to 75 Tg CO and 0.42 Tg BC for 1998. The annual burned area was generally coupled with the fire emissions. This is consistent with the result of van der Werf et al. [2006], but our averaged BA from 1998 to 2004 (2.3  106 km2) was significantly larger than their estimate (0.81  106 km2) from Giglio et al. [2006a]. Since both of BA were based on the calibration methods, the major difference in the total amounts might be associated with the burned area data sets of Roy et al. [2002] and Giglio et al. [2006a], both of which used the MODIS observations but different techniques. Appropriate validation is needed to assess the uncertainties associated with satellite-data-based products [Roy et al., 2005a]. [37] Barbosa et al. [1999] showed significantly smaller biomass burned after an El Nin˜o in 1986/1987. Their possible explanation was that the drought might have lead to a lower biomass production with a consequent reduction in the number of fires, burned areas and biomass burned. During the 1997/1998 episode, the typical drought that occurs during the ENSO warm events was not as severe and extensive in southern Africa [Anyamba et al., 2002]. Therefore the biomass production did not significantly change during the growing season for 1997/1998 in southern Africa. Our estimate of CO emission for 1998 was larger (24%) than the average, which was consistent with the result of van der Werf et al. [2006] (25%). 3.4. Carbon Monoxide Emissions [38] The annual amounts of total area burned, the effective CO emissions per area burned and CO emissions for southern Africa are compared among various modeling studies in Table 2. Our CO emission was in good agreement with those of van der Werf et al. [2006], although our annually averaged burned area was significantly larger by a factor of 2.8. Our annually averaged CO emission per area burned was 27 (g CO m2), which was in good agreement with the regional estimates of Sinha et al. [2004] (34) and Ito and Penner [2005b] (35). That of van der Werf et al. [2006] (104) was in better agreement with that of global emission model of Ito and Penner [2004] (130). The main reason is that in the global models, woody vegetation contributes substantially to emissions in savanna grid cells, while the fine fuel is a dominant component in the regional

Figure 7. Spatial distributions of the peaks in (a) burned areas, (b) combustion factors, and (c) CO emissions for the year 2000.

Table 1. Annual Amounts of Burned Areas, CO, and BC Emissions From Open Vegetation Fires in Southern Africa Year

[BA]a

[CO]b

[BC]c

(106 km2 yr1) decreased from 0.36 to 0.24. The fuel loads and the emission factors did not change significantly. As a result, CO emission (Tg CO yr1) increased from 4.9 to 11.9.

1998 1999 2000 2001 2002 2003 2004 2005

3.0 2.4 2.4 2.3 1.8 2.0 1.9 2.2

75 67 68 64 45 51 51 60

0.49 0.43 0.43 0.41 0.29 0.33 0.33 0.39

3.3. Interannual Variability [36] The interannual variations of BA, CO and BC emissions are presented in Table 1. The annual estimates of BA from open vegetation burning in southern Africa

8 of 13

Burned areas (106 km2 yr1). CO emissions (Tg CO yr1). c BC emissions (Tg BC yr1). a

b

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

Table 2. Comparison of Annual Amounts of Areas Burned, CO Emissions per Area Burned and CO Emissions From Open Vegetation Fires in Southern Africa Among Various Modeling Studies Study

[BA]a

[CO]/[BA]b

Scholes et al. [1996a, 1996b] Hoelzemann et al. [2004] Ito and Penner [2004] Palacios-Orueta et al. [2004] Sinha et al. [2004] Korontzi [2005] van der Werf et al. [2006]d This worke Ito and Penner [2005b] Arellano et al. [2006] Pe´tron et al. (manuscript in preparation, 2007)

1.68 0.58 1.00 1.08 0.77 1.04 0.79 2.26

9 119 130 202 34 11 104 27 35

[CO]c 15 69 130 218 26 12 82 60 124 116 – 150f

Burned areas (106 km2 yr1). Effective CO emissions per area burned (g CO m2). c CO emissions (Tg CO yr1). d The 1997 – 2004 average. e The 1998 – 2005 average. f Range during 2000 – 2003. a

b

models. Thus the average fuel loads in the global models are higher than the values in regional models, which are based on the regional measurements [Shea et al., 1996; Hoffa et al., 1999]. The fuel wood collection is accounted for the global models in the coarse grid cells, but this human

GB2011

practice must be treated as sub-grid-scale phenomena, as Shea et al. [1996] and Hoffa et al. [1999] reported that the majority of coarse woody debris (CWD) had been previously removed for use as fuels prior to the fires by local farmers. Satellite measurements of vegetation changes together with the burned areas may provide better space-andtime-resolved data in the future. [39] In order to examine differences in bottom-up estimates of CO emissions from different data sets with those from top-down estimates, we compare the estimates of CO to the ranges deduced from the inverse modeling studies [Arellano et al., 2006] (also G. Pe´tron et al., Multiyear inversion of CO surface sources using the MOPITT satellite data, manuscript in preparation, 2007) (hereinafter referred to as Pe´tron et al., manuscript in preparation, 2007). In the Arellano et al. [2006] and Pe´tron et al. (manuscript in preparation, 2007) studies, the monthly amounts of the prescribed CO emissions were optimized, so that the predicted CO mixing ratios from a chemical transport model fit the measured monthly average CO mixing ratios derived from the MOPITT instrument. Biofuel emissions were separated from the open vegetation fire category in these studies and thus were also not included here. [40] The comparisons of monthly variations of CO emissions are shown in Figure 8. The seasonal variations of our CO emissions, which were derived from the burned areas using VIRS fire counts and time-dependent effective CO

Figure 8. Monthly estimates of CO emissions for (a) 2000, (b) 2001, (c) 2002, and (d) 2003. The solid with circle line denotes our estimates and the solid with triangle from van der Werf et al. [2006], while the dashed line represents those from Arellano et al. [2006] and the dashed line with square from Pe´tron et al. (manuscript in preparation, 2007). 9 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

GB2011

Africa for 2003 and 2004. Accurate measurements of the burned areas are required to improve the emission estimates in southern Africa.

Figure 9. Seasonal variations derived from different sources for the year 2000. (a) The solid line with circle denotes our estimates of BC emissions and the solid with triangle from van der Werf et al. [2006], while the dashed line with square represents those from Penner et al. [2004]. (b) The dashed line with square denotes TOMS AI, while the solid with triangle represents MODIS AOD. emissions per area burned, were better agreement with those from previous top-down estimates using the MOPITT data for 2000 and 2001. Although our annual emissions were smaller than the top-down estimates, the smaller peaks derived from the observations from January to March were likely influenced by biomass burning in the Northern Hemisphere Africa. Thus sum of monthly emissions from May to September was calculated for the comparison. Better agreement was obtained among the four estimates of Arellano et al. [2006], Pe´tron et al. (manuscript in preparation, 2007), van der Werf et al. [2006], and this work (74 ± 14 Tg CO for 2000). In later years, however, the agreement was worse (75 ± 38 Tg CO for 2003). Our calibration methods for the burned area estimates are based on the data sets solely for 2000. Thus the extrapolation to different years causes additional uncertainties in later years. van der Werf et al. [2006] used the Giglio et al. [2006a] data set, which was extrapolated from the selected areas and months of the prototype MODIS burned areas between January 2001 and December 2004 to fourteen geographic regions during the period. Gloudemans et al. [2006] compared the SCIAMACHY CO column with the model CO estimates using the van der Werf et al. [2006] data set and found significant underestimates of the model results in southern

3.5. Black Carbon Emissions [41] We compare our estimates of BC emissions from open biomass burning for 2000 to those deduced from the inverse modeling study [Penner et al., 2004] in regions of southern Africa defined by Zhang et al. [2005]. We note that the values shown in Figure 9a is smaller (about 10%) than those presented in section 3.1 (Table 1) because the land area in the domain defined by Zhang et al. [2005] is slightly smaller than the area in the region defined in this paper as the Southern Hemisphere Africa. In the Penner et al. [2004] study, the monthly amounts of the prescribed BC emissions from open vegetation fires were optimized using the method developed by Zhang et al. [2005], based on the Total Ozone Mapping Spectrometer (TOMS) Aerosol Index (AI), a 3-D chemical transport model (IMPACT), and the Bayesian inversion approach [Tarantola and Valette, 1982a, 1982b]. The TOMS AI is a measure of the change of spectral contrast in the near ultraviolet due to the radiative transfer effects of aerosols in a Rayleigh scattering atmosphere [Herman et al., 1997]. The IMPACT model was developed at the Lawrence Livermore National Laboratories [Rotman et al., 2004] and extended to treat aerosols and detailed chemical reactions for a wider set of volatile organic compounds at the University of Michigan [Liu et al., 2005; Feng and Penner, 2007; Ito et al., 2007]. [42] The monthly optimized emissions using a priori data set of Ito and Penner [2005a] are shown in Figure 9a. Sum of monthly emissions from May to September was in good agreement among the three estimates of Penner et al [2004], van der Werf et al. [2006], and this work (0.46 ± 0.08 Tg BC yr1). Although our BC emissions were in good agreement with those of van der Werf et al. [2006], our annually averaged BC emission per area burned was 0.17 (g BC m2), which was in good agreement with the regional estimate of Sinha et al. [2004] (0.22). [43] Here we show the monthly averaged TOMS AI data in Figure 9b, corresponding to absorbing aerosols (0.2 < AI) [Torres et al., 2002]. For each month, we used the total number of days of available data to calculate the weighted monthly averages. We also show the monthly averaged MODIS Aerosol Optical Depth (AOD) data [Remer et al., 2005]. The peak in our BC emissions during open biomass burning season for 2000 was identical to those from previous top-down estimate and the TOMS AI data. [44] Our estimated amount of peak emission of BC was in good agreement with the inverse modeling estimate, but we underestimated CO (about 30%). Although this difference is close to the uncertain range of inverse estimates [Arellano et al., 2006; Pe´tron et al., manuscript in preparation, 2007], this systematic bias may be associated with the assumption of no coarse fuels burning in this work. Generally, the coarse fuels are consumed during the smoldering process in savanna [Ward et al., 1996; Bertschi et al., 2003]. Because the emissions from the prolonged, smoldering fires associated with coarse woody debris are injected near the surface, the transport of these emissions from the surface to the

10 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

altitude where the satellite instrument has enough sensitivity to measure the mixing ratios associated with these emissions is particularly important [Ito and Penner, 2005a; Arellano and Hess, 2006]. Especially, the TOMS AI may not detect smoke plumes from low-intensity fires [Duncan et al., 2003].

4. Conclusions [45] We estimated the emissions of CO and BC from open vegetation fires in southern Africa over the 1998 – 2005 period, using satellite derived information and the biogeochemical model. Overall, our estimates of relative greenness indices, fuel consumption and emission factors from June to September were in good agreement with those of the measurements. [46] This study implies that process-based carbon cycle models may reasonably retrieve the distributions of fuels in natural land ecosystems. However, these models were not sufficiently validated for the spatial and temporal variability in fuel load, which is actually affected by fire history and human activities. Incorporating the impacts of disturbance is a challenging issue for carbon cycle model study in conjunction with observational study; especially, wild fire regime is interactive with the carbon cycle through fuel accumulation and consumption. [47] The monthly averaged burned area from VIRS fire counts peaked earlier in June and July, than those in CO emissions in August and September. This characteristic delay between burned areas and emissions was mainly explained by significant changes in combustion factors for woodlands. Even though more fire affected areas were observed in the early burning season by VIRS, the averaged combustion factor was substantially low owing to the wet fuels. [48] This is the first bottom-up emission model which shows that the peaks in CO and BC emissions during open biomass burning season for 2000 were identical to those from previous top-down estimates using the MOPITT and the TOMS AI data [Penner et al., 2004; Arellano et al., 2006; Pe´tron et al., manuscript in preparation, 2007]. Sum of monthly emissions during burning season from May to September in 2000 was in good agreement between our bottom-up estimates and previous top-down estimates. In later years, however, the bottom-up estimates were lower than the top-down estimates. Clearly, further work is required to reduce the differences between these two approaches.

[49] Acknowledgments. The SPOT VEGETATION data sets were distributed by Spot Image and produced by VITO, Belgium. We wish to thank the Oak Ridge National Laboratory Distributed Active Archive Center at Goddard Space Flight Center, for producing the SAFARI 2000 data in its present format and distributing them.

References Adler, R. F., et al. (2003), The Version-2 Global Precipitation Climatology Project (GPCP) monthly precipitation analysis (1979 – present), J. Hydrometeorol., 4, 1147 – 1167. Andreae, M. O., and P. J. Crutzen (1997), Atmospheric aerosols: Biogeochemical sources and role in atmospheric chemistry, Science, 276, 1052 – 1058. Andreae, M. O., and P. Merlet (2001), Emission of trace gases and aerosols from biomass burning, Global Biogeochem. Cycles, 15, 955 – 966.

GB2011

Anyamba, A., C. J. Tucker, and R. Mahoney (2002), From El Nin˜o to La Nin˜a: Vegetation response patterns over east and southern Africa during the 1997 – 2000 period, J. Clim., 15, 3096 – 3103. Arellano, A. F., Jr., and P. G. Hess (2006), Sensitivity of top-down estimates of CO sources to GCTM transport, Geophys. Res. Lett., 33, L21807, doi:10.1029/2006GL027371. Arellano, A. F., Jr., P. S. Kasibhatla, L. Giglio, G. R. van der Werf, and J. T. Randerson (2004), Top-down estimates of global CO sources using MOPITT measurements, Geophys. Res. Lett., 31, L01104, doi:10.1029/ 2003GL018609. Arellano, A. F., Jr., P. S. Kasibhatla, L. Giglio, G. R. van der Werf, J. T. Randerson, and G. J. Collatz (2006), Time-dependent inversion estimates of global biomass-burning CO emissions using Measurement of Pollution in the Troposphere (MOPITT) measurements, J. Geophys. Res., 111, D09303, doi:10.1029/2005JD006613. Arino, O., and S. Plummer (2001), Along Track Scanning Radiometer World Fire Atlas: Validation of the 1997 – 98 Active Fire Product, Eur. Space Agency, Eur. Space Res. Inst., Frascati, Italy. Barbosa, P. M., D. Stroppiana, J. Gre´goire, and J. M. C. Pereira (1999), An assessment of vegetation fire in Africa (1981 – 1991): Burned areas, burned biomass, and atmospheric emissions, Global Biogeochem. Cycles, 13, 933 – 950. Bertschi, I., R. J. Yokelson, D. E. Ward, R. E. Babbitt, R. A. Susott, J. G. Goode, and W. M. Hao (2003), Trace gas and particle emissions from fires in large diameter and belowground biomass fuels, J. Geophys. Res., 108(D13), 8472, doi:10.1029/2002JD002100. Boschetti, L., H. D. Eva, P. A. Brivio, and J. M. Gre´goire (2004), Lessons to be learned from the comparison of three satellite-derived biomass burning products, Geophys. Res. Lett., 31, L21501, doi:10.1029/ 2004GL021229. Burgan, R. E., R. W. Klaver, and J. M. Klaver (1998), Fuel models and fire potential from satellite and surface observations, Int. J. Wildland Fires, 8, 159 – 170. Cahoon, J., B. Stocks, J. Levine, W. Cofer, and K. O’Neill (1992), Seasonal distribution of African savannah fires, Nature, 359, 812 – 815. Carmona-Moreno, C., A. Belward, J.-P. Malingreau, A. Hartley, M. GarciaAlegre, M. Antonovskiy, V. Buchshtaber, and V. Pivovarov (2005), Characterizing interannual variations in global fire calendar using data from Earth observing satellites, Global Change Biol., 11, 1537 – 1555. Chuvieco, E., and E. S. Kasischke (2007), Remote sensing information for fire management and fire effects assessment, J. Geophys. Res., 112, G01S90, doi:10.1029/2006JG000230. Crutzen, P. J., and M. O. Andreae (1990), Biomass burning in the tropics: Impact on atmospheric chemistry and biogeochemical cycles, Science, 250, 1669 – 1678. Crutzen, P. J., L. E. Heidt, J. P. Krasnec, W. H. Pollock, and W. Seiler (1979), Biomass burning as a source of atmospheric gases CO, H2, N2O, NO, CH3Cl and COS, Nature, 282, 253 – 256. Desanker, P. V., P. G. H. Frost, C. O. Justice, and R. J. Schloes (Eds.) (1997), The Miombo Network: Framework for a terrestrial transect study of land-use and land-cover change in the Miombo ecosystems of Central Africa, conclusions of the Miombo Network Workshop, Rep. 41, 109 pp., Int. Geosphere-Biosphere Programme, Stockholm. Duncan, B. N., R. V. Martin, A. C. Staudt, R. Yevich, and J. A. Logan (2003), Interannual and seasonal variability of biomass burning emissions constrained by satellite observations, J. Geophys. Res., 108(D2), 4100, doi:10.1029/2002JD002378. Edwards, D. P., et al. (2006), Satellite-observed pollution from Southern Hemisphere biomass burning, J. Geophys. Res., 111, D14312, doi:10.1029/2005JD006655. Feng, Y., and J. E. Penner (2007), Global modeling of nitrate and ammonium: Interaction of aerosols and tropospheric chemistry, J. Geophys. Res., 112, D01304, doi:10.1029/2005JD006404. Giglio, L., J. D. Kendall, and R. Mack (2003), A multi-year active fire dataset for the tropics derived from the TRMM VIRS, Int. J. Remote Sens., 24, 4505 – 4525. Giglio, L., G. R. van der Werf, J. T. Randerson, G. J. Collatz, and P. Kasibhatla (2006a), Global estimation of burned area using MODIS active fire observations, Atmos. Chem. Phys., 6, 957 – 974. Giglio, L., I. Csiszar, and C. O. Justice (2006b), Global distribution and seasonality of active fires as observed with the Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) sensors, J. Geophys. Res., 111, G02016, doi:10.1029/2005JG000142. Gloudemans, A. M. S., M. C. Krol, J. F. Meirink, A. T. J. de Laat, G. R. van der Werf, H. Schrijver, M. M. P. van den Broek, and I. Aben (2006), Evidence for long-range transport of carbon monoxide in the Southern Hemisphere from SCIAMACHY observations, Geophys. Res. Lett., 33, L16807, doi:10.1029/2006GL026804.

11 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

Goode, J. G., R. J. Yokelson, D. E. Ward, R. A. Susott, R. E. Babbitt, M. A. Davies, and W. M. Hao (2000), Measurements of excess O3, CO2, CO, CH4, C2H4, C2H2, HCN, NO, NH3, HCOOH, CH3COOH, HCHO, and CH3OH in 1997 Alaskan biomass burning plumes by airborne Fourier transform infrared spectroscopy (AFTIR), J. Geophys. Res., 105, 22,147 – 22,166. Hansen, M. C., R. S. DeFries, J. R. G. Townshend, R. Sohlberg, C. Dimiceli, and M. Carroll (2002), Towards an operational MODIS continuous field of percent tree cover algorithm: Examples using AVHRR and MODIS data, Remote Sens. Environ., 83, 303 – 319. Hansen, M. C., R. S. DeFries, J. R. G. Townshend, M. Carroll, C. Dimiceli, and R. A. Sohlberg (2003), Global percent tree cover at a spatial resolution of 500 meters: First results of the MODIS vegetation continuous fields algorithm, Earth Interact., 7(10), doi:10.1175/1087-3562 (2003)0072.0.CO;2. Hansen, M. C., J. R. G. Townshend, R. S. DeFries, and M. Carroll (2006), Estimation of tree cover using MODIS data at global, continental and regional/local scales, Int. J. Remote Sens., 26, 4359 – 4380. He´ly, C., S. Alleaume, R. J. Swap, H. H. Shugart, and C. O. Justice (2003), SAFARI-2000 characterization of fuels, fire behavior, combustion completeness, and emissions from experimental burns in infertile grass savannas in western Zambia, J. Arid Environ., 54, 381 – 394. Herman, J. R., P. K. Bhartia, O. Torres, C. Hsu, C. Seftor, and E. Celarier (1997), Global distribution of UV-absorbing aerosols from Nimbus-7/ TOMS data, J. Geophys. Res., 102, 16,911 – 16,922. Hoelzemann, J. J., M. G. Schultz, G. P. Brasseur, C. Granier, and M. Simon (2004), The Global Wildland Fire Emission Model GWEM: Evaluating the use of global area burnt data, J. Geophys. Res., 109, D14S04, doi:10.1029/2003JD003666. Hoffa, E. A., D. E. Ward, W. M. Hao, R. A. Susott, and R. H. Wakimoto (1999), Seasonality of carbon emissions from biomass burning in a Zambian savanna, J. Geophys. Res., 104, 13,841 – 13,853. Hurst, D. F., D. W. T. Griffith, and G. D. Cook (1994a), Trace gas emissions and biomass burning in tropical Australian savannas, J. Geophys. Res., 99, 16,441 – 16,456. Hurst, D. F., D. W. T. Griffith, J. N. Carras, D. J. Williams, and P. J. Fraser (1994b), Measurements of trace gas emitted by Australian savanna fires during the 1990 dry season, J. Atmos. Chem., 18, 33 – 56. Ito, A. H., and T. Oikawa (2002), A simulation model of the carbon cycle in land ecosystems (Sim-CYCLE): A description based on dry-matter production theory and plot-scale validation source, Ecol. Modell., 151, 143 – 176. Ito, A., and J. E. Penner (2004), Global estimates of biomass burning emissions based on satellite imagery for the year 2000, J. Geophys. Res., 109, D14S05, doi:10.1029/2003JD004423. Ito, A., and J. E. Penner (2005a), Historical emissions of carbonaceous aerosols from biomass and fossil fuel burning for the period 1870 – 2000, Global Biogeochem. Cycles, 19, GB2028, doi:10.1029/ 2004GB002374. Ito, A., and J. E. Penner (2005b), Estimates of CO emissions from open biomass burning in southern Africa for the year 2000, J. Geophys. Res., 110, D19306, doi:10.1029/2004JD005347. Ito, A. H., H. Muraoka, H. Koizumi, N. Saigusa, S. Murayama, and S. Yamamoto (2006), Seasonal variation in leaf properties and ecosystem carbon budget in a cool-temperate deciduous broad-leaved forest: Simulation analysis at Takayama site, Japan, Ecol. Res., 21, 137 – 149. Ito, A., S. Sillman, and J. E. Penner (2007), Effects of additional nonmethane volatile organic compounds, organic nitrates and direct emissions of oxygenated organic species on global tropospheric chemistry, J. Geophys. Res., 112, D06309, doi:10.1029/2005JD006556. Jain, A. K., Z. Tao, X. Yang, and C. Gillespie (2006), Estimates of global biomass burning emissions for reactive greenhouse gases (CO, NMHCs, and NOx) and CO2 , J. Geophys. Res., 111, D06304, doi:10.1029/ 2005JD006237. Justice, C. O., L. Giglio, S. Korontzi, J. Owens, J. T. Morisette, D. Roy, J. Descloitres, S. Alleaume, F. Petitcolin, and Y. Kaufman (2002), The MODIS fire products, Remote Sens. Environ., 83, 244 – 262. Kasischke, E. S., J. H. Hewson, B. Stocks, G. van der Werf, and J. Randerson (2003), The use of ATSR active fire counts for estimating relative patterns of biomass burning—A study from the boreal forest region, Geophys. Res. Lett., 30(18), 1969, doi:10.1029/2003GL017859. Kempeneers, P., G. Lissens, F. Fierens, and J. van Rensbergen (2000), Development of a cloud, snow and cloud shadow mask for VEGETATION imagery, in Proceedings of the Symposium VEGETATION 2000: Two Years of Operation for Preparation for the Future, edited by G. Saint, pp. 303 – 306, Joint Res. Cent., Eur. Comm., Ispera, Italy. Kistler, R., et al. (2001), The NCEP-NCAR 50-year reanalysis: Monthly means CD-ROM and documentation, Bull. Am. Meteorol. Soc., 82, 247 – 267.

GB2011

Kogan, F. N. (1997), Global drought watch from space, Bull. Am. Meteorol. Soc., 78, 621 – 636. Korontzi, S. (2005), Seasonal patterns in biomass burning emissions from southern African vegetation fires for the year 2000, Global Change Biol., 11, 1680 – 1700. Korontzi, S., D. E. Ward, R. A. Susott, R. J. Yokelson, C. O. Justice, P. V. Hobbs, E. A. H. Smithwick, and W. M. Hao (2003), Seasonal variation and ecosystem dependence of emission factors for selected trace gases and PM2.5 for southern African savanna fires, J. Geophys. Res., 108(D24), 4758, doi:10.1029/2003JD003730. Kummerow, C., W. Barnes, T. Kozu, J. Shiue, and J. Simpson (1998), The Tropical Rainfall Measuring Mission (TRMM) sensor package, J. Atmos. Oceanic Technol., 15, 809 – 817. Lindesay, J. A., M. O. Andreae, J. G. Goldammer, G. Harris, H. J. Annegarn, M. Garstang, R. J. Scholes, and B. W. van Wilgen (1996), International Geosphere-Biosphere Programme/International Global Atmospheric Chemistry SAFARI-92 field experiment: Background and overview, J. Geophys. Res., 101, 23,521 – 23,530. Liu, X., J. E. Penner, and M. Herzog (2005), Global modeling of aerosol dynamics: Model description, evaluation, and interactions between sulfate and nonsulfate aerosols, J. Geophys. Res., 110, D18206, doi:10.1029/ 2004JD005674. Maisongrande, P., B. Duchemin, and G. Dedieu (2004), VEGETATION/ SPOT: An operational mission for the Earth monitoring; presentation of new standard products, Int. J. Remote Sens., 25, 9 – 14. Menaut, J. C., L. Abbadie, F. Lavenu, P. Loudjani, and A. Podaire (1991), Biomass burning in west African savannas, in Global Biomass Burning, Atmospheric, Climatic, and Biospheric Implications, edited by J. S. Levine, pp. 133 – 142, MIT Press, Cambridge, Mass. Palacios-Orueta, A., A. Parra, E. Chuvieco, and C. Carmona-Moreno (2004), Remote sensing and geographic information systems methods for global spatiotemporal modeling of biomass burning emissions: Assessment in the African continent, J. Geophys. Res., 109, D14S09, doi:10.1029/2004JD004734. Penner, J. E., S. Zhang, and A. Ito (2004), Estimates of black carbon emissions from open biomass burning, Eos Trans. AGU, 85(47), Fall Meet. Suppl., Abstract A12B-01. Pe´tron, G., C. Granier, B. Khattatov, V. Yudin, J. Lamarque, L. Emmons, J. Gille, and D. P. Edwards (2004), Monthly CO surface sources inventory based on the 2000 – 2001 MOPITT satellite data, Geophys. Res. Lett., 31, L21107, doi:10.1029/2004GL020560. Plummer, S., O. Arino, M. Simon, and W. Steffen (2006), Establishing an Earth observation product service for the terrestrial carbon community: The GLOBCARBON initiative, Mitigation Adapt. Strat. Global Change, 11, 97 – 111. Remer, L. A., et al. (2005), The MODIS aerosol algorithm, products and validation, J. Atmos. Sci., 62, 947 – 973. Rotman, D. A., et al. (2004), IMAPCT, the LLNL 3-D global atmospheric chemical transport model for the combined troposphere and stratosphere: Model description and analysis of ozone and other trace gases, J. Geophys. Res., 109, D04303, doi:10.1029/2002JD003155. Roy, D. P. (2003), SAFARI 2000 July and September MODIS 500 m burned area products for Southern Africa, in SAFARI 2000 CD-ROM Series [CD-ROM], vol. 3, edited by J. Nickeson et al., Natl. Aeronaut. and Space Admin., Goddard Space Flight Cent., Greenbelt, Md. Roy, D. P., P. E. Lewis, and C. O. Justice (2002), Burned area mapping using multi-temporal moderate spatial resolution data—A bi-directional reflectance model-based expectation approach, Remote Sens. Environ., 83, 263 – 286. Roy, D. P., et al. (2005a), The Southern Africa Fire Network (SAFNet) regional burned area product validation protocol, Int. J. Remote Sens., 26, 4265 – 4292. Roy, D. P., Y. Jin, P. E. Lewis, and C. O. Justice (2005b), Prototyping a global algorithm for systematic fire-affected area mapping using MODIS time series data, Remote Sens. Environ., 97, 137 – 162. Scholes, R. J., J. Kendall, and C. O. Justice (1996a), The quantity of biomass burned in southern Africa, J. Geophys. Res., 101, 23,677 – 23,676. Scholes, R. J., D. E. Ward, and C. O. Justice (1996b), Emissions of trace gases and aerosol particles due to vegetation burning in Southern Hemisphere Africa, J. Geophys. Res., 101, 23,677 – 23,682. Seiler, W., and P. J. Crutzen (1980), Estimates of gross and net fluxes of carbon between the biosphere and the atmosphere from biomass burning, Clim. Change, 2, 207 – 247. Shea, R. W., B. W. Shea, J. B. Kauffman, D. E. Ward, C. I. Haskins, and M. C. Scholes (1996), Fuel biomass and combustion factors associated with fires in savanna ecosystems of South Africa and Zambia, J. Geophys. Res., 101, 23,551 – 23,568.

12 of 13

GB2011

ITO ET AL.: SEASONAL CO/BC EMISSIONS

Sheffield, J., G. Goteti, and E. F. Wood (2006), Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling, J. Clim., 19, 3088 – 3111. Shirai, T., et al. (2003), Emission estimates of selected volatile organic compounds from tropical savanna burning in northern Australia, J. Geophys. Res., 108(D3), 8406, doi:10.1029/2001JD000841. Silva, J. M. N., J. M. C. Pereira, A. I. Cabral, A. C. L. Sa´, M. J. P. Vasconcelos, B. Mota, and J. Gre´goire (2003), An estimate of the area burned in southern Africa during the 2000 dry season using SPOTVEGETATION satellite data, J. Geophys. Res., 108(D13), 8498, doi:10.1029/2002JD002320. Simon, M., S. Plummer, F. F. Fierens, J. J. Hoelzemann, and O. Arino (2004), Burnt area detection at global scale using ATSR-2; The GLOBSCAR products and their qualifications, J. Geophys. Res., 109, D14S02, doi:10.1029/2003JD003622. Sinha, P., P. V. Hobbs, R. J. Yokelson, I. T. Bertschi, D. R. Blake, I. J. Simpson, S. Gao, T. W. Kirchstetter, and T. Novakov (2003), Emissions of trace gases and particles from savanna fires in southern Africa, J. Geophys. Res., 108(D13), 8487, doi:10.1029/2002JD002325. Sinha, P., P. V. Hobbs, R. J. Yokelson, D. R. Blake, S. Gao, and T. W. Kirchstetter (2004), Emissions from miombo woodland and dambo grassland savanna fires, J. Geophys. Res., 109, D11305, doi:10.1029/ 2004JD004521. Swap, R. J., H. J. Annegarn, J. T. Suttles, M. D. King, S. Platnick, J. L. Privette, and R. J. Scholes (2003), Africa burning: A thematic analysis of the Southern African Regional Science Initiative (SAFARI 2000), J. Geophys. Res., 108(D13), 8465, doi:10.1029/2003JD003747. Tansey, K., et al. (2004), Vegetation burning in the year 2000: Global burned area estimates from SPOT VEGETATION data, J. Geophys. Res., 109, D14S03, doi:10.1029/2003JD003598. Tarantola, A., and B. Valette (1982a), Generalized nonlinear inverse problems solved using the least square criterion, Rev. Geophys., 20, 219 – 232. Tarantola, A., and B. Valette (1982b), Inverse problems—Quest of information, J. Geophys., 50, 159 – 170. Torres, O., P. K. Bhartia, J. R. Herman, A. Sinyuk, P. Ginoux, and B. Holben (2002), A long-term record of aerosol optical depth from TOMS observations and comparison to AERONET measurements, J. Atmos. Sci., 59, 398 – 413. Trenberth, K. E., and C. J. Guillemot (1998), Evaluation of the atmospheric moisture and hydrological cycle in the NCEP/NCAR reanalyses, Clim. Dyn., 14, 213 – 231.

GB2011

Tucker, C. J. (1979), Red and photographic infrared linear combinations for monitoring, Remote Sens. Environ., 8, 127 – 150. van der Werf, G. R., J. T. Randerson, G. J. Collatz, and L. Giglio (2003), Carbon emissions from fires in tropical and subtropical ecosystems, Global Change Biol., 9, 547 – 562. van der Werf, G. R., J. T. Randerson, L. Giglio, G. J. Collatz, P. S. Kasibhatla, and A. F. Arellano Jr. (2006), Interannual variability in global biomass burning emissions from 1997 to 2004, Atmos. Chem. Phys., 6, 3423 – 3441. Verbesselt, J., B. Somers, J. van Aardt, I. Jonckheere, and P. Coppin (2006), Monitoring herbaceous biomass and water content with SPOT VEGETATION time-series to improve fire risk assessment in savanna ecosystems, Remote Sens. Environ., 101, 399 – 414. Ward, D. E., and L. F. Radke (1993), Emissions measurements from vegetation fires: A comparative evaluation of methods and results, in Fire in the Environment: The Ecological, Atmospheric and Climatic Importance of Vegetation Fires, edited by P. J. Crutzen and J. G. Goldammer, pp. 53 – 76, John Wiley, Hoboken, N. J. Ward, D. E., R. A. Susott, J. B. Kauffman, R. E. Babbitt, D. L. Cummings, B. Dias, B. N. Holden, Y. J. Kaufman, R. A. Rasmussen, and A. W. Setzer (1992), Smoke and fire characteristics for cerrado and deforestation burns in Brazil: BASE-B experiment, J. Geophys. Res., 97, 14,601 – 14,619. Ward, D. E., W. M. Hao, R. A. Susott, R. E. Babbitt, R. W. Shea, J. B. Kauffman, and C. O. Justice (1996), Effect of fuel composition on combustion efficiency and emission factors for African savanna ecosystems, J. Geophys. Res., 101, 23,569 – 23,576. Yokelson, R. J., I. T. Bertschi, T. J. Christian, P. V. Hobbs, D. E. Ward, and W. M. Hao (2003), Trace gas measurements in nascent, aged, and cloudprocessed smoke from African savanna fires by airborne Fourier transform infrared spectroscopy (AFTIR), J. Geophys. Res., 108(D13), 8478, doi:10.1029/2002JD002322. Zhang, S., J. E. Penner, and O. Torres (2005), Inverse modeling of biomass burning emissions using Total Ozone Mapping Spectrometer aerosol index for 1997, J. Geophys. Res., 110, D21306, doi:10.1029/ 2004JD005738.



H. Akimoto, A. Ito, and A. Ito, Frontier Research Center for Global Change, JAMSTEC, 3173-25 Showa-machi, Kanazawa-ku, Yokohama, Kanagawa 236-0001, Japan. ([email protected])

13 of 13

Suggest Documents