Search NASA⌕ Search

SEARCH · Search NASA

Results for “geologic model”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Applications and Results of MODIS Lunar Observations

Both Terra and Aqua MODIS have been making periodic lunar observations since their launch in December 1999 and May 2002, respectively. The primary objective of the MODIS lunar observations is to monitor on-orbit response changes for the reflective solar bands (RSB) and to track long-term radiometric calibration stability. The Moon is an extremely stable surface reflectance reference and, together with the lunar radiometric model developed recently by the US Geological Survey (USGS), the applications of lunar observations for Earth remote sensing missions have been greatly enhanced and expanded. In particular, calibration inter-comparisons among sensors on different platforms and operated in different time frames have been enabled in the VIS/NIR/SWIR spectral regions. In this paper, we discuss various applications of MODIS lunar observations, such as monitoring the RSB radiometric stability and calibration consistency between the two sensors, evaluating the calibration differences among detectors within the same spectral band, and tracking the on-orbit changes of band-to-band registration (BBR). Other applications designed for the thermal emissive bands (TEB) and special sensor characterization purposes are also illustrated. Results from multi-year Terra and Aqua MODIS lunar observations demonstrate clearly that the two sensors have been calibrated consistently to within 1% for most RSB compared to the 2% radiometric calibration requirement. In addition the spatial characterization results derived from lunar observations are in good agreement with those derived from the instrument's on-board calibrator. The applications and results of MODIS lunar observations presented in this paper provide useful references for other sensors that also make use of lunar surface observations.

Xiong, X.↗
NASA NTRS2024

Combining Multi-Faceted Laboratory Studies of 74001-2 and Regional Remote Sensing to Address How Pyroclastic Eruptions Record and Affect the Lunar Volatile Budget

Basaltic magmatism is an efficient process for bringing volatiles from a planetary interior to its surface, with the possibility of generation of a transient lunar atmosphere as abundant volcanic materials de-gassed. However, pyroclastic deposits are locations where trapped gases may be studied [e.g., 2,3]. Volatile-rich pyroclastic deposits occur over a wide surface area of the Moon, indicating that the transport of volatiles and associated pyroclastic materials from the Moon’s mantle to the surface was a wide-spread phenomenon. Numerous studies analyzing remotely sensed data and using empirical modeling have demonstrated that various stages of pyroclastic eruptions significantly influence gas release patterns, morphology, and mineralogy of the deposit. Many observations based on mare basalts and pyroclastic deposits have identified potential histories of gas release [e.g., 2-10] and their influence on volatiles and their stable isotopes [e.g., 11-13]. The best representation of these pyroclastic deposits in the sample collection is core sample 74001-74002 that was collected during the Apollo 17 mission to the Taurus-Littrow Valley (TLV). The double drive tube penetrated a part of a regional-scale pyroclastic deposit and sampled approximately 68.1 cm of that deposit in the TLV. Remnants of this and other pyroclastic depos-its are represented throughout and beyond the TLV [e.g., 14-16]. The stratigraphy of this core has been investigated and defined by numerous studies. The CASA Moon SSERVI research team is conducting a multi-faceted analytical study of this deposit. Data generated from revisiting the stratigraphy of 74001-74002 will be used to place stable isotopes (H, B, Cl, S, Zn, Cu, Rb, Ga, Pb), Ar-Ar and U-Pb chronology, geochemistry, nanometer-scale observations of mineral surfaces, orbital observations, and experiments and modeling within a stratigraphic, eruptive, and geologic context. It is important to place these data into such a context. For example, recent S isotope measurements reported by Dottin et al. show differences within this core that may be related to either changes in source or eruptive process over the course of the eruption (vs. multiple eruptions). A fuller understanding of the stratigraphy is fundamental to resolving this interpretation. This comprehensive approach can only be achieved within the context of a program such as SSERVI. In addition, imaging produced in this project will be incorporated into a citizen scientist program to further identify many of the textural features of this double drive tube.

SSERVI↗
DOE OSTI2020

Investigation and Analysis of Dolomite Dissolution in Variable Strength Systems Relevant to the WIPP

