Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface waters”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Self-Assembly and Phase Separation for Transport: A brief argument for the continued exploration of liquid crystal flows and electrodeposition in micro-gravity

Phase separation, interfacial tension gradients, and molecular self-assembly are all phenomena frequently observed when manipulating soft condensed matter systems on Earth. Depending on the system these processes can also be quite predictable. For instance, on Earth water droplets placed on hydrophobic surfaces will have contact angles > 90° with limited spread, contrary to their wetting behavior on hydrophilic surfaces. However, in an environment that is in a constant state of free-fall, where the gravitational field is reduced by 10%, such as on the International Space Station (ISS), buoyancy and its effect on fluid convection also become negligible. Fluid flow and phase behavior become less predictable. Thus, continuing the prior example, even if water droplets still remain in contact with a surface in microgravity (µg), the interface between the air and droplets will still remain significantly curved, appearing almost as if they were placed over completely nonwetting surfaces (Fig. 1a). Moreover, the negligible opposing buoyant force would also mean that if immiscible liquid droplets of a less dense character, or air bubbles even (Fig. 1b), were to be suspended within a larger water drop, these immiscible droplets or air bubbles would not migrate toward any one direction unless either artificially propelled by an external force or internal reaction inducing flow.

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

Methanol at Water–TiO 2 Interfaces: Free Energies of Water and Methanol Dissociation

Methanol adsorption on TiO 2 surfaces has long been studied due to its role in enhancing photocatalytic hydrogen evolution, yet how it modulates surface chemistry under aqueous conditions remains little understood. Using molecular dynamics with an ab initio-based deep neural network potential, we find that methanol adsorption induces markedly different effects on the aqueous surfaces of anatase and rutile, the two common phases of TiO 2 . In anatase, methanol adsorption significantly enhances water dissociation, which is otherwise rare at the neat water interface. This enhancement arises from an alternative dissociation pathway mediated by surface-bound methoxyl groups. In contrast, methanol adsorption tends to suppress water dissociation on rutile, replacing it with thermodynamically favored methanol dissociation. Overall, methanol adsorption in an aqueous environment alters not only the availability of key reactive intermediates involved in hydrogen evolution but also the hydrogen source, which turns out to be primarily methanol on rutile, whereas both water and methanol are consumed on anatase. These results provide mechanistic insights into the coupled roles of organic adsorbates and water at photocatalytic interfaces, with implications on how methanol enhances the activity of H 2 evolution.

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

Constraining the Hydration of Clay Minerals and Abundances of Amorphous Phases in Gale Crater, Mars

Both water and organic matter are required for the development and persistence of life. Phyllosilicates (clay minerals) have high surface areas that easily sorb water and organic matter. The Curiosity rover has investigated several hundred meters of stratigraphy in Gale crater, including where clays were detected from orbit. Previous results have suggested that subsurface hydration is greatest in units with the most abundant clays, suggesting that these minerals may be hydrated. Organics have also been found throughout Gale crater. Smectites are the most common and abundant phyllosilicates in Gale crater samples and can expand and sorb water and organics in interlayer sites. The most common organic sorption processes on Earth typically involve water or hydroxyl, so hydrated phyllosilicates are good candidates for organic preservation. Using newly derived subsurface hydration results with previously published mineralogy and geochemistry, we derived modeled constraints on the abundances of hydrated amorphous phases, “excess” water, and “excess” cations. These “excess” phases are not accounted for by published crystalline phase abundances or by amorphous phases constrained here. We found correlations between smectites and both “excess” water and “excess” cation abundances, indicating that smectites in Gale crater are hydrated and that cation bridging could be a mechanism for sorption of organics. Our results also show the persistence of amorphous sulfates, opal-A, and volcanic or impact glass, which indicate low water-rock interactions. Increased abundances of sulfates and glass in stratigraphically higher samples may indicate lower water availability and environmental aridification during the time these units were being deposited.

58 GEOSCIENCES↗
DOE OSTI2021

Performance of III–V Solar Cells Grown on Reformed Mesoporous Ge Templates

