Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vadose zone”

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 19 records

Analyzing Recharge Dynamics and Storage in a Thick, Karstic Vadose Zone

Karst vadose zone heterogeneity creates complex transmission and storage dynamics that affect the timing and magnitude of aquifer recharge. Young, high-matrix permeability eogenetic karst aquifers may have significantly higher matrix storage than older, lower matrix permeability counterparts. In vulnerable and water-limited karst regions, the timescales of storage may be important for seasonal and sub-seasonal water resource management. We create a framework to quantify storage dynamics in karst aquifers using high-resolution precipitation and groundwater levels from the Northern Guam Lens Aquifer in the US territory of Guam. We estimate recharge using the Water Table Fluctuation method, and then develop transfer functions between precipitation and recharge to quantify storage and release of water from the vadose zone. The transfer functions are partitioned into different flow pathways including conduit, conduit/matrix, and slowly draining matrix. Probability distributions are fit to each pathway to determine the average travel times of infiltrated waters. The results show that aquifer recharge through secondary porosity features typically occurs within a few hours of a rainfall event, but this rapid recharge accounts for only 12%–28% of total recharge. The majority of aquifer recharge (>70%) occurs within a month of a contributing storm event. An additional 10% of recharge, on average, took longer than a month to reach the water table. The framework established can help improve the hydrological modeling and freshwater management for karst aquifers.

Patricia Spellman↗

Composite Analysis Vadose Zone Transport of Selected Hanford Site Composite Analysis: Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Waste Sites

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Vadose Zone Model for SALDS for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP) is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the State-Approved Land Disposal Site (SALDS) model. The scope of this ECF is to document the development and results of the SALDS vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Model for B-3A/B Pond for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP) simulator officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the B-3A/B Ponds model. The scope of this ECF is to document the development and results of the B-3 A/B Ponds vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Transport of Tritium and Nitrate under Ponded Water Conditions

Vadose zone transport of tritium and nitrate can be important considerations at radioactive waste sites, landfills, or areas with industrial impacts. These contaminants are of particular concern because they typically have a relatively higher mobility in the subsurface compared to other compounds. Here, we describe a semiarid site with tritium and nitrate contamination involving a manmade ponded water source above a thick unsaturated zone at Los Alamos National Laboratory in New Mexico, USA. This study demonstrates the value of vadose zone flow and transport modeling for the development of field investigation plans (i.e., identifying optimal borehole locations and depths for contaminant characterization), and how a combination of modeling with isotope and geochemical measurements can provide insight into how tritium and nitrate transport in the vadose zone in semiarid environments. Modeling results suggest that at this location, tritium transport is well predicted by classical multiphase theory. Our work expands the demonstrated usefulness of a standard tritium conceptual model to sites with ponded surface conditions and agrees with previous results where a standard model was able to explain the evolution of a tritium plume at an arid waste disposal site. In addition, depth-based analyses of δ 18 O and δ 2 H of pore waters helped confirm the extent of pond infiltration into the vadose zone, and the δ 15 N of nitrate showed that the contaminant release history of the site was different than originally assumed.

58 GEOSCIENCES↗

Field Test Bed for Vadose-Zone Monitoring Approaches - 20404

