Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aquifer”

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 667 records · Page 37

Shallow Geothermal Potential of the Snake River Plain

The Snake River Plain (SRP) terrestrial heat flow and subsurface thermal regime are not well understood but are important for assessing the local geothermal resource potential, both for conventional and for Enhanced Geothermal Systems (EGS) development in the region. Resource evaluation for the SRP is complicated by the disparate data density, along with the known lateral advection of heat in the Eastern Snake River Plain Aquifer and vertical heat transport by fluids in the bounding faults, primarily in the southwestern section. Fortunately, recent studies, e.g., the Snake River Plain Play Fairway Analysis, the Idaho FORGE site, and site-specific investigations, which included drilling within the Camas Prairie and on the Mountain Home Air Force Base, near Twin Falls, and in the Eastern Snake River Plain as part of the HOTSPOT Project, add both additional drilling and geophysical data. The SMU Geothermal Laboratory has conducted detailed studies of SRP tectonics and heat flow since the 1970's and used this knowledge as part of the EGS geothermal potential estimation for the conterminous United States in 2006 and again in 2011, calculating geothermal potential from 3.5 km to 10 km depth. Recent temperature modeling refined the calculation methodology to estimate shallow (1 km to 4 km) resource potential using an improved thermal conductivity model and incorporation of shallow groundwater flow. By incorporating the new SRP geology, geophysics, and 206 thermal data sites into the SMU thermal modeling methodology, this project updates the resource estimate for the SRP, and generates new temperature-at-depth maps for the shallow subsurface (1 km to 4 km). The project results highlight the EGS potential resource areas (=150°C) and areas with more exploration risks based on minimal and/or low-quality data. The newest temperature modeling results suggest EGS potential is near five times greater in the SRP than previously estimated.

EGS↗

Impact of sub-core scale heterogeneity on CO2/brine multiphase flow for geological carbon storage in the Minnelusa sandstone

CO2 geological storage in deep saline aquifers is a mitigation option for CO2 emissions due to its large storage capacity and immediate accessibility. Accurately determining the CO2-brine relative permeability curves is key to the evaluation of CO2 injectivity and sweep efficiency in reservoir simulation as well as the CO2 injection in the field. This study highlights the remarkable effects of sub-core scale heterogeneity on the CO2-brine multiphase flow properties of the Minnelusa Sandstone in Wyoming. Two unsteady state CO2-brine drainage experiments were performed on the two samples. The first sample exhibits slanted laminated structure, while the second one represents a more homogeneous sandstone system. The CO2 saturation distributions during drainage reveals that the main variation in multiphase flow properties of two core samples were attributed to the porosity distribution that leads to the capillary pressure heterogeneity. Assisted history matching was used to obtain the respective relative permeability curves, which suggests heterogeneity-dependent behavior. In addition, sensitive and uncertainty analyses indicate that physical and petro-physical properties of low-porosity and low-permeability bedding layers exert marked effects on CO2 front breakthrough time and brine production. The results presented in this study help to gain insight into the CO2-brine multiphase flow properties in heterogeneous sandstones and can pave the way for the upscaling of CO2 migration and field-scale simulation accurately. This work is funded under the Department of Energy CarbonSAFE program (awards DE-FE0031624 DE0031891).

Kou, Zuhao↗

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Response to Community-Scale Geothermal Energy Production

Geothermal heat, as renewable energy, shows great advantage with respect to its environmental impact due to its significantly lower CO2 emissions than conventional fossil fuel. Open and closed-loop geothermal heat pumps, which utilize shallow geothermal systems, are an efficient technology for cooling and heating buildings, especially in urban areas. Integrated use of geothermal energy technologies for district heating, cooling, and thermal energy storage can be applied to optimize the subsurface for communities to provide them with multiple sustainable energy and community resilience benefits. The utilization of the subsurface resources may lead to a variation in the underground environment, which might further impact local environmental conditions. However, very few simulators can handle such a highly complex set of coupled computations on a regional or city scale. We have developed high-performance computing (HPC) based hydrothermal finite element (FE) simulator that can simulate the subsurface and its hydrothermal conditions at a scale of tens of km. The HPC simulator enables us to investigate the subsurface thermal and hydrologic response to the built underground environment (such as basements and subways) at the community scale. In this study, a coupled hydrothermal simulator is developed based on the open-source finite element library deal.II. The HPC simulator was validated by comparing the results of a benchmark case study against COMSOL Multiphysics, in which Aquifer Thermal Energy Storage (ATES) is modeled and a process of heat injection into ATES is simulated. The use of an energy pile system at the Treasure Island redevelopment site (San Francisco, CA, USA) was selected as a case study to demonstrate the HPC capability of the developed simulator. The simulator is capable of modeling multiple city-scale geothermal scenarios in a reasonable amount of time.

