Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4

A ModEx Framework for Watershed Subsurface Investigation With Limited Geophysical Data Using Machine Learning and Hydrologic Modeling

Abstract Subsurface heterogeneity influences watershed hydrology strongly but remains difficult to characterize at catchment scales with sparse and costly field data. Geophysical surveys such as electromagnetic induction (EMI) provide local spatial subsurface images yet scaling them to watershed scales and converting EMI‐derived resistivity into hydraulic properties remains a challenge. We present a Model–Experiment (ModEx) framework that integrates limited EMI data with machine learning (ML) and hydrologic modeling to improve process representation and guide field investigations. Sparse EMI surveys were scaled to the catchment scale using a Random Forest model, and the resulting resistivity fields were combined with nearby borehole constraints to parameterize a hydrologic model. The EMI‐informed hydrological simulations improved predictions of streamflow sustained by subsurface flow and shallow saturation patterns. By combining EMI data and ML with hydrologic modeling, the ModEx framework guides future subsurface surveys, providing a transferable and efficient strategy for data–model integration across diverse watersheds. Plain Language Summary Mapping the underground network of soil and rock that controls water is essential for predicting floods and droughts, but seeing underground is difficult and expensive. We cannot drill everywhere, so scientists use geophysical tools to scan broad areas. There are two key challenges: these geophysical scans are often sparse across the whole watershed, and the geophysical data is hard to translate into water‐related properties. We used artificial intelligence to solve these problems. We taught a computer to find patterns linking the limited geophysical data to the land surface properties. This allowed it to fill in the gaps and create a complete, useful subsurface map for the entire watershed. This new map improves hydrologic simulations, leading to more accurate predictions of water movement in the watershed. It also helps scientists build better models with less data and generates a priority map showing where to measure next, making future investigations more efficient. Key Points Limited EMI scaled with ML improves catchment‐scale subsurface parameterization for hydrologic models The framework integrates hydrologic modeling with limited geophysical data to support subsurface investigation design ModEx framework offers a transferable data–model integration strategy that quantifies and reduces uncertainty guiding watershed studies

Chen, Hang↗

The Role of Topography in Controlling Evapotranspiration Age

Abstract Evapotranspiration (ET) age is a key metric of water sustainability but a major unknown partly due to the extreme difficulty in modeling it. Groundwater is found to be important in ET age variations in small‐scale studies, yet our understanding is insufficient because groundwater systems are nested across scales. Here, we conducted GPU‐accelerated particle tracking with integrated hydrologic modeling to quantify the variations in ET age at a regional scale of ∼0.4 M km 2 . Simulation results reveal topography‐driven flow paths shaping the spatial and temporal patterns of ET age variations. On ridges, where root zone decoupling with deep subsurface storage, ET age is generally young, with seasonal variations dominated by meteorological conditions. In the valley bottom, ET age is generally old, with significant subseasonal variations caused by the convergence of subsurface flow paths. On hillslopes with water table depths ranging from 1 to 10 m, ET age shows strong seasonal variations caused by the connections with lateral groundwater regulated by ET demand. Our modeling approach provides insights into the basic linkages between ET age and topography at large scale. Our work highlights the perspective of multiscale studies of ET age, suggesting new field experiments to test these process connections and to determine if such linkages warrant inclusion in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗

Impacts into a Strength‐Layered Target: The Time Evolution of the Excavation‐Stage Flow

Impact cratering has dominated the evolution and modification of planetary surfaces throughout the history of the solar system. Impact craters can serve as probes to understanding the details of a planetary subsurface; for example, Oberbeck and Quaide suggested that crater morphology can be used to estimate the thickness of a regolith layer on top of a more competent unit. Lunar craters show a morphological progression from a simple bowl shape to flat‐floored and concentric craters as diameter increases for a given regolith thickness. The final shape of the impact crater is a result of the subsurface flow‐field initiated as the projectile transfers its energy and momentum to the target surface at the moment of impact. Therefore, when a regolith layer is present over a stronger substrate, such as is the case on the lunar mare, the substrate modifies the flow‐field resulting in the distinctive final crater morphology. Here we report on a series of experimental impacts into targets composed of a thin (2 cm) layer of loose sand on top of a stronger substrate. We use the Ejection‐Velocity Measurement System (EVMS, Fig.1) developed by [3] to examine the ejecta kinematics at specific times during the formation of these craters in strength‐layered targets.

J. L. B. Anderson↗

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

Interflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

hillslope hydrology↗

Non-linear hydrologic organization

We revisit three variants of the well-known Stommel diagrams that have been used to summarize knowledge of characteristic scales in time and space of some important hydrologic phenomena and modified these diagrams focusing on spatiotemporal scaling analyses of the underlying hydrologic processes. In the present paper we focus on soil formation, vegetation growth, and drainage network organization. We use existing scaling relationships for vegetation growth and soil formation, both of which refer to the same fundamental length and timescales defining flow rates at the pore scale but different powers of the power law relating time and space. The principle of a hierarchical organization of optimal subsurface flow paths could underlie both root lateral spread (RLS) of vegetation and drainage basin organization. To assess the applicability of scaling, and to extend the Stommel diagrams, data for soil depth, vegetation root lateral spread, and drainage basin length have been accessed. The new data considered here include timescales out to 150 Myr that correspond to depths of up to 240 m and horizontal length scales up to 6400 km and probe the limits of drainage basin development in time, depth, and horizontal extent.

