Search NASA⌕ Search

SEARCH · Search NASA

Results for “Uncertainty visualization”

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 163 records · Page 9

Assessment of Local Observation of Atomic Ordering in Alloys via the Radial Distribution Function: A Computational and Experimental Approach

As a powerful analytical technique, atom probe tomography (APT) has the capacity to acquire the spatial distribution of millions of atoms from a complex sample. However, extracting information at the Ångstrom-scale on atomic ordering remains a challenge due to the limits of the APT experiment and data analysis algorithms. The development of new computational tools enable visualization of the data and aid understanding of the physical phenomena such as disorder of complex crystalline structures. Here, we report progress towards this goal using two steps. We describe a computational approach to evaluate atomic ordering in the crystal structure by generating radial distribution functions (RDF). Atomic ordering is rendered as the Fractional Cumulative Radial Distribution Function (FCRDF) which allows for greater visibility of local compositions at short range in the structure. Further, we accommodate in the analysis additional parameters such as uncertainty in the atomic coordinates and the atomic abundance to ascertain short-range ordering in APT data sets. We applied the FCRDF analysis to synthetic and experimental APT data sets for Ni 3 Al. The ability to observe a signal of atomic ordering consistent with the known L1 2 crystal structure is heavily dependent on spatial uncertainty, irrespective of abundance. Detection of atomic ordering is subject to an upper limit of spatial uncertainty of atoms described with Gaussian distributions with a standard deviation of 1.3 Å. The FCRDF analysis was also applied to the APT data set for a six-component alloy, Al 1.3 CoCrCuFeNi. In this case, we are currently able to visualize elemental segregation at the nanoscale, though unambiguous identification of atomic ordering at the Ångstrom (nearest-neighbor) scale remains a goal.

36 MATERIALS SCIENCE↗

Assessing the Sensitivity of Pourbaix Diagrams to Computational Protocols: Electrochemical Stability of Ni Oxides as a Case Study

Pourbaix diagrams stand as a useful tool in assessing and visualizing materials’ electrochemical stability and are widely used for electrocatalyst design. However, their reliability hinges on the accuracy of the chemical potentials of involved phases, which may bear uncertainties and can be significantly impacted by decision-making steps in the computational protocol. Here, this study introduces a robust sensitivity analysis framework, exemplified through a detailed examination of the computational Pourbaix diagram of Ni, the oxides of which are used as high-activity and cost-friendly catalysts for many electrochemical reactions. Quantities of interest derived from the Pourbaix diagram include the appearance and stability domain of the catalytically active Ni oxide phases along with the onset electrochemical potentials of phase transitions. These metrics can guide the design of operational conditions for Ni oxide electrocatalysts. We find that the employed DFT exchange-correlation functional has the most significant influence on the computed Pourbaix diagram. Uncertainties on crystal structures, along with their related ab initio energetics, are also found to affect the size of the phase stability domain. Higher-order coupling among input parameters is found to play a crucial role in influencing the appearance and distribution of Ni phases in the diagram. Our findings suggest a need to consider variations and uncertainties associated with the computational procedures on predicted Pourbaix diagrams for materials design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗

SU-F-T-566: Absolute Film Dosimetry for Stereotactic Radiosurgery and Stereotactic Body Radiotherapy Quality Assurance Using Gafchromic EBT3 Films

Purpose: To evaluate the dosimetric uncertainty associated with Gafchromic (EBT3) films and establish an absolute dosimetry protocol for Stereotactic Radiosurgery (SRS) and Stereotactic Body Radiotherapy (SBRT). Methods: EBT3 films were irradiated at each of seven different dose levels between 1 and 15 Gy with open fields, and standard deviations of dose maps were calculated at each color channel for evaluation. A scanner non-uniform response correction map was built by registering and comparing film doses to the reference diode array-based dose map delivered with the same doses. To determine the temporal dependence of EBT3 films, the average correction factors of different dose levels as a function of time were evaluated up to four days after irradiation. An integrated film dosimetry protocol was developed for dose calibration, calibration curve fitting, dose mapping, and profile/gamma analysis. Patient specific quality assurance (PSQA) was performed for 93 SRS/SBRT treatment plans. Results: The scanner response varied within 1% for the field sizes less than 5 × 5 cm{sup 2}, and up to 5% for the field sizes of 10 × 10 cm{sup 2}. The scanner correction method was able to remove visually evident, irregular detector responses found for larger field sizes. The dose response of the film changed rapidly (∼10%) in the first two hours and plateaued afterwards, ∼3% change between 2 and 24 hours. The mean uncertainties (mean of the standard deviations) were <0.5% over the dose range 1∼15Gy for all color channels for the OD response curves. The percentage of points passing the 3%/1mm gamma criteria based on absolute dose analysis, averaged over all tests, was 95.0 ± 4.2. Conclusion: We have developed an absolute film dose dosimetry protocol using EBT3 films. The overall uncertainty has been established to be approximately 1% for SRS and SBRT PSQA. The work was supported by a Research Scholar Grant, RSG-15-137-01-CCE from the American Cancer Society.

