Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inference mechanisms”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗
DOE OSTI2020

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2–3 × 10 −4 (Δ E / E ). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (Δ E / E = 3 × 10 −3 ) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Analysis and Integration of the Hydraulic Fracturing Test Site-2 (HFTS-2) Comprehensive Dataset

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Permian (Delaware) Basin. The unique aspect of this program was the acquisition of a unique, comprehensive, diagnostic dataset. Additionally, shorter parent wells drilled three years before the child wells offered clear distinction between the stages influenced by parent-child effects and the stages without any effects. The goal of this study was to analyze and integrate this comprehensive diagnostic dataset to understand the areal and vertical extent of hydraulic fractures (HF). The paper also provides insights on the effects of parent wells’ depletion on child well HF geometry based on various monitoring methods and subsurface models. Areal and vertical coverage for all HFTS-2 wells during stimulation and depletion was estimated based on analysis and interpretation of diagnostics and advanced modeling results. HFTS-2 diagnostics included microseismic (MS), pre- and post-stimulation logs and cores, bottomhole gauges, and fiber optic (FO) data. The diagnostics results (MS, FO, image logs) were integrated and used to calibrate subsurface models. Additional field tests were designed and implemented for depletion monitoring. The tailored program for monitoring depletion included vertical and slant well pressures, interference testing, and a vertical strain depletion trial. Areal Coverage: Conventional MS (and FO MS) were used to compute HF dimensions, which were compared with diagnostics (FO strain, gauge, image logs) observations and calibrated subsurface models. A post-production interference test did not show offset well communication. Vertical Coverage: Vertical coverage during stimulation was monitored using a vertical monitoring well. The stronger mechanical strain signals showed good correlation with MS event intensities, geomechanical properties, and gauge inferences. Vertical depletion was estimated based on vertical/slant well gauges and strain depletion tests. Parent-Child Effects: Diagnostics and calibrated subsurface models show asymmetry in child well fracture geometries for stages that overlap parent wells. Child well image logs serve as a good indicator for parent Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R004/2477423/urtec-2021-5241-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5241 well HF tracking. Child well MS events had an eastward bias, in line with pre-stimulation image logs, and was confirmed by parent well frac hits. Novel/Additive Information: The dataset presents a unique, over-constrained problem space to compare independent techniques to arrive at HF metrics (i.e., stimulation height and/or half-length), unlike a single source dataset, in which calibration is done using available data to guide predictions. Here, the asymmetry in HF geometry seen in the stages influenced by parent-child effects offers unique insights into well spacing and landing, which are key capital decisions the unconventional resources industry is seeking to optimize.

58 GEOSCIENCES↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2024

Identifying Heterogeneous Micromechanical Properties of Biological Tissues via Physics–Informed Neural Networks

The heterogeneous micromechanical properties of biological tissues have profound implications across diverse medical and engineering domains. However, identifying full-field heterogeneous elastic properties of soft materials using traditional engineering approaches is fundamentally challenging due to difficulties in estimating local stress fields. Recently, there has been a growing interest in data-driven models for learning full-field mechanical responses, such as displacement and strain, from experimental or synthetic data. However, research studies on inferring full-field elastic properties of materials, a more challenging problem, are scarce, particularly for large deformation, hyperelastic materials. Here, a physics-informed machine learning approach is proposed to identify the elasticity map in nonlinear, large deformation hyperelastic materials. This study reports the prediction accuracies and computational efficiency of physics-informed neural networks (PINNs) in inferring the heterogeneous elasticity maps across materials with structural complexity that closely resemble real tissue microstructure, such as brain, tricuspid valve, and breast cancer tissues. Further, the improved architecture is applied to three hyperelastic constitutive models: Neo-Hookean, Mooney Rivlin, and Gent. Furthermore, the improved network architecture consistently produces accurate estimations of heterogeneous elasticity maps, even when there is up to 10% noise present in the training data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

From remotely sensed solar-induced chlorophyll fluorescence to ecosystem structure, function, and service: Part I—Harnessing theory

