Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Set”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.

benchmark↗

Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways and Discovering Blind Geothermal Systems in the Great Basin Region, Western USA: A Progress Report

We have undertaken an integrated geologic, geochemical, and geophysical study of a broad 240-km-wide, 400-km-long transect stretching from west-central to eastern Nevada in the Great Basin region of the western USA. The main goal of this study is to produce a comprehensive geothermal potential map that incorporates up to 11 parameters and identifies geothermal play fairways that represent potential blind or hidden geothermal systems. Our new geothermal potential map incorporates: 1) heat flow; 2) geochemistry from springs and wells; 3) structural setting; 4) recency of faulting; 5) slip rates on Quaternary faults; 6) regional strain rate; 7) slip and dilation tendency on Quaternary faults; 8) seismologic data; 9) gravity data; 10) magnetotelluric data (where available); and 11) seismic reflection data (primarily from the Carson Sink and Steptoe basins). The transect is respectively anchored on its western and eastern ends by regional 3D modeling of the Carson Sink and Steptoe basins, which will provide more detailed geothermal potential maps of these two promising areas. To date, geological, geochemical, and geophysical data sets have been assembled into an ArcGIS platform and combined into a preliminary predictive geothermal play fairway model using various statistical techniques. The fairway model consists of the following components, each of which are represented in grid-cell format in ArcGIS and combined using specified weights and mathematical operators: 1) structural component of permeability; 2) regional-scale component of permeability; 3) combined permeability, and 4) heat source model. The preliminary model demonstrates that the multiple data sets can be successfully combined into a comprehensive favorability map. An initial evaluation using known geothermal systems as benchmarks to test interpretations indicates that the preliminary modeling has done a good job assigning relative ranks of geothermal potential. However, a major challenge is defining logical relative rankings of each parameter and how best to combine the multiple data sets into the geothermal potential/ permeability map. Ongoing feedback and data analysis are in use to revise the grouping and weighting of some parameters in order to develop a more robust, optimized, final model. The final product will incorporate more parameters into a geothermal potential map than any previous effort in the region and may serve as a prototype to develop comprehensive geothermal potential maps for other regions.

15 GEOTHERMAL ENERGY↗

Internal respiratory surrogate in multislice 4D CT using a combination of Fourier transform and anatomical features

Purpose: The purpose of this study was to develop a novel algorithm to create a robust internal respiratory signal (IRS) for retrospective sorting of four-dimensional (4D) computed tomography (CT) images. Methods: The proposed algorithm combines information from the Fourier transform of the CT images and from internal anatomical features to form the IRS. The algorithm first extracts potential respiratory signals from low-frequency components in the Fourier space and selected anatomical features in the image space. A clustering algorithm then constructs groups of potential respiratory signals with similar temporal oscillation patterns. The clustered group with the largest number of similar signals is chosen to form the final IRS. To evaluate the performance of the proposed algorithm, the IRS was computed and compared with the external respiratory signal from the real-time position management (RPM) system on 80 patients. Results: In 72 (90%) of the 4D CT data sets tested, the IRS computed by the authors’ proposed algorithm matched with the RPM signal based on their normalized cross correlation. For these data sets with matching respiratory signals, the average difference between the end inspiration times (Δt{sub ins}) in the IRS and RPM signal was 0.11 s, and only 2.1% of Δt{sub ins} were more than 0.5 s apart. In the eight (10%) 4D CT data sets in which the IRS and the RPM signal did not match, the average Δt{sub ins} was 0.73 s in the nonmatching couch positions, and 35.4% of them had a Δt{sub ins} greater than 0.5 s. At couch positions in which IRS did not match the RPM signal, a correlation-based metric indicated poorer matching of neighboring couch positions in the RPM-sorted images. This implied that, when IRS did not match the RPM signal, the images sorted using the IRS showed fewer artifacts than the clinical images sorted using the RPM signal. Conclusions: The authors’ proposed algorithm can generate robust IRSs that can be used for retrospective sorting of 4D CT data. The algorithm is completely automatic and requires very little processing time. The algorithm is cost efficient and can be easily adopted for everyday clinical use.

