Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability study”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Decade-long time-monitoring of candidate luminous blue variable stars in the two very metal-deficient star-forming galaxies DDO 68 and PHL 293B

We have studied the spectral time variations of candidate luminous blue variable (cLBV) stars in two low-metallicity star-forming galaxies, DDO 68 and PHL 293B. The LBV in DDO 68, located in H II region #3, shows an outburst, with an increase of more than 1000 times in H α luminosity during the period 2008–2010. The broad emission of the H I and He I lines display a P Cygni profile, with a relatively constant terminal velocity of ~800 km s -1 , reaching a maximum luminosity L(H α ) of ~2 × 10 38 erg s -1 , with a full width at half-maximum (FWHM) of ~1000–1200 km s -1 . On the other hand, since the discovery of a cLBV in 2001 in PHL 293B, the fluxes of the broad components and the broad-to-narrow flux ratios of the H I and He I emission lines in this galaxy have remained nearly constant over 16 yr, with small variations. The luminosity of the broad H α component varies between ~2 × 10 38 erg s -1 and ~10 39 erg s -1 , with the FWHM varying in the range ~500–1500 km s -1 . Unusually persistent P Cygni features are clearly visible until the end of 2020 despite a decrease of the broad-to-narrow flux ratio in the most recent years. A terminal velocity of ~800 km s -1 is measured from the P Cygni profile, similar to the one in DDO 68, although the latter is 3.7 more metal-deficient than PHL 293B. In conclusion, the relative constancy of the broad H α luminosity in PHL 293B suggests that it is due to a long-lived stellar transient of type LBV/SN IIn.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Experimental test of model predictive control in a variable air volume system

Model predictive control (MPC) has been widely studied as a promising approach for improving energy efficiency and operational flexibility in buildings, yet its real-world performance for commercial variable air volume (VAV) systems remains insufficiently characterized. In particular, the impacts of model mismatch on control robustness, real-time computational burden, and device-level operation are rarely evaluated using long-term field data. Here, this study presents a comprehensive experimental evaluation of MPC applied to a full-scale VAV system in Oak Ridge National Laboratory’s Flexible Research Platform-2 building with constant cooling/heating temperature setpoints and no occupancy. The study offers three key advantages over existing work: (1) it uses a representative building in a full-scale experimental test, capturing realistic system dynamics and complexity; (2) it evaluates a relatively sophisticated MPC formulation using two different optimization solvers (Gurobi and PSO), fully accounting for computational complexity and methodological diversity; and (3) it systematically assesses potential negative impacts on various building devices, benchmark against a well-established baseline, ASHRAE Guideline 36 (G36). To isolate zone- and air-handling-unit–level supervisory control effects, the supply fan was operated with a fixed static pressure setpoint under all strategies, and the trim-and-response static pressure reset in G36 was not enabled. Results show that MPC maintained thermal comfort while improving energy efficiency. Abrupt solar radiation variations degraded performance. Computation times ranged from ∼1 s (Gurobi) to ∼ 70 s (PSO). Compared with G36, MPC achieves 33% energy savings and reduces median reheat coil output by approximately a factor of 5–10 for a representative cooling day under matched weather conditions. However, it increases the maximum discomfort deviation from 0.5 to 1°C and results in a 32% increase in staging frequency. In addition, PSO-based MPC introduced damper oscillations, also affecting actuator longevity.

ASHRAE guideline 36↗
DOE OSTI2024

Data driven discovery and quantification of hyperspectral leaf reflectance phenotypes across a maize diversity panel

