Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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 397 records · Page 22

The Completion of a Geosynchronous Earth Orbit Survey with the Eugene Stansbery-Meter Class Autonomous Telescope

The Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT) is the primary optical sensor used by the NASA Orbital Debris Program Office (ODPO) to statistically characterize the geosynchronous Earth orbit (GEO) debris environment and support future Orbital Debris Engineering Model (ORDEM) releases. The ES-MCAT completed its first optical survey of the GEO region from 2020 to 2022. The primary goal of this survey was to autonomously collect and process GEO data with calculated photometric and astrometric uncertainties. A pointing plan was developed to provide uniform sampling within the region of interest (ROI) while accounting for predicted downtime due to insufficient observing conditions. Detections are autonomously correlated to the Space Surveillance Network (SSN) catalog to determine if objects are correlated targets (CTs) or uncorrelated targets (UCTs), the latter of which are of interest for modeling the GEO orbital debris environment. To assess the size detection sensitivity over time and monitor the general performance of the telescope’s optics and software, the optical throughput of the system and limiting magnitudes are evaluated on a routine basis. While the telescope’s ability to operate autonomously and remotely allowed for the GEO survey to continue throughout the COVID-19 pandemic, travel restrictions hampered routine cleaning of the optics during this time, and the primary mirror degraded enough to require recoating. The mirror was removed in 2022, concluding the first GEO survey. The primary mirror received a new coating designed to be more robust against the harsh environment surrounding the ES-MCAT’s location on Ascension Island, accounting for experience gained during operations over the first GEO survey. In early 2023, the recoated primary mirror was reinstalled, and the second GEO survey was initiated. The primary goal of the second GEO survey is to characterize the evolving GEO debris environment with updated optics, software, and pointing strategies while allowing for the inclusion of non-GEO regimes or those that are outside of the ROI. While the pointing method implemented in the first survey allowed for adequate coverage of the ROI over two years, it has been improved to include pointings that avoid the Moon’s position and the galactic plane to reduce software processing time and maximize the detection capabilities of fainter objects. This method also accounts for the changing weather patterns throughout the year and reduces coverage gaps in the ROI. Provided the success of the first two-year GEO survey using autonomous operations, the ODPO is actively collaborating with the United States Space Force (USSF) to make the ES-MCAT a contributing sensor to the SSN. This paper presents results from the first GEO survey including magnitude distributions and orbital parameters for CTs and UCTs. Details are provided for the automated processing pipeline and the optical system throughput for the previous and current primary mirror coatings. In addition, an updated strategy for the second GEO survey to optimize coverage over the ROI is discussed, as are preliminary results from the ongoing second survey.

Corbin Cruz↗

The Completion of a Geosynchronous Earth Orbit Survey with the Eugene Stansbery-Meter Class Autonomous Telescope

The Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT) is the primary optical sensor used by the NASA Orbital Debris Program Office (ODPO) to statistically characterize the geosynchronous Earth orbit (GEO) debris environment and support future Orbital Debris Engineering Model (ORDEM) releases. The ES-MCAT completed its first optical survey of the GEO region from 2020 to 2022. The primary goal of this survey was to autonomously collect and process GEO data with calculated photometric and astrometric uncertainties. A pointing plan was developed to provide uniform sampling within the region of interest (ROI) while accounting for predicted downtime due to insufficient observing conditions. Detections are autonomously correlated to the Space Surveillance Network (SSN) catalog to determine if objects are correlated targets (CTs) or uncorrelated targets (UCTs), the latter of which are of interest for modeling the GEO orbital debris environment. To assess the size detection sensitivity over time and monitor the general performance of the telescope’s optics and software, the optical throughput of the system and limiting magnitudes are evaluated on a routine basis. While the telescope’s ability to operate autonomously and remotely allowed for the GEO survey to continue throughout the COVID-19 pandemic, travel restrictions hampered routine cleaning of the optics during this time, and the primary mirror degraded enough to require recoating. The mirror was removed in 2022, concluding the first GEO survey. The primary mirror received a new coating designed to be more robust against the harsh environment surrounding the ES-MCAT’s location on Ascension Island, accounting for experience gained during operations over the first GEO survey. In early 2023, the recoated primary mirror was reinstalled, and the second GEO survey was initiated. The primary goal of the second GEO survey is to characterize the evolving GEO debris environment with updated optics, software, and pointing strategies while allowing for the inclusion of non-GEO regimes or those that are outside of the ROI. While the pointing method implemented in the first survey allowed for adequate coverage of the ROI over two years, it has been improved to include pointings that avoid the Moon’s position and the galactic plane to reduce software processing time and maximize the detection capabilities of fainter objects. This method also accounts for the changing weather patterns throughout the year and reduces coverage gaps in the ROI. Provided the success of the first two-year GEO survey using autonomous operations, the ODPO is actively collaborating with the United States Space Force (USSF) to make the ES-MCAT a contributing sensor to the SSN. This paper presents results from the first GEO survey including magnitude distributions and orbital parameters for CTs and UCTs. Details are provided for the automated processing pipeline and the optical system throughput for the previous and current primary mirror coatings. In addition, an updated strategy for the second GEO survey to optimize coverage over the ROI is discussed, as are preliminary results from the ongoing second survey.

