Articles | Volume 23, issue 13
Technical note
10 Jul 2023
Technical note |  | 10 Jul 2023

Technical note: Isolating methane emissions from animal feeding operations in an interfering location

Megan E. McCabe, Ilana B. Pollack, Emily V. Fischer, Kathryn M. Steinmann, and Dana R. Caulton

Agricultural emissions, including those from concentrated animal feeding operations (CAFOs) for beef and dairy cattle, make up a large portion of the United States' total greenhouse gas (GHG) emissions. However, many CAFOs reside in areas where methane (CH4) from oil and natural gas (ONG) complicates the quantification of CAFO emissions. Traditional approaches to quantify emissions in such regions often relied on inventory subtraction of other known sources. We compare the results of two approaches to attribute the CAFO CH4 emission rate from the total CH4 emission rate derived from an aircraft mass balance technique. These methods make use of the mixing ratio data of CH4, ethane (C2H6), and ammonia (NH3) that were collected simultaneously in-flight downwind of CAFOs in northeastern Colorado. The first approach, the subtraction method (SM), is similar to inventory subtraction, except the amount to be removed is derived from the observed C2H6 to CH4 ratio rather than an inventory estimate. The results from this approach showed high uncertainty, primarily due to how error propagates through subtraction. Alternatively, multivariate regression (MVR) can be used to estimate CAFO CH4 emissions using the NH3 emission rate and an NH3 to CH4 ratio. These results showed significantly less uncertainty. We identified criteria to determine the best attribution method; these criteria can support attribution in other regions. The final emission estimates for the CAFOs presented here were 13 ± 3 g of CH4 per head per hour and 13 ± 2 g of NH3 per head per hour. These estimates are higher than the inventory of the US Environmental Protection Agency (EPA) and previous studies highlighting the need for more measurements of CH4 and NH3 emission rates.

1 Introduction

Livestock produce large amounts of greenhouse gases (GHGs) and reactive nitrogen species, including methane (CH4) and ammonia (NH3), through enteric fermentation and waste generation. Ruminant animals (e.g., cattle, buffalo, sheep, goats, and camels) constitute a significant source of CH4 as their digestive systems break down coarse plant material through microbial fermentation in their rumen stomach (large frontal stomach) and subsequently release the produced gas (CH4). From 1990 to 2019, CH4 emissions from enteric fermentation grew 8.4 % in the United States, making agriculture the largest source of US anthropogenic CH4 emissions in 2020 (EPA2022). Waste and manure management are also significant emission sources of CH4 and NH3. Together, enteric fermentation and manure management account for more than 30 % of US anthropogenic CH4 emissions (Maasakkers et al.2016; EPA2022); CH4 is a crucial GHG due to its high global warming potential (Myhre et al.2013; Moumen et al.2016; Smith et al.2021).

Large uncertainty remains around the magnitude of livestock emissions. Beef and dairy cattle are confined to feedlots where they are fed and kept in tight areas to process the animals efficiently. These feedlots are known as concentrated animal feeding operations (CAFOs). The compactness of CAFOs has been shown to create significant emissions (Golston et al.2020; Hacker et al.2016; Eilerman et al.2016; Staebler et al.2009). However, observations indicated large variability in these CH4 emissions, which creates uncertainty in cumulative estimates of agricultural emissions (Golston et al.2020). There are many different factors determining the amount of CH4 released from enteric fermentation and manure management, and practices may vary from farm to farm (EPA2022; Maasakkers et al.2016). However, interfering sources of CH4 and NH3, such as oil and natural gas (ONG), waste pools, and landfills, may also complicate measurements from individual CAFOs. Improving the methodology for isolating CAFO signals from other interfering sources will allow more accurate measurements and provide new information to constrain greenhouse gas emissions from the agriculture sector.

When CH4 concentration data are used alone for the quantification of emissions in regions with multiple CH4 sources, there is not enough information to distinguish the contributing sources. However, previous studies have attributed CH4 emissions in complicated regions using a variety of methods. A few examples include attribution by subtracting inventory data (Caulton et al.2014; Peischl et al.2015, 2018), collecting ground-based isotope data to attribute CH4 signals (Townsend-Small et al.2016), using C2H6 as a tracer to subtract or attribute the ONG fraction (Mielke-Maday et al.2019), and using multivariate regression (MVR) with independent gas tracers to attribute sectors (Kille et al.2019; Pollack et al.2022).

There are many considerations when determining which attribution method is most appropriate, and often attribution is done based on the data available rather than the ideal methodology. Attribution using the subtraction of a value determined by inventory or calculated ratio is best done for larger regions due to the variability of emissions from individual sources. For ONG, individual sites can release different amounts of CH4 and ethane (C2H6) compared to other wells and compressors (Yacovitch et al.2014; Zimmerle et al.2022). When quantifying emissions in a small region (as in this work), an estimate of emissions from each type of ONG source or a local C2H6 : CH4 ratio near the CAFO is required to accurately separate the CH4 emissions into contributions from the CAFO and nearby ONG activities.

Another concern with using a tracer like C2H6 in isolation to estimate the contribution of ONG is that in complicated regions, individual CH4 signals become mixed, making it possible that the ratio observed by the aircraft is not representative of the original ratio at the ground. This is the theory behind tracer release, for example, where a tracer gas is released at a known rate near a source of interest and used to back out the source emission rate (Roscioli et al.2015). The added gas does not have to be introduced exactly at the source emission point (which may be unknown), provided sampling occurs far enough downwind where the species are well mixed. Typically, this means that the tracer gas must be released within 100 m of the source and measured > 500 m downwind (Roscioli et al.2015). For airborne data, there are many situations where we would expect signals to be mixed, complicating the use of airborne ratio analysis. Townsend-Small et al. (2016) circumvented this by combining ground-based isotope ratios (which show distinct ratios for particular sources) with aircraft data (which showed only one ratio). However, because they used a single isotope ratio, there is large uncertainty in their results stemming from the single isotope ratio that must be attributed to multiple contribution signals (one equation, multiple unknowns). This is similar to the difficulties that Smith et al. (2015) encountered using C2H6 as a tracer in a complicated region with multiple C2H6 : CH4 ratios.

On the other hand, an approach like MVR, which makes use of multiple tracer gases, requires sufficient data and is subject to its own sensitivities (Kille et al.2019). MVR is best used when there are multiple tracer gases that can be treated as independent variables and one dependent variable (i.e., CH4). The more independent variables that are included, the more data are needed to produce statistically significant results. These independent variables can be thought of as independent source terms. Although MVR is generally an effective method, it is important to note that it may not be appropriate in all situations. Specifically, using MVR with too many independent variables in a region where there are only a few sources may result in misleading or inaccurate results. However, using MVR with a limited number of independent variables that correspond to known sources can yield reliable results for those specific sources. Even in cases where there is a possibility of additional sources, fewer independent variables will produce more robust results and other sources can be treated in the extra term or background.

Here, we demonstrate a methodology to isolate and quantify emission rates for CAFOs in the northeastern Colorado Front Range (NCFR) using airborne measurements of CH4, C2H6, and NH3. We investigated two methods for CAFO CH4 emission isolation: (1) a subtraction method (SM) using the C2H6 : CH4 ratio and (2) an MVR method using CH4, C2H6, and NH3. This study focuses on the NCFR, where there is a high density of large CAFOs. Figure 1 shows a map of the NCFR with CAFOs for beef cattle, dairy cattle, chickens, sheep, and swine in the area. The area is dominated by CAFOs for beef cattle and dairy cattle. However, the NCFR contains a large mixture of CH4 emissions due to the high production of ONG. In this region, prior estimates indicate that natural gas accounts for 38.5 % of the state-wide CH4 emissions, while agriculture accounts for 22.3 % of the state-wide CH4 emissions (Arnold et al.2014). The NCFR is home to the Denver–Julesburg Basin (DJB), with over 52 000 ONG wells, and it has an abundance of compressors and processing plants, many of which are in close proximity to CAFOs (Higley and Cox2007).

