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. Partial lower tropospheric ozone columns provided by the IASI (Infrared Atmospheric Sounding Interferometer) instrument have been assimilated into a chemistrytransport model at continental scale (CHIMERE) using an Ensemble Square Root Kalman Filter (EnSRF). Analyses are made for the month of July 2007 over the European domain. Launched in 2006, aboard the MetOp-A satellite, IASI shows high sensitivity for ozone in the free troposphere and low sensitivity at the ground; therefore it is important to evaluate if assimilation of these observations can improve free tropospheric ozone, and possibly surface ozone. The analyses are validated against independent ozone observations from sondes, MOZAIC 1 aircraft and ground based stations (AIR-BASE -the European Air quality dataBase) and compared with respect to the free run of CHIMERE. These comparisons show a decrease in error of 6 parts-per-billion (ppb) in the free troposphere over the Frankfurt area, and also a reduction of the root mean square error (respectively bias) at the surface of 19 % (33 %) for more than 90 % of existing ground stations. This provides evidence of the potential of data assimilation of tropospheric IASI columns to better describe the tropospheric ozone distribution, including surface ozone, despite the lower sensitivity.The changes in concentration resulting from the observational constraints were quantified and several geophysical explanations for the findings of this study were drawn. The corrections were most pronounced over Italy and the Mediter-1 Measurements of OZone, water vapour, carbon monoxide and nitrogen oxides by in-service AIrbus airCraft (http://mozaic.aero. obs-mip.fr/web/). ranean region, we noted an average reduction of 8-9 ppb in the free troposphere with respect to the free run, and still a reduction of 5.5 ppb at ground, likely due to a longer residence time of air masses in this part associated to the general circulation pattern (i.e. dominant western circulation) and to persistent anticyclonic conditions over the Mediterranean basin. This is an important geophysical result, since the ozone burden is large over this area, with impact on the radiative balance and air quality.

Abstract. An ensemble Kalman filter (EnKF) has been coupled to the CHIMERE chemical transport model in order to assimilate ozone ground-based measurements on a regional scale. The number of ensembles is reduced to 20, which allows for future operational use of the system for air quality analysis and forecast. Observation sites of the European ozone monitoring network have been classified using criteria on ozone temporal variability, based on previous work by Flemming et al. (2005). This leads to the choice of specific subsets of suburban, rural and remote sites for data assimilation and for evaluation of the reference run and the assimilation system. For a 10-day experiment during an ozone pollution event over Western Europe, data assimilation allows for a significant improvement in ozone fields: the RMSE is reduced by about a third with respect to the reference run, and the hourly correlation coefficient is increased from 0.75 to 0.87. Several sensitivity tests focus on an a posteriori diagnostic estimation of errors associated with the background estimate and with the spatial representativeness of observations. A strong diurnal cycle of both these errors with an amplitude up to a factor of 2 is made evident. Therefore, the hourly ozone background error and the observation error variances are corrected online in separate assimilation experiments. These adjusted background and observational error variances provide a better uncertainty estimate, as verified by using statistics based on the reduced centered random variable. Over the studied 10-day period the overall EnKF performance over evaluation stations is found relatively unaffected by different formulations of observation and simulation errors, probably due to the large density of observation sites. From these sensitivity tests, an optimal configuration was chosen for an assimilation experiment extended over a three-month summer period. It shows a similarly good performance as the 10-day experiment.