Results for “Drainage” — Page 7

DOE OSTI2020

Field-Scale Lysimeter Studies of Glass and Cementitious Waste Forms at the Hanford Site - 20392

The Hanford site Integrated Disposal Facility (IDF) will receive waste forms from vitrification activities at the Hanford Waste Treatment and Immobilization Plant (WTP). The waste form inventory to be disposed of at the IDF will consist of vitrified low-activity waste (LAW) in glass forms and solidified (or encapsulated) secondary wastes in cementitious forms. The IDF is a near-surface burial facility located near Hanford's Waste Treatment and Immobilization Plant in the 200 East Area of the site's Central Plateau. An extensive set of laboratory experimental data has been collected to support IDF performance assessment calculations of the eventual degradation of waste forms and mobility of contaminants. This paper describes the start of a field experiment to generate waste-form performance data on a larger scale (tens of centimeters), over a longer duration (five years or more), and under field conditions representative of the IDF. Results from this study are expected to improve model descriptions of contaminant mobility, reduce uncertainties about the representativeness of laboratory results in the IDF performance assessment, improve stakeholder confidence in the safe disposal of treated waste at the IDF, and facilitate the adoption of informed facility designs with the potential to reduce operational costs. A lysimeter test facility near Hanford's 200 West Area has been repurposed to carry out long-term assessments of waste forms buried in sediments excavated from the IDF. The waste forms are buried in large drainage lysimeters, open caissons that are two meters in diameter and three meters deep. Waste forms are buried in up to six locations in each lysimeter. Samplers collect pore water and drainage water beneath the waste forms and air from the pore space adjacent to waste forms. Sensors measure temperature, soil water content, soil water tension, and drainage flux. One lysimeter containing cementitious waste forms and one containing glass waste forms have been completed, along with a control lysimeter (without waste forms). The cementitious waste forms are Cast Stone grout formulations using a liquid secondary waste simulant and a generic Hanford high-salt simulant, as well as a hydrated lime-based grout formulation using a liquid secondary waste simulant. Technetium-99 and iodine-127 were added as tracers to monitor contaminant mobility. Glass waste forms were fabricated in two formulations (LAWA44 and ORLEC28), with rhenium and molybdenum tracers to monitor waste form degradation. The predicted performance of the waste forms in the lysimeters was modeled prior to installation to determine waste form size, tracer concentration/activity, and locations of sensors and samplers. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Revisiting the definition of field capacity as a functional parameter in a layered agronomic soil profile beneath irrigated maize

The soil water content at the condition of field capacity (θ FC ) is a key parameter in irrigation scheduling and has been suggested to be determined by running a synthetic drainage experiment until the flux rate (q) at the bottom of the soil profile achieves a predefined negligible value (q FC ). We question the impact of q FC on the assessment of field capacity. Moreover, calculating θ FC as the integral mean of the water content profile when q is equal to q FC is strictly valid only for uniform soil profiles. By contrast, this practice is ambiguous and biased for stratified soil profiles due to the soil water content discontinuity at the layer interfaces. In this study, the concept of field capacity was revisited and adapted to practical agronomic heuristics. By resorting to the assessment of root-zone water storage capacity (W), we envision field capacity as a functional hydraulic parameter derived from synthetic irrigation scheduling scenarios to minimize drought stress, drainage, and nitrate leachate below the root zone. A functional analysis was carried out on a 135-cm-thick layered soil profile beneath maize in eastern Nebraska. On-farm irrigation scheduling applications and agricultural practices were recorded for 20 years (2001–2020) at a daily time step. Hydrus-1D was calibrated and validated with direct measurements of the soil water retention curve and soil water content data, respectively, in each soil layer. A set of functional field capacity values was derived from 24 irrigation scheduling scenarios, and the optimal water storage capacity at field capacity (W FC ) was approximately 50 cm (corresponding to about 80% saturation in the soil profile). An average irrigation amount of 217.5 mm distributed over 21 events was obtained by using optimal irrigation scheduling, which was initiated when the matric pressure head took on a value of –700 cm and the irrigation rate was set at 1.0 cm d –1 . This irrigation practice ensured water storage at approximately the same level (ideally at W FC ) by sustaining only evapotranspiration fluxes in the uppermost portion of the root zone and by limiting excessive drainage. This protocol can be transferred to other agricultural fields.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗

DOE OSTI2026

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

DOE OSTI2020

Nevada National Security Site 2019 Waste Management Monitoring Report Area 3 and Area 5 Radioactive Waste Management Sites

