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

Evaluating Remediation Options for a Perched-Water Unit at the Hanford Site - 20378

At the U.S. Department of Energy's Hanford Site, liquid-phase waste was discharged using surface disposal features such as ponds, trenches, and cribs. These practices distributed contaminants in the subsurface, some of which are currently in the vadose zone above the groundwater. At the 200-DV-1 operable unit (OU) at the Hanford Site, a key feature is a perched-water zone. This perched-water zone contains several million liters of water contaminated with high concentrations of uranium, technetium-99 (Tc-99), and nitrate that serves as a continuing source of contamination to the underlying aquifer. Currently, a pumping system is used to remove contaminant mass from the perched-water zone. However, water extraction will become more challenging as the perched unit dewaters. As a result, residual contaminants remaining after pumping is no longer effective needs to be addressed with post-extraction mitigation approaches, appropriate monitoring, and selection of suitable remediation endpoints. Thus, a remedy evaluation needs to include the unique conditions of this perched-water zone. As part of the initial evaluation, a set of technologies have been preliminarily identified as potentially applicable to the perched-water zone, including in situ remediation technologies to sequester uranium and Tc-99 and to degrade nitrate. While there is sufficient information available for direct consideration of some of the reviewed technologies in a feasibility study, five in situ remediation technologies were identified as needing additional information on their effectiveness under the geochemical and contaminant conditions of the perched-water zone. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Applications of Hyporheic Zone Sampling for Assessing VOC Plume Engineered Remediation and Natural Attenuation Effectiveness - 20414

Determining the efficacy of engineered remedial actions and/or natural attenuation as an acceptable EPA Record of Decision (ROD) must address long-term temporal and spatial variability. Since final plume downslope destinations are aquifers and related surface streams and rivers, hyporheic zone sampling is essential. Comprehensive hyporheic sampling documents plume groundwater remediation prior to upwelling into stream benthic strata and natural attenuation within this dynamic zone itself. Our intensive sampling of the Savannah River Site (SRS) Chemical, Metals, and Pesticides (CMP) Pits VOC plume-fringe from 2005 to 2018 documented PCE natural attenuation and effectiveness of In Situ Thermal Treatment (ISTT). Most plume hyporheic upwelling was delimited to a 20 m inflow reach along Pen Branch. Our study utilized corroborative sampling methodologies: ambient hole-water grab samples: passive diffusion bags (PDB); and hole-core sediment samples to confirm station to station spatial and temporal trends. VOC hot spot detections were confirmed from both groundwater and soil core samples. VOC pathways through the Pen Branch hyporheic zone were assessed at hole-installation using sequential soil samples from augered cores and grab-samples of hole-bottom water. Longer term, composite water samples were collected from PDBs deployed for a minimum two-week equilibration period. Soil samples were collected from the 15 and 65 cm core levels in December 2017 and June 2018 at six stations. Additional core depths of 40 and 85 cm were sampled in June 2018 at the highest VOC station (5DB80). All samples were analyzed in an EPA certified lab by purge and trap GC-MS. Effectiveness of previous ISTT and soil vapor extraction (SVE) for this CMP Pits PCE plume, was clearly documented by significant downward trends in VOC at hyporheic zone stations prior to December 2015. Achieving final stage conversions of PCE daughter products to ethene and CO{sub 2} may be limited depending upon suitable interactions between hyporheic zone sediment morphology, hydrology, and redox conditions. Because declines in VOC daughter products at our hyporheic stations may also be caused by physical actions of dilution and dispersion rather than true biochemical dechlorination, full quantification of natural attenuation requires the use of compound-specific isotope ratio analysis (CSIA) and microbial 16S rRNA with qPCR. To assist these next-phase efforts, a natural attenuation index (NAI) was used to help produce quicker, more cost-effective decision making for locating CSIA and microbial sampling efforts. The NAI used natural logarithm transformed GC-MS results for station VOC degradation compounds versus VOC parent compounds. Wide differences were observed between stations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

235-F GoldSim Fate and Transport Model: Uncertainty Quantification

