Robert Alexander Field

Methane mitigation is planned within administrative boundaries, but atmospheric transport redistributes methane across jurisdictions. Here we show that the methane-column response to Shanxi coal-mine emissions over the Beijing-Tianjin-Hebei megacity region is episodically amplified, while a key terrain gate is persistently undersampled by TROPOMI. A 4.5-year TROPOMI record, supported by independent Xianghe TCCON observations, demonstrates that the associated atmospheric organization and receptor relationships recur significantly beyond the originally simulated summer, including after excluding the entire 2022 calendar year. Controlled WRF-GHG experiments show that observation-constrained refinement of the Shanxi inventory increases mean downwind methane column mass by 15.2×10 ³ tons relative to EDGAR, equivalent to 5.1% of the monthly regional emission benchmark (EBE), and that this response is amplified 3.36-fold under strong southwesterly flow. Valid TROPOMI retrievals cover only 0.13% of the gate-grid days, versus 8.6% regionally, with matched methane-carbon monoxide diagnostics linking the deficit to terrain-dependent screening. These findings reframe cross-boundary methane attribution as an episode resolved source-gate-receptor problem, linking source representation, atmospheric state, and observational coverage. Effective verification therefore requires coordinated source, pathway and receptor observations across the coupled airshed.

Other Test Method 33A (OTM 33A) is a near-source flux measurement method developed by the Environmental Protection Agency (EPA) primarily used to locate and estimate emission fluxes of methane from oil and gas (O&G) production facilities without requiring site access. A recent national estimate of methane emissions from O&G production included a large number of flux measurements of upstream O&G facilities made using OTM 33A and concluded the EPA National Emission Inventory underestimates this sector by a factor of ∼2.1 (Alvarez et al., 2018). The study presented here investigates the accuracy of OTM 33A through a series of test releases performed at the Methane Emissions Technology Evaluation Center (METEC), a facility designed to allow quantified amounts of natural gas to be released from decommissioned O&G equipment to simulate emissions from real facilities (Fig. 1). This study includes test releases from single and multiple points, from equipment locations at different heights, and spanned methane release rates ranging from 0.16 to 2.15 kg h⁻¹. Approximately 95 % of individual measurements (N=45) fell within ±70 % of the known release rate. A simple linear regression of OTM 33A versus known release rates at the METEC site gives an average slope of 0.96 with 95 % CI (0.66,1.28), suggesting that an ensemble of OTM 33A measurements may have a small but statistically insignificant low bias.

Flux estimates of volatile organic compounds (VOCs) from oil and gas (O&G) production facilities are fundamental in understanding hazardous air pollutant concentrations and ozone formation. Previous off-site emission estimates derive fluxes by ratioing VOCs measured in canisters to methane fluxes measured in the field. This study uses the Environmental Protection Agency's Other Test Method 33A (OTM 33A) and a fast-response proton transfer reaction mass spectrometer to make direct measurements of VOC emissions from O&G facilities in the Upper Green River Basin, Wyoming. We report the first off-site direct flux estimates of benzene, toluene, ethylbenzene and xylenes from upstream O&G production facilities and find that these estimates can vary significantly from flux estimates derived using both the canister ratio technique and from the emission inventory. The 32 OTM 33A flux estimates had arithmetic mean (and 95% CI) as follows: benzene 17.83 (0.22, 98.05) g/h, toluene 34.43 (1.01, 126.76) g/h, C8 aromatics 37.38 (1.06, 225.34) g/h, and methane 2.3 (1.7, 3.1) kg h-1. 20% of facilities measured account for ~67% of total BTEX emissions. While this heavy tail is less dramatic than previous observations of methane in other basins, it is more prominent than predicted by the emission inventory.

Atmospheric methane emissions from active natural gas production sites in normal operation were quantified using an inverse Gaussian method (EPA’s OTM 33a) in four major U.S. basins/plays: Upper Green River (UGR, Wyoming), Denver-Julesburg (DJ, Colorado), Uintah (Utah), and Fayetteville (FV, Arkansas). In DJ, Uintah, and FV, 72–83% of total measured emissions were from 20% of the well pads, while in UGR the highest 20% of emitting well pads only contributed 54% of total emissions. The total mass of methane emitted as a percent of gross methane produced, termed throughput-normalized methane average (TNMA) and determined by bootstrapping measurements from each basin, varied widely between basins and was (95% CI): 0.09% (0.05–0.15%) in FV, 0.18% (0.12–0.29%) in UGR, 2.1% (1.1–3.9%) in DJ, and 2.8% (1.0–8.6%) in Uintah. Overall, wet-gas basins (UGR, DJ, Uintah) had higher TNMA emissions than the dry-gas FV at all ranges of production per well pad. Among wet basins, TNMA emissions had a strong negative correlation with average gas production per well pad, suggesting that consolidation of operations onto single pads may reduce normalized emissions (average number of wells per pad is 5.3 in UGR versus 1.3 in Uintah and 2.8 in DJ).

This study presents the results of a campaign that estimated methane emissions at 268 gas production facilities in the Fayetteville shale gas play using onsite measurements (261 facilities) and two downwind methods – the dual tracer flux ratio method (Tracer Facility Estimate – TFE, 17 facilities) and the EPA Other Test Method 33a (OTM33A Facility Estimate – OFE, 50 facilities). A study onsite estimate (SOE) for each facility was developed by combining direct measurements and simulation of unmeasured emission sources, using operator activity data and emission data from literature. The SOE spans 0–403 kg/h and simulated methane emissions from liquid unloadings account for 88% of total emissions estimated by the SOE, with 76% (95% CI [51%–92%]) contributed by liquid unloading at two facilities. TFE and SOE show overlapping 95% CI between individual estimates at 15 of 16 (94%) facilities where the measurements were paired, while OFE and SOE show overlapping 95% CI between individual estimates at 28 of 43 (65%) facilities. However, variance-weighted least-squares (VWLS) regressions performed on sets of paired estimates indicate statistically significant differences between methods. The SOE represents a lower bound of emissions at facilities where onsite direct measurements of continuously emitting sources are the primary contributor to the SOE, a sub-selection of facilities which minimizes expected inter-method differences for intermittent pneumatic controllers and the impact of episodically-emitting unloadings. At 9 such facilities, VWLS indicates that TFE estimates systematically higher emissions than SOE (TFE-to-SOE ratio = 1.6, 95% CI [1.2 to 2.1]). At 20 such facilities, VWLS indicates that OFE estimates systematically lower emissions than SOE (OFE-to-SOE ratio of 0.41 [0.26 to 0.90]). Given that SOE at these facilities is a lower limit on emissions, these results indicate that OFE is likely a less accurate method than SOE or TFE for this type of facility.