Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

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 397 records · Page 22

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

Utah FORGE 4-2492: Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery - 2024 Annual Workshop Presentation

This is a presentation on the Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery Efficiency by The University of Texas at Austin, presented by Mukul M. Sharma. This video slide presentation discusses the following objectives: (1) to place fractures uniformly in a horizontal well to ensure a uniform distribution of flow into the fractures; (2) to maximize the area of the created fracture; (3) to ensure connectivity of the created fracture network; and (4) to ensure fracture size is optimized. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 14, 2024.

15 GEOTHERMAL ENERGY↗

Geothermal Deep Direct Use for Turbine Inlet Cooling in East Texas

The National Renewable Energy Laboratory (NREL), the Southern Methodist University Geothermal Laboratory (SMU), Eastman Chemical (Longview, TX), and TAS (Houston, TX) evaluated the feasibility of using geothermal heat to improve the performance of a natural-gas power plant in East Texas. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of a geologic region known as the Sabine Uplift. The feasibility study focused on determining the potential for accessing a subsurface hot-water geothermal resource within a 10-km radius of the site to provide thermal energy for absorption chillers. Wells within a 20-km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, field porosity and permeability. The lithologies of most interest are the Lower Cretaceous Trinity Group and Upper Jurassic Cotton Valley Group. The deeper Cotton Valley formations are hotter (averaging 117 to 130°C), yet permeability and porosity are low. The shallower Trinity Group contains more variability in permeability and porosity and lower temperatures averaging about 98 to 117°C. The shallower formations are considered despite the lower temperature because of increased ability to produce larger volumes of water and extract enough heat before reinjection. The complete SMU analysis is available in the National Geothermal Data System (NGDS). Tapping such deep geothermal sources for direct heating (as opposed to power generation) is known as geothermal deep direct use (DDU). Geothermal DDU has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for what are typically small-scale, variable-demand projects. This project examines the feasibility of geothermal energy integration in a natural-gas combined cycle power station in East Texas. The DDU resource is tapped to drive absorption chillers (24/7) for production of chilled water at 5-10°C (41-50°F). This chilled water is stored until needed, which allows for continuous operation with a relatively small-capacity geothermal/absorption chiller system. When conditions are favorable, the chilled water is dispatched to cool the air entering the compressor stage of a gas combustion turbine. This process, known as turbine inlet cooling (TIC), boosts power production during periods of high temperature and high-power demand. Such systems can enhance grid reliability and reduce the cost for peak-demand power. A simulation model of the power plant was developed in IPSEpro software and validated against operational data from the plant. This model allowed the team to estimate the additional power that could be produced by applying TIC under different operating and ambient conditions. Absorption chiller performance was estimated from vendor sources to determine the production rate of chilled water from the geothermal resource. Geothermal drilling and development costs were estimated using NREL's GEOPHIRES 2.0. The expected lower drilling costs in this region led to an estimated cost of geothermal heat of about $4/MMBtu (1.4 cents/kWh t ). The estimated cost for the absorption chillers and TIC hardware were obtained from literature sources and project partners. Hourly data were obtained for weather, natural gas and electricity prices, and plant operating state for 2017, which served as a representative year. NREL estimated the capital cost, operating cost, and additional electricity production and revenue for different combinations of geothermal capacity, chiller capacity, and water storage-tank size. The analysis drove toward smaller geothermal and chiller systems to reduce equipment cost. A relatively low-cost water storage tank accumulated the near-continuous chilled water output for later use when TIC was most valued.

15 GEOTHERMAL ENERGY↗

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗

GEOPHIRES Simulations for Deep Direct Use (DDU) Projects