Corbin Cruz↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗

Detection of Vegetation Trends in Highly Variable Environments After Grazing Exclusion in Mongolia

Aims: Environmental variability (e.g. in precipitation) has a large effect on vegetation dynamics, and this often makes it difficult to assess the recovery of vegetation after a disturbance. In this study, we assessed vegetation recovery trends in response to grazing exclusion while considering for the annual environmental variability. Location: Two regions with highly variable precipitation: a steppe near Mandalgobi, in Mongolia's Central Gobi province, and a desert steppe near Bulgan, in Mongolia's South Gobi province. Methods: Changes in vegetation were observed along grazing gradients at the above two sites, and vegetation thresholds were identified. We established reference plots in pre-threshold areas along the grazing gradients. We removed the impact of livestock grazing at various locations along the grazing gradients by establishing exclosures, and investigated vegetation from 2005 to 2013. We developed a smoothed hierarchical model within a Bayesian framework, and examined the effect of grazing exclusion on vegetation, focusing especially on the extent of grass cover recovery. In addition, we compared soil nutrient conditions in the reference plots and inside and outside each exclosure along the grazing gradients. Results: Temporal trends in the cover of perennial grass in each plot inside and outside of the exclosures largely coincided, irrespective of grazing intensity, and exclosure had no effect (Bulgan) or a negative effect (Mandalgobi) on vegetation recovery. Soil nutrient content was not significantly affected by exclosure at a given distance from the grazing source, but decreased significantly with decreasing grazing intensity. Thus, recovery of the land from a post-threshold state may not be apparent even after 9 yr of grazing exclusion in environments with highly variable precipitation. Conclusion: The effect of exclosure duration on perennial grass cover was limited, even after controlling for environmental variability. Once a vegetation threshold has been crossed, merely removing livestock from the landscape may not be sufficient for that area to recover.

Bayesian statistics;Drought;Environmental variabil↗

Artificial Intelligence Application to D and D - 20492

As aging facilities across the DOE complex await decommissioning, there is an ongoing need to understand any changes in the structural conditions. Many of these facilities were built over 50 years ago and, in some cases, these facilities have gone beyond the expected operational lifetime. Many facilities have been placed in a state of 'cold and dark,' sitting unused and awaiting decommissioning. Especially challenging are the aging facilities that provide unique operational/production capabilities to support critical DOE missions and cannot be shut down. In any of these scenarios, the structural integrity of these facilities may become compromised as time passes. It is critical that adequate inspections be performed on a continual basis and that the data collected undergoes sufficient analysis to support timely identification of any new or worsening structural issues as well as prompt needed maintenance and repairs to maintain the facilities in a safe condition. In recent days, Artificial Intelligence (AI) [1] and its application to various domains are growing at fast speed. FIU is performing research in this area and exploring the associated technologies to solve nuclear decommissioning problems. Artificial intelligence refers to the capability of a program to autonomously act, react and adapt to the working environment. AI enables the machine to behave like humans and perform the cognitive functions such as 'learning' and 'problem solving'. AI systems gradually moving from traditional approaches (algorithms and expert systems) towards more efficient and advanced technologies (machine learning [1] and deep learning [2] [3]). AI is the study of algorithms and statistical models that is being used by computers to perform specific tasks without using explicit instructions. FIU is working to develop a pilot-scale infrastructure to implement structural health monitoring using AI technologies with focus on machine learning, deep learning. This research is focused on Computer Vision/Image Classification area of AI applications. This can also be expanded to other areas of AI related to Object Recognition and Character Recognition in images. In addition to utilizing existing data sets, FIU will collect and investigate image and video data using FIU test-bed mockups to monitor structural health of the facility. Resulting data will be processed and analyzed using machine learning/deep learning technologies. The proposed pilot system is intended to serve as a starting point to engage the DOE field sites on related data sets and their decision making needs. It is anticipated that proposed machine learning/deep learning technologies can be effectively employed using anomaly detection to solve EM challenges in surveillance and maintenance of the D and D facilities. FIU will work with research stakeholders to identify applications at various sites and other DOE facilities. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of the 𝐶⁢𝑃 asymmetry in 𝐷 + → 𝜋 + ⁢𝜋 0 decays at Belle II