60 APPLIED LIFE SCIENCES↗

A cross-study analysis of drug response prediction in cancer cell lines

Abstract To enable personalized cancer treatment, machine learning models have been developed to predict drug response as a function of tumor and drug features. However, most algorithm development efforts have relied on cross-validation within a single study to assess model accuracy. While an essential first step, cross-validation within a biological data set typically provides an overly optimistic estimate of the prediction performance on independent test sets. To provide a more rigorous assessment of model generalizability between different studies, we use machine learning to analyze five publicly available cell line-based data sets: National Cancer Institute 60, ancer Therapeutics Response Portal (CTRP), Genomics of Drug Sensitivity in Cancer, Cancer Cell Line Encyclopedia and Genentech Cell Line Screening Initiative (gCSI). Based on observed experimental variability across studies, we explore estimates of prediction upper bounds. We report performance results of a variety of machine learning models, with a multitasking deep neural network achieving the best cross-study generalizability. By multiple measures, models trained on CTRP yield the most accurate predictions on the remaining testing data, and gCSI is the most predictable among the cell line data sets included in this study. With these experiments and further simulations on partial data, two lessons emerge: (1) differences in viability assays can limit model generalizability across studies and (2) drug diversity, more than tumor diversity, is crucial for raising model generalizability in preclinical screening.

59 BASIC BIOLOGICAL SCIENCES↗

A Bayesian Approach to Period Searching in Solar Coronal Loops

We have applied a Bayesian generalized Lomb–Scargle period searching algorithm to movies of coronal loop images obtained with the Hinode X-ray Telescope (XRT) to search for evidence of periodicities that would indicate resonant heating of the loops. The algorithm makes as its only assumption that there is a single sinusoidal signal within each light curve of the data. Both the amplitudes and noise are taken as free parameters. It is argued that this procedure should be used alongside Fourier and wavelet analyses to more accurately extract periodic intensity modulations in coronal loops. The data analyzed are from XRT Observation Program 129C: “MHD Wave Heating (Thin Filters),” which occurred during 2006 November 13 and focused on active region 10293, which included coronal loops. The first data set spans approximately 10 min with an average cadence of 2 s, 2″ per pixel resolution, and used the Al-mesh analysis filter. The second data set spans approximately 4 min with a 3 s average cadence, 1″ per pixel resolution, and used the Al-poly analysis filter. The final data set spans approximately 22 min at a 6 s average cadence, and used the Al-poly analysis filter. In total, 55 periods of sinusoidal coronal loop oscillations between 5.5 and 59.6 s are discussed, supporting proposals in the literature that resonant absorption of magnetic waves is a viable mechanism for depositing energy in the corona.

79 ASTRONOMY AND ASTROPHYSICS↗

MetaboTools: A comprehensive toolbox for analysis of genome-scale metabolic models

Metabolomic data sets provide a direct read-out of cellular phenotypes and are increasingly generated to study biological questions. Previous work, by us and others, revealed the potential of analyzing extracellular metabolomic data in the context of the metabolic model using constraint-based modeling. With the MetaboTools, we make our methods available to the broader scientific community. The MetaboTools consist of a protocol, a toolbox, and tutorials of two use cases. The protocol describes, in a step-wise manner, the workflow of data integration, and computational analysis. The MetaboTools comprise the Matlab code required to complete the workflow described in the protocol. Tutorials explain the computational steps for integration of two different data sets and demonstrate a comprehensive set of methods for the computational analysis of metabolic models and stratification thereof into different phenotypes. The presented workflow supports integrative analysis of multiple omics data sets. Importantly, all analysis tools can be applied to metabolic models without performing the entire workflow. Taken together, the MetaboTools constitute a comprehensive guide to the intra-model analysis of extracellular metabolomic data from microbial, plant, or human cells. In conclusion, this computational modeling resource offers a broad set of computational analysis tools for a wide biomedical and non-biomedical research community.

