Results for “sampling variability” — Page 46

DOE OSTI2020

Predicting Biological Cleanliness: An Empirical Bayes Approach for Spacecraft Bioburden Accounting

To comply with the international planetary protection policy set forth by the Committee on Space Research and NASA Agency level requirements, spacecraft destined to biologically sensitive planetary bodies have to minimize terrestrial biological contamination. Analysis, testing and inspection are the standard forward verification activities that are used to demonstrate compliance with the biological contamination requirements. For testing of spacecraft surface areas, a swab or wipe sample is collected from surfaces prior to last access and subsequently processed in the lab using NASA Approved Planetary Protection Methods for Culture Based Assays. Raw data resulting from this assay is then statistically treated employing a mathematical paradigm stemming from the 1970’s Viking Lander Project to generate the bioburden density and total microbial bioburden present. This standard approach arbitrarily accounts for error and provides an upper conservative bound as it reports the maximum number of spores estimated to be present on flight hardware surfaces. A bioburden density estimate factors in the following variables: the observed bioburden count, representative volume processed, sampling efficiencies. Notably, to account for error in the approach, a 0 observed count is arbitrarily changed to a count of 1 for each hardware grouping. The data generated by spacecraft bioburden verification campaigns in the past have resulted in <80% of wipes and <90% of swabs containing a bioburden count of 0. As such, having a robust and well documented statistical approach for dealing with the probability of low incident rates is necessary to be able to estimate spacecraft bioburden. Being able to statistically describe the bioburden distribution and associated confidence level is a gamechanger for the development of bioburden allocations during mission design and will allow for tighter management of risk throughout spacecraft build. Thus, Empirical Bayes statistical approach was evaluated to estimate the microbial bioburden on spacecraft to mitigate the aforementioned mathematical concerns and provide a probabilistic bioburden distribution of the flight hardware surface. For application of this approach to performing bioburden calculations, a range of non-informative prior assumptions on hardware surfaces are explored for Bayesian analyses while informative priors using posterior distributions from prior assays are utilized for Empirical Bayes analyses. Several non-informative priors are currently under investigation to assess fitness including use of these priors to serve as a foundation to build off of NASA specification values or a basis of risk to account for unknowns during the integration and testing process. Informative priors under consideration are generated using sampled bioburden values from hardware originating within like processing environments (e.g. vendor cleaning process or similar assembly process), temporal spacecraft status events as a prediction for hardware cleanliness of future samples, and heritage system bioburden actuals to predict allocation for subsequent missions. Informative priors and probabilistic bioburden distributions are then validated using data sets from the Mars Exploration Rover, Mars Science Laboratory, and InSight missions. Using Empirical Bayes approach to generate a probabilistic bioburden distribution as demonstrated through mission use cases provides a valid approach for use in the end-to-end requirements verification process.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Bayesian Framework for Bioburden Density Estimation in Planetary Protection