Abstract Estimates of plant traits derived from hyperspectral reflectance data have the potential to efficiently substitute for traits, which are time or labor intensive to manually score. Typical workflows for estimating plant traits from hyperspectral reflectance data employ supervised classification models that can require substantial ground truth datasets for training. We explore the potential of an unsupervised approach, autoencoders, to extract meaningful traits from plant hyperspectral reflectance data using measurements of the reflectance of 2151 individual wavelengths of light from the leaves of maize ( Zea mays ) plants harvested from 1658 field plots in a replicated field trial. A subset of autoencoder‐derived variables exhibited significant repeatability, indicating that a substantial proportion of the total variance in these variables was explained by difference between maize genotypes, while other autoencoder variables appear to capture variation resulting from changes in leaf reflectance between different batches of data collection. Several of the repeatable latent variables were significantly correlated with other traits scored from the same maize field experiment, including one autoencoder‐derived latent variable (LV8) that predicted plant chlorophyll content modestly better than a supervised model trained on the same data. In at least one case, genome‐wide association study hits for variation in autoencoder‐derived variables were proximal to genes with known or plausible links to leaf phenotypes expected to alter hyperspectral reflectance. In aggregate, these results suggest that an unsupervised, autoencoder‐based approach can identify meaningful and genetically controlled variation in high‐dimensional, high‐throughput phenotyping data and link identified variables back to known plant traits of interest.

Tross, Michael C.↗
DOE OSTI2021

A search for time-dependent astrophysical neutrino emission with icecube data from 2012 to 2017

High-energy neutrinos are unique messengers of the high-energy universe, tracing the processes of cosmic ray acceleration. This work presents analyses focusing on time-dependent neutrino point-source searches. A scan of the whole sky, making no prior assumption about source candidates, is performed, looking for a space and time clustering of high-energy neutrinos in data collected by the IceCube Neutrino Observatory between 2012 and 2017. No statistically significant evidence for a time-dependent neutrino signal is found with this search during this period, as all results are consistent with the background expectation. Within this study period, the blazar 3C 279, showed strong variability, inducing a very prominent gamma-ray flare observed in 2015 June. This event motivated a dedicated study of the blazar, which consists of searching for a time-dependent neutrino signal correlated with the gamma-ray emission. No evidence for a time-dependent signal is found. Hence, an upper limit on the neutrino fluence is derived, allowing us to constrain a hadronic emission model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Advantages of a variable‐resolution global climate model in reproducing the seasonal evolution of East Asian summer monsoon

Abstract The East Asian summer monsoon (EASM) is unique among monsoon systems that it features meridional evolution of the summer monsoon. In this study, we evaluate the performances of a Variable‐Resolution Community Earth System Model (VR‐CESM) regionally refined over eastern China (14 km) in reproducing the seasonal evolution of EASM precipitation over China. Compared with reference datasets, VR‐CESM shows better performance than the corresponding globally uniform coarse‐resolution model CESM (quasi‐uniform 1°), especially over western China where complex local topography exists. The northward monsoon migration is closely related to low‐level southerly flows and vertical moisture advection, which are more reasonably simulated in VR‐CESM. The four critical timings of the EASM (monsoon onset, withdrawal, peak, and duration) are also better captured in VR‐CESM than in CESM. The corresponding spatial Pearson correlation coefficients of the four critical timings with respect to reference datasets are about 0.1 higher in VR‐CESM than those in CESM. Both models are most accurate in simulating monsoon onset and least accurate at simulating the monsoon peak. The overestimated zonal thermal contrast in CESM is responsible for the earlier monsoon onset and excessive precipitation in September over the Yangtze River valley. Finer resolution in VR‐CESM, especially over the Tibetan Plateau (TP), appears to be a main factor in simulating better zonal thermal contrast and seasonal evolution of the EASM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2025

Evapotranspiration Partitioning Using Flux Tower Data in a Semi-Arid Ecosystem

