Abstract. A representation of atmospheric chemistry has been included in the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The new chemistry modules complement the aerosol modules of the IFS for atmospheric composition, which is named C-IFS. C-IFS for chemistry supersedes a coupled system in which chemical transport model (CTM) Model for OZone and Related chemical Tracers 3 was two-way coupled to the IFS (IFS-MOZART). This paper contains a description of the new on-line implementation, an evaluation with observations and a comparison of the performance of C-IFS with MOZART and with a re-analysis of atmospheric composition produced by IFS-MOZART within the Monitoring Atmospheric Composition and Climate (MACC) project. The chemical mechanism of C-IFS is an extended version of the Carbon Bond 2005 (CB05) chemical mechanism as implemented in CTM Transport Model 5 (TM5). CB05 describes tropospheric chemistry with 54 species and 126 reactions. Wet deposition and lightning nitrogen monoxide (NO) emissions are modelled in C-IFS using the detailed input of the IFS physics package. A 1 year simulation by C-IFS, MOZART and the MACC re-analysis is evaluated against ozonesondes, carbon monoxide (CO) aircraft profiles, European surface observations of ozone (O3), CO, sulfur dioxide (SO2) and nitrogen dioxide (NO2) as well as satellite retrievals of CO, tropospheric NO2 and formaldehyde. Anthropogenic emissions from the MACC/CityZen (MACCity) inventory and biomass burning emissions from the Global Fire Assimilation System (GFAS) data set were used in the simulations by both C-IFS and MOZART. C-IFS (CB05) showed an improved performance with respect to MOZART for CO, upper tropospheric O3, and wintertime SO2, and was of a similar accuracy for other evaluated species. C-IFS (CB05) is about 10 times more computationally efficient than IFS-MOZART.