Figure 1Map of northeastern Colorado. CAFOs are colored by animal type and sized by maximum animal units. Animal units are equivalent to 1 live beef cattle such that 1 head of beef cattle = 0.7 dairy cattle = 2.5 swine = 10 sheep = 100 poultry (CDPHE2017). Note that the maximum animal units represent the maximum animal capacity of a given facility and not necessarily the actual number of animals present at that facility at the time of sampling. CAFO data as of 2017 are registered with the Colorado Department of Public Health and Environment (CDPHE2017).

2 Materials and methods

2.1 Data collection

This study was conducted in the NCFR near Greeley, CO; the analysis includes data collected over Weld, Morgan, Logan, Larimer, and Washington counties. There are many CAFOs within these five counties with an area-wide maximum capacity of > 1 million heads of cattle (United States Department of Agriculture [USDA]2018). This study used the University of Wyoming King Air (UWKA), which is a relatively small research aircraft capable of flying at low altitudes (100 m a.g.l.) and slow flight speeds (95 m s−1). The UWKA is a national aircraft research facility owned and operated by the University of Wyoming. Flights departed from and returned to the Laramie Airport in Laramie, Wyoming (KLAR). Data outside of the NCFR were removed to ensure regional data only. This was done by screening out everything north and west of 41 latitude, 105.25 longitude. Three flights were performed in November 2019, departing around 12:00 MST (Mountain Standard Time) and lasting 2–4 h.

The UWKA was instrumented to measure CH4, NH3, and C2H6 mixing ratios. A Picarro G2401-m flight-ready analyzer measured CH4, carbon monoxide (CO), carbon dioxide (CO2), and water vapor (H2O) at 0.25 Hz through infrared cavity ring-down spectroscopy (Crosson2008). This Picarro and other Picarro models were tested previously and found to be stable and suitable for airborne field measurements (Richardson et al.2012). Two separate Aerodyne commercial quantum–cascade tunable infrared laser direct absorption spectrometers (QC-TILDAS) measured NH3 and C2H6. These instruments are described in detail in Yacovitch et al. (2014), Pollack et al. (2019), and Pollack et al. (2022). Measurements of NH3 are collected at 10 Hz and averaged to 1 Hz for reporting and this analysis; C2H6 measurements are collected and reported at 1 Hz. All chemical data were adjusted for time lag between instruments and further averaged to 0.25 Hz for emission calculations. Instruments were calibrated on the ground before and after flights. Instrument zeros were routinely measured in flight by overblowing the instrument inlets with a bottled source of synthetic “zero” air. Other in situ measurements from the UWKA standard instrument package included pressure, temperature, three-dimensional winds, GPS position, aircraft altitude, and heading. The wind speed is measured at 25 Hz then averaged to 0.25 Hz. The precision on the reported wind speed is 0.14 m s−1 with an expected wind direction precision of  5 (Strauss et al.2015).

2.2 Flight patterns

Flights were designed to identify the best flight patterns to simultaneously quantify CH4 and NH3 emission fluxes, CH4 to NH3 ratios, and NH3 deposition downwind of CAFOs. Prior to each flight, forecast meteorology was used to identify the ideal CAFOs to sample based on prevailing wind direction, isolation from other CAFO plumes, and other logistical constraints (e.g., proximity to urban areas, towers, and airports). Once airborne, the pilot conducted a vertical profile to characterize the mixed boundary layer (MBL). Selected CAFOs were located by coordinates and, once close enough, by sight. Once a selected CAFO was identified by sight, the UWKA pilot would perform a visual safety inspection of the area and then fly up or down in a spiral pattern centered on the selected CAFO. The aircraft proceeded to circle the target CAFO at a low altitude to confirm in situ enhancements of CH4 and NH3 mixing ratios and the direction of the outflow plume. Flight altitudes ranged from 0.1 to 3 km a.g.l., and near-CAFO flight altitudes depended on safety constraints.

Ammonia (NH3) deposition calculations are not the focus of this work; however, observations of deposition require multiple downwind observations. NH3 : CH4 ratios should be calculated near the source, as documented in Pollack et al. (2022), which used the same dataset as this work. We investigated both spiral (which can be completed quickly) and horizontal (which require much more flying time) transect flight patterns, shown in Fig. 2. The spiral patterns can be executed in quick succession, provide information as to the upwind background, and provide multiple downwind distances for analysis of NH3 deposition. However, the spiral transects were found to be undesirable for quantifying emissions because the aircraft could not get far enough outside the plume to characterize background conditions. For example, the further-downwind spiral transects had high enhancements of CH4 and NH3 for the full width of the spiral, indicating that the UWKA did not leave the plume. The horizontal transects, while ideal for emission calculation, did not provide much information about the evolution of downwind NH3 deposition. In order to sample the full plume efficiently and provide multiple downwind observations, racetrack patterns or boxes were later identified as a preferred approach for future sampling.

Figure 2Observed CAFO flight path color-filled with CH4 (ppbv). Animal operations are indicated by the different colored circles as in Fig. 1. Green dots represent ONG wells; data of ONG are as of 2015 (Colorado Department of Natural Resources Oil & Gas Conservation Commission2016).

2.3 Emission calculations

CAFOs suited for emission calculations were identified based on the following requirements: (1) enhancements of CH4 and NH3 above background conditions, (2) flight path includes multiple transects downwind at different altitudes within the MBL, and (3) enough data near the CAFO, but outside the CAFO plume, to characterize background mixing ratios. Only one CAFO sampled on 13 November 2019 satisfied the stated requirements to be suitable for emission quantification and is used in the remainder of this work. The flight on 13 November 2019 (denoted as F2 from hereon) occurred during a period with strong and steady winds (average wind speed of 8.4 m s−1) from the north–northeast (average wind direction 32). The MBL was well mixed with a top at 1200 ± 150 m a.g.l. (Fig. S2 in the Supplement). The flight track was carefully planned to target a specific feedlot. The aircraft performed spiral transects at distances ranging from 4–14 km downwind of the feedlot, and stacked horizontal transects were conducted 12 km downwind. The horizontal transects were extended to account for high NH3 and CH4 values observed during the flight. However, it was found that this extended area included more than one feedlot, as determined by NOAA HYSPLIT back trajectories performed on the two peaks of the plume (Fig. S1). The analysis showed that the entire width of the plume contained emissions from multiple feedlots. The combined maximum capacity of all the feedlots within the plume was estimated to be 173 800 maximum heads based on the 2017 data from the Colorado Department of Public Health and Environment (Colorado Department of Natural Resources Oil & Gas Conservation Commission2016). The sampling region during F2 is surrounded by many ONG wells and is located 16 km southeast of Greeley, CO, an urban area with a population > 100 000.

We used a mass balance approach to calculate emissions. Aircraft mass balance has been used to quantify emissions from a variety of source types (examples include Cambaliza et al.2014; Caulton et al.2014; Karion et al.2015; Peischl et al.2015, 2018). Briefly, the mass flow rate of a species through a crosswind plane downwind of the source is approximated by the integration of enhancement above a background concentration over the width and height of the plume. The emissions are derived using Eq. (1) shown here:

(1) M ( u ) = 0 Z MBL - x x ( C u - C b ) × U × ρ ( z ) d x d z .

In Eq. (1), M represents the molar flux (mol s−1) of a gas downwind of the source. To find the enhancement, the local background concentration, Cb (ppbv), is subtracted from the measured concentration, Cu. The ideal gas law is used to calculate the air density (ρ) at every data point, using the universal gas constant at  8.31 J mol−1 K−1. The values ±x are the horizontal limits of the plume width from the center point, and ZMBL is the top of the MBL. The full plume's width was identified by consistent NH3 values above 5 ppbv. The plume's height is constrained to be from ground level to the top of the MBL. A vertical profile near the source was used to identify the MBL level using H2O vapor and calculated potential temperature (Fig. S2). We used a constant value to represent the regional background. The background region was determined from the edges of the horizontal transects for the observed CAFO and used to calculate average values for CH4, NH3, and C2H6.