Sun, Xiang↗

Accurate and Timely Forecasts of Geologic Carbon Storage using Machine Learning Methods

Carbon capture and storage is one strategy to reduce greenhouse gas emissions. One approach to storing the captured CO2 is to inject it into deep saline aquifers. However, dynamics of the injected CO2 plume is uncertain and the potential for leakage back to the atmosphere must be assessed. Thus, accurate and timely forecasts of CO2 storage via real-time measurements integration becomes very crucial. This study proposes a learning-based, inverse-free prediction method that can accurately and rapidly forecast CO2 movement and distribution with uncertainty quantification based on limited simulation and observation data. The machine learning techniques include dimension reduction, multivariate data analysis, and Bayesian learning. The outcome is expected to provide CO2 storage site operators with an effective tool for real-time decision making.

Lu, Dan↗

Geological Thermal Energy Storage Using Solar Thermal and Carnot Batteries: Techno-Economic Analysis: Preprint

Energy storage is increasingly necessary as Variable Renewable Energy (VRE) technologies are deployed. Seasonal energy storage can shift energy generation from the summer to the winter, but these technologies must have extremely large energy capacities and low costs. Geological Thermal energy storage (GeoTES) is proposed as a solution for long-term energy storage. Excess thermal energy can be stored in permeable reservoirs such as aquifers and depleted hydrocarbon reservoirs for several months. In this article, we describe a techno-economic model that has been developed to evaluate GeoTES systems. The models are developed by combining the output of specialist models which enables the performance and cost of both the subsurface and surface systems to be captured. Off-design models are developed so that the performance can be evaluated at each hour of the year. GeoTES can be charged with two different energy sources: (1) concentrating solar thermal and (2) renewable electricity using heat pumps (henceforth known as a "Carnot Battery"). The stored thermal energy can be used to generate electricity and - uniquely - also directly produce heat that can be used by industrial processes. Furthermore, Carnot-Battery-GeoTES can also be used to form a cold storage reservoir. Preliminary results that quantify the technical and economic performance of these two GeoTES systems are presented.

carnot battery↗

SHASTA-HELP: Hydrogen Estimator for Logistical Planning Tool