Building 235-F was configured with two missions in mind: Actinide Billet Line (ABL) and the fabrication of Pu-238 oxide for space program applications. ABL produced Np-237 billets for use in SRS reactors, whereas the design process, fabrication, and examination of Pu-238 oxide powder occurred within the following areas, respectively: Plutonium Experimental Facility (PEF), Plutonium Fuel Form (PuFF), and Old Metallography Lab (OML). By 1990 production ceased and by 2006 de-inventory occurred; however, assays have shown significant holdup remains within ABL and PuFF. As a result, 235-F is a Category 2 nuclear facility, with plans to undergo deactivation and decommission (D and D) via In-Situ Disposal (ISD). The purpose of this project is to ensure United States Environmental Protection Agency (USEPA) groundwater radiation maximum contaminant level (MCL) and dosage standards are met during the D and D of 235-F by quantifying uncertainty through probabilistic modeling and evaluation of various ISD alternatives. GoldSim is a dynamic modeling software package with a graphical, object-oriented interface capable of capturing the influence of complex system input variability on probabilistic system outcomes. A GoldSim stochastic fate and transport model for 235-F was developed and matched with a PORFLOW deterministic model to simulate probabilistic release and flow of radionuclides from ABL and PuFF into the vadose zone, the Upper Three Runs (UTR) Aquifer, and UTR Creek. The GoldSim model was used to probabilistically evaluate four ISD scenarios against USEPA groundwater radiation MCLs and dosage standards. The deterministic 235-F GoldSim fate and transport model continues to be refined to match the results of the PORFLOW deterministic model to ensure the model accurately represents radionuclide movement through the groundwater system. The stochastic variables that are utilized within the GoldSim model are founded on the most current data; a conservative perspective is taken where needed. Alignment with the PORFLOW deterministic model, coupled with input stochastic variability, allows the probabilistic 235-F GoldSim model to capture the conservative breadth of possible outcomes for radionuclide fate within this particular system. This ensures that the USEPA MCLs and dosage limits hold even in the worst case scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Deep Isolation-Development of the Safety Case for Disposal of Radioactive Wastes in Horizontal Boreholes - 20028

Deep Isolation has developed a safe, secure, and permanent geological disposal method for high-level waste, including spent nuclear fuel as well as sealed sources and other highly radioactive materials. The method leverages well established directional drilling technology to create horizontal repositories deep underground. The combination of great depth (1-3 km) and the ability to precisely position repositories in a horizontal orientation provides access to geologic strata that are inaccessible for typical mined and deep vertical borehole repositories. A number of potential direct and indirect safety benefits accrue including:1) greater depth below regional freshwater aquifers; 2) increased flexibility and more geologic options for siting; 3) access to deep formations with persistent and sustained reducing conditions. In addition there are a number of safety elements related specifically to the horizontal repository geometry including: 1) passive direction of thermally driven fluid and radionuclide movement away from the vertical access hole and toward the 'dead end' portion of the repository; 2) mitigation of seismic hazards by orienting repositories parallel to local and regional fault structures. In this paper we explore and discuss some of the key technologic, geologic and hydrologic elements that support the deep horizontal borehole safety case. The stalemate seen across the globe on the disposal of spent nuclear fuel and high-level waste can be broken. Deep Isolation offers a novel option for safe, secure, and permanent deep geological disposal of nuclear waste that can be developed as centralized repositories or adapted to smaller regional or site-specific repositories located near waste sources. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Development of a Conditioned System-Level Groundwater Model to Evaluate Long-Term Groundwater Impacts within a Performance Assessment - 20115

