Simulating the Impacts of Marine Organic Emissions on Global Atmospheric Chemistry and Aerosols Using an Online-Coupled Meteorology and Chemistry Model

To realistically simulate the impacts of marine isoprene and primary organic aerosols (POA) on atmospheric chemistry, a unified model framework with online emissions, comprehensive treatment of gas-phase chemistry, and advanced aerosol microphysics is required. In this work, the global-through-urban WRF/Chem model (GU-WRF/Chem) implemented with the online emissions of marine isoprene and size-resolved marine POA is applied to examine such impacts. The net effect of these emissions was increased surface concentrations of isoprene and organic aerosols and decreased surfaced concentrations of hydroxyl radical and ozone over most marine regions. With the inclusion of these emissions, GU-WRF/Chem better predicted the surface concentrations of isoprene and organic aerosols and the aerosol number size distribution when compared to measurements in clean marine conditions.


Introduction
The chemical composition and aerosols of the natural background atmosphere are important to model predictions of air quality and climate [1] [2].Most models, however, have an incomplete emission inventory of natural marine aerosols and precursor gases beyond sea-salt aerosol and dimethyl sulfide.Two major emission types neglected from the default setting of many current models are biogenic volatile organic compounds (BVOCs) of marine origin and primary organic aerosol (POA) associated with sea spray.Isoprene is the most widely-reported BVOC species with a marine emission source [3] and has been shown to have a small impact on surface ozone (O 3 ) and secondary organic aerosol (SOA) formation in coastal regions [4].Marine POA has been suggested to affect the surface aerosol mass concentration and cloud condensation nuclei (CCN) in remote marine regions [5]- [8].However, large uncertainties in chemical composition [9] [10] and magnitude [11] [12] of sea spray aerosol have been shown to play an important role in determining the direct and indirect radiative forcing of marine aerosols [13].Previous modeling efforts have included marine isoprene and POA emissions in several atmospheric chemistry [14]- [19] and online-coupled global climate models [6] [13] [20] [21].
In this study, we use the online-coupled global-through-urban weather and forecasting model with chemistry (GU-WRF/Chem) to quantify the impact of online emissions of marine isoprene and POA on atmospheric chemistry and aerosol abundance in remote marine and coastal regions.GU-WRF/Chem is capable of simulating these impacts across a wide range of spatial scales while reducing uncertainties from the use of offline-coupled model systems with inconsistent model treatments [22].Comprehensive treatments of gas-phase chemistry and aerosol microphysics, in conjunction with online emissions of marine BVOC and POA, are expected to improve model predictions of surface concentrations of BVOCs, O 3 , and organic aerosols.The overall objective of this study is to estimate the impact of marine isoprene and POA emissions on remote marine atmospheric chemistry through comprehensive evaluation of model predictions against observations from marine regions.

