Interplay between microorganisms and geochemistry in geological carbon storage
Explore the source record for details and available documents.
nuclear (including radiation effects), carbon sequ↗
SEARCH · Search NASA
Results for “Geological 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.
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↗
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.↗
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.↗
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↗
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.↗
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↗
Three-Dimensional Finite Difference Simulation of Ground Motions from the August 24, 2014 South Napa Earthquake
We performed three-dimensional (3D) anelastic ground motion simulations of the South Napa earthquake to investigate the performance of different finite rupture models and the effects of 3D structure on the observed wavefield. We considered rupture models reported by Dreger et al. (2015), Ji et al., (2015), Wei et al. (2015) and Melgar et al. (2015). We used the SW4 anelastic finite difference code developed at Lawrence Livermore National Laboratory (Petersson and Sjogreen, 2013) and distributed by the Computational Infrastructure for Geodynamics. This code can compute the seismic response for fully 3D sub-surface models, including surface topography and linear anelasticity. We use the 3D geologic/seismic model of the San Francisco Bay Area developed by the United States Geological Survey (Aagaard et al., 2008, 2010). Evaluation of earlier versions of this model indicated that the structure can reproduce main features of observed waveforms from moderate earthquakes (Rodgers et al., 2008; Kim et al., 2010). Simulations were performed for a domain covering local distances (< 25 km) and resolution providing simulated ground motions valid to 1 Hz.
58 GEOSCIENCES↗
The geophysics, geology and mechanics of slow fault slip
We report that modern geodetic and seismologic observations describe the behavior of fault slip over a vast range of spatial and temporal scales. Slip at sub-seismogenic speeds is evident from top to bottom of lithospheric faults and plays an important role throughout the earthquake cycle. Where earthquakes and tremor accompany slow slip, they help illuminate the spatiotemporal evolution of fault slip. Geophysical subsurface imaging and geologic field studies provide information about suitable environments of slow slip. In particular, exhumed fault and shear zones from various depths reveal the importance of multiple deformation processes and fault-zone structures. Most geologic examples feature frictionally weak and velocity-strengthening materials, well-developed mineral fabrics, and abundant veining indicative of near-lithostatic fluid pressure. To produce transient slow slip events and tremor, in addition to the presence of high-pressure fluids a heterogeneous fault-zone structure, composition, and/or metamorphic assemblage may be needed. Laboratory and computational models suggest that velocity-weakening slip patches smaller than a critical dimension needed for earthquake nucleation will also fail in slow slip events. Changes in fluid pressure or slip rate can cause a fault to transition between stable and unstable fault slip behavior. Lastly, future interdisciplinary investigations of slow fault slip, directly integrating geophysical, geological and modeling investigations, will further improve our understanding of the dynamics of slow slip and aid in providing more accurate earthquake hazard characterizations.
58 GEOSCIENCES↗
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↗
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↗
Impact of Three-Phase Relative Permeability and Hysteresis Models on Forecasts of Storage Associated With CO 2 -EOR
Geological CO 2 sequestration in conjunction with enhanced oil recovery (CO 2 -EOR) includes complex multiphase flow processes compared to CO 2 storage in deep saline aquifers. Two of the most critical factors affecting multiphase flow in CO 2 -EOR are three-phase relative permeability and associated hysteresis, both of which are difficult to measure and are usually represented by numerical interpolation models. The purpose of this research is to improve understanding of (1) the relative impacts of different three-phase relative permeability models and hysteresis models on CO 2 trapping mechanisms, and (2) uncertainty associated with these two factors. Four different three-phase relative permeability models and three hysteresis models were applied to simulations of an active CO 2 -EOR site, the SACROC unit located in western Texas. To eliminate potential bias of deterministic parameters, we utilized a sequential Gaussian simulation technique to generate 50 realizations to describe heterogeneity of porosity and permeability, based on data obtained from well logs and seismic survey. Simulation results of forecasted CO 2 storage suggested that (1) the choice of three-phase relative permeability model and hysteresis model led to noticeable impacts on forecasted CO 2 sequestration capacity; (2) impacts of three-phase relative permeability models and hysteresis models on CO 2 trapping are small during the CO 2 -EOR injection period, and increase during the post-EOR CO 2 injection period; (3) the specific choice of hysteresis model is more important relative to the choice of three-phase relative permeability model; and (4) using the recommended three-phase WAG (Water-Alternating-Gas) hysteresis model may increase the impact of three-phase relative permeability models and uncertainty due to heterogeneity.
42 ENGINEERING↗
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.↗
Modeling Background Radiation Using Geochemical Data
Aerial gamma ray surveys have many applications in geology and science in general, such as locating mining prospects, defining radioactive plumes, and detecting nuclear weapons. Unfortunately there is currently no simple way to separate the natural gamma radiation of soil and rocks from that of contaminants such as radioactive plumes. This project used geochemical data (uranium, potassium and thorium concentrations) collected from national databases, private companies, and the NURE (National Uranium Resource Evaluation) Survey, to create forward models of exposure rates measured by aerial gamma ray surveys. We developed these techniques using an area in north central Arizona known as the Navajo Mines area, chosen for its optimal conditions for aerial gamma ray surveys and readily available survey data. Models based on geochemical analyses from databases were not found to be successful, in part due to lack of data for some units. Models based on NURE data sorted by geologic unit were successful at replicating the aerial gamma ray survey, though units with Uranium mineralization proved difficult. ASTER visualizations were effectively used to create subunits of similar exposure rate within the Chinle Formation which contains multiple lithologies. For alluvial units, models based on drainage basin were attempted with success. With these models we are developing techniques to estimate the natural radiation generated by the rocks and soils of an area, making aerial ray surveys more effective.
58 GEOSCIENCES↗
Performance Assessment of a Generic Nuclear Waste Repository in Shale - 18392
A performance assessment (PA) is an important component of a comprehensive safety analysis for a nuclear waste repository. In a PA, probabilistic simulations of the total repository system are performed, and results are evaluated against performance metrics. Uncertainty and sensitivity analyses help prioritize additional research and model development. The United States Department of Energy (DOE) has been developing a state of the art PA simulation software tool-kit, Geologic Disposal Safety Assessment (GDSA) Framework, that couples increasingly higher fidelity models of subsystem processes into total system PA simulations. Over the past several years, PAs of generic repositories in three geologic media (salt, shale, and crystalline rock) have demonstrated ongoing developments in capability. The current PA of a nuclear waste repository in a generic shale formation showcases GDSA Framework, including capabilities in domain discretization (Cubit), multi-physics simulations (PFLOTRAN), uncertainty and sensitivity analysis (Dakota), and visualization (Paraview). The generic shale reference case considers the disposal of 22,000 metric tons heavy metal of commercial spent nuclear fuel (SNF) in a generic shale formation. PA simulations account for the thermal load and radionuclide inventory of the waste form, components of the engineered barrier system including waste package and buffer, and components of the natural barrier system including the host rock shale and underlying and overlying stratigraphic units. Two repository layouts are considered, one for emplacement of waste packages containing 12 pressurized water reactor (PWR) assemblies, and one for 4-PWR waste package emplacement. Model domains are half-symmetry and contain between 7 and 22 million grid cells. Grid refinement captures the detail of individual waste packages, emplacement drifts, access drifts, and shafts. Simulations are run in a high-performance computing (HPC) environment on 512 or 2048 processes. The governing equations describing coupled heat and fluid flow and reactive transport are solved with PFLOTRAN, an open-source, massively parallel multiphase flow and reactive transport code. Additional simulated processes include waste package degradation; waste form dissolution; radioactive decay and ingrowth in aqueous, solid, and adsorbed phases; and a simple biosphere based on a pumping well. Simulations are run to 106 y; dose and radionuclide concentrations are observed within aquifers at a point approximately 5 km downgradient of the repository. Dakota is used to sample likely ranges of input parameters including waste form and waste package degradation rates and properties of engineered and natural materials to quantify uncertainty in predicted concentrations and sensitivity to input parameters. Given the assumptions of the reference case, results of 12-PWR and 4-PWR simulations are very similar. Dose at the pumping well 5-km downgradient of the repository does not exceed 10-9 Sv/y. I-129 is the largest contributor to dose. The dose due to Cl-36, the next largest contributor, is orders of magnitude smaller. At the pumping well, the concentration of I-129 5-km downgradient does not exceed 10-14 mol/L. Rank correlation coefficients indicate that the maximum concentration of I-129 within the aquifer is sensitive to shale (repository host rock) porosity and aquifer permeability. These results affirm that HPC-capable codes can be used to simulate important multi-physics couplings directly in a total system performance assessment of a geologic repository, and can be used in prioritization of future research and model development. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
International Collaboration on Spent Fuel Disposition in Crystalline Media: FY17 Progress Report
Active participation in international R&D is crucial for achieving the Spent Fuel Waste Science & Technology (SFWST) long-term goals of conducting “experiments to fill data needs and confirm advanced modeling approaches” and of having a “robust modeling and experimental basis for evaluation of multiple disposal system options” (by 2020). DOE’s Office of Nuclear Energy (NE) has developed a strategic plan to advance cooperation with international partners. The international collaboration on the evaluation of crystalline disposal media at Sandia National Laboratories (SNL) in FY17 focused on the collaboration through the Development of Coupled Models and their Validation against Experiments (DECOVALEX-2019) project. The DECOVALEX project is an international research and model comparison collaboration, initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. SNL has been participating in three tasks of the DECOVALEX project: Task A. Modeling gas injection experiments (ENGINEER), Task C. Modeling groundwater recovery experiment in tunnel (GREET), and Task F. Fluid inclusion and movement in the tight rock (FINITO).
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Final Report on Update and Enhancement of Shale Gas Outlooks
The study performed by an interdisciplinary team of the Bureau of Economic Geology (Jackson School of Geosciences, The University of Texas at Austin) presents an improved production outlook assessments for the major U.S. shale gas plays, including the Barnett, Fayetteville, Haynesville, and Marcellus. The team of geoscientists, engineers, statisticians, and economists has reviewed and analyzed developments in the understanding of geologic, production-technology, and economic factors and has developed a model to run various scenarios with respect to energy prices, costs, technological advances, changes in geologic attributes, and physical and regulatory constraints.
03 NATURAL GAS↗
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↗