Subscriber access provided by University of Sussex Library
Environmental Modeling
Parameterized yields of semi-volatile products from isoprene oxidation under different NOx levels: impacts of chemical aging and wall-loss of reactive gases Li Xing, ManishKumar Shrivastava, T. -M. Fu, Pontus Roldin, Yun Qian, Lu Xu, Nga Lee Ng, John E. Shilling, Alla Zelenyuk, and Christopher David Cappa Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b00373 • Publication Date (Web): 20 Jul 2018 Downloaded from http://pubs.acs.org on July 21, 2018
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 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 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.
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 38
Environmental Science & Technology
1
Parameterized yields of semi-volatile products from
2
isoprene oxidation under different NOx levels:
3
impacts of chemical aging and wall-loss of reactive
4
gases
5
Li Xing,†,‡,§ Manish Shrivastava,*, ‡ Tzung-May Fu,*, † Pontus Roldin,|| Yun Qian, ‡ Lu Xu,⊥,# Nga
6
L. Ng,⊥,▽ John Shilling,‡ Alla Zelenyuk,‡ and Christopher D. Cappa○
7
†
8
Atmosphere Studies, School of Physics, Peking University, Beijing, 100871, China
9
‡
Pacific Northwest National Laboratory, Richland, Washington, 99352, USA
10
§
Key Lab of Aerosol Chemistry & Physics, State Key Laboratory of Loess and Quaternary
11
Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi’an, 710061, China
12
||
13
⊥
14
Department of Atmospheric and Oceanic Sciences and Laboratory for Climate and Ocean-
Division of Nuclear Physics, Lund University, P.O. Box 118, 221 00 Lund, Sweden School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
15
#
16
California 91125, USA
Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena,
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 38
17
▽
18
30332, USA
19
○
20
California, 95616, USA
21
* Corresponding authors: M. Shrivastava (
[email protected]) and T.-M. Fu
22
(
[email protected])
School of Earth and Atmospheric Sciences, Georgia Institute of Technology, Atlanta, Georgia
Department of Civil and Environmental Engineering, University of California, Davis,
23 24
ACS Paragon Plus Environment
2
Page 3 of 38
Environmental Science & Technology
25
ABSTRACT:
26
We developed a parameterizable box model to empirically derive the yields of semi-volatile
27
products from VOC oxidation using chamber measurements, while explicitly accounting for the
28
multigenerational chemical aging processes (such as the gas-phase fragmentation and
29
functionalization and aerosol-phase oligomerization and photolysis) under different NOx levels
30
and the loss of particles and gases to chamber walls. Using the oxidation of isoprene as an
31
example, we showed that the assumptions regarding the NOx-sensitive, multigenerational aging
32
processes of VOC oxidation products have large impacts on the parameterized product yields
33
and SOA formation. We derived sets of semi-volatile product yields from isoprene oxidation
34
under different NOx levels. However, we stress that these product yields must be used in
35
conjunction with the corresponding multigenerational aging schemes in chemical transport
36
models. As more mechanistic insights regarding SOA formation from VOC oxidation emerge,
37
our box model can be expanded to include more explicit chemical aging processes and help
38
ultimately bridge the gap between the process-based understanding of SOA formation from VOC
39
oxidation and the bulk-yield parameterizations used in chemical transport models.
ACS Paragon Plus Environment
3
Environmental Science & Technology
40
Page 4 of 38
Abstract Art:
41 42 43
ACS Paragon Plus Environment
4
Page 5 of 38
44
Environmental Science & Technology
INTRODUCTION
45
Many volatile organic compounds (VOCs) oxidize in the atmosphere to produce lower
46
volatility products that form secondary organic aerosols (SOA).1-4 Laboratory studies have
47
shown that the oxidation pathways of VOCs are sensitive to ambient NOx levels, and that the
48
molecular complexity and aging of VOC oxidation products have large impacts on the yields and
49
formation timescales of SOA.1-3,5,6 Many current chemical transport models simulated the SOA
50
formation from VOC oxidation using parameterized yields derived from simple empirical
51
theories.7,8 At the same time, an increasing number of those same models are implementing
52
complex chemical aging processes based on new mechanistic insights on SOA formation gained
53
from laboratory studies.9-11 This discrepancy between the implementation of detailed chemical
54
mechanism and the use of simplified parameterized yields in models have so far been overlooked.
55
Here, we presented a new way to parameterize the product yields from VOC oxidation using
56
chamber measurements, taking isoprene oxidation as a specific example. We showed that the
57
explicit accounting of chemical aging and chamber wall-loss processes changed the
58
parameterized yields and volatility distribution of products, which will in turn affect the SOA
59
simulation in chemical transport models.
60
The volatility basis set (VBS) framework, which expands on the two-product model,12 is a
61
widely used empirical approach for modeling SOA formation.13-16 Under the VBS framework, a
62
VOC precursor oxidizes to produce semi-volatile and intermediate volatility (SV/IV) products,
63
which are lumped into n (typically between 4 and 8) bins of effective volatility (represented by
64
bins of effective saturation vapor concentration, C*i, i = 1 to n).8,9,13 These SV/IV products then
65
partition into the aerosol phase according to their respective volatility. To date, most studies
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 38
66
parameterized the effective stoichiometric mass yields (αi) of the SV/IV products by fitting to
67
smog chamber measurements using Eq (1),7,8
68
≡
OA
ROG
= ∑ ∙
∗
Eq (1)
69
where ξ was the aerosol mass fraction (AMF), which was the bulk yield of SOA mass from a
70
reacted VOC precursor measured in the chamber. ∆ROG was the reacted VOC precursor mass,
71
and ∆COA was the total SOA mass formed (adjusted for particle wall-loss). Ci* was the effective
72
saturation vapor concentration of products in the volatility bin i, and COA was the total SOA
73
mass. Typically, a series of experiments were performed by injecting different initial amounts of
74
the VOC precursor into a chamber to be oxidized, and the values of ξ were calculated using the
75
measured time series of the SOA mass and the reacted VOC mass. The product yields in
76
different volatility bins (αi, i = 1 to n) were then determined by fitting to the resulting series of ξ
77
values.8,9,17
78
A frequently overlooked fact is that the yields empirically obtained from fitting Eq (1) were
79
theory-specific. It was assumed that (1) all wall-losses of particles during the measurement were
80
accounted for prior to the fitting; (2) no SV/IV gases were lost to the walls; (3) the volatility of
81
SV/IV products can be lumped into one single, static set of VBS products, which allowed no
82
evolution of their volatility by chemical aging; and (4) the gas/aerosol partitioning of the SV/IV
83
products followed the absorptive partitioning theory. Recent laboratory results have shown that
84
these assumptions pertaining to Eq (1) may not apply.1,18-20 For example, SV/IV gases might be
85
lost to chamber walls during the experiments, leading to a factor of 2 to 10 underestimation of
86
the product yields, especially when the seed-to-chamber surface area ratio was relatively low,21-23
ACS Paragon Plus Environment
6
Page 7 of 38
Environmental Science & Technology
87
although SOA mass yields may not be affected by vapor wall-loss if the SOA formation was
88
governed by fast quasi-equilibrium growth on seed particles.24
89
More importantly, the product yields parameterized from Eq (1) were based only on the total
90
SOA mass formed, with little relevance to the changes of volatility distribution due to aging
91
either within or beyond the chamber experiment timescale (typically a few hours). In reality, the
92
SV/IV products from VOC oxidation undergo mutigenerational aging processes that were highly
93
sensitive to NOx, within the timescale of the chamber experiments and also beyond,1,19 leading to
94
significant changes in volatility and thus SOA yields. For example, under low-NOx (HO2-
95
dominant) conditions, chamber experiments of isoprene oxidation showed that SOA mass began
96
to decline once all isoprene was consumed.1,5 This decline of SOA mass may be due to the
97
photolysis of hydroperoxides in the aerosol phase, or the evaporation of SOA mass once gas-
98
phase compounds were reacted away.5,6 At the same time, the mean SOA oxidation state
99
(represented by 2O/C-H/C)25 continued to increase during an 18-hour experiment under low-NOx
100
conditions1, indicating the formation of highly oxidized (and potentially less volatile) molecules
101
containing one or more peroxy, hydroxyl, or carbonyl function groups.1,26-28 On the other hand,
102
in chamber experiments under high-NOx conditions, SOA mass continued to increase even after
103
all the isoprene has reacted. The mean oxidation state of SOA was higher than that of low-NOx
104
experiment and remained nearly constant during an 18-hour experiment under high-NOx
105
conditions.1
106
Several model studies have tried to incorporate these new mechanistic insights gained from
107
chamber experiments to improve the simple VBS n-product scheme.15,16,23,29 One way to do this
108
was to allow the generation-one VBS products to further oxidize and undergo
109
functionalization/fragmentation at assumed branching ratios.15,16 This approach allowed the
ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 38
110
formation of higher-generation oxidation products with evolving volatility profiles.15,29 In so
111
doing, the great complexity of the multigenerational aging of thousands of oxidation products
112
was reduced and represented with a small set of lumped VBS products. To date, models using
113
such VBS-plus-aging schemes have used products yields empirically obtained from Eq (1) and
114
added subsequent multigenerational aging processes.15,16 However, if multigenerational aging
115
was not included during the parameterization of product yields, it could bias the product yields
116
that serve as initial concentrations for subsequent chemistry. Consequently, there is strong
117
potential for the VBS-plus-aging schemes to overestimate SOA.29,30
118
In this study, we constructed a box model that included multigenerational chemistry of
119
isoprene oxidation and the losses of gas and particles to chamber walls. We then ran the box
120
model with different chemical aging scenarios to fit the measured time series of SOA mass
121
concentrations from chamber experiments under different NOx level. In this way, we investigate
122
the effects of chemical aging and gas wall losses on the parameterized yields and SOA formation
123
from isoprene.
124 125
METHODOLOGY
126
Chamber experiment results from Xu et al.1 We used the time series of SOA mass
127
concentrations in the isoprene oxidation experiments by Xu et al.1, which were conducted in the
128
Pacific Northwest National Laboratory (PNNL) dual 10.6 m3 Teflon environmental chambers
129
with different initial concentrations of isoprene and NOx. The chambers were flushed with pure
130
air prior to each experiment and no seed particles were used. UV lamps initiated the
131
photochemical reactions. A proton transfer reaction mass spectrometer (PTR-MS) measured the
ACS Paragon Plus Environment
8
Page 9 of 38
Environmental Science & Technology
132
concentrations of isoprene and two of its major oxidation products, methacrolein (MACR) and
133
methyl vinyl ketone (MVK). A scanning mobility particle sizer (SMPS) measured the aerosol
134
size distribution between 14.1 and 710.5 nm every five minutes. SOA mass concentrations were
135
calculated from the aerosol volume concentrations measured by the SMPS, assuming particle
136
densities of 1.3 g cm-3 (experiment 2) and 1.4 g cm-3 (experiments 6 and 8). We analyzed results
137
from three experiments with different initial NO/isoprene ratios: ~0, 3.0, and 7.3 (experiments 2,
138
6, and 8 in Xu et al.1) to examine the impact of NOx. In experiment 2, the NOx concentration was
139
below the detection limit (1 ppb) throughout the experiment, such that organic peroxy radicals
140
(RO2) mainly reacted with HO2; this was referred to as an “HO2-dominant” experiment. In
141
experiments 6 and 8, NO was injected into the chambers and RO2 radicals may react with HO2,
142
NO, and NO2. We referred to experiments 6 and 8 as “intermediate-NOx mixed” and “high-NOx
143
mixed” experiments, respectively. Measurements from the other five experiments in Xu et al.1
144
were used for validation (supporting information).
145 146
Box model for SOA formation from isoprene oxidation. We constructed a parameterizable
147
box model to simulate the oxidation of isoprene and the subsequent multigenerational chemistry
148
of the gas and aerosol products. Our box model was built on the framework of the Model for
149
Simulating Aerosol Interactions and Chemistry (MOSAIC),31 which was a discretized-size-bin
150
model originally developed to describe the chemical and microphysical evolution of inorganic
151
aerosols. Here we included a multigenerational aging scheme for isoprene oxidation using the
152
modified VBS aging scheme developed by Shrivastava et al.15 Fig. 1 and Table S1 show the gas-
153
and aerosol-phase reactions and gas-particle mass transfer pathways in the model. Isoprene
154
reacts with OH (reaction R1, rate constant k1=1.0×10-10 cm3 molecule-1 s-1)32 to produce
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 38
155
generation-one semi-volatile gas products (G1,i, i=1 to 4) in four effective saturation vapor
156
concentration bins (C*i, i=1 to 4): 0.1, 1, 10, and 100 µg m-3 at yields αi (i=1 to 4), respectively.
157
Some intermediate-volatility products with saturation vapor concentrations between 103-105 µg
158
m-3 may potentially oxidize to form SOA.33 However, the maximum SOA concentration
159
measured during the experiments by Xu et al.1 was