Spatial and temporal variability in the ratio of trace gases emitted from biomass burning

Fires are a major source of trace gases and aerosols to the atmosphere. The amount of biomass burned is becoming better known, most importantly due to improved burned area datasets and a better representation of fuel consumption. The spatial and temporal variability in the partitioning of biomass burned into emitted trace gases and aerosols, however, has received relatively little attention. To convert estimates of biomass burned to trace gas and aerosol emissions, most studies have used emission ratios (or emission factors (EFs)) based on the arithmetic mean of field measurement outcomes, stratified by biome. However, EFs vary substantially in time and space, even within a single biome. In addition, it is unknown whether the available field measurement locations provide a representative sample for the various biomes. Here we used the available body of EF literature in combination with satellite-derived information on vegetation characteristics and climatic conditions to better understand the spatio-temporal variability in EFs. While focusing on CO, CH<sub>4</sub>, and CO<sub>2</sub>, our findings are also applicable to other trace gases and aerosols. We explored relations between EFs and different measurements of environmental variables that may correlate with part of the variability in EFs (tree cover density, vegetation greenness, temperature, precipitation, and the length of the dry season). Although reasonable correlations were found for specific case studies, correlations based on the full suite of available measurements were lower and explained about 33%, 38%, 19%, and 34% of the variability for respectively CO, CH<sub>4</sub>, CO<sub>2</sub>, and the Modified Combustion Efficiency (MCE). This may be partly due to uncertainties in the environmental variables, differences in measurement techniques for EFs, assumptions on the ratio between flaming and smoldering combustion, and incomplete information on the location and timing of EF measurements. We derived new mean EFs, using the relative importance of each measurement location with regard to fire emissions. These weighted averages were relatively similar to the arithmetic mean. When using relations between the environmental variables and EFs to extrapolate to regional and global scales, we found substantial differences, with for savannas 13% and 22% higher CO and CH<sub>4</sub> EFs than the arithmetic mean of the field studies, possibly linked to an underrepresentation of woodland fires in EF measurement locations. We argue that from a global modeling perspective, future measurement campaigns could be more beneficial if measurements are made over the full fire season, and if relations between ambient conditions and EFs receive more attention

Global estimation of burned area using MODIS active fire observations

We present a method for estimating monthly burned area globally at 1° spatial resolution using Terra MODIS data and ancillary vegetation cover information. Using regression trees constructed for 14 different global regions, MODIS active fire observations were calibrated to burned area estimates derived from 500-m MODIS imagery based on the assumption that burned area is proportional to counts of fire pixels. Unlike earlier methods, we allow the constant of proportionality to vary as a function of tree and herbaceous vegetation cover, and the mean size of monthly cumulative fire-pixel clusters. In areas undergoing active deforestation, we implemented a subsequent correction based on tree cover information and a simple measure of fire persistence. Regions showing good agreement between predicted and observed burned area included Boreal Asia, Central Asia, Europe, and Temperate North America, where the estimates produced by the regression trees were relatively accurate and precise. Poorest agreement was found for southern-hemisphere South America, where predicted values of burned area are both inaccurate and imprecise; this is most likely a consequence of multiple factors that include extremely persistent cloud cover, and lower quality of the 500-m burned area maps used for calibration. Application of our approach to the nine remaining regions yielded comparatively accurate, but less precise, estimates of monthly burned area. We applied the regional regression trees to the entire archive of Terra MODIS fire data to produce a monthly global burned area data set spanning late 2000 through mid-2005. Annual totals derived from this approach showed good agreement with independent annual estimates available for nine Canadian provinces, the United States, and Russia. With our data set we estimate the global annual burned area for the years 2001-2004 to vary between 2.97 million and 3.74 million km², with the maximum occurring in 2001. These coarse-resolution burned area estimates may serve as a useful interim product until long-term burned area data sets from multiple sensors and retrieval approaches become available

Time-dependent inversion estimates of global biomass-burning CO emissions using Measurement of Pollution in the Troposphere (MOPITT) measurements