Environmental monitoring data are collected at and around the Area 3 and Area 5 Radioactive Waste Management Sites (RWMSs) within the Nevada National Security Site (NNSS). This report summarizes the 2019 environmental data to provide an overall evaluation of RWMS performance and to support environmental compliance and performance assessment (PA) activities. Some of these data (e.g., radiation exposure, air, and groundwater) are presented in other reports (Mission Support and Test Services, LLC 2019, 2020a, 2020b). Direct radiation monitoring data indicate exposure levels at the Area 3 and Area 5 RWMSs are within the range of background levels measured at the NNSS. Slightly elevated exposure levels outside the Area 3 RWMS are attributed to nearby historical aboveground nuclear weapons tests. Air monitoring data at the Area 3 and Area 5 RWMSs show that tritium concentrations in water vapor and americium and plutonium concentrations in air particles are below Derived Concentration Standards for these radionuclides. Groundwater monitoring data indicate the groundwater in the uppermost aquifer beneath the Area 5 RWMS is not impacted by RWMS operations. Results of groundwater analysis from wells around the Area 5 RWMS are below established investigation levels. Leachate samples collected from the leachate collection systems at the Area 5 mixed low-level waste disposal unit are below established contaminant regulatory limits. During 2019, precipitation at the Area 3 RWMS was 77 percent above average, and precipitation at the Area 5 RWMS was 69 percent above average. Water balance measurements indicate that evapotranspiration from the vegetated weighing lysimeter at the Area 5 RWMS dries the soil and prevents downward percolation of precipitation more effectively than evaporation as measured from the bare-soil weighing lysimeter. Vadose zone monitoring in the Area 3 and Area 5 RWMS soil covers shows no evidence of precipitation percolating through the covers to the waste. Moisture from precipitation did not percolate below 120 centimeters (3.9 feet [ft]) in the vegetated final cover on the U-3ax/bl disposal unit at the Area 3 RWMS during 2019. There was no drainage through 2.4 meters (8 ft) of soil as indicated from the Area 3 drainage lysimeters that received only natural precipitation. At the Area 3 RWMS, which received three times the natural precipitation, 57 percent of the applied precipitation and irrigation drained from the bare-soil drainage lysimeter. All 2019 monitoring data indicate that the Area 3 and Area 5 RWMSs are performing within expectations of the model and parameter assumptions for the facilities’ PAs.

2019↗

DOE OSTI2022

Clean Condensing Gas Furnace

Natural gas furnaces are the most common space heating equipment in the U.S. residential and commercial building markets. However, current residential natural gas condensing furnaces generate substantial acidic condensate as well as significant emissions of sulfur oxides (SOx), nitrogen oxides (NOx), carbon monoxide (CO), hydrocarbons (HC), and methane (CH 4 ) contributing to environmental degradation of air, water, and soil. This report describes a novel solution to reduce the environmental impact of natural gas condensing furnaces based on the technology of a monolithic acidic gas reduction (AGR) catalyst for SOx trapping, NOx redox to nitrogen, and oxidation of formic acid, CO, HC, and CH 4 . The AGR technology offers the following benefits: (1) a neutral furnace condensate with a pH of ~7, allowing its safe release into the sewer system thus eliminating a second drainage system; (2) trapping and removing nearly all SOx emissions; (3) NOx emissions nearly at nearly 1-2 ng/J, more than 95% lower than new emissions standards in California; (4) the use of a low-cost heat exchanger as a condensing heat exchanger (HX) since the condensate is not acidic, avoiding the need for expensive stainless steel alloys; and (5) unburnt fuel energy recovery to boost efficiency.The AGR component and AGR-enabled furnace performance were broadly tested to determine their effects on long-term reliability and durability, as well as SOx storage and regeneration activity. The AGR regeneration does not impair the performance in achieving neutral condensate and ultra-low NOx emissions, and the AGR catalyst subjected to regeneration activities continued to function well and achieved slightly better annual fuel utilization efficiency (AFUE). The 400-hour reliability and durability test of the retrofitted condensing furnace with the AGR component shows that the furnace unit still achieves a neutral furnace condensate with a pH of ~7 and enables 0~3 ng/J of NOx emissions. However, the 400-hour operation slightly degraded the AFUE because of soot particle accumulation caused by frequent incomplete combustion owing to inappropriate condensate drainage during testing. Thus, proper condensate drainage is critical for AGR-enabled furnaces. Furthermore, neutron computed tomography was employed to survey the aged AGR component and demonstrate high-resolution 2D and 3D representations for the nondestructive diagnosis of the AGR component. The tomography showed that the AGR component did not deform or suffer broken AGR channels. A new AGR catalyst with low precious metal loading was preliminarily explored to identify a pathway of optimizing AGR material loading and maximizing acidic gas reduction at low cost. The new AGR component can reduce precious metal loading by 38% and still achieve neutral condensate and ultralow NOx emissions. The furnace with the AGR component of low Pt/Rh loading enables a maximum AFUE of 97%, which is meaningfully higher than the original furnace. Long-duration testing for the furnace enabled with the low precious metal loading AGR component will be vital in future research. Although the current work demonstrates a proof of concept for the AGR-enabled furnace, the AGR assembly needs to be optimized and integrated into the design of new OEM furnace products. Furthermore, the AGR technology can be applied not only for residential gas furnaces, but also for commercial rooftop units, gas heat pumps, gas-fired water heaters, combustion boilers, and other systems.