Information about evapotranspiration (ET) and its components, that is, evaporation and transpiration, is crucial for a wide range of water and ecosystem management applications. However, partitioning ET into its two components is often challenging because of their spatiotemporal variabilities and lack of process understanding. This study developed a machine learning (ML) framework to shed light on ET processes and assess the relative importance of different drivers by incorporating hydrometeorology and biomass productivity variables. The Shapley Additive Explanations (SHAP) approach was applied to enhance explainability and rank the importance of ET drivers and their components. A total of 62 variables covering hydrometeorological and biomass productivity dimensions were considered from the Reynolds Creek Critical Zone Observatory (CZO) station in Idaho. The variable importance assessment identified the leading drivers individually for evaporation, transpiration and ET (soil water content for evaporation, vapour pressure deficit for transpiration and soil water content for ET). The results further highlighted the value of combining hydrometeorological and biomass productivity variables to achieve better predictability of ET processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Strong Red Noise Ocean Forcing on Atlantic Multidecadal Variability Assessed from Surface Heat Flux: Theory and Application

The role of ocean forcing on Atlantic multidecadal variability (AMV) is assessed from the (downward) heat flux–SST relation in the framework of a new stochastic climate theory forced by red noise ocean forcing. Previous studies suggested that atmospheric forcing drives SST variability from monthly to interannual time scales, with a positive heat flux–SST correlation, while heat flux induced by ocean processes can drive SST variability at decadal and longer time scales, with a negative heat flux–SST correlation. Here, first, we develop a theory to show how the sign of heat flux–SST correlation is affected by atmospheric and oceanic forcing with time scale. In particular, a red noise ocean forcing is necessary for the sign reversal of heat flux–SST correlation. Furthermore, this sign reversal can be detected equivalently in three approaches: the low-pass correlation at lag zero, the unfiltered correlation at long (heat flux) lead, and the real part of the heat flux–SST coherence. Second, we develop a new scheme in combination with the theory to assess the magnitude and time scale of the red noise ocean forcing for AMV in the GFDL SPEAR model (Seamless System for Prediction and Earth System Research) and observations. In both the model and observations, the ocean forcing on AMV is in general comparable with the atmospheric forcing, with a 90% probability greater than the atmospheric forcing in observations. In contrast to the white noise atmospheric forcing, the ocean forcing has a persistence time comparable or longer than a year, much longer than the SST persistence of ~3 months. This slow ocean forcing is associated implicitly with slow subsurface ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ground surface temperature derived Snow Cover Properties, Seward Peninsula, Alaska, 2019-2023

Snow-ground interface temperatures have been collected at the Teller mile marker 27 and Kougarok mile marker 64 field sites on the Seward Peninsula, Alaska from 2019 through 2023 (with data missing from Fall 2020 through Summer 2021 due to COVID). Temperatures were measured using iButton Link DS1921G-F5# Thermochron miniature temperature sensors and Tinytag TGP-4017 internal sensors deployed across the Kougarok 64 and Teller 27 field sites. These sensors are a cost-efficient way to collect snow-ground interface temperatures at a high spatial resolution, and when paired with air temperature data these measurements can provide insight into fine-scale variability in snowpack characteristics across the study sites. From this data, snow process metrics were calculated at each sensor location based on the methods outlined in Staub and Delaloye, 2017. Metrics are calculated daily for each sensor as well as over the entire season. These metrics include ground surface temperature (°C), the number of days under snow cover (number of days), the insulation effect of snow (unitless), the length of the transitional snow periods (number of days), as well as intermediaries such as temperature variability. Calculating these snow processes relies on the assumption that when snow covers a temperature sensor, it is buffered from diurnal fluctuations in air temperature by the insulating snow layer. More information on the calculated metrics can be found in the User Guide of this dataset, as well as in Staub and Delaloye’s 2017 publication Using Near-Surface Ground Temperature Data to Derive Snow Insulation and Melt Indices for Mountain Permafrost Applications. This dataset includes one daily and one seasonal *.csv file of metrics for every year of data, a daily and a seasonal *.csv data dictionary, and one User Guide document (*.pdf) describing data collection and processing.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Metallothionein expression in Eastern mosquitofish ( Gambusia holbrooki ) exposed to 137 Cs and lead: implications for using metallothionein as a metal exposure biomarker

