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 55 records · Page 3

Deep Vadose Zone Monitoring Test Bed (FY23 Status Report)

A subsurface air injection at the Hanford Site’s Deep Vadose Zone Monitoring Test Bed was completed to realize a change in subsurface hydrologic conditions in accordance with a soil desiccation remedy. The injection mimicked a previous injection at the site that was performed in accordance with a vadose zone treatability study. Unlike the previous test which relied on electrical methods only, the change in hydrologic conditions during the recent test was also monitored using cross-hole seismic sensing methods to assess the ability of the seismic methods to evaluate changes in moisture conditions of desiccated sediment at the field scale. Data from in situ neutron probes indicates a reduction of soil moisture in the vicinity of the injection well due to the air injection. Similar changes were observed in the time-lapse electrical resistivity and seismic data, which indicates a loss of soil moisture over time. Tomographic inversions of the time-lapse geophysical data illustrate the 2D and 3D features of the soil moisture distribution over time. Time-lapse electrical resistivity tomography (ERT) results show a reduction in the electrical conductivity of the subsurface in the vicinity of the injection well, with most changes occurring within the screened interval. Similar patterns are observed in the seismic tomography results, with both methods illustrating two lobe-shaped features of reduced soil moisture. Use of seismic and ERT technologies in tandem takes advantage of two complementary geophysical monitoring technologies, providing increased sensitivity to specific hydrologic conditions. The multiphysics approach, therefore, has the potential to improve the ability to estimate subsurface moisture conditions from sensor-based and remotely sensed geophysical data that will ultimately improve the ability of remediation contractors to evaluate remedy performance.

58 GEOSCIENCES↗

Hanford Site Composite Analysis: LLBG-200-W B Vadose Zone Model

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 to the CA saturated zone model (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 1 ) 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 saturated zone model. The rates will be summed over the 100 by 100 m saturated zone model grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Flow and Transport Parameters Data Package for the Hanford Site Composite Analysis

This report provides a description of the basis for the development and implementation of a conceptual model for vadose zone flow and transport for the Hanford Site composite analysis (CA) groundwater pathway analysis. The parameterization for a numerical model is intimately linked to the conceptual model framework. This report describes the basis for the selection of hydraulic and transport parameters for the hydrostratigraphic units (HSUs) identified in the 200 East and 200 West Areas. The following three-step process was used to develop the hydraulic properties and transport parameters for the CA modeling: 1. Breakdown of existing database on hydraulic properties by HSUs 2. Constitutive model parameters for laboratory, core-scale hydraulic properties 3. Upscaling for macroscopic, field-scale flow, and transport parameters This document provides data and information on the following items: A brief description and technical basis for the conceptual model selected to represent vadose zone flow and transport for the CA modeling. Variation in sediment particle size distribution for representative sites in the 200 East and West Areas. A summary of effective (upscaled) moisture retention, and saturated and unsaturated hydraulic conductivity estimates for the HSUs in the 200 East and West Areas. A summary of effective transport parameters, including sediment bulk density and macrodispersivity estimates for various HSUs. The 200 East Area B Complex perched water aquifer hydraulic properties

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗

Vadose Zone Soil Flushing for Chromium Remediation: A Laboratory Investigation to Support Field‐scale Application

Cr(VI) flushing from the vadose zone to the groundwater (with subsequent Cr(VI) removal in groundwater by pump-and-treat system) is a promising remedial technique that has recently been used at field scale. This laboratory study was conducted to provide the technical basis to design a field soil flushing strategy. The objectives were to (1) quantify the relationship between sediment Cr(VI) and Cr(III) mass and release rates and subsequent Cr(VI) leaching; (2) investigate different methodologies to maximize Cr(VI) leaching, and (3) investigate methods to minimize leaching of remaining residual Cr. Characterization of Cr-contaminated sediments (Hanford Site, WA) exhibited Cr(VI) showed that leach rates that were correlated to different Cr surface phases. Sediments with low leachable Cr(VI) (<2 μg/g) leached Cr rapidly, so slow infiltration of water in a single pulse was sufficient to leach most Cr. In contrast, sediments with high Cr (2 to 200 μg/g) released some Cr(VI) quickly but 10 to 50% Cr(VI) slowly (tens to hundreds of hours). Efficient unsaturated leaching of these sediments required a different infiltration strategy that includes: multiple slow leach pulses with time between flushing cycles; the use of a surfactant to increase Cr leaching from low-permeability zones, and the use of a reductant (Na-dithionite or Ca-polysulfide) in the final leach water was highly effective at decreasing residual Cr leaching. This study clearly demonstrated that the methodology of basing laboratory Cr flushing on parameters such as Cr release mass and rates could be used to improve the efficiency of soil flushing at field scale.

