Search NASA⌕ Search

SEARCH · Search NASA

Results for “chain distribution”

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 181 records · Page 10

Population balance models for polymer upcycling: signatures of the mechanism in the molecular weight evolution

Chemical and catalytic upcycling processes could help realize a circular plastics economy, but current models for testing mechanistic hypotheses and designing catalysts remain primitive. This work shows how proposed catalytic mechanisms can be incorporated into population balance models to predict the time evolution of molecular weight distributions. We develop models for homogeneous and heterogeneous catalysts, including catalysts that cut at chain ends and catalysts that cut at random locations. For heterogeneous catalysts, we illustrate the effect of adsorption constants that depend on polymer chain length. Here, we discuss ongoing efforts and challenges in measuring and modeling the time evolving molecular weight distributions in polymer upcycling processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of the feedstock supply for siting single- and multiple-feedstock biorefineries in the USA and identification of prevalent feedstocks

An integrated multi-feedstock bioenergy (i.e., biofuel, biopower, or bioproduct) supply system has potential to reduce biomass supply system uncertainties and costs. This study identifies optimal configurations of multi-feedstock biomass-to-biorefinery supply chains and pertinent feedstock combinations based on spatial distribution of feedstock and lowest delivered cost to the biorefinery. We used the Supply Characterization Model (SCM) to allocate feedstock supplies to candidate biorefinery facilities. Model runs were performed for herbaceous energy crops, agriculture residue, and woody biomass available in 2017, 2022, 2025, and 2030 as estimated by the Policy Analysis System (POLYSYS) and Forest Sustainable and Economic Analysis Model (ForSEAM) models. Three feedstock supply scenarios were compared: (a) an herbaceous scenario: switchgrass, miscanthus, biosorghum, and corn stover; (b) a woody scenario: coppice wood, noncoppice wood, whole trees, and forestry residues, and (c) a mixed scenario: a combination of all feedstocks in herbaceous and woody scenarios. By 2030 the analyses predicted that 323, 168, and 473 biorefineries were sited in the herbaceous, woody, and mixed scenario, respectively, in the conterminous USA. Feedstock mixes supplied to the biorefineries were mostly dominated by a single feedstock. The most prominent feedstock mixes identified were: (1) switchgrass and miscanthus; (2) coppice and noncoppice wood; and (3) coppice wood, noncoppice wood, switchgrass and miscanthus. Biorefineries using multi-feedstock would be beneficial for growth of bioeconomy, however flexible and cost-effective conversion platforms should be developed to efficiently utilize multiple feedstocks. This analysis identifies biorefinery locations and feedstock supply mixes while minimizing delivered feedstock costs based on spatial and temporal feedstock availability.

09 BIOMASS FUELS↗

Direct Structural Evidence for Interfacial Gradients in Asymmetric Polymer Nanocomposite Blends

Understanding the complex structure of polymer blends filled with nanoparticles (NPs) is key to design their macroscopic properties. In this work, the spatial distribution of hydrogenated (H) and deuterated (D) polymer chains asymmetric in mass is studied by small-angle neutron scattering. Depending on the chain mass, a qualitatively new large-scale organization of poly(vinyl acetate) chains beyond the random-phase approximation is evidenced in nanocomposites with attractive polymer–silica interactions. The silica is found to systematically induce bulk segregation. Only with long H-chains, a strong scattering signature is observed in the q range of the NP size: it is the sign of interfacial isotopic enrichment, that is, of contrasted polymer shells close to the NP surface. A quantitative model describing both the bulk segregation and the interfacial gradient (over ca. 10–20 nm depending on the NP size) is developed, showing that both are of comparable strength. In all cases, NP surfaces trap the polymer blend in a non-equilibrium state, with preferential adsorption around NPs only if the chain length and isotopic preference toward the surface combine their entropic and enthalpic driving forces. This structural evidence for interfacial polymer gradients will open the road for quantitative understanding of the dynamics of many-chain nanocomposite systems.

36 MATERIALS SCIENCE↗

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. Furthermore, the combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

