SEARCH · Engineering Papers

Results for “Groundwater Flux”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Borehole Disposal of Radioactive Waste in Israel - 20388

Israel is assessing borehole disposal of radioactive waste. With limited geological options for disposal, intermediate-depth borehole disposal is being considered in the arid Yamin Plain region of the northeastern Negev desert at depths of several hundred meters below ground surface in the vadose zone. Unlike deep borehole disposal of several kilometers, which relies on emplacement below the depth of recirculating groundwater, the safety case for intermediate-depth borehole disposal relies more on the aridity of the vadose zone, the robustness of the waste package, and the other engineered barriers (e.g., seals, backfill materials) in the disposal borehole. As a tradeoff, the shallower depth may allow the use of a larger diameter borehole than would be possible for a deeper borehole. To study the suitability of the Yamin Plain region for borehole disposal, a small-diameter characterization borehole is being planned to retrieve core samples and to better understand the vadose zone geomechanical and hydrogeochemical properties and percolation flux. The information from the characterization borehole will inform the safety case, and together with performance assessment analyses will help to identify key areas of uncertainty and guide future research and development activities aimed at demonstrating the feasibility of the intermediate-depth borehole disposal concept in Israel. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Data for: Spatial access and resource limitations control carbon mineralization in soils

This dataset contains data and code used for the paper "Spatial access and resource limitations control carbon mineralization in soils", https://doi.org/10.1016/j.soilbio.2021.108427. Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale water-soluble soil organic matter (SOM) mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble carbon (acetate) and nitrogen (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Soil respiration was measured during the incubation, after which pore water was extracted from the saturated soils and analyzed for water soluble organic carbon concentrations and characterization. Our results showed that carbon (C) amendments increased the cumulative carbon dioxide (CO2) evolved from the soil cores, suggesting that the soils were C-limited. Drought and rewetting increased soil respiration, and there was a greater abundance of complex aromatic molecules in pore waters sampled from these soils. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We had hypothesized that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (pore water samples), drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules, suggesting that fresh inputs stimulated the C-and N-limited microbial populations present deeper in the soil profile. Our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/TES_spatial_access_2021 (DOI: 10.5281/zenodo.5522938).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Composite Analysis Solid Waste Release Data Reduction of Activity Flux from Waste Sites to the Vadose Zone for Baseline Assessment

This environmental calculation file (ECF) describes the data reduction methodology for the solid waste release datasets generated by the Composite Analysis Solid Waste Release (CASWR) model (CP-62766, Model Package Report: Composite Analysis Solid Waste Release Model [CASWR Model]) and documents the resulting reduced datasets for the updated Composite Analysis (CA) baseline assessment. This ECF is in support of the updated Hanford Site CA. The purpose of the updated Hanford Site CA is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The vadose zone modeling activities associated with the CA includes solid waste forms as contaminant sources. To support the modeling of solid waste transport from the vadose zone to the groundwater at the Hanford Site, the CASWR model was designed to generate deterministic radionuclide release rates for the Hanford Site Central Plateau solid waste disposal sites. The complete CASWR model-generated dataset for the CA Baseline Assessment consists of 2,378,831 sets of time and release rates (1,073,869 sets for the 200 East Area and 1,304,962 sets for the 200 West Area), representing a course of the simulation from year 1953 to year 12,069, 16 radiological contaminants of potential concern (COPCs), and 237 solid waste sources as documented in ECFHANFORD- 19-0112, Solid Waste Release Calculations for the Composite Analysis Baseline Assessment. For several of the STOMP vadose zone models, the number of release rate data pairs required for release at all waste sites would exceed the limitations imposed by the STOMP modelling software package used to simulate contaminant transport through the vadose zone. As a result, a data reduction of the CASWR dataset is required prior to its use as a source input for the CA vadose zone modeling activities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

SGD-OD: investigating the potential oxygen demand of submarine groundwater discharge in coastal systems

Submarine groundwater discharge (SGD) supplies nutrients, carbon, metals, and radionuclide tracers to estuarine and coastal waters. One aspect of SGD that is poorly recognized is its direct effect on dissolved oxygen (DO) demand in receiving waters, denoted here as SGD-OD. Sulfate-mediated oxidation of organic matter in salty coastal aquifers produces numerous reduced byproducts including sulfide, ammonia, dissolved organic carbon and nitrogen, methane, and reduced metals. When these byproducts are introduced to estuarine and coastal systems by SGD and are oxidized, they may substantially reduce the DO concentration in receiving waters and impact organisms living there. We consider six estuarine and coastal sites where SGD derived fluxes of reduced byproducts are well documented. Using data from these sites we present a semiquantitative model to estimate the effect of these byproducts on DO in the receiving waters. Without continued aeration with atmospheric oxygen, the study sites would have experienced periodic hypoxic conditions due to SGD-OD. The presence of H2S supplied by SGD could also impact organisms. This process is likely prevalent in other systems worldwide.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Anion Data for the East River Watershed, Colorado (2014-2025)

