Next Article in Journal
Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics
Previous Article in Journal
Citric Acid Inhibits Cd Absorption and Transportation by Improving the Antagonism of Essential Elements in Rice Organs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using the Multicomponent Aerosol FORmation Model (MAFOR) to Determine Improved VOC Emission Factors in Ship Plumes

by
Lea Fink
1,
Matthias Karl
1,*,
Volker Matthias
1,
Andreas Weigelt
2,
Matti Irjala
3 and
Pauli Simonen
4
1
Helmholtz-Zentrum Hereon, Department of Coastal Environmental Chemistry, 21052 Geesthacht, Germany
2
Bundesamt für Seeschifffahrt und Hydrographie, 20359 Hamburg, Germany
3
Aeromon Oy, FI-00240 Helsinki, Finland
4
Faculty of Engineering and Natural Sciences, Tampere University, FI-33720 Tampere, Finland
*
Author to whom correspondence should be addressed.
Toxics 2024, 12(6), 432; https://doi.org/10.3390/toxics12060432
Submission received: 21 May 2024 / Revised: 11 June 2024 / Accepted: 12 June 2024 / Published: 14 June 2024
(This article belongs to the Section Air Pollution and Health)

Abstract

:
International shipping’s particulate matter primary emissions have a share in global anthropogenic emissions of between 3% and 4%. Ship emissions of volatile organic compounds (VOCs) can play an important role in the formation of fine particulate matter. Using an aerosol box model for the near-plume scale, this study investigated how the changing VOC emission factor (EF) for ship engines impacts the formation of secondary PM2.5 in ship exhaust plumes that were detected during a measurement campaign. The agreement between measured and modeled particle number size distribution was improved by adjusting VOC emissions, in particular of intermediate-, low-, and extremely low-volatility compounds. The scaling of the VOC emission factor showed that the initial emission factor, based on literature data, had to be multiplied by 3.6 for all VOCs. Information obtained from the box model was integrated into a regional-scale chemistry transport model (CTM) to study the influence of changed VOC ship emissions over the Mediterranean Sea. The regional-scale CTM run with adjusted ship emissions indicated a change in PM2.5 of up to 5% at the main shipping routes and harbor cities in summer. Nevertheless, overall changes due to a change in the VOC EF were rather small, indicating that the size of grid cells in CTMs leads to a fast dilution.

1. Introduction

Marine transport emissions contribute considerably to global air pollution, e.g., [1,2,3]. Ships have the potential to decline the air quality in coastal regions because roughly 70% of ship emissions occur within a distance of 400 km from land [4]. According to Klimont et al. [5], the proportion of international shipping’s particulate matter primary emissions to global anthropogenic emissions is between 3% and 4%, which is comparable to road traffic. Through a combination of chemical reactions and microphysical processes involving precursor pollutants, secondary pollutants are formed in the atmosphere. Secondary atmospheric particulate matter is formed through chemical reactions in the atmosphere involving gaseous precursors.
There is evidence linking PM2.5 (particulate matter < 2.5 µm) exposure to the development of certain lung conditions, cancer, or type 2 diabetes [6,7,8]. These small particles can penetrate deeply into lung alveoli, by inhalation, and cause important adverse effects on human health [9,10,11]. The smaller the particle, the easier and deeper it can enter into pulmonary systems [12]. Long-term exposure to PM2.5 was identified to lead to decreased life expectancy, early death, and morbidity [13,14,15]. Further, there is no safe threshold for PM2.5 concentrations according to the WHO [16].
Volatile organic compounds (VOCs) play an important role as precursors for the formation of fine particulate matter in the atmosphere [17,18]. Most of the particle growth is caused by organic compounds in the gas phase [19,20]. VOCs react in the atmosphere with atmospheric oxidants such as the hydroxyl (OH) radical, the nitrate (NO3) radical, ozone (O3), and chlorine (Cl) atoms [17]. The initial oxidation step produces alkyl or substituted alkyl radicals by the H-atom abstraction from C–H bonds or the addition to C=C bonds. The alkyl radicals react with molecular oxygen (O2) to give organic peroxy radicals, which subsequently can react with nitrogen oxides (NOx), originating from various emission sources, to produce ozone under sunlight. In the fresh exhaust of diesel engines, most primary organic particulate emissions are semi-volatile [21]; thus, to a certain extent, they evaporate accompanied by atmospheric dilution, generating large amounts of low-volatility gas-phase material that is, after several oxidation steps, available for later condensation to particles [22]. Photo-oxidation initiates the formation of particles, as it results in the generation of abundant precursors for the nucleation and growth of fine particles in air [23,24]. Experiments in the laboratory indicate that photo-oxidation of diesel exhaust rapidly produces organic aerosols [22]. Oxidation of gas-phase organic species and chemical aging leads to the production of secondary organic aerosol (SOA) [25,26,27,28]. The SOA is, in large part, responsible for the formation of fine particulate matter [25,29]. SOA contributes 20% to 60% globally and, during periods with high pollution, up to 30% to 77% to total PM2.5 [30,31,32,33]. Particularly in summer, the aerosol mixing ratios in the Mediterranean area are higher than those in most continental European regions [34].
Previous modeling studies have shown the relative potential ship impact on PM2.5. In Karl et al. [35], a regional study on the impacts of ship emissions in the Baltic Sea region was done using three regional-scale chemical transport models. According to that study, the average contribution of ships to PM2.5 levels varied from 4.15% to 6.5% over the whole Baltic Sea region and from 3.15% to 5.7% along the shore. Studies regarding the Mediterranean region found a relative ship impact on the PM2.5 concentrations of a share of 5% to 20% of the total PM2.5 concentrations, e.g., [36,37,38].
Regional-scale chemistry transport model (CTM) systems tend to underestimate the measured PM (particulate matter) [38,39]. This can be traced back to the relatively coarse grid resolution in regional-scale models that do not account for chemical transformation mechanisms in sub-grid pollution plumes. Although regional CTM systems generally contain chemically speciated PM2.5, they do not divide the PM into smaller particles nor consider ultrafine particles (UFPs) with aerodynamic diameters of <0.1 µm [40]. While UFPs have a small contribution to particle mass, the physical processes transforming them are highly dynamic and of great importance to a detailed investigation of PM formation. The physical transformation of ultrafine particles rapidly changes the ambient PNC (particle number concentration) due to nucleation, condensation, evaporation, coagulation, and dry deposition [40]. Through particle growth, a sub-set of these particles (larger than about 50–80 nm) can produce cloud condensation nuclei (CCN), influencing cloud formation and the radiative balance [41].
Gonzales et al. [42] found inland transport of UFPs from ship plumes in Santa Cruz de Tenerife, with ship emissions having a share of 65% to 70% of the UFP concentrations. Pirjola et al. [43] also reported an increase in UFP concentrations due to ship emissions in European coastal cities. Secondary volatile particles make up a large fraction of the UFPs that form in the nucleation of gas-phase species and grow fast in size by condensation of organic vapors and sulfuric acid. Shipping is a likely source of secondary particle formation because of the high sulfur content of the fuel oil they burn [44]. Particles coming from ship exhaust contain several chemical compounds such as particulate sulfate, black carbon, organics, ash, and heavy metals, each having an impact on the properties of particle emission [42,45]. Due to this, particle size distribution in ship plumes can range from 5 nm to >3 µm [3]. Still, few studies have been conducted to investigate the change in particle characteristics in the exhaust plume of ships directly in the stack and in parallel in aged plumes on land [46]. Karl et al. [47] estimated the exposure to UFPs from ship exhaust in major harbor areas in two cities in Finland by measuring and modeling the dispersion and dynamics of aerosols. They found increased UFP concentrations due to ship emissions up to a distance of 3600 m. According to a study conducted in the port area of Marseille, the mean potential impact from ships in July 2020 was 6% to 9% for SOAs and 27% to 51% for total particle number concentrations [48].
Ship plume parametrizations on plume dispersion and chemistry in plumes have been evaluated in previous studies, e.g., [49,50,51], but usually not for aerosols.
To describe the development of the particle number and mass size distribution from the point of emission to the site of interest, the modeling of particle transformations in parallel with plume dispersion is required. It is crucial to evaluate the evolution close to the source with high spatial and temporal resolution since the particle size and composition changes fast after exhaust release [52]. Lagrangian size-resolved aerosol models can be used for investigations when higher resolutions in time and space are necessary.
VOCs are considered to be important for the exhaust plume development and the size distribution, since they might condense shortly after they enter the atmosphere from the emission source, depending on their volatility. Despite the fast condensation during the dilution of the hot exhaust plume, they can still reach gas-phase saturation in the expanding plume, due to the high amounts of VOCs present in ship exhaust. Only a few studies focus on VOCs in ship exhaust. The shift from high-sulfur residual fuel oil to low-sulfur diesel or heavy oils that are often richer in short-chain hydrocarbons was shown to have a major impact on the rise in VOC emissions [53]. The VOC mixture emitted from ships tends to produce more secondary organic aerosol than diesel and gasoline vehicles, while the ozone-forming potential is comparable. The different contents of organic components with high molecular weight is the major reason for this. Ship exhaust contains more heavy VOCs than automotive exhaust (i.e., alkanes and aromatics). Aromatic hydrocarbons play a significant role in SOA formation [54].
The goal of the present study was to simulate particle transformation within aging ship plumes that were observed in a measurement campaign within the SCIPPER project. The aim of this project was to obtain further insights in the capability of several monitoring approaches for shipping emissions (https://www.scipper-project.eu/, accessed on 12 May 2024). The in-plume particle evolution was simulated with an aerosol box model, ship emissions were adjusted, and simulations with the adjusted emissions were performed using a larger-scale CTM. VOC emission factors (EFs) for ships were used to establish a link between the box model and the regional-scale model. Two related questions to be answered in this study are (1) how the changing VOC emission factors for ships impact the formation of particulate matter and (2) how the output of a box model for the near-plume scale can be integrated into the regional-scale CTM. The sectional aerosol box model MAFOR v2.1.0 (Multicomponent Aerosol FORmation model) was applied in the present study [52,55] to simulate the particle size distribution between the ship stack and the surroundings for up to 1 h travel time downwind [47]. Furthermore, in the applied model, gas-phase chemistry is coupled with aerosol dynamics. For the model input, measurement data of a ship exhaust plume measured on board a passenger ferry operated in the Baltic Sea, and air pollutant concentrations and meteorological data, both measured 500 m distant from the shipping lane, were used. The simulations of the field campaign data were used to adjust the VOC emission factors in the ship emission dataset STEAM (version 3.3.0.) of the Finnish Meteorological Institute (FMI) [56,57,58,59]. STEAM ship emissions were further used as input for regional-scale model simulations with the Community Multiscale Air Quality Model (CMAQ) [60]. CMAQ was applied as described in Fink et al. [38,61] on a model domain covering most of the Mediterranean Sea. The Mediterranean Sea was in the focus due to the higher insolation and SOA formation potential compared to the Baltic Sea, leading to a stronger response and better possibility for validation.

2. Materials and Methods

2.1. Measurements on Board and on Land

The MAFOR model runs were initialized by using measurement data as input data. The measurements of ship emissions were done on board a passenger ferry, which is a Passenger/Ro-Ro Cargo Ship powered by two four-stroke Wärtsilä 8ZAL40S engines with a maximum continuous rating (MCR) of 6000 kW and a nominal speed of 510 rev/min. All primary engines have SCR catalysts and are converted to dual-fuel engines that can run either on methanol or regular marine fuels. The used fuel from plumes investigated in the current study was marine gas oil (MGO). Carbon dioxide (CO2), nitrogen monoxide (NO), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), and ammonia (NH3) as initial gases in the ship plume were measured in the ship stack by an Aeromon BH-12 measurement device [62]. The uncertainty for NOx measured by the Aeromon BH-12 is in a range of ±19.7%. The measurement device was located in the stack after the catalyst [63].
Particle number and mass size distributions from the passenger ferry were measured at a resolution of one second by an Engine Exhaust Particle Sizer (EEPS; size distribution: 5.6 nm–560 nm) and an Electrical Low Pressure Impactor (ELPI; size distribution of 7 nm–9.9 µm). The measurement point for these instruments was in the stack after the Selective Catalytic Reduction (SCR) and it was measured downstream of a combination of a porous tube diluter and ejector diluter [63].

Ship Plume Event

