Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface hydrology”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stability of Floodplain Subsurface Microbial Communities Through Seasonal Hydrological and Geochemical Cycles

Riparian floodplains represent an interaction zone between the terrestrial subsurface and rivers, where regional groundwater flows, infiltration, and evapotranspiration drive mixing of water and import/export of nutrients and contaminants. These dynamics create seasonally transient redox conditions that drive biogeochemical transformations, which strongly modify groundwater quality. Microbial responses to changing hydrological conditions are perhaps the critical step connecting hydrology to geochemical transformations and groundwater quality, yet are not well understood. We aimed to address this knowledge gap by monitoring seasonal transitions at the U.S. Department of Energy legacy uranium ore processing site in Riverton, WY, through spring-summer-fall hydrological transitions. Our goal was to characterize the microbial community throughout the soil profile, down to the saturated aquifer, and observe its response to wet-dry transitions across a full season and compare to changes in geochemistry and hydrology. Next-generation sequencing was employed to identify biogeochemically-relevant microbial taxa based on the 16S rRNA gene; we found a broad diversity of microbial clades including taxa involved in sulfur and metal cycling, as well as nitrification. These data were paired with measurements of soil moisture, major nutrients and cations, and trace elements. Overall microbial community composition was dependent on soil depth or type, with seasonal effects only observed in the topsoil or subsurface aquifer. This finding indicates that microbial communities in the transiently-reduced center of the soil profile at the Riverton, WY site are remarkably stable, despite moisture and redox inversions. In addition, these communities likely impact the communities in surrounding soil horizons through export of metabolites and solutes as the water table rises and falls throughout the season.

16S rRNA↗
DOE OSTI2025

Model data for a watershed-scale study in the Portage River Basin (OH) examining the effects of subsurface drainage on the hydrologic response of an agricultural watershed.

This study builds on Rathore et al. (2024, WRR) and investigates the role of artificial tile-drainage on various aspects of watershed hydrological response, with a particular focus on peakflow. The model-data for the original modeling-focused paper (Rathore et al., 2024, WRR) is archived at Rathore et al. (2024, ESS-DIVE). Hence, this model-data archive provides scripts that are specific to this study that includes model updates, processing and analysis scripts. For details and models files of original model, readers are referred to Rathore et al. (2024, ESS-DIVE). The key difference between the model configuration in this study and Rathore et al. (2024, WRR) is that the tile drains are applied to the entire domain, to study the impact of tile-drains on different aspects of hydrological response. Additional scenario considering intensified precipitation after a dry period was also simulated. The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2025

Coupled Hydrological-Thermal-Biogeochemical Modeling for Predicting Arctic Carbon Emissions (CH4PACE)

The Arctic, with nearly 1,700 billion metric tons of carbon stored in its soils, is experiencing significant warming trends that raise concerns about the potential release of greenhouse gases, particularly methane. This report investigates the conversion of carbon stored in Arctic soils and permafrost into methane and its subsequent atmospheric release over the coming decades. Recent data indicate a rise in global methane concentrations, highlighting the importance of understanding the mechanisms of methane generation under anaerobic conditions prevalent in thawing permafrost. Our study employs a multi-faceted approach that integrates fieldwork, benchtop experiments, and numerical simulations to explore the subsurface hydrological-thermal-biogeochemical systems influencing methane emissions. We focus on the relationship between microbial activity in decomposing organic matter under oxygen-limited conditions and the transport of methane out of these saturated wetland soils. Key hydrological parameters, including permeability, porosity, and thermal conductivity, were measured on soil samples from Fairbanks, Alaska, and used to parameterize PFLOTRAN, a coupled hydrological-thermal-reactive transport simulator. Our model enables prediction of the evolution of the active layer and its impact on methane production and release. Our findings underscore the critical need for a mechanistic understanding of methane dynamics in permafrost environments to improve predictions of methane emissions in Earth systems models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The hidden influence of terrestrial groundwater on salt marsh function and resilience

Salt marshes are hotspots of nutrient processing and carbon sequestration. So far, studies addressing spatiotemporal variability in and drivers of salt marsh biogeochemical function, carbon storage and resilience have focused on ocean-driven surface hydrologic influences, neglecting effects of terrestrial hydrology through subsurface connections. Here we evaluate drivers of salt marsh redox potential, a proxy for biogeochemical state, through wavelet analyses and information theory using data from seven marshes. The results point to terrestrial groundwater level as a dominant control on redox variability across all sites. Because redox is a key driver of biogeochemical processes, and specifically oxidation of organic matter that sequesters carbon and maintains marsh elevation, these terrestrial influences are critical to understanding marsh function and evolution. Further, the newly identified links between onshore groundwater levels and marsh redox conditions shift the traditional paradigm and suggest that terrestrial hydrology is a primary control on salt marsh carbon sequestration potential and resilience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sapflow and xylem water isotopes from Snodgrass Mountain, East River Watershed, Colorado USA

