ASU Scholarship Showcase

The ‘Hestia Project’ uses a bottom-up approach to quantify fossil fuel CO₂(FFCO₂) emissions spatially at the building/street level and temporally at the hourly level. Hestia FFCO₂ emissions are provided in the form of a group of sector-specific vector layers with point, line, and polygon sources to support carbon cycle science and climate policy. Application to carbon cycle science, in particular, requires regular gridded data in order to link surface carbon fluxes to atmospheric transport models. However, the heterogeneity and complexity of FFCO₂ sources within regular grids is sensitive to spatial resolution. From the perspective of a data provider, we need to find a balance between resolution and data volume so that the gridded data product retains the maximum amount of information content while maintaining an efficient data volume.

The Shannon entropy determines the minimum bits that are needed to encode an information source and can serve as a metric for the effective information content. In this paper, we present an analysis of the Shannon entropy of gridded FFCO₂ emissions with varying resolutions in four Hestia study areas, and find: (1) the Shannon entropy increases with smaller grid resolution until it reaches a maximum value (the max-entropy resolution); (2) total emissions (the sum of several sector-specific emission fields) show a finer max-entropy resolution than each of the sector-specific fields; (3) the residential emissions show a finer max-entropy resolution than the commercial emissions; (4) the max-entropy resolution of the onroad emissions grid is closely correlated to the density of the road network. These findings suggest that the Shannon entropy can detect the information effectiveness of the spatial resolution of gridded FFCO₂ emissions. Hence, the resolution-entropy relationship can be used to assist in determining an appropriate spatial resolution for urban CO₂ flux studies. We conclude that the optimal spatial resolution for providing Hestia total FFCO₂ emissions products is centered around 100 m, at which the FFCO₂ emissions data can not only fully meet the requirement of urban flux integration, but also be effectively used in understanding the relationships between FFCO₂ emissions and various social-economic variables at the U.S. census block group level.

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.

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.

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.

Recent advances in fossil fuel CO₂ (FFCO₂) emission inventories enable sensitivity tests of simulated atmospheric CO₂ concentrations to sub-annual variations in FFCO₂ emissions and what this implies for the interpretation of observed CO₂. Six experiments are conducted to investigate the potential impact of three cycles of FFCO₂ emission variability (diurnal, weekly and monthly) using a global tracer transport model. Results show an annual FFCO₂ rectification varying from −1.35 to +0.13 ppm from the combination of all three cycles. This rectification is driven by a large negative diurnal FFCO₂ rectification due to the covariation of diurnal FFCO₂ emissions and diurnal vertical mixing, as well as a smaller positive seasonal FFCO₂ rectification driven by the covariation of monthly FFCO₂ emissions and monthly atmospheric transport. The diurnal FFCO₂ emissions are responsible for a diurnal FFCO₂ concentration amplitude of up to 9.12 ppm at the grid cell scale. Similarly, the monthly FFCO₂ emissions are responsible for a simulated seasonal CO₂ amplitude of up to 6.11 ppm at the grid cell scale. The impact of the diurnal FFCO₂ emissions, when only sampled in the local afternoon, is also important, causing an increase of +1.13 ppmv at the grid cell scale. The simulated CO₂ concentration impacts from the diurnally and seasonally varying FFCO₂ emissions are centered over large source regions in the Northern Hemisphere, extending to downwind regions. This study demonstrates the influence of sub-annual variations in FFCO₂ emissions on simulated CO₂ concentration and suggests that inversion studies must take account of these variations in the affected regions.