In the Central Plateau at the U.S. Department of Energy Hanford Site, a large inventory of contaminants resides in unsaturated sediments within the approximately 100-meter-thick vadose zone, posing a potential continuing risk to groundwater. Vadose zone remedies used to address these contaminants will require performance monitoring to provide feedback during implementation and for long-term verification that remedial action objectives have been met. Passive approaches may also need long-term monitoring to demonstrate that the flux of contaminants from the vadose zone to the groundwater are below thresholds established to meet groundwater protection goals. Collection of physical (e.g., groundwater or sediment) samples is a common method for identifying contaminant concentration distributions, but this approach is limited by the number of locations and the frequency with which data can be collected. In situ vadose zone measurements have evolved over the past few years to include key measurements of water content, soil water pressure, temperature, and chemical concentration. However, the current generation of sensors is designed for relatively short-term use in near-surface soils or sediments. Geophysical methods have been evolving but are also limited in that they have not been designed for the specific long-term vadose zone monitoring needs at the Hanford Site. Overall, monitoring under unsaturated conditions can be difficult due to the need to install and maintain instrumentation over a large area and depth and the need to identify preferential flow pathways due to geologic and chemical heterogeneities over long periods. Thus, a vadose zone monitoring test bed was initiated to address these challenges and identify cost-effective approaches for implementation and postclosure monitoring of the deep vadose zone. The monitoring test bed is expected to provide valuable field-scale information for the design of vadose zone monitoring systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Central Plateau Vadose Zone Geoframework

The objective of this environmental calculation file (ECF) is to document the revision of the three dimensional (3D) geoframework model (GFM) of the Central Plateau Vadose Zone (CPVZ) on the Hanford Site. CP-60925, Model Package Report: Central Plateau Vadose Zone Geoframework, documents the original version of the CPVZGFM (herein after called the CPVZ-Rev.0) and is the underlying basis for the CPVZ conceptual model. The CPVZ GFM provides a more detailed interpretation of the vadose zone hydrostratigraphic units (HSUs) in the Central Plateau area of the greater Hanford South geoframework model (HSGFM) (ECF-HANFORD-13-0029, Development of the Hanford South Geologic Framework Model, Hanford Site, Washington). The development of the CPVZ GFM provides a consistent vadose zone geologic framework that is needed to support the 3D, multiple-waste-site modeling approach proposed in CP-60405, Hanford Site Composite Analysis Technical Approach Description: Vadose Zone, while incorporating an appropriate level of geologic data and interpretations for the size of the model using current site-standard HSUs. Future more detailed site-specific investigations within the Central Plateau can then use the CPVZ GFM as a starting point instead of the larger scale HSGFM. This GFM was developed using the site standard and approved for use IHS Kingdom software. This ECF is a revision to the original model described in CP-60925, but the underlying assumptions and conceptual model are unchanged. Therefore, this revision is a new calculation of version 1.0 of the CPVZ GFM and not a new version of the GFM. Throughout this report, the original model documented in CP-60925 is called CPVZ-Rev.0and the calculation for this revision of the model is called CPVZ-Rev.1.

58 GEOSCIENCES↗

Estimating Vadose Zone Flow Properties at the 100 K-East Soil Flushing Site Using ERT Monitoring Data: 2023 Interim Report - 100 KE Soil Property Estimation

In situ soil flushing is being using at the Hanford 100 K-East (100 KE) area to transport mobile chromium contamination in the vadose zone to the water table, where it can be collected and treat through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. In practice, it is infeasible to comprehensively monitor which regions of the vadose zone are being infiltrated through direct sampling of pore water. Consequently, there can be significant uncertainty about which regions of the subsurface have been treated, especially if hydrogeologic conditions are favorable for the development of unstable flows and preferred flow pathways through the vadose zone (Jarvis, Koestel, and Larsbo 2016). Current approaches for quantitative monitoring of soil flushing performance rely on contaminant concentration measurements collected from extractions wells. There is no quantitative information on the volume of flush water delivered to targeted regions of the vadose zone at the Hanford Site, leading to significant uncertainty regarding source term removal and long-term impacts to groundwater. If the subsurface hydrogeologic properties at the 100 KE Area were adequately known, qualitative metrics of soil flushing performance could be simulated, thereby negating expenses required to obtain quantitative performance information through borehole drilling/sampling. However, estimating in situ hydrogeologic properties has long proven elusive, due primarily to a lack of sufficient information to constrain heterogeneous property estimates to a useful degree of certainty. Estimating vadose zone hydrogeologic properties is particularly challenging due the dependence of hydraulic conductivity on saturation. This report describes progress toward a first-of-its-kind demonstration using surface time-lapse 3D electrical resistivity tomography (ERT) monitoring data to estimate the hydrogeologic properties that control flush water transport at the 100 KE soil flushing site. The ultimate objective is (1) to verify sufficient information exists in the ERT monitoring data to adequately resolve vadose zone hydraulic properties, and (2) generate a “digital twin” (i.e., a numerical simulator) that can be used to simulate the amount of flush water that has been delivered to each targeted region of the vadose zone, and thereby assess the efficacy of flush water delivery. Resulting performance estimates can be used in leu of comprehensive borehole drilling and direct sampling (or wellbore logging) that would otherwise be required to obtain the same information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Using Strain in the Vadose Zone During Sinusoidal-Rate Well Tests to Characterize an Aquifer

Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.

