Parameterized Yields of Semivolatile Products from Isoprene

Jul 20, 2018 - As more mechanistic insights regarding SOA formation from VOC oxidation emerge, our box model can be expanded to include more explicit ...
0 downloads 0 Views 1MB Size
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