We demonstrate a solar cell on reformed porous Ge with an efficiency of 7.7%. We generate mesopores in (100) Ge by bipolar electrochemical etching and anneal them at high temperature. The pores coalesce deep in the structure rather than at the surface as desired, although resulting in coarse superficial morphology unsuitable for device growth. To combat this issue, we developed a surface treatment involving an HBr dip, annealing at 415 degrees C, and a postannealing ultrasonic De-ionized water dip to improve the surface structure, resulting in a smoother reformed surface on which we grow a GaInAs solar cell. The structure retains embedded pores after growth and the transitions between Ge and III-V layers are distinct. The solar cell fabricated using the improved coalescence has an efficiency of 4.5%. The efficiency improves to 7.7% by isolating the rest of the device from three limiting localized shunt areas. Protruding defects in the porous Ge and III-V layers still limit the performance, but this work establishes a step toward the technical viability of this exfoliation approach, showing decent efficiency if protruding defects can be removed or reduced.

14 SOLAR ENERGY↗
DOE OSTI2024

Patch-level CO2 and CH4 fluxes and porewater concentrations in experimental wetlands, 2 PPT saltwater intrusion simulations, Aug-Oct 2022: Louisiana

This dataset contains carbon dioxide (CO2) and methane (CH4) flux measurements from patches of wetland vegetation dominated by Typha domingensis and Panicum hemitomon, which were conducted to assess flux responses to acute saltwater intrusion. The measurements occurred before, during, and after simulated acute saltwater intrusion events of low concentrations of ~ 2 PPT. The measurements comprise gas fluxes from the wetland surface (i.e., soil-water column and vegetation) and fluxes from the soil-water column exclusively. These two sets of fluxes are separated into two files and are complemented with four more files containing porewater concentrations of CO2 and CH4 collected at 0-5 cm, 10-15 cm, and 20-25 cm depth increments, spectral indices measurements, biomass, and sediment elevation table measurements. The files can be opened with regular text editors or spreadsheet programs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

TEAMER - Extreme Events Modeling for the MARMOK-OWC Wave Energy Converter

Through the TEAMER program, Sandia National Laboratories (SNL) collaborated with IDOM Incorporated to study their MARMOK-Oscillating Water Column (MARMOK-OWC) wave energy conversion device. The study yielded a quantitative understanding of hydrodynamic pressures on the oscillating water column (OWC) device surfaces, the mooring tensions, and the dynamic performance of the device under extreme ocean wave conditions. This project utilized a comprehensive multi-phase Navier-Stokes flow solver with an overset body-fit mesh to predict fluid velocities and hydrodynamic forces on the MARMOK-OWC device. Computational Fluid Dynamics (CFD) analysis were conducted using OpenFOAM. This data includes the OpenFOAM cases (setup and data) to run the extreme events developed during the project. This project is part of the TEAMER RFTS 4 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

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

The Oceanic Barrier Layer in the Eastern Indian Ocean as a Predictor for Rainfall Over Indonesia and Australia

Barrier layers in the tropics trap heat in a shallow and stable near-surface layer and limit entrainment of cooler water from below. Both processes act to increase sea surface temperature and enhance atmospheric convection. The high resolution fully coupled pre-industrial Energy Exascale Earth System Model version 0 (E3SMv0) is used to investigate the relationship between barrier layers in the eastern Indian Ocean during the wet season with local atmospheric convection and remote rainfall. A partial least squares regression reveals a significant relationship between Australasian rainfall and the barrier layer thickness (BLT) west of Sumatra, occurring one month earlier. Furthermore, the largest positive regression coefficients are over northern Australia. The region west of Sumatra is strategically located where the East-Asian monsoon moisture flows toward northern Australia. Thickening of the west Sumatra BLT intensifies evaporation and local convection and amplifies the moisture transported to Australia acting to increase the terrestrial rainfall.

58 GEOSCIENCES↗
DOE OSTI2024

Perovskite/Perovskite Tandem Photoelectrodes for Low-Cost Unassisted Photoelectrochemical Water Splitting

In this project, we aim to address the challenges of achieving efficient and cost-effective unassisted photoelectrochemical water splitting using perovskite/perovskite tandem photoelectrodes. We utilized the following unique approaches to advance our innovation. (1) We developed stable and efficient low-E g (1.2 – 1.4 eV) perovskites as the bottom electrodes. We demonstrated approaches to improve the photovoltaic performance and photothermal stability of low-E g Sn-Pb iodide perovskite solar cells. (2) We developed a low-temperature synthesis route to fabricate wide-E g (>1.8 eV) perovskite top electrodes. By combining theoretical and experimental investigation, we found methods to reduce the formation of defects and dislocations in the wide-E g mixed halide perovskites and suppress halide segregation. (3) We developed a robust metal oxide-based interconnecting layer to integrate two perovskite layers into tandem photoelectrodes. Our monolithically integrated tandem devices feature a high V OC of more than 2 V and a high J SC of more than 15 mA/cm 2 . We showed the new design of the tandem photoelectrodes that is critical to the stable operation of tandem photoelectrodes for unassisted PV/PEC water splitting. (4) We demonstrated a water-impermeable barrier of carbon paste/epoxy/metal foil composite, which can prevent photocorrosion and water ingress of perovskite active layers. This surface protection enabled the operation of perovskite photoelectrodes in water and significantly enhanced the long-term stability of our tandem devices. (5) We conducted standardized PEC characterization in collaboration with NREL and reported accurate determination of solar-to-hydrogen conversion efficiencies of perovskite/perovskite tandem photoelectrodes. At the end of the project, we demonstrated perovskite/perovskite tandem photoelectrodes with STH efficiencies of up to 18% and less than 20% efficiency loss after continuous operation for more than 500 hours in water. Our results demonstrate the potential to develop a low-cost, durable, and efficient water-splitting system that meets the DOE 2026 and 2031 cost targets for hydrogen production.