well test, sinusoidal, periodic pumping, strain↗

Composite Analysis Solid Waste Release Data Reduction of Activity Flux from Waste Sites to the Vadose Zone for Baseline Assessment

This environmental calculation file (ECF) describes the data reduction methodology for the solid waste release datasets generated by the Composite Analysis Solid Waste Release (CASWR) model (CP-62766, Model Package Report: Composite Analysis Solid Waste Release Model [CASWR Model]) and documents the resulting reduced datasets for the updated Composite Analysis (CA) baseline assessment. This ECF is in support of the updated Hanford Site CA. The purpose of the updated Hanford Site CA is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The vadose zone modeling activities associated with the CA includes solid waste forms as contaminant sources. To support the modeling of solid waste transport from the vadose zone to the groundwater at the Hanford Site, the CASWR model was designed to generate deterministic radionuclide release rates for the Hanford Site Central Plateau solid waste disposal sites. The complete CASWR model-generated dataset for the CA Baseline Assessment consists of 2,378,831 sets of time and release rates (1,073,869 sets for the 200 East Area and 1,304,962 sets for the 200 West Area), representing a course of the simulation from year 1953 to year 12,069, 16 radiological contaminants of potential concern (COPCs), and 237 solid waste sources as documented in ECFHANFORD- 19-0112, Solid Waste Release Calculations for the Composite Analysis Baseline Assessment. For several of the STOMP vadose zone models, the number of release rate data pairs required for release at all waste sites would exceed the limitations imposed by the STOMP modelling software package used to simulate contaminant transport through the vadose zone. As a result, a data reduction of the CASWR dataset is required prior to its use as a source input for the CA vadose zone modeling activities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Rapid Decision Support Tool for Estimating Impacts of a Vadose Zone Volatile Organic Compound Source on Groundwater and Soil Gas

Diminishing rates of subsurface volatile contaminant removal by soil vapor extraction (SVE) oftentimes warrants an in‐depth performance assessment to guide remedy decision‐making processes. Such a performance assessment must include quantitative approaches to better understand the impact of remaining vadose zone contamination on soil gas and groundwater concentrations. The spreadsheet‐based Soil Vapor Extraction Endstate Tool (SVEET) software functionality has recently been expanded to facilitate quantitative performance assessments. The updated version, referred to as SVEET2, includes expansion of the input parameter ranges for describing a site (site geometry, source characteristics, etc.), an expanded list of contaminants, and incorporation of elements of the Vapor Intrusion Estimation Tool for Unsaturated‐zone Sources software to provide soil gas concentration estimates for use in vapor intrusion evaluation. As part of the update, SVEET2 was used to estimate the impact of a tetrachloroethene (PCE) vadose zone source on groundwater concentrations, comparing SVEET2 results to field‐observed values at an undisclosed site where SVE was recently terminated. PCE concentrations from three separate monitoring wells were estimated by SVEET2 to be within the range of 6.0–6.7 μg/L, as compared to actual field concentrations that ranged from 3 to 11 μg/L PCE. These data demonstrate that SVEET2 can rapidly provide representative quantitative estimates of impacts from a vadose zone contaminant source at field sites. Finally, in the context of the SVE performance assessment, such quantitative estimates provide a basis to support remedial and/or regulatory decisions regarding the continued need for vadose zone volatile organic compound remediation or technical justification for SVE termination, which can significantly reduce the cost to complete for a site.

