SEARCH · Engineering Papers

Results for “High resolution data analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗

DOE OSTI2024

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗

DOE OSTI2026

Statistical White-Line Analysis in High-Throughput TXM-XANES for Chemical State Quantification

The transmission X-ray microscopy (TXM) based X-ray absorption near-edge structure (XANES) technique provides three-dimensional mapping of element-specific chemical states at nanometer-scale spatial resolution and micrometer-scale fields of view. However, compared to conventional volume-averaged XANES (VA-XANES) measurements, the inherently small voxel size in TXM-XANES leads to a lower signal-to-noise ratio, making full-spectrum analysis computationally demanding and less robust. Here, we present the structural and compositional conditions for a statistical white-line analysis framework under which chemical state information can be directly extracted from the white-line peak position in voxel spectra without the need for voxel-wise background subtraction or normalization, under well-defined structural and compositional conditions. The method is validated on layered oxide cathode materials, where low-order polynomial fitting accurately reproduces white-line features, and the extracted energy distributions correlate strongly with VA-XANES results. This statistical approach enables high-throughput, dose-efficient, and noise-robust chemical state quantification in TXM-XANES, offering broad applicability to functional materials requiring nanoscale oxidation-state mapping.

TXM↗

DOE OSTI2020

Climate-Related Trends of Within-Storm Intensities Using Dimensionless Temporal-Storm Distributions

Huff curves are probabilistic time distributions of rainfall expressed as dimensionless cumulative percentages of storm depth and duration. Previous studies have documented development factors, spatial robustness, and the utility of Huff curves in practical applications. However, the effects of trending rainfall on Huff curve intensity patterns have not yet been studied. As such, the goal of this paper is to fill this gap by studying Huff curve patterns in a watershed with demonstrated increasing trends of temperature and precipitation, with the intention that it can be generalized to other areas in the US and the world. To achieve this goal, the high temporal resolution precipitation data collected from a high spatial density, 72-year precipitation-gauge network on the 4.25-km 2 North Appalachian Experimental Watershed in east-central Ohio were used. Seasonal storm pattern trends from 1939 to 2010 were investigated using dimensionless depth (with the frequency of 50%, $d_{50}$) and the curve variability $(V = d_{80} - d_{20})$ at three dimensionless within-storms time periods (three verticals). The Spearman rank correlation procedure (correlation coefficient, ρ and significance probability, p) was used to statistically determine trends over time using 8 periods of 4-season sets of Huff curves over the 72 years. Two averages of ρ and p were computed: (1) by averaging the individual ρ and p obtained from 10 gauges (AvgI); and (2) by grouped averaging of all individual gauge values of $d_{50}$ and V and then computing ρ and p (AvgG). The test results of individual gauges showed that 4 cases for $d_{50}$ and 23 cases for V were significant for all seasons and verticals (total of 120 cases for each variable). The test results of AvgI for $d_{50}$ and V and AvgG for $d_{50}$ showed no significant trends in all seasons and verticals. Only the AvgG for V led to a significant trend for V in spring and fall at different times within storm patterns. The data do not provide sufficient evidence at the p=0.05 significance level to reject the null hypothesis of unchanging position of the dimensionless depth of the 50% Huff curves for individual or averages for all seasons and verticals. Also, there is insufficient evidence to reject the null hypothesis of unchanging variability, V, using the average ρ of individual gauges (AvgI for V) for all seasons and verticals. AvgG results showed a significant trend in V; however, this analysis may be affected by the nonindependence of storm data. The results suggest that it is likely that there is little if any effect of trending climate over approximately 70 years on Huff curve patterns. The results of this study add to the robustness characteristics of Huff curves and to their potential use in hydrological practice as design storms, as the foundation of stochastic storm generation, and for storm disaggregation, and they deserve further investigation. These different forms of inputs to watershed models have the potential to improve runoff estimation. The results suggest that, if verified in other studies, they have applicability to provide useful stationary precipitation patterns across the US and other areas of the world in areas of nonstationary climate. Also, individual rain gauge data may not be representative of trends even over small areas, and seasonal differences were noticeable as found in other studies. Recommendations are provided.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

The Influence of Soil Properties on Sea-Breeze Circulations in the Southeast U.S.