We measure the 𝐶⁢𝑃 asymmetry in 𝐷 + → 𝜋 + ⁢𝜋 0 decays reconstructed in 𝑒 + ⁢𝑒 − collisions at the Belle II experiment using a dataset corresponding to an integrated luminosity of 428 fb −1 . A control sample of 𝐷 + → 𝜋 + ⁢𝐾$^0_𝑆$ decays is used to correct for detection and production asymmetries. The result 𝐴 𝐶⁢𝑃 ⁡(𝐷 + → 𝜋 + ⁢𝜋 0 ) = (−1.8 ± 0.9 ± 0.1)%, where the first uncertainty is statistical and the second systematic, is the most precise determination to date. It agrees with the prediction of 𝐶⁢𝑃 symmetry from the standard model and with results of previous measurements.

CP violation↗

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗

Measurement of the Born cross section for 𝑒 +⁢ 𝑒 − → 𝜂⁢ℎ 𝑐 at center-of-mass energies between 4.1 and 4.6 GeV

We measure the Born cross section for the reaction 𝑒 + ⁢𝑒 − → 𝜂⁢ℎ 𝑐 from $\sqrt{𝑠}$ = 4.129 to 4.600 GeV using datasets collected by the BESIII detector running at the BEPCII collider. A resonant structure in the cross-section line shape near 4.200 GeV is observed with a statistical significance of 7⁢𝜎. The parameters of this resonance are measured to be 𝑀 = 4188.8 ± 4.7 ± 8.0 MeV/𝑐 2 and Γ = 49 ± 16 ± 19 MeV, where the first uncertainties are statistical and the second systematic.

exotic mesons↗

Cosmological Constraints from Combining Photometric Galaxy Surveys and Gravitational Wave Observatories

Spatial variations in survey properties due to selection effects generate substantial systematic errors in large-scale structure measurements in optical galaxy surveys on very large scales. On such scales, the statistical sensitivity of optical surveys is also limited by their finite sky coverage. By contrast, gravitational wave (GW) sources appear to be relatively free of these issues, provided the angular sensitivity of GW experiments can be accurately characterized. We quantify the expected cosmological information gain from combining the forecast LSST 3$\times$2pt analysis (combination of three 2-point correlations of galaxy density and weak lensing shear fields) with the large-scale auto-correlation of GW sources from proposed next-generation GW experiments. We find that in $\Lambda$CDM and $w$CDM models, there is no significant improvement in cosmological constraints from combining GW with LSST 3$\times$2pt over LSST alone, due to the large shot noise for the former; however, this combination does enable a $\sim6\%$ constraint on the linear galaxy bias of GW sources. More interestingly, the optical-GW data combination provides tight constraints on models with primordial non-Gaussianity (PNG), due to the predicted scale-dependent bias in PNG models on large scales. Assuming that the largest angular scales that LSST will probe are comparable to those in Stage III surveys ($\ell_{\rm min}\sim50$), the inclusion of next-generation GW measurements could improve constraints on the PNG parameter $f_{\rm NL}$ by up to a factor of $\simeq6.6$ compared to LSST alone, yielding $\sigma(f_{\rm NL})=8.5$. These results assume the expected capability of a network of Einstein Telescope-like GW observatories, with a detection rate of $10^6$ events/year. We investigate the sensitivity of our results to different assumptions about future GW detectors as well as different LSST analysis choices.

79 ASTRONOMY AND ASTROPHYSICS↗

The Cognition Battery: Developing a Normative Database for Spaceflight and Examining the Impact of Prolonged Isolation and Confinement on Cognitive Performance

