We present a high-resolution atmospheric inversion system combining a Lagrangian Particle Dispersion Model (LPDM) and the Weather Research and Forecasting model (WRF), and test the impact of assimilating meteorological observation on transport accuracy. A Four Dimensional Data Assimilation (FDDA) technique continuously assimilates meteorological observations from various observing systems into the transport modeling system, and is coupled to the high resolution CO₂ emission product Hestia to simulate the atmospheric mole fractions of CO₂. For the Indianapolis Flux Experiment (INFLUX) project, we evaluated the impact of assimilating different meteorological observation systems on the linearized adjoint solutions and the CO₂ inverse fluxes estimated using observed CO₂ mole fractions from 11 out of 12 communications towers over Indianapolis for the Sep.-Nov. 2013 period. While assimilating WMO surface measurements improved the simulated wind speed and direction, their impact on the planetary boundary layer (PBL) was limited. Simulated PBL wind statistics improved significantly when assimilating upper-air observations from the commercial airline program Aircraft Communications Addressing and Reporting System (ACARS) and continuous ground-based Doppler lidar wind observations. Wind direction mean absolute error (MAE) decreased from 26 to 14 degrees and the wind speed MAE decreased from 2.0 to 1.2 m s_-1, while the bias remains small in all configurations (< 6 degrees and 0.2 m s_-1). Wind speed MAE and ME are larger in daytime than in nighttime. PBL depth MAE is reduced by ~10%, with little bias reduction. The inverse results indicate that the spatial distribution of CO₂ inverse fluxes were affected by the model performance while the overall flux estimates changed little across WRF simulations when aggregated over the entire domain. Our results show that PBL wind observations are a potent tool for increasing the precision of urban meteorological reanalyses, but that the impact on inverse flux estimates is dependent on the specific urban environment.

The INFLUX experiment has taken multiple approaches to estimate the carbon dioxide (CO₂) flux in a domain centered on the city of Indianapolis, Indiana. One approach, Hestia, uses a bottom-up technique relying on a mixture of activity data, fuel statistics, direct flux measurement and modeling algorithms. A second uses a Bayesian atmospheric inverse approach constrained by atmospheric CO₂ measurements and the Hestia emissions estimate as a prior CO₂ flux. The difference in the central estimate of the two approaches comes to 0.94 MtC (an 18.7% difference) over the eight-month period between September 1, 2012 and April 30, 2013, a statistically significant difference at the 2-sigma level. Here we explore possible explanations for this apparent discrepancy in an attempt to reconcile the flux estimates. We focus on two broad categories: 1) biases in the largest of bottom-up flux contributions and 2) missing CO₂ sources. Though there is some evidence for small biases in the Hestia fossil fuel carbon dioxide (FFCO2) flux estimate as an explanation for the calculated difference, we find more support for missing CO₂ fluxes, with biological respiration the largest of these. Incorporation of these differences bring the Hestia bottom-up and the INFLUX inversion flux estimates into statistical agreement and are additionally consistent with wintertime measurements of atmospheric ¹⁴CO₂. We conclude that comparison of bottom-up and top-down approaches must consider all flux contributions and highlight the important contribution to urban carbon budgets of animal and biotic respiration. Incorporation of missing CO₂ fluxes reconciles the bottom-up and inverse-based approach in the INFLUX domain.

We assess the detectability of city emissions via a tower-based greenhouse gas (GHG) network, as part of the Indianapolis Flux (INFLUX) experiment. By examining afternoon-averaged results from a network of carbon dioxide (CO₂), methane (CH₄), and carbon monoxide (CO) mole fraction measurements in Indianapolis, Indiana for 2011–2013, we quantify spatial and temporal patterns in urban atmospheric GHG dry mole fractions. The platform for these measurements is twelve communications towers spread across the metropolitan region, ranging in height from 39 to 136 m above ground level, and instrumented with cavity ring-down spectrometers. Nine of the sites were deployed as of January 2013 and data from these sites are the focus of this paper. A background site, chosen such that it is on the predominantly upwind side of the city, is utilized to quantify enhancements caused by urban emissions. Afternoon averaged mole fractions are studied because this is the time of day during which the height of the boundary layer is most steady in time and the area that influences the tower measurements is likely to be largest. Additionally, atmospheric transport models have better performance in simulating the daytime convective boundary layer compared to the nighttime boundary layer. Averaged from January through April of 2013, the mean urban dormant-season enhancements range from 0.3 ppm CO₂ at the site 24 km typically downwind of the edge of the city (Site 09) to 1.4 ppm at the site at the downwind edge of the city (Site 02) to 2.9 ppm at the downtown site (Site 03). When the wind is aligned such that the sites are downwind of the urban area, the enhancements are increased, to 1.6 ppm at Site 09, and 3.3 ppm at Site 02. Differences in sampling height affect the reported urban enhancement by up to 50%, but the overall spatial pattern remains similar. The time interval over which the afternoon data are averaged alters the calculated urban enhancement by an average of 0.4 ppm. The CO₂ observations are compared to CO₂ mole fractions simulated using a mesoscale atmospheric model and an emissions inventory for Indianapolis. The observed and modeled CO₂ enhancements are highly correlated (r² = 0.94), but the modeled enhancements prior to inversion average 53% of those measured at the towers. Following the inversion, the enhancements follow the observations closely, as expected. The CH₄ urban enhancement ranges from 5 ppb at the site 10 km predominantly downwind of the city (Site 13) to 21 ppb at the site near the landfill (Site 10), and for CO ranges from 6 ppb at the site 24 km downwind of the edge of the city (Site 09) to 29 ppb at the downtown site (Site 03). Overall, these observations show that a dense network of urban GHG measurements yield a detectable urban signal, well-suited as input to an urban inversion system given appropriate attention to sampling time, sampling altitude and quantification of background conditions.

