Methane emissions from the U.S. oil and natural gas supply chain were estimated by using ground-based, facility-scale measurements and validated with aircraft observations in areas accounting for ~30% of U.S. gas production. When scaled up nationally, our facility-based estimate of 2015 supply chain emissions is 13 ± 2 teragrams per year, equivalent to 2.3% of gross U.S. gas production. This value is ~60% higher than the U.S. Environmental Protection Agency inventory estimate, likely because existing inventory methods miss emissions released during abnormal operating conditions. Methane emissions of this magnitude, per unit of natural gas consumed, produce radiative forcing over a 20-year time horizon comparable to the CO from natural gas combustion. Substantial emission reductions are feasible through rapid detection of the root causes of high emissions and deployment of less failure-prone systems.

Based on a uniquely dense network of surface towers measuring continuously the atmospheric concentrations of greenhouse gases (GHGs), we developed the first comprehensive monitoring systems of CO 2 emissions at high resolution over the city of Indianapolis. The urban inversion evaluated over the 2012-2013 dormant season showed a statistically significant increase of about 20% (from 4.5 to 5.7 MtC ± 0.23 MtC) compared to the Hestia CO 2 emission estimate, a state-of-the-art building-level emission product. Spatial structures in prior emission errors, mostly undetermined, appeared to affect the spatial pattern in the inverse solution and the total carbon budget over the entire area by up to 15%, while the inverse solution remains fairly insensitive to the CO 2 boundary inflow and to the different prior emissions (i.e., ODIAC). Preceding the surface emission optimization, we improved the atmospheric simulations using a meteorological data assimilation system also informing our Bayesian inversion system through updated observations error variances. Finally, we estimated the uncertainties associated with undetermined parameters using an ensemble of inversions. The total CO 2 emissions based on the ensemble mean and quartiles (5.26-5.91 MtC) were statistically different compared to the prior total emissions (4.1 to 4.5 MtC). Considering the relatively small sensitivity to the different parameters, we conclude that atmospheric inversions are potentially able to constrain the carbon budget of the city, assuming sufficient data to measure the inflow of GHG over the city, but additional information on prior emission error structures are required to determine the spatial structures of urban emissions at high resolution.

The Indianapolis Flux Experiment (INFLUX) aims to develop and assess methods for quantifying urban greenhouse gas emissions. Here we use CO 2 , 14 CO 2 , and CO measurements from tall towers around Indianapolis, USA, to determine urban total CO 2 , the fossil fuel derived CO 2 component (CO 2 ff), and CO enhancements relative to background measurements. When a local background directly upwind of the urban area is used, the wintertime total CO 2 enhancement over Indianapolis can be entirely explained by urban CO 2 ff emissions. Conversely, when a continental background is used, CO 2 ff enhancements are larger and account for only half the total CO 2 enhancement, effectively representing the combined CO 2 ff enhancement from Indianapolis and the wider region. In summer, we find that diurnal variability in both background CO 2 mole fraction and covarying vertical mixing makes it difficult to use a simple upwind-downwind difference for a reliable determination of total CO 2 urban enhancement. We use characteristic CO 2 ff source sector CO:CO 2 ff emission ratios to examine the contribution of the CO 2 ff source sectors to total CO 2 ff emissions. This method is strongly sensitive to the mobile sector, which produces most CO. We show that the inventory-based emission product ("bottom up") and atmospheric observations ("top down") can be directly compared throughout the diurnal cycle using this ratio method. For Indianapolis, the top-down observations are consistent with the bottom-up Hestia data product emission sector patterns for most of the diurnal cycle but disagree during the nighttime hours. Further examination of both the top-down and bottom-up assumptions is needed to assess the exact cause of the discrepancy.