Solar-induced chlorophyll fluorescence (SIF) is a remotely sensed optical signal emitted during the light reactions of photosynthesis. The past two decades have witnessed an explosion in availability of SIF data at increasingly higher spatial and temporal resolutions, sparking applications in diverse research sectors (e.g., ecology, agriculture, hydrology, climate, and socioeconomics). These applications must deal with complexities caused by tremendous variations in scale and the impacts of interacting and superimposing plant physiology and three-dimensional vegetation structure on the emission and scattering of SIF. At present, these complexities have not been overcome. To advance future research, the two companion reviews aim to (1) develop an analytical framework for inferring terrestrial vegetation structures and function that are tied to SIF emission, (2) synthesize progress and identify challenges in SIF research via the lens of multi-sector applications, and (3) map out actionable solutions to tackle these challenges and offer our vision for research priorities over the next 5-10 years based on the proposed analytical framework. This paper is the first of the two companion reviews, and theory oriented. It introduces a theoretically rigorous yet practically applicable analytical framework. Guided by this framework, we offer theoretical perspectives on three overarching questions: (1) The forward (mechanism) question—How are the dynamics of SIF affected by terrestrial ecosystem structure and function? (2) The inference question: What aspects of terrestrial ecosystem structure, function, and service can be reliably inferred from remotely sensed SIF and how? (3) The innovation question: What innovations are needed to realize the full potential of SIF remote sensing for real-world applications under climate change? In conclusion, the analytical framework elucidates that process complexity must be appreciated in inferring ecosystem structure and function from the observed SIF; this framework can serve as a diagnosis and inference tool for versatile applications across diverse spatial and temporal scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Vapor-phase self-assembly for generating thermally stable single-atom catalysts

Preparation of thermally stable metal single-atom catalysts (SACs) is a challenge in heterogeneous catalysis, especially on conventional supports that provide a weak metal-support interaction. In this work, we report that a modified support MgAl 2 O 4 can stabilize Pt single atoms by a mechanism of vapor-phase self-assembly in a high-temperature treatment (800°C, air). The experimental results on the formation mechanism and the structure are validated by DFT and ab initio molecular dynamics simulations. We infer that stable triangular K 3 O 3 structures help stabilize Pt single atoms at high temperatures in oxidizing conditions, exhibiting excellent reactivity for methane oxidation. The obtained Pt/K/MgAl 2 O 4 SAC presents excellent stability in methane oxidation after steam treatment at elevated temperatures, whereas the Pt/MgAl 2 O 4 nanocatalyst suffers from rapid deactivation due to Pt nanoparticle growth. Now this work paves the way for preparing thermally stable and highly active SACs using conventional high-surface-area supports, despite the weak metal-support interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗
DOE OSTI2024

Codiscovering graphical structure and functional relationships within data: A Gaussian Process framework for connecting the dots

Most problems within and beyond the scientific domain can be framed into one of the following three levels of complexity of function approximation. Type 1: Approximate an unknown function given input/output data. Type 2: Consider a collection of variables and functions, some of which are unknown, indexed by the nodes and hyperedges of a hypergraph (a generalized graph where edges can connect more than two vertices). Given partial observations of the variables of the hypergraph (satisfying the functional dependencies imposed by its structure), approximate all the unobserved variables and unknown functions. Type 3: Expanding on Type 2, if the hypergraph structure itself is unknown, use partial observations of the variables of the hypergraph to discover its structure and approximate its unknown functions. These hypergraphs offer a natural platform for organizing, communicating, and processing computational knowledge. While most scientific problems can be framed as the data-driven discovery of unknown functions in a computational hypergraph whose structure is known (Type 2), many require the data-driven discovery of the structure (connectivity) of the hypergraph itself (Type 3). We introduce an interpretable Gaussian Process (GP) framework for such (Type 3) problems that does not require randomization of the data, access to or control over its sampling, or sparsity of the unknown functions in a known or learned basis. Its polynomial complexity, which contrasts sharply with the super-exponential complexity of causal inference methods, is enabled by the nonlinear ANOVA capabilities of GPs used as a sensing mechanism.

Science & Technology - Other Topics↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2023

A comparison of histopathology imaging comprehension algorithms based on multiple instance learning

