Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty modeling”

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 397 records · Page 22

IceCube's convex all-sky neutrino spectrum consistent with the magnetically powered corona scenario for active galactic nuclei

High-energy multimessenger background analyses over the past decade have provided evidence for a population of hidden neutrino sources that are opaque to GeV-TeV gamma rays, a picture bolstered by recent observations of the nearby active galaxy NGC 1068. The coronal regions in the hearts of active galactic nuclei (AGNs) have been proposed as the most promising sites for such hidden nonthermal particle production, and NGC 1068 is expected to be the most neutrino-active galaxy for IceCube. We demonstrate that the latest all-sky neutrino spectrum, exhibiting a spectral bend around 3-30 TeV, is consistent with predictions of the magnetically powered corona scenario, and the models for the all-sky neutrino flux can simultaneously explain the multimessenger data from NGC 1068 within observational and modeling uncertainties. We further show, in a largely model-independent way, that the contribution from NGC 1068-like sources does not overshoot the observed medium-energy neutrino flux. Finally, we highlight the key role of the Eddington ratio, which can drive substantial variations in the predicted neutrino fluxes of nearby AGNs, and we encourage systematic multimessenger searches for the neutrino-brightest AGNs.

Murase, Kohta [Penn State U.; Penn State U., Unive↗

Quantifying subsurface parameter and transport uncertainty using surrogate modelling and environmental tracers

Here, we combine physics-based groundwater reactive transport modelling with machine-learning techniques to quantify hydrogeological model and solute transport predictive uncertainties. We train an artificial neural network (ANN) on a dataset of groundwater hydraulic heads and 3 H concentrations generated using a high-fidelity groundwater reactive transport model. Using the trained ANN as a surrogate model to reproduce the input–output response of the high-fidelity reactive transport model, we quantify the posterior distributions of hydrogeological parameters and hydraulic forcing conditions using Markov chain Monte Carlo calibration against field observations of groundwater hydraulic heads and 3 H concentrations. We demonstrate the methodology with a model application that predicts Chlorofluorocarbon-12 (CFC-12) solute transport at a contaminated field site in Wyoming, United States. Our results show that including 3 H observations in the calibration dataset reduced the uncertainty in the estimated permeability field and infiltration rates, compared to calibration against hydraulic heads alone. However, predictive uncertainty quantification shows that CFC-12 transport predictions conditioned to the parameter posterior distributions cannot reproduce the field measurements. We found that calibrating the model to hydraulic head and 3 H observations results in groundwater mean ages that are too large to explain the observed CFC-12 concentrations. The coupling of the physics-based reactive transport model with the machine-learning surrogate model allows us to efficiently quantify model parameter and predictive uncertainties, which is typically computationally intractable using reactive transport models alone.

58 GEOSCIENCES↗

Adjusting neutrino interaction models and evaluating uncertainties using NOvA near detector data

The two-detector design of the NOvA neutrino oscillation experiment, in which two functionally identical detectors are exposed to an intense neutrino beam, aids in canceling leading order effects of cross-section uncertainties. However, limited knowledge of neutrino interaction cross sections still gives rise to some of the largest systematic uncertainties in current oscillation measurements. We show contemporary models of neutrino interactions to be discrepant with data from NOvA, consistent with discrepancies seen in other experiments. Adjustments to neutrino interaction models in GENIE are presented, creating an effective model that improves agreement with our data. We also describe systematic uncertainties on these models, including uncertainties on multi-nucleon interactions from a newly developed procedure using NOvA near detector data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uncertainty Quantification using Deep Ensembles for Decision Making in Cyber-Physical-Human Systems

In this paper and its companion, Differential Equation Approximation Using Gradient-Boosted Quantile Regression, Robison et al., we examine an approach to quantifying model uncertainty with the aim of increasing the trustworthiness of computational models in human-machine interactions. In Differential Equation Approximation Using Gradient-Boosted Quantile Regression, we focus on gradient-boosted decision trees, while in this one, we give more details about deep ensembles. Uncertainty quantification is crucial for building trustworthy autonomous decision-making agents in human-machine teams. There are two types of uncertainties: aleatoric and epistemic. The former is related to the inherent stochasticity (noise) of the process, whereas the latter is associated with the lack of knowledge or representation capability of models, such as neural networks. By lack of knowledge, we mean the model’s inability to accurately predict outputs for all possible inputs. The aleatory uncertainty can be estimated fairly easily with, for example, filters, whereas epistemic uncertainty is challenging to compute. This paper uses deep ensembles to quantify both aleatory and epistemic uncertainty. It can act as an uncertainty-aware surrogate transition model for decision-making frameworks. "Uncertainty-aware" means that the surrogate transition model should make predictions along with confidence in those predictions. In the context of decision-making, the transition models are ordinary differential equations (ODEs). Since ODEs can be simulated to make one-step or multi-step predictions, a good surrogate model for them should perform reasonably well in both modes. In a multi-step approach, the trajectory sampling method TS∞ was used to propagate uncertainty over multiple steps. The cartpole dynamical system was selected to demonstrate the ability of deep ensembles as good surrogate transition models for decision-making frameworks. The deep ensembles modeled the dynamics of cartpole ODEs and made uncertainty-aware predictions in single-step and multi-step transition modes.

CPH systems↗

First Measurement of Inclusive Muon Neutrino Charged Current Triple Differential Cross Section on Argon

The field of accelerator neutrino experiments is entering an era of precision oscillation measurements where the remaining unknown neutrino measurements will be determined. The upcoming DUNE and Hyper-K experiments aim to determine the neutrino mass hierarchy and degree of Charge-Parity (CP) violation in the neutrino sector, providing potential insight on the matter-antimatter imbalance observed in the universe. However, these experiments require highly accurate measurements, and neutrino cross section modeling uncertainties may limit their capabilities. Cross section measurements at current- generation experiments can aid the development of neutrino interaction models to reduce these uncertainties. This is especially true for measurements of neutrino energy, as it drives neutrino oscillations and is of key importance to oscillation experiments. The MicroBooNE experiment uses a Liquid Argon Time Projection Chamber (LArTPC) to produce neutrino-argon cross sections as one of its physics goals. The MicroBooNE detector’s fully active volume, precision reconstruction, and calorimetry information are leveraged in the Wire-Cell analysis to produce a muon neutrino selection that is 92% pure while maintaining 68% efficiency. A reconstruction chain featuring a fully 3D charge reconstruction and a graph-based particle trajectory fit are used to produce accurate measurements of lepton kinematics as well as visible hadronic energy produced in a neutrino interaction. This thesis presents the first neutrino-argon triple-differential cross section measurement, targeting inclusive charged-current final states. Wiener SVD unfolding is used to produce a measurement over neutrino energy, muon momentum, and muon scattering angle. A series of constrained goodness of fit tests are used to demonstrate the validity of MicroBooNE’s model in describing the distribution of reconstructed kinematics seen in data to ensure the accuracy of unfolding. The validated unfolding to neutrino energy represents a step forward in the field of neutrino cross sections, and demonstrates the capabilities of the LArTPC detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗

Characterizing Performance of Freshwater Wetland Methane Models Across Time Scales at FLUXNET-CH 4 Sites Using Wavelet Analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

58 GEOSCIENCES↗

Projected runoff declines from plant physiological effects on precipitation

Abstract The impact of plants on runoff under high atmospheric CO 2 is a major uncertainty for future water resources. Theory and Earth system models (ESMs) suggest that stricter plant stomatal regulation under high CO 2 will reduce transpiration, potentially boosting runoff. Yet, across a 12-member ensemble of idealized ESM simulations that isolate plant responses to CO 2 , we show that lower transpiration robustly enhances runoff over only 5% of modelled global land area. Precipitation changes are five times more important than transpiration changes in driving runoff responses and are a significant signal of CO 2 physiological forcing over 31–57% of land areas across models. Crucially, ESMs largely disagree on where physiologically forced precipitation changes occur but agree that plant responses in most locations are as likely to reduce runoff as increase it. These results imply that large model uncertainties in precipitation responses, rather than transpiration responses, explain why ESMs disagree on plant physiologically driven runoff changes.

Lesk, Corey S. (ORCID:0000000308684040)↗

Analysis and Optimization of Seismo-Acoustic Monitoring Networks with Bayesian Optimal Experimental Design

The Bayesian optimal experimental design (OED) problem seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty. The goal of OED is to find an experiment that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge about expected data. Therefore, within the context of seismic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize the ability to locate seismic events from arrival time data of detected seismic phases. In order to do utilize Bayesian OED we must develop four elements:1. A likelihood function that describes the uncertainty of detection and travel times; 2. A Bayesian solver that takes a prior and likelihood to identify the posterior; 3. An algorithm to compute EIG; and, 4. An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we can explore many relevant questions to monitoring such as: how and what multiphenomenology data can be used to optimally reduce uncertainty, how to trade off sensor fidelity and earth model uncertainty, and how sensor types, number, and locations influence uncertainty

47 OTHER INSTRUMENTATION↗

Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6

Abstract. The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.

54 ENVIRONMENTAL SCIENCES↗

A layered approach to structural control system design for the JPL phase B testbed

This paper describes our efforts in structural control experiments for a flexible precision structure using both passive and active control. Specifically, a layered structural control approach utilizing passive viscous dampers, colocated and noncolocated active control is used. The passive dampers and active colocated control loops are used to enhance stability and robustness and the noncolocated multivariable controller is designed with respect to a partially controlled structure to further improve the performance. Our approach to the multivariable design problem requires not only a state space model description of the plant, but also an accompanying uncertainty model. To support the multivariable robust control design development, a significant effort in system identification was undertaken. The process of identification and identified results are discussed. Noncolocated controllers are designed using the H(infinity)/mu-synthesis methodology. The order of the controller is typically large, hence, model reduction is performed for practical digital implementation on a real-time control system developed at JPL. Finally, experimental results are presented and some of lessons are discussed.

Chu, Cheng-Chih↗

Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6

The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 °C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 °C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 °C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 °C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.d-2090s.

Scenario Model Intercomparison Project (ScenarioMI↗

Uncertainty quantification in machine learning for engineering design and health prognostics: A tutorial

On top of machine learning (ML) models, uncertainty quantification (UQ) functions as an essential layer of safety assurance that could lead to more principled decision making by enabling sound risk assessment and management. The safety and reliability improvement of ML models empowered by UQ has the potential to significantly facilitate the broad adoption of ML solutions in high-stakes decision settings, such as healthcare, manufacturing, and aviation, to name a few. In this tutorial, we aim to provide a holistic lens on emerging UQ methods for ML models with a particular focus on neural networks and the applications of these UQ methods in tackling engineering design as well as prognostics and health management problems. Towards this goal, we start with a comprehensive classification of uncertainty types, sources, and causes pertaining to UQ of ML models. Next, we provide a tutorial-style description of several state-of-the-art UQ methods: Gaussian process regression, Bayesian neural network, neural network ensemble, and deterministic UQ methods focusing on spectral-normalized neural Gaussian process. Established upon the mathematical formulations, we subsequently examine the soundness of these UQ methods quantitatively and qualitatively (by a toy regression example) to examine their strengths and shortcomings from different dimensions. Then, we review quantitative metrics commonly used to assess the quality of predictive uncertainty in classification and regression problems. Afterward, we discuss the increasingly important role of UQ of ML models in solving challenging problems in engineering design and health prognostics. In conclusion, two case studies with source codes available on GitHub are used to demonstrate these UQ methods and compare their performance in the life prediction of lithium-ion batteries at the early stage (case study 1) and the remaining useful life prediction of turbofan engines (case study 2).

97 MATHEMATICS AND COMPUTING↗

Uncertainty Evaluation of Computational Model Used to Support the Integrated Powerhead Demonstration Project

NASA and the U.S. Air Force are working on a joint project to develop a new hydrogen-fueled, full-flow, staged combustion rocket engine. The initial testing and modeling work for the Integrated Powerhead Demonstrator (IPD) project is being performed by NASA Marshall and Stennis Space Centers. A key factor in the testing of this engine is the ability to predict and measure the transient fluid flow during engine start and shutdown phases of operation. A model built by NASA Marshall in the ROCket Engine Transient Simulation (ROCETS) program is used to predict transient engine fluid flows. The model is initially calibrated to data from previous tests on the Stennis E1 test stand. The model is then used to predict the next run. Data from this run can then be used to recalibrate the model providing a tool to guide the test program in incremental steps to reduce the risk to the prototype engine. In this paper, they define this type of model as a calibrated model. This paper proposes a method to estimate the uncertainty of a model calibrated to a set of experimental test data. The method is similar to that used in the calibration of experiment instrumentation. For the IPD example used in this paper, the model uncertainty is determined for both LOX and LH flow rates using previous data. The successful use of this model is then demonstrated to predict another similar test run within the uncertainty bounds. The paper summarizes the uncertainty methodology when a model is continually recalibrated with new test data. The methodology is general and can be applied to other calibrated models.

Steele, W. G.↗

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral shortwave instruments produce high-resolution and information-dense spectra, thus offering the opportunity to reduce uncertainties in retrieval products. We are in the process of developing a retrieval that uses the very fast Principal Component Radiative Transfer Model in the solar spectral region (PCRTM-Solar) in the forward model calculations. This retrieval will use measured reflectances from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder) instruments. In this manuscript we present progress towards a reference retrieval employing a widely used, verified, accurate, yet computationally slower radiative transfer modeling technique. The reference retrieval, while too slow for using the complete hyperspectral measurement, will allow us to study the behavior of retrieval products and help us verify results from our in-development fast retrieval. Both retrievals will use the optimal estimation retrieval framework. In this manuscript we present results from an uncertainty analysis considering three uncertainty sources for a cirrus cloud retrieval in the form of error covariance matrices: reflectance uncertainty due to water vapor, the reflectance uncertainty due to ice crystal scattering assumptions, and the instrument measurement uncertainty. Results show that the uncertainty due to habit selection is the largest, while that due to water vapor is at most 0.6% relative to channel reflectance. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water.

Cirrus cloud↗

Creating high-precision reference gas standards of 85Kr for groundwater age-dating

Absolute gas counting (AGC) was applied to two gas blends of 85Kr in argon-methane (P10) counting gas to establish a high-precision specific activity (Bq/cm3) reference value for characterizing 85Kr detection efficiency for groundwater age dating measurements. The AGC or length-compensated technique has been utilized by the metrology community for decades and is an accepted method for developing radioactive gas standards. The AGC capability at Pacific Northwest National Laboratory (PNNL) uses a set of nine unequal-length proportional counters with precisely-measured internal volumes, and a gas loading system with high-precision pressure and temperature sensors. A series of AGC measurements were collected at multiple pressures to determine the inverse pressure relationship (1/P) for 85Kr and define a wall-effect correction that accounts for events decaying into the detector wall and not depositing sufficient energy in the gas to be detected. In addition to the wall-effect, two additional corrections were evaluated and are discussed in detail. Specifically, the threshold effect which accounts for events deposited below the analysis threshold and a detection efficiency as a function of detector volume effect that was observed during analysis. A robust uncertainty model was developed using the Guide to the expression of Uncertainty in Measurements (GUM) approach. The combination of carefully scrutinized correction factors, precise measurements of pressure, temperature and detector volume, and robust counting statistics resulted in the determination of high-precision specific activity values with 0.50% or less total combined uncertainty for two Kr-in-P10 reference gas standards (KP10) that will enable new groundwater age-dating measurements at PNNL.

Kr-85↗

High–Resolution Maps of Near–Surface Permafrost for Three Watersheds on the Seward Peninsula, Alaska Derived From Machine Learning

Permafrost soils are a critical component of the global carbon cycle and are locally important because they regulate the hydrologic flux from uplands to rivers. Furthermore, degradation of permafrost soils causes land surface subsidence, damaging infrastructure that is crucial for local communities. Regional and hemispherical maps of permafrost are too coarse to resolve distributions at a scale relevant to assessments of infrastructure stability or to illuminate geomorphic impacts of permafrost thaw. Here we train machine learning models to generate meter–scale maps of near–surface permafrost for three watersheds in the discontinuous permafrost region. The models were trained using ground truth determinations of near–surface permafrost presence from measurements of soil temperature and electrical resistivity. We trained three classifiers: extremely randomized trees (ERTr), support vector machines (SVM), and an artificial neural network (ANN). Model uncertainty was determined using k–fold cross validation, and the modeled extents of near–surface permafrost were compared to the observed extents at each site. At–a–site near–surface permafrost distributions predicted by the ERTr produced the highest accuracy (70%–90%). However, the transferability of the ERTr to the sites outside of the training data set was poor, with accuracies ranging from 50% to 77%. The SVM and ANN models had lower accuracies for at–a–site prediction (70%–83%), yet they had greater accuracy when transferred to the non–training site (62%–78%). These models demonstrate the potential for integrating high–resolution spatial data and machine learning models to develop maps of near–surface permafrost extent at resolutions fine enough to assess infrastructure vulnerability and landscape morphology influenced by permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗

FY2022 Status Update: A Probabilistic Model for Stress Corrosion Cracking of SNF Dry Storage Canisters

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