Search

An updated and expanded representation of organics in the chemistry general circulation model EMAC (ECHAM5/MESSy for Atmospheric Chemistry) has been evaluated. First, the comprehensive Mainz Organic Mechanism (MOM) in the submodel MECCA (Module Efficiently Calculating the Chemistry of the Atmosphere) was activated with explicit degradation of organic species up to five carbon atoms and a simplified mechanism for larger molecules. Second, the ORACLE submodel (version 1.0) now considers condensation on aerosols for all organics in the mechanism. Parameterizations for aerosol yields are used only for the lumped species that are not included in the explicit mechanism. The simultaneous usage of MOM and ORACLE allows an efficient estimation of not only the chemical degradation of the simulated volatile organic compounds but also the contribution of organics to the growth and fate of (organic) aerosol, with the complexity of the mechanism largely increased compared to EMAC simulations with more simplified chemistry. The model evaluation presented here reveals that the OH concentration is reproduced well globally, whereas significant biases for observed oxygenated organics are present. We also investigate the general properties of the aerosols and their composition, showing that the more sophisticated and process-oriented secondary aerosol formation does not degrade the good agreement of previous model configurations with observations at the surface, allowing further research in the field of gas–aerosol interactions.

We developed a high-resolution surface flux inversion system based on the global Eulerian–Lagrangian coupled tracer transport model composed of the National Institute for Environmental Studies (NIES) transport model (TM; collectively NIES-TM) and the FLEXible PARTicle dispersion model (FLEXPART). The inversion system is named NTFVAR (NIES-TM–FLEXPART-variational) as it applies a variational optimization to estimate surface fluxes. We tested the system by estimating optimized corrections to natural surface CO2 fluxes to achieve the best fit to atmospheric CO2 data collected by the global in situ network as a necessary step towards the capability of estimating anthropogenic CO2 emissions. We employed the Lagrangian particle dispersion model (LPDM) FLEXPART to calculate surface flux footprints of CO2 observations at a spatial resolution of 0.1∘×0.1∘. The LPDM is coupled with a global atmospheric tracer transport model (NIES-TM). Our inversion technique uses an adjoint of the coupled transport model in an iterative optimization procedure. The flux error covariance operator was implemented via implicit diffusion. Biweekly flux corrections to prior flux fields were estimated for the years 2010–2012 from in situ CO2 data included in the Observation Package (ObsPack) data set. High-resolution prior flux fields were prepared using the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) for fossil fuel combustion, the Global Fire Assimilation System (GFAS) for biomass burning, the Vegetation Integrative SImulator for Trace gases (VISIT) model for terrestrial biosphere exchange, and the Ocean Tracer Transport Model (OTTM) for oceanic exchange. The terrestrial biospheric flux field was constructed using a vegetation mosaic map and a separate simulation of CO2 fluxes at a daily time step by the VISIT model for each vegetation type. The prior flux uncertainty for the terrestrial biosphere was scaled proportionally to the monthly mean gross primary production (GPP) by the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD17 product. The inverse system calculates flux corrections to the prior fluxes in the form of a relatively smooth field multiplied by high-resolution patterns of the prior flux uncertainties for land and ocean, following the coastlines and fine-scale vegetation productivity gradients. The resulting flux estimates improved the fit to the observations taken at continuous observation sites, reproducing both the seasonal and short-term concentration variabilities including high CO2 concentration events associated with anthropogenic emissions. The use of a high-resolution atmospheric transport in global CO2 flux inversions has the advantage of better resolving the transported mixed signals from the anthropogenic and biospheric sources in densely populated continental regions. Thus, it has the potential to achieve better separation between fluxes from terrestrial ecosystems and strong localized sources, such as anthropogenic emissions and forest fires. Further improvements in the modelling system are needed as our posterior fit was better than that of the National Oceanic and Atmospheric Administration (NOAA)'s CarbonTracker for only a fraction of the monitoring sites, i.e. mostly at coastal and island locations where background and local flux signals are mixed.

Gas flares are a regionally and globally significant source of atmospheric pollutants. They can be detected by satellite remote sensing. We calculate the global flared gas volume and black carbon emissions in 2017 by applying (1) a previously developed hot spot detection and characterisation algorithm to all observations of the Sea and Land Surface Temperature Radiometer (SLSTR) instrument on board the Copernicus satellite Sentinel-3A and (2) newly developed filters for identifying gas flares and corrections for calculating both flared gas volumes (billion cubic metres, BCM) and black carbon (BC) emissions (g). The filter to discriminate gas flares from other hot spots uses the observed hot spot characteristics in terms of temperature and persistence. A regression function is used to correct for the variability of detection opportunities. A total of 6232 flaring sites are identified worldwide. The best estimates of the annual flared gas volume and the BC emissions are 129 BCM with a confidence interval of [35, 419 BCM] and 73 Gg with a confidence interval of [20, 239 Gg], respectively. Comparison of our activity (i.e. BCM) results with those of the Visible Infrared Imaging Radiometer Suite (VIIRS) Nightfire data set and SWIR-based calculations show general agreement but distinct differences in several details. The calculation of black carbon emissions using our gas flaring data set with a newly developed dynamic assignment of emission factors lie in the range of recently published black carbon inventories, albeit towards the lower end. The data presented here can therefore be used e.g. in atmospheric dispersion simulations. The advantage of using our algorithm with Sentinel-3 data lies in the previously demonstrated ability to detect and quantify small flares, the long-term data availability from the Copernicus programme, and the increased detection opportunity of global gas flare monitoring when used in conjunction with the VIIRS instruments. The flaring activity and related black carbon emissions are available as “GFlaringS3” on the Emissions of atmospheric Compounds and Compilation of Ancillary Data (ECCAD) website (https://doi.org/10.25326/19, Caseiro and Kaiser, 2019).

Black carbon (BC) aerosols influence the Earth's atmosphere and climate, but their microphysical properties, spatiotemporal distribution, and long-range transport are not well constrained. This study presents airborne observations of the transatlantic transport of BC-rich African biomass burning (BB) smoke into the Amazon Basin using a Single Particle Soot Photometer (SP2) as well as several complementary techniques. We base our results on observations of aerosols and trace gases off the Brazilian coast onboard the HALO (High Altitude and LOng range) research aircraft during the ACRIDICON-CHUVA campaign in September 2014. During flight AC19 over land and ocean at the northeastern coastline of the Amazon Basin, we observed a BC-rich layer at ∼3.5 km altitude with a vertical extension of ∼0.3 km. Backward trajectories suggest that fires in African grasslands, savannas, and shrublands were the main source of this pollution layer and that the observed BB smoke had undergone more than 10 d of atmospheric transport and aging over the South Atlantic before reaching the Amazon Basin. The aged smoke is characterized by a dominant accumulation mode, centered at about 130 nm, with a particle concentration of Nacc=850±330 cm−3. The rBC particles account for ∼15 % of the submicrometer aerosol mass and ∼40 % of the total aerosol number concentration. This corresponds to a mass concentration range from 0.5 to 2 µg m−3 (1st to 99th percentiles) and a number concentration range from 90 to 530 cm−3. Along with rBC, high cCO (150±30 ppb) and cO3 (56±9 ppb) mixing ratios support the biomass burning origin and pronounced photochemical aging of this layer. Upon reaching the Amazon Basin, it started to broaden and to subside, due to convective mixing and entrainment of the BB aerosol into the boundary layer. Satellite observations show that the transatlantic transport of pollution layers is a frequently occurring process, seasonally peaking in August/September. By analyzing the aircraft observations together with the long-term data from the Amazon Tall Tower Observatory (ATTO), we found that the transatlantic transport of African BB smoke layers has a strong impact on the northern and central Amazonian aerosol population during the BB-influenced season (July to December). In fact, the early BB season (July to September) in this part of the Amazon appears to be dominated by African smoke, whereas the later BB season (October to December) appears to be dominated by South American fires. This dichotomy is reflected in pronounced changes in aerosol optical properties such as the single scattering albedo (increasing from 0.85 in August to 0.90 in November) and the BC-to-CO enhancement ratio (decreasing from 11 to 6 ng m−3 ppb−1). Our results suggest that, despite the high fraction of BC particles, the African BB aerosol acts as efficient cloud condensation nuclei (CCN), with potentially important implications for aerosol–cloud interactions and the hydrological cycle in the Amazon.

Fires and the aerosols that they emit impact air quality, health, and climate, but the abundance and properties of carbonaceous aerosol (both black carbon and organic carbon) from biomass burning (BB) remain uncertain and poorly constrained. We aim to explore the uncertainties associated with fire emissions and their air quality and radiative impacts from underlying dry matter consumed and emissions factors. To investigate this, we compare model simulations from a global chemical transport model, GEOS-Chem, driven by a variety of fire emission inventories with surface and airborne observations of black carbon (BC) and organic aerosol (OA) concentrations and satellite-derived aerosol optical depth (AOD). We focus on two fire-detection-based and/or burned-area-based (FD-BA) inventories using burned area and active fire counts, respectively, i.e., the Global Fire Emissions Database version 4 (GFED4s) with small fires and the Fire INventory from NCAR version 1.5 (FINN1.5), and two fire radiative power (FRP)-based approaches, i.e., the Quick Fire Emission Dataset version 2.4 (QFED2.4) and the Global Fire Assimilation System version 1.2 (GFAS1.2). We show that, across the inventories, emissions of BB aerosol (BBA) differ by a factor of 4 to 7 over North America and that dry matter differences, not emissions factors, drive this spread. We find that simulations driven by QFED2.4 generally overestimate BC and, to a lesser extent, OA concentrations observations from two fire-influenced aircraft campaigns in North America (ARCTAS and DC3) and from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network, while simulations driven by FINN1.5 substantially underestimate concentrations. The GFED4s and GFAS1.2-driven simulations provide the best agreement with OA and BC mass concentrations at the surface (IMPROVE), BC observed aloft (DC3 and ARCTAS), and AOD observed by MODIS over North America. We also show that a sensitivity simulation including an enhanced source of secondary organic aerosol (SOA) from fires, based on the NOAA Fire Lab 2016 experiments, produces substantial additional OA; however, the spread in the primary emissions estimates implies that this magnitude of SOA can be neither confirmed nor ruled out when comparing the simulations against the observations explored here. Given the substantial uncertainty in fire emissions, as represented by these four emission inventories, we find a sizeable range in 2012 annual BBA PM2.5 population-weighted exposure over Canada and the contiguous US (0.5 to 1.6 µg m−3). We also show that the range in the estimated global direct radiative effect of carbonaceous aerosol from fires (−0.11 to −0.048 W m−2) is large and comparable to the direct radiative forcing of OA (−0.09 W m−2) estimated in the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC). Our analysis suggests that fire emissions uncertainty challenges our ability to accurately characterize the impact of smoke on air quality and climate.

