skip to main content


Title: A new model mechanism for atmospheric oxidation of isoprene: global effects on oxidants, nitrogen oxides, organic products, and secondary organic aerosol
Abstract. Atmospheric oxidation of isoprene, the most abundantly emitted non-methane hydrocarbon, affects the abundances of ozone (O3), the hydroxyl radical (OH), nitrogen oxide radicals (NOx), carbon monoxide (CO), oxygenated and nitrated organic compounds, and secondary organic aerosol (SOA). We analyze these effects in box models and in the global GEOS-Chem chemical transport model using the new reduced Caltech isoprene mechanism (RCIM) condensed from a recently developed explicit isoprene oxidation mechanism. We find many similarities with previous global models of isoprene chemistry along with a number of important differences. Proper accounting of the isomer distribution of peroxy radicals following the addition of OH and O2 to isoprene influences the subsequent distribution of products, decreasing in particular the yield of methacrolein and increasing the capacity of intramolecular hydrogen shifts to promptly regenerate OH. Hydrogen shift reactions throughout the mechanism lead to increased OH recycling, resulting in less depletion of OH under low-NO conditions than in previous mechanisms. Higher organonitrate yields and faster tertiary nitrate hydrolysis lead to more efficient NOx removal by isoprene and conversion to inorganic nitrate. Only 20 % of isoprene-derived organonitrates (excluding peroxyacyl nitrates) are chemically recycled to NOx. The global yield of formaldehyde from isoprene is 22 % per carbon and less sensitive to NO than in previous mechanisms. The global molar yield of glyoxal is 2 %, much lower than in previous mechanisms because of deposition and aerosol uptake of glyoxal precursors. Global production of isoprene SOA is about one-third from each of the following: isoprene epoxydiols (IEPOX), organonitrates, and tetrafunctional compounds. We find a SOA yield from isoprene of 13 % per carbon, much higher than commonly assumed in models and likely offset by SOA chemical loss. We use the results of our simulations to further condense RCIM into a mini Caltech isoprene mechanism (Mini-CIM) for less expensive implementation in atmospheric models, with a total size (108 species, 345 reactions) comparable to currently used mechanisms.  more » « less
Award ID(s):
1643217
NSF-PAR ID:
10197175
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Atmospheric Chemistry and Physics
Volume:
19
Issue:
14
ISSN:
1680-7324
Page Range / eLocation ID:
9613 to 9640
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The daytime oxidation of biogenic hydrocarbons is attributed to both OH radicals and O3, while nighttime chemistry is dominated by the reaction with O3 and NO3 radicals. Here, the diurnal pattern of Secondary Organic Aerosol (SOA) originating from biogenic hydrocarbons was intensively evaluated under varying environmental conditions (temperature, humidity, sunlight intensity, NOx levels, and seed conditions) by using the UNIfied Partitioning Aerosol phase Reaction (UNIPAR) model, which comprises multiphase gas-particle partitioning and in-particle chemistry. The oxidized products of three different hydrocarbons (isoprene, α-pinene, and β-caryophyllene) were predicted by using near explicit gas mechanisms for four different oxidation paths (OH, O3, NO3, and O(3P)) during day and night. The gas mechanisms implemented the Master Chemical Mechanism (MCM v3.3.1), the reactions that formed low volatility products via peroxy radical (RO2) autoxidation, and self- and cross-reactions of nitrate-origin RO2. In the model, oxygenated products were then classified into volatility-reactivity base lumping species, which were dynamically constructed under varying NOx levels and aging scales. To increase feasibility, the UNIPAR model that equipped mathematical equations for stoichiometric coefficients and physicochemical parameters of lumping species was integrated with the SAPRC gas mechanism. The predictability of the UNIPAR model was demonstrated by simulating chamber-generated SOA data under varying environments day and night. Overall, the SOA simulation decoupled to each oxidation path indicated that the nighttime isoprene SOA formation was dominated by the NO3-driven oxidation, regardless of NOx levels. However, the oxidation path to produce the nighttime α-pinene SOA gradually transited from the NO3-initiated reaction to ozonolysis as NOx levels decreased. For daytime SOA formation, both isoprene and α-pinene were dominated by the OH-radical initiated oxidation. The contribution of the O(3P) path to all biogenic SOA formation was negligible in daytime. Sunlight during daytime promotes the decomposition of oxidized products via photolysis and thus, reduces SOA yields. Nighttime α-pinene SOA yields were significantly higher than daytime SOA yields, although the nighttime α-pinene SOA yields gradually decreased with decreasing NOx levels. For isoprene, nighttime chemistry yielded higher SOA mass than daytime at the higher NOx level (isoprene/NOx > 5 ppbC/ppb). The daytime isoprene oxidation at the low NOx level formed epoxy-diols that significantly contributed SOA formation via heterogeneous chemistry. For isoprene and α-pinene, daytime SOA yields gradually increased with decreasing NOx levels. The daytime SOA produced more highly oxidized multifunctional products and thus, it was generally more sensitive to the aqueous reactions than the nighttime SOA. β-Caryophyllene, which rapidly oxidized and produced SOA with high yields, showed a relatively small variation in SOA yields from changes in environmental conditions (i.e., NOx levels, seed conditions, and diurnal pattern), and its SOA formation was mainly attributed to ozonolysis day and night. To mimic the nighttime α-pinene SOA formation under the polluted urban atmosphere, α-pinene SOA formation was simulated in the presence of gasoline fuel. The simulation suggested the growth of α-pinene SOA in the presence of gasoline fuel gas by the enhancement of the ozonolysis path under the excess amount of ozone, which is typical in urban air. We concluded that the oxidation of the biogenic hydrocarbon with O3 or NO3 radicals is a source to produce a sizable amount of nocturnal SOA, despite of the low emission at night. 
    more » « less
  2. null (Ed.)
    Abstract. Reactions of the hydroxyl (OH) and peroxy (HO2 and RO2) radicals playa central role in the chemistry of the atmosphere. In addition to controlling the lifetimes ofmany trace gases important to issues of global climate change, OH radical reactionsinitiate the oxidation of volatile organic compounds (VOCs) which can lead to the production ofozone and secondary organic aerosols in the atmosphere. Previous measurements of these radicalsin forest environments characterized by high mixing ratios of isoprene and low mixing ratios ofnitrogen oxides (NOx) (typically less than 1–2 ppb) have shown seriousdiscrepancies with modeled concentrations. These results bring into question our understanding ofthe atmospheric chemistry of isoprene and other biogenic VOCs under low NOxconditions. During the summer of 2015, OH and HO2 radical concentrations, as well as totalOH reactivity, were measured using laser-induced fluorescence–fluorescence assay by gasexpansion (LIF-FAGE) techniques as part of the Indiana Radical Reactivity and Ozone productioN InterComparison (IRRONIC). This campaign took place in a forested area near Indiana University's Bloomington campus which is characterized by high mixing ratios of isoprene (average daily maximum ofapproximately 4 ppb at 28 ∘C) and low mixing ratios of NO (diurnal averageof approximately 170 ppt). Supporting measurements of photolysis rates, VOCs,NOx, and other species were used to constrain a zero-dimensional box model basedon the Regional Atmospheric Chemistry Mechanism (RACM2) and the Master Chemical Mechanism (MCM 3.2),including versions of the Leuven isoprene mechanism (LIM1) for HOx regeneration(RACM2-LIM1 and MCM 3.3.1). Using an OH chemical scavenger technique, the study revealed thepresence of an interference with the LIF-FAGE measurements of OH that increased with bothambient concentrations of ozone and temperature with an average daytime maximum equivalentOH concentration of approximately 5×106 cm−3. Subtraction of theinterference resulted in measured OH concentrations of approximately4×106 cm−3 (average daytime maximum) that were in better agreement with modelpredictions although the models underestimated the measurements in the evening. The addition ofversions of the LIM1 mechanism increased the base RACM2 and MCM 3.2 modeled OH concentrationsby approximately 20 % and 13 %, respectively, with the RACM2-LIM1 mechanism providing thebest agreement with the measured concentrations, predicting maximum daily OH concentrationsto within 30 % of the measured concentrations. Measurements of HO2 concentrationsduring the campaign (approximately a 1×109 cm−3 average daytime maximum)included a fraction of isoprene-based peroxy radicals(HO2*=HO2+αRO2) and were found to agree with modelpredictions to within 10 %–30 %. On average, the measured reactivity was consistent with thatcalculated from measured OH sinks to within 20 %, with modeled oxidation productsaccounting for the missing reactivity, however significant missing reactivity (approximately40 % of the total measured reactivity) was observed on some days. 
    more » « less
  3. We use a 0-D photochemical box model and a 3-D global chemistry-climate model, combined with observations from the NOAA Southeast Nexus (SENEX) aircraft campaign, to understand the sources and sinks of glyoxal over the Southeast United States. Box model simulations suggest a large difference in glyoxal production among three isoprene oxidation mechanisms (AM3ST, AM3B, and MCM v3.3.1). These mechanisms are then implemented into a 3-D global chemistry-climate model. Comparison with field observations shows that the average vertical profile of glyoxal is best reproduced by AM3ST with an effective reactive uptake coefficient γglyx of 2 × 10-3, and AM3B without heterogeneous loss of glyoxal. The two mechanisms lead to 0-0.8 µg m-3 secondary organic aerosol (SOA) from glyoxal in the boundary layer of the Southeast U.S. in summer. We consider this to be the lower limit for the contribution of glyoxal to SOA, as other sources of glyoxal other than isoprene are not included in our model. In addition, we find that AM3B shows better agreement on both formaldehyde and the correlation between glyoxal and formaldehyde (RGF = [GLYX]/[HCHO]), resulting from the suppression of δ-isoprene peroxy radicals (δ-ISOPO2). We also find that MCM v3.3.1 may underestimate glyoxal production from isoprene oxidation, in part due to an underestimated yield from the reaction of IEPOX peroxy radicals (IEPOXOO) with HO2. Our work highlights that the gas-phase production of glyoxal represents a large uncertainty in quantifying its contribution to SOA. 
    more » « less
  4. Isoprene is one of the most common biogenic volatile organic compounds (BVOC) in the atmosphere, produced by many plants. Isoprene undergoes oxidation to form gaseous isoprene epoxydiols (IEPOX) under low-NOx conditions, which can lead to the formation of secondary organic aerosol (SOA) particles. SOA-containing particles affect climate by scattering and absorbing solar radiation or acting as cloud condensation nuclei (CCN). High concentrations of SOA are also associated with adverse health impacts in people. While in the atmosphere, IEPOX SOA particles continue to undergo reactions with atmospheric oxidants, including hydroxyl radical (OH). To isolate and probe this process, we studied atmospheric chemical processes in an aerosol chamber to better understand the evolution of heterogeneous OH oxidation of IEPOX-derived SOA particles. Since very little is understood about the structural and spectroscopic properties because of the complexity of their many sources and atmospheric processing, individual particle measurements are necessary to provide better understanding of the composition of IEPOX SOA. We injected particles composed of mixtures of ammonium sulfate and sulfuric acid across a range of acidities(PH = 0.5 – 2.5) and gas-phase IEPOX into the chamber to generate SOA. The SOA particles were then sent to an oxidation flow reactor, and exposed to different OH concentrations representative of aging of a number of days. We kept relative humidity (RH) constant at ~65%, the temperature was ~23 °C, and levels of oxidation were controlled by adjusting lamp intensity. After oxidized SOA was impacted on quartz substrates, we used single-particle Raman microspectroscopy to identify their functional group compositions. From the Raman vibrational spectra of submicron particles (~500-1000 nm aerodynamic diameter), we observed a distinct difference in core-shell morphology and composition: an organic outer layer and an aqueous-inorganic core. The core also has significantly more CH-stretch than the shell. Small changes were also observed with increasing oxidation, which are important to consider when predicting SOA particle evolution in the atmosphere. 
    more » « less
  5. Abstract. While camphene is one of the dominant monoterpenesmeasured in biogenic and pyrogenic emission samples, oxidation of camphenehas not been well-studied in environmental chambers and very little is knownabout its potential to form secondary organic aerosol (SOA). The lack ofchamber-derived SOA data for camphene may lead to significant uncertaintiesin predictions of SOA from oxidation of monoterpenes using existingparameterizations when camphene is a significant contributor to totalmonoterpenes. Therefore, to advance the understanding of camphene oxidationand SOA formation and to improve representation of camphene in air qualitymodels, a series of experiments was performed in the University ofCalifornia Riverside environmental chamber to explore camphene SOA massyields and properties across a range of chemical conditions atatmospherically relevant OH concentrations. The experimental results werecompared with modeling simulations obtained using two chemically detailedbox models: Statewide Air Pollution Research Center (SAPRC) and Generatorfor Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A).SOA parameterizations were derived from the chamber data using both thetwo-product and volatility basis set (VBS) approaches. Experiments performedwith added nitrogen oxides (NOx) resulted in higher SOA mass yields (upto 64 %) than experiments performed without added NOx (up to 28 %).In addition, camphene SOA mass yields increased with SOA mass (Mo) atlower mass loadings, but a threshold was reached at higher mass loadings inwhich the SOA mass yields no longer increased with Mo. SAPRC modelingof the chamber studies suggested that the higher SOA mass yields at higherinitial NOx levels were primarily due to higher production of peroxyradicals (RO2) and the generation of highly oxygenated organicmolecules (HOMs) formed through unimolecular RO2 reactions. SAPRCpredicted that in the presence of NOx, camphene RO2 reacts with NOand the resultant RO2 undergoes hydrogen (H)-shift isomerizationreactions; as has been documented previously, such reactions rapidly addoxygen and lead to products with very low volatility (i.e., HOMs). The endproducts formed in the presence of NOx have significantly lowervolatilities, and higher O : C ratios, than those formed by initial campheneRO2 reacting with hydroperoxyl radicals (HO2) or other RO2.Further analysis reveals the existence of an extreme NOx regime, whereinthe SOA mass yield can be suppressed again due to high NO / HO2 ratios.Moreover, particle densities were found to decrease from 1.47 to 1.30 g cm−3 as [HC]0 / [NOx]0 increased and O : C decreased. Theobserved differences in SOA mass yields were largely explained by thegas-phase RO2 chemistry and the competition between RO2+HO2, RO2+ NO, RO2+ RO2, and RO2 autoxidationreactions. 
    more » « less