Underground reservoirs offer a promising option for medium- and long-term hydrogen storage that could potentially be safer and more cost-effective than surface alternatives. Currently, there are three underground hydrogen storage (UHS) facilities in the United States, which operate solely in salt caverns and store hydrogen primarily for industrial purposes. Salt formations favorable for geologic storage are geographically limited in the U.S. If a hydrogen economy is to develop, expansion of UHS to other geologic formations such as depleted gas fields and saline aquifers will be necessary to increase the availability of hydrogen storage. While research into the feasibility of storing hydrogen in these formations is ongoing, tools and methodologies are needed for characterizing the hydrogen storage potential in these reservoirs. We introduce the Hydrogen Estimator for Logistical Planning (HELP) tool from the Subsurface Hydrogen Assessment, Storage, and Technology Acceleration (SHASTA) project (https://edx.netl.doe.gov/shasta). SHASTA-HELP is a suite of web-based tools that have numerous capabilities for UHS project planning.

Haeri, Foad↗

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↗

Geological Thermal Energy Storage Using Solar Thermal and Carnot Batteries: Techno-Economic Analysis

Energy storage is increasingly necessary as variable renewable energy technologies are deployed. Seasonal energy storage can shift energy generation from the summer to the winter, but these technologies must have extremely large energy capacities and low costs. Geological thermal energy storage (GeoTES) is proposed as a solution for long-term energy storage. Excess thermal energy can be stored in permeable reservoirs such as aquifers and depleted hydrocarbon reservoirs for several months. In this article, we describe a techno-economic model that has been developed to evaluate GeoTES systems. The models are developed by combining the output of specialist models, which enables the performance and cost of both the subsurface and surface systems to be captured. Off-design models are developed so that the performance can be evaluated at each hour of the year. GeoTES can be charged with two different energy sources: (1) concentrating solar thermal and (2) renewable electricity using heat pumps (henceforth known as a "Carnot Battery"). The stored thermal energy can be used to generate electricity and, uniquely, also directly produce heat that can be used by industrial processes. Furthermore, Carnot Battery GeoTES can also be used to form a cold storage reservoir. Preliminary results that quantify the technical and economic performance of these two GeoTES systems are presented.

carnot battery↗

A Comparison of Groundwater Sampling Technologies, Including Passive Diffusion Sampling, for Radionuclide Contamination - 20205

Using traditional high-flow purge methods for long-term water quality monitoring of deep groundwater wells can be expensive, affect contaminant migration, and produce excessive volumes of discharge water that can be difficult to manage. The use of low-flow pumping methods and depth discrete bailers (DDBs) can reduce the cost of sampling deep groundwater wells. In general, using different pumping methods to obtain reproducible and representative groundwater can be challenging. Passive diffusion samplers (PDSs) have successfully been used in long-term monitoring for volatile organic compounds, major ions, and trace-elements, but application has been limited for stable and radioactive isotopes. Beginning in 2018, the United States Geological Survey (USGS) conducted three sampling events to test the ability of PDSs and DDBs to obtain reproducible and representative groundwater samples. The sampled well is completed in a regional, permeable carbonate aquifer and is one in a set of wells that have historically been used for tritium tracer testing. All samples were obtained at 180 m (590 ft) below land surface (bls) in a 13.97 cm (5.5 inch) uncased well that has a total depth of 202 m (662 ft) bls. The first sampling event deployed a regenerated cellulose dialysis membrane (RCDM) PDS for 14 days with deionized water as the blank. The second sampling event deployed a RCDM for 27 days also with deionized water as the blank. The third sampling event deployed a Dual Membrane (DM) PDS for 65 days using a blank of tritium-dead carbonate water. The DM PDS was used to assess the effect of longer term deployment on tritium concentrations and address whether or not the PDSs reached equilibrium with ambient groundwater. The day after each of the three passive samplers were retrieved a DDB was used to obtain discrete non-integrated groundwater samples. For each DDB sampling day, the bailer was lowered into the well 10 consecutive times to determine if the water chemistry changed from the first to the last bailed sample. Quality assurance samples including blanks and duplicates were obtained during all three sampling events. All blank waters had tritium concentrations less than 21±33 pCi/L. Major ion (e.g. calcium, chloride, sodium, and sulfate) results were compared between all samples obtained with RCDM and DDB. Major ion concentrations showed a coefficient of variation of less than 6% between all RCDM and DDB samples; however, the coefficient of variation between the different deployment times and the two different methods for trace-element concentrations was much larger, particularly for manganese (82%), lead (41%), and zinc (33%). Stable isotope values were compared between the RCDM and DDB samples. The DDB sample results all fell within analytical uncertainty and were considered representative of the formation groundwater. The stable isotope values from the RCDM samples indicated that a longer deployment time was necessary to gain equilibrium and to obtain representative groundwater samples. Tritium results from groundwater samples obtained from the RCDM, DM, and DDB indicate that groundwaters obtained with PDSs produced tritium concentrations 2 to 3 times higher (between 1,426 and 3,060±87 pCi/L) than groundwaters obtained with a DDB (479 to 1,219±51 pCi/L). The longer the passive diffusion samplers were deployed, the higher the tritium concentration, suggesting that equilibrium with tritium was not reached within a 27-day deployment. DDB samples showed tritium results declining from the first to the last bailed sample for all three sampling events. This research suggests that tritium results from groundwater samples obtained from PDSs are more reproducible than samples obtained from DDBs. Also, PDSs likely do not accumulate isotopes of water but rather equilibrate with the ambient groundwater. On the other hand, both PDSs and DDBs were able to provide representative groundwater samples for major ions and have the potential to produce representative groundwater samples for stable isotopes of water (δD and δ{sup 18}O). Utilizing the most effective method for sampling groundwater is important for long-term groundwater monitoring. Quantifying the ability of passive diffusion groundwater samplers to obtain representative radioactive chemistry samples will help determine if these methods can be a useful alternative to pumping and purging wells for long-term monitoring. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SUBTASK 1.6 – BASIN ELECTRIC CARBON STORAGE RESEARCH PROJECT: NOVEL MONITORING TECHNIQUES

The Energy & Environmental Research Center (EERC) conducted baseline activities associated with an applied research project at Basin Electric Power Cooperative’s (Basin’s) carbon capture and storage (CCS) site in Beulah, North Dakota, to establish novel carbon storage-monitoring techniques as commercial methods under Cooperative Agreement No. DE-FE0024233, Subtask 1.6. The following report summarizes the baseline activities performed and briefly describes the subsequent (operational monitoring) activities that have been proposed to the U.S. Department of Energy (DOE) as part of the overall project to develop and demonstrate novel monitoring techniques at North America’s largest permitted CCS operation. Dakota Gasification Company (DGC), a wholly owned subsidiary of Basin, owns and operates the Great Plains Synfuels Plant (GPSP) approximately 5 miles northwest of the town of Beulah, North Dakota (Figure 1). In 2023, DGC received approval from the North Dakota Industrial Commission (NDIC) to develop a storage facility on-site for injecting a stream of carbon dioxide (CO2) captured from GPSP. DGC will transport the captured CO2 stream with approximately 6.8 miles of transmission lines that extend north of GPSP and inject >1 million tonnes (MMt) of CO2 annually (>1 MMt/yr) over a 12-year period with up to six underground injection control (UIC) Class VI-compliant injection wells completed in the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer underlying GPSP. The Broom Creek Formation lies approximately 5900 feet (ft) below ground surface (bgs) at GPSP. The commercial scale (i.e., >1 MMt/yr) of DGC’s permitted carbon storage project is ideal for developing and testing the novel monitoring techniques included within Subtask 1.6. The goals of this project are to demonstrate 1) the cost-effectiveness of novel monitoring technologies included as part of this research, 2) technology capability for tracking the CO2 plume and/or associated pressure response in the subsurface and monitoring out-of-zone migration, and 3) compliance with UIC Class VI program requirements. The research activities proposed for the overall project include 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; 5) advanced wellbore-monitoring methods; 6) deployment of an AIM monitoring network; 7) EM monitoring of CO2 with real-time data processing; 8) continued seasonal drone-based surveillance studies; 9) seismic monitoring with passive and active surveys; and 10) wellbore monitoring with nuclear magnetic resonance (NMR) for near-surface characterization. Completion of Activities 1.0–5.0 (baseline activities) are described in this report. Upon authorization of funding by DOE, the EERC will initiate Activities 6.0– 10.0 (operational monitoring activities). Current state-of-the-art (SOA) carbon storage-monitoring techniques require countless labor hours dedicated to the acquisition of data. Once data are gathered, these SOA techniques often rely on commercial facilities to process raw data from the field. However, it is anticipated that next-generation monitoring techniques, such as those being demonstrated, will lower acquisition footprints, be less operationally intensive, and improve data acquisition efficiencies. These new techniques are more conducive to the application of machine learning, artificial intelligence, and automation, thus providing a pathway for integration into active control systems, informing site operability, and improving the integration of data for future CCS projects across the United States. Additionally, reclaimed and active mining lands are present within the project site, creating a unique opportunity to demonstrate the effectiveness of remote sensing and surface-based geophysics monitoring techniques at similar project sites that may include disturbed, unconsolidated, or actively excavated near-surface environments. The efforts included in the overall project will produce necessary designs, learnings, and data acquired during the baseline and operational monitoring periods that are necessary for time-lapse demonstration and validation of the described monitoring techniques. In addition, it is anticipated that the monitoring technologies included in this study will be compliant with UIC Class VI requirements to enable the potential for implementation at other CCS sites across the United States.

