Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling statistics”

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 559 records · Page 31

A Future Moon Mission: Curatorial Statistics on Regolith Fragments Applicable to Sample Collection by Raking

The strategy of raking rock fragments from the lunar regolith as a means of acquiring representative samples has wide support due to science return, spacecraft simplicity (reliability) and economy [3, 4, 5]. While there exists widespread agreement that raking or sieving the bulk regolith is good strategy, there is lively discussion about the minimum sample size. Advocates of consor-tium studies desire fragments large enough to support petrologic and isotopic studies. Fragments from 5 to 10 mm are thought adequate [4, 5]. Yet, Jolliff et al. [6] demonstrated use of 2-4 mm fragments as repre-sentative of larger rocks. Here we make use of cura-torial records and sample catalogs to give a different perspective on minimum sample size for a robotic sample collector.

Allton, J. H.↗

Regulator Interface Strategies Implemented at the Y-12 National Security Complex Old Salvage Yard Soils Remediation Project, Oak Ridge, TN - 12162

The Oak Ridge Y-12 National Security Complex housed an area known as the Old Salvage Yard (OSY) that was approximately 7 acres. The OSY was used as an area for the accumulation, processing and storage of scrap metal and equipment from Y-12 operations extending from 1968 until 2009. Areas in the northern sections of OSY also have been used for the storage of used oils containing solvents and the accumulation and recycling or de-heading and crushing of 55-gal metal drums. Scrap metal operations historically involved the accumulation, sorting, storage, public sale or disposal of scrap metal and equipment. Non-containerised storage of scrap metal was routine until 1995 when scrap metal received at OSY was placed in B-24 and B-25 boxes. Under the American Recovery and Reinvestment Act (ARRA), approximately 26,759 cubic meters of scrap metal and debris were removed and disposed at both on and off-site disposal facilities including the on-site, Oak Ridge Comprehensive Environmental Response, Compensation and Liability Act (CERCLA) landfill in 2010 and 2011. This removal action was performed in accordance with a CERCLA Record of Decision (ROD) and a close working relationship with both the U.S. Environmental Protection Agency (EPA) Region IV and Tennessee Department of Environment and Conservation (TDEC). Due to efficiencies and the excellent cooperative relationship forged with EPA Region IV and TDEC for Y-12 ARRA Cleanup Projects, a surplus of funding was available for additional remediation work that was completed in fiscal year (FY) 2011. The underlying OSY soils were targeted for characterization and potential remediation. To expedite these important activities, the U.S. Department of Energy Oak Ridge Environmental Management partnered with the regulators during detailed planning sessions through a variety of means to quickly and efficiently characterize and pinpoint areas requiring remediation according to previous ROD commitments. Data Quality Objectives (DQOs), data-sharing, real-time characterization reporting, surface and groundwater modeling and other interface planning activities were utilized to help facilitate and complete characterization and remediation activities. As a result of these strategies, the surgical extraction of one contiguous area of soil approximately 354 cubic meters is planned for FY12. The strategies discussed resulted in a major reduction of footprint remediation (i.e., 2.8% of the original estimate) which was originally estimated at over 26,759 cubic meters. The original estimate was developed using historical data collected at various times over the period of 20 years. By leveraging a hybrid sampling approach that involved both statistically-based and biased sampling locations, the area of contamination was significantly reduced resulting in both a compliant remedial design that is cost effective while mitigating a principle threat sources to surface and groundwater at the Y-12 plant. One remedial action boundary of 354 cubic meters was verified in the northern section of the Western OSY area known as the old drum de-header station for VOCs. The original estimate for disposal was in excess of 26,759 cubic meters. This area is scheduled for waste characterization and profile development in the first half of fiscal year 2012. The anticipated disposal facility is an on-site Oak Ridge CERCLA disposal landfill known as the Environmental Management Waste Management Facility (EMWMF). By utilizing the careful strategic planning, field-based screening and close cooperation of regulatory stakeholders as detailed in this paper, the total area of soil requiring remedial action within the Y-12 OSY footprint was 354 cubic meters or 2.8% of the original planned estimate. A potential waste reduction of 97.2% was realized over the original planned estimate for OSY Soils. Significant cost savings were achieved by - Minimizing the footprint of the remedial action; - Confirmatory analysis of soils instead of use of historical sampling results for waste profile development; - Targeting traditional laboratory analysis for soils that failed field screen protocols for radionuclides, heavy metals and VOCs; - Providing accurate and precise definition of RABs; and - Reduction of target contaminants for focused laboratory analysis instead of broad-range contaminants. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Composite Sampling of a Bacillus anthracis Surrogate with Cellulose Sponge Surface Samplers from a Nonporous Surface

