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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Chapman Conference on the Hydrologic Aspects of Global Climate Change, Lake Chelan, WA, June 12-14, 1990, Selected Papers

The present conference on the hydrological aspects of global climate change discusses land-surface schemes for future climate models, modeling of the land-surface boundary in climate models as a composite of independent vegetation, a land-surface hydrology parameterizaton with subgrid variability for general circulation models, and conceptual aspects of a statistical-dynamical approach to represent landscape subgrid-scale heterogeneities in atmospheric models. Attention is given to the impact of global warming on river runoff, the influence of atmospheric moisture transport on the fresh water balance of the Atlantic drainage basin, a comparison of observations and model simulations of tropospheric water vapor, and the use of weather types to disaggregate the prediction of general circulation models. Topics addressed include the potential response of an Arctic watershed during a period of global warming and the sensitivity of groundwater recharge estimates to climate variability and change.

Lettenmaier, Dennis P.↗

NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗

NASA NTRS1998

Magmatic Intrusions and a Hydrothermal Origin for Fluvial Valleys on Mars

Numerical models of Martian hydrothermal systems demonstrate that systems associated with magmatic intrusions greater than several hundred cubic kilometers can provide sufficient groundwater outflow to form the observed fluvial valleys, if subsurface permeability exceeds about 1.0 darcy. Groundwater outflow increases with increasing intrusion volume and subsurface permeability and is relatively insensitive to intrusion depth and subsurface porosity within the range considered here. Hydrothermally-derived fluids can melt through 1 to 2 km thick ice-rich permafrost layers in several thousand years. Hydrothermal systems thus provide a viable alternative to rainfall for providing surface water for valley formation. This mechanism can form fluvial valleys not only during the postulated early warm, wet climatic epoch, but also during more recent epochs when atmospheric conditions did not favor atmospheric cycling of water. The clustered distribution of the valley networks on a given geologic surface or terrain unit of Mars may also be more compatible with localized, hydrothermally-driven groundwater outflow than regional rainfall. Hydrothermal centers on Mars may have provided appropriate environments for the initiation of life or final oases for the long-term persistence of life.

Gulick, Virginia C↗

DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗

DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Exploration of groundwater and geothermal characteristics of Tohatchi Hot Springs aquifer at local and regional scales

Objective: Aquifer sustainability for geothermal resources exploration of Tohatchi Hotsprings, New Mexico. Tasks: (1) Generate 3D numerical models using PFLOTRAN to explore the geothermal resources in the Tohatchi spring aquifer and its contiguous areas and (2) estimate the sustainability of the aquifer for future geothermal resource development.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1976

International Conference on Environmental Sensing and Assessment, Las Vegas, Nev., September 14-19, 1975, Proceedings. Volumes 1 & 2

The papers deal with the detection of hazardous environmental pollutants, the development of emission control plans, and the design of compliance monitoring systems. Topics include remote sensing techniques in environmental pollution monitoring, monitoring of atmospheric particulate matter, air pollution due to sulfur dioxide and other inorganic compounds, marine pollution, atmospheric aerosols, industrial pollution, and legal aspects of pollution monitoring. Other papers examine the toxic effects of heavy metals and halogenated hydrocarbons, pollution associated with waste-disposal processes, pesticide residues in soil and groundwater, evaluations of groundwater quality, and monitoring of nuclear wastes. The interaction of climate and pollution is also discussed along with global pollutant transport, environmental modeling, ambient environmental air quality, aircraft and ground-vehicle emissions, and pollution associated with energy extraction and utilization processes. Individual items are announced in this issue.

Source record↗

NASA NTRS2011

Use of GRACE Terrestrial Water Storage Retrievals to Evaluate Model Estimates by the Australian Water Resources Assessment System

Terrestrial water storage (TWS) estimates retrievals from the Gravity Recovery and Climate Experiment (GRACE) satellite mission were compared to TWS modeled by the Australian Water Resources Assessment (AWRA) system. The aim was to test whether differences could be attributed and used to identify model deficiencies. Data for 2003 2010 were decomposed into the seasonal cycle, linear trends and the remaining de-trended anomalies before comparing. AWRA tended to have smaller seasonal amplitude than GRACE. GRACE showed a strong (greater than 15 millimeter per year) drying trend in northwest Australia that was associated with a preceding period of unusually wet conditions, whereas weaker drying trends in the southern Murray Basin and southwest Western Australia were associated with relatively dry conditions. AWRA estimated trends were less negative for these regions, while a more positive trend was estimated for areas affected by cyclone Charlotte in 2009. For 2003-2009, a decrease of 7-8 millimeter per year (50-60 cubic kilometers per year) was estimated from GRACE, enough to explain 6-7% of the contemporary rate of global sea level rise. This trend was not reproduced by the model. Agreement between model and data suggested that the GRACE retrieval error estimates are biased high. A scaling coefficient applied to GRACE TWS to reduce the effect of signal leakage appeared to degrade quantitative agreement for some regions. Model aspects identified for improvement included a need for better estimation of rainfall in northwest Australia, and more sophisticated treatment of diffuse groundwater discharge processes and surface-groundwater connectivity for some regions.

van Dijk, A. I. J. M.↗

DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗

DOE OSTI2026

Corrosion of 316 SS in direct contact with Zircaloy-4 in repository relevant conditions

The effects of interfacial contact between 316 stainless steel (316 SS) waste canister and spent fuel cladding material (Zircaloy-4) in case of water intrusion into a deep geological repository are investigated. Hydrothermal corrosion experiments of 316 SS and Zircaloy-4 in direct contact were performed at elevated temperatures in MilliQ water and simulated Mont Terri groundwater. Hematite, magnetite, chromite, and baddeleyite form as secondary phases. pH of the Mont Terri solution increases with corrosion duration and magnetite forms preferentially instead of hematite. Thermodynamic modeling shows that iron oxidation in 316 SS controls the pH and redox potential of the aqueous environment and is in good agreement with experimental Mont Terri results, showing that experimental Mont Terri results trend to equilibrium. All oxides except hematite are predicted to exist at thermodynamic equilibrium based on calculations, implying that experimentally formed hematite will be reduced or transformed to magnetite at equilibrium. The thermodynamically predicted chromite was observed to form only at the interface where 316 SS was in contact with Zircaloy-4 at 250°C, showing that contact with Zircaloy-4 seems to accelerate the 316 SS corrosion to thermodynamic equilibrium. Finally, Mont Terri groundwater readily trends towards thermodynamic equilibrium, leading to a pH and redox buffered system that is dominated by iron in solution coming from the accelerated corrosion of 316 SS caused by Zircaloy-4 contact.

316 stainless steel↗

DOE OSTI2022

Aggregation in bottom-up vulnerability assessments and equity implications: The case of Jordanian households’ water supply

Bottom-up methods for water resources modeling rely on acceptability thresholds to find, through a response surface, which deeply uncertain futures lead to system failure. They commonly treat water users as aggregate actors, which may preclude analysis of the equity impacts of interventions. This paper explores how aggregation choices for large groups of water users lead to different policy recommendations in response surface assessments. Herein, two aggregation methods with varying parameters are considered: percentile satisfaction targets and generalized mean. A 2-dimensional stress-test assessment across groundwater availability and population is applied to household water supply in Jordan. The study compares six different policies covering supply enhancement and rebalancing, using a country-wide multi-agent model that characterizes households across socioeconomic strata. For different aggregation levels, policies are ordered by their associated robustness index. Results show that aggregation choices may modify response surfaces as much as policy changes and strongly determine policy preference. Modifying allocation rules can substantially reduce the disparity in household vulnerability. Preferences defined by aggregation intervals provide a finer understanding of trade-offs among water users and may improve deliberation over equity under deep uncertainty.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

NEWTS USGS Brackish Water Database

Data from the USGS Brackish Water Database on brackish water compositions from across the US. Data has been filtered for high salinity and completeness of data attributes. Data has been transposed and converted for ease of input into aqueous chemistry modeling software including OLI Systems and Geochemist's Workbench. Original data from: Qi, S.L., and Harris, A.C., 2017, Geochemical Database for the Brackish Groundwater Assessment of the United States: U.S. Geological Survey data release, https://doi.org/10.5066/F72F7KK1.

Brackish Groundwater↗

DOE OSTI2021

PORFLOW Modeling of Vadose Zone Flow and Transport for the E- Area Low Activity Waste Vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Low Activity Waste Vault (LAWV) and local vadose zone has been developed using the PORFLOW TM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the LAWV 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 accessed by members of the public following site closure.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

PORFLOW Modeling of Vadose Zone Flow and Transport for 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 PORFLOW TM 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 is accessible to members of the public following site closure.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Alternative Conceptual Models of the Subsurface at the Hanford Site

This report describes the use of a multipoint geostatistics (MPS) framework with multiple data types for developing alternative conceptual-mathematical models of the subsurface at the Hanford Site. Developing multiple conceptual-mathematical models of the site is important for assessing potential uncertainties related to the features, events, and processes that control the subsurface flow of groundwater and the transport of contaminants of concern (COCs). Evaluating the effects of these uncertainties can result in better site characterization, monitoring, and remedial decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Modelling yield cap evolution in sandstone based on brittle creep experiments

The Zenifim Formation is being considered as a potential disposal formation for a deep borehole nuclear repository concept in Israel. Site selection and repository construction are intended to ensure that waste is separated from circulating groundwater, but long-term deformation of the wellbore could potentially create fluid flow pathways. To understand how time-dependent rock strength could affect wellbore stability, here we conducted creep tests under low to moderate confining pressures on retrieved core from the Zenifim formation. During creep, samples strain slowly as gradual damage accumulation progressively weakens the samples. Failure eventually occurred through the near-instantaneous formation of a shear fracture. Experimental results were used to calibrate a continuum damage poro-elastic model for sandstones. The calibrated damage-poro-elastic model successfully simulates different types of loading experiments including quasi-static and creep. The state of strain in experiments is close to yield during loading as the yield cap continuously evolves with damage accumulation. For creep tests, most damage occurs during triaxial loading. Minor damage accumulation occurs under constant load until the final stage of creep, where damage accelerates and promotes unstable fracturing.

58 GEOSCIENCES↗