Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 217 records · Page 12

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

Nevada Higher Education Benefits from US Department of Energy (DOE) Environmental Management (EM) Nevada Program's Transfer of Geologic Samples - 20442

The Nevada National Security Site (NNSS), formerly the Nevada Test Site, was the location of 100 historic atmospheric and 828 historic underground nuclear tests from 1951 to 1992. Related to this historic nuclear testing, geologic and hydrologic studies of the site were conducted utilizing the skills and the expertise of the U.S. Geological Survey (USGS) and national laboratories (Los Alamos National Laboratory, Lawrence Livermore National Laboratory), who were principal in leading the development of weapons and were responsible for specific underground testing programs. The subsurface samples and technical/scientific data associated with NNSS geologic studies (past and present) are preserved and stored at the USGS Mercury Core Library and Data Center located at the NNSS. Currently, the facility stores over 2,000,000 linear feet of cores and cuttings from more than 2,600 drill holes that can be accessed for study. Most of the samples and historic geologic work was focused directly on subsurface geologic settings that relate directly to historic underground nuclear testing. Since 1992, the United States has observed a unilateral moratorium on full-scale nuclear testing, and the U.S. Department of Energy (DOE) Environmental Management (EM) Nevada Program is now responsible for hydrogeologic characterization of the potential impacts to the natural groundwater systems that may have resulted from historic underground nuclear tests. As part of this effort, the EM Nevada Program Underground Test Area (UGTA) Activity has drilled and completed over 58 deep (2,000 - 7,000 ft.) characterization wells, totaling in excess of 170,000 linear feet of cuttings and core samples. Cutting samples were collected as triplicate samples for each respective depth interval, to account for potential later nondestructive/destructive analysis and to preserve samples for regulatory purposes. Recently, it was recognized by the EM Nevada Program that opportunities may exist to reduce the cost and floor space required for the storage of geologic samples at the Mercury Core Library without impacting the integrity and representative nature of the samples necessary for project execution. An initiative was sponsored by EM Nevada Program to evaluate several options: 1) disposal of a portion of cuttings and cores in a land fill setting; 2) reduction in sample volume through skeletonizing core and cuttings from wells to a representative but much smaller number of samples; and 3) solicit potential interest in the academic community where other geoscientists could freely access the samples for studies. The Nevada state university system through the Nevada Bureau of Mines and Geology responded positively to the opportunity to receive these NNSS samples from the EM Nevada Program. In June 2019, approximately 17,000 geologic samples, representing greater than 170,000 linear feet of drilling and weighing over 20,000 pounds, were shipped from the NNSS to the Great Basin Science Sample and Records Library located in Reno, Nevada. The benefits resulting from this transfer were realized by the EM Nevada Program in terms of cost and space savings for geologic sample storage at the NNSS. The transfer did not impact the EM Nevada Program mission as access to representative geologic samples for regulatory and scientific purposes was preserved. Moreover, the Nevada Bureau of Mines and Geology acquired a significant resource of subsurface geologic samples and supporting technical data to support academic and scientific studies in a complex volcanic setting in southwest Nevada. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Surface Atmosphere Integrated Field Laboratory (SAIL) Campaign

The science of mountainous hydrology spans the atmosphere through the bedrock and inherently crosses physical and disciplinary boundaries: land-atmosphere interactions in complex terrain enhance clouds and precipitation, while watersheds retain and release water over a large range of spatial and temporal scales. Limited observations in complex terrain challenge efforts to improve predictive models of the hydrology in the face of rapid changes. The Upper Colorado River exemplifies these challenges, especially with ongoing mismatches between precipitation, snowpack, and discharge. Consequently, the U.S. Department of Energy’s (DOE) Atmospheric Radiation Measurement (ARM) user facility has deployed an observatory to the East River Watershed near Crested Butte, Colorado between September 2021 and June 2023 to measure the main atmospheric drivers of water resources, including precipitation, clouds, winds, aerosols, radiation, temperature and humidity. This effort, called the Surface Atmosphere Integrated Field Laboratory (SAIL), is also working in tandem with DOE-sponsored surface and subsurface hydrologists and other federal, state, and local partners. SAIL data can be benchmarks for model development by producing a wide range of observational information on precipitation and its associated processes, including those processes that impact snowpack sublimation and redistribution, aerosol direct radiative effects in the atmosphere and in the snowpack, aerosol impacts on clouds and precipitation, and processes controlling surface fluxes of energy and mass. Preliminary data from SAIL’s first year showcase the rich information content in SAIL’s many data-streams and support testing hypotheses that will ultimately improve scientific understanding and predictability of Upper Colorado River hydrology in 2023 and beyond.

