Nine years of global hydrocarbon emissions based on source inversion of OMI formaldehyde observations
Authors: Bauwens Maite , Stavrakou Trissevgeni , Müller Jean-François , De Smedt Isabelle , Van Roozendael Michel , van der Werf Guido R. , Wiedinmyer Christine , Kaiser Johannes W. , Sindelarova Katerina , Guenther Alex 
Published: August 2016

Abstract. As formaldehyde (HCHO) is a high-yield product in the oxidation of most volatile organic compounds (VOCs) emitted by fires, vegetation, and anthropogenic activities, satellite observations of HCHO are well-suited to inform us on the spatial and temporal variability of the underlying VOC sources. The long record of space-based HCHO column observations from the Ozone Monitoring Instrument (OMI) is used to infer emission flux estimates from pyrogenic and biogenic volatile organic compounds (VOCs) on the global scale over 2005–2013. This is realized through the method of source inverse modeling, which consists in the optimization of emissions in a chemistry-transport model (CTM) in order to minimize the discrepancy between the observed and modeled HCHO columns. The top–down fluxes are derived in the global CTM IMAGESv2 by an iterative minimization algorithm based on the full adjoint of IMAGESv2, starting from a priori emission estimates provided by the newly released GFED4s (Global Fire Emission Database, version 4s) inventory for fires, and by the MEGAN-MOHYCAN inventory for isoprene emissions. The top–down fluxes are compared to two independent inventories for fire (GFAS and FINNv1.5) and isoprene emissions (MEGAN-MACC and GUESS-ES). The inversion indicates a moderate decrease (ca. 20 %) in the average annual global fire and isoprene emissions, from 2028 Tg C in the a priori to 1653 Tg C for burned biomass, and from 343 to 272 Tg for isoprene fluxes. Those estimates are acknowledged to depend on the accuracy of formaldehyde data, as well as on the assumed fire emission factors and the oxidation mechanisms leading to HCHO production. Strongly decreased top–down fire fluxes (30–50 %) are inferred in the peak fire season in Africa and during years with strong a priori fluxes associated with forest fires in Amazonia (in 2005, 2007, and 2010), bushfires in Australia (in 2006 and 2011), and peat burning in Indonesia (in 2006 and 2009), whereas generally increased fluxes are suggested in Indochina and during the 2007 fires in southern Europe. Moreover, changes in fire seasonal patterns are suggested; e.g., the seasonal amplitude is reduced over southeast Asia. In Africa, the inversion indicates increased fluxes due to agricultural fires and decreased maxima when natural fires are dominant. The top–down fire emissions are much better correlated with MODIS fire counts than the a priori inventory in regions with small and agricultural fires, indicating that the OMI-based inversion is well-suited to assess the associated emissions. Regarding biogenic sources, significant reductions in isoprene fluxes are inferred in tropical ecosystems (30–40 %), suggesting overestimated basal emission rates in those areas in the bottom–up inventory, whereas strongly positive isoprene emission updates are derived over semiarid and desert areas, especially in southern Africa and Australia. This finding suggests that the parameterization of the soil moisture stress used in MEGAN greatly exaggerates the flux reduction due to drought in those regions. The isoprene emission trends over 2005–2013 are often enhanced after optimization, with positive top–down trends in Siberia (4.2 % year−1) and eastern Europe (3.9 % year−1), likely reflecting forest expansion and warming temperatures, and negative trends in Amazonia (−2.1 % year−1), south China (−1 % year−1), the United States (−3.7 % year−1), and western Europe (−3.3 % year−1), which are generally corroborated by independent studies, yet their interpretation warrants further investigation.