Fire emissions are a critical component of carbon and nutrient cycles and strongly affect climate and air quality. Dynamic global vegetation models (DGVMs) with interactive fire modeling provide important estimates for long-term and large-scale changes in fire emissions. Here we present the first multi-model estimates of global gridded historical fire emissions for 1700–2012, including carbon and 33 species of trace gases and aerosols. The dataset is based on simulations of nine DGVMs with different state-of-the-art global fire models that participated in the Fire Modeling Intercomparison Project (FireMIP), using the same and standardized protocols and forcing data, and the most up-to-date fire emission factor table based on field and laboratory studies in various land cover types. We evaluate the simulations of present-day fire emissions by comparing them with satellite-based products. The evaluation results show that most DGVMs simulate present-day global fire emission totals within the range of satellite-based products. They can capture the high emissions over the tropical savannas and low emissions over the arid and sparsely vegetated regions, and the main features of seasonality. However, most models fail to simulate the interannual variability, partly due to a lack of modeling peat fires and tropical deforestation fires. Before the 1850s, all models show only a weak trend in global fire emissions, which is consistent with the multi-source merged historical reconstructions used as input data for CMIP6. On the other hand, the trends are quite different among DGVMs for the 20th century, with some models showing an increase and others a decrease in fire emissions, mainly as a result of the discrepancy in their simulated responses to human population density change and land use and land cover change (LULCC). Our study provides an important dataset for further development of regional and global multi-source merged historical reconstructions, analyses of the historical changes in fire emissions and their uncertainties, and quantification of the role of fire emissions in the Earth system. It also highlights the importance of accurately modeling the responses of fire emissions to LULCC and population density change in reducing uncertainties in historical reconstructions of fire emissions and providing more reliable future projections.

The Amazon rain forest experiences the combined pressures from human-made deforestation and progressing climate change, causing severe and potentially disruptive perturbations of the ecosystem's integrity and stability. To intensify research on critical aspects of Amazonian biosphere–atmosphere exchange, the Amazon Tall Tower Observatory (ATTO) has been established in the central Amazon Basin. Here we present a multi-year analysis of backward trajectories to derive an effective footprint region of the observatory, which spans large parts of the particularly vulnerable eastern basin. Further, we characterize geospatial properties of the footprint regions, such as climatic conditions, distribution of ecoregions, land cover categories, deforestation dynamics, agricultural expansion, fire regimes, infrastructural development, protected areas, and future deforestation scenarios. This study is meant to be a resource and reference work, helping to embed the ATTO observations into the larger context of human-caused transformations of Amazonia. We conclude that the chances to observe an unperturbed rain forest–atmosphere exchange at the ATTO site will likely decrease in the future, whereas the atmospheric signals from human-made and climate-change-related forest perturbations will increase in frequency and intensity.

Lagrangian particle dispersion models (LPDMs) in backward mode are widely used to quantify the impact of transboundary pollution on downwind sites. Most LPDM applications count particles with a technique that introduces a so-called footprint layer (FL) with constant height, in which passing air tracer particles are assumed to be affected by surface emissions. The mixing layer dynamics are represented by the underlying meteorological model. This particle counting technique implicitly assumes that the atmosphere is well mixed in the FL. We have performed backward trajectory simulations with the FLEXPART model starting at Cyprus to calculate the sensitivity to emissions of upwind pollution sources. The emission sensitivity is used to quantify source contributions at the receptor and support the interpretation of ground measurements carried out during the CYPHEX campaign in July 2014. Here we analyse the effects of different constant and dynamic FL height assumptions. The results show that calculations with FL heights of 100 and 300 m yield similar but still discernible results. Comparison of calculations with FL heights constant at 300 m and dynamically following the planetary boundary layer (PBL) height exhibits systematic differences, with daytime and night-time sensitivity differences compensating for each other. The differences at daytime when a well-mixed PBL can be assumed indicate that residual inaccuracies in the representation of the mixing layer dynamics in the trajectories may introduce errors in the impact assessment on downwind sites. Emissions from vegetation fires are mixed up by pyrogenic convection which is not represented in FLEXPART. Neglecting this convection may lead to severe over- or underestimations of the downwind smoke concentrations. Introducing an extreme fire source from a different year in our study period and using fire-observation-based plume heights as reference, we find an overestimation of more than 60  % by the constant FL height assumptions used for surface emissions. Assuming a FL that follows the PBL may reproduce the peak of the smoke plume passing through but erroneously elevates the background for shallow stable PBL heights. It might thus be a reasonable assumption for open biomass burning emissions wherever observation-based injection heights are not available.