We present estimates of regional methane (CH4) emissions from oil and natural gas operations in the Barnett Shale, Texas, using airborne atmospheric measurements. Using a mass balance approach on eight different flight days in March and October 2013, the total CH4 emissions for the region are estimated to be 76 ± 13 × 10(3) kg hr(-1) (equivalent to 0.66 ± 0.11 Tg CH4 yr(-1); 95% confidence interval (CI)). We estimate that 60 ± 11 × 10(3) kg CH4 hr(-1) (95% CI) are emitted by natural gas and oil operations, including production, processing, and distribution in the urban areas of Dallas and Fort Worth. This estimate agrees with the U.S. Environmental Protection Agency (EPA) estimate for nationwide CH4 emissions from the natural gas sector when scaled by natural gas production, but it is higher than emissions reported by the EDGAR inventory or by industry to EPA's Greenhouse Gas Reporting Program. This study is the first to show consistency between mass balance results on so many different days and in two different seasons, enabling better quantification of the related uncertainty. The Barnett is one of the largest production basins in the United States, with 8% of total U.S. natural gas production, and thus, our results represent a crucial step toward determining the greenhouse gas footprint of U.S. onshore natural gas production.

Urban environments are the primary contributors to global anthropogenic carbon emissions. Because much of the growth in CO 2 emissions will originate from cities, there is a need to develop, assess, and improve measurement and modeling strategies for quantifying and monitoring greenhouse gas emissions from large urban centers. In this study the uncertainties in an aircraft-based mass balance approach for quantifying carbon dioxide and methane emissions from an urban environment, focusing on Indianapolis, IN, USA, are described. The relatively level terrain of Indianapolis facilitated the application of mean wind fields in the mass balance approach. We investigate the uncertainties in our aircraft-based mass balance approach by (1) assessing the sensitivity of the measured flux to important measurement and analysis parameters including wind speed, background CO 2 and CH 4 , boundary layer depth, and interpolation technique, and (2) determining the flux at two or more downwind distances from a point or area source (with relatively large source strengths such as solid waste facilities and a power generating station) in rapid succession, assuming that the emission flux is constant. When we quantify the precision in the approach by comparing the estimated emis-sions derived from measurements at two or more downwind distances from an area or point source, we find that the minimum and maximum repeatability were 12 and 52 %, with an average of 31 %. We suggest that improvements in the experimental design can be achieved by careful determination of the background concentration, monitoring the evolution of the boundary layer through the measurement period, and increasing the number of downwind horizontal transect measurements at multiple altitudes within the boundary layer.

Abstract. We performed an atmospheric inversion of the CO 2 fluxes over Iowa and the surrounding states, from June to December 2007, at 20 km resolution and weekly timescale. Eight concentration towers were used to constrain the carbon balance in a 1000×1000 km 2 domain in this agricultural region of the US upper midwest. The CO 2 concentrations of the boundaries derived from CarbonTracker were adjusted to match direct observations from aircraft profiles around the domain. The regional carbon balance ends up with a sink of 183 Tg C±35 Tg C over the area for the period June-December, 2007. Potential bias from incorrect boundary conditions of about 0.55 ppm over the 7 months was corrected using mixing ratios from four different aircraft profile sites operated at a weekly time scale, acting as an additional source of uncertainty of 24 Tg C. We used two different prior flux estimates, the SiBCrop model and the inverse flux product from the CarbonTracker system. We show that inverse flux estimates using both priors converge to similar posterior estimates (20 Tg C difference), in our reference inversion, but some spatial structures from the prior fluxes remain in the posterior fluxes, revealing the importance of the prior flux resolution and distribution despite the large amount of atmospheric data available. The retrieved fluxes were compared to eddy flux towers in the corn and grassland areas, revealing an improvement in the seasonal cycles between the two compared to the prior fluxes, despite large absolute differences due to representation errors. The uncertainty of 34 Tg C (or 34 g C m 2 ) was derived from the posterior uncertainty obtained with our reference inversion of about 25 to 30 Tg C and from sensitivity tests of the assumptions made in the inverse system, for a mean carbon balance over the region of −183 Tg C, slightly weaker than the reference. Because of the potential large bias (∼24 Tg C in this case) due to choice of background conditions, proportional to the surface but not to the regional flux, this methodology seems limited to regions with a large signal (sink or source), unless additional observations can be used to constrain the boundary inflow.