The anion data for the East River Watershed, Colorado, consist of fluoride, chloride, sulfate, nitrate, and phosphate concentrations collected at multiple, long-term monitoring sites that include stream, groundwater, and spring sampling locations. These locations represent important and/or unique end-member locations for which solute concentrations can be diagnostic of the connection between terrestrial and aquatic systems. Such locations include drainages underlined entirely or largely by shale bedrock, land covered dominated by conifers, aspens, or meadows, and drainages impacted by historic mining activity and the presence of naturally mineralized rock. Developing a long-term record of solute concentrations from a diversity of environments is a critical component of quantifying the impacts of both climate change and discrete climate perturbations, such as drought, forest mortality, and wildfire, on the riverine export of multiple anionic species. Such data may be combined with stream gauging stations co-located at each monitoring site to directly quantify the seasonal and annual mass flux of these anionic species out of the watershed. This data package contains (1) a zip file (anion_data_2014_2025.zip) containing a total of 386 files: 387 data files of anion data from across the Lawrence Berkeley National Laboratory (LBNL) Watershed Function Scientific Focus Area (SFA) which is reported in .csv files per location and a locations.csv (1 file) with latitude and longitude for each location; (2) a file-level metadata (v7_20260901_flmd.csv) file that lists each file contained in the dataset with associated metadata; (3) a data dictionary (v7_20260901_dd.csv) file that contains terms/column_headers used throughout the files along with a definition, units, and data type; and (4) a anion MDL fact sheet (anion_MDLs_202608 in PDF and docx formats). Missing values within the anion data files are noted as either "-9999" or "0.0" for not detectable (N.D.) data. There are a total of 47 locations containing anion data. Update on 2022-06-10: versioned updates to this dataset was made along with these changes: (1) updated anion data for all locations up to 2021-12-31, (2) removal of units from column headers in datafiles, (3) added row underneath headers to contain units of variables, (4) restructure of units to comply with CSV reporting format requirements, and (5) the addition of the file-level metadata (flmd.csv) and data dictionary (dd.csv) were added to comply with the File-Level Metadata Reporting Format. Update on 2022-09-09: Updates were made to reporting format specific files (file-level metadata and data dictionary) to correct swapped file names, add additional details on metadata descriptions on both files, add a header_row column to enable parsing, and add version number and date to file names (v2_20220909_flmd.csv and v2_20220909_dd.csv). Update on 2022-12-20: Updates were made to both the data files and reporting format specific files. Conversion issues affecting ER-PLM locations for anion data was resolved for the data files. Additionally, the flmd and dd files were updated to reflect the updated versions of these files. Available data was added up until 2022-03-14. Update on 2023-08-08: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-05-19. The file level metadata and data dictionary files were updated to reflect the additional data added. Update on 2024-03-11: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-09-11. Further, revisions to the data files were made to remove incorrect data points (from 1970 and 2001). The reporting format specific files were updated to reflect the additional data added. Update on 2025-05-15: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until the end of WY2024 (September 30, 2024). International Generic Sample Numbers (IGSNs), when registered, were added to the data files. The reporting format specific files were updated to reflect the additional data added. Update on 2026-09-01: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until the end of WY2025 (September 30, 2025). An anion MDL document was included in this update.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1991

The role of thermal vapor diffusion in the subsurface hydrologic evolution of Mars

The hydrologic response of groundwater to the thermal evolution of the early martian crust is considered. When a temperature gradient is present in a moist porous medium, it gives rise to a vapor-pressure gradient that drives the diffusion of water vapor from regions of high to low temperature. By this process, a geothermal gradient as small as 15 K/km could drive the vertical transport of 1 km of water to the freezing front at the base of the martian crysophere every 10 exp 6-10 exp 7 years, or the equivalent of about 100-1000 km of water over the course of martian geologic history. Models of the thermal history of Mars suggest that this thermally-driven vapor flux may have been as much as 3-5 times greater in the past. The magnitude of this transport suggests that the process of geothermally-induced vapor diffusion may have played a critical role in the initial emplacement of ground ice and the subsequent geomorphic and geochemical evolution of the martian crust.

Clifford, Stephen M.↗

DOE OSTI2022

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗

NASA NTRS1985

Paleohydrologic Implications of Valley Networks of Mars