54 ENVIRONMENTAL SCIENCES↗

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

On the similarity of hillslope hydrologic function: a clustering approach based on groundwater changes

Abstract. Hillslope similarity is an active topic in hydrology because of its importance in improving our understanding of hydrologic processes and enabling comparisons and paired studies. In this study, we propose a holistic bottom-up hillslope clustering based on a region's integrative hydrodynamic response quantified by the seasonal changes in groundwater levels ΔP. The main advantage of the ΔP clustering is its ability to capture recharge and discharge processes. We test the performance of the ΔP clustering by comparing it to seven other common hillslope clustering approaches. These include clustering approaches based on the aridity index, topographic wetness index, elevation, land cover, and machine-learning that jointly integrate multiple data. We assess the ability of these clustering approaches to identify and categorize hillslopes with similar static characteristics, hydroclimate, land surface processes, and subsurface dynamics in a mountainous watershed – the East River – located in the headwaters of the Upper Colorado River Basin. The ΔP clustering performs very well in identifying hillslopes with six out of the nine characteristics studied. The variability among clusters as quantified by the coefficient of variation (0.2) is less in the ΔP and the machine learning approaches than in the others (> 0.3 for TWI, elevation, and land cover). We further demonstrate the robustness of the ΔP clustering by testing its ability to predict hillslope responses to wet and dry hydrologic conditions, of which it performs well when based on average conditions.

54 ENVIRONMENTAL SCIENCES↗

Brief communication: Decadal changes in topography, surface water and subsurface structure across an Arctic coastal tundra site

Abstract. In ice-rich polygonal tundra, spatiotemporal heterogeneity in ground-ice melt reshapes topography, impacting infrastructure, water and carbon cycles. This study evaluates changes in topography and subsurface structure at a coastal Arctic site by comparing data from two surveys conducted a decade apart. Each survey includes electrical resistivity tomography, active layer thickness, photogrammetry, and topographic data. Results reveal subsidence and decrease in permafrost table elevation with varying intensity and spatial distribution across polygons, alongside diverse thermal-hydrological responses, such as thermokarst pool formation in high-centered-polygons and more even subsidence in flat-centered-polygons. The study also underscores the value and limitations of sporadic surveys.

Bachman, Jonathan A↗

Propyne: Determination of Physical Properties and Unit Cell Parameters under Titan-Relevant Conditions

With its large size, dense atmosphere, methane-based hydrological-like cycle, and diverse surface features, the Saturnian moon Titan is one of the most unique of the outer Solar System satellites. Study of the photochemically produced molecules in Titan’s atmosphere is critical in order to understand the mechanics of the atmosphere and, by extension, the interactions between atmosphere, surface, and subsurface water ocean. One example is propyne vapor, a photochemically produced species in Titan’s upper atmosphere expected to condense in Titan’s stratosphere at lower altitudes. Propyne may also be a trace species in Titan’s stratospheric co-condensed ice clouds detected by the Cassini Composite InfraRed Spectrometer. Bulk structural characterization of propyne ice is currently incomplete and is lacking in published laboratory Raman spectra and X-ray diffraction data. Here, we present a laboratory characterization of propyne ice, including the first published X-ray diffraction and Raman spectroscopy results for propyne ice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (Final Scientific/Technical Report)

The ENIGMA science focus area (SFA) is a multi-disciplinary, multi-institutional research effort focused on addressing foundational knowledge gaps in environmental microbial communities by studying groundwater and sediment microbiomes in the shallow subsurface at the contaminated Oak Ridge Field Research Site (FRC). We seek to discover and characterize the reciprocal interactions between the microbial communities and the geochemical and geophysical parameters of the shallow subsurface within the contamination plume. The goal of the Fields lab is to help manage and conduct laboratory experiments informed by field observations using ecological and physiological approaches to study microbial populations in situ as well as ex situ (in the laboratory). In particular, we aim to study and characterize the impact of hydrological constraints on free-living and biofilm biomass and activity with increasing spatial and temporal resolution under static and flow conditions. We also aim to help characterize novel microbial groups that are present and active under relevant field conditions, including the development of molecular techniques for in situ detection as well as metabolic interactions that underlie pertinent physiology and ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗

Machine learning-enabled model-data integration for predicting subsurface water storage

Subsurface water storage (SWS) is a key variable of the climate system and a storage component for precipitation and radiation anomalies, inducing persistence in the climate system. It plays a critical role in climate-change projections and can mitigate the impacts of climate change on ecosystems. However, because of the difficult accessibility of the underground, hydrologic properties and dynamics of SWS are poorly known. Direct observations of SWS are limited, and accurate incorporation of SWS dynamics into Earth system land models remains challenging. We propose a machine learning-enabled model-data integration framework to improve the SWS prediction at local to conus scales in a changing climate by leveraging all the available observation and simulation resources, as well as to inform the model development and guide the observation collection. The accurate prediction will enable an optimal decision of water management and land use and improve the ecosystem's resilience to the climate change.

