ASU Scholarship Showcase

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.

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.

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.

Atmospheric radiocarbon (¹⁴C) represents an important observational constraint on emissions of fossil-fuel derived carbon into the atmosphere due to the absence of ¹⁴C in fossil fuel reservoirs. The high sensitivity and precision that accelerator mass spectrometry (AMS) affords in atmospheric ¹⁴C analysis has greatly increased the potential for using such measurements to evaluate bottom-up emissions inventories of fossil fuel CO₂(CO₂ff), as well as those for other co-emitted species. Here we use observations of ¹⁴CO₂ and a series of primary hydrocarbons and combustion tracers from discrete air samples collected between June 2009 and September 2010 at the National Oceanic and Atmospheric Administration Boulder Atmospheric Observatory (BAO; Lat: 40.050° N, Lon: 105.004° W) to derive emission ratios of each species with respect to CO₂ff. The BAO tower is situated at the boundary of the Denver metropolitan area to the south and a large industrial and agricultural region to the north and east, making it an ideal location to study the contrasting mix of emissions from the activities in each region. The species considered in this analysis are carbon monoxide (CO), methane (CH₄), acetylene (C₂H₂), benzene (C₆H₆), and C₃–C₅ alkanes. We estimate emissions for a subset of these species by using the Vulcan high resolution CO2ff emission data product as a reference. We find that CO is overestimated in the 2008 National Emissions Inventory (NEI08) by a factor of ~2. A close evaluation of the inventory suggests that the ratio of CO emitted per unit fuel burned from on-road gasoline vehicles is likely over-estimated by a factor of 2.5. Using a wind-directional analysis of the data, we find enhanced concentrations of CH₄, relative to CO₂ff, in air influenced by emissions to the north and east of the BAO tower when compared to air influenced by emissions in the Denver metro region to the south. Along with enhanced CH₄, the strongest enhancements of the C₃–C₅ alkanes are also found in the north and east wind sector, suggesting that both the alkane and CH₄ enhancements are sourced from oil and gas fields located to the northeast, though it was not possible to rule out the contribution of non oil and gas CH₄ sources.

Urban environments are the primary contributors to global anthropogenic carbon emissions. Because much of the growth in CO₂ 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₂ and CH₄, 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 emissions 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.

Large urban emissions of greenhouse gases result in large atmospheric enhancements relative to background that are easily measured. Using CO₂ mole fractions and Δ¹⁴C and δ¹³C values of CO₂ in the Los Angeles megacity observed in inland Pasadena (2006–2013) and coastal Palos Verdes peninsula (autumn 2009–2013), we have determined time series for CO₂ contributions from fossil fuel combustion (C_ff) for both sites and broken those down into contributions from petroleum and/or gasoline and natural gas burning for Pasadena. We find a 10 % reduction in Pasadena C_ff during the Great Recession of 2008–2010, which is consistent with the bottom-up inventory determined by the California Air Resources Board. The isotopic variations and total atmospheric CO₂ from our observations are used to infer seasonality of natural gas and petroleum combustion. The trend of CO₂ contributions to the atmosphere from natural gas combustion is out of phase with the seasonal cycle of total natural gas combustion seasonal patterns in bottom-up inventories but is consistent with the seasonality of natural gas usage by the area's electricity generating power plants. For petroleum, the inferred seasonality of CO₂ contributions from burning petroleum is delayed by several months relative to usage indicated by statewide gasoline taxes. Using the high-resolution Hestia-LA data product to compare C_ff from parts of the basin sampled by winds at different times of year, we find that variations in observed fossil fuel CO₂ reflect seasonal variations in wind direction. The seasonality of the local CO₂ excess from fossil fuel combustion along the coast, on Palos Verdes peninsula, is higher in autumn and winter than spring and summer, almost completely out of phase with that from Pasadena, also because of the annual variations of winds in the region. Variations in fossil fuel CO₂ signals are consistent with sampling the bottom-up Hestia-LA fossil CO₂ emissions product for sub-city source regions in the LA megacity domain when wind directions are considered.

Streamer biofilm communities (SBC) are often observed within chemosynthetic zones of Yellowstone hot spring outflow channels, where temperatures exceed those conducive to photosynthesis. Nearest the hydrothermal source (75–88°C) SBC comprise thermophilic Archaea and Bacteria, often mixed communities including Desulfurococcales and uncultured Crenarchaeota, as well as Aquificae and Thermus, each carrying diagnostic membrane lipid biomarkers. We tested the hypothesis that SBC can alternate their metabolism between autotrophy and heterotrophy depending on substrate availability. Feeding experiments were performed at two alkaline hot springs in Yellowstone National Park: Octopus Spring and “Bison Pool,” using various ¹³C-labeled substrates (bicarbonate, formate, acetate, and glucose) to determine the relative uptake of these different carbon sources. Highest ¹³C uptake, at both sites, was from acetate into almost all bacterial fatty acids, particularly into methyl-branched C₁₅, C₁₇ and C₁₉ fatty acids that are diagnostic for Thermus/Meiothermus, and some Firmicutes as well as into universally common C_16:0 and C_18:0 fatty acids. ¹³C-glucose showed a similar, but a 10–30 times lower uptake across most fatty acids. ¹³C-bicarbonate uptake, signifying the presence of autotrophic communities was only significant at “Bison Pool” and was observed predominantly in non-specific saturated C₁₆, C₁₈, C₂₀, and C₂₂ fatty acids. Incorporation of ¹³C-formate occurred only at very low rates at “Bison Pool” and was almost undetectable at Octopus Spring, suggesting that formate is not an important carbon source for SBC. ¹³C-uptake into archaeal lipids occurred predominantly with ¹³C-acetate, suggesting also that archaeal communities at both springs have primarily heterotrophic carbon assimilation pathways. We hypothesize that these communities are energy-limited and predominantly nurtured by input of exogenous organic material, with only a small fraction being sustained by autotrophic growth.

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.

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.