We transformed all observations surrounding the target CAFO sampled during F2 (Fig. 2) onto a polar coordinate system (r,θ) using the center of the target CAFO as the origin, following the process described in Nathan et al. (2015). The location of the data point on the polar coordinate system (θ) is perpendicular to the theoretical flux surface. Thus, U (m s−1 ) is the corrected perpendicular wind found by taking the cosine of the location of the data point on the polar coordinate system (θ), subtracted by wind direction (ϕ) and multiplied by the wind magnitude (V):

(2) U = cos ( θ - ϕ ) × V .

2.4 Uncertainty analysis

We conducted an uncertainty analysis using a Monte Carlo approach, creating a pseudo-distribution of the data using observed means and standard deviations and recalculating the emission rates. Due to the nature of the plume, not all variables may be represented using a pseudo-distribution. We define the final emission uncertainty as the magnitude of the change in the emission from a combined function of the pseudo-distributions of the following five parameters: background value, perpendicular wind speed, density, MBL depth, and instrumental uncertainty. The Monte Carlo recalculations were first done for individual parameters, then as a combined pseudo-distribution of perpendicular winds, density, background values, and MBL depth to calculate the final uncertainty. The uncertainty for the attribution methods followed the same approach but is addressed separately.

For the uncertainty analysis, we used a Gaussian distribution based on the mean and standard deviation of the original background value to select a new background value randomly. We found the density to have a pattern of increased values to one side of the horizontal transects. Therefore, to account for the pattern, we created a Gaussian distribution using the standard deviation found in each transect and added that onto a moving mean of 1 min. The pseudo-distribution of the MBL height was formed using a uniform randomly selected value of 2600 ± 150 m a.m.s.l. Throughout the day, the MBL height changed by an average of  150 m. The perpendicular winds consist of two variables: wind speed and wind direction. For each separate transect, we created a pseudo-Gaussian distribution of wind speed and wind direction based on the transects' mean and standard deviations for each variable. These Gaussian distributions of wind speed and direction were recalculated through Eq. (2) to produce the new perpendicular wind speeds. The emission recalculations were done 1000 times using all 4 parameters for CH4, C2H6, and NH3. The uncertainty of final emission estimates were calculated by using the 95 % confidence interval (CI) from the Monte Carlo approach divided by the average emission.

When considering the effect of instrumental uncertainty, we only examine the factors that actually impact the measurements used in the mass balance equation: the delta value (Cu-Cb). There are two factors that impact this value: (1) the intercept (bias) and (2) the slope (calibration factor) of the instrument. The bias, which can be affected by drift in the zero reading, is quantified for NH3 and C2H6. To further clarify why these are the only instrumental uncertainty factors we evaluate and why total uncertainty does not affect our results, we consider the following hypothetical situation: two otherwise identical instruments are used to quantify the same enhancement. One instrument reads a background value of 1×106± 0.1 and peak value of 1.000001×106. The second reads a background value of 1±0.1 and an enhancement of 2. The resulting delta values have the same absolute value and uncertainty because of how error propagates through subtraction regardless of which instrument was accurate. We assume that the variability in the bias is dwarfed by real variability in the background, or if the bias is actually large, it similarly affects the background reading. Thus, we expect our analysis of background variability to be the appropriate metric to account for bias.

The uncertainty of the calibration factor is also possible to analyze and include. As far as we are aware, this source of uncertainty does not appear to be routinely reported or incorporated into uncertainty analysis in the mass balance literature. We assumed that the applied calibration factor can vary randomly and applied a randomly picked factor to the delta value to create a pseudo-distribution of possible enhancements. Note that the accuracy of the sensor does not matter. If the uncertainty of the delta value is calculated in the fashion described, there would be equivalent uncertainty of an accurate or inaccurate reading provided the variability of the background and calibration factor was the same. The variability of the calibration factor was <0.1 % for CH4 and 9 % for C2H6. For NH3, Pollack et al. (2019) found the variability to be 2 %.

The effects of instrumental precision, which affect both background values and enhancements, are neglected because we average a large area to calculate the background and bin average 5 s of data for the enhancements. Averaging has the effect of decreasing the random error. The 1 Hz in-flight precision is already low at 1 ppbv for CH4 (Picarro2023), 200 pptv for C2H6, and 60 pptv for NH3 (Pollack et al.2019). Comparatively, the observed variability in the background values dwarfs the error from precision at 33 ppbv, 3 ppbv, and 4 pbbv for CH4, C2H6, and NH3, respectively. Thus, we would expect the error introduced from precision to be negligible.

2.5 CH4 attribution

The two methods described here were applied to either all the data within the MBL during the flight (abbreviated F2) or only the downwind transects used for emission calculation (abbreviated Transect). This provides a total of four scenarios that will be analyzed. Ratios are reported in percentage (ratio ×100 %) across all methods.

2.5.1 Subtraction method

The first attribution approach, referred to as the subtraction method (SM), removes the CH4 emissions related to ONG and attributes the excess CH4 emission to the CAFOs. Rather than using inventory estimates for interfering sources, we calculated observation-based C2H6 : CH4 ratios. The F2 ratio (11 ± 0.02 %) is from Pollack et al. (2022). Briefly, NH3 values >5 ppbv were used to screen out data points associated with CAFOs. The remaining data points were assumed to be associated with ONG sources, which is consistent with observations in this region that CAFOs and ONG are the dominant sources (Kille et al.2019). The slope was calculated using least squares orthogonal distance regressions (ODRs) and is equivalent to the C2H6 : CH4 enhancement ratio. Such a ratio can also be calculated from the CAFO transects as there was a region of CH4 and C2H6 signals that did not include NH3 interference located on the northwest end of the transect. The resulting ratio for the Transect data following the identical calculation is 14.7 ± 0.6 %. The ONG CH4 emissions were removed through Eq. (3):

(3) M CH 4 ag = M CH 4 - CH 4 : C 2 H 6 × M C 2 H 6 .

The inverse of the C2H6 : CH4 ratio is used with the emission rate of C2H6 (mol s−1) to retrieve the portion of CH4 associated with ONG; that value is subtracted from the total CH4 emission rate to result in the CH4 emission associated with the CAFOs. Finally, the molar emission rate (mol s−1) is multiplied by the molar mass to return the emission estimate in grams per second (g s−1).

2.5.2 Multivariate regression

An alternative to the SM is to directly calculate the CAFO CH4 emissions from an NH3 : CH4 ratio and the NH3 emission rate using MVR. Because of the widespread ONG activity in this region, it is not always possible to have clear regions with which to calculate a ratio using a traditional regression approach (as in the SM). Indeed, there is also concern that the ratios calculated in the SM may not be accurate due to the influence of diffuse CAFO or agricultural signals. Figure 3 presents C2H6 and NH3 mixing ratios plotted against CH4. In both plots, there are elevated regions of the other species (i.e., regions of elevated C2H6 in the NH3 versus CH4 plot). Instead, MVR using CH4, C2H6, and NH3 data can be used to calculate NH3 : CH4 and C2H6 : CH4 ratios, as described in Pollack et al. (2022). Briefly, C2H6 and NH3 are assumed to be independent tracers (associated with ONG and CAFOs, respectively), and CH4 is the dependent variable, as shown in Eq. (4).

(4) CH 4 = a + b × NH 3 + c × C 2 H 6

In this equation, a is the background CH4 mixing ratio, b is the inverse effective NH3 : CH4 ratio, and c is the inverse SM C2H6 : CH4 ratio. Unlike Kille et al. (2019) and Pollack et al. (2022), we did not subtract a background mixing ratio from the observed CH4, NH3, or C2H6 mixing ratios, thus the a variable actually represents the local background, and we can compare its value to the observed background that we identified at the edges of the transect. Kille et al. (2019) performed sensitivity analyses on their MVR results. Following the guidance of Kille et al. (2019), we only use the MVR results when all three variables are positive, R2>0.5, and all variables are statistically different from 0. We also tested scenarios with background-subtracted CH4, NH3, and C2H6 to compare with the approaches used in Kille et al. (2019) and Pollack et al. (2022). We performed this analysis on the entire F2 dataset versus the Transect only data. The results of these sensitivity analyses are reported in Tables S1–S2 in the Supplement. The choice of background made no difference to the ratios and only affected value a. This is consistent with Kille et al. (2019) and Pollack et al. (2022). Generally, scenarios where NH3 and C2H6 were not background-subtracted produced a values more consistent with the observed CH4 background and were used for the remainder of the analysis. Slight differences in the MVR results were observed from Pollack et al. (2022). This is attributed to the differences in the area used for MVR between these studies. Pollack et al. (2022) isolated specific source regions for MVR analyses while this study uses all of the data in the study region.