60 APPLIED LIFE SCIENCES↗

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗

Observations of Binary Stars with the Differential Speckle Survey Instrument. IX. Observations of Known and Suspected Binaries, and a Partial Survey of Be Stars

We report 370 measures of 170 components of binary and multiple-star systems, obtained from speckle imaging observations made with the Differential Speckle Survey Instrument at Lowell Observatory’s Discovery Channel Telescope in 2015 through 2017. Of the systems studied, 147 are binary stars, 10 are seen as triple systems, and 1 quadruple system is measured. Seventy-six high-quality nondetections and 15 newly resolved components are presented in our observations. The uncertainty in relative astrometry appears to be similar to our previous work at Lowell, namely, linear measurement uncertainties of approximately 2 mas, and the relative photometry appears to be uncertain at the 0.1–0.15 mag level. Using these measures and those in the literature, we calculate six new visual orbits, including one for the Be star 66 Oph and two combined spectroscopic–visual orbits. The latter two orbits, which are for HD 22451 (YSC 127) and HD 185501 (YSC 135), yield individual masses of the components at the level of 2% or better, and independent distance measures that in one case agrees with the value found in the Gaia DR2 and in the other disagrees at the 2σ level. We find that HD 22451 consists of an F6V+F7V pair with orbital period of 2401.1 ± 3.2 days and masses of 1.342 ± 0.029 and 1.236±0.026 M{sub ⊙}. For HD 185501, both stars are G5 dwarfs that orbit one another with a period of 433.94 ± 0.15 days, and the masses are 0.898 ± 0.012 and 0.876±0.012 M{sub ⊙}. We discuss the details of both the new discoveries and the orbit objects.

79 ASTRONOMY AND ASTROPHYSICS↗

Cesium Batch Contact Equilibrium Testing of Crystalline Silicotitanate (CST) Sieve Cuts in SRS Average Simulant and Examination of CST Samples Aged in Caustic and Simulant in Support of Tank Closure Cesium Removal 1A (TCCR-1A)

Batch contact testing to determine cesium equilibrium loading on Crystalline Silicotitanate ion exchange media in Savannah River Site Average Waste Simulant at 25 °C indicated that smaller diameter particles isolated by sieving pretreated CST media may load slightly higher (<10%) amounts of cesium, though the differences are within analytical uncertainty. In addition, ion exchange media sub-samples stored in 2-4 M NaOH and caustic simulant solutions for ~2.5 years were examined by optical microscopy and the 3 M NaOH sample was also analyzed to determine whether changes in the particle size distribution occurred during storage. No visual indications of particle attrition or agglomeration were observed for any sample. Particle size analysis indicated that a slight decrease occurred in the average particle diameter following contact with 3 M NaOH (541 µm average diameter versus 566 µm for the pretreated CST prior to contact). A small increase (from 0 to <0.5 wt. %) in the number of particles ranging from 271 and 322 µm was also observed for the CST sample contacted with 3 M NaOH relative to a sample of the original pretreated material. However, this small change could be due to sub-sampling differences or analytical uncertainty. It does not appear that small particles are formed to a significant degree during CST caustic contact or that small particles which do form (presumably from attrition of larger particles during pretreatment) load significantly more cesium than the bulk material. Minimal other negative consequences were observed associated with CST extended caustic or simulant contact, except for the tendency for more concentrated salt solutions to form some salt crystals which deposit on the media over time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-E-J-85: The Effect of Different Imaging Modalities On the Delineation of the True Spinal Cord for Spinal Stereotactic Body Radiation Therapy