42 ENGINEERING↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa–Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat–Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

54 ENVIRONMENTAL SCIENCES↗

Groundwater level elevation and temperature data, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, CO

This data package includes a time series of water level and temperature measurements from October 2018 to December 2021 in groundwater and surface water from the Slate River floodplain outside Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The data were recorded by pressure transducers were installed in four types of piezometers: 1) a network of gravel bed ("GB") drive-point piezometers with a 6-inch screen interval installed all at ~330 cm below ground surface. 2) a network of piezometers screened across the water table, used to measure evapotranspiration ("ET") using the White method. Each of these piezometers is screened along almost its entire length.3) a suite nested piezometers used the measure the vertical hydraulic gradient ("VHG") across the fines-cobble interface. Each of these piezometers uses a 6-inch screen length.4) a group of piezometers scattered across the boundaries of the floodplain, used to monitor boundary conditions ("BC") flowing into and out of the floodplain. With the exception of "SR-BD-WT", each of these piezometers is screened along its entire length. Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 11 "*data.csv" files, one for each piezometer type for each calendar year. Because piezometers have been added over time, not every sensor has data dating back to Oct 2018. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, elevation, screen depth, sediment type and sensor metadata for each piezometer and pressure transducer.

54 ENVIRONMENTAL SCIENCES↗

Estimating irrigation water use from remotely sensed evapotranspiration data: Accuracy and uncertainties at field, water right, and regional scales