Figure 3Scatter plots of (a–b) C2H6 versus CH4 colored by the mixing ratio of NH3 and (c–d) NH3 versus CH4 colored by the C2H6 mixing ratio sampled during the horizontal transects on 13 November 2019.


3 Results and discussions

3.1 CAFO emissions

Figure 4 shows curtain plots of the measured CH4, NH3, and C2H6 mixing ratios from the horizontal transects. The curtain plots display the data from the flight paths in color-filled boxes representing the area used for the calculation in Eq. (1). In general, elevated NH3 coincides with elevated CH4, indicating the presence of an agricultural plume. However, a high C2H6 enhancement appears to be embedded in the agricultural plume, verifying mixed sources of CH4 in this region. To isolate the agricultural plume and minimize the influence of other sources of CH4, we created a mask around the NH3 signal and screened out points that were <5 ppbv of NH3 (Fig. 4e). Further, the plume was limited to regions with contiguous NH3 values above 5 ppbv. This removed some isolated areas of NH3 enhancement at the lowest altitude. Because the lower transects showed a consistent drop in NH3, we extrapolated that boundary upward to the highest transect as indicated by the vertical dashed line on the left side of Fig. 4. This removes some signal that is primarily downwind of Greeley, CO. The horizontal dotted line represents the boundary layer, and points above this layer were excluded. This threshold was identified by Pollack et al. (2022) and was calculated as 2 times the maximum NH3 mixing ratio observed outside of the plumes of individual CAFOs and within the MBL. The resulting mask was used to limit the area integrated for emission estimates according to Eq. (1). The calculated emission rates are 3330, 609, and 542 g s−1 for CH4, NH3, and C2H6, respectively.

Figure 4Vertical curtain plots of horizontal transects for (a) CH4, (b) NH3, and (c) C2H6 12 km downwind of the observed CAFO. For reference, the surface is located at  1400 m.s.l. Vertical curtain plots of (d) CH4, (e) NH3, and (f) C2H6, where NH3>5 ppbv downwind of the observed CAFO. Points above the boundary layer were removed for calculation. Additionally, points on the left side of the vertical dashed line, above 5 NH3 ppbv, were removed due to their likely origin from Greeley, CO.


We present an in-depth look at the uncertainty contributed by each parameter for CH4 emissions only. Perpendicular winds, density, and MBL depth are consistent in the calculations of NH3 and C2H6 emissions; thus, we only discuss total uncertainty for these. Total uncertainty distribution plots are shown in Fig. S3.

Variability in density had little impact (<1 %) on the final CH4 emission rate. The winds during the plume had a standard deviation of ±1.3, ±1.9, and ±1.8 m s−1 for lowest to highest transect altitude, respectively, with small changes in the wind direction with increased altitude. The wind speed was recalculated into the perpendicular wind and showed minimal change to the final CH4 emissions by  4 %. The location of the background in this study had a relatively consistent CH4 mixing ratio (mean = 1933 ppbv, standard deviation = 1 ppbv). The background is similar to the regional background (1990 ppbv) identified by Pollack et al. (2022). Background variation affects the final CH4 emissions by 5 %. NH3 had higher variations in background values leading to increased uncertainty overall. MBL height was the largest driver of uncertainty and was associated with an 8 % change in CH4 emissions. The significant uncertainty due to MBL depth is expected as changes in MBL depth would result in the interpolation of a different area without a response in concentration measurements. Thus, if the MBL depth were to increase, there would likely be a corresponding decrease in concentrations. The absence of data between the lowest altitude transects and the surface may affect the accuracy of the results, but the associated uncertainty cannot be quantified because we do not have data at the surface. Instrumental uncertainty had little impact, changing the emission rates by <1 % for all: CH4 (0.3 %), NH3 (0.6 %), and C2H6 (0.5 %).

Cambaliza et al. (2014) carried out a detailed analysis on uncertainty from aircraft mass balance calculations. They found that MBL depth and background mixing ratios may have uncertainties up to 19 % and 21 %, respectively, on the overall emission rate estimates. Other studies have confirmed that uncertainty associated with the MBL depth makes a large contribution to overall uncertainty for these types of calculations (Karion et al.2013; Peischl et al.2015). The winds in the mass balance are assumed to be consistent from the release of the emissions to the location of the measurement. Therefore, the natural variability of winds may contribute to the uncertainty more with less consistent winds (including direction and speed) and may have a large effect on the final uncertainty similar to MBL and background value (Karion et al.2013; Peischl et al.2015). During F2, winds were consistent and made a smaller contribution to the overall uncertainty on this particular day than noted by these other studies.

The total uncertainty for CH4 emissions calculated from horizontal transects was ±10 %. For the C2H6 and NH3 emission estimates, the total uncertainty was ±14 % and ±17 %, respectively. Other factors that may influence the accuracy of the emissions include smaller-scale variations in the mixing ratios and regions of the plume that the flights did not sample. As it is not possible to sample the entire MBL from top to bottom, the vertical spacing between the horizontal transects may result in errors. Errors associated with interpolation were not explored here but are expected to be small (Cambaliza et al.2014). Further, times when the bulk of the mixing ratio is near the ground are seen to have the most uncertainty (Gordon et al.2015). The elevated mixing ratio values of CH4, NH3, and C2H6 were seen throughout the MBL. The uncertainty estimates are specific to this flight. Different meteorological conditions can lead to varying uncertainty estimates. These uncertainties are specific to each flight and can be influenced by factors such as smaller-scale variations in mixing ratios, unexplored errors associated with interpolation, and the proximity of the bulk mixing ratio to the ground. In particular, MBL uncertainty or growth may significantly affect the magnitude of the uncertainty. The uncertainty analysis methodology presented here can be applied to any future mass balance emission estimates.

3.2 CAFO CH4 attribution: comparison of methods

Results from the SM and MVR CAFO CH4 attribution and associated uncertainty are presented in Table 1. Reported uncertainties in Table 1 represent results from combining the ratio uncertainty to the rest of the pseudo-distributions as described in Sect. 2.4. Both the magnitude of the attributed CAFO CH4 emission and its uncertainty vary between approaches. These methods represent typical approaches that can be undertaken to isolate the emissions for a given facility. In order to identify the optimal method, we identify four criteria: (1) the relationships should predict the total CH4 well; in other words, traditional goodness-of-fit (GOF) values should be optimal (the average residual of the fit should be near 0 and R2 values should be high); (2) the relationships for NH3 : CH4 and C2H6 : CH4 should be consistent with observations; (3) the method should err on the side of being conservative, meaning it should be more likely to under-attribute than over-attribute the CAFO emissions; and (4) uncertainty of the result should be low. As currently presented, these criteria are qualitative because each analysis is unique. We will discuss the implications and further refinements of these criteria later on. These criteria are presented in order of their importance to the final recommendation. This approach ensures that low uncertainty is not the primary deciding factor.

Table 1CH4 Attribution sensitivity.

a SE = standard error. b Pollack et al. (2022). c Criteria are defined in Sect. 3.2.

Download Print Version | Download XLSX

To investigate criterion 1, we calculated total fits, residuals, and R2 values for all four methods and assessed how well the NH3 : CH4 and C2H6 : CH4 ratios represented the Transect data. The central assumption is that regardless of the attribution method, all methods should predict the observed CH4 values well as evaluated by traditional metrics of goodness of fit (e.g., R2, residuals). MVR directly produces both ratios such that predicted CH4 can be calculated from the NH3 and C2H6 time series. To construct a CH4 prediction for the SM, which does not require calculation of the NH3 : CH4 ratios, we inferred the effective NH3 : CH4 ratio from the ratio of the calculated NH3 emissions to the CH4 attributed to the CAFO by the SM. The resulting ratio is then used in Eq. (4) to produce an effective MVR prediction. Table S3 shows the full range of variables for each scenario. GOF statistics are calculated for the Transect data only.