03 NATURAL GAS↗

DOE OSTI2024

Computational Analysis of Hydraulic Capacity of Ohio DOT Catch Basin No. 6 in On-Grade and Sag Locations

Stormwater runoff from streets and highways is typically captured by drainage structures strategically placed when roadways use curb or curb and gutter systems. These structures include catch basins with grates, inlets, or combination grates/inlets that collect and discharge storm water runoff to buried pipe conveyance systems. The performance of these drainage structures is measured in terms of hydraulic efficiency, which is defined as the percentage of flow captured by the basin as compared to the total flow draining to the structure. Understanding the performance of these drainage structures allows for a safe and economical design which prevents flooding along Ohio roadways.

42 ENGINEERING↗

DOE OSTI2021

Computational Study of Hydraulic Performance of South Carolina DOT Catch Basins Type CB 1, CB 16, and CB 17

Proper design of surface drainage of roadways is essential to minimize flooding and provide for traffic safety. Inlets collect the excess storm water from the drainage area of a roadway and discharge it to storm drains. Knowing the hydraulic efficiency of inlets, defined as the percentage of intercepted flow to the total street flow, is necessary in drainage design to determine inlet spacing such that the system can transport all or the majority of the road surface flow during rain events off of the road into the catch basins.

42 ENGINEERING↗

DOE OSTI2025

Prototype Hanford Barrier Performance Monitoring Report: Fiscal Year 2024

This report provides an annual update on performance monitoring of the Prototype Hanford Barrier (PHB) at the Hanford Site. The PHB is part of a long-term study that serves as the scientific basis for many of the evapotranspiration-capillary barrier designs currently used globally and for future engineered barrier designs planned for remedial actions over waste and demolition sites at the Hanford Site. The PHB allows us to identify potential future issues and develop better monitoring techniques before these engineered barriers are constructed. Surface barriers like the PHB are essential for preventing water infiltration and curbing the spread of contaminants to groundwater. Effective monitoring of soil moisture levels above and below these barriers is crucial given that performance metrics could span up to 1,000 years. From July 2023 to June 2024, the water flux (as measured by tipping buckets) through the 2-m-thick silt loam layer of the PHB – a fine silty material that stores water under high capillary tension – remained well below the 0.5-mm-per-year performance threshold, demonstrating its effectiveness in preventing water penetration. In 2024, degraded tipping bucket gauges were replaced to ensure the accuracy of future data. Additionally, neutron probes revealed that the wetting front from the rainy season only penetrated to a maximum depth of 1.2 m into the silt loam. This further demonstrated the barrier's efficiency, as the wetting front did not fully penetrate the 2-m-thick capillary barrier. The western gravel slope of the PHB, composed of Hanford Site pit gravel of varying sizes, demonstrated low flux rates. In contrast, the eastern riprap slope exhibited higher flux rates. Zhang (2017) found that the riprap side slope of the PHB had the highest drainage rates in January and the lowest in late summer or early fall, indicating significant summer evaporation. Drainage rates increased significantly under enhanced precipitation conditions, far exceeding the design criterion, which could lead to water infiltration into the waste zone. The study introduced the “edge effect,” where elevated drainage rates from the riprap may migrate laterally beneath the barrier, compromising its ability to isolate waste, and recommended expanding the barrier and conducting further research to mitigate this issue. This report provides a review on surface barrier edge effects and their potential impact to subsurface contaminant migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.

Earth Systems↗

DOE OSTI2023

MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.

Earth Systems↗

DOE OSTI2023

MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.

Earth Systems↗

DOE OSTI2023

MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.

Earth Systems↗