To comply with the international planetary protection policy set forth by the Committee on Space Research and NASA Agency level requirements, spacecraft destined to biologically sensitive planetary bodies have to minimize terrestrial biological contamination. Analysis, testing and inspection are the standard forward verification activities that are used to demonstrate compliance with the biological contamination requirements. For testing of spacecraft surface areas, a swab or wipe sample is collected from surfaces prior to last access and subsequently processed in the lab using NASA Approved Planetary Protection Methods for Culture Based Assays. Raw data resulting from this assay is then statistically treated employing a mathematical paradigm stemming from the 1970’s Viking Lander Project to generate the bioburden density and total microbial bioburden present. This standard approach arbitrarily accounts for error and provides an upper conservative bound as it reports the maximum number of spores estimated to be present on flight hardware surfaces. A bioburden density estimate factors in the following variables: the observed bioburden count, representative volume processed, sampling efficiencies. Notably, to account for error in the approach, a 0 observed count is arbitrarily changed to a count of 1 for each hardware grouping. The data generated by spacecraft bioburden verification campaigns in the past have resulted in <80% of wipes and <90% of swabs containing a bioburden count of 0. As such, having a robust and well documented statistical approach for dealing with the probability of low incident rates is necessary to be able to estimate spacecraft bioburden. Being able to statistically describe the bioburden distribution and associated confidence level is a gamechanger for the development of bioburden allocations during mission design and will allow for tighter management of risk throughout spacecraft build. Thus, Empirical Bayes statistical approach was evaluated to estimate the microbial bioburden on spacecraft to mitigate the aforementioned mathematical concerns and provide a probabilistic bioburden distribution of the flight hardware surface. For application of this approach to performing bioburden calculations, a range of non-informative prior assumptions on hardware surfaces are explored for Bayesian analyses while informative priors using posterior distributions from prior assays are utilized for Empirical Bayes analyses. Several non-informative priors are currently under investigation to assess fitness including use of these priors to serve as a foundation to build off of NASA specification values or a basis of risk to account for unknowns during the integration and testing process. Informative priors under consideration are generated using sampled bioburden values from hardware originating within like processing environments (e.g. vendor cleaning process or similar assembly process), temporal spacecraft status events as a prediction for hardware cleanliness of future samples, and heritage system bioburden actuals to predict allocation for subsequent missions. Informative priors and probabilistic bioburden distributions are then validated using data sets from the Mars Exploration Rover, Mars Science Laboratory, and InSight missions. Using Empirical Bayes approach to generate a probabilistic bioburden distribution as demonstrated through mission use cases provides a valid approach for use in the end-to-end requirements verification process.

97 - MATHEMATICS AND COMPUTING↗

NASA NTRS1989

Rain attenuation measurements: Variability and data quality assessment

Year to year variations in the cumulative distributions of rain rate or rain attenuation are evident in any of the published measurements for a single propagation path that span a period of several years of observation. These variations must be described by models for the prediction of rain attenuation statistics. Now that a large measurement data base has been assembled by the International Radio Consultative Committee, the information needed to assess variability is available. On the basis of 252 sample cumulative distribution functions for the occurrence of attenuation by rain, the expected year to year variation in attenuation at a fixed probability level in the 0.1 to 0.001 percent of a year range is estimated to be 27 percent. The expected deviation from an attenuation model prediction for a single year of observations is estimated to exceed 33 percent when any of the available global rain climate model are employed to estimate the rain rate statistics. The probability distribution for the variation in attenuation or rain rate at a fixed fraction of a year is lognormal. The lognormal behavior of the variate was used to compile the statistics for variability.

Crane, Robert K.↗

NASA NTRS1997

Evaluation of Tropospheric Water Vapor Simulations from the Atmospheric Model Intercomparison Project

Simulations of humidity from 28 general circulation models for the period 1979-88 from the Atmospheric Model Intercomparison Project are compared with observations from radiosondes over North America and the globe and with satellite microwave observations over the Pacific basin. The simulations of decadal mean values of precipitable water (W) integrated over each of these regions tend to be less moist than the real atmosphere in all three cases; the median model values are approximately 5% less than the observed values. The spread among the simulations is larger over regions of high terrain, which suggests that differences in methods of resolving topographic features are important. The mean elevation of the North American continent is substantially higher in the models than is observed, which may contribute to the overall dry bias of the models over that area. The authors do not find a clear association between the mean topography of a model and its mean W simulation, however, which suggests that the bias over land is not purely a matter of orography. The seasonal cycle of W is reasonably well simulated by the models, although over North America they have a tendency to become moister more quickly in the spring than is observed. The interannual component of the variability of W is not well captured by the models over North America. Globally, the simulated W values show a signal correlated with the Southern Oscillation index but the observations do not. This discrepancy may be related to deficiencies in the radiosonde network, which does not sample the tropical ocean regions well. Overall, the interannual variability of W, as well as its climatology and mean seasonal cycle, are better described by the median of the 28 simulations than by individual members of the ensemble. Tests to learn whether simulated precipitable water, evaporation, and precipitation values may be related to aspects of model formulation yield few clear signals, although the authors find, for example, a tendency for the few models that predict boundary layer depth to have large values of evaporation and precipitation. Controlled experiments, in which aspects of model architecture are systematically varied within individual models, may be necessary to elucidate whether and how model characteristics influence simulations.

