attribution in a U.S. dry gas basin using spatial patterns of ground and airborne ethane and methane measurements

,

estimates between different sources of CH 4 . The mass balance method used with aircraft data provides total CH 4 emissions for an area, which then must be apportioned between O&NG and other sources. To do so, previous top-down regional and national studies have often relied completely or partly on bottom-up information for CH 4 emissions apportionment or non-O&NG emissions estimation (Caulton et al., 2014;Karion et al., 2013;Pétron et al., 2012), even in cases where C 2 H 6 was used to distinguish between sources (Peischl et al., 2015. Atmospheric measurements of CH 4 alone, unless directly downwind of a point source Frankenberg et al., 2016), are usually not enough to distinguish and separate contributions from different CH 4 sources. A few apportionment studies in O&NG fields have used alkane ratios or CH 4 isotopic signatures to distinguish between emissions from various CH 4 sources (Hopkins et al., 2016;Smith et al., 2015Smith et al., , 2017Townsend-Small et al., 2015). However, results from these studies can be limited by a relatively small number of discrete samples or days of measurements. Additionally, mass balance and attribution methods using alkane or isotope enhancement ratios require subtracting local background values, which have not always been sufficiently documented .
The work presented here is part of a larger coordinated intensive study conducted in the Fayetteville Shale dry gas play in north-central Arkansas in September and October 2015 to compare state of the science bottom-up and topdown regional NG CH 4 emission estimates. Emissions quantification at the facility or equipment level is presented by , Conley et al. (2017), Robertson et al. (2017), Vaughn et al. (2017), Yacovitch et al. (2017), and Zimmerle et al. (2017). An aircraft mass balance estimate of spatially resolved CH 4 emissions from the study area on two consecutive days is described in Schwietzke et al. (2017), and a reconciliation between this top-down estimate and a bottom-up CH 4 emissions estimates is provided by Vaughn et al. (2018). The bottom-up model presented by Vaughn et al. revealed a strong diurnal variability in total NG CH 4 emissions, with the peak occurring during the mid-afternoon due to manual liquid unloading events at 107 wells over the two days modeled . Modeled bottom-up emissions estimates captured the same spatial variability in CH 4 emissions throughout the study area with a reasonable agreement in emissions estimates between the two approaches for the total study area, eastern half, and western half. NG CH 4 emissions normalized by natural gas production were found to be two times higher in the western half than in the eastern half .
Presented here are atmospheric in situ CH 4 and C 2 H 6 measurements from the ground and aircraft to apportion CH 4 emissions between NG operations and other sources in the Fayetteville dry shale gas field. C 2 H 6 to CH 4 enhancement ratios in observed emission plumes (ERs, the molar ratio of C 2 H 6 and CH 4 relative to their respective local background values) are compiled downwind of NG facilities in the Fayetteville Shale to characterize the spatial variability of NG ERs throughout the study area. These source ERs are then used together with total emission ERs from aircraft measurements to attribute CH 4 emissions between NG and other sources.
High frequency in situ measurements permitted the determination of enhancement ratios at more locations than would have been possible with discrete flask samples. The ability to quickly make measurements at a given location eliminated the need to correct for the species background, which can vary over time and space. The use of data collected on multiple days allows an assessment of the variability in ERs at different NG facilities and variability among ERs on flights across different days. This method is atmospheric measurement-based with incorporation of publicly available data and does not require inventory or operator-provided data. The unique combination of ground and aircraft data incorporates information about both source signatures and the mix of sources in the area into the attribution. Our work is an effort to further refine top-down measurement and attribution methods to help improve independent constraints on regional long-lived greenhouse gas emission attribution.

