Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collective variables”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Deep X-Ray Observation of NGC 3221: Everything Everywhere All at Once

We present a comprehensive analysis of 475 ks (438 ks unpublished and 37 ks archival) XMM-Newton/EPIC-pn observations of a nearby, highly inclined, star-forming, luminous infrared galaxy NGC 3221 through spatial, temporal, and spectral information. We confirm the presence of a low-luminosity (presumably Compton-thick) active galactic nucleus (AGN). The 0.4–12 keV luminosity and the hardness ratio of the six ultraluminous X-ray sources previously identified in Chandra data exhibit diverse variability on day scales. The collective emission from unresolved sources exhibits a different day-scale variability. We have also discovered two new predominantly soft (<1 keV) sources. One of these has an enigmatic spectral shape featuring a soft component, which we interpret as a superbubble in NGC 3221, and a variable hard component from a compact object, unresolved from the superbubble. We do not confidently detect any X-ray emission from SN 1961L. The hot gas in the interstellar medium (ISM, out to ±6 kpc from the disk plane) and that in the extraplanar region (6–12 kpc) both require two thermal phases at ∼0.15 keV and ∼0.55 keV. The ∼0.55 keV component is fainter in the ISM than the ∼0.15 keV component, but the emission from the latter falls off more steeply with disk height than the emission from the former. This makes the extraplanar region hotter and less dense than the ISM. The proximity of NGC 3221 and the occurrence of the underluminous AGN offer a unique observing opportunity to study the hot diffuse medium along with nuclear and diskwide point sources.

Das, Sanskriti [Stanford Univ., CA (United States)↗
DOE OSTI2024

Timeseries Unlabeled and Labeled Photos, Modeled Stream Elevation, and (Meta)Data of Variably Inundated Streams Across The Yakima River Basin, Washington, United States (v2)

This dataset is associated with the “River Monitoring Photos” (RMP) study and subsequent manuscript (Bao et al. 2025. Monitoring river flow status using low-cost wildlife camera and image segmentation artificial intelligence doi: 10.1016/j.envsoft.2025.106715). Game camera timeseries photos were collected to evaluate stream variable inundation via changes in width. A subset of photos was labeled for training the YOLOv8 and Mask2Former models and used to segment water surface fractions from all the game camera photos.This data package was originally published in March 2024. It was updated in October 2025 (v2) to add additional photos and files associated with the manuscript (i.e., processed data, labeled photos, and trained models). For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.In addition to a readme, this data package also includes two file-level metadata (FLMD) files that describes each file and two data dictionaries (DD) that describe all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) readme; (4) field metadata; and (5) folders containing game camera photos and manuscript-associated files. Each Yakima River Basin site has a folder that contains subfolders for each month photos were collected. There is also a folder for files associated with the manuscript which has subfolders for labeled data, trained models, Yakima River Basin site water surface fractions, and USGS site water surface fractions. All files are .csv, .json, .txt, .yaml, .pth, .pt, or .pdf. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Learning Planar Ising Models Software

Learning Planar Ising Models is a software package written in Matlab for learning relationships among variable in a dataset using graphical models. The software package implements a generally-applicable algorithm for learning planar Ising models from any multivariate dataset. The code provides an algorithm for learning the best planar Ising model to approximate an arbitrary collection of binary random variables (possibly from sample data). Given the set of all pairwise correlations among variables, we select a planar graph and optimal planar Ising model defined on this graph to best approximate that set of correlations. The software includes demonstrations of the algorithm in simulations and for applications on publicly available datasets. Details of the algorithm, demonstration simulations, and applications are given in Johnson, et al; 2016. Reference: Johnson, J. K., Oyen, D., Chertkov, M., and Netrapalli, P. (2016). Learning planar Ising models. Journal of Machine Learning Research.

Oyen, Diane↗
DOE OSTI2022

A Framework for Closed-Loop Optimization of an Automated Mechanical Serial-Sectioning System via Run-to-Run Control as Applied to a Robo-Met.3D

Optimization of automated data collection is gaining increased interest for the purposes of enabling closed-loop self-correcting systems that inherently maximize operational efficiencies and reduce waste. Many data collection systems have several variables which influence data accuracy or consistency and which can require frequent user interaction to be monitored and maintained. Operating upon a Robo-MET.3D™ automated mechanical serial-sectioning system, a run-to-run control algorithm has been developed to accelerate data collection and reduce data inconsistency. Here, using historical data amassed over a decade of experiments, a linear regression model of the deterministic system dynamics is created and used to employ a run-to-run control algorithm that optimizes selected system inputs to reduce operator intervention and increase efficacy while reducing variance of system output.

42 ENGINEERING↗
DOE OSTI2022