Given the range of threats facing aquatic ecosystems, it is critical we have monitoring tools to quickly and efficiently identify stressors within ecosystems. Metallothionein’s (MTs) are a group of proteins widely used as biomarkers of metal exposure; however, their suitability has been questioned due to their influenceable nature and ability to be induced by various stressors. Through our study, we sought to (1) evaluate induction of MT in eastern mosquitofish (Gambusia holbrooki) exposed to a radionuclide ( 137 Cs) and metal (Pb), (2) assess the influence of in situ 137 Cs pre-exposure on Pb-induced MT expression, and (3) quantify changes in MT concentrations during Pb depuration. We conducted our experiment in three phases using mosquitofish as a model species: Phase I—pre-exposure to 137 Cs in the field, Phase II—exposure to Pb in the lab immediately following 137 Cs exposure, and Phase III—depuration in the laboratory. We quantified MT concentrations over time and assessed differences among treatments. We did not find in situ exposure of mosquitofish to 137 Cs to induce MT, nor did pre-exposure to 137 Cs appear to influence induction of MT following Pb exposure. Metallothionein was rapidly induced in mosquitofish following Pb exposure; however, MT induction was also observed in control fish. Additionally, there was no apparent relationship between MT and Pb body burdens. During depuration, we found MT concentrations to be highly variable among fish demonstrating no clear response. Furthermore, our study demonstrates the high variability exhibited by MT, the need for further validation of MT as a metal exposure biomarker, and the value in additional research on these proteins under varying exposure scenarios, particularly environmentally realistic concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Integrating biomechanics in evolutionary studies, with examples from the amphidromous goby model system

The functional capacities of animals are a primary factor determining survival in nature. In this context, understanding the biomechanical performance of animals can provide insight into diverse aspects of their biology, ranging from ecological distributions across habitat gradients to the evolutionary diversification of lineages. To survive and reproduce in the face of environmental pressures, animals must perform a wide range of tasks, some of which entail tradeoffs between competing demands. Moreover, the demands encountered by animals can change through ontogeny as they grow, sexually mature or migrate across environmental gradients. To understand how mechanisms that underlie functional performance contribute to survival and diversification across challenging and variable habitats, we have pursued diverse studies of the comparative biomechanics of amphidromous goby fishes across functional requirements ranging from prey capture and fast-start swimming to adhesion and waterfall climbing. The pan-tropical distribution of these fishes has provided opportunities for repeated testing of evolutionary hypotheses. By synthesizing data from the lab and field, across approaches spanning high-speed kinematics, selection trials, suction pressure recordings, mechanical property testing, muscle fiber-type measurements and physical modeling of bioinspired designs, we have clarified how multiple axes of variation in biomechanical performance associate with the ecological and evolutionary diversity of these fishes. Here our studies of how these fishes meet both common and extreme functional demands add new, complementary perspectives to frameworks developed from other systems, and illustrate how integrating knowledge of the mechanical underpinnings of diverse aspects of performance can give critical insights into ecological and evolutionary questions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Single rhenium atoms on nanomagnetite: Probing the recharge process that controls the fate of rhenium in the environment

Understanding the redox transitions that control rhenium geochemistry is central to paleoredox and geochronology studies, as well as predicting the fate of chemically similar hazardous oxyanions in the environment such as pertechnetate. However, detailed mechanistic information regarding rhenium redox transitions in anoxic systems is scarce. Here, we performed a comprehensive laboratory study of rhenium redox transitions on variably oxidized magnetite nanoparticle surfaces. Through high-end spectroscopic and microscopic tools, we propose an abiotic transition pathway in which aqueous iron(II) ions in the presence of pure or preoxidized magnetite serve as an electron source to reduce rhenium(VII) to individual rhenium(IV) atoms or small polynuclear species on nanoparticle surfaces. Notably, iron(II) ions recharged preoxidized magnetite nanoparticles exhibit a maghemite core and a magnetite shell, challenging the traditional core-shell magnetite-maghemite model. This study provides a fundamental understanding of redox processes governing rhenium fate and transport in the environment and enables an improved basis for predicting its speciation in geochemical systems.