A preliminary performance assessment (PA) of single-shell tank Waste Management Area (WMA) A-AX located at the U.S. Department of Energy (DOE)'s Hanford Site in southeastern Washington is being conducted to satisfy the requirements of DOE Order 435.1 [1] as it relates to closure of the single-shell radioactive waste tanks in the WMA A-AX tank farms. A PA assesses the fate, transport, and impacts of radionuclides within a low-level radioactive waste disposal facility in its assumed closure configuration and the subsequent potential doses to humans over a 1,000-year compliance period and a 10,000-year performance evaluation period. The WMA A-AX preliminary PA evaluation is structured around the complementary use of process-level and system-level models to calculate the facility performance against established DOE Order 435.1 [1] performance objectives. Process-level models are those that represent a detailed phenomenological representation of processes of concern in the PA. Process models typically only represent one or a few of the components of the PA, such as groundwater flow and transport, and must be integrated with other modeling elements to perform PA calculations. System-level models are those that are abstracted from the process models, retaining the essential features of the process model, while allowing integration of all aspects of the PA in a single modeling framework. System-level models are often characterized by coarser numerical discretization, lower dimensionality, or other similar simplifications compared to the process-level model. Traditionally, a three-dimensional (3-D) process-level model is utilized primarily to evaluate the long-term impact on groundwater and the potential doses to individuals who consume contaminated groundwater. System-level models are also utilized to evaluate the groundwater pathway in PAs. These models typically have reduced dimensionality (1-D) and are conditioned utilizing flow fields (Darcy fluxes) and moisture content distributions that are abstracted from the process level models. The abstraction approach assures that the flow field in both models is consistent for a specific set of input parameters for flow, differing only in the discretization and dimensionality of the two models. The preliminary WMA A-AX PA incorporates a detailed representation of the geological system and hydraulic properties within the 3-D model STOMP{sup C} numerical code so that the effects of relevant features and processes on water flow and radionuclide transport in the subsurface can be evaluated. The complementary system-level model is developed utilizing the GoldSim{sup C} code to implement a simplified, 1-D equivalent model to represent the groundwater pathway. Contaminant transport through the vadose zone and unconfined aquifer for Tc-99 and I-129 were evaluated in each of the models. Adjustments to the saturated portion of the GoldSim{sup C} 1-D model were required to mimic the dispersion effect captured with the 3-D STOMP{sup C} model. Once this conditioning was conducted, highly similar results for the transport of Tc-99 and I-129 were achieved at the point of calculation, located 100 m downgradient from the WMA A-AX fenceline. These radionuclides represent elements that are regarded as primary dose drivers in the PA groundwater pathway analysis. The high degree of conformance between the two models suggests that the use of the equivalent 1-D system model is suitable for evaluating the full suite of radionuclides that will be released from the tank sources within WMA A-AX over the 1,000-year compliance period and over the 10,000-year evaluation period. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Study of an Unrefined Humate Solution as a Possible Attenuation-based Remedy for Uranium Contamination in Acidic Groundwater - 20351

The Savannah River Site (SRS), joined the environmental cleanup program in 1981 after the Cold War, and was designated as a hazardous waste management facility. About 1.8 billion gallons of acid wastes were disposed into the F/H Area Seepage Basin that included many radionuclides and dissolved metals, resulting in highly contaminated groundwater plumes with pH of 3 - 5.5. The acidity of the plumes contributes to the mobility of several constituents of concern (COC) such as H-3, U-238, I-129, and Sr-90 for the F-Area plume and H-3, Sr-90 and mercury for the H-Area plume. An active treatment (pump-and-treat and re-injection) system was constructed and put in place in 1997 to address the removal of metals and radionuclides from the contaminated groundwater. As remediation projects advanced, active treatments were transitioned to more passive or enhanced-passive approaches such as the subsurface barrier with base injection system, which has been controlling the mobility of contaminants. However, new and more efficient attenuation-based remedies are always sought. Several studies have proposed that humic substances can be used to remediate sites contaminated with heavy metals by creating a permeable reactive barrier. Humic substances (HS) are major components of soil organic matter. HS are polyfunctional organic macromolecules formed by the chemo-microbiological decomposition of biomass or dead organic matter. These substances are usually divided into three main fractions: humin (insoluble at all pH values), humic acid (soluble at pH greater than 3.5), and fulvic acid (soluble at all pH values). Humic substances are major components of the soil at SRS and, in addition, are also helpful in the removal process of uranium. This research focuses on uranium (VI), which is a key contaminant of concern in the F-Area groundwater plume. The interaction of uranium with sediment in the presence or absence of humic substances involves complex mechanisms that are not yet well understood. The interaction of U(VI) with humic substances can affect the adsorption of U on sediment, altering its mobility in the subsurface. The objective of this study is to study the mobility of uranium in the presence of humic acid and to determine if sediments amended by humic acid can enhance the sorption of heavy metals onto sediments. In this research, experiments were conducted using background (clean) F-Area aquifer sediments to understand and predict uranium mobility. A blended material, humate, containing both humic and fulvic acids and a chemically modified humate (KW15 modified humic or mod-HA) was studied as a possible amendment for uranium remediation in SRS groundwater. Batch experiments were conducted in triplicates using 20 ml of aqueous suspension in DI water containing SRS sediment and mod-HA and studied the effect of pH, kinetics of humate sorption onto SRS sediment and also studied the sorption of uranium onto mod-HA amended SRS sediment. Results indicate that mod-HA amended sediment increases uranium removal significantly. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Uranium Mill Tailings Radiation Control Act Title II DOE Due Diligence and Lessons Learned from a Previous Site Transfer - 20352

