Search NASA⌕ Search

SEARCH · Search NASA

Results for “Total variance”

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

Compositional effects on cure kinetics, mechanical properties and printability of dual-cure epoxy/acrylate resins for DIW additive manufacturing

Interest in 3D printing of thermoset resins has increased significantly in recent years. One approach to additive manufacturing of thermoset resins is printing dual-cure resins with direct ink write (DIW). Dual-cure resins are multi-component resins which employ an in situ curable constituent to enable net-shape fabrication while a second constituent and cure mechanism contribute to the final mechanical properties of the printed materials. In this work, the cure kinetics, green strength, printability, and print fidelity of dual-cure epoxy/acrylate thermoset resins are investigated. Resin properties are evaluated as a function of acrylate concentration and in situ UV exposure conditions. The acrylate cure kinetics are probed using photo-differential scanning calorimetry and the impacts of resin composition and UV cure profile on the acrylate extent of conversion are presented. Continuous and pulsed UV cure profiles are shown to affect total conversion due to variances in radical efficiency at different UV intensities and acrylate concentrations. The effects of acrylate concentration on the kinetics of the epoxy thermal cure and the final mechanical properties are also investigated using dynamic mechanical analysis and three-point bend measurements. The glass transition temperature is dependent on formulation, with increasing acrylate content decreasing the T g . However, the room temperature shear moduli, flexural moduli, strength, strain-to-failure, and toughness values are relatively independent of resin composition. The similarity of the final properties allows for greater flexibility in resin formulation and in situ cure parameters, which can enable the printing of complex parts that require high green strength. We found that the in situ UV print intensities and exposure profiles that are necessary to achieve the best print quality are not, in most cases, the conditions that maximize conversion of the acrylate network. Lastly, this highlights the importance of developing optimized resin compositions which enable complete cure of the acrylate network by promoting acrylate dark cure or thermal cure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Variance Preserving Spectral Subsampling

Generating statistically faithful short-duration gamma-ray spectra from a single long measurement is essential in nuclear safeguards, supporting tasks such as algorithm development and machine-learning applications, especially when list-mode data are unavailable. Existing subsampling methods often distort the statistical characteristics of genuine short-duration measurements, leading to biased or unreliable analytical outcomes and thereby undermining downstream tasks. In this work, we compare five subsampling approaches using a benchmark set of 156 genuine replicate spectra collected with a high-purity germanium detector. We evaluate each method with respect to run-to-run variance, channel-to-channel variance, and preservation of total counts (losslessness). Across a wide range of subsampling ratios, only binomial subsampling without replacement consistently reproduces the statistical properties of genuine short-duration spectra, maintaining proper dispersion even in sparse spectral regions and perfectly preserving total counts. These results provide a mathematically principled and practically validated framework for generating synthetically shortened spectra when true short-duration measurements are unavailable.

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

Genomic selection in algae with biphasic lifecycles: A Saccharina latissima (sugar kelp) case study

Introduction Sugar kelp ( Saccharina latissima ) has a biphasic life cycle, allowing selection on both thediploid sporophytes (SPs) and haploid gametophytes (GPs). Methods We trained a genomic selection (GS) model from farm-tested SP phenotypic data and used a mixed-ploidy additive relationship matrix to predict GP breeding values. Topranked GPs were used to make crosses for further farm evaluation. The relationship matrix included 866 individuals: a) founder SPs sampled from the wild; b) progeny GPs from founders; c) Farm-tested SPs crossed from b); and d) progeny GPs from farm-tested SPs. The complete pedigree-based relationship matrix was estimated for all individuals. A subset of founder SPs ( n = 58) and GPs ( n = 276) were genotyped with Diversity Array Technology and whole genome sequencing, respectively. We evaluated GS prediction accuracy via cross validation for SPs tested on farm in 2019 and 2020 using a basic GBLUP model. We also estimated the general combining ability (GCA) and specific combining ability (SCA) variances of parental GPs. A total of 11 yield-related and morphology traits were evaluated. Results The cross validation accuracies for dry weight per meter ( r ranged from 0.16 to 0.35) and wet weight per meter ( r ranged 0.19 to 0.35) were comparable to GS accuracy for yield traits in terrestrial crops. For morphology traits, cross validation accuracy exceeded 0.18 in all scenarios except for blade thickness in the second year. Accuracy in a third validation year (2021) was 0.31 for dry weight per meter over a confirmation set of 87 individuals. Discussion Our findings indicate that progress can be made in sugar kelp breeding by using genomic selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

SDSS-IV MaNGA: Exploring the Local Scaling Relations for N/O