INTRODUCTION: Astronauts on future long duration space exploration will be required to execute complex tasks in which even minor errors could have devastating consequences. Intact cognition is critical to maintain exceptional performance standards and it is possible that variability even at the highest ends of the performance spectrum will impact operational tasks. Exposure to spaceflight hazards could compromise cognitive performance, and decrements have been documented under conditions of altered gravity and radiation exposure. The cognitive impacts of prolonged isolation and confinement remain relatively unknown. Traditional neuropsychological assessments cannot support earth independent monitoring, do not have normative comparisons for high functioning individuals, and most are not sensitive enough to detect small performance decrements. The Cognition Battery was developed to address these limitations but currently lacks comprehensive normative data. In this project, we aimed to 1) develop a preliminary normative database for astronauts and astronaut surrogates using the Cognition Battery, and 2) characterize differences in cognitive performance after short and long duration analog missions. PARTICIPANTS & METHODS: We assessed baseline cognitive performance by administering the Cognition Battery to 97 astronaut and astronaut surrogates recruited to approximate astronaut demographics (mean age: 39.38, SD=7.62; 35.1% female; 91.7% advanced degrees). For aim 1, we calculated speed and accuracy outcomes for each subtest and summarized the data with descriptive statistics. We examined the relationship between age and performance with Pearsons’s correlations and the relationship between gender and performance with independent samples t-tests. A subset of individuals on short duration (n=48; mean age: 38.18, SD=7.14; 37.5% female; 91.7% advanced degrees) missions of 45 days and long duration (n=34; mean age: 42.12, SD=8.39; 38.2% female; 94.1% advanced degrees) missions ranging from 4-8 months were administered the Cognition Battery during and after their respective missions. For aim 2, we accounted for practice effects using published corrections, and z-transformed post-mission scores using the full sample’s baseline scores. One-way analysis of covariance tests determined the main effect of analog duration on post-mission performance, after accounting for age and gender. The False Discovery Rate method was applied to adjust for multiple comparisons. RESULTS: In the full normative sample, older age was associated with slower reaction times on a processing speed task (r=.39, p=.02), and men were more accurate on processing speed (t=-2.19, p=.03) and faster on sustained attention (t=2.02, p=.049) and risk-taking tasks (t=2.49, p=.02). We observed a main effect of duration on performance on tasks of visual memory (f(1)=5.62, p=.04), processing speed (f(1)=5.48, p=.03), and sensorimotor functioning (f(1)=39.52, p < .001), such that slower performance was observed after long duration missions relative to short. Only sensorimotor functioning was significant after adjustments for multiple comparisons (adjusted p<.001). DISCUSSION: The full sample represents the largest dataset of cognitive performance assessed by the Cognition Battery available and can be used for further research in spaceflight and high-performance populations. Relationships between performance and key demographic variables suggest future research that include age and gender stratification is needed. Our results also show minimal changes in cognitive performance between longer and shorter periods of isolation and confinement.

S I Dev↗

Process Anomaly Detection for Sparsely Labeled Events in Nuclear Power Plants

An essential aspect of online monitoring, subtle anomaly detection increases the detection lead time for equipment failure and enables a nuclear power plant (NPP) to mitigate unexpected partial or full outages, resulting in significant cost saving to the plant. Once an anomaly is detected by plant staff, its cause and severity are investigated. Because the vast majority of anomalies require some level of investigation, including some that require time-consuming examination, before they are passed over to the engineering organization for further analysis, plants are often equipped with tools to assist the staff in performing anomaly detection. Those tools operate as a black box and are often based on statistical methods that establish sensor correlations using preconfigured mathematical models and flag correlation deviations as anomalies. Due to the number of anomalies detected at a given NPP on a daily basis, a significant number of flagged anomalies usually await examination for days or weeks. A primary cause of this backlog is that the methods used by the tools generate many false positives. Though this is usually attributed to oversensitive model settings due to very narrow normal operation bands, it can also be associated with the model development being inadequate for the process being monitored, or with missing model inputs that could have explained misclassified positives. The performance of anomaly detection tools impacts their plant acceptance and utilization, especially when the effort to address false positives generated by the tool depletes the value or cost saved by using that tool. Thus, means to advance anomaly detection performance have been investigated by the Department of Energy’s Light Water Reactor Sustainability program. Previous and ongoing efforts have targeted unsupervised machine-learning (ML) methods, which do not require the labeling of any data fed into the ML model. By contrast, in supervised anomaly detection methods, every data point is labeled as either a normal or abnormal process condition, and the model is trained to replicate the classification process. Supervised methods usually outperform unsupervised methods, due to the added value in differentiating normal from anomalous states of the monitored process. An NPP’s corrective action program requires it to track and document, via a dedicated report, the resolution of any issues that occur within the plant. Once created, each report is reviewed by a plant screening committee, and several classifications and decisions are made. Recently, a collaborating NPP developed an artificial intelligence and ML-based classifier to categorize a condition report (CR) into classes that can serve to label the data as normal or anomalous. Applying CRs as labels represents a semi-supervised use case. Semi-supervised ML assumes that labels exist for some data points (i.e., labeled anomalies, in this case) but not for the rest. In this effort, semi-supervised ML methods were used to fuse data from CRs with anomaly detection methods in order to test the hypothesis that partially labeled anomalies would improve the accuracy of the anomaly detection methods. Specifically, two methods were used. The first is the deep Semi-supervised Anomaly Detection (deep SAD) method, which can handle labels ranging from fully unsupervised to fully supervised cases. The second is a newly designed ML method developed specifically for this effort and referred to as the high-order feature (HOF)-based method. To evaluate these two methods in controlled environments, synthetic data generators were developed and used. The first datasets used a spring-mass-damper (SMD) system simulator commonly found in mechanical engineering references. This was used to create two use cases: a one- and a three-mass system. Anomalies were introduced by changing the spring and damper coefficients while the system was actuated by random forces. The second datasets used the commercial Dymola-Modelica software to build a simplified nuclear reactor model. Anomalies were added in the form of corrupted sensor readings and/or control commands. The deep SAD method was tested using the SMD system, while the HOF method was tested using both datasets. Application of the deep SAD semi-supervised ML method demonstrated that labels can generate increased confidence in detecting true anomalies. This helped increase the number of true positives and decrease the number of false negatives—something that would aid in addressing the backlog of possible anomalies. Application of the HOF method demonstrated that labels can aid in down selecting from a candidate set of features to a more optimal subset in order to better differentiate between normal and anomalous conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analysis of Piezoelectric Structural Sensors with Emergent Computing Techniques

