Baseflow Age Distributions and Depth of Active Groundwater Flow in a Snow‐Dominated Mountain Headwater Basin
Not Available
SEARCH · Search NASA
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.
Not Available
In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The updated Hanford CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within and downgradient of the Central Plateau.
Abstract Water resources, including groundwater and prominent rivers worldwide, are under duress because of excessive contaminant and nutrient loads. To help mitigate this problem, the United States Department of Energy (DOE) has supported research since the late 1980s to improve our fundamental knowledge of processes that could be used to help clean up challenging subsurface problems. Problems of interest have included subsurface radioactive waste, heavy metals, and metalloids (e.g. uranium, mercury, arsenic). Research efforts have provided insights into detailed groundwater biogeochemical process coupling and the resulting geochemical exports of metals and nutrients to surrounding environments. Recently, an increased focus has been placed on constraining the exchanges and fates of carbon and nitrogen within and across bedrock to canopy compartments of a watershed and in river–floodplain settings, because of their important role in driving biogeochemical interactions with contaminants and the potential of increased fluxes under changing precipitation regimes, including extreme events. While reviewing the extensive research that has been conducted at DOE’s representative sites and testbeds (such as the Oyster Site in Virginia, Savannah River Site in South Carolina, Oak Ridge Reservation in Tennessee, Hanford in Washington, Nevada National Security Site in Nevada, Riverton in Wyoming, and Rifle and East River in Colorado), this review paper explores the nature and distribution of contaminants in the surface and shallow subsurface (i.e. the critical zone) and their interactions with carbon and nitrogen dynamics. We also describe state-of-the-art, scale-aware characterization approaches and models developed to predict contaminant fate and transport. The models take advantage of DOE leadership-class high-performance computers and are beginning to incorporate artificial intelligence approaches to tackle the extreme diversity of hydro-biogeochemical processes and measurements. Recognizing that the insights and capability developments are potentially transferable to many other sites, we also explore the scientific implications of these advances and recommend future research directions.
This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.
In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)
Abstract Groundwater pumping can reduce streamflow in nearby waterways (‘streamflow depletion’), a process which must be accounted for in integrated management of surface and groundwater resources. However, causal identification of streamflow depletion from hydrographs alone is challenging because pumping impacts are masked by other drivers of hydrologic variability. To identify potential indicators of streamflow depletion, we used synthetic hydrographs and an analytical streamflow depletion model to assess potential pumping impacts on specific hydrograph characteristics (‘hydrologic signatures’) for 215 streamgages spanning the conterminous United States (CONUS). We found that streamflow depletion commonly impacts signatures associated with seasonal and annual low flows and low flow recessions. The largest impacts occurred during dry years, suggesting streamflow depletion may be evident in dry years even where impacts are unmeasurable in wet years. Random forest models indicated that streamflow depletion could significantly impact Annual, Summer, and Fall signatures in most streams. Our finding that multiple hydrologic signatures are consistently responsive to streamflow depletion across CONUS suggests that the underlying hydrological processes linking pumping to streamflow reductions are consistent across diverse settings, information that will aid in identifying indicators of streamflow depletion from streamflow hydrographs.
Abstract A substantial fraction of nitrogen (N) fertilizer applied in agricultural systems is not incorporated into crops and moves below the rooting zone as nitrate (NO 3 − ). Understanding mechanisms for soil N retention below the rooting zone and leaching to groundwater is essential for our ability to track the fate of added N. We used dual stable isotopes of nitrate (δ 15 N–NO 3 − and δ 18 O–NO 3 − ) and water (δ 18 O–H 2 O and δ 2 H–H 2 O) to understand the mechanisms driving nitrate leaching at three depths (0.8, 1.5, and 3.0 m) of an irrigated corn field sampled every 2 weeks from 2016 to 2020 in the southern Willamette Valley, Oregon, USA. Distinct periods of high nitrate concentrations with lower δ 15 N–NO 3 − values indicated that a portion of that nitrate was from recent fertilizer applications. We used a mixing model to quantify nitrate fluxes associated with recently added fertilizer N versus older, legacy soil N during these “fertilizer signal periods.” Nitrate leached below 3.0 m in these periods made up a larger proportion of the total N leached at that depth (∼52%) versus the two shallower depths (∼13%–16%), indicating preferential movement of recently applied fertilizer N through the deep soil into groundwater. Further, N associated with recent fertilizer additions leached more easily when compared to remobilized legacy N. A high volume of fall and winter precipitation may push residual fertilizer N to depth, potentially posing a larger threat to groundwater than legacy N. Optimizing fertilizer N additions could minimize fertilizer losses and reduce nitrate leaching to groundwater.
Increasing competition for water resources in the United States could create future challenges for allocating and using thermal cooling water in the U.S. electric power sector. While thermal power plant retirements and the growth of wind and solar technologies can reduce national aggregate power sector cooling water use, local water constraints and growing demand for agricultural or municipal supply could create greater needs for higher-cost alternative water supplies such as groundwater or recycled wastewater. For some regions, these incentives could change future electricity planning and operational decisions. These relationships and impacts are studied here using the National Renewable Energy Laboratory (NREL) Regional Energy Deployment System (ReEDS), a national electric sector planning model that has recently been upgraded to include a highly detailed representation of thermal cooling water demand and supply. Thermal power technologies are differentiated by both cooling technology and water source type to track and constrain thermal cooling water use in a way that incorporates both physical and legal water considerations in the United States. This capability is exercised under a range of electricity sector futures with alternative technology costs, fuel prices, and water constraints to illustrate ways that U.S. electric sector water use could evolve under uncertain future electric sector drivers. In exploring changes to regional generation and transmission planning and operation, water requirements, and cost, we highlight the environmental and economic impacts of future power sector water decisions.
This document presents an evaluation of a potential groundwater tracer test in the uppermost RUM (Ringold Formation member of Wooded Island – upper mud unit) aquifer in the 100-H Area of the Hanford Site. The results from field electrical resistivity tomography (ERT) investigations, numerical flow and transport (F&T) modeling, coupled ERT-F&T model simulations, and laboratory resin experiments provide the technical basis for evaluating tracer test scenarios that minimize impacts to the 100-HX pump and treat (HX P&T) operations, and maximize hydrologic data acquisitions including the use of advanced surface geophysical techniques for monitoring tracer transport.
Among the most common amendments added to groundwater during site remediation are compounds used to adjust or maintain the pH. This research describes an approach to encapsulate mineral particles (MgO and CaCO 3 ) within oil droplets suspended within an aqueous phase for the purpose of delivery to the subsurface environment. A series of batch experiments was combined with mathematical modeling to illustrate the encapsulation and understand the influence of particle encapsulation on rates and extents of alkalinity release. The encapsulation of the alkalinity-releasing particles results in slower rates of amendment release as compared to rates obtained using suspensions of bare mineral particles, allowing for the possibility of control as a function of the pH. The results indicate that the alkalinity release from particle suspensions followed a mineral dissolution mechanism that could not explain the rate of the alkalinity release of the encapsulated particles. The reduction in mineral dissolution rates observed with the encapsulated particles was found to result from a mass transfer limitation. This limitation was well described using a linear driving force expression to account for the resistance to mass transfer at the oil–water interface.
A series of batch experiments were performed to assess the uranium sorption capacity of four mineralogically distinct lithologies from the Negev Desert, Israel, to evaluate the suitability of a potential site for subsurface radioactive waste disposal. The rock specimens consisted of an organic-rich phosphorite, a bituminous marl, a chalk, and a sandstone. The sorption data for each lithology were fitted using a general composite surface complexation model (GC SCM) implemented in PHREEQC. Sorption data were also fitted by a non-mechanistic Langmuir sorption isotherm, which can be used as an alternative to the GC SCM to provide a more computationally efficient method for uranium sorption. This is because all the rocks tested have high pH/alkalinity/calcium buffering capacities that restrict groundwater chemistry variations, so that the use of a GC SCM is not advantageous. The mineralogy of the rocks points to several dominant sorption phases for uranyl (UO 2 2+ ), including apatite, organic carbon, clays, and iron-bearing phases. The surface complexation parameters based on literature values for the minerals identified overestimate the uranium sorption capacities, so that for our application, an empirical approach that makes direct use of the experimental data to estimate mineral-specific sorption parameters appears to be more practical for predicting uranium sorption.
Overview and How-to Tutorial Videos for Using NEWTS Data Video 1: Overview of NEWTS Database Video 2: How-to tutorial for EPA Flue Gas Desulfurization (FGD) Effluent NEWTS dataset Video 3: How-to tutorial for USGS Produced Waters NEWTS dataset Video 4: How-to tutorial for EPA Ash NEWTS dataset Video 5: How-to tutorial for Quillinan, et al 2018 DOE Geothermal Technology Office REE dataset Video 6: Tutorial video on navigating the NEWTS Dashboard, with an overview of NEWTS and navigating between the NEWTS Dashboard and Datasets (https://netl-doe.maps.arcgis.com/apps/dashboards/a5fa4192f7c6478dab3d6180d9c30b84) Video 7: Additional tutorial video on navigating the NEWTS Dashboard and investigating specific data points in the Dashboard and Datasets Video 8: Re-record of recent webinar giving an overview of the NEWTS Database and Dashboard, including interacting with the NEWTS Dashboard, locating specific data points, and finding the relevant streams in the NEWTS Database and datasets on EDX. Includes overview of the datasets, case studies, and steps for taking stream data from the database and modeling stream data in OLI Studio and Geochemist's Workbench. Note: Video 3 tutorial is also applicable to the USGS Brackish Water NEWTS dataset.
Site operations at the Los Alamos National Laboratory (LANL) resulted in the release of oxidized chromium, Cr(VI), into Sandia Canyon from cooling tower effluent from 1956 until 1972. The chromium traveled with the surface water approximately 3 miles downstream before migrating below ground surface. Chromium concentrations exceed 50 μg/L in the upper portion of the aquifer. Another LANL groundwater plume of concern is associated with RDX (Royal Demolition Explosives, 1,3,5-trinitro-1,3,5-triazine). Between 1951 and 1996, RDX was released to the mesa-top facilities' process water outfall, adjacent and underlying soils, and alluvial sediments, along with surface water in Cañon de Valle. Between 2000 and 2010, two remedial actions were deployed, removing much of the near-surface RDX, however, recharge due to precipitation has transported RDX into the perched-intermediate zone and into the regional aquifer. The report documents an independent technical review by scientists from the Department of Energy (DOE) Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) to provide recommendations for potential near term actions to address and optimize remediation for both the Cr(VI) and RDX plumes. The proposed near-term remedial actions include design of pump and treat systems for Cr(VI) and monitoring and study for natural attenuation for RDX. The review assesses existing data, conceptual and numerical modeling, and it recommends a technical integration process to support identifying and implementing strategic, effective and efficient remedies. The DOE Environmental Management Los Alamos Field Office (EM-LA) and their cleanup contractor Newport News Nuclear-BWTX, LLC Los Alamos (N3B) provided the information required for the review. Interviews were also conducted with regulators to obtain the full spectrum of technical, regulatory and scientific perspectives. The independent review team was impressed by the capabilities, experiences, innovativeness, and insightfulness of the technical representatives from both the regulator, the New Mexico Environment Department (NMED) and N3B. Incorporation of vadose zone flow pathways in the conceptual site model (CSM) and configuring the numerical modeling for the site was generally state-of-the-practice (or better). This could be considered state-of-the-art by addressing uncertainties related to spatial extent of hydraulic windows. The reviews from the regulators were thorough and often provided useful concepts for consideration and future/study resolutions. The overarching consensus recommendation of independent review team is that the LANL groundwater plumes should be addressed in context of the emerging "management of complex sites" paradigm.
Using a projected end-state date of 2065 (SRNS 2015b), the Savannah River Site (SRS) Composite Analysis (CA) modeling for each facility and waste site began on the inventory year assigned to it so that source depletion and radionuclide transport out of the system could be appropriately captured. Some SRS waste sites that have already achieved their end states (i.e. end-state inventories and end-state configuration) are currently contributing to the potential off-site public dose through source release, groundwater transport, discharge to on-site surface streams, and stream transport to the CA point of assessments (POAs). The inventory year assigned to these waste sites is 2002 or before. This means that SRS CA results from 2002 and beyond are a reasonable representation for these waste sites that have already achieved their end states and are currently contributing to the potential off-site public dose. The SRS Annual Environmental Report (AER) monitoring can differentiate and separate liquid pathway data allowing the data representing only waste sites at their end state to be produced. Because the SRS CA has projected reasonable end-state impacts from 2002 and beyond, and the AER monitoring can differentiate and separate operating and end-state contributions to annual liquid pathway release, an opportunity exists to use the AER monitoring data to validate the SRS CA model. The CA model validation program uses a graded and systematic approach for taking corrective action, starting with an SRS established administrative dose limit of 15 mrem/yr, below which no action is required. Based on the location of the 2010 SRS CA POAs, the only potential exposure pathway for the public is through surface water. The completion of the FY2019 CA model validation indicates that the SRS CA projected dose, while generally conservative, provides a reasonable representation of the maximum annual doses. These doses are well below the administrative limit; therefore, no additional action is required.
Using a projected end-state date of 2065 (SRNS 2015b), the Savannah River Site (SRS) Composite Analysis (CA) modeling for each facility and waste site began on the inventory year assigned to it so that source depletion and radionuclide transport out of the system could be appropriately captured. Some SRS waste sites that have already achieved their end states (i.e., end-state inventories and end-state configuration) are currently contributing to the potential off-site public dose through source release, groundwater transport, discharge to on-site surface streams, and stream transport to the CA point of assessments (POAs). The inventory year assigned to these waste sites is 2002 or before. This means that SRS CA results from 2002 and beyond are a reasonable representation for these waste sites that have already achieved their end states and are currently contributing to the potential off-site public dose. The SRS Annual Environmental Report (AER) monitoring can differentiate and separate liquid pathway data allowing the data representing only waste sites at their end state to be produced. Because the SRS CA has projected reasonable end-state impacts from 2002 and beyond, and the AER monitoring can differentiate and separate operating and end-state contributions to annual liquid pathway release, an opportunity exists to use the AER monitoring data to validate the SRS CA model.
Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.
Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.
Abstract Electricity and water systems in the Western US (WUS) are closely connected, with hydropower comprising 20% of total annual WUS generation, and electricity related to water comprising about 7% of total WUS electricity use. Because of these interdependencies, the threat of climate change to WUS resources will likely have compounding electricity impacts on the Western Interconnect grid. This study describes a WUS‐wide water system model with a particular emphasis on estimating climate impacts on hydropower generation and water‐related electricity use, which can be linked with a grid expansion model to support climate‐resilient electricity planning. The water system model combines climatically‐driven physical hydrology and management of both water supply and demand allocation, and is applied to an ensemble of 15 climate scenarios out to 2050. Model results show decreasing streamflow in key basins of the WUS under most scenarios. Annual water‐related electricity use increases up to 4%, and by up to 6% during the summer months, driven by growing agricultural demands met increasingly through a shift toward energy‐intensive groundwater to replace declining surface water. Total annual hydropower generation changes by +5% to −20% by mid‐century but declines in most scenarios, with decreases in summer generation by up to nearly −30%. Water‐related electricity use increases tend to coincide with hydropower generation declines, annually and seasonally, demonstrating the importance of concurrently evaluating the climate signal on both water‐for‐energy and energy‐for‐water to inform planning for grid reliability and decarbonization goals.