The Waste Isolation Pilot Plant (WIPP) was opened in 1999 as a solution to the long-term disposal of the nation's nuclear defense-generated transuranic (TRU) waste. The salt beds were recommended as the ideal location because they are free from flowing fresh water, easily mined, impermeable, and geologically stable. However, for performance assessment modeling, release scenarios must be considered. The most likely release scenario was determined to be human intrusion leading to direct and/or long-term brine release (US DOE 1995, US DOE 1996, Perkins et al., 1999). The focus of this research is the Rustler formation because it is located above the WIPP and is the most transmissive layer (Perkins et al., 1999). Objective: To determine the effect of chelating ligands ( e.g. EDTA) on dissolution of dolomite in high ionic strength systems. Do brines impact dolomite dissolution and how is it affected by EDTA binding with cations? By varying the ionic strengths, different brine conditions can be simulated similar to those in and around the Culebra member and evaluate trends of dissolved Ca and Mg in the aqueous phase from the quantified data. Conclusions: <2% of dolomite dissolved over one week at pH 8.5. Greater dissolution occurred at higher ionic strength and in the presence of NaCl with approximately 0.4 vs. 0.9% of Mg removed from dolomite in 0.1 and 1.0 M Na. Further, 0.3 vs. 0.9% of Mg was removed from dolomite in 1.0 M IS CaCl{sub 2} versus NaCl. EDTA increased dissolution of dolomite at low ionic strength due to its strong complexation affinity for Ca/Mg with slightly greater removal of Ca than Mg. Significant differences were not observed at high ionic strength due to the presence of large concentrations of competing cations (Na, Ca, Mg). Slightly higher concentrations of Ca were removed from dolomite than Mg after one week, suggesting either incongruent dissolution or release through ion exchange of Ca from trace mineral impurities. Further, when EDTA was present, it may have increased Ca removal due to its greater complexation for Ca than Mg. Environmental Implications: High ionic strength brines and low ionic strength solutions with organic ligands increase dissolution of dolomite. Due to the long half-lives of radionuclides in the TRU waste, it may be important to consider the long term impact of dolomite dissolution for potential release scenarios.

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

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗
NASA NTRS1982

The geology of Europa

The map units and lineations of Europa are detailed, and the geologic processes, and history, and thick and thin ice models of the satellite are discussed. It is concluded that Europa lacks evidence of a horizontally stratified crust; the geology appears characterized by disruption of the crust and intrusions into an icy shell. The surface consists of plains and mottled terrain, the former being older. Numerous straight and curved lineations, streaks, stripes, and bands cross EuropA's surface on a global and surface scale. Most lineations appear related to fractures in the crust. Five fresh craters in the 10 to 30 km diameter range are visible. The dark spots, stripes, and bands that appear to have replaced sections of the crust suggest that material was transported to the surface from the subjacent silicate lithosphere. The apparent low density of craters superposed on Europa's surface suggests that the surface is about 100 million years old.

Lucchitta, B. K.↗
NASA NTRS1985

Impact of GRM: New evidence from the Soviet Union

Gravity information released by the Soviet Union allows the quantitative assessment of how the geopotential research mission (GRM) mission might effect the ability to use global gravity data for continental tectonic interpretation. The information is of an isostatic response spectra for eight individual tectonic units in the USSR. The regions examined include the Caroathians, Caucasus, Urals, Pamirs, Tien-Shan, Altal, Chersky Ridge, and East Siberian Platform. The 1 deg x 1 deg gravity data are used to calculate the admittances are used in two different sorts of tectonic studies of mountain belts in the USSR: (1) interpretation of isostatic responses in terms of plate models of compensation for mountainous terrain. Using geologic information concerning time of the orogeny, lithospheric plates involved, and polarity of subduction in collision zones, they convert the best-fitting flexural rigidity to an elastic plate thickness for the lithospheric plate inferred to underlie the mountains; the isostatic admittance functions is an attempt to directly model gravity and topography data for a few select regions in the Soviet Union. By knowing the value of the expected correlation between topography and gravity from the admittances, the Artemjev's map in mountainous areas can be calibrated, and the maps are converted back to Bouguer gravity. This procedure is applied to the Caucasus and southern Urals.

Mcnutt, M.↗
DOE OSTI2020

Optimal carbon storage reservoir management through deep reinforcement learning

Model-based optimization plays a central role in energy system design and management. The complexity and high-dimensionality of many process-level models, especially those used for geosystem energy exploration and utilization, often lead to formidable computational costs when the dimension of decision space is also large. This work adopts elements of recently advanced deep learning techniques to solve a sequential decision-making problem in applied geosystem management. Specifically, a deep reinforcement learning framework was formed for optimal multiperiod planning, in which a deep Q-learning network (DQN) agent was trained to maximize rewards by learning from high-dimensional inputs and from exploitation of its past experiences. To expedite computation, deep multitask learning was used to approximate high-dimensional, multistate transition functions. Both DQN and deep multitask learning are pattern based. As a demonstration, the framework was applied to optimal carbon sequestration reservoir planning using two different types of management strategies: monitoring only and brine extraction. Both strategies are designed to mitigate potential risks due to pressure buildup. Results show that the DQN agent can identify the optimal policies to maximize the reward for given risk and cost constraints. Finally, experiments also show that knowledge the agent gained from interacting with one environment is largely preserved when deploying the same agent in other similar environments.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Geoanalytical Evaluation of Saline Storage (GEESS) Geodatabase v2.0