A series of experiments was conducted to explore the utility of composite-based collection of surface samples for the detection of a Bacillus anthracis surrogate using cellulose sponge samplers on a nonporous stainless steel surface. Two composite-based collection approaches were evaluated over a surface area of 3716 cm 2 (four separate 929 cm 2 areas), larger than the 645 cm 2 prescribed by the standard Centers for Disease Control (CDC) and Prevention cellulose sponge sampling protocol for use on nonporous surfaces. The CDC method was also compared to a modified protocol where only one surface of the sponge sampler was used for each of the four areas composited. Differences in collection efficiency compared to positive controls and the potential for contaminant transfer for each protocol were assessed. The impact of the loss of wetting buffer from the sponge sampler onto additional surface areas sampled was evaluated. Statistical tests of the results using ANOVA indicate that the collection of composite samples using the modified sampling protocol is comparable to the collection of composite samples using the standard CDC protocol (p = 0.261). Most of the surface-bound spores are collected on the first sampling pass, suggesting that multiple passes with the sponge sampler over the same surface may be unnecessary. The effect of moisture loss from the sponge sampler on collection efficiency was not significant (p = 0.720) for both methods. Contaminant transfer occurs with both sampling protocols, but the magnitude of transfer is significantly greater when using the standard protocol than when the modified protocol is used (p<0.001). The results of this study suggest that composite surface sampling, by either method presented here, could successfully be used to increase the surface area sampled per sponge sampler, resulting in reduced sampling times in the field and decreased laboratory processing cost and turn-around times.

59 BASIC BIOLOGICAL SCIENCES↗

A BASIS FOR MODIFYING THE TANK 12 COMPOSITE SAMPLING DESIGN

The SRR sampling campaign to obtain residual solids material from the Savannah River Site (SRS) Tank Farm Tank 12 primary vessel resulted in obtaining appreciable material in all 6 planned source samples from the mound strata but only in 5 of the 6 planned source samples from the floor stratum. Consequently, the design of the compositing scheme presented in the Tank 12 Sampling and Analysis Plan, Pavletich (2014a), must be revised. Analytical Development of SRNL statistically evaluated the sampling uncertainty associated with using various compositing arrays and splitting one or more samples for compositing. The variance of the simple mean of composite sample concentrations is a reasonable standard to investigate the impact of the following sampling options. Composite Sample Design Option (a). Assign only 1 source sample from the floor stratum and 1 source sample from each of the mound strata to each of the composite samples. Each source sample contributes material to only 1 composite sample. Two source samples from the floor stratum would not be used. Composite Sample Design Option (b). Assign 2 source samples from the floor stratum and 1 source sample from each of the mound strata to each composite sample. This infers that one source sample from the floor must be used twice, with 2 composite samples sharing material from this particular source sample. All five source samples from the floor would be used. Composite Sample Design Option (c). Assign 3 source samples from the floor stratum and 1 source sample from each of the mound strata to each composite sample. This infers that several of the source samples from the floor stratum must be assigned to more than one composite sample. All 5 source samples from the floor would be used. Using fewer than 12 source samples will increase the sampling variability over that of the Basic Composite Sample Design, Pavletich (2013). Considering the impact to the variance of the simple mean of the composite sample concentrations, the recommendation is to construct each sample composite using four or five source samples. Although the variance using 5 source samples per composite sample (Composite Sample Design Option (c)) was slightly less than the variance using 4 source samples per composite sample (Composite Sample Design Option (b)), there is no practical difference between those variances. This does not consider that the measurement error variance, which is the same for all composite sample design options considered in this report, will further dilute any differences. Composite Sample Design Option (a) had the largest variance for the mean concentration in the three composite samples and should be avoided. These results are consistent with Pavletich (2014b) which utilizes a low elevation and a high elevation mound source sample and two floor source samples for each composite sample. Utilizing the four source samples per composite design, Pavletich (2014b) utilizes aliquots of Floor Sample 4 for two composite samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multivariate Pointwise Information-Driven Data Sampling and Visualization