Model Description and Setup
GU-WRF/Chem is developed upon the global WRF (GWRF) version 3.0 [23]- [26] and the mesoscale WRF/ Chem version 3.0 [27], and described in detail by Zhang et al. [22].Simulations in this study are conducted with the 2005 Carbon Bond chemical mechanism with global extension (CB05_GE) [28].The Model of Aerosol Dynamics, Reaction, Ionization, and Dissolution (MADRID) aerosol microphysics module [29] [30] used in GU-WRF/Chem treats all major aerosol chemical and microphysical processes including inorganic aerosol thermodynamic equilibrium, SOA formation, nucleation, gas/particle mass transfer, condensation, and coagulation with eight size sections to represent the particle size distribution.MADRID was selected from among the three aerosol modules in this version of GU-WRF/Chem because it is the only one that treats isoprene as an SOA precursor.32 aerosol types were considered in GU-WRF/Chem with anthropogenic and terrestrial emissions described in detail in Zhang et al. [22].Model configurations not pertaining to chemistry were consistent with the baseline GWRF simulation described in Zhang et al. [26].The model was reinitialized from the NCEP Final Global Data Assimilation System every 24 hours to allow meteorological feedbacks while constraining the simulation with re-initialized meteorological conditions based on observations.Marine isoprene emissions implemented into GU-WRF/Chem were based on the parameterization developed by Gantt et al. [31], which experimentally determined that marine isoprene production is a function of phytoplankton abundance, speciation, and photosynthetically active radiation.Model implementation of marine isoprene emissions of Gantt et al. [31] has been previously described in detail in Gantt et al. [4], which uses the isoprene production rate for diatoms, model-predicted solar radiation, and Level 3 monthly average chlorophylla concentration ([Chl-a]) and diffuse attenuation coefficient derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) satellite to drive the emissions.Despite evidence for marine production of monoterpenes [32], this emission source was not included in GU-WRF/Chem due to the very low production rates relative to marine isoprene emissions [33].Marine isoprene emissions were calculated online within GU-WRF/Chem and added to the existing terrestrial isoprene emissions of Guenther et al. [34].
Marine POA emissions implemented into GU-WRF/Chem were based on the parameterization developed by Gantt et al. [35].These size-dependent aerosol emissions were based on the calculation of the organic mass fraction of sea spray aerosol (OM SSA ), which Gantt et al. [35] found to be positively related to phytoplankton abundance and aerosol diameter and negatively related to surface wind speed.Within GU-WRF/Chem, OM SSA- was calculated online using monthly-mean SeaWiFS-derived [Chl-a] and model-predicted surface wind speed and applied to the existing sea spray emission module of Gong et al. [36] and O'Dowd et al. [37].Model implementation of these emissions performed here was similar to the method described in detail by Meskhidze et al. [6].The size-dependent mass and number emissions of marine POA were added online to the existing terrestrial POA based on the sectional size bins of MADRID in GU-WRF/Chem, with the atmospheric evolution of marine/terrestrial POA and SOA is simulated explicitly.
Two sets of GU-WRF/Chem simulations were performed: 1) the Default simulation without marine organic emissions and 2) the MarOrg simulation with marine isoprene and POA emissions described above.Each set of simulations was conducted over a 12-month period for two years (2001 and 2009) with a one-month spin-up for each year at a spatial resolution of 4˚ × 5˚ and 27 vertical layers (1 st vertical layer height: 35 m) up to 50 mb.These years were chosen for direct evaluation with field campaigns [38]- [40] measuring marine isoprene emissions and surface concentrations of isoprene and organic aerosols in remote marine environments.Other than the marine isoprene and POA emissions, all other chemical and meteorological model configurations for the Default and MarOrg simulations were kept constant.