Errors in the specification or utilization of fossil fuel CO₂ emissions within carbon budget or atmospheric CO₂ inverse studies can alias the estimation of biospheric and oceanic carbon exchange. A key component in the simulation of CO₂ concentrations arising from fossil fuel emissions is the spatial distribution of the emission near coastlines. Regridding of fossil fuel CO₂ emissions (FFCO₂) from fine to coarse grids to enable atmospheric transport simulations can give rise to mismatches between the emissions and simulated atmospheric dynamics which differ over land or water. For example, emissions originally emanating from the land are emitted from a grid cell for which the vertical mixing reflects the roughness and/or surface energy exchange of an ocean surface. We test this potential "dynamical inconsistency" by examining simulated global atmospheric CO₂ concentration driven by two different approaches to regridding fossil fuel CO₂ emissions. The two approaches are as follows: (1) a commonly used method that allocates emissions to grid cells with no attempt to ensure dynamical consistency with atmospheric transport and (2) an improved method that reallocates emissions to grid cells to ensure dynamically consistent results. Results show large spatial and temporal differences in the simulated CO₂ concentration when comparing these two approaches. The emissions difference ranges from −30.3 TgC grid cell^-1 yr^-1 (−3.39 kgC m^-2 yr^-1) to +30.0 TgC grid cell^-1 yr^-1 (+2.6 kgC m^-2 yr^-1) along coastal margins. Maximum simulated annual mean CO₂ concentration differences at the surface exceed ±6 ppm at various locations and times. Examination of the current CO₂ monitoring locations during the local afternoon, consistent with inversion modeling system sampling and measurement protocols, finds maximum hourly differences at 38 stations exceed ±0.10 ppm with individual station differences exceeding −32 ppm. The differences implied by not accounting for this dynamical consistency problem are largest at monitoring sites proximal to large coastal urban areas and point sources. These results suggest that studies comparing simulated to observed atmospheric CO₂ concentration, such as atmospheric CO₂ inversions, must take measures to correct for this potential problem and ensure flux and dynamical consistency.

A globally integrated carbon observation and analysis system is needed to improve the fundamental understanding of the global carbon cycle, to improve our ability to project future changes, and to verify the effectiveness of policies aiming to reduce greenhouse gas emissions and increase carbon sequestration. Building an integrated carbon observation system requires transformational advances from the existing sparse, exploratory framework towards a dense, robust, and sustained system in all components: anthropogenic emissions, the atmosphere, the ocean, and the terrestrial biosphere. The paper is addressed to scientists, policymakers, and funding agencies who need to have a global picture of the current state of the (diverse) carbon observations.

We identify the current state of carbon observations, and the needs and notional requirements for a global integrated carbon observation system that can be built in the next decade. A key conclusion is the substantial expansion of the ground-based observation networks required to reach the high spatial resolution for CO₂ and CH₄ fluxes, and for carbon stocks for addressing policy-relevant objectives, and attributing flux changes to underlying processes in each region. In order to establish flux and stock diagnostics over areas such as the southern oceans, tropical forests, and the Arctic, in situ observations will have to be complemented with remote-sensing measurements. Remote sensing offers the advantage of dense spatial coverage and frequent revisit. A key challenge is to bring remote-sensing measurements to a level of long-term consistency and accuracy so that they can be efficiently combined in models to reduce uncertainties, in synergy with ground-based data.

Bringing tight observational constraints on fossil fuel and land use change emissions will be the biggest challenge for deployment of a policy-relevant integrated carbon observation system. This will require in situ and remotely sensed data at much higher resolution and density than currently achieved for natural fluxes, although over a small land area (cities, industrial sites, power plants), as well as the inclusion of fossil fuel CO₂ proxy measurements such as radiocarbon in CO₂ and carbon-fuel combustion tracers. Additionally, a policy-relevant carbon monitoring system should also provide mechanisms for reconciling regional top-down (atmosphere-based) and bottom-up (surface-based) flux estimates across the range of spatial and temporal scales relevant to mitigation policies. In addition, uncertainties for each observation data-stream should be assessed. The success of the system will rely on long-term commitments to monitoring, on improved international collaboration to fill gaps in the current observations, on sustained efforts to improve access to the different data streams and make databases interoperable, and on the calibration of each component of the system to agreed-upon international scales.