With increasing computing capabilities of modern supercomputers, the size of the data generated from the scientific simulations is growing rapidly. As a result, application scientists need effective data summarization techniques that can reduce large-scale multivariate spatiotemporal data sets while preserving the important data properties so that the reduced data can answer domain-specific queries involving multiple variables with sufficient accuracy. While analyzing complex scientific events, domain experts often analyze and visualize two or more variables together to obtain a better understanding of the characteristics of the data features. Therefore, data summarization techniques are required to analyze multi-variable relationships in detail and then perform data reduction such that the important features involving multiple variables are preserved in the reduced data. To achieve this, in this work, we propose a data sub-sampling algorithm for performing statistical data summarization that leverages pointwise information theoretic measures to quantify the statistical association of data points considering multiple variables and generates a sub-sampled data that preserves the statistical association among multi-variables. Using such reduced sampled data, we show that multivariate feature query and analysis can be done effectively. The efficacy of the proposed multivariate association driven sampling algorithm is presented by applying it on several scientific data sets.

97 MATHEMATICS AND COMPUTING↗

Metabolome patterns identify active dechlorination in bioaugmentation consortium SDC-9™

Ultra-high performance liquid chromatography–high-resolution mass spectrometry (UPHLC–HRMS) is used to discover and monitor single or sets of biomarkers informing about metabolic processes of interest. The technique can detect 1000’s of molecules (i.e., metabolites) in a single instrument run and provide a measurement of the global metabolome, which could be a fingerprint of activity. Despite the power of this approach, technical challenges have hindered the effective use of metabolomics to interrogate microbial communities implicated in the removal of priority contaminants. Herein, our efforts to circumvent these challenges and apply this emerging systems biology technique to microbiomes relevant for contaminant biodegradation will be discussed. Chlorinated ethenes impact many contaminated sites, and detoxification can be achieved by organohalide-respiring bacteria, a process currently assessed by quantitative gene-centric tools (e.g., quantitative PCR). This laboratory study monitored the metabolome of the SDC-9™ bioaugmentation consortium during cis-1,2-dichloroethene (cDCE) conversion to vinyl chloride (VC) and nontoxic ethene. Untargeted metabolomics using an UHPLC-Orbitrap mass spectrometer and performed on SDC-9™ cultures at different stages of the reductive dechlorination process detected ~10,000 spectral features per sample arising from water-soluble molecules with both known and unknown structures. Multivariate statistical techniques including partial least squares-discriminate analysis (PLSDA) identified patterns of measurable spectral features (peak patterns) that correlated with dechlorination (in)activity, and ANOVA analyses identified 18 potential biomarkers for this process. Statistical clustering of samples with these 18 features identified dechlorination activity more reliably than clustering of samples based only on chlorinated ethene concentration and Dhc 16S rRNA gene abundance data, highlighting the potential value of metabolomic workflows as an innovative site assessment and bioremediation monitoring tool.

environmental monitoring↗

Calculating uncertainty on K-effective with MONK10 - 14583

