Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater Modeling”

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
NASA NTRS2021

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land-surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snow pack evolution now support a range of drought indicators that better reflect coupled water, energy and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Drought↗
DOE OSTI2020

Component-In-Grout Model Implementation for the E-Area Low-Level Waste Facility's Performance Assessment

The component-in-grout (CIG) disposal units are below grade earthen trenches that contain grout encapsulated waste components. Components disposed of within the CIG segments consist of large radioactively contaminated equipment and smaller waste forms (e.g., B-25 boxes and SeaLand containers) to fill the space around and above the large equipment. In the 2022 revision of the E-Area Low-Level Waste Facility’s (ELLWF) performance assessment, groundwater radionuclide contaminant transport through the vadose zone will be modeled using the PORFLOW software package for nine existing CIG segments located within the Slit Trench (ST) 23 footprint - no additional CIG segments are planned for at this time. The nine CIG segments have been placed within two of the nominally 20-foot-wide by 656-foot-long ST segments and are surrounded by no less than 1 foot of grout on all sides and up to 4 feet of backfill material. The two most recently placed CIG segments have an additional reinforced concrete mat to improve structural stability.

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

Permafrost Promotes Shallow Groundwater Flow and Warmer Headwater Streams

Abstract The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from 11 headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7% of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5% moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Policy Innovation and Governance for Irrigation Sustainability in the Arid, Saline San Joaquin River Basin

This paper provides a chronology and overview of events and policy initiatives aimed at addressing irrigation sustainability issues in the San Joaquin River Basin (SJRB) of California. Although the SJRB was selected in this case study, many of the same resource management issues are being played out in arid, agricultural regions around the world. The first part of this paper provides an introduction to some of the early issues impacting the expansion of irrigated agriculture primarily on the west side of the San Joaquin Valley and the policy and capital investments that were used to address salinity impairments to the use of the San Joaquin River (SJR) as an irrigation water supply. Irrigated agriculture requires large quantities of water if it is to be sustained, as well as supply water of adequate quality for the crop being grown. The second part of the paper addresses these supply issues and a period of excessive groundwater pumping that resulted in widespread land subsidence. A joint federal and state policy response that resulted in the facilities to import Delta water provided a remedy that lasted almost 50 years until the Sustainable Groundwater Management Act of 2014 was passed in the legislature to address a recurrence of the same issue. The paper describes the current state of basin-scale simulation modeling that many areas, including California, are using to craft a future sustainable groundwater resource management policy. The third section of the paper deals with unique water quality issues that arose in connection with the selenium crisis at Kesterson Reservoir and the significant threats to irrigation sustainability on the west side of the San Joaquin Valley that followed. The eventual policy response to this crisis was incremental, spanning two decades of University of California-led research programs focused on finding permanent solutions to the salt and selenium contamination problems constraining irrigated agriculture, primarily on the west side. Arid-zone agricultural drainage-induced water quality problems are becoming more ubiquitous worldwide. One policy approach that found traction in California is an innovative variant on the traditional Total Maximum Daily Load (TMDL) approach to salinity regulation, which has features in common with a scheme in Australia’s Hunter River Basin. The paper describes the real-time salinity management (RTSM) concept, which is geared to improving coordination of west side agricultural and wetland exports of salt load with east side tributary reservoir release flows to improve compliance with river salinity objectives. RTSM is a concept that requires access to continuous flow and electrical conductivity data from sensor networks located along the San Joaquin River and its major tributaries and a simulation model-based decision support designed to make salt load assimilative capacity forecasts. Web-based information dissemination and data sharing innovations are described with an emphasis on experience with stakeholder engagement and participation. The last decade has seen wide-scale, global deployment of similar technologies for enhancing irrigation agriculture productivity and protecting environmental resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hydrological Control on Soil Redox Condition and Carbon Loss of Coastal Wetland Under Sea-Level Rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Chen, Kewei [ORNL] (ORCID:000000032580514X)↗
DOE OSTI2020

Characterizing Physical Properties of Streambed Interface Sediments Using In Situ Complex Electrical Conductivity Measurements

Streambed sediment physical properties such as surface area, are difficult to quantify in situ but exert a high-level control on a wide range of biogeochemical processes and sorption of contaminants. We introduce the use of complex electrical conductivity (CC) methods (also known as spectral-induced polarization (SIP)) that measure both real and imaginary conductivity to noninvasively and efficiently characterize shallow streambed sediments. We explore the method through synthetic modeling, laboratory, and field measurements to demonstrate the sensitivity of imaginary conductivity to sediment surface area, controlled in part by fine-grained iron oxides produced by anoxic groundwater discharge. Laboratory measurements verify expected relationships between CC parameters and sediment properties. Synthetic modeling using a 1D analytical model illustrates the influence of water layer depth and conductivity on the field CC measurements made at the streambed-stream water interface. Specifically, the inverted sediment imaginary conductivity is less impacted by uncertainty in the water layer depth and conductivity relative to the real conductivity and phase shift. Field CC measurements along a landfill-impacted river reveal discrete streambed zones with enhanced bulk surface area generally corresponding to anoxic groundwater discharge zones with high concentrations of fine-grained iron oxide precipitates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
NASA NTRS2017

Applications of Satellite Data to Support Improvements in Irrigation and Groundwater Management in California

