Juliane L. Fry

Biomass burning dominates brown carbon absorption in the rural southeastern United States: Biomass burning dominates brown carbon

Brown carbon aerosol consists of light-absorbing organic particulate matter with wavelength-dependent absorption. Aerosol optical extinction, absorption, size distributions, and chemical composition were measured in rural Alabama during summer 2013. The field site was well located to examine sources of brown carbon aerosol, with influence by high biogenic organic aerosol concentrations, pollution from two nearby cities, and biomass burning aerosol. We report the optical closure between measured dry aerosol extinction at 365 nm and calculated extinction from composition and size distribution, showing agreement within experiment uncertainties. We find that aerosol optical extinction is dominated by scattering, with single-scattering albedo values of 0.94 ± 0.02. Black carbon aerosol accounts for 91 ± 9% of the total carbonaceous aerosol absorption at 365 nm, while organic aerosol accounts for 9 ± 9%. The majority of brown carbon aerosol mass is associated with biomass burning, with smaller contributions from biogenically derived secondary organic aerosol.

Secondary Organic Aerosol Formation and Organic Nitrate Yield from NO3 Oxidation of Biogenic Hydrocarbons

The secondary organic aerosol (SOA) mass yields from NO3 oxidation of a series of biogenic volatile organic compounds (BVOCs), consisting of five monoterpenes and one sesquiterpene (α-pinene, β-pinene, Δ-3-carene, limonene, sabinene, and β-caryophyllene), were investigated in a series of continuous flow experiments in a 10 m3 indoor Teflon chamber. By making in situ measurements of the nitrate radical and employing a kinetics box model, we generate time-dependent yield curves as a function of reacted BVOC. SOA yields varied dramatically among the different BVOCs, from zero for α-pinene to 38–65% for Δ-3-carene and 86% for β-caryophyllene at mass loading of 10 μg m–3, suggesting that model mechanisms that treat all NO3 + monoterpene reactions equally will lead to errors in predicted SOA depending on each location’s mix of BVOC emissions. In most cases, organonitrate is a dominant component of the aerosol produced, but in the case of α-pinene, little organonitrate and no aerosol is formed.

Nitrate radicals and biogenic volatile organic compounds: oxidation, mechanisms, and organic aerosol

Oxidation of biogenic volatile organic compounds (BVOC) by the nitrate radical (NO3) represents one of the important interactions between anthropogenic emissions related to combustion and natural emissions from the biosphere. This interaction has been recognized for more than 3 decades, during which time a large body of research has emerged from laboratory, field, and modeling studies. NO3-BVOC reactions influence air quality, climate and visibility through regional and global budgets for reactive nitrogen (particularly organic nitrates), ozone, and organic aerosol. Despite its long history of research and the significance of this topic in atmospheric chemistry, a number of important uncertainties remain. These include an incomplete understanding of the rates, mechanisms, and organic aerosol yields for NO3-BVOC reactions, lack of constraints on the role of heterogeneous oxidative processes associated with the NO3 radical, the difficulty of characterizing the spatial distributions of BVOC and NO3 within the poorly mixed nocturnal atmosphere, and the challenge of constructing appropriate boundary layer schemes and non-photochemical mechanisms for use in state-of-the-art chemical transport and chemistry–climate models. This review is the result of a workshop of the same title held at the Georgia Institute of Technology in June 2015. The first half of the review summarizes the current literature on NO3-BVOC chemistry, with a particular focus on recent advances in instrumentation and models, and in organic nitrate and secondary organic aerosol (SOA) formation chemistry. Building on this current understanding, the second half of the review outlines impacts of NO3-BVOC chemistry on air quality and climate, and suggests critical research needs to better constrain this interaction to improve the predictive capabilities of atmospheric models.

Modeling the Current and Future Roles of Particulate Organic Nitrates in the Southeastern United States

Organic nitrates are an important aerosol constituent in locations where biogenic hydrocarbon emissions mix with anthropogenic NOx sources. While regional and global chemical transport models may include a representation of organic aerosol from monoterpene reactions with nitrate radicals (the primary source of particle-phase organic nitrates in the Southeast United States), secondary organic aerosol (SOA) models can underestimate yields. Furthermore, SOA parametrizations do not explicitly take into account organic nitrate compounds produced in the gas phase. In this work, we developed a coupled gas and aerosol system to describe the formation and subsequent aerosol-phase partitioning of organic nitrates from isoprene and monoterpenes with a focus on the Southeast United States. The concentrations of organic aerosol and gas-phase organic nitrates were improved when particulate organic nitrates were assumed to undergo rapid (τ = 3 h) pseudohydrolysis resulting in nitric acid and nonvolatile secondary organic aerosol. In addition, up to 60% of less oxidized-oxygenated organic aerosol (LO-OOA) could be accounted for via organic nitrate mediated chemistry during the Southern Oxidants and Aerosol Study (SOAS). A 25% reduction in nitrogen oxide (NO + NO2) emissions was predicted to cause a 9% reduction in organic aerosol for June 2013 SOAS conditions at Centreville, Alabama.