This folder contains the GEOPHIRES codes and input files for running the base case scenarios for the six deep direct-use (DDU) projects. The six DDU projects took place during 2017-2020 and were funded by the U.S. Department of Energy Geothermal Technologies Office. They investigated the potential of geothermal deep direct-use at six locations across the country. The projects were conducted by Cornell University, West Virginia University (WVU), University of Illinois (U of IL), Sandia National Laboratory (SNL), Portland State University (PSU), and National Renewable Energy Laboratory (NREL). Four projects (Cornell, WVU, U of IL, SNL) investigated geothermal for direct heating of a local campus or community, the project by PSU considered seasonal subsurface storage of solar heating, and the NREL project investigated geothermal heating for turbine inlet cooling using absorption chillers. To allow comparison of techno-economic results across the six DDU projects, GEOPHIRES simulations were set up and conducted for each project. The GEOPHIRES code was modified for each project to simulate the local application and incorporate project-specific assumptions and results such as reservoir production temperature or financing conditions. The base case input file is included which simulates the base case conditions assumed by each project team. The levelized cost of heat (LCOH) is calculated and matches the base case LCOH reported by the project teams.

15 GEOTHERMAL ENERGY↗

Model Development for Thermal-Hydrology Simulations of a Full-Scale Heater Experiment in Opalinus Clay

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. Here, this study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

This data package contains input files for TOUGHREACT for a modeling study examining the effects of managed aquifer recharge on agricultural lands on nitrate cycling and transport in the Central Valley of CA near Modesto. The files contain all the geochemical species, reactions, and hydrological parameters for the model. The files are text files used for the TOUGH family of code created by LBNL. To use the files a license is required. https://tough.lbl.gov/licensing-download/toughreact-licensing-download/Accompanying Paper Abstract: Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO3-) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying water all-at-once rather than in increments for a fixed volume of recharge transports higher concentrations of NO3- deeper into the profile, which has the potential to exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as riparian corridors, floodplains, wetlands, and other managed landscapes like dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗

Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets

Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.

02 PETROLEUM↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated, including a single-, double-, and multi-lateral design. The subsurface within the area of interest is characterized by deep, hot (> 200 degrees C) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lie(s) within a target depth of 6 km, where temperatures are estimated to approach 300 degrees C. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot, dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down techno-economic analysis model. The results show that the performance- and cost-optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60 degrees C/km. By injecting 20 degrees C-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth (i.e., 2.2 MWe net plant output) can be achieved, resulting in an LCOE and LCOH of $136/MWhe and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

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

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

54 ENVIRONMENTAL SCIENCES↗

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

Agricultural managed aquifer recharge (AgMAR) is a strategy whereby surface water is used to intentionally flood croplands to recharge underlying aquifers. However, nitrate (NO 3 - ) contamination in agriculturally intensive regions poses a threat to groundwater resources under AgMAR. We use a reactive transport model to understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on NO 3 - leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. We examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and with preferential flow paths characteristic of agricultural fields. Simulations indicate finer textured sediments, alone or embedded within/adjacent to high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results suggest that applying water all-at-once rather than in increments transports higher concentrations of NO 3 - deeper into the profile, which may exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR can be designed to enhance denitrification in the subsurface and reduce NO 3 - leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as floodplains and dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗

Exploring Plasma Heating in the Current Sheet Region in a Three-dimensional Coronal Mass Ejection Simulation

We simulate a coronal mass ejection using a three-dimensional magnetohydrodynamic code that includes coronal heating, thermal conduction, and radiative cooling in the energy equation. The magnetic flux distribution at 1 Rs is produced by a localized subsurface dipole superimposed on a global dipole field, mimicking the presence of an active region within the global corona. Transverse electric fields are applied near the polarity inversion line to introduce a transverse magnetic field, followed by the imposition of a converging flow to form and destabilize a flux rope, producing an eruption. We examine the quantities responsible for plasma heating and cooling during the eruption, including thermal conduction, radiation, adiabatic effects, coronal heating, and ohmic heating. We find that ohmic heating is an important contributor to hot temperatures in the current sheet region early in the eruption, but in the late phase, adiabatic compression plays an important role in heating the plasma there. Thermal conduction also plays an important role in the transport of thermal energy away from the current sheet region throughout the reconnection process, producing a “thermal halo” and widening the region of high temperatures. We simulate emission from solar telescopes for this eruption and find that there is evidence for emission from heated plasma above the flare loops late in the eruption, when the adiabatic heating is the dominant heating term. These results provide an explanation for hot supra-arcade plasma sheets that are often observed in X-rays and extreme ultraviolet wavelengths during the decay phase of large flares.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