Enhancements in Ammonia and Methane from Agricultural Sources in the Northeastern Colorado Front Range Using Observations from a Small Research Aircraft.

Quantifying ammonia (NH3) to methane (CH4) enhancement ratios from agricultural sources is important for understanding air pollution and nitrogen deposition. The northeastern Colorado Front Range is home to concentrated animal feeding operations (CAFOs) that produce large emissions of NH3 and CH4. Isolating enhancements of NH3 and CH4 in this region due to agriculture is complicated because CAFOs are often located within regions of oil and natural gas (O&NG) extraction that are a major source of CH4 and other alkanes. Here, we utilize a small research aircraft to collect in situ 1 Hz measurements of gas-phase NH3, CH4, and ethane (C2H6) downwind of CAFOs during three flights conducted in November 2019. Enhancements in NH3 and CH4 are distinguishable up to 10 km downwind of CAFOs with the most concentrated portions of the plumes typically below 0.25 km AGL. We demonstrate that NH3 and C2H6 can be jointly used to separate near-source enhancements in CH4 from agriculture and O&NG. Molar enhancement ratios of NH3 to CH4 are quantified for individual CAFOs in this region, and they range from 0.8 to 2.7 ppbv ppbv-1. A multivariate regression model produces enhancement ratios and quantitative regional source contributions that are consistent with prior studies.

Cows as canaries: The effects of ambient air pollution exposure on milk production and somatic cell count in dairy cows.

Exposure to air pollution, including criteria pollutants such as fine particulate matter (PM2.5) and ozone (O3), has been associated with morbidity and mortality in mammals. As a genetically homogenous population that is closely monitored for health, dairy cattle present a unique opportunity to assess the association between changes in air pollution and mammalian health. Milk yield decreases in the summer if temperature and humidity, measured by the Temperature Humidity Index (THI). As O3 levels increase with warmer temperatures, and summer PM2.5 may increase with wildfire smoke, dairy cows may serve as a useful sentinel species to evaluate subacute markers of inflammation and metabolic output and ambient pollution. Over two years, we assessed summertime O3 and PM2.5 concentrations from local US EPA air quality monitors into an auto-regressive mixed model of the association between THI and daily milk production data and bulk tank somatic cell count (SCC). In unadjusted models, a 10 unit increase THI was associated with 28,700 cells/mL (95% CI: 17,700, 39,690) increase in SCC. After controlling for ambient air pollutants, THI was associated with a 14,500 SCC increase (95% CI: 3,400, 25,680), a 48% decrease in effect compared to the crude model. Further, in fully adjusted models, PM2.5 was associated with a 105,500 cells/mL (95% CI: 90,030, 121,050) increase in SCC. Similar results were found for milk production. Results were amplified when high PM2.5 days (95th percentile of observed values) associated with wildfire smoke were removed from the analyses. Our results support the hypothesis that PM2.5 confounds the relationships between THI and milk yield and somatic cell count. The results of this study can be used to inform strategies for intervention to mitigate these impacts at the dairy level and potentially contribute to a model where production animals can act as air quality sentinels.

Anthropogenic enhancements to production of highly oxygenated molecules from autoxidation.

Atmospheric oxidation of natural and anthropogenic volatile organic compounds (VOCs) leads to secondary organic aerosol (SOA), which constitutes a major and often dominant component of atmospheric fine particulate matter (PM2.5). Recent work demonstrates that rapid autoxidation of organic peroxy radicals (RO2) formed during VOC oxidation results in highly oxygenated organic molecules (HOM) that efficiently form SOA. As NOx emissions decrease, the chemical regime of the atmosphere changes to one in which RO2 autoxidation becomes increasingly important, potentially increasing PM2.5, while oxidant availability driving RO2 formation rates simultaneously declines, possibly slowing regional PM2.5 formation. Using a suite of in situ aircraft observations and laboratory studies of HOM, together with a detailed molecular mechanism, we show that although autoxidation in an archetypal biogenic VOC system becomes more competitive as NOx decreases, absolute HOM production rates decrease due to oxidant reductions, leading to an overall positive coupling between anthropogenic NOx and localized biogenic SOA from autoxidation. This effect is observed in the Atlanta, Georgia, urban plume where HOM is enhanced in the presence of elevated NO, and predictions for Guangzhou, China, where increasing HOM-RO2 production coincides with increases in NO from 1990 to 2010. These results suggest added benefits to PM2.5 abatement strategies come with NOx emission reductions and have implications for aerosol-climate interactions due to changes in global SOA resulting from NOx interactions since the preindustrial era.

Quantifying Methane and Ozone Precursor Emissions from Oil and Gas Production Regions across the Contiguous US.