Accelerated Cr Flushing↗

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗

Nitrate, ammonium, and DON mass time series output for East River stream, vadose zone and groundwater subwatersheds from HAN-SoMo model

There is a growing understanding of the role that bedrock weathering can play as a source of nitrogen (N) to soils, groundwater and river systems. The significance is particularly apparent in mountainous environments where weathering fluxes can be large. However, our understanding of the relative contributions of rock-derived, or geogenic, N to the total N supply of mountainous watersheds remains poorly understood. This dataset contains model output from the High-Altitude Nitrogen Suite of Models (HAN-SoMo), a watershed-scale ensemble of process-based models to quantify the relative sources, transformations, and sinks of geogenic and atmospheric N through a mountain watershed. Our study is based in the East River Watershed (ERW) in the Upper Colorado River Basin. The East River is a near-pristine headwater watershed underlain primarily by an N-rich Mancos Shale bedrock, enabling the timing and magnitude of geogenic and atmospheric contributions to watershed scale dissolved N-exports to be quantified. Six calibration scenarios were developed to explore equifinality using >1600 N concentration measurements from streams, groundwater, and vadose zone samples collected over the course of four years across the watershed, from 2014-2018. The files in this dataset are named according to which calibration scenario they represent (1, 2, 3, NM1, NM2, and NC), and include daily time series model output of nitrate, ammonium and dissolved organic nitrogen concentrations for all six calibration scenarios described in Maavara et al., 2021, "Modeling geogenic and atmospheric nitrogen through the East River Watershed, Colorado Rocky Mountains," PLOS One (in press), for the stream, vadose zone, and groundwater in the five discretized subwatersheds: Rustlers Gulch, Copper Creek, East above Quigley (EAQ), the Middle East River (ME), and Lower Triangle (LT) from Gothic to the Pumphouse. The READ_ME.txt file includes details on units and naming conventions in each of these files. Model files are in .mat format, to be opened in MATLAB.

54 ENVIRONMENTAL SCIENCES↗

Vadose Zone Model for B Complex 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B Complex vadose zone model. The B Complex model encompasses the region around single-shell tank farms 241-B, 241-BX, and 241-BY, and adjacent cribs, trenches, and burial grounds in the northwestern part of the 200 East Area. Sources of waste to this area were primarily from B Plant (221-B) and U Plant (221-U).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hanford Site Composite Analysis: Recharge and Inventory Sensitivity Analyses for the A Trenches, BC Cribs and Trenches, and PUREX Area Vadose Zone Models

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 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↗

Vadose Zone Model for A Trenches Area 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. 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 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 scope of this ECF is to document the development and results of the A Trenches Area vadose zone model. The A Trenches Area model covers the region between the 200 East Area and southwest of the 216-B-3 Pond. The primary waste sites in this area are cribs associated with PUREX Plant (202-A).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Model for B Plant Area 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B Plant Area vadose zone model. The B Plant Area model encompasses the region around B Plant (221-B) in the western part of 200 East Area. The plant was used to separate plutonium from irradiated fuel using the bismuth-phosphate process from 1945 to 1956, and for strontium and cesium recovery from 1968 to 1985.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Model for B-3 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3 Pond vadose zone model. The B-3 Pond model covers the main lobe of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Model for B-3C 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3C Pond vadose zone model. The B-3C Pond model covers the 216-B-3C Expansion Pond of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bedrock Vadose Zone Storage Dynamics Under Extreme Drought: Consequences for Plant Water Availability, Recharge, and Runoff