58 GEOSCIENCES↗

Determining the biogeochemical transformations of organic matter composition in rivers using molecular signatures

Inland waters are hotspots for biogeochemical activity, but the environmental and biological factors that govern the transformation of organic matter (OM) flowing through them are still poorly constrained. Here we evaluate data from a crowdsourced sampling campaign led by the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to investigate broad continental-scale trends in OM composition compared to localized events that influence biogeochemical transformations. Samples from two different OM compartments, sediments and surface water, were collected from 97 streams throughout the Northern Hemisphere and analyzed to identify differences in biogeochemical processes involved in OM transformations. By using dimensional reduction techniques, we identified that putative biogeochemical transformations and microbial respiration rates vary across sediment and surface water along river continua independent of latitude (18°N–68°N). In contrast, we reveal small- and large-scale patterns in OM composition related to local (sediment vs. water column) and reach (stream order, latitude) characteristics. These patterns lay the foundation to modeling the linkage between ecological processes and biogeochemical signals. We further showed how spatial, physical, and biogeochemical factors influence the reactivity of the two OM pools in local reaches yet find emergent broad-scale patterns between OM concentrations and stream order. OM processing will likely change as hydrologic flow regimes shift and vertical mixing occurs on different spatial and temporal scales. As our planet continues to warm and the timing and magnitude of surface and subsurface flows shift, understanding changes in OM cycling across hydrologic systems is critical, given the unknown broad-scale responses and consequences for riverine OM.

rivers↗

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗

Tracer Gas Model Development and Verification in PFLOTRAN

Tracer gases, whether they are chemical or isotopic in nature, are useful tools in examining the flow and transport of gaseous or volatile species in the underground. One application is using detection of short-lived argon and xenon radionuclides to monitor for underground nuclear explosions. However, even chemically inert species, such as the noble gases, have bene observed to exhibit non-conservative behavior when flowing through porous media containing certain materials, such as zeolites, due to gas adsorption processes. This report details the model developed, implemented, and tested in the open source and massively parallel subsurface flow and transport simulator PFLOTRAN for future use in modeling the transport of adsorbing tracer gases.

07 ISOTOPE AND RADIATION SOURCES↗

The Lakshmi Plateau structure as an indicator of asthenosphere horizontal flows on Venus

The structure of Lakshmi Planum in the western part of Ishtar Terra in a fold-fault setting which conforms to the basic massif of the plateau with eruptive centers is constructed concentrically and is interpreted from the point of view of the subsurface flow of materials in the form of horizontally diverging asthenospheric flows and gravitational creep. The surrounding structures are formed by the deformation of the more rigid lithosphere as it breaks away from the asthenospheric flow.

Pronin, A. A.↗

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗

Geophysical Methods for Stratigraphic Identification

The goal of this work is to identify geophysical methods that can be used to delineate (in 3D) a large-scale paleochannel that has been inferred to exist between the 200 West and 200 East Areas of the Hanford Site’s Central Plateau . This paleochannel is of particular importance because it will significantly impact flow and contaminant transport from 200 West to 200 East. Its high permeability can both accelerate contaminant transport and lower contaminant concentrations through spreading. Information about the spatial extent and structure of such features is needed to inform subsurface flow and contaminant transport models that support remedy decisions at Hanford.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparison of Ring Diagrams Based on the Doppler Shifts of Synthetic Data Obtained with Bisector Method and SDO/HMI Pipeline

Ring diagrams are cross-sections of three-dimensional spatiotemporal power spectra of solar oscillations. The rings reveal information about sub-surface flows and represent an important tool for helioseismology. Ring diagrams can be constructed using Doppler velocity or intensity maps of spectral lines. How the velocities are computed is an important factor for accuracy of information we can retrieve on subsurface flows. In our work, the ring diagrams are generated from Doppler shift data of synthesized Fe I 6173 Å line. We compare ring diagrams computed by two methods–HMI line-of-sight pipeline and the bisector of Fe line. Fe I line is synthesized for StellarBox 3D Radiative hydrodynamic simulations under LTE assumption. We aim to answer the following questions: 1.How do power spectra obtained from velocities computed with the HMI pipeline and bisector of the Fe I 6173Å compare? 2.How do the power spectral density retrieved with each method vary with heliocentric angle? 3.What is the effect of changing resolution on the power spectral density in ring diagrams obtained with the two methods?

SMD↗

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Understanding the relative importance of vertical and horizontal flow in ice-wedge polygons

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

54 ENVIRONMENTAL SCIENCES↗