Gaffen, Dian J.↗

NASA NTRS2016

Discordant K-Ar and Young Exposure Dates for the Windjana Sandstone, Kimberley, Gale Crater, Mars

K-Ar and noble gas surface exposure age measurements were carried out on the Windjana sandstone, Kimberley region, Gale Crater, Mars, by using the Sample Analysis at Mars instrument on the Curiosity rover. The sandstone is unusually rich in sanidine, as determined by CheMin X-ray diffraction, contributing to the high K2O concentration of 3.09 +/- 0.20 wt % measured by Alpha-Particle X-ray Spectrometer analysis. A sandstone aliquot heated to approximately 915 C yielded a K-Ar age of 627 +/- 50 Ma. Reheating this aliquot yielded no additional Ar. A second aliquot heated in the same way yielded a much higher K-Ar age of 1710 +/- 110 Ma. These data suggest incomplete Ar extraction from a rock with a K-Ar age older than 1710 Ma. Incomplete extraction at approximately 900 C is not surprising for a rock with a large fraction of K carried by Ar-retentive K-feldspar. Likely, variability in the exact temperature achieved by the sample from run to run, uncertainties in sample mass estimation, and possible mineral fractionation during transport and storage prior to analysis may contribute to these discrepant data. Cosmic ray exposure ages from He-3 and Ne-21 in the two aliquots are minimum values given the possibility of incomplete extraction. However, the general similarity between the He-3 (57 +/- 49 and 18 +/- 32 Ma, mean 30 Ma) and Ne-21 (2 +/- 32 and 83 +/- 24 Ma, mean 54 Ma) exposure ages provides no evidence for underextraction. The implied erosion rate at the Kimberley location is similar to that reported at the nearby Yellowknife Bay outcrop.

K-Ar↗

DOE OSTI2023

Effects of finite sampling on fatigue damage estimation of wind turbine components: A statistical study

The variability of the wind turbine loads complicates fatigue assessment in the design phase, as performing simulations covering the entire lifetime is computationally expensive. The current work provides important information for assessing the uncertainty in fatigue damage estimation due to finite data. Here we study the sample size effect on mean, variance, and skewness of damage in each wind bin, identify the important wind bins, and study the uncertainty propagation from each wind bin to the lifetime damage using 3600 aeroelastic simulations and bootstrapping. To achieve less than 1% error in the damage estimation across all load channels in the current case study, at least 100 turbulence seeds are needed. Damage in different wind bins follows a lognormal distribution when using the conventional approach of six seeds. The provided insights and information allow the designer to achieve a specific level of accuracy for a given computational cost using strategic bin sampling.

17 WIND ENERGY↗

DOE OSTI2021

Low-Activity Waste Glass Standards Preparation and Characterization for Calibration of Analytical Instruments