59 BASIC BIOLOGICAL SCIENCES↗

Temperature and precipitation extremes in century-long gridded observations, reanalyses, and atmospheric model simulations

Knowledge about long-term changes in climate extremes is vital to better understand multidecadal climate variability and long-term changes and to place today's extreme events in a historical context. While global changes in temperature and precipitation extremes since the midtwentieth century are well studied, knowledge about century-scale changes is limited. This paper analyses a range of largely independent observations-based data sets covering 1901–2010 for long-term changes and interannual variability in daily scale temperature and precipitation extremes. We compare across data sets for consistency to ascertain our confidence in century-scale changes in extremes. We find consistent warming trends in temperature extremes globally and in most land areas over the past century. For precipitation extremes we find global tendencies toward more intense rainfall throughout much of the twentieth century; however, local changes are spatially more variable. While global time series of the different data sets agree well after about 1950, they often show different changes during the first half of the twentieth century. In regions with good observational coverage, gridded observations and reanalyses agree well throughout the entire past century. Finally, simulations with an atmospheric model suggest that ocean temperatures and sea ice may explain up to about 50% of interannual variability in the global average of temperature extremes, and about 15% in the global average of moderate precipitation extremes, but local correlations are mostly significant only in low latitudes.

54 ENVIRONMENTAL SCIENCES↗

Multivariate physiological recordings in an experimental hemorrhage model

In this paper we describe a data set of multivariate physiological measurements recorded from conscious sheep ( N = 8; 37.4 ± 1.1 kg) during hemorrhage. Hemorrhage was experimentally induced in each animal by withdrawing blood from a femoral artery at two different rates (fast: 1.25 mL/kg/min; and slow: 0.25 mL/kg/min). Data, including physiological waveforms and continuous/intermittent measurements, were transformed to digital file formats (European Data Format [EDF] for waveforms and Comma-Separated Values [CSV] for continuous and intermittent measurements) as a comprehensive data set and stored and publicly shared here (Appendix A). The data set comprises experimental information (e.g., hemorrhage rate, animal weight, event times), physiological waveforms (arterial and central venous blood pressure, electrocardiogram), time-series records of non-invasive physiological measurements (SpO 2 , tissue oximetry), intermittent arterial and venous blood gas analyses (e.g., hemoglobin, lactate, SaO 2 , SvO 2 ) and intermittent thermodilution cardiac output measurements. A detailed explanation of the hemodynamic and pulmonary changes during hemorrhage is available in a previous publication (Scully et al., 2016)

Yaghouby, Farid↗

Tensor Minkowski Functionals: first application to the CMB

Tensor Minkowski Functionals (TMFs) are tensor generalizations of the usual Minkowski Functionals which are scalar quantities. We introduce them here for use in cosmological analysis, in particular to analyze the Cosmic Microwave Background (CMB) radiation. They encapsulate information about the shapes of structures and the orientation of distributions of structures. We focus on one of the TMFs, namely W {sub 2}{sup 1,1}, which is the (1,1) rank tensor generalization of the genus. The ratio of the eigenvalues of the average of W {sub 2}{sup 1,1} over all structures, α, encodes the net orientation of the structures; and the average of the ratios of the eigenvalues of W {sub 2}{sup 1,1} for each structure, β, encodes the net intrinsic anisotropy of the structures. We have developed a code that computes W {sub 2}{sup 1,1}, and from it α and β, for a set of structures on the 2-dimensional Euclidean plane. We use it to compute α and β as functions of chosen threshold levels for simulated Gaussian and isotropic CMB temperature and E mode fields. We obtain the value of α to be one for both temperature and E mode, which means that we recover the statistical isotropy of density fluctuations that we input in the simulations. We find that the standard ΛCDM model predicts a charateristic shape of β for temperature and E mode as a function of the threshold, and the average over thresholds is β∼ 0.62 for temperature and β∼ 0.63 for E mode. Accurate measurements of α and β can be used to test the standard model of cosmology and to search for deviations from it. For this purpose we compute α and β for temperature and E mode data of various data sets from PLANCK mission. We compare the values measured from observed data with those obtained from simulations to which instrument beam and noise characteristics of the 44GHz frequency channel have been added (which are provided as part of the PLANCK data release). We find very good agreement of β and α between all PLANCK temperature data sets with ΛCDM expectations. E mode data show good agreement for β but α for all data sets deviate from ΛCDM predictions higher than 3−σ. It is most likely that the deviations are probing the anisotropy of the noise field and beam characteristics of the detector rather than the true E mode signal since for 44GHz the signal-to-noise ratio is well below one. This will be further investigated after the full PLANCK data becomes publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SPRUCE Whole Ecosystems Warming (WEW) Environmental Data Beginning August 2015