This dataset includes sapflux and stable water isotopes of soil water and xylem water for aspen, fir and spruce trees along the Snodgrass Mountain transect in the East River Watershed, Colorado USA. The purpose of generating this dataset was to understand: (1) the total flux of water being used by trees in the watershed, (2) separate the component of transpiration that was derived from recent precipitation vs. older water such as winter precipitation of groundwater and (3) understand how total water use and water sources for the trees varies between species and position on a hillslope. The data were collected from May 2019 until October 2022. The sap flux data were collected using ICT SFM1 sensors and are presented in both units of cm h^-1 and as kg h^-1 by multiplying the sap flux by the sapwood area of the tree. All sap flow data has been been corrected using estimates of wounding diameter, water content of wood and sap wood depth. The xylem water isotope data were collected approximately weekly from each of the trees instrumented with sap flux. The stems were collected and the water extracted using classic cryogenic methods. Isotope measurements were done using a Picarro 2140i analyzer. We provide an estimate of the Seasonal Origin Index for each measurements following Allen et al., 2019 (10.5194/hess-23-1199-2019) where values of -1 equate to trees relying on winter precipitation and +1 tree relying on summer precipitation. We also provide the simultaneous sap flux for each isotope measurement when this data was available. This is an update to an earlier data repository with the same name that only included 2019 data. The new dataset was posted in January 2023 and is inclusive of the original 2019 data but now includes 2021 and 2022 data. Please note there was an error in the units of transpiration in the original dataset. It was incorrectly listed as mm h-1 when it should have been kg h-1. This change was made in May 2024.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2020

Predicting Water Cycle Characteristics from Percolation Theory and Observational Data

The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth’s surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET 0 /P, with P being the precipitation and ET 0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko’s model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Low-Power, Flexible Sensor Arrays with Solderless Board-to-Board Connectors for Monitoring Soil Deformation and Temperature

Landslides are a global and frequent natural hazard, affecting many communities and infrastructure networks. Technological solutions are needed for long-term, large-scale condition monitoring of infrastructure earthworks or natural slopes. However, current instruments for slope stability monitoring are often costly, require a complex installation process and/or data processing schemes, or have poor resolution. Wireless sensor networks comprising low-power, low-cost sensors have been shown to be a crucial part of landslide early warning systems. Here, we present the development of a novel sensing approach that uses linear arrays of three-axis accelerometers for monitoring changes in sensor inclination, and thus the surrounding soil’s deformation. By combining these deformation measurements with depth-resolved temperature measurements, we can link our data to subsurface thermal–hydrological regimes where relevant. In this research, we present a configuration of cascaded I2C sensors that (i) have ultra-low power consumption and (ii) enable an adjustable probe length. From an electromechanical perspective, we developed a novel board-to-board connection method that enables narrow, semi-flexible sensor arrays and a streamlined assembly process. The low-cost connection method relies on a specific FR4 printed circuit board design that allows board-to-board press fitting without using electromechanical components or solder connections. The sensor assembly is placed in a thin, semi-flexible tube (inner diameter 6.35 mm) that is filled with an epoxy compound. The resulting sensor probe is connected to an AA-battery-powered data logger with wireless connectivity. We characterize the system’s electromechanical properties and investigate the accuracy of deformation measurements. Our experiments, performed with probes up to 1.8 m long, demonstrate long-term connector stability, as well as probe mechanical flexibility. Furthermore, our accuracy analysis indicates that deformation measurements can be performed with a 0.390 mm resolution and a 95% confidence interval of ±0.73 mm per meter of probe length. This research shows the suitability of low-cost accelerometer arrays for distributed soil stability monitoring. In comparison with emerging low-cost measurements of surface displacement, our approach provides depth-resolved deformation, which can inform about shallow sliding surfaces.

deformation monitoring↗
NASA NTRS2005

An Ice Thickness Study Utilizing Ground Penetrating Radar on the Lower Jamapa Glacier of Citlaltepetl (El Pico de Orizaba), Mexico

Citlalt6petl (Pico de Orizaba) is a dormant stratovolcano located at the eastern end of the trans-Mexican Volcanic Belt at approximately 19 degrees of latitude. It is one of the largest stratovolcanos in the world and at 5,630 meters above sea level, the highest mountain in Mexico and the third highest in North America. Situated on the summit cone and north face of the volcano is a permanent ice cap known as the Jamapa Glacier. Recent and historical studies of Citlaltepetl have been based primarily on volcanic risk assessment, in particular stability assessments of the summit cone. Relatively little work has been directed toward the glacial environment of the mountain, possibly due in part to its high altitude, steep slopes, and general inaccessibility. In addition to this glacier's potential to contribute to a better understanding of climate change, the Jamapa glacier and its environmental, cryologic and geologic setting could also serve as a valuable terrestrial analog to studies of Martian geology, hydrology, and subsurface ice.