Trace Gases
Figure 1 shows that marine isoprene emissions in the MarOrg simulation resulted in a change in surface isoprene concentrations from the simulation with default treatments over marine regions.The strong summertime solar radiation and high [Chl-a] levels in northern Atlantic and Pacific Oceans contributed to the high emission rates in those regions, while moderate ocean biological activity throughout the year over much of the Southern Ocean also resulted in high marine isoprene emission rates.Over remote marine regions, marine isoprene emissions increased surface isoprene concentrations by 5 -10 ppt (see Figure 1).This increase was within the range of most isoprene concentrations (<1 to 300 ppt, typically < 50 ppt) observed over remote marine regions [3].The large emission rates occurring near coastal regions (and some inland lakes) were due to the combination of isolated regions of very high [Chl-a] and coarse (4˚ × 5˚) grid resolution.The greatest increases in surface isoprene concentrations were predicted over the northern Atlantic/Pacific and Southern Oceans where marine isoprene emission rates are high.These areas also had low oxidant levels and therefore long lifetimes of isoprene emitted from marine sources and transported from terrestrial regions.Over terrestrial areas such as the rainforests of South America, Southeast Asia, and Africa, a mixture of positive and negative changes in the surface isoprene concentrations is evident.Such changes were likely due to the fact that terrestrial emissions were computed online and were sensitive to the slight variability in meteorology occurring due to model reinitializations.Although the changes in surface isoprene concentrations over terrestrial areas were of similar magnitude as those of marine regions, they were much smaller when calculated in terms of percentage changes (typically < 5%).
Predictions of other atmospheric chemistry-relevant gases in remote marine regions were also affected by emissions of marine isoprene.The surface concentrations of hydroxyl radical (OH) decreased by 10 −3 ppt or more in parts of the northern Pacific and Southern Oceans(see Figure 1) because of its reaction with the additional isoprene from marine sources.Small changes in surface O 3 concentration due to marine isoprene emissions(both increases and decreases of over 0.1 ppb) were also predicted; oceanic regions immediately downwind of anthropogenic pollution centers (China, U.S., etc.) experienced increased concentrations while remote marine regions typically have decreased concentrations (Figure 1).Decreases in predicted surface O 3 concentrations over most of the eastern North Pacific Ocean were consistent with Gantt et al. [16], who found that the emission of isoprene and monoterpenes from marine sources caused a small (−0.1%) decrease in CMAQ-predicted O 3 mixing ratios off the Pacific coast of the U.S.
Although coincident measurements of isoprene emissions and concentrations over remote marine regions are very limited [3], the May 2001 simulation period allowed for model evaluation of predicted isoprene emissions/ concentrations in the northwest Pacific Ocean [38].Compared to the Matsunaga et al. [38] observations, the model overpredicted marine isoprene emissions (especially during the daytime) and underpredicted surface concentrations (Figure 2).Concurrent over predictions of emissions and under predictions of surface concentrations suggested that the simulated isoprene lifetime by GU-WRF/Chem over the marine boundary layer was too short, potentially due to the coarse grid resolution.The low temporal resolution of the observed emissions in Matsunaga et al. [38] made a more quantitative analysis of the diurnal cycle of emissions difficult, but the few observations available were inconsistent with the large increase in daytime emissions predicted by the model.The very low surface isoprene concentrations predicted by the Default simulation do not support the suggestion from Matsunaga et al. [38] that terrestrial isoprene sources influence the surface concentrations in this region; although we cannot definitely rule out the possibility of terrestrial isoprene sources affecting surface concentrations because of the use of a coarse model grid resolution in our simulations, proximity to the Japanese coast (~500 km) and the short isoprene lifetime predicted by GU-WRF/Chem in the marine boundary layer.The chemical mechanism used in GU-WRF/Chem, CB05-GE, is part of series developed for urban photochemistry [41] that may have difficulty representing remote marine conditions despite recent updates [38].Nonetheless, according to Figure 2, the inclusion of marine isoprene emissions in the MarOrg simulation brought the predicted concentrations closer to observations.

Aerosols
Globally, the implementation of marine isoprene and POA emissions resulted in increased surface PM 2.5 organic aerosol concentrations (> 0.5 µg•m −3 ) over biologically-active remote marine regions in the northern Atlantic/Pacific and Southern Oceans (Figure 3).While the predicted increases in organic aerosol concentration occurred year-round in the Southern Ocean due to the persistent strong winds and moderate [Chl-a], increases in northern Atlantic and Pacific were the highest in the spring and lowest in the winter months.Surface organic aerosol concentrations in tropical oceanic regions were largely unchanged between the Default and Mar Org simulations due to low biological activity and sea spray aerosol emissions relative to mid-and high-latitude oceanic regions.Comparison of the predicted concentrations with observations in clean marine conditions compiled from the years 1973 to 2009 [9] indicated a latitudinal model bias, with over predictions in mid-and highlatitude regions and underpredictions in the tropics.Although only a few sites have measurements throughout the year (e.g., Mace Head, Amsterdam Island), comparison with model predictions in Figure 3 suggests that the MarOrg simulation did not replicate the observed seasonal cycle of concentrations.This result is consistent with Gantt et al. [42] who found that none of the existing sea spray aerosol-based marine POA emission parameterizations are able to replicate the observed seasonality of organic aerosol concentrations.
Comparison of GU-WRF/Chem and aerosol mass spectrometer observations from Mace Head, Ireland [39] for the year 2009 revealed that the inclusion of marine isoprene and POA emissions in the MarOrg simulation improved the prediction of surface submicron organic aerosol mass concentrations (see Figure 4).Without marine isoprene and POA emissions, the Default simulation frequently predicted surface organic aerosol concentrations  at Mace Head < 0.025 µg•m −3 during clean marine conditions (defined as periods having aerosol number concentration < 700 cm −3 , onshore wind direction, and black carbon concentrations < 0.05 µg•m −3 ) while the observed concentrations typically ranged between 0.05 and1 µg•m −3 [39].Organic aerosol concentrations rarely fell below 0.05 µg•m −3 in the MarOrg simulation.Model-measurement discrepancies persisted in the MarOrg simulation, especially during the marine organic aerosol plume period (16 August 2009) observed by Ovadnevaite et al. [40].During this plume event, the MarOrg simulation predicted organic aerosol concentrations (0.2 µg•m −3 ) that were an order of magnitude lower than the observed maximum concentrations (3.8 µg•m −3 ).The change of annual average submicron organic (secondary and primary) aerosol concentrations at Mace Head between the Default and MarOrg simulations was 0.15 µg•m −3 , a value similar to the observed annual concentrations (0.25 µg•m −3 ) at Mace Head during clean marine conditions.
During the marine organic aerosol plume (16 August 2009) observed by Ovadnevaite et al. [40] (see Figure 5), the MarOrg simulation had a slightly higher aerosol number concentration than the Default simulation in the size range over which both simulations had the largest underpredictions (0.05 to 0.5 µm in diameter) relative to observations.Much of this difference is attributable to the increased organic aerosol number concentrations for diameters > 0.05 µm.In other clean marine periods where the MarOrg and Default simulations predicted substantially different submicron organic aerosol mass concentrations, the increase in the aerosol number concentration for the MarOrg simulation in the 0.05 to 0.5 µm diametersize range is more evident.These model results are consistent with the observations of multiple marine organic aerosol plumes from Ovadnevaite et al. [40], who found that periods with the highest mass fraction of organics in submicron aerosols also had aerosol number size distributions peaking between 0.05 and 0.2 µm in diameter.It's clear, however, that the MarOrg simulation did not accurately simulate this marine organic aerosol plume both in terms of magnitude and size distribution.