54 ENVIRONMENTAL SCIENCES↗

Multisystem feedbacks from a changing climate: Do altered hydrological dynamics control vadose zone carbon nutrient cycling and storage in shallow aquifer systems?

The vadose zone, which extends from upper soils to the subsurface water table, consists of many distinct habitats (including the critical zone), each with its own physical characteristics. Upper soils are typically richer in organic carbon chemical diversity and concentration, while deeper portions near the water table have less labile carbon and a greater percentage of humic acids and other long-lived organics. The availability of carbon and oxygen constrain the habitability of these zones. Typically, microorganisms (bacteria, archaea, and fungi) extend throughout the vadose zone and potentially deeper into the bedrock, while higher eukaryotes (i.e., arthropods and plants) are limited to the surficial soils. An exception to this is deep taproots of some tree species that can extend tens of meters into the subsurface. In subsurface systems, microbial metabolisms are constrained by the availability of carbon (organic and inorganic) and electron acceptors.

54 ENVIRONMENTAL SCIENCES↗

Use of an Alternative Conceptual Model of Vadose Zone Heterogeneity to Evaluate Past Tank Leaks and Other Unplanned Releases within a Tank Farm at the Hanford Site - 20061

Washington State Department of Ecology (Ecology) requested that U.S. Department of Energy (DOE)-Office of River Protection (ORP) consider an evaluation of effects of fine-grained thin sediment layers on transport with a separate alternative conceptual model in its evaluation of the potential impact of tank leaks and other unplanned releases within the tank farm at the Hanford Site. For this alternative model, Ecology recommended that the model be developed based on the general framework of fine-grained units identified by a stakeholder group in their interpretations of variability in moisture content data collected in the vicinity of Waste Management Area (WMA) C. In discussions with DOE-ORP, Ecology acknowledged that the underlying data and interpretations of the occurrence and lateral continuity of the fine-grained thin layers identified by a stakeholder group are uncertain. However, DOE-ORP agreed to Ecology's recommendation and has provided support for the requested evaluation that involved development of an alternative model based on the general framework of the unpublished report by the stakeholder. DOE-ORP considered this evaluation to be a hypothetical evaluation of vadose zone heterogeneities at WMA C. General observations from the range of simulation cases examined in the evaluation of the effects of hypothetical vadose zone heterogeneities at WMA C are as follows. - The movement of the center of mass of the simulated plumes was generally vertically downward below the source for all simulations, including those that incorporated the hypothetical fine-grained units. - All simulations that incorporated hypothetical heterogeneity produced additional plume spreading over what was produced in simulations using Equivalent Homogeneous Media (EHM) model(s)a. The spreading resulted in a broadening of the fringes of the plume, resulting in a wider region of low concentration, but lower peak concentrations associated with the center of mass of the plume. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced similar spreading and slightly lower peak mass flux at the water table when compared to the EHM modeling results. a The EHM-based models used at WMA C do not explicitly include small-scale fine-grained heterogeneities used in the model advocated by Ecology. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced less spreading and an earlier arrival of mass flux at the water table when compared to the use of another set of hydraulic properties from a silty-sand sample collected at a nearby disposal facility. - The EHM representation of the vadose zone generally produced higher peak mass flux and an earlier occurrence of peak fluxes at the water table compared to all analyses incorporating additional hypothetical heterogeneity. Results of this alternative model evaluation of past leaks provided some insight into the transport effects of vadose zone heterogeneities at WMA C that helped resolve Ecology's comments and issues related to the effects of vadose zone heterogeneities on past leaks and losses from the WMA C tank farms area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Water reactive polyurethane grouting for deep vadose zone contaminant immobilization