Abstract. Megacities are major sources of anthropogenic fossil fuel CO 2 (FFCO 2 ) emissions. The spatial extents of these large urban systems cover areas of 10 000 km 2 or more with complex topography and changing landscapes. We present a high-resolution land-atmosphere modelling system for urban CO 2 emissions over the Los Angeles (LA) megacity area. The Weather Research and Forecasting (WRF)-Chem model was coupled to a very high-resolution FFCO 2 emission product, Hestia-LA, to simulate atmospheric CO 2 concentrations across the LA megacity at spatial resolutions as fine as ∼ 1 km. We evaluated multiple WRF configurations, selecting one that minimized errors in wind speed, wind direction, and boundary layer height as evaluated by its performance against meteorological data collected during the CalNex-LA campaign (May-June 2010). Our results show no significant difference between moderate-resolution (4 km) and high-resolution (1.3 km) simulations when evaluated against surface meteorological data, but the highresolution configurations better resolved planetary boundary layer heights and vertical gradients in the horizontal mean winds. We coupled our WRF configuration with the Vulcan 2.2 (10 km resolution) and Hestia-LA (1.3 km resolution) fossil fuel CO 2 emission products to evaluate the impact of the spatial resolution of the CO 2 emission products and the meteorological transport model on the representation of spatiotemporal variability in simulated atmospheric CO 2 concentrations. We find that high spatial resolution in the fossil fuel CO 2 emissions is more important than in the atmospheric model to capture CO 2 concentration variability across the LA megacity. Finally, we present a novel approach that employs simultaneous correlations of the simulated atmospheric CO 2 fields to qualitatively evaluate the greenhouse gas measurement network over the LA megacity. Spatial correlations in the atmospheric CO 2 fields reflect the coverage of individual measurement sites when a statistically significant number of sites observe emissions from a specific source or location. We conclude that elevated atmospheric CO 2 concentrations over the LA megacity are composed of multiple fine-scale plumes rather than a single homogenous urban dome. Furthermore, we conclude that FFCO 2 emissions monitoring in the LA megacity requires FFCO 2 emissions modelling with ∼ 1 km resolution becausePublished by Copernicus Publications on behalf of the European Geosciences Union. 9020 S. Feng et al.: LA megacity GHG modelling system coarser-resolution emissions modelling tends to overestimate the observational constraints on the emissions estimates.

This paper describes process-based estimation of CH4 emissions from sources in Indianapolis, IN and compares these with atmospheric inferences of whole city emissions. Emissions from the natural gas distribution system were estimated from measurements at metering and regulating stations and from pipeline leaks. Tracer methods and inverse plume modeling were used to estimate emissions from the major landfill and wastewater treatment plant. These direct source measurements informed the compilation of a methane emission inventory for the city equal to 29 Gg/yr (5% to 95% confidence limits, 15 to 54 Gg/yr). Emission estimates for the whole city based on an aircraft mass balance method and from inverse modeling of CH4 tower observations were 41 ± 12 Gg/yr and 81 ± 11 Gg/yr, respectively. Footprint modeling using 11 days of ethane/methane tower data indicated that landfills, wastewater treatment, wetlands, and other biological sources contribute 48% while natural gas usage and other fossil fuel sources contribute 52% of the city total. With the biogenic CH4 emissions omitted, the top-down estimates are 3.5-6.9 times the nonbiogenic city inventory. Mobile mapping of CH4 concentrations showed low level enhancement of CH4 throughout the city reflecting diffuse natural gas leakage and downstream usage as possible sources for the missing residual in the inventory.