Conclusion
In this study, we implemented marine isoprene and size-resolved marine POA emissions into GU-WRF/Chem, an online-coupled meteorology and chemistry model.The resulting increases in the surface concentrations of isoprene and organic aerosols over remote marine regions led to small decreases in the surface concentration of OH and O 3 and brought the model predictions into better agreement with isoprene emissions and surface concentrations in the northwest Pacific Ocean.Predictions of the aerosol number size distribution during an observed marine organic aerosol plume were also improved when these emissions were included.Although the inclusion of these emission sources did not bring the model into complete agreement with observed concentrations of isoprene and organic aerosols in remote marine regions, the improved predictions suggested that these emissions are important processes affecting atmospheric chemistry in these areas.As many coastal cities experience periods when clean marine conditions impact the background concentrations of pollutants, the inclusion of these marine emissions would benefit future air quality studies.The improved representation of aerosol number size distributions in clean marine conditions may also benefit the prediction of cloud condensation nuclei concentrations and associated cloud properties for future climate studies.For reference, the organic aerosol size distribution for both simulations is given by the dotted lines.Note that the sectional size distribution predicted by MADRID in GU-WRF/Chem has been mapped to a modal size distribution through smoothing to better compare with observations.

Figure 1 .
Figure 1.Multi-year (2001 and 2009) annual average marine isoprene emissions and differences in the predicted surface concentrations of isoprene, OH, and O 3 between the Default and MarOrg simulations.

Figure 2 .
Figure 2. Time series comparison of observed [38] and modeled isoprene emissions (top) and surface concentrations (bottom) from the northwest Pacific Ocean in May 2001.

Figure 3 .
Figure 3. Multi-year (2001 and 2009) seasonally-averaged changes in the surface PM 2.5 organic aerosol concentrations between the Default and MarOrg simulations superimposed with observed surface organic aerosol concentrations in clean marine conditions compiled by Gantt and Meskhidze [9].

Figure 4 .
Figure 4. Time series comparison of modeled organic aerosol mass concentrations by GU-WRF/Chem with those observed at Mace Head, Ireland for the year 2009.

Figure 5 .
Figure 5. Observed[39] and predicted aerosol size distribution for the Default and MarOrg simulations on 16 August 2009.For reference, the organic aerosol size distribution for both simulations is given by the dotted lines.Note that the sectional size distribution predicted by MADRID in GU-WRF/Chem has been mapped to a modal size distribution through smoothing to better compare with observations.