Journal topic
Atmos. Chem. Phys., 19, 7019–7053, 2019
https://doi.org/10.5194/acp-19-7019-2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
Atmos. Chem. Phys., 19, 7019–7053, 2019
https://doi.org/10.5194/acp-19-7019-2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 24 May 2019

Research article | 24 May 2019

# Effects of ship emissions on air quality in the Baltic Sea region simulated with three different chemistry transport models

Effects of ship emissions on air quality in the Baltic Sea region simulated with three different chemistry transport models
Matthias Karl1, Jan Eiof Jonson2, Andreas Uppstu3, Armin Aulinger1, Marje Prank3,a, Mikhail Sofiev3, Jukka-Pekka Jalkanen3, Lasse Johansson3, Markus Quante1, and Volker Matthias1 Matthias Karl et al.
• 1Institute of Coastal Research, Helmholtz-Zentrum Geesthacht, 21502 Geesthacht, Germany
• 2Norwegian Meteorological Institute, Oslo, Norway
• 3Atmospheric Composition Research, Finnish Meteorological Institute, P.O. Box 503, 00101 Helsinki, Finland
• anow at: Department of Earth and Atmospheric Sciences, Cornell University, Ithaka, NY, USA

Correspondence: Matthias Karl (matthias.karl@hzg.de)

Abstract
Back to toptop

The Baltic Sea is a highly frequented shipping area with busy shipping lanes close to densely populated regions. Exhaust emissions from ship traffic into the atmosphere do not only enhance air pollution, they also affect the Baltic Sea environment through acidification and eutrophication of marine waters and surrounding terrestrial ecosystems. As part of the European BONUS project SHEBA (Sustainable Shipping and Environment of the Baltic Sea region), the transport, chemical transformation and fate of atmospheric pollutants in the Baltic Sea region were simulated with three regional chemistry transport model (CTM) systems, CMAQ, EMEP/MSC-W and SILAM, with grid resolutions between 4 and 11 km. The main goal was to quantify the effect that shipping emissions have on the regional air quality in the Baltic Sea region when the same shipping emission dataset but different CTMs are used in their typical set-ups. The performance of these models and the shipping contribution to the results of the individual models were evaluated for sulfur dioxide (SO2), nitrogen dioxide (NO2), ozone (O3) and particulate matter (PM2.5). Model results from the three CTMs for total air pollutant concentrations were compared to observations from rural and urban background stations of the AirBase monitoring network in the coastal areas of the Baltic Sea region. Observed PM2.5 in summer was underestimated strongly by CMAQ and to some extent by EMEP/MSC-W. Observed PM2.5 in winter was underestimated by SILAM. In autumn all models were in better agreement with observed PM2.5. The spatial average of the annual mean O3 in the EMEP/MSC-W simulation was ca. 20 % higher compared to the other two simulations, which is mainly the consequence of using a different set of boundary conditions for the European model domain. There are significant differences in the calculated ship contributions to the levels of air pollutants among the three models. EMEP/MSC-W, with the coarsest grid, predicted weaker ozone depletion through NO emissions in the proximity of the main shipping routes than the other two models. The average contribution of ships to PM2.5 levels in coastal land areas is in the range of 3.1 %–5.7 % for the three CTMs. Differences in ship-related PM2.5 between the models are mainly attributed to differences in the schemes for inorganic aerosol formation. Differences in the ship-related elemental carbon (EC) among the CTMs can be explained by differences in the meteorological conditions, atmospheric transport processes and the applied wet-scavenging parameterizations. Overall, results from the present study show the sensitivity of the ship contribution to combined uncertainties in boundary conditions, meteorological data and aerosol formation and deposition schemes. This is an important step towards a more reliable evaluation of policy options regarding emission regulations for ship traffic and the planned introduction of a nitrogen emission control area (NECA) in the Baltic Sea and the North Sea in 2021.

1 Introduction
Back to toptop

International shipping is important for the economic exchange in Europe: almost 90 % of the European Union (EU) import and export freight trade is seaborne. Compared to other modes of transport such as trucks and air freight, shipping is far more energy efficient per ton of cargo. The Baltic Sea is one of the most densely trafficked sea regions in the world. Roughly 407 500 ship crossings in the Baltic Sea were recorded in 2012 (HELCOM2014), with ship types including passenger, cargo and tanker. Maritime transport of goods between the main EU ports and ports located in the Baltic Sea has a share of 22 % (in 2016) of the total shipping tonnages within European seas (EUROSTAT2018).

Ship traffic is associated with exhaust emissions of a wide range of air pollutants, among them nitrogen oxides (NOx=NO+NO2), black carbon (BC), sulfur dioxide (SO2), non-methane volatile organic compounds (NMVOC) and particulate matter, as well as greenhouse gases (mainly carbon dioxide, CO2). The emitted amounts and size spectrum of particulate matter depends on the type of fuel and its sulfur content and the ship engine type (e.g. ). Particulate matter with aerodynamic diameters of less than 10 µm (PM10) and less than 2.5 µm (PM2.5) from ship exhausts has been associated with adverse health effects (e.g. Corbett et al.2007). A global model study by demonstrated the health benefits from reducing the ship-related fine-particulate matter with low-sulfur ship fuels in densely populated, major trading nations.

The atmospheric transformation of primary emitted gases from shipping is especially relevant for the formation of ozone (O3) and the deposition of sulfur and nitrogen compounds distant from the shipping lanes . Emissions from ships are transported in the atmosphere over several hundreds of kilometres . Exhaust emissions from ship traffic in the Baltic Sea has the potential to degrade air quality in the coastal areas and to significantly affect the Baltic Sea environment through acidification and eutrophication of marine waters and surrounding terrestrial ecosystems (). Acidification is a major challenge in the Baltic Sea region today where the critical load (CL) for acidification is exceeded, especially in the southern part . Despite the considerable improvement concerning critical loads with respect to acidification, there are still regions in the Baltic Sea catchment, for which further reductions in acidification are desirable . Atmospheric deposition of nitrogen compounds plays a role in the eutrophication of the coastal marine environment (e.g. Paerl1995) and threatens biodiversity in forests, semi-natural vegetation and freshwater catchments through excessive nitrogen input . Even though exceedances of CL for eutrophication has decreased over the past decades, critical loads are still exceeded in about 65 % of the European ecosystems .

Air pollution from ships is increasingly controlled worldwide by the International Maritime Organization (IMO) through the Marine Pollution Convention (MARPOL) Annex VI – Regulations for the Prevention of Air Pollution from Ships (IMO2008). The Baltic Sea has been a sulfur emission control area (SECA) since May 2006, with stepwise reductions in the sulfur content in ship fuels; from 2015 onwards the sulfur content of any fuel oil used on board ships within the Baltic Sea had to be 0.1 % or less . The effect of regulation on nitrogen emissions has been small until now, as these are only enforced for the new built ships. A nitrogen emission control area (NECA) for the Baltic Sea, North Sea and English Channel will become effective in 2021, but only newly built ships have to comply with the new regulation. The consequences of establishing the NECA on the future air quality in the Baltic Sea region are investigated in the companion paper by .

The transport, chemical transformation and fate of atmospheric pollutants in the Baltic Sea region can be simulated with 3-D chemistry transport model (CTM) systems. Previous air quality model studies related to shipping in the North Sea and Baltic Sea () used CTM systems to investigate the effect of implementation of MARPOL regulations on sulfur emissions by ships, the effect of establishing the NECA and other ship emission control scenarios. The studies quantified the contributions from shipping to the total air concentrations, deposition of nitrogen and sulfur, as well as air quality and health indicators.

studied the effects of shipping on the global scale, including the effects of shipping in the Baltic Sea and the North Sea compared to total anthropogenic emissions in a global CTM with 0.5 × 0.5 resolution. They found that a significant fraction, ranging from 5 % to more than 10 %, of the PM2.5, and the deposition of nitrogen of anthropogenic origin in bordering countries can be attributed to ship emissions in the two sea areas.

Using the EMEP/MSC-W model with a horizontal resolution of 14 × 14 km2 on the regional scale, assessed the effect of reduced sulfur content (2015 value of 0.1 %) and regulation of NECAs on the air quality, deposition of nitrogen and related impacts on human health in the Baltic and North seas. used the Community Multiscale Air Quality (CMAQ) model v4.7.1 with a horizontal resolution of 24 × 24 km2 to investigate the effects of different future developments of shipping emissions in the North Sea area on air quality in the North Sea region. used the SILAM CTM system with spatial resolution of ca. 8 × 8 km2 in a cost-benefit analysis of the sulfur reduction policy in the Baltic Sea but did not investigate the effects of shipping emissions on ozone concentration or nitrogen deposition. The study by used the EMEP/MSC-W model with a much coarser resolution (50 × 50 km2). They find the highest contribution of international shipping in the Baltic Sea and North Sea to ambient levels of air pollutants near large harbour cities and along the main shipping lanes.

The use of relatively coarse model grids (coarser than 10 km resolution) in some of the previous CTM simulations raises concerns about non-linear chemical effects, particularly for O3, since a high source strength from shipping in the proximity to large land-based emissions (inside the same grid cell of the model) often results in very high levels of NOx and excessive ozone titration . Moreover, shipping releases large amounts of NOx from a moving point source within the relatively clean maritime atmosphere. In regional CTMs, these NOx emissions are diluted into large grid volumes, which can lead to a systematic overestimation of the ozone production and artificially increases the lifetime of NOx ().

Despite previous model based work on the effects of shipping on air quality and human health in the Baltic Sea region, there is a need for more localised studies building on a much higher level of details, i.e. concerning shipping activity, for the quantification of regional ship-related air pollution. Knowledge on air quality impacts of shipping with a finer spatial resolution than in previous model studies is required for the identification of best suited sustainable development options for the shipping sector, especially if a varying suite of competing environmental and economic drivers is to be considered in different subregions.

With the goal to support EU policies on environmental and economic aspect of the shipping sector the BONUS project SHEBA (Sustainable Shipping and Environment of the Baltic Sea Region) was established in 2015. The overarching aim of BONUS SHEBA was an integrated and in-depth analysis of the ecological, economic and social impacts of shipping in the Baltic Sea. As part of the SHEBA project, the transport, chemical transformation and fate of atmospheric pollutants in the Baltic Sea region was simulated with three different regional CTM systems (CMAQ, EMEP/MSC-W and SILAM) to investigate the effect of ship emissions on the regional air quality in the Baltic Sea region. The EMEP/MSC-W model (MET Norway) is also included with the same model configuration in the Interreg Baltic Sea region project EnviSuM (Environmental Impact of Low Emission Shipping: Measurements and Modelling Strategies). The main focus of the EnviSuM project is on sulfur emissions from shipping. EnviSuM investigates the effects of the implementation of the stricter SECA from 2015 onwards, combining measurements and modelling. evaluated the skill of air quality models including SILAM, EMEP and CMAQ to reproduce the particulate matter concentration and composition on a European scale. The chosen CTM systems are well established in Europe and have been extensively tested in several multi-model assessment studies (). All three models have been used previously in the North Sea and Baltic Sea region for estimating the effect of shipping (). The model set-up with CMAQ used in has been evaluated for the larger North Sea region .

This study takes a multi-model approach using three CTM systems to assess the uncertainties connected with the atmospheric transport and transformation of air pollutants. The comparison of air concentration of regulatory pollutants between the models is the primary focus of this study. Modelled total air concentrations of NOx, O3, SO2 and PM2.5 from the three CTMs are compared to observations from rural and urban background stations in the coastal areas of the Baltic Sea region. Statistical performance analysis of the comparison of modelled against observation data of total concentrations was carried out for all CTM systems and the performance of the models was intercompared based on several statistical indicators. Specifically, we want to evaluate the contribution of shipping emissions to modelled surface concentrations of important air pollutants. The significance of the ship contribution to ambient NO2 observations at coastal monitoring stations is evaluated for the different models. The use of three CTM systems, together with comparison to ground-based observations, provides a comprehensive view on the current air quality situation of the Baltic Sea region and how it is affected by emissions from shipping. The combination of three CTMs also provides a more robust estimate of the ship-related contribution to ambient atmospheric concentrations.

2 Description of the CTM systems and set-up of the model comparison
Back to toptop

The set-up of the three CTM systems for this study with respect to drivers for meteorology, boundary conditions and emissions was specific for each model system. The models were set up in the same way as they are typically used in air quality studies in European regions. However, the applied CTMs use a much higher spatial and temporal resolution than previous modelling of the air quality in the Baltic Sea region. Shipping activities are considered with a high degree of detail in the simulations, using automatic identification system (AIS) position data and up-to-date load-dependent emission factors for all important air pollutants. The dynamic ship emission inventory, Ship Traffic Emission Assessment Model (STEAM; ), was applied in all CTMs. STEAM takes into account the emission control areas and regulations, emission abatement equipment on board the ships as well as fuel sulfur content modelling separately for main and auxiliary engines (). All three regional air quality models implement state-of-the-art formulations of atmospheric transport, atmospheric chemistry and aerosol formation, which are updated compared to the model versions used in the previous studies. Partly the same or similar drivers for anthropogenic emissions were used in the CTMs. Ship exhaust emissions from the North Sea are handled in the same manner as the Baltic Sea emissions since they affect the western part of the Baltic Sea region. By this procedure it is ensured that all interactions of shipping emissions with pollutants in the regional background and from land-based emission sources are correctly considered. With all models a reference run for the current air quality situation was performed including all emissions (“base”) and one run without the emissions from shipping (“noship”). The difference between the run with all emissions and the run without shipping emissions is used to determine the contribution of ships to the ambient pollutant concentration. Previous calculations have shown that the assumption of linearity, by adding the contributions from different emission sources, is reasonable for ozone and other pollutants and that the associated error is within a few percent ().

## 2.1 Description of the models

### 2.1.1 CMAQ model

The CMAQ model v5.0.1 () computes the air concentration and deposition fluxes of atmospheric gases and aerosols as a consequence of emission, transport and chemical transformation. The atmospheric chemistry is treated by the modified Carbon Bond V mechanism cb05tucl with updated toluene chemistry and chlorine radical chemistry (). The aerosol scheme AERO5 is used for the formation of secondary inorganic aerosol (SIA). The gas phase–aerosol equilibrium partitioning of sulfuric acid (H2SO4), nitric acid (HNO3), hydrochloric acid (HCl) and ammonia (NH3) is solved with the ISORROPIA v1.7 mechanism (). The formation of secondary organic aerosol (SOA) from isoprene, monoterpenes, sesquiterpenes, benzene, toluene, xylene and alkanes () is included.

The dry deposition parameterization is presented in and . Wet deposition of gases and particles is computed by the resolved cloud model of CMAQ, which estimates how much certain vertical model layers contributed to the precipitation . Sea salt emissions were calculated in line by the parameterization of Gong (2003) (as described in ). Sea salt surf zone emissions were deactivated because of considerable overestimations in some coastal regions . Biogenic emissions (NMVOC from vegetation and soil NO) were calculated offline with the biogenic Emission Inventory System BEIS v3.4 (). Emissions of wind-blown dust were not considered.

The Meteorology-Chemistry Interface Processor (MCIP; ) processes meteorological model output into the input format required for CMAQ. The vertical dimension of the model extends up to 100 hPa in a sigma hybrid pressure coordinate system with 30 layers. Twenty of these layers are below approximately 2 km; the lowest layer extends to ca. 42 m above the ground. A spin-up period of 1 month (December 2011) was used for the initialization of the model runs, sufficiently long to prevent initial conditions having an effect on the simulated atmospheric concentrations of the investigated period (year 2012).