The ship plumes were from a regular passenger ferry that travels between Kiel and Gothenburg. All captured plumes originated from that vessel. The ship plume used as a reference for the aerosol formation model MAFOR was detected during the SCIPPER measurement campaign (from 31 August 2021 to 6 September 2021) on 3 September 2021 in Laboe at the Kiel Fjord (“Kieler Förde”) (54.392809, 10.209186; Figure 1). During that event, both onboard and on-shore observations were available. Thus, this plume is used as a reference and all adjustments in the configuration of MAFOR simulations were done based on this occurrence. The plume arrived at the mobile laboratory at 17:41:37 UTC. The evolution of the plume and the decrease in particle number over time was used to determine the time of the plume peak (i.e., when the total concentration at Laboe was highest; Figure A1) and the time when only background particles were measured. The wind direction was 284° (±15° in one hour). If not indicated otherwise, the term ship plume hereafter refers to the reference plume.
For comparison with this reference plume, an additional 12 on-shore plume events from the same ferry between 2021 and 2023 detected at the mobile lab were evaluated in the present study. The plume events were chosen when (i) measured in the afternoon between 16:00 UTC and 18:00 UTC, (ii) departing from Kiel, and (iii) measured in the time frame of mid-July to mid-October. A detailed overview is given in Table A1. The ship emissions were assumed to be the same as in the base case because the ferry passes the measurement location always at the same speed and engine load.
On-shore measurements of particles, gases, and meteorological parameters were carried out at the mobile lab at Laboe close to Kiel. For measuring the particles on shore, a Fast Mobility Particle Sizer (FMPS spectrometer Model 3091; size distribution: 5.6 nm–560 nm, manufacturer: TSI Incorporated) was used for small particles. Larger particles were measured by an OPS (Optical Particle Sizer Model 3330; size distribution: 300 nm–10 µm, manufacturer: TSI Incorporated). The time resolutions for FMPS and OPS were one second. The used air intake system for both particle measurement devices was the Sampling System for Atmospheric Particles from TSI for particles up to 10 µm. The plumes were detected automatically.
Gases and meteorological parameters were measured with an airpointer 4D measurement system (manufacturer: JCT NextGen AQMS mlu-recordum) for NOx, O3, SO2, and a weather station. The lower detection limits are 0.4 ppb for NOx and 0.5 ppb for O3 and SO2. The uncertainty for NOx measurements is given by the manufacturer and is 1 ppb (2 µg/m3) for measured concentrations <500 ppb. The peaks of measured NOx concentrations after the passing of the ships are shown in Figure A2.

2.2. MAFOR

2.2.1. Input Data and Model Configuration

The open-source model MAFOR v2.1.0 was used in the present study. It is a zero-dimensional Lagrangian-type sectional aerosol box model [52,55] available on Github: https://github.com/mafor2/mafor, accessed on 15 May 2024. The model calculates the development of particles in terms of particle number and mass as well as changes in particle average density. The model was run in plume dispersion mode to display a single exhaust plume along one dimension in space. All simulated particles are assumed to be spherical. The model is linked to the MECCA [64,65] chemistry module, which contains up-to-date photolysis rates for VOC chemistry (Figure 2).
The model relies on several input data files that include general data (meteorological data, time, and location), configuration (processes that are used in the model), initial aerosol size distribution (in terms of mass concentration), gas-phase emission (concentration) and background aerosols (particle mass concentration), organics (properties of the organic vapors), and dispersion (plume dispersion and deposition parameterization). The gas measurements and meteorology data obtained at the mobile laboratory on shore were used in the model and defined the air parcels’ trajectory.
The following processes were considered in the MAFOR run: dry deposition over the water surface [66], nucleation, condensation of organics and of water, and coagulation. The coagulation of particles by Brownian diffusion was considered. The dilution of the plume through expansion into the background air with a low particle concentration was done with the dispersion parameterization for ship plumes under convective conditions in the open sea.
The plume dilution after MAFOR model initiation followed the parametrization of Chosson et al. [67] for the convective boundary layer to match the modeled concentrations in the ship plume. The parameters for plume dispersion were not adjusted and used as provided in the MAFOR model. The parameters for plume dispersion in the MAFOR model and the original parameters of Chosson et al. [67] are listed in Table A2. The plume height approximation follows a formulation of von Glasow et al. [68] (Table A3).
The nucleation mechanism of homogeneous H2SO4–H2O [69,70] was applied in the MAFOR simulation. For this, the consideration of H2SO4 is of importance [71]. Starting concentrations for H2SO4 in the stack were set according the concentrations received from Karl et al. [47], with a value of 0.1 × 1011 cm−3.
The EF of NOx was measured on board by the Aeromon BH-12 measurement device [61] in the exhaust stack and was used to calculate the EF for other gases. The initial guess for NOx EF was based on Fridell et al. [72].
The measurements were too close to the emission source for the MAFOR model to start (it starts usually 1 s after the plume left the stack); therefore, a dilution ratio for the gases of 1:8 was used as described in Karl et al. [47]. Karl et al. [47] determined that the plume reached the background temperature after 1 s, when the nucleation was complete and was, therefore, no longer taken into account. Therefore, they only considered the error due to neglecting coagulation during exhaust cooling. In the present study, the simulation of the plume also started after the first second from the point of release; therefore, the initial condensation effect during cooling was not considered. The particle instruments measuring the size distribution in the hot exhaust did not detect significant numbers of particles <7 nm diameter, probably due to the lower size cut-off at 5.6 nm. The SOA concentrations were simulated using a hybrid method of condensation/evaporation and absorptive partitioning into an organic liquid (according to Kerminen et al. [73]) and using a 2D VBS set [74] consisting of nine organic compounds with different volatility. The simulations included the basic gas-phase reactions in the troposphere and the Mainz Organic Mechanism (MOM) as an oxidation scheme for VOCs. The chemistry of the following VOCs were considered in the simulations: ethane (C2H6), ethylene (C2H4), propane (C3H8), trimethylbenzene (TMB), toluene (TOL), xylene (XYL), and primary emitted organic vapors, such as long-chained n-alkanes, which are represented by PIOV (intermediate volatility), PSOV (semi-volatile), and PELV (extremely low-volatile). Reactions were taken from MECCA [64,65].
To start the MAFOR simulation, the particle mass initial at the stack was calculated from the measured particle number. For this, the particles in sizes from 6.04 nm to 523.3 nm measured by EEPS (32 bins) and larger particles with sizes from 764.4 nm to 6.3 µm (4 bins) measured by ELPI were used. The resulting 36 bins were assigned to the four modes in MAFOR (NU: nucleation mode, AI: Aitken mode, AS1: accumulation mode 1, AS2: accumulation mode 2), as described in the following, and the mass was calculated as a sum within those by assuming a density of 1.0 g/cm3. The fourth mode of the measured size distribution at the ship stack is a second accumulation mode rather than a coarse mode; therefore, it was divided into a first and second accumulation mode. The median of in total 366 measurements of particle number and mass on board the ship on 03.09.2021 between 17:49:00 UTC and 22:54:30 UTC was calculated for each bin to determine the input data for the ship plume aerosol. This time window for measurements on the ferry was chosen according to the engine load and fuel used. At this time, the exhaust plume characteristics were expected to be the same as in the plume that came out of the stack at the time when the ship passed the measuring station. Measurements at the exact time of the passing were not available, but they were a few minutes thereafter. Therefore, the median from onboard measurements was used. The received measured aerosol mass was transferred into a MAFOR input file that provides the initial aerosol size distribution of the ship plume. The total aerosol mass concentration was divided into the four modes (NU, AI, AS1, AS2). The aerosol initial mass composition was subsequently converted to particle number based on the material densities of the different aerosol components, assuming spherical particles to ensure consistency in terms of mass and number [52]. The calculation of particle number and mass concentrations in the model was then done using a sectional size representation of the aerosol. The initial chemical composition in the ship plume was based on data from the literature [47,53,72] (Table A4).
Within MAFOR, the total mass was divided into H2SO4, OC, NH4, NO3, MSAp (methane sulfonate), salt, POA (primary organic material), EC, and ash for the ship aerosol and the background aerosol (Table 1). The background composition was taken from Karl et al. [46]. The emissions measured on the ship and their chemical compositions for the reference plume were also used for the simulation of the 12 other ship plumes.
For background aerosols, the particle mass was calculated from the particle number measured at the on-shore mobile lab (from FMPS and OPS devices). To avoid any influence of ship emissions, the median of the measured particle number from half an hour before the arrival of the plume was calculated for the background aerosol data. The distribution of the plume aerosol and background aerosol mass concentrations in the four modes is displayed in Table 1. The mass was distributed to fit the particle number distribution of the respective measurement.
Because of multicomponent condensation and particle coagulation, the composition of particles in the size bins may vary over time. MAFOR calculates rapid changes in the aerosol size distribution every 0.1 s. The simulation of the ship plume starts one second after release from the ship’s stack. The distance from the stack to the mobile laboratory was 530 m (corresponding to 150 s travel time of the ship plume).

2.2.2. Comparison Simulated against Measured Data

The comparison between modeled and measured data was done by comparing particle number concentrations and mean diameter for different size classes on which the measured and the modeled bins were mapped. These classes were set following the size class distribution in Karl et al. [52] by adapting the modeled size distribution with the size classes (bins) of the measured size distribution. The resulting size classes have the diameter ranges S1: 1–10 nm; S2: 10–20 nm; S3: 20–50 nm; S4: 50–100 nm; S5: 100–300 nm; S6: 300–600 nm; S7: >600 nm. The mean diameter of each size class was calculated as follows:
d n ¯ = i ( N i d i ) i N i
where dn denotes the mean diameter per size class, Ni is the number of particles, and di is the midpoint of the size class; i is running over the size bins, n is running over the size classes.

2.3. CMAQ and Ship Emission Dataset