In agricultural regions around the world, threats to water supplies from drought and groundwater depletion are driving increased demand for tools to advance agricultural water use efficiency and support sustainable groundwater management. Satellite mapping of evapotranspiration (ET) from irrigated agricultural lands can provide agricultural producers and water resource managers with information that can be used to both optimize ag water use and improve estimates of groundwater withdrawals for irrigation. We describe the development of two remote sensing-based tools for ET mapping in California, including important lessons in terms of system design, partnership development, and transition to operations. For irrigation management, the integration of satellite data and surface sensor networks to provide timely delivery of information on crop water requirements can make irrigation scheduling more practical, convenient, and accurate.Developed through a partnership between NASA and the CA Department of Water Resources, the Satellite Irrigation Management Support (SIMS) framework integrates satellite data with information from agricultural weather networks to map crop canopy development and crop water requirements at the scale of individual fields. Information is distributed to agricultural producers and water managers via a web-based interface and web data services. SIMS also provides an API that facilitates integration with other irrigation decision support tools, such as CropManage and IrriQuest. Field trials using these integrated tools have shown that they can be used to sustain yields while improving water use efficiency and nutrient management. For sustainable groundwater management, the combination of satellite-derived estimates of ET and data on surface water deliveries for irrigation can increase the accuracy of estimates of groundwater pumping. We are developing an OpenET platform to facilitate access to ET data from multiple models and accelerate operational use of ET data in support of a range of water management applications, including implementation of the Sustainable Groundwater Management Act in CA. By providing a shared basis for decision making, we anticipate that the OpenET platform will accelerate implementation of solutions for sustainable groundwater management.

Satellite↗
DOE OSTI2025

Distinct Fracture Mineralogy That is Out of Equilibrium With Modern Groundwaters Provides Important Context for Subsurface Life

Rock fracture surfaces in the crust are essential habitat for microorganisms. Fracture‐groundwater interfaces provide physical substrates for biofilm growth and are sources of carbon, nutrients, and electron donors and acceptors. To better understand geochemical processes impacting fracture surfaces and the subsurface microbiome, we identified fractures in archived rock cores from the Soudan formation, which is known to host saline groundwaters and isolated microbial communities dependent on rock‐water interactions. Cores with open fractures were thin sectioned and studied via electron microprobe and synchrotron X‐ray fluorescence microprobe. Most fracture surfaces had mineralogy distinct from that of the bulk rock. Chlorite minerals were abundant on fracture surfaces and had elemental compositions suggesting deposition during late‐stage hydrothermal alteration. Fracture‐lining chlorites likely limit access to iron oxide and sulfide minerals that are active in subsurface biogeochemical cycles. Calcium‐rich rinds were also observed along fracture edges. These rinds were too thin and poorly ordered to be identified via light microscopy or X‐ray diffraction; however, Ca K‐edge micro‐X‐ray absorption near‐edge structure spectroscopy identified them as carbonates, minerals not observed in the bulk rock. Thermodynamic modeling shows that carbonate precipitation is largely unfavorable in Soudan groundwaters, indicating that fracture edge conditions differed from those in modern water samples. Because of the low carbon concentrations in Soudan groundwaters, carbonate rinds likely play an important role in subsurface carbon cycling and may mark fracture surfaces that once hosted biofilms. Overall, this study suggests that fracture alteration can both play an active role in and suppress rock‐water interactions essential to subsurface life.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

3D characterisation and quantification of an offshore freshened groundwater system in the Canterbury Bight

Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Differentiable, Learnable, Regionalized Process-Based Models With Multiphysical Outputs can Approach State-Of-The-Art Hydrologic Prediction Accuracy

Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2018

Benefits and Pitfalls of GRACE Terrestrial Water Storage Data Assimilation

Satellite observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) mission have a coarse resolution in time (monthly) and space (roughly 150,000 sq km at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Nonetheless, data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This presentation illustrates some of the benefits and drawbacks of assimilating TWS observations from GRACE into a land surface model over the continental United States and India. The assimilation scheme yields improved skill metrics for groundwater compared to the no-assimilation simulations. A smaller impact is seen for surface and root-zone soil moisture. Further, GRACE observes TWS depletion associated with anthropogenic groundwater extraction. Results from the assimilation emphasize the importance of representing anthropogenic processes in land surface modeling and data assimilation systems.

TWS↗
DOE OSTI2022

Predictive Flow Simulation with the P2R Model for the 200-IA-1 Preliminary Remediation Goal Saturated Zone Abstraction

In order to develop and evaluate alternatives for remediation of waste sites constituting the 200-IA-1 Operable Unit (OU), a focused feasibility study (FFS) is currently underway. The 200-IA-1 OU FFS requires preliminary remediation goals (PRGs) protective of groundwater be evaluated which includes estimates of fate and transport of contaminants in the groundwater from multiple sources within the Central Plateau Inner Area. This document details the application of the Plateau-to-River (P2R) Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the 200-IA-1 OU PRG calculation. The simulated flow field will support the simulation of fate and transport of contaminants for use in estimating peak unit-length concentration as part of the 200-IA-1 OU PRG calculation presented in a separate environmental calculation file (ECF) (ECF-200IA1-20-0108, Determination of Preliminary Remediation Goals for 200-IA-1 Source Operable Unit).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