Quantifying greenhouse gas (GHG) emissions from cities is a key challenge towards effective emissions management. An inversion analysis from the INdianapolis FLUX experiment (INFLUX) project, as the first of its kind, has achieved a top-down emission estimate for a single city using CO₂ data collected by the dense tower network deployed across the city. However, city-level emission data, used as a priori emissions, are also a key component in the atmospheric inversion framework. Currently, fine-grained emission inventories (EIs) able to resolve GHG city emissions at high spatial resolution, are only available for few major cities across the globe. Following the INFLUX inversion case with a global 1.1 km ODIAC fossil fuel CO₂ emission dataset, we further improved the ODIAC emission field and examined its utility as a prior for the city scale inversion. We disaggregated the 1.1 km ODIAC non-point source emissions using geospatial datasets such as the global road network data and satellite-data driven surface imperviousness data to a 30 . 30 m resolution. We assessed the impact of the improved emission field on the inversion result, relative to priors in previous studies (Hestia and ODIAC). The posterior total emission estimate (5.1 MtC/yr) remains statistically similar to the previous estimate with ODIAC (5.3 MtC/yr). However, the distribution of the flux corrections was very close to those of Hestia inversion and the model-observation mismatches were significantly reduced both in forward and inverse runs, even without hourly temporal changes in emissions. EIs reported by cities often do not have estimates of spatial extents. Thus, emission disaggregation is a required step when verifying those reported emissions using atmospheric models. Our approach offers gridded emission estimates for global cities that could serves as a prior for inversion, even without locally reported EIs in a systematic way to support city-level Measuring, Reporting and Verification (MRV) practice implementation.

The objective of the Indianapolis Flux Experiment (INFLUX) is to develop, evaluate and improve methods for measuring greenhouse gas (GHG) emissions from cities. INFLUX’s scientific objectives are to quantify CO₂ and CH₄ emission rates at 1 km₂ resolution with a 10% or better accuracy and precision, to determine whole-city emissions with similar skill, and to achieve high (weekly or finer) temporal resolution at both spatial resolutions. The experiment employs atmospheric GHG measurements from both towers and aircraft, atmospheric transport observations and models, and activity-based inventory products to quantify urban GHG emissions. Multiple, independent methods for estimating urban emissions are a central facet of our experimental design. INFLUX was initiated in 2010 and measurements and analyses are ongoing. To date we have quantified urban atmospheric GHG enhancements using aircraft and towers with measurements collected over multiple years, and have estimated whole-city CO₂ and CH₄ emissions using aircraft and tower GHG measurements, and inventory methods. Significant differences exist across methods; these differences have not yet been resolved; research to reduce uncertainties and reconcile these differences is underway. Sectorally- and spatially-resolved flux estimates, and detection of changes of fluxes over time, are also active research topics. Major challenges include developing methods for distinguishing anthropogenic from biogenic CO₂ fluxes, improving our ability to interpret atmospheric GHG measurements close to urban GHG sources and across a broader range of atmospheric stability conditions, and quantifying uncertainties in inventory data products. INFLUX data and tools are intended to serve as an open resource and test bed for future investigations. Well-documented, public archival of data and methods is under development in support of this objective.

Megacities are major sources of anthropogenic fossil fuel CO₂(FFCO₂) emissions. The spatial extents of these large urban systems cover areas of 10 000 km² or more with complex topography and changing landscapes. We present a high-resolution land–atmosphere modelling system for urban CO₂ emissions over the Los Angeles (LA) megacity area. The Weather Research and Forecasting (WRF)-Chem model was coupled to a very high-resolution FFCO₂ emission product, Hestia-LA, to simulate atmospheric CO₂ 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 high-resolution 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₂ emission products to evaluate the impact of the spatial resolution of the CO₂ emission products and the meteorological transport model on the representation of spatiotemporal variability in simulated atmospheric CO₂ concentrations. We find that high spatial resolution in the fossil fuel CO₂ emissions is more important than in the atmospheric model to capture CO₂ concentration variability across the LA megacity. Finally, we present a novel approach that employs simultaneous correlations of the simulated atmospheric CO₂ fields to qualitatively evaluate the greenhouse gas measurement network over the LA megacity. Spatial correlations in the atmospheric CO₂ 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₂ concentrations over the LA megacity are composed of multiple fine-scale plumes rather than a single homogenous urban dome. Furthermore, we conclude that FFCO₂ emissions monitoring in the LA megacity requires FFCO₂ emissions modelling with  ∼1 km resolution because coarser-resolution emissions modelling tends to overestimate the observational constraints on the emissions estimates.