A matrix of seven low-activity waste (LAW) glasses was fabricated and characterized to be used by the Hanford Waste Treatment and Immobilization Plant (WTP) as analytical instrument calibration and matrix interference standards when measuring LAW glass compositions. The matrix of glasses was designed to represent the range of LAW glass compositions that will be generated by the WTP. Samples from each of the seven glasses were measured with electron probe microanalysis (EPMA), ion chromatography (IC), and bulk inductively coupled plasma mass spectroscopy (ICP-MS) to determine the variability and accuracy across batches and within individually poured glass bar samples. High relative percent differences and standard deviations were calculated for components with low concentrations due to instrument detection limits. IC and ICP-MS percent differences and standard deviations were higher at low component concentrations. Otherwise, variability in measured compositions from IC and ICP-MS was generally lower than measurements from EPMA, with the exclusion of Cl, F, Si, and B. Overall, the two methods proved comparable and complementary. From EPMA, IC, and ICP-MS data, volatile elements were underrepresented compared to batched compositions, which suggested loss of these elements during melting. EPMA data was statistically analyzed to evaluate homogeneity within single bars and whole compositions for each of the seven glasses. Most of the variability within glass compositions occurred at the bar-to-bar level; however, no clear systematic trends were observed. Based on these analyses, when using these glasses as reference material, it is recommended that the overall mean and variance of each composition be used for performing instrument calibrations and analytical corrections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Time variability of ultra fast BAL outflows using SALT: C IV equivalent width analysis

We study the time variability (over ≤7.3 yr) of ultra fast outflows (UFOs) detected in a sample of 64 C IV broad absorption line (BAL) quasars (with 80 distinct BAL components) monitored using the Southern African Large Telescope. By comparing the properties of the quasar in our sample with those of a control sample of non-BAL quasars, we show that the distributions of black hole mass are different and the bolometric luminosities and optical photometric variations of UFO BAL quasars are slightly smaller compared to that of non-BAL quasars. The detection fraction of C IV equivalent width (W) variability (∼95 per cent), the fractional variability amplitude $\left(\frac{\Delta W}{W}\right)$ and the fraction of ‘highly variable’ BAL (i.e. $\big|$$\frac{\Delta W}{W}$$\big|$ > 0.67) components (∼33 per cent) are higher in our sample compared to the general BAL population. The scatter in $\frac{\Delta W}{W}$ and the fraction of ‘highly variable’ BALs increase with the time-scale probed. The $\frac{\Delta W}{W}$ distribution is asymmetric at large time scales. We attribute this to the BAL strengthening time-scales being shorter than the weakening time-scales. The BAL variability amplitude correlates strongly with the BAL properties compared to the quasar properties. BALs with low W, high-velocity, shallow profiles, and low-velocity width tend to show more variability. When multiple BAL components are present, a correlated variability is seen between low- and high-velocity components with the latter showing a larger amplitude variations. We find an anticorrelation between the fractional variations in the continuum flux and W. While this suggests photoionization induced variability, the scatter in continuum flux is much smaller than that of W.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

Soil moisture variation patterns observed in Hand County, South Dakota

Soil moisture data were taken during 1976 (April, June, October), 1977 (April, May, June), and 1978 (May, June, July) Hand County, South Dakota as part of the ground truth used in NASA's aircraft experiments to study the use of microwave radiometers for the remote sensing of soil moisture. The spatial variability observed on the ground during each of the sampling events was studied. The data reported are the mean gravimetric soil moisture contained in three surface horizon depths: 0 to 2.5, 0 to 5 and 0 to 10 cm. The overall moisture levels ranged from extremely dry conditions in June 1976 to very wet in May 1978, with a relatively even distribution of values within that range. It is indicated that well drained sites have to be partitioned from imperfectly drained areas when attempting to characterize the general moisture profile throughout an area of varying soil and cover type conditions. It is also found that the variability in moisture content is greatest in the 0 to 2.5 cm measurements and decreases as the measurements are integrated over a greater depth. It is also determined that the sampling intensity of 10 measurements per km is adequate to estimate the mean moisture with an uncertainty of + or - 3 percent under average moisture conditions in areas of moderate to good drainage.

Jones, E. B.↗

DOE OSTI2018

Integrated Component-Specific Measurements to Develop Emission Factors for Compressors and Gas Gathering Lines