Criticality safety assessments require a demonstration that a particular configuration of fissile material has an adequate sub-critical margin (k-effective sufficiently below unity) to ensure that the risk of criticality under normal operation and accident conditions is acceptable. The required sub-critical margin depends upon the uncertainty in the estimated value of k-effective. The uncertainty in the calculated value of k-effective arises from a number of sources, including: manufacturing tolerances on input data to the code (affecting geometry, compositions and densities), uncertainty in the nuclear data used by the code, stochastic uncertainty resulting from Monte Carlo simulation and modelling approximations/errors, including the inevitable bugs in the software. The ANSWERS Software Service, in collaboration with industrial partners, is developing a number of techniques to better understand and quantify uncertainty on predicted values of k-effective, using MONK. The SPRUCE utility code has been developed to allow uncertainty to be estimated using sampling methods. This can include the sampling of input parameters (including dimensions, compositions and densities) from statistical distributions. It can also include sampling different nuclear data libraries. A set of nuclear data libraries has been generated for this purpose by sampling from statistical distributions that represent the uncertainties in the published nuclear data evaluated files; a set of libraries has been produced for Latin Hypercube Sampling. By varying the input data and nuclear data, separate and combined uncertainties due to manufacturing tolerances and nuclear data can be derived. By performing least squares fitting on the results it is also possible to estimate the contribution of each of the uncertain inputs and a sensitivity method in MONK can break down the nuclear data uncertainty. (authors)

42 ENGINEERING↗

Statistical searches for microlensing events in large, non-uniformly sampled time-domain surveys: A test using palomar transient factory data

Many photometric time-domain surveys are driven by specific goals, such as searches for supernovae or transiting exoplanets, which set the cadence with which fields are re-imaged. In the case of the Palomar Transient Factory (PTF), several sub-surveys are conducted in parallel, leading to non-uniform sampling over its ∼20,000 deg{sup 2} footprint. While the median 7.26 deg{sup 2} PTF field has been imaged ∼40 times in the R band, ∼2300 deg{sup 2} have been observed >100 times. We use PTF data to study the trade off between searching for microlensing events in a survey whose footprint is much larger than that of typical microlensing searches, but with far-from-optimal time sampling. To examine the probability that microlensing events can be recovered in these data, we test statistics used on uniformly sampled data to identify variables and transients. We find that the von Neumann ratio performs best for identifying simulated microlensing events in our data. We develop a selection method using this statistic and apply it to data from fields with >10 R-band observations, 1.1 × 10{sup 9} light curves, uncovering three candidate microlensing events. We lack simultaneous, multi-color photometry to confirm these as microlensing events. However, their number is consistent with predictions for the event rate in the PTF footprint over the survey's three years of operations, as estimated from near-field microlensing models. This work can help constrain all-sky event rate predictions and tests microlensing signal recovery in large data sets, which will be useful to future time-domain surveys, such as that planned with the Large Synoptic Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

The cosmology dependence of galaxy clustering and lensing from a hybrid N -body–perturbation theory model

