ASU Scholarship Showcase

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.

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.

High-resolution, global quantification of fossil fuel CO[subscript 2] emissions is emerging as a critical need in carbon cycle science and climate policy. We build upon a previously developed fossil fuel data assimilation system (FFDAS) for estimating global high-resolution fossil fuel CO[subscript 2] emissions. We have improved the underlying observationally based data sources, expanded the approach through treatment of separate emitting sectors including a new pointwise database of global power plants, and extended the results to cover a 1997 to 2010 time series at a spatial resolution of 0.1°. Long-term trend analysis of the resulting global emissions shows subnational spatial structure in large active economies such as the United States, China, and India. These three countries, in particular, show different long-term trends and exploration of the trends in nighttime lights, and population reveal a decoupling of population and emissions at the subnational level. Analysis of shorter-term variations reveals the impact of the 2008–2009 global financial crisis with widespread negative emission anomalies across the U.S. and Europe. We have used a center of mass (CM) calculation as a compact metric to express the time evolution of spatial patterns in fossil fuel CO[subscript 2] emissions. The global emission CM has moved toward the east and somewhat south between 1997 and 2010, driven by the increase in emissions in China and South Asia over this time period. Analysis at the level of individual countries reveals per capita CO[subscript 2] emission migration in both Russia and India. The per capita emission CM holds potential as a way to succinctly analyze subnational shifts in carbon intensity over time. Uncertainties are generally lower than the previous version of FFDAS due mainly to an improved nightlight data set.

The sensitivity of Earth’s wetlands to observed shifts in global precipitation and temperature patterns and their ability to produce large quantities of methane gas are key global change questions. We present a microwave satellite-based approach for mapping fractional surface water (FW) globally at 25-km resolution. The approach employs a land cover-supported, atmospherically-corrected dynamic mixture model applied to 20+ years (1992–2013) of combined, daily, passive/active microwave remote sensing data. The resulting product, known as Surface Water Microwave Product Series (SWAMPS), shows strong microwave sensitivity to sub-grid scale open water and inundated wetlands comprising open plant canopies. SWAMPS’ FW compares favorably (R² = 91%–94%) with higher-resolution, global-scale maps of open water from MODIS and SRTM-MOD44W. Correspondence of SWAMPS with open water and wetland products from satellite SAR in Alaska and the Amazon deteriorates when exposed wetlands or inundated forests captured by the SAR products were added to the open water fraction reflecting SWAMPS’ inability to detect water underneath the soil surface or beneath closed forest canopies. Except for a brief period of drying during the first 4 years of observation, the inundation extent for the global domain excluding the coast was largely stable. Regionally, inundation in North America is advancing while inundation is on the retreat in Tropical Africa and North Eurasia. SWAMPS provides a consistent and long-term global record of daily FW dynamics, with documented accuracies suitable for hydrologic assessment and global change-related investigations.

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.

A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO₂) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960–2009 at 0.5° resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO₂ fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data.

The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g C m^-2 yr^-2, equivalent to 3 to 340 % of the respective model means, over the analysis period. For the multimodel average the increase is 135 % of the mean from the first to last 10 years of record (1960–1969 vs. 2000–2009), with a weakening CO₂ sink over the latter decades. Vegetation net primary productivity increased by 8 to 30 % from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO₂ sink strength.

The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO₂ exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

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.

Soil temperature (T_s) change is a key indicator of the dynamics of permafrost. On seasonal and interannual timescales, the variability of T_s determines the active-layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T_s not only drives permafrost thaw/retreat but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of T_s in permafrost regions. There is a large spread of T_s trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 °C yr^-1. Most models show smaller increase in T_s with increasing depth. Air temperature (Tsub>a) and longwave downward radiation (LWDR) are the main drivers of T_s trends, but their relative contributions differ amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in T_s trends, while trends of T_a only explain 5 % of the differences in T_s trends. Uncertain climate forcing contributes a larger uncertainty in T_s trends (0.021 ± 0.008 °C yr^-1, mean ± standard deviation) than the uncertainty of model structure (0.012 ± 0.001 °C yr^-1), diagnosed from the range of response between different models, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active-layer thickness (ALT) is less than 3 m loss rate, is found to be significantly correlated with the magnitude of the trends of T_s at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T_s at 1 m is estimated to be of −2.80 ± 0.67 million km²°C^-1. Finally, by using two long-term LWDR data sets and relationships between trends of LWDR and T_s across models, we infer an observation-constrained total boreal near-surface permafrost area decrease comprising between 39 ± 14  ×  10³ and 75 ± 14  ×  10³km²yr^-1 from 1960 to 2000. This corresponds to 9–18 % degradation of the current permafrost area.

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.