Abstract. 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 multimodel ensemble. The MACC-II system provides daily 96 h forecasts with hourly outputs of 10 chemical species/aerosols (O 3 , NO 2 , SO 2 , CO, PM 10 , PM 2.5 , NO, NH 3 , total NMVOCs (non-methane volatile organic compounds) and PAN+PAN Published by Copernicus Publications on behalf of the European Geosciences Union. V. Marécal et al.:A regional air quality forecasting system over Europe 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, NO 2 and PM 10 . 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 PM 10 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 PM 10 (∼ 0.52) than for ozone and the correlation is lower (∼ 0.35 for PM 10 and ∼ 0.54 for ...

Abstract. This paper deals with recent improvements to the global chemical transport model of Météo-France MOCAGE (Modèle de Chimie Atmosphérique à Grande Echelle) that consists of updates to different aerosol parameterizations. MOCAGE only contains primary aerosol species: desert dust, sea salt, black carbon, organic carbon, and also volcanic ash in the case of large volcanic eruptions. We introduced important changes to the aerosol parameterization concerning emissions, wet deposition and sedimentation. For the emissions, size distribution and wind calculations are modified for desert dust aerosols, and a surface sea temperature dependant source function is introduced for sea salt aerosols. Wet deposition is modified toward a more physically realistic representation by introducing re-evaporation of falling rain and snowfall scavenging and by changing the in-cloud scavenging scheme along with calculations of precipitation cloud cover and rain properties. The sedimentation scheme update includes changes regarding the stability and viscosity calculations. Independent data from satellites (MODIS, SEVIRI), the ground (AERONET, EMEP), and a model inter-comparison project (AeroCom) are compared with MOCAGE simulations and show that the introduced changes brought a significant improvement on aerosol representation, properties and global distribution. Emitted quantities of desert dust and sea salt, as well their lifetimes, moved closer towards values of AeroCom estimates and the multimodel average. When comparing the model simulations with MODIS aerosol optical depth (AOD) observations over the oceans, the updated model configuration shows a decrease in the modified normalized mean bias (MNMB; from 0.42 to 0.10) and a better correlation (from 0.06 to 0.32) in terms of the geographical distribution and the temporal variability. The updates corrected a strong positive MNMB in the sea salt representation at high latitudes (from 0.65 to 0.16), and a negative MNMB in the desert dust representation in the African dust outflow region (from − 1.01 to −0.22). The updates in sedimentation produced a modest difference; the MNMB with MODIS data from 0.10 in the updated configuration went to 0.11 in the updated configuration only without the sedimentation updates. Yet, the updates in the emissions and the wet deposition made a stronger impact on the results; the MNMB was 0.27 and 0.21 in updated configurations only without emission, and only without wet deposition updates, respectively. Also, the lifetime, the extent, and the strength of the episodic aerosol events are better reproduced in the updated configuration. The wet deposition processes and the differences between the various configurations that were tested greatly influence the representation of the episodic events. However, wet deposition is not a continuous process; it has a local and episodic signature and its representation depends strongly on the precipitation regime in the model.

Abstract. We report on an evaluation of tropospheric ozone and its precursor gases in three atmospheric chemistry versions as implemented in the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), referred to as IFS(CB05BASCOE), IFS(MOZART) and IFS(MOCAGE). While the model versions were forced with the same overall meteorology, emissions, transport and deposition schemes, they vary largely in their parameterisations describing atmospheric chemistry, including the organics degradation, heterogeneous chemistry and photolysis, as well as chemical solver. The model results from the three chemistry versions are compared against a range of aircraft field campaigns, surface observations, ozone-sondes and satellite observations, which provides quantification of the overall model uncertainty driven by the chemistry parameterisations. We find that they produce similar patterns and magnitudes for carbon monoxide (CO) and ozone (O3), as well as a range of non-methane hydrocarbons (NMHCs), with averaged differences for O3 (CO) within 10 % (20 %) throughout the troposphere. Most of the divergence in the magnitude of CO and NMHCs can be explained by differences in OH concentrations, which can reach up to 50 %, particularly at high latitudes. There are also comparatively large discrepancies between model versions for NO2, SO2 and HNO3, which are strongly influenced by secondary chemical production and loss. Other common biases in CO and NMHCs are mainly attributed to uncertainties in their emissions. This configuration of having various chemistry versions within IFS provides a quantification of uncertainties induced by chemistry modelling in the main CAMS global trace gas products beyond those that are constrained by data assimilation.

Abstract. A representation of atmospheric chemistry has been included in the Integrated Forecasting System (IFS) of the European Centre for Medium-range Weather Forecasts (ECMWF). The new chemistry modules complement the aerosol modules of the IFS for atmospheric composition, which is named C-IFS. C-IFS for chemistry supersedes a coupled system, in which the Chemical Transport Model (CTM) Model for OZone and Related chemical Tracers 3 was two-way coupled to the IFS (IFS-MOZART). This paper contains a description of the new on-line implementation, an evaluation with observations and a comparison of the performance of C-IFS with MOZART and with a re-analysis of atmospheric composition produced by IFS-MOZART within the Monitoring Atmospheric Composition and Climate (MACC) project. The chemical mechanism of C-IFS is an extended version of the Carbon Bond 2005 (CB05) chemical mechanism as implemented in the CTM Transport Model 5 (TM5). CB05 describes tropospheric chemistry with 54 species and 126 reactions. Wet deposition and lightning nitrogen monoxide (NO) emissions are modelled in C-IFS using the detailed input of the IFS physics package. A one-year simulation by C-IFS, MOZART and the MACC re-analysis is evaluated against ozonesondes, carbon monoxide (CO) aircraft profiles, European surface observations of ozone (O3), CO, sulphur dioxide (SO2) and nitrogen dioxide (NO2) as well as satellite retrievals of CO, tropospheric NO2 and formaldehyde. Anthropogenic emissions from the MACC/CityZen (MACCity) inventory and biomass burning emissions from the Global Fire Assimilation System (GFAS) data set were used in the simulations by both C-IFS and MOZART. C-IFS (CB05) showed an improved performance with respect to MOZART for CO, upper tropospheric O3, winter time SO2 and was of a similar accuracy for other evaluated species. C-IFS (CB05) is about ten times more computationally efficient than IFS-MOZART.

Abstract. The tropical transport processes of 14 different models or model versions were compared, within the framework of the SCOUT-O3 (Stratospheric-Climate Links with Emphasis on the Upper Troposphere and Lower Stratosphere) project. The tested models range from the regional to the global scale, and include numerical weather prediction (NWP), chemical transport, and chemistry-climate models. Idealised tracers were used in order to prevent the model's chemistry schemes from influencing the results substantially, so that the effects of modelled transport could be isolated. We find large differences in the vertical transport of very short-lived tracers (with a lifetime of 6 h) within the tropical troposphere. Peak convective outflow altitudes rangeCorrespondence to: C. R. Hoyle (christopher.hoyle@env.ethz.ch) from around 300 hPa to almost 100 hPa among the different models, and the upper tropospheric tracer mixing ratios differ by up to an order of magnitude. The timing of convective events is found to be different between the models, even among those which source their forcing data from the same NWP model (ECMWF). The differences are less pronounced for longer lived tracers, however they could have implications for modelling the halogen burden of the lowermost stratosphere through transport of species such as bromoform, or short-lived hydrocarbons into the lowermost stratosphere. The modelled tracer profiles are strongly influenced by the convective transport parameterisations, and different boundary layer mixing parameterisations also have a large impact on the modelled tracer profiles. Preferential locations for rapid transport from the surface into the upper troposphere are similar in all models, and are mostly concentrated over the western Pacific, the Maritime Continent and the Indian Published by Copernicus Publications on behalf of the European Geosciences Union. Ocean. In contrast, models do not indicate that upward transport is highest over western Africa.

Abstract. This paper presents the first ensemble modelling experiment in relation to birch pollen in Europe. The sevenmodel European ensemble of MACC-ENS, tested in trial simulations over the flowering season of 2010, was run through the flowering season of 2013. The simulations have been compared with observations in 11 countries, all members of the European Aeroallergen Network, for both individual models and the ensemble mean and median. It is shown that the models successfully reproduced the timing of the very late season of 2013, generally within a couple of days from the observed start of the season. The end of the season was generally predicted later than observed, by 5 days or more, which is a known feature of the source term used in the study. Absolute pollen concentrations during the seaPublished by Copernicus Publications on behalf of the European Geosciences Union. M. Sofiev et al.: MACC regional multi-model ensemble simulations of birch pollen dispersionson were somewhat underestimated in the southern part of the birch habitat. In the northern part of Europe, a recordlow pollen season was strongly overestimated by all models. The median of the multi-model ensemble demonstrated robust performance, successfully eliminating the impact of outliers, which was particularly useful since for most models this was the first experience of pollen forecasting.

Abstract. Fast convective transport in the tropics can efficiently redistribute water vapour and pollutants up to the upper troposphere. In this study we compare tropical convection characteristics for the year 2005 in a range of atmospheric models, including numerical weather prediction (NWP) models, chemistry transport models (CTMs), and chemistry-climate models (CCMs). The model runs have been performed within the framework of the SCOUT-O3 (Stratospheric-Climate Links with Emphasis on the Upper Troposphere and Lower Stratosphere) project. The characteristics of tropical convection, such as seasonal cycle, land/sea contrast and vertical extent, are analysed using satellite observations as a benchmark for model simulations. The observational datasets used in this work comprise precipitation rates, outgoing longwave radiation, cloud-top pressure, and water vapour from a number of independent sources, including ERA-Interim analyses. Most models are generally able to reproduce the seasonal cycle and strength of precipitation for continental regions but show larger discrepancies with observations for the Maritime Continent region. The frequency distribution of high clouds from models and observations is calculated using highly temporally-resolved (up to 3-hourly) cloud top data. The percentage of clouds above 15 km varies significantly between the models. Vertical profiles of water vapour in the upper troposphere-lower stratosphere (UTLS)Correspondence to: M. Russo (maria.russo@atm.ch.cam.ac.uk) show large differences between the models which can only be partly attributed to temperature differences. If a convective plume reaches above the level of zero net radiative heating, which is estimated to be ∼15 km in the tropics, the air detrained from it can be transported upwards by radiative heating into the lower stratosphere. In this context, we discuss the role of tropical convection as a precursor for the transport of short-lived species into the lower stratosphere.