Methane and nonmethane volatile organic compounds (VOCs) were monitored near Boulder in the Northern Colorado Front Range to investigate their spatial distribution and sources as a part of the Front Range Air Pollution and Photochemistry Experiment (FRAPPE) and the Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) campaign, in summer 2014. A particular emphasis was the study of the contribution of emissions from oil and natural gas (O&NG) operations on the regional air quality. One network extended along an elevation gradient from the City of Boulder (elevation ≈1,600 m) to the University of Colorado Mountain Research Station (≈2900 m) on the eastern slopes of the Rocky Mountains. Light alkane petroleum hydrocarbons had the highest mole fraction of the VOCs that could be analyzed with the applied techniques. The longer lived VOCs ethane and propane decreased with increasing elevation, suggesting that Boulder and the surrounding plains were a source of these anthropogenic compounds. VOC diurnal time series showed a few events with elevated mole fractions at the mountain sites, which were likely the result of the upslope transport of plumes with elevated VOCs from the plains. Within the other site network, which extended into suburban East Boulder County (EBC), VOCs were monitored at 5 sites increasingly close to O&NG development in the Denver Julesburg Basin. Mean mole fractions and variability of primarily O&NG-associated VOCs (ethane, propane, butane isomers) increased by a factor of 2.4–5.2 with closer proximity to the O&NG producing region. Median mole fractions of C2–C5 n-alkanes and of imuch-butane at the EBC sites were higher than those previously reported from 28 larger urban areas in the United States. Among the VOCs that could be quantified with the gas chromatography methods, VOCs most clearly associated to O&NG-related emissions (C2–C5 alkanes) accounted for 52%–79% of the VOC hydroxyl radical reactivity (OHR). The horizontal gradient in OHR of the considered VOCs, with ≈3 times higher values at the furthest eastern sites, points toward higher chemical reactivity and ozone production potential from these ozone precursors in the eastern area of the county than within the City of Boulder.

