We present an overview of state-of-the-art chemistry–climate and chemistry transport models that are used within phase 1 of the Chemistry–Climate Model Initiative (CCMI-1). The CCMI aims to conduct a detailed evaluation of participating models using process-oriented diagnostics derived from observations in order to gain confidence in the models' projections of the stratospheric ozone layer, tropospheric composition, air quality, where applicable global climate change, and the interactions between them. Interpretation of these diagnostics requires detailed knowledge of the radiative, chemical, dynamical, and physical processes incorporated in the models. Also an understanding of the degree to which CCMI-1 recommendations for simulations have been followed is necessary to understand model responses to anthropogenic and natural forcing and also to explain inter-model differences. This becomes even more important given the ongoing development and the ever-growing complexity of these models. This paper also provides an overview of the available CCMI-1 simulations with the aim of informing CCMI data users.

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 (O₃), CO, sulfur dioxide (SO₂) and nitrogen dioxide (NO₂) as well as satellite retrievals of CO, tropospheric NO₂ 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 O₃, and wintertime SO₂, and was of a similar accuracy for other evaluated species. C-IFS (CB05) is about 10 times more computationally efficient than IFS-MOZART.

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 (O₃, NO₂, SO₂, CO, PM₁₀, PM_2.5, NO, NH₃, 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, NO₂ and PM₁₀. 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⁻³ on average. Locally, during the studied ozone episodes, the maxima from the ensemble median are often lower than observations by 30–50 μg m⁻³. 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₁₀ for winter 2013–1014. There is an underestimation of most models leading the ensemble median to a mean bias of −4.5 μg m⁻³. The ensemble median fractional gross error is larger for PM₁₀ (~ 0.52) than for ozone and the correlation is lower (~ 0.35 for PM₁₀ and ~ 0.54 for ozone). This is related to a larger spread of the seven model scores for PM₁₀ 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.

>We analyse simulations performed for the Chemistry-Climate Model Initiative (CCMI) to estimate the return dates of the stratospheric ozone layer from depletion caused by anthropogenic stratospheric chlorine and bromine. We consider a total of 155 simulations from 20 models, including a range of sensitivity studies which examine the impact of climate change on ozone recovery. For the control simulations (unconstrained by nudging towards analysed meteorology) there is a large spread (±20 DU in the global average) in the predictions of the absolute ozone column. Therefore, the model results need to be adjusted for biases against historical data. Also, the interannual variability in the model results need to be smoothed in order to provide a reasonably narrow estimate of the range of ozone return dates. Consistent with previous studies, but here for a Representative Concentration Pathway (RCP) of 6.0, these new CCMI simulations project that global total column ozone will return to 1980 values in 2049 (with a 1σ uncertainty of 2043–2055). At Southern Hemisphere mid-latitudes column ozone is projected to return to 1980 values in 2045 (2039–2050), and at Northern Hemisphere mid-latitudes in 2032 (2020–2044). In the polar regions, the return dates are 2060 (2055–2066) in the Antarctic in October and 2034 (2025–2043) in the Arctic in March. The earlier return dates in the Northern Hemisphere reflect the larger sensitivity to dynamical changes. Our estimates of return dates are later than those presented in the 2014 Ozone Assessment by approximately 5–17 years, depending on the region, with the previous best estimates often falling outside of our uncertainty range. In the tropics only around half the models predict a return of ozone to 1980 values, around 2040, while the other half do not reach the 1980 value. All models show a negative trend in tropical total column ozone towards the end of the 21st century. The CCMI models generally agree in their simulation of the time evolution of stratospheric chlorine and bromine, which are the main drivers of ozone loss and recovery. However, there are a few outliers which show that the multi-model mean results for ozone recovery are not as tightly constrained as possible. Throughout the stratosphere the spread of ozone return dates to 1980 values between models tends to correlate with the spread of the return of inorganic chlorine to 1980 values. In the upper stratosphere, greenhouse gas-induced cooling speeds up the return by about 10–20 years. In the lower stratosphere, and for the column, there is a more direct link in the timing of the return dates of ozone and chlorine, especially for the large Antarctic depletion. Comparisons of total column ozone between the models is affected by different predictions of the evolution of tropospheric ozone within the same scenario, presumably due to differing treatment of tropospheric chemistry. Therefore, for many scenarios, clear conclusions can only be drawn for stratospheric ozone columns rather than the total column. As noted by previous studies, the timing of ozone recovery is affected by the evolution of N₂O and CH₄. However, quantifying the effect in the simulations analysed here is limited by the few realisations available for these experiments compared to internal model variability. The large increase in N₂O given in RCP 6.0 extends the ozone return globally by ∼ 15 years relative to N₂O fixed at 1960 abundances, mainly because it allows tropical column ozone to be depleted. The effect in extratropical latitudes is much smaller. The large increase in CH₄ given in the RCP 8.5 scenario compared to RCP 6.0 also lengthens ozone return by ∼ 15 years, again mainly through its impact in the tropics. Overall, our estimates of ozone return dates are uncertain due to both uncertainties in future scenarios, in particular those of greenhouse gases, and uncertainties in models. The scenario uncertainty is small in the short term but increases with time, and becomes large by the end of the century. There are still some model–model differences related to well-known processes which affect ozone recovery. Efforts need to continue to ensure that models used for assessment purposes accurately represent stratospheric chemistry and the prescribed scenarios of ozone-depleting substances, and only those models are used to calculate return dates. For future assessments of single forcing or combined effects of CO₂, CH₄, and N₂O on the stratospheric column ozone return dates, this work suggests that it is more important to have multi-member (at least three) ensembles for each scenario from every established participating model, rather than a large number of individual models.

The preprocessor PREP-CHEM-SRC presented in the paper is a comprehensive tool aiming at preparing emission fields of trace gases and aerosols for use in atmospheric-chemistry transport models. The considered emissions are from the most recent databases of urban/industrial, biogenic, biomass burning, volcanic, biofuel use and burning from agricultural waste sources. For biomass burning, emissions can be also estimated directly from satellite fire detections using a fire emission model included in the tool. The preprocessor provides emission fields interpolated onto the transport model grid. Several map projections can be chosen. The inclusion of these emissions in transport models is also presented. The preprocessor is coded using Fortran90 and C and is driven by a namelist allowing the user to choose the type of emissions and the databases.

The modeling study presented here aims to estimate how uncertainties in global hydroxyl radical (OH) distributions, variability, and trends may contribute to resolving discrepancies between simulated and observed methane (CH₄) changes since 2000. A multi-model ensemble of 14 OH fields was analyzed and aggregated into 64 scenarios to force the offline atmospheric chemistry transport model LMDz (Laboratoire de Meteorologie Dynamique) with a standard CH₄ emission scenario over the period 2000–2016. The multi-model simulated global volume-weighted tropospheric mean OH concentration ([OH]) averaged over 2000–2010 ranges between 8.7×10⁵ and 12.8×10⁵ molec cm⁻³. The inter-model differences in tropospheric OH burden and vertical distributions are mainly determined by the differences in the nitrogen oxide (NO) distributions, while the spatial discrepancies between OH fields are mostly due to differences in natural emissions and volatile organic compound (VOC) chemistry. From 2000 to 2010, most simulated OH fields show an increase of 0.1–0.3×10⁵ molec cm⁻³ in the tropospheric mean [OH], with year-to-year variations much smaller than during the historical period 1960–2000. Once ingested into the LMDz model, these OH changes translated into a 5 to 15 ppbv reduction in the CH₄ mixing ratio in 2010, which represents 7 %–20 % of the model-simulated CH₄ increase due to surface emissions. Between 2010 and 2016, the ensemble of simulations showed that OH changes could lead to a CH₄ mixing ratio uncertainty of $>\pm \mathrm{30}$ ppbv. Over the full 2000–2016 time period, using a common state-of-the-art but nonoptimized emission scenario, the impact of [OH] changes tested here can explain up to 54 % of the gap between model simulations and observations. This result emphasizes the importance of better representing OH abundance and variations in CH₄ forward simulations and emission optimizations performed by atmospheric inversions.