We present an updated fuel-based oil and gas (FOG) inventory with estimates of nitrogen oxide (NOx) emissions from oil and natural gas production in the contiguous US (CONUS). We compare the FOG inventory with aircraft-derived ("top-down") emissions for NOx over footprints that account for ∼25% of US oil and natural gas production. Across CONUS, we find that the bottom-up FOG inventory combined with other anthropogenic emissions is on average within ∼10% of top-down aircraft-derived NOx emissions. We also find good agreement in the trends of NOx from drilling- and production-phase activities, as inferred by satellites and in the bottom-up inventory. Leveraging tracer-tracer relationships derived from aircraft observations, methane (CH4) and non-methane volatile organic compound (NMVOC) emissions have been added to the inventory. Our total CONUS emission estimates for 2015 of oil and natural gas are 0.45 ± 0.14 Tg NOx/yr, 15.2 ± 3.0 Tg CH4/yr, and 5.7 ± 1.7 Tg NMVOC/yr. Compared to the US National Emissions Inventory and Greenhouse Gas Inventory, FOG NOx emissions are ∼40% lower, while inferred CH4 and NMVOC emissions are up to a factor of ∼2 higher. This suggests that NMVOC/NOx emissions from oil and gas basins are ∼3 times higher than current estimates and will likely affect how air quality models represent ozone formation downwind of oil and gas fields.

HONO Emissions from Western U.S. Wildfires Provide Dominant Radical Source in Fresh Wildfire Smoke.

Wildfires are an important source of nitrous acid (HONO), a photolabile radical precursor, yet in situ measurements and quantification of primary HONO emissions from open wildfires have been scarce. We present airborne observations of HONO within wildfire plumes sampled during the Western Wildfire Experiment for Cloud chemistry, Aerosol absorption and Nitrogen (WE-CAN) campaign. ΔHONO/ΔCO close to the fire locations ranged from 0.7 to 17 pptv ppbv-1 using a maximum enhancement method, with the median similar to previous observations of temperate forest fire plumes. Measured HONO to NOx enhancement ratios were generally factors of 2, or higher, at early plume ages than previous studies. Enhancement ratios scale with modified combustion efficiency and certain nitrogenous trace gases, which may be useful to estimate HONO release when HONO observations are lacking or plumes have photochemical exposures exceeding an hour as emitted HONO is rapidly photolyzed. We find that HONO photolysis is the dominant contributor to hydrogen oxide radicals (HOx = OH + HO2) in early stage (<3 h) wildfire plume evolution. These results highlight the role of HONO as a major component of reactive nitrogen emissions from wildfires and the main driver of initial photochemical oxidation.

Development of a Fuel-Based Oil and Gas Inventory of Nitrogen Oxides Emissions.

In this study, we develop an alternative Fuel-based Oil and Gas inventory (FOG) of nitrogen oxides (NO x) from oil and gas production using publicly available fuel use records and emission factors reported in the literature. FOG is compared with the Environmental Protection Agency's 2014 National Emissions Inventory (NEI) and with new top-down estimates of NO x emissions derived from aircraft and ground-based field measurement campaigns. Compared to our top-down estimates derived in four oil and gas basins (Uinta, UT, Haynesville, TX/LA, Marcellus, PA, and Fayetteville, AR), the NEI overestimates NO x by over a factor of 2 in three out of four basins, while FOG is generally consistent with atmospheric observations. Challenges in estimating oil and gas engine activity, rather than uncertainties in NO x emission factors, may explain gaps between the NEI and top-down emission estimates. Lastly, we find a consistent relationship between reactive odd nitrogen species (NO y) and ambient methane (CH4) across basins with different geological characteristics and in different stages of production. Future work could leverage this relationship as an additional constraint on CH4 emissions from oil and gas basins.

Modeling Ozone in the Eastern U.S. using a Fuel-Based Mobile Source Emissions Inventory.

Recent studies suggest overestimates in current U.S. emission inventories of nitrogen oxides (NO x = NO + NO2). Here, we expand a previously developed fuel-based inventory of motor-vehicle emissions (FIVE) to the continental U.S. for the year 2013, and evaluate our estimates of mobile source emissions with the U.S. Environmental Protection Agency's National Emissions Inventory (NEI) interpolated to 2013. We find that mobile source emissions of NO x and carbon monoxide (CO) in the NEI are higher than FIVE by 28% and 90%, respectively. Using a chemical transport model, we model mobile source emissions from FIVE, and find consistent levels of urban NO x and CO as measured during the Southeast Nexus (SENEX) Study in 2013. Lastly, we assess the sensitivity of ozone (O3) over the Eastern U.S. to uncertainties in mobile source NO x emissions and biogenic volatile organic compound (VOC) emissions. The ground-level O3 is sensitive to reductions in mobile source NO x emissions, most notably in the Southeastern U.S. and during O3 exceedance events, under the revised standard proposed in 2015 (>70 ppb, 8 h maximum). This suggests that decreasing mobile source NO x emissions could help in meeting more stringent O3 standards in the future.