The CMAQ model v5.2 was applied with configurations as described in Fink et al. [38,61] for four months (March, June, September, December) of 2015 on a model domain covering most of the Mediterranean Sea with a resolution of 12 × 12 km2. CMAQ calculates atmospheric concentrations as well as deposition fluxes of gases and aerosols based on emission input data [75,76]. The Carbon Bond 5 (CB05) chemical mechanism with updated toluene chemistry cb05tucl [77] including the chlorine chemistry extension (CB05-TUCL; https://www.airqualitymodeling.org/index.php/CMAQv5.0_Chemistry_Notes, accessed 21 June 2023) was used for the gas-phase chemistry. Secondary inorganic aerosols were formed using the AERO6 aerosol mechanism. The ISORROPIA model was applied for solving the gas-phase–aerosol partition equilibrium of sulfuric acid (H2SO4), nitric acid (HNO3), hydrochloric acid (HCl), and ammonia (NH3) [78,79]. Secondary organic aerosol (SOA) is formed from isoprene, terpenes, benzene, toluene, xylene, and alkanes [80,81]. CMAQ allows for dynamic mass transfer of semi-volatile inorganic gases to coarse mode particles, which facilitates the replacement of chloride by NO3 in sea salt aerosols [82].
Kelly et al. [83] described the method for calculating sea salt emissions as used in this model configuration. Biogenic emissions (NMVOC from plant and soil NO) were previously calculated using the MEGAN Model v3 [84] and then combined with the land-based emissions. Windblown dust emissions and NOx lightning treatment were not taken into account. For land-based anthropogenic emissions, the gridded emissions obtained from the CAMS-REG v2.2 emission inventory were used as input. Gridded emission files contain GNFR (Gridded Nomenclature for Reporting) emission sectors for each country for the air pollutants NOx, SO2, NMVOC, NH3, CO, PM10, PM2.5, and CH4. The detailed distribution and splitting of land-based emissions used for the CMAQ run is described in Fink et al. [38,61].
The ship emission dataset received from STEAM (version 3.3.0) [56,57,58,59] was used as input data for CMAQ. Mineral ash, CO, CO2, elemental carbon (EC), NOx, organic carbon (OC; assumed to be non-volatile), PM2.5, particle number count (PNC), sulfate (SO4), SOx (containing SO2 and SO3), and VOC emissions were provided in two vertical layers (0 m to 36 m; 36 m to 1000 m). VOC emissions were split by FMI into four groups based on their characteristics as a function of engine load to limit the number of produced emissions and the computing resources required to run the STEAM model. The VOC groups contain reactive (volatile) VOCs as well as organic compounds of different volatility. VOC EFs in STEAM were calculated using average values from previous publications [85,86,87,88].
The COSMO model, version COSMO5-CLM16, was used to simulate the meteorological data for CMAQ [89,90]. The MCIP (Meteorology–Chemistry Interface Processor) converted meteorological model output into the CMAQ input format. The meteorological model has a vertical resolution of 40 terrain-following geometric height levels up to 22 km. IFS-CAMS cycle45r1 (Integrated Forecasting System—Copernicus Atmosphere Monitoring Service) [91] was applied as the Boundary Condition driver, with a vertical resolution of 60 sigma levels up to 65 km. CMAQ simulates 30 vertical levels, with the lowest layer ranging from 0 m to 42 m and the second layer ranging from 42 to 85 m.
Two simulations were performed with CMAQ: one run with the ship emission data as provided by STEAM and used in Fink et al. [38,61] and a second run with the ship emissions of the VOC groups adjusted by the adjustment factor derived from MAFOR, as described in Section 2.4, next. The simulated output was afterwards compared with observations at 28 monitoring stations distributed over the domain and as described in Fink et al. [61].

2.4. VOC Emission Factor Scaling

In the following, the adjustment procedure for one VOC EF, independently from volatility, for the ship emissions will be explained. This is done based on the reference plume on 03 September 2021 and by running different sensitivity cases with MAFOR. Toluene, xylene, trimethylbenzene, propane, ethane, and ethylene are VOCs used in MAFOR. The MAFOR model, furthermore, contains three volatility classes for primary delayed organic aerosols represented by PIOV, PSOV, and PELV [52]. Concentrations for PIOV, PSOV, and PELV at 530 m are 0.026 µg/m3, 0.036 µg/m3, and 0.25 µg/m3. These classes contain the condensing VOCs with low vapor pressure. Sulfate in MAFOR is condensed and the emissions of aromatics are oxidized and also contribute to the growth of particles. The formation of secondary aromatic aerosol in the particle phase is neglectable and is not displayed in the investigated time scale. The semi-volatile organic aromatics remained in the gas phase with a concentration at 530 m of 1.64 × 10−4 ng/m3. The input to initialize the MAFOR run is shown in Table A4.
The uncertainties of VOC emissions in ship emission inventories are large. They can range from 0.11 g/kg fuel to 6.93 g/kg fuel [53,56,57,58,59,85,86,87,88].
The approach to receive a new VOC EF is presented in Figure 3. The NOx emission factor NOx EFnew was the only EF that had been measured on board the ship by the Aeromon device. Therefore, all emission factors (EFall) were normalized to the initial NOx EF (Equation (2)). The initial VOC EFinit received from the literature was multiplied by the normalizing factor NF_NOx (Equation (3)):
N F _ N O x = N O x   E F n e w N O x   E F i n i t
N F _ N O x V O C _ E F i n i t = V O C _ E F n e w
The particle size distribution obtained from MAFOR simulations for the plume at the on-shore site for different VOC concentrations was evaluated by comparing against the measured particle size distribution and mean diameter, as described in Equation (1).
To investigate changes in the model performance in connection with the change in VOC emission factor, several sensitivity cases were conducted (Section 3.3, Table 2). This was done because no VOC EF was available for the present plume, and it cannot be taken for granted that the changes in the EF for VOC (compared to the engine testing conditions) are related to changes in the NOx emission factor. The sensitivity cases allowed us to derive a new (real-world) emissions factor for VOCs, VOC EFnew, by multiplication with a correcting factor CF:
V O C _ E F n e w C F = V O C _ E F n e w ,   c o r r
The sensitivity cases are named after the used correction factor (i.e., 0.5 means VOC EFnew ∗ 0.5) within a range from 0.5 to 5.0. With a changing VOC EF, the size distribution of particles changes in the model. With a higher VOC EF, there will be larger particles, which are formed primarily by the growth of smaller particles through condensation. The particle number size distribution computed by MAFOR was compared to the measurement results. The modeled PNSD (particle number size distribution) was divided into seven size classes (described in Section 3.3). The best-fitting CF value was chosen for the size classes S2–S5 by calculating the sum of absolute differences between MAFOR simulations to measurements (Section 3.3, Table 2). These size classes (S2–S5) were chosen since they contain the most sensitive values. The optimal CF is the value where the sum of absolute relative differences is the lowest. The VOC EFnew,corr with the lowest deviation from simulations to measurements serves as a base for comparison with the EF used in STEAM.
The emission factors for VOCs from MAFOR were obtained from measurements on the ship, operating at 50% engine load. Since the aim was to derive a new VOC emission factor for the STEAM emission dataset, the EF used in STEAM for VOCs for vessels at 50% engine load had to be related to the EFs for VOCs that were found in the optimization procedure using MAFOR (i.e., VOC EFnew,corr). It is further assumed that the ratio of the VOC EFnew,corr to the VOC EFSTEAM remains the same for all engine loads. The ratio between the two emission factors was calculated and the STEAM ship emissions for VOCs were adjusted by this factor, f_STEAM, for the CMAQ run:
f _ S T E A M = V O C   E F n e w ,   c o r r V O C   E F S T E A M
The STEAM data were used as ship emissions and for the comparison run with the regional-scale chemical transport model CMAQ, which is described in detail in Section 2.3.

3. Results and Discussion

3.1. Dilution of Gases and Particles between Ship and Remote Monitoring Station

To obtain more representative ship emissions for use in regional air quality models, the emission process from the stack to the background air needs to be considered in detail. Therefore, the emission process, e.g., in-plume chemistry and aerosol dynamics, should be addressed due to the extremely nonlinear reactions among gases and particle coagulation [47].
In addition to particle nucleation, condensation, and coagulation, dry deposition of particles, and gas-phase chemistry inside the plume, the model considered the mixing of air parcels with gases and particles in the background air.
NOx dilution simulated by MAFOR suitably agrees with the observations at the sniffer location and the background concentrations (Figure 4). The dilution curve of NOx serves as a base for the dilution of other chemicals. Particles in the MAFOR model evolve after the model’s starting time of one second after the exhaust leaves the stack.
As described in Section 2.1, the peaks in the NOx concentration and the particle number concentration at the shore-based sniffer station were measured at 17:41:00 UTC (Appendix A). The peak of NOx was followed by a second peak at 17:42:00 UTC. A second peak was also observed in the particle number concentration, with the first peak appearing at 17:41:37 and the second one at 17:42:42 UTC. There is a delay in instrumentation and the instrument times were synchronized using the timing of the peaks. The aerosol measurements were obtained at a high temporal resolution, and the aerosol instrument response time was shorter than that of the NOx sensor. Thus, the exact peak time could be slightly shifted and recorded. The explanation of the two peaks is the meandering behavior of the ship plume. In addition, the ship was moving, and the plume released one minute after the first plume may follow a different path.

3.2. Aerosol Size Distribution

MAFOR was used to calculate the changes in the exhaust particle size distribution over time with increasing downwind distance from the ship stack. MAFOR was adopted to simulate the particle size distribution starting with the input mass concentrations (Table 1). The measured particle background concentrations were observed at the same station as the plume measurements but at a different time. The station was located 530 m distant from the shipping lane. As shown in Figure 5, there is a number concentration peak in the Aitken mode for the measured background particles, indicating a certain influence of shipping-related emissions or other sources of UFPs.
In this study, one peak was found at 15 nm (Figure 5). Regarding the plume at the sniffer location (on-shore mobile lab), only the measurement with the highest concentration was considered; thus, no percentiles are displayed.
The EF for VOCs was adjusted as described in Section 2.2, and the amount of converted H2SO4 from SO2 was set according to the concentrations retrieved from Karl et al. (2020) [47]. With these settings, the simulated plume particle size distribution at a 530 m distance was optimized to narrow the discrepancy with the measured data at the sniffer location (Figure 5, red line). The optimized curve after using the correction factor is displayed in Figure 5, magenta line. This was done to obtain a correction factor by comparing several scenarios of the modeled output against the measured data from the sniffer station in regard to the number concentration and size distribution (Table 2). Only slight variations in the measurements on the ferry and at the mobile laboratory were found, as revealed by the 25th and 75th percentiles of the median of each bin for the measured data at the ship location and at the sniffer location without a plume (Figure 5).
It has been previously shown that the particle number size distribution of fresh diesel exhaust exhibits a bimodal character [47]. The measured size distribution at the sniffer station (Figure 5) exhibits a bimodal structure, indicating that the plume concentration has not yet reached the background concentration. Similar to this study, Pirjola et al. [43] found one dominating peak for the number size distribution in the Aitken mode, with one predominant size distribution in two modes and the dominating mode peaking from 20 nm to 30 nm. Nevertheless, the position of the maximum varied between 30 and 50 nm depending on the type of vessel; thus, comparing the results of this work with those of previous studies is hardly possible.
Karl et al. [47] applied the emission factor for the particle number concentration obtained by Moldanova et al. [45] directly to exhaust flow. However, no data were obtained from measurements on the ship at the same time as those on land, as was done in the present investigation.

3.3. Comparison between Simulated and Measured Data

In the common size class distribution as in Karl et al. [47,52], the integral number concentration and mean diameter of the modeled data were compared to those of the measured data (Table 2).
The best fit was selected based on the differences in size classes S2 to S5. Larger size classes were not used because of the high uncertainty caused by the fluctuations in the effective density of soot particles that could generate a systematic error of approximately 20% in PM estimation using the ELPI [92]. Regarding the small size classes (<15 nm), particle measurements suffer higher uncertainty due to the lower particle charging efficiency or higher diffusion loss [23].
The smallest total difference between the measured and modeled values was found when multiplying the initial VOC EF by a factor of 3.6 (Table 2). The best match for the number concentration and mean diameter could yield different CFs. Thus, the smallest total difference was based on the sum of the absolute difference of the relative deviation for the number concentration plus the sum of the difference of the relative deviation for the mean diameter in size classes S2 to S5 under the respective sensitivity case.
The evolution of the particle size distribution in the ship plume was found to be very sensitive to changes in the VOC EF. The relative differences between the simulated and measured size distribution data could either be positive (simulated values < measured values) or negative (measured values < simulated values).
If there was under- or overestimation of the VOCs, there could be a shift in the peak and a change in the distribution of the particles in the various size ranges. This occurred because all VOC EFs were adjusted, i.e., all volatility classes were adjusted. This also resulted in different deviations (positive or negative) in size classes S1 to S7.

3.4. Validation of New VOC EFnew,corr

The new VOC EF new,corr obtained for the reference plume was tested against 12 additional on-land measurements of plumes from the same vessel, with the assumption that the ship emissions and dilution conditions would remain the same (Table A5, Figures S1–S13). This leads to the possibility to add the confidence interval of the VOC EF, based on the range of the measured and simulated number concentration of each plume. This resulted in a VOC EFnew,corr of 0.323 0.1 + 0.6 g/kg. The difference in number concentration shows that large size classes (S6 and S7) are overestimated by the model, whereas most middle (S3–S5) and small size classes (S1 and S2) are underestimated. The underestimation of nucleation mode particles indicates that the binary homogeneous nucleation process may not efficiently explain new particle formation in ship plumes. Non-volatile particles below 7 nm in size [93] forming in the hot exhaust, so-called core particles, may further grow by the condensation of volatile hydrocarbons and sulfur compounds during cooling [94]. Heterogeneous nucleation onto the solid core particles may play a role [95]. Arnold et al. (2012) [96] have observed condensable dicarboxylic acids in the exhaust of heavy-duty diesel vehicles equipped with modern exhaust after-treatment systems, which can contribute to the atmospheric nucleation processes. Pirjola et al. (2015) [97] investigated the potential of several nucleation mechanisms to predict particle formation in diesel exhaust and found that the best fit with particle size distribution measurements was predicted by a heteromolecular nucleation mechanism in which both sulfuric acid and semi-volatile organic acid molecules participate. However, nucleation is extremely dependent upon dilution conditions in the ship exhaust plume, and the measurement of small particles may be strongly influenced by diffusive losses and sampling artifacts, which makes it impossible to draw firm conclusions on the nucleation processes.
The pattern slightly differs for the mean diameter, where there are more cases of overestimations by the model. In general, the deviations between simulation and measurement for mean diameter are smaller than for number concentration.
It is generally assumed that wind speed has a large influence on dilution, so it was considered in detail for the 13 plumes.
At higher wind speeds (increase of 1 m/s and 2 m/s), the underestimation of the number concentration decreased, but only to a small extent (Figures S1–S13). This indicates that a change in wind speed has only a small impact on the aging in the current investigations. Celik et al. [98] found in a measurement study in the Mediterranean Sea that particle mass emission factors decrease (less condensation) as wind speed increases (stronger dilution), while particle number emission factors increase (less coagulation). They showed that the particle number in the plume was significantly influenced by wind speed. In the presented simulations of this study, changes in wind speed do not translate into different plume dispersion. Therefore, the apparent dependence of particle number concentration on wind speed in [98] is rather caused by different prevailing dilution regimes in the plume dispersion (slow dilution at low wind speed, fast dilution at high wind speed).

3.5. Application of Adjusted VOC Ship Emissions in a Regional-Scale CTM

In the following, the results of the ship emissions used in a regional-scale model are displayed to see whether there is a general improvement in the correlation coefficient as well as normalized mean bias (NMB) between simulated and measured PM2.5 data, when the aerosol modeling and the plume evolution are considered on a small scale.
By utilizing the adjusted VOC emission factor derived from the smallest-difference fit between the modeled and observed particle size distributions (Section 3.3), the ship emissions for the CTM run were modified.
The same adjustment factor (Equation (5)) was used to adjust the VOC emissions in all volatility classes. VOC emissions that have lower volatility and are partly in the condensed phase at standard conditions enter the CMAQ model as particulate organic carbon (POC) emissions. The volatilities of organic compounds in the CMAQ model are calculated internally based on the volatility basis set (VBS) approach [22,74], providing a framework for gas aerosol partitioning and chemical aging of both POAs and SOAs. The CMAQ model contains two basis sets for emitted primary organic aerosols, with each set containing five volatility bins: the first is for the directly emitted POA, while the second is for their multigenerational oxidation products [99,100]. Hence, the model redistributes the POC emission to species of the first volatility base set.
Although the CMAQ model includes many more (nonmethane) VOCs in the CB05 chemical mechanism, regarding the connection with the MAFOR model, only the POCs were modified, including semi- and low-volatile VOCs.
During the preprocessing of CMAQ emissions, the derived STEAM ship emissions were added to the land-based emissions, and the pollutants were assigned to the chemical species of CB05 in the CMAQ model. In this process, the STEAM VOC emissions were multiplied by the adjustment factor f_STEAM with a value of 18.13. This factor was received by dividing the VOC EFnew,corr of 0.323 0.1 + 0.6 g/kg by the VOC EFSTEAM of 0.01779 g/kg fuel. The adjustment factor f_STEAM was applied to the four VOC groups in STEAM. The ship emissions of other substances derived from STEAM, as listed in Section 2.3, remained unchanged.
The CMAQ simulations with adjusted semi-volatile VOC ship emissions were compared to simulations without semi-volatile VOC emissions adjusted to study the effect on PM2.5 concentrations. The increase in semi-volatile VOC ship emissions by a factor of 18.13 caused an increase in the PM2.5 concentration along the major shipping routes by up to 5% (Figure 6). A detailed investigation on the individual aerosol species considered in the CMAQ aerosol module revealed that all showed changing values. The highest increase was found in anthropogenic semi-volatile primary organic compounds, particulate nitrate, and particulate chloride, with an over 40% increase in June. Thus, the scaling of POC emissions in CMAQ leads to particle growth and production of SOAs.
In this study, the highest increase was found in June (Figure 6c), with an increase in the total PM2.5 mass ranging from 3% to 5% over the sea. In the colder months of December and March with less solar radiation, the increase in the PM2.5 concentration varied between 0% and 2% (Figure 6a,d). There were no negative values. This indicates that the impact of VOCs on the background concentration of PM2.5 is relatively limited in the winter and early spring, which is similar to the results of the sensitivity study of Lee et al. [101] and can explained by decreased speed of chemical reactions at low temperatures. Also, emissions from ship traffic might be lower in winter. Absolute values in PM2.5 concentrations were mainly below 10 µg/m3 over water for both CMAQ simulations and indicate only slight differences in the run with adjusted and with original ship emissions (Figures S14 and S15).
The change in PM2.5 daily mean concentrations remained well below 0.3 µg/m3, as revealed by the time series based on data at 28 monitoring stations (Figure A3). The correlation between measured and simulated data from the previous and the new run show that the change in VOC EFs at monitoring stations on land do not cause an increase in correlation coefficients with observed PM2.5. Also, the comparison of the NMB does not show any larger changes (not more than 5%); in all cases, the PM2.5 concentration was underestimated by the model. This is the case for all observed months (March, June, September, December; Figure A3). One might conclude from this that the distance to shipping lanes is either too large, so that the changes in PM2.5 are not observed at measurement stations, or that the effect of an adjusted VOC EF in ship emissions on total PM2.5 concentrations is rather small. Similar results were found in a study carried out in harbor city Marseille, where the ship impact on PM2.5 concentration was almost neglectable, whereas it was drastic for UFP numbers [48].
The overall change in PM2.5 concentration was small along the main shipping routes as well as at several harbor cities (i.e., Barcelona, Genoa, Nador, Oran, Algiers). Benchrif et al. [102] investigated PM sources in Moroccan harbor cities and found that the largest influence of PM in the Western Mediterranean Sea on air quality comes from maritime traffic emissions. In areas with increased PM2.5 concentration, they can have effects on radiative forcing and cloud formation.

4. Conclusions

The aerosol box model MAFOR was used to calculate the particle number and mass evolution in fresh ship plumes based on observed data received from a measurement campaign targeting particle emissions from a regular passenger ship. The model was validated by conducting several sensitivity runs with changes in the VOC emission factor to get the best approximation from modeled to observed particle size distribution data. The best-fitting approximation was used for correcting the emission factor of VOCs for ship engines. For the validation of the corrected VOC EF, 12 ship plumes monitored at the mobile lab on land were evaluated. This was done based on the assumption that each plume was homogenous and weather conditions were comparable during all 12 plumes that were measured on shore. With this, a confidence interval for the new VOC EF was calculated that indicates the range of the VOC EF and the variability that can occur due to the different conditions in the plumes and the environmental conditions. The newly determined VOC EF is confined by a confidence interval from −44% to 192%. The integration of the output obtained with the MAFOR model into the regional-scale model CMAQ was done via the ship emission dataset from STEAM by using the new VOC emission factor. Two runs with CMAQ were simulated, one with changes in ship emissions and one with the initial ship emissions. Results indicate that there exist effects of changes in VOC ship emissions on ambient PM2.5 concentration, but the overall changes in the concentration were rather small relative to the increase in VOC ship emissions. This low relative increase might be due to low levels of pre-existing particulate organic matter, which were present in the CMAQ simulations. This leads to generally lower particle formation. In the present study, only the lowest atmospheric layer was investigated. The concentration and amount of organic mass could be different in higher layers and deliver changed results.
Uncertainties in meteorology can also have an impact on the CMAQ simulations and the simulated concentrations, as already discussed in Fink et al. [38,61]. Even small changes in meteorological conditions (e.g., wind speed, humidity) might impact the particle formation processes.
Direct dilution of ship plume particles within the grid box volume of a regional-scale chemistry transport model leads to the underestimation of ambient particulate matter concentrations at coastal sites. This problem should be addressed in the future by integrating a sub-grid aerosol box model for the treatment of ship plumes in regional-scale chemistry transport models.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/toxics12060432/s1, Figure S1: Ship plume at 21.07.2021 17:31:10 UTC with changing wind speed. Figure S2: Ship plume at 02.08.2021 17:22:46 UTC with changing wind speed. Figure S3: Ship plume at 30.08.2021 17:34:00 UTC with changing wind speed. Figure S4: Ship plume at 01.09.2021 17:30:10 UTC with changing wind speed. Figure S5: Ship plume at 03.09.2021 17:41:40 UTC with changing wind speed. Figure S6: Ship plume at 08.07.2022 17:22:10 UTC with changing wind speed. Figure S7: Ship plume at 24.07.2022 16:28:20 UTC with changing wind speed. Figure S8: Ship plume at 26.07.2022 17:31:04 UTC with changing wind speed. Figure S9: Ship plume at 08.10.2022 16:29:49 UTC with changing wind speed. Figure S10: Ship plume at 19.07.2023 17:32:54 UTC with changing wind speed. Figure S11: Ship plume at 29.07.2023 16:40:50 UTC with changing wind speed. Figure S12: Ship plume at 26.08.2023 16:25:00 UTC with changing wind speed. Figure S13: Ship plume at 28.08.2023 17:28:10 UTC with changing wind speed. Figure S14: PM2.5 old absolute concentration, (a) mean value in March, (b) mean value in June, (c) mean value in September, (d) mean value in December. Figure S15: PM2.5 new absolute concentration, (a) mean value in March, (b) mean value in June, (c) mean value in September, (d) mean value in December.

Author Contributions

Conceptualization, L.F., M.K. and V.M.; methodology, L.F. and M.K.; software, M.K.; validation, L.F. and M.K.; formal analysis, L.F., M.K., V.M. and A.W.; resources, A.W., P.S. and M.I.; data curation, L.F.; writing—original draft preparation, L.F.; writing—review and editing, M.K., V.M., A.W., P.S. and M.I.; visualization, L.F.; supervision, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union’s Horizon 2020 research and innovation program under grant agreement No. 814893 (the project “Shipping Contributions to Inland Pollution Push for the Enforcement of Regulations”, SCIPPER).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The model output data and the measurement data used in this study will be made available upon request.

Acknowledgments

The authors are very grateful to Joni Kalliokoski and Lassi Markkula from Tampere University, who participated in installing and operating the instruments. Furthermore, we would like to thank everyone who participated in the measurement campaign, set up instrumentation, and provided data for the model runs.

Conflicts of Interest

Author Matti Irjala was employed by the company Aeromon Oy. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Figure A1. Total particle count measured during the SCIPPER measurement campaign on 03.09.2021 in Laboe at the Kieler Förde (54.392809, 10.209186). The plume arrived at the measurement station at 17:41:37 UTC and the plume peak was determined when the total concentration was highest. The uncertainty for the total particle number is approximately ±10 %, as reported for the German Ultrafine Aerosol Network [103].
Figure A1. Total particle count measured during the SCIPPER measurement campaign on 03.09.2021 in Laboe at the Kieler Förde (54.392809, 10.209186). The plume arrived at the measurement station at 17:41:37 UTC and the plume peak was determined when the total concentration was highest. The uncertainty for the total particle number is approximately ±10 %, as reported for the German Ultrafine Aerosol Network [103].
Toxics 12 00432 g0a1
Figure A2. NOx concentrations after the Germanica passed the measurement station. Gases and meteorological parameters were measured by an airpointer 4D measurement system for NOx, O3, SO2, and a weather station during the SCIPPER measurement campaign on 03.09.2021 in Laboe at the Kieler Förde (54.392809, 10.209186). The uncertainty for NOx measurements is given by the manufacturer and is 2 µg/m3 (1 ppb) for measured concentrations <500 ppb.
Figure A2. NOx concentrations after the Germanica passed the measurement station. Gases and meteorological parameters were measured by an airpointer 4D measurement system for NOx, O3, SO2, and a weather station during the SCIPPER measurement campaign on 03.09.2021 in Laboe at the Kieler Förde (54.392809, 10.209186). The uncertainty for NOx measurements is given by the manufacturer and is 2 µg/m3 (1 ppb) for measured concentrations <500 ppb.
Toxics 12 00432 g0a2

Appendix B

Table A1. Overview of the evaluated plumes. WD = wind direction (in degree), SD = shipping direction (in degree), WSP = wind speed (in m/s), SSP = ship speed (in m/s), Vector = calculated wind speed vector over the deck (in m/s), Travel = travel time of the plume from the ship to the sniffer (in s), RH = relative humidity (in %), Max. Meas. dN/dlogDp = maximum of measured dN/dlogDp, Max. Sim. dN/dlogDp = maximum of simulated dN/dlogDp in 530 m.
Table A1. Overview of the evaluated plumes. WD = wind direction (in degree), SD = shipping direction (in degree), WSP = wind speed (in m/s), SSP = ship speed (in m/s), Vector = calculated wind speed vector over the deck (in m/s), Travel = travel time of the plume from the ship to the sniffer (in s), RH = relative humidity (in %), Max. Meas. dN/dlogDp = maximum of measured dN/dlogDp, Max. Sim. dN/dlogDp = maximum of simulated dN/dlogDp in 530 m.
DateTimeWDSDWDPSSPVectorTravel RHMax. Meas. dN/dlogDpMax. Sim. dN/dlogDp
21.07.202117:31:10263 ± 1614.8 ± 0.23.0 ± 0.66.3 ± 1.96.8177.8562.6 ± 1.64757.848196.38
02.08.202117:22:40268. ± 1217.0 ± 0.24.9 ± 1.16.7 ± 1.48.1964.7172.1 ± 1.43595.692540.21
30.08.202117:34:0013. ± 1937.9 ± 0.15.1 ± 0.86.6 ± 1.18.2064.6183.3 ± 1.610,054.932009.83
01.09.202117:30:10290 ± 1218.6 ± 0.23.2 ± 0.56.7 ± 1.87.3771.9072.2 ± 2.319,417.267751.90
03.09.202117:41:40247 ± 1018.3 ± 0.13.3 ± 0.65.4 ± 0.86.1885.7074.0 ± 1.310,283.3812,648.08
08.07.202217:22:10253 ± 1316.1 ± 0.35.6 ± 0.96.6 ± 1.18.6861.0568.9 ± 1.25539.442194.99
24.07.202216:28:20277 ± 2217.7 ± 0.12.7 ± 0.96.2 ± 1.26.8477.5159.6 ± 2.87250.3211,694.63
26.07.202217:30:30276 ± 1013.1 ± 0.15.8 ± 1.26.1 ± 2.28.3363.6057.3 ± 1.27350.022191.53
08.10.202216:29:30258 ± 1315.2 ± 0.24.8 ± 0.96.4 ± 2.18.0366.0471.7 ± 2.310,077.736174.34
19.07.202317:32:10254 ± 1517.4 ± 0.12.5 ± 0.55.3 ± 2.15.6993.0778.8 ± 1.82941.2911,420.58
29.07.202316:40:50267 ± 1018.8 ± 0.25.1 ± 0.94.5 ± 2.46.8976.9170.1 ± 1.911,034.411453.65
26.08.202316:25:30232 ± 914.6 ± 0.26.9 ± 0.86.0 ± 1.19.2457.3375.2 ± 1.46366.785407.44
28.08.202317:28:30293 ± 1317.5 ± 0.12.9 ± 0.75.9 ± 1.26.4682.0965.9 ± 1.69805.028832.01

Appendix C

Table A2. Comparison of plume dispersion in MAFOR and Chosson et al. [67]. Unit for dil_rateM is s−1 and for dil_rateC is min−1.
Table A2. Comparison of plume dispersion in MAFOR and Chosson et al. [67]. Unit for dil_rateM is s−1 and for dil_rateC is min−1.
MAFORChosson et al. (2008) [67]
CalculationParameters (1)CalculationParameters (2)
d i l _ r a t e M = ( a d i l _ t i m e ) b a′ = 1.659 d i l _ r a t e C = a ( t * t ) b a = 0.051 [min−1]
b = 1.133b = 1.08
dil_time = input:
time passed in plume [s]
t* = turn-over
time scale [min]
t = time passed in plume [min]
(1) Based on adjustments from Moldanova, J. (personal communication). (2) Parameters are for a buoyancy flux of 250 m4/s3.
Table A3. Plume height and width approximation following the formulation of von Glasow et al. [68].
Table A3. Plume height and width approximation following the formulation of von Glasow et al. [68].
Plume WidthPlume Height
w p l t = w 0 ( t t 0 ) α wpl = plume width h p l t = h 0 ( t t 0 ) β hpl = plume height
t = plume time
t0 = reference time after plume release
w0 = reference dimension of the plume at time t0
α = plume expansion rates in the horizontal
t = plume time
t0 = reference time after plume release
h0 = reference dimension of the plume at time t0
β = plume expansion rates in the vertical

Appendix D

Table A4. Procedure and intermediate steps to obtain the final input and starting concentrations for MAFOR from the emission factors of the individual substances. First-guess emission factors were received from the literature. The emission factor was multiplied by a changing factor of 0.67 to normalize to the emission factor of NOx, which was measured by Aeromon measurement device on board. Substance split was taken from the literature; for NOx it was assumed that it consists of 90% NO.
Table A4. Procedure and intermediate steps to obtain the final input and starting concentrations for MAFOR from the emission factors of the individual substances. First-guess emission factors were received from the literature. The emission factor was multiplied by a changing factor of 0.67 to normalize to the emission factor of NOx, which was measured by Aeromon measurement device on board. Substance split was taken from the literature; for NOx it was assumed that it consists of 90% NO.
SubstanceFirst-Guess Emission Factor (g/kg_fuel)Emission Factor (g/kg_fuel) Received with Changing Factor 0.67Diluted Gases (g/m3) (1)Substance SplitRecalculation to cm3MAFOR Input after Dilution of 1/8 (in cm3)MAFOR Input VOCs without Dilution of 1/8 for VOCs (in cm3) (6)
NOx40.91 (1)27.480.11NO2.041 × 10152.55 × 1014
NO22.268 × 10142.83 × 1013
SOx0.34 (2)0.230.00093SO28.74 × 10121.09 × 1012
H2SO41 × 1010 (4)1.25 × 109
SO38.74 × 108 (5)1.09 × 108
CO23.64 (1)15.880.065 1.40 × 10151.75 × 1014
HCHO0.0022 (1)0.00156.15 × 10−6 1.23 × 10111.54 × 1010
TOL0.0078 (3)0.00522.15 × 10−5 1.41 × 10111.76 × 1010
XYL0.00295 (3)0.0028.15 × 10−6 4.63 × 10105.78 × 109
TMB0.0028 (3)0.00197.68 × 10−6 3.85 × 10104.82 × 109
IVOC0.0076 (4)0.0187.56 × 10−5 2.53 × 1011 4.21 × 1011
SVOC0.0105 (4)0.0250.0001 2.51 × 10114.19 × 1011
ELVOC0.105 (4)0.250.001 1.96 × 10123.27 × 1012
(1) Based on Fridell et al., 2021 [72] (corresponding fuel mass flow: 0.00411 kg_fuel/m3). (2) Based on FSC limits. (3) Based on Wu et al., 2020 [53]. (4) Based on Karl et al., 2020 [47]. (5) Based on the assumption of 1% SO3. (6) Based on the concentration in cm−3 from Karl et al. 2020, which is already diluted; initial value is multiplied by 6 according to Wu et al., 2020 [53].

Appendix E

Table A5. Size ranges of the size classes; S1: 1–10 nm; S2: 10–20 nm; S3: 20–50 nm; S4: 50–100 nm; S5: 100–300 nm; S6: 300–600 nm; S7: >600 nm; differences are referring to the deviations in the modeled and measured data. A CF (correcting factor) value of 3.6 was used in all simulations.
Table A5. Size ranges of the size classes; S1: 1–10 nm; S2: 10–20 nm; S3: 20–50 nm; S4: 50–100 nm; S5: 100–300 nm; S6: 300–600 nm; S7: >600 nm; differences are referring to the deviations in the modeled and measured data. A CF (correcting factor) value of 3.6 was used in all simulations.
Difference in Number Concentrations Compared to Measurements (in %)
DateS1S2S3S4S5S6S7|S2| + |S3| +
|S4| + |S5|
21.07.202169.9562.1666.19−23.82−1.1147.2495.05153
02.08.202172.6252.49−77.02−174.3−71.1870.1−29.35375
30.08.2021−9.8813.91−192.51−2.88−3.775.371.62213
01.09.2021−0.3318.670.25−185.7−60.7856.999.2265.4
03.09.2021−149.343.714.915.670.0578.3983.5314.35
08.07.202246.0326.9−234.71−398.34−156.8942.8590.65816.83
24.07.202210.3347.6766.42−25.36−83.8380.998.2223
26.07.2022−42.65−51.75−177.08−92.6−166.390.4782.84487.8
08.10.202254.2244.05−114.64−144.76−54.5225.6393.64357.96
19.07.202390.8570.4686.889.0818.3323.397184.8
29.07.2023−74.36−167.9−704.88−146.31−40.31−46.0363.71059.4
26.08.202358.8955.01−27.51−25.24−57.7648.3291.98165.52
28.08.2023−4.94−0.9845.46−19.47−47.3356.6294.03113.25
Difference in Mean Diameter Compared to Measurements (in %)
DateS1S2S3S4S5S6S7|S2| + |S3| +
|S4| + |S5|
21.07.2021−25.448.07−14.7714.15−13.547.223.3150.5
02.08.2021−11.88−2.3−22.0114.71−4.434.461.0243.43
30.08.2021−42.84−4.47−16.1314.25−4.84−5.87−13.6439.7
01.09.2021−27.026.17−15.814.35−7.577.7−2.5543.89
03.09.20219.19−5.9610.25−12.257.63−3.426.1113.93
08.07.2022−14.22−2.37−22.8818.31−13.4516.75−3.0257.01
24.07.2022−10.48.86−123.54−1.314.43−8.9925.7
26.07.2022−25.941.31−14.2612.43−10.764.41−10.2638.8
08.10.2022−16.65−2.54−20.7317.50−9.252.294.5750.02
19.07.2023−30.1912.03−13.416.74−10.21−2.912.6842.38
29.07.2023−16.6−3.00−14.7314.15−4.77−3.873.2236.63
26.08.2023−25.28−0.95−16.172.84−3.819.603.7223.77
28.08.2023−25.039.46−9.569.42−7.289.873.9235.72

Appendix F

Figure A3. Time series at 28 monitoring stations over the whole domain; stations are described in Fink et al. [51]. (a) = March 2015, (b) = June 2015, (c) = September 2015, (d) = December 2015. Rold is the correlation for the PM2.5 concentration between measured and simulated concentrations for the run without the adjusted VOC emissions; Rnew is for the run with the adjusted VOC emissions. NMBold is the normalized mean bias for the PM2.5 concentration between measured and simulated concentrations for the run without the adjusted VOC emissions; NMBnew is for the run with the adjusted VOC emissions.
Figure A3. Time series at 28 monitoring stations over the whole domain; stations are described in Fink et al. [51]. (a) = March 2015, (b) = June 2015, (c) = September 2015, (d) = December 2015. Rold is the correlation for the PM2.5 concentration between measured and simulated concentrations for the run without the adjusted VOC emissions; Rnew is for the run with the adjusted VOC emissions. NMBold is the normalized mean bias for the PM2.5 concentration between measured and simulated concentrations for the run without the adjusted VOC emissions; NMBnew is for the run with the adjusted VOC emissions.
Toxics 12 00432 g0a3

References

  1. Wang, C.; Corbett, J.J.; Firestone, J. Improving spatial representation of global ship emissions inventories. Environ. Sci. Technol. 2008, 42, 193–199. [Google Scholar] [CrossRef]
  2. Sofiev, M.; Winebrake, J.J.; Johansson, L.; Carr, E.W.; Prank, M.; Soares, J.; Vira, J.; Kouznetsov, R.; Jalkanen, J.-P.; Corbett, J.J. Cleaner fuels for ships provide public health benefits with climate tradeoffs. Nat. Commun. 2018, 9, 406. [Google Scholar] [CrossRef] [PubMed]
  3. Ausmeel, S.; Eriksson, A.; Ahlberg, E.; Kristensson, A. Methods for identifying aged ship plumes and estimating contribution to aerosol exposure downwind of shipping lanes. Atmos. Meas. Tech. 2019, 12, 4479–4493. [Google Scholar] [CrossRef]
  4. Endresen, Ø.; Sørgård, E.; Sundet, J.K.; Dalsøren, S.B.; Isaksen, I.S.A.; Berglen, T.F.; Gravir, G. Emission from international sea transportation and environmental impact. J. Geophys. Res. 2003, 108, D17. [Google Scholar] [CrossRef]
  5. Klimont, Z.; Kupiainen, K.; Heyes, C.; Purohit, P.; Cofala, J.; Rafaj, P.; Borken-Kleefeld, J.; Schöpp, W. Global anthropogenic emissions of particulate matter including black carbon. Atmos. Chem. Phys. 2017, 17, 8681–8723. [Google Scholar] [CrossRef]
  6. Heusinkveld, H.J.; Wahle, T.; Campbell, A.; Westerink, R.H.S.; Tran, L.; Johnston, H.; Stone, V.; Cassee, F.R.; Schins, R.P.F. Neurodegenerative and neurological disorders by small inhaled particles. Neurotoxicology 2016, 56 (Suppl. C), 94–106. [Google Scholar] [CrossRef]
  7. Chen, R.; Hu, B.; Liu, Y.; Xu, J.; Yang, G.; Xu, D.; Chen, C. Beyond PM2.5: The role of ultrafine particles on adverse health effects of air pollution. Biochim. Biophys. Acta (BBA)-Gen. Subj. 2016, 1860, 2844–2855. [Google Scholar] [CrossRef] [PubMed]
  8. Gao, R.; Sang, N. Quasi-ultrafine particles promote cell metastasis via HMGB1-mediated cancer cell adhesion. Environ. Pollut. 2020, 256, 113390. [Google Scholar] [CrossRef]
  9. Callén, M.S.; Iturmendi, A.; López, J.M. Source apportionment of atmospheric PM2.5-bound polycyclic aromatic hydrocarbons by a PMF receptor model. Assessment of potential risk for human health. Environ. Pollut. 2014, 195, 167–177. [Google Scholar] [CrossRef] [PubMed]
  10. Khaniabadi, Y.O.; Polosa, R.; Chuturkova, R.Z.; Daryanoosh, M.; Goudarzi, G.; Borgini, A.; Tittarelli, A.; Basiri, H.; Armin, H.; Nourmoradi, H.; et al. Human health risk assessment due to ambient PM10 and SO2 by an air quality modeling technique. Process Saf. Environ. Prot. 2017, 111, 346–354. [Google Scholar] [CrossRef]
  11. Amoatey, P.; Omidvarborna, H.; Baawain, M. The modeling and health risk assessment of PM2.5 from Tema Oil Refinery. Hum. Ecol. Risk Assess. Int. J. 2018, 24, 1181–1196. [Google Scholar] [CrossRef]
  12. Canu, I.G.; Bateson, T.F.; Bouvard, V.; Debia, M.; Dion, C.; Savolainen, K.; Yu, I.-J. Human exposure to carbon-based fibrous nanomaterials: A review. Int. J. Hyg. Environ. Health 2016, 219, 166–175. [Google Scholar] [CrossRef] [PubMed]
  13. Pope, C.A.; Burnett, R.T.; Turner, M.C.; Cohen, A.; Krewski, D.; Jerrett, M.; Gapstur, S.M.; Thun, M.J. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: Shape of the exposure–response relationships. Environ. Health Perspect. 2011, 119, 1616–1621. [Google Scholar] [CrossRef] [PubMed]
  14. Burnett, R.T.; Pope, C.A.; Ezzati, M.; Olives, C.; Lim, S.S.; Mehta, S.; Shin, H.H.; Singh, G.; Hubbell, B.; Brauer, M.; et al. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ. Health Perspect. 2014, 122, 397–403. [Google Scholar] [CrossRef] [PubMed]
  15. Apte, J.S.; Marshall, J.D.; Cohen, A.J.; Brauer, M. Addressing Global Mortality from Ambient PM2.5. Environ. Sci. Technol. 2015, 49, 8057–8066. [Google Scholar] [CrossRef] [PubMed]
  16. Karamfilova, E. BRIEFING Implementation Appraisal, Revision of the EU Ambient Air Quality Directives, EPRS European Parliamentary Research Service. 2022. Available online: https://www.europarl.europa.eu/RegData/etudes/BRIE/2022/734679/EPRS_BRI(2022)734679_EN.pdf (accessed on 1 August 2023).
  17. Atkinson, R. Atmospheric chemistry of VOCs and NOx. Atmos. Environ. 2000, 34, 2063–2101. [Google Scholar] [CrossRef]
  18. Seinfeld, J.H.; Pandis, S.N. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change; John Wiley & Sons.: Hoboken, NJ, USA, 2006; ISBN 9780471720188. [Google Scholar]
  19. Yli-Juuti, T.; Mohr, C.; Riipinen, I. Open questions on atmospheric nanoparticle growth. Commun. Chem. 2020, 3, 106. [Google Scholar] [CrossRef] [PubMed]
  20. Gómez-Moreno, F.J.; Alonso-Blanco, E.; Díaz, E.; Coz, E.; Molero, F.; Núñez, L.; Palacios, M.; Barreiro, M.; Fernández, J.; Salvador, P.; et al. On the influence of VOCs on new particle growth in a Continental-Mediterranean region. Environ. Res. Commun. 2022, 4, 125010. [Google Scholar] [CrossRef]
  21. Lu, Q.; Zhao, Y.; Robinson, A.L. Comprehensive organic emission profiles for gasoline, diesel, and gas-turbine engines including intermediate and semi-volatile organic compound emissions. Atmos. Chem. Phys. 2018, 18, 17637–17654. [Google Scholar] [CrossRef]
  22. Robinson, A.L.; Donahue, N.M.; Shrivastava, M.K.; Weitkamp, E.A.; Sage, A.M.; Grieshop, A.P.; Lane, T.E.; Pierce, J.R.; Pandis, S.N. Rethinking organic aerosols: Semivolatile emissions and photochemical aging. Science 2007, 315, 1259–1262. [Google Scholar] [CrossRef]
  23. Wang, X.; Grose, M.A.; Caldow, R.; Osmondson, B.L.; Swanson, J.J.; Chow, J.C.; Watson, J.G.; Kittelson, D.B.; Li, Y.; Xue, J.; et al. Improvement of Engine Exhaust Particle Sizer (EEPS) size distribution measurement—II. Engine exhaust particles. J. Aerosol Sci. 2016, 92, 83–94. [Google Scholar] [CrossRef]
  24. Guo, S.; Hu, M.; Peng, J.; Wu, Z.; Zamora, M.L.; Shang, D.; Du, Z.; Zheng, J.; Fang, X.; Tang, R.; et al. Remarkable nucleation and growth of ultrafine particles from vehicular exhaust. Proc. Natl. Acad. Sci. USA 2020, 117, 3427–3432. [Google Scholar] [CrossRef]
  25. Hallquist, M.; Wenger, J.C.; Baltensperger, U.; Rudich, Y.; Simpson, D.; Claeys, M.; Dommen, J.; Donahue, N.M.; George, C.; Goldstein, A.H.; et al. The formation, properties and impact of secondary organic aerosol: Current and emerging issues. Atmos. Chem. Phys. 2009, 9, 5155–5236. [Google Scholar] [CrossRef]
  26. Almeida, J.; Schobesberger, S.; Kürten, A.; Ortega, I.K.; Kupiainen-Määttä, O.; Praplan, A.P.; Adamov, A.; Amorim, A.; Bianchi, F.; Breitenlechner, M.; et al. Molecular understanding of sulphuric acid–amine particle nucleation in the atmosphere. Nature 2013, 502, 359–363. [Google Scholar] [CrossRef] [PubMed]
  27. Kulmala, M.; Kontkanen, J.; Junninen, H.; Lehtipalo, K.; Manninen, H.E.; Nieminen, T.; Petäjä, T.; Sipilä, M.; Schobesberger, S.; Rantala, P.; et al. Direct observations of atmospheric aerosol nucleation. Science 2013, 339, 943–946. [Google Scholar] [CrossRef]
  28. Pusfitasari, E.D.; Ruiz-Jimenez, J.; Tiusanen, A.; Suuronen, M.; Haataja, J.; Wu, Y.; Kangasluoma, J.; Luoma, K.; Petäjä, T.; Jussila, M.; et al. Vertical profiles of volatile organic compounds and fine particles in atmospheric air by using an aerial drone with miniaturized samplers and portable devices. Atmos. Chem. Phys. 2023, 23, 5885–5904. [Google Scholar] [CrossRef]
  29. Ge, X.; Wexler, A.S.; Clegg, S.L. Atmospheric amines—Part I. A review. Atmos. Environ. 2011, 45, 524–546. [Google Scholar] [CrossRef]
  30. Jimenez, J.L.; Canagaratna, M.R.; Donahue, N.M.; Prevot, A.S.H.; Zhang, Q.; Kroll, J.H.; Decarlo, P.F.; Allan, J.D.; Coe, H.; Ng, N.L.; et al. Evolution of organic aerosols in the atmosphere. Science 2009, 326, 1525–1529. [Google Scholar] [CrossRef]
  31. Guo, S.; Hu, M.; Zamora, M.L.; Peng, J.; Shang, D.; Zheng, J.; Du, Z.; Wu, Z.; Shao, M.; Zeng, L.; et al. Elucidating severe urban haze formation in China. Proc. Natl. Acad. Sci. USA 2014, 111, 17373–17378. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, T.; Liu, J.; Ma, Q.; Chu, B.; Zhang, P.; Ma, J.; Liu, Y.; Zhong, C.; Liu, P.; Wang, Y.; et al. Measurement report: Effects of photochemical aging on the formation and evolution of summertime secondary aerosol in Beijing. Atmos. Chem. Phys. 2021, 21, 1341–1356. [Google Scholar] [CrossRef]
  33. Lee, P.K.H.; Brook, J.R.; Dabek-Zlotorzynska, E.; Mabury, S.A. Identification of the major sources contributing to PM2.5 observed in Toronto. Environ. Sci. Technol. 2003, 37, 4831–4840. [Google Scholar] [CrossRef] [PubMed]
  34. Debevec, C.; Sauvage, S.; Gros, V.; Sellegri, K.; Sciare, J.; Pikridas, M.; Stavroulas, I.; Leonardis, T.; Gaudion, V.; Depelchin, L.; et al. Driving parameters of biogenic volatile organic compounds and consequences on new particle formation observed at an eastern Mediterranean background site. Atmos. Chem. Phys. 2018, 18, 14297–14325. [Google Scholar] [CrossRef]
  35. Karl, M.; Bieser, J.; Geyer, B.; Matthias, V.; Jalkanen, J.-P.; Johansson, L.; 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. 2019, 19, 1721–1752. [Google Scholar] [CrossRef]
  36. Aksoyoglu, S.; Baltensperger, U.; Prévôt, A.S.H. Contribution of ship emissions to the concentration and deposition of air pollutants in Europe. Atmos. Chem. Phys. 2016, 16, 1895–1906. [Google Scholar] [CrossRef]
  37. Nunes, R.A.O.; Alvim-Ferraz, M.C.M.; Martins, F.G.; Calderay-Cayetano, F.; Durán-Grados, V.; Moreno-Gutiérrez, J.; Jalkanen, J.-P.; Hannuniemi, H.; Sousa, S.I.V. Shipping emissions in the Iberian Peninsula and the impacts on air quality. Atmos. Chem. Phys. 2020, 20, 9473–9489. [Google Scholar] [CrossRef]
  38. Fink, L.; Karl, M.; Matthias, V.; Oppo, S.; Kranenburg, R.; Kuenen, J.; Jutterström, S.; Moldanova, J.; Majamäki, E.; Jalkanen, J.-P. A multimodel evaluation of the potential impact of shipping on particle species in the Mediterranean Sea. Atmos. Chem. Phys. 2023, 23, 10163–10189. [Google Scholar] [CrossRef]
  39. Im, U.; Bianconi, R.; Solazzo, E.; Kioutsioukis, I.; Badia, A.; Balzarini, A.; Baró, R.; Bellasio, R.; Brunner, D.; Chemel, C.; et al. Evaluation of operational online-coupled regional air quality models over Europe and North America in the context of AQMEII phase 2. Part II: Particulate matter. Atmos. Environ. 2015, 115, 421–441. [Google Scholar] [CrossRef]
  40. Frohn, L.M.; Ketzel, M.; Christensen, J.H.; Brandt, J.; Im, U.; Massling, A.; Andersen, C.; Plejdrup, M.S.; Nielsen, O.-K.; van der Gon, H.D.; et al. Modelling ultrafine particle number concentrations at address resolution in Denmark from 1979–2018—Part 1: Regional and urban scale modelling and evaluation. Atmos. Environ. 2021, 264, 118631. [Google Scholar] [CrossRef]
  41. Intergovernmental Panel on Climate Change. Climate Change 2013—The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change: Clouds and Aerosols; Cambridge University Press: Cambridge, UK, 2014; ISBN 9781107415324. [Google Scholar] [CrossRef]
  42. González, Y.; Rodríguez, S.; García, J.C.G.; Trujillo, J.L.; García, R. Ultrafine particles pollution in urban coastal air due to ship emissions. Atmos. Environ. 2011, 45, 4907–4914. [Google Scholar] [CrossRef]
  43. Pirjola, L.; Pajunoja, A.; Walden, J.; Jalkanen, J.-P.; Rönkkö, T.; Kousa, A.; Koskentalo, T. Mobile measurements of ship emissions in two harbour areas in Finland. Atmos. Meas. Tech. 2014, 7, 149–161. [Google Scholar] [CrossRef]
  44. Vouitsis, I.; Portugal, J.; Kontses, A.; Karlsson, H.L.; Faria, M.; Elihn, K.; Juárez-Facio, A.T.; Amato, F.; Piña, B.; Samaras, Z. Transport-related airborne nanoparticles: Sources, different aerosol modes, and their toxicity. Atmos. Environ. 2023, 301, 119698. [Google Scholar] [CrossRef]
  45. Moldanova, J.; Fridell, E.; Winnes, H.; Holmin-Fridell, S.; Boman, J.; Jedynska, A.; Tishkova, V.; Demirdjian, B.; Joulie, S.; Bladt, H.; et al. Physical and chemical characterisation of PM emissions from two ships operating in European Emission Control Areas. Atmos. Meas. Tech. 2013, 6, 3577–3596. [Google Scholar] [CrossRef]
  46. Jonsson, Å.M.; Westerlund, J.; Hallquist, M. Size-resolved particle emission factors for individual ships. Geophys. Res. Lett. 2011, 38, L13809. [Google Scholar] [CrossRef]
  47. Karl, M.; Pirjola, L.; Karppinen, A.; Jalkanen, J.-P.; Ramacher, M.O.P.; Kukkonen, J. Modeling of the Concentrations of Ultrafine Particles in the Plumes of Ships in the Vicinity of Major Harbors. Int. J. Environ. Res. Public Health 2020, 17, 777. [Google Scholar] [CrossRef] [PubMed]
  48. Karl, M.; Ramacher, M.O.P.; Oppo, S.; Lanzi, L.; Majamäki, E.; Jalkanen, J.-P.; Lanzafame, G.M.; Temime-Roussel, B.; Le Berre, L.; D’anna, B. Measurement and Modeling of Ship-Related Ultrafine Particles and Secondary Organic Aerosols in a Mediterranean Port City. Toxics 2023, 11, 771. [Google Scholar] [CrossRef]
  49. Huszar, P.; Cariolle, D.; Paoli, R.; Halenka, T.; Belda, M.; Schlager, H.; Miksovsky, J.; Pisoft, P. Modeling the regional impact of ship emissions on NOx and ozone levels over the Eastern Atlantic and Western Europe using ship plume parameterization. Atmos. Chem. Phys. 2010, 10, 6645–6660. [Google Scholar] [CrossRef]
  50. Kim, H.S.; Kim, Y.H.; Han, K.M.; Kim, J.; Song, C.H. Ozone production efficiency of a ship-plume: ITCT 2K2 case study. Chemosphere 2016, 143, 17–23. [Google Scholar] [CrossRef] [PubMed]
  51. Murena, F.; Mocerino, L.; Quaranta, F.; Toscano, D. Impact on air quality of cruise ship emissions in Naples, Italy. Atmos. Environ. 2018, 187, 70–83. [Google Scholar] [CrossRef]
  52. Karl, M.; Pirjola, L.; Grönholm, T.; Kurppa, M.; Anand, S.; Zhang, X.; Held, A.; Sander, R.; Dal Maso, M.; Topping, D.; et al. Description and evaluation of the community aerosol dynamics model MAFOR v2.0. Geosci. Model Dev. 2022, 15, 3969–4026. [Google Scholar] [CrossRef]
  53. Wu, Z.; Zhang, Y.; He, J.; Chen, H.; Huang, X.; Wang, Y.; Yu, X.; Yang, W.; Zhang, R.; Zhu, M.; et al. Dramatic increase in reactive volatile organic compound (VOC) emissions from ships at berth after implementing the fuel switch policy in the Pearl River Delta Emission Control Area. Atmos. Chem. Phys. 2020, 20, 1887–1900. [Google Scholar] [CrossRef]
  54. Xiao, Q.; Li, M.; Liu, H.; Fu, M.; Deng, F.; Lv, Z.; Man, H.; Jin, X.; Liu, S.; He, K. Characteristics of marine shipping emissions at berth: Profiles for particulate matter and volatile organic compounds. Atmos. Chem. Phys. 2018, 18, 9527–9545. [Google Scholar] [CrossRef]
  55. Karl, M.; Gross, A.; Pirjola, L.; Leck, C. A new flexible multicomponent model for the study of aerosol dynamics in the marine boundary layer. Tellus B Chem. Phys. Meteorol. 2011, 63, 1001–1025. [Google Scholar] [CrossRef]
  56. Jalkanen, J.-P.; Brink, A.; Kalli, J.; Pettersson, H.; Kukkonen, J.; Stipa, T. A modelling system for the exhaust emissions of marine traffic and its application in the Baltic Sea area. Atmos. Chem. Phys. 2009, 9, 9209–9223. [Google Scholar] [CrossRef]
  57. Jalkanen, J.-P.; Johansson, L.; Kukkonen, J.; Brink, A.; Kalli, J.; Stipa, T. Extension of an assessment model of ship traffic exhaust emissions for particulate matter and carbon monoxide. Atmos. Chem. Phys. 2012, 12, 2641–2659. [Google Scholar] [CrossRef]
  58. Johansson, L.; Jalkanen, J.-P.; Kalli, J.; Kukkonen, J. The evolution of shipping emissions and the costs of regulation changes in the northern EU area. Atmos. Chem. Phys. 2013, 13, 11375–11389. [Google Scholar] [CrossRef]
  59. Johansson, L.; Jalkanen, J.-P.; Kukkonen, J. Global assessment of shipping emissions in 2015 on a high spatial and temporal resolution. Atmos. Environ. 2017, 167, 403–415. [Google Scholar] [CrossRef]
  60. US EPA Office of Research and Development. CMAQv5.0.2 (5.0.2); Zenodo: Genève, Switzerland, 2014. [CrossRef]
  61. Fink, L.; Karl, M.; Matthias, V.; Oppo, S.; Kranenburg, R.; Kuenen, J.; Moldanova, J.; Jutterström, S.; Jalkanen, J.-P.; Majamäki, E. Potential impact of shipping on air pollution in the Mediterranean region—A multimodel evaluation: Comparison of photooxidants NO2 and O3. Atmos. Chem. Phys. 2023, 23, 1825–1862. [Google Scholar] [CrossRef]
  62. Nordic Drones. AEROMON BH-12, Real-Time Data from the Air. 2023. Available online: https://nordicdrones.fi/en/products/product-information/aeromon-bh-12/ (accessed on 20 September 2023).
  63. Kangasniemi, O.; Simonen, P.; Moldanová, J.; Timonen, H.; Barreira, L.M.F.; Hellén, H.; Jalkanen, J.-P.; Majamäki, E.; D’anna, B.; Lanzafame, G.; et al. Volatility of a Ship’s Emissions in the Baltic Sea Using Modelling and Measurements in Real-World Conditions. Atmosphere 2023, 14, 1175. [Google Scholar] [CrossRef]
  64. Sander, R.; Kerkweg, A.; Jöckel, P.; Lelieveld, J. Technical note: The new comprehensive atmospheric chemistry module MECCA. Atmos. Chem. Phys. 2005, 5, 445–450. [Google Scholar] [CrossRef]
  65. Sander, R.; Baumgaertner, A.; Gromov, S.; Harder, H.; Jöckel, P.; Kerkweg, A.; Kubistin, D.; Regelin, E.; Riede, H.; Sandu, A.; et al. The atmospheric chemistry box model CAABA/MECCA-3.0. Geosci. Model Dev. 2011, 4, 373–380. [Google Scholar] [CrossRef]
  66. Schack, C.J.; Pratsinis, S.E.; Friedlander, S. A general correlation for deposition of suspended particles from turbulent gases to completely rough surfaces. Atmos. Environ. (1967) 1985, 19, 953–960. [Google Scholar] [CrossRef]
  67. Chosson, F.; Paoli, R.; Cuenot, B. Ship plume dispersion rates in convective boundary layers for chemistry models. Atmos. Chem. Phys. 2008, 8, 4841–4853. [Google Scholar] [CrossRef]
  68. von Glasow, R.; Lawrence, M.G.; Sander, R.; Crutzen, P.J. Modeling the chemical effects of ship exhaust in the cloud-free marine boundary layer. Atmos. Chem. Phys. 2003, 3, 233–250. [Google Scholar] [CrossRef]
  69. Vehkamäki, H.; Kulmala, M.; Napari, I.; Lehtinen, K.E.J.; Timmreck, C.; Noppel, M.; Laaksonen, A. An improved parameterization for sulfuric acid–water nucleation rates for tropospheric and stratospheric conditions. J. Geophys. Res. Atmos. 2002, 107, AAC 3-1–AAC 3-10. [Google Scholar] [CrossRef]
  70. Vehkamäki, H.; Kulmala, M.; Lehtinen, K.E.J.; Noppel, M. Modelling binary homogeneous nucleation of water−sulfuric acid vapours: Parameterisation for high temperature emissions. Environ. Sci. Technol. 2003, 37, 3392–3398. [Google Scholar] [CrossRef] [PubMed]
  71. Kulmala, M.; Vehkamäki, H.; Petäjä, T.; Maso, M.D.; Lauri, A.; Kerminen, V.-M.; Birmili, W.; McMurry, P. Formation and growth rates of ultrafine atmospheric particles: A review of observations. J. Aerosol Sci. 2004, 35, 143–176. [Google Scholar] [CrossRef]
  72. Fridell, E.; Salberg, H.; Salo, K. Measurements of Emissions to Air from a Marine Engine Fueled by Methanol. J. Mar. Sci. Appl. 2021, 20, 138–143. [Google Scholar] [CrossRef]
  73. Kerminen, V.-M.; Pirjola, L.; Boy, M.; Eskola, A.; Teinilä, K.; Laakso, L.; Asmi, A.; Hienola, J.; Lauri, A.; Vainio, V.; et al. Interaction between SO2 and submicron atmospheric particles. Atmos. Res. 2000, 54, 41–57. [Google Scholar] [CrossRef]
  74. Donahue, N.M.; Robinson, A.L.; Stanier, C.O.; Pandis, S.N. Coupled partitioning, dilution, and chemical aging of semivolatile organics. Environ. Sci. Technol. 2006, 40, 2635–2643. [Google Scholar] [CrossRef]
  75. Byun, D.; Schere, K.L. Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System. Appl. Mech. Rev. 2006, 59, 51–77. [Google Scholar] [CrossRef]
  76. 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.; et al. Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1. Geosci. Model Dev. 2017, 10, 1703–1732. [Google Scholar] [CrossRef] [PubMed]
  77. Whitten, G.Z.; Heo, G.; Kimura, Y.; McDonald-Buller, E.; Allen, D.T.; Carter, W.P.; Yarwood, G. A new condensed toluene mechanism for Carbon Bond: CB05-TU☆. Atmos. Environ. 2010, 44, 5346–5355. [Google Scholar] [CrossRef]
  78. Nenes, A.; Pandis, S.N.; Pilinis, C. ISORROPIA: A New Thermodynamic Equilibrium Model for Multiphase Multicomponent Inorganic Aerosols. Aquat. Geochem. 1998, 4, 123–152. [Google Scholar] [CrossRef]
  79. Fountoukis, C.; Nenes, A. ISORROPIA II: A computationally efficient thermodynamic equilibrium model for K+–Ca2+–Mg2+–NH4+–Na+–SO42−–NO3–Cl–H2O aerosols. Atmos. Chem. Phys. 2007, 7, 4639–4659. [Google Scholar] [CrossRef]
  80. Carlton, A.G.; Bhave, P.V.; Napelenok, S.L.; Edney, E.O.; Sarwar, G.; Pinder, R.W.; Pouliot, G.A.; Houyoux, M. Model representation of secondary organic aerosol in CMAQv4.7. Environ. Sci. Technol. 2010, 44, 8553–8560. [Google Scholar] [CrossRef] [PubMed]
  81. Pye, H.O.T.; Pouliot, G.A. Modeling the role of alkanes, polycyclic aromatic hydrocarbons, and their oligomers in secondary organic aerosol formation. Environ. Sci. Technol. 2012, 46, 6041–6047. [Google Scholar] [CrossRef] [PubMed]
  82. 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.; et al. Incremental testing of the Community Multiscale Air Quality (CMAQ) modeling system version 4.7. Geosci. Model Dev. 2010, 3, 205–226. [Google Scholar] [CrossRef]
  83. Kelly, J.T.; Bhave, P.V.; Nolte, C.G.; Shankar, U.; Foley, K.M. Simulating emission and chemical evolution of coarse sea-salt particles in the Community Multiscale Air Quality (CMAQ) model. Geosci. Model Dev. 2010, 3, 257–273. [Google Scholar] [CrossRef]
  84. Guenther, A.B.; Jiang, X.; Heald, C.L.; Sakulyanontvittaya, T.; Duhl, T.; Emmons, L.K.; Wang, X. The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): An extended and updated framework for modeling biogenic emissions. Geosci. Model Dev. 2012, 5, 1471–1492. [Google Scholar] [CrossRef]
  85. Agrawal, H.; Welch, W.A.; Henningsen, S.; Miller, J.W.; Cocker, D.R. Emissions from main propulsion engine on container ship at sea. J. Geophys. Res. 2010, 115, D23205. [Google Scholar] [CrossRef]
  86. Agrawal, H.; Welch, W.A.; Miller, J.W.; Cocker, D.R. Emission measurements from a crude oil tanker at sea. Environ. Sci. Technol. 2008, 42, 7098–7103. [Google Scholar] [CrossRef] [PubMed]
  87. Sippula, O.; Stengel, B.; Sklorz, M.; Streibel, T.; Rabe, R.; Orasche, J.; Lintelmann, J.; Michalke, B.; Abbaszade, G.; Radischat, C.; et al. Particle emissions from a marine engine: Chemical composition and aromatic emission profiles under various operating conditions. Environ. Sci. Technol. 2014, 48, 11721–11729. [Google Scholar] [CrossRef]
  88. Reichle, L.J.; Cook, R.; Yanca, C.A.; Sonntag, D.B. Development of organic gas exhaust speciation profiles for nonroad spark-ignition and compression-ignition engines and equipment. J. Air Waste Manag. Assoc. 2015, 65, 1185–1193. [Google Scholar] [CrossRef] [PubMed]
  89. Schultze, M.; Rockel, B. Direct and semi-direct effects of aerosol climatologies on long-term climate simulations over Europe. Clim. Dyn. 2018, 50, 3331–3354. [Google Scholar] [CrossRef]
  90. Petrik, R.; Geyer, B.; Rockel, B. On the diurnal cycle and variability of winds in the lower planetary boundary layer: Evaluation of regional reanalyses and hindcasts. Tellus A Dyn. Meteorol. Oceanogr. 2022, 73, 1804294. [Google Scholar] [CrossRef]
  91. Inness, A.; Ades, M.; Agustí-Panareda, A.; Barré, J.; Benedictow, A.; Blechschmidt, A.-M.; Dominguez, J.J.; Engelen, R.; Eskes, H.; Flemming, J.; et al. The CAMS reanalysis of atmospheric composition. Atmos. Chem. Phys. 2019, 19, 3515–3556. [Google Scholar] [CrossRef]
  92. Maricq, M.M.; Xu, N.; Chase, R.E. Measuring Particulate Mass Emissions with the Electrical Low Pressure Impactor. Aerosol Sci. Technol. 2006, 40, 68–79. [Google Scholar] [CrossRef]
  93. Lähde, T.; Rönkkö, T.; Virtanen, A.; Schuck, T.J.; Pirjola, L.; Hämeri, K.; Kulmala, M.; Arnold, F.; Rothe, D.; Keskinen, J. Heavy Duty diesel engine exhaust aerosol particle and ion measurements. Environ. Sci. Technol. 2009, 43, 163–168. [Google Scholar] [CrossRef] [PubMed]
  94. Rönkkö, T.; Lähde, T.; Heikkilä, J.; Pirjola, L.; Bauschke, U.; Arnold, F.; Schlager, H.; Rothe, D.; Yli-Ojanperä, J.; Keskinen, J. Effects of gaseous sulphuric acid on diesel exhaust nanoparticle formation and characteristics. Environ. Sci. Technol. 2013, 47, 11882–11889. [Google Scholar] [CrossRef]
  95. Rönkkö, T.; Virtanen, A.; Kannosto, J.; Keskinen, J.; Lappi, M.; Pirjola, L. Nucleation mode particles with a nonvolatile core in the exhaust of a heavy duty diesel vehicle. Environ. Sci. Technol. 2007, 41, 6384–6389. [Google Scholar] [CrossRef]
  96. Arnold, F.; Pirjola, L.; Rönkkö, T.; Reichl, U.; Schlager, H.; Lähde, T.; Heikkilä, J.; Keskinen, J. First online measurements of sulfuric acid gas in modern heavy-duty diesel engine exhaust: Implications for nanoparticle formation. Environ. Sci. Technol. 2012, 46, 11227–11234. [Google Scholar] [CrossRef] [PubMed]
  97. Pirjola, L.; Karl, M.; Rönkkö, T.; Arnold, F. Model studies of volatile diesel exhaust particle formation: Are organic vapours involved in nucleation and growth? Atmos. Chem. Phys. 2015, 15, 10435–10452. [Google Scholar] [CrossRef]
  98. Celik, S.; Drewnick, F.; Fachinger, F.; Brooks, J.; Darbyshire, E.; Coe, H.; Paris, J.-D.; Eger, P.G.; Schuladen, J.; Tadic, I.; et al. Influence of vessel characteristics and atmospheric processes on the gas and particle phase of ship emission plumes: In situ measurements in the Mediterranean Sea and around the Arabian Peninsula. Atmos. Chem. Phys. 2020, 20, 4713–4734. [Google Scholar] [CrossRef]
  99. Murphy, B.N.; Woody, M.C.; Jimenez, J.L.; Carlton, A.M.G.; Hayes, P.L.; Liu, S.; Ng, N.L.; Russell, L.M.; Setyan, A.; Xu, L.; et al. Semivolatile POA and parameterized total combustion SOA in CMAQv5.2: Impacts on source strength and partitioning. Atmos. Chem. Phys. 2017, 17, 11107–11133. [Google Scholar] [CrossRef] [PubMed]
  100. Koo, B.; Knipping, E.; Yarwood, G. 1.5-Dimensional volatility basis set approach for modeling organic aerosol in CAMx and CMAQ. Atmos. Environ. 2014, 95, 158–164. [Google Scholar] [CrossRef]
  101. Lee, B.-K.; Choi, S.-D.; Shin, B.; Kim, S.-J.; Lee, S.-J.; Kim, D.-G.; Lee, G.; Kang, H.-J.; Kim, H.-S.; Park, D.-Y. Sensitivity analysis of volatile organic compounds to PM2.5 concentrations in a representative industrial city of Korea. Asian J. Atmos. Environ. 2023, 17, 3. [Google Scholar] [CrossRef]
  102. Benchrif, A.; Tahri, M.; Oujidi, B.; Bounouira, H.; Zahry, F. Comprehensive Analysis of PM10 and PM2.5 Composition in four Moroccan Cities: Emission Sources and Geographical Variations. Presented at the IGAC-iCACGP ECR Conference, Online, 17 November 2023. [Google Scholar]
  103. Birmili, W.; Weinhold, K.; Rasch, F.; Sonntag, A.; Sun, J.; Merkel, M.; Wiedensohler, A.; Bastian, S.; Schladitz, A.; Löschau, G.; et al. Long-term observations of tropospheric particle number size distributions and equivalent black carbon mass concentrations in the German Ultrafine Aerosol Network (GUAN). Earth Syst. Sci. Data 2016, 8, 355–382. [Google Scholar] [CrossRef]