DOE OSTI2021

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↗

DOE OSTI2023

Genetic monitoring of steelhead in the Klickitat River to estimate productivity, straying, and migration timing

Abstract Objective Salmonids with a complex life history variation present challenges for conservation management, but genetic approaches alongside fisheries monitoring can address questions regarding the viability of the natural populations. Methods We genotyped adult (n = 3108) and juvenile (n = 2624) samples of steelhead Oncorhynchus mykiss that were collected in the Klickitat River, Washington, USA at traps in the lower drainage to examine tributary level productivity, straying from outside sources, and variation in adult migration timing. Result Genetic assignment of steelhead from this system indicated that the majority were produced within or near tributaries of the middle Klickitat River (juvenile mean = 72.8%; adult mean = 87.3%). Analyses with parentage-based tagging identified that most hatchery-origin adults were assigned to the Skamania Hatchery (80.8%), as expected, since this has been the release stock for decades within the Klickitat River drainage. Hatchery-origin adults were also identified from programs operating outside the Klickitat River, which were primarily strays from Snake River hatcheries. Most natural-origin steelhead were assigned to the Klickitat River, but there were also natural-origin fish identified as strays from other regions of the Columbia River (22.3% of natural returns). We also examined genes known to be associated with migration timing in adult steelhead observed at the trap and observed a strong relationship between migration date and alleles for early and late migration, but individual outliers were detected across seasons. Conclusion Our results indicate that genetic variation of steelhead in the Klickitat River has been influenced by hatchery programs as well as natural-origin straying from other subbasins, but genetic diversity remains high throughout the subbasin, and both early and late migration alleles are maintained. The genetic diversity present in Klickitat River steelhead may enable this Endangered Species Act listed (threatened) species to better adapt to stochastic environmental conditions compared to less diverse populations.

Collins, Erin E. (ORCID:0000000260978479)↗

DOE OSTI2023

Watersheds and stream networks viewed longitudinally: Example insights from novel spatial portrayals of watershed characteristics

Abstract Longitudinal depictions of watershed structure and characteristics, including topography, stream networks, wetlands, ground water levels, and land use, can provide watershed knowledge and understanding unavailable from standard plan view maps. Three case studies provide examples of knowledge gained by applying longitudinal views of stream networks, watershed hydrologic behavior, and land use distributions. Longitudinal views of mountain stream networks show extreme variability in the slope‐area relationships of low Strahler order streams, large discontinuities in drainage area (large parts of drainage area space are absent in networks), and large variations in network curvature. Longitudinal views of a groundwater‐dominated headwater watershed increase the inference available from limited groundwater observations and clearly reveal how groundwater connections affect the permanence of surface water features and the distribution of vadose zone storage in the landscape. Plotting land uses longitudinally illuminates and allows a quantitative analysis of how land uses are distributed relative to topographic position. Viewing watersheds and stream networks longitudinally can provide new insights into watershed forms and processes and motivate new questions and research.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Interdisciplinary Approaches Improve Understanding of Cryptogenic Species: A Historical Case Study of Crayfish in Montana, USA

ABSTRACT Cryptogenic species are those that are not yet “demonstrably native or introduced” in a given area, such as crayfish in western Montana, USA. Delving into evidence from diverse fields can help clarify the status of cryptogenic species. Primary historical sources and indigenous knowledge have informed various ecological questions but are seldom used to clarify the status of cryptogenic species, especially aquatic taxa. We summarize types of evidence used to illuminate species status and offer a case study applying historical sources to cryptogenic crayfish in Montana. We searched for crayfish mentions in dictionaries of indigenous languages and primary historical documents (e.g., travel journals, newspaper accounts, early biological surveys). Early explorer accounts we examined did not mention crayfish in Montana, though they noted them in the nearby Snake River drainage of Idaho and Wyoming. The first mention of crayfish in Montana that we found was from an 1868 newspaper account, probably referencing the upper Missouri River. The species was not noted, but an 1887 article suggested signal crayfish Pacifastacus leniusculus introductions near Bozeman, where non‐native signal crayfish still persist. The earliest crayfish mention we found west of the Great Divide in Montana was from 1934 in Ninepipes and probably referred to non‐native virile crayfish ( Faxonius virilis ). Signal crayfish were not mentioned west of the Divide until a 1944 newspaper account noted their 1937 introduction to the Bitterroot River drainage. The case study demonstrates the value of historical sources in identifying the absence of species and documenting their presence or introduction during periods predating scientific data.

Kennedy, Hampton L. [USDA Forest Service, Southern↗