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 PM₁₀ 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.

This study presents the new aerosol assimilation system, developed at the European Centre for Medium‐Range Weather Forecasts, for the Global and regional Earth‐system Monitoring using Satellite and in‐situ data (GEMS) project. The aerosol modeling and analysis system is fully integrated in the operational four‐dimensional assimilation apparatus. Its purpose is to produce aerosol forecasts and reanalyses of aerosol fields using optical depth data from satellite sensors. This paper is the second of a series which describes the GEMS aerosol effort. It focuses on the theoretical architecture and practical implementation of the aerosol assimilation system. It also provides a discussion of the background errors and observations errors for the aerosol fields, and presents a subset of results from the 2‐year reanalysis which has been run for 2003 and 2004 using data from the Moderate Resolution Imaging Spectroradiometer on the Aqua and Terra satellites. Independent data sets are used to show that despite some compromises that have been made for feasibility reasons in regards to the choice of control variable and error characteristics, the analysis is very skillful in drawing to the observations and in improving the forecasts of aerosol optical depth.

This paper presents the aerosol modeling now part of the ECMWF Integrated Forecasting System (IFS). It includes new prognostic variables for the mass of sea salt, dust, organic matter and black carbon, and sulphate aerosols, interactive with both the dynamics and the physics of the model. It details the various parameterizations used in the IFS to account for the presence of tropospheric aerosols. Details are given of the various formulations and data sets for the sources of the different aerosols and of the parameterizations describing their sinks. Comparisons of monthly mean and daily aerosol quantities like optical depths against satellite and surface observations are presented. The capability of the forecast model to simulate aerosol events is illustrated through comparisons of dust plume events. The ECMWF IFS provides a good description of the horizontal distribution and temporal variability of the main aerosol types. The forecast‐only model described here generally gives the total aerosol optical depth within 0.12 of the relevant observations and can therefore provide the background trajectory information for the aerosol assimilation system described in part 2 of this paper.

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. NO₂ fields from the reanalysis show the right seasonality over polluted urban areas of the NH and over tropical biomass burning areas, but underestimate wintertime NO₂ maxima over anthropogenic pollution regions and overestimate NO₂ 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.

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⁻¹ (range 34–144 Tg a⁻¹) and the median SOA source strength (natural and anthropogenic) is 19 Tg a⁻¹ (range 13–121 Tg a⁻¹). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a⁻¹ (range 16–121 Tg a⁻¹), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a⁻¹; range 13–20 Tg a⁻¹, with one model at 37 Tg a⁻¹). 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⁻¹ (range 28–209 Tg a⁻¹), 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.

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 10⁵ 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⁻². The RMSE ranges from 20 W m⁻² (3% of the mean observed irradiance) to 36 W m⁻² (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⁻². The root mean square error (RMSE) ranges from 33 W m⁻² (5%) to 64 W m⁻² (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.