Valley networks in the heavily cratered terrains of Mars represent an ancient epoch of hydrologic conditions greatly different from those of today. Available crater counts on the valley networks indicate formation during the high flux of impacting bodies charaterizing the early heavy bombardment phase of Martian history. Two populations of valleys are recognized in the equatorial regions of Mars: pristine and degraded. The latter probably formed at the very end of the heavy bombardment phase, extending into the post-heavy bombardment by formation in the intercrater plains. Pristine valleys generally form segments of larger networks with degraded components. This suggests that valley formation was a prolonged process coeval with the heavy bombardment period and extending just beyond that period in martian history. The pristine networks and pristine portions of compound networks on Mars show morphological attributes consistent with an origin by headward growth through spring sapping. On Earth spring sapping occurs where groundwater out-streamflow can be generated by insolation changes associated with orbital parameters or with geothermal effects, such as might be associated with impact or with the volcanic emplacement of the intercrater plains. Thus, it is appropriate to specify the most conservation deviation from modern hydrologic conditions on Mars that could account for the ancient epoch of valley formation.

Baker, V. R.↗

DOE OSTI2020

The Use of Silver Chloride Injection in Remediation of Iodine-129 by In Situ Capture as Silver Iodide at the F-Area Seepage Basin - 20225

The Savannah River Site (SRS) produced tritium, plutonium, and special nuclear materials for national defense, medicine, and the space programs. As part of operations, the F-Area Seepage Basins operated until 1988 for the disposition of deionized acidic wastewater from the F Separations Facility. The wastewater contained dilute nitric acid and low concentrations of non-radioactive metals, and radionuclides, with the major isotopes being Cs-137, Sr-90, U-235, U-238, Pu-239, Tc-99, I-129, and tritium. The seepage basins were closed in 1988 and backfilled and capped by 1991. The groundwater emanating from beneath the closed and capped seepage basins is acidic and contains elevated levels of both chemical and radiological contaminants. Releases from the groundwater plumes sourced from the F-Area Seepage Basins have impacted the water quality of Fourmile Branch, which is a small tributary to the Savannah River, a regional water source. A large pump-and-treat system was constructed in 1997 and operated until 2003 in an attempt to capture the releases to Fourmile Branch. The system in F Area and a similar system in H Area were expensive (∼$1.3 M/month) to operate and produced large quantities of radioactive waste, with concentrations of I- 129 too high to be disposed of at the SRS. In 2004, SRS replaced pump-and-treat with a funnel and gate system that along with operation of a base injection system at the gates reduces the flux of contaminants to the wetlands adjacent to Fourmile Branch. The alkaline solution injected into the aquifer neutralizes the acidic plume and immobilizes many of the cationic constituents. However, base injection is not effective in managing the release of iodine-129, an anionic contaminant. To address iodine-129 SRS and the Savannah River National Laboratory developed an in situ technology that uses ultra-fine ground silver chloride (AgCl) as an injectable capture medium for the sequestration of iodine-129. The AgCl amendment has a very small particle size and is designed to be injected into the contaminated aquifer to capture iodine-129. Dissolved iodine-129 forms a stable and highly insoluble solid (silver iodide) upon contact with AgCl. Laboratory studies, a field scale pilot test (2009), and three deployments (2011, 2015 and 2019) of AgCl have been successfully performed at the F-Area Seepage Basins. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

The Challenge of Achieving Site Closure of a Complex Groundwater Plume - 20531

Over the last three decades, DOE has made substantial progress in nearly every area of nuclear waste cleanup and has closed and transitioned 91 of the 107 sites in the EM complex. Several of the remaining groundwater plumes at Hanford, Los Alamos, Paducah, Oak Ridge, and Savannah River are complex and tend to have unique site or contaminant characteristics which provide an almost impossible challenge to remediate to regulatory levels. As a result, these complex sites will require a large investment to remediate. In order for DOE to achieve regulatory closure, a technically sound case must be developed and presented to regulators and stakeholders that clearly demonstrates that remedial goals have been or will be irreversibly achieved or that there is an alternate acceptable end-state. Five lines of evidence are required: - A conceptual site model that is as complete as possible, with the knowledge that it may continue to evolve with time; - Demonstration that trends in contaminant concentrations, or flux, versus time indicate that remedial goals will be achieved on schedule or within a reasonable time thereafter and that trends that deviate from this path are explainable by the conceptual model; - Demonstration that remobilization of attenuated contaminants is unlikely to occur in the future; - Explicit consideration of remedial contingencies, should deviations in the path toward remedial goals occur; - Development of an effective long-term monitoring plan to detect any conditions that might alter the path toward remedial goals. SRS has made significant progress in the remediation of contamination associated with the F-area and H-Area Seepage Basins, and it is currently developing a strategic closure path that addresses each of these lines of evidence and identifies both basic and applied scientific issues that will need to be addressed in order to obtain regulatory closure. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Hydraulic redistribution by deeply rooted grasses and its ecohydrologic implications in the southern Great Plains of North America