This data set provides a record of the half-hourly averages of automated environmental data collected for 12 SPRUCE plots (4, 6, 7, 8, 10, 11, 13, 16, 17, 19, 20, and 21) beginning during deep peat heating (DPH) and continuing throughout the whole ecosystem warming (WEW) manipulations for the SPRUCE Project (Hanson et al. 2017). In August 2015, WEW was initiated at 5 warming levels (+0, +2.25 +4.5, +6.75 and +9 °C) with 2 plots per warming level. DPH measurements were underway before the initiation of WEW heating treatments and both are expected to operate through 2025. This current version includes data from 2014 through 2024. This data set includes 15 data files provided in comma separated (*.csv) format: 12 individual data files for each of the monitored SPRUCE plots, two data files corresponding to figures from Hanson et al. (2017) that support key analyses of the performance of the WEW systems, and one data file containing snow and ice depths.

54 ENVIRONMENTAL SCIENCES↗

Full-sky Gravitational Lensing Simulation for Large-area Galaxy Surveys and Cosmic Microwave Background Experiments

We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological N-body simulations. The data sets include full-sky convergence and shear maps from redshifts z = 0.05 to 5.3 at intervals of $150\,{h}^{-1}\mathrm{Mpc}$ comoving radial distance (corresponding to a redshift interval of ${\rm{\Delta }}z\simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to z = 5.3. The dark matter halos are identified from the same N-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and luminous red galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy–galaxy and cluster–galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy–galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within 5% up to ${\ell }=3000$ (or at an angular scale $\theta \gt 0.5$ arcmin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey, but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing/.

79 ASTRONOMY AND ASTROPHYSICS↗

Measures of soil, sediment, and groundwater background used by Los Alamos National Laboratory, 2012

This document summarizes approaches used by Los Alamos National Laboratory (LANL) to estimate background concentrations of constituents in soils, sediment, and groundwater. For the sake of presentation and interpretation, LANL frequently reports a “background value” that statistically represents an upper bound for a data set of background concentrations. However, actual background concentrations vary by soil and bedrock type, disturbance regimes, precipitation levels, vegetative cover, and a wide variety of other factors. The LANL reported background values are estimates for comparison, not absolute limits of possible background concentrations. In evaluating whether a sample result represents a background concentration, LANL staff evaluate the results statistically (if the sample size is large enough) or graphically against the complete distribution of the relevant background data set. For analytes with no established background data set, LANL site results are sometimes compared to upgradient or geographically distant results on a study-specific basis.

Environmental Protection↗

Variations in the Contouring of Organs at Risk: Test Case From a Patient With Oropharyngeal Cancer

Purpose: Anatomy contouring is critical in radiation therapy. Inaccuracy and variation in defining critical volumes will affect everything downstream: treatment planning, dose-volume histogram analysis, and contour-based visual guidance used in image-guided radiation therapy. This study quantified: (1) variation in the contouring of organs at risk (OAR) in a clinical test case and (2) corresponding effects on dosimetric metrics of highly conformal plans. Methods and Materials: A common CT data set with predefined targets from a patient with oropharyngeal cancer was provided to a population of clinics, which were asked to (1) contour OARs and (2) design an intensity-modulated radiation therapy plan. Thirty-two acceptable plans were submitted as DICOM RT data sets, each generated by a different clinical team. Using those data sets, we quantified: (1) the OAR contouring variation and (2) the impact this variation has on dosimetric metrics. New technologies were employed, including a software tool to quantify three-dimensional structure comparisons. Results: There was significant interclinician variation in OAR contouring. The degree of variation is organ-dependent. We found substantial dose differences resulting strictly from contouring variation (differences ranging from -289% to 56% for mean OAR dose; -22% to 35% for maximum dose). However, there appears to be a threshold in the OAR comparison metric beyond which the dose differences stabilize. Conclusions: The effects of interclinician variation in contouring organs-at-risk in the head and neck can be large and are organ-specific. Physicians need to be aware of the effect that variation in OAR contouring can play on the final treatment plan and not restrict their focus only to the target volumes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Prediction of Distributed River Sediment Respiration Rates Using Community-Generated Data and Machine Learning

River sediment microbial respiration is a key indicator of ecosystem functioning and the biogeochemical fluxes across this critical zone link surface and subsurface waters. As such, there is tremendous interest in measuring and mapping these respiration rates. Respiration observations are expensive and labor intensive; there is limited data available to the community. An open science, collaborative initiative is collecting samples for respiration rate analysis and multi-scale metadata; this evolving data set is being used for making machine learning (ML) predictions at unsampled sites to help inform continued community engagement. However, it is a challenge to find an optimum configuration for ML models to work with this feature-rich (i.e., 100+ possible input variables) data set. Here, we present results from a two-tiered approach to managing the analysis of this complex data set: (a) a stacked ensemble of models that automatically optimizes hyperparameters and manages the training of many models and (b) feature permutation importance to detect the most important features in the models. The major elements of this workflow are modular, portable, open, and cloud-based thus making this implementation a potential template for other applications. The models developed here predict that sediment organic matter chemistry is one of the most important features for predicting sediment respiration rate. Other larger-scale, important features fall into the categories of climatic, ecological, geological, and fluvial settings. Leveraging these larger-scale features to generate data-driven estimates of river sediment respiration rates reveals spatially consistent but heterogeneous patterns across the river network of the Columbia River Basin.

54 ENVIRONMENTAL SCIENCES↗

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗

A nonvoxel-based dose convolution/superposition algorithm optimized for scalable GPU architectures

Purpose: Real-time adaptive planning and treatment has been infeasible due in part to its high computational complexity. There have been many recent efforts to utilize graphics processing units (GPUs) to accelerate the computational performance and dose accuracy in radiation therapy. Data structure and memory access patterns are the key GPU factors that determine the computational performance and accuracy. In this paper, the authors present a nonvoxel-based (NVB) approach to maximize computational and memory access efficiency and throughput on the GPU. Methods: The proposed algorithm employs a ray-tracing mechanism to restructure the 3D data sets computed from the CT anatomy into a nonvoxel-based framework. In a process that takes only a few milliseconds of computing time, the algorithm restructured the data sets by ray-tracing through precalculated CT volumes to realign the coordinate system along the convolution direction, as defined by zenithal and azimuthal angles. During the ray-tracing step, the data were resampled according to radial sampling and parallel ray-spacing parameters making the algorithm independent of the original CT resolution. The nonvoxel-based algorithm presented in this paper also demonstrated a trade-off in computational performance and dose accuracy for different coordinate system configurations. In order to find the best balance between the computed speedup and the accuracy, the authors employed an exhaustive parameter search on all sampling parameters that defined the coordinate system configuration: zenithal, azimuthal, and radial sampling of the convolution algorithm, as well as the parallel ray spacing during ray tracing. The angular sampling parameters were varied between 4 and 48 discrete angles, while both radial sampling and parallel ray spacing were varied from 0.5 to 10 mm. The gamma distribution analysis method (γ) was used to compare the dose distributions using 2% and 2 mm dose difference and distance-to-agreement criteria, respectively. Accuracy was investigated using three distinct phantoms with varied geometries and heterogeneities and on a series of 14 segmented lung CT data sets. Performance gains were calculated using three 256 mm cube homogenous water phantoms, with isotropic voxel dimensions of 1, 2, and 4 mm. Results: The nonvoxel-based GPU algorithm was independent of the data size and provided significant computational gains over the CPU algorithm for large CT data sizes. The parameter search analysis also showed that the ray combination of 8 zenithal and 8 azimuthal angles along with 1 mm radial sampling and 2 mm parallel ray spacing maintained dose accuracy with greater than 99% of voxels passing the γ test. Combining the acceleration obtained from GPU parallelization with the sampling optimization, the authors achieved a total performance improvement factor of >175 000 when compared to our voxel-based ground truth CPU benchmark and a factor of 20 compared with a voxel-based GPU dose convolution method. Conclusions: The nonvoxel-based convolution method yielded substantial performance improvements over a generic GPU implementation, while maintaining accuracy as compared to a CPU computed ground truth dose distribution. Such an algorithm can be a key contribution toward developing tools for adaptive radiation therapy systems.

60 APPLIED LIFE SCIENCES↗

Global Carbon Budget 2017

Here an accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the global carbon budget – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink ( S OCEAN ) and terrestrial CO 2 sink ( S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance ( B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1 σ . For the last decade available (2007–2016), E FF was 9.4 ± 0.5 GtC yr –1 , E LUC 1.3 ± 0.7 GtC yr –1 , G ATM 4.7 ± 0.1 GtC yr –1 , S OCEAN 2.4 ± 0.5 GtC yr –1 , and S LAND 3.0 ± 0.8 GtC yr –1 , with a budget imbalance B IM of 0.6 GtC yr –1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ±\ 0.5 GtC yr –1 . Also for 2016, E LUC was 1.3 ± .7 GtC yr –1 , G ATM was 6.1 ± 0.2 GtC yr –1 , S OCEAN was 2.6 ± 0.5 GtC yr –1 , and S LAND was 2.7 ± 1.0 GtC yr –1 , with a small B IM of –0.3 GtC. G ATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small S LAND consistent with El Niño conditions. The global atmospheric CO 2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in E FF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set.

54 ENVIRONMENTAL SCIENCES↗

Prostate implant reconstruction from C-arm images with motion-compensated tomosynthesis

Purpose: Accurate localization of prostate implants from several C-arm images is necessary for ultrasound-fluoroscopy fusion and intraoperative dosimetry. The authors propose a computational motion compensation method for tomosynthesis-based reconstruction that enables 3D localization of prostate implants from C-arm images despite C-arm oscillation and sagging. Methods: Five C-arm images are captured by rotating the C-arm around its primary axis, while measuring its rotation angle using a protractor or the C-arm joint encoder. The C-arm images are processed to obtain binary seed-only images from which a volume of interest is reconstructed. The motion compensation algorithm, iteratively, compensates for 2D translational motion of the C-arm by maximizing the number of voxels that project on a seed projection in all of the images. This obviates the need for C-arm full pose tracking traditionally implemented using radio-opaque fiducials or external trackers. The proposed reconstruction method is tested in simulations, in a phantom study and on ten patient data sets. Results: In a phantom implanted with 136 dummy seeds, the seed detection rate was 100% with a localization error of 0.86 {+-} 0.44 mm (Mean {+-} STD) compared to CT. For patient data sets, a detection rate of 99.5% was achieved in approximately 1 min per patient. The reconstruction results for patient data sets were compared against an available matching-based reconstruction method and showed relative localization difference of 0.5 {+-} 0.4 mm. Conclusions: The motion compensation method can successfully compensate for large C-arm motion without using radio-opaque fiducial or external trackers. Considering the efficacy of the algorithm, its successful reconstruction rate and low computational burden, the algorithm is feasible for clinical use.

62 RADIOLOGY AND NUCLEAR MEDICINE↗