We present, for the first time, the relationship between local stellar mass surface density, Σ * , and N/O derived from SDSS-IV MaNGA data, using a sample of 792,765 high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modeling and partial correlation analysis, we find that Σ * alone is insufficient to predict the N/O in MaNGA spaxels and that there is an additional dependence on the local star formation rate surface density, Σ SFR . This effect is a factor of 3 stronger than the dependence of 12+log(O/H) on Σ SFR . Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed Σ * and the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when Σ * and Σ SFR are controlled for. We show that ~50% of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the buildup of oxygen in galaxies. The precise form of the N/O–O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2023

Efficient QAOA Optimization using Directed Restarts and Graph Lookup

Variational Quantum Algorithms (VQA) aim to enhance the capabilities of Noisy Intermediate-Scale Quantum (NISQ) devices. These algorithms utilize parameterized circuits and classical optimizers to iteratively execute circuits with varying parameters. However, VQA faces computational overheads due to repeated iterations and random restarts. Prior work suggests using basic sub-graphs to transfer parameters for the input graph, reducing optimizer overheads but limiting applicability to structured regular graphs. In real-world applications, random irregular graphs are common, and existing methods are not scalable or practical for such graphs. This paper presents a framework that aims to improve random irregular graphs in VQA. The framework uses graph similarity and important features like total edge counts, average edge counts, and variance. It follows an iterative process to choose basis sub-graphs from a small database and adjust parameters accordingly. Classical optimizers then utilize these parameters to determine when to restart and perform gradient descent. This approach increases the chances of reaching global maximum points.

Wang, Meng↗
DOE OSTI2021

Heterogeneous Composition of Oxygen-Evolving Complexes in Crystal Structures of Dark-Adapted Photosystem II

Photosystem II (PSII) is a homodimeric protein complex that catalyzes water oxidation at the oxygen-evolving complex (OEC), a heterocubanoid calcium-tetra-manganese cluster. Here, we analyze the omit electron-density peaks of the OEC’s metal ions in five X-ray free-electron laser PSII structures at a resolution between 2.15 Å and 1.95 Å. The omit peaks can be described by the total number of electrons and approximated by the variance of electron density distribution when the distributions are spherically symmetric. We show that the number of electrons of metal centers are different in the two OECs of PSII dimers, implying that either the oxidation states and/or occupancies of individual metal ions are different in the two monomers. In either case, we find that the two OECs of dark-adapted PSII dimers in crystals are not fully synchronized in the same expected S 1 state. Differences in redox states of the OEC in PSII only partially account for the observation that the electron densities integrate to a smaller number of electrons than expected. Differences between the determined and expected relative electron numbers are much larger than the estimated errors, indicating heterogeneity in the OEC composition. Here, the reported findings suggest that two OECs of dark-adapted PSII crystals are typically not fully synchronized in the S 1 state and likely generate heterogeneity in the crystallographical data of PSII intermediate structures.

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

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗
DOE OSTI2024

Precipitation hardening of laser powder bed fusion Ti-6Al-4V

Here, the laser powder bed fusion (PBF-L) additive manufacturing (AM) community has dedicated significant efforts into process optimization and control for defect-free Ti-6Al-4V. As defects become less of an issue for PBF-L Ti-6Al-4V, the processing-structure-properties (PSP) relationships between AM microstructures can now be explored to optimize mechanical properties. Lower temperature aging treatments around 550 °C in wrought Ti-6A-4V have been historically understood to precipitation harden the α phase with an ordered, hexagonal close packed (HCP) α 2 phase. The α 2 phase existed as nanoscale Ti 3 Al precipitates coherent with the parent α phase. The goal of the present investigation was to implement a vacuum heat treatment of 545 °C for 100 hours on PBF-L Ti-6Al-4V. This vacuum heat treatment took place after an initial hot isostatic pressure (HIP) treatment that decomposed the as-built, martensitic microstructure. The vacuum aging successfully produced nanoscale precipitates of α 2 phase within α-laths, confirmed via atom probe tomography (APT). Microstructural-length scale and quasi-static mechanical properties were investigated by nanoindentation and uniaxial tensile tests of miniaturized test specimens. Given the sensitivity of Ti-6Al-4V mechanical properties to small changes in chemistry, all test specimens originated from the same build. The α 2 precipitation resulted in significantly harder α-laths (≈ 2 GPa) as measured via nanoindentation, as well as a relative yield and ultimate tensile strength increase of 60 MPa and 38 MPa, respectively. Analysis of Variance (ANOVA) of the datasets revealed no statistical differences in total elongation between the HIPed and HIPed + aged specimens, indicative of the aging treatment producing a net benefit of strength with no loss in ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pyrolysis Molecular Beam Mass Spectrometry_Analysis_of_Natural_Variants_of_Poplulus_Trichocarpa_Leaves