Observational Constraints on the Oxidation of NOx in the Upper Troposphere.

NOx (NOx ≡ NO + NO2) regulates O3 and HOx (HOx ≡ OH + HO2) concentrations in the upper troposphere. In the laboratory, it is difficult to measure rates and branching ratios of the chemical reactions affecting NOx at the low temperatures and pressures characteristic of the upper troposphere, making direct measurements in the atmosphere especially useful. We report quasi-Lagrangian observations of the chemical evolution of an air parcel following a lightning event that results in high NOx concentrations. These quasi-Lagrangian measurements obtained during the Deep Convective Clouds and Chemistry experiment are used to characterize the daytime rates for conversion of NOx to different peroxy nitrates, the sum of alkyl and multifunctional nitrates, and HNO3. We infer the following production rate constants [in (cm(3)/molecule)/s] at 225 K and 230 hPa: 7.2(±5.7) × 10(-12) (CH3O2NO2), 5.1(±3.1) × 10(-13) (HO2NO2), 1.3(±0.8) × 10(-11) (PAN), 7.3(±3.4) × 10(-12) (PPN), and 6.2(±2.9) × 10(-12) (HNO3). The HNO3 and HO2NO2 rates are ∼ 30-50% lower than currently recommended whereas the other rates are consistent with current recommendations to within ±30%. The analysis indicates that HNO3 production from the HO2 and NO reaction (if any) must be accompanied by a slower rate for the reaction of OH with NO2, keeping the total combined rate for the two processes at the rate reported for HNO3 production above.

Air quality implications of the Deepwater Horizon oil spill.

During the Deepwater Horizon (DWH) oil spill, a wide range of gas and aerosol species were measured from an aircraft around, downwind, and away from the DWH site. Additional hydrocarbon measurements were made from ships in the vicinity. Aerosol particles of respirable sizes were on occasions a significant air quality issue for populated areas along the Gulf Coast. Yields of organic aerosol particles and emission factors for other atmospheric pollutants were derived for the sources from the spill, recovery, and cleanup efforts. Evaporation and subsequent secondary chemistry produced organic particulate matter with a mass yield of 8 ± 4% of the oil mixture reaching the water surface. Approximately 4% by mass of oil burned on the surface was emitted as soot particles. These yields can be used to estimate the effects on air quality for similar events as well as for this spill at other times without these data. Whereas emission of soot from burning surface oil was large during the episodic burns, the mass flux of secondary organic aerosol to the atmosphere was substantially larger overall. We use a regional air quality model to show that some observed enhancements in organic aerosol concentration along the Gulf Coast were likely due to the DWH spill. In the presence of evaporating hydrocarbons from the oil, NO(x) emissions from the recovery and cleanup operations produced ozone.

Observations of ozone, acyl peroxy nitrates, and their precursors during summer 2019 at Carlsbad Caverns National Park, New Mexico.

Carlsbad Caverns National Park (CAVE) is located in southeastern New Mexico and is adjacent to the Permian Basin, one of the most productive oil and natural gas (O&G) production regions in the United States. Since 2018, ozone (O3) at CAVE has frequently exceeded the 70 ppbv 8-hour National Ambient Air Quality Standard. We examine the influence of regional emissions on O3 formation using observations of O3, nitrogen oxides (NOx = NO + NO2), a suite of volatile organic compounds (VOCs), peroxyacetyl nitrate (PAN), and peroxypropionyl nitrate (PPN). Elevated O3 and its precursors are observed when the wind is from the southeast, the direction of the Permian Basin. We identify 13 days during the July 25 to September 5, 2019 study period when the maximum daily 8-hour average (MDA8) O3 exceeded 65 ppbv; MDA8 O3 exceeded 70 ppbv on 5 of these days. The results of a positive matrix factorization (PMF) analysis are used to identify and attribute source contributions of VOCs and NOx. On days when the winds are from the southeast, there are larger contributions from factors associated with primary O&G emissions; and, on high O3 days, there is more contribution from factors associated with secondary photochemical processing of O&G emissions. The observed ratio of VOCs to NOx is consistently high throughout the study period, consistent with NOx-limited O3 production. Finally, all high O3 days coincide with elevated acyl peroxy nitrate abundances with PPN to PAN ratios > 0.15 ppbv ppbv-1 indicating that anthropogenic VOC precursors, and often alkanes specifically, dominate the photochemistry.Implications: The results above strongly indicate NOx-sensitive photochemistry at Carlsbad Caverns National Park indicating that reductions in NOx emissions should drive reductions in O3. However, the NOx-sensitivity is largely driven by emissions of NOx into a VOC-rich environment, and a high PPN:PAN ratio and its relationship to O3 indicate substantial influence from alkanes in the regional photochemistry. Thus, simultaneous reductions in emissions of NOx and non-methane VOCs from the oil and gas sector should be considered for reducing O3 at Carlsbad Caverns National Park. Reductions in non-methane VOCs will have the added benefit of reducing formation of other secondary pollutants and air toxics.