The purpose of this project was to try to interpret the results of some tests that were performed earlier this year and to demonstrate a possible use of emergence in computing to solve IVHM problems. The test data used was collected with piezoelectric sensors to detect mechanical changes in structures. This project team was included of Dr. Doug Ramers and Dr. Abdul Jallob of the Summer Faculty Fellowship Program, Arnaldo Colon-Lopez - a student intern from the University of Puerto Rico of Turabo, and John Lassister and Bob Engberg of the Structural and Dynamics Test Group. The tests were performed by Bob Engberg to compare the performance two types of piezoelectric (piezo) sensors, Pb(Zr(sub 1-1)Ti(sub x))O3, which we will label PZT, and Pb(Zn(sub 1/3)Nb(sub 2/3))O3-PbTiO, which we will label SCP. The tests were conducted under varying temperature and pressure conditions. One set of tests was done by varying water pressure inside an aluminum liner covered with carbon-fiber composite layers (a cylindrical "bottle" with domed ends) and the other by varying temperatures down to cryogenic levels on some specially prepared composite panels. This report discusses the data from the pressure study. The study of the temperature results was not completed in time for this report. The particular sensing done with these piezo sensors is accomplished by the sensor generating an controlled vibration that is transmitted into the structure to which the sensor is attached, and the same sensor then responding to the induced vibration of the structure. There is a relationship between the mechanical impedance of the structure and the resulting electrical impedance produced in the in the piezo sensor. The impedance is also a function of the excitation frequency. Changes in the real part of impendance signature relative to an original reference signature indicate a change in the coupled structure that could be the results of damage or strain. The water pressure tests were conducted by pressurizing the bottle on a test stand, and running sweeps of excitations frequencies for each of the piezo sensors and recording the resulting impedance. The sweeps were limited to 401 points by the available analyzer, and it was decided to perform individual sweeps at five different excitation frequency ranges. The frequency ranges used for the PZTs were different in two of the five ranges from the ranges used for the SCP. The bottles were pressurized to empty (no water), 0psig, 77 psig, 155 psig, 227 psig in nearly uniform increments of about 77psi. One of each of the two types of piezo sensors was fastened on to the bottle surface at two locations: about midway between the ends on cylindrical portion of the bottle and at the very edge of one of the end domes. The data was collected in files by sensor type (2 cases), by location (2 cases), by frequency range (5 cases), and pressure (5cases) to produce 100 data sets of 401 impedances. After familiarization with the piezo sensing technology and obtaining the data, the team developed a set of questions to try to answer regarding the data and made assignments of responsibilities. The next section lists the questions, and the remainder of the report describes the data analysis work performed by Dr. Ramers. This includes a discussion of the data, the approach to answering the question using statistical techniques, the use of an emergent system to investigate the data where statistical techniques were not usable, conclusions regarding the data, and recommendations.

Ramers, Douglas L.↗

Low Count and Background Radionuclides Analysis - 20488