### 2.1.2 SILAM model

The SILAM (System for Integrated modeLling of Atmospheric coMposition) model v5.5 (; http://silam.fmi.fi/, last access: 17 April 2019) was used as second CTM in this study. The gas phase chemistry was simulated with the carbon bond (CB) mechanism, CBM-IV, with reaction rates updated according to the recommendations of the International Union of Pure and Applied Chemistry (IUPAC, http://iupac.pole-ether.fr, last access: 17 April 2019) and the NASA Jet Propulsion Laboratory (JPL; http://jpldataeval.jpl.nasa.gov, last access: 17 April 2019) and the terpenes oxidation added from CB05 reaction list (Yarwood et al., 2005). The sulfur chemistry and secondary inorganic aerosol formation is computed with an updated version of the DMAT scheme (Sofiev2000) and secondary organic aerosol formation with the volatility basis set (VBS, ), which is the volatility distribution of anthropogenic organic carbon taken from . Organic aerosol in SILAM was evaluated in a recent study by .

The dry deposition scheme is described in and the wet deposition in . Natural emissions included in the simulations are sea salt emissions as in , biogenic NMVOC emissions as in , wildfire emissions as in and wind-blown desert dust.

SILAM includes a meteorological pre-processor for diagnosing the basic features of the boundary layer and the free troposphere from the meteorological fields provided by various meteorological models . In total 10 vertical layers, extending up to 2000 m above the surface, are included. The lowest layer extends to 20 m above the surface. No spin-up period was applied.

### 2.1.3 EMEP/MSC-W model

The third CTM used in this study is the EMEP/MSC-W chemical transport model, version rv4.8. This model, which is open source (https://github.com/metno/emep-ctm, last access: 17 April 2019), has been described in detail in , with various updates, see and references within. The chemistry scheme of the gas phase in the model is EmChem09, having 72 chemical compounds including 10 “surrogate” VOCs, of which isoprene represents BVOCs and 137 reactions. This scheme is an update of previous chemical schemes (e.g. ). The EMEP scheme involves relatively more details on peroxy radical (RO2) chemistry than, e.g. CB schemes. SOA is calculated using a VBS scheme, which tracks the semi-volatile products of VOC oxidation and dynamically partitions these between the gas and aerosol phases (e.g. Robinson et al.2007). A number of schemes were tested in , but here the standard NPAS scheme as described in is used.

For Europe the model is regularly evaluated against measurements in the EMEP annual reports (see https://www.emep.int, last access: 17 April 2019). In addition the EMEP model has been included in model intercomparisons and model validations in a number of peer-reviewed publications (). Biogenic emissions (NMVOC, soil NO), emissions of dimethyl sulfide (DMS) from oceans, sea salt, dust, road dust, emissions from aviation on cruising altitude, lightning, volcanic emissions and emissions from forest fires are included as separate databases or calculated within the model (EMEP2015).

EMEP is driven by the meteorological data of the European Centre for Medium-Range Weather Forecasts (ECMWF) based on the Cy40r1 version. An important addition to the forecast ensemble in cycle Cy40r1 has been the introduction of ocean–atmosphere coupling from day 0 instead of from day 10 as in the previous cycles. Vertically, the meteorological fields from ECMWF are interpolated onto 20 EMEP sigma levels between the surface and 100 hPa. Initial and boundary concentrations are based on long-term observations and some model data. No spin-up period was applied.

## 2.2 Set-up of the CTM systems

The CTMs were coupled offline with different meteorological models (COSMO-CLM, ECMWF-IFS and WRF). CMAQ and SILAM were operated with high horizontal resolution (4 km) on the inner nest, representing the Baltic Sea region starting from simulations of a coarser European domain. EMEP MSC-W was operated on 0.1 resolution for the whole of Europe. Ship emissions from the STEAM model () were gridded to the respective model's grid resolution. Land-based emissions were from SMOKE-EU or ECLIPSE () databases; annual totals were comparable.

Figure 1Geographic map of the study domain for the CTM comparison, spanning from latitude 53.40 N (south) to 65.80 N (north) and longitude 9.00 E (west) to 31.10 E (east). The area of the Baltic Sea used in this study is shown in blue.

Download

### 2.2.1 Model domains and nesting

The spatial extent for the intercomparison study covers the Baltic Sea region is from latitude 53.50 N (south) to 66.00 N (north) and longitude 9.00 E (west) to 31.00 E (east). Parts of the Kattegat and a small part of the Norwegian Sea are covered by this area but not considered in the comparison. The extent of the geographic domain is displayed in Fig. 1. Nested simulations were done with CMAQ and SILAM models, using the output of the finer resolved inner nest, whereas the simulation with the EMEP/MSC-W model covered the European domain. The SILAM model was operated on rotated grids centred on the respective domain. The horizontal grid resolution of the output was 4 km for CMAQ, 0.04 (∼4km) for SILAM and 0.1 (∼11km) for EMEP/MSC-W. Note that the grid distance in the x direction becomes smaller with increasing latitude (for instance, 0.1 in longitude corresponds to 6.2 km at 56 N). Anthropogenic emissions from the continent and the shipping emissions in the North and Baltic seas were identical (CMAQ and SILAM) or similar in the spatial distribution and magnitude (EMEP/MSC-W). The EMEP/MSC-W model used monthly averaged gridded ship emissions, while the other two models used hourly emissions. Daily or hourly emissions reflect ship traffic pattern changes due to meteorological conditions or due to sea ice. Using a coarser time resolution for shipping thus mainly neglects the influence of weather and ice on ship operations . Table 1 gives an overview of the model set-ups of the three CTM systems for use in the intercomparison study.

Table 1Description of the model set-up of the three CTM systems.

Nested simulations with CMAQ were performed with a coarse outer domain for the whole of Europe with a grid cell size of 64×64km2, an intermediate domain with 16×16km2 for northern Europe and an inner domain with a horizontal resolution of 4×4km2 for the entire Baltic Sea. Model results for the intercomparison were taken from the inner domain for the coastal regions and from the intermediate domain for parts of Sweden, Finland and the Baltic states. For details on the high-resolution output from CMAQ and an evaluation of the model set-up with a limited number of regional background stations, refer to .

For the SILAM model, the grid cell size was roughly 70×70km2 for the outer domain, roughly 18×18km2 for the central domain and roughly 4×4km2 for the inner domain. The simulation time steps were 20, 10 and 4 min. Model results for the intercomparison were mostly from the inner domain, with parts of Finland and eastern Europe from the central domain. SILAM took part in AQMEII 1 and 3 intercomparisons showing comparable performance with other European state-of-the-art air quality models (). The EMEP/MSC-W model was run with a 0.1×0.1 resolution for the whole of Europe. A comprehensive description, including model evaluations, of the model results with the 0.1×0.1 application of the EMEP model for 2013 can be found in .

### 2.2.2 Meteorology

The SILAM model is run with meteorological input from a simulation with the Weather Research and Forecast (WRF) model v3.7.1 using original resolutions of 4.0, 16.0 and 64.0 km, for inner, central and outer domains, respectively. WRF was driven with large-scale meteorological forcing data taken from the ERA-Interim reanalysis . In general, linear interpolation was applied for the simulation, but conservation of mass was used where applicable. The high-resolution inner domain extended up to 2000 m in height and was therefore less influenced by upper-tropospheric dynamics of WRF. , using WRF in a similar configuration, evaluated the WRF meteorological fields against station observations in Poland. The 2 m air temperature (T2) was underestimated in winter (bias smaller than −0.6 K), while temperature in the warm season was overestimated (bias up to +1.0 K). The largest errors for the 10 m wind speed (WS10) occurred in late summer and autumn and the largest errors for wind direction were in spring and summer. The error of wind direction was very small in winter. The spatial distribution of meteorological variables obtained from WRF were in close agreement with the station measurements, but the model performance was found to be worse for the seashore and mountain areas than for other inland areas .

High-resolution meteorological fields for CMAQ were obtained from the COSMO-CLM model v5.0. The meteorological fields were converted to the extension, resolution and projection of the CMAQ nested grids, using an in-house modified version of MCIP. More details on the meteorological forcing data and the evaluation of precipitation can be found in . Here we include an evaluation of T2 and WS10 in the southern part of the Baltic Sea region. Temperature was compared against gridded observational dataset E-OBS v.16 . Wind speed was compared against observational data from MiKlip DecReg of the German Weather Service (DWD). Monthly mean T2 in Denmark and southern Sweden was underestimated in winter (bias smaller than −1.4 K) and overestimated in summer. The warm bias in summer was higher in Sweden (+1.4 K) than in Denmark (+0.4 K). The spatial correlation of T2 in the southern Baltic Sea region based on 3 daily averages was remarkably good. Monthly mean WS10 was slightly overestimated in most parts of the region. The largest errors in wind speed occurred in Denmark and northern Poland during May and June.

EMEP/MSC-W was driven by meteorological data from the Integrated Forecasting System (IFS) of the ECMWF, version IFS38r2, with t1279 resolution (about 0.16 resolution) interpolated to 0.1×0.1. The ECMWF forecasting system of weather parameters is regularly validated by comparing them against European synoptic observation data available on the Global Telecommunication System (GTS). The evaluation of the weather forecast from cycle Cy40r1 is summarized as follows . The frequency of light precipitation is overestimated with a bias of 1.2–1.4 mm d−1 (for precipitation amounts >1 mm d−1). T2 has a negative night-time temperature bias over Europe in winter and early spring. For total cloudiness, bias and standard deviation are small in 2012. For WS10, the standard deviation is low and the night-time bias is very small.

The use of different meteorological datasets introduces additional variability which is, on one hand, wanted to achieve a wider range of possible results for estimating the effect of shipping on air quality but, on the other hand, complicates the interpretation of differences between the models.

### 2.2.3 Boundary conditions

The initial conditions (ICONs) for the simulation and the lateral boundary conditions (BCONs) for the outer European domain are taken from FMI APTA global reanalysis . The global boundary conditions results have been interpolated in time and space to provide hourly boundary conditions for the respective outer domains of the CMAQ and SILAM simulations. The set-up for initial and boundary concentrations for EMEP/MSC-W is described in . ICONs and BCONs are based on long-term observations. For ozone, 3-D fields for the whole domain are specified from climatological ozone sonde datasets, modified monthly against clean-air surface observation. For most other chemical compounds they are defined by simple functions based on measurements and/or model calculations, prescribing concentrations in terms of latitude and time of year or time of day.

### 2.2.4 Anthropogenic land-based emissions

Anthropogenic land-based emissions at hourly resolution obtained from the SMOKE-EU emission inventory were provided for CMAQ and SILAM. These emissions are based on officially reported EMEP emissions which are then distributed in time and space using appropriate surrogates like population density maps, street maps or land use maps. Point sources from the European point source emission register are considered. Vertical distribution of point source emissions is based on real-world stack information and calculated within SMOKE-EU . Dynamic emissions from agricultural activity and animal husbandry depending on meteorological variability are considered . The EMEP/MSC-W model uses anthropogenic emissions from the ECLIPSE gridded emission inventory (http://www.iiasa.ac.at/web/home/research/researchPrograms/air/ECLIPSEv5.html, last access: 17 April 2019). These emissions differ slightly from the reported national total EMEP emissions for 2012; see . For the countries bordering the Baltic Sea (excluding Russia), the national total sulfur emissions from ECLIPSE are about 6 % higher and the NO2 emissions about 10 % lower than the corresponding EMEP emissions.

### 2.2.5 Shipping emissions

The STEAM inventory for the Baltic Sea shipping emissions used in the SHEBA project consists of hourly updated 2×2km2 gridded data for NOx, SOx, carbon monoxide (CO) and particulate matter, which are further divided into elemental carbon, organic carbon, sulfate (${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$) and mineral ash. For the North Sea and other European seas the STEAM data for 2011 were used. Ship emissions were used with hourly time resolution in CMAQ and SILAM, whereas they were used with monthly resolution in EMEP/MSC-W. The use of monthly aggregated ship emissions in EMEP/MSC-W is justified by the fact that the same set of ship emissions from FMI is applied for different meteorological years in the routine application of EMEP modelling and that ship emissions from other seas were only available for 2011. Previous tests with daily and monthly aggregated ship emissions showed that the differences in results are very small. The use of North Sea ship emissions from 2011 on an hourly basis in CMAQ and SILAM causes some inconsistency because meteorological data from 2012 are used in the CTM simulations. Because we are mainly interested in the seasonal variability of pollutant concentrations based on daily averages, the outcome of this study will be less affected by the inconsistency between the timing of ship emissions and the meteorological conditions. STEAM emission data were provided for two vertical layers (below 36 m, from 36–100 m). In CMAQ and SILAM, emissions below 36 m were attributed to the vertical model layers below 42 m, while emissions above 36 m were attributed to the model layers between 42 and 84 m. In EMEP/MSC-W all ship emissions were attributed to the lowest vertical model layer, which typically has a height of 92 m.

## 2.3 Statistical analysis

### 2.3.1 Evaluation method for the total air pollutant concentrations

Model results for total surface concentrations of NO2, O3, SO2 and PM2.5 from the three CTMs are evaluated against available measurements of the air quality monitoring network from the AirBase v8 database (Simoens2014). AirBase is the air quality information system maintained by the European Environmental Agency (EEA) through the European topic centre on Air Pollution and Climate Change Mitigation. Table S1 in the Supplement gives a list of all rural and regional background monitoring stations. Concentrations of NO2 are monitored at 17 stations, O3 at 35, SO2 at 11 and PM2.5 at 8 rural and regional background stations. Table S2 gives all urban and suburban background monitoring stations included in the statistical evaluation of the models. Concentrations of NO2, O3, SO2 and PM2.5 are monitored at 52, 46, 37 and 10 stations of the urban background, respectively. Figure 2a shows locations of stations with NO2 and with O3 measurements. Figure 2b shows locations of the stations with SO2 and with PM2.5 measurements.

Figure 2Map of the Baltic Sea region with the location of background monitoring stations used in the statistical performance analysis with observations of (a) NO2 (filled red circles), O3 (filled green circles), (b) SO2 (filled dark green circles) and PM2.5 (filled yellow circles). Same domain as in Fig. 1.

Download

The model output of surface concentration fields of each CTM is used with its original horizontal resolution to calculate daily mean concentrations. The modelled concentrations are extracted from the respective model grid cell, where the selected monitoring stations are located. The evaluation was done for the entire year of 2012 based on daily means. The model output for PM2.5 was taken from the modelled PM2.5 containing aerosol water at 50 % relative humidity.

The performance of each model is quantified in terms of mean values (μMod and μObs), normalized mean bias (NMB), Spearman's correlation coefficient (R), root mean square error of the modelled values (RMSE) and fraction of model values within a factor of 2 of the observations (FAC2). Definitions of NMB, R, RMSE and FAC2 are given in Appendix A. The model performance analysis is discussed separately for rural background stations and urban background stations. In order to better highlight model differences in terms of urban areas and station types (i.e. rural, suburban, urban background sites), groups of stations (rural vs. urban) are generated in which statistical performance indicators are averaged. In the rural group, rural background and regional background stations are included, while in the urban group, urban background and suburban background stations are included. Monitoring stations classified as traffic stations and industrial stations were not included in the comparison, since the regional CTM systems applied here do not handle the local-scale dispersion near emission sources.

In the context of this evaluation of predicted air pollutant concentrations, we consider a correlation coefficient of more than 0.5 to indicate a correlation between modelled and observed time series, while values of 0.7 and above are considered to be good correlations. define certain acceptance criteria for model performance based on their experience in conducting a large number of model evaluation exercises. For rural stations, FAC2 values >0.5 and for urban stations FAC2 values >0.3 indicate acceptable performance. We adopt these bounds in the present study to characterize the predictive strength of the models with respect to the pollutant concentrations.

We compare the performance between models with the help of a graphical comparison in the form of box plots. Box plots of the correlation coefficient, NMB and RMSE, including either all rural or all urban monitoring stations, were prepared. The box plots show the median as a line dividing the box into two parts, the upper and lower quartiles as the ends of the box and the minimum and maximum values of the data and outliers.

In addition to the model performance for temporal correlation we evaluated the spatial correlation of the total air pollutant concentrations with the observations of the AirBase network for the three CTMs.

### 2.3.2 Significance of the ship contribution

The method described in was used to assess the significance of the ship influence on ambient NO2 at the monitoring stations. The ship influence at a station was positively confirmed in the tests if (1) the concentrations increased and (2) the temporal correlation improved when shipping emissions were included in the CTM simulation.

By means of a paired t test it was first tested whether the modelled NO2 concentrations at the monitoring stations with available NO2 observations (Table S1) significantly increased if shipping emissions were considered. This test estimated whether the mean concentration difference between the noship run and the base run (noship – base) is significantly equal to or greater than zero, indicated by the probability pbias. If the value of pbias was greater than the level of significance of 0.05, then this hypothesis was confirmed, which means that no difference between the base run and the noship run could be statistically proven. In case pbias was less than 0.05, it was decided that the model run without shipping emissions led to lower concentrations, confirming the ship influence.

The significance of the improvement in the correlation between simulations and observations was tested by calculating the Fisher z transformation of the two correlation coefficients for the two model runs (noship and base) and testing the hypothesis “greater than”. Correlation coefficients were calculated with Spearman's method for consistency with the statistical evaluation. The probability pcorr was calculated for the hypothesis that the correlation between the base run and observations is greater than the correlation between the noship run and observations. We accepted this hypothesis if the probability was higher than 0.9. Therefore, in the following, a station i with ${p}_{\text{corr},i}>\mathrm{0.9}$ for a specific CTM simulation is termed “ship influenced”.

3 Results
Back to toptop

## 3.1 Statistical evaluation of air pollutant concentrations

### 3.1.1 Rural vs. urban sites

A statistical performance analysis for each of the three CTMs was undertaken using the available observation data form AirBase for 2012 based on daily mean total concentrations. The evaluation results for the temporal variation of air pollutants are summarized in Table S3 for daily mean NO2, in Table S4 for daily mean O3, in Table S5 for daily mean SO2 and in Table S6 for daily mean PM2.5. In the following, the performances of the models that simulate air pollutant concentrations are compared and discussed separately for the group of rural stations and for the group of urban stations in order to highlight differences in the predictive capability of the models for rural vs. urban sites.

Figure 3Comparison of statistical indicators for NO2 daily means (in the order R, NMB and RMSE) between three CTMs at (a) rural background stations and (b) urban background stations. Outliers shown as small circles.

Download

Figure 4Comparison of statistical indicators for O3 daily means (in the order R, NMB and RMSE) between three CTMs at (a) rural background stations and (b) urban background stations. Outliers shown as small circles.

Download

The atmospheric lifetime of NO2 is relatively short: a few hours in summer and up to 1 day in winter ; hence differences between rural and urban sites are expected due to the higher emission density in urban or industrial areas. The urban station average of observed annual mean NO2 is more than twice the concentration average at rural sites. The three CTMs underestimate the annual means at the urban sites. The overall correlation of NO2 for rural stations is good for all models. The overall correlation of NO2 for urban sites is lower than at the rural sites (Fig. 3). At most urban stations, models underestimate the observed daily mean NO2 by ca. 40 % (Table S3). The general underestimation of NO2 at urban sites has been evident in other multi-model air quality studies in Europe (e.g. Giordano et al.2015). The finer horizontal resolution of CMAQ and SILAM (4 km) compared to EMEP/MSC-W (11 km) does not result in a significant improvement in the urban bias and urban temporal correlation. This result was expected based on the study by , who found no further improvement in the urban signal, i.e. the concentration difference between high-emission areas and their surroundings, when increasing the resolution from 14 to 7 km in the same model. Moreover, increasing the spatial resolution in the model does not help to significantly improve the performance in time because the temporal variability of pollutants is affected by the meteorological conditions and pollution levels upwind .

Figure 5Comparison of statistical indicators for SO2 daily means (in the order R, NMB and RMSE) between three CTMs at (a) rural background stations and (b) urban background stations. Outliers shown as small circles.

Download

Tropospheric ozone is largely controlled by the atmospheric transport from regions outside the study area, by stratosphere–troposphere exchange and by the photochemical production through the oxidation of VOCs and carbon monoxide (CO) in the presence of NOx and sunlight. The higher density of NOx emissions in urban areas is expected to lead to a larger titration effect of NO on ozone, which results in lower average O3 at the urban sites compared to rural sites. The models slightly overestimate the O3 measurements at urban sites, with CMAQ having the smallest bias. CMAQ and SILAM predict similar annual mean concentrations as observed for both rural and urban sites, whereas EMEP/MSC-W predicts higher annual mean ozone (Fig. 4). The ozone bias might be linked to boundary conditions : the EMEP model uses ozone boundary conditions from long-term observations, whereas CMAQ and SILAM models use boundary conditions from the FMI APTA global reanalysis. The models slightly overestimate the O3 measurements at urban sites, with CMAQ having the smallest bias (average NMB =0.08). The average RMSE values for the rural sites and the urban sites are similar for the three models (Fig. 4), indicating comparable model performance for the CTMs with respect to daily mean O3 concentrations.

Figure 6Comparison of statistical indicators for PM2.5 daily means (in the order R, NMB and RMSE) between three CTMs at (a) rural background stations and (b) urban background stations. Outliers shown as small circles.

Download

Another major air pollutant is SO2. It is primarily emitted from anthropogenic emission sources such as coal power plants, residential heating, waste incineration and shipping activities. SO2 acts as a precursor to sulfates, which are one of the main components of particulate matter in the atmosphere. The atmospheric lifetime of SO2 is on the order of a few days . SO2 can still be considered to be relatively short-lived and thus less influenced by transport from regions outside the study area. Most emission sources of SO2 are located in urban areas. In the case of power plants, the emissions of SO2 are, however, injected at elevated height and therefore do not directly impact the surface concentrations in the urban area. The urban station average of observed annual mean SO2 is 3 times higher than at the rural stations. At rural and urban sites, the modelled daily mean SO2 from CMAQ and SILAM has a positive bias, whereas the modelled daily mean SO2 from EMEP/MSC-W has a slight negative bias (Fig. 5). For EMEP/MSC-W some urban stations have a FAC2 below 0.3 due to the underestimation of observed SO2. At urban stations, the temporal correlation between model data and observed SO2 shows a mixed performance by the models, with good correlation at some stations and poor correlation at others. The model data and the observations for daily mean SO2 at the rural sites are correlated, but with only 10 stations, the rural station group for SO2 is rather small, limiting the conclusions that can be drawn from the statistical analysis. The weaker performance of the models for SO2 at the rural sites is related to uncertainties in local residential heating emissions, as the timing of use and the sulfur content of burned fuels are difficult to predict.

Ambient PM2.5 is a widespread pollutant, which is directly emitted by biomass and fossil fuel combustion in domestic and industrial activities, and it is also formed from gaseous precursors such as NOx, SO2, NH3 and NMVOC in the atmosphere. The atmospheric lifetime of PM2.5 is on the order of days or weeks and thus PM2.5 can be subject to long-range transport. For PM2.5, smaller differences between rural and urban stations are expected than for NO2 and SO2 because PM2.5 has a large secondary component, which is generally more homogeneously distributed over rural and urban areas. SILAM is able to reproduce annual mean PM2.5 concentrations for urban stations, whereas CMAQ and EMEP/MSC-W give lower annual mean values than observed. For urban stations, the temporal correlation of daily mean PM2.5 is good for all models, whereas for rural stations the temporal correlation is slightly better for EMEP/MSC-W than for CMAQ and SILAM (Fig. 6). At both rural and urban sites, the modelled daily mean PM2.5 from CMAQ and EMEP/MSC-W has a slightly negative bias, whereas modelled daily mean PM2.5 from SILAM has no bias.

### 3.1.2 Spatial correlation

The spatial correlation between modelled and observed annual mean total pollutant concentrations for the three CTMs is presented in Fig. 7. Because NOx is mainly emitted near the ground, the spatial distribution of NO2 is expected to be highly correlated with emissions. The improvement in spatial resolution of emissions should lead to improved spatial correlation between modelled and observed concentrations. In contrast with this expectation, EMEP/MSC-W shows the best correlation and the lowest bias. Observed annual mean NO2 at urban stations is strongly underestimated by SILAM and in CMAQ. The positive bias indicates that observed NO2 at rural stations tends to be overestimated by these two models. The annual mean O3 is closely linked to the annual mean NO2 through the local titration effect. Hence, lower than observed O3 at rural stations for SILAM and CMAQ is related to the higher than observed NO2. The models are capable of representing the seasonal variation of ozone, with highest average concentrations in spring, followed by summer, and lowest concentrations in winter and autumn.

Figure 7Spatial correlation of annual mean concentrations (µg m−3) and the seasonal averages from CMAQ (left column), SILAM (middle column) and EMEP (right column) in the Baltic Sea region for (a) O3, (b) NO2, (c) SO2 and (d) PM2.5.

Download

The spatial correlation between modelled and observed SO2 is weaker than that for NO2, which is probably due to the fact that most SO2 sources are emitting into higher vertical layers. CMAQ and SILAM overestimate SO2 in autumn and winter at many stations by a factor of 2 or more, which is likely related to the uncertainty in the residential heating emissions. EMEP/MSC-W underestimates SO2 in summer, which might be connected to uncertainties in the vertical emission distribution. EMEP-MSC/W shows the best spatial correlation for PM2.5 with almost no bias. However, annual average PM2.5 is underestimated by 23 % on average. For CMAQ and SILAM the spatial correlation for annual mean PM2.5 has a positive bias due to overestimation at the rural stations. At almost all stations, CMAQ underestimates PM2.5 in summer. This has also been evident for regional background stations of the EMEP monitoring network and can partly be attributed to the underestimation of secondary organic aerosols in the CMAQ simulation . SILAM underestimates PM2.5 in winter. Since PM2.5 in winter is mainly from anthropogenic sources and SILAM uses the same emissions as CMAQ, the discrepancy in winter PM2.5 is attributed to problems with simulating stagnant meteorological conditions.

## 3.2 Evaluation of ship-related concentration contributions

A direct comparison of the modelled ship-related concentration contribution to measurements of the shipping signal (in exceedance of the background air) is hampered by the fact that measured concentration increases due to individual ship plumes do not reflect the entire contribution of shipping at sea. In order to evaluate the modelled ship contributions, a statistical method (Sect. 2.3.2) was applied to decide whether the modelled concentration as well as the correlation with observed daily mean NO2 concentration at a specific station increases significantly when ship emissions are included in the CTM simulation. The results of the significance test are summarized in Table S7.

A significant concentration increase was found at all 69 stations for the three CTMs. However, the significance of the concentration increase only shows that the modelled concentrations at a station are sensitive to ship emissions. The correlation increases significantly (on 0.9 or 0.95 level) at 10, 7 and 8 stations for CMAQ, SILAM and EMEP/MSC-W, respectively (Table S7).

Four ship-influenced stations were identified by all models: Vilsandi (EE0011R, Estonia), Utö (FI00349, Finland) at the shoreline, Norr Malma (SE0066A), 20 km inland at the east coast of Sweden and Lübeck-St Jürgen, (DESH023, Germany) close to a port (time series plots in Appendix B). At Norr Malma, NO2 is largely overestimated by the models in summer when ship emissions are included. We suggest that, due to the high spatial resolution of the model, ship exhaust plumes are resolved but not adequately dispersed in the models. Since the models do not specifically treat the plume dispersion of individual ships, the spreading of the plume might not be sufficiently large or the plume height from ship exhausts is not properly considered with the applied vertical profile of ship emissions.

Ship-influenced stations found by any of the models included shoreline stations (Virolahti, FI00351, Finland; Lahemaa, EE0009R, Estonia), stations in harbour cities (Rostock Warnemünde, DEMV021, Germany; Kiel, DESH033, Germany; Aalborg/8158, DK0053A, Denmark, Södermalm, SE0022A, Sweden; Gdansk Pm.a09aN, PL0053A; Poland) and one urban inland station (Pm.63.wDSAa, PL0171A, Poland). The corresponding time series plots of daily mean NO2 are shown in Appendix B. Observed daily mean NO2 at the two urban stations in Poland is underestimated by all models, indicating missing local emissions from other sectors. The ship influence at station Rostock Warnemünde, located close to a harbour, was significant in EMEP/MSC-W but not in the other two models (Fig. B1i). This could indicate that differences in the meteorological data, in particular wind flow fields, are responsible for the different ship influences. Although, the timing and location of ship exhaust plumes – based on AIS data – should be accurate during the port stays, the emission fluxes at berth are more challenging to estimate, because this involves an estimation of electrical power usage during the port stays. Evaluation of the ship contribution in Rostock using an urban air quality model with a high degree of detail on ship emissions and other urban emissions showed that shipping significantly impacts on annual averaged NO2 in the city domain .

For all ship-influenced stations, time series plots of daily mean O3 are compiled in Fig. S1 in the Supplement. Including the shipping emissions in the model simulations affected ozone concentrations mainly in the summer months. The change in modelled O3 due to shipping was below 6 % on average in summer at the ship-influenced stations. At most of the stations, including ship emissions increased the modelled O3 concentration as a consequence of photochemical ozone production. The seasonal variation of NO2 at the ship-influenced sites with peak concentrations in winter is in general reproduced by the models. Including ship emissions improved the agreement between modelled and measured NO2 daily mean concentrations at about 50 % of the stations. For more than 70 % of all stations the observations of total NO2 concentrations are reproduced by the models within a NMB range from −0.5 to 0.5 in the base simulation (Fig. S2).

## 3.3 Comparison of the spatial distribution of air quality indicators

### 3.3.1 Spatial distribution of annual mean NO2

A strong south–north gradient for annual mean NO2 concentrations is found for the Baltic Sea region in the three base simulations, with 4–5 times higher NO2 concentrations in the southwestern part than in the northern part of the region (Fig. 8a). High modelled NO2 concentrations are predicted in Denmark, northern Germany and Poland as well as over the Danish Straits and in the urbanized areas of the region. Modelled annual mean NO2 concentrations in the proximity of the main shipping routes several times exceed the concentrations in the regional background. EMEP/MSC-W shows the strongest concentration gradients between urban and rural areas and between shipping lanes and the surrounding sea. The simulations with the other two models result in a wider spread of the NOx emissions from the shipping routes and the urban centres, indicating stronger horizontal transport by advection and diffusion in CMAQ and SILAM. This finding is counterintuitive as the NOx emissions are initially less diluted than in the EMEP/MSC-W simulation because of the smaller volume of the grid boxes and therefore result in higher NOx concentrations near the emission sources. Atmospheric transport by diffusion processes are subgrid mixing processes, which are not resolved by the given resolution of the applied models. For large grid cells, e.g. 50×50km2, the numerical diffusion will usually be much larger than the physical diffusion in the horizontal direction. However, at finer-resolution scales, the physical diffusion will gradually become more important than numerical diffusion and becomes greater than numerical diffusion for a cell size of 5×5km2 or below . The wider spread of elevated NO2 concentrations is also indicative of a longer atmospheric lifetime of NO2 in CMAQ and SILAM compared to the simulation with EMEP/MSC-W. NO2 is removed relatively quickly in the lower troposphere through the reaction with hydroxyl (OH) radicals to form HNO3. The rate coefficient for this reaction, k(NO2+OH), is similar in the three models ((1.1–1.2)$×{\mathrm{10}}^{-\mathrm{11}}$cm3 s−1 at 298 K). Thus, differences in the NO2 lifetime are mainly due to different abundances of OH radicals in the simulations. High NO2 concentrations in Belarus and Russia in the SILAM simulation are an artefact from merging with the output of the coarser central model grid (Sect. 2.2.1).

Table 2Spatial averages of the annual mean concentrations of NO2, O3, SO2, PM2.5 and EC in µg m−3 for the study domain (Baltic Sea region as in Fig. 1).

Figure 8Comparison of the spatial distribution of annual mean concentrations (µg m−3) from CMAQ (left column), SILAM (middle column) and EMEP (right column) in the Baltic Sea region for (a) O3, (b)  NO2, (c) SO2 and (d) PM2.5. Empty areas correspond to concentrations between zero and the lowest value in the legend.

Download

### 3.3.2 Spatial distribution of annual mean O3

Modelled annual mean O3 concentrations over the Baltic Sea are 15 %–25 % higher than over land but are reduced along the shipping lanes due to the titration effect caused by the ship-emitted NOx (Fig. 8b). Lowest ozone concentrations are seen for St Petersburg (< 32 µg m−3) in the three simulations, although we note that the city is outside of the high-resolution grid in the case of SILAM. The spatial average of annual mean O3 is clearly higher for the EMEP/MSC-W simulation, by ca. 20 % compared to the other two simulations (Table 2). The most probable reason for the difference is the application of different sets of boundary conditions for the European model domains, as discussed in Sect. 3.1.1. Model simulations for Europe have shown a high sensitivity of ozone changes to the dry deposition to vegetation . Thus differences in the deposition schemes may partly explain the different O3 levels over the continent, e.g. when comparing ozone over Sweden and Finland between CMAQ and SILAM.

### 3.3.3 Spatial distribution of annual mean SO2

Clear differences in the spatial distribution of the annual mean SO2 concentrations are found between CMAQ and SILAM on one hand and EMEP/MSC-W on the other hand (Fig. 8c). The simulation with CMAQ and SILAM shows a southeast–northwest gradient with elevated SO2 over large parts of the southern Baltic Sea region, Poland, Belarus, Russia and the Baltic States with annual mean concentrations in the range of 1.3–3.0 µg m−3. Residential heating emissions and power plant emissions for district heating in the urban centres and rural areas strongly contribute to the high SO2 concentrations in this subregion. In the EMEP/MSC-W simulation, elevated SO2 concentrations are present along the main shipping routes in urban areas and in Poland, whereas the levels of SO2 outside of these areas are much lower. The concentration gradients between urban and rural areas and between shipping lanes and surrounding sea are up to 2.5 µg m−3 for EMEP/MSC-W, while it is only up to 0.7 µg m−3 for CMAQ and SILAM. Factors contributing to the different gradients are differences in the representation of horizontal transport (see Sect. 3.3.1), differences in the meteorological conditions and differences in the atmospheric lifetime of SO2. The atmospheric lifetime of SO2 is determined by its reaction with the OH radical and by its removal via dry deposition. In EMEP/MSC-W, the canopy uptake of SO2 is strongly controlled by NH3 levels, and the implemented deposition parameterization accounts for co-deposition effects on the dry deposition of SO2 . Co-deposition effects are not considered in the other two models.

### 3.3.4 Spatial distribution of annual mean PM2.5

Modelled annual mean PM2.5 is higher in the southern part, both over land and sea, than in the northern part of the Baltic Sea region (Fig. 8d). On annual average, PM2.5 concentrations are not elevated along the shipping routes. The seasonal differences between summer and winter will be discussed below (Sect. 3.5) and will help to understand differences between the models. High PM2.5 levels (8 µg m−3 and higher) are simulated in the urban areas of major cities like Copenhagen, Oslo, Helsinki, Riga, Tallinn and St Petersburg. The high PM2.5 levels over the continent in the southern part of the Baltic Sea region presumably result from a combination of land-based primary emissions, long-range transported particles and the formation of secondary particulate matter.

Table 3Relative ship contribution to the spatial average of annual mean NO2, O3, SO2, PM2.5 and EC in percent for the study domain (Baltic Sea region as in Fig. 1). Values in brackets denote the average relative ship contribution in the coastal land areas of the domain.

### 3.3.5 Recommendations from the comparison between the CTM systems

The applied CTM modelling systems originate from different institutes and represent independent lines of development. Their operations require varying degrees of the user experience, input data requirements and computational demand. All three systems are open source, installed and used in a number of countries and possess long records of operational and research applications. The EMEP and SILAM models are usually less demanding than CMAQ in terms of computational resources and input data. Yearly totals of anthropogenic emissions can be used as input to the two models, which perform the temporal disaggregation of emissions in line with the computation. CMAQ probably has the most-extensive user community with support provided by the developers from the US Environmental Protection Agency .

Comparison of the model performances does not give an unequivocal answer: all model skills are within the uncertainty in the corresponding parameters (Figs. 36). One should, however, bear in mind that the EMEP model was run with a lower resolution than two other models. There are certain systematic differences between model results from SILAM and CMAQ on one hand and EMEP on the other: e.g. higher NO2 concentrations at rural stations and lower at urban stations (compare the spatial correlation of annual station averages in Fig. 7). To a large part, this mismatch can be attributed to the spatial distribution of anthropogenic emissions in the SMOKE inventory, applied in the two CTM systems in comparison with ECLIPSE emissions used by EMEP.

The EMEP/MSC-W model is routinely used for multi-year calculations, facilitating its use for the HELCOM (Baltic Marine Environment Protection Commission – Helsinki Commission) evaluation of trends in the deposition of nitrogen and sulfur in the Baltic Sea region. CMAQ is being used for a variety of environmental modelling problems including regulatory applications and evaluation of emission control strategies . CMAQ is continuously updated to remain a state-of-the-art regional CTM. A specific advantage of SILAM is the online computation of wildfire emissions and operational input of hourly STEAM ship emissions. The evaluation of modelled daily mean PM2.5 showed that RMSE station values for SILAM are within a smaller range (between lower and upper quartiles) than the other two models (Fig. 6). The model was also recently applied to 35-year-long global-to-local reanalysis of air quality by . SILAM can therefore be recommended for use in advanced research applications, specifically addressing the abundance and composition of particulate matter.

## 3.4 Comparison of the ship contribution in the three CTMs

The influence of shipping emissions on the air quality was evaluated for the annual mean concentrations of the three CTMs. The results for the impact of shipping emissions were calculated as differences between the base and the noship simulations. Results for the absolute ship-related concentrations of O3, NO2, SO2 and PM2.5 are shown in Fig. 9, the resulting relative ship contribution to annual mean concentrations is shown in Fig. S3 and the spatial average of the relative ship contribution is given in Table 3.

### 3.4.1 Ship contribution to annual mean NO2

The ship-related annual mean NO2 concentrations from the three CTMs are in the range of 3–5 µg m−3 along the main shipping routes. The NO2 ship contribution decreases to Baltic Sea background values (about 1 µg m−3) within a few hundred kilometres from the centre of the shipping routes. Ships emit NOx mainly in the form of NO, which is, however, quickly converted to NO2. Thus atmospheric NOx is mainly in the form of NO2. The relative contribution of ship emissions to annual mean NO2 is more than 40 % over the Baltic Sea (Fig. S3), 22 %–28 % for the entire Baltic Sea region and 16 %–20 % in the coastal land areas (Table 3). In particular, NOx emissions from shipping affect the harbour cities of the region and coastal areas in southern Sweden. Local differences between the models might be due to the different meteorological drivers or differences in the titration efficiency for ozone.

Figure 9Comparison of the spatial distribution of annual mean ship-related concentrations (absolute ship contributions in µg m−3) of the CMAQ (left column), SILAM (middle column) and EMEP (right column) models in the Baltic Sea region for (a) O3, (b) NO2, (c) SO2 and (d) PM2.5.

Download

### 3.4.2 Ship contribution to annual mean O3

In the proximity of the main shipping routes, negative concentration differences for the modelled annual mean O3 between the base and the noship simulation are obtained as a result of the titration effect by the NOx emissions from shipping. The highest ozone reduction due to shipping is found in the western part of the Baltic Sea. In the CMAQ simulation the depletion of ozone is stronger than in the other two models, with O3 reduction of 6–12 µg m−3 in the Kattegat and in the Danish Straits. The hourly variation of ship emissions is represented in the simulations with CMAQ and SILAM, whereas monthly averaged ship emissions are used in the EMEP/MSC-W simulation. Emission peaks of NOx from ships that are present in the hourly data can result in occasional stronger ozone titration, leading to overall higher reduction of ozone than is the case for monthly averaged ship emissions. Over the coastal land areas, the average impact on annual mean of O3 is very small, with ozone increases between 0.1 % and 0.4 % for the models.

### 3.4.3 Ship contribution to annual mean SO2

Ship emissions of SO2 have a high contribution to annual mean SO2 concentrations over the Baltic Sea. The ship contribution to SO2 is 0.5–0.7 µg m−3 in a wide corridor around the main shipping routes of the Baltic Sea. While the absolute ship contribution of the three CTMs is similar, the relative ship contribution in the EMEP/MSC-W simulation is higher in most areas of the Baltic Sea and in Sweden, because the background atmospheric SO2 levels in this simulation are lower than in CMAQ and SILAM.

### 3.4.4 Ship contribution to annual mean PM2.5

The ship contribution to annual mean PM2.5 shows a gradient from southwest to north with highest concentrations over Denmark, the west coast of Sweden, the Belt Sea/Kattegat and over the sea south of Sweden, with maximum values up to 0.9 µg m−3. The relative contribution in these ship-impacted areas is up to 10 %. The average ship contribution for the three CTMs is in the range of 4.1 %–6.5 % in the entire Baltic Sea region and 3.1 %–5.7 % in the coastal land areas. The absolute ship contribution in SILAM is slightly smaller than for the other two models, in particular in the southwestern part of the Baltic Sea region (Fig. 9d). The ship-related PM2.5 affects the coastal areas in the Baltic Sea region, as its influence extends further inland than is the case for ship-related NO2 or SO2. This can be attributed to the formation of secondary particulate matter in the ship exhaust plume during its transport away from the main shipping routes.

## 3.5 Comparison of PM2.5 in summer and autumn

CMAQ and EMEP/MSC-W simulations predict higher seasonal mean concentrations of PM2.5 in autumn (average of September, October and November, SON) than in summer (average of June, July and August, JJA), whereas the SILAM simulation predicts higher PM2.5 in summer (Fig. 10a, c; Table S10). The temporal correlation between model data and observations of daily mean PM2.5 for the average of the AirBase stations is slightly better in autumn than in summer (Tables S8 and S9). Observed PM2.5 in summer is underestimated strongly by CMAQ (at almost all stations; see Sect. 3.1.2) and to some extent by EMEP/MSC-W. In autumn all models are in better agreement with observed PM2.5.

The SOA formation mechanism in the applied version of CMAQ (i.e. v5.0.1) is probably not adequate for reproducing the summertime aerosol. Primary organic aerosol (POA), SOA and organic vapours in the atmosphere should be considered a dynamic system that constantly evolves due to multigeneration oxidation . We note that multigenerational ageing chemistry for the semi-volatile POA was introduced in CMAQ v5.2, based on the approach by , which considers the functionalization and fragmentation of organic vapours upon oxidation. In addition, wildfire emissions have not been considered in the simulation with CMAQ. Wildfires emit large quantities of organic material and are associated with high biogenic VOC emissions due to high temperatures, leading to increased SOA formation .

In summer, modelled mean PM2.5 in the region is much higher in the SILAM simulation (5–8 µg m−3 in most parts, 7.4 µg m−3 on regional average, 5.4 µg m−3 on average in coastal land areas) than for the other two models (<4µg m−3, except for the urban areas). The higher summertime PM2.5 in SILAM is most likely due to more efficient SOA production and primary emission from wildfires and/or mineral dust. A previous comparison of the models to PM2.5 observations from the EMEP station network in Europe reported similar seasonal mean concentrations of the SIA components, i.e. nitrate (${\mathrm{NO}}_{\mathrm{3}}^{-}$), ammonium (${\mathrm{NH}}_{\mathrm{4}}^{+}$) and ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, for the three CTMs in summer .

Figure 10Comparison of PM2.5 in summer and autumn from CMAQ (left column), SILAM (middle column) and EMEP (right column) in the Baltic Sea region for (a) JJA mean concentration (µg m−3), (b) JJA mean ship contribution (µg m−3), (c) SON mean concentration and (d) SON mean ship contribution. Empty areas correspond to concentrations between zero and the lowest value in the legend.

Download

Table 4Spatial averages of seasonal mean concentrations of EC (µg m−3) from the base simulation and ship contributions to EC levels (µg m−3) for the study domain (Baltic Sea region as in Fig. 1). Mean values are given for spring (March to May; MAM), summer (JJA), autumn (September to November; SON) and winter (January, February and December 2012; DJF). Values in brackets denote the seasonal mean concentrations in the coastal land areas of the domain.

The calculated ship contribution from all models is higher in summer than in autumn (Table S10). The simulations reflect the greater importance of shipping activities during summer and their influence on PM2.5 levels over the entire Baltic Sea and the coastal areas (Fig. 10b). In particular Denmark and the Swedish west coast is highly impacted in summer, with a ship contribution of 0.5–0.9 µg m−3 to ambient PM2.5 levels.

In autumn, all CTMs predict high levels of PM2.5 in the southern part of the Baltic Sea region, exceeding 6 µg m−3. High PM2.5 in autumn is typically attributed to stagnant meteorological conditions and higher emissions of primary particulate matter from residential heating and energy production. Modelled PM2.5 in Sweden and Finland is higher in SILAM than in the other two models. SILAM overestimates observed PM2.5 at the stations in Sweden, Lithuania and Finland in summer (NMB: 0.54 on average) and autumn (NMB: 0.31 on average). In an earlier model comparison, all three models were shown to overestimate ${\mathrm{NO}}_{\mathrm{3}}^{-}$ and ${\mathrm{NH}}_{\mathrm{4}}^{+}$ in autumn, while SILAM also overestimated ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$ in autumn .

The ship contribution in autumn in the southwestern part of the region is higher in EMEP/MSC-W compared to the other two models (Fig. 10d), which is obviously a result of larger secondary formation of particulate matter, as mainly the coastal regions are impacted. The formation of SIA in autumn is favoured by lower temperature and higher humidity compared to summer. The higher autumn ship contribution in the EMEP model can be attributed to differences in land-based anthropogenic emissions of NH3 and NO2 or to differences in the schemes for inorganic aerosol formation. The investigation into differences between the SIA formation schemes is, however, beyond the scope of this study.

## 3.6 Comparison of elemental carbon related to ship emissions

Primary carbonaceous particles emitted from ships are the product of incomplete fuel combustion and consist of a mixture of elemental carbon and non-polar organic carbon. In the STEAM ship emission inventory these are separated into emissions of elemental carbon (EC) and organic carbon (OC). The terms EC and BC are used interchangeably in the models; however both can only be regarded as proxies for the concentration of soot particles . Here we are mainly interested in the atmospheric fate of EC from ship emissions, as simulated by the models. EC particles are associated with adverse human health effects and contribute to regional haze and poor visibility (e.g. Odman et al.2007). The atmospheric lifetime of EC is relatively long, around 6 d in the continental outflow and 4–8 d on a global scale , with a large uncertainty due to soot ageing processes and wet deposition .

The spatial averages of the EC concentrations (base simulation) and the ship-contributed EC concentrations are given in Table 4. The seasonality of ship-related EC predicted by the three CTMs is shown in Fig. 11. The levels of ship-related EC are higher in spring and summer than in autumn and winter due to more intense shipping activities. Therefore the ship contribution peaks in the seasons when ambient EC concentrations are lowest (Table 4; Fig. S4). The highest levels of ship-related EC, in the range of 0.03–0.04 µg m−3, occur along the main shipping routes and in the main ports of the region. The average ship contribution to annual mean EC is 4 %–5 % over coastal land regions (Table 3).

Measurements of the ship contribution to equivalent black carbon (eBC) concentrations at a shoreline location in southern Sweden (Falsterbo [55.3843 N, 12.8164 E] downwind of main shipping lanes, based on 113 individual plumes, reported a value of 0.0035 µg m−3 as an average of the winter campaign in 2016 . Wintertime average modelled ship-related EC at this location is factor of 4 to 6 higher than the measured value (CMAQ: 0.0207 µg m−3; SILAM: 0.0144 µg m−3, EMEP/MSC-W: 0.0149 µg m−3). The discrepancy might arise from comparison with a different year than used in the model simulations. Another reason for the higher model values is that the CTMs consider all ships within a radius of 50 km upwind, whereas measurements considered individual ships passing by in a limited sea area.

SILAM predicts a stronger seasonal variability of the ship-related EC than the other models. In particular, modelled EC ship contribution in winter is lower than for CMAQ and EMEP/MSC-W. Shipping emissions of EC are identical in the three CTMs on a monthly basis. Differences between the models are therefore explained by differences in the meteorological conditions, deposition schemes and the treatment of atmospheric transport in the models.

Figure 11Spatial distribution of the seasonal mean EC ship contribution (µg m−3) from CMAQ (left column), SILAM (middle column) and EMEP (right column) in the Baltic Sea region for (a) spring, (b) summer, (c) autumn, and (d) winter. Empty areas correspond to concentrations between zero and the lowest value in the legend.

Download

In stable conditions, the boundary layer (BL) height over the Baltic Sea is often at or below 500 m (). Climatological simulations over the Baltic Sea show that there is a strong seasonality in the atmospheric stability over the sea with more than 50 % stable conditions in spring, whereas during the other seasons unstable conditions dominate together with occasionally neutral conditions . The lower surface concentrations of ship-related EC in SILAM and EMEP/MSC-W during spring compared to CMAQ might be attributed to different atmospheric stability in the simulations, i.e. more frequent occurrence of neutral conditions, which dilutes the concentration of the emitted pollutant. The fact that ship-related EC in EMEP/MSC-W is more confined to the shipping routes and shows a limited spatial spreading compared to the other models might also indicate less efficient horizontal diffusion. In winter, SILAM and EMEP/MSC-W ship-related EC is very low and the impacted area has a smaller extent than in the other months (white areas in Fig. 11), indicating faster removal of EC particles than in the CMAQ simulation. Different treatment of the hygroscopicity and ageing processes of EC particles, affecting their wet scavenging, could have contributed to the differences among the models.

4 Summary and conclusions
Back to toptop

The effect of ship emissions on the regional air quality in the Baltic Sea region was investigated with three regional CTM systems (CMAQ, EMEP/MSC-W and SILAM) that simulate the transport, chemical transformation and fate of atmospheric pollutants. The models were applied with their typical set-up for air quality studies in European regions. The same ship emission dataset from the STEAM model based on ship movements from AIS records, detailed ship characteristics and up-to-date load-dependent emission factors were used in all CTMs. The models were set up with a finer grid resolution (4 to 11 km grid length) than was the case for previous air quality studies in the Baltic Sea region, potentially enabling better treatment of the dispersion and photochemistry in exhaust plumes from shipping along the major ship lanes of the Baltic Sea.

The comparison of total concentrations of regulatory air pollutants among the models is the primary focus of this study. Results from the three CTMs were compared to observations from rural and urban background stations of the AirBase monitoring network in the coastal areas of the Baltic Sea region. The finer resolution of CMAQ and SILAM (4 km) compared to EMEP/MSC-W (11 km) did not lead to a significant improvement in the urban bias and urban temporal correlation for daily mean NO2 concentrations. The benefit from using high-resolution grids depends on the availability of accurate urban emission data with high spatial resolution and realistic temporal profiles . While the STEAM inventory provides this data for shipping, the compilation of urban emission inventories is more challenging because they are based on specific information for each sector, such as housing units for domestic heating or number of vehicles . Observed PM2.5 in summer is underestimated strongly by CMAQ at all stations and to some extent by EMEP/MSC-W. Observed PM2.5 in winter is underestimated by SILAM. In autumn all models are in better agreement with observed PM2.5. The low summer PM2.5 in CMAQ has been attributed to the underestimation of secondary organic aerosols and to the missing emissions of wind-blown soil dust particles . Particulate matter emissions from wind-blown dust and forest fires were included in EMEP/MSC-W and SILAM. Ship-related concentrations of NO2 were evaluated at coastal monitoring stations by testing if the agreement between predicted and observed total concentrations improves significantly when ship emissions are included in the simulation. Including ship emissions improved the agreement between modelled and measured NO2 daily mean concentrations at about 50 % of the stations. The change of modelled O3 due to shipping was below 6 % on summer average at the shoreline stations and mainly reflects additional ozone production due to ship emissions. Ship-influenced stations identified by the models are mainly located close to the shoreline or close to a port, with a maximum distance of 20 km from the sea. However, modelled peaks of high daily mean NO2 at coastal rural sites during summer that are not present in the measurements indicate that the models often did not properly resolve the ship plumes due to the subgrid variability of the plume dispersion of individual ships.

The spatial average of annual mean O3 concentrations in the EMEP/MSC-W simulation is ca. 20 % higher compared to the other two simulations. EMEP/MSC-W overestimates the measurements of daily mean O3 concentrations at rural stations by 17 % on average. The higher ozone concentrations in the EMEP model are mainly the consequence of using a different set of boundary conditions for the European model domain. The concentration gradients of NO2 and SO2 between urban and rural areas and between shipping lanes and the surrounding sea are larger in EMEP/MSC-W than in the other models. Factors contributing to the different gradients are differences in the representation of horizontal transport, differences in the meteorological driving data and differences in the atmospheric lifetime of NO2 and SO2 in the models.

There are significant differences in the calculated ship contributions to the levels of air pollutants among the three models. In the proximity of the main shipping routes, ozone is depleted as a result of the titration effect by NO emissions from shipping. Ozone depletion in EMEP/MSC-W is weaker than in the other two models, due to a combination of the larger grid cell volume causing higher initial dilution of ship emissions and of the use of monthly averages for ship emissions, averaging out hourly emission peaks.

The ship-related PM2.5 affects the coastal areas in the Baltic Sea region, as its influence extends further inland than is the case for ship-related NO2 and SO2. The average contribution of ships to levels of PM2.5 is calculated by the three CTMs ranges between 4.15 % and 6.5 % in the entire Baltic Sea region and between 3.15 % and 5.7 % in the coastal land areas. Differences in ship-related PM2.5 among the models are mainly attributed to differences in the schemes for inorganic aerosol formation. Since shipping emissions of elemental carbon are identical in the three CTMs on a monthly basis, differences for ship-related EC can be explained by differences in the meteorological conditions, the treatment of atmospheric transport and the wet-scavenging parameterizations in the models.

Results obtained from the use of three CTMs give a more robust estimate of the ship contribution to atmospheric pollutant concentrations than a single model. By using several models the sensitivity of the ship contribution to uncertainties in boundary conditions, meteorological data as well as aerosol formation and deposition schemes is taken into account. This is an important step towards a more reliable evaluation of policy options regarding emission regulations for ship traffic and the introduction of a nitrogen emissions control area (NECA) in the Baltic Sea and the North Sea in 2021.

Data availability
Back to toptop
Data availability.

Data from the simulations with CMAQ, SILAM and EMEP/MSC-W on air pollutant concentrations and nitrogen deposition are available upon request.

Appendix A: Statistical indicators
Back to toptop

The performance of each model is quantified in terms of mean values (μMod and μObs), normalized mean bias (NMB), Spearman's correlation coefficient (R), root mean square error (RMSE) and FAC2.

The normalized mean bias is given by the following:

$\begin{array}{}\text{(A1)}& \text{NMB}=\frac{\stackrel{\mathrm{‾}}{M}-\stackrel{\mathrm{‾}}{O}}{\stackrel{\mathrm{‾}}{O}}\phantom{\rule{0.33em}{0ex}},\end{array}$

where M and O stand for model and observation results, respectively. The overbars indicate the time average over N time intervals (number of observations). The time average is done for 1 year; hence $\stackrel{\mathrm{‾}}{M}$ corresponds to μMod and $\stackrel{\mathrm{‾}}{O}$ corresponds to μObs.

The Spearman's correlation coefficient is defined as the Pearson's correlation coefficient between the ranked variables. For a sample of size N, the N raw scores Mi and Oi for each time step i are converted to ranks rg(Mi) and rg(Oi). The Spearman's correlation coefficient is then computed from

$\begin{array}{}\text{(A2)}& \mathrm{R}=\frac{\mathrm{cov}\left(\mathrm{rg}\left(M\right),\mathrm{rg}\left(O\right)\right)}{{\text{SDR}}_{M}\cdot {\text{SDR}}_{O}},\end{array}$

where cov(rg(M),rg(O)) is the covariance of the rank variables. SDRM and SDRO are the standard deviations of the ranks of the model and observation data. The Spearman correlation between two variables is equal to the Pearson correlation between the rank values of those two variables; while Pearson's correlation assesses linear relationships, Spearman's correlation assesses monotonic relationships (whether linear or not). If there are no repeated data values, a perfect Spearman correlation of +1 or −1 occurs when each of the variables is a perfect monotone function of the other.

The root mean squared error combines the magnitudes of the errors in predictions for various times into a single measure and is defined as follows:

$\begin{array}{}\text{(A3)}& \text{RMSE}=\sqrt{\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}{\left({M}_{i}-{O}_{i}\right)}^{\mathrm{2}}}\phantom{\rule{0.33em}{0ex}}.\end{array}$

RMSE is a measure of accuracy and allows prediction errors of different models to be compared for a particular dataset.

FAC2 is the fraction of model values within a factor of 2 of the observations :

Appendix B: Time series plots
Back to toptop

The appendix contains the time series plots of NO2 concentrations at selected ship-influenced monitoring stations for the meteorological year 2012. Time series plots of O3 concentrations at the corresponding sites are shown in Fig. S1.

Figure B1Ship influence on daily mean NO2 concentration (µg m−3): (a) EE0011A, (b) FI00349, (c) SE0066A, (d) DESH023, (e) FI00351, (f) EE0009R, (g) PL0053A, (h) PL0171A, (i) DEMV021, (j) DESH033, (k) DK0053A, and (l) SE0022A. Model data from CMAQ (dark red line), SILAM (blue line), EMEP/MSC-W (green line), the respective noship run (dashed red line) and observations (black line with open circles).

Supplement
Back to toptop
Supplement.

The supplement related to this article is available online at: https://doi.org/10.5194/acp-19-7019-2019-supplement.

Author contributions
Back to toptop
Author contributions.

MK was responsible for the overall structure and did most of the writing, performed simulations with CMAQ, developed a framework for statistical processing, visualization and plotting and evaluated the air concentration data from all models. JEJ performed simulations with the EMEP/MSC-W model and provided model data for the comparison. AU and MP performed SILAM simulations and provided model output. AA carried out the statistical evaluation of the ship contributions and programmed R-scripts for visualization and data processing. MS gave advice on performing the intercomparison and SILAM simulations. JPJ and LJ provided ship emissions with STEAM. MQ developed research questions, discussed the manuscript structure and assisted with the introduction, conclusions and abstract. VM developed research questions and coordinated the intercomparison. All contributed to the writing of the paper and discussion of the model results.

Competing interests
Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Back to toptop
Special issue statement.

This article is part of the special issue “Shipping and the Environment – From Regional to Global Perspectives (ACP/OS inter-journal SI)”. It is a result of the Shipping and the Environment – From Regional to Global Perspectives, Gothenburg, Sweden, 23–24 October 2017.

Acknowledgements
Back to toptop
Acknowledgements.

This work resulted from the BONUS SHEBA project and was supported by BONUS (Art 185), which is jointly funded by the EU, the Academy of Finland and Forschungszentrum Jülich Beteiligungsgesellschaft mbH. The contributions from the EMEP model were funded by the European Union (European Regional Development Fund) project EnviSuM and partially by EMEP under UNECE. The Baltic Sea Regional Development Fund is acknowledged for financial support. We are grateful to the HELCOM member states for allowing the use of HELCOM AIS data in this research. We thank Jana Moldanová (IVL) and Malin Gustafsson (IVL) for providing WRF data that were used in the SILAM simulation. Johannes Bieser (HZG) is thanked for providing SMOKE-EU emission datasets and Ronny Petrik (HZG) for evaluation of COSMO-CLM data using the HZG Evaluation Suite. The computational time for EMEP model runs was supported by the Research Council of Norway through the NOTUR project EMEP (NN2890K) for CPU and NorStore project European Monitoring and Evaluation Programme (NS9005K) for storage of data. The air quality model CMAQ is developed and maintained by the U.S. Environmental Protection Agency. COSMO-CLM is the community model of the German climate research. The simulations with COSMO-CLM and CMAQ were performed at the German Climate Computing Centre (DKRZ) within the project “Regional Atmospheric Modelling” (project ID 0302). The European Environmental Agency (EEA) is thanked for the AirBase air quality database system. We acknowledge the E-OBS dataset from the EU-FP6 project UERRA (http://www.uerra.eu, last access: 17 April 2019) and the Copernicus Climate Change Service and the data providers in the ECA&D project (https://www.ecad.eu, last access: 17 April 2019).

Financial support
Back to toptop
Financial support.

The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.

Review statement
Back to toptop
Review statement.

This paper was edited by Huan Liu and reviewed by Huan Liu and one anonymous referee.

References
Back to toptop

Andersson, C. and Engardt, M.: European ozone in a future climate: Importance of changes in dry deposition and isoprene emissions, J. Geophys. Res.-Atmos., 115, D02303, https://doi.org/10.1029/2008JD011690, 2010. a

Andersson-Sköld, Y. and Simpson, D.: Comparison of the chemical schemes of the EMEP MSC-W and the IVL photochemical trajectory models, Atmos. Environ., 33, 1111–1129, 1999. a

Angelbratt, J., Mellqvist, J., Simpson, D., Jonson, J. E., Blumenstock, T., Borsdorff, T., Duchatelet, P., Forster, F., Hase, F., Mahieu, E., De Mazière, M., Notholt, J., Petersen, A. K., Raffalski, U., Servais, C., Sussmann, R., Warneke, T., and Vigouroux, C.: Carbon monoxide (CO) and ethane (C2H6) trends from ground-based solar FTIR measurements at six European stations, comparison and sensitivity analysis with the EMEP model, Atmos. Chem. Phys., 11, 9253–9269, https://doi.org/10.5194/acp-11-9253-2011, 2011. a

Antturi, J., Hänninen, O., Jalkanen, J.-P., Johansson, L., Prank, M., Sofiev, M., and Ollikainen, M.: Costs and benefits of low-sulphur fuel standard for Baltic Sea shipping, J. Environ. Manage., 184, 431–440, 10.1016/j.jenvman.2016.09.064, 2016. a, b, c

Amann, M., Borken-Kleefeld, J., Cofala, J., Heyes, C., Klimont, Z., Rafaj, P., Purohit, P., Schoepp, W. and Winiwarter, W.: Future emissions of air pollutants in Europe – Current legislation baseline and the scope for further reductions, TSAP Report #1, International Institute for Applied Systems Analysis, Laxenburg, Austria, 2012. a

Amann, M., Klimont, Z., and Wagner, F.: Regional and Global Emissions of Air Pollutants: Recent Trends and Future Scenarios, Annu. Rev. Env. Resour., 38, 31–55, 2013. a

Appel, K. W., Pouliot, G. A., Simon, H., Sarwar, G., Pye, H. O. T., Napelenok, S. L., Akhtar, F., and Roselle, S. J.: Evaluation of dust and trace metal estimates from the Community Multiscale Air Quality (CMAQ) model version 5.0, Geosci. Model Dev., 6, 883–899, https://doi.org/10.5194/gmd-6-883-2013, 2013. a

Appel, K. W., Napelenok, S. L., Foley, K. M., Pye, H. O. T., Hogrefe, C., Luecken, D. J., Bash, J. O., Roselle, S. J., Pleim, J. E., Foroutan, H., Hutzell, W. T., Pouliot, G. A., Sarwar, G., Fahey, K. M., Gantt, B., Gilliam, R. C., Heath, N. K., Kang, D., Mathur, R., Schwede, D. B., Spero, T. L., Wong, D. C., and Young, J. O.: Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1, Geosci. Model Dev., 10, 1703–1732, https://doi.org/10.5194/gmd-10-1703-2017, 2017. a

Aulinger, A., Matthias, V., Zeretzke, M., Bieser, J., Quante, M., and Backes, A.: The impact of shipping emissions on air pollution in the greater North Sea region – Part 1: Current emissions and concentrations, Atmos. Chem. Phys., 16, 739–758, https://doi.org/10.5194/acp-16-739-2016, 2016. a, b

Ausmeel, S., Eriksson, A., Ahlberg, E., and Kristensson, A.: Methods for identifying aged ship plumes and estimating contribution to aerosol exposure downwind of shipping lanes, Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2018-445, in review, 2019. a

Backes, A., Aulinger, A., Bieser, J., Matthias, V., and Quante, M.: Ammonia emissions in Europe, part I: Development of a dynamical ammonia emission inventory, Atmos. Environ., 131, 55–66, 2016. a

Bergström, R., Denier van der Gon, H. A. C., Prévôt, A. S. H., Yttri, K. E., and Simpson, D.: Modelling of organic aerosols over Europe (2002–2007) using a volatility basis set (VBS) framework: application of different assumptions regarding the formation of secondary organic aerosol, Atmos. Chem. Phys., 12, 8499–8527, https://doi.org/10.5194/acp-12-8499-2012, 2012. a

Bieser, J., Aulinger, A., Matthias, V., Quante, M., and Builtjes, P.: SMOKE for Europe – adaptation, modification and evaluation of a comprehensive emission model for Europe, Geosci. Model Dev., 4, 47–68, https://doi.org/10.5194/gmd-4-47-2011, 2011a. a, b

Bieser, J., Aulinger, A., Matthias, V., Quante, M., and Denier van der Gon, H. A. C.: Vertical emission profiles for Europe based on plume rise calculations, Environ. Pollut., 159, 2935–2946, https://doi.org/10.1016/j.envpol.2011.04.030, 2011b. a

Binkowski, F. S. and Roselle, S. J.: Models-3 Community Multiscale Air Quality (CMAQ) model aerosol component 1. Model description, J. Geophys. Res.-Atmos., 108, 4183, https://doi.org/10.1029/2001JD001409, 2003. a

Binkowski, F. S. and Shankar, U.: The Regional Particulate Matter Model: 1. Model description and preliminary results, J. Geophys. Res.-Atmos., 100, 26191–26209, https://doi.org/10.1029/95JD02093, 1995. a

Byun, D. and Schere, K.: Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System, Appl. Mech. Rev., 59, 51–77, 2006. a, b

Carlton, A. G., Bhave, P. V., Napelenok, S. L., Edney, E. D., Sarwar, G., Pinder, R. W., Pouliot, G. A., and Houyoux, M.: Model representation of secondary organic aerosol in CMAQ v4.7, Environ. Sci. Technol., 44, 8553–8560, https://doi.org/10.1021/es100636q, 2010. a

Cofala, J., Amann, M., Heyes, C., Klimont, Z., Posch, M., Schöpp, W., Tarasson, L., Jonson, J., Whall, C., and Stavrakaki, A.: Final Report: Analysis of Policy Measures to Reduce Ship Emissions in the Context of the Revision of the National Emissions Ceilings Directive, March 2007, International Institute for Applied Systems Analysis, Laxenburg, Austria, 74 pp., 2007. a

Colette, A., Granier, C., Hodnebrog, Ø., Jakobs, H., Maurizi, A., Nyiri, A., Bessagnet, B., D'Angiola, A., D'Isidoro, M., Gauss, M., Meleux, F., Memmesheimer, M., Mieville, A., Rouïl, L., Russo, F., Solberg, S., Stordal, F., and Tampieri, F.: Air quality trends in Europe over the past decade: a first multi-model assessment, Atmos. Chem. Phys., 11, 11657–11678, https://doi.org/10.5194/acp-11-11657-2011, 2011. a, b

Colette, A., Granier, C., Hodnebrog, Ø., Jakobs, H., Maurizi, A., Nyiri, A., Rao, S., Amann, M., Bessagnet, B., D'Angiola, A., Gauss, M., Heyes, C., Klimont, Z., Meleux, F., Memmesheimer, M., Mieville, A., Rouïl, L., Russo, F., Schucht, S., Simpson, D., Stordal, F., Tampieri, F., and Vrac, M.: Future air quality in Europe: a multi-model assessment of projected exposure to ozone, Atmos. Chem. Phys., 12, 10613–10630, https://doi.org/10.5194/acp-12-10613-2012, 2012. a, b

Corbett, J. J., Winebrake, J. J., Green, E. H., Kasibhatla, P., Eyring, V., and Lauer, A.: Mortality from ship emissions: a global assessment, Environ. Sci. Technol., 41, 8512–8518, 2007. a

Cornes, R., van der Schrier, G., van den Besselaar, E. J. M., and Jones, P. D.: An ensemble version of the E-OBS temperature and precipitation datasets, J. Geophys. Res.-Atmos., 123, 9391–9409, https://doi.org/10.1029/2017JD028200, 2018. a

Claremar, B., Haglund, K., and Rutgersson, A.: Ship emissions and the use of current air cleaning technology: contributions to air pollution and acidification in the Baltic Sea, Earth Syst. Dynam., 8, 901–919, https://doi.org/10.5194/esd-8-901-2017, 2017. a, b

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, 2011. a

Dockery, D. W., Pope III, C. A., Xu, X., Spengler, J. D., Ware, J. H., Fay, M. E., Ferris Jr., B. G., and Speizer, F. E.: An association between air pollution and mortality in six U.S. cities, New Engl. J. Med., 329, 1753–1759, https://doi.org/10.1056/NEJM199312093292401, 1993. a

Donahue, N. M., Epstein, S. A., Pandis, S. N., and Robinson, A. L.: A two-dimensional volatility basis set: 1. organic-aerosol mixing thermodynamics, Atmos. Chem. Phys., 11, 3303–3318, https://doi.org/10.5194/acp-11-3303-2011, 2011. a

Donahue, N. M., Kroll, J. H., Pandis, S. N., and Robinson, A. L.: A two-dimensional volatility basis set – Part 2: Diagnostics of organic-aerosol evolution, Atmos. Chem. Phys., 12, 615–634, https://doi.org/10.5194/acp-12-615-2012, 2012. a

Dore, A. J., Carslaw, D. C., Braban, C., Cain, M., Chemel, C., Conolly, C., Derwent, R. G., Griffiths, S. J., Hall, J., Hayman, G., Lawrence, S., Metcalfe, S. E., Redington, A., Simpson, D., Sutton, M. A., Sutton, P., Tang, Y. S., Vieno, M., Werner, M., and Whyatt, J. D.: Evaluation of the performance of different atmospheric chemical transport models and intercomparison of nitrogen and sulphur deposition estimates for the UK, Atmos. Environ., 119, 131–143, 2015. a

EMEP: The EMEP/MSC-W Model – User's Guide, version rv4.8, October 2015, available at: https://github.com/metno/emep-ctm/releases/tag/rv4_8 (last access: 1 November 2018), 2015. a

Endresen, Ø., Sørgård, E., Sundet, J. K., Dalsøren, S. B., Isaksen, I. S., Berglen, T. F., and Gravir, G.: Emission from international sea transportation and environmental impact, J. Geophys. Res.-Atmos., 108, 4560, https://doi.org/10.1029/2002JD002898, 2003. a

EUROSTAT, Maritime transport statistics – short sea shipping of goods, available at: https://ec.europa.eu/eurostat/statistics-explained/index.php/Maritime_transport_statistics_-_short_sea_shipping_of_goods, last access: 19 October 2018. a

Eyring, V., Isaksen, I., Berntsen, T., Collins, W., Corbett, J., Endresen, Ø., Grainger, R., Moldanová, J., Schlager, H., and Stevenson, D.: Transport impacts on atmosphere and climate: shipping, Atmos. Environ., 44, 4735–4771, 2010. a

Foley, K. M., Roselle, S. J., Appel, K. W., Bhave, P. V., Pleim, J. E., Otte, T. L., Mathur, R., Sarwar, G., Young, J. O., Gilliam, R. C., Nolte, C. G., Kelly, J. T., Gilliland, A. B., and Bash, J. O.: Incremental testing of the Community Multiscale Air Quality (CMAQ) modeling system version 4.7, Geosci. Model Dev., 3, 205–226, https://doi.org/10.5194/gmd-3-205-2010, 2010. a

Fountoukis, C. and Nenes, A.: ISORROPIA II: a computationally efficient thermodynamic equilibrium model for K+-Ca2+-Mg2+-${\mathrm{NH}}_{\mathrm{4}}^{+}$-Na+-${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$-${\mathrm{NO}}_{\mathrm{3}}^{-}$-Cl-H2O aerosols, Atmos. Chem. Phys., 7, 4639–4659, https://doi.org/10.5194/acp-7-4639-2007, 2007. a

Fridell, E., Steen, E., and Peterson, K.: Primary particles in ship emissions, Atmos. Environ., 42, 1160–1168, 2008. a

Giordano, L., Brunner, D., Flemming, J., Hogrefe, C., Im, U., Bianconi, R., Badia, A., Balzarini, A., Baró, R., Chemel, C., Curci, G., Forkel, R., Jiménez-Guerrero, P., X Hirtl, M., Hodzic, A., Honzak, L., Jorba, O., Knote, C., Kuenen, J. J. P., Makar, P. A., Manders-Groot, A., Neal, L., Pérez, J. L., Pirovano, G., Pouliot, G., San José, R., Savage, N., Schröder, W., Sokhi, R S., Syrakov, D., Torian, A., Tuccella, P., Werhahn, J., Wolke, R., Yahya, K., Žabkar, R., Zhang, Y., and Galmarini, S., Assessment of the MACC reanalysis and its influence as chemical boundary conditions for regional air quality modeling in AQMEII-2, Atmos. Environ., 115, 371–388, https://doi.org/10.1016/j.atmosenv.2015.02.034, 2015. a, b

Gong, S. L.: A parameterization of sea-salt aerosol source function for sub- and super-micron particles, Global Biogeochem. Cy., 17, 1097, https://doi.org/10.1029/2003GB002079, 2003. a

Gryning, S.-E. and Batchvarova, E.: Marine boundary layer and turbulent fluxes over the Baltic Sea: measurements and modelling, Bound.-Lay. Meteorol., 103, 29–47, 2002. a

Guevara, M., Lopez-Aparicio, S., Cuvelier, C., Tarrason, L., Clappier, A., and Thunis, P.: A benchmarking tool to screen and compare bottom up and top-down atmospheric emission inventories, Air Qual. Atmos. Hlth., 10, 627–642, https://doi.org/10.1007/s11869-016-0456-6, 2016. a

Haiden, T., Janousek, M., Bauer, P., Bidlot, J., Ferranti, L., Hewson, T., Prates, F., Richardson, D. S., and Vitart, F.: Evaluation of ECMWF forecasts, including 2013–2014 upgrades, ECMWF Technical Memorandum No. 742, December 2014, European Centre for Medium-Range Weather Forecasts, Reading, UK, 2014. a

Hanna, S. and Chang, J.: Acceptance criteria for urban dispersion model evaluation, Meteorol. Atmos. Phys., 116, 133–146, https://doi.org/10.1007/s00703-011-0177-1, 2012. a, b

HELCOM: Eutrophication in the Baltic Sea – An integrated thematic assessment of the effects of nutrient enrichment and eutrophication in the Baltic Sea region, Balt. Sea Environ. Proc., 115B, 1–152, 2009. a

HELCOM: Annual report on shipping accidents in the Baltic Sea area during 2012, HELCOM – Baltic Marine Environment Protection Commission, Helsinki, Finland, 43 pp., 2014. a

Hettelingh, J.-P., Posch, M., and Slootweg, J.: European critical loads: database, biodiversity and ecosystems at risk, CCE Final Report 2017, Coordination Centre for Effects, RIVM Report 2017-0155, National Institute for Public Health and the Environment, Bilthoven, The Netherlands, 2017. a

Hongisto, M.: Impact of the emissions of international sea traffic on airborne deposition to the Baltic Sea and concentrations at the coastline, Oceanologia, 56, 349–372, https://doi.org/10.5697/oc.56-2.349, 2014. a

Hunter, K. A., Liss, P. S., Surapipith, V., Dentener, F., Duce, R., Kanakidou, M., Kubilay, N., Mahowald, N., Okin, G., Sarin, M., Uematsu, M., and Zhu, T.: Impacts of anthropogenic SOx, NOx and NH3 on acidification of coastal waters and shipping lanes, Geophys. Res. Lett., 38, L13602, https://doi.org/10.1029/2011GL047720, 2011. a

IMO: The International Convention for the Prevention of Marine Pollution from Ships, 1973 as modified by the Protocol of 1978 relating thereto (MARPOL 73/78), Annex VI: Prevention of Air Pollution from Ships, 2008 Amendments, International Maritime Organization, 2008a. a

Jalkanen, J.-P., Brink, A., Kalli, J., Pettersson, H., Kukkonen, J., and Stipa, T.: A modelling system for the exhaust emissions of marine traffic and its application in the Baltic Sea area, Atmos. Chem. Phys., 9, 9209–9223, https://doi.org/10.5194/acp-9-9209-2009, 2009. a, b

Jalkanen, J.-P., Johansson, L., Kukkonen, J., Brink, A., Kalli, J., and Stipa, T.: Extension of an assessment model of ship traffic exhaust emissions for particulate matter and carbon monoxide, Atmos. Chem. Phys., 12, 2641–2659, https://doi.org/10.5194/acp-12-2641-2012, 2012. a, b, c

Johansson, L., Jalkanen, J.-P., Kalli, J., and Kukkonen, J.: The evolution of shipping emissions and the costs of regulation changes in the northern EU area, Atmos. Chem. Phys., 13, 11375–11389, https://doi.org/10.5194/acp-13-11375-2013, 2013. a

Johansson, L., Jalkanen, J.-P., and Kukkonen, J.: Global assessment of shipping emissions in 2015 on a high spatial and temporal resolution. Atmos. Environ., 167, 403–415, https://doi.org/10.1016/j.atmosenv.2017.08.042, 2017. a

Jonson, J. E., Simpson, D., Fagerli, H., and Solberg, S.: Can we explain the trends in European ozone levels?, Atmos. Chem. Phys., 6, 51–66, https://doi.org/10.5194/acp-6-51-2006, 2006. a

Jonson, J. E., Tarrason, L., Klein, H., Vestreng, V., Cofala, J., and Whall, C.: Effects of ship emissions on European ground level ozone in 2020, Int. J. Remote Sens., 30, 4099–4110, https://doi.org/10.1080/01431160902821858, 2009. a

Jonson, J. E., Stohl, A., Fiore, A. M., Hess, P., Szopa, S., Wild, O., Zeng, G., Dentener, F. J., Lupu, A., Schultz, M. G., Duncan, B. N., Sudo, K., Wind, P., Schulz, M., Marmer, E., Cuvelier, C., Keating, T., Zuber, A., Valdebenito, A., Dorokhov, V., De Backer, H., Davies, J., Chen, G. H., Johnson, B., Tarasick, D. W., Stübi, R., Newchurch, M. J., von der Gathen, P., Steinbrecht, W., and Claude, H.: A multi-model analysis of vertical ozone profiles, Atmos. Chem. Phys., 10, 5759–5783, https://doi.org/10.5194/acp-10-5759-2010, 2010. a

Jonson, J. E., Jalkanen, J. P., Johansson, L., Gauss, M., and Denier van der Gon, H. A. C.: Model calculations of the effects of present and future emissions of air pollutants from shipping in the Baltic Sea and the North Sea, Atmos. Chem. Phys., 15, 783–798, https://doi.org/10.5194/acp-15-783-2015, 2015. a, b, c, d, e

Jonson, J. E., Gauss, M., Schulz, M. and Nyíri, A.: Effects of international shipping, in: Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components. EMEP Status Report 1/2018, Norwegian Meteorological Institute (EMEP/MSC-W), Oslo, Norway, 1–204, 2018a. a, b

Jonson, J. E., Schulz, M., Emmons, L., Flemming, J., Henze, D., Sudo, K., Tronstad Lund, M., Lin, M., Benedictow, A., Koffi, B., Dentener, F., Keating, T., Kivi, R., and Davila, Y.: The effects of intercontinental emission sources on European air pollution levels, Atmos. Chem. Phys., 18, 13655–13672, https://doi.org/10.5194/acp-18-13655-2018, 2018b. a, b

Karl, M., Castell, N., Simpson, D., Solberg, S., Starrfelt, J., Svendby, T., Walker, S.-E., and Wright, R. F.: Uncertainties in assessing the environmental impact of amine emissions from a CO2 capture plant, Atmos. Chem. Phys., 14, 8533–8557, https://doi.org/10.5194/acp-14-8533-2014, 2014. a

Karl, M., Bieser, J., Geyer, B., Matthias, V., Jalkanen, J.-P., Johansson, L., and Fridell, E.: Impact of a nitrogen emission control area (NECA) on the future air quality and nitrogen deposition to seawater in the Baltic Sea region, Atmos. Chem. Phys., 19, 1721–1752, https://doi.org/10.5194/acp-19-1721-2019, 2019. a, b, c, d, e, f

Kelly, J. T., Bhave, P. V., Nolte, C. G., Shankar, U., and Foley, K. M.: Simulating emission and chemical evolution of coarse sea-salt particles in the Community Multiscale Air Quality (CMAQ) model, Geosci. Model Dev., 3, 257–273, https://doi.org/10.5194/gmd-3-257-2010, 2010. a

Kouznetsov, R. and Sofiev, M.: A methodology for evaluation of vertical dispersion and dry deposition of atmospheric aerosols, J. Geophys. Res.-Atmos., 117, D01202, https://doi.org/10.1029/2011JD016366, 2012. a

Kouznetsov, R. and Sofiev, M.: A New Wet Deposition Module in SILAM Chemical Transport Model, American Geophysical Union, Fall Meeting 2013, abstract A51G-0117, 2013. a

Kryza, M., Wałaszek, K., Ojrzyńska, H., Szymanowski, M., Werner, M., and Dore, A. J.: High-resolution dynamical downscaling of ERA-Interim using the WRF regional climate model for the area of Poland. Part 1: Model configuration and statistical evaluation for the 1981-2010 period, Pure Appl. Geophys., 174, 511–526, https://doi.org/10.1007/s00024-016-1272-5, 2017. a, b

Kukkonen, J., Olsson, T., Schultz, D. M., Baklanov, A., Klein, T., Miranda, A. I., Monteiro, A., Hirtl, M., Tarvainen, V., Boy, M., Peuch, V.-H., Poupkou, A., Kioutsioukis, I., Finardi, S., Sofiev, M., Sokhi, R., Lehtinen, K. E. J., Karatzas, K., San José, R., Astitha, M., Kallos, G., Schaap, M., Reimer, E., Jakobs, H., and Eben, K.: A review of operational, regional-scale, chemical weather forecasting models in Europe, Atmos. Chem. Phys., 12, 1–87, https://doi.org/10.5194/acp-12-1-2012, 2012. a

Kukkonen, J., Kangas, L., Kauhaniemi, M., Sofiev, M., Aarnio, M., Jaakkola, J. J. K., Kousa, A., and Karppinen, A.: Modelling of the urban concentrations of PM2.5 on a high resolution for a period of 35 years, for the assessment of lifetime exposure and health effects, Atmos. Chem. Phys., 18, 8041–8064, https://doi.org/10.5194/acp-18-8041-2018, 2018. a

Langner, J., Engardt, M., Baklanov, A., Christensen, J. H., Gauss, M., Geels, C., Hedegaard, G. B., Nuterman, R., Simpson, D., Soares, J., Sofiev, M., Wind, P., and Zakey, A.: A multi-model study of impacts of climate change on surface ozone in Europe, Atmos. Chem. Phys., 12, 10423–10440, https://doi.org/10.5194/acp-12-10423-2012, 2012. a

Lee, S., Kim, H. K., Yan, B., Cobb, C. E., Hennigan, C., Nichols, S., Chamber, M., Edgerton, E. S., Jansen, J. J., Hu, Y., Zheng, M., Weber, R. J., and Russell, A. G.: Diagnosis of aged prescribed burning plumes impacting an urban area, Environ. Sci. Technol., 42, 1438–1444, https://doi.org/10.1021/es7023059, 2008. a

Lee, C., Martin, R. V., van Donkelaar, A., Lee, H., Dickerson, R. R., Hains, J. C., Krotkov, N., Richter, A., Vinnikov, K., and Schwab, J. J.: SO2 emissions and lifetimes: estimates from inverse modeling using in situ and global, space-based (SCIAMACHY and OMI) observations, J. Geophys. Res.-Atmos., 116, D06304, https://doi.org/10.1029/2010JD014758, 2011. a

Marécal, V., Peuch, V.-H., Andersson, C., Andersson, S., Arteta, J., Beekmann, M., Benedictow, A., Bergström, R., Bessagnet, B., Cansado, A., Chéroux, F., Colette, A., Coman, A., Curier, R. L., Denier van der Gon, H. A. C., Drouin, A., Elbern, H., Emili, E., Engelen, R. J., Eskes, H. J., Foret, G., Friese, E., Gauss, M., Giannaros, C., Guth, J., Joly, M., Jaumouillé, E., Josse, B., Kadygrov, N., Kaiser, J. W., Krajsek, K., Kuenen, J., Kumar, U., Liora, N., Lopez, E., Malherbe, L., Martinez, I., Melas, D., Meleux, F., Menut, L., Moinat, P., Morales, T., Parmentier, J., Piacentini, A., Plu, M., Poupkou, A., Queguiner, S., Robertson, L., Rouïl, L., Schaap, M., Segers, A., Sofiev, M., Tarasson, L., Thomas, M., Timmermans, R., Valdebenito, Á., van Velthoven, P., van Versendaal, R., Vira, J., and Ung, A.: A regional air quality forecasting system over Europe: the MACC-II daily ensemble production, Geosci. Model Dev., 8, 2777–2813, https://doi.org/10.5194/gmd-8-2777-2015, 2015. a

Matthias, V., Bewersdorff, I., Aulinger, A., and Quante, M.: The contribution of ship emissions to air pollution in the North Sea regions, Environ. Pollut., 158, 2241–2250, https://doi.org/10.1016/j.envpol.2010.02.013, 2010. a, b

Matthias, V., Aulinger, A., Backes, A., Bieser, J., Geyer, B., Quante, M., and Zeretzke, M.: The impact of shipping emissions on air pollution in the greater North Sea region – Part 2: Scenarios for 2030, Atmos. Chem. Phys., 16, 759–776, https://doi.org/10.5194/acp-16-759-2016, 2016. a, b, c, d

Moldanová, J., Fridell, E., Popovicheva, O., Demirdjian, B., Tishkova, V., Faccinetto, A., and Focsa, C.: Characterisation of particulate matter and gaseous emissions from a large ship diesel engine, Atmos. Environ., 43, 2632–2641, https://doi.org/10.1016/j.atmosenv.2009.02.008, 2009. a

Myers, L. and Sirois, M. J.: Spearman Correlation Coefficients, Differences between, in: Encyclopedia of Statistical Sciences, edited by: Kotz, S., Read, C. B., Balakrishnan, N., Vidakovic, B., and Johnson, N. L., https://doi.org/10.1002/0471667196.ess5050.pub2, 2006. a, b

Nenes, A., Pandis, S. N., and Pilinis, C.: Continued development and testing of a new thermodynamic aerosol module for urban and regional air quality models, Atmos. Environ., 33, 1553–1560, https://doi.org/10.1016/S1352-2310(98)00352-5, 1999. a

Neumann, D., Matthias, V., Bieser, J., Aulinger, A., and Quante, M.: A comparison of sea salt emission parameterizations in northwestern Europe using a chemistry transport model setup, Atmos. Chem. Phys., 16, 9905–9933, https://doi.org/10.5194/acp-16-9905-2016, 2016. a

Neumann, D., Friedland, R., Karl, M., Radtke, H., Matthias, V., and Neumann, T.: Importance of high resolution nitrogen deposition data for biogeochemical modeling in the western Baltic Sea and the contribution of the shipping sector, Ocean Sci. Discuss., https://doi.org/10.5194/os-2018-71, 2018. a

Odman, M. T., Hu, Y., Unal, A., Russell, A. G., and Boylan, J. W.: Determining the sources of regional haze in the southeastern U.S. using the CMAQ model, J. Appl. Meteorol. Clim., 46, 1731–1743, https://doi.org/10.1175/2007JAMC1430.1, 2007. a

Otte, T. L. and Pleim, J. E.: The Meteorology-Chemistry Interface Processor (MCIP) for the CMAQ modeling system: updates through MCIPv3.4.1, Geosci. Model Dev., 3, 243–256, https://doi.org/10.5194/gmd-3-243-2010, 2010. a, b, c

Paerl, H. W.: Coastal Eutrophication in relation to atmospheric nitrogen deposition: Current perspectives, Ophelia, 41, 237–259, 1995. a

Park, R. J., Jacob, D. J., Palmer, P. I., Clarke, A. D., Weber, R. J., Zondlo, M. A., Eisele, F. L., Bandy, A. R., Thornton, D. C., Sachse, G. W., and Bond, T. C.: Export efficiency of black carbon aerosol in continental outflow: Global implications, J. Geophys. Res., 110, D11205, https://doi.org/10.1029/2004JD005432, 2005. a

Poupkou, A., Giannaros, T., Markakis, K., Kioutsioukis, I., Curci, G., Melas, D., and Zerefos, C.: A model for European Biogenic Volatile Organic Compound emissions: Software development and first validation, Environ. Modell. Softw., 25, 1845–1856, 2010. a

Prank, M., Sofiev, M., Tsyro, S., Hendriks, C., Semeena, V., Vazhappilly Francis, X., Butler, T., Denier van der Gon, H., Friedrich, R., Hendricks, J., Kong, X., Lawrence, M., Righi, M., Samaras, Z., Sausen, R., Kukkonen, J., and Sokhi, R.: Evaluation of the performance of four chemical transport models in predicting the aerosol chemical composition in Europe in 2005, Atmos. Chem. Phys., 16, 6041–6070, https://doi.org/10.5194/acp-16-6041-2016, 2016. a, b, c

Prank, M., Vira, J., Ots, R., and Sofiev, M.: Evaluation of Organic Aerosol and Its Precursors in the SILAM Model, in: Air Pollution Modeling and its Application XXV, edited by: Mensink C. and Kallos G., ITM 2016, Springer Proceedings in Complexity, Springer, Cham, 2018. a

Pye, H. O. T. and Pouliot, G. A.: Modeling the role of alkanes, polycyclic aromatic hydrocarbons, and their oligomers in secondary organic aerosol formation, Environ. Sci. Technol., 46, 6041–6047, https://doi.org/10.1021/es300409w, 2012. a

Ramacher, M. O. P., Karl, M., Bieser, J., Jalkanen, J.-P., and Johansson, L.: Urban population exposure to NOx emissions from local shipping in three Baltic Sea harbour cities – a generic approach, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2019-127, in review, 2019. a

Raudsepp, U., Laanemets, J., Maljutenko, I., Hongisto, M., and Jalkanen, J.-P.: Impact of ship-borne nitrogen deposition on the Gulf of Finland ecosystem: an evaluation, Oceanologia, 55, 837–857, https://doi.org/10.5697/oc.55-4.837, 2013. a

Robinson, A. L., Donahue, N. M., Shrivastava, M. K., Weitkamp, E. A., Sage, A. M., Grieshop, A. P., Lane, T. E., Pierce, J. R., and Pandis, S. N.: Rethinking organic aerosols: semivolatile emissions and photochemical aging, Science, 315, 1259–1262, https://doi.org/10.1126/science.1133061, 2007. a, b

Rockel, B., Will, A., and Hense, A.: The Regional Climate Model COSMO-CLM (CCLM), Meteorol. Z., 17, 347–348, 2008. a

Sarwar, G., Simon, H., Bhave, P., and Yarwood, G.: Examining the impact of heterogeneous nitryl chloride production on air quality across the United States, Atmos. Chem. Phys., 12, 6455–6473, https://doi.org/10.5194/acp-12-6455-2012, 2012. a

Schaap, M., Cuvelier, C., Hendriks, C., Bessagnet, B., Baldasano, J. M., Colette, A., Thunis, P., Karam, D., Fagerli, H., Graff, A., Kranenburg, R., Nyíri, A., Pay, M. T., Rouïl, L., Schulz, M., Simpson, D., Stern, R., Terrenoire, E., and Wind, P.: Performance of European chemistry transport models as function of horizontal resolution, Atmos. Environ., 112, 90–105, https://doi.org/10.1016/j.atmosenv.2015.04.003, 2015. a, b, c

Schaub, D., Brunner, D., Boersma, K. F., Keller, J., Folini, D., Buchmann, B., Berresheim, H., and Staehelin, J.: SCIAMACHY tropospheric NO2 over Switzerland: estimates of NOx lifetimes and impact of the complex Alpine topography on the retrieval, Atmos. Chem. Phys., 7, 5971–5987, https://doi.org/10.5194/acp-7-5971-2007, 2007. a

Schwede, D., Pouliot, G., and Pierce, T.: Changes to the Biogenic Emissions Inventory System Version 3 (BEIS3), Proceedings of the 4th CMAS Models-3 Users' Conference 26–28 September 2005, Chapel Hill, NC, available at: http://www.cmascenter.org/conference/2005/abstracts/2_7.pdf (last access: 24 October 2018), 2005. a

Shrivastava, M., Fast, J., Easter, R., Gustafson Jr., W. I., Zaveri, R. A., Jimenez, J. L., Saide, P., and Hodzic, A.: Modeling organic aerosols in a megacity: comparison of simple and complex representations of the volatility basis set approach, Atmos. Chem. Phys., 11, 6639–6662, https://doi.org/10.5194/acp-11-6639-2011, 2011. a

Simoens, D.: AirBase version 8 data products on EEA data service, European Environment Agency, available at: https://www.eea.europa.eu/data-and-maps/data/airbase-the-european-air-quality-database-8 (last access: 24 October 2018), 2014. a

Simpson, D.: Long-period modelling of photochemical oxidants in Europe. Model calculations for July 1985, Atmos. Environ., 26A, 1609–1634, 1992. a

Simpson, D., Butterbach-Bahl, K., Fagerli, H., Kesik, M., Skiba, U., and Tang, S.: Deposition and Emissions of Reactive Nitrogen over European Forests: A Modelling Study, Atmos. Environ., 40, 5712–5726, 2006a. a

Simpson, D., Fagerli, H., Hellsten, S., Knulst, J. C., and Westling, O.: Comparison of modelled and monitored deposition fluxes of sulphur and nitrogen to ICP-forest sites in Europe, Biogeosciences, 3, 337–355, https://doi.org/10.5194/bg-3-337-2006, 2006b. a

Simpson, D., Benedictow, A., Berge, H., Bergström, R., Emberson, L. D., Fagerli, H., Flechard, C. R., Hayman, G. D., Gauss, M., Jonson, J. E., Jenkin, M. E., Nyíri, A., Richter, C., Semeena, V. S., Tsyro, S., Tuovinen, J.-P., Valdebenito, Á., and Wind, P.: The EMEP MSC-W chemical transport model – technical description, Atmos. Chem. Phys., 12, 7825–7865, https://doi.org/10.5194/acp-12-7825-2012, 2012. a, b, c, d, e

Simpson, D., Nyíri, A., Tsyro, S., Valdebenito, Á., and Wind, P.: Updates to the EMEP/MSC-W model, 2015–2016, in: Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components. EMEP Status Report 1/2016, 133–139, The Norwegian Meteorological Institute, Oslo, Norway, 2016. a

Soares, J., Sofiev, M., Hakkarainen, J.: Uncertainties of wild-land fires emission in AQMEII phase 2 case study, Atmos. Environ., 115, 361–370, 2015. a

Sofiev, M.: A model for the evaluation of long-term airborne pollution transport at regional and continental scales, Atmos. Environ., 34, 2481–2493, 2000. a

Sofiev, M., Genikhovich, E., Keronen, P., Vesala, T.: Diagnosing the surface layer parameters for dispersion models within the meteorological-to-dispersion modeling interface, J. Appl. Meteorol. Clim., 49, 221–233, https://doi.org/10.1175/2009JAMC2210.1, 2010. a

Sofiev, M., Soares, J., Prank, M., de Leeuw, G., Kukkonen, J.: A regional-to-global model of emission and transport of sea salt particles in the atmosphere, J. Geophys. Res.-Atmos., 116, D21302, https://doi.org/10.1029/2010JD014713, 2011. a

Sofiev, M., Vira, J., Kouznetsov, R., Prank, M., Soares, J., and Genikhovich, E.: Construction of the SILAM Eulerian atmospheric dispersion model based on the advection algorithm of Michael Galperin, Geosci. Model Dev., 8, 3497–3522, https://doi.org/10.5194/gmd-8-3497-2015, 2015. a, b

Sofiev, M., Kouznetsov, R., Prank, M., Soares, J., Vira, J., Tarvainen, V., and Sofieva, V.: A Long-Term Re-Analysis of Atmospheric Composition and Air Quality, in: Air Pollution Modeling and its Application XXV, edited by: Mensink, C. and Kallos, G., 55–59, Springer International Publishing, Cham, https://doi.org/10.1007/978-3-319-57645-9_9, 2018a. a

Sofiev, M. Winebrake, J. J., Johansson, L., Carr, E. W., Prank, M., Soares, J., Vira, J., Kouznetsov, R., Jalkanen, J.-P., and Corbett, J. J.: Cleaner fuels for ships provide public health benefits with climate tradeoffs, Nat. Commun., 9, 1–12, https://doi.org/10.1038/s41467-017-02774-9, 2018b. a

Solazzo, E., Bianconi, R., Pirovano, G., Matthias, V., Vautard, R., Moran, M. D., Appel, K. W., Bessagnet, B., Brandt, J., Christensen, J. H., Chemel, C., Coll, I., Ferreira, J., Forkel, R., Francis, X. V., Grell, G., Grossi, P., Hansen, A. B., Miranda, A. I., Nopmongcol, U., Prank, M., Sartelet, K. N., Schaap, M., Silver, J. D., Sokhi, R. S., Vira, J., Werhahn, J., Wolke, R., Yarwood, G., Zhang, J., Rao, S. T., and Galmarini, S.: Operational model evaluation for particulate matter in Europe and North America in the context of AQMEII, Atmos. Environ., 53, 75–92, https://doi.org/10.1016/j.atmosenv.2012.02.045, 2012a. a, b

Solazzo, E., Bianconi, R., Vautard, R., Appel K. W., Moran, M. D., Hogrefe, C., Bessagnet, B., Brandt, J., Christensen, J. H., Chemel, C., Coll, I., Denier van der Gon, H., Ferreira, J., Forkel, R., Francis, X. V., Grell, G., Grossi, P., Hansen, A. B., Jeričević, A., Kraljević, L., Miranda, A. I., Nopmongcol, U., Pirovano, G., Prank, M., Riccio, A., Sartelet, K. N., Schaap, M., Silver, J. D., Sokhi, R. S., Vira, J., Werhahn, J., Wolke, R., Yarwood, G., Zhang, J., Rao, S. T., and Galmarini, S.: Model evaluation and ensemble modelling for surface-level ozone in Europe and North America in the context of AQMEII, Atmos. Environ., 53, 60–74, https://doi.org/10.1016/j.atmosenv.2012.01.003, 2012b. a, b

Solazzo, E., Bianconi, R., Pirovano, G., Moran, M. D., Vautard, R., Hogrefe, C., Appel, K. W., Matthias, V., Grossi, P., Bessagnet, B., Brandt, J., Chemel, C., Christensen, J. H., Forkel, R., Francis, X. V., Hansen, A. B., McKeen, S., Nopmongcol, U., Prank, M., Sartelet, K. N., Segers, A., Silver, J. D., Yarwood, G., Werhahn, J., Zhang, J., Rao, S. T., and Galmarini, S.: Evaluating the capability of regional-scale air quality models to capture the vertical distribution of pollutants, Geosci. Model Dev., 6, 791–818, https://doi.org/10.5194/gmd-6-791-2013, 2013. a, b

Solazzo, E., Bianconi, R., Hogrefe, C., Curci, G., Tuccella, P., Alyuz, U., Balzarini, A., Baró, R., Bellasio, R., Bieser, J., Brandt, J., Christensen, J. H., Colette, A., Francis, X., Fraser, A., Vivanco, M. G., Jiménez-Guerrero, P., Im, U., Manders, A., Nopmongcol, U., Kitwiroon, N., Pirovano, G., Pozzoli, L., Prank, M., Sokhi, R. S., Unal, A., Yarwood, G., and Galmarini, S.: Evaluation and error apportionment of an ensemble of atmospheric chemistry transport modeling systems: multivariable temporal and spatial breakdown, Atmos. Chem. Phys., 17, 3001–3054, https://doi.org/10.5194/acp-17-3001-2017, 2017. a, b

Song, C., Chen, G., Hanna, S., Crawford, J., and Davis, D.: Dispersion and chemical evolution of ship plumes in the marine boundary layer: Investigation of O3/NOy/HOX chemistry, J. Geophys. Res., 108, 4143, https://doi.org/10.1029/2002JD002216, 2003. a

Stjern, C. W., Samset, B. H., Myhre, G., Bian, H., Chin, M., Davila, Y., Dentener, F., Emmons, L., Flemming, J., Haslerud, A. S., Henze, D., Jonson, J. E., Kucsera, T., Lund, M. T., Schulz, M., Sudo, K., Takemura, T., and Tilmes, S.: Global and regional radiative forcing from 20 % reductions in BC, OC and SO4 – an HTAP2 multi-model study, Atmos. Chem. Phys., 16, 13579–13599, https://doi.org/10.5194/acp-16-13579-2016, 2016. a

Svensson, N., Bergström, H., Sahlée, E., and Rutgersson, A.: Stable atmospheric conditions over the Baltic Sea: model evaluation and climatology, Boreal Environ. Res., 21, 387–404, 2016. a, b

Textor, C., Schulz, M., Guibert, S., Kinne, S., Balkanski, Y., Bauer, S., Berntsen, T., Berglen, T., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Feichter, H., Fillmore, D., Ghan, S., Ginoux, P., Gong, S., Grini, A., Hendricks, J., Horowitz, L., Huang, P., Isaksen, I., Iversen, I., Kloster, S., Koch, D., Kirkevåg, A., Kristjansson, J. E., Krol, M., Lauer, A., Lamarque, J. F., Liu, X., Montanaro, V., Myhre, G., Penner, J., Pitari, G., Reddy, S., Seland, Ø., Stier, P., Takemura, T., and Tie, X.: Analysis and quantification of the diversities of aerosol life cycles within AeroCom, Atmos. Chem. Phys., 6, 1777–1813, https://doi.org/10.5194/acp-6-1777-2006, 2006. a

Tsyro, S., Nyíri, A., and Wind, P.: Model calculations for 2013 in the new EMEP grid, in: Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components, EMEP/MSC-W Status Report 1/2015, The Norwegian Meteorological Institute, Oslo, Norway, 2015. a

Tsyro, S., Aas, W., Solberg, S., Benedictow, A., Fagerli, H., and Posch, M.: Status of transboundary air pollution in 2016, in: Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components, EMEP/MSC-W Status Report 1/2018, Norwegian Meteorological Institute (EMEP/MSC-W), Oslo, Norway, 204 pp., 2018. a, b

van Aardenne, J., Colette, A., Degraeuwe, B., Hammingh, P., Viana, M., and Vlieger, I.: The impact of international shipping on European air quality and climate forcing, Technical Report No. 4/2013, EEA (European Environment Agency), Copenhagen, Denmark, 1–84, 2013. a

Vautard, R., Moran, M. D., Solazzo, E., Gilliam, R. C., Matthias, V., Bianconi, R., Chemel, C., Ferreira, J., Geyer, B., Hansen, A. B., Jericevic, A., Prank, M., Segers, A., Silver, J. D., Werhahn, J., Wolke, R., Rao, S. T., and Galmarini, S.: Evaluation of the meteorological forcing used for the Air Quality Model Evaluation International Initiative (AQMEII) air quality simulations, Atmos. Environ., 53, 15–37, https://doi.org/10.1016/j.atmosenv.2011.10.065, 2012. a

Vignati, E., Karl, M., Krol, M., Wilson, J., Stier, P., and Cavalli, F.: Sources of uncertainties in modelling black carbon at the global scale, Atmos. Chem. Phys., 10, 2595–2611, https://doi.org/10.5194/acp-10-2595-2010, 2010. a, b

Vinken, G. C. M., Boersma, K. F., Jacob, D. J., and Meijer, E. W.: Accounting for non-linear chemistry of ship plumes in the GEOS-Chem global chemistry transport model, Atmos. Chem. Phys., 11, 11707–11722, https://doi.org/10.5194/acp-11-11707-2011, 2011. a

Vivanco, M. G., Theobald, M. R., García-Gómez, H., Garrido, J. L., Prank, M., Aas, W., Adani, M., Alyuz, U., Andersson, C., Bellasio, R., Bessagnet, B., Bianconi, R., Bieser, J., Brandt, J., Briganti, G., Cappelletti, A., Curci, G., Christensen, J. H., Colette, A., Couvidat, F., Cuvelier, C., D'Isidoro, M., Flemming, J., Fraser, A., Geels, C., Hansen, K. M., Hogrefe, C., Im, U., Jorba, O., Kitwiroon, N., Manders, A., Mircea, M., Otero, N., Pay, M.-T., Pozzoli, L., Solazzo, E., Tsyro, S., Unal, A., Wind, P., and Galmarini, S.: Modeled deposition of nitrogen and sulfur in Europe estimated by 14 air quality model systems: evaluation, effects of changes in emissions and implications for habitat protection, Atmos. Chem. Phys., 18, 10199–10218, https://doi.org/10.5194/acp-18-10199-2018, 2018.  a, b

von Glasow, R., Lawrence, M. G., Sander, R., and Crutzen, P. J.: Modeling the chemical effects of ship exhaust in the cloud-free marine boundary layer, Atmos. Chem. Phys., 3, 233–250, https://doi.org/10.5194/acp-3-233-2003, 2003. a

Vukovich, J. and Pierce, T.: The Implementation of BEIS3 within the SMOKE Modeling Framework, in: Proceedings of the 11th International Emissions Inventory Conference, Atlanta, Georgia, 15–18 April 2002, available at: http://www.epa.gov/ttn/chief/conference/ei11/modeling/vukovich.pdf (last access: 24 October 2018), 2002. a

Wankmüller, R. and Mareckova, K.: Emissions for 2012, in: Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components, EMEP Status Report 1/2014, Norwegian Meteorological Institute (EMEP/MSC-W), Oslo, Norway, 2014 a

Whitten, G. Z., Heo, G., Kimura, Y., McDonald-Buller, E., Allen, D. T., Carter, W. P. L., and Yarwood, G.: A new condensed toluene mechanism for Carbon Bond: CB05-TU, Atmos. Environ., 44, 5346–5355, 2010. a

Yarwood, G., Rao, S., Yocke, M., and Whitten, G. Z.: Updates to the Carbon Bond Mechanism: CB05, Report to the U.S. Environmental Protection Agency, RT-04-00675, available at: http://www.camx.com/publications/default.aspx (last access: 24 October 2018), 2005. a