The Geoanalytical Economic Evaluation of Saline Storage (GEESS) geodatabase was developed to support the United States Department of Energy (DOE) and National Energy Technology Laboratory (NETL) in their geologic carbon storage efforts by characterizing saline geologic formations present in the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM) [1]. Using publicly available literature and data, the GEESS geodatabase characterizes 57 geologic formations across the lower-48 U.S. states in what are called Fully Integrated Geodatabases (FIGs). The FIG is a vector polygon feature containing thousands or tens of thousands of individual polygons, which each contain discrete geologic parameter values. A list of the critical geologic parameters that are characterized in the GEESS geodatabase are described in the “Processing Steps and Workflow” part of the ReadMe file, as well as the Data Catalog accompanying the GEESS geodatabase. The FIG is the basis of the GEESS geodatabase and is the direct representation of the collected geologic data. In addition to the FIG, the GEESS system contains grid files. Due to the complexity of the FIGs, grids are used to sample the geologic data so they can be exercised within CO2_S_COM. The grid files contain the geologic data sampled from the FIG, as well as estimates of “Plume Uncertainty Diameter” and “First-year Break-even Price of CO2” derived from CO2_S_COM based on the GEESS grid data.

carbon↗
NASA NTRS1994

Venus climate stability and volcanic resurfacing rates

The climate of Venus is to a large degree controlled by the radiative properties of its massive atmosphere. In addition, outgassing due to volcanic activity, exospheric escape processes, and surface/atmosphere interactions may all be important in moderating the abundances of atmospheric CO2 and other volatiles. We have developed an evolutionary climate model for Venus using a systems approach that emphasizes feedbacks between elements in the climate system. Modules for atmospheric radiative transfer, surface/atmosphere interactions, tropospheric chemistry, and exospheric escape processes have so far been developed. Climate feedback loops result from interconnections between modules, in the form of the environmental parameters pressure, temperature, and atmospheric mixing ratios. The radiative transfer module has been implemented by using Rosseland mean opacities in a one dimensional grey radiative-convective model. The model has been solved for the static (time independent) case to determine climate equilibrium points. The dynamics of the model have also been explored by employing reaction/diffusion kinetics for possible surface atmosphere heterogeneous reactions over geologic timescales. It was found that under current conditions, the model predicts that the climate of Venus is at or near an unstable equilibrium point. The effects of constant rate volcanism and corresponding exsolution of volatiles on the stability of the climate model were also explored.

Bullock, M. A.↗
NASA NTRS2003

The Formation of Life-sustaining Planets in Extrasolar Systems

The spatial exploration is providing us a large quantity of information about the composition of the planets and satellites crusts. However, most of the experiences that are proposed in the guides of activities in Planetary Geology are based exclusively on the images utilization: photographs, maps, models or artistic reconstructions [1,2]. That things help us to recognize shapes and to deduce geological processes, but they says us little about the materials that they are implicated. In order to avoid this dicotomy between shapes and materials, we have designed an experience in the one which, employing of rocks and landscapes of our geological environment more next, the pupils be able to do an exercise of compared planetology analyzing shapes, processes and material of several planetary bodies of the Solar System.

Chambers, J. E.↗
DOE OSTI2026

Data for The Value of Reversible Carbon Storage in a Zero-Emissions World

Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.

Carbon↗
NASA NTRS1985

Dislocation model for aseismic fault slip in the transverse ranges of Southern California

Geodetic data at a plate boundary can reveal the pattern of subsurface displacements that accompany plate motion. These displacements are modelled as the sum of rigid block motion and the elastic effects of frictional interaction between blocks. The frictional interactions are represented by uniform dislocation on each of several rectangular fault patches. The block velocities and fault parameters are then estimated from geodetic data. Bayesian inversion procedure employs prior estimates based on geological and seismological data. The method is applied to the Transverse Ranges, using prior geological and seismological data and geodetic data from the USGS trilateration networks. Geodetic data imply a displacement rate of about 20 mm/yr across the San Andreas Fault, while the geologic estimates exceed 30 mm/yr. The prior model and the final estimates both imply about 10 mm/yr crustal shortening normal to the trend of the San Andreas Fault. Aseismic fault motion is a major contributor to plate motion. The geodetic data can help to identify faults that are suffering rapid stress accumulation; in the Transverse Ranges those faults are the San Andreas and the Santa Susana.

Cheng, A.↗
DOE OSTI2025

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗
DOE OSTI2022

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

The excavation stage of basin formation - A qualitative model

One of the most complex problems in planetary geology and geophysics is the determination of the nature of the impact cratering processes at scales of tens to thousands of kilometers that produce the complex morphological structures of multiring basins. The cratering process is frequently considered to be divided into three stages, including a short high-pressure stage of initial contact between the projectile and the planetary crust, a longer excavation or cratering flow stage culminating in the formation of a transient crater, and a still longer modification stage during which the transient crater is modified into the observed final geologic form. The transient crater may be considered as the initial boundary condition of the modification stage. In the present investigation, the nature of the transient crater is indicated by the cratering flow field determined from numerical simulations of the excavation stage. Attention is given to empirical and theoretical scaling.