The US Department of Energy (DOE) is often faced with the need to evaluate radionuclides at low concentrations. When site sample data are likely to be close to threshold activity concentrations of interest, then the means by which the radiochemical analysis is performed and reported is critical. This situation can occur when differentiating from zero (presence/absence) for radionuclides that do not occur naturally, close comparison with environmental background for naturally occurring radionuclides, close comparison with a risk- or dose-based threshold concentrations of interest, or even comparisons across studies. There are several analytical issues that are of concern, but the two that appear to cause incorrect decisions to be made most often involve establishing detection limits and subtracting ambient background conditions in the laboratory. These issues are not critical when radionuclide activity concentrations are large relative to thresholds of concern, but they seem to be poorly understood when it matters. When the comparisons are important and are likely to be close to a threshold of interest, then the general contract with the analytical laboratories needs to be changed so that the right or appropriate data are obtained. The concern is that important decisions are made incorrectly more often as greater scrutiny is placed on DoE's radionuclide cleanup or monitoring decisions by the public and other stakeholders. Examples are presented of problems that have been observed for different projects, both within and outside the realm of DOE and NRC remediation and radioactive waste disposal problems, and solutions are offered that should lead to better data from which important decisions need to be made. The first example is from Los Alamos National Laboratory (LANL) and involves radionuclide concentrations in soil and rock beneath LANL's Material Disposal Area (MDA) G. An initial review of the data led to a conclusion that americium and plutonium are a long way present beneath MDA G. A more thorough review of the data that accounted properly for ambient background and the detection limits that had been established led to the opposite conclusion. Another example is from the Nevada National Security Site where tritium results from one of the wells were unexpectedly high. Proper understanding and analysis of ambient background led to the conclusion that the increased concentrations were not so obvious, and that a different contract with the analytical laboratory was needed to provide more appropriate data to support a better determination. Other examples are used from regulatory review of projects in Nevada, where background levels and secular equilibrium for naturally occurring radionuclides are not established correctly because of analytical issues. The same basic issues have also been found to create difficulties analyzing historical data from the West Valley Demonstration Project. There is evidence in the data that the apparent lack of secular equilibrium where it is expected to exist is related to ambient background subtraction or other analytical issues. A final example is presented for analysis of Tc-99 in samples of depleted uranium. In this case, two different studies that were performed only three months apart provide quite different results. The US Environmental Protection Agency (EPA) established the data quality objectives (DQO) process in the mid-1980's to establish decision performance criteria for data collection. EPA guidance (EPA G-4, for example) clearly distinguishes between DQOs and measurement performance objectives (MQOs) that should be addressed for laboratory analysis of samples. The language of DQOs and MQOs has become confused over time it seems, and the subsequent effects seem to include a lack of attention to decision performance and a routine approach to measurement quality. In order to better address radionuclide sample analysis when the concentrations are close to thresholds of concern, which might be zero for some radionuclides, background for others, and risk-based thresholds for yet others, it is important that routine laboratory analysis methods are adjusted, and that the project team and the laboratory work closely together to ensure that the data meets the MQO requirements of laboratory analysis and reporting of results, and that the MQOs effectively support project-specific DQOs. This basic approach will be applied in Los Alamos in the coming year to the collection of moisture data from underneath MDA T that will be analyzed for americium and neptunium isotopes. Proper understanding of the radiochemistry methods and reporting, and of appropriate statistical methods is critical to the success of such projects, ensuring that the right decisions are made. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Partial wave analysis of 𝑒 + ⁢𝑒 − → 𝜋 + ⁢𝜋 − ⁢𝐽/𝜓 and cross section measurement of 𝑒 + ⁢𝑒 − → 𝜋 ± ⁢𝑍 𝑐 ⁢(3900) ∓ from 4.1271 to 4.3583 GeV

Based on 12.0 fb −1 of 𝑒 + ⁢𝑒 − collision data samples collected by the BESIII detector at center-of-mass energies from 4.1271 to 4.3583 GeV, a partial wave analysis is performed for the process 𝑒 + ⁢𝑒 − → 𝜋 + ⁢𝜋 − ⁢𝐽/𝜓. The cross sections for the subprocesses 𝑒 + ⁢𝑒 − → 𝜋 + ⁢𝑍 𝑐 ⁢(3900) − + c.c. → 𝜋 + ⁢𝜋 − ⁢𝐽/𝜓, 𝑓 0 ⁡(980)⁢(→ 𝜋 + ⁢𝜋 − )⁢𝐽/𝜓, and (𝜋 + ⁢𝜋 − ) S−wave⁢ 𝐽/𝜓 are measured for the first time. The mass and width of the 𝑍 𝑐 ⁢(3900) ± are determined to be 3884.6 ± 0.7 ± 3.3 MeV/𝑐 2 and 37.2 ± 1.3 ± 6.6 MeV, respectively. The first errors are statistical and the second systematic. The final state (𝜋 + ⁢𝜋 − ) S−wave ⁢𝐽/𝜓 dominates the process 𝑒 + ⁢𝑒 − → 𝜋 + ⁢𝜋 − ⁢𝐽/𝜓. By analyzing the cross sections of 𝜋 ±⁢ 𝑍 𝑐 ⁢(3900) ∓ and 𝑓 0 ⁡(980)⁢𝐽/𝜓, 𝑌⁡(4220) has been observed. Its mass and width are determined to be 4225.7 ± 4.1 ± 3.4 MeV/𝑐 2 and 57.5 ± 9.4 ± 12.1 MeV, respectively.

lepton colliders↗

Photosynthetic Efficiency of Northern Forest Ecosystems Using a MODIS-Derived Photochemical Reflectance Index (PRI)

