ALBERT

Description: The 2008 Antarctic ozone hole was one of the largest and most long-lived in recent years. Predictions of the ozone hole were made in near-real time (NRT) and hindcast mode with the Integrated Forecast System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The forecasts were carried out both with and without assimilation of satellite observations from multiple instruments to provide more realistic initial conditions. Three different chemistry schemes were applied for the description of stratospheric ozone chemistry: (i) a linearization of the ozone chemistry, (ii) the stratospheric chemical mechanism of the Model of Ozone and Related Chemical Tracers, version 3, (MOZART-3) and (iii) the relaxation to climatology as implemented in the Transport Model, version 5, (TM5). The IFS uses the latter two schemes by means of a two-way coupled system. Without assimilation, the forecasts showed model-specific shortcomings in predicting start time, extent and duration of the ozone hole. The assimilation of satellite observations from the Microwave Limb Sounder (MLS), the Ozone Monitoring Instrument (OMI), the Solar Backscattering Ultraviolet radiometer (SBUV-2) and the SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) led to a significant improvement of the forecasts when compared with total columns and vertical profiles from ozone sondes. The combined assimilation of observations from multiple instruments helped to overcome limitations of the ultraviolet (UV) sensors at low solar elevation over Antarctica. The assimilation of data from MLS was crucial to obtain a good agreement with the observed ozone profiles both in the polar stratosphere and troposphere. The ozone analyses by the three model configurations were very similar despite the different underlying chemistry schemes. Using ozone analyses as initial conditions had a very beneficial but variable effect on the predictability of the ozone hole over 15 days. The initialized forecasts with the MOZART-3 chemistry produced the best predictions of the increasing ozone hole whereas the linear scheme showed the best results during the ozonehole closure.

Description: We present a comparison of tropospheric NO2 from OMI measurements to the median of an ensemble of Regional Air Quality (RAQ) models, and an intercomparison of the contributing RAQ models and two global models for the period July 2008–June 2009 over Europe. The model forecasts were produced routinely on a daily basis in the context of the European GEMS ("Global and regional Earth-system (atmosphere) Monitoring using Satellite and in-situ data") project. The tropospheric vertical column of the RAQ ensemble median shows a spatial distribution which agrees well with the OMI NO2 observations, with a correlation r=0.8. This is higher than the correlations from any one of the individual RAQ models, which supports the use of a model ensemble approach for regional air pollution forecasting. The global models show high correlations compared to OMI, but with significantly less spatial detail, due to their coarser resolution. Deviations in the tropospheric NO2 columns of individual RAQ models from the mean were in the range of 20–34% in winter and 40–62% in summer, suggesting that the RAQ ensemble prediction is relatively more uncertain in the summer months. The ensemble median shows a stronger seasonal cycle of NO2 columns than OMI, and the ensemble is on average 50% below the OMI observations in summer, whereas in winter the bias is small. On the other hand the ensemble median shows a somewhat weaker seasonal cycle than NO2 surface observations from the Dutch Air Quality Network, and on average a negative bias of 14%. Full profile information was available for two RAQ models and for the global models. For these models the retrieval averaging kernel was applied. Minor differences are found for area-averaged model columns with and without applying the kernel, which shows that the impact of replacing the a priori profiles by the RAQ model profiles is on average small. However, the contrast between major hotspots and rural areas is stronger for the direct modeled vertical columns than the columns where the averaging kernels are applied, related to a larger relative contribution of the free troposphere and the coarse horizontal resolution in the a priori profiles compared to the RAQ models. In line with validation results reported in the literature, summertime concentrations in the lowermost boundary layer in the a priori profiles from the DOMINO product are significantly larger than the RAQ model concentrations and surface observations over the Netherlands. This affects the profile shape, and contributes to a high bias in OMI tropospheric columns over polluted regions. The global models indicate that the upper troposphere may contribute significantly to the total column and it is important to account for this in comparisons with RAQ models. A combination of upper troposphere model biases, the a priori profile effects and DOMINO product retrieval issues could explain the discrepancy observed between the OMI observations and the ensemble median in summer.

Description: 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.

Description: This paper presents the aerosol budget over Europe in 2006 calculated with the global transport model TM5 coupled to the size-resolved aerosol module M7. Comparison with ground observations indicates that the model reproduces the observed concentrations quite well with an expected slight underestimation of PM10 due to missing emissions (e.g. resuspension). We model that a little less than half of the anthropogenic aerosols emitted in Europe are exported and the rest is removed by deposition. The anthropogenic aerosols are removed mostly by rain (95%) and only 5% is removed by dry deposition. For the larger natural aerosols, especially sea salt, a larger fraction is removed by dry processes (sea salt: 70%, mineral dust: 35%). We model transport of aerosols in the jet stream in the higher atmosphere and an import of Sahara dust from the south at high altitudes. Comparison with optical measurements shows that the model reproduces the Ångström parameter very well, which indicates a correct simulation of the aerosol size distribution. However, we underestimate the aerosol optical depth. Because the surface concentrations are close to the observations, the shortage of aerosol in the model is probably at higher altitudes. We show that the discrepancies are mainly caused by an overestimation of wet-removal rates. To match the observations, the wet-removal rates have to be scaled down by a factor of about 5. In that case the modelled ground-level concentrations of sulphate and sea salt increase by 50% (which deteriorates the match), while other components stay roughly the same. Finally, it is shown that in particular events, improved fire emission estimates may significantly improve the ability of the model to simulate the aerosol optical depth. We stress that discrepancies in aerosol models can be adequately analysed if all models would provide (regional) aerosol budgets, as presented in the current study.

