Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “importance sampling”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mycorrhiza Better Predict Soil Fungal Community Composition and Function than Aboveground Traits in Temperate Forest Ecosystems

Forests in the northeastern US are experiencing shifts in community composition due to the northward migration of warm-adapted tree species and certain species’ declines (for example, white ash and eastern hemlock) due to invasive insects. Changes in belowground fungal communities and associated functions will inevitably follow. Therefore, we sought to investigate the relative importance of two important tree characteristics—mycorrhizal type [ectomycorrhizal (EcM) or arbuscular mycorrhizal (AM)] and leaf habit (deciduous or evergreen) on soil fungal community composition and organic matter cycling. We sampled soil in the organic and mineral horizons beneath two AM-associated (Fraxinus americana and Thuja occidentalis) and two ECM-associated tree species (Betula alleghaniensis and Tsuga canadensis), with an evergreen and deciduous species in each mycorrhizal group. To characterize fungal communities and organic matter decomposition beneath each tree species, we sequenced the ITS1 region of fungal DNA and measured the potential activity of carbon- and nitrogen-targeting extracellular enzymes. Each tree species harbored distinct fungal communities, supporting the need to consider both mycorrhizal type and leaf habit. However, between tree characteristics, mycorrhizal type better predicted fungal communities. Across fungal guilds, saprotrophic fungi were the most important group in shaping fungal community differences in soils beneath all tree species. The effect of leaf habit on carbon- and nitrogen-targeting hydrolytic enzymes depended on tree mycorrhizal association in the organic horizon, while oxidative enzyme activities were higher beneath EcM-associated trees across both soil horizons and leaf habits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A compact dication source for Ba 2+ tagging and heavy metal ion sensor development

Here, we present a tunable metal ion beam that delivers controllable ion currents in the picoamp range for testing of dry-phase ion sensors. Ion beams are formed by sequential atomic evaporation and single or multiple electron impact ionization, followed by acceleration into a sensing region. Controllability of the ionic charge state is achieved through tuning of electrode potentials that influence the retention time in the ionization region. Barium, lead, and cadmium samples have been used to test the system, with ion currents identified and quantified using a quadrupole mass analyzer. Realization of a clean Ba 2+ ion beam within a bench-top system represents an important technical advance toward the development and characterization of barium tagging systems for neutrinoless double beta decay searches in xenon gas. This system also provides a testbed for investigation of novel ion sensing methodologies for environmental assay applications, with dication beams of Pb 2+ and Cd 2+ also demonstrated for this purpose.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Point-of-use filtration units as drinking water distribution system sentinels

Abstract Municipal drinking water distribution systems (DWDSs) and associated premise plumbing (PP) systems are vulnerable to proliferation of opportunistic pathogens, even when chemical disinfection residuals are present, thus presenting a public health risk. Monitoring the structure of microbial communities of drinking water is challenging because of limited continuous access to faucets, pipes, and storage tanks. We propose a scalable household sampling method, which uses spent activated carbon and reverse osmosis (RO) membrane point-of-use (POU) filters to evaluate mid- to long-term occurrence of microorganisms in PP systems that are relevant to consumer exposure. As a proof of concept, POU filter microbiomes were collected from four different locations and analyzed with 16S rRNA gene amplicon sequencing. The analyses revealed distinct microbial communities, with occasional detection of potential pathogens. The findings highlight the importance of local, and if possible, continuous monitoring within and across distribution systems. The continuous operation of POU filters offers an advantage in capturing species that may be missed by instantaneous sampling methods. We suggest that water utilities, public institutions, and regulatory agencies take advantage of end-of-life POU filters for microbial monitoring. This approach can be easily implemented to ensure drinking water safety, especially from microbes of emerging concerns; e.g., pathogenic Legionella and Mycobacterium species.