Select leaves from natural variants of Poplar (Populus Trichocarpa) grown in a greenhouse at Oak Ridge National Laboratory were analyzed by Pyrolysis-Molecular Beam Mass Spectrometry (Py-MBMS). Leaves were harvested, cryomilled and kept frozen until analysis. Py-MBMS analysis was conducted using approximately 4 mg of biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗
DOE OSTI2026

Pyrolysis_Molecular_Beam_Mass_Spectrometry_Analysis_of_Specific_Switchgrass_Genotypes

Select natural variant switchgrass genotypes grown in Tifton, GA were analyzed by Pyrolysis-Molecular Beam Mass Spectrometry (Py-MBMS). Biomass was harvested, milled, several genotypes were analyzed with and without being destarched and extracted with ethanol prior to analysis (indicated with -DE if destarched and extracted). Py-MBMS analysis was conducted using approximately 4 mg of biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗
DOE OSTI2026

Pyrolysis_Molecular_Beam_Mass_Spectrometry_Analysis_of_hybrid_cross_of_Populus_tremula_x_P_alba_717-1B4_and_overexpression_of_a_lectin_receptor-like_kinase_(PtLecRLK1)

Stem tissues from the hybrid poplarPopulus tremula × P. albaclone 717-1B4 and from lectin receptor-like kinase overexpression lines PP7 and PP19 were individually colonized with the ectomycorrhizal fungiLaccaria bicolorstrain S238N,Hyaloscypha finlandicastrain PMI746, orUmbelopsis vinaceastrain PMI3018, as well as with a mixed fungal inoculum; non-inoculated plants served as controls. Plants were grown in a greenhouse at Oak Ridge National Laboratory and harvested in January 2025. Stem samples were analyzed using Pyrolysis–Molecular Beam Mass Spectrometry (Py-MBMS). Stems were harvested, debarked, dried, milled, destarched and ethanol extracted prior to analysis. Py-MBMS analysis was conducted using approximately 4 mg of wood from biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗
DOE OSTI2020

Fisher for complements: extracting cosmology and neutrino mass from the counts-in-cells PDF

In this work, we comprehensively analyse the cosmology dependence of counts-in-cells statistics. Furthermore, we focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ 8 . Additionally, using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ 8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ω m by a factor of 5 and σ 8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h -1 ) 3 . We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Temporal deposition of copper and zinc in the sediments of metal removal constructed wetlands

The objective of this study was to explore the effects of time, seasons, and total carbon (TC) on Copper (Cu) and Zinc (Zn) deposition in the surface sediments. This study was performed at the H-02 constructed wetland on the Savannah River Site (Aiken, SC, USA). Covering both warm (April-September) and cool (October-March) seasons, several sediment cores were collected twice a year from the H-02 constructed wetland cells from 2007 to 2013. Total concentrations of Cu and Zn were measured in the sediments. Concentrations of Cu and Zn (mean ± standard deviation) in the surface sediments over 7 years of operation increased from 6.0 ± 2.8 and 14.6 ± 4.5 mg kg -1 to 139.6 ± 87.7 and 279.3 ± 202.9 mg kg -1 dry weight, respectively. The linear regression model explained the behavior and the variability of Cu deposition in the sediments. On the other hand, using the generalized least squares extension with the linear regression model allowed for unequal variance and thus produced a model that explained the variance properly, and as a result, was more successful in explaining the pattern of Zn deposition. Total carbon significantly affected both Cu ( p = 0.047) and Zn ( p < 0.001). Time effect on Cu deposition was statistically significant ( p = 0.013), whereas Zn was significantly affected by the season ( p = 0.009).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Unbalanced Nested Random Effects Estimation of Variance Components

To investigate the contributing factors of variance in the measurement of an iodine 127 (127I) sample, we implement a nested random effects analysis of variance (ANOVA). Historically, the reported uncertainty on a measurement of 127I has been obtained by methods of forward uncertainty propagation because there is typically only one replicate of a given sample for which to estimate the uncertainty. When samples are processed there are several types of quality control (QC) standards analyzed along-side the unknown samples with two to five replicates for each. Assuming the variance observed in these replicate QC standards is representative of that of the unknown samples, we use these data in a nested random effects ANOVA to estimate the total uncertainty of a measurement. We demonstrate this approach with two sets of measurements from Idaho National Laboratory and compare the results with the forward uncertainty propagation approach. Variance component estimates for the coarsest level of the nesting structure were most imprecise because replicates were most limited at these levels. We find that the results largely agree between forward propagation and ANOVA, and the greatest contributors of variance are due to instrument variation and chemical processing, with human processing being among the smallest contributors. This analysis provides reassurance that the reported uncertainties using forward propagation are reasonable and the process is well controlled. We propose a future designed experiment to increase replicates at the coarsest level of the hierarchy to improve estimates of these variance components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