Purpose: SBRT allows the delivery of high dose radiation treatments to localized tumors while minimizing dose to surrounding tissues. Due to the large doses delivered, accurate contouring of organs at risk is essential. In this study, differences between the true spinal cord as seen using MRI and CT myelogram (CTM) have been assessed in patients with spinal metastases treated using SBRT. Methods: Ten patients were identified that have both a CTM and a MRI. Using rigid registration tools, the MRI was fused to the CTM. The thecal sac and true cord were contoured using each imaging modality. Images were exported and analyzed for similarity by computing the Dice similarity coefficient and the modified Hausdorff distance (greatest distance from a point in one set to the closest point in the other set). Results: The Dice coefficient was calculated for the thecal sac (0.81 ±0.06) and true cord (0.63 ±0.13). These two measures are correlated; however, some points show a low true cord overlap despite a high overlap for the thecal sac. The Hausdorff distance for structure comparisons was also calculated. For thecal sac structures, the average value, 1.6mm (±1.1), indicates good overlap. For true cord comparison, the average value, 0.3mm (±0.16), indicates very good overlap. The minimum Hausdorff distance between the true cord and thecal sac was on average 1.6mm (±0.9) Conclusion: The true cord position as seen in MRI and CTM is fairly constant, although care should be taken as large differences can be seen in individual patients. Avoidning the true cord in spine SBRT is critical, so the ability to visualize the true cord before performing SBRT to the vertebrae is essential. Here, CT myelogram appears an excellent, robust option, that can be obtained the day of treatment planning and is unaffected by uncertainties in image fusion.

60 APPLIED LIFE SCIENCES↗

Microreactor Optimization Using Simulation And Economics (mouse)

Microreactor Optimization Using Simulation and Economics (MOUSE) is a tool that integrates both nuclear microreactor design and reactor economics to provide comprehensive evaluations and optimizations. This tool enables stakeholders to explore the interplay between technical and economic variables, guiding them towards effective and competitive microreactor solutions. For the reactor core simulations, MOUSE leverages the OpenMC Monte Carlo Particle Transport Code to perform detailed core simulations for various microreactor designs. The included OpenMC models are 2D core designs of a Liquid Metal Thermal Microreactor (LMTR), a Gas-Cooled TRISO-Fueled Microreactor (GCMR), and a Heat Pipe Microreactor. Beyond core design, MOUSE includes simplified calculations for: - Calculating the masses of heat exchangers within the system. - Mechanical power of pumps. - Estimating the area occupied by various buildings within the nuclear plant. For the economic analysis, MOUSE provides detailed bottom-up cost estimates, encompassing a wide range of costs including preconstruction costs, direct costs, indirect costs, training costs, financial costs, operation & maintenance (O&M) costs, and fuel costs. These cost estimations are developed using data from the MARVEL project and additional literature sources, enabling the calculation of total capital costs and levelized cost of energy for both first-of-a-kind and nth-of-a-kind microreactors. MOUSE also enables analysis of the cost drivers and competitiveness in the electricity market. MOUSE allows users to modify a wide array of technical and economic parameters to evaluate different scenarios and their impacts. Examples of these parameters include: Fuels, coolants, or reflector materials Enrichment levels Control drum materials and geometry Fuel pin geometry and materials Moderator pin geometry and materials Reactor core and reflector dimensions Packing factor for the TRISO particles Nuclear reactor power and reactor burnup Number of sensors Shielding thickness Reactor vessel and guard vessel dimensions Operational staff requirements Number of emergency shutdowns Levelization period Interest rate Construction duration Since MOUSE is powered by the WATTS toolkit, it supports optimization studies, parametric analyses, and uncertainty calculations/propagation. The optimization techniques enable users to identify optimal design and economic configurations. The parametric analysis tools allow users to explore the sensitivity of various parameters, while uncertainty propagation helps quantify the impact of uncertainties on overall performance and cost. User Interface and Workflow: Currently, MOUSE is a command-line-based tool. Users can input various reactor design or economic parameters, modify the designs, run simulations, and visualize results through comprehensive data visualization and reporting capabilities. The typical workflow involves setting up the reactor model, defining economic parameters, running simulations, and analyzing the results to make informed decisions. By combining advanced design calculations with detailed economic modeling, MOUSE provides a robust framework for optimizing nuclear microreactor technologies, enhancing their competitiveness, and guiding stakeholders towards innovative and cost-effective solutions.