We present an inverse-modeling analysis of CO emissions using column CO retrievals from the Measurement of Pollution in the Troposphere (MOPITT) instrument and a global chemical transport model (GEOS-CHEM). We first focus on the information content of MOPITT CO column retrievals in terms of constraining CO emissions associated with biomass burning and fossil fuel/biofuel use. Our analysis shows that seasonal variation of biomass-burning CO emissions in Africa, South America, and Southeast Asia can be characterized using monthly mean MOPITT CO columns. For the fossil fuel/biofuel source category the derived monthly mean emission estimates are noisy even when the error statistics are accurately known, precluding a characterization of seasonal variations of regional CO emissions for this source category. The derived estimate of CO emissions from biomass burning in southern Africa during the June-July 2000 period is significantly higher than the prior estimate (prior, 34 Tg; posterior, 13 Tg). We also estimate that emissions are higher relative to the prior estimate in northern Africa during December 2000 to January 2001 and lower relative to the prior estimate in Central America and Oceania/Indonesia during April-May and September-October 2000, respectively. While these adjustments provide better agreement of the model with MOPITT CO column fields and with independent measurements of surface CO from National Oceanic and Atmospheric Administration Climate Monitoring and Diagnostics Laboratory at background sites in the Northern Hemisphere, some systematic differences between modeled and measured CO fields persist, including model overestimation of background surface CO in the Southern Hemisphere. Characterizing and accounting for underlying biases in the measurement model system are needed to improve the robustness of the top-down estimate. Copyright 2006 by the American Geophysical Union

Estimates of fire emissions from an active deforestation region in the southern Amazon based on satellite data and biogeochemical modelling

Tropical deforestation contributes to the build-up of atmospheric carbon dioxide in the atmosphere. Within the deforestation process, fire is frequently used to eliminate biomass in preparation for agricultural use. Quantifying these deforestation-induced fire emissions represents a challenge, and current estimates are only available at coarse spatial resolution with large uncertainty. Here we developed a biogeochemical model using remote sensing observations of plant productivity, fire activity, and deforestation rates to estimate emissions for the Brazilian state of Mato Grosso during 2001–2005. Our model of DEforestation CArbon Fluxes (DECAF) runs at 250-m spatial resolution with a monthly time step to capture spatial and temporal heterogeneity in fire dynamics in our study area within the ''arc of deforestation'', the southern and eastern fringe of the Amazon tropical forest where agricultural expansion is most concentrated. Fire emissions estimates from our modelling framework were on average 90 Tg C year<sup>−1</sup>, mostly stemming from fires associated with deforestation (74%) with smaller contributions from fires from conversions of Cerrado or pastures to cropland (19%) and pasture fires (7%). In terms of carbon dynamics, about 80% of the aboveground living biomass and litter was combusted when forests were converted to pasture, and 89% when converted to cropland because of the highly mechanized nature of the deforestation process in Mato Grosso. The trajectory of land use change from forest to other land uses often takes more than one year, and part of the biomass that was not burned in the dry season following deforestation burned in consecutive years. This led to a partial decoupling of annual deforestation rates and fire emissions, and lowered interannual variability in fire emissions. Interannual variability in the region was somewhat dampened as well because annual emissions from fires following deforestation and from maintenance fires did not covary, although the effect was small due to the minor contribution of maintenance fires. Our results demonstrate how the DECAF model can be used to model deforestation fire emissions at relatively high spatial and temporal resolutions. Detailed model output is suitable for policy applications concerned with annual emissions estimates distributed among post-clearing land uses and science applications in combination with atmospheric emissions modelling to provide constrained global deforestation fire emissions estimates. DECAF currently estimates emissions from fire; future efforts can incorporate other aspects of net carbon emissions from deforestation including soil respiration and regrowth

Global emissions of non-methane hydrocarbons deduced from SCIAMACHY formaldehyde columns through 2003-2006