Title II of the Uranium Mill Tailings Radiation Control Act (UMTRCA) established that a government agency will provide perpetual care for closed uranium and thorium ore-processing sites that were operating under an NRC source material license in 1978 or were licensed thereafter. Commercial owners (licensees) operating under an NRC or agreement state specific license when UMTRCA was passed are responsible for conducting reclamation of any byproduct material remaining from uranium-ore processing operations in accordance with an NRC or agreement state approved reclamation plan. Reclamation includes both surface and groundwater remedies. Upon completion of reclamation and approval by NRC, the site is required to be transferred to either the host state or the DOE for long-term surveillance and maintenance. Since UMTRCA's enactment, six Title II sites have been transferred to DOE; an additional 24 Title II sites are anticipated to be transferred before 2050. DoE's role mandated under UMTRCA Title II as the long-term care custodian is to perform 'monitoring, maintenance, and emergency measures necessary to protect the public health and safety.' UMTRCA requires that the licensee pay a long-term surveillance charge 'sufficient to cover the annual costs of site surveillance.' However, at some sites such as the Bluewater, New Mexico, Disposal Site, this mandate has required additional effort and expense by DOE, beyond the originally anticipated and intended scope within UMTRCA, but within the authority of DOE under UMTRCA. In 1997, the Bluewater site became the second UMTRCA Title II site to be transferred to DOE. The site was the location of a uranium mill operated from 1953 until 1982. The specific licensee began site reclamation in 1991, and by 1995 all tailings and contaminated materials were encapsulated in two tailings disposal cells and other disposal areas. In addition to surface contamination, milling activities impacted groundwater in the two upper aquifers. In 1989, the specific licensee attempted active groundwater remediation; however, no significant reduction in contaminant concentrations was observed. As a result, the specific licensee applied to NRC for alternate concentration limits (ACLs) in 1990, which were approved in 1996 as being protective, after additional corrective actions were performed. Since transfer of the Bluewater site to DOE, unforeseen challenges have occurred, requiring additional actions. The first challenge is the occurrence of surface depressions located on the northern section of the main tailings disposal cell. The depressions were first observed during DoE's initial inspection in 1998; however, evidence of these can be observed on satellite images taken prior to transfer. Since being first observed, the depressions have continued to grow both in depth and areal extent. Due to the design of the main tailings disposal cell, the depressions impede storm water from being effectively shed off the 101- hectare (250-acre) top slope of the main tailings disposal cell. Instead, storm water accumulates in the depressions, forming a large ephemeral pond that has stored up to 16.3 x 10{sup 6} liters (4.3 million gallons) of stormwater. The ponding poses a potential risk to the integrity of the main tailings disposal cell in the case of a large storm event, with the potential to cause the pond to overtop and erode the cover material and underlying waste. DOE has taken a number of short-term actions to monitor, measure, and reduce the ponding and is currently working with the US Army Corps of Engineers to design and construct a repair. Additional challenges are associated with groundwater at the Bluewater site. Nine wells were present on the 1335-hectare (3300-acre) site upon transfer. Groundwater compliance was called into question after the State of New Mexico reduced its uranium groundwater standard from 5.0 to 0.03 milligrams per liter in 2004, and when an ACL for uranium was exceeded in a site monitoring well in 2010. Acquiring historical groundwater data and subsequent evaluations as well as additional DOE groundwater monitoring led to installing 10 new monitoring wells and performing additional site hydrogeology recharacterization. DOE continues to evaluate groundwater conditions at the site and works with the NRC to determine regulatory requirements and a path forward. As a result of lessons learned at the Bluewater site and other Title II sites, improved processes have been implemented at a programmatic level to increase due diligence before site transfer and prevent similar issues from occurring at other UMTRCA Title II sites under long-term management. DoE's due diligence process is documented in the Process for Transition of UMTRCA Title II Disposal Sites to DOE for Long-Term Surveillance and Maintenance and is designed to ensure that DOE has no technical or compliance concerns with regulatory decisions that might compromise protectiveness following site transfer to DOE. Implementation of the due diligence process has increased DoE's role prior to site transfer and has been effective in identifying potential issues. Actions by NRC and specific licensees, in response to enhanced due diligence efforts by DOE, are expected to minimize, if not totally prevent, the need for unanticipated actions by DOE pertaining to the surface and groundwater remedies after site transfer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