Fires have influenced atmospheric composition and climate since the rise of vascular plants, and satellite data have shown the overall global extent of fires. Our knowledge of historic fire emissions has progressively improved over the past decades due mostly to the development of new proxies and the improvement of fire models. Currently, there is a suite of proxies including sedimentary charcoal records, measurements of fire-emitted trace gases and black carbon stored in ice and firn, and visibility observations. These proxies provide opportunities to extrapolate emission estimates back in time based on satellite data starting in 1997, but each proxy has strengths and weaknesses regarding, for example, the spatial and temporal extents over which they are representative. We developed a new historic biomass burning emissions dataset starting in 1750 that merges the satellite record with several existing proxies and uses the average of six models from the Fire Model Intercomparison Project (FireMIP) protocol to estimate emissions when the available proxies had limited coverage. According to our approach, global biomass burning emissions were relatively constant, with 10-year averages varying between 1.8 and 2.3 Pg C yr−1. Carbon emissions increased only slightly over the full time period and peaked during the 1990s after which they decreased gradually. There is substantial uncertainty in these estimates, and patterns varied depending on choices regarding data representation, especially on regional scales. The observed pattern in fire carbon emissions is for a large part driven by African fires, which accounted for 58 % of global fire carbon emissions. African fire emissions declined since about 1950 due to conversion of savanna to cropland, and this decrease is partially compensated for by increasing emissions in deforestation zones of South America and Asia. These global fire emission estimates are mostly suited for global analyses and will be used in the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations.

The Global Fire Assimilation System (GFAS) assimilates fire radiative power (FRP) observations from satellite-based sensors to produce daily estimates of biomass burning emissions. It has been extended to include information about injection heights derived from fire observations and meteorological information from the operational weather forecasts of ECMWF. Injection heights are provided by two distinct methods: the Integrated Monitoring and Modelling System for wildland fires (IS4FIRES) parameterisation and the one-dimensional plume rise model (PRM). A global database of daily biomass burning emissions and injection heights at 0.1° resolution has been produced for 2003–2015 and is continuously extended in near-real time with the operational GFAS service of the Copernicus Atmospheric Monitoring Service (CAMS). In this study, the two injection height data sets were compared with the new MPHP2 (MISR Plume Height Project 2) satellite-based plume height retrievals. The IS4FIRES parameterisation showed a better overall agreement than the observations, while the PRM was better at capturing the variability of injection heights. The performance of both parameterisations is also dependent on the type of vegetation. Furthermore, the use of biomass burning emission heights from GFAS in atmospheric composition forecasts was assessed in two case studies: the South AMerican Biomass Burning Analysis (SAMBBA) campaign which took place in September 2012 in Brazil, and a series of large fire events in the western USA in August 2013. For these case studies, forecasts of biomass burning aerosol species by the Composition Integrated Forecasting System (C-IFS) of CAMS were found to better reproduce the observed vertical distribution when using PRM injection heights from GFAS compared to aerosols emissions being prescribed at the surface. The globally available GFAS injection heights introduced and evaluated in this study provide a comprehensive data set for future fire and atmospheric composition modelling studies.

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.

Landscape fires occur on a large scale in (sub)tropical savannas and grasslands, affecting ecosystem dynamics, regional air quality and concentrations of atmospheric trace gasses. Fuel consumption per unit of area burned is an important but poorly constrained parameter in fire emission modelling. We combined satellite-derived burned area with fire radiative power (FRP) data to derive fuel consumption estimates for land cover types with low tree cover in South America, Sub-Saharan Africa, and Australia. We developed a new approach to estimate fuel consumption, based on FRP data from the polar-orbiting Moderate Resolution Imaging Spectroradiometer (MODIS) and the geostationary Spinning Enhanced Visible and Infrared Imager (SEVIRI) in combination with MODIS burned-area estimates. The fuel consumption estimates based on the geostationary and polar-orbiting instruments showed good agreement in terms of spatial patterns. We used field measurements of fuel consumption to constrain our results, but the large variation in fuel consumption in both space and time complicated this comparison and absolute fuel consumption estimates remained more uncertain. Spatial patterns in fuel consumption could be partly explained by vegetation productivity and fire return periods. In South America, most fires occurred in savannas with relatively long fire return periods, resulting in comparatively high fuel consumption as opposed to the more frequently burning savannas in Sub-Saharan Africa. Strikingly, we found the infrequently burning interior of Australia to have higher fuel consumption than the more productive but frequently burning savannas in northern Australia. Vegetation type also played an important role in explaining the distribution of fuel consumption, by affecting both fuel build-up rates and fire return periods. Hummock grasslands, which were responsible for a large share of Australian biomass burning, showed larger fuel build-up rates than equally productive grasslands in Africa, although this effect might have been partially driven by the presence of grazers in Africa or differences in landscape management. Finally, land management in the form of deforestation and agriculture also considerably affected fuel consumption regionally. We conclude that combining FRP and burned-area estimates, calibrated against field measurements, is a promising approach in deriving quantitative estimates of fuel consumption. Satellite-derived fuel consumption estimates may both challenge our current understanding of spatiotemporal fuel consumption dynamics and serve as reference datasets to improve biogeochemical modelling approaches. Future field studies especially designed to validate satellite-based products, or airborne remote sensing, may further improve confidence in the absolute fuel consumption estimates which are quickly becoming the weakest link in fire emission estimates.

Fires associated with land use and land cover changes release large amounts of aerosols and trace gases into the atmosphere. Although several inventories of biomass burning emissions cover Brazil, there are still considerable uncertainties and differences among them. While most fire emission inventories utilize the parameters of burned area, vegetation fuel load, emission factors, and other parameters to estimate the biomass burned and its associated emissions, several more recent inventories apply an alternative method based on fire radiative power (FRP) observations to estimate the amount of biomass burned and the corresponding emissions of trace gases and aerosols. The Brazilian Biomass Burning Emission Model (3BEM) and the Fire Inventory from NCAR (FINN) are examples of the first, while the Brazilian Biomass Burning Emission Model with FRP assimilation (3BEM_FRP) and the Global Fire Assimilation System (GFAS) are examples of the latter. These four biomass burning emission inventories were used during the South American Biomass Burning Analysis (SAMBBA) field campaign. This paper analyzes and inter-compared them, focusing on eight regions in Brazil and the time period of 1 September–31 October 2012. Aerosol optical thickness (AOT550 nm) derived from measurements made by the Moderate Resolution Imaging Spectroradiometer (MODIS) operating on board the Terra and Aqua satellites is also applied to assess the inventories' consistency. The daily area-averaged pyrogenic carbon monoxide (CO) emission estimates exhibit significant linear correlations (r, p > 0.05 level, Student t test) between 3BEM and FINN and between 3BEM_ FRP and GFAS, with values of 0.86 and 0.85, respectively. These results indicate that emission estimates in this region derived via similar methods tend to agree with one other. However, they differ more from the estimates derived via the alternative approach. The evaluation of MODIS AOT550 nm indicates that model simulation driven by 3BEM and FINN typically underestimate the smoke particle loading in the eastern region of Amazon forest, while 3BEM_FRP estimations to the area tend to overestimate fire emissions. The daily regional CO emission fluxes from 3BEM and FINN have linear correlation coefficients of 0.75–0.92, with typically 20–30 % higher emission fluxes in FINN. The daily regional CO emission fluxes from 3BEM_FRP and GFAS show linear correlation coefficients between 0.82 and 0.90, with a particularly strong correlation near the arc of deforestation in the Amazon rainforest. In this region, GFAS has a tendency to present higher CO emissions than 3BEM_FRP, while 3BEM_FRP yields more emissions in the area of soybean expansion east of the Amazon forest. Atmospheric aerosol optical thickness is simulated by using the emission inventories with two operational atmospheric chemistry transport models: the IFS from Monitoring Atmospheric Composition and Climate (MACC) and the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modelling System (CCATT-BRAMS). Evaluation against MODIS observations shows a good representation of the general patterns of the AOT550 nm time series. However, the aerosol emissions from fires with particularly high biomass consumption still lead to an underestimation of the atmospheric aerosol load in both models.