54 ENVIRONMENTAL SCIENCES↗

Detecting outbreaks using a spatial latent field

In this paper, we present a method for estimating the infection-rate of a disease as a spatial-temporal field. Our data comprises time-series case-counts of symptomatic patients in various areal units of a region. We extend an epidemiological model, originally designed for a single areal unit, to accommodate multiple units. The field estimation is framed within a Bayesian context, utilizing a parameterized Gaussian random field as a spatial prior. We apply an adaptive Markov chain Monte Carlo method to sample the posterior distribution of the model parameters condition on COVID-19 case-count data from three adjacent counties in New Mexico, USA. Our results suggest that the correlation between epidemiological dynamics in neighboring regions helps regularize estimations in areas with high variance (i.e., poor quality) data. Using the calibrated epidemic model, we forecast the infection-rate over each areal unit and develop a simple anomaly detector to signal new epidemic waves. Our findings show that anomaly detector based on estimated infection-rates outperforms a conventional algorithm that relies solely on case-counts.

Safta, Cosmin [Sandia National Laboratories (SNL-C↗

Scalable Algorithms for Inverse Problems With High-Dimensional Parameter Spaces

Inverse problems, which involve inferring unknown parameters from observed data, present significant computational challenges, especially in large-scale settings with high-dimensional unknown parameters and nonlinear relationships between the unknowns and observations. Bayesian inference provides an approach for addressing these problems, often relying on sequential sampling methods like Markov chain Monte Carlo (MCMC) to approximate the posterior distribution of the parameters. However, MCMC methods become computationally demanding as the dimensionality of the problem increases, particularly in large-scale systems where likelihood evaluations rely on solving partial differential equations (PDEs) on large spatial domains with finely resolved meshes. To overcome these limitations, recent advancements have focused on designing scalable computa tional techniques – for both PDE simulations and sampling strategies – to make Bayesian methods feasible for high-dimensional problems.

97 MATHEMATICS AND COMPUTING↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Insight into the Competitive Adsorption Behavior of Polymer Chains in Silica Nanopores by Small-Angle Neutron Scattering

Processive hydrogenolysis catalysts, in which a metal nanoparticle (e.g., Pt) embedded at the bottom of a cylindrical silica nanopore can repeatedly cleave polymer chains and produce value-added hydrocarbon products, offer a potential solution for billions of tons of waste plastics. As the chain stays longer near the Pt catalysts, it would have a higher chance of getting cut, and therefore the molecular weight distribution of the product could be affected by the adsorption behavior of virgin chains (long polymers) versus cleaved chains (short polymers) into the nanopores. Further, this work reports a model study to understand the competitive adsorption behavior of the two different molecular weight polymers that are mixed, mimicking the reaction medium in the intermediate stage of the catalytic reaction. This study employs small-angle neutron scattering (SANS), which takes advantage of contrast differences between hydrogenous and deuterated polystyrenes to experimentally observe the relative composition of the two polymers in the silica nanopores. Our results reveal preferential adsorption of longer chains in the silica nanopores, which is consistent with the theoretical prediction in the literature for the case of the enthalpic attraction between polymers and pore walls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Life Cycle Analysis of Natural Gas Supply Chain and End Use Applications in the United States

Natural gas (NG) plays a crucial role in current and future energy systems in the United States due to its abundance and affordability. In this study a life cycle analysis of the NG supply chain in the United States was conducted using Argonne's R&D GREET model, examining stages from recovery to distribution using reported field data processed and documented by National Energy Technology Laboratory. Supply chain emissions were evaluated across multiple spatial scales, including national average, overall regional production, region-to-region, and basin-to-region scenarios. The GHG intensity of the U.S. average NG supply chain was estimated at 10.3 kg CO 2 e/MMBtu (lower heating value), with a range across regions from 7.8 kg CO 2 e/MMBtu (Northeast) to 15.1 kg CO 2 e/MMBtu (Pacific). The analysis further assessed how upstream NG emissions influence the life cycle GHG emissions of key end-use applications, including electricity generation (0.044–0.086 kg CO 2 e/kWh from upstream NG in combined cycle facilities), hydrogen production (1.04–2.20 kg CO 2 e/kg H 2 for steam methane reforming [SMR] and 1.06–2.23 kg CO 2 e/kg for autothermal reforming [ATR]), and transit bus operation utilizing compressed natural gas fuel (0.19–0.37 kg CO 2 e/mile) and hydrogen fuel (0.12–0.25 kg CO 2 e/mile for hydrogen produced in SMR and ATR).

compression↗

Many-body Nuclear Dynamics

The principal goal of this work was to better understand the fusion of neutron-rich nuclei a topic relevant to the fields of both nuclear physics and nuclear astrophysics. The work provides insight into the structure and reactions of neutron-rich nuclei namely the extent of their neutron density distribution and its polarizability. By comparing the fusion excitation functions for a chain of isotopes with a common target nucleus changes in the attractive nuclear potential are assessed. This change in the attractive potential is related to changes in the neutron density distribution with increasing number of neutrons or changes in fusion dynamics with increasing neutron number. The impact of the pairing of valence neutrons and protons on the fusion cross-section is also examined. Measurement of an isotopic chain is a powerful tool to address this topic. The experimental program made use of several different accelerator facilities. The core of the experimental work involved experiments at the ReA3 accelerator and the Facility for Rare Isotope Beams (FRIB) situated at Michigan State University, at GANIL, the French national nuclear physics laboratory, and the University of Notre Dame. At ReA3 the degree to which α-clusters associated with fusion of 28,30,32 Si + 28 Si result from the collision dynamics or reflect an initial α-cluster structure was explored. Alpha clusters observed in fusion reactions exceed the predictions of the standard statistical model. This experiment provides a measure of how clusterization is impacted by the increasing neutron-richess of the system. In the experiment at GANIL we investigated fusion in 19 O + 12 C and 20 O + 12 C. This experiment utilized the recently developed active-target detector MuSIC@Indiana. Using this proven, highly efficient, active-target detector we measured the fusion excitation function for these neutron-rich systems. Grounded by these experimental measurements, theoretical models were used to examine the role of unpaired valence neutrons on the fusion cross-section at energies just above the fusion barrier.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting

Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.

36 MATERIALS SCIENCE↗

In Pursuit of Abundance Anomalies and the FIP Effect in Late-Type Stellar Coronae

Key EUVE and ASCA data have been retrieved from their respective archives. New software has been written in the IDL language to carry out data analysis and to interface with the relevant atomic physics databases. During the analysis of ASCA spectra, it was found that the abundances of elements other than Fe could not be constrained very well, and Fe abundances were not constrained unless the underlying emission measure distribution model was reasonably well-determined. Consequently, the study has concentrated on the quantity Fe/H. A method has been developed as a means of deriving Fe/H based on fitting the continuum to EUVE spectra, thereby using the Fe lines to determine the Fe abundance. A Monte Carlo Markov Chain algorithm was also developed to determine the emission measure distribution based on observed spectral lines. This is the first application of this type of monte carlo approach to this scientific problem. This work has resulted in three scientific publications, one of which published, one of which is now ready for submission and the other of which is still in preparation:

Drake, Jeremy↗

Viscoelastic Response of Dispersed Entangled Polymer Melts

Any polymer synthesis routes result in molecular weight dispersity DAI, depending on the polymerization mechanism. The lowest dispersity polymers are those made by anionic and atom-transfer radical polymerization, which exhibit relatively narrow distributions DM=g,„1A172 — 1.02-1.04. This small divergence from monodispersity results in significant number of molecules that differ in their molecular weight from the average and their impact on the viscoelastic response remains an open question. Here the effects of low dispersity on stress relaxation and shear viscosity of entangled melts are studied using a coarse-grained model for polyethylene. Polymer melts with chain length dispersity set to follow a Schulz-Zimm distribution of Dm = 1.0 — 1.16 for an entangled polyethylene melt using molecular dynamic simulations. The systems were studied to times up to 800 ,us which is beyond the terminal time. These systems are compared to those with binary and ternary distributions. The stress relaxation functions were extracted from the Green-Kubo relation and from stress-strain relaxation following a uniaxial extension. We find on the entanglement time scale, both the mean squared displacement and the stress relaxation are independent of D m . The tube diameter and the rubbery plateau for entangled dynamics do not depend on D m . At longer times, the terminal relaxation time decreases as the faster motion of the shorter chains results in constraint release for the longer chains. Probing this low dispersity regime opens the way to evaluate the degree of dispersity that affects mechanical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

[Using Molecular Simulations to Understand Complex Nanoscale Dynamic Phenomena in Polymer Solutions]

The first half of the project concentrated on molecular simulation studies of the translocation of model molecules for single-stranded DNA through a nanosized pore. This has resulted in the publication, Translocation of a polymer chain across a nanopore: A Brownian dynamics simulation study, by Pu Tian and Grant D. Smith, JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 21 1 DECEMBER 2003, which is attached to this report. In this work we carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient) designed to mimic an electrostatic field. The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient). We focused on the latter case in our studies. Calculation of radius of gyration of the translocating chain at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tube-like pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied. Attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.