Brown, S. B.↗
NASA NTRS2017

Enhanced Groundwater Recharge Rates and Altered Recharge Sensitivity to Climate Variability Through Subsurface Heterogeneity

Our environment is heterogeneous. In hydrological sciences, the heterogeneity of subsurface properties, such as hydraulic conductivities or porosities, exerts an important control on water balance. This notably includes groundwater recharge, which is an important variable for efficient and sustainable groundwater resources management. Current large-scale hydrological models do not adequately consider this subsurface heterogeneity. Here we show that regions with strong subsurface heterogeneity have enhanced present and future recharge rates due to a different sensitivity of recharge to climate variability compared with regions with homogeneous subsurface properties. Our study domain comprises the carbonate rock regions of Europe, Northern Africa, and the Middle East, which cover 25 of the total land area. We compare the simulations of two large-scale hydrological models, one of them accounting for subsurface heterogeneity. Carbonate rock regions strongly exhibit karstification, which is known to produce particularly strong subsurface heterogeneity. Aquifers from these regions contribute up to half of the drinking water supply for some European countries. Our results suggest that water management for these regions cannot rely on most of the presently available projections of groundwater recharge because spatially variable storages and spatial concentration of recharge result in actual recharge rates that are up to four times larger for present conditions and changes up to five times larger for potential future conditions than previously estimated. These differences in recharge rates for strongly heterogeneous regions suggest a need for groundwater management strategies that are adapted to the fast transit of water from the surface to the aquifers.

water resources↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗
DOE OSTI2022

Subsurface redox potential and water level at the Elkhorn Slough NERR

This resource contains various hydrological, and geochemical data from Elkhorn Slough National Estuarine Research Reserve from the years 2020 and 2021. These data have been used to assess how continuous measurements of environmental variables. The data can be used to understand processes at timescales over which biochemical transformations can happen. Especially, the data were used to explain the local subsurface hydrology, and its implication, in an experimental transect in a coasta estuary. Water level and water temperature were measured in-situ with Solinst pressure transducer loggers (Ontario, Canada). Redox potential was collected using in-situ, redox sensors (Paleoterra, The Netherlands) connected to CR1000X Campbell data loggers (Logan, Utah). Meteorological data was gathered from the Elkhorn Slough meteorological station.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Real-time and post-hoc compression for data from Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) is an emerging sensing technology that records the strain-rate along fiber optic cables at high spatial and temporal resolution. This technique is becoming a popular tool in seismology, hydrology, and other subsurface monitoring applications. However, due to the large coverage (10’s of km) and high density of measurements (1m spacing at 100’s of Hz), a DAS installation could produce terabytes of data records per day. Because many DAS instruments are deployed in remote locations, this large data size poses significant challenges to its transfer and storage. In this paper, we explore lossless compression methods to reduce the storage requirement in both real-time and post-hoc scenarios. Here we propose a two-stage compression method to improve the compression ratio and compression speed. This two-stage compression method could reduce the storage requirement by 40%, which is 20% more than other lossless methods, such as ZSTD. We demonstrate that the compression method could complete its operation well before the DAS instrument needs to output the next file, making it suitable for real-time DAS acquisition. We also implement a parallel compression method for a post-hoc scenario and demonstrate that our method could effectively utilize a parallel computer. With 256 CPU cores, our parallel compression method achieves the speed of 26GB/second.

58 GEOSCIENCES↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

WHONDRS River Corridor Surface Water Metabolites and Geochemistry from Global Sites

This dataset supports a broader study examining the character of organic matter that may be delivered to subsurface sediments via hydrologic exchange. To implement the global survey, free stream sampling kits were provided to interested volunteers throughout the world. Samples were collected with minimal constraints in terms of location, but following strict protocols, and shipped for metabolomic analysis via Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS). In addition, basic geochemistry analyses (e.g., dissolved organic matter concentration) were conducted, standardized photos of each field system were taken, and extensive metadata were captured. Sampling began in 2018 and is ongoing as of 2025. This dataset is comprised of one folders of field photos, one folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data, and one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) international generic sample number (IGSN) mapping file; (6) field protocol; and (7) a subfolder with sample data. The sample data subfolder contains (1) surface water dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) methods codes; (3) surface water FTICR methods; and (4) a subfolder of 12 Tesla (12T) FTICR-MS data. This folder contains three subfolders, one containing the.xml files, one containing the CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, or .png. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.

Biogeochemistry↗