Overview
The study area (Figure 1) comprised the core of the Fayetteville Shale production basin in Arkansas, USA, and spanned Cleburne, Conway, Faulkner, Independence, Jackson, Van Buren, and White counties to form an area ~6 5 km (north to south) by ~1 50 km (west to east). This area is a dry NG production region (Arkansas Geological Survey, 2015) with no oil production or processing of NG , and accounted for ~3 .1% of national dry NG production in September 2015 (EIA, 2017). NG facilities in the area in 2015 totaled ~5 500 active NG wells, 125 gathering stations, and six transmission stations. Variability in the C 2 H 6 to CH 4 emission ratio for NG sources is expected to be due to variability in raw gas composition, as there is no gas processing at production facilities or at centralized locations. In addition to NG operations, potential CH 4 sources included agriculture (cattle and poultry farms), wetlands, geologic seeps, stationary combustion, and one small landfill in Faulkner County . In the absence of biomass burning, only NG related sources in the study area are expected to co-emit C 2 H 6 with CH 4 . A fire ban was in effect in much of the study area during the study, so biomass burning is not expected to have occurred in the majority of the study area.
Fast response in situ CH 4 and C 2 H 6 measurements were obtained using both ground-based and airborne platforms. An instrumented van (NOAA Global Monitoring Division (GMD) mobile laboratory) measured CH 4 and C 2 H 6 in CH 4 emission plumes from public roads downwind of NG facilities in order to gather data needed to calculate representative ERs. C 2 H 6 and CH 4 data were also collected by two mobile laboratories operated by Aerodyne Research, Inc. performing dual tracer flux ratio measurements, a downwind method used to determine facility-level CH 4 emissions (Roscioli et al., 2015;Yacovitch et al., 2017). A single-engine instrumented aircraft operated by Scientific Aviation was used to obtain the C 2 H 6 and CH 4 data (Schwietzke et al., 2017) for regional-scale ERs that reflected the mix of CH 4 emissions in the study area.
Mielke-Maday et al: Methane source attribution in a U.S. dry gas basin using spatial patterns of ground and airborne ethane and methane measurements Art. 13,page 3 of 16 Each platform's instrumentation and sampling strategy are described below.