Source characterization of volatile organic compounds at Carlsbad Caverns National Park.

Carlsbad Caverns National Park (CAVE), located in southeastern New Mexico, experiences elevated ground-level ozone (O3) exceeding the National Ambient Air Quality Standard (NAAQS) of 70 ppbv. It is situated adjacent to the Permian Basin, one of the largest oil and gas (O&G) producing regions in the US. In 2019, the Carlsbad Caverns Air Quality Study (CarCavAQS) was conducted to examine impacts of different sources on ozone precursors, including nitrogen oxides (NOx) and volatile organic compounds (VOCs). Here, we use positive matrix factorization (PMF) analysis of speciated VOCs to characterize VOC sources at CAVE during the study. Seven factors were identified. Three factors composed largely of alkanes and aromatics with different lifetimes were attributed to O&G development and production activities. VOCs in these factors were typical of those emitted by O&G operations. Associated residence time analyses (RTA) indicated their contributions increased in the park during periods of transport from the Permian Basin. These O&G factors were the largest contributor to VOC reactivity with hydroxyl radicals (62%). Two PMF factors were rich in photochemically generated secondary VOCs; one factor contained species with shorter atmospheric lifetimes and one with species with longer lifetimes. RTA of the secondary factors suggested impacts of O&G emissions from regions farther upwind, such as Eagle Ford Shale and Barnett Shale formations. The last two factors were attributed to alkenes likely emitted from vehicles or other combustion sources in the Permian Basin and regional background VOCs, respectively.Implications: Carlsbad Caverns National Park experiences ground-level ozone exceeding the National Ambient Air Quality Standard. Volatile organic compounds are critical precursors to ozone formation. Measurements in the Park identify oil and gas production and development activities as the major contributors to volatile organic compounds. Emissions from the adjacent Permian Basin contributed to increases in primary species that enhanced local ozone formation. Observations of photochemically generated compounds indicate that ozone was also transported from shale formations and basins farther upwind. Therefore, emission reductions of volatile organic compounds from oil and gas activities are important for mitigating elevated O3 in the region.

Evaluating the Impact of Chemical Complexity and Horizontal Resolution on Tropospheric Ozone Over the Conterminous US With a Global Variable Resolution Chemistry Model.

A new configuration of the Community Earth System Model (CESM)/Community Atmosphere Model with full chemistry (CAM-chem) supporting the capability of horizontal mesh refinement through the use of the spectral element (SE) dynamical core is developed and called CESM/CAM-chem-SE. Horizontal mesh refinement in CESM/CAM-chem-SE is unique and novel in that pollutants such as ozone are accurately represented at human exposure relevant scales while also directly including global feedbacks. CESM/CAM-chem-SE with mesh refinement down to ∼14 km over the conterminous US (CONUS) is the beginning of the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICAv0). Here, MUSICAv0 is evaluated and used to better understand how horizontal resolution and chemical complexity impact ozone and ozone precursors over CONUS as compared to measurements from five aircraft campaigns, which occurred in 2013. This field campaign analysis demonstrates the importance of using finer horizontal resolution to accurately simulate ozone precursors such as nitrogen oxides and carbon monoxide. In general, the impact of using more complex chemistry on ozone and other oxidation products is more pronounced when using finer horizontal resolution where a larger number of chemical regimes are resolved. Large model biases for ozone near the surface remain in the Southeast US as compared to the aircraft observations even with updated chemistry and finer horizontal resolution. This suggests a need for adding the capability of replacing sections of global emission inventories with regional inventories, increasing the vertical resolution in the planetary boundary layer, and reducing model biases in meteorological variables such as temperature and clouds.

Observational constraints on glyoxal production from isoprene oxidation and its contribution to organic aerosol over the Southeast United States.

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.