Characterising the dynamics of landscape-scale wildfires at very high temporal resolutions is best achieved using observations from Earth Observation (EO) sensors mounted onboard geostationary satellites. As a result, a number of operational active fire products have been developed from the data of such sensors. An example of which are the Fire Radiative Power (FRP) products, the FRP-PIXEL and FRP-GRID products, generated by the Land Surface Analysis Satellite Applications Facility (LSA SAF) from imagery collected by the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard the Meteosat Second Generation (MSG) series of geostationary EO satellites. The processing chain developed to deliver these FRP products detects SEVIRI pixels containing actively burning fires and characterises their FRP output across four geographic regions covering Europe, part of South America and Northern and Southern Africa. The FRP-PIXEL product contains the highest spatial and temporal resolution FRP data set, whilst the FRP-GRID product contains a spatio-temporal summary that includes bias adjustments for cloud cover and the non-detection of low FRP fire pixels. Here we evaluate these two products against active fire data collected by the Moderate Resolution Imaging Spectroradiometer (MODIS) and compare the results to those for three alternative active fire products derived from SEVIRI imagery. The FRP-PIXEL product is shown to detect a substantially greater number of active fire pixels than do alternative SEVIRI-based products, and comparison to MODIS on a per-fire basis indicates a strong agreement and low bias in terms of FRP values. However, low FRP fire pixels remain undetected by SEVIRI, with errors of active fire pixel detection commission and omission compared to MODIS ranging between 9–13 % and 65–77 % respectively in Africa. Higher errors of omission result in greater underestimation of regional FRP totals relative to those derived from simultaneously collected MODIS data, ranging from 35 % over the Northern Africa region to 89 % over the European region. High errors of active fire omission and FRP underestimation are found over Europe and South America and result from SEVIRI's larger pixel area over these regions. An advantage of using FRP for characterising wildfire emissions is the ability to do so very frequently and in near-real time (NRT). To illustrate the potential of this approach, wildfire fuel consumption rates derived from the SEVIRI FRP-PIXEL product are used to characterise smoke emissions of the 2007 "mega-fire" event focused on Peloponnese (Greece) and used within the European Centre for Medium-Range Weather Forecasting (ECMWF) Integrated Forecasting System (IFS) as a demonstration of what can be achieved when using geostationary active fire data within the Copernicus Atmosphere Monitoring Service (CAMS). Qualitative comparison of the modelled smoke plumes with MODIS optical imagery illustrates that the model captures the temporal and spatial dynamics of the plume very well, and that high temporal resolution emissions estimates such as those available from a geostationary orbit are important for capturing the sub-daily variability in smoke plume parameters such as aerosol optical depth (AOD), which are increasingly less well resolved using daily or coarser temporal resolution emissions data sets. Quantitative comparison of modelled AOD with coincident MODIS and AERONET (Aerosol Robotic Network) AOD indicates that the former is overestimated by ~ 20–30 %, but captures the observed AOD dynamics with a high degree of fidelity. The case study highlights the potential of using geostationary FRP data to drive fire emissions estimates for use within atmospheric transport models such as those implemented in the Monitoring Atmospheric Composition and Climate (MACC) series of projects for the CAMS.

Detecting the optical properties of aerosols using passive satellite-borne measurements alone is a difficult task due to the broadband effect of aerosols on the measured spectra and the influences of surface and cloud reflection. We present another approach to determine aerosol type, namely by studying the relationship of aerosol optical depth (AOD) with trace gas abundance, aerosol absorption, and mean aerosol size. Our new Global Aerosol Classification Algorithm, GACA, examines relationships between aerosol properties (AOD and extinction Ångström exponent from the Moderate Resolution Imaging Spectroradiometer (MODIS), UV Aerosol Index from the second Global Ozone Monitoring Experiment, GOME-2) and trace gas column densities (NO2, HCHO, SO2 from GOME-2, and CO from MOPITT, the Measurements of Pollution in the Troposphere instrument) on a monthly mean basis. First, aerosol types are separated based on size (Ångström exponent) and absorption (UV Aerosol Index), then the dominating sources are identified based on mean trace gas columns and their correlation with AOD. In this way, global maps of dominant aerosol type and main source type are constructed for each season and compared with maps of aerosol composition from the global MACC (Monitoring Atmospheric Composition and Climate) model. Although GACA cannot correctly characterize transported or mixed aerosols, GACA and MACC show good agreement regarding the global seasonal cycle, particularly for urban/industrial aerosols. The seasonal cycles of both aerosol type and source are also studied in more detail for selected 5° × 5° regions. Again, good agreement between GACA and MACC is found for all regions, but some systematic differences become apparent: the variability of aerosol composition (yearly and/or seasonal) is often not well captured by MACC, the amount of mineral dust outside of the dust belt appears to be overestimated, and the abundance of secondary organic aerosols is underestimated in comparison with GACA. Whereas the presented study is of exploratory nature, we show that the developed algorithm is well suited to evaluate climate and atmospheric composition models by including aerosol type and source obtained from measurements into the comparison, instead of focusing on a single parameter, e.g., AOD. The approach could be adapted to constrain the mix of aerosol types during the process of a combined data assimilation of aerosol and trace gas observations.

This paper describes the pre-operational analysis and forecasting system developed during MACC (Monitoring Atmospheric Composition and Climate) and continued in the MACC-II (Monitoring Atmospheric Composition and Climate: Interim Implementation) European projects to provide air quality services for the European continent. This system is based on seven state-of-the art models developed and run in Europe (CHIMERE, EMEP, EURAD-IM, LOTOS-EUROS, MATCH, MOCAGE and SILAM). These models are used to calculate multi-model ensemble products. The paper gives an overall picture of its status at the end of MACC-II (summer 2014) and analyses the performance of the multi-model ensemble. The MACC-II system provides daily 96 h forecasts with hourly outputs of 10 chemical species/aerosols (O3, NO2, SO2, CO, PM10, PM2.5, NO, NH3, total NMVOCs (non-methane volatile organic compounds) and PAN+PAN precursors) over eight vertical levels from the surface to 5 km height. The hourly analysis at the surface is done a posteriori for the past day using a selection of representative air quality data from European monitoring stations. The performance of the system is assessed daily, weekly and every 3 months (seasonally) through statistical indicators calculated using the available representative air quality data from European monitoring stations. Results for a case study show the ability of the ensemble median to forecast regional ozone pollution events. The seasonal performances of the individual models and of the multi-model ensemble have been monitored since September 2009 for ozone, NO2 and PM10. The statistical indicators for ozone in summer 2014 show that the ensemble median gives on average the best performances compared to the seven models. There is very little degradation of the scores with the forecast day but there is a marked diurnal cycle, similarly to the individual models, that can be related partly to the prescribed diurnal variations of anthropogenic emissions in the models. During summer 2014, the diurnal ozone maximum is underestimated by the ensemble median by about 4 μg m−3 on average. Locally, during the studied ozone episodes, the maxima from the ensemble median are often lower than observations by 30–50 μg m−3. Overall, ozone scores are generally good with average values for the normalised indicators of 0.14 for the modified normalised mean bias and of 0.30 for the fractional gross error. Tests have also shown that the ensemble median is robust to reduction of ensemble size by one, that is, if predictions are unavailable from one model. Scores are also discussed for PM10 for winter 2013–1014. There is an underestimation of most models leading the ensemble median to a mean bias of −4.5 μg m−3. The ensemble median fractional gross error is larger for PM10 (~ 0.52) than for ozone and the correlation is lower (~ 0.35 for PM10 and ~ 0.54 for ozone). This is related to a larger spread of the seven model scores for PM10 than for ozone linked to different levels of complexity of aerosol representation in the individual models. In parallel, a scientific analysis of the results of the seven models and of the ensemble is also done over the Mediterranean area because of the specificity of its meteorology and emissions. The system is robust in terms of the production availability. Major efforts have been done in MACC-II towards the operationalisation of all its components. Foreseen developments and research for improving its performances are discussed in the conclusion.