Formaldehyde columns retrieved from the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography/Chemistry (SCIAMACHY) instrument onboard ENVISAT satellite through 2003 to 2006 are used as top-down constraints to derive updated global biogenic and biomass burning flux estimates for the non-methane volatile organic compounds (NMVOCs) precursors of formaldehyde. Our interest is centered over regions experiencing strong emissions, and hence exhibiting a high signal-to-noise ratio and lower measurement uncertainties. The formaldehyde dataset used in this study has been recently made available to the community and complements the long record of formaldehyde measurements from the Global Ozone Monitoring Experiment (GOME). We use the IMAGESv2 global chemistry-transport model driven by the Global Fire Emissions Database (GFED) version 1 or 2 for biomass burning, and from the newly developed MEGAN-ECMWF isoprene emission database. The adjoint of the model is implemented in a grid-based framework within which emission fluxes are derived at the model resolution, together with a differentiation of the sources in a grid cell. Two inversion studies are conducted using either the GFEDv1 or GFEDv2 as a priori for the pyrogenic fluxes. Although on the global scale the inferred emissions from the two categories exhibit only weak deviations from the corresponding a priori estimates, the regional updates often present large departures from their a priori values. The posterior isoprene emissions over North America, amounting to about 34 Tg C/yr, are estimated to be on average by 25% lower than the a priori over 2003–2006, whereas a strong increase (55%) is deduced over the south African continent, the optimized emission being estimated at 57 Tg C/yr. Over Indonesia the biogenic emissions appear to be overestimated by 20–30%, whereas over Indochina and the Amazon basin during the wet season the a priori inventory captures both the seasonality and the magnitude of the observed columns. Although neither biomass burning inventory seems to be consistent with the data over all regions, pyrogenic estimates inferred from the two inversions are reasonably similar, despite their a priori deviations. A number of sensitivity experiments are conducted in order to assess the impact of uncertainties related to the inversion setup and the chemical mechanism. Whereas changes in the background error covariance matrix have only a limited impact on the posterior fluxes, the use of an alternative isoprene mechanism characterized by lower HCHO yields (the GEOS-Chem mechanism) increases the posterior isoprene source estimate by 11% over northern America, and by up to 40% in tropical regions

Evaluating the performance of pyrogenic and biogenic emission inventories against one decade of space-based formaldehyde columns

A new one-decade (1997–2006) dataset of formaldehyde (HCHO) columns retrieved from GOME and SCIAMACHY is compared with HCHO columns simulated by an updated version of the IMAGES global chemical transport model. This model version includes an optimized chemical scheme with respect to HCHO production, where the short-term and final HCHO yields from pyrogenically emitted non-methane volatile organic compounds (NMVOCs) are estimated from the Master Chemical Mechanism (MCM) and an explicit speciation profile of pyrogenic emissions. The model is driven by the Global Fire Emissions Database (GFED) version 1 or 2 for biomass burning, whereas biogenic emissions are provided either by the Global Emissions Inventory Activity (GEIA), or by a newly developed inventory based on the Model of Emissions of Gases and Aerosols from Nature (MEGAN) algorithms driven by meteorological fields from the European Centre for Medium-Range Weather Forecasts (ECMWF). The comparisons focus on tropical ecosystems, North America and China, which experience strong biogenic and biomass burning NMVOC emissions reflected in the enhanced measured HCHO columns. These comparisons aim at testing the ability of the model to reproduce the observed features of the HCHO distribution on the global scale and at providing a first assessment of the performance of the current emission inventories. The high correlation coefficients (<i>r</i>>0.7) between the observed and simulated columns over most regions indicate a good consistency between the model, the implemented inventories and the HCHO dataset. The use of the MEGAN-ECMWF inventory improves the model/data agreement in almost all regions, but biases persist over parts of Africa and Australia. Although neither GFED version is consistent with the data over all regions, a better agreement is achieved over Indonesia and Southern Africa when GFEDv2 is used, but GFEDv1 succeeds better in getting the correct seasonal patterns and intensities of the fire episodes over the Amazon basin, as reflected in the significantly higher correlations calculated in this region. Although the uncertainties in the HCHO retrievals, especially over fire scenes, can be quite large, this study provides a first assessment about whether the improved methodologies and input data implemented in GFEDv2 and MEGAN-ECMWF lead to better results in the comparisons of modelled with observed HCHO column measurements

Scanning Imaging Absorption Spectrometer for Atmospheric Chartography carbon monoxide total columns: Statistical evaluation and comparison with chemistry transport model results

