Subscriber access provided by UNIV OF TASMANIA
Article
Simulating Aqueous-Phase Isoprene-Epoxydiol (IEPOX) Secondary Organic Aerosol Production During the 2013 Southern Oxidant and Aerosol Study (SOAS) Sri Hapsari Budisulistiorini, Athanasios Nenes, Ann Marie Grover Carlton, Jason Douglas Surratt, V. Faye McNeill, and Havala Olson Taylor Pye Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05750 • Publication Date (Web): 10 Apr 2017 Downloaded from http://pubs.acs.org on April 10, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 31
Environmental Science & Technology
1
Simulating Aqueous-Phase Isoprene-Epoxydiol (IEPOX) Secondary Organic
2
Aerosol Production During the 2013 Southern Oxidant and Aerosol Study
3
(SOAS)
4 5
Sri Hapsari Budisulistiorini†,∞,#, Athanasios Nenes‡,§,¶,ƒ, Annmarie G. Carlton∫, Jason D. Surratt∞, V. Faye McNeill◊,*, Havala O. T. Pye†,*
6 7
†
National Exposure Research Laboratory, US EPA, Research Triangle Park, North Carolina 27711, United States
8 9
∞
Department of Environmental Sciences and Engineering, Gillings School of Global Public Health, University of North Carolina, Chapel Hill, North Carolina 27599, United States
10 11
#
12 13
‡
14 15
§
16 17
¶
18 19
ƒ
Institute of Environmental Research and Sustainable Development, National Observatory of Athens, Palea Penteli, GR-15236, Greece
20
∫
21 22
◊
Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
23
*
Earth Observatory of Singapore, Nanyang Technological University, Singapore 639798, Singapore School of Earth and Atmospheric Sciences, Georgia Institute of Technology, Atlanta, Georgia 30322, United States School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30322, United States Institute of Chemical Engineering Sciences, Foundation for Research and Technology - Hellas, GR-26504 Patras, Greece
Department of Chemistry, University of California, Irvine, California, 92617, United States
To whom correspondence should be addressed
24 25
For submission to Environmental Science and Technology
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 31
26
Abstract
27
The lack of statistically robust relationships between IEPOX (isoprene epoxydiol)-derived SOA
28
(IEPOX SOA) and aerosol liquid water and pH observed during the 2013 Southern Oxidant and
29
Aerosol Study (SOAS) emphasizes the importance of modeling the whole system to understand
30
the controlling factors governing IEPOX SOA formation. We present a mechanistic modeling
31
investigation predicting IEPOX SOA based on Community Multiscale Air Quality (CMAQ)
32
model algorithms and a recently introduced photochemical box model, simpleGAMMA. We aim
33
to 1) simulate IEPOX SOA tracers from the SOAS Look Rock ground site, 2) compare the two
34
model formulations, 3) determine the limiting factors in IEPOX SOA formation, and 4) test the
35
impact of a hypothetical sulfate reduction scenario on IEPOX SOA. The estimated IEPOX SOA
36
mass variability is in similar agreement (r2 ~ 0.6) with measurements. Correlations of the
37
estimated and measured IEPOX SOA tracers with observed aerosol surface area (r2 ~ 0.5–0.7),
38
rate of particle-phase reaction (r2 ~ 0.4–0.7), and sulfate (r2 ~ 0.4–0.5) suggest an important role
39
of sulfate in tracer formation via both physical and chemical mechanisms. A hypothetical 25%
40
reduction of sulfate results in ~70% reduction of IEPOX SOA formation, reaffirming the
41
importance of aqueous phase chemistry in IEPOX SOA production.
42
ACS Paragon Plus Environment
2
Page 3 of 31
43
Environmental Science & Technology
Introduction
44
Particles with aerodynamic diameters less than 2.5 µm (PM2.5) are known to affect human
45
health and regional climate forcing. 1 Organic aerosol (OA) constitutes a large fraction of PM2.5
46
mass.
47
plays an important role in the spatial and temporal distribution of aerosol optical thickness,
48
potentially causing regional cooling over the southeastern United States. 3-5
2
In summer, biogenic secondary organic aerosol (SOA) is a major fraction of OA and
49
Isoprene (2-methyl-1,3-butadiene) emission is estimated to be 600 Tg yr-1, 6 which makes it the
50
most abundant non-methane volatile organic compound (VOC) emitted globally. It also plays a
51
substantial role in SOA production. 7 On average, isoprene-derived SOA tracers account for up to
52
20% of total OA observed in the southeastern U.S. (SE US), and compounds derived from
53
isoprene epoxydiol (IEPOX) chemical pathways can account for >95% of isoprene SOA.
54
IEPOX-OA, a positive matrix factor linked to IEPOX-derived SOA (IEPOX SOA) tracers,
55
comprises one-third of ambient organic aerosol mass in urban and rural areas in the SE US, 9-11
56
suggesting a regional source. Among isoprene-derived SOA tracers, 2-methyltetrols (tetrols) and
57
their sulfate ester derivatives (organosulfates) contribute ~ 33% and ~ 34% of IEPOX-OA,
58
respectively, at the Look Rock (LRK), Tennessee site. 9
8,9
59
Inclusion of biogenic and anthropogenic SOA precursors and mechanisms in the Community
60
Multiscale Air Quality (CMAQ) model improves the temporal correlation between predicted and
61
semi-empirical secondary organic carbon (OC). 12 Addition of the acid-catalyzed reactive uptake
62
mechanism of IEPOX into CMAQ further improves predictions of tetrols against field
63
measurements.
64
reported a strong correlation between predicted and observed IEPOX SOA, although the
65
simulations underestimate the measurements by an order of magnitude.
13
Using the updated SOA module in the CMAQ model, Karambelas et al.
14
ACS Paragon Plus Environment
3
Environmental Science & Technology
66
McNeill et al.
15
Page 4 of 31
introduced GAMMA (Gas-Aerosol Model for Mechanism Analysis), a
67
detailed photochemical box model of aqueous-aerosol SOA (aaSOA) formation. A version of
68
GAMMA with simplified aqueous reaction mechanism, simpleGAMMA, has been shown to be
69
computationally efficient and suitable for coupling with large-scale atmospheric chemistry
70
models. 16 The aqueous aerosol phase mechanism in simpleGAMMA treats SOA formation from
71
glyoxal and IEPOX.
72
IEPOX-derived organosulfates (IEPOXOS) and tetrols, which have been quantified in 17,18
73
laboratory experiments,
have been shown to contribute up to half of all known IEPOX-
74
derived aerosol-phase species in the southeastern U.S.
75
IEPOX SOA is through acid-catalyzed ring-opening reactions of isoprene epoxides in concert
76
with nucleophilic addition in aqueous particles. 18-20 An airborne study over the southeastern US
77
found a strong dependency between IEPOXOS and acidic sulfate aerosol. 21 However, ground-
78
based studies in the eastern US suggest weak or no correlation between IEPOX SOA and aerosol
79
pH or liquid water content (LWC), but a positive correlation with sulfate.
80
preliminary modeling study employing simpleGAMMA found good correlation (r2 ~ 0.6)
81
between estimated IEPOX SOA and measurements made during the 2013 Southern Oxidant
82
Aerosol Study (SOAS) campaign. 9
8,9
The accepted formation pathway for
9,11
Additionally, a
83
This study aims to investigate the formation of two IEPOX SOA tracers, IEPOXOS and
84
tetrols, using CMAQ and simpleGAMMA algorithms, and evaluate the model response to
85
perturbations that replicate atmospheric conditions. Predictions from these two different
86
algorithms are compared with measurement results from the 2013 SOAS campaign. Only two
87
IEPOX SOA tracers are predicted and investigated by these algorithms since they are the most
ACS Paragon Plus Environment
4
Page 5 of 31
Environmental Science & Technology
88
abundant tracers during the entire field observation. 9 Limiting factors in IEPOX SOA formation
89
are investigated and the IEPOX SOA response due to changes in SOx emission were estimated.
90 91
Method
92
Two different representations of IEPOX SOA formation were implemented in box-model
93
form, namely CMAQ-box and simpleGAMMA, using conditions representative of the SOAS
94
Look Rock (LRK) site.
95
CMAQ-box. CMAQ-box is based on the algorithms in the CMAQ model, which are detailed 13
96
in the work of Pye et al.
97
uptake of gaseous IEPOX into the aerosol phase as a first order reaction:
98 99 100 101 102
Briefly, CMAQ treats IEPOX SOA formation as heterogeneous
→
(R1)
Governed by a first order heterogeneous rate constant, ,22 such that:
= −
(1)
where
=
!
(2)
" #$
103
SA is the aerosol surface area (µm2 cm-3) upon which uptake occurs, % is the mean molecular
104
speed (m s-1) in the gas phase, and &' is the effective particle radius (cm). ( is the diffusivity of
105
IEPOX in the gas phase, which is estimated as ( = 1.9 ,-. / 01 23
106
molecular weight (g mol-1).
107
uptake coefficient (6) calculated with the following. 23,24
108
0.5 4,
1 7
1
= 8 + ;< ∗ >?
:
15
with MW as
The heterogeneous reaction is parameterized using a reactive
1
@AB CBDEFGH IJKL M015M
(3)
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 31
109
Here, N is the accommodation coefficient (0.02) 15, O ∗ is the Henry’s Law coefficient (M atm-1),
110
R is the gas constant, T is temperature (K), (P is diffusivity in the aerosol phase (1 × 10-9 m2 s-1),
111
25
112
rate constant (s-1) for reaction of IEPOX in the aerosol phase. This formulation allows for the
113
model to predict whether the particle-phase reaction will occur throughout the particle, or as a
114
surface process based on the relative timescales for reaction Q1⁄ 'PRSTU W vs. aerosol-phase
115
diffusion Q&' . ⁄(P W. The Henry’s law coefficient was set to 3 × 107 M atm-1 based on the value
116
measured by Nguyen et al. and used in CMAQ starting with v5.226 and simpleGAMMA v1.0. 16
117
q is the diffuso-reactive parameter Q&' @ 'PRSTU ⁄(P W, and 'PRSTU is the pseudo-first order
simpleGAMMA. A detailed description of simpleGAMMA is available in Woo and McNeill.
118
16
119
aqueous-phase reaction23 to form IEPOX SOA consisting of tetrols and IEPOXOS: → PM →
120 121
Briefly, simpleGAMMA represents IEPOX SOA in terms of aqueous uptake followed by
(R2)
Formation of the aaSOA of IEPOX SOA (mol L-1 of aerosol) is given by:
122
=
CXD,Z[\]^ >?
−
CXD,Z[\]^ < ∗ >?
_PM ` + ∑C &C,PM
(4)
123
where, is the gas-phase partial pressure of species IEPOX. &C,PM are the rates of formation
124
or
125
(&C,PM = 'PRSTU _PM `). b, is the gas-to-aerosol mass transfer coefficient of
126
IEPOX, which can be described as follows: 16
127
decomposition
b, =
c d
of
IEPOX
due
to
chemical
reactions
1 !
in
aerosol
phase
(5)
" d#e
128
Particle-phase Reaction Parameters. The particle-phase reaction rate constant for IEPOX is
129
calculated assuming protonation of the epoxide oxygen and nucleophilic addition in both CMAQ
130
and simpleGAMMA following an A2 mechanism by Eddingsaas et al. 19 While CMAQ does not
ACS Paragon Plus Environment
6
Page 7 of 31
Environmental Science & Technology
131
treat NH4+ as an acid for epoxide uptake purposes, we implemented 'PRSTU such that it was the
132
same in both simpleGAMMA and CMAQ-box calculations using the most recent values from
133
literature which included updating the organosulfate formation rate constant following Riedel et
134
al. 27 Thus, the particle-phase rate constant for both CMAQ and simpleGAMMA simulations in
135
this study is:
136 137
'PRSTU = < f ,gPR h< f i3j +