Perennial bioenergy crops with deep (>1 m) rooting systems, such as switchgrass (Panicum virgatum L.), are hypothesized to increase carbon storage in deep soil. Deeply rooted plants may also affect soil hydrology by accessing deep soil water for transpiration, which can affect soil water content and infiltration in deep soil layers, thereby affecting groundwater recharge. Using stable H and O isotope (δ 2 H and δ 18 O) and 3 H values, we studied the soil water conditions at 20–30 cm intervals to depths of 2.4–3.6 m in paired fields of switchgrass and shallow rooted crops at three sites in the southern Great Plains of North America. We found that soil under switchgrass had consistently higher soil water content than nearby soil under shallow-rooted annual crops by a margin of 15%–100%. Soil water content and isotopic depth profiles indicated that hydraulic redistribution of deep soil water by switchgrass roots explained these observed soil water differences. To our knowledge, these are the first observations of hydraulic redistribution in deeply rooted grasses, and complement earlier observations of dynamic soil water fluxes under shallow-rooted grasses. Hydraulic redistribution by switchgrass may be a strategy for drought avoidance, wherein the plant may actively prevent water limitation. Furthermore, this raises the possibility that deeply rooted grasses may be used to passively subsidize soil water to more shallow-rooted species in inter-cropping arrangements.

3H↗

NASA NTRS1989

Microbial mats in playa lakes and other saline habitats: Early Mars analog?

Microbial mats are cohesive benthic microbial communities which inhabit various Terra (Earth-based) environments including the marine littoral and both permanent and ephemeral (playa) saline lakes. Certain geomorphological features of Mars, such as the Margaritifer Sinus, were interpreted as ancient, dried playa lakes, presumably formed before or during the transition to the present Mars climate. Studies of modern Terran examples suggest that microbial mats on early Mars would have had the capacity to survive and propagate under environmental constraints that would have included irregularly fluctuating regimes of water activity and high ultraviolet flux. Assuming that such microbial communities did indeed inhabit early Mars, their detection during the Mars Rover Sample Return (MRSR) mission depends upon the presence of features diagnostic of the prior existence of these communities or their component microbes or, as an aid to choosing suitable landing, local exploration or sampling sites, geomorphological, sedimentological or chemical features characteristic of their playa lake habitats. Examination of modern Terran playas (e.g., the Lake Eyre basin) shows that these features span several orders of magnitude in size. While stromatolites are commonly centimeter-meter scale features, bioherms or fields of individuals may extend to larger scales. Preservation of organic matter (mats and microbes) would be favored in topographic lows such as channels or ponds of high salinity, particularly those receiving silica-rich groundwaters. These areas are likely to be located near former zones of groundwater emergence and/or where flood channels entered the paleo-playa. Fossil playa systems which may aid in assessing the applicability of this particular Mars analog include the Cambrian Observatory Hill Beds of the Officer Basin and the Eocene Wilkins Peak Member of the Green River Formation.

Bauld, John↗

DOE OSTI2024

Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗

NASA NTRS2025

Integrating the Interconnections Between Groundwater and Land Surface Processes Through the Coupled NASA Land Information System and ParFlow Environment

Understanding the interactions between the atmosphere, the land, and the subsurface is fundamental to hydrology and is critical for a better assessment of the impacts of climate change and human management on hydrological systems. However, many land surface models simplify the subsurface hydrology and thereby these interactions. In this study, we couple the land surface model Noah-MP included in the NASA Land Information System (LIS) with the integrated hydrologic model ParFlow (ParFlow-LIS) using the Earth System Modeling Framework (ESMF) and the National United Operational Prediction Capability (NUOPC). This coupling improves the simulation of water and energy cycle processes by adding the three-dimensional variably saturated and heterogeneous flow in the subsurface using sophisticated and nonlinear physics-based equations as well as the advances in satellite remote sensing-based data assimilation of the land surface, thereby benefiting the integrated hydrologic modeling and data assimilation community. We use the High Plains aquifer, located in the central United States, as a testbed to evaluate the coupled ParFlow-LIS system. The new ParFlow-LIS system accounts for the effects of topographically driven flows on the land surface, producing more fine-scale patterns of land surface states and fluxes than standalone LIS. In addition, ParFlow-LIS enables the consideration of the effect of subsurface water storage on evapotranspiration. This is particularly important in areas and times with dry soils, such as during drought conditions or in the presence of a cone of depression due to pumping.

land surface modeling↗

NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗

DOE OSTI2026

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