The city of Boulder, CO, is located in the Northern Colorado Front Range (NCFR) east of the Rocky Mountains. Local wind patterns expose Boulder to air plumes from many directions. Relatively clean air reaches Boulder from over the mountains to the west. This downslope flow occurs mostly at night, when cooler air from the mountains descends and flows to the east (Johnson and Toth, 1982). Upslope flow that causes air flow from east to west is especially dominant during daytime in summer. This flow regime can transport polluted air from the plains up the mountain slope from midmorning to late afternoon. This is reflected in diurnal cycles of pollutants in data from elevation sites exhibiting distinctly increasing concentrations during midday to late afternoon. The onset of the increase is becoming later in the day with increasing elevation, reflecting the east-to-west and upward (in elevation) progression of the polluted air masses (Brodin et al., 2010; Benedict et al., 2018; Darby et al., 2021). There have been a few reported cases that demonstrated that this flow can make it all across the high mountain tops of the Continental Divide and descend into the Frazier Valley to the west of the Divide (Pfister et al., 2017b; Flocke et al., 2020). Midday to afternoon concentration increases resulting from upslope flow have been observed both for primary emissions (i.e., volatile organic compounds (VOCs; Benedict et al., 2019), oxidized nitrogen compounds (Benedict et al., 2018; Benedict et al., 2019), and secondary pollutants (i.e., ozone (Brodin et al., 2010; Darby et al., 2021), peroxy acetyl nitrate and alkyl nitrates (Benedict et al., 2019), and the OH radical (Mount and Harder, 1995).

In 2016, the annual average daily vehicle count from the 6 most heavily trafficked access points to Boulder was ≈160,000 vehicles (Colorado-DOT, 2016), making vehicle traffic into and within the city a major contributor to air emissions. Being ≈50 km northwest of Denver, Boulder is subjected to transport of urban and industrial emissions from the Denver Metro Area when winds are from the southeast. Weld County is ≈50 km northeast of Boulder and is home to most of the oil and natural gas (O&NG) drilling and extraction in the Denver Julesburg Basin (DJB). Besides O&NG emissions, transportation sector emissions and agricultural operations are another considerable contributing source from this sector. There has been a remarkable growth of active wells and O&NG production in Weld County. From 2010 to 2014, annual oil production approximately quadrupled (21–81 million barrels), and annual gas production nearly doubled (219–392 billion cubic feet; COGCC, 2014, 2016).

Enhanced levels of VOCs have been observed in ground-level air near O&NG wells, including the carcinogens benzene and toluene (Verma et al., 2000; Huff, 2007; Macey et al., 2014; Sovacool, 2014; McKenzie et al., 2018). Exposure to emissions related to O&NG development has demonstrated deleterious effects on human health; a relation has been observed between maternal residence within 10 miles of O&NG wells and congenital heart defects in their children (McKenzie et al., 2014). Other studies have demonstrated that health risks such as birth defects, low birth weight, preterm birth, and cancer are increasingly likely when exposed to air near O&NG sites (McKenzie et al., 2014; McKenzie et al., 2018). Further, VOCs contribute to the photochemical formation of ozone and secondary aerosol (Chameides et al., 1992; Bowman and Seinfeld, 1994), which pose additional health concerns. Fann et al. (2018) forecasted that ozone and particulate matter pollution resulting from O&NG emissions will be responsible for 70 premature deaths in Colorado by 2025.

The NCFR was classified as a “moderate” nonattainment area (NAA) for the 8-h 75 ppb ozone National Ambient Air Quality (NAAQS) standard in 2012 (Bien and Helmig, 2018). After continuing exceedances of the NAAQS in subsequent years, the NAA classification was bumped up to a “serious” status in 2019 (EPA, 2019). This places the area under an obligation to reduce ozone precursor compounds. Identification and characterization of ozone precursors, including sources of nitrogen oxides and biogenic and anthropogenic VOC, and their contribution to the regional ozone production is an imperative in that process.

Previous research reported that O&G operations in the DJB are an important and potentially dominant source of regional VOC emissions, resulting in highly elevated ambient air concentrations of VOCs (Pétron et al., 2012; Brown et al., 2013; Gilman et al., 2013; Swarthout et al., 2013; Pétron et al., 2014; Thompson et al., 2014) and deteriorated air quality (Helmig, 2020). While numerous studies have pointed toward a probable important contribution of O&NG emissions to the regional ozone production and NAAQS exceedances, these mostly single-point and short-term observations have been insufficient for quantitatively assessing the O&NG contribution to local ozone production in the NCFR under the diverse spatial distribution of sources categories, land cover, topography, and meteorological conditions.

These questions motivated the Front Range Air Pollution and Photochemistry Experiment (FRAPPE) and the Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) study that took place in parallel in the wider NCFR during the summer of 2014 (Flocke et al., 2020). The research reported here contributed to these campaigns. The emphasis of our study was the investigation of spatial surface gradients of VOCs within and in the area surrounding Boulder, in particular of O&NG-related VOCs and their relation to O&NG well locations.

VOC and methane data were collected from 2 networks of sites. One network followed an elevation gradient from Boulder (≈1,600 m above sea level [asl]) to the University of Colorado Research Station (≈2,900 m asl) near Ward, CO. The other network included 5 sites in the eastern portion of Boulder County (EBC).

2.1. Elevation gradient measurements

Locations, elevations, site abbreviations, and sampling methods used at sites in the elevation gradient network are available in Table 1; locations are also marked on the map in Figure 1. Stainless steel (Summa) canisters were deployed to each site for whole air sample collection on 8 days in July and August 2014. A programmable flask package (PFPs) was deployed at the Lost Angels (LA) site on 4 days in July and August 2014. Continuous ozone observations from these sites during June–August, 2014, are reported by Oltmans et al. (2019).

Table 1.

Monitoring sites with site name abbreviation, location, elevation, and sampling methods. DOI: https://doi.org/10.1525/elementa.2019.00036.t1

Site NameSite Abbr.Lat.Lon.Elevation (m)Sampling Methods EmployedDistance to Center of O&NG Region (km)Population Density (km–2)
Soddie Lab (MRS) SO 40° 2′ 52″ N 105° 34′ 15″ W 3,345 3-h canisters — — 
Coughlin CO 40° 0′ 13″ N 105° 28′ 43″ W 2,523 3-h canisters — — 
Lost Angels Road LA 40° 1′ 8″ N 105° 24′ 18″ W 2,392 PFP, 3-h canisters — — 
Sugarloaf Road SL 40° 1′ 1″ N 105° 21′ 28″ W 1,978 PFP, 3-h canisters — — 
INSTAAR, CU-Boulder BO 40° 0′ 48” 105° 15′ 11″ W 1,607 in situ GC — — 
Boulder County Public Health BCPH 40° 2′ 14″ N 105° 16′ 53″ W 1,569 3-day canisters 48 1,568 
Lafayette Fire Department LFD 40° 0′ 6″ N 105° 5′ 40″ W 1,600 3-day canisters 44 1,095 
Dawson School DAW 40° 3′ 50″ N 105° 6′ 37″ W 1,562 3-day canisters 39 77 
St. Luke’s Church SLC 40° 1′ 56″ N 105° 3′ 24″ W 1,550 3-day canisters 38 427 
Stephen Day Park SDP 40° 10′ 57” 105° 3′ 49″ W 1,522 3-day canisters 30 1,210 
Site NameSite Abbr.Lat.Lon.Elevation (m)Sampling Methods EmployedDistance to Center of O&NG Region (km)Population Density (km–2)
Soddie Lab (MRS) SO 40° 2′ 52″ N 105° 34′ 15″ W 3,345 3-h canisters — — 
Coughlin CO 40° 0′ 13″ N 105° 28′ 43″ W 2,523 3-h canisters — — 
Lost Angels Road LA 40° 1′ 8″ N 105° 24′ 18″ W 2,392 PFP, 3-h canisters — — 
Sugarloaf Road SL 40° 1′ 1″ N 105° 21′ 28″ W 1,978 PFP, 3-h canisters — — 
INSTAAR, CU-Boulder BO 40° 0′ 48” 105° 15′ 11″ W 1,607 in situ GC — — 
Boulder County Public Health BCPH 40° 2′ 14″ N 105° 16′ 53″ W 1,569 3-day canisters 48 1,568 
Lafayette Fire Department LFD 40° 0′ 6″ N 105° 5′ 40″ W 1,600 3-day canisters 44 1,095 
Dawson School DAW 40° 3′ 50″ N 105° 6′ 37″ W 1,562 3-day canisters 39 77 
St. Luke’s Church SLC 40° 1′ 56″ N 105° 3′ 24″ W 1,550 3-day canisters 38 427 
Stephen Day Park SDP 40° 10′ 57” 105° 3′ 49″ W 1,522 3-day canisters 30 1,210 

MRS stands for University of Colorado’s Mountain Research Station. Distances to the O&NG producing region are included for East Boulder County (EBC) sites. These were determined using the distance function in Google Earth from each of the sites to the center of Weld County. Population density data were found at Bureau (2010) and Population (2016). EBC sites are listed in order of decreasing distance to the center of Weld County. GC = gas chromatography.

Figure 1.

Map displaying sites where VOC samples and methane measurements were collected. Blue stars indicate the elevation gradient sites, red stars the east Boulder County (EBC) sites, and purple dots indicate the locations of O&NG wells, respectively. The map with O&G well locations was obtained from the Colorado Oil and Gas Commission website (COGCC, 2016). DOI: https://doi.org/10.1525/elementa.2019.00036.f1

Figure 1.

Map displaying sites where VOC samples and methane measurements were collected. Blue stars indicate the elevation gradient sites, red stars the east Boulder County (EBC) sites, and purple dots indicate the locations of O&NG wells, respectively. The map with O&G well locations was obtained from the Colorado Oil and Gas Commission website (COGCC, 2016). DOI: https://doi.org/10.1525/elementa.2019.00036.f1

Close modal

The 6-L Summa canisters (part #27405, Restek, Bellefonte, PA, USA) were sampled at the elevation gradient sites for 3-h periods on 8 days when FRAPPE/DISOVER-AQ aircraft missions took place. The canisters were conditioned prior to deployment by 3 cycles of filling with humidified zero air and evacuation to a final pressure of ≈1–2 in Hg. Canisters were placed within temperature-controlled buildings, and air was sampled from the outside via 0.64 cm o.d. sampling lines with an in-line ozone scrubber (glass fiber filter soaked in sodium thiosulfate solution [preparation details available in Helmig, 1997]) from inlets that were mounted 3–5 m above the surface. The sampling line was purged continuously by a pump at 1.2 L min–1. Collection of air into the canister was through a tee in the sampling line and controlled by a programmable timer and Veriflo passive flow devices (part #24232, Restek) that maintained a constant flow of ≈25 mL min–1 into the canister. The sampling was conducted over the 3-h period from noon to 3 PM. The final pressure in the canisters was ≈19 in Hg. Canister samples were analyzed via gas chromatography (GC) within 1–3 days of sampling with the column flow split for compound identification by mass spectrometry (MS) and quantification by flame ionization detection (FID; more detail on the analysis method is provided in Section 2.4).

The PFP (Atmospheric Observing Systems, INC., Boulder, CO, USA) was deployed at the Lost Angels Road fire station site (LA) to investigate diurnal cycles of VOCs. The PFP contains twelve ≈0.8 L whole air sampling glass flasks, a set of shutoff valves, and an external compressor, to allow purging and then pressurizing the flasks sequentially with ambient air. Air was collected through an inlet line equipped with a sodium thiosulfate-coated inlet filter mounted at 2 m height outside of the building. Air samples were collected at 2-h intervals during 5 deployments (July 23, July 31, August 1, August 3–4, August 10–11), yielding a total of 42 samples. These samples were analyzed with the same analytical system as the canisters. The analysis procedure is detailed in Section 2.4.

2.2. East Boulder measurements

These samples were collected to investigate the influence of the nearby O&NG drilling region on the surface horizontal gradients of atmospheric methane and VOCs in the plains. The 5 sampling sites locations, elevations, sampling methods, and distances from the center of Weld County are listed in Table 1. Figure 1 shows the locations of these sites in relation to locations of O&NG wells.

The Boulder County Public Health (BCPH) building was chosen as the “urban reference” site for contrasting the sites in EBC. The BCPH site is located in the City of Boulder and is subjected to typical urban emission profiles emanating from mobile and industrial sources and is the furthest of the 5 sites from Boulder County’s border with Weld County (≈25 km), which marks the southwest corner of the region of concentrated O&NG operations. The population densities displayed in Table 1 were calculated based on the population and area of the nearest city, town, or incorporated settlement to the monitoring site (Bureau, 2010; Population, 2016).

Samples at the EBC sites were also collected in 6 L-Summa stainless steel canisters (Restek #27420, Bellefonte, PA). Canisters were deployed on 12 occasions for 3-day periods from mid-May to late August. Because these samples were collected over multi-day timescales, they are more representative of average conditions at a given site than shorter sampling techniques such as grab or aircraft sampling. The sampling schedule is provided in Figure S1. An example of the sampling setup is presented in Figure S2. In total, 36 days of canister sampling were conducted at each site. Consequently, the data obtained from this sampling represent the average of a total of 36 sampling days spread over the summer of 2014. To minimize temperature variations during the sampling, canisters were placed in ≈45 L insulated ice chest boxes with ≈6 L of water contained in plastic bottles. The sampling line, made of 0.16 cm o.d. silcosteel tubing (part #29191, Restek, Bellefonte, PA, USA), was fed through a hole in the box, and the inlet was fit with an ozone scrubber (same as described above). Effective sampling heights were ≈3 m above ground level. Canister sampling was controlled using Veriflo passive flow devices (part #24232, Restek) that maintained a constant flow of ≈1 mL min–1 into the canister over the 3-day sampling window. At the end of the sampling period, the final pressures in the canisters were ≈19 in Hg. Canisters were capped and transported to the lab for analysis (details in Section 2.4). Air samples were analyzed for methane and C2–C8 VOCs (including benzene, toluene, ethylbenzene, and xylenes [BTEX compounds]) using GC-FID/MS (Section 2.4).

2.3. Continuous GC data

Air was sampled from an inlet outside a second-story window at the former INSTAAR building (Rose Litman Laboratory 1) on the University of Colorado’s east campus (BO site, Figure 1). The sample preconcentration and injection system was similar in its configuration as the instrument described in Tanner et al. (2006) using an SRI 310 GC/FID (SRI Instruments #0310-6003, Torrance, CA) with a 0.32 mm i.d., 30 m length GasPro PLOT column (Agilent Technologies #113-4332, Santa Clara, CA). Ethane was not quantified because the GC ethane peak was not adequately separated. VOCs were analyzed on FRAPPE/DISCOVER-AQ flight days continuously at 2-h time resolution.

2.4. Canister and PFP sample analysis

First, methane was analyzed by drawing sample from the canisters (or PFPs) through a 3 mL sample loop at a flow rate of 50 mL min–1 for 3 min. The loop content was then injected onto a Porapak Q packed column (0.31 cm, 0.21 cm I.D., 3.6 m long; Agilent Technologies #113-4362, Santa Clara, CA). GC separation was isothermal (at 50°C) with FID. The FID was calibrated using a whole air standard made of compressed air (breathing air); its methane mole fraction of 1,898.73 ppb was determined by reference analysis against the NOAA Global Monitoring Division methane scale. The methane standard was run before and after each set of samples (no more than 7 samples analyzed at a time).

Following methane analysis, VOC sample volumes of 500 mL were first scrubbed of water by passing through a cryogenic trap (–30°C). VOCs were then preconcentrated on a Peltier-cooled micro-adsorbent trap (–30°C; Pollmann et al., 2006) and then injected into a GC-FID/MS equipped with a 60 m, 0.32 mm i.d. GAS-PRO PLOT column (Agilent Technologies #80468-800, Santa Clara, CA), with the column flow split between a second FID (i.e., in addition to the FID used for methane analysis) and an MS. The GC oven temperature began at 50°C for 2 min, was then ramped at a rate of 8°C min–1, and held at 225°C for 15 min.

Blanks and calibration standards were analyzed prior to and immediately after each set of 7 canisters, and canister analysis results were blank-subtracted. Blank sample air was prepared by a custom zero-air generator that used a heated PtO catalyst to remove hydrocarbons. Blank experiments were conducted by direct injection of zero air into the GC-MS system and by retrieving zero air from canisters after filling them with the blank air. Calibrations were performed using a series of gravimetrically prepared VOC standards (Table S1). Our laboratory VOC calibration scale has been cross-referenced against other national and international laboratories (Rhoderick et al., 2014), including 3 audits by the World Calibration Centre for VOCs (WCC-VOC, 2016). Measurement uncertainties were between 2% and 5% for the C2–C6 alkanes, 9% for acetylene, and 25% for benzene (due to benzene background issues with Summa canisters and a higher benzene background from the preconcentration system, as well as a slight overlap with a neighboring peak in the chromatogram, the benzene uncertainty is greater than for other VOCs). Detection limits were mostly in the range of 0.01–0.05 ppb (Table S2). Ethene and higher alkene results were omitted from the data because of uncertain accuracy of the data, as was indicated by inconsistencies between results from Summa canisters and in situ monitoring at sub-ppb mole fractions for ethene. An assessment of the relative significance of not quantified alkenes to the total VOCs spectrum can be derived from currently ongoing in situ VOCs monitoring at the Boulder Reservoir (Boulder_AIR, 2021). It needs to be emphasized that these data are from different, albeit close (5–25 km) locations and not contemporaneous, with the Boulder Reservoir monitoring having commenced in 2017. Nonetheless, most likely, these data comparison will allow for a reasonable order of magnitude evaluation. In the Boulder Reservoir data, ethene was determined to constitute 1.6%, propene 0.3%, and other minor VOCs a total of 3.1% of the total VOC carbon (excluding isoprene; Table S4). We therefore estimate that the VOCs species considered here amount to ≈95% of the total VOCs that are measureable with the applied sampling and GC method. Isoprene results were omitted from this estimation as ambient air isoprene mole fractions have been shown to be highly variable in the region, likely from the influence of the nearby vegetation type and density at the particular site. In addition, there are oxygenated VOCs present in the atmosphere, for example, methanol, formaldehyde, ethanol, acetaldehyde, and so on, that are not captured by the applied analytical sampling method and were therefore not included in the data analyses.

PFP blanks were determined by filling the PFP flasks with VOC-free air (zero air), letting it sit for several days, and then conducting the same analysis as for the ambient samples. Blanks for the VOCs considered here were nondetects or below 5% of typical ambient mole fractions and were therefore not considered in the quantification calculation of ambient samples. Data are visualized using box whisker plots. Data comparisons for days when PFPs and canisters were sampled along the elevation gradient on the same days are presented in Table S3.

3.1. Elevation gradient data

Figure 2 shows VOC mole fractions from the canister measurements along the elevation gradient, with sites arranged by decreasing elevation from left to right. Ethane mole fractions and the variability in the data increase as elevation decreases. For all sites but Soddie (SO), the ethane mean mole fraction is at least 2 ppb higher than the median, indicating that means are biased upward by selected occurrences of high mole fractions. Propane shows a similar behavior, except that mole fractions are highest and most variable at the Lost Angels (LA) site. Mountain homes in the area are primarily fueled by propane gas with storage tanks that are located on the premises. This accounts to a high number of potential point sources within the area that may explain why observations at LA deviate from the elevation gradient seen for ethane.

Figure 2.

VOC mole fractions from canister measurements along the elevation gradient. Sites on the x-axis are in order of decreasing elevation. Boxes show the range of 25th–75th percentile data. The middle line represents the median. The diamond represents the mean. Whiskers show 10th and 90th percentile ranges. DOI: https://doi.org/10.1525/elementa.2019.00036.f2

Figure 2.

VOC mole fractions from canister measurements along the elevation gradient. Sites on the x-axis are in order of decreasing elevation. Boxes show the range of 25th–75th percentile data. The middle line represents the median. The diamond represents the mean. Whiskers show 10th and 90th percentile ranges. DOI: https://doi.org/10.1525/elementa.2019.00036.f2

Close modal

For the Coughlin (CO) site, 2 samples showed anomalously high (more than 10 times higher than the median of other data) mole fractions for VOCs > C3. This mostly concerned gasoline-type compounds. Those samples also showed an i/n-pentane ratio that was larger than 1.5, which is characteristic for gasoline and urban emissions, while data from the other, further eastern sites, had lower i/n-pentane ratios. We suspect that these sample values resulted from a nearby petroleum source, possibly a gasoline leak from a nearby residence. It is possible that the mole fractions observed at the CO site are exacerbated by the site’s location in a valley defined by Bald Mountain, Klondike Mountain, and Hurricane Hill. Those data were therefore excluded from further consideration.

For the other samples, mole fractions for >C3 compounds also generally increase as elevation decreases, but higher median mole fractions are observed at LA than at the lower elevation Sugarloaf (SL) site for i-butane, the pentane isomers, hexane, benzene, and toluene. It is noteworthy that overall, the concentration gradient is most evident for the natural gas VOCs (i.e., ethane, propane, butane isomers) and becomes less pronounced with increasing molecular weight of the hydrocarbon. Scatterplots that show individual canister measurements as a function of elevation are available in Figure S3.

3.2. East Boulder County observations

3.2.1. Methane

Figure 3 shows box whisker plots of the statistical distribution of methane values at sites in the EBC network. The regional background level of methane (estimated from the 10th percentile value from NOAA monitoring at the Mountain Research Station Niwot Ridge site in the summer of 2014) is shown as a red horizontal line in the figure. Methane was lowest at the BCPH site and increased toward the east of the county; the overall highest levels were observed at Stephen Day Park (SDP) in Longmont. The methane enhancement (values above background) at SDP was approximately twice that seen in Boulder. There was a noticeable increase in variability with increasing absolute levels, indicated by the wider spread of the data at the more easterly EBC monitoring sites. These 2 features (higher absolute levels as well as higher variability) are a clear indication of higher abundance and stronger strength of methane sources in the eastern part of the county.

Figure 3.

Spatial gradient of methane in 3-day integrated canister samples collected at the 5 study sites. Each box was derived from 9 to 11 individual sample collections, with each individual sample being the mean for a 3-day sampling interval. The red line shows the approximate background level of methane (1.87 ppm) during summer 2014. Sites on the x-axis are ordered west to east, that is, by decreasing distance to Weld County. DOI: https://doi.org/10.1525/elementa.2019.00036.f3

Figure 3.

Spatial gradient of methane in 3-day integrated canister samples collected at the 5 study sites. Each box was derived from 9 to 11 individual sample collections, with each individual sample being the mean for a 3-day sampling interval. The red line shows the approximate background level of methane (1.87 ppm) during summer 2014. Sites on the x-axis are ordered west to east, that is, by decreasing distance to Weld County. DOI: https://doi.org/10.1525/elementa.2019.00036.f3

Close modal

Methane is emitted from a variety of sources. Correlations of methane with characteristic VOCs can provide insight into the contributing local sources. Pearson correlation coefficients (PCC) and the statistical significance (P-values) of orthogonal linear regressions of mole fractions of methane and VOCs in all data combined are reported in Table 2. Figure 4 visualizes these relationships. Data points are color coded by site; PCC are available in the title of each plot.

Table 2.

Correlation results of orthogonal linear regression (OR) analyses of volatile organic compound (VOC) with methane at EBC sites. DOI: https://doi.org/10.1525/elementa.2019.00036.t2

CompoundSlope of OR (Compound [ppb]/Methane [ppb])PCCP-Value
Ethane 9.12 × 10–2 0.75 0.000 
Propane 4.57 × 10–2 0.54 0.000 
i-Butane 1.22 × 10–2 0.41 0.006 
n-Butane 1.22 × 10–2 0.46 0.002 
i-Pentane 5.68 × 10–3 0.37 0.014 
n-Pentane 5.75 × 10–3 0.37 0.014 
CompoundSlope of OR (Compound [ppb]/Methane [ppb])PCCP-Value
Ethane 9.12 × 10–2 0.75 0.000 
Propane 4.57 × 10–2 0.54 0.000 
i-Butane 1.22 × 10–2 0.41 0.006 
n-Butane 1.22 × 10–2 0.46 0.002 
i-Pentane 5.68 × 10–3 0.37 0.014 
n-Pentane 5.75 × 10–3 0.37 0.014 

Slopes of the orthogonal linear regression (OR) fit, Pearson correlation coefficients (PCCs), and P-values are tabulated for several VOCs. P-values < 0.001 are reported as 0.000.

Figure 4.

Comparison of VOC and methane mole fractions at EBC sites. The data were fit with an orthogonal regression (OR). The Pearson Correlation Coefficients (PCCs) are reported in the title of selected graphs and in Table 2; slopes of the fits are also included in Table 2. Colored markers indicate the site at which the measurement was taken. Note that y-axis scales vary between plots. OR and PCC are not shown for compounds with PCCs fits <0.4. DOI: https://doi.org/10.1525/elementa.2019.00036.f4

Figure 4.

Comparison of VOC and methane mole fractions at EBC sites. The data were fit with an orthogonal regression (OR). The Pearson Correlation Coefficients (PCCs) are reported in the title of selected graphs and in Table 2; slopes of the fits are also included in Table 2. Colored markers indicate the site at which the measurement was taken. Note that y-axis scales vary between plots. OR and PCC are not shown for compounds with PCCs fits <0.4. DOI: https://doi.org/10.1525/elementa.2019.00036.f4

Close modal

The C2-C5 alkanes have statistically significant correlations with methane at >95% confidence (Table 2), with ethane and propane having the highest correlations (PCC = 0.75 and 0.54, respectively). The butane and pentane isomers have correlations with methane that range from PCC = 0.37 to 0.46. Poorer correlations were found for other compounds. Controlled and fugitive emissions, and incomplete flaring of natural gas is the globally dominant source of atmospheric ethane (Tzompa-Sosa et al., 2017). Secondary sources are biofuels and biomass burning. These 2 latter source categories were of negligible importance for the NCFR during the summer of 2014; therefore, it is highly likely that emissions from natural gas extraction industry are the dominant regional source of ethane. The tight correlation between methane and ethane in our data indicates that a major portion of methane and the light n-alkanes (especially ethane and propane) resulted from common sources. This has been pointed out previously: Pétron et al. (2014) estimated from their 2012 experiment that 75%; Peischl et al. (2018) estimated that during a campaign in 2014, 75 ± 37%, and Kille et al. (2019) estimated that in 2015, 63 ± 17% of observed methane emissions in the DJB were attributable to O&NG sources. Our mean ethane/methane ratio (EMR) result for the linear regression analysis (9.12% molar) for all 4 site data is ≈50% lower than mean slope values reported by Peischl et al. (2018) and Kille et al. (2019; that span a range of 13%–16%), which can possibly be attributed to a relatively lower O&NG influence in our sample population compared to these 2 comparison studies (i.e., methane contributions from non-O&NG methane sources will pull the EMR to lower values) that relied more heavily on sampling within the central region of the DJB. Similarly, our propane/methane molar ratio of 4.6% is smaller than a ratio of ≈15% that was found from mobile lab and stationary surface sampling in the central DJB by Pétron et al. (2014). The fact that VOC/methane correlations become weaker for the heavier VOCs (i.e., toluene and ethylbenzene) is likely due to 2 reasons: (1) The atmospheric lifetimes of VOCs from reaction with OH become shorter with increasing molecular size. Therefore, the fraction of the VOC that was quantified in each sample in relation to its fraction in the emissions will become smaller with increasing molecule size, and the variability in its oxidation rate becomes an increasing factor in determining absolute levels and variability in mole fraction. This is expected to weaken the correlation with the much longer-lived methane. (2) The heavier VOCs are emitted by more diverse sources, with a relatively smaller fraction of these VOC potentially being co-emitted with methane.

3.2.2. VOCs

The statistical distribution of VOC results from the sites in EBC (Figure 5) shows increasing mole fractions of the C2–C5 alkanes as the distance to the O&NG region decreases. Ethane and propane increase by > 10 ppb going from BCPH to SDP. The variability in mole fractions, as shown by the relative range of the box-whiskers, also increases as the distance to Weld County decreases. These data features are indicative of O&NG emission being a strong regional emission source for these VOCs. The spatial gradient becomes weaker with increasing alkane molecular weight. Other VOC, for example BTEX compounds, show more similar mole fractions and variability across all sites. This suggests that multiple sources are contributing to their emissions, which is consistent with findings from other previous NCFR studies (Pétron et al., 2012; Pétron et al., 2014).

Figure 5.

Box whisker plots of VOCs quantified in canister samples from the City of Boulder (BCPH) and the 4 Eastern Boulder County sites. The bottom right panel shows the (i/n) pentane ratio of the data. DOI: https://doi.org/10.1525/elementa.2019.00036.f5

Figure 5.

Box whisker plots of VOCs quantified in canister samples from the City of Boulder (BCPH) and the 4 Eastern Boulder County sites. The bottom right panel shows the (i/n) pentane ratio of the data. DOI: https://doi.org/10.1525/elementa.2019.00036.f5

Close modal

The observed gradients do not show an obvious correlation with population and traffic density, which is highest at the BCPH site (Table 1). As a matter of fact, median mole fractions are lowest at BCPH for all of the C2–C5 alkanes. None of the compounds in Figure 5 exhibit a gradient that correlates with population density. This behavior suggests that population and traffic density have a relatively small influence on the composition of the petroleum hydrocarbons at these sites. At all of the EBC network sites, median results for ethane and propane are all higher than those reported in data from 28 much larger urban U.S. cities (Baker et al., 2008; except for 1 case, i.e., BCPH propane compared to Los Angeles). Since petroleum hydrocarbons constitute the bulk of the overall VOC carbon, O&NG sources have a strong role in determining the regional VOC gradients. C2–C5 alkanes are the most abundantly emitted VOCs from O&NG operations (Pétron et al., 2012; Swarthout et al., 2013; Helmig et al., 2014; Pétron et al., 2014; Thompson et al., 2014); their highly elevated concentrations, and the clear gradients of these compounds as distance to the O&NG region decreases, demonstrate how selective these compounds are as O&NG emissions tracers, and how strongly these O&G industry emissions influence VOC composition and absolute levels not just in the immediate surrounding of well sites, but on large geographical scales in the NCFR.

The variability of the observed mole fractions for the C2–C5 alkanes, shown by the range of the box-whiskers in Figure 5, increased with decreasing distance to Weld County. High variability has indicated proximity to a source in other previous work (Russo et al., 2011; Sharma et al., 2016; Zheng et al., 2017). It is remarkable that this behavior is this pronounced despite that these results are from 3-day integrated sampling. Emission plumes are most compact at their source and gain volume during transport from dispersion and entrainment of the surrounding air. Their relatively small initial volume makes it more likely that factors such as wind direction will result in a nearby monitoring site “missing” the plume. For the same reason, if the center of the plume passes through the monitoring site, then high mole fractions will be observed because this portion of the plume has undergone the least mixing with surrounding air. Therefore, sampling at a site near a source will produce a dataset with a large amount of variation. The high degree of variability seen at the eastern network sites (especially at SDP) for the petroleum hydrocarbons can be deemed as an indicator of the closer proximity and stronger emission sources of these compounds at these sites.

The bottom right panel of Figure 5 shows the i/n-pentane isomeric ratios observed at each of the EBC sites. Gilman et al. (2013) investigated i/n-pentane isomeric ratios in different urban and O&NG influenced settings and defined an overall range of 0.9–3, with the lower end values being characteristic of air that is strongly influenced by O&NG emissions, and the higher end values being typical urban signatures. The highest i/n-pentane ratios were observed at BCPH. Median ratios of i/n-pentane at LFD, DAW, SLC, and SDP are 0.09–0.19 lower than at BCPH, and other percentile values of the i/n-pentane ratio similarly show a general trend towards smaller values from west-to-east and closer proximity to the DJB, which likely is indicative of the increasing influence of O&NG emissions on the pentanes composition along this transect.

The atmospheric photooxidation of VOCs by the hydroxyl (OH) radical fuels the diurnal ozone production. The contribution of the identified VOCs to the OH reactivity (OHR) and the spatial gradients of OHR from the identified VOCs was estimated for the C2–C8 alkanes, BTEX compounds, acetylene, and propyne. Please note that these calculations represent only a portion of the total OHR. There are many other atmospheric gases, that is, nitric oxide, carbon monoxide, alkenes, unsaturated biogenic VOCs, and partially oxidized VOCs (i.e., alcohols, aldehydes, ketones), that also contribute to the total OHR and that were not included in our assessment of the OH reactivity because of the lack of their observations. Therefore, analyses and interpretations presented here reflect the relative contribution among the identified VOCs species only and not of the overall reactivity of all VOCs present in ambient air.

OHR were determined by multiplying the mole fraction of each compound quantified during deployment periods by its OH reaction rate constants (Hatakeyama et al., 1986; Atkinson, 2003). OHRs of considered compounds were then totaled for each deployment period (ΣOHRs). The statistical distribution of ΣOHRs from all deployments at each site and the spatial gradient of the results from the 5 considered sites is shown in Figure 6. Figure 7 shows pie charts of the campaign median mole fractions by compound class and scaled to the number of carbons for each considered VOC (ppbC). The bottom row of graphs depicts the OHRs calculated from these median values. (Please note that the median mole fractions for each compound were not necessarily observed during the same deployment period; therefore, the sum of OHRs in the pie charts may not be exactly identical to the median value in the box plots in Figure 6).

Figure 6.

Box whisker plots of the sum of OH reactivity (OHR) of the quantified VOCs in the EBC network. Considered compounds were C2–C8 alkanes, BTEX, acetylene, and propyne. DOI: https://doi.org/10.1525/elementa.2019.00036.f6

Figure 6.

Box whisker plots of the sum of OH reactivity (OHR) of the quantified VOCs in the EBC network. Considered compounds were C2–C8 alkanes, BTEX, acetylene, and propyne. DOI: https://doi.org/10.1525/elementa.2019.00036.f6

Close modal
Figure 7.

VOC classes and their OH reactivity at EBC network sites. Top row: Contribution of measured compounds classified into C2–C5 alkanes, C6–C10 alkanes, BTEX, and alkenes + alkynes to the total, with each individual VOCs mole fraction multiplied by its number of carbons (ppbC) at the EBC sites. The sum of all compound class mole fractions (in ppbC) is given underneath the pie graphs. Bottom row: Median OH reactivity (OHR) of the same VOCs classes as in the upper row (OHR, lower row [s–1]). DOI: https://doi.org/10.1525/elementa.2019.00036.f7

Figure 7.

VOC classes and their OH reactivity at EBC network sites. Top row: Contribution of measured compounds classified into C2–C5 alkanes, C6–C10 alkanes, BTEX, and alkenes + alkynes to the total, with each individual VOCs mole fraction multiplied by its number of carbons (ppbC) at the EBC sites. The sum of all compound class mole fractions (in ppbC) is given underneath the pie graphs. Bottom row: Median OH reactivity (OHR) of the same VOCs classes as in the upper row (OHR, lower row [s–1]). DOI: https://doi.org/10.1525/elementa.2019.00036.f7

Close modal

Both analyses result in BCPH having the lowest total OHR and the lowest OHR variability. There is a general tendency of OHR variability and absolute percentile values increasing with the site being located further to the east. Figure 7 shows that the C2–C5 alkanes made up 72%–88% of the mole fractions of carbon (ppbC) of the quantified compounds and 52%–79% of ΣOHR at each site. The median total amount of VOC carbon of quantified species almost triples going from BCPH to SDP. Despite constituting only 4%–12% of the ppbC observed at each site, the alkynes and BTEX considered by this analysis accounted for 8%–22% of the ΣOHR at the EBC sites due to the relatively fast reaction rate constants with OH, demonstrating the reactive potential of the unsaturated VOCs to the ΣOHR.

3.3. Diurnal cycle of VOC

Figure 8 shows 7 diurnal cycles of VOCs measured by the in situ GC at INSTAAR. A striking feature is the high variability in the data and the lack of a consistent diurnal cycle. During most times, individual VOC species behave quite similarly in their dynamic behavior, suggesting that they mostly arise from common source emissions that are transported to the site. Mole fractions of some VOCs peaked during the day on 3 of the sampling dates. Propane and isomers of butane and pentane are the most frequently and most highly elevated compounds. Despite this sampling being from an inner-city site, and within ≈100 m of a very busy traffic intersection and a gas station, benzene and toluene are rarely above ≈0.2 ppb.

Figure 8.

VOC diurnal cycles measured by in-situ GC/FID at INSTAAR on several dates in July and August 2014. The x-axis displays the time (MDT) in 24-h format. The legend applies to all plots. Note that both the left and right y-axes scales vary by plot. DOI: https://doi.org/10.1525/elementa.2019.00036.f8

Figure 8.

VOC diurnal cycles measured by in-situ GC/FID at INSTAAR on several dates in July and August 2014. The x-axis displays the time (MDT) in 24-h format. The legend applies to all plots. Note that both the left and right y-axes scales vary by plot. DOI: https://doi.org/10.1525/elementa.2019.00036.f8

Close modal

The high variability seen in these data suggests that these light alkane VOCs likely arise from upwind emission sources that are variable in strength and/or that transport regimes vary such that polluted air from these sources is only transported to the site with particular air flow conditions. The absence of a distinct diurnal cycle for these compounds indicates that nearby emissions (within ≈1 day of transport time distance) have a stronger influence than daily photochemical oxidation and boundary layer dynamics, which both would result in concentration minima during midday hours. During the events on July 20 and July 31, the i/n-pentane ratio rises, which suggests urban or industrial sources (Gilman et al., 2013). This is supported by HYSPLIT back trajectories (Figures S4 and S6), which show that the air parcels that reached INSTAAR at the time of the event had passed through rural and urban areas in the 8 h preceding the event. In contrast, the VOC peak seen on July 26 is accompanied by a drop in the i/n-pentane ratio. The HYSPLIT back trajectory for this event (Figure S5) shows that the air parcel had circulated northeast of the city of Boulder and had passed through the southwestern part of Weld County. East to west transport as seen on July 26 is fairly common in this region, especially during mid to afternoon hours in the summer months, as elucidated in the introduction section.

The diurnal cycle data also do not show an obvious correlation with the diurnal traffic behavior (there are no peaks in VOCs concentrations seen during traffic rush hours). With mixed boundary layer depths being the highest during midday to afternoon hours, particularly during the midsummer, the large enhancements that were seen during midday hours in this surface sampling can be deemed representative of a large atmospheric column and therefore an indication of the strength of these sources.

Diurnal cycles were collected via PFP sampling at the LA site and are available in Figure S7. A comparison of the 3-h integrated samples with the PFP results is provided in Table S3. The PFP measurements support the canister measurements (Figure 2) in that propane is often higher than ethane in these elevation samples. This also supports the conclusion that the LA site, or potentially a larger area in this mountainous region, is exposed to local propane sources, possibly related to domestic gas usage. Expected diurnal cycles, that is, minimum mole fractions during the day and maximum mole fractions at night, are not clearly present in these diurnal cycles (Figure S7), nor are any obvious pollution transport events during these sample dates. The lack of consistent behavior suggests daily transport patterns are variable, at least in this limited number of samples. As for the canister data, the i/n-pentane isomeric ratios in the PFP data are close to the O&NG signature, which shows that these high elevation sites exhibit VOC signatures that are predominately influenced by VOC sources and the upslope transport of polluted air from the plains.

Observations from 2 monitoring networks showed clear elevation gradients for VOCs and spatial gradients of methane and VOCs in the plains between the City of Boulder and the eastern Boulder County line. Mole fractions of C2–C5 alkanes decreased with increasing elevation and increased with decreasing distance to the O&NG producing areas in Weld County to the east. The increase in absolute mole fractions was paralleled by an increase in the variability of mole fractions. Among the identified VOCs, the total VOC carbon was dominated by light alkane hydrocarbons throughout the networks. C2–C8 alkanes made up 88%–95% of the total VOCs carbon (ppbC) that was quantified by this study. This high contribution of O&NG petroleum hydrocarbons to total VOCs indicates that O&NG VOCs are abundant throughout the plains east of the Front Range and not only in the areas immediately surrounding O&NG development. This conclusion is further supported by the i/n-pentane ratios seen in the data. Medians of i/n-pentane ratios at the EBC sites ranged from 1.01 to 1.20, which is characteristic of O&NG influence. These ratios confirm the previously reported ratios of 1.1 at observation sites in Erie, CO, and Longmont, CO (Gilman et al., 2013; Thompson et al., 2014).

The spatial gradients, with increasingly elevated VOC toward the center of Weld County, as well as the predominance of petroleum hydrocarbons in the VOCs spectrum, are in accord with other previous VOC monitoring results. Unusually high levels of VOCs have been reported in measurements taken at the Boulder Atmospheric Observatory (BAO), which was located near the town of Erie in the southwest corner of Weld County (Pétron et al., 2012; Brown et al., 2013; Gilman et al., 2013; Swarthout et al., 2013; Oltmans et al., 2021). VOCs associated to O&NG sector emissions were the single largest VOC-reactivity contributing sector (40%–60% of total; with others being traffic, background, and secondary chemical production) throughout the morning hours (a biogenic factor also became important during 2–10 PM; Abeleira et al., 2017). Thompson et al. (2014) reported that O&NG petroleum hydrocarbon C2–C5 alkanes were enhanced by a factor of 18–77 relative to the regional background at residences in southwest Weld County. Unexpectedly high benzene mole fractions were reported in measurements taken near Platteville, approximately 28 km to the northeast of BAO, and further toward the center of the DJB (Halliday et al., 2016). The maximum observed value was 29.3 ppb. The authors implicated emissions from O&NG facilities to the southwest of the site. Surface automobile transects data collected inside Weld County showed tight n-alkane-methane and alkane-alkane correlations in the proximity of O&NG facilities (Pétron et al., 2014). These geographic patterns are not only present at the surface but also in the mixed boundary layer, as seen in aircraft measurements: Mixed boundary layer mole fractions of methane and the O&NG tracer ethane were highest within and downwind of the DJB (Pétron et al., 2014; Richter et al., 2015; Peischl et al., 2018; Flocke et al., 2020). The data and analyses presented here are consistent with this larger body of research that shows that the DJB O&NG operations constitute a strong source of atmospheric petroleum hydrocarbons that may dominate the VOCs spectrum and total carbon fraction from VOCs.

The higher time resolved data from the in situ monitoring within the City of Boulder showed evidence for pollution transport from the center of the O&NG basin into Boulder County to the west. Elevated values of the O&NG tracer ethane, low i/n-pentane values, and back trajectory analyses associated an event with the most enhanced levels of VOCs to O&NG emissions transported from the DJB. The observations presented here further demonstrate that the O&NG influence extends at the surface well beyond the boundaries of Weld County toward adjacent counties and into the eastern foothills of the Rocky Mountains: The influence of VOC emissions and transport was evident at the CO site, some 1,000 m in elevation above, and ≈20 km to the west of Boulder, demonstrating that the spatial and vertical extent of pollution transport from the NCFR plains extends into the foothills on the eastern flanks of the Rocky Mountains.

The OHR contribution from petroleum hydrocarbons generally increased with decreasing distance to the DJB. C2–C5 alkanes contributed 0.36–1.36 s–1, or 52%–79% to the OHR of the compounds that were quantified at the EBC sites. This is a sizable contribution given the relatively slow OH rate constants for these compounds. OHR from petroleum hydrocarbons was lowest at the urban/most densely populated City of Boulder site (BCPH). Our results for C2–C8 alkane mole fractions are within the range of the alkane portion of OHR determined by McDuffie et al. (2016), who found that petroleum alkanes were responsible for 50% of the OHR at the BAO. Using atmospheric modeling, these authors went a step further and estimated that these alkanes contributed ≈20% to the regional photochemical ozone production. The horizontal gradient in OHR with approximately 3 times higher values at the furthest eastern site points toward faster and more chemical reactivity and ozone production from these VOC species in the area represented by the eastern network sites than in the City of Boulder.

VOCs from industrial and mobile sources have decreased in most areas in the United States and Europe (Warneke et al., 2012; Monks et al., 2015; Rossabi and Helmig, 2018). Reductions in automobile VOCs emission have made a major contribution to these trends (Von Schneidemesser et al., 2010; Warneke et al., 2012). Declining automobile-related emissions have caused other emission categories, such as household and personal care products, to increase in their relative importance (Coggon et al., 2018; McDonald et al., 2018). Stricter regulations targeting O&NG emissions appear to have been successful in reducing average emissions per well unit; however, these reductions have likely been offset by the increase in the number of wells, with the overall result being only a modest overall decline in emissions (Peischl et al., 2018). A recent trend analysis of 8 years of BAO data concluded that actual total emission reductions are likely approximately only 1 quarter of regulatory inventory projections (Oltmans et al., 2021).

The results from our and these referenced research studies show that O&NG VOC emissions in the NCFR are a strong source of total VOC carbon and may contribute as much as half of the OHR. The VOC mixture in the NCFR contains a higher fraction and abundance of alkanes from O&NG sources compared to other urban regions (Zaragoza et al., 2017). These O&NG VOC emissions contribute to photochemical ozone production, and their atmospheric oxidation appears to be particularly crucial in the occurrences of peak ozone events in the NCFR (Zaragoza et al., 2017), contributing 30%–40% to the regional ozone production (Pfister et al., 2017a; Oltmans et al., 2019; Flocke et al., 2020). Cheadle et al. (2017) estimated that O&NG emissions contribute locally up to 30 ppb on high ozone days. While improvements in surface ozone have been achieved in many regions of the United States (Chang et al., 2017), surface ozone in the NCFR thus far has not shown clear signs of improvement (Bien and Helmig, 2018) despite notable reductions in nitrogen oxides emissions (Abeleira and Farmer, 2017; Bien and Helmig, 2018). Nine counties within the NCFR remain in nonattainment for the ozone NAAQS. Flocke et al. (2020), in their campaign overview publication, conclude about O&NG emissions that, “This emission sector is also identified as carrying the largest uncertainties in currently available emission inventories, both with regard to magnitude and emission factors.” The findings of our study illustrate another complexity to this evaluation: Ambient air concentrations of methane and VOCs show spatial gradients that reflect the distribution and heterogeneity of the location of emissions from the O&NG extraction operations and their transport under the diverse surface flow and meteorological conditions in the NCFR. This and other FRAPPE and DISCOVER-AQ studies reemphasize that O&NG sector emissions exert an important influence on the regional air quality not just within the DJB and Weld County but also into neighboring counties and into the Rocky Mountain Foothills.

The data from the surface sites monitoring are available in the FRAPPÉ/DISCOVER-AQ data archive (https://www-air.larc.nasa.gov/cgi-bin/ArcView/discover-aq.co-2014?GROUND-BOULDER=1).

The supplemental files for this article can be found as follows:

Figures S1–S7. Tables S1–S4. Docx.

Assistance from these INSTAAR/University of Colorado colleagues is acknowledged: John Ortega helped with the sampler installation and Brendan Blanchard with code development for data processing. Katie Smith and Reed Terrell assisted with the sampling and analysis work. Jana Milford and Mike Hannigan, Mechanical Engineering, University of Colorado, contributed to the planning of the study.

The East Boulder County sampling program and DH’s work on the manuscript were sponsored through contracts from Boulder County Public Health. The State of Colorado Department of Public Health and Environment (CDPHE) provided funding for the elevation gradient project. The National Science Foundation (NSF) funded the facilities and research at the Mountain Research Station site/Niwot Ridge under the Long-term Ecological Research program (award # DEB 1027341). FRAPPE was funded through a federal–state partnership between CDPHE and the National Science Foundation. DISCOVER-AQ was a NASA-led and funded campaign.

DH is the Editor-in-Chief of the Atmospheric Science domain of Elementa. He was not involved in any aspects of the editorial review process and evaluation of this manuscript.

Analyzed data, interpreted data, prepared graphics, and the manuscript: SR.

Prepared sampling equipment, collected and analyzed samples, and performed data quality control: JH.

Contributed to the design of the study: PM.

Set up GC system at INSTAAR, analyzed data, and contributed to the manuscript preparation: WW.

Contributed to the planning, design, sites selection, and execution of the study; analyzed data, prepared, and revised the manuscript: DH.

Abeleira
,
A
,
Pollack
,
IB
,
Sive
,
B
,
Zhou
,
Y
,
Fischer
,
EV
,
Farmer
,
DK
.
2017
.
Source characterization of volatile organic compounds in the Colorado Northern Front Range Metropolitan Area during spring and summer 2015
.
Journal of Geophysical Research
122
:
3595
3613
. DOI: http://dx.doi.org/10.1002/2016jd026227.
Abeleira
,
AJ
,
Farmer
,
DK
.
2017
.
Summer ozone in the Northern Front Range Metropolitan Area: Weekend-weekday effects, temperature dependences, and the impact of drought
.
Atmospheric Chemistry and Physics
17
:
6517
6529
. DOI: http://dx.doi.org/10.5194/acp-17-6517-2017.
Atkinson
,
R
.
2003
.
Kinetics of the gas-phase reactions of OH radicals with alkanes and cycloalkanes
.
Atmospheric Chemistry and Physics
3
:
2233
2307
.
Baker
,
AK
,
Beyersdorf
,
AJ
,
Doezema
,
LA
,
Katzenstein
,
A
,
Meinardi
,
S
,
Simpson
,
IJ
,
Blake
,
DR
,
Rowland
,
FS
.
2008
.
Measurements of nonmethane hydrocarbons in 28 United States cities
.
Atmospheric Environment
42
:
170
182
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2007.09.007.
Benedict
,
KB
,
Prenni
,
AJ
,
Sullivan
,
AP
,
Evanoski-Cole
,
AR
,
Fischer
,
EV
,
Callahan
,
S
,
Sive
,
BC
,
Zhou
,
Y
,
Schichtel
,
BA
,
Collett
Jr
JL
.
2018
.
Impact of Front Range sources on reactive nitrogen concentrations and deposition in Rocky Mountain National Park
.
PeerJ
6
:
e4759
. DOI: http://dx.doi.org/10.7717/peerj.4759.
Benedict
,
KB
,
Zhou
,
Y
,
Sive
,
BC
,
Prenni
,
AJ
,
Gebhart
,
KA
,
Fischer
,
EV
,
Evanoski-Cole
,
A
,
Sullivan
,
AP
,
Callahan
,
S
,
Schichtel
,
BA
,
Mao
,
H
.
2019
.
Volatile organic compounds and ozone in Rocky Mountain National Park during FRAPPE
.
Atmospheric Chemistry and Physics
19
:
499
521
. DOI: http://dx.doi.org/10.5194/acp-19-499-2019.
Bien
,
T
,
Helmig
,
D
.
2018
.
Changes in summertime ozone in Colorado during 2000-2015
.
Elementa: Science of the Anthropocene
6
:
1
25
. DOI: http://dx.doi.org/10.1525/elementa.300.
Boulder_AIR
.
2021
.
Current weather and atmospheric chemical conditions in Boulder, Colorado
.
Available at
https://www.bouldair.com/boulder.htm.
Accessed 7 December 2021
.
Bowman
,
FM
,
Seinfeld
,
JH
.
1994
.
Ozone productivity of atmospheric organics
.
Journal of Geophysical Research: Atmospheres
99
:
5309
5324
. DOI: http://dx.doi.org/10.1029/93jd03400.
Brodin
,
M
,
Helmig
,
D
,
Oltmans
,
S
.
2010
.
Seasonal ozone behavior along an elevation gradient in the Colorado front range mountains
.
Atmospheric Environment
44
:
5305
5315
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2010.06.033.
Brown
,
SS
,
Thornton
,
JA
,
Keene
,
WC
,
Pszenny
,
AA
,
Sive
,
BC
,
Dubé
,
WP
,
Wagner
,
NL
,
Young
,
CJ
,
Riedel
,
TP
,
Roberts
,
JM
,
VandenBoer
,
TC
,
Bahreini
,
R
,
Ozturk
,
F
,
Middlebrook
,
AM
,
Kim
,
S
,
Hubler
,
G
,
Wolfe
,
DE
.
2013
.
Nitrogen, aerosol composition, and halogens on a tall tower NACHTT: Overview of a wintertime air chemistry field study in the Front Range urban corridor of Colorado
.
Journal of Geophysical Research
118
:
8067
8085
. DOI: http://dx.doi.org/10.1002/jgrd.50537.
Bureau, U. S. C
.
2010
.
Boulder County, Colorado, edited, United States Census Bureau
,
Available at
https://www.census.gov/quickfacts/bouldercountycolorado.
Accessed 7 December 2021
.
Chameides
,
WL
,
Fehsenfeld
,
F
,
Rodgers
,
MO
,
Cardelino
,
C
,
Martinez
,
J
,
Parrish
,
D
,
Lonneman
,
W
,
Lawson
,
DR
,
Rasmussen
,
RA
,
Zimmerman
,
P
,
Greenberg
,
J
.
1992
.
Ozone precursor relationships in the ambient atmosphere
.
Journal of Geophysical Research
97
:
6037
6055
.
Chang
,
KL
,
Petropavlovskikh
,
I
,
Cooper
,
OR
,
Schultz
,
MG
,
Wang
,
T
,
Helmig
,
D
,
Lewis
,
A
.
2017
.
Regional trend analysis of surface ozone observations from monitoring networks in eastern North America, Europe and East Asia
.
Elementa: Science of the Anthropocene
5
. DOI: http://dx.doi.org/10.1525/elementa.243.
Cheadle
,
LC
,
Oltmans
,
SJ
,
Pétron
,
G
,
Schnell
,
RC
,
Mattson
,
EJ
,
Herndon
,
SC
,
Thompson
,
AM
,
Blake
,
DR
,
McClure-Begley
,
A
,
Helmig
,
D
,
Lamb
,
B
.
2017
.
Surface ozone in the Colorado Northern Front Range and the influence of oil and gas development during FRAPPE/DISCOVER-AQ in summer 2014
.
Elementa: Science of the Anthropocene
:
1
23
. DOI: http://dx.doi.org/10.1525/elementa.254.
COGCC
.
2014
.
County level production data
.
Available in the “Data: COGIS DataBase” section
.
Colorado Oil & Gas Conservation Commission
.
Available at
http://cogcc.state.co.us/.
Accessed 7 December 2021
.
COGCC
.
2016
.
Oil and natural gas wells
.
Available in the “Export of Data” section
.
Colorado Oil & Gas Conservation Commission
.
Available at
https://cogcc.state.co.us/dashboard.html#/dashboard.
Accessed 7 December 2021
.
Coggon
,
MM
,
McDonald
,
BC
,
Vlasenko
,
A
,
Veres
,
PR
,
Bernard
,
F
,
Koss
,
AR
,
Yuan
,
B
,
Gilman
,
JB
,
Peischl
,
J
,
Aikin
,
KC
,
DuRant
,
J
,
Warneke
,
C
,
Li
,
SM
,
Gouw
,
JAD
.
2018
.
Diurnal variability and emission pattern of decamethylcyclopentasiloxane D5 from the application of personal care products in two North American cities
.
Environmental Science & Technology
52
:
5610
5618
.
Colorado-DOT
.
2016
.
Online transportation information system, edited, Colorado Department of Transportation
.
Available at
http://dtdapps.coloradodot.info/otis/TrafficData#ui/0/0/0/criteria/7850//true/true/.
Accessed 7 December 2021
.
Darby
,
LS
,
Senff
,
CJ
,
Alvarez
,
RJ
II
,
Banta
,
RM
,
Bianco
,
L
,
Helmig
,
D
,
White
,
AB
.
2021
.
Spatial and temporal variability of ozone along the Colorado Front Range occurring over two days with contrasting wind flow
.
Elementa: Science of the Anthropocene
9
(
1
):
00146
. DOI: https://doi.org/10.1525/elementa.2020.00146.
EPA
.
2019
.
News releases from region 08. EPA reclassifies Denver area to “Serious” non-attainment for ozone
.
Available at
https://www.epa.gov/newsreleases/epa-reclassifies-denver-area-serious-nonattainment-ozone.
Accessed 7 December 2021
.
Fann
,
N
,
Baker
,
KR
,
Chan
,
EA
,
Eyth
,
A
,
Macpherson
,
A
,
Miller
,
E
,
Snyder
,
J
.
2018
.
Assessing human health PM2.5 and ozone impacts from US oil and natural gas sector emissions in 2025
.
Environmental Science & Technology
52
:
8095
8103
. DOI: http://dx.doi.org/10.1021/acs.est.8b02050.
Flocke
,
F
,
Pfister
,
G
,
Crawford
,
JH
,
Pickering
,
KE
,
Pierce
,
G
,
Bon
,
D
,
Reddy
,
P
.
2020
.
Air quality in the Northern Colorado Front Range metroarea: The Front Range air pollution andphotochemistry éxperiment FRAPPÉ
.
Journal of Geophysical Research
125
:
1
22
. DOI: http://dx.doi.org/e2019JD031197.
Gilman
,
JB
,
Lerner
,
BM
,
Kuster
,
WC
,
De Gouw
,
JA
.
2013
.
Source signature of volatile organic compounds from oil and natural gas operations in northeastern Colorado
.
Environmental Science & Technology
47
:
1297
1305
. DOI: http://dx.doi.org/10.1021/es304119a.
Halliday
,
HS
,
Thompson
,
AM
,
Wisthaler
,
A
,
Blake
,
DR
,
Hornbrook
,
RS
,
Mikoviny
,
T
,
Müller
,
M
,
Eichler
,
P
,
Apel
,
EC
,
Hills
,
AJ
.
2016
.
Atmospheric benzene observations from oil and gas production in the Denver-Julesburg Basin in July and August 2014
.
Journal of Geophysical Research
121
:
11055
11074
. DOI: http://dx.doi.org/10.1002/2016jd025327.
Hatakeyama
,
S
,
Washida
,
N
,
Akimoto
,
H
.
1986
.
Rate constants and mechanisms for the reaction of OH (OD) radicals with acetylene, propyne, and 2-butyne in air at 297 +/- 2 K
.
Journal of Physical Chemistry
90
:
173
178
. DOI: http://dx.doi.org/10.1021/j100273a039.
Helmig
,
D
.
1997
.
Ozone removal techniques in the sampling of atmospheric volatile organic trace gases
.
Atmospheric Environment
31
:
3635
3651
.
Helmig
,
D
.
2020
.
Policy bridge: Air quality impacts from oil and natural gas development in Colorado
.
Elementa: Science of the Anthropocene
8
:
1
33
. DOI: http://dx.doi.org/10.1525/elementa.398.
Helmig
,
D
,
Thompson
,
CR
,
Evans
,
J
,
Boylan
,
P
,
Hueber
,
J
,
Park
,
JH
.
2014
.
Highly elevated atmospheric levels of volatile organic compounds in the Uintah Basin, Utah
.
Environmental Science & Technology
48
:
4707
4715
. DOI: http://dx.doi.org/10.1021/es405046r.
Huff
,
J
.
2007
.
Benzene-induced cancers: Abridged history and occupational health impact
.
International Journal of Occupational and Environmental Health
13
:
213
221
.
Johnson
,
R
,
Toth
,
JJ
.
1982
.
A climatology of the July 1981 surface flow over northeast Colorado
.
Department of Atmospheric Science, Colorado State University
,
Paper No. 342
.
Available at
https://mountainscholar.org/bitstream/handle/10217/335/0342_Bluebook.pdf;sequence=1.
Kille
,
N
,
Chiu
,
R
,
Frey
,
M
,
Hase
,
F
,
Sha
,
MK
,
Blumenstock
,
T
,
Hannigan
,
JW
,
Orphal
,
J
,
Bon
,
D
,
Volkamer
,
R
.
2019
.
Separation of methane emissions from agricultural and natural gas sources in the Colorado Front Range
.
Geophysical Research Letters
46
:
3990
3998
, DOI: http://dx.doi.org/10.1029/2019gl082132.
Macey
,
GP
,
Breech
,
R
,
Chernaik
,
M
,
Cox
,
C
,
Larson
,
D
,
Thomas
,
D
,
Carpenter
,
DO
.
2014
.
Air concentrations of volatile compounds near oil and gas production: a community-based exploratory study
.
Environmental Health
13
: DOI: http://dx.doi.org/10.1186/1476-069x-13-82.
McDonald
,
BC
,
De Gouw
,
JA
,
Gilman
,
JB
,
Jathar
,
SH
,
Akherati
,
A
,
Cappa
,
CD
,
Jimenez
,
JL
,
Lee-Taylor
,
J
,
Hayes
,
PL
,
McKeen
,
SA
,
Cui
,
YY
,
Kim
,
SW
,
Gentner
,
DR
,
Isaacman-VanWertz
,
G
,
Goldstein
,
AH
,
Harley
,
RA
,
Frost
,
GJ
,
Roberts
,
JM
,
Ryerson
,
TB
,
Trainer
,
M
.
2018
.
Volatile chemical products emerging as largest petrochemical source of urban organic emissions
.
Science
359
:
760
764
. DOI: http://dx.doi.org/10.1126/science.aaq0524.
McDuffie
,
EE
,
Edwards
,
PM
,
Gilman
,
JB
,
Lerner
,
BM
,
Dubé
,
WP
,
Trainer
,
M
,
Wolfe
,
DE
,
Angevine
,
WM
,
deGouw
,
J
,
Williams
,
EJ
,
Tevlin
,
AG
,
Murphy
,
JG
,
Fischer
,
EV
,
McKeen
,
S
,
Ryerson
,
TB
,
Peischl
,
J
,
Holloway
,
JS
,
Aikin
,
K
,
Langford
,
AO
,
Senff
,
CJ
,
Alvarez
,
RJ
,
Hall
,
SR
,
Ullmann
,
K
,
Lantz
,
KO
,
Brown
,
SS
.
2016
.
Influence of oil and gas emissions on summertime ozone in the Colorado Northern Front Range
.
Journal of Geophysical Research: Atmospheres
121
:
8712
8729
. DOI: http://dx.doi.org/10.1002/2016jd025265.
McKenzie
,
LM
,
Blair
,
B
,
Hughes
,
J
,
Allshouse
,
WB
,
Blake
,
NJ
,
Helmig
,
D
,
Milmoe
,
P
,
Halliday
,
H
,
Blake
,
DR
,
Adgate
,
JL
.
2018
.
Ambient nonmethane hydrocarbon levels along Colorado’s Northern Front Range: Acute and chronic health risks
.
Environmental Science & Technology
52
:
4514
4525
. DOI: http://dx.doi.org/10.1021/acs.est.7b05983.
McKenzie
,
LM
,
Guo
,
R
,
Witter
,
RZ
,
Savitz
,
DA
,
Newman
,
LS
,
Adgate
,
JL
.
2014
.
Birth outcomes and maternal residential proximity to natural gas development in rural Colorado
.
Environmental Health Perspectives
122
:
412
417
. DOI: http://dx.doi.org/10.1289/ehp.1306722.
Monks
,
PS
,
Archibald
,
AT
,
Colette
,
A
,
Cooper
,
O
,
Coyle
,
M
,
Derwent
,
R
,
Fowler
,
D
,
Granier
,
C
,
Law
,
KS
,
Mills
,
GE
,
Stevenson
,
DS
,
Tarasova
,
O
,
Thouret
,
V
,
von Schneidemesser
,
E
,
Sommariva
,
R
,
Wild
,
O
,
Williams
,
ML
.
2015
.
Tropospheric ozone and its precursors from the urban to the global scale from air quality to short-lived climate forcer
.
Atmospheric Chemistry and Physics
15
:
8889
8973
. DOI: http://dx.doi.org/10.5194/acp-15-8889-2015.
Mount
,
GH
,
Harder
,
JW
.
1995
.
The measurement of tropospheric trace gases at Fritz Peak Observatory, Colorado, by long-path absorption: OH and ancillary gases
.
Journal of Geophysical Research
52
:
3342
3353
.
Oltmans
,
SJ
,
Cheadle
,
LC
,
Helmig
,
D
,
Angot
,
H
,
Pétron
,
G
,
Montzka
,
SA
,
Dlugokencky
,
EJ
,
Miller
,
B
,
Hall
B
,
Schnell
,
RC
,
Kofler
,
J
,
Wolter
,
S
,
Crotwell
,
M
,
Siso
,
C
,
Tans
,
P
,
Andrews
,
A
.
2021
.
Atmospheric oil and natural gas hydrocarbon trends in the Northern Colorado Front Range are notably smaller than inventory emissions reductions
.
Elementa: Science of the Anthropocene
9
:
1
16
. DOI: https://doi.org/10.1525/elementa.2020.00136.
Oltmans
,
SJ
,
Cheadle
,
LC
,
Johnson
,
BJ
,
Schnell
,
RC
,
Helmig
,
D
,
Thompson
,
AM
,
Cullis
,
P
,
Hall
,
E
,
Jordan
A
,
Sterling
,
C
,
McClure-Begley
,
A
.
2019
.
Boundary layer ozone in the Northern Colorado Front Range in July-August 2014 during FRAPPE and DISCOVER-AQ from vertical profile measurements
.
Elementa: Science of the Anthropocene
7
:
1
14
. DOI: https://doi.org/10.1525/elementa.345.
Peischl
,
J
,
Eilerman
,
SJ
,
Neuman
,
JA
,
Aikin
,
KC
,
de Gouw
,
J
,
Gilman
,
JB
,
Herndon
,
SC
,
Nadkarni
,
R
,
Trainer
,
M
,
Warneke
,
C
,
Ryerson
,
TB
.
2018
.
Quantifying methane and ethane emissions to the atmosphere from Central and Western U.S. oil and natural gas production regions
.
Journal of Geophysical Research
123
:
7725
7740
. DOI: http://dx.doi.org/10.1029/2018JD028622.
Pétron
,
G
,
Frost
,
G
,
Miller
,
BR
,
Hirsch
,
AI
,
Montzka
,
SA
,
Karion
,
A
,
Trainer
,
M
,
Sweeney
,
C
,
Andrews
,
AE
,
Miller
,
L
,
Kofler
,
J
,
Bar-Ilan
,
A
,
Dlugokencky
,
EJ
,
Patrick
,
L
,
Moore
,
CT
,
Ryerson
,
TB
,
Siso
,
C
,
Kolodzey
,
W
,
Lang
,
PM
,
Conway
,
T
,
Novelli
,
P
,
Masarie
,
K
,
Hall
,
B
,
Guenther
,
D
,
Kitzis
,
D
,
Miller
,
J
,
Welsh
,
D
,
Wolfe
,
D
,
Neff
,
W
,
Tans
,
P
.
2012
.
Hydrocarbon emissions characterization in the Colorado Front Range: A pilot study
.
Journal of Geophysical Research
117
:
1
19
. DOI: http://dx.doi.org/10.1029/2011jd016360.
Pétron
,
G
,
Karion
,
A
,
Sweeney
,
C
,
Miller
,
BR
,
Montzka
,
SA
,
Frost
,
GJ
,
Trainer
,
M
,
Tans
,
P
,
Andrews
,
A
,
Kofler
,
J
,
Helmig
,
D
,
Guenther
,
D
,
Dlugokencky
,
E
,
Lang
,
P
,
Newberger
,
T
,
Wolter
,
S
,
Hall
,
B
,
Novelli
,
P
,
Brewer
,
A
,
Conley
,
S
,
Hardesty
,
M
,
Banta
,
R
,
White
,
A
,
Noone
,
D
,
Wolfe
,
D
,
Schnell
,
R
.
2014
.
A new look at methane and nonmethane hydrocarbon emissions from oil and natural gas operations in the Colorado Denver-Julesburg Basin
.
Journal of Geophysical Research
119
:
6836
6852
. DOI: http://dx.doi.org/10.1002/2013jd021272.
Pfister
,
G
,
Flocke
,
F
,
Hornbrook
,
R
,
Orlando
,
J
,
Lee
,
S
,
Schroeder
,
J
.
2017
a
31
Jul
.
Process-based and regional source impact analysis for FRAPPE and DISCOVER-AQ 2014
.
Final Report to the Colorado Department of Public Health and Environment
.
Available at
https://www.colorado.gov/airquality/tech_doc_repository.aspx?action=open&file=FRAPPE-NCAR_Final_Report_July2017.pdf.
Pfister
,
GG
,
Reddy
,
PJ
,
Barth
,
MC
,
Flocke
,
FF
,
Fried
,
A
,
Herndon
,
SC
,
Sive
,
BC
,
Sullivan
,
JT
,
Thompson
,
AM
,
Yacovitch
,
TI
,
Weinheimer
,
AJ
.
2017
b.
Using observations and source-specific model tracers to characterize pollutant transport during FRAPPÉ and DISCOVER-AQ
.
Journal of Geophysical Research
122
:
10,510-510,538
. DOI: http://dx.doi.org/10.1002/2017JD027257.
Pollmann
,
J
,
Helmig
,
D
,
Hueber
,
J
,
Tanner
,
D
,
Tans
,
PP
.
2006
.
Evaluation of solid adsorbent materials for cryogen-free trapping - gas chromatographic analysis of atmospheric C2-C6 non-methane hydrocarbons
.
Journal of Chromatography A
1134
:
1
15
.
Population
.
2016
.
Boulder County population
.
Available at
https://worldpopulationreview.com/us-counties/co/boulder-county-population.
Accessed 7 December 2021
.
Rhoderick
,
GC
,
Duewer
,
DL
,
Apel
,
E
,
Baldan
,
A
,
Hall
,
B
,
Harling
,
A
,
Helmig
,
D
,
Heo
,
GS
,
Hueber
,
J
,
Kim
,
ME
,
Kim
,
YD
.
2014
.
International comparison of a hydrocarbon gas standard at the picomol per mol level
.
Analytical Chemistry
86
:
2580
2589
. DOI: http://dx.doi.org/10.1021/ac403761u.
Richter
,
D
,
Weibring
,
P
,
Walega
,
JG
,
Fried
,
A
,
Spuler
,
SM
,
Taubman
,
MS
.
2015
.
Compact highly sensitive multi-species airborne mid-IR spectrometer
.
Applied Physics B
119
:
119
131
. DOI: http://dx.doi.org/10.1007/s00340-015-6038-8.
Rossabi
,
S
,
Helmig
,
D
.
2018
.
Changes in atmospheric butanes and pentanes and their isomeric ratios in the Continental United States
.
Journal of Geophysical Research
123
:
3772
3790
. DOI: http://dx.doi.org/10.1002/2017jd027709.
Russo
,
RS
,
White
,
ML
,
Zhou
,
Y
,
Haase
,
KB
,
Ambrose
,
JL
,
Conway
,
L
,
Mentis
,
E
,
Talbot
,
R
,
Sive
,
BC
.
2011
. Spatial variation, sources, and emission rates of volatile organic compounds over the Northeastern U.S., in
Mazzeo
,
N
ed.,
Air quality-models and applications
.
Durham, NH
:
Climate Change Research Center, Institute for the Study of Earth, Oceans, and Space, University of New Hampshire
:
233
260
.
Sharma
,
S
,
Giri
,
B
,
Patel
,
KS
.
2016
.
Ambient volatile organic compounds in the atmosphere of industrial central India
.
Journal of Atmospheric Chemistry
73
:
381
395
. DOI: http://dx.doi.org/10.1007/s10874-016-9329-5.
Sovacool
,
BK
.
2014
,
Cornucopia or curse? Reviewing the costs and benefits of shale gas hydraulic fracturing fracking
,
Renewable & Sustainable Energy Reviews
37
:
249
-
264
, DOI: http://dx.doi.org/10.1016/j.rser.2014.04.068.
Swarthout
,
RF
,
Russo
,
RS
,
Zhou
,
Y
,
Hart
,
AH
,
Sive
,
BC
.
2013
.
Volatile organic compound distributions during the NACHTT campaign at the Boulder Atmospheric Observatory: Influence of urban and natural gas sources
.
Journal of Geophysical Research
118
:
10614
10637
.
Tanner
,
D
,
Helmig
,
D
,
Hueber
,
J
,
Goldan
,
P
.
2006
.
Gas chromatography system for the automated, unattended, and cryogen-free monitoring of C2 to C6 non-methane hydrocarbons in the remote troposphere
.
Journal of Chromatography A
1111
:
76
88
. DOI: http://dx.doi.org/10.1016/j.chroma.2006.01.100.
Thompson
,
CR
,
Hueber
,
J
,
Helmig
,
D
.
2014
.
Influence of oil and gas emissions on ambient atmospheric non-methane hydrocarbons in residential areas of Northeastern Colorado
.
Elementa: Science of the Anthropocene
3
:
1
17
. DOI: http://dx.doi.org/10.12952/journal.elementa.000035.
Tzompa-Sosa
,
ZA
,
Mahieu
,
E
,
Franco
,
B
,
Keller
,
CA
,
Turner
,
AJ
,
Helmig
,
D
,
Fried
,
A
,
Richter
,
D
,
Weibring
,
P
,
Walega
,
J
,
Yacovitch
,
TI
.
2017
.
Revisiting global fossil fuel and biofuel emissions of ethane
.
Journal of Geophysical Research-Atmospheres
122
:
2493
2512
. DOI: http://dx.doi.org/10.1002/2016jd025767.
Verma
,
DK
,
Johnson
,
DM
,
McLean
,
JD
.
2000
.
Benzene and total hydrocarbon exposures in the upstream petroleum oil and gas industry
.
American Industrial Hygiene Association Journal
61
:
255
263
, DOI: http://dx.doi.org/10.1080/15298660008984534.
Von Schneidemesser
,
E
,
Monks
,
PS
,
Plass-Duelmer
,
C
.
2010
.
Global comparison of VOC and CO observations in urban areas
.
Atmospheric Environment
44
:
5053
5064
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2010.09.010.
Warneke
,
C
,
De Gouw
,
JA
,
Holloway
,
JS
,
Peischl
,
J
,
Ryerson
,
TB
,
Atlas
,
E
,
Blake
,
D
,
Trainer
,
M
,
Parrish
,
DD
.
2012
.
Multiyear trends in volatile organic compounds in Los Angeles, California: Five decades of decreasing emissions
.
Journal of Geophysical Research-Atmospheres
117
. DOI: http://dx.doi.org/D00v1710.1029/2012jd017899.
WCC-VOC
.
2016
.
World Calibration Centre for Volatile Organic Compounds
.
Karlsruhe Institute of Technology
.
Available at
http://www.imk-ifu.kit.edu/wcc-voc/.
Accessed 7 December 2021
.
Zaragoza
,
J
,
Callahan
,
S
,
McDuffie
,
EE
,
Kirkland
,
J
,
Brophy
,
P
,
Durrett
,
L
,
Farmer
,
DK
,
Zhou
,
Y
,
Sive
,
B
,
Flocke
,
F
,
Pfister
,
G
.
2017
.
Observations of acyl peroxy nitrates during the Front Range air pollution and photochemistry experiment FRAPPE
.
Journal of Geophysical Research-Atmospheres
122
:
12416
12432
. DOI: http://dx.doi.org/10.1002/2017jd027337.
Zheng
,
C
,
Shen
,
J
,
Zhang
,
Y
,
Huang
,
W
,
Zhu
,
X
,
Wu
,
X
,
Chen
,
L
,
Gao
,
X
,
Cen
,
K
.
2017
.
Quantitative assessment of industrial VOC emissions in China: Historical trend, spatial distribution, uncertainties, and projection
.
Atmospheric Environment
150
:
116
125
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2016.11.023.

How to cite this article: Rossabi, S, Hueber, J, Wang, W, Milmoe, P, Helmig, D. 2021. Spatial distribution of atmospheric oil and natural gas volatile organic compounds in the Northern Colorado Front Range. Elementa: Science of the Anthropocene 9(1). DOI: https://doi.org/10.1525/elementa.2019.00036

Associate Editor: Michael E. Chang, Georgia Institute of Technology, Atlanta, GA, USA

Knowledge Domain: Atmospheric Science

Part of an Elementa Forum: Oil and Natural Gas Development: Air Quality, Climate Science, and Policy

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.

Supplementary Material