Figure 1. Shipping lane (blue dots), measurement station = mobile lab (single blue square) and the wind direction (284° ± 15°, blue arrow). The distance from the ship lane to the mobile laboratory is 530 m.
Figure 1. Shipping lane (blue dots), measurement station = mobile lab (single blue square) and the wind direction (284° ± 15°, blue arrow). The distance from the ship lane to the mobile laboratory is 530 m.
Toxics 12 00432 g001
Figure 2. Illustration of the model structure, as presented in Karl et al. [52]. Dashed outline contains the MAFOR model. External modules are MECCA v4.0 and MOSAIC solver (the models are not part of MAFOR but included in the model distribution).
Figure 2. Illustration of the model structure, as presented in Karl et al. [52]. Dashed outline contains the MAFOR model. External modules are MECCA v4.0 and MOSAIC solver (the models are not part of MAFOR but included in the model distribution).
Toxics 12 00432 g002
Figure 3. Approach for receiving ship emissions with adjusted VOC EFs.
Figure 3. Approach for receiving ship emissions with adjusted VOC EFs.
Toxics 12 00432 g003
Figure 4. Dilution of NOx on logarithmic scale; dark red dots are the measured NOx values. Stack, measured by Aeromon BH-12 with an uncertainty of the measurement device of 19.7%, stack with 1:8 dilution, mobile lab, measured by airpointer 4D measurement system, background, measured by airpointer 4D measurement system with an uncertainty of 1 ppb (2 µg/m3).
Figure 4. Dilution of NOx on logarithmic scale; dark red dots are the measured NOx values. Stack, measured by Aeromon BH-12 with an uncertainty of the measurement device of 19.7%, stack with 1:8 dilution, mobile lab, measured by airpointer 4D measurement system, background, measured by airpointer 4D measurement system with an uncertainty of 1 ppb (2 µg/m3).
Toxics 12 00432 g004
Figure 5. Comparison of the modeled particle size distributions with the observed particle size distributions. “ferry” = measured particle size distribution at the stack, with the shaded blue area indicating the 25th and 75th percentile. “ferry mod” = simulated particle size distribution at the stack. “SNIFFER (plume)” = particle size distributions recorded at the sniffer at the plume arrival with the shaded red area indicating the uncertainty range of ±40%. “SNIFFER (no plume)” = particle size distributions recorded at the sniffer without a plume event, with shaded grey area indicating the 25th and 75th percentile. “backgr mod (530 m)” = simulated particle size distribution at the sniffer without a plume event. “plume mod (530 m)” = particle size distributions simulated at the sniffer at the plume arrival.
Figure 5. Comparison of the modeled particle size distributions with the observed particle size distributions. “ferry” = measured particle size distribution at the stack, with the shaded blue area indicating the 25th and 75th percentile. “ferry mod” = simulated particle size distribution at the stack. “SNIFFER (plume)” = particle size distributions recorded at the sniffer at the plume arrival with the shaded red area indicating the uncertainty range of ±40%. “SNIFFER (no plume)” = particle size distributions recorded at the sniffer without a plume event, with shaded grey area indicating the 25th and 75th percentile. “backgr mod (530 m)” = simulated particle size distribution at the sniffer without a plume event. “plume mod (530 m)” = particle size distributions simulated at the sniffer at the plume arrival.
Toxics 12 00432 g005
Figure 6. Changes in PM2.5 for the CTM run with adjusted ship emissions to initial ship emissions, displayed as mean values for (a) = March 2015; (b) = June 2015; (c) = September 2015; (d) = December 2015.
Figure 6. Changes in PM2.5 for the CTM run with adjusted ship emissions to initial ship emissions, displayed as mean values for (a) = March 2015; (b) = June 2015; (c) = September 2015; (d) = December 2015.
Toxics 12 00432 g006
Table 1. Mass distribution for background and ship aerosols; mass concentration is displayed in ng/m3; MSAp = methane sulfonate; PBA = primary biological material; v = volatile; nv = non-volatile; modes: NU = nucleation (Dp < 10 nm); AI = Aitken (10 nm < Dp < 50 nm); AS1 = accumulation 1 (50 nm < Dp < 500 nm); AS2 = accumulation 2 (>500 nm).
Table 1. Mass distribution for background and ship aerosols; mass concentration is displayed in ng/m3; MSAp = methane sulfonate; PBA = primary biological material; v = volatile; nv = non-volatile; modes: NU = nucleation (Dp < 10 nm); AI = Aitken (10 nm < Dp < 50 nm); AS1 = accumulation 1 (50 nm < Dp < 500 nm); AS2 = accumulation 2 (>500 nm).
MODEDp (m)SIGH2SO4OCNH4NO3MSApSALTPBAECASHMTOT
vvvvvnvnvnvnv
BACKGROUND AEROSOL
NU9.90 × 10−91.700.0010.00200000000.003
AI2.3 × 10−81.520.320.250.150.1500.0700.270.171.38
AS11.0 × 10−71.7554330104323
AS27.0 × 10−72.1014117703012761
SHIP AEROSOL
NU9.6 × 10−91.6500.02600000.011000.05
AI1.5 × 10−81.35283989.9573800424.2494.301886.14
AS11.3 × 10−71.7512063657241161001567.475636438038.47
AS27.5 × 10−72.1010,72632,535000013,943.78357510,72671,505.78
Table 2. Size ranges of the size classes. S1: 1 nm–10 nm; S2: 10 nm–20 nm; S3: 20 nm–50 nm; S4: 50 nm–100 nm; S5: 100 nm–300 nm; S6: 300 nm–600 nm; S7: >600 nm; differences are referring to the relative deviations in the modeled and measured data. CF = correcting factor., Because the CF is 1, 1.0 means no changes. The CF with the smallest deviation is shown in bold.
Table 2. Size ranges of the size classes. S1: 1 nm–10 nm; S2: 10 nm–20 nm; S3: 20 nm–50 nm; S4: 50 nm–100 nm; S5: 100 nm–300 nm; S6: 300 nm–600 nm; S7: >600 nm; differences are referring to the relative deviations in the modeled and measured data. CF = correcting factor., Because the CF is 1, 1.0 means no changes. The CF with the smallest deviation is shown in bold.
Difference in Number Concentrations Compared to Measurements (in %)
CFS1S2S3S4S5S6S7|S2| + |S3| +
|S4| + |S5|
0.5038.6637.72−206.513.63−7.7977.7082.83−172.94
1.00−11.2737.30−117.903.91−6.4877.8182.95−83.17
1.50−55.6133.30−69.654.16−5.2277.9383.07−37.41
2.00−87.1427.74−40.204.41−3.9778.0483.18−12.02
2.50−110.0521.19−20.724.68−2.7178.1583.292.44
3.00−128.9813.76−6.945.05−1.4678.2683.4010.41
3.10−132.1412.34−4.905.13−1.2378.2883.4211.33
3.20−136.1610.46−2.425.24−0.9478.3183.4512.33
3.30−139.288.95−0.575.33−0.7178.3283.4713.00
3.40−142.357.391.195.42−0.4878.3583.4913.52
3.50−146.315.373.315.56−0.1878.3783.5114.05
3.60−149.343.714.915.670.0578.3983.5314.35
3.70−152.302.086.395.800.2678.4183.5514.53
3.80−155.240.427.805.930.4978.4383.5714.63
3.90−158.95−1.719.486.100.7578.4683.6014.63
4.00−162.00−3.5310.826.260.9878.4883.6214.52
4.50−177.29−13.1316.627.162.0278.5883.7212.67
5.00−193.18−23.9721.428.373.1078.6983.838.93
Difference in Mean Diameter Compared to Measurements (in %)
CFS1S2S3S4S5S6S7|S2| + |S3| +
|S4| + |S5|
0.508.556.6610.53−12.237.91−3.326.23162.50
1.0010.353.1012.30−12.287.88−3.346.21137.99
1.5011.580.2912.63−12.337.84−3.356.19117.65
2.0011.49−1.9112.40−12.377.80−3.376.17100.36
2.5010.82−3.5911.88−12.397.75−3.386.1585.21
3.0010.04−4.8611.20−12.367.70−3.406.1372.03
3.109.90−5.0511.07−12.357.69−3.406.1360.29
3.209.73−5.2910.89−12.337.68−3.416.1249.45
3.309.60−5.4610.75−12.317.67−3.416.1239.75
3.409.47−5.6210.60−12.307.66−3.416.1230.62
3.509.31−5.8110.40−12.277.64−3.426.1121.99
3.609.19−5.9610.25−12.257.63−3.426.1113.93
3.709.07−6.0910.09−12.227.62−3.426.1013.93
3.808.96−6.219.93−12.197.60−3.436.1014.39
3.908.82−6.359.73−12.167.58−3.436.1014.39
4.008.70−6.469.56−12.127.57−3.436.0914.39
4.508.16−6.888.67−11.927.48−3.456.0720.30
5.007.63−7.127.71−11.647.38−3.476.0527.89
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fink, L.; Karl, M.; Matthias, V.; Weigelt, A.; Irjala, M.; Simonen, P. Using the Multicomponent Aerosol FORmation Model (MAFOR) to Determine Improved VOC Emission Factors in Ship Plumes. Toxics 2024, 12, 432. https://doi.org/10.3390/toxics12060432

AMA Style

Fink L, Karl M, Matthias V, Weigelt A, Irjala M, Simonen P. Using the Multicomponent Aerosol FORmation Model (MAFOR) to Determine Improved VOC Emission Factors in Ship Plumes. Toxics. 2024; 12(6):432. https://doi.org/10.3390/toxics12060432

Chicago/Turabian Style

Fink, Lea, Matthias Karl, Volker Matthias, Andreas Weigelt, Matti Irjala, and Pauli Simonen. 2024. "Using the Multicomponent Aerosol FORmation Model (MAFOR) to Determine Improved VOC Emission Factors in Ship Plumes" Toxics 12, no. 6: 432. https://doi.org/10.3390/toxics12060432

APA Style

Fink, L., Karl, M., Matthias, V., Weigelt, A., Irjala, M., & Simonen, P. (2024). Using the Multicomponent Aerosol FORmation Model (MAFOR) to Determine Improved VOC Emission Factors in Ship Plumes. Toxics, 12(6), 432. https://doi.org/10.3390/toxics12060432

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop