We evaluate public health and climate impacts of low-sulphur fuels in global shipping. Using high-resolution emissions inventories, integrated atmospheric models, and health risk functions, we assess ship-related PM2.5 pollution impacts in 2020 with and without the use of low-sulphur fuels. Cleaner marine fuels will reduce ship-related premature mortality and morbidity by 34 and 54%, respectively, representing a ~ 2.6% global reduction in PM2.5 cardiovascular and lung cancer deaths and a ~3.6% global reduction in childhood asthma. Despite these reductions, low-sulphur marine fuels will still account for ~250k deaths and ~6.4 M childhood asthma cases annually, and more stringent standards beyond 2020 may provide additional health benefits. Lower sulphur fuels also reduce radiative cooling from ship aerosols by ~80%, equating to a ~3% increase in current estimates of total anthropogenic forcing. Therefore, stronger international shipping policies may need to achieve climate and health targets by jointly reducing greenhouse gases and air pollution.

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. We present a comparison of tropospheric NO 2 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 NO 2 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 Correspondence to: V. Huijnen (huijnen@knmi.nl) to OMI, but with significantly less spatial detail, due to their coarser resolution. Deviations in the tropospheric NO 2 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 NO 2 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 NO 2 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.Published by Copernicus Publications on behalf of the European Geosciences Union. V. Huijnen et al.:Comparison of NO 2 in regional and global models to OMI 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 effec...

More than ten state-of-the-art regional air quality models have been applied as part of the Air Quality Model Evaluation International Initiative (AQMEII). These models were run by twenty independent groups in Europe and North America. Standardised modelling outputs over a full year (2006) from each group have been shared on the web-distributed ENSEMBLE system, which allows for statistical and ensemble analyses to be performed by each group. The estimated ground-level ozone mixing ratios from the models are collectively examined in an ensemble fashion and evaluated against a large set of observations from both continents. The scale of the exercise is unprecedented and offers a unique opportunity to investigate methodologies for generating skilful ensembles of regional air quality models outputs. Despite the remarkable progress of ensemble air quality modelling over the past decade, there are still outstanding questions regarding this technique. Among them, what is the best and most beneficial way to build an ensemble of members? And how should the optimum size of the ensemble be determined in order to capture data variability as well as keeping the error low? These questions are addressed here by looking at optimal ensemble size and quality of the members. The analysis carried out is based on systematic minimization of the model error and is important for performing diagnostic/probabilistic model evaluation. It is shown that the most commonly used multi-model approach, namely the average over all available members, can be outperformed by subsets of members optimally selected in terms of bias, error, and correlation. More importantly, this result does not strictly depend on the skill of the individual members, but may require the inclusion of low-ranking skill-score members. A clustering methodology is applied to discern among members and to build a skilful ensemble based on model association and data clustering, which makes no use of priori knowledge of model skill. Results show that, while the methodology needs further refinement, by optimally selecting the cluster distance and association criteria, this approach can be useful for model applications beyond those strictly related to model evaluation, such as air quality forecasting

Abstract. Data assimilation is used in atmospheric chemistry models to improve air quality forecasts, construct re-analyses of three-dimensional chemical (including aerosol) concentrations and perform inverse modeling of input variables or model parameters (e.g., emissions). Coupled chemistry meteorology models (CCMM) are atmospheric chemistry models that simulate meteorological processes and chemical transformations jointly. They offer the possibility to assimilate both meteorological and chemical data; however, because CCMM are fairly recent, data assimilation in CCMM has been limited to date. We review here the current status of data assimilation in atmospheric chemistry models with a particular focus on future prospects for data assimilation in CCMM. We first review the methods available for data assimilation in atmospheric models, including variational methods, ensemble Kalman filters, and hybrid methods. Next, we review past applications that have included chemical data assimilation in chemical transport models

Ten state-of-the-science regional air quality (AQ) modeling systems have been applied to continental-scale domains in North America and Europe for full-year simulations of 2006 in the context of Air Quality Model Evaluation International Initiative (AQMEII), whose main goals are model intercomparison and evaluation. Standardised modeling outputs from each group have been shared on the web-distributed ENSEMBLE system, which allows statistical and ensemble analyses to be performed. In this study, the one-year model simulations are inter-compared and evaluated with a large set of observations for ground-level particulate matter (PK10 and PM2.5) and its chemical components. Modeled concentrations of gaseous PM precursors, SO2 and NO2, have also been evaluated against observational data for both continents. Furthermore, modeled deposition (dry and wet) and emissions of several species relevant to PM are also inter-compared. The unprecedented scale of the exercise (two continents, one full year, fifteen modeling groups) allows for a detailed description of AQ model skill and uncertainty with respect to PM. Analyses of PM10 yearly time series and mean diurnal cycle show a large underestimation throughout the year for the AQ models included in AQMEII. The possible causes of PM bias, including errors in the emissions and meteorological inputs (e.g., wind speed and precipitation), and the calculated deposition are investigated. Further analysis of the coarse PM components, PM2.5 and its major components (SO4, NH4, NO3, elemental carbon), have also been performed, and the model performance for each component evaluated against measurements. Finally, the ability of the models to capture high PM concentrations has been evaluated by examining two separate PM2.5 episodes in Europe and North America. A large variability among models in predicting emissions, deposition, and concentration of PM and its precursors during the episodes has been found. Major challenges still remain with regards to identifying and eliminating the sources of PM bias in the models. Although PM2.5 was found to be much better estimated by the models than PM10, no model was found to consistently match the observations for all locations throughout the entire yea

Abstract. The paper presents the transport module of the System for Integrated modeLling of Atmospheric coMposition SILAM v.5 based on the advection algorithm of Michael Galperin. This advection routine, so far weakly presented in the international literature, is positively defined, stable at any Courant number, and efficient computationally. We present the rigorous description of its original version, along with several updates that improve its monotonicity and shape preservation, allowing for applications to long-living species in conditions of complex atmospheric flows. The scheme is connected with other parts of the model in a way that preserves the sub-grid mass distribution information that is a cornerstone of the advection algorithm. The other parts include the previously developed vertical diffusion algorithm combined with dry deposition, a meteorological preprocessor, and chemical transformation modules.The quality of the advection routine is evaluated using a large set of tests. The original approach has been previously compared with several classic algorithms widely used in operational dispersion models. The basic tests were repeated for the updated scheme and extended with real-wind simulations and demanding global 2-D tests recently suggested in the literature, which allowed one to position the scheme with regard to sophisticated state-of-the-art approaches. The advection scheme performance was fully comparable with other algorithms, with a modest computational cost.

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.