California Trees Seasonally Use Augmented Water Sources: Water Isotope Tracking in a Groundwater‐Dependent Ecosystem

Sustainable groundwater management must account for the needs of groundwater dependent ecosystems. To understand the relationship of ecosystems and seasonal water use, we studied the stable isotope composition (δ 18 O and δ 2 H) of water in streamside trees in a semi-arid streamside environment (Livermore, California, USA). We sampled seven trees at two sites every other month from April 2024 through April 2025 for tree xylem stable water isotope signatures. These data were compared to potential source waters: precipitation, imported surface water, soil water and regional groundwaters. Large daily precipitation events were found to be isotopically similar to regional groundwater and were thus treated as one water source. A Bayesian mixing model using stable water isotopes was used to determine the ratios of these three potential source waters (small daily precipitation events, groundwater/large daily precipitation events and imported water) present in tree xylem water. On average, tree water sources include 32% imported water (SD = 9%), 35% small daily precipitation events (SD = 10%) and 33% groundwater (SD = 3%), with significant seasonal variation (t summer-winter = 30.8, p < 0.01), particularly drawing more imported water (more than 55%) in the summer. While small daily precipitation events contribute only 10% of the total precipitation in our dataset, it represents a third of water used by trees. In addition, while these ecosystems are designated as groundwater dependent ecosystems, the trees use approximately one third imported water and even more during dry summer months. This approach provides water managers with a practical tool for quantifying ecosystem water needs, supporting data-driven decisions and regulatory compliance. While California's Sustainable Groundwater Management Act emphasizes supporting groundwater dependent ecosystems, it allows flexibility in demonstrating benefits. In conclusion, our methodology provides a way to document that management actions (such as managed aquifer recharge with imported water) deliver measurable co-benefits to GDEs.

Geosciences↗
DOE OSTI2021

The stable isotopes of natural waters at the Marcell Experimental Forest

Here, we present a new data set from the Marcell Experimental Forest (MEF) that compiles water isotope measurements from multiple research catchments, some of which have been studied since the 1960s. The MEF is located in northern Minnesota, USA, and is home to heavily studied and monitored forests, streams, bogs, and fens. Peat-forming systems (bogs and fens) are an important component of the MEF landscape and have a profound impact on the water cycle in these catchments. Within the last decade, analysis of stable isotopes of water (expressed as δD and δ 18 O) has been implemented to characterize the different components of the water budget, and to allow researchers to look at catchment and peatland-specific hydrologic effects in the watershed. This δD and δ 18 O data set of natural waters from MEF catchments is primarily composed of measurements from three peatlands (S1, S2, S6) during an 11-year period. More recently collection and analysis were expanded to also include samples from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project in the S1 bog, peatlands S3, S4, S5, as well as nearby lakes. We establish a local meteoric water line by analyzing the isotopic composition of precipitation, which fills a void in regional meteoric water lines for Minnesota. Furthermore, we establish baseline isotopic composition for bog outlet streams, bog porewater, aquifer groundwater, overland flow, subsurface stormflow, and snowpack, as well as runoff from the SPRUCE experimental chambers. These data are publicly available and will be expanded upon in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Thermal insulation versus capacitance: A simulation experiment comparing effects of shade and hyporheic exchange on daily and seasonal stream temperature cycles