Croft, S. K.↗
DOE OSTI2021

Vadose and Saturated Zone Flow and Transport Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

This model package report (MPR) documents the development of the integrated vadose and saturated zone flow and transport model developed for the performance assessment of Trenches 31 and 34 in the low-level burial grounds (LLBGs). This modeling capability is intended for use in addressing the analysis requirements outlined in DOE O 435.1, Chg 1, Radioactive Waste Management1. The overall objective of the modeling effort is to provide a basis for making informed disposal decisions pertinent to Trenches 31 and 34. The purpose of the MPR is to document the development of the three-dimensional numerical vadose zone and saturated flow and transport model test case and evaluate its adequacy to support the LLBG performance assessment. The purpose is not to present results for DOE 435.1 decision making. The use of the model to perform base case and sensitivity analysis for the LLBG performance assessment, including inputs and results, is documented separately in subsequent environmental calculation files. This report discusses the development and translation of the conceptual model for flow and contaminant transport into the LLBG performance assessment three-dimensional numerical flow and transport model evaluated using the Subsurface Transport Over Multiple Phases (STOMP©) simulator. The development of representative geologic framework is described along with the implementation of waste release models used to represent contaminant release from waste disposed in the trenches. The report also provides the technical basis for specific model parameters and boundary conditions, along with description of modeling assumptions. This MPR includes certain calculations that are necessary to demonstrate the soundness of the model. Results provided by the model include vadose zone and saturated zone flow fields, and estimates of the possible future concentration in groundwater of technetium-99 and iodine-129 as example test cases. As an evaluation of the test cases, the model estimates of current vadose conditions are compared to available and analogous field and laboratory data, and compared to the results from the original performance assessment model (WHC-EP-0645, Performance Assessment for the Disposal of Low-Level Waste in the 200 West Area Burial Grounds 2 ). The features, events, and processes applicable to vadose zone and saturated zone flow and transport model are identified, and representative initial estimates for various parameters are documented. Note that the parameter estimates presented in this MPR are for illustration purposes, and may or may not reflect values selected for eventual performance assessment. Several key topical discussions (i.e., basis for recharge estimates, basis for vadose zone modeling, basis for saturated zone model development, and calibration) are included, which serve as the groundwork for confidence building for the groundwater pathway modeling and results. Numerical simulation results based on an example test case are included to illustrate the use of the combined saturated-unsaturated model in performance assessment calculations. Sensitivity and uncertainty results are not included in this MPR. Those results will be included in future environmental calculation files. The inclusion of the trapezoidal trench geometry and construction details in the finite difference grid introduces some gross simplifications regarding the trench and liner systems. The model test case presented in this MPR includes the assumption that the trench liner system does not affect flow through or around the trenches after the liner system is assumed to fail. This MPR also includes results of other test cases that involve alternate assumptions about how to incorporate the hydraulic effects of the trenches into the vadose zone of the model. These alternate cases do not attempt to account for the presence of the liner system after its assumed failure either. None of these cases is considered to be the base case at this time. Analysis and alternate cases that attempt to account for the presence of the liner system in greater detail, and its effect on flow through and around it, are to be documented in subsequent environmental calculation files. Depending on the results of those analyses and alternate cases, the eventual base case may involve more detailed inclusion of the effects of the liner system hydraulics.

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

Well–based monitoring of CO 2 geological sequestration operations in saline aquifers: Critical insights into key questions

Geological carbon sequestration in saline aquifers is one of the most promising strategies to help mitigate emissions of CO 2 to the atmosphere. Significant challenges in ensuring the security of the sequestration process rest in the evolution and expansion of the CO 2 plume in the subsurface. The ability to track the movement of the injected CO 2 poses another challenge. Critical questions related to the integrity of the sequestration operations in saline aquifers relate to plume characteristics that can we monitor using well-based variables. We addressed this and related questions using an integrated modeling framework through a numerical investigation of carbon sequestration in saline aquifers during long-term and post-injection periods. This modeling paradigm incorporates the effect of structural, geological, and petrophysical characteristics. That way, we can account for critical physicochemical processes, rock-fluid interactions, and lithology dependencies. The well fluid variables investigated include fluid composition, pH, fluid density, and ion activity. We learned that fluid property analytics and diagnostics can be powerful tools to estimate the movement of CO 2 and its storage in different trapping mechanisms. These analytics can help optimize operational aspects and simplify reservoir-scale models while still reflecting the complex nature of the CO 2 interactions underground and offering insights into plume evolution.

58 GEOSCIENCES↗