The El Niño–Southern Oscillation (ENSO) not only affects meteorological fields but also has a large impact on atmospheric composition. Atmospheric composition fields from the Monitoring Atmospheric Composition and Climate (MACC) reanalysis are used to identify the ENSO signal in tropospheric ozone, carbon monoxide, nitrogen oxide and smoke aerosols, concentrating on the months October to December. During El Niño years, all of these fields have increased concentrations over maritime South East Asia in October. The MACC Composition Integrated Forecasting System (C-IFS) model is used to quantify the relative magnitude of dynamically induced and emission-driven changes in the atmospheric composition fields. While changes in tropospheric ozone are a combination of dynamically induced and emission-driven changes, the changes in carbon monoxide, nitrogen oxides and smoke aerosols are almost entirely emission-driven in the MACC model. The ozone changes continue into December, i.e. after the end of the Indonesian fire season while changes in the other fields are confined to the fire season.

Accurate near real time fire emissions estimates are required for air quality forecasts. To date, most approaches are based on satellite-derived estimates of fire radiative power (FRP), which can be converted to fire radiative energy (FRE) which is directly related to fire emissions. Uncertainties in these FRE estimates are often substantial. This is for a large part because the most often used low-Earth orbit satellite-based instruments such as the Moderate Resolution Imaging Spectroradiometer (MODIS) have a relatively poor sampling of the usually pronounced fire diurnal cycle. In this paper we explore the spatial variation of this fire diurnal cycle and its drivers using data from the geostationary Meteosat Spinning Enhanced Visible and Infrared Imager (SEVIRI). In addition, we sampled data from the SEVIRI instrument at MODIS detection opportunities to develop two approaches to estimate hourly FRE based on MODIS active fire detections. The first approach ignored the fire diurnal cycle, assuming persistent fire activity between two MODIS observations, while the second approach combined knowledge on the climatology of the fire diurnal cycle with active fire detections to estimate hourly FRE. The full SEVIRI time series, providing full coverage of the fire diurnal cycle, were used to evaluate the results. Our study period comprised of 3 years (2010–2012), and we focused on Africa and the Mediterranean basin to avoid the use of potentially lower quality SEVIRI data obtained at very far off-nadir view angles. We found that the fire diurnal cycle varies substantially over the study region, and depends on both fuel and weather conditions. For example, more "intense" fires characterized by a fire diurnal cycle with high peak fire activity, long duration over the day, and with nighttime fire activity are most common in areas of large fire size (i.e., large burned area per fire event). These areas are most prevalent in relatively arid regions. Ignoring the fire diurnal cycle generally resulted in an overestimation of FRE, while including information on the climatology of the fire diurnal cycle improved FRE estimates. The approach based on knowledge of the climatology of the fire diurnal cycle also improved distribution of FRE over the day, although only when aggregating model results to coarser spatial and/or temporal scale good correlation was found with the full SEVIRI hourly reference data set. We recommend the use of regionally varying fire diurnal cycle information within the Global Fire Assimilation System (GFAS) used in the Copernicus Atmosphere Monitoring Services, which will improve FRE estimates and may allow for further reconciliation of biomass burning emission estimates from different inventories.

Wildfires represent a major source for aerosols impacting atmospheric radiation, atmospheric chemistry and cloud micro-physical properties. Previous case studies indicated that the height of the aerosol–radiation interaction may crucially affect atmospheric radiation, but the sensitivity to emission heights has been examined with only a few models and is still uncertain. In this study we use the general circulation model ECHAM6 extended by the aerosol module HAM2 to investigate the impact of wildfire emission heights on atmospheric long-range transport, black carbon (BC) concentrations and atmospheric radiation. We simulate the wildfire aerosol release using either various versions of a semi-empirical plume height parametrization or prescribed standard emission heights in ECHAM6-HAM2. Extreme scenarios of near-surface or free-tropospheric-only injections provide lower and upper constraints on the emission height climate impact. We find relative changes in mean global atmospheric BC burden of up to 7.9±4.4 % caused by average changes in emission heights of 1.5–3.5 km. Regionally, changes in BC burden exceed 30–40 % in the major biomass burning regions. The model evaluation of aerosol optical thickness (AOT) against Moderate Resolution Imaging Spectroradiometer (MODIS), AErosol RObotic NETwork (AERONET) and Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) observations indicates that the implementation of a plume height parametrization slightly reduces the ECHAM6-HAM2 biases regionally, but on the global scale these improvements in model performance are small. For prescribed emission release at the surface, wildfire emissions entail a total sky top-of-atmosphere (TOA) radiative forcing (RF) of −0.16±0.06 W m−2. The application of a plume height parametrization which agrees reasonably well with observations introduces a slightly stronger negative TOA RF of −0.20±0.07 W m−2. The standard ECHAM6-HAM2 model in which 25 % of the wildfire emissions are injected into the free troposphere (FT) and 75 % into the planetary boundary layer (PBL), leads to a TOA RF of −0.24±0.06 W m−2. Overall, we conclude that simple plume height parametrizations provide sufficient representations of emission heights for global climate modeling. Significant improvements in aerosol wildfire modeling likely depend on better emission inventories and aerosol process modeling rather than on improved emission height parametrizations.

Among the atmospheric emission sources, wildfires are episodic events characterized by large spatial and temporal variability. Therefore, accurate information on gaseous and aerosol emissions from fires for specific regions and seasons is critical for air quality forecasts. The Spinning Enhanced Visible and Infrared Imager (SEVIRI) in geostationary orbit provides fire observations over Africa and the Mediterranean with a temporal resolution of 15 min. It thus resolves the complete fire life cycle and captures the fires' peak intensities, which is not possible in Moderate Resolution Imaging Spectroradiometer (MODIS) fire emission inventories like the Global Fire Assimilation System (GFAS). We evaluate two different operational fire radiative power (FRP) products derived from SEVIRI, by studying a large forest fire in Antalya, Turkey, in July–August 2008. The EUMETSAT Land Surface Analysis Satellite Applications Facility (LSA SAF) has higher FRP values during the fire episode than the Wildfire Automated Biomass Burning Algorithm (WF_ABBA). It is also in better agreement with the co-located, gridded MODIS FRP. Both products miss small fires that frequently occur in the region and are detected by MODIS. Emissions are derived from the FRP products. They are used along-side GFAS emissions in smoke plume simulations with the Weather Research and Forecasting (WRF) model and the Community Multiscale Air Quality (CMAQ) model. In comparisons with MODIS aerosol optical thickness (AOT) and Infrared Atmospheric Sounding Interferometer (IASI), CO and NH3 observations show that including the diurnal variability of fire emissions improves the spatial distribution and peak concentrations of the simulated smoke plumes associated with this large fire. They also show a large discrepancy between the currently available operational FRP products, with the LSA SAF being the most appropriate.

