Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAILS”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Wastewater hydrogen nexus (WwHeN): Greening the wastewater industry via integration with the hydrogen economy

As industries rally toward achieving net zero emissions, hydrogen and hydrogen-based fuels are emerging as key players in decarbonization efforts. Although the primary technology for producing renewable natural gas has been well implemented in the wastewater industry, the “decarbonization based goal setting” is trailing. This perspective assimilates existing literature presented in other contexts to highlight the need for framing the decarbonization dialog by using green hydrogen as a potential pathway for the wastewater industry. Specifically, we note the importance of (a) developing the decarbonization or net zero focus in the wastewater industry, and (b) colocating the wastewater industry with hydrogen production facilities. Finally, we also delve into technological, cost, and operational considerations to understand the readiness level of key stakeholders to identify future research and development opportunities for the wastewater hydrogen nexus.

08 HYDROGEN↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2025

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗
DOE OSTI2025

Catalytic Electrolyte Additive for High-Loading and Lean Electrolyte Li–S Batteries

The cycle life of high-energy Li–S cells is largely constrained by the quick electrolyte depletion. LiNO 3 has been a well-established additive known for protecting the Li metal anode and stabilizing the battery from polysulfide “shuttling”. However, it can be depleted prematurely and can pose safety risks when exposed to carbon, sulfur, or Li metal under harsh conditions. Here, in this study, LiPO 2 F 2 was explored as a safe and durable alternative additive in ether-based electrolytes. LiPO 2 F 2 demonstrates superior performance in Li/S batteries, especially under high sulfur loading (∼4 mg/cm 2 ) and lean electrolyte conditions (E/S = 4), achieving a long-term cycling stability of 40%, compared to 14.7% with LiNO 3 . This additive facilitates the disproportionation of polysulfides, reducing their dissolution and mitigating the shuttle effect. Additionally, LiPO 2 F 2 promotes the formation of a stable solid-electrolyte interphase (SEI) composed of inorganic anion-derived species, improving the battery’s overall stability and functionality. These findings blaze a trail in the design of safer and more durable electrolytes for Li–S batteries.

Li-S batteries↗
DOE OSTI2025

Cu-, Co-, and Zn-Based Metal–Organic Framework-Derived Nanoporous Ion Emitters for Picogram Level Analysis of Actinides

Thermal ionization mass spectrometry (TIMS) is often regarded as the preferred technique for trace-level isotopic analysis of actinides owing to its high sensitivity and absence of carry-over effects. However, actinide sample utilization efficiency (SUE) is typically low (<0.05%) without the use of activators or specialized loading approaches which can yield SUEs of >5%. To this effect, we investigate a series of metal–organic framework (MOF)-based nanoporous ion emitters (nano-PIEs) that show enhanced ionization of actinides when used in TIMS loading. We study the impact of physical and chemical properties of MOFs on TIMS SUEs using two families of MOFs that can be synthesized under similar reaction conditions. The structural and chemical properties of these MOFs can be systematically modified one at a time while keeping other features the same. This allows us to strategically investigate their impact on SUEs. The first family of MOFs considered in this study is Zeolitic Imidazole Frameworks (ZIFs) which are made using 2 methyl-imidazole as an organic linker with zinc (ZIF-8) and cobalt (ZIF-67) as metal centers. Additionally, the effect of morphology was also studied using Zn-based ZIF-L with a 2-dimensional structure. The second family of MOFs was synthesized using benzene-1,3,5-tricarboxylate (BTC) with copper (Cu-BTC) and zinc (Zn-BTC) as the metal center. Among the MOFs tested, Cu-BTC showed the highest SUE with an average SUE of 0.27 ± 0.15%, followed closely by Zn-BTC (0.24 ± 0.10%), ZIF-8 (0.17 ± 0.10%), and trailed by the other MOFs. When the MOFs were pyrolyzed in N 2 before loading, an apparent increase in the SUE was observed with ZIF-8 and Cu-BTC showing average SUEs of 0.25 ± 0.08% and 0.34 ± 0.13%, respectively, with the highest measured SUE of 0.53% for pyrolyzed Cu-BTC. In conclusion, this observed increase in SUE by up to an order of magnitude compared to bare filaments demonstrates the efficacy and potential of MOF-derived nano-PIEs for TIMS application.

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

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Shrub Expansion Can Counteract Carbon Losses From Warming Tundra