In streams where water temperatures stress native biota, management of riparian shade or hyporheic exchange are both considered viable management strategies for reducing the peaks of daily and seasonal stream channel temperature cycles. Although shade and hyporheic exchange may have similar effects on stream temperatures, their mechanisms differ. Improved understanding of the heat-exchange mechanisms influenced by shade and hyporheic exchange will aid in the appropriate application of either stream temperature management strategy. To illustrate a conceptual model highlighting shade as ‘thermal insulation’ and hyporheic exchange imparting ‘thermal capacitance’ to a stream reach, we conducted an in-silico simulation modelling experiment increasing shade or hyporheic exchange parameters on an idealized, hypothetical stream. We assessed the potential effects of increasing shade or hyporheic exchange on a stream reach using an established process-based heat-energy budget model of stream-atmosphere heat exchange and incorporated an advection-driven hyporheic heat exchange routine. The model tracked heat transport through the hyporheic zone and exchange with the stream channel, while including the effects of hyporheic water age distribution on upwelling hyporheic temperatures. Results showed that shade and hyporheic exchange similarly damped diurnal temperature cycles and differentially altered seasonal cycles of our theoretical stream. In winter, hyporheic exchange warmed simulated channel temperatures whereas shade had little effect. In summer, both shade and hyporheic exchange cooled channel temperatures, though the effects of shade were more pronounced. Our simple-to-grasp analogies of ‘thermal insulation’ for shade effects and ‘thermal capacitance’ for hyporheic exchange effects on stream temperature encourage more accurate conceptualization of complex, dynamic heat exchange processes among the atmosphere, stream channel, and alluvial aquifer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Supercell precipitation contribution to the United States hydroclimate

Abstract This research seeks to understand simulated supercell precipitation characteristics across the conterminous United States (CONUS) using high‐resolution, convection‐permitting, dynamically downscaled simulations for three 15‐year epochs. Epochs include a historical end‐of‐20th‐century period (1990–2005) and two end‐of‐21st‐century (2085–2100) scenarios for intermediate and pessimistic greenhouse gas concentration trajectories. Simulated updraft helicity, which measures the corkscrew flow within a storm's updraft, is used as a proxy for supercells. An algorithm tracks and catalogues updraft helicity swaths that, when buffered, are used to acquire simulated precipitation from supercells. The historical epoch provides a baseline climatology of supercell precipitation for a contemporary climate, which is then compared against the two future epochs to assess how supercell precipitation may change during the 21st century. Despite their relatively small size, supercells provide critical precipitation to the Wheat and Corn Belts, large expanses of CONUS pasture and rangeland, regional aquifers and several large river basins. Many areas in the central CONUS receive upwards of 3%–6% of their annual and 5%–8% of their warm‐season precipitation from these storms. Results suggest that precipitation contribution from supercells will decrease in the future across most of the High Plains and Central and Northern Great Plains with robust increases likely across the south‐central and Southeast regions. Supercell precipitation rates are expected to increase for large portions of the CONUS by the end‐of‐the‐21st‐century, suggesting a growing threat for flash floods from these storms as they become more efficient precipitation producers. This research provides an initial perspective on the magnitude of supercell precipitation and potential changes to this important hydrologic input to assist water‐sensitive industries, private and public insurance markets, agriculture entities, as well as inform plans to mitigate and build resilience to rapid environmental and societal change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Effects of structure and volcanic stratigraphy on groundwater and surface water flow: Hat Creek basin, California, USA