Sea-breeze circulations (SBCs) are common weather phenomena at and near coastal regions. They form because of a thermal gradient between the land surface at the coast and the sea surface. In a mid-day regime, a “thermal low” generated at the warm coast will lead to rising air motion, creating a wind shift coming from the sea near the surface displacing the coastal air. A “return flow” moving back towards the sea is generated by upper-level divergence because of the rising motion from the thermal low. SBCs propagate and serve as a method of urban pollutant dispersion in the Los Angeles region of California and are constrained to the coast due to the topography of surrounding mountains serving as a boundary for further inland propagation. Within the northeast U.S. SBCs are seen in the warm season but tend to remain coastally bound due to Coriolis distortion over long distances. Within the southeast U.S. (SEUS), the paradigmatic example of SBCs occurs over the Florida peninsula, where thunderstorms form on a nearly daily occurrence due to the convergence of SBCs from the east and west sides of the peninsula. However, there are further examples of sea-breezes in the SEUS that warrant study. Within the region bordering the SEUS and the Mid-Atlantic, just east of the southern Appalachian Mountains, warm-season SBCs form at the coast of Georgia and the Carolinas. Relatively flat topography ~150-200km inland allows for mostly unimpeded inland SBC propagation. Through visual analysis, Viner et al. catalogued several SBCs that propagated as far inland as the Central Savannah River Area surrounding Augusta, Georgia. Wermter et al. found that while the land-sea thermal gradient at the coast can influence coastal SBC genesis, the inland-coastal thermal gradient over the land is the primary influencer on the speed and depth of inland propagation of SBCs in this region. Additionally, soil moisture itself is a known correlative factor to sea-breeze formation, as it influences the soil temperature and the thermal gradient needed for SBC genesis and inland penetration. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Conversely, drier soil enhances the thermal gradient and promotes SBC formation. However, while there is an inverse relationship between soil moisture and SBCs, higher soil moisture can actually promote more convective rainfall following a SBC if it does not significantly impact the thermal gradient. While the relationship between soil moisture and SBC formation has been conceptually explored and modeled numerically, there is a research gap in observed connections. The Soil Moisture Active Passive (SMAP) satellite mission has been operational since 2015 and has been used to create high-resolution re-analytical Level 4 (L4) datasets of soil moisture and soil temperature at different soil depths: the surface (0-5cm) and rootzone (0-1m). The surface soil temperature effectively acts as the “skin temperature” of the surface at these levels, and a spatial map of the land-sea as well as the coastal-inland thermal gradients can be represented. SMAP data are also assimilated in some atmospheric models such at the High Resolution Rapid Refresh (HRRR) mesoscale model. We propose leveraging the use of SMAP products to fill in spatial gaps left by weather and mesonet stations within the SEUS region, as well as assessing the effectiveness of utilizing SMAP products towards SBC forecasting in both deterministic and machine learning (ML) models.

58 GEOSCIENCES↗

DOE OSTI2026

Final DOE-ASR Report for the Project “Using LASSO to bridge the gap between model and observations and to learn about atmospheric convection”

Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Fermi-LAT View on Three Ultra-high-energy 1LHAASO Sources in the 52° < l < 55° Region

Using more than 17 yr of Fermi-LAT data, we performed a detailed investigation of the complex 52° < l < 55° region, which encompasses the three ultra-high-energy sources 1LHAASO J1928+1746u, 1LHAASO J1928+1813u, and 1LHAASO J1929+1846u. This region hosts multiple supernova remnants (SNRs), pulsars, GeV and TeV sources. Our analysis resolves the GeV emission into three pointlike sources (J1925+1729P, J1930+1851P, and J1932+1916P) and two extended sources (J1929+1732E and J1930+1826E), and improves significantly on the description based on the 4FGL-DR4 catalog. Source J1932+1916P is identified as the known gamma-ray pulsar PSR J1932+1916, while J1925+1729P may be a new gamma-ray pulsar candidate distinct from the known gamma-ray pulsar PSR J1925+1720. This warrants future investigation and a search for pulsations. Source J1930+1851P coincides with the TeV source PWN/SNR G54.1+0.3, and its GeV−TeV spectrum is consistent with both leptonic and hadronic interpretations, although a leptonic origin in relation to the known pulsar wind nebula (PWN) is more likely. The GeV−TeV spectrum of J1929+1732E is consistent with a hybrid lepto-hadronic scenario in which the TeV emission traces the PWN powered by the pulsar PSR J1928+1746, while the GeV emission may result from interactions between particles escaped from the parent SNR and illuminating the gas environment. Similarly, J1930+1826E is likely connected to PWN/SNR G54.1+0.3 under a hadronic scenario involving escaped particles in their early propagation stage. Owing to spectral and/or morphological mismatches, the connection of these five GeV sources to the three LHAASO sources is not clear. This warrants deeper observations with HAWC and LHAASO, and a dedicated study of the modeling of the Galactic diffuse emission. Future CTAO observations with higher angular resolution are expected to deliver crucial information for the study of this region.