Science & Technology - Other Topics↗
DOE OSTI2025

SMBH binary candidate PKS J2134−0153: possible multi-band periodic variability and inter-band time lags

ABSTRACT Studying the periodic flux-variation behaviour of blazars is vital for probing supermassive black hole binaries and the kinematics of relativistic jets. In this work, we report the detection of multi-band possible periodic variations of the blazar PKS J2134−0153, including infrared ($1.6(\pm 0.4)\times 10^3$ day) and optical ($1.8(\pm 1)\times 10^3$ day) wavelengths. The periods in the infrared and optical bands are statistically consistent with the period in the radio band ($P_{\mathrm{Radio}} = 1760\pm 33$ days, obtained from our previous work). Moreover, flux variations in different bands are correlated with evident inter-band time delays, and the time lags of infrared and optical emission with respect to radio emission are $(3.3\pm 2.3)\times 10^{2}$ days and $(3.0\pm 2.3)\times 10^{2}$ days, respectively. The cross-correlations indicate a common origin of radio, infrared, and optical emission. The relative positions between emission regions of infrared and optical emission and radio emission are estimated according to the time lags, i.e. $0.37\pm 0.26$ and $0.33\pm 0.26$ pc. The relative distances seem to be quantitatively consistent with the theoretical prediction.

Ren, Guo-Wei (ORCID:0000000214978371)↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗
DOE OSTI2022

Variable wildfire impacts on the seasonal water temperatures of western US streams: A retrospective study

Recent increases in the burn area and severity of wildfires in the western US have raised concerns about the impact on stream water temperature–a key determinant of cold-water fish habitats. However, the effect on seasonal water temperatures of concern, including winter and summer, are not fully understood. In this study, we assessed the impact of wildfire burns at Boulder Creek (Oregon), Elk Creek (Oregon), and Gibbon River (Wyoming) watersheds on the downstream winter and summer water temperatures for the first three post-fire years. To obtain results independent of the choice of the analytical method, we evaluated the consequence of each burn using three different statistical approaches that utilize local water temperature data. Our results from the three approaches indicated that the response of water temperatures to wildfire burns varied across seasons and sites. Wildfire burns were associated with a median increase of up to 0.56°C (Standard Error; S.E. < 0.23°C) in the summer mean water temperatures (MWT) and 62 degree-day Celsius (DDC; S.E. < 20.7 DDC) in the summer accumulated degree days (ADD) for the three subsequent years across studied stream sites. Interestingly, these burns also corresponded to a median decrease of up to 0.49°C (S.E. < 0.45°C) in the winter MWT and 39 DDC (S.E. < 40.5 DDC) in the winter ADD for the same period across sites. Wildfire effects on the downstream water temperatures diminished with increasing site distance from the burn perimeter. Our analyses demonstrated that analytical methods that utilize local watershed data could be applied to evaluate fire effects on downstream water temperatures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uncertainty-aware mixed-variable machine learning for materials design

Abstract Data-driven design shows the promise of accelerating materials discovery but is challenging due to the prohibitive cost of searching the vast design space of chemistry, structure, and synthesis methods. Bayesian optimization (BO) employs uncertainty-aware machine learning models to select promising designs to evaluate, hence reducing the cost. However, BO with mixed numerical and categorical variables, which is of particular interest in materials design, has not been well studied. In this work, we survey frequentist and Bayesian approaches to uncertainty quantification of machine learning with mixed variables. We then conduct a systematic comparative study of their performances in BO using a popular representative model from each group, the random forest-based Lolo model (frequentist) and the latent variable Gaussian process model (Bayesian). We examine the efficacy of the two models in the optimization of mathematical functions, as well as properties of structural and functional materials, where we observe performance differences as related to problem dimensionality and complexity. By investigating the machine learning models’ predictive and uncertainty estimation capabilities, we provide interpretations of the observed performance differences. Our results provide practical guidance on choosing between frequentist and Bayesian uncertainty-aware machine learning models for mixed-variable BO in materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