Hourly water-carbon interactions modulate decadal water-use efficiency trends inferred from ecosystem-scale measurements

Plant stomatal conductance regulates photosynthesis and transpiration. This physiological link affects ecosystem responses to microclimate and harmonizes carbon, energy, and water exchanges between the biosphere and atmosphere. The relationship between water losses via transpiration and carbon gains via photosynthesis can be quantified by plant water-use efficiency (WUE). While leaf- and ecosystem-scale observations both suggest rising WUE in recent decades, WUE trends inferred from the ecosystem scale are much larger than those inferred from the leaf scale or implied by theory. The unexpectedly large ecosystem-scale WUE trends complicate interpretation of ecophysiological responses to changing environmental conditions. Here, we analyze ecosystem-scale WUE inferred from 40 FLUXNET sites, each with at least 10 years of measurements. Furthermore, our results demonstrate that observed ecosystem-scale WUE trends are more sensitive to hourly weather conditions than longer-term changes in atmospheric carbon dioxide or vapor pressure deficit. Our analysis shows that Earth System Models participating in CMIP6 did not capture the observed WUE sensitivity to inter-site variability and microclimatic conditions. Collectively, our findings suggest that ecosystem-scale WUE trends reflect water-carbon interactions across multiple temporal scales, and disentangling factors contributing to emergent ecosystem responses is needed to infer ecophysiological relationships and model structures from observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Coercive Fields Exceeding 30 T in the Mixed-Valence Single-Molecule Magnet (Cp iPr5 ) 2 Ho 2 I 3

Mixed-valence dilanthanide complexes of the type (Cp iPr5 ) 2 Ln 2 I 3 (Cp iPr5 = pentaisopropylcyclopentadienyl; Ln = Gd, Tb, Dy) featuring a direct Ln–Ln σ-bonding interaction have been shown to exhibit well-isolated high-spin ground states and, in the case of the Tb and Dy variants, a strong axial magnetic anisotropy that gives rise to a large magnetic coercivity. Here, we report the synthesis and characterization of two new mixed-valence dilanthanide compounds in this series, (Cp iPr5 ) 2 Ln 2 I 3 (1-Ln; Ln = Ho, Er). Both compounds feature a Ln–Ln bonding interaction, the first such interaction in any molecular compounds of Ho or Er. Like the Tb and Dy congeners, both complexes exhibit high-spin ground states arising from strong spin–spin coupling between the lanthanide 4f electrons and a single σ-type lanthanide–lanthanide bonding electron. Beyond these similarities, however, the magnetic properties of the two compounds diverge. In particular, 1-Er does not exhibit observable magnetic blocking or slow magnetic relaxation, while 1-Ho exhibits magnetic blocking below 28 K, which is the highest temperature among Ho-based single-molecule magnets, and a spin reversal barrier of 556(4) cm –1 . Additionally, variable-field magnetization data collected for 1-Ho reveal a coercive field of greater than 32 T below 8 K, more than 6-fold higher than observed for the bulk magnets SmCo 5 and Nd 2 Fe 14 B, and the highest coercive field reported to date for any single-molecule magnet or molecule-based magnetic material. Multiconfigurational calculations, supported by far-infrared magnetospectroscopy data, reveal that the stark differences in magnetic properties of 1-Ho and 1-Er arise from differences in the local magnetic anisotropy of the lanthanide centers.

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

BCFW tilings and cluster adjacency for the amplituhedron

In 2005, Britto, Cachazo, Feng, and Witten gave a recurrence (now known as the BCFW recurrence) for computing scattering amplitudes in N = 4 super Yang–Mills theory. Arkani-Hamed and Trnka subsequently introduced the amplituhedron to give a geometric interpretation of the BCFW recurrence. Arkani-Hamed and Trnka conjectured that each way of iterating the BCFW recurrence gives a “triangulation” or “tiling” of the m=4 amplituhedron. In this article, we prove the BCFW tiling conjecture of Arkani-Hamed and Trnka. We also prove the cluster adjacency conjecture for BCFW tiles of the amplituhedron, which says that facets of tiles are cut out by collections of compatible cluster variables for the Grassmannian Gr4,n. Moreover we show that each BCFW tile is the subset of the Grassmannian where certain cluster variables have particular signs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Aboveground Rather Than Belowground Productivity Drives Variability in Miscanthus × giganteus Net Primary Productivity