ABSTRACT We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from N-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of N-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these non-linear spectra from redshifts z = 0 to z = 2. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h \,{\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift z ≃ 0.4 to scales of $k_{\rm max} \approx 0.6\, h\,\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy–galaxy lensing analyses with data drawn from an independent N-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming data sets from galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

How to avoid multiple scattering in strongly scattering SANS and USANS samples

Small Angle Neutron Scattering (SANS) and Ultra Small Angle Neutron Scattering (USANS) are the only available experimental techniques to provide seamless non-destructive measurements of the geometry of the accessible and inaccessible pore structure of rocks from sub-nanopore size to the scale of macropores. They have therefore become the measurement of choice for tight reservoir rocks such as organic rich shales. A simplifying assumption in the analysis is, however, that during the path of neutrons through the sample each neutron is only scattered once. Shales are samples with a high scattering power and Multiple Scattering (MS) may occur which requires special modelling for deconvolution of the results. The approach to avoid MS is to simply reduce the sample thickness to <0.15–0.5 mm. Here, in this work, we present a systematic method on wavelength selection and preparation of samples to optimise extraction of microstructural data and minimise parasitic errors. Experimentally measured SAS transmission (TSAS) values are used as a practical criterion for estimation of the extent of MS. Generous beamtime allocations allowed robust testing revealing that sample thicknesses can be twice as thick as predicted using the standard protocol. Analysing thicker samples is particularly beneficial for statistically relevant characterisation of heterogeneous samples making the new protocol the method of choice for such samples.

(U)SANS↗

Fast-neutron Multiplicity Counter for Passive and Active Measurements of Special Nuclear Material

We developed a fast-neutron multiplicity counter (FNMC) based on organic scintillators (liquid and stilbene). The system is able to detect correlated photon and neutron multiplets emitted by fission reactions. We tested the system both in passive and active operation mode, and present the results of the experiments, both carried out at the Zero Power Physics Reactor of Idaho National Laboratory. The system was operated in coincidence mode (with a gate time of 50 nanoseconds) and two calibration-based approaches were used to estimate the 240Pu effective (240Pueff) and 235U mass of the samples. In passive mode, we assayed two types of plutonium metal plates, i.e. PAHN and PANN, with a 240Pueff mass of 4.7 g and 25.05 g per plate, respectively. The 239Pu mass percentage was 74% for the PAHN and 95% for the PANN series. We measured neutrons emitted in coincidence by several plate assemblies, having an overall 240Pueff mass in the 4.7 g to 500 g range. In active mode, we characterized several uranium samples. We used two AmLi sources (~50000 neutrons emitted per second) to induce fission in the samples. The first six measured samples were certified reference materials (CRM, set 149), having a constant enrichment, i.e. 93.2 wt%, and mass in the range between 0.5 and 4 kg. Additionally, we assayed three CRM-146 samples, with constant total mass of 500 g and an enrichment of 20.1 wt%, 52.8 wt% and 93.2 wt%. In both configurations, the system succeeded in estimating the fissile mass of an unknown sample with a statistical uncertainty lower than 1% (4-minute assay) and 7% (10-minute assay), in passive and active mode, respectively. In passive mode, the statistical uncertainty is lower than the one achieved by standard 3He counters for plutonium metal samples of similar mass. The difference between the estimated and the declared 240Pueff mass, i.e., the bias error, measured by the fast neutron multiplicity counter was consistently lower compared to the bias error of a standard thermal system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Interpretable Machine Learning for Molecular Biosignatures: a Novel Single-Sample Feature Importance Method That Is Sensitive To Statistical Interactions

Isotope ratio mass spectrometry (IRMS) of volatiles (e.g., CO 2 ) promises to be a powerful tool for potential biosignature detection for future missions to ocean worlds (OW) such as Europa and Enceladus. Machine learning (ML) methods for IRMS data could enable science autonomy by onboard prediction of seawater chemistry and biosignature presence. However, ML models are likely to be complex and involve statistical interactions between features (variables), which can make predictions seem opaque and enigmatic. For ML predictions as significant as extraterrestrial biosignatures, we must place extraordinary confidence in models. It is therefore essential that these models make interpretable predictions (i.e., human-understandable) and include false-prediction diagnostics. We achieve high accuracy and interpretability in ML biosignature and seawater chemistry models for OW through a nearest-neighbors feature selection tool that detects statistical interactions between predictors, constructs interaction networks for visualization of selected features working together to make a prediction, and reports single-sample feature importance scores for false-detection diagnostics. Here we develop a novel single-sample nearest-neighbors projected distance regression(ssNPDR) feature selection method that improves upon existing single-sample algorithms through the inclusion of statistical interactions while providing false-prediction diagnostics for ML models.

geochemistry↗

Innovative Extension of the EDF Inventory Methodology for Irradiated Graphite Sleeves - 16084

About 17,000 t of EDF irradiated graphite wastes will be produced from the decommissioning of the six French gas-cooled nuclear reactors. Determining the radionuclide content of this graphite is of capital importance for safety purposes and in order to determine the best way to manage them. A method using an identification calculation-measurement process has been proposed in order to assess a radiological inventory as precise as possible and has been validated by ANDRA, the French governmental radioactive waste management agency and by IRSN, the technical support of the French nuclear safety authority. As presented at WM2013, this method is based on (a) radionuclide measurements from sufficient samples from a statistical point of view and (b) activation coefficients calculated from the neutron flux histories of the measured samples. For that purpose, it is important to know the precise location of each measured graphite sample, which is easy with pile sampling. About 2,000 t of EDF irradiated graphite is composed of sleeves which contained fuel cartridges and they have been stored in silos of St-Laurent site. The general method needs to know the precise history including neutron flux of any sampled sleeve. As there is no possibility to individually identify sleeves, such a requirement on sleeves in a silo is impossible. The paper presents both the general method for piles and an innovative way to replace that knowledge for long lived radionuclide of irradiated sleeves. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bootstrap performance profiles in stochastic algorithms assessment

Optimization with stochastic algorithms has become a relevant research field. Due to its stochastic nature, its assessment is not straightforward and involves integrating accuracy and precision. Performance profiles for the mean do not show the trade-off between accuracy and precision, and parametric stochastic profiles require strong distributional assumptions and are limited to the mean performance for a large number of runs. In this work, bootstrap performance profiles are used to compare stochastic algorithms for different statistics. This technique allows the estimation of the sampling distribution of almost any statistic even with small samples. Multiple comparison profiles are presented for more than two algorithms. The advantages and drawbacks of each assessment methodology are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE OPTICALLY UNBIASED GAMMA-RAY BURST HOST (TOUGH) SURVEY. I. SURVEY DESIGN AND CATALOGS

Long-duration gamma-ray bursts (GRBs) are powerful tracers of star-forming galaxies. We have defined a homogeneous subsample of 69 Swift GRB-selected galaxies spanning a very wide redshift range. Special attention has been devoted to making the sample optically unbiased through simple and well-defined selection criteria based on the high-energy properties of the bursts and their positions on the sky. Thanks to our extensive follow-up observations, this sample has now achieved a comparatively high degree of redshift completeness, and thus provides a legacy sample, useful for statistical studies of GRBs and their host galaxies. In this paper, we present the survey design and summarize the results of our observing program conducted at the ESO Very Large Telescope (VLT) aimed at obtaining the most basic properties of galaxies in this sample, including a catalog of R and K{sub s} magnitudes and redshifts. We detect the host galaxies for 80% of the GRBs in the sample, although only 42% have K{sub s} -band detections, which confirms that GRB-selected host galaxies are generally blue. The sample is not uniformly blue, however, with two extremely red objects detected. Moreover, galaxies hosting GRBs with no optical/NIR afterglows, whose identification therefore relies on X-ray localizations, are significantly brighter and redder than those with an optical/NIR afterglow. This supports a scenario where GRBs occurring in more massive and dusty galaxies frequently suffer high optical obscuration. Our spectroscopic campaign has resulted in 77% now having redshift measurements, with a median redshift of 2.14 {+-} 0.18. TOUGH alone includes 17 detected z > 2 Swift GRB host galaxies suitable for individual and statistical studies-a substantial increase over previous samples. Seven hosts have detections of the Ly{alpha} emission line and we can exclude an early indication that Ly{alpha} emission is ubiquitous among GRB hosts, but confirm that Ly{alpha} is stronger in GRB-selected galaxies than in flux-limited samples of Lyman break galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Multilevel Monte Carlo Estimation of Unbiased Expectation via Sample Reuse and the Low Variance Estimation of Asymptotic Rates

A new variant of the multilevel Monte Carlo estimator [5, 3, 9, 12] is presented for the estimation of expectation statistics that utilizes sample reuse in specified levels, explicitly removes approximation error bias associated with numerically computed output quantities of interest that have an asymptotic limit behavior, and permits a low variance estimate of the asymptotic rate of convergence to that limit. In addition, it is shown that this new multilevel Monte Carlo variant can yield a computational cost savings. A review of Monte Carlo and multilevel Monte Carlo estimators is presented that includes analysis of expected value, expected mean squared error, and the calculation of optimized multilevel sample size parameters. The multilevel Monte Carlo estimator produces estimates of expectation for numerically approximated output quantities of interest that are biased by approximation error. When the quantity of interest can be modeled as the asymptotic limit of numerically approximated output quantities of interest, it is theoretically possible to remove this approximation error bias in the multilevel Monte Carlo estimator. In actual implementations, however, this procedure is unreliable due to statistical variability and inaccuracy in estimating the needed asymptotic limit. Analysis and numerical experiment show that the proposed variant of the multilevel Monte Carlo method greatly reduces (in some cases eliminates) the statistical variability in this limit estimation.

Barth, Timothy↗

ACTINIDE AND ULTRA-HEAVY ABUNDANCES IN THE LOCAL GALACTIC COSMIC RAYS: AN ANALYSIS OF THE RESULTS FROM THE LDEF ULTRA-HEAVY COSMIC-RAY EXPERIMENT

The LDEF Ultra-Heavy Cosmic-Ray Experiment (UHCRE) detected Galactic cosmic rays (GCRs) of charge Z {>=} 70 in Earth orbit with an exposure factor of 170 m{sup 2} sr yr, much larger than any other experiment. The major results include the first statistically significant uniform sample of GCR actinides with 35 events passing quality cuts, evidence for the existence of transuranic nuclei in the GCR with one {sub 96}Cm candidate event, and a low {sub 82}Pb/{sub 78}Pt ratio consistent with other experiments. The probability of the existence of a transuranic component is estimated as 96%, while the most likely {sub 92}U/{sub 90}Th ratio is found to be 0.4 within a wide 70% confidence interval ranging from 0 to 0.96. Overall, the results are consistent with a volatility-based acceleration bias and source material which is mainly ordinary interstellar medium material with some recent contamination by freshly synthesized material. Uncertainty in the key {sub 92}U/{sub 90}Th ratio is dominated by statistical errors resulting from the small sample size and any improved determination will thus require an experiment with a substantially larger exposure factor than the UHCRE.

79 ASTRONOMY AND ASTROPHYSICS↗

Does Collection Time Bias the Ecology of Cleanroom Air Samples?

Microbial monitoring of astromaterials collections has taken on increased importance with the return of biologically sensitive samples from the asteroids Ryugu and Bennu and the initiation of the Mars Sample Return Program. Terrestrial bacteria and fungi can alter the mineralogy and organic composition of our collections causing irreversible contamination of pristine samples and increasing the risk of false positives for life detection measurements. NASA has conducted routine microbial monitoring of its existing collections since 20181. Initial monitoring focused on surface samples collected with foam swabs. Although, airborne microbiology is often decoupled from surface microbiology in the built environment2 culture-based air sampling techniques like impactors were not compliant with existing contamination control requirements. Bringing organic rich media, gelatin or liquids into curation cleanrooms presents an unacceptable risk to pristine samples. In 2022 NASA purchased a materials complaint air sampler and began collecting air samples from the cleanrooms in addition to surface samples3. The new instrument uses an electret filter to collect samples that are suitable for cultivating organisms or for direct DNA sequencing. Preliminary DNA sequencing results appeared to indicate that longer sampling times biased the microbial community in favor of hearty, spore-forming bacteria3. We present the results of a study comparing overnight sampling (17 hours) to short (1 hour) sampling of unoccupied curation cleanrooms. The results will help us optimize our monitoring protocols and develop a more detailed inventory of the ecology of astromaterials curation cleanrooms. Methods: We analyzed 72 paired air samples from six different cleanrooms including the meteorite processing lab (ISO 7 equivalent, 16 samples), the lunar lab (ISO 6 equivalent, 10 samples), the stardust lab (ISO 5 equivalent 14 samples), the OSIRIS-REx lab (ISO 5 equivalent, 12 samples), the Hayabusa2 lab (ISO 5 equivalent, 14 samples), and the Genesis lab (ISO 4 equivalent, 6 samples). All the samples were collected with an InnovaPrep Bobcat air sampler operating at a sampling rate of 200 L/min. The sampler operates for 5 minutes out of every 20 minute period. Half of the samples were collected by filtering 3,000L (15 min. of active sampling) of air across an electret filter for one hour. The rest of the samples were collected by filtering approximately 51,000 L air across the filter overnight (~17 hours, 255 min. of active sampling). Cells were eluted from the filter using 6-7 ml of pressurized 0.15% tween 20 in PBS (phosphate buffered saline). This liquid was used to cultivate bacteria according to previously published methods1,4,5 and for DNA extraction and next generation sequencing. DNA was extracted with a Qiagen MagAttract PowerMicrobiome kit6. To identify bacteria and archaea, the 16S rRNA gene was amplified using Earth Microbiome primers for the V4 region 7. The amplified DNA was sequenced on an Illumina MiSeq using a V3 reagent kit. The resulting sequences were processed using DADA2 and QIIME2 as implemented on the EDGE bioinformatics platform8–10. Results: Only two of the 72 samples had no amplifiable DNA. Amplified DNA concentrations ranged from 2.67 – 0.272 ng/µl. The median concentration of amplified DNA for the 1 hour samples was 0.770 ± 0.368 ng/µl. The median concentration of amplified DNA for the overnight samples was 0.877 ± 0.434 ng/µl. On average the overnight samples had slightly more sequences (58,960 vs. 59,456) and ASV’s (amplicon sequence variants) (60 vs 64.5) than the one hour samples, but these differences are not statistically significant. The most abundant ASV in every sample mapped to the genus Cupravidus. ASV’s mapping to the genuses Bacillus, Schlegelella, Thermus, and Staphylococcus were also common. Discussion and Future Work: Alpha diversity statistics like Shannon Entropy and Faith Phylogenetic Diversity are used to describe the diversity of organisms in a single sample. If a longer sampling time was biasing the data, we would expect to see a change in these diversity statistics vs. sample time. However, we did not observe this in our data. The median Shannon entropy was slightly higher for the overnight samples (3.773 vs 3.611) as was the Faith Phylogenetic Diversity (4.042 vs 3.596), but both values were within a standard deviation of each other for the two sampling times (Fig. 1). It is unlikely, that the longer sampling time is introducing bias into our data. We do observe a significant decrease in diversity when comparing the air samples by lab. The Genesis lab (ISO 4 equivalent) has a lower median number of ASV’s (45.5) than the other labs (62). Median values for Shannon Entropy (3.717 vs. 3.430) and Faith Phylogenetic Diversity (3.796 vs. 3.548) are also lower for Genesis, but those values are with one standard deviation of each other for the different sampling times. This is consistent with previous culture-based results suggesting that the environment in cleanrooms tends to select for a core group of organisms capable of surviving under dry, low nutrient, conditions. The presence of the ASV’s mapping to Cupravidus and Thermus in our sequencing blanks and controls suggests that several of the most common organisms in our samples represent contaminants from the reagents used to perform the DNA extractions and sequencing. Further work is needed to identify these contaminants, remove them from our data and recalculate the diversity statistics. This is a systematic error. Therefore, we do not expect removing the sequencing contaminants to change our conclusions. Longer air sample collection times appear to result in slightly higher diversity and do not bias the results towards “hardy” bacteria like spore-formers. Based on these preliminary results we conclude that sampling at least 3,000 liters of air is sufficient to capture the microbial diversity of cleanrooms, and that air samples can also be collected overnight without negatively impacting diversity. These results allow us to be flexible when designing microbial monitoring plans so that they do not interfere with routine lab activity. References: 1. Regberg, A. B. et al. 49th Lunar and Planetary Science Conference (2018). 2. The United States Pharmacopeial Convention. USP General Chapter <1116> (2013). 3. Regberg, A. B., et al. 54th Lunar and Planetary Science Conference (2023). 4. Regberg, A. B. et al. 53rd Lunar and Planetary Science Conference ( 2022). 5. Davis, R. E.,et al. 50th Lunar and Planetary Science Conference (2019). 6. Qiagen. MagAttract® PowerMicrobiome® DNA/RNA EP Kit Handbook. (2018). 7. Walters, W. et al. mSystems 1, (2015). 8. Callahan, B. J. et al. Nat. Methods 13, 581–583 (2016). 9. Hall, M. & Beiko, R. G. Microbiome Analysis: Methods and Protocols113–129 (Springer, 2018). 10. Philipson, C. et al. Bio-Protoc. 7, e2622 (2017).

A. B. Regberg↗