The Department of Energy (DOE) National Energy Technology Laboratory (NETL) and other critical stakeholders (e.g., U.S. Environmental Protection Agency (USEPA), state regulators, and industry) seek to more accurately characterize and quantify methane emissions from compressors and gathering lines within the Gathering and Boosting (G&B) segment of the natural gas system network. These efforts will reduce uncertainties in emissions estimates in USEPA’s current Green House Gas Inventory (GHGI) under their Greenhouse Gas Reporting Program (GHGRP). For this study, four G&B stations in the Gulf Coast Area were visited four times over the course of one year. Components were classified according to equipment and component type and emissions from select components were measured using multiple measurement technologies. The following major equipment types were considered: i) compressors, ii) separators, iii) dehydrators, iv) coalescers, v) slug catchers, vi) yard piping, vii) tanks, and viii) gathering lines. Individual components on each equipment type were further divided into the following main categories: i) connectors, ii) valves, iii) pressure relief valves (PRV), iv) meters, v) gauges, vi) regulators, vii) pneumatic device vents, and viii) compressor vents (e.g., distance piece vents [DPVs] and rod packing vents [RPVs]). Components were screened for the presence of leaks and venting emissions with an optical gas imaging (OGI) camera, and methane emission rates were measured using a microdilution high volume sampler. Downwind emissions were measured using an Open Path Fourier Transform Infrared (OP-FTIR) Spectrometer. Data were evaluated to assess correlation with site, equipment, time, and operational variables. Key findings and recommendations from the study include: A) Different counting methodologies can result in vastly different component counts, thus influencing total emissions calculations using Subpart W population EFs. A standard methodology should be adopted for classifying and counting components in order to minimize inconsistencies across studies and instruct stakeholders responsible for GHGI reporting. B) Results support current component subcategories listed in USEPA Subpart W, including the subdivision of connectors (flanged and “other”) and pneumatic device vents (intermittent, continuous low bleed, and continuous high bleed). C) Emission rates vary across sites and a relationship between higher emissions and compressor age was observed for some components. D) Component emission rates varied between field campaigns. The type of components contributing to overall measured emissions also varied among field campaigns. However, the emission rate variability was larger between different components (of the same category) than the variability in one component over the four field campaigns. Therefore, when calculating emission factors, it is not necessary to take repeat measurements to capture component variability. E) Measuring the emission rate of a sample for an extended period of time is not necessary to capture the average component emission rate. F) For reporting purposes, excluding liquid lines from site-wide screening may have a negligible effect on component-related emissions. However, from a site management perspective, malfunctioning components were not uncommon on liquid lines and contributed substantially to fugitive emissions at some sites. G) Population EFs, whether including or excluding liquid lines, were a conservative estimate of fugitive emissions. Leaker EFs consistently underestimated emissions at the field sites. Measurements at more sites, particularity sites in other areas of the country, need to be completed before determining the accuracy of EFs.

03 NATURAL GAS↗

DOE OSTI2016

THE MIRA–TITAN UNIVERSE: PRECISION PREDICTIONS FOR DARK ENERGY SURVEYS

Large-scale simulations of cosmic structure formation play an important role in interpreting cosmological observations at high precision. The simulations must cover a parameter range beyond the standard six cosmological parameters and need to be run at high mass and force resolution. A key simulation-based task is the generation of accurate theoretical predictions for observables using a finite number of simulation runs, via the method of emulation. Using a new sampling technique, we explore an eight-dimensional parameter space including massive neutrinos and a variable equation of state of dark energy. We construct trial emulators using two surrogate models (the linear power spectrum and an approximate halo mass function). The new sampling method allows us to build precision emulators from just 26 cosmological models and to systematically increase the emulator accuracy by adding new sets of simulations in a prescribed way. Emulator fidelity can now be continuously improved as new observational data sets become available and higher accuracy is required. Finally, using one ΛCDM cosmology as an example, we study the demands imposed on a simulation campaign to achieve the required statistics and accuracy when building emulators for investigations of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2000

Ethylene monitoring and control system