42 ENGINEERING↗
DOE OSTI2025

Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation (emulator) used in the mocks. Importantly, we find a close agreement (≤ 8% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, making it particularly relevant to the broad current and future analyses of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Machine-to-machine variability of roughness and corrosion in additively manufactured 316L stainless steel

Numerous studies on the corrosion response of metal AM have been conducted. Nonetheless, the specimens being tested are commonly ground or polished to remove the outer as-built surface. If metal AM is truly going to be employed as a transformative technology that can produce complex shapes that do not require traditional machining, then the material needs to be evaluated in the as-built state. The reality is that AM alloys have shown significant inconsistencies regarding as-built surface texture, topology, and residual stress. One metric that has shown significant unpredictability in the literature is the susceptibility to localized corrosion of typically passive alloys, such as stainless steel (SS). There are a large number of studies that have attempted to understand the corrosion response of metal AM materials, but these studies are typically performed on materials that have been printed on a single machine and often mechanically polished to a smooth finish. This study compares the corrosion response of as-built AM, laser-beam powder bed fusion (LB-PBF), 316L SS parts that have been fabricated on five different machines. The majority of this work focused on understanding the susceptibility to localized corrosion of AM metals with respect to machine-dependent variables, namely surface roughness, and build angle. Surface roughness data was collected using scanning white light triangulation, laser scanning confocal microscopy, and coherence scanning interferometry. The results show that there is significant variability (p two-tail < 0.05) in the susceptibility to local corrosion initiation of LB-PBF 316L SS samples built on different machines. Surface oxides were probed with electron dispersive spectroscopy and revealed that variations in local corrosion susceptibility likely arise from differences in the stability of the passive film caused by chemical segregation, unique microstructure, and tortuous roughness features at the as-built surfaces. The variability of roughness and corrosion properties from test samples printed on different machines was corroborated by property measurements performed at five different testing sites, proving reproducibility of the data. Most importantly, this study shows that if the as-built surface layer is removed through grinding or electropolishing the machine-to-machine variation observed in the corrosion susceptibility is reduced.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Benchmark to Test Generalization Capabilities of Deep Learning Methods to Classify Severe Convective Storms in a Changing Climate

Abstract This is a test case study assessing the ability of deep learning methods to generalize to a future climate (end of 21st century) when trained to classify thunderstorms in model output representative of the present‐day climate. A convolutional neural network (CNN) was trained to classify strongly rotating thunderstorms from a current climate created using the Weather Research and Forecasting model at high‐resolution, then evaluated against thunderstorms from a future climate and found to perform with skill and comparatively in both climates. Despite training with labels derived from a threshold value of a severe thunderstorm diagnostic (updraft helicity), which was not used as an input attribute, the CNN learned physical characteristics of organized convection and environments that are not captured by the diagnostic heuristic. Physical features were not prescribed but rather learned from the data, such as the importance of dry air at mid‐levels for intense thunderstorm development when low‐level moisture is present (i.e., convective available potential energy). Explanation techniques also revealed that thunderstorms classified as strongly rotating are associated with learned rotation signatures. Results show that the creation of synthetic data with ground truth is a viable alternative to human‐labeled data and that a CNN is able to generalize a target using learned features that would be difficult to encode due to spatial complexity. Most importantly, results from this study show that deep learning is capable of generalizing to future climate extremes and can exhibit out‐of‐sample robustness with hyperparameter tuning in certain applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Better Method to Calculate Fuel Burnup in Pebble Bed Reactors Using Machine Learning

Burnup measurement is an important step in material control and accountancy (MC&A) at nuclear reactors, and may be done by examining gamma spectra of fuel samples. Traditional approaches rely on known correlations to specific photopeaks (e.g. 137 Cs) and operate via a standard linear regression method. However, the quality of these regression methods is limited even in the best case, and is significantly poorer at short fuel cool-down times, due to the elevated radiation background by short life-time isotopes, and self-shielding effect of the fuel. For practical operation of pebble bed reactors (PBRs), quick measurements (in minutes) and short cooling times (in hours) are required from a safety and security perspective. We investigated the efficacy and performance of machine learning (ML) methods to predict the burnup of the pebble fuel from full gamma spectra (rather than specific discrete photopeaks) and found a full-spectrum ML approach to far outperform baseline regression predictions in all measurement and cooling conditions - including in operational-like measurement conditions. We also performed model and data ablation experiments to determine the relative performance impact of our ML methods' capacity to model data nonlinearities and the inherent additional information in full spectra. Applying our ML methods, we found a number of surprising results, including improved accuracy at shorter fuel cooling times (the opposite of the norm), remarkable robustness to spectrum compression (via rebinning), and competitive burnup predictions even when using background signal only (i.e. explicitly omitting known isotope photopeaks).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Interplay between magnetism and superconductivity in UTe 2

Time-reversal symmetry breaking (TRSB) in UTe 2 was inferred from observations of a spontaneous Kerr response in the superconducting state after cooling in zero magnetic field, while a finite c-axis magnetic field training was further used to determine the nature of the nonunitary composite order parameter of this material. Here, we present an extensive study of the magnetic-field-trained Kerr effect, which unveils a unique critical state of pinned “ferromagnetic vortices.” We show that a remanent Kerr signal appears following the removal of a training magnetic field, which reflects the response to the TRSB order parameter and the external magnetic field through the paramagnetic susceptibility. Furthermore, this unambiguously demonstrates the importance of the magnetic fluctuations and their intimate relation to the composite order parameter. Focusing the beam on the center of the sample, we are able to accurately determine the maximum field that is screened by the critical state and the respective critical current. Measurements in the presence of magnetic field show the tendency of the superconductor to produce shielding currents that oppose the increase in vortex-induced magnetization due to the diverging paramagnetic susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Optimal Twirling Depth for Classical Shadows in the Presence of Noise

The classical shadows protocol is an efficient strategy for estimating properties of an unknown state p using a small number of state copies and measurements. In its original form, it involves twirling the state with unitaries from some ensemble and measuring the twirled state in a fixed basis. It was recently shown that for computing local properties, optimal sample complexity (copies of the state required) is remarkably achieved for unitaries drawn from shallow depth circuits composed of local entangling gates, as opposed to purely local (zero depth) or global twirling (infinite depth) ensembles. Here, we consider the sample complexity as a function of the depth of the circuit, in the presence of noise. We find that this noise has important implications for determining the optimal twirling ensemble. Under fairly general conditions, we (i) show that any single-site noise can be accounted for using a depolarizing noise channel with an appropriate damping parameter f, (ii) compute thresholds f th at which optimal twirling reduces to local twirling for Pauli operators, (iii) nth order Renyi entropies (n ≥2), and (iv) provide a meaningful upper bound t max on the optimal circuit depth for any finite noise strength f, which applies to observables and entanglement entropy measurements. In conclusion, these thresholds strongly constrain the search for optimal strategies to implement shadow tomography and are easily tailored to the experimental system at hand.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Bottomonium suppression in pNRQCD and open quantum system approach

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sample Preparation for 3D Characterization of Irradiated Fuel

Thorough characterization of nuclear fuel is essential when striving to understand its microstructural evolution in response to varying irradiation conditions. A detailed understanding of the microstructural evolution of irradiated fuel can be useful when predicting the fuel’s response to off-normal conditions and related phenomena such as high burnup fuel fragmentation (HBFF) and fission gas release. This sort of characterization and the study of complex phenomena like HBFF that are related to fuel performance are an important part of the US Department of Energy Office of Nuclear Energy’s Advanced Fuels Campaign (AFC). Multiple characterization techniques are commonly used to perform post-irradiation examination (PIE) analysis on an irradiated fuel sample, but many of these techniques are limited to the 2D surface analysis of the fuel pellet. Although there is a lot of useful information to be gained by analyzing the fuel pellet surface, a 2D analysis provides only an approximation of the nature of the material’s features, such as porosity and grain structure. Therefore, it is critical to analyze the 3D structure of the material, as materials exist and respond to conditions in all three dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Evaluating Limits of Machine Learning-Assisted Raman Spectroscopy in Classification of Biological Samples

Machine learning (ML)-assisted Raman spectroscopy has become a powerful analytical tool for the classification and identification of analytes; however, technical challenges impacting its detection accuracy have not been thoroughly investigated. This study explores experimental factors affecting classification performance. Among the evaluated ML models, ML algorithms show minimal impact on classification accuracy. Instead, experimental factors, including spectral similarity between tested samples and data quality, dominate detection performance. Increases in spectral noise and spectral similarity significantly reduce classification accuracy. In well-controlled samples with low experimental noise, ML-assisted Raman spectroscopy can discriminate lipid mixtures with a composition difference of 1.85 mol %. To assess the effect of biological heterogeneity, we analyzed single-cell Raman spectra from Saccharomyces cerevisiae strains carrying single, double, or triple gene mutations. Intrinsic cell-to-cell variability introduced substantial spectral differences, severely reducing the accuracy of multiclass classification of these genetically similar strains at the single-cell level. Averaging Raman spectra across multiple cells improved classification accuracy by reducing this spectral variability. We also assess the effectiveness of transfer learning across different Raman spectrometers, specifically by applying an ML model trained on one instrument to another Raman spectrometer. Transfer learning can be improved with proper instrument calibration, highlighting the importance of instrument standardization. Overall, our results demonstrate that data quality and spectral similarity are the primary bottlenecks in ML-assisted Raman spectroscopy. Careful attention to sample preparation, data acquisition, measurement conditions, and instrument calibration is critical to achieving robust and reliable classification performance.

Fungi↗
DOE OSTI2023

Bootstrapping outperforms community‐weighted approaches for estimating the shapes of phenotypic distributions

Abstract Estimating phenotypic distributions of populations and communities is central to many questions in ecology and evolution. These distributions can be characterized by their moments (mean, variance, skewness and kurtosis) or diversity metrics (e.g. functional richness). Typically, such moments and metrics are calculated using community‐weighted approaches (e.g. abundance‐weighted mean). We propose an alternative bootstrapping approach that allows flexibility in trait sampling and explicit incorporation of intraspecific variation, and show that this approach significantly improves estimation while allowing us to quantify uncertainty. We assess the performance of different approaches for estimating the moments of trait distributions across various sampling scenarios, taxa and datasets by comparing estimates derived from simulated samples with the true values calculated from full datasets. Simulations differ in sampling intensity (individuals per species), sampling biases (abundance, size), trait data source (local vs. global) and estimation method (two types of community‐weighting, two types of bootstrapping). We introduce the traitstrap R package, which contains a modular and extensible set of bootstrapping and weighted‐averaging functions that use community composition and trait data to estimate the moments of community trait distributions with their uncertainty. Importantly, the first function in the workflow, trait_fill , allows the user to specify hierarchical structures (e.g. plot within site, experiment vs. control, species within genus) to assign trait values to each taxon in each community sample. Across all taxa, simulations and metrics, bootstrapping approaches were more accurate and less biased than community‐weighted approaches. With bootstrapping, a sample size of 9 or more measurements per species per trait generally included the true mean within the 95% CI. It reduced average percent errors by 26%–74% relative to community‐weighting. Random sampling across all species outperformed both size‐ and abundance‐biased sampling. Our results suggest randomly sampling ~9 individuals per sampling unit and species, covering all species in the community and analysing the data using nonparametric bootstrapping generally enable reliable inference on trait distributions, including the central moments, of communities. By providing better estimates of community trait distributions, bootstrapping approaches can improve our ability to link traits to both the processes that generate them and their effects on ecosystems.

Maitner, Brian S.↗
DOE OSTI2023

Arctic Permafrost Thawing Enhances Sulfide Oxidation

Permafrost degradation is altering biogeochemical processes throughout the Arctic. Thaw-induced changes in organic matter transformations and mineral weathering reactions are impacting fluxes of inorganic carbon (IC) and alkalinity (ALK) in Arctic rivers. However, the net impact of these changing fluxes on the concentration of carbon dioxide in the atmosphere (pCO 2 ) is relatively unconstrained. Resolving this uncertainty is important as thaw-driven changes in the fluxes of IC and ALK could produce feedbacks in the global carbon cycle. Enhanced production of sulfuric acid through sulfide oxidation is particularly poorly quantified despite its potential to remove ALK from the ocean-atmosphere system and increase pCO 2 , producing a positive feedback leading to more warming and permafrost degradation. In this work, we quantified weathering in the Koyukuk River, a major tributary of the Yukon River draining discontinuous permafrost in central Alaska, based on water and sediment samples collected near the village of Huslia in summer 2018. Using measurements of major ion abundances and sulfate (SO 4 2- ) sulfur ( 34 S/ 32 S) and oxygen ( 18 O/ 16 O) isotope ratios, we employed the MEANDIR inversion model to quantify the relative importance of a suite of weathering processes and their net impact on pCO 2 . Calculations found that approximately 80% of SO 4 2- in mainstem samples derived from sulfide oxidation with the remainder from evaporite dissolution. Moreover, 34 S/ 32 S ratios, 13 C/ 12 C ratios of dissolved IC, and sulfur X-ray absorption spectra of mainstem, secondary channel, and floodplain pore fluid and sediment samples revealed modest degrees of microbial sulfate reduction within the floodplain. Weathering fluxes of ALK and IC result in lower values of pCO 2 over timescales shorter than carbonate compensation (~10 4 yr) and, for mainstem samples, higher values of pCO 2 over timescales longer than carbonate compensation but shorter than the residence time of marine (~10 7 yr). Furthermore, the absolute concentrations of and Mg 2+ in the Koyukuk River, as well as the ratios of SO 4 2- and Mg 2+ to other dissolved weathering products, have increased over the past 50 years. Through analogy to similar trends in the Yukon River, we interpret these changes as reflecting enhanced sulfide oxidation due to ongoing exposure of previously frozen sediment and changes in the contributions of shallow and deep flow paths to the active channel. Overall, these findings confirm that sulfide oxidation is a substantial outcome of permafrost degradation and that the sulfur cycle responds to permafrost thaw with a timescale-dependent feedback on warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Elucidating Processes Controlling Arctic Atmospheric Aerosol Sources, Aging, and Mixing States (Final Report)

Atmospheric aerosols play critical roles in the Earth’s energy budget, directly by scattering or absorbing solar and terrestrial radiation and indirectly by serving as seeds (nuclei) for cloud droplet and ice crystal formation and by depositing on snow and ice surface, thereby changing the surface albedo. These effects are dependent on aerosol particle size and chemical composition and impact the hydrological cycle as well. This project provided single-particle size and chemical composition measurements across the entire annual cycle in the high Arctic and in the Alaskan Arctic during fall – winter, addressing the most significant gaps in Arctic aerosol observational data. These needs were based on recent rapid sea ice loss across the entire Arctic, as well as the major annual delays in sea ice freeze-up during fall in the Chukchi Sea and increased wintertime sea ice fracturing in the Beaufort Sea, both off the North Slope of Alaska. Two DOE Atmospheric Radiation Measurement (ARM) field campaigns were conducted for atmospheric aerosol sampling. The Aerosols during the Polar Utqiagvik Night (APUN – ‘snow on ground’ in Iñupiaq) ARM field campaign at Utqagivik, Alaska was conducted from Oct. 28 – Dec. 22, 2018. Aerosol sizing instrumentation and a single-particle mass spectrometer were successfully deployed for size-resolved number concentration measurements and measurements of individual particle size and chemical composition, respectively. These results show the influence of locally-produced sea spray aerosol, with high cloud-forming potential, due to delayed sea ice freeze-up in the fall. During the 2019‐2020 international Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, daily atmospheric aerosol particles were collected aboard the German icebreaker Polarstern in the Central Arctic from Nov. 2019 – Oct. 2020. Sea salt aerosol and marine organics were observed year-round during MOSAiC with varying morphologies and sources. These findings are important because most Arctic models do not include a sea spray aerosol source, despite this source increasing with declining sea ice extent. In addition to collecting new samples and data, this project also conducted further analysis of previously collected single-particle chemical composition measurements within the North Slope of Alaska oil fields and at Utqiaġvik, AK, during Aug. – Sep. 2015 and 2016 field campaigns. This work resulted in the discovery of chemical reactions of oil field combustion emissions occurring within fog droplets across the North Slope of Alaska and forming secondary aerosol, showing the impact of Arctic oil field emissions beyond black carbon aerosol and greenhouse gases. In addition, the distribution of chemical species across the aerosol population within the oil fields was quantified, using these data and a previously development framework. We also presented the first ambient evidence of the collision of two atmospheric particles resulting in formation of an organic-coated ammonium sulfate particle of marine origin, which has implications for cloud formation with declining sea ice extent. Overall, this project has elucidated connections between seawater biogeochemistry, resource extraction activities, atmospheric composition, clouds, and the energy budget of the Arctic region. The results of this project are expected to improve weather and sea ice forecasting for security and development in the Arctic and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sampling Size Optimization for Bioburden Density Estimation in Planetary Protection

Planetary protection (PP) is a discipline that focuses on minimizing the biological contamination of spacecraft to ensure compliance with international policy. Precise estimation of bioburden - the total number of microbes in or on spacecraft hardware – and the bioburden density are of utmost importance for PP. Such estimation is the way concordance with requirements is demonstrated, and it is critical for quantifying the potential risk of inadvertently contaminating other planetary bodies. Although a suite of molecular techniques have been used to thoroughly characterize and profile the microbiome of various cleanroom environments and spacecraft, the gold standard remains the physical enumeration of microbes via culturing of samples directly taken from spacecraft and associated surfaces. However, due to technical, budgetary, and programmatic constraints, only a manageable portion (around 10%) of the entire spacecraft surface is directly sampled with cotton swabs or wipes. To generate the bioburden current best estimate (CBE) for components not directly verifiable, the accepted approach is to apply a NASA-defined bioburden estimate based on the components’ manufacturing or assembly environment. This approach utilizes a prespecified bioburden density estimation that applies a maximum value across the total surface area of the specified component. For hardware components that underwent similar assembly processes, an implied bioburden is adopted for all components, based on a direct verification of a representative component within the same lot. Once all components have a CBE, the bioburden estimates are generated. In previous publication [ 1], we have shown that statistical risks quantifying the accuracy of the estimates for sampled, prespecified, and implied components can be derived and ranked. For mean squared error (MSE) function, the risks are available analytically and hence a cost function can be obtained to optimize the risks with respect to the sampling area and sampling cost. Since the sampling area and sampling cost are two complimentary variables, their sum will have a well-defined minimum. This paper presents the multivariate optimization of the integrated risk of an empirical Bayes estimator to determine the optimal sampling schedule for a given number of components. It is assumed that given a number of components, N, the bioburden density for each component can either be sampled, implied, or prespecified. The multivariate optimization searches through different options to sample, imply or prespecify the bioburden density for a component, and account for the component’s surface area and cost of sampling. The idea of the optimization is based on the observation that the statistical risk of using an estimator is a monotonically decreasing function of the sampled area. The larger the sampled area, the lower the risk of using the estimator as the estimator becomes more and more accurate as the sampling area increases. On the other hand, the cost of sampling is monotonically increasing as the sampled surface grows. This makes the risk and total cost of sampling complimentary variables which can be counterbalanced to achieve an optimal overall value with respect to the sampled surface. In this paper, the integrated risk has been used to quantify the accuracy of the estimator. This risk has been selected because it depends on neither the true value of the parameter nor on the collected data. The cost of each sample was also available to obtain the total cost of sampling of N components. The paper will present the results based on computer-simulated data as well as the data collected during the InSight mission. The computer-simulated data have N components with randomly generated total areas and each component assigned to one of the three categories according to the method of estimating of bioburden density: sampled, implied, or prespecified. The cost of sampling is also available. The cost of sampling is estimated based on a cost model provided by the planetary protection group at JPL. For this paper, the overall cost was assumed to be a linear function of exposure. The optimization process finds the allocation of the components to the three categories that minimizes the tradeoff between integrated risk and total cost. For the InSight data, a set of components is selected representing all three categories, and optimization is performed to determine if the performed allocation was optimal or if a better allocation could have been obtained. To the best of our knowledge, this work is the first attempt not only perform an accurate estimation of bioburden density but also do it in an optimal way.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Geochemical Assessment for Carbon Sequestration in the Conasauga Group, Northwest Georgia, USA

Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.

42 ENGINEERING↗