The results of this comparison are shown in Fig. 5, with the prediction limited to the Transect data, because these are the data that are actually used to calculate the CH4 emission rate. The Transect MVR approach was the method with the best GOF statistics including an average residual closest to 0 and the highest R2 value; it is also the only approach that passes criterion 1. In some ways, this is unsurprising because fits are calculated by minimizing the sum of the residual squared, and this is the only method that directly fits the Transect data (Skoog et al.2004). Still, this analysis does provide the means to compare the other methods and can be useful when site-specific MVR is not possible. The background values were consistent in all methods and consistent with observations. The three other methods all overpredicted CH4 and had average negative residuals and lower R2 values. The overprediction suggests that these methods produce an incorrect relationship between CH4 and one (or both) of the tracer gases, but it is not possible to identify which relationship is incorrect. The method that best simulates the data after the Transect MVR is the F2 MVR.

Figure 5(a) Residual CH4 (observation–prediction) versus predicted CH4 (ppbv). The ratios corresponding to the transect MVR results are shown in red, the F2 MVR results are shown in green, the transect SM results are shown in blue, and the F2 SM results are shown in cyan. The black line in both plots represents the zero line. (b) Violin plots of residuals for the four different methods: the black lines are the calculated mean residuals, the red lines are the median residuals, and the closed black dots are the sum of the residuals squared (normalized by the minimum value ×10). The violin plots are ordered by increasing absolute average residual.


The effective NH3 : CH4 ratios produced for this method can also be compared directly in support of criterion 2. Criterion 2 is difficult to evaluate because the actual ratios would ideally be known. We can compare these to literature values, compare values between scenarios, and look at the underlying assumptions of the calculations to evaluate this criterion. The C2H6 : CH4 ratios calculated for the SM and MVR varied by 40 %, which substantially affects how much CH4 is attributed to ONG. However, the range of ONG ratios is broadly consistent with other observations in this region (Kille et al.2017, 2019; Yacovitch et al.2017; Peischl et al.2018).

As a plume moves downwind from its source, it is expected that the plume will disperse over the MBL. NH3 has a short lifetime ( hours) and may be removed through dry deposition and transition into the particle phase, even at relatively short distances downwind from sources (Staebler et al.2009; Miller et al.2015). Thus, there may be lower concentrations downwind of gaseous NH3 than expected from dilution alone. The range of reported NH3 : CH4 ratios in the literature is quite variable and includes observations from 1 %–200 % (Eilerman et al.2016; Golston et al.2020). We used Pollack et al. (2022) as our primary reference for comparison. Comparison to this study revealed that the transects located  12 km downwind of the observed target CAFO do not have the same NH3 : CH4 ratio as transects collected closer to the CAFO. There is also a wide variety of ratios reported in Pollack et al. (2022) ranging from 80 % to 270 % across different facilities which illustrate why the MVR results from the full F2 data may not be appropriate to represent an individual site. In general, the SM produced ratios that were lower than MVR and observations near the source (Pollack et al.2022). The SM Transect ratio is particularly low (45 %), about half the value of any other ratio produced from the other attribution methods or results from Pollack et al. (2022). This indicates that this approach actually fails criterion 2.

On a theoretical basis, the MVR ratios calculated from the F2 data include locations near other CAFOs and locations near ONG sources. Therefore, the F2 MVR results may not accurately represent the area around this particular plume. The SM using F2 data has similar challenges. The SM assumes that C2H6 : CH4 is constant throughout the area and that there are no other sources of CH4 other than the CAFO and ONG. However, we know there may be other trace amounts of CH4 from vehicles, waste areas, and other small sources. Generally, we would expect data nearest the source to produce the most accurate ratios and pass criterion 2. This would generally exclude the F2 results and for this reason, we assume that these results fail criterion 2.

For criterion 3, the two methods are different in approach and it is quite easy to identify which approach is more conservative. The SM attributes anything outside of ONG to the CAFOs plume; thus, there is no unattributed CH4 for these methods. The SM is, therefore, expected to be less conservative and an upper limit of the CAFOs CH4 emissions due to the possibility of other sources in the region.

In the MVR approach, the R2 value of the fit gives an estimate how much of the CH4 signal is explained by the chosen tracers. The MVR fits had R2 values of 0.72 for F2 and 0.74 for Transects, which leaves an amount of CH4 that is not well correlated to NH3 or C2H6. In this scenario, this excess CH4 is left unattributed in the MVR approach (56 % for F2 and 24 % for Transect). A visual representation of the unattributed portion of the CH4 signal for the transect MVR is presented in Fig. S4. The unattributed CH4 is broadly distributed and cannot be attributed to any particular source. The MVR approach avoids over-attributing the CH4 signal to the CAFOs.

In general, we expect the Transect MVR approach, which uses the data closest to the source, to provide the best estimate of these CAFO emissions in accordance with the first three criteria. This approach produces local ratios on par with the literature, has the closest average residual to 0, and is conservative in its approach to attribution. Transect MVR has slightly greater uncertainty than F2 MVR and both have lower uncertainty than SM because they do not require subtraction. The impact of subtraction on uncertainty is evident in Table 1 and is the basis for including criterion 4. While the ONG ratios calculated either by SM or MVR using the F2 data differ by  40 %, the uncertainty of the results using MVR attribution is a factor of 3 lower. This is despite the fact that the SM ratio uncertainty is actually a factor of 5 lower than the MVR ratio. If we had used a comparable error on the C2H6 : CH4 ratio to MVR with SM, the error would approach 100 %. The primary reason for the difference in uncertainty between the SM and MVR approaches relates to how error is propagated through subtraction versus multiplication or division. For subtraction, absolute errors add in quadrature, while for multiplication or division, relative errors add in quadrature (Skoog et al.2004). The net effect is that error on the SM will be very high when the absolute amount to be subtracted is large. The ONG ratio is not directly used for attribution in the MVR approach, which is why the CAFO-attributed emissions by the F2 MVR and SM approaches vary by  30 %. However, the relative uncertainty (standard error or ratio) of the NH3 : CH4 MVR ratio is comparable to the relative uncertainty of the MVR ONG ratio.

The subtraction method is very similar to inventory subtraction in principle. Error estimates are often not provided for inventory data, but the error introduced to emission estimates by subtracting inventories may be substantial. Zimmerle et al. (2022) recently showed how ONG sources could vary and produce uncertainties in excess of estimates relying on traditional inventory data. The MVR approach is generally attractive to avoid such errors when interfering signals are high and there are enough data to produce robust MVR results. The plume from the observed CAFOs was broad, allowing for multiple data points within the plume, and produced results with low uncertainty (20 %). For sites with narrower and lower enhancements, local MVR may not be possible. Our criteria would suggest that using MVR on the full flight data may be an appropriate proxy when local ratios with MVR cannot be calculated. However, the results from the MVR Transect and F2 data are statistically different in this case. Additional analysis may need to be done to ensure that the MVR results from larger datasets are applicable to the site of interest. Repeat measurements or ground observations would help to constrain such results.

3.3 Additional case studies

To test the methodology we have developed so far for methane attribution in complicated regions more thoroughly, we have applied our strategy to a case study from a larger dataset of observations in the region. These data were collected as part of the Transport and Transformation of Ammonia (Trans2Am) campaign during the summer of 2021 and 2022, where updated flight plans were performed for CAFOs in the same region as this work. While full results from that campaign are outside the scope of this work, we compare our results to one research flight (research flight 13, denoted as RF13 in the remaining text) in the same vicinity as the facilities measured during the original F2. The updated flight plans, shown in Fig. 6, included multiple downwind transects to allow computation of emissions at consecutive distances over the course of  1.5 h. We calculated CH4 and NH3 emission rates and uncertainties identically to the previously described methodology. We would expect the CH4 emission rate to be conserved at different downwind distances. The NH3 emission may be affected by deposition to the surface and thus may remain the same or decrease downwind from a facility.