08 HYDROGEN↗
DOE OSTI2021

Reversible electrowetting transitions on superhydrophobic surfaces

An electric field applied across the interface has been shown to enable transitions from the Cassie to the Wenzel state on superhydrophobic surfaces with miniature corrugations. Molecular dynamics (MD) simulations manifest the possibility of reversible cycling between the two states when narrow surface wells support spontaneous expulsion of water in the absence of the field. With approximately 1 nm sized wells between the surface asperities, the response times to changes in the electric field are of O(0.1) ns, allowing up to GHz frequency of the cycle. Because of the orientation preferences of interfacial water in contact with the solid, the phenomenon depends on the polarity of the field normal to the interface. The threshold field strength for the Cassie-to-Wenzel transition is significantly lower for the field pointing from the aqueous phase to the surface; however, once in the Wenzel state, the opposite field direction secures tighter filling of the wells. Considerable hysteresis revealed by the delayed water retraction at decreasing field strength indicates the presence of moderate kinetic barriers to expulsion. Known to scale approximately with the square of the length scale of the corrugations, these barriers preclude the use of increased corrugation sizes while the reduction of the well diameter necessitates stronger electric fields. Field-controlled Cassie-to-Wenzel transitions are therefore optimized by using superhydrophobic surfaces with nanosized corrugations. Abrupt changes indicate a high degree of cooperativity reflecting the correlations between the wetting states of interconnected wells on the textured surface.

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

Panel Session 4 and 15: Japan Fukushima Daiichi D and D Update and Technological Challenges at Japan Fukushima Daiichi D and D and Update on Nuclear Overview and Development in Japan - Nuclear Fuel Cycle

Sessions 4 and 15 represent a two-part panel series that discusses progress and challenges associated with cleanup at Fukushima. Severe limitations on availability of original panelists from Japan due to strict restrictions put in place to alleviate the spread of coronavirus necessitated changes to both panels. The result was a significantly modified panel for Session 04 (shown above) and the elimination of all panelists for Session 15. The 4 and 15 Panel Sessions provide an overview of activities related to both the progress and challenges of cleanup and decommissioning of the Fukushima Daiichi Nuclear Power Station (NPS) in Japan. Five panelists discussed perspectives of the cleanup following a Tokyo Electric Power Company (TEPCO) video showing the progress on site since the devastating Great East Earthquake and tsunami that caused the explosions at three of the six reactors on the site. Three of the five panelists discussed on-going work being performed for the effort, while the other two provided expert perspectives of on strategic efforts at the site. The panel was attended by over 80 technologists and policy makers spanning the globe and was opened by Dr. Monica Regalbuto of Idaho National Laboratory and a short video provided by TEPCO. The video described changes at the site that spanned the cleanup efforts from stabilizing water intrusion into the contaminated reactor buildings to construction of new administrative facilities. The video explained the processes underway to retrieve spent fuel rods and challenges in retrieval of the compromised fuel debris. The video highlighted working condition improvements that included establishment of rest housing and a small convenience store on the site, and the rollback of protective equipment around the site due to decreases contamination. Panelists with presentations: Revision of 'the Mid-and-Long-Term Road-map towards the Decommissioning of TEPCO's Fukushima Daiichi Nuclear Power Station' (Paul Dickman); Sharing UK experience at Fukushima Daiichi (Adrian Simper); SRNL Japan (Andrew Fellinger); ABLE's Initiative to Dismantle the Exhaust Stack (Daniel Walter); JAEA R and D in Fukushima (Tokio Fukahori); TEPCO - Overview and Update of the Fukushima Decommissioning Process (Monica Regalbuto); Toshiba's Involvement in the Decommissioning of the Fukushima Daiichi Nuclear Power Plant (Yasuhiro Yuguchi); Remote Dismantling of the Exhaust Stack At Fukushima Dai-ichi NPS (Takashi Okutsu)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Characterization of Lanthanum Monazite Surface Chemistry and Crystal Morphology through Density Functional Theory and Experimental Approaches