This paper presents a detailed statistical analysis of one year (September 2003 to August 2004) of global Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) carbon monoxide (CO) total column retrievals from the Iterative Maximum Likelihood Method (IMLM) algorithm, version 6.3. SCIAMACHY provides the first solar reflectance measurements of CO and is uniquely sensitive down to the boundary layer. SCIAMACHY measurements and chemistry transport model (CTM) results are compared and jointly evaluated. Significant improvements in agreement occur, especially close to biomass burning emission regions, when the new Global Fire Emissions Database version 2 (GFEDv2) is used with the CTM. Globally, the seasonal variation of the model is very similar to that of the SCIAMACHY measurements. For certain locations, significant differences were found, which are likely related to modeling errors due to CO emission uncertainties. Statistical analysis shows that differences between single SCIAMACHY CO total column measurements and corresponding model results are primarily explained by random instrument noise errors. This strongly suggests that the random instrument noise errors are a good diagnostic for the precision of the measurements. The analysis also indicates that noise in single SCIAMACHY CO measurements is generally greater than actual variations in total columns. It is thus required to average SCIAMACHY data over larger temporal and spatial scales to obtain valuable information. Analyses of monthly averaged SCIAMACHY measurements over 3° × 2° geographical regions indicates that they are of sufficient accuracy to reveal valuable information about spatial and temporal variations in CO columns and provide an important tool for model validation. A large spatial and temporal variation in instrument noise errors exists which shows a close correspondence with the spatial distribution of surface albedo and cloud cover. This large spatial variability is important for the use of monthly and annual mean SCIAMACHY CO total column measurements. The smallest instrument noise errors of monthly mean 3° × 2° SCIAMACHY CO total columns measurements are 0.01 × 1018 molecules/cm2 for high surface albedo areas over the Sahara. Errors in SCIAMACHY CO total column retrievals due to errors other than instrument noise, like cloud cover, calibration, retrieval uncertainties and averaging kernels are estimated to be about 0.05–0.1 × 1018 molecules/cm2 in total. The bias found between model and observations is around 0.05–0.1 1018 molecules/cm2 (or about 5%) which also includes model errors. This thus provides a best estimate of the currently achievable measurement accuracy for SCIAMACHY CO monthly mean averages

Towards multi-tracer data-assimilation: biomass burning and carbon isotope exchange in SiBCASA

We present an enhanced version of the SiBCASA photosynthetic/biogeochemical model for a future integration with a multi-tracer data-assimilation system. We extended the model with (a) biomass burning emissions from the SiBCASA carbon pools using remotely sensed burned area from Global Fire Emissions Database (GFED) version 3.1, (b) a new set of 13C pools that cycle consistently through the biosphere, and (c), a modified isotopic discrimination scheme to estimate variations in 13C exchange as a~response to stomatal conductance. Previous studies suggest that the observed variations of atmospheric 13C/12C are driven by processes specifically in the terrestrial biosphere rather than in the oceans. Therefore, we quantify in this study the terrestrial exchange of CO2 and 13CO2 as a function of environmental changes in humidity and biomass burning. Based on an assessment of observed respiration signatures we conclude that SiBCASA does well in simulating global to regional plant discrimination. The global mean discrimination value is 15.2‰, and ranges between 4 and 20‰ depending on the regional plant phenology. The biomass burning emissions (annually and seasonally) compare favorably to other published values. However, the observed short-term changes in discrimination and the respiration 13C signature are more difficult to capture. We see a too weak drought response in SiBCASA and too slow return of anomalies in respiration. We demonstrate possible ways to improve this, and discuss the implications for our current capacity to interpret atmospheric 13C observation

Optimizing global CO emission estimates using a four-dimensional variational data assimilation system and surface network observations

We apply a four-dimensional variational (4D-VAR) data assimilation system to optimize carbon monoxide (CO) emissions for 2003 and 2004 and to reduce the uncertainty of emission estimates from individual sources using the chemistry transport model TM5. The system is designed to assimilate large (satellite) datasets, but in the current study only a limited amount of surface network observations from the National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL) Global Monitoring Division (GMD) is used to test the 4D-VAR system. By design, the system is capable to adjust the emissions in such a way that the posterior simulation reproduces background CO mixing ratios and large-scale pollution events at background stations. Uncertainty reduction up to 60 % in yearly emissions is observed over well-constrained regions and the inferred emissions compare well with recent studies for 2004. However, with the limited amount of data from the surface network, the system becomes data sparse resulting in a large solution space. Sensitivity studies have shown that model uncertainties (e.g., vertical distribution of biomass burning emissions and the OH field) and the prior inventories used, influence the inferred emission estimates. Also, since the observations only constrain total CO emissions, the 4D-VAR system has difficulties in separating anthropogenic and biogenic sources in particular. The inferred emissions are validated with NOAA aircraft data over North America and the agreement is significantly improved from the prior to posterior simulation. Validation with the Measurements Of Pollution In The Troposphere (MOPITT) instrument version 4 (V4) shows a slight improved agreement over the well-constrained Northern Hemisphere and in the tropics (except for the African continent). However, the model simulation with posterior emissions underestimates MOPITT CO total columns on the remote Southern Hemisphere (SH) by about 10 %. This is caused by a reduction in SH CO sources mainly due to surface stations on the high southern latitudes