Airborne measurements of organosulfates over the continental U.S.

Organosulfates are important secondary organic aerosol (SOA) components and good tracers for aerosol heterogeneous reactions. However, the knowledge of their spatial distribution, formation conditions, and environmental impact is limited. In this study, we report two organosulfates, an isoprene-derived isoprene epoxydiols (IEPOX) (2,3-epoxy-2-methyl-1,4-butanediol) sulfate and a glycolic acid (GA) sulfate, measured using the NOAA Particle Analysis Laser Mass Spectrometer (PALMS) on board the NASA DC8 aircraft over the continental U.S. during the Deep Convective Clouds and Chemistry Experiment (DC3) and the Studies of Emissions and Atmospheric Composition, Clouds, and Climate Coupling by Regional Surveys (SEAC4RS). During these campaigns, IEPOX sulfate was estimated to account for 1.4% of submicron aerosol mass (or 2.2% of organic aerosol mass) on average near the ground in the southeast U.S., with lower concentrations in the western U.S. (0.2-0.4%) and at high altitudes (<0.2%). Compared to IEPOX sulfate, GA sulfate was more uniformly distributed, accounting for about 0.5% aerosol mass on average, and may be more abundant globally. A number of other organosulfates were detected; none were as abundant as these two. Ambient measurements confirmed that IEPOX sulfate is formed from isoprene oxidation and is a tracer for isoprene SOA formation. The organic precursors of GA sulfate may include glycolic acid and likely have both biogenic and anthropogenic sources. Higher aerosol acidity as measured by PALMS and relative humidity tend to promote IEPOX sulfate formation, and aerosol acidity largely drives in situ GA sulfate formation at high altitudes. This study suggests that the formation of aerosol organosulfates depends not only on the appropriate organic precursors but also on emissions of anthropogenic sulfur dioxide (SO2), which contributes to aerosol acidity. KEY POINTS: IEPOX sulfate is an isoprene SOA tracer at acidic and low NO conditions Glycolic acid sulfate may be more abundant than IEPOX sulfate globally SO2 impacts IEPOX sulfate by increasing aerosol acidity and water uptake.

Infrared overtone spectroscopy and unimolecular decay dynamics of peroxynitrous acid.

Peroxynitrous acid (HOONO) is generated in a pulsed supersonic expansion through recombination of photolytically generated OH and NO(2) radicals. A rotationally resolved infrared action spectrum of HOONO is obtained in the OH overtone region at 6971.351(4) cm(-1) (origin), providing definitive spectroscopic identification of the trans-perp (tp) conformer of HOONO. Analysis of the rotational band structure yields rotational constants for the near prolate asymmetric top, the ratio of the a-type to c-type components of the transition dipole moment for the hybrid band, and a homogeneous linewidth arising from intramolecular vibrational energy redistribution and/or dissociation. The quantum state distribution of the OH (nu=0,J(OH)) products from dissociation is well characterized by a microcanonical statistical distribution constrained only by the energy available to products, 1304+/-38 cm(-1). This yields a 5667+/-38 cm(-1) [16.2(1) kcal mol(-1)] binding energy for tp-HOONO. An equivalent available energy and corresponding binding energy are obtained from the highest observed OH product state. Complementary high level ab initio calculations are carried out in conjunction with second-order vibrational perturbation theory to predict the spectroscopic observables associated with the OH overtone transition of tp-HOONO including its vibrational frequency, rotational constants, and transition dipole moment. The same approach is used to compute frequencies and intensities of multiple quantum transitions that aid in the assignment of weaker features observed in the OH overtone region, in particular, a combination band of tp-HOONO involving the HOON torsional mode.

Fluorescence-dip infrared spectroscopy and predissociation dynamics of OH A 2Sigma+ (v = 4) radicals.

Fluorescence-dip infrared spectroscopy, an UV-IR double-resonance technique, is employed to characterize the line positions, linewidths, and corresponding lifetimes of highly predissociative rovibrational levels of the excited A (2)Sigma(+) electronic state of the OH radical. Various lines of the 4 <--2 overtone transition in the excited A (2)Sigma(+) state are observed, from which the rotational, centrifugal distortion, and spin-rotation constants for the A (2)Sigma(+) (v = 4) state are determined, along with the vibrational frequency for the overtone transition. Homogeneous linewidths of 0.23-0.31 cm(-1) full width at half maximum are extracted from the line profiles, demonstrating that the N = 0 to 7 rotational levels of the OH A (2)Sigma(+) (v = 4) state undergo rapid predissociation with lifetimes of < or =23 ps. The experimental linewidths are in near quantitative agreement with first-principles theoretical predictions.