Whole slide imaging (WSI), also called digital virtual microscopy, is a new imaging modality. It allows for the application of AI and machine learning methods to cancer pathology to help establish a means for the automatic diagnosis of cancer cases. However, designing machine-learning models for WSI is computationally challenging due to its required ultra-high resolution. The current state-of-the-art models use multiple instance learning (MIL). MIL is a weakly-supervised learning method in which the model uses an array of inferences from many smaller instances to make a final classification about the entire set. In the context of WSI, researchers divide the ultra-high-resolution image into many patches. The model then classifies the slide based on an array of inferences from the patches. Among several ways of making the final classification, attention-based mechanisms have resulted in superb accuracy scores. The Transformer, one attention-based algorithm, has reported substantial improvements for WSI comprehension tasks. In this project, we studied and compared several WSI comprehension algorithms. We used the following three datasets: CAMELYON16+17, TCGALung, and TCGA-Kidney. We found that attention-based MIL algorithms performed better than standard MIL algorithms for classifying WSI images, achieving a higher mean accuracy and AUC. However, none of the attention-based algorithms performed significantly better than the others, reporting accuracy scores that varied widely. Presumably, it is due to the limited availability of training samples in the data corpus. Since it is not easy to increase the samples from human subjects, some machine learning techniques like transfer learning could help mitigate this issue.

Saunders, Adam↗
DOE OSTI2022

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Identifying crack tip position and stress intensity factors from displacement data

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, in this work, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗
DOE OSTI2021

Long-Term 13 C Uptake by 12 C-Enriched Calcite

Knowledge of the exchange of carbon isotopes between dissolved inorganic carbon and calcite minerals is of long-standing importance for the interpretation of sedimentary paleoclimate records and 14 C transport in the geosphere. To assess the mechanism and rates of carbon isotope exchange, we equilibrated 12 C-pure synthetic calcite particles in water, first in a glovebox and then in contact with atmospheric P CO 2 and 13 C/ 12 C ratio, at two different temperatures. Cavity ring-down infrared spectroscopy δ 13 C measurements of the solid revealed sustained 13C incorporation for over a period of 500 days (21 °C) and 125 days (50 °C). We developed a quantitative model for recrystallization and isotope exchange, assuming that the interfacial free energy provides a thermodynamic driving force for the growth of larger particles at the expense of smaller ones. Furthermore, this Ostwald ripening model did not reproduce the kinetics of 13C uptake and required greater coarsening than observed. Rather, the data were best explained by a mechanism involving surface exchange and solid-phase diffusion of 13C into the particles with an inferred effective diffusion constant at 21 °C of about 10 –25 m 2 /s. Although this work cannot rule out the possibility that structural or chemical aspects of the synthetic particles enabled faster 13 C uptake than could be observed in natural systems, this study adds to the body of the recent work, suggesting that fast exchange processes are possible, likely through grain boundaries and other defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Informed Multi‐Scale Approach Applied to the British Columbia Fires of Late Summer 2017

Abstract Pyrocumulonimbus (PyroCb) clouds have a complex origin dependent on fire dynamics and meteorological conditions. When a pyrocumulonimbus cloud develops and is maintained over a period of time, it can inject significant aerosol into the troposphere and lower stratosphere, resulting in a longer‐term (months to years) occurrence of aerosol in the stratosphere. In this work, we investigate the British Columbia wildfires on 12–13 August 2017 using a multi‐scale simulation framework. We use the output of a physics‐based wildfire model (FIRETEC) with parameterized energy, particle, and gas emissions to drive the upper atmospheric aerosol mass injection within a regional cloud resolving model (HIGRAD). We demonstrate that vertical motions produced by latent heat release of the condensation of ice and cloud particles within the PyroCbs induce another 5 km of lifting of the simulated aerosol plume. Primary black carbon and organic aerosols (OAs) alone may not be enough to explain the observed aerosol burden, thus we show that secondary OA produced via condensation of gases by the fires, ash, and possibly dust can enhance lofted aerosol mass. A simulation with all emission mechanisms active, driven by the observed fuel load and environmental conditions, reasonably reproduces an aerosol profile inferred from observational data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

58 GEOSCIENCES↗