Bedrock vadose zone water storage (i.e., rock moisture) dynamics are rarely observed but potentially key to understanding drought responses. Exploiting a borehole network at a Mediterranean blue oak savanna site—Rancho Venada—we document how water storage capacity in deeply weathered bedrock profiles regulates woody plant water availability and groundwater recharge. The site is in the Northern California Coast Range within steeply dipping turbidites. In a wet year (water year 2019; 647 mm of precipitation), rock moisture was quickly replenished to a characteristic storage capacity, recharging groundwater that emerged at springs to generate streamflow. In the subsequent rainless summer growing season, rock moisture was depleted by about 93 mm. In two drought years that followed (212 and 121 mm of precipitation) the total amount of rock moisture gained each winter was about 54 and 20 mm, respectively, and declines were documented exceeding these amounts, resulting in progressively lower rock moisture content. Oaks, which are rooted into bedrock, demonstrated signs of water stress in drought, including reduced transpiration rates and extremely low water potentials. In the 2020–2021 drought, precipitation did not exceed storage capacity, resulting in variable belowground water storage, increased plant water stress, and no recharge or runoff. Rock moisture deficits (rather than soil moisture deficits) explain these responses.

54 ENVIRONMENTAL SCIENCES↗

Deep Roots Supply Reactivity and Enhance Silicate Weathering in the Bedrock Vadose Zone

In upland environments, roots commonly extend deep below soil into partially saturated bedrock. This Bedrock Vadose Zone (BVZ) has been shown to store and circulate water, host organic carbon respiration and serve as a critical source of rock-derived nutrients. However, the extent to which deep roots influence chemical weathering rates remains poorly understood. Here, we report 4 years of depth-resolved major ion chemistry over a 16-m thick BVZ hosting a deep rhizosphere in a catchment subject to a Mediterranean climate. These data allow development and validation of a reactive transport model (RTM), revealing that the timescales of water storage and drainage in the BVZ are sufficient to facilitate substantial chemical weathering of the shale bedrock. However, observed solute concentrations are only reproduced by the RTM when we explicitly include measured rates of CO 2(g) production meters below soil driven by the deeply rooted forest. By combining direct observations and a process-based RTM we conclude that the carbon respiration promoted by deep roots significantly enhances chemical weathering rates in the BVZ, constituting 43% ± 3% of total solute flux from the base of the BVZ to the water table.

Osorio‐Leon, Ivan D. [University of Illinois Urban↗

Selection of Vadose Zone Flow and Transport Properties with Gravel Fraction Corrections for the Hanford Site Composite Analysis and Cumulative Impact Evaluation

This environmental calculation file (ECF) is a compilation of hydrologic and radiological/chemical properties to be used for the updated Hanford Site Composite Analysis (CA) and Cumulative Impact Evaluation (CIE) vadose zone (VZ) modeling. The CA will provide an all-pathways dose projection to a hypothetical future member of the public from all planned low-level radioactive waste disposal facilities and potential contributions from all other projected end-state sources of radioactive material left at the Hanford Site following site closure. Its primary purpose is to support the decision-making process of the U.S. Department of Energy (DOE) under DOE O 435.1 Chg 1, Radioactive Waste Management, related to managing low-level waste disposal facilities at the Hanford Site. The CIE evaluates the effects of cleanup decisions regarding groundwater quality in the Hanford Site Central Plateau (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document). Due to the complexity and large number of waste sites in source operable units (OUs), the computational tools used for the CIE must be capable of representing a range of site conditions and source terms in the VZ while also efficiently computing the impact that cleanup decisions have on the underlying aquifer. Waste-site proximity between and within source OUs has resulted in contaminants commingling in the vadose and saturated zones in complex ways. Plume commingling requires cleanup decisions to be evaluated considering the surrounding waste sites and existing groundwater contamination, therefore demonstrating the need to evaluate cumulative impacts from the VZ to groundwater.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Composite Analysis Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Sources.

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 (eSTOMP1) 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↗

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