NOAA mobile laboratory
The NOAA mobile laboratory (ML), a 2005 Chevrolet 15 passenger van, was equipped with a four species cavity ring down spectrometer (CRDS) (Picarro model #G2401-m) that provided sequential measurements of CH 4 , CO 2 , and CO mole fractions, and measurements of water vapor every three seconds. A tunable infrared laser direct absorption spectroscopy (TILDAS) analyzer (Aerodyne Research, Inc.) with a 2990 cm -1 laser measured C 2 H 6 , CH 4 , and water vapor at 1 Hz. Stationary 1-sigma noise over 30 seconds under the deployment conditions was approximately 0.30 ppb on CH 4 for the CRDS, 0.05 ppb on C 2 H 6 , and 0.40 ppb on CH 4 for the TILDAS. Prior to and following ML drives each day, three gas standards prepared by NOAA GMD in Boulder, Colorado were measured by the CRDS and TILDAS instruments to identify potential instrument drift (see Supplemental Material). The ML was also equipped with a global positioning system (GPS) and NOAA GMD flask air sampling system. Further details of the ML set-up can be found in the Supplemental Material.
Discrete air samples, used to compare with in situ measurements, were collected using a programmable compressor package (PCP) and programmable flask packages (PFPs). The PCP contains the pump, batteries, and a microprocessor to control sampling, and the PFP consists of a manifold and twelve cylindrical, 0.7 L borosilicate glass flasks that are filled one by one (Aircraft Program Methods, 2018). These discrete air samples were analyzed by NOAA GMD using the same systems used for measurements in the Global Greenhouse Gas Reference Network (ESRL, 2017). CH 4 was analyzed by gas chromatography with flame ionization detection (Dlugokencky et al., 1994), and C 2 H 6 , by gas chromatography mass spectrometry (see Supplemental Material). C 2 H 6 and CH 4 dry air mole fractions have 1-sigma uncertainties of 11 ppt (NOAA 2016 scale) and 1.1 ppb (NOAA X2004A scale), respectively.

Aerodyne mobile laboratory
Aerodyne Research, Inc. operated two mobile laboratories during the time when the NOAA ML was also deployed. Each housed TILDAS analyzers used for measuring CH 4 , C 2 H 6 and the trace gases used in their tracer release method for facility-level quantification of CH 4 emissions. Details of the mobile laboratory set-ups can be found in the supplemental material of Yacovitch et al. (2017).

Scientific aviation aircraft
Instrumentation onboard the Mooney Ovation aircraft included a NOAA GMD flask air sampling system, a cavity ring down spectrometer (Picarro model # 2301-f) with 0.5 Hz water-corrected measurements of CH 4 (<2 ppb uncertainty (Karion et al., 2013) CO 2 , and water vapor, and a TILDAS C 2 H 6 mini trace gas monitor (Aerodyne Research, Inc.). This instrument, with a 2997 cm -1 laser, (~200 ppt noise in the configuration used in this study) was used to obtain 1 Hz measurements of C 2 H 6 . Though C 2 H 6 values were not water-corrected (not removing water during calculations ignores dilution and results in a lower mole fraction than the dry air mole fraction), the highest water content observed on any flight (<2%) would necessitate a data correction that is small relative to the C 2 H 6 signal observed on flights. Additional instrumentation included a Vaisala relative humidity and temperature probe and GPS (Hemisphere VS101). Details of the aircraft set-up can be found in Conley et al. (2017).

Methane emission source attribution model
A top-down attribution of CH 4 emissions between NG and other CH 4 sources was accomplished using the approach described below. Atmospheric ER measurements are assumed to be equivalent to emission ratios from CH 4 sources since the average global atmospheric lifetimes of C 2 H 6 and CH 4 (approximately 2 months and 9 years, respectively (Rudolph and Enhalt, 1981;Nisbet et al., 2016)) make chemical losses negligible when examining plumes and air masses downwind of a facility or source region over timescales on the order of minutes to a few hours. It is assumed that the ER from a given source in the study area does not change from hour-to-hour or day-today due to the fact that the produced gas is dry and there are no gas treatment facilities or liquid condensate collection activities in this basin. Results presented in a later section support this assumption.
For an area of interest, the total CH 4 and C 2 H 6 emissions (E) can be separated into the contributions from NG and other sources (denoted "other").
The relative contribution of CH 4 from NG related sources in the study area can thus be calculated using Eq. 7, the ratio of Eq. 6 to Eq. 5.
The result is a mean estimate of the fraction of the total CH 4 emissions that originate from NG sources in the area of interest during the time of the campaign. This model is only applicable in cases in which there are only two distinct CH 4 source categories with a small range of ERs, as another CH 4 source with its own characteristic ER would introduce another mode in the ER distribution and would require another independent trace gas marker and more equations to constrain the attribution. As previously stated, all C 2 H 6 enhancements are assumed to originate from NG related sources, and the ER distribution for NG point source emissions and regional area source plumes is assumed constant over the time period of the study.

NOAA mobile laboratory
Twelve daytime drives were conducted with the NOAA ML, ranging from six to ten hours each (typically from about 10:00-18:00 local time), to compile ERs in CH 4 plumes throughout the study area. Drives were made on public roads throughout extensive portions of the study area, as shown in Figure 1. Distances from the ML on the road to a facility were typically 0.05 to 0.30 km.
When a CH 4 emission plume (at least 50 ppb above local background) was observed downwind of a potential emission source, the ML was intentionally positioned to acquire 30 seconds to a few minutes of in situ data to measure multiple CH 4 plumes to calculate the ER, check for temporal stability in the ER at the source, and in some instances, collect flask air samples in the local background air and in the facility emission plume. When collecting in situ data, the ML was driven slowly out of the CH 4 plume to allow CH 4 to return to its background level. In this way, a range of CH 4 and C 2 H 6 mole fractions could be measured and a slope calculated to determine an ER. Section 3.1 provides an assessment of the temporal consistency of ERs when measurements at the same facility were repeated on multiple days.
To collect flask air samples, the ML engine was first turned off. Then, the PFP sampling manifold and flask were first flushed (at ~1 0 and ~1 7 SLPM, respectively) for 45 seconds and, when complete, air sample collection was manually triggered using the real-time display of in situ measurements to target different mole fractions. Fill times were ~1 0 seconds. The core GMD chemical analysis instrumentation focuses on tracking background air composition, and to avoid contamination of the GC-FID columns used to measure samples in the Global Greenhouse Gas Reference Network, we intentionally only collected PFP air samples with CH 4 mole fractions below 5 ppm.

Aerodyne mobile laboratory
Included in our analysis are ERs from Aerodyne Research, Inc. for which: (i) the CH 4 emissions at the NG facility could be determined using the dual tracer release method, (ii) there was a coefficient of determination of 0.50 or greater, and (iii) manual inspection of the CH 4 and C 2 H 6 plumes did not indicate any problems with the Aerodyne dual tracer release at the target site, as described in Section 6 of the Supplemental Material in Yacovitch et al. (2017).
Here we use only ERs for which the reported 95% confidence interval is less than or equal to 30% of the mean ER, leaving 46 of the original 49 facilities at which ERs were obtained by Yacovitch et al. (2017).

Scientific aviation aircraft
Fifteen flights were conducted between September 21 and October 14, 2015. Data from October 1 and October 2 were used by Schwietzke et al. (2017) to calculate a total CH 4 emission estimate for the study area using a mass balance approach. Flights with box, raster, or spiral sampling patterns were conducted over the entirety and subsections of the study area to gather data for CH 4 emission quantification, CH 4 source attribution, and individual facility CH 4 emission quantification, respectively (see Supplemental Material). Aircraft measurements from five raster flights, the paths of which are shown in Figure 2, were used to investigate mean ERs for the western and eastern halves of the study area.

Analysis of ground measurements in facilitylevel plumes
Mobile laboratory data were used to determine representative ERs for NG sources in the study area. CH 4 plumes measured by the ML were identified in CH 4 data from the CRDS using an instantaneous 50 ppb CH 4 enhancement threshold above the local CH 4 background. This threshold was chosen to distinguish enhancements caused by local emissions from instrument noise and background variability without biasing plume selection toward relatively large CH 4 plumes. The start time of each plume was determined by the time immediately before CH 4 increased by 50 ppb and the end time was determined by the time CH 4 returned to its start time CH 4 value, with low variability in CH 4 prior to and following the plume. Plumes were identified by manually screening the time series multiple times to prevent omissions and errors. If data were acquired when the ML was moving, only plumes less than one minute in duration were included in the analysis to avoid including multiple facility plumes and potentially varying background levels. Using code described by Yacovitch et al. (2014), a C 2 H 6 vs. CH 4 scatterplot was created for each plume with data from the TILDAS, and the ER for each plume was determined using the slope calculated by orthogonal distance regression. Because the C 2 H 6 and CH 4 data were acquired from the same absorption spectra on the same instrument, no sampling time offset optimization was required. Additionally, no assignment of C 2 H 6 and CH 4 background was needed for each plume because in the method described above, the background is accounted for in the intercept of the scatterplot and the slope determined by enhancements alone, reflecting the C 2 H 6 to CH 4 ratio in isolated plumes.
An ER was assigned to a source only if the coefficient of determination, R 2 , was at least 0.65. For CH 4 sources where more than one ER was determined, a time-weighted average was calculated to assign a single mean ER to the source. Only results for which a facility location could be readily identified as the emission plume source are presented. CH 4 plumes that could not be attributed to a specific source are not included in the analysis. Plumes from 71 NG facilities measured by the NOAA ML passed the filtering criteria and though the minimum R 2 required was 0.65, 85% of the NG plumes had an R 2 associated with the ER that was 0.95 or greater.

Analysis of aircraft measurements for areascale plumes
Aircraft measurements in the boundary layer over the study area were used to determine area-scale plume ERs (see Eq. 5). The portion of the data for which CH 4 was enhanced above background was identified on downwind transects (for box pattern flights) and legs (for raster pattern flights) and an ER was calculated using orthogonal distance regression weighted by the airborne instrumentation noise of 1.3 ppb for CH 4  and 200 ppt for C 2 H 6 . We divide the study area into western and eastern halves using -92.1° longitude according to observed spatial differences in aircraft-based CH 4 emission estimates normalized by NG production . In the event that multiple transects or legs from a single flight were downwind of the study area, the ERs for all transects or legs were averaged. Typical ranges observed while flying transects were approximately 1940 to 2300 ppb for CH 4 and 1.6 to 4.0 ppb for C 2 H 6 . Mole fractions observed in the aircraft data are lower than those observed in the ML, where CH 4 and C 2 H 6 enhancements very close to sources could be up to several thousand ppb CH 4 and several hundred ppb C 2 H 6 . As the aircraft measured downwind of different regions of the study area, background CH 4 levels sometimes changed. To distinguish background variability from emission sources, ERs were determined only for those subsections of a leg or transect that included CH 4 enhancements of at least 50 ppb. Criteria for acceptance of ERs were that measurements were made in the planetary boundary layer, a CH 4 enhancement in area plumes of at least 50 ppb above background and a C 2 H 6 enhancement of at least 1 ppb to ensure a high enough signal to noise ratio, and a coefficient of determination of at least 0.65. After filtering using the above criteria, flight pattern, and wind direction, five flights could be used in the determination of an area ER for the western half of the study area. Several transects in these flight and others had low signal to noise and low enhancements, disqualifying them for use in ER calculations. Low signal to noise occurred due to low CH 4 and C 2 H 6 enhancements in the planetary boundary layer and was observed more often over the eastern half of the study area, where total CH 4 emissions were found to be lower than in the western half .

Mobile laboratory drives
Facility plume C 2 H 6 to CH 4 enhancement ratios on different days Seven of eight facilities at which plume measurements were repeated over multiple days show consistent (difference <20%) ERs across days, as shown in Table 1. This excludes 106 facilities with measurements on only one day (see next subsection). ERs shown in Table 1 are for each plume measured at the facility rather than the timeweighted average for the facility. All ERs have an R 2 of 0.65 or greater, and 85% of the NG plumes had an R 2 associated with the ER that was 0.95 or greater. The repeatability of ERs on different days at these facilities suggests that an ER measured at a facility by the ML on any given day can be used as a representative value for that facility even without repeated measurements on different days and that emission ratios at NG facilities were likely constant over the course of the study. The spread in ERs observed at Production Pad 7 on September 24 suggests another emission source in addition to vented gas. CH 4 plumes at this pad that are enhanced in C 2 H 6 also show a strong correlation with the combustion product, CO 2 (see Supplemental Material), though no correlation between CO and CH 4 was observed. The most likely sources of combustion at production pads in this dry gas basin are on-site NG-powered compressor engines, used at some production pads to increase flow from a well. It is possible that the natural gas-fueled compressor was running at the time of the ML measurements, though the activity data needed for confirmation are not available. Using ERs only from plumes for which the CH 4 enhancement was at least 50 ppb, the time-averaged ER for production pad 7 is 1.7%. Low temporal variability between ERs at a NG facility is expected because the ER from mostly fugitive (in this region) emissions is dependent upon the composition of the raw natural gas being extracted, which remains constant over the timescale at which this study was conducted. However, variability in natural gas composition due to geologic differences throughout a region can lead to spatial variability in ERs.

Range and spatial distribution of C 2 H 6 to CH 4 ERs observed at different NG facility types
In this section, we present ER results from the 114 NG facilities measured by the NOAA ML and the Aerodyne mobile laboratories (AML) . For the three facilities that were measured by both NOAA and the AML, the ERs at both gathering stations are identical (1.0% and 1.4% for the two gathering stations) and within 10% for the production pad (1.1% measured by NOAA and 1.0% measured by the AML). The dataset presented from this point forward consists of results from both mobile laboratories.
ERs in NG facility plumes, including production pads, gathering stations, transmission stations, pipelines and dehydrator stations, display some spatial variability throughout the study area, especially on the eastern side, as shown in Figure 3. Table 2 shows the range of ERs observed at each type of NG facility measured on the ground in the western and eastern halves of the study area. The sample size of three of the NG source types (pipeline junctions, dehydrator stations, and transmission stations) is small, making it  difficult to draw conclusions about the observed ER ranges at these types of facilities. More measurements were made throughout the study area at production pads and gathering stations, which represent the majority (~98%) of the active facilities and about 97% of estimated CH 4 emissions  in the study area. Production pads and gathering stations display similar ER ranges of 0.8 to 7.8% and 1.0 to 7.9%, respectively. The similar ER ranges between these two types of facilities are to be expected given the lack of gas processing in this dry gas field, which has very low water and non-methane hydrocarbon content and no co-produced natural gas condensate or oil. The ranges of ERs observed at production pads are similar to gathering stations within the western half and within the eastern half of the study area. In the western half, production pad ERs range from 0.8 to 2.1% and gathering stations from 1.0 to 1.7%. In the eastern half, production pad ERs range from 0.9 to 7.8% and gathering stations from 1.1 to 7.9%. Facilities in the eastern half have a greater range of ERs that tend toward higher values. Only one of 13 facilities with an ER of 2.1% or greater is located in the western half. Of the remaining facilities with relatively high ERs, all of which are located in the eastern half, eight are located in White county (2.1-5.7%), and four are located in Independence county (6.2-7.9%). Though the range of ERs observed at NG facilities in the eastern half is larger, facilities with relatively high ERs tend to be in the same geographic subregion of the study area (see Figure 3).
At the gathering station at which a 7.9% ER was observed, combustion was evident, as CH 4 plumes correlated with both CO 2 and CO. At gathering stations, which have on-site compressors, combustion processes are expected but CH 4 correlation with combustion tracers was observed at only one of the eight gathering stations surveyed by the NOAA ML. It is possible that at the other stations, hot combustion plumes rose too high and that the ML was not far enough downwind to be able to measure the plume where it reached the surface. However, even though combustion occurs at compressor stations, higher ERs at these facilities are not necessarily expected because combustion slip (unburned fuel in the compressor engine exhaust) has been shown to be the main source of CH 4 emissions at gathering stations in the study area . Figure 4 shows the distribution of ERs measured at NG facilities by the MLs. The majority (64%) of NG ERs observed in the western half range from 1.2 to 1.5%, with a mode of 1.3%. The eastern half has a similar range of frequent ERs of 1.2 to 1.5% with a peak in the distribution at 1.3% also. However, there is significant variability in ER in the eastern half, where 34% of measured facilities have ERs of 2% or greater. Higher ERs in the eastern half are confirmed by comparison with C 2 H 6 to CH 4 ratios from gas composition data provided by operators for a subset of production pads and gathering stations (Figure 5). In cases where results from multiple gas composition samples were provided for a facility, the mean C 2 H 6 to CH 4 ratio was used. Gas composition data exhibit modes for both the western and eastern halves of 1.3 to 1.5% respectively, consistent with ML plume data. Composition data also reveal a portion of facilities in the eastern half with C 2 H 6 to CH 4 ratios greater than 2%. This portion is 28% compared with the 34% determined by ML measurements. Like the ML results, all facilities in the eastern half with C 2 H 6 to CH 4 ratios of at least 2% are located in White and Independence counties. Higher ERs in the eastern half are likely a result of the thermal maturity of the geologic formation being drilled. Thermal maturity refers to the extent to which reactions involving heat that convert organic matter to oil and natural gas have occurred. Generally, natural gas is expected to be wetter (higher C 2 H 6 to CH 4 ratio) as the thermal maturity of its source rock decreases (Visschedijk et al., 2018).
In the Barnett Shale play in Texas, where natural gas is produced throughout the entirety of the study area, the distribution of 172 NG ERs measured by Yacovitch et al. (2015) showed a peak at approximately 1.5%. The highest ERs were observed in the northwest, where oil is co-produced.
The aircraft performed spiral flight patterns around several potential CH 4 sources, including non-NG sources, and ERs were calculated. Although results are not used in the attribution, they can be found in the Supplemental Material. Mielke-Maday et al: Methane source attribution in a U.S. dry gas basin using spatial patterns of ground and airborne ethane and methane measurements Art. 13, page 9 of 16

Measurements at poultry farms
Arkansas is among the top states for the production of poultry in the United States, ranking third and fourth for broiler production in terms of head of chickens and pounds produced, respectively, in 2015 (USDA, 2015). While chickens are not known to be a CH 4 source as a result of enteric fermentation like cattle, CH 4 could be emitted from their manure (Monteny et al., 2001). Chicken litter, a combination of manure and the bedding in the chicken house, is often temporarily stored on-site until it can be removed and applied as a fertilizer for crops. However, best management practices recommend covered (dry) storage, which serves to minimize the potential for CH 4 and ammonia production (Ogejo et al., 2009).
Commercial poultry farms in the study area were long enclosed buildings with fans for ventilation and no visible liquid manure storage ponds. Measurements were made downwind of nine poultry farms (see Supplemental Material for locations), identified by the strong downwind localized smell of ammonia. At the time of the study, 292 poultry farms were known to be in the study area. Using the same 50 ppb threshold above the local background as was used to identify NG CH 4 plumes, poultry farms exhibited no substantial CH 4 enhancements. For a poultry farm measurement shown in Figure 6, CH 4 enhancements were approximately 10 ppb during each ML transect. For comparison, background variability for surface CH 4 on ML drives had a range of approximately 300 ppb. Based on these data, it is not expected that poultry contributed significantly to total CH 4 emissions in the study area. The bottom-up CH 4 emissions estimate for poultry farms provided by  of <200 kg CH 4 /hr for the entire study area also suggests that poultry farms are a small source of CH 4 emissions when compared to the total (30.6 tons CH 4 /hr).

Determination of mean NG C 2 H 6 to CH 4 ERs
To derive mean ERs for the full study area as well as the western and eastern halves, we calculate spatially distinct means of the ML NG facility plume ERs weighted by NG production volumes using a Delaunay Triangulation (see Supplemental Material). This method accounts for spatial variability in the ER across each half and for sections of the study area with less measurement coverage than others. The Delaunay Triangulation algorithm in MATLAB connects the ER measurement locations into triangles (using their GPS coordinates) so that no other measurement location lies in a circumcircle of a given triangle and the value of the minimum angle of each triangle is maximized (See Figure S5). The varying sizes of the resulting triangles account for disparity in measurement density in the study area and ensure that there are no triangles without associated ER measurements. Each vertex of a triangle represents a measured ER. For each triangle, the three ERs are averaged and the mean value is assigned as the mean ER for all NG emissions in that triangle. Using September 2015 NG production data from the Arkansas Oil and Gas Commission (AOGC) (AOGC, 2015), the NG production of all production pads contained within each triangle was summed. With this production data and the mean ER for each triangle, a production weighted mean NG ER was calculated for the study area, western half, and eastern half. A simulation was performed in which a simple random sample of 80% of the NOAA and AML facility ER data points was selected without replacement and the Delaunay Triangulation method was applied. Performing this simulation 30 times generated a mean ER and standard deviation, which was used to assign an uncertainty to the mean NG ERs generated by Delaunay Triangulation. The resulting mean NG ERs are 1.68 ± 0.10% for the full study area, 1.34 ± 0.03% for the western half, and 2.10 ± 0.22% for the eastern half.

Determination of area-scale C 2 H 6 to CH 4 ERs from aircraft data
While ML data were used to determine ERs characteristic of NG operations in the study area, aircraft data were used to calculate area-scale ERs that include CH 4 emissions from all sources in the study area. CH 4 plumes measured by the aircraft are from a combination of multiple point sources after emissions have been transported from their sources and mixing has occurred downwind. NG source ERs and area-scale ERs were then combined (Eq. 7) to estimate the portion of the total CH 4 emissions in the study area originating from NG operations.
Area-scale ERs were calculated by fitting orthogonal distance regression slopes through the C 2 H 6 vs. CH 4 in situ aircraft measurements for multiple flight legs in the boundary layer downwind of or over the study area and calculating a mean for each flight. The resulting aircraft ERs are shown in Table 3. ERs for each transect or leg of each flight used for attribution can be found in the Supplemental Material. Figure 7A shows an example raster flight track for the determination of an ER for one leg (leg 4). For this flight on October 5, 2015 with wind from the North, CH 4 enhancements were observed over the western half of the study area, beginning with leg 4 (starting upwind in the northwestern corner and moving further southeast). The portion of leg 4 for which CH 4 is enhanced is shown with a box in Figure 7B and the corresponding C 2 H 6 to CH 4 scatterplot in Figure 7C. The slope is 1.25 ± 0.04% (R 2 = 0.88).
Representative area-scale ERs were calculated for the western half using five raster flights conducted on September 22, October 2, October 5, October 6 and October 7, 2015 (Table 3), and a mean C 2 H 6 to CH 4 slope for the enhancements observed over and/or downwind of the western half of 1.21 ± 0.02% was calculated.
While this result is higher than the 0.6% determined by Peischl et al. (2015) from two flights over the Fayetteville Shale, the ER reported in this study is only for the western half of the study area, whereas the 0.6% reported by Peischl et al. is for the entire region. The mean areascale ER of 1.21 ± 0.02% from this study falls within the range of ~1 -3% from dry gas sources with a mode of 1.8% determined from flights in the Barnett Shale .
Larger enhancements over the western half permitted calculation of area-scale ERs whereas in the eastern half, C 2 H 6 signal with the same magnitude as the measurement noise made calculating ERs impossible. CH 4 emissions from NG operations normalized by NG production were found to be about one half in the eastern half relative to the western half on two days . Since C 2 H 6 is co-emitted with CH 4 , it is likely that C 2 H 6 emissions are also lower in the eastern half, and a spatially-resolved bottom-up CH 4 emissions inventory of the total study area showed no substantial differences in non-NG related CH 4 emissions between the western and eastern half. Assuming similar atmospheric mixing and dispersion conditions in the eastern half as in the western half, lower emissions would produce smaller enhancements, which matches observations. Though C 2 H 6 mixing ratios over the eastern half of the study area were generally lower than those observed over the western half, a raster flight on October 6 revealed C 2 H 6 enhancements in the northeastern section of the study area (Figure 8). ERs on the four legs of this flight with C 2 H 6 mixing ratios above 3.5 ppb (shown in light blue on Figure 8) were 1.5 ± 0.1%, 2.0 ± 0.1%, 2.5 ± 0.2%, and 7.1 ± 0.2%. The plume with an ER of 7.1% (shown in red on Figure 8) represents a localized enhancement in both CH 4 and C 2 H 6 .

Attribution of emissions
The mean fraction of CH 4 emissions attributable to NG sources in the western half of the study area was calculated using Eq. 7. With [C 2 H 6 /CH 4 ] area = 1.21 ± 0.02% and [C 2 H 6 /CH 4 ] NG = 1.34 ± 0.03%, the mean estimate of the portion of the total CH 4 emissions originating from NG related sources in the western half of the study area was 90 ± 5% (1-sigma), and the portion of the total CH 4 emissions originating from other sources was estimated to be 10 ± 5%.
Our estimate for the portion of total CH 4 emissions originating from NG sources compares well with bottom-up estimates made for the western half of the study area by Vaughn et al. (2018). Bottom-up emission estimates were made for the two afternoons during which the mass balance flights used by Schwietzke et al. (2017) to calculate total CH 4 emissions occurred. Vaughn et al. (2018) estimated 89.9% and 87.8% of total CH 4 emissions in the western half to have originated from NG sources on October 1 and October 2, 2015, respectively, yielding a mean of 88.8%. The highest estimated CH 4 emissions from non-NG sources were from beef cattle and geologic seeps .
Several factors prevented attribution of CH 4 emissions in the eastern half. First, CH 4 and C 2 H 6 enhancements measured by the aircraft in the planetary boundary layer over the eastern half were low compared to those in the western half (as shown in Figure 8 and by Schwietzke et al. (2017). In particular, the signal to noise ratio for the C 2 H 6 enhancements onboard the aircraft was too low over the eastern half to allow estimation of ERs for aerial plumes mapped over the subregion. The C 2 H 6 enhancements were also too low along the downwind transects of the two mass balance flights analyzed by Schwietzke et al. (2017) to determine a mass balance estimation of the C 2 H 6 emissions for the entire study region, with C 2 H 6 ranges on the downwind transects on October 1 st and 2 nd of 0.8 to 2.4 ppb and 1.0 to 2.0 ppb, respectively.

Conclusions
Ground-based, atmospheric measurements were used to determine a representative range of C 2 H 6 to CH 4 facilitylevel plume ERs, confirmed by produced gas composition data. Using airborne and mobile surface atmospheric measurements of C 2 H 6 and CH 4 and a simple two source mixing model, we show that 85-95% of total CH 4 emissions in the western half of the study area originate from NG operations. The western half represents 76% of total study area CH 4 emissions based on aircraft measurements . Our result agrees with the NG CH 4 emissions estimate from a bottom-up model combining field measurements and spatiotemporal NG operations data by Vaughn et al. (2018). Attribution of CH 4 emissions in the eastern half of the study area was not possible due to low boundary layer enhancements in CH 4 and C 2 H 6 . Low signal to noise for C 2 H 6 also prohibited the use of the aircraft mass balance model for obtaining a C 2 H 6 flux in both halves of the study area. Additionally, the assumption of a single mode NG ER distribution is not valid for the eastern half due to the more complex NG ER distribution.
The spatial variability in ER across the study area, as demonstrated by the difference in ER distributions in the western half compared with the eastern half, is important to account for when developing aggregated source emission signatures in a NG producing region. Although the majority of NG ERs observed in both the western and eastern halves fall in the same range of 1.2 to 1.5%, the eastern half has a second mode in its distribution, in which 34% of measured facilities have ERs of 2% or greater.
The Fayetteville Shale was chosen for the field campaign partly because CH 4 emissions from sources other than NG systems were expected to be low, making CH 4 source attribution less complicated than in areas with many potential CH 4 sources. Results from this study confirm that NG operations are the primary source of CH 4 in the western portion of the study area. Poultry farms were not found to be a significant source of CH 4 .
Use of a ML facilitated the collection of an extensive amount of data on specific types of NG facilities and their C 2 H 6 to CH 4 source signatures. For facility emission sampling, the ML was positioned directly adjacent to and downwind of the target facility such that ERs could be assigned to specific facilities. Fast response C 2 H 6 and CH 4 in situ measurements made it possible to determine ERs in multiple plumes downwind of individual emission sources, yielding robust facility-level statistics compared to those from a simple linear regression using a small number of discrete flask samples collected in sources plumes or background conditions at different locations and times. Since ML routes were restricted only by public road access, it was possible to perform repeat measurements downwind of several facilities and examine day-today ER variability. With the exception of one production pad at which some combustion was likely occurring, ERs at all seven facilities with repeated measurements over several days were consistent over time (<10% variation), suggesting that measurements of C 2 H 6 and CH 4 by ML surveys represent daily ERs at a given facility during the campaign.
Aircraft measurements in the atmospheric boundary layer were critical for obtaining information about areaintegrated emission signals. By flying over and downwind of the entire western half of the study area, where CH 4 and C 2 H 6 enhancements were greatest, the aircraft sampled area plumes in the western half and provided a representative ER for the mix of emission sources. The attribution model used in this study assumes that all C 2 H 6 emissions originate from NG sources. Additionally, mean ERs for NG and total emissions are assumed to remain constant from day to day. The two-source model can only be applied in a simple case with two source categories with a single mode each, and would not be appropriate for use in an area with more complicated mix of sources and source signatures, including wetter gas and oil production, which could introduce other modes in the ER distribution. A complex mix of co-located sources with several different source signatures can present a challenge in that multiple tracers or constraints may be required for attribution. In such a case, a representative sample population of facility plumes with sufficient signal to noise may be used as a proxy for the mix of emissions as in Smith et al. (2015), preferably for multiple days and different flight paths. To further advance the value of top-down methods, it would be beneficial to test how stable the emission estimates and attribution solution are for different flight patterns, wind configurations, and dates/seasons, and to have access to composition data for key emission sources in a study region.

Data Accessibility Statement
Enhancement ratios obtained by the mobile laboratories are uploaded in Microsoft Excel .xlsx format. NOAA mobile laboratory C 2 H 6 and CH 4 enhancements relative to the local background are provided for ten different facilities (Sites 1-10) as text files. Files for the remaining 104 facilities can be provided upon request. Due to a governance agreement with project sponsors, locations of high frequency in situ ground measurements cannot be published. Aircraft data for the flights used for attribution and all discrete flask data are provided. Data are available at the following address: ftp://aftp.cmdl.noaa.gov/data/ campaign/mls/RPSEA2015/.

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