Contaminants in the deep vadose zone (DVZ) pose a long-term threat to groundwater, human health, and the environment. Polyurethane grouting is a type of chemical grouting that is potentially advantageous in immobilizing contaminants in DVZs. Polyurethane resin has low viscosity that makes it feasible to penetrate and fill the pore space between fine particles. Here, in this study a series of laboratory pressure grouting and leaching tests were conducted to simulate and assess the effectiveness of polyurethane grouting for immobilizing contaminants in DVZs. Soil with fine particles was prepared with 127 I (as iodide) that served as a non-radioactive surrogate for radioactive 129 I. After grouting and curing, leaching tests were used to measure and compare contaminant diffusivity and leachability index values. Additionally, changes in porosity and saturated hydraulic conductivity of grouted soil were measured. X-Ray Computer Tomography (XCT) results showed that the cured polyurethane was distributed nearly homogeneously and approximately half of the voids were filled with cured polyurethane. Grouting reduced the saturated hydraulic conductivity of soil by 37 %. The effective diffusivity decreased by more than 80 % as compared with the ungrouted soil. The leachability index of the grouted soil was 6.5; meeting the criteria established by the U.S. Nuclear Regulatory Commission (NRC) standard. The results obtained in this study provide a valuable assessment of polyurethane grouting for iodide immobilization in the DVZ and indicate this approach may be a viable method for contaminant remediation in DVZ soils.

Contaminant Diffusion↗

SX Tank Farm Interim Surface Barrier Vadose Zone Monitoring

The Hanford Site, in southeastern Washington State, has 149 underground single-shell tanks to store hazardous radioactive waste. A number of these tanks and their associated infrastructure (e.g., pipelines, diversion boxes) have leaked. While some of the leaked waste has entered the groundwater, the majority still resides in the vadose zone. Corrective measures are being taken to mitigate the groundwater impacts from past leaks and spills from contaminants that still reside within the vadose zone. Washington River Protection Solutions, LLC, constructed two interim surface barriers (identified as south and north ISBs) in 2018 over the southern portion of the 241-SX Tank Farm to minimize the migration of the residual contaminants in the vadose zone. Two monitoring instrument nests were installed in 2018 and the data acquisition system was put into function in September 2019. This document summarizes the monitoring methods, the instrument calibration and installation, and the vadose zone monitoring plan for the interim barriers in the SX tank farm.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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

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

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of a vadose zone advanced monitoring system: Tools to assess groundwater vulnerability

Performing repeat pore-fluid sampling over long time-scales can provide valuable information on unsaturated zone contaminants and their potential flux to ground water. This information can be used to manage groundwater remedies and identify contaminants that need to be sequestered in the vadose zone to minimize flux to ground water. Pore-water samples are commonly used to obtain contaminant concentrations within the vadose zone, but existing methods are limited as they only provide a single sample at one location and time. The vadose zone advanced monitoring system (VZAMS) has been designed to integrate multiple technologies into a single down-borehole system that allows for sampling of pore fluids (liquid and gas) to provide information about contamination and hydraulic conditions at multiple depths (~0.3-m intervals) within a cased borehole. Testing has been completed at the laboratory scale to verify the sampling elements of VZAMS, including geochemical testing for representative contaminants known to exist at the Hanford Site, located in southeastern Washington State. Physical tests focused on the ability of the sampler to draw fluid under unsaturated conditions. Initial geochemical testing showed that the stainless steel material used with the porous cuff may affect the sampled concentrations of redox-sensitive contaminants under very dry conditions. Additional laboratory testing demonstrated that the VZAMS components are able to collect representative samples for substances of interest under expected field conditions. In this paper, the design and functionality of a novel instrument are demonstrated in support of subsequent testing in the field.

54 ENVIRONMENTAL SCIENCES↗