Smith, Grant↗

Translocation of a Polymer Chain across a Nanopore: A Brownian Dynamics Simulation Study

We carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient). The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient across the wall), we focused on the latter case in our studies. Calculation of radius of gyrations at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tubelike pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied, attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.

Tian, Pu↗

Advanced Precipitation Radar Antenna to Measure Rainfall From Space

To support NASA s planned 20-year mission to provide sustained global precipitation measurement (EOS-9 Global Precipitation Measurement (GPM)), a deployable antenna has been explored with an inflatable thin-membrane structure. This design uses a 5.3 5.3-m inflatable parabolic reflector with the electronically scanned, dual-frequency phased array feeds to provide improved rainfall measurements at 2.0-km horizontal resolution over a cross-track scan range of up to 37 , necessary for resolving intense, isolated storm cells and for reducing the beam-filling and spatial sampling errors. The two matched radar beams at the two frequencies (Ku and Ka bands) will allow unambiguous retrieval of the parameters in raindrop size distribution. The antenna is inflatable, using rigidizable booms, deployable chain-link supports with prescribed curvatures, a smooth, thin-membrane reflecting surface, and an offset feed technique to achieve the precision surface tolerance (0.2 mm RMS) for meeting the low-sidelobe requirement. The cylindrical parabolic offset-feed reflector augmented with two linear phased array feeds achieves dual-frequency shared-aperture with wide-angle beam scanning and very low sidelobe level of -30 dB. Very long Ku and Ka band microstrip feed arrays incorporating a combination of parallel and series power divider lines with cosine-over-pedestal distribution also augment the sidelobe level and beam scan. This design reduces antenna mass and launch vehicle stowage volume. The Ku and Ka band feed arrays are needed to achieve the required cross-track beam scanning. To demonstrate the inflatable cylindrical reflector with two linear polarizations (V and H), and two beam directions (0deg and 30deg), each frequency band has four individual microstrip array designs. The Ku-band array has a total of 166x2 elements and the Ka-band has 166x4 elements with both bands having element spacing about 0.65 lambda(sub 0). The cylindrical reflector with offset linear array feeds reduces the complexity from "NxN" transmit/receive (T/R) modules of a conventional planar-phased array to just "N" T/R modules. The antenna uses T/R modules with electronic phase-shifters for beam steering. The offset reflector does not provide poor cross-polarization like a double- curved offset reflector would, and it allows the wide scan angle in one plane required by the mission. Also, the cylindrical reflector with two linear array feeds provides dual-frequency performance with a single, shared aperture. The aperture comprises a reflective surface with a focal length of 1.89 m and is made from aluminized Kapton film. The reflective surface is of uniform thickness in the range of a few thousandths of an inch and is attached to the chain-link support structure via an adjustable suspension system. The film aperture rolls up, together with the chain-link structure, for launch and can be deployed in space by the deployment of the chain-link structure.

Rahmat-Samii, Yahya↗

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↗