ASU Scholarship Showcase

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.

Mild Cognitive Impairment (MCI) is a transitional stage between normal aging and dementia and people with MCI are at high risk of progression to dementia. MCI is attracting increasing attention, as it offers an opportunity to target the disease process during an early symptomatic stage. Structural magnetic resonance imaging (MRI) measures have been the mainstay of Alzheimer's disease (AD) imaging research, however, ventricular morphometry analysis remains challenging because of its complicated topological structure. Here we describe a novel ventricular morphometry system based on the hyperbolic Ricci flow method and tensor-based morphometry (TBM) statistics. Unlike prior ventricular surface parameterization methods, hyperbolic conformal parameterization is angle-preserving and does not have any singularities. Our system generates a one-to-one diffeomorphic mapping between ventricular surfaces with consistent boundary matching conditions. The TBM statistics encode a great deal of surface deformation information that could be inaccessible or overlooked by other methods. We applied our system to the baseline MRI scans of a set of MCI subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI: 71 MCI converters vs. 62 MCI stable). Although the combined ventricular area and volume features did not differ between the two groups, our fine-grained surface analysis revealed significant differences in the ventricular regions close to the temporal lobe and posterior cingulate, structures that are affected early in AD. Significant correlations were also detected between ventricular morphometry, neuropsychological measures, and a previously described imaging index based on fluorodeoxyglucose positron emission tomography (FDG-PET) scans. This novel ventricular morphometry method may offer a new and more sensitive approach to study preclinical and early symptomatic stage AD.

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.

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.

Alzheimer's disease (AD) is a progressive brain disease. Accurate detection of AD and its prodromal stage, mild cognitive impairment (MCI), are crucial. There is also a growing interest in identifying brain imaging biomarkers that help to automatically differentiate stages of Alzheimer's disease. Here, we focused on brain structural networks computed from diffusion MRI and proposed a new feature extraction and classification framework based on higher order singular value decomposition and sparse logistic regression. In tests on publicly available data from the Alzheimer's Disease Neuroimaging Initiative, our proposed framework showed promise in detecting brain network differences that help in classifying different stages of Alzheimer's disease.

Alzheimer’s disease (AD) involves a gradual breakdown of brain connectivity, and network analyses offer a promising new approach to track and understand disease progression. Even so, our ability to detect degenerative changes in brain networks depends on the methods used. Here we compared several tractography and feature extraction methods to see which ones gave best diagnostic classification for 202 people with AD, mild cognitive impairment or normal cognition, scanned with 41-gradient diffusion-weighted magnetic resonance imaging as part of the Alzheimer’s Disease Neuroimaging Initiative (ADNI) project. We computed brain networks based on whole brain tractography with nine different methods – four of them tensor-based deterministic (FACT, RK2, SL, and TL), two orientation distribution function (ODF)-based deterministic (FACT, RK2), two ODF-based probabilistic approaches (Hough and PICo), and one “ball-and-stick” approach (Probtrackx). Brain networks derived from different tractography algorithms did not differ in terms of classification performance on ADNI, but performing principal components analysis on networks helped classification in some cases. Small differences may still be detectable in a truly vast cohort, but these experiments help assess the relative advantages of different tractography algorithms, and different post-processing choices, when used for classification.

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.

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.

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.

Recent neuroimaging findings have highlighted the impact of premature birth on subcortical development and morphological changes in the deep grey nuclei and ventricular system. To help characterize subcortical microstructural changes in preterm neonates, we recently implemented a multivariate tensor-based method (mTBM). This method allows to precisely measure local surface deformation of brain structures in infants. Here, we investigated ventricular abnormalities and their spatial relationships with surrounding subcortical structures in preterm neonates. We performed regional group comparisons on the surface morphometry and relative position of the lateral ventricles between 19 full-term and 17 preterm born neonates at term-equivalent age. Furthermore, a relative pose analysis was used to detect individual differences in translation, rotation, and scale of a given brain structure with respect to an average. Our mTBM results revealed broad areas of alterations on the frontal horn and body of the left ventricle, and narrower areas of differences on the temporal horn of the right ventricle. A significant shift in the rotation of the left ventricle was also found in preterm neonates. Furthermore, we located significant correlations between morphology and pose parameters of the lateral ventricles and that of the putamen and thalamus. These results show that regional abnormalities on the surface and pose of the ventricles are also associated with alterations on the putamen and thalamus. The complementarity of the information provided by the surface and pose analysis may help to identify abnormal white and grey matter growth, hinting toward a pattern of neural and cellular dysmaturation.