Figure 6Path of Trans2Am RF13, colored by CH4 (ppbv). Black numbers represent the corresponding transect number. Animal operations are indicated by the different colored circles as in Fig. 1. Green dots represent ONG wells; data of ONG are as of 2015 (Colorado Department of Natural Resources Oil & Gas Conservation Commission2016).

Table 2 shows results from RF13, which was collected on 23 August 2021 near Greeley, CO. The target CAFO was the same target feedlot in F2. The top of the MBL was located at 2000 ± 200 m a.g.l. An average temperature of 24 C and southwesterly winds with an average wind speed of 4.1 m s−1 were observed during the flight, and a total of five downwind transects were performed. Background values and uncertainties for CH4 were obtained from the MVR fit and varied for each transect, ranging from 1886 to 1927 ppbv. The background NH3 value (3.023 ppbv) was calculated by constraining observations to those associated with the target CAFO (Fig. 6) and identifying the minimum NH3 concentration. Enhanced concentrations of CH4 and NH3 were observed during all transects, spanning about 17 km downwind of the target facility. CAFO CH4 was attributed using the transect MVR approach for all five transects described previously.

Table 2CH4 Attribution for Trans2Am RF13.

* SE = standard error.

Download Print Version | Download XLSX

Attributed CH4 ranged between 189 ± 55 to 533 ±162 g s−1 with a total average of 365 ± 89 g s−1. The average emission rates for this site convert to 12 ± 3 g of CH4 per head per hour and 13 ± 3 g of NH3 per head per hour. There is a decreasing trend in the NH3 : CH4 ratio observed as the transect distance increases, indicating that NH3 depletion is occurring, likely by dry deposition. The per head values use the combined maximum number of cattle from all facilities within the plume (109 500 cattle), according to the CAFO Permit Database 2021 (CDPHE2022). The total average uncertainty range (95 % CI) for this site was calculated from the distribution of the sum of all five transect Monte Carlo simulations. The total average relative uncertainty for this case (24 %) is higher than the uncertainty for the previous case (20 %) due to different environmental conditions. While there is considerable variability in the CH4 emissions from this site, we note that the NH3 emissions show similar variability. There is a factor of 2.7 between the minimum and maximum CH4 emission observations and a factor of 2.4 between the minimum and maximum NH3 emission observations. Since the NH3 emission rate is not affected by any attribution method, it is assumed that the variability in CH4 is primarily caused by on-site or environmental conditions and not the MVR attribution method. Comparisons between the MVR and SM attribution methods for this facility show (Table S4) that the SM method produces more variable results (factor of 3.9 min–max range) and overpredicts the total average methane by about 49 %, while the total average relative uncertainty associated with the SM method for this site is almost double the uncertainty associated with the MVR method. These results are in agreement with our previous results. The average per head CH4 and NH3 results for this site are not statistically different from the November 2019 results despite the different wind direction, background values, etc.

3.4 Comparison to literature

A few recent studies have reported CH4 emission rates from CAFOs. The emission factors calculated from the Transect MVR method in this work are 13 ± 3 g of CH4 per head per hour and 13 ± 2 g of NH3 per head per hour in November 2019 (F2) and 12 ± 3 g of CH4 per head per hour, and 13 ± 3 g of NH3 per head per hour August 2021 (RF13). The calculation of emission rates as a value per head is based on the maximum amount of cattle allowed in the CAFOs, although the actual number of cattle in each CAFO at the time of the measurement is unknown. Our maximum number of cattle identified for the F2 data consisted of mostly beef cattle at 161 500 heads, and only 12 300 heads were dairy cattle (total 173 800 heads). The Environmental Protection Agency (EPA) estimated CH4 emissions for beef cattle is 7.2 g per head per hour and previous studies range from 7–9 g per head per hour (EPA Annex2021; Golston et al.2020; Hacker et al.2016). Dairy cattle have high rates of CH4 emissions due to higher exertion on the animal with rates ranging from 14–39 g per head per hour and the EPA national average is 18 g per head per hour (EPA Annex2021; Leytem et al.2011; Bjorneberg et al.2009; Griffith et al.2008). Golston et al. (2020) conducted a ground-based study in the NCFR and found the highest emission rates in dairy cattle at 39.32 ± 3.92 g of CH4 per head per hour. In contrast, their reported emissions from beef cattle were much lower and similar to other studies (9.48 ± 0.93 g of CH4 per head per hour). However, they noted that there was a significant difference in emission rates of repeat observation for an individual CAFO.

The difference between the CH4 emissions per head of cattle in this work compared to prior works could be due to the management of the CAFOs. Enteric fermentation can be altered through the composition of the food given to the animals. For instance, the type and maturity of diets provided to the animals may modify the nutrients and digestibility of the food (Archimède et al.2011). CH4 emissions from cattle are also known to depend on exertion on the animals including exercise and stress (EPA2022). In this study, the food source and feeding schedule were unknown. We also note that one CAFO was extremely large with 98 000–100 000 reported maximum cattle heads during the time of study and one of the largest in Colorado. Most emission estimates are based on considerably smaller CAFOs.

The NH3 emissions per head reported here are also slightly higher than previous estimates, which ranged from 2–12 g per head per hour (Staebler et al.2009; Golston et al.2020; Hacker et al.2016). The NH3 emission per head calculation in our work does not require removal of interfering signals like the CH4 emissions. The NH3 emissions are highly affected by meteorological conditions including seasonality and time of day. Waste that is stored outside may exacerbate NH3 emissions in warmer temperatures (Montes et al.2013). There is also evidence of a diurnal pattern of NH3 emissions from CAFOs with a peak during noon local time, near the time of sampling in this work (Shonkwiler and Ham2017). NH3 : CH4 ratios follow a similar pattern with higher ratios during the midday and lower ratios at night (Eilerman et al.2016).

4 Conclusions

We demonstrate an approach to isolate and quantify CH4 plumes of agricultural CAFO sources with interfering sources by using the SM and MVR. The SM uses NH3 as a tracer to identify a CAFO plume before using a C2H6 : CH4 ratio to remove any ONG CH4 emissions interfering with the CAFO source. It is an optimal method when there is little ONG influence to keep the error introduced from subtraction small. The MVR method uses C2H6 and NH3 as tracers and provides lower uncertainty of emission estimates. This approach is appropriate when tracer data are available and there is enough signal to produce statistically significant relationships and thus site-specific ratios. The criteria to identify the best approach may be useful for isolating and attributing emissions from specific sources in other regions. Overall, our best estimates of emissions from the observed CAFOs are 13 (±2) g per head per hour for CH4 and 13 (±2) g per head per hour for NH3. These findings are significantly higher than EPA inventory estimates and previous studies highlighting the need for more observations and estimates of CAFO CH4 and NH3 emission rates.

Data availability

The merged data file for the F2 data are available from WyoScholar (Caulton2023). This long-term repository is open-access and provides freely available data stewarded by the University of Wyoming Libraries. The RF13 data collected as part of the Trans2Am campaign are publicly available via the National Center for Atmospheric Research Earth Observing Laboratory.


The supplement related to this article is available online at:

Author contributions

MEM: writing – original draft, formal analysis, visualization, data acquisition and curation. IBP: writing – reviewing and editing, data acquisition and curation, funding acquisition. EVF: writing – reviewing and editing, funding acquisition. KMS: formal analysis, data acquisition and curation, writing – reviewing and editing. DRC: conceptualization, supervision, writing – reviewing and editing, funding acquisition, data collection.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to thank the University of Wyoming King Air facility for their support of this project. We would also like to thank Kristen Pozsonyi for help with instrument installation and data collection.

Financial support

This research has been supported by the National Science Foundation (grant no. 2020151).

Review statement

This paper was edited by Eleanor Browne and reviewed by two anonymous referees.