Large-scale fires occur frequently across Indonesia, particularly in the southern region of Kalimantan and eastern Sumatra. They have considerable impacts on carbon emissions, haze production, biodiversity, health, and economic activities. In this study, we demonstrate that severe fire and haze events in Indonesia can generally be predicted months in advance using predictions of seasonal rainfall from the ECMWF System 4 coupled ocean–atmosphere model. Based on analyses of long, up-to-date series observations on burnt area, rainfall, and tree cover, we demonstrate that fire activity is negatively correlated with rainfall and is positively associated with deforestation in Indonesia. There is a contrast between the southern region of Kalimantan (high fire activity, high tree cover loss, and strong non-linear correlation between observed rainfall and fire) and the central region of Kalimantan (low fire activity, low tree cover loss, and weak, non-linear correlation between observed rainfall and fire). The ECMWF seasonal forecast provides skilled forecasts of burnt and fire-affected area with several months lead time explaining at least 70% of the variance between rainfall and burnt and fire-affected area. Results are strongly influenced by El Niño years which show a consistent positive bias. Overall, our findings point to a high potential for using a more physical-based method for predicting fires with several months lead time in the tropics rather than one based on indexes only. We argue that seasonal precipitation forecasts should be central to Indonesia's evolving fire management policy.

The South American Biomass Burning Analysis (SAMBBA) field campaign took detailed in situ flight measurements of aerosol during the 2012 dry season to characterise biomass burning aerosol and improve understanding of its impacts on weather and climate. Developments have been made to the Weather Research and Forecast model with chemistry (WRF-Chem) model to improve the representation of biomass burning aerosol in the region, by coupling a sectional aerosol scheme to the plume-rise parameterisation. Brazilian Biomass Burning Emissions Model (3BEM) fire emissions are used, prepared using PREP-CHEM-SRC, and mapped to CBM-Z and MOSAIC species. Model results have been evaluated against remote sensing products, AERONET sites, and four case studies of flight measurements from the SAMBBA campaign. WRF-Chem predicted layers of elevated aerosol loadings (5–20 μg sm−3) of particulate organic matter at high altitude (6–8 km) over tropical forest regions, while flight measurements showed a sharp decrease above 2–4 km altitude. This difference was attributed to the plume-rise parameterisation overestimating injection height. The 3BEM emissions product was modified using estimates of active fire size and burned area for the 2012 fire season, which reduced the fire size. The enhancement factor for fire emissions was increased from 1.3 to 5 to retain reasonable aerosol optical depths (AODs). The smaller fire size lowered the injection height of the emissions, but WRF-Chem still showed elevated aerosol loadings between 4–5 km altitude. Over eastern cerrado (savannah-like) regions, both modelled and measured aerosol loadings decreased above approximately 4 km altitude. Compared with MODIS satellite data and AERONET sites, WRF-Chem represented AOD magnitude well (between 0.3–1.5) over western tropical forest fire regions in the first half of the campaign, but tended to over-predict them in the second half, when precipitation was more significant. Over eastern cerrado regions, WRF-Chem tended to under-predict AODs. Modelled aerosol loadings in the east were higher in the modified emission scenario. The primary organic matter to black carbon ratio was typically between 8–10 in WRF-Chem. This was lower than the western flight measurements (interquartile range of 11.6–15.7 in B734, 14.7–24.0 in B739), but similar to the eastern flight B742 (8.1–10.4). However, single scattering albedo was close to measured over the western flights (0.87–0.89 in model; 0.86–0.91 in flight B734, and 0.81–0.95 in flight B739 measurements) but too high over the eastern flight B742 (0.86–0.87 in model, 0.79–0.82 in measurements). This suggests that improvements are needed to both modelled aerosol composition and optical properties calculations in WRF-Chem.

Fires are important emitters of aerosol and trace gases and as such need to be taken into account in any atmospheric composition modelling enterprise. One method to estimate these emissions is to convert fire radiative power (FRP) analysis into dry matter burnt and emissions of smoke constituents using land-cover-dependent conversion factors. Inventories like the Global Fire Assimilation System (GFAS) follow this approach by calculating daily global smoke emissions from FRP observed by the MODIS instruments onboard the Terra and Aqua satellites. Observations with different overpass times systematically sample fires at different stages in the strong diurnal fire cycle. For some time periods, observations are available from only one instrument, which leads to a bias in the observed average FRP. We develop a method to correct this bias in daily FRP observations from any low Earth orbit (LEO) satellite, so that the budget of daily smoke emissions remains independent of the number of satellites from which FRP observations are taken into account. This ensures the possibility of running, for example, GFAS in case of failure of one of the MODIS instruments. It also enables the extension GFAS to 2000–2002 and the inclusion of FRP observations from upcoming satellite missions. The correction combines linear and non-linear regressions and uses an adaptive regionalization algorithm. It decreases the bias in daily average FRP from Terra and Aqua by more than 95%, and RMSE by 75% for Aqua and 55% for Terra. The correction algorithm is applied to Terra observations from 25 February 2000 to 31 December 2002, when Aqua observations were not available. The database of fire emissions GFASv1.0 is extended correspondingly.

A method to constrain carbon dioxide (CO2) emissions from open biomass burning by using satellite observations of co-emitted species and a chemistry-transport model (CTM) is proposed and applied to the case of wildfires in Siberia. CO2 emissions are assessed by means of an emission model assuming a direct relationship between the biomass burning rate (BBR) and the fire radiative power (FRP) derived from MODIS measurements. The key features of the method are (1) estimating the FRP-to-BBR conversion factors (α) for different vegetative land cover types by assimilating the satellite observations of co-emitted species into the CTM, (2) optimal combination of the estimates of α derived independently from satellite observations of different species (CO and aerosol in this study), and (3) estimation of the diurnal cycle of the fire emissions directly from the FRP measurements. Values of α for forest and grassland fires in Siberia and their uncertainties are estimated using the Infrared Atmospheric Sounding Interferometer (IASI) carbon monoxide (CO) retrievals and MODIS aerosol optical depth (AOD) measurements combined with outputs from the CHIMERE mesoscale chemistry-transport model. The constrained CO emissions are validated through comparison of the respective simulations with independent data of ground-based CO measurements at the ZOTTO site. Using our optimal regional-scale estimates of the conversion factors (which are found to be in agreement with earlier published estimates obtained from local measurements of experimental fires), the total CO2 emissions from wildfires in Siberia in 2012 are estimated to be in the range from 280 to 550 Tg C, with the optimal (maximum likelihood) value of 392 Tg C. Sensitivity test cases featuring different assumptions regarding the injection height and diurnal variations of emissions indicate that the derived estimates of the total CO2 emissions in Siberia are robust with respect to the modeling options (the different estimates vary within less than 15% of their magnitude). The CO2 emission estimates obtained for several years are compared with independent estimates provided by the GFED3.1 and GFASv1.0 global emission inventories. It is found that our "top-down" estimates for the total annual biomass burning CO2 emissions in the period from 2007 to 2011 in Siberia are by factors of 2.5 and 1.8 larger than the respective bottom-up estimates; these discrepancies cannot be fully explained by uncertainties in our estimates. There are also considerable differences in the spatial distribution of the different emission estimates; some of those differences have a systematic character and require further analysis.

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

Nitrogen oxides (NOx) play key roles in atmospheric chemistry, air pollution, and climate. While the largest fraction of these reactive gases is released by anthropogenic emission sources, a significant amount can be attributed to vegetation fires. In this study, NO2 from GOME-2 on board EUMETSAT's MetOp-A and OMI on board NASA's Aura as well as fire radiative power (FRP) from the measurements of MODIS on board NASA's Terra and Aqua satellites are used to derive fire emission rates (FERs) of NOx for different types of vegetation using a simple statistical approach. Monthly means of tropospheric NO2 vertical columns (TVC NO2) have been analyzed for their temporal correlation with the monthly means of FRP for five consecutive years from 2007 to 2011 on a horizontal 1° × 1° grid. The strongest correlation is found to be largely confined to tropical and subtropical regions, which account for more than 80% of yearly burned area, on average, globally. In these regions, the seasonal variation of fire intensity, expressed by the FRP data, is similar to the pattern of TVC NO2. As chemical models typically require values for the amount of NOx being released as a function of time, we have converted the retrieved TVC NO2 into production rates of NOx from fire (Pf) by assuming a constant lifetime of NOx. The comparison between Pf and NOx emissions from the Global Fire Emissions Database (GFEDv3.1) over 5 characteristic biomass burning regions in the tropics and subtropics shows good agreement. By separating the monthly means of Pf and FRP according to land cover type, FERs of NOx could be derived for different biomes. The estimated FERs for the dominating types of vegetation burned are lowest for open shrublands and savannas (0.28–1.03 g NOx s−1 MW−1) and highest for croplands and woody savannas (0.82–1.56 g NOx s−1 MW−1). This analysis demonstrates that the strong empirical relationship between TVC NO2 and FRP and the following simplified assumptions are a useful tool for the characterization of NOx emission rates from vegetation fires in the tropics and subtropics. Possible factors affecting the magnitude of the obtained values are discussed.