Abstract Hydrogeologic systems in the southern Cascade Range in California (USA) develop in volcanic rocks where morphology, stratigraphy, extensional structures, and attendant basin geometry play a central role in groundwater flow paths, groundwater/surface-water interactions, and spring discharge locations. High-volume springs (greater than 3 m 3 /s) flow from basin-filling (<800 ka) volcanic rocks in the Hat Creek and Fall River tributaries and contribute approximately half of the average annual flow of the Pit River, the largest tributary to Shasta Lake. A hydrogeologic conceptual framework is constructed for the Hat Creek basin combining new geologic mapping, water-well lithologic logs, a database of active faults, LiDAR mapping of faults and volcanic landforms, streamflow measurements and airborne thermal infrared remote sensing of stream temperature. These data are used to integrate the geologic structure and the volcanic and volcaniclastic stratigraphy to create a three-dimensional interpretation of the hydrogeology in the basin. Two large streamflow gains from focused groundwater discharge near Big Spring and north of Sugarloaf Peak result from geologic barriers that restrict lateral groundwater flow and force water into Hat Creek. The inferred groundwater-flow barriers divide the aquifer system into at least three leaky compartments. The two downstream compartments lose streamflow in the upstream reaches (immediately downstream of the groundwater-flow barriers) and gain in downstream reaches with the greatest inflows immediately upstream of the barriers.

Geology↗
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2022

Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock

Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.

58 GEOSCIENCES↗
DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
DOE OSTI2022

GW-PINN: A deep learning algorithm for solving groundwater flow equations

Machine learning methods provide new perspective for more convenient and efficient prediction of groundwater flow. In this study, a deep learning method “GW-PINN” without labeled data for solving groundwater flow equations with wells was proposed. GW-PINN takes the physics inform neural network (PINN) as the backbone and uses either the hard or soft constraint in the loss function for training. A locally refined sampling strategy (LRS) is adopted to generate the consistent spatial sampling points for problems with strong hydraulic head change, and then combined with an appropriate temporal sampling scheme to obtain the final spatial-temporal sampling points. A snowball-style two-stage training strategy by dividing the temporal domain into two subdomains is designed to decrease the sampling points. Five cases were designed to test the training performance of GW-PINN under different sampling strategies and two constraints. The predicted results of GW-PINN were compared with MODFLOW and the analytical solution. The results demonstrate that GW-PINN possesses strong ability in capturing the hydraulic head change for both confined and un-confined aquifers. The hard constraint owns more robust learning ability than the soft constraint. The LRS strategy can generate more accurate results with much fewer sampling points than traditional sampling strategies, and the snowball-style two-stage training strategy is significantly efficient for problems with the drastic change of hydraulic head. Additionally, the application of GW-PINN as a surrogate model for parameterized groundwater flow equations is illustrated. This study provides an option tool for efficient groundwater flow simulation, especially for those with local refinements are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Numerical investigations to identify environmental factors for field-scale reactive transport of pathogens at riverbank filtration sites

While induced bank filtration is a proven method for facilitating sustainable drinking water production, it is at risk from surface water contaminations (e.g., pathogens). Induced bank filtration and pathogen transport in groundwater have been studied extensively. However, long-term studies that consider real-world conditions are missing. These conditions include seasonal changes to environmental conditions and waterworks operations. Therefore, to analyze the effect of seasonal changes on the transport of human pathogenic viruses and their indicators in induced bank filtration, concentrations of adenoviruses and pathogen indicators were monitored over 16 months at an active bank filtration plant at the Rhine River, in Düsseldorf (Germany). Based on this data, a 2D groundwater model was created in PFLOTRAN that simulated flow, heat transport, conservative transport of chloride and the resulting electrical conductivity, reactive transport of oxygen and nitrate, and colloid-based transport of coliforms, somatic coliphages, and adenoviruses. The results show that reduced travel time was the key factor determining periods with a low removal of coliforms and somatic coliphages in the aquifer. Travel time was controlled by river level variations during rainy seasons, and the waterworks extraction rates during dry seasons. Further, for adenovirus transport, travel distance in the subsurface appeared to be the key factor, while travel time had no significant impact. Coliform removal increased when the colmation layer permeability decreased, while coliphage and adenovirus removal was unaffected by the colmation layer permeability. Seasonal changes in temperature and oxygen content did not significantly impact the removal of coliphages and adenoviruses in groundwater. Denitrifying conditions correlated with a lowered coliform removal, but the modelling could not establish a connection between denitrifying conditions and coliform removal. Our study showed that removal of pathogens and pathogen indicators at induced bank filtration plants varies greatly in time and space (e.g., for coliforms from 1 to 4 log-levels at 20m travel distance), and that adenovirus transport differs considerably from transport of coliforms and somatic coliphages.

59 BASIC BIOLOGICAL SCIENCES↗