Faculty and Staff

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.

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.

Atmospheric CO₂ inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO₂ measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr^-1 (±0.5 Pg C yr^-1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr^-1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr^-1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr[superscript −1] for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr^-1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr^-1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr^-1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.

This paper presents an analysis of methane emissions from the Los Angeles Basin at monthly timescales across a 4-year time period – from September 2011 to August 2015. Using observations acquired by a ground-based near-infrared remote sensing instrument on Mount Wilson, California, combined with atmospheric CH₄–CO₂ tracer–tracer correlations, we observed −18 to +22 % monthly variability in CH₄ : CO₂ from the annual mean in the Los Angeles Basin. Top-down estimates of methane emissions for the basin also exhibit significant monthly variability (−19 to +31 % from annual mean and a maximum month-to-month change of 47 %). During this period, methane emissions consistently peaked in the late summer/early fall and winter. The estimated annual methane emissions did not show a statistically significant trend over the 2011 to 2015 time period.

Oppia Quadrangle Av-10 (288–360°E, ±22°) is a junction of key geologic features that preserve a rough history of Asteroid (4) Vesta and serves as a case study of using geologic mapping to define a relative geologic timescale. Clear filter images, stereo-derived topography, slope maps, and multispectral color-ratio images from the Framing Camera on NASA’s Dawn spacecraft served as basemaps to create a geologic map and investigate the spatial and temporal relationships of the local stratigraphy. Geologic mapping reveals the oldest map unit within Av-10 is the cratered highlands terrain which possibly represents original crustal material on Vesta that was then excavated by one or more impacts to form the basin Feralia Planitia. Saturnalia Fossae and Divalia Fossae ridge and trough terrains intersect the wall of Feralia Planitia indicating that this impact basin is older than both the Veneneia and Rheasilvia impact structures, representing Pre-Veneneian crustal material. Two of the youngest geologic features in Av-10 are Lepida (∼45 km diameter) and Oppia (∼40 km diameter) impact craters that formed on the northern and southern wall of Feralia Planitia and each cross-cuts a trough terrain. The ejecta blanket of Oppia is mapped as ‘dark mantle’ material because it appears dark orange in the Framing Camera ‘Clementine-type’ color-ratio image and has a diffuse, gradational contact distributed to the south across the rim of Rheasilvia. Mapping of surface material that appears light orange in color in the Framing Camera ‘Clementine-type’ color-ratio image as ‘light mantle material’ supports previous interpretations of an impact ejecta origin. Some light mantle deposits are easily traced to nearby source craters, but other deposits may represent distal ejecta deposits (emplaced >5 crater radii away) in a microgravity environment.

We produced two 1:250,000 scale geologic maps of the adjacent quadrangles Av-6 Gegania and Av-7 Lucaria, located in the equatorial region of (4) Vesta (0–144°E, 22°S to 22°N). The mapping is based on clear and color filter images of the Framing Camera (FC) onboard the Dawn spacecraft, which has captured the entire illuminated surface of Vesta with high spatial resolution (up to ∼20 m/pixel), and on a digital terrain model derived from FC imagery. Besides the geologic mapping itself, a secondary purpose of this work is to investigate one of the most prominent morphological features on Vesta, namely the aggregation of several giant equatorial troughs termed the Divalia Fossae, most probably formed during the Rheasilvia impact near Vesta’s south pole. The up to 465 km long and 22 km wide troughs show height differences of up to 5 km between adjacent troughs and ridges. Another imprint of the Rheasilvia impact is the >350 km long and ∼250 km wide swath of ejecta crossing quadrangle Av-6 Gegania. This lobe shows a distinct appearance in FC color ratios and a high albedo in FC images, indicating a mineralogical similarity to material typically found within the Rheasilvia basin, in particular composed of diogenite-rich howardites. Almost the entire northern half of the mapping area shows the oldest surface, being dominated by upper crustal basaltic material. To the south, increasingly younger formations related to the Rheasilvia impact occur, either indicated by the troughs formed by Rheasilvia or by the Rheasilvia ejecta itself. Only medium sized impact craters with diameters less than 22 km occur within the two mapped quadrangles. Some of the craters exhibit ejecta blankets and/or distinctly dark or bright ejecta material in ejecta rays outside and exposures within the crater, and mass-wasting deposits down crater slopes, forming the youngest surfaces.

We used Dawn spacecraft data to identify and delineate geological units and landforms in the Marcia quadrangle of Vesta as a means to assess the role of the large, relatively young impact craters Marcia (∼63 km diam.) and Calpurnia (∼53 km diam.) and their surrounding ejecta field on the local geology. We also investigated a local topographic high with a dark-rayed crater named Aricia Tholus, and the impact crater Octavia that is surrounded by a distinctive diffuse mantle. Crater counts and stratigraphic relations suggest that Marcia is the youngest large crater on Vesta, in which a putative impact melt on the crater floor ranges in age between ∼40 and 60 Ma (depending upon choice of chronology system), and Marcia’s ejecta blanket ranges in age between ∼120 and 390 Ma (depending upon choice of chronology system).

We interpret the geologic units in and around Marcia crater to mark a major vestan time-stratigraphic event, and that the Marcia Formation is one of the geologically youngest formations on Vesta. Marcia crater reveals pristine bright and dark material in its walls and smooth and pitted terrains on its floor. The smooth unit we interpret as evidence of flow of impact melts and (for the pitted terrain) release of volatiles during or after the impact process. The distinctive dark ejecta surrounding craters Marcia and Calpurnia is enriched in OH- or H-bearing phases and has a variable morphology, suggestive of a complex mixture of impact ejecta and impact melts including dark materials possibly derived from carbonaceous chondrite-rich material. Aricia Tholus, which was originally interpreted as a putative vestan volcanic edifice based on lower resolution observations, appears to be a fragment of an ancient impact basin rim topped by a dark-rayed impact crater. Octavia crater has a cratering model formation age of ∼280–990 Ma based on counts of its ejecta field (depending upon choice of chronology system), and its ejecta field is the second oldest unit in this quadrangle. The relatively young craters and their related ejecta materials in this quadrangle are in stark contrast to the surrounding heavily cratered units that are related to the billion years old or older Rheasilvia and Veneneia impact basins and Vesta’s ancient crust preserved on Vestalia Terra.

The purpose of this paper is to introduce the Geologic Mapping of Vesta Special Issue/Section of Icarus, which includes several papers containing geologic maps of the surface of Vesta made to support data analysis conducted by the Dawn Science Team during the Vesta Encounter (July 2011–September 2012). In this paper we briefly discuss pre-Dawn knowledge of Vesta, provide the goals of our geologic mapping campaign, discuss the methodologies and materials used for geologic mapping, review the global geologic context of Vesta, discuss the challenges of mapping the geology of Vesta as a small airless body, and describe the content of the papers in this Special Issue/Section. We conclude with a discussion of lessons learned from our quadrangle-based mapping effort and provide recommendations for conducting mapping campaigns as part of planetary spacecraft nominal missions.

Attributing observed CO2 variations to human or natural cause is critical to deducing and tracking emissions from observations. We have used in situ CO2, CO, and planetary boundary layer height (PBLH) measurements recorded during the CalNex-LA (CARB et al., 2008) ground campaign of 15 May-15 June 2010, in Pasadena, CA, to deduce the diurnally varying anthropogenic component of observed CO2 in the megacity of Los Angeles (LA). This affordable and simple technique, validated by carbon isotope observations and WRF-STILT (Weather Research and Forecasting model - Stochastic Time-Inverted Lagrangian Transport model) predictions, is shown to robustly attribute observed CO2 variation to anthropogenic or biogenic origin over the entire diurnal cycle. During CalNex-LA, local fossil fuel combustion contributed up to similar to 50% of the observed CO2 enhancement overnight, and similar to 100% of the enhancement near midday. This suggests that sufficiently accurate total column CO2 observations recorded near midday, such as those from the GOSAT or OCO-2 satellites, can potentially be used to track anthropogenic emissions from the LA megacity.

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.