A new fast clear-sky model called McClear was developed to estimate the downwelling shortwave direct and global irradiances received at ground level under clear skies. It is a fully physical model replacing empirical relations or simpler models used before. It exploits the recent results on aerosol properties, and total column content in water vapour and ozone produced by the MACC project (Monitoring Atmosphere Composition and Climate). It accurately reproduces the irradiance computed by the libRadtran reference radiative transfer model with a computational speed approximately 105 times greater by adopting the abaci, or look-up table, approach combined with interpolation functions. It is therefore suited for geostationary satellite retrievals or numerical weather prediction schemes with many pixels or grid points, respectively. McClear irradiances were compared to 1 min measurements made in clear-sky conditions at several stations within the Baseline Surface Radiation Network in various climates. The bias for global irradiance comprises between −6 and 25 W m−2. The RMSE ranges from 20 W m−2 (3% of the mean observed irradiance) to 36 W m−2 (5%) and the correlation coefficient ranges between 0.95 and 0.99. The bias for the direct irradiance comprises between −48 and +33 W m−2. The root mean square error (RMSE) ranges from 33 W m−2 (5%) to 64 W m−2 (10%). The correlation coefficient ranges between 0.84 and 0.98. This work demonstrates the quality of the McClear model combined with MACC products, and indirectly the quality of the aerosol properties modelled by the MACC reanalysis.

The fires around Moscow in July and August 2010 emitted a large amount of pollutants to the atmosphere. Here we estimate the carbon monoxide (CO) source strength of the Moscow fires in July and August by using the TM5-4DVAR system in combination with CO column observations of the Infrared Atmospheric Sounding Interferometer (IASI). It is shown that the IASI observations provide a strong constraint on the total emissions needed in the model. Irrespective of the prior emissions used, the optimised CO fire emission estimates from mid-July to mid-August 2010 amount to approximately 24 Tg CO. This estimate depends only weakly (< 15%) on the assumed diurnal variations and injection height of the emissions. However, the estimated emissions might depend on unaccounted model uncertainties such as vertical transport. Our emission estimate of 22–27 Tg CO during roughly one month of intense burning is less than suggested by another recent study, but substantially larger than predicted by the bottom-up inventories. This latter discrepancy suggests that bottom-up emission estimates for extreme peat burning events require improvements.

An eight-year long reanalysis of atmospheric composition data covering the period 2003–2010 was constructed as part of the FP7-funded Monitoring Atmospheric Composition and Climate project by assimilating satellite data into a global model and data assimilation system. This reanalysis provides fields of chemically reactive gases, namely carbon monoxide, ozone, nitrogen oxides, and formaldehyde, as well as aerosols and greenhouse gases globally at a horizontal resolution of about 80 km for both the troposphere and the stratosphere. This paper describes the assimilation system for the reactive gases and presents validation results for the reactive gas analysis fields to document the data set and to give a first indication of its quality. Tropospheric CO values from the MACC reanalysis are on average 10–20% lower than routine observations from commercial aircrafts over airports through most of the troposphere, and have larger negative biases in the boundary layer at urban sites affected by air pollution, possibly due to an underestimation of CO or precursor emissions. Stratospheric ozone fields from the MACC reanalysis agree with ozonesondes and ACE-FTS data to within ±10% in most seasons and regions. In the troposphere the reanalysis shows biases of −5% to +10% with respect to ozonesondes and aircraft data in the extratropics, but has larger negative biases in the tropics. Area-averaged total column ozone agrees with ozone fields from a multi-sensor reanalysis data set to within a few percent. NO2 fields from the reanalysis show the right seasonality over polluted urban areas of the NH and over tropical biomass burning areas, but underestimate wintertime NO2 maxima over anthropogenic pollution regions and overestimate NO2 in northern and southern Africa during the tropical biomass burning seasons. Tropospheric HCHO is well simulated in the MACC reanalysis even though no satellite data are assimilated. It shows good agreement with independent SCIAMACHY retrievals over regions dominated by biogenic emissions with some anthropogenic input, such as the eastern US and China, and also over African regions influenced by biogenic sources and biomass burning.

The severe wildfires in western Russia during July–August 2010 coincided with a strong heat wave and led to large emissions of aerosols and trace gases such as carbon monoxide (CO), hydrocarbons and nitrogen oxides into the troposphere. This extreme event is used to evaluate the ability of the global MACC (Monitoring Atmospheric Composition and Climate) atmospheric composition forecasting system to provide analyses of large-scale pollution episodes and to test the respective influence of a priori emission information and data assimilation on the results. Daily 4-day hindcasts were conducted using assimilated aerosol optical depth (AOD), CO, nitrogen dioxide (NO2) and ozone (O3) data from a range of satellite instruments. Daily fire emissions were used from the Global Fire Assimilation System (GFAS) version 1.0, derived from satellite fire radiative power retrievals. The impact of accurate wildfire emissions is dominant on the composition in the boundary layer, whereas the assimilation system influences concentrations throughout the troposphere, reflecting the vertical sensitivity of the satellite instruments. The application of the daily fire emissions reduces the area-average mean bias by 63% (for CO), 60% (O3) and 75% (NO2) during the first 24 h with respect to independent satellite observations, compared to a reference simulation with a multi-annual mean climatology of biomass burning emissions. When initial tracer concentrations are further constrained by data assimilation, biases are reduced by 87, 67 and 90%. The forecast accuracy, quantified by the mean bias up to 96 h lead time, was best for all compounds when using both the GFAS emissions and assimilation. The model simulations suggest an indirect positive impact of O3 and CO assimilation on hindcasts of NO2 via changes in the oxidizing capacity. However, the quality of local hindcasts was strongly dependent on the assumptions made for forecasted fire emissions. This was well visible from a relatively poor forecast accuracy quantified by the root mean square error, as well as the temporal correlation with respect to ground-based CO total column data and AOD. This calls for a more advanced method to forecast fire emissions than the currently adopted persistency approach. The combined analysis of fire radiative power observations, multiple trace gas and aerosol satellite observations, as provided by the MACC system, results in a detailed quantitative description of the impact of major fires on atmospheric composition, and demonstrate the capabilities for the real-time analysis and forecasts of large-scale fire events.

