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):
Author(s) / Creator(s):
Date Published:
Journal Name:
Atmospheric Chemistry and Physics
Page Range / eLocation ID:
9613 to 9640
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. 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
  3. Isoprene emitted by vegetation is an important precursor of secondary organic aerosol (SOA), but the mechanism and yields are uncertain. Aerosol is prevailingly aqueous under the humid conditions typical of isoprene-emitting regions. Here we develop an aqueous-phase mechanism for isoprene SOA formation coupled to a detailed gas-phase isoprene oxidation scheme. The mechanism is based on aerosol reactive uptake coefficients (γ) for water-soluble isoprene oxidation products, including sensitivity to aerosol acidity and nucleophile concentrations. We apply this mechanism to simulation of aircraft (SEAC4RS) and ground-based (SOAS) observations over the southeast US in summer 2013 using the GEOS-Chem chemical transport model. Emissions of nitrogen oxides (NOx  ≡  NO + NO2) over the southeast US are such that the peroxy radicals produced from isoprene oxidation (ISOPO2) react significantly with both NO (high-NOx pathway) and HO2 (low-NOx pathway), leading to different suites of isoprene SOA precursors. We find a mean SOA mass yield of 3.3 % from isoprene oxidation, consistent with the observed relationship of total fine organic aerosol (OA) and formaldehyde (a product of isoprene oxidation). Isoprene SOA production is mainly contributed by two immediate gas-phase precursors, isoprene epoxydiols (IEPOX, 58 % of isoprene SOA) from the low-NOx pathway and glyoxal (28 %) from both low- and high-NOx pathways. This speciation is consistent with observations of IEPOX SOA from SOAS and SEAC4RS. Observations show a strong relationship between IEPOX SOA and sulfate aerosol that we explain as due to the effect of sulfate on aerosol acidity and volume. Isoprene SOA concentrations increase as NOx emissions decrease (favoring the low-NOx pathway for isoprene oxidation), but decrease more strongly as SO2 emissions decrease (due to the effect of sulfate on aerosol acidity and volume). The US Environmental Protection Agency (EPA) projects 2013–2025 decreases in anthropogenic emissions of 34 % for NOx (leading to a 7 % increase in isoprene SOA) and 48 % for SO2 (35 % decrease in isoprene SOA). Reducing SO2 emissions decreases sulfate and isoprene SOA by a similar magnitude, representing a factor of 2 co-benefit for PM2.5 from SO2 emission controls. 
    more » « less
  4. The UNIfied Partitioning-Aerosol phase Reaction (UNIPAR) model was established on the Comprehensive Air quality Model with extensions (CAMx) to process Secondary Organic Aerosol (SOA) formation by capturing multiphase reactions of hydrocarbons (HCs) in regional scales. SOA growth was simulated using a wide range of anthropogenic HCs including ten aromatics and linear alkanes with different carbon-lengths. The atmospheric processes of biogenic HCs (isoprene, terpenes, and sesquiterpene) were simulated for the major oxidation paths (ozone, OH radicals, and nitrate radicals) to predict day and night SOA formation. The UNIPAR model streamlined the multiphase partitioning of the lumping species originating from semi-explicitly predicted gas products and their heterogeneous chemistry to form non-volatile oligomeric species in both organic aerosol and inorganic aqueous phase. The CAMx-UNIPAR model predicted SOA formation at four ground urban sites (San Jose, Sacramento, Fresno, and Bakersfield) in California, United States during wintertime 2018. Overall, the simulated mass concentrations of the total organic matter, consisting of primary OA (POA) and SOA, showed a good agreement with the observations. The simulated SOA mass in the urban areas of California was predominated by alkane and terpene. During the daytime, low-volatile products originating from the autoxidation of long-chain alkanes considerably contributed to the SOA mass. In contrast, a significant amount of nighttime SOA was produced by the reaction of terpene with ozone or nitrate radicals. The spatial distributions of anthropogenic SOA associated with aromatic and alkane HCs were noticeably affected by the southward wind direction owing to the relatively long lifetime of their atmospheric oxidation, whereas those of biogenic SOA were nearly insensitive to wind direction. During wintertime 2018, the impact of inorganic aerosol hygroscopicity on the total SOA budget was not evident because of the small contribution of aromatic and isoprene products that are hydrophilic and reactive in the inorganic aqueous phase. However, an increased isoprene SOA mass was predicted during the wet periods, although its contribution to the total SOA was little. 
    more » « less
  5. Chemical mechanisms play an important role in simulating the atmospheric chemistry of volatile organic compound oxidation. Comparison of mechanism simulations with laboratory chamber data tests our level of understanding of the prevailing chemistry as well as the dynamic processes occurring in the chamber itself. α-Pinene photooxidation is a well-studied system experimentally, for which detailed chemical mechanisms have been formulated. Here, we present the results of simulating low-NO α-pinene photooxidation experiments conducted in the Caltech chamber with the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) under varying concentrations of seed particles and OH levels. Unexpectedly, experiments conducted at low and high OH levels yield the same secondary organic aerosol (SOA) growth, whereas GECKO-A predicts greater SOA growth under high OH levels. SOA formation in the chamber is a result of a competition among the rates of gas-phase oxidation to low-volatility products, wall deposition of these products, and condensation into the aerosol phase. Various processes – such as photolysis of condensed-phase products, particle-phase dimerization, and peroxy radical autoxidation – are explored to rationalize the observations. In order to explain the observed similar SOA growth at different OH levels, we conclude that vapor wall loss in the Caltech chamber is likely of order 10−5 s−1, consistent with previous experimental measurements in that chamber. We find that GECKO-A tends to overpredict the contribution to SOA of later-generation oxidation products under high-OH conditions. Moreover, we propose that autoxidation may alternatively resolve some or all of the measurement–model discrepancy, but this hypothesis cannot be confirmed until more explicit mechanisms are established for α-pinene autoxidation. The key role of the interplay among oxidation rate, product volatility, and vapor–wall deposition in chamber experiments is illustrated. 
    more » « less