Monazite is a rare earth element (REE)-containing mineral that consists of (REE)PO 4 formal units and is one of the most important sources of these critical materials. The concentration of REEs from mined monazite ore often involves froth flotation, which is a beneficiation process that enhances the efficiency of downstream processing. The effectiveness of froth flotation is largely governed by the ability of collector agents to selectively bind to monazite particles. Thus, a molecular-level understanding of monazite interfacial chemistry is integral to the design of effective collector agents. To address this need, we performed density functional theory (DFT) calculations and a variety of experimental techniques to characterize La-monazite and elucidate its crystal morphology. Interestingly, we find minimal differences in the predicted morphologies of La-monazite for hydrous and anhydrous environments, which are largely dominated by low-index facets (e.g., {110}, {100}, and {010}). Indexing of synthesized La-monazite crystals via X-ray diffraction also uncovers {110} and {100} as the predominant facets. The average surface energies of 0% and 100% water coverage La-monazite crystals were predicted to be 0.87 and 0.76 J/m 2 , respectively, while calorimetry suggests values of 1.30 and 1.15 J/m 2 , respectively. The apparent discrepancies between the theoretical and experimental values are expected and attributed to defects present in physical crystals, in contrast to the perfect mineral surfaces in simulations. The difference in surface energy between the 0% and 100% water coverage morphologies predicted by theory is consistent with the value measured via calorimetry. DFT reveals a wide range of adsorption energies for water across the studied facets, but in all cases, water is predicted to strongly bind to monazite surfaces with an average adsorption energy of -92.7 kJ/mol for a La-monazite single crystal. Finally, this study provides the groundwork necessary for the rational design of froth flotation collector agents by granting molecular-level insight into the predominant facets of monazite.

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

The role of surface hydroxyls on the radiolysis of gibbsite and boehmite nanoplatelets

Understanding mechanistic pathways to radiolytic hydrogen generation by metal oxyhydroxide nanomaterials is challenging because of the difficulties of distinguishing key locations of OH bond scission, from structural interiors to hydroxylated surfaces to physi-sorbed water molecules. In this work we exploited the interface-selectivity of vibrational sum frequency generation (VSFG) to isolate surface versus bulk hydroxyl groups for gibbsite and boehmite nanoplatelets before and after 60 Co irradiation at dose levels of approximately 7.0 and 29.6 Mrad. While high-resolution microscopy revealed no effect on particle bulk and surface structures, VSFG results clearly indicated up to 83% and 94% radiation-induced surface OH bond scission for gibbsite and boehmite, respectively, a substantially higher proportion than observed for interior OH groups by IR and Raman spectroscopy. Electron paramagnetic spectroscopy revealed that the major radiolysis products bound in the mineral structures are trapped electrons, O , O 2 – and possibly F-centers in gibbsite, and H , O and O 3 – in boehmite, which persist on the time frame of several months. The entrapped radiolysis products appear to be highly stable, enduring re-hydration of particle surfaces, and likely reflect a permanent adjustment in the thermodynamic stabilities of these nanomaterials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Effects of Surface Terminations of 2D Bi 2 WO 6 on Photocatalytic Hydrogen Evolution from Water Splitting

Two-dimensional (2D)-structured photocatalysts with atomically thin layers not only have the potential to enhance hydrogen generation efficiency but also allow more direct investigations of the effects of surface terminations on photocatalytic activity. In this work, taking 2D Bi 2 WO 6 as a model, we found that the configuration of bilayer Bi 2 O 2 sandwiched by alternating WO 4 layers enabled the thermodynamic driving potential for photocatalytic hydrogen evolution. Without Pt deposition, the H 2 generation efficiency can reach to 56.9 μmol/g/h by 2D Bi 2 WO 6 as compared with no activity of Bi 2 WO 6 nanocrystals under simulated solar light. This configuration is easily functionalized by adsorption of Cl – /Br – to form Bi–Cl/Bi–Br bonds, which leads to the decrease of recombination in photogenerated charge carriers and narrower band gaps. This work highlights an effective way to design photocatalysts with efficient hydrogen evolution by tuning the surface terminations.

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

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