Lu, Dan↗

Improving collaborations between empiricists and modelers to advance grassland community dynamics in ecosystem models

Climate change, increasing atmospheric CO 2 , and land use change have altered biogeochemical and hydrologic cycles world-wide, with grassland systems being particularly vulnerable to resulting vegetation shifts (Komatsu et al., 2019). Therefore, incorporating plant community dynamics into ecosystem models is critical for accurate forecasting of ecosystem responses to global change (Levine, 2016). Process-based ecosystem models, which simulate the biogeochemical transfers of mass and energy among biota, the subsurface, and atmosphere, require representation of dynamic composition of organisms within ecosystems. For example, these models simulate leaf and plant-level characteristics, such as electron transport rate and allometry of carbon (C) allocation, to predict how net primary productivity and other ecosystem processes respond to abiotic drivers. These models are particularly useful in scaling from organismal to ecosystem levels but are still underdeveloped in their ability to capture community change, especially in grassland ecosystems. To represent compositional changes, these models must simulate competition, mortality, establishment, and reproduction of plant populations within communities. Yet, current ecosystem modeling approaches to forecast plant community change have derived from studies of forested systems and are either too coarse to capture fine-scale community dynamics (e.g. dynamic global vegetation models (DGVMs)) or too complex to be used at large spatial scales (e.g. forest gap models).

54 ENVIRONMENTAL SCIENCES↗

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

60 APPLIED LIFE SCIENCES↗

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES↗

Model output from simulations of manganese-carbon interactions in temperate forest soil profiles

This archive contains model output, code, and scripts for simulations of coupled manganese-carbon cycling in temperate forest soil profiles. These model results were generated as part of a study investigating how manganese availability influences soil organic carbon stocks and demonstrating a new model framework for coupling carbon and manganese cycling. The simulations were in support of a manuscript: "Modeling interactive effects of manganese bioavailability, nitrogen deposition, and warming on soil carbon storage." The study addresses the research questions: How does Mn bioavailability, as driven by subsurface mineral properties, pH, and redox status, affect temperate forest soil organic carbon and litter carbon stocks?How is the relationship between Mn bioavailability and carbon cycling affected by changes in temperature and nitrogen deposition?"Model simulations were conducted in a reactive transport modeling framework using PFLOTRAN coupled to python. Multiple model simulations testing different Mn-bearing mineral solubilities, hydrological patterns, nitrogen deposition rates, and temperatures are included. Soil properties, including total and exchangeable Mn concentrations, are based on values reported for the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), a temperate forested watershed in central Pennsylvania, U.S.A where Mn cycling through vegetation has been documented.File formats include netCDF (.nc), python script (.py), shell script (.sh), plain text PFLOTRAN input file (.in), and plain text PFLOTRAN database file (.dat), and gzipped tar archive (tar.gz).

54 ENVIRONMENTAL SCIENCES↗

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗

Decovalex-2019 (Executive Summary)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. The overall aim of DECOVALEX-2019 was to increase the understanding of various thermo-hydro-mechanical-chemical processes of importance for radionuclide release and transport from a repository to the biosphere and how they can be described and modelled using mathematical models. The scientific and technical objectives are: to increase the basic understanding of T-H-M-C coupled processes in fractured rocks (crystalline, sedimentary, argillaceous) and buffer materials; to investigate the predictive capabilities of different codes to field experiments and to perform verification of codes; to exchange experimental data, and improve the understanding of the constitutive behavior of crystalline and argillaceous rock masses and buffer materials; and to perform THMC calculations in a performance/safety assessment context.

58 GEOSCIENCES↗

DECOVALEX-2019 (Task E Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. This document is the final report of Task E which was proposed and coordinated by Andra, the National Radioactive Waste Management Agency in France, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. The purpose of Task E of the DECOVALEX-2019 project is to investigate upscaling THM modelling from small-scale experiments (some cubic meters) to full-scale experiments (some ten cubic meters) and finally to the scale of the waste repository (cubic kilometers). To achieve this aim, the data of two in-situ heating experiments performed by Andra (the French National Radioactive Waste Management Agency) in the Meuse/Haute-Marne Underground Research Laboratory (MHM URL) have formed the basis for the understanding of the THM behavior of the COx at different scales. The first experiment provided the reference values of the THM parameters by means of a calibration exercise and they were used for a blind prediction and an interpretative analysis of the second one.

58 GEOSCIENCES↗