Archimède, H., Eugène, M., Marie Magdeleine, C., Boval, M., Martin, C., Morgavi, D. P., Lecomte, P., and Doreau, M.: Comparison of methane production between C3 and C4 grasses and legumes, Anim. Feed Sci. Tech., 166–167, 59–64,, 2011. a

Arnold, S., Dileo, J., and Takushi, T.: Colorado Greenhouse Gas Inventory — 2014 Update Including Projections to 2020 & 2030, Department of Public Heath & Environment, (last access: 26 April 2021), 2014. a

Bjorneberg, D. L., Leytem, A. B., Westermann, D. T., Griffiths, P. R., Shao, L., and Pollard, M. J.: Measurement of Atmospheric Ammonia, Methane, and Nitrous Oxide at a Concentrated Dairy Production Facility in Southern Idaho Using Open-Path FTIR Spectrometry, T. ASABE, 52, 1749–1756, 2009. a

Cambaliza, M. O. L., Shepson, P. B., Caulton, D. R., Stirm, B., Samarov, D., Gurney, K. R., Turnbull, J., Davis, K. J., Possolo, A., Karion, A., Sweeney, C., Moser, B., Hendricks, A., Lauvaux, T., Mays, K., Whetstone, J., Huang, J., Razlivanov, I., Miles, N. L., and Richardson, S. J.: Assessment of uncertainties of an aircraft-based mass balance approach for quantifying urban greenhouse gas emissions, Atmos. Chem. Phys., 14, 9029–9050,, 2014. a, b, c

Caulton, D.: Ammonia phase PARtitioning and Transport-lite (APART-lite), University of Wyoming [data set],, 2023. a

Caulton, D. R., Shepson, P. B., Santoro, R. L., Sparks, J. P., Howarth, R. W., Ingraffea, A. R., Cambaliza, M. O., Sweeney, C., Karion, A., Davis, K. J., Stirm, B. H., Montzka, S. A., and Miller, B. R.: Toward a better understanding and quantification of methane emissions from shale gas development, P. Natl. Acad. Sci. USA, 111, 6237–6242,, 2014. a, b

CDPHE: Colorado Department of Public Health and Environment CAFO Permit Database, (last access: 27 January 2021), 2017. a, b

CDPHE: Colorado Department of Public Health and Environment CAFO Permit Database 2021, (last access: March 2022), 2022. a

Colorado Department of Natural Resources Oil & Gas Conservation Commission: Colorado Oil and Gas Information System (COGIS), (last access: April 2016), 2016. a, b, c

Crosson, E. R.: A cavity ring-down analyzer for measuring atmospheric levels of methane, carbon dioxide, and water vapor, Appl. Phys. B-Lasers O., 92, 403–408,, 2008. a

Eilerman, S. J., Peischl, J., Neuman, J. A., Ryerson, T. B., Aikin, K. C., Holloway, M. W., Zondlo, M. A., Golston, L. M., Pan, D., Floerchinger, C., and Herndon, S.: Characterization of Ammonia, Methane, and Nitrous Oxide Emissions from Concentrated Animal Feeding Operations in Northeastern Colorado, Environ. Sci. Technol., 50, 10885–10893,, 2016. a, b, c

EPA: Inventory of U.S. Greenhouse Gas Emissions and Sinks: 1990–2020, U.S. Environmental Protection Agency, Epa 430-R-18-003, 1–205, (last access: 9 May 2022), 2022. a, b, c, d

EPA Annex: U.S. Green House Gas Emissions and Sinks 1990–2019 Annex 3 Part B, 4, 1–71, (last access: 29 February 2021), 2021. a, b

Golston, L. M., Pan, D., Sun, K., Tao, L., Zondlo, M. A., Eilerman, S. J., Peischl, J., Neuman, J. A., and Floerchinger, C.: Variability of Ammonia and Methane Emissions from Animal Feeding Operations in Northeastern Colorado, Environ. Sci. Technol., 54, 11015–11024,, 2020. a, b, c, d, e, f

Gordon, M., Li, S.-M., Staebler, R., Darlington, A., Hayden, K., O'Brien, J., and Wolde, M.: Determining air pollutant emission rates based on mass balance using airborne measurement data over the Alberta oil sands operations, Atmos. Meas. Tech., 8, 3745–3765,, 2015. a

Griffith, D. W., Bryant, G. R., Hsu, D., and Reisinger, A. R.: Methane emissions from free-ranging cattle: Comparison of tracer and integrated horizontal flux techniques, Aust. J. Exp. Agr., 48, 582–591,, 2008. a

Hacker, J. M., Chen, D., Bai, M., Ewenz, C., Junkermann, W., Lieff, W., McManus, B., Neininger, B., Sun, J., Coates, T., Denmead, T., Flesch, T., McGinn, S., and Hill, J.: Using airborne technology to quantify and apportion emissions of CH4 and NH3 from feedlots, Anim. Prod. Sci., 56, 190–203,, 2016. a, b, c

Higley, B. D. K. and Cox, D. O.: Oil and Gas Exploration and Development along the Front Range in the Denver Basin of Colorado, Nebraska, and Wyoming, USGS Province 39: U.S. Geological Survey Digital Data Series, Ch. 2, (last access: 11 April 2021), 2007. a

Karion, A., Sweeney, C., Pétron, G., Frost, G., Michael Hardesty, R., Kofler, J., Miller, B. R., Newberger, T., Wolter, S., Banta, R., Brewer, A., Dlugokencky, E., Lang, P., Montzka, S. A., Schnell, R., Tans, P., Trainer, M., Zamora, R., and Conley, S.: Methane emissions estimate from airborne measurements over a western United States natural gas field, Geophys. Res. Lett., 40, 4393–4397,, 2013. a, b

Karion, A., Sweeney, C., Kort, E. A., Shepson, P. B., Brewer, A., Cambaliza, M., Conley, S. A., Davis, K., Deng, A., Hardesty, M., Herndon, S. C., Lauvaux, T., Lavoie, T., Lyon, D., Newberger, T., Pétron, G., Rella, C., Smith, M., Wolter, S., Yacovitch, T. I., and Tans, P.: Aircraft-Based Estimate of Total Methane Emissions from the Barnett Shale Region, Environ. Sci. Technol., 49, 8124–8131,, 2015. a

Kille, N., Baidar, S., Handley, P., Ortega, I., Sinreich, R., Cooper, O. R., Hase, F., Hannigan, J. W., Pfister, G., and Volkamer, R.: The CU mobile Solar Occultation Flux instrument: structure functions and emission rates of NH3, NO2 and C2H6, Atmos. Meas. Tech., 10, 373–392,, 2017. a

Kille, N., Chiu, R., Frey, M., Hase, F., Sha, M. K., Blumenstock, T., Hannigan, J. W., Orphal, J., Bon, D., and Volkamer, R.: Separation of Methane Emissions From Agricultural and Natural Gas Sources in the Colorado Front Range, Geophys. Res. Lett., 46, 3990–3998,, 2019. a, b, c, d, e, f, g, h, i

Leytem, A. B., Dungan, R. S., Bjorneberg, D. L., and Koehn, A. C.: Emissions of Ammonia, Methane, Carbon Dioxide, and Nitrous Oxide from Dairy Cattle Housing and Manure Management Systems, J. Environ. Qual., 40, 1383–1394,, 2011. a

Maasakkers, J. D., Jacob, D. J., Sulprizio, M. P., Turner, A. J., Weitz, M., Wirth, T., Hight, C., DeFigueiredo, M., Desai, M., Schmeltz, R., Hockstad, L., Bloom, A. A., Bowman, K. W., Jeong, S., and Fischer, M. L.: Gridded National Inventory of U.S. Methane Emissions, Environ. Sci. Technol., 50, 13123–13133,, 2016. a, b

Mielke-Maday, I., Schwietzke, S., Yacovitch, T. I., Miller, B., Conley, S., Kofler, J., Handley, P., Thorley, E., Herndon, S. C., Hall, B., Dlugokencky, E., Lang, P., Wolter, S., Moglia, E., Crotwell, M., Crotwell, A., Rhodes, M., Kitzis, D., Vaughn, T., Bell, C., Zimmerle, D., Schnell, R., and Pétron, G.: Methane source attribution in a U.S. dry gas basin using spatial patterns of ground and airborne ethane and methane measurements, Elementa, 7, 13,, 2019. a