Cis-cis and trans-perp HOONO: Action spectroscopy and isomerization kinetics

The weakly bound HOONO product of the OH+NO2+M reaction is studied using the vibrational predissociation that follows excitation of the first OH overtone (2nu1). We observe formation of both cis-cis and trans-perp conformers of HOONO. The trans-perp HOONO 2nu1 band is observed under thermal (223-238 K) conditions at 6971 cm(-1). We assign the previously published (warmer temperature) HOONO spectrum to the 2nu1 band at 6365 cm(-1) and 2nu1-containing combination bands of the cis-cis conformer of HOONO. The band shape of the trans-perp HOONO spectrum is in excellent agreement with the predicted rotational contour based on previous experimental and theoretical results, but the apparent origin of the cis-cis HOONO spectrum at 6365 cm(-1) is featureless and significantly broader, suggesting more rapid intramolecular vibrational redistribution or predissociation in the latter isomer. The thermally less stable trans-perp HOONO isomerizes rapidly to cis-cis HOONO with an experimentally determined lifetime of 39 ms at 233 K at 13 hPa (in a buffer gas of predominantly Ar). The temperature dependence of the trans-perp HOONO lifetime in the range 223-238 K yields an isomerization barrier of 33+/-12 kJ/mol. New ab initio calculations of the structure and vibrational mode frequencies of the transition state perp-perp HOONO are performed using the coupled cluster singles and doubles with perturbative triples [CCSD(T)] model, using a correlation consistent polarized triple zeta basis set (cc-pVTZ). The energetics of cis-cis, trans-perp, and perp-perp HOONO are also calculated at this level [CCSD(T)/cc-pVTZ] and with a quadruple zeta basis set using the structure determined at the triple zeta basis set [CCSD(T)/cc-pVQZ//CCSD(T)/cc-pVTZ]. These calculations predict that the anti form of perp-perp HOONO has an energy of DeltaE0=42.4 kJ/mol above trans-perp HOONO, corresponding to an activation enthalpy of DeltaH298 (double dagger 0)=41.1 kJ/mol. These results are in good agreement with statistical simulations based on a model developed by Golden, Barker, and Lohr. The simulated isomerization rates match the observed decay rates when modeled with a trans-perp to cis-cis HOONO isomerization barrier of 40.8 kJ/mol and a strong collision model. The quantum yield of cis-cis HOONO dissociation to OH and NO2 is also calculated as a function of photon excitation energy in the range 3500-7500 cm(-1), assuming D0=83 kJ/mol. The quantum yield is predicted to vary from 0.15 to 1 over the observed spectrum at 298 K, leading to band intensities in the action spectrum that are highly temperature dependent; however, the observed relative band strengths in the cis-cis HOONO spectrum do not change substantially with temperature over the range 193-273 K. Semiempirical calculations of the oscillator strengths for 2nu1(cis-cis HOONO) and 2nu1(trans-perp HOONO) are performed using (1) a one-dimensional anharmonic model and (2) a Morse oscillator model for the OH stretch, and ab initio dipole moment functions calculated using Becke, Lee, Yang, and Parr density functional theory (B3LYP), Møller-Plesset pertubation theory truncated at the second and third order (MP2 and MP3), and quadratic configuration interaction theory using single and double excitations (QCISD). The QCISD level calculated ratio of 2nu1 oscillator strengths of trans-perp to cis-cis HOONO is 3.7:1. The observed intensities indicate that the concentration of trans-perp HOONO early in the OH+NO2 reaction is significantly greater than predicted by a Boltzmann distribution, consistent with statistical predictions of high initial yields of trans-perp HOONO from the OH+NO2+M reaction. In the atmosphere, trans-perp HOONO will isomerize nearly instantaneously to cis-cis HOONO. Loss of HOONO via photodissociation in the near-IR limits the lifetime of cis-cis HOONO during daylight to less than 45 h, other loss mechanisms will reduce the lifetime further.

Role of OH-stretch/torsion coupling and quantum yield effects in the first OH overtone spectrum of cis-cis HOONO

A joint theoretical and experimental investigation is undertaken to study the effects of OH-stretch/HOON torsion coupling and of quantum yield on the previously reported first overtone action spectrum of cis-cis HOONO (peroxynitrous acid). The minimum energy path along the HOON dihedral angle is computed at the coupled cluster singles and doubles with perturbative triples level with correlation consistent polarized quadruple zeta basis set, at the structure optimized using the triple zeta basis set (CCSD(T)/cc-pVQZ//CCSD(T)/cc-pVTZ). The two-dimensional ab initio potential energy and dipole moment surfaces for cis-cis HOONO are calculated as functions of the HOON torsion and OH bond length about the minimum energy path at the CCSD(T)/cc-pVTZ and QCISD/AUG-cc-pVTZ (QCISD-quadratic configuration interaction with single and double excitation and AUG-augmented with diffuse functions) level of theory/basis, respectively. The OH-stretch vibration depends strongly on the torsional angle, and the torsional potential possesses a broad shelf at approximately 90 degrees , the cis-perp conformation. The calculated electronic energies and dipoles are fit to simple functional forms and absorption spectra in the region of the OH fundamental and first overtone are calculated from these surfaces. While the experimental and calculated spectra of the OH fundamental band are in good agreement, significant differences in the intensity patterns are observed between the calculated absorption spectrum and the measured action spectrum in the 2nu(OH) region. These differences are attributed to the fact that several of the experimentally accessible states do not have sufficient energy to dissociate to OH+NO(2) and therefore are not detectable in an action spectrum. Scaling of the intensities of transitions to these states, assuming D(0)=82.0 kJ/mol, is shown to produce a spectrum that is in good agreement with the measured action spectrum. Based on this agreement, we assign two of the features in the spectrum to Deltan=0 transitions (where n is the HOON torsion quantum number) that are blue shifted relative to the origin band, while the large peak near 7000 cm(-1) is assigned to a series of Deltan=+1 transitions, with predominant contributions from torsionally excited states with substantial cis-perp character. The direct absorption spectrum of cis-cis HOONO (6300-6850 cm(-1)) is recorded by cavity ringdown spectroscopy in a discharge flow cell. A single band of HOONO is observed at 6370 cm(-1) and is assigned as the origin of the first OH overtone of cis-cis HOONO. These results imply that the origin band is suppressed by over an order of magnitude in the action spectrum, due to a reduced quantum yield. The striking differences between absorption and action spectra are correctly predicted by the calculations.

Testing Atmospheric Oxidation in an Alabama Forest

The chemical species emitted by forests create complex atmospheric oxidation chemistry and influence global atmospheric oxidation capacity and climate. The Southern Oxidant and Aerosol Study (SOAS) provided an opportunity to test the oxidation chemistry in a forest where isoprene is the dominant biogenic volatile organic compound. Hydroxyl (OH) and hydroperoxyl (HO_2) radicals were two of the hundreds of atmospheric chemical species measured, as was OH reactivity (the inverse of the OH lifetime). OH was measured by laser-induced fluorescence (LIF) and by taking the difference in signals without and with an OH scavenger that was added just outside the instrument’s pinhole inlet. To test whether the chemistry at SOAS can be simulated by current model mechanisms, OH and HO_2 were evaluated with a box model using two chemical mechanisms: Master Chemical Mechanism, version 3.2 (MCMv3.2), augmented with explicit isoprene chemistry and MCMv3.3.1. Measured and modeled OH peak at about 10^6 cm^(−3) and agree well within combined uncertainties. Measured and modeled HO_2 peak at about 27 pptv and also agree well within combined uncertainties. Median OH reactivity cycled between about 11 s^(−1) at dawn and about 26 s^(−1) during midafternoon. A good test of the oxidation chemistry is the balance between OH production and loss rates using measurements; this balance was observed to within uncertainties. These SOAS results provide strong evidence that the current isoprene mechanisms are consistent with measured OH and HO_2 and, thus, capture significant aspects of the atmospheric oxidation chemistry in this isoprene-rich forest.

Vibrational overtone initiated unimolecular dissociation of HOCH2OOH and HOCD2OOH : Evidence for mode selective behavior

The vibrational overtone induced unimolecular dissociation of HMHP (HOCH(2)OOH) and HMHP-d(2) (HOCD(2)OOH) into OH and HOCH(2)O (HOCD(2)O) fragments is investigated in the region of the 4nu(OH) and 5nu(OH) bands. The unimolecular dissociation rates in the threshold region, corresponding to the 4nu(OH) band, exhibit measurable differences associated with excitation of the OH stretch of the alcohol versus the peroxide functional group, with the higher energy alcohol OH stretching state exhibiting a slower dissociation rate compared to the lower energy peroxide OH stretch in both HMHP and HMHP-d(2). Predictions using the Rice-Ramsperger-Kassel-Marcus theory give rates that are in reasonably good agreement with the measured dissociation rate for the alcohol OH stretch but considerably differ from the measured rates for the peroxide OH stretch in both isotopomers. The present results are interpreted as suggesting that the extent of intramolecular vibrational energy redistribution (IVR) is different for the two OH stretching states associated with the two functional groups in HMHP, with IVR being substantially less complete for the peroxide OH stretch. Analysis of the OH fragment product state distributions in conjunction with phase-space theory simulation gives a D(0) value of 38+/-0.7 kcal/mole for breaking the peroxide bond in HMHP.

Coupling of organic and inorganic aerosol systems and the effect on gas–particle partitioning in the southeastern US

Several models were used to describe the partitioning of ammonia, water, and organic compounds between the gas and particle phases for conditions in the southeastern US during summer 2013. Existing equilibrium models and frameworks were found to be sufficient, although additional improvements in terms of estimating pure-species vapor pressures are needed. Thermodynamic model predictions were consistent, to first order, with a molar ratio of ammonium to sulfate of approximately 1.6 to 1.8 (ratio of ammonium to 2× sulfate, RN/2S ≈ 0.8 to 0.9) with approximately 70% of total ammonia and ammonium (NHx) in the particle. Southeastern Aerosol Research and Characterization Network (SEARCH) gas and aerosol and Southern Oxidant and Aerosol Study (SOAS) Monitor for AeRosols and Gases in Ambient air (MARGA) aerosol measurements were consistent with these conditions. CMAQv5.2 regional chemical transport model predictions did not reflect these conditions due to a factor of 3 overestimate of the nonvolatile cations. In addition, gas-phase ammonia was overestimated in the CMAQ model leading to an even lower fraction of total ammonia in the particle. Chemical Speciation Network (CSN) and aerosol mass spectrometer (AMS) measurements indicated less ammonium per sulfate than SEARCH and MARGA measurements and were inconsistent with thermodynamic model predictions. Organic compounds were predicted to be present to some extent in the same phase as inorganic constituents, modifying their activity and resulting in a decrease in [H+]air (H+ in μgm−3 air), increase in ammonia partitioning to the gas phase, and increase in pH compared to complete organic vs. inorganic liquid–liquid phase separation. In addition, accounting for nonideal mixing modified the pH such that a fully interactive inorganic–organic system had a pH roughly 0.7 units higher than predicted using traditional methods (pH = 1.5 vs. 0.7). Particle-phase interactions of organic and inorganic compounds were found to increase partitioning towards the particle phase (vs. gas phase) for highly oxygenated (O : C≥0.6) compounds including several isoprene-derived tracers as well as levoglu-cosan but decrease particle-phase partitioning for low O: C, monoterpene-derived species.

Alkoxy Radical Bond Scissions Explain the Anomalously Low Secondary Organic Aerosol and Organonitrate Yields From α-Pinene + NO3

Oxidation of monoterpenes (C10H16) by nitrate radicals (NO3) constitutes an important source of atmospheric secondary organic aerosol (SOA) and organonitrates. However, knowledge of the mechanisms of their formation is incomplete and differences in yields between similar monoterpenes are poorly understood. In particular, yields of SOA and organonitrates from α-pinene + NO3 are low, while those from Δ3-carene + NO3 are high. Using computational methods, we suggest that bond scission of the nitrooxy alkoxy radicals from Δ3-carene lead to the formation of reactive keto-nitrooxy-alkyl radicals, which retain the nitrooxy moiety and can undergo further reactions to form SOA. By contrast, bond scissions of the nitrooxy alkoxy radicals from α-pinene lead almost exclusively to the formation of the relatively unreactive and volatile product pinonaldehyde (C10H16O2), thereby limiting organonitrate and SOA formation. This hypothesis is supported by laboratory experiments that quantify products of the reaction of α-pinene + NO3 under atmospherically relevant conditions.