Description: Three global Chemistry Transport Models – MOZART, MOCAGE, and TM5 – as well as MOZART coupled to the IFS meteorological model including assimilation of ozone (O3) and carbon monoxide (CO) satellite column retrievals, have been compared to surface measurements and MOZAIC vertical profiles in the troposphere over Western/Central Europe for summer 2003. The models reproduce the meteorological features and enhancement of pollution during the period 2–14 August, but not fully the ozone and CO mixing ratios measured during that episode. Modified normalised mean biases are around −25% (except ~5% for MOCAGE) in the case of ozone and from −80% to −30% for CO in the boundary layer above Frankfurt. The coupling and assimilation of CO columns from MOPITT overcomes some of the deficiencies in the treatment of transport, chemistry and emissions in MOZART, reducing the negative biases to around 20%. The high reactivity and small dry deposition velocities in MOCAGE seem to be responsible for the overestimation of O3 in this model. Results from sensitivity simulations indicate that an increase of the horizontal resolution to around 1°×1° and potential uncertainties in European anthropogenic emissions or in long-range transport of pollution cannot completely account for the underestimation of CO and O3 found for most models. A process-oriented TM5 sensitivity simulation where soil wetness was reduced results in a decrease in dry deposition fluxes and a subsequent ozone increase larger than the ozone changes due to the previous sensitivity runs. However this latest simulation still underestimates ozone during the heat wave and overestimates it outside that period. Most probably, a combination of the mentioned factors together with underrepresented biogenic emissions in the models, uncertainties in the modelling of vertical/horizontal transport processes in the proximity of the boundary layer as well as limitations of the chemistry schemes are responsible for the underestimation of ozone (overestimation in the case of MOCAGE) and CO found in the models during this extreme pollution event.

Description: Despite the developments in the global modelling of chemistry and of the parameterization of the physical processes, carbon monoxide (CO) concentrations remain underestimated during Northern Hemisphere (NH) winter by most state-of-the-art chemistry transport models. The consequential model bias can in principle originate from either an underestimation of CO sources or an overestimation of its sinks. We address both the role of surface sources and sinks with a series of MOZART (Model for Ozone And Related Tracers) model sensitivity studies for the year 2008 and compare our results to observational data from ground-based stations, satellite observations, and vertical profiles from measurements on passenger aircraft. In our base case simulation using MACCity (Monitoring Atmospheric Composition and Climate project) anthropogenic emissions, the near-surface CO mixing ratios are underestimated in the Northern Hemisphere by more than 20 ppb from December to April, with the largest bias of up to 75 ppb over Europe in January. An increase in global biomass burning or biogenic emissions of CO or volatile organic compounds (VOCs) is not able to reduce the annual course of the model bias and yields concentrations over the Southern Hemisphere which are too high. Raising global annual anthropogenic emissions with a simple scaling factor results in overestimations of surface mixing ratios in most regions all year round. Instead, our results indicate that anthropogenic CO and, possibly, VOC emissions in the MACCity inventory are too low for the industrialized countries only during winter and spring. Reasonable agreement with observations can only be achieved if the CO emissions are adjusted seasonally with regionally varying scaling factors. A part of the model bias could also be eliminated by exchanging the original resistance-type dry deposition scheme with a parameterization for CO uptake by oxidation from soil bacteria and microbes, which reduces the boreal winter dry deposition fluxes. The best match to surface observations, satellite retrievals, and aircraft observations was achieved when the modified dry deposition scheme was combined with increased wintertime road traffic emissions over Europe and North America (factors up to 4.5 and 2, respectively). One reason for the apparent underestimation of emissions may be an exaggerated downward trend in the Representative Concentration Pathway (RCP) 8.5 scenario in these regions between 2000 and 2010, as this scenario was used to extrapolate the MACCity emissions from their base year 2000. This factor is potentially amplified by a lack of knowledge about the seasonality of emissions. A methane lifetime of 9.7 yr for our basic model and 9.8 yr for the optimized simulation agrees well with current estimates of global OH, but we cannot fully exclude a potential effect from errors in the geographical and seasonal distribution of OH concentrations on the modelled CO.

Description: 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.

Description: 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 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 gases analysis fields to document the dataset 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 ozone sondes and ACE-FTS data to within ±10% in most situations. In the troposphere the reanalysis shows biases of −5% to +10% with respect to ozone sondes 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.

Description: 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 analyze 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), 38% (O3) and 64% (NO2) during the first 24 h, 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, 38 and 80%. 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 depending on the assumptions made for forecasted fire emissions. This was well visible from a relatively rapid increase by the root mean square error with respect to ground-based data for AOD, and satellite based NO2. 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.