A system that can accurately monitor and control low concentrations of ethylene gas includes a test chamber configured to receive sample gas potentially containing an ethylene concentration and ozone, a detector configured to receive light produced during a reaction between the ethylene and ozone and to produce signals related thereto, and a computer connected to the detector to process the signals to determine therefrom a value of the concentration of ethylene in the sample gas. The supply for the system can include a four way valve configured to receive pressurized gas at one input and a test chamber. A piston is journaled in the test chamber with a drive end disposed in a drive chamber and a reaction end defining with walls of the test chamber a variable volume reaction chamber. The drive end of the piston is pneumatically connected to two ports of the four way valve to provide motive force to the piston. A manifold is connected to the variable volume reaction chamber, and is configured to receive sample gasses from at least one of a plurality of ports connectable to degreening rooms and to supply the sample gas to the reactive chamber for reaction with ozone. The apparatus can be used to monitor and control the ethylene concentration in multiple degreening rooms.

Nelson, Bruce N.↗

NASA NTRS2001

Ethylene monitoring and control system

A system that can accurately monitor and control low concentrations of ethylene gas includes a test chamber configured to receive sample gas potentially containing an ethylene concentration and ozone, a detector configured to receive light produced during a reaction between the ethylene and ozone and to produce signals related thereto, and a computer connected to the detector to process the signals to determine therefrom a value of the concentration of ethylene in the sample gas. The supply for the system can include a four way valve configured to receive pressurized gas at one input and a test chamber. A piston is journaled in the test chamber with a drive end disposed in a drive chamber and a reaction end defining with walls of the test chamber a variable volume reaction chamber. The drive end of the piston is pneumatically connected to two ports of the four way valve to provide motive force to the piston. A manifold is connected to the variable volume reaction chamber, and is configured to receive sample gasses from at least one of a plurality of ports connectable to degreening rooms and to supply the sample gas to the reactive chamber for reaction with ozone. The apparatus can be used to monitor and control the ethylene concentration in multiple degreening rooms.

Nelson, Bruce N.↗

DOE OSTI2019

Polynomial chaos expansions for dependent random variables

Polynomial chaos expansions (PCE) are well-suited to quantifying uncertainty in models parameterized by independent random variables. The assumption of independence leads to simple strategies for building multivariate orthonormal bases and for sampling strategies to evaluate PCE coefficients. In contrast, the application of PCE to models of dependent variables is much more challenging. Three approaches can be used to construct PCE of models of dependent variables. The first approach uses mapping methods where measure transformations, such as the Nataf and Rosenblatt transformation, can be used to map dependent random variables to independent ones; however we show that this can significantly degrade performance since the Jacobian of the map must be approximated. A second strategy is the class of dominating support methods. In these approaches a PCE is built using independent random variables whose distributional support dominates the support of the true dependent joint density; we provide evidence that this approach appears to produce approximations with suboptimal accuracy. A third approach, the novel method proposed here, uses Gram-Schmidt orthogonalization (GSO) to numerically compute orthonormal polynomials for the dependent random variables. This approach has been used successfully when solving differential equations using the intrusive stochastic Galerkin method, and in this paper we use GSO to build PCE using a non-intrusive stochastic collocation method. The stochastic collocation method treats the model as a black box and builds approximations of the input-output map from a set of samples. Building PCE from samples can introduce ill-conditioning which does not plague stochastic Galerkin methods. To mitigate this ill-conditioning we generate weighted Leja sequences, which are nested sample sets, to build accurate polynomial interpolants. We show that our proposed approach, GSO with weighted Leja sequences, produces PCE which are orders of magnitude more accurate than PCE constructed using mapping or dominating support methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

The Influence of Demographic Variables on the Pooled Rideshare Acceptance Model Multigroup Analyses (PRAMMA)