This study evaluates a direct remote sensing approach from space for the determination of ecosystem photosynthetic light use efficiency (LUE), through measurement of vegetation reflectance changes expressed with the Photochemical Reflectance Index (PRI). The PRI is a normalized difference index based on spectral changes at a physiologically active wavelength (approximately 531 nanometers) as compared to a reference waveband, and is only available from a very few satellites. These include the two Moderate-Resolution Imaging Spectroradiometers (MODIS) on the Aqua and Terra satellites each of which have a narrow (10-nanometer) ocean band centered at 531 nanometers. We examined several PRI variations computed with candidate reference bands, since MODIS lacks the traditional 570-nanometer reference band. The PRI computed using MODIS land band 1 (620-670 nanometers) gave the best performance for daily LUE estimation. Through rigorous statistical analyses over a large image collection (n equals 420), the success of relating in situ daily tower-derived LUE to MODIS observations for northern forests was strongly influenced by satellite viewing geometry. LUE was calculated from CO2 fluxes (moles per moles of carbon absorbed quanta) measured at instrumented Canadian Carbon Program flux towers in four Canadian forests: a mature fir site in British Columbia, mature aspen and black spruce sites in Saskatchewan, and a mixed deciduous/coniferous forest site in Ontario. All aspects of the viewing geometry had significant effects on the MODIS-PRI, including the view zenith angle (VZA), the view azimuth angle, and the displacement of the view azimuth relative to the solar principal plane, in addition to illumination related variables.Nevertheless, we show that forward scatter sector views (VZA, 16 degrees-45 degrees) provided the strongest relationships to daily LUE, especially those collected in the early afternoon by Aqua (r squared = 0.83, RMSE (root mean square error) equals 0.003 moles per moles of carbon absorbed quanta). Nadir (VZA, 0 degrees plus or minus 15 degrees) and backscatter views (VZA, −16 degrees to −45 degrees) had lower performance in estimating LUE (nadir: r squared approximately equal to 0.62-0.67; backscatter: r squared approximately equal to 0.54-0.59) and similar estimation error (RMSE equals 0.004-0.005).When directional effects were not considered, only a moderately successful MODIS-PRI vs. LUE relationship (r squared equals 0.34, RMSE equals 0.007) was obtained in the full dataset (all views & sites, both satellites), but site-specific relationships were able to discriminate between coniferous and deciduous forests. Overall, MODIS-PRI values from Terra (late morning) were higher than those from Aqua (early afternoon), before/after the onset of diurnal stress responses expressed spectrally. Therefore, we identified ninety-two Terra-Aqua "same day" pairs, for which the sum of Terra morning and Aqua afternoon MODIS-PRI values (PRI (sub sum) using all available directional observations was linearly correlated with daily tower LUE (r squared equals 0.622, RMSE equals 0.013) and independent of site differences or meteorological information. Our study highlights the value of off-nadir directional reflectance observations, and the value of pairing morning and afternoon satellite observations to monitor stress responses that inhibit carbon uptake in Canadian forest ecosystems. In addition, we show that MODIS-PRI values, when derived from either: (i) forward views only, or (ii) Terra/Aqua same day (any view) combined observations, provided more accurate estimates of tower-measured daily LUE than those derived from either nadir or backscatter views or those calculated by the widely used semi-operational MODIS GPP model (MOD17) which is based on a theoretical maximum LUE and environmental data. Consequently, we demonstrate the importance of diurnal as well as off-nadir satellite observations for detecting vegetation physiological processes.

Middleton, E. M.↗

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↗

Influence of Coronal Abundance Variations

The PI of this project was Jeff Scargle of NASA/Ames. Co-I's were Alma Connors of Eureka Scientific/Wellesley, and myself. Part of the work was subcontracted to Eureka Scientific via SAO, with Vinay Kashyap as PI. This project was originally assigned grant number NCC2-1206, and was later changed to NCC2-1350 for administrative reasons. The goal of the project was to obtain, derive, and develop statistical and data analysis tools that would be of use in the analyses of high-resolution, high-sensitivity data that are becoming available with new instruments. This is envisioned as a cross-disciplinary effort with a number of "collaborators" including some at SA0 (Aneta Siemiginowska, Peter Freeman) and at the Harvard Statistics department (David van Dyk, Rostislav Protassov, Xiao-li Meng, Epaminondas Sourlas, et al). We have developed a new tool to reliably measure the metallicities of thermal plasma. It is unfeasible to obtain high-resolution grating spectra for most stars, and one must make the best possible determination based on lower-resolution, CCD-type spectra. It has been noticed that most analyses of such spectra have resulted in measured metallicities that were significantly lower than when compared with analyses of high- resolution grating data where available (see, e.g., Brickhouse et al., 2000, ApJ 530,387). Such results have led to the proposal of the existence of so-called Metal Abundance Deficient, or "MAD" stars (e.g., Drake, J.J., 1996, Cool Stars 9, ASP Conf.Ser. 109, 203). We however find that much of these analyses may be systematically underestimating the metallicities, and using a newly developed method to correctly treat the low-counts regime at the high-energy tail of the stellar spectra (van Dyk et al. 2001, ApJ 548,224), have found that the metallicities of these stars are generally comparable to their photospheric values. The results were reported at the AAS (Sourlas, Yu, van Dyk, Kashyap, and Drake, 2000, BAAS 196, v32, #54.02), and at the conference on Statistical Challenges in Modem Astronomy (Sourlas, van Dyk, Kashyap, Drake, and Pease, 2003, SCMA 111, Eds. E.D.Feigelson, G.J.Babu, New York:Springer, p489-490). We also described the limitations of one of the most egregiously misused and misapplied statistical tests in astrophysical literature, the F-test for verifying model components (Protassov, van Dyk, Connors, Kashyap, and Siemiginowska, 2002, ApJ, 571,545). Indeed, a search through the ApJ archives turned up 170 papers in the 5 previous years that used the F-test explicitly in some form or the other, and with the vast majority of them not using it correctly! Indeed, looking at just 4 issues of the ApJ in 2001, we found 13 instances of its use, of which nine were demonstrably incorrect. Clearly, it is difficult to understate the importance of this issue. We also worked on speeding up Bayes Blocks and Sparse Bayes Blocks algorithms to make them more tractable for large searches. We also supported staistics students and postdocs in both explicit physics- model-based (spectra with tens of thousands of atomic lines) and "model-free" -- i.e. non-parametric or semi-parametric -- algorithms. Work on using more of the latter is just beginning; while using multi-scale methods for Poisson imaging has come to hition. In fact, "An Image Restoration Technique with Error Estimates", by D. Esch, A. Connors, M. Karovska, and D. van Dyk, was published by ApJ (Esch et a1.2004, ApJ, 610, 1213). The code has been delivered to M. Karovska for CXC; and is available for beta-testing upon request. The other large project we worked on was on the self-consistent modeling of logN-logs curves in the Poisson limit. logN-logs curves are a fundamental tool in the study of source populations, luminosity functions, and cosmological parameters. However, their determination is hampered by statistical effects such as the Eddington bias, incompleteness due to detection efficiency, faint source flux fluctuations, etc. We have develed a new and powerful method using the full Poisson machinery that allows us to model the logN-logs distribution of X-ray sources in a self-consistent manner. Because we properly account for all the above statistical effects, our modeling is valid over the full range of the data, and not just for strong sources, as is normally done. Using a Bayesian approach and modeling the fluxes with known functional forms such as simple or broken power-laws, and conditioning the expected photon counts on the fluxes, the background contamination, effective area, detector vignetting, and detection probability, we can delve deeply into the low counts regime and extend the usefulness of medium sensitivity surveys such as ChAMP by orders of magnitude. The built-in flexibility of the algorithm also allows a simultaneous analysis of multiple datasets. We have applied this analysis to a set a Chandra observations (Sourlas, Kashyap, Zezas, van Dyk, 2004, HEAD #8, #16.32)

Scargle, Jeffrey D.↗

Millimeter-wave Point Sources from the 2500 Square Degree SPT-SZ Survey: Catalog and Population Statistics

In this work, we present a catalog of emissive point sources detected in the SPT-SZ survey, a contiguous 2530- square-degree area surveyed with the South Pole Telescope (SPT) from 2008 - 2011 in three bands centered at 95, 150, and 220GHz. The catalog contains 4845 sources measured at a significance of 4.5σ or greater in at least one band, corresponding to detections above approximately 9.8, 5.8, and 20.4 mJy in 95, 150, and 220 GHz, respectively. Spectral behavior in the SPT bands is used for source classification into two populations based on the underlying physical mechanisms of compact, emissive sources that are bright at millimeter wavelengths: synchrotron radiation from active galactic nuclei and thermal emission from dust. The latter population includes a component of high-redshift sources often referred to as submillimeter galaxies (SMGs). In the relatively bright flux ranges probed by the survey, these sources are expected to be magnified by strong gravitational lensing. The survey also contains sources consistent with protoclusters, groups of dusty galaxies at high redshift undergoing collapse. We cross-match the SPT-SZ catalog with external catalogs at radio, infrared, and X-ray wavelengths and identify available redshift information. The catalog splits into 3980 synchrotron-dominated and 865 dust-dominated sources and we determine a list of 506 SMGs. Ten sources in the catalog are identified as stars. We calculate number counts for the full catalog, and synchrotron and dusty components, using a bootstrap method and compare our measured counts with models. This paper represents the third and final catalog of point sources in the SPT-SZ survey.

79 ASTRONOMY AND ASTROPHYSICS↗