Quantifying the carbon (C) uptake of Miscanthus × giganteus ( M × g ) in both aboveground and belowground structures (e.g., net primary productivity (NPP)) and differences among methodological approaches is crucial. Our objectives were to directly measure Mxg NPP and evaluate the effects of nitrogen application, location, and belowground biomass sampling methods. We hypothesize that increased nitrogen application increases the overall NPP of M × g and that quantifying rhizome biomass using excavations will produce the lowest variability between replicates. We collected biomass from mature M × g stands from three locations in Iowa with three nitrogen application rates and one site in Illinois. We destructively sampled at two time points, when rhizome mass is anticipated to be at a minimum (initial) and anticipated to be at its maximum (peak). Biomass was collected from 1 × 1 m quadrats in which one in-clump and one beside-clump cores were collected and then excavated to 30 cm depth to extract all rhizomes. We found that aboveground M × g NPP ranged from 15.4 Mg DM ha –1 year –1 to 36.4 Mg DM ha –1 year–1 and belowground M × g NPP ranged from 4.4 Mg DM ha –1 year –1 to 19.6 Mg DM ha –1 year –1 . M × g NPP varied across sites, fertilization, and calculation assumptions. Aboveground NPP (yield) was on average 68.7% of the total NPP. Root-to-shoot ratios at peak biomass decreased with nitrogen application rate, from an average of 1.9 for 0 N plots to 0.89 for 224 N fertilized plots. There was more variation in core data than from excavations; however, when in-clump and beside-clump cores were averaged together, core and excavation averages were not different. Overall, these results show that the range of mature M × g NPP is driven by aboveground productivity, influenced by nitrogen application and site. Our results provide useful data to constrain agro-ecosystem models and provide crucial insights for future perennial belowground sampling.