The Global Fire Assimilation System (GFASv1.0) calculates biomass burning emissions by assimilating Fire Radiative Power (FRP) observations from the MODIS instruments onboard the Terra and Aqua satellites. It corrects for gaps in the observations, which are mostly due to cloud cover, and filters spurious FRP observations of volcanoes, gas flares and other industrial activity. The combustion rate is subsequently calculated with land cover-specific conversion factors. Emission factors for 40 gas-phase and aerosol trace species have been compiled from a literature survey. The corresponding daily emissions have been calculated on a global 0.5° × 0.5° grid from 2003 to the present. General consistency with the Global Fire Emission Database version 3.1 (GFED3.1) within its accuracy is achieved while maintaining the advantages of an FRP-based approach: GFASv1.0 makes use of the quantitative information on the combustion rate that is contained in the FRP observations, and it detects fires in real time at high spatial and temporal resolution. GFASv1.0 indicates omission errors in GFED3.1 due to undetected small fires. It also exhibits slightly longer fire seasons in South America and North Africa and a slightly shorter fire season in Southeast Asia. GFASv1.0 has already been used for atmospheric reactive gas simulations in an independent study, which found good agreement with atmospheric observations. We have performed simulations of the atmospheric aerosol distribution with and without the assimilation of MODIS aerosol optical depth (AOD). They indicate that the emissions of particulate matter need to be boosted by a factor of 2–4 to reproduce the global distribution of organic matter and black carbon. This discrepancy is also evident in the comparison of previously published top-down and bottom-up estimates. For the time being, a global enhancement of the particulate matter emissions by 3.4 is recommended. Validation with independent AOD and PM10 observations recorded during the Russian fires in summer 2010 show that the global Monitoring Atmospheric Composition and Change (MACC) aerosol model with GFASv1.0 aerosol emissions captures the smoke plume evolution well when organic matter and black carbon are enhanced by the recommended factor. In conjunction with the assimilation of MODIS AOD, the use of GFASv1.0 with enhanced emission factors quantitatively improves the forecast of the aerosol load near the surface sufficiently to allow air quality warnings with a lead time of up to four days.

The Preliminary VOCALS Model Assessment (PreVOCA) aims to assess contemporary atmospheric modeling of the subtropical South East Pacific, with a particular focus on the clouds and the marine boundary layer (MBL). Models results from fourteen modeling centers were collected including operational forecast models, regional models, and global climate models for the month of October 2006. Forecast models and global climate models produced daily forecasts, while most regional models were run continuously during the study period, initialized and forced at the boundaries with global model analyses. Results are compared in the region from 40° S to the equator and from 110° W to 70° W, corresponding to the Pacific coast of South America. Mean-monthly model surface winds agree well with QuikSCAT observed winds and models agree fairly well on mean weak large-scale subsidence in the region next to the coast. However they have greatly differing geographic patterns of mean cloud fraction with only a few models agreeing well with MODIS observations. Most models also underestimate the MBL depth by several hundred meters in the eastern part of the study region. The diurnal cycle of liquid water path is underestimated by most models at the 85° W 20° S stratus buoy site compared with satellite, consistent with previous modeling studies. The low cloud fraction is also underestimated during all parts of the diurnal cycle compared to surface-based climatologies. Most models qualitatively capture the MBL deepening around 15 October 2006 at the stratus buoy, associated with colder air at 700 hPa.

SCIAMACHY limb scatter radiance measurements at selected wavelengths in the HARTLEY bands have been used to retrieve ozone profiles in the upper stratosphere and lower mesosphere. Comparisons with profiles measured by a ground based radiometer in Norway, MIPAS on board ENVISAT, HALOE on UARS and MLS on AURA indicate an agreement within 15% between 40 and 55 km and show that the retrieval provides reliable ozone profiles at these altitudes. Above 55 km, an increasing overestimation is observed. Beside the profile comparisons, further retrieval features of the current retrieval (version 1.26) are described.

Nitrogen dioxide (NO2) vertical tropospheric column densities (VTCs) retrieved from the Global Ozone Monitoring Experiment (GOME) are compared to coincident ground-based tropospheric NO2 columns. The ground-based columns are deduced from in situ measurements at different altitudes in the Alps for 1997 to June 2003, yielding a unique long-term comparison of GOME NO2 VTC data retrieved by a collaboration of KNMI (Royal Netherlands Meteorological Institute) and BIRA/IASB (Belgian Institute for Space Aeronomy) with independently derived tropospheric NO2 profiles. A first comparison relates the GOME retrieved tropospheric columns to the tropospheric columns obtained by integrating the ground-based NO2 measurements. For a second comparison, the tropospheric profiles constructed from the ground-based measurements are first multiplied with the averaging kernel (AK) of the GOME retrieval. The second approach makes the comparison independent from the a priori NO2 profile used in the GOME retrieval. This allows splitting the total difference between the column data sets into two contributions: one that is due to differences between the a priori and the ground-based NO2 profile shapes, and another that can be attributed to uncertainties in both the remaining retrieval parameters (such as, e.g., surface albedo or aerosol concentration) and the ground-based in situ NO2 profiles. For anticyclonic clear sky conditions the comparison indicates a good agreement between the columns (n=157, R=0.70/0.74 for the first/second comparison approach, respectively). The mean relative difference (with respect to the ground-based columns) is −7% with a standard deviation of 40% and GOME on average slightly underestimating the ground-based columns. Both data sets show a similar seasonal behaviour with a distinct maximum of spring NO2 VTCs. Further analysis indicates small GOME columns being systematically smaller than the ground-based ones. The influence of different shapes in the a priori and the ground-based NO2 profile is analysed by considering AK information. It is moderate and indicates similar shapes of the profiles for clear sky conditions. Only for large GOME columns, differences between the profile shapes explain the larger part of the relative difference. In contrast, the other error sources give rise to the larger relative differences found towards smaller columns. Further, for the clear sky cases, errors from different sources are found to compensate each other partially. The comparison for cloudy cases indicates a poorer agreement between the columns (n=60, R=0.61). The mean relative difference between the columns is 60% with a standard deviation of 118% and GOME on average overestimating the ground-based columns. The clear improvement after inclusion of AK information (n=60, R=0.87) suggests larger errors in the a priori NO2 profiles under cloudy conditions and demonstrates the importance of using accurate profile information for (partially) clouded scenes.

High tropospheric NO2 amounts are occasionally detected by space-borne spectrometers above cloudy scenes. For monitoring of near-ground air pollution such data are not directly applicable because clouds shield the highly polluted planetary boundary layer (PBL). We present a method based on trajectories which implicitly estimates the additional sub-cloud NO2 distribution in order to model concentrations at ground stations. The method is applied to a transboundary pollution transport episode which led to high NO2 vertical tropospheric column densities (VTCs) over middle Europe observed by the Global Ozone Monitoring Experiment (GOME) instrument above clouds on 17 February 2001. The case study shows that pollution originally residing near the ground in central Germany, the Ruhr area and adjacent parts of the Netherlands and Belgium has been advected to higher tropospheric levels by a passing weather front. Combining the above-cloud NO2 VTCs with trajectory information covering the GOME columns and including their sub-cloud part yields an estimate of the total NO2 distribution within the tropospheric columns. The highly polluted air masses are then traced by forward trajectories starting from the GOME columns to move further to the Alpine region and their impact there is assessed. Considering ground-based in-situ measurements in the Alpine region, we conclude that for this episode, at least 50% of the NO2 concentration recorded at the sites can be attributed to transboundary transport during the frontal passage. This study demonstrates the potential of using NO2 VTCs from GOME detected above clouds when combined with transport modelling.

Limb scattering retrievals of atmospheric minor constituent profiles require highly accurate knowledge of the tangent heights during the measurements. The limb scattering measurements of the Scanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) on Envisat are affected by tangent height errors of up to 2 km. This contribution provides a summary of the temporal and spatial variation of the SCIAMACHY limb pointing errors during the first three years of the SCIAMACHY mission. The tangent height errors are retrieved from the limb measurements in the UV-B spectral range. A seasonal modulation of the monthly mean tangent height offsets is identified with amplitudes of 800m (220m) before (after) the improvement of the Envisat orbit propagator model in December 2003. Even after the December 2003 orbit model improvement a constant offset component of about 1km is present. Furthermore, pointing discontinuities are identified that coincide with the daily updates of the on-board orbit propagator model. In order to reduce the errors in ozone profile retrievals caused by pointing errors to less than 5%, the tangent heights have to be known to within 250m.