Arctic warming is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and tree-line expansion and densification. However, predicting the carbon trajectories of the changing Arctic is challenging due to interacting feedbacks between vegetation composition and structure, and surface characteristics. We conduct a sensitivity analysis of the current-date to 2100 projected surface energy fluxes, soil carbon pools, and CO 2 fluxes to different shrub expansion rates under future emission scenarios (intermediate—RCP4.5, and high—RCP8.5) using the Arctic-focused configuration of E3SM Land Model (ELM). We focus on Trail Valley Creek (TVC), an upland tundra site in the western Canadian Arctic, which is experiencing shrub densification and expansion. We find that shrub expansion did not significantly alter the modeled surface energy and water budgets. However, the carbon balance was sensitive to shrub expansion, which drove higher rates of carbon sequestration as a consequence of higher shrubification rates. Thus, at low shrub expansion rates, the site would become a carbon source, especially under RCP8.5, due to higher temperatures, which deepen the active layer and enhance soil respiration. At higher shrub expansion rates, TVC would become a net CO 2 sink under both Representative Concentration Pathway scenarios due to higher shrub productivity outweighing temperature-driven respiration increase. Our simulations highlight the effect of shrub expansion on Arctic ecosystem carbon fluxes and stocks. We predict that at TVC, shrubification rate would interact with climate change intensity to determine whether the site would become a carbon sink or source under projected future climate.