Hartman, Theodore [Univ. of Illinois at Urbana-Cha↗
DOE OSTI2026

Synoptic NO3 in Slate River Watershed, Colorado (2022)

The synoptic nitrate (NO3) dataset in the Slate River Watershed, Colorado consists of NO3 data collected at 19 locations three times during the summer of 2022. Stream samples were collected in early summer (early July), mid summer (late August), and late summer (late September). The samples include mainstem, tributary, and point source input water samples. These data were collected to evaluate spatiotemporal variability in stream NO3 during the summer, and evaluate anthropogenic controls on stream NO3 dynamics. This data package contains: (1) a csv of all NO3 samples and (2) a csv of locations for each sampling site. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > BIOSPHERE > ECOSYSTEMS > FRESHWATE↗
DOE OSTI2025

Assessing Heterogeneity of Surface Water Temperature Following Stream Restoration and a High-Intensity Fire from Thermal Imagery

Thermal heterogeneity of rivers is essential to support freshwater biodiversity. Salmon behaviorally thermoregulate by moving from patches of warm water to cold water. When implementing river restoration projects, it is essential to monitor changes in temperature and thermal heterogeneity through time to assess the impacts to a river’s thermal regime. Lightweight sensors that record both thermal infrared (TIR) and multispectral data carried via unoccupied aircraft systems (UASs) present an opportunity to monitor temperature variations at high spatial (<0.5 m) and temporal resolution, facilitating the detection of the small patches of varying temperatures salmon require. Here, we present methods to classify and filter visible wetted area, including a novel procedure to measure canopy cover, and extract and correct radiant surface water temperature to evaluate changes in the variability of stream temperature pre- and post-restoration followed by a high-intensity fire in a section of the river corridor of the South Fork McKenzie River, Oregon. We used a simple linear model to correct the TIR data by imaging a water bath where the temperature increased from 9.5 to 33.4 °C. The resulting model reduced the mean absolute error from 1.62 to 0.35 °C. We applied this correction to TIR-measured temperatures of wetted cells classified using NDWI imagery acquired in the field. We found warmer conditions (+2.6 °C) after restoration (p < 0.001) and median absolute deviation for pre-restoration (0.30) to be less than both that of post-restoration (0.85) and post-fire (0.79) orthomosaics. In addition, there was statistically significant evidence to support the hypothesis of shifts in temperature distributions pre- and post-restoration (KS test 2009 vs. 2019, p < 0.001, D = 0.99; KS test 2019 vs. 2021, p < 0.001, D = 0.10). Moreover, we used a Generalized Additive Model (GAM) that included spatial and environmental predictors (i.e., canopy cover calculated from multispectral NDVI and photogrammetrically derived digital elevation model) to model TIR temperature from a transect along the main river channel. This model explained 89% of the deviance, and the predictor variables showed statistical significance. Collectively, our study underscored the potential of a multispectral/TIR sensor to assess thermal heterogeneity in large and complex river systems.

Barker, Matthew I. (ORCID:0000000252864930)↗
DOE OSTI2020

Adaptive Neuro Fuzzy Inference System for Cyber-Intrusion Detection in a Smart Grid

The evolution of the power grid has brought increasing deployment of advance metering infrastructure, penetration of intelligent electronic devices, and integration of physical power system components with information and communications technologies. With the fast-expanding connectivity, cyber vulnerabilities arise due to the use of internet-based communication systems. These systems are targets of cyber-intrusions which attempt to disturb the normal power system functions. Traditional intrusion detection algorithms have been developed without an explicit model of the cyber components. In this paper, an algorithm to detect false data injections in the power system is proposed considering both cyber and physical models of the power system. The algorithm is based on an Adaptive Neuro Fuzzy Inference System (ANFIS) which collects information from state variables of the cyber-physical system to meet the performance requirements of the grid. Simulations of the proposed approach using the IEEE 13-bus test system validate the effectiveness of this artificial intelligence-based algorithm.

Adaptive Neuro Fuzzy Inference System (ANFIS), cyb↗
DOE OSTI2024

Identifying Controlling Variables for Mercury Vapors in Alpha-4 at Y-12: Two Year Data Collection Update

Multiple sensor packages were deployed at Alpha-4 by SRNL, in collaboration with United Cleanup Oak Ridge LLC (UCOR), to monitor mercury vapor concentrations and meteorological parameters. These sensors collected data, inside and outside of the legacy-use facility, for approximately two years. Though data gaps still exist, particularly in colder months, several controlling variables were identified that govern mercury vapor concentrations within Alpha-4. Temperature, barometric pressure gradients, humidity, and wind speed have been identified as controlling variables. A strong positive correlation was seen between mercury vapor concentrations and temperature which generally followed diurnal fluctuations. Temperatures below approximately 10 degrees Celsius did not show any spikes above the PEL, indicating more work can be performed at any time during the winter months – more data should be collected to confirm consistency in this finding. Additionally, spikes in mercury vapor concentrations above that of the permissible exposure limit (PEL; 100 µg/m 3 ) occurred primarily in late afternoon or evening/overnight hours (between 3 PM and 6 AM), which suggests D&D operations might be best scheduled during morning or daytime hours prior to the late afternoon. However, a limited number of spikes did occur outside of the identified window, although this may be attributed to disturbances in air flow and mercury vapor release from work activities performed inside of the Alpha-4 building. The analysis conducted allows for a strong predictive capability for estimating mercury vapor concentrations based upon accurate meteorological parameters. Still, additional data collection, particularly in the winter months, could help to strengthen the predictive power and validate the identified data trends. Further, increased temporal resolution could also help to better characterize the incipient stages of the increases and decreases in the mercury vapor concentration. Continued monitoring support by SRNL at Y-12 is underway at Alpha-4 to further close remaining data gaps and support deactivation and decommissioning work. Within a collaborative effort with UCOR, the SRNL team is collecting mercury vapor data to study the efficacy of a novel mercury suppressant, FerroBlack® which was recently deployed at Alpha-4. In addition, modifications to the current monitoring setup to increase measurement resolution is also being investigated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

NGEE Arctic, Canopy Spectral Reflectance, Utqiagvik (Barrow), Alaska, 2014-2016

Measurements of full-spectrum (i.e. 350-2500 nm) canopy spectral reflectance of Arctic plant species within the BEO, Barrow, Alaska. Spectra were collected using Spectra Vista Corporation (SVC) HR-2014i and Spectral Evolution (SE) PSR+ instruments mounted on a tripod or monopod together with a Spectralon white plate to calibrate each measurement under variable illumination conditions. Data were collected in Utqiagvik (Barrow), Alaska during the 2014 to 2016 period. Data package includes processed spectra (*.csv) files and original instrument output and processing files (*.zip).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) 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 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 OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Interplay of Crystal Structure and Magnetic Properties of the Eu 5.08-x Sr x Al 3 Sb 6 Solid Solution

Zintl phases containing rare-earth metals have gained attention due to their magnetic, electronic, and thermoelectric properties. Eu 5.08 Al 3 Sb 6 is a new structure type (monoclinic space group C2/m) that can be described as a pseudorock-salt EuSb motif with the Eu-centered Sb octahedra at the origin of the unit cell, and on the C-face center, containing either Eu (8%) or an Al 4 tetrahedron modeled as a dual tetrahedron (37.5%). The complete solid solution of Eu 5.08–x Sr x Al 3 Sb 6 can be prepared; however, the cation totals vacillate from 5 to 5.24 depending on the Al content. Al K-edge XANES shows a shift to higher energy relative to the Al metal but at slightly lower energy relative to AlSb, indicating an intermediate oxidation state closer to +3 than 0. The lack of an Al K-edge shift with the incorporation of Sr suggests that changes in Sr content do not have a meaningful impact on the electronics of the Al tetrahedra. Investigation of the solid solution structures provides evidence for classifying this structure type as a polar intermetallic phase with variable composition. Magnetization measurements were collected for the solid solution and show complex magnetic ordering with competing ferromagnetic and antiferromagnetic interactions as the Sr content increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