Irrigated agriculture is the dominant user of water globally, but most water withdrawals are not monitored or reported. As a result, it is largely unknown when, where, and how much water is used for irrigation. Here, we evaluated the ability of remotely sensed evapotranspiration (ET) data, integrated with other datasets, to calculate irrigation water withdrawals and applications in an intensively irrigated portion of the United States. We compared irrigation calculations based on an ensemble of satellite-driven ET models from OpenET with reported groundwater withdrawals from hundreds of farmer irrigation application records and a statewide flowmeter database at three spatial scales (field, water right group, and management area). At the field scale, we found that ET-based calculations of irrigation agreed best with reported irrigation when the OpenET ensemble mean was aggregated to the growing season timescale (bias = 1.6–4.9%, R 2 = 0.53–0.74), and agreement between calculated and reported irrigation was better for multi-year averages than for individual years. At the water right group scale, linking pumping wells to specific irrigated fields was the primary source of uncertainty. At the management area scale, calculated irrigation exhibited similar temporal patterns as flowmeter data but tended to be positively biased with more interannual variability. Disagreement between calculated and reported irrigation was strongly correlated with annual precipitation, and calculated and reported irrigation agreed more closely after statistically adjusting for annual precipitation. The selection of an ET model was also an important consideration, as variability across ET models was larger than the potential impacts of conservation measures employed in the region. From these results, we suggest key practices for working with ET-based irrigation data that include accurately accounting for changes in soil moisture, deep percolation, and runoff; careful verification of irrigated area and well-field linkages; and conducting application-specific evaluations of uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

Geothermal district energy systems coupled with seasonal underground thermal energy storage: a U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

25 ENERGY STORAGE↗

Short-Term Groundwater Level Fluctuations Drive Subsurface Redox Variability

As global change processes modify the extent and functions of terrestrial-aquatic interfaces, the variability of critical and dynamic transitional zones between wetlands and uplands increases. However, it is still unclear how fluctuating water levels at these dynamic boundaries alter groundwater biogeochemical cycling. Here, we used high-temporal resolution data along gradients from wetlands to uplands and during fluctuating water levels at freshwater coastal areas to capture spatiotemporal patterns of groundwater redox potential (E h ). We observed that topography influences groundwater E h that is higher in uplands than wetlands; however, the high variability within TAI zones challenged the establishment of distinct redox zonation. Declining water levels generally decreased E h , but most locations exhibited significant E h variability, which is associated with rare instances of short-term water level fluctuations introducing oxygen. The E h -oxygen relationship showed distinct hysteresis patterns, reflecting redox poising capacity at higher E h , maintaining more oxidizing states longer than dissolved oxygen presence. Surprisingly, we observed more frequent oxidizing states in transitional areas and wetlands than in uplands. We infer that occasional oxygen entering specific wetland-upland boundaries acts as critical biogeochemical control points. Finally, high-resolution data can capture such rare yet significant biogeochemical instances, supporting redox-informed models and advancing the predictability of climate change feedback.

54 ENVIRONMENTAL SCIENCES↗