Yazbeck, Theresia [The Ohio State Univ., Columbus,↗
DOE OSTI2026

Deep Convection‐Driven Downward Transport of Trace Gases and Aerosols From the Free Troposphere to the Boundary Layer

Deep convective clouds regulate Earth's energy and moisture budgets, yet their impact on the atmospheric boundary layer (BL) composition remains underexplored. Using long-term observations from three mid-latitude sites, we show that deep convection (DC) consistently enhances nighttime surface ozone and is often accompanied by modest increases in ultrafine particle concentrations. Within the BL, the condensational growth of these transported ultrafine particles may contribute up to 60% of total cloud condensation nuclei (CCN). Mass flux calculations suggest that short-lived convective cores (∼30 min) account for ∼2% of total vertical air mass transport relative to steady entrainment, increasing to ∼13% when the trailing stratiform regions are included. These results show that DC provides an episodic but efficient pathway linking the free troposphere and BL, influencing oxidant budgets, CCN variability, and climate forcing. Accurately representing this process in climate models may help reduce uncertainties in climate projections, under both preindustrial and present-day conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

X-ray phase-contrast imaging of strong shocks on OMEGA EP

The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. Furthermore, the interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗
DOE OSTI2024

Analytical noise bias correction for precise weak lensing shear inference

Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the ‘renoising’ noise bias correction approach introduced by METACALIBRATION. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin Legacy Survey of Space and Time (LSST) imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a per cent, meeting LSST Dark Energy Science Collaboration requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Direct Evidence of Multispecies Hydrodynamics in Ignition-Scale Hohlraums

A targeted experiment at the National Ignition Facility (NIF) confirms the presence of multispecies hydrodynamics in inertial confinement fusion hohlraums relevant to ignition. Here, the effects are identified by filling the gold hohlraum with a deuterium-tritium (DT) gas mixture instead of helium. As the hohlraum is heated by the NIF lasers, it implodes inward, compressing and heating the DT, which leads to fusion. The resulting DT-fusion neutrons are measured in space, time, yield, angle, and energy. A distinct, peaked, triangular shape in the radial neutron emission profile provides evidence of a “leaky piston” effect caused by the interpenetration of DT into the expanding gold. This process reduces the reversibility of DT compression and decreases neutron generation on the temporal trailing edge compared to the leading edge. These results are well described by multispecies hydrodynamics simulations, which reproduce the observed spatial and temporal features, as well as the total neutron yield and angularly resolved energy spectra. In contrast, conventional simulations that use only single-species physics fail to match the experimental data. They overpredict the neutron yield, produce a flatter-than-observed spatial profile, and show excessive emission at late times compared to the experimental measurements.

Higginson, Drew P. [Lawrence Livermore National La↗
DOE OSTI2024

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems

This work presents methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices play a critical role in many power system applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

large scale↗
DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗
DOE OSTI2026

Hydrometeor Transport and the Stratiform Water Budget in Amazonian Mesoscale Convective Systems

A significant portion of precipitation falls within the stratiform regions of mesoscale convective systems (MCSs). These regions are formed from convective clouds that have either weakened and formed a horizontally stratified structure or have advected hydrometeors to form a trailing-stratiform region, with further precipitation growth promoted by in situ dynamical and microphysical processes. Green Ocean Amazon 2014/5 (GoAmazon2014/5) field campaign observations are used to quantify the stratiform water budget of six MCSs (three each from the wet and dry seasons) utilizing the framework from Leary and Houze, with an emphasis on the transfer of hydrometeors (C A ) from the convective to the stratiform regions. The storms were predominantly leading-line/trailing-stratiform MCSs, but two had embedded convection in the stratiform region, which can complicate water budget calculations. A pulse-like pattern with a mean interval of 21 min occurs in C A during each storm resulting from a lag between the initial convective updrafts and the subsequent hydrometeor transport and sedimentation in the stratiform region. Here, the peaks in C A were followed by peaks in hydrometeor growth associated with the mesoscale updraft (C mu ) in the stratiform region. Averaged over the six storms, the positive contributions of C A and C mu to the stratiform water budget were roughly equal, although storm-to-storm variations existed. Negative contributions from evaporation below cloud base in the mesoscale downdraft and detrainment from the stratiform region into the nonprecipitating anvil cloud and environment accounted for 53% and 21% of the sum of C A and C mu , respectively, resulting in a mean stratiform region precipitation efficiency of 26%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Moisture Recycling through Pumping by Mesoscale Convective Systems

Abstract Moisture recycling, the contribution of local evapotranspiration (ET) to precipitation, has been studied using bulk models assuming a well-mixed atmosphere. The latter is inconsistent with a climatologically stratified atmosphere that slants across latitudes. Reconciling the two views requires an understanding of overturning associated with different weather systems. In this study, we aim to better understand moisture recycling associated with mesoscale convective systems (MCSs). Using a convection-permitting WRF simulation equipped with water vapor tracers (WRF-WVT), we tag moisture from terrestrial ET in the U.S. Southern Great Plains during May 2015, when more than 20 MCS events occurred and produced significant precipitation and flooding. Water budget analysis reveals that approximately 76% of terrestrial ET is advected away from the region while the remaining 24% of terrestrial ET is “pumped” upward within the region, accounting for 12% of precipitation. Moisture recycling peaks during early night hours (1800–2400 LT) due to the mixing of the daytime accumulated ET by active convection. By focusing on five “diurnally driven” MCSs with less large-scale circulation influence than other MCSs during the same period, we find an upright pumping of terrestrial ET at the MCS initiation and development stages, which diverges into two branches during the MCS mature and decaying stages. One branch in the upper level advects the ET-sourced moisture downstream, while the other branch in the mid-to-upper level contributes to the trailing precipitation upstream. Overall, our analysis depicts a pumping mechanism associated with MCSs that mixes local ET vertically, highlighting its specific contributions to enhancing convective precipitation processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Stream discharge and temperature data collected within the East and Taylor Watershed, Colorado for the Lawrence Berkeley National Laboratory Watershed Function Science Focus Area (water years 2019 to 2025)

This dataset contains stream discharge and temperature data for water years 2019 to 2025 from the East and Taylor Watersheds in Colorado, United States. This data was collected to understand hydrological processes occurring in the East River and Taylor River Watersheds, Colorado, which is part of the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area. Data includes instantaneous observed discharge using salt dilution and acoustic doppler velocimeter techniques, raw pressure transducer downloaded data, sub-hourly temperature as well as corrected water level and associated stream discharge and mean daily values. Notes on water level corrections, rating curve development and metadata provided. A rating curve is the translation of depth to streamflow. The rating curve can be used as a quantitative measure of the “quality of the data.” Data within this dataset is formatted using ESS-DIVE’s Hydrological Monitoring Reporting Format. This data package contains (1) a zip file (Stream_Discharge_Data_WY19-WY25.zip) containing stream discharge and temperature data organized by location; (2) an InstallationMethods file (InstallationMethods.csv) describing metadata about the installation; (3) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; (4) a data dictionary (dd.csv) file that contains terms/column_headers used throughout the files along with a definition, units, and data type; (5) a locations metadata file (locations.csv); (6) and a sensor metadata file (sensors.csv). All data files are in non-proprietary formats (csv, png, or pdf formats). Please contact Rosemary Carroll, Curtis Beutler, or Austin Shirley for any support in accessing the files. Update on 2023-05-12: Additional data from WYs 2021 and 2022 were added. Additionally, the dataset was converted using ESS-DIVE’s Hydrological Monitoring Reporting Format. Data files were reformatted to match reporting format guidance, new metadata files were added, and files were converted from excel to CSV. Update on 2025-05-16: Additional data from WYs 2022 (for locations not previously included), 2023, and 2024 were added. An additional descriptive PDF (WFSFA_Streamflow_Hydrograph_Disclaimer.pdf) was added. Metadata files were updated to reflect the addition of new data and locations. Update on 2026-05-18: Additional data from WY 2025 were added, including a new location Upper Trail Creek (TR-TCG2). Metadata files were updated to reflect the addition of new data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Integrated Hourly Meteorological Database of 20 Meteorological Stations (1981-2022) for Watershed Function SFA Hydrological Modeling

This dataset contains (a) a script “R_met_integrated_for_modeling.R”, and (b) associated input CSV files: 3 CSV files per location to create a 5-variable integrated meteorological dataset file (air temperature, precipitation, wind speed, relative humidity, and solar radiation) for 19 meteorological stations and 1 location within Trail Creek from the modeling team within the East River Community Observatory as part of the Watershed Function Scientific Focus Area (SFA). As meteorological forcings varied across the watershed, a high-frequency database is needed to ensure consistency in the data analysis and modeling. We evaluated several data sources, including gridded meteorological products and field data from meteorological stations. We determined that our modeling efforts required multiple data sources to meet all their needs. As output, this dataset contains (c) a single CSV data file (*_1981-2022.csv) for each location (20 CSV output files total) containing hourly time series data for 1981 to 2022 and (d) five PNG files of time series and density plots for each variable per location (100 PNG files). Detailed location metadata is contained within the Integrated_Met_Database_Locations.csv file for each point location included within this dataset, obtained from Varadharajan et al., 2023 doi:10.15485/1660962. This dataset also includes (e) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata and (f) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. Review the (g) ReadMe_Integrated_Met_Database.pdf file for additional details on the script, methods, and structure of the dataset.The script integrates Northwest Alliance for Computational Science and Engineering’s PRISM gridded data product, National Oceanic and Atmospheric Administration’s NCEP-NCAR Reanalysis 1 gridded data product (through the `RCNEP` R package, Kemp et al., doi:10.32614/CRAN.package.RNCEP), and analytical-based calculations. Further, this script downscales the input data into hourly frequency, which is necessary for the modeling efforts.

54 ENVIRONMENTAL SCIENCES↗