Building on our prior research with a national survey sample of 5385 US participants, the Pooled Rideshare Acceptance Model (PRAM) was built upon two factor analyses. This exploratory study extends the PRAM framework using the Pooled Rideshare Acceptance Model Multigroup Analyses (PRAMMA) to examine how 16 demographic variables influence and interact with the acceptance of Pooled Rideshare (PR), filling a gap in understanding user segmentation and personalization. Using a national sample of 5385 US participants, this methodological approach allowed for the evaluation of how PRAM variables such as safety, privacy, service experience, and environmental impact vary across diverse groups, including gender, generation, driver’s license, rideshare experience, education level, employment status, household size, number of children, income, vehicle ownership, and typical commuting practices. Factors such as convenience, comfort, and passenger safety did not show significant differences across the moderators, suggesting their universal importance across all demographics. Furthermore, geographical differences did not significantly impact the relationships within the model, suggesting consistent relationships across different regions. The findings highlight the need to move beyond a “one size fits all” approach, demonstrating that tailored strategies may be crucial for enhancing the adoption and satisfaction of PR services among various demographic groups. The analyses provide valuable insight for policymakers and rideshare companies looking to optimize their services and increase user engagement in PR.

moderator↗

DOE OSTI2018

Automated Classification and Analysis of Non-metallic Inclusion Data Sets

The aim of this study is to utilize principal component analysis (PCA), clustering methods, and correlation analysis to condense and examine large, multivariate data sets produced from automated analysis of non-metallic inclusions. Non-metallic inclusions play a major role in defining the properties of steel and their examination has been greatly aided by automated analysis in scanning electron microscopes equipped with energy dispersive X-ray spectroscopy. The methods were applied to analyze inclusions on two sets of samples: two laboratory-scale samples and four industrial samples from a near-finished 4140 alloy steel components with varying machinability. The laboratory samples had well-defined inclusions chemistries, composed of MgO-Al{sub 2}O{sub 3}-CaO, spinel (MgO-Al{sub 2}O{sub 3}), and calcium aluminate inclusions. The industrial samples contained MnS inclusions as well as (Ca,Mn)S + calcium aluminate oxide inclusions. PCA could be used to reduce inclusion chemistry variables to a 2D plot, which revealed inclusion chemistry groupings in the samples. Clustering methods were used to automatically classify inclusion chemistry measurements into groups, i.e., no user-defined rules were required.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Variation in peatland porewater chemistry over time and space along a bog to fen gradient

Porewater chemistry is an integrative measure of the physical, chemical, and biological processes occurring within peatland ecosystems, and therefore some chemistry measures (e.g., pH, calcium concentrations) have been used to classify bog vs fen peatlands. However, porewater sampling is often limited in spatial and temporal resolution, highlighting the need for a more comprehensive analysis of spatiotemporal variation in porewater chemistry. We examined depth profiles of porewater chemistry in four nearby peatlands that fall along a bog to rich fen gradient in northcentral Minnesota, USA. Porewater was sampled ~monthly during one ice-free season from three replicate piezometer nests per peatland to quantify temporal and spatial variability of those depth profiles. Porewater depth profiles of pH, calcium and total organic carbon concentrations, and δ 18 O-H 2 O varied along the bog to fen gradient, but total nitrogen and total phosphorus concentrations did not. Porewater chemistry was similar in the bogs and poor fen which were all quite different from the rich fen. In contrast, temporal and spatial variation in porewater chemistry, quantified using coefficients of variation, did not differ between bogs and fens despite the hypothesis that variation in porewater chemistry would be lower in fens than bogs due to the perennial throughput of large volumes of discharging groundwater in the rich fen. Spatial and temporal variability in porewater chemistry across all peatland types highlights the importance of collecting porewater samples in multiples over time in both near-surface and deeper peats. Furthermore, this variation can be important when scaling findings to the peatland scale, assessing the representativeness of peatlands within a larger landscape, and understanding variability in solute export to downstream ecosystems.

54 ENVIRONMENTAL SCIENCES↗