Hanna, Botros [Idaho National Laboratory (INL), Id↗

Orbits and Masses of Binaries from Speckle Interferometry at SOAR

We present results from Speckle inteferometric observations of 15 visual binaries and one double-line spectroscopic binary, carried out with the HRCam Speckle camera of the SOAR 4.1 m telescope. These systems were observed as a part of an on-going survey to characterize the binary population in the solar vicinity, out to a distance of 250 pc. We obtained orbital elements and mass sums for our sample of visual binaries. The orbits were computed using a Markov Chain Monte Carlo algorithm that delivers maximum likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate their uncertainty. Their periods cover a range from 5 yr to more than 500 yr; and their spectral types go from early A to mid M, implying total system masses from slightly more than 4M{sub ⊙} down to 0.2M {sub ⊙}. They are located at distances between approximately 12 and 200 pc, mostly at low Galactic latitude. For the double-line spectroscopic binary YSC8, we present the first combined astrometric/radial-velocity orbit resulting from a self-consistent fit, leading to individual component masses of 0.897 ± 0.027 M {sub ⊙} and 0.857 ± 0.026 M {sub ⊙}; and an orbital parallax of 26.61 ± 0.29 mas, which compares very well with the Gaia DR2 trigonometric parallax (26.55 ± 0.27 mas). In combination with published photometry and trigonometric parallaxes, we place our objects on an H-R diagram and discuss their evolutionary status. We also present a thorough analysis of the precision and consistency of the photometry available for them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗

Model-Data for Joint Estimation of Biogeochemical Model Parameters from Multiple Experiments: A Bayesian Approach Applied to Mercury Methylation

This modeling archive supports the manuscript submitted for publication in the Environmental Modeling and Software. This study is supported by ORNL-SFA and IDEAS-Watershed. This study aims to improve calibration of complex biogeochemical models using datasets from multiple experiments targeting specific subprocesses. The proposed Bayesian joint-fitting scheme calibrates the entire biogeochemical model in one go using all the available datasets and estimate parameter uncertainties using Markov Chain Monte Carlo (MCMC). This allows for complete propagation of uncertainties and utilization of the information shared between different datasets. Mapping joint distribution of parameters guides model improvement by identifying null spaces in the parameter space. This archive contains files used to perform MCMC, post-process outputs and visualize results.

East Fork Poplar Creek↗

Visual Exploration of Semantic Relationships in Neural Word Embeddings

Constructing distributed representations for words through neural language models and using the resulting vector spaces for analysis has become a crucial component of natural language processing (NLP). But, despite their widespread application, little is known about the structure and properties of these spaces. To gain insights into the relationship between words, the NLP community has begun to adapt high-dimensional visualization techniques. Particularly, researchers commonly use t-distributed stochastic neighbor embeddings (t-SNE) and principal component analysis (PCA) to create two-dimensional embeddings for assessing the overall structure and exploring linear relationships (e.g., word analogies), respectively. Unfortunately, these techniques often produce mediocre or even misleading results and cannot address domain-specific visualization challenges that are crucial for understanding semantic relationships in word embeddings. We introduce new embedding techniques for visualizing semantic and syntactic analogies, and the corresponding tests to determine whether the resulting views capture salient structures. Additionally, we introduce two novel views for a comprehensive study of analogy relationships. Finally, we augment t-SNE embeddings to convey uncertainty information in order to allow a reliable interpretation. Combined, the different views address a number of domain-specific tasks difficult to solve with existing tools.

97 MATHEMATICS AND COMPUTING↗

Semi-inclusive Neutral Current Neutral Pion Production Selection at the NOvA (NuMI Off-axis Electron Neutrino Appearance) Near Detector Using Prong Level Convolutional Neural Networks

The NONA neutrino experiment based in Fermilab is designed to measure ν µ →ν e neutrino oscillations. This experiment will give us insight into the properties of massive neutrinos. Neutral current (NC) ν µ,e neutral pion production events can mimic the ν µ →ν e oscillation signal and therefore are an important background for NONA to understand. Neutral pions decay into two photons which can fake a single electron shower (νe appearance signal) in two ways: either the 2 photons can merge together or one of them may escape detection. In order to constrain this background, NONA utilizes the Near Detector to measure neutral current π 0 neutrino interactions. In this analysis, neutrino-Nucleus (ν µ →N) NC π 0 interactions with total pion energy greater then 0.3 GeV are studied by selecting two prong events with two final state photons as determined by prong based Convolutional Visual Networks (CVN). The analysis is performed on 3.54x10 21 Protons On Target (POT) of NONA Near Detector simulated data and compared to 8.09x10 20 POT of data. Optimization of the selection based on fractional cross-section uncertainty and an initial energy resolution study of the final sample are presented. The final 2 prong selection using prong based CVN gave a purity of 74%, selection efficiency of ~ 1.8%, and an expected (NC) ν µ,e neutral pion cross section of ~ 15.2%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Bayesian Machine Learning Model for Estimating Building Occupancy from Open Source Data

Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian based informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000ft2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian model and informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000 ft 2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