Liu, Linjie [Chinese Academy of Sciences (CAS), Ku↗

DOE OSTI2020

Spectral Variability of VHS J1256–1257b from 1 to 5 μm

Multiwavelength time-resolved observations of rotationally modulated variability from brown dwarfs and giant exoplanets are the most effective method for constraining their heterogeneous atmospheric structures. In a companion paper, we reported the discovery of strong near-infrared variability in HST/WFC3/G141 light curves of the very red L-dwarf companion VHS J1256–1257b. In this paper, we present a follow-up 36 hr Spitzer/IRAC Channel 2 light curve together with an in-depth analysis of the Hubble space telescope (HST) and the Spitzer data. The combined data set provides time-resolved light curves of VHS1256b sampling 1.1 to 4.5 μm. The Spitzer light curve is best fit with a single sine wave with a period of 22.04 ± 0.05 hr and a peak-to-peak amplitude of 5.76 ± 0.04%. Combining the period with a previously measured projected rotational velocity (v sin i), we find that VHS1256b is most consistent with equatorial viewing geometry. The HST/G141+Spitzer spectral energy distribution favors a model with a T{sub eff} of 1000 K and low surface gravity with disequilibrium chemistry. The spectral variability of VHS1256b is consistent with predictions from partly cloudy models, suggesting that heterogeneous clouds are the dominant source of the observed modulations. We find evidence at the 3.3σ level for amplitude variations within the 1.67 μm CH{sub 4} band, which is the first such detection for a variable L-dwarf. We compare the HST/G141 time-resolved spectra of three red L-dwarfs with high-amplitude near-infrared rotational modulations and find that although their time-averaged spectra are similar, their spectroscopic variabilities exhibit notable differences. This diversity reinforces the advantage of time-resolved spectroscopic observations for understanding the atmospheres of brown dwarfs and directly imaged exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

An Indirect Search for Weakly Interacting Massive Particles in the Sun Using Upward-going Muons in NOvA

I present the first Dark Matter search results using the full data set collected with the upward-going muon trigger in NOvA. Weakly Interactive Massive Particles (WIMPs) are a theoretical non-baryonic form of Dark Matter. The nature of Dark Matter is one of the most exciting open questions in modern physics. Though its existence can be inferred by astrophysical evidence, its properties are not yet understood. If we assume that Dark Matter particles can produce Standard Model particles through their interactions, an indirect search can help shed light on this mystery.The NOvA collaboration has built a 14 kton, fine-grained, low-Z, total absorption tracking calorimeter at an off-axis angle to the NuMI neutrino beam. Even though the detector is optimized to observe electron neutrino appearance from a muon neutrino beam, it has a unique potential for more exotic searches given its excellent granularity and energy resolution and relatively low-energy neutrino thresholds. In fact, with an efficient upward-going muon trigger and sufficient background suppression offline, NOvA is capable of a competitive indirect Dark Matter search for low-mass WIMPs.The idea of the upward-going muon trigger is first to select high-quality muon tracks, then use the timing information of all of the hits of each track to estimate directionality. In this way, the background flux is suppressed by more than a factor of 10^5 at trigger level to a rate of approximately 1 Hz. To further optimize this search, we use only upward-going muons that point to the Sun, so our search occurs at night when the Sun is on the other side of the Earth. This strategy also allows us to use the time when the Sun is above the horizon as a control region to estimate the background. Ultimately, implementation of a cut based and maximum likelihood analysis provides a powerful tool for rejecting background and selecting a sample of neutrino-induced upward-going muons. The overall background rejection power achieved by the analysis is substantial and impressive. Starting with approximately 150,000 events per second, we reduced it to 40 events per year. Since no statistically significant excess was found, a 90\% C.L. upper limit on the expected muon flux of upward-going muons has been set using the upper limit on the number of events given the number of observed events in the signal region. Lastly, by assuming the theory behind the upward-going muon flux, a limit on the WIMP-nucleon spin-dependent cross-section in the Sun was estimated. Although the limits on the spin-dependent cross-section do not appear to be competitive with previous indirect Dark Matter searches, the upward-going muon flux limits are promising. The upward-going muon flux limits could extend these results to a broader class of models that are not specific to the dark matter theory but produce upward-going muons, leading to competitive results.

Principato, Cristiana↗

DOE OSTI2023

Advanced Analytical Methodologies Enable Feedstock Screening and Correlation to Pyrolysis Yields and Catalytic Upgrading

Pyrolysis of lignocellulosic biomass can serve as a powerful pathway for the generation of renewable fuels and industrial products. However, procurement of some biomass may be costly, especially when there is competition for alternative uses. Forestry and agricultural byproducts may serve as valuable pyrolysis feedstocks due to their high availability and low cost. Although waste products may fill demand for affordable pyrolysis feedstocks, they are often high in ash and extractives which contribute to unfavorable changes in pyrolysis vapor composition, as well as end products. The work presented here investigates the correlations between various pine anatomical fractions to understand their effect on pyrolysis oil composition for forest residues. Analytical pyrolysis coupled to direct analysis by molecular beam mass spectrometry (MBMS) provides real-time monitoring of pyrolysis vapor composition, even for compounds that are not amenable to gas chromatography. Statistical analysis of MBMS data from the fast pyrolysis (FP) of clean pine, cambium, twigs and branches, needles, bark, and a forest residue samples revealed that bark and needles made the greatest differences in the vapor composition of forest residues during pyrolysis and were enriched in resin acids, furfural derivatives, and ions that have been previously observed from catalysis corresponding to ash content. These results were also compared to catalytic fast pyrolysis (CFP) using a platinum on titania catalyst, and correlations between FP and CFP products were observed for the anatomical fractions. The forest residues and corresponding anatomical fractions from a stand of 23-year-old pine will be converted by fast pyrolysis in a bench-scale 2-inch fluidized-bed reactor, and statistical correlations between the bench-scale products and analytical-scale pyrolysis vapors will be shown. This work highlights the power of advanced analytical methodologies in determining the contributions of vapor components to catalytic pyrolysis outcomes for mixed feedstocks and demonstrates the use of analytical methodologies to better screen feedstocks prior to pyrolysis.

analytical↗

DOE OSTI2022

Synchrotron-based techniques for characterizing STCH water-splitting materials

Understanding the role of oxygen vacancy–induced atomic and electronic structural changes to complex metal oxides during water-splitting processes is paramount to advancing the field of solar thermochemical hydrogen production (STCH). The formulation and confirmation of a mechanism for these types of chemical reactions necessitate a multifaceted experimental approach, featuring advanced structural characterization methods. Synchrotron X-ray techniques are essential to the rapidly advancing field of STCH in part due to properties such as high brilliance, high coherence, and variable energy that provide sensitivity, resolution, and rapid data acquisition times required for the characterization of complex metal oxides during water-splitting cycles. X-ray diffraction (XRD) is commonly used for determining the structures and phase purity of new materials synthesized by solid-state techniques and monitoring the structural integrity of oxides during water-splitting processes (e.g., oxygen vacancy–induced lattice expansion). X-ray absorption spectroscopy (XAS) is an element-specific technique and is sensitive to local atomic and electronic changes encountered around metal coordination centers during redox. While in operando measurements are desirable, the experimental conditions required for such measurements (high temperatures, controlled oxygen partial pressures, and H 2 O) practically necessitate in situ measurements that do not meet all operating conditions or ex situ measurements. Here, we highlight the application of synchrotron X-ray scattering and spectroscopic techniques using both in situ and ex situ measurements, emphasizing the advantages and limitations of each method as they relate to water-splitting processes. The best practices are discussed for preparing quenched states of reduction and performing synchrotron measurements, which focus on XRD and XAS at soft (e.g., oxygen K-edge, transition metal L-edges, and lanthanide M-edges) and hard (e.g., transition metal K-edges and lanthanide L-edges) X-ray energies. The X-ray absorption spectra of these complex oxides are a convolution of multiple contributions with accurate interpretation being contingent on computational methods. The state-of-the-art methods are discussed that enable peak positions and intensities to be related to material electronic and structural properties. Through careful experimental design, these studies can elucidate complex structure–property relationships as they pertain to nonstoichiometric water splitting. A survey of modern approaches for the evaluation of water-splitting materials at synchrotron sources under various experimental conditions is provided, and available software for data analysis is discussed.

08 HYDROGEN↗

DOE OSTI2026

Hydrodynamic Analysis and Optimization of Aquantis Marine Turbine: Cooperative Research and Development (Final Report)

The primary aim of this proposal is to improve the accurate prediction of hydrodynamic performance and dynamic load responses of the AQ10 floating axial-flow tidal turbine with a tri-cat mooring configuration. The validation of reduced-order modeling approaches with high-fidelity model will be implemented. Additionally, the frequency response domain, Response Amplitude Floating Wind (RAFT) toolbox plus an optimizer expanded for marine hydrokinetic turbines under the Submarine Hydrokinetic And Riverine Kilo-megawatt. Systems (SHARKS) program will be used for designing and exploring different key design parameters (platform dimension, mooring layout and its parameters) of next marine hydrokinetic (MHK) turbine generation.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

A multi-scale temperature-based strategy to map hydrologic exchange flows in highly dynamic systems

Mapping and quantifying hydrologic exchange flows (HEFs) is critical to environmental monitoring and remediation at contaminated sites; however, these objectives are challenging in highly dynamic systems, e.g., along dam-regulated rivers, where HEFs vary rapidly. Direct seepage measurements are labor-intensive and difficult to automate, whereas indirect (e.g., thermal) and remote sensing methods have potential to allow continuous monitoring with limited field effort. We present a preliminary assessment of a multi-scale temperature-based strategy for monitoring HEFs along the Hanford Reach of the Columbia River, in eastern WA, United States. Five thermal methods were assessed. First, a vertical temperature profile (VTP) was installed into the streambed. The VTP data were analyzed using a data assimilation algorithm designed for automated real-time estimation in dynamic systems. Second, a thermal infrared (TIR) camera was used in roving surveys to identify seeps. Third, a TIR camera was stationed at the VTP site to collect images at 1-h intervals. Together, the two TIR datasets provided a basis to assess the potential for drone-based TIR. Fourth, temperature was measured at the sediment/water interface to assess fiber-optic distributed temperature sensing. Fifth, imagery from the ECOSTRESS satellite mission was acquired to assess the potential of spaceborne thermal monitoring. Based on our preliminary assessment, VTP, TIR, and bed temperature measurements provide complementary spatial coverage, temporal sampling, and resolution; these methods have potential for long-term, automated monitoring of HEFs. The publicly available spaceborne imagery, however, proved inadequate because of insufficient spatial resolution and data gaps resulting from cloud cover and revisit frequency.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

UHT-CAMANCHE: Ultra-High Temperature Ceramic Additively Manufactured Compact Heat Exchangers

The conceptual basis for this project is the convergence of advanced ultra-high temperature ceramic materials and additive manufacturing technologies to produce compact ceramic heat exchangers with complex internal flow path geometries. Task areas were broadly divided into materials and manufacturing development, heat exchanger design, component testing, and techno-economic analysis. Technical challenges included the design and commissioning of new test facilities, improving feature resolution and deposition rate of ceramic additive manufacturing techniques, establishing process-structure-property relationships in additively manufactured ultra-high temperature ceramics, and assessing high temperature materials compatibility in CO 2 environments. The primary candidate material evaluated in this work is a composite comprising zirconium diboride (ZrB2) with 30 vol. % silicon carbide (SiC) which was selected based on its desirable combination of high temperature mechanical properties, high thermal and electrical conductivities, and oxidation resistance. High solids loaded ZrB2-SiC pastes suitable for extrusion-based additive manufacturing were developed for the first time as part of this work. Materials compatibility studies indicate this material oxidizes in CO 2 to form a protective borosilicate scale which transforms to pure silica above 1000°C. Parts made by additive manufacturing displayed enlarged grain sizes produced by pressureless sintering as compared to hot-press sintering. Increases in microstructural coarseness have outsized effect on oxidation performance up to 1400°C due to incomplete oxidation of coarse large diameter SiC particles resulting in lower amounts of silica that apparently inhibit protective scale formation. Additive manufacturing as a forming technique did not appear to significantly affect thermal conductivity, hardness, or elastic modulus, though flexural strength was reduced by half or more as compared to traditionally hot-pressed materials. This effect was attributed to the presence of strength-limiting flaws (ca. 40 microns in size) originating from extrudate inhomogeneities that could potentially be eliminated with further process improvements. Attempts to attain economies of scale for production of multi-kilowatt scale heat exchangers by ceramic additive manufacturing proved difficult. Lack of automation and a modest extrusion rate while retaining fine feature resolution made the overall process labor intensive and limited experimental throughput. A number of full-scale components were taken through post-process heat treatments including drying, binder burnout, and sintering; however, none survived without significant flaws or cracks. Therefore, no operational data from a newly installed heat exchanger test loop were able to be obtained during the performance period. Continued research and development is recommended to improve economic feasibility of ceramic additive manufacturing by standardizing the use of advanced sensors, artificial intelligence, and automation tools to reduce associated labor costs and accelerate production rates. The materials and manufacturing techniques demonstrated in this work are likely to find applications in defense and energy applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗

DOE OSTI2020

COMBLE-ISLAS water isotopes in precipitation (COMBLEISLASISO)

During the measurement period, precipitation was sampled on daily to sub-daily basis with a sampling kit consiting of a clear plastic box (60x40x40 cm) mounted to a ground structure. Daily sampling lasted from 22 Feb to 24 Mar 2020. At a fixed time of day (10 UTC), samples were collected and the sampling box cleaned for the next sample. Additionally, during intense observations periods (IOPs) announced by ISLAS, ARMS scientists increased sampling frequency to sub-daily sampling. During the campaign period, 2 IOPs were conducted, with in total 12 samples taken. IOP1 lasted from 27 Feb 10 to 02 Mar 2020, and IOP2 lasted from 12 Mar to 14 Mar 2020. Solid precipitation in the box was then melted, transferred to an 8 ml glass vial, sealed and shipped to FARLAB (University of Bergen, Norway) for stable water isotope analysis. Measured stable isotope composition of precipitation samples will be related to conditions at other sampling locations, and at the evaporation site to inform about conservation of water isotope quantities. In the period from 22 Feb to 24 Mar, 23 daily samples were collected at 10 UTC for stable water isotope analysis. During the entire sampling period, a high variability of precipitation phase was encountered. Most samples were recovered in the form of graupel/snow (12 samples), as rain (7 samples), and the remaining samples fell as graupel, snow, rain or frozen rain. Precipitation amounts varied between <1 mm (26 Feb) and 140 mm (25 Feb) within the sampling box. Furthermore, higher-resolution sampling was conducted during 2 intense observation periods (IOPs). During IOP1, lasting from 27 Feb 10 UTC to 02 Mar 10 UTC, 7 samples were collected, all of them as solid phase precipitation. IOP2 lasted from 12 Mar 10 UTC to 14 Mar 09 UTC, and 5 samples (all as graupel/snow) were collected. Samples were processed according to FARLAB standard measurement procedures. In short, samples were transferred to 1.5 ml glass vials with rubber/PTFE septa (part #548-0907, VWR, USA). An autosampler (A0325, Picarro Inc) transferred ca. 2µl per injection into a high-precision vapourizer (A0211, Picarro Inc, USA) heated to 110°C. After blending with dry N2 (< 5 ppm H2O) the gas mixture was directed into the measurement cavity of a Cavity-Ring Down Spectrometer (L2140-i, Picarro Inc) for about 7 min with a typical water concentration of 20 000 ppm. Memory effects were reduced by two times measuring a vapour mixture at a mixing ratio of 50 000 ppm, obtained from 2 injections of 2 µl for 5 min at the beginning of each new sample vial. Thereafter, another 6 injections of 2 µl per sample were measured individually as described above, and averages of the last 5 injections were used for further processing. Three standards were measured at the beginning and end of each batch, including a drift standard DI2 (δD: -50.72±0.73 permil, δ18O: -7.63±0.10 permil), and for calibration the laboratory standards GLW (δD: -307.79±0.75 permil, δ18O: -40.02±0.07 permil) and EVAP2 (δD: 9.52±0.65 permil, δ18O: 1.81±0.13 permil). A detailed calibration report is included with the final uploaded data set.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