Miller, D. J., Sun, K., Tao, L., Pan, D., Zondlo, M. A., Nowak, J. B., Liu, Z., Diskin, G., Sachse, G., Beyersdorf, A., Ferrare, R., and Scarino, A. J.: Ammonia and methane dairy emission plumes in the San Joaquin Valley of California from individual feedlot to regional scales, J. Geophys. Res.-Atmos., 120, 9718–9738,, 2015. a

Montes, F., Meinen, R., Dell, C., Rotz, A., Hristov, A. N., Oh, J., Waghorn, G., Gerber, P. J., Henderson, B., Makkar, H. P., and Dijkstra, J.: SPECIAL TOPICS-Mitigation of methane and nitrous oxide emissions from animal operations: II. A review of manure management mitigation options, J. Anim. Sci., 91, 5070–5094,, 2013. a

Moumen, A., Azizi, G., ChekrounK, B., and Baghour, M.: The effects of livestock methane emission on the global warming: A review, Int. J. Global Warming, 9, 229–253,, 2016. a

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., Zhan, H., and Zhang, H.: Anthropogenic and Natural Radiative Forcing: Supplementary Material, Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., 1–44, (last access: 4 September 2022), 2013. a

Nathan, B. J., Golston, L. M., Brien, A. S. O., Ross, K., Harrison, W. A., Tao, L., Lary, D. J., Johnson, D. R., Covington, A. N., Clark, N. N., and Zondlo, M. A.: Near-Field Characterization of Methane Emission Variability from a Compressor Station Using a Model Aircraft, Environ. Sci. Technol., 49, 7896–7903,, 2015. a

Peischl, J., Ryerson, T. B., Aikin, K. C., Gouw, J. A., Gilman, J. B., Holloway, J. S., Lerner, B. M., Nadkarni, R., Neuman, J. A., Nowak, J. B., Trainer, M., Warneke, C., and Parrish, D. D.: Quantifying Methane and Ethane Emissions to the Atmosphere From Central and Western U.S. Oil and Natural Gas Production Regions, J. Geophys. Res., 120, 2119–2139,, 2015. a, b, c, d

Peischl, J., Eilerman, S. J., Neuman, J. A., Aikin, K. C., de Gouw, J., Gilman, J. B., Herndon, S. C., Nadkarni, R., Trainer, M., Warneke, C., and Ryerson, T. B.: Quantifying Methane and Ethane Emissions to the Atmosphere From Central and Western U.S. Oil and Natural Gas Production Regions, J. Geophys. Res.-Atmos., 123, 7725–7740,, 2018. a, b, c

Picarro: G2401 Analyzer Datasheet, (last access: 12 August 2022), 2023. a

Pollack, I. B., Lindaas, J., Roscioli, J. R., Agnese, M., Permar, W., Hu, L., and Fischer, E. V.: Evaluation of ambient ammonia measurements from a research aircraft using a closed-path QC-TILDAS operated with active continuous passivation, Atmos. Meas. Tech., 12, 3717–3742,, 2019. a, b, c

Pollack, I. B., McCabe, M. E., Caulton, D. R., and Fischer, E. V.: Enhancements in Ammonia and Methane from Agricultural Sources in the Northeastern Colorado Front Range Using Observations from a Small Research Aircraft, Environ. Sci. Technol., 56, 2236–2247,, 2022. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q

Richardson, S. J., Miles, N. L., Davis, K. J., Crosson, E. R., Rella, C. W., and Andrews, A. E.: Field testing of cavity ring-down spectroscopy analyzers measuring carbon dioxide and water vapor, J. Atmos. Ocean Tech., 29, 397–406,, 2012. a

Roscioli, J. R., Zahniser, M. S., Nelson, D. D., Herndon, S. C., and Kolb, C. E.: New Approaches to Measuring Sticky Molecules: Improvement of Instrumental Response Times Using Active Passivation, J. Phys. Chem. A, 120, 1347–1357,, 2015. a, b

Shonkwiler, K. B. and Ham, J. M.: Agr. Forest Meteorol. Ammonia emissions from a beef feedlot: Comparison of inverse modeling techniques using long-path and point measurements of fenceline NH3, Agr. Forest Meteorol., 258, 29–42,, 2017. a

Skoog, D. A., West, D. M., Holler, F. J., and Crouch, S. R.: Fundamentals of Analytical Chemistry, 8th Edn., Thomson-Brooks/Cole, ark:/13960/t6vx0xq3j, 2004. a, b

Smith, C., Nicholls, Z., Armour, K., Collins, W., Forster, P., Meinshausen, M., Palmer, M., and Watanabe, M.: The Earth’s Energy Budget, Climate Feedbacks, and Climate Sensitivity Supplementary Material, (last access: 2 November 2021), 2021. a

Smith, M. L., Kort, E. A., Karion, A., Sweeney, C., Herndon, S. C., and Yacovitch, T. I.: Airborne Ethane Observations in the Barnett Shale: Quantification of Ethane Flux and Attribution of Methane Emissions, Environ. Sci. Technol., 49, 8158–8166,, 2015. a

Staebler, R. M., McGinn, S. M., Crenna, B. P., Flesch, T. K., Hayden, K. L., and Li, S. M.: Three-dimensional characterization of the ammonia plume from a beef cattle feedlot, Atmos. Environ., 43, 6091–6099,, 2009. a, b, c

Strauss, L., Serafin, S., Haimov, S., and Grubišić, V.: Turbulence in breaking mountain waves and atmospheric rotors estimated from airborne in situ and Doppler radar measurements, Q. J. Roy. Meteor. Soc., 141, 3207–3225,, 2015. a

Townsend-Small, A. E., Botner, C., Jimenez, K. L., Schroeder, J. R., Blake, N. J., Meinardi, S., Blake, D. R., Sive, B. C., Bon, D., Crawford, J. H., Pfist, G., and Frank, F. M.: Using stable isotopes of hydrogen to quantify biogenic and thermogenic atmospheric methane sources: A case study from the Colorado Front Range Amy, Geophys. Res. Lett., 31, 1–8,, 2016. a, b

United States Department of Agriculture [USDA]: Colorado Annual Bulletin Colorado Agricultural Statistics, National Agricultural Statistics Service 1, 1–58, (last access: 5 January 2021), 2018.  a

Yacovitch, T. I., Herndon, S. C., Roscioli, J. R., Floerchinger, C., McGovern, R. M., Agnese, M., Pétron, G., Kofler, J., Sweeney, C., Karion, A., Conley, S. A., Kort, E. A., Nähle, L., Fischer, M., Hildebrandt, L., Koeth, J., McManus, J. B., Nelson, D. D., Zahniser, M. S., and Kolb, C. E.: Demonstration of an Ethane Spectrometer for Methane Source Identification, Environ. Sci. Technol., 48, 8028–8034,, 2014. a, b

Yacovitch, T. I., Daube, C., Vaughn, T. L., Bell, C. S., Roscioli, J. R., Knighton, W. B., Nelson, D. D., Zimmerle, D., Pétron, G., and Herndon, S. C.: Natural gas facility methane emissions: measurements by tracer flux ratio in two US natural gas producing basins, Elementa, 5, 69,, 2017. a

Zimmerle, D., Duggan, G., Vaughn, T., Bell, C., Lute, C., Bennett, K., Kimura, Y., Cardoso-Saldaña, F. J., and Allen, D. T.: Modeling air emissions from complex facilities at detailed temporal and spatial resolution: The Methane Emission Estimation Tool (MEET), Sci. Total Environ., 824, 153653,, 2022. a, b

Short summary
Agriculture emissions, including those from beef and dairy cattle feeding operations, make up a large portion of the United States’ total greenhouse gas emissions, but many of these operations reside in areas where methane from oil and natural gas is prevalent, making it difficult to attribute methane in these areas. This work investigates two approaches to emission attribution for cattle feeding operations and provides guidance for emission attribution in other complicated regions.
Final-revised paper