Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electrical conductivity of soil”

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 73 records · Page 4

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↗

3D hydrogeophysical characterization of managed aquifer recharge basins

Aquifers are increasingly stressed. Managed aquifer recharge provides a potential solution to mitigate this stress and provide sustainable groundwater resources. Subsurface properties are known to have a strong control on the infiltration rates that can be achieved. However, these properties are often highly heterogeneous and difficult to assess with conventional probing techniques. Here, we show the application of 3D geophysical imaging to assess the recharge potential and its variation across several basins used for managed aquifer recharge. We link in-situ measurements of saturated hydraulic conductivity with the electrical resistivity of the subsurface to establish petrophysical relationships and use those relationships to estimate the distribution of hydraulic conductivity throughout the five recharge basins. Our results show a considerable variability in the hydraulic properties, i.e., soil texture and saturated hydraulic conductivity, that have a direct impact on potential infiltration rates. We use the 3D hydraulic property distributions to model groundwater recharge and provide estimates for infiltration rates and volumes, and use this approach to assess the impact of management activities on groundwater recharge performance. Having such data not only enables us to predict infiltration rates, but also provides means for optimizing such water infrastructure.

54 ENVIRONMENTAL SCIENCES↗

Cytochrome Nanowire XRD

Microbial metabolism plays an important role in the global cycling of carbon, nutrients, and metals. Recent discovery of microorganisms sharing energy using direct electrical connections might represent a prominent strategy by which anaerobic microorganisms interact with the diversity of environments. The purpose of the proposed research is to apply these newly developed principles of electron exchange in microbial communities to the study of biogeochemical cycling in climate and subsurface systems. Previous studies revealed that soil bacteria Geobacter sulfurreducens produce conductive pili nanofilaments that facilitate long-range electron transport to extracellular electron acceptors and other cells. This discovery has already transformed our understanding of the function of microbial communities in diverse environments, making basic contributions to microbial ecology and helping to improve practical applications such as increasing current output of microbial fuel cells and enhancing the conversion of organic waste to methane. Microbial electron exchange via pili may also contribute to methane production in terrestrial environments that are major sources of atmospheric methane. However, the mechanism of electron transport through pili is not well understood. Initial studies have suggested that pili exhibit metallic-like conductivity similar to conducting polymer polyaniline. However, these studies were limited to electrical measurements. X-ray diffraction studies have suggested that aromatic amino acids tyrosine and phenylalanine are closely packed in pili. But these studies used homology modeling to predict the pilus assembly and information about the number of monomers in each pilus assembly unit and their organization are still lacking. The central hypothesis of this proposal is aromatics in pili enable intermolecular electron delocalization due to pi-stacking that can give rise to metallic-like conductivity. By analyzing pili structure using infrared scattering-scanning near field optical microscope (IR s-SNOM) and combining with solution NMR, Focused Ion Beam (FIB and computational capabilities, this work will provide comprehensive understanding of the conduction mechanism from a structure-function perspective. The specific aims are: 1) With infrared scattering scanning near field optical microscope (IR s-SNOM), image the organization and density of pilin monomers and pH-induced conformational changes in a pilus filament. 2) With solution NMR, resolve pH-induced structural changes in the pilin monomer. 3) Using Focused Ion Beam (FIB), draw electrodes on single pili to measure their conductivity. 4) With computational modeling using the NWChem, simulate the effect of environment on the electronic properties of pili as a function of pH. The proposed work will gain the nanoscale insight into the structure of the pili, as the basis for understanding interspecies electrical communications. The improved understanding will be harnessed to accelerate bioenergy production and bioremediation using pili as well as for predictive modeling of carbon cycling in diverse environments to reach the central goals of DOE-BER to develop sustainable energy sources, regulate the contaminants in the subsurface and understand the effects of greenhouse gas emissions on biosphere. The proposed studies will provide design principles for engineering of electrical interactions among these communities by manipulating the amino acid composition of pili. Moreover, these studies will identify the forces that maintain the function of these communities in the face of environmental perturbation. Therefore, these studies will help to address the need "to incorporate process understanding of biogeochemical cycling into a scalable hierarchy of predictive capabilities" specified in EMSL's Science Themes.

59 BASIC BIOLOGICAL SCIENCES↗

Life Cycle Emissions and Health Cost Impacts of Producing Ethanol and Electricity from Willow and Switchgrass in the Riparian Buffers of Mid-Atlantic United States

The United States (US Mid-Atlantic Region (MAR) has the potential to grow a variety of perennial feedstocks such as switchgrass and shrub willow to increase domestic energy production. These cellulosic feedstocks have also shown improved ecosystem services, such as soil carbon sequestration, nitrate leaching reduction, and flood mitigation along rivers and streams as partially harvested riparian buffers. To examine the effects on greenhouse gases (GHGs) and criteria air pollutants (CAPs) from using these feedstocks to produce ethanol or electricity, we conducted a comprehensive life cycle assessment (LCA) and estimated the impact on human health costs when land use is changed from corn production for ethanol. Results indicate up to 54% reduced GHG per hectare from using willow and switchgrass feedstock sources to produce ethanol instead of corn. However, there was a trade-off in terms of CAP emission, as grass-based energy emitted more NO x and SO x compared to the corn ethanol pathway, except for SO X emissions from willow-based electricity. Electricity from cellulosic biomass had higher particulate matter (PM) emission compared to that from corn ethanol. Estimates for health cost to society ranged from $2498 ha –1 for electricity from switchgrass to a net benefit of $448 ha –1 for ethanol production from willow, depending on varying biomass yield under different market scenarios. Although using cellulosic feedstocks to produce bioenergy has great potential to reduce GHG emissions, CAP control measures are needed to manage CAP-induced health costs.

Bioenergy↗

Spread Spectrum Time Domain Reflectivity for String Monitoring in PV Power Plants (Final Technical Report)

This final report describes the methods and results of applying Spread Spectrum Time Domain Reflectivity (SSTDR) for String Monitoring in PV Power Plants for DE-EE0008169. The project created a new system for both detecting and locating electrical faults in photovoltaic systems. In this work, we address photovoltaic electric faults that are both common and costly. Based on interviews with photovoltaic power plant owners, operators, and maintainers, three types of faults are common and of significant interest: disconnects, ground faults, and arc faults. Disconnects can originate from many sources. They are often due to everyday events, such as lawnmowing (accidentally running over a cable), animals eating through the cables, or degradation that occurs over time due to corrosion or general degradation. Ground faults occur when the cables (for example, due to frayed insolation) connect to the ground, relaying current into the ground. These faults are particularly problematic since the ground faults are often intermittent. That is, ground faults commonly appear during rain storms due to a change in soil conductivity and then disappear when the rain ends. This makes the ground fault difficult to find because while current systems can detect the overall change in voltage and current associated with a ground fault, technicians are necessary to locate the fault. As a result, ground faults may disappear before the technician arrives at the power plant. Hence, locating and fixing ground faults often require multiple trips. We also study arc faults, which can result when nearby conductors create an arc of electrical current through the air. While less common, arc faults can be extremely dangerous. The energetic electrical arc can cause fires and destroy equipment, costing significant damage. Overall all three types of faults cost owners and operators money, either from the destruction of equipment or from technician time. Furthermore, while devices exist for detecting ground faults (ground fault circuit interrupters) and arc faults (arc fault circuit interrupter), these systems only search patterns of electrical current that correspond to each fault. This information cannot be used to locate the fault. In addition, these protection systems experience nuisance trips due to nearby electromagnetic interference, such as from a lawn mower or other motors that produce significant amounts of electromagnetic radiation. Hence, the overall goal of this project is to create an SSTDR tool that provides photovoltaic power plants with more reliable fault detection in addition to the localization of faults. SSTDR works by transmitting electrical signals into the photovoltaic string. Those signals reflect from impedance discontinuities (i.e., disconnects, ground faults, and arc faults). These faults are then detected by measuring the presence of a reflection at the SSTDR and can be located by identifying the location of that reflection in time. In addition, unlike current protection systems, these systems do not experience nuisance trips since their low amplitude, high frequency, and coded signal can by analyzed without interference from the regular operational voltage on the photovoltaic string.

14 SOLAR ENERGY↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗

An integrated Imaging and Modeling Toolbox for Accelerated Development of Root-focused Crops at Field Scales

This project, led by Yuxin Wu at the LBNL, was proposed to develop an integrated imaging-modeling toolbox allowing for accelerated development of root-focused crops at field scales. Our approach is based on a novel root phenotyping method we termed Tomographic Electrical Rhizosphere Imaging (TERI), which provides advanced phenotyping of key root traits by sensing the electrical impedance response of roots and soil to external electrical excitations, and mapping these responses to root and soil properties. The role of Noble Research Institute in this project is to assist TERI model establishment and validation by collecting various above- and below-ground traits data from container- and field-grown plants (mainly wheat). Over the last three years, Noble Research Institute has conducted 3 phases of experiments as proposed including: (Phase 1) greenhouse and hoop-house experiments in wheat and pecan plants grown under controlled container-grown conditions; (Phase 2) space-planted individual plant field experiment using different wheat varieties; and (Phase 3) small-plot planted field experiment using different wheat varieties, as well as using different small grains species including wheat, rye, triticale, oat and barley for increasing plant variation.

60 APPLIED LIFE SCIENCES↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗

NETL Plastic Pipes Project - June 2022 Progress Report

Plastic or composite pipelines have been the bane of the utility locating industry because conduits are neither conductive nor magnetic, which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipeline where the resistive cover allows for adequate penetration of radar energy. However, GPR cannot be used in areas where the soil cover is conductive due to significant clay and/or salt content. This study takes a comprehensive look at near-surface geophysical methods that potentially are useful for locating buried plastic/composite pipelines, either singly or in combination with other geophysical methods. Specifically, this modeling study uses computational numerical methods to forward model the response of GPR, resistivity, seismic, and gravity gradiometry methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft to 4 ft; (3) various degrees in contrast in physical properties (i.e., permittivity, elasticity, electrical resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar frequencies, seismic source frequency, electrode spacing).

47 OTHER INSTRUMENTATION↗

Defining a compact dry cooler design to reduce LCOE contribution in a CSP facility

Concentrating solar power (CSP), when coupled with a supercritical carbon dioxide (sCO 2 ) power cycle and sensible heat storage, presents a renewable and clean alternative for utility-scale power generation. However, in order to be competitive in the current and future markets, CSP facilities must limit their levelized cost of electricity (LCOE) by minimizing capital costs and reducing operating costs over the lifetime of the plant. Targeting this goal, this study investigates the LCOE impact of the power cycle pre-cooler. This study considers a compact dry cooler with micro-channel technology on the CO 2 side and formed fin geometry on the air side, using directly-coupled centrifugal fans and a transition duct to improve air distribution across the fins as well as protect the fins from contaminants which may cause blockage, soiling, fouling, and damage. In an effort to better understand the dry cooler impact on LCOE, a sensitivity study was conducted using various combinations of end-to-end approach temperatures, air-side pressure drop values, CO 2 -side pressure drop values, fan types, cooler turndown control schemes, cooler module sizes, and design-point ambient temperatures. Furthemore, off-design cycle performance data was calculated for each dry cooler design using NPSS simulation software; cycle performance data were then input to System Advisor Model (SAM) along with the associated capital costs for LCOE prediction of a 100 MW system over a 30 year plant lifetime. Results of this study show the LCOE is most sensitive to air-side performance, followed by heat transfer effectiveness and capital cost. It was found that a power cycle with a mid- to high-performance dry cooler will produce the most competitive power-production costs. Designing at the extreme ends for approach temperature (or effectiveness), design-point ambient temperature, and compactness (footprint) produce higher LCOE values; mid-range values for these parameters balance performance, operating costs, and associated capital cost to optimize LCOE.

14 SOLAR ENERGY↗

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗

Microbial potentiometric sensor array measurements in unsaturated soils

The overarching goal of this study is to demonstrate a novel technology for monitoring changes in electrical potential of unsaturated soils using biofilm-populated electrodes. The novelty of the study stems from the fact that it demonstrates a method for measuring open-circuit potentials (OCP) in environments without the presence of an electrolyte solution. This study also reveals that using a biofilm-populated electrode as a reference in stable environments could successfully be employed to assess and monitor the electrochemical potential generated by plants and microorganisms. The findings imply that long-term (months to years) and real-time measurements of the open-circuit potential in unsaturated soils are possible. Because MPS arrays can directly measure open-circuit potential from the biofilm, the challenges related to locally induced electrochemical changes caused by microorganisms in the soil to achieve optimum physiological levels are eliminated. The simplicity of the technology, which allows for multiple indicator electrodes to be referenced against an “internal” reference electrode, enables spatial-temporal monitoring of the changes in the soil and the generation of 2D- or 3D-signal patterns. Once a signal pattern, generated by an array of sensors, develops (usually after 30 to 90 days), it does not significantly change unless the soil is exposed to external stimuli. The observed OCP phenomena, however, suggests that the change in OCP signal is independent of changes in soil conductivity resulting from the addition of water. In brief, findings suggest that the proposed technology can enable multidimensional profiling and long-term monitoring of changes occurring in unsaturated soils without direct implications of presence of water. The changes in the 2D or 3-D signal patterns, however, can be correlated to other important parameters that characterize soil health.

54 ENVIRONMENTAL SCIENCES↗

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗

Subtask 1.3 – Integrated Carbon Capture and Storage for North Dakota Ethanol Production

The Energy & Environmental Research Center (EERC), in partnership with Red Trail Energy, LLC (RTE), a North Dakota ethanol producer; the North Dakota Industrial Commission (NDIC) Renewable Energy Program (REP); and the U.S. Department of Energy (DOE), conducted a feasibility and implementation study for a commercial carbon capture and storage (CCS) effort. This subtask provided technical support and developed recommended practices of how small-scale industrial CO2 emitters (<1,000,000 tonnes of CO2 emitted annually) may economically deploy CCS. The 64-million-gallon dry mill RTE ethanol facility, emitting an average 180,000 tonnes of CO2 annually, was the subject of the case study. The positive outcome of this research, which shows technical and economic potential for ethanol–CCS in North Dakota, has resulted in RTE acquiring an approved Permit to Drill on December 2, 2019, for a stratigraphic test well in early 2020, a necessary step toward a North Dakota CO2 Storage Facility Permit (SFP) for the RTE CCS effort. Outcomes include 1) validating the Broom Creek Formation as a regional target for CCS, 2) determining the full carbon life cycle of an industrial fuel production facility with CCS, 3) developing a field implementation plan (FIP) for small-scale CCS, and 4) determining the validity and pathway for using CCS to meet low-carbon fuel (LCF) standards. The technical team included the EERC, RTE, Trimeric Corporation, Schlumberger Carbon Services, and Computer Modelling Group (CMG). Findings from activities conducted November 2016 – May 2020 are summarized as follows. Several key steps have been accomplished toward validating the Broom Creek Formation as a CO2 injection and storage target at the RTE CCS site, including verification of the presence and structure of sandstone layers that may comprise the potential CO2 storage reservoir and the several thousand feet of overlying confining zone. Work included site characterization using existing data, interpretation of seismic data within the study area, and geologic modeling simulation of CO2 injection. Interpreted results estimate 3000 feet of confining zone between the Broom Creek Formation (storage target) and the lowermost underground source of drinking water (e.g., the Fox Hills Formation). The thickness of the Broom Creek injection target varies 230–420 ft within the survey area. Results were used to inform the location of a stratigraphic test well and associated characterization test program. No impediments were identified within the targeted CO2 storage complex that would prevent the project from moving forward. A stratigraphic test well is the next step to validate these results and to acquire remaining data necessary to develop a North Dakota CO2 SFP application. The EERC generated full carbon life cycle estimates for CCS integration with the RTE ethanol facility. Results indicated that an average 40% reduction in CO2 emissions is possible through CCS implementation. Approximately half of the carbon in the overall life cycle for a dry mill ethanol plant is generated through the fermentation process and emitted to the atmosphere; this is the CO2 stream targeted for CCS. The remaining carbon is attributed to corn feedstock farming (diesel, fertilizer), energy for fuel processing (natural gas, electricity), and transportation (diesel) of the fuel product. Life cycle carbon estimates are also affected by the anticipated energy consumption of a potential capture facility, which depend on the type of CO2 product generated. For example, a 30%–40% net CO2 emission reduction is estimated if a liquefied CO2 facility were incorporated compared to a 40%–50% net CO2 reduction if a supercritical “injection-grade” CO2 product is generated; i.e., more energy is required to further refine the CO2 stream, affecting the full carbon life cycle estimates. A CCS FIP was developed and initiated at the RTE CCS site, resulting in the development of several guidance documents: a CO2 Capture Process Design Package, a North Dakota CO2 geologic Storage Permits Template, and a Public Outreach Package for CCS in North Dakota. General FIP components include CO2 capture system, pipeline and well designs; monitoring, verification, and accounting (MVA) plans; geologic characterization and testing programs; and permitting and outreach plans. Vendor bids were also acquired for the CO2 liquefaction facility. Near-surface characterization (groundwater and soil gas sampling) and geologic characterization (seismic survey) were initiated to inform development of a UIC Class VI-compliant MVA plan compliant with a North Dakota CO2 SFP. Designs (well and geologic testing) were completed for a stratigraphic test well compliant with a North Dakota CO2 SFP. In addition, the outreach plan was executed, including community open houses, meetings with city/county/state officials, and development of public materials. Although other entities continued to mature incentive programs in 2019–2020, California and the Internal Revenue Service (IRS) currently provide the most advanced economic opportunities for CCS integrated with fuel production. The California Low-Carbon Fuel Standard (LCFS) adopted a CCS Protocol in January 2019, allowing submittal of a design-based pathway (DBP) application for an approved temporary (not certified) carbon intensity value for a fully engineered facility. An ethanol–CCS DBP application to the California LCFS Program (officially approved February 28, 2020) was developed to show that the RTE CCS effort meets LCFS requirements. The approved DBP provides confidence to advance the project and supports potential investments. Other entities continue to mature incentive programs. The IRS issued guidance in February 2020 that addresses the definition of beginning of construction and revenue procedure on partnerships for the Enhancement of Carbon Dioxide Sequestration Credit (a.k.a. Section 45Q) CCS tax credit program; the IRS anticipates issuing further guidance on issues such as secure geologic storage, utilization qualifications, and recapture of claimed credits. Maturing incentive programs coupled with workable permitting regulations provide confidence to advance CCS projects in North Dakota and support financial investment to proceed with designing, constructing, and implementing CCS projects at small-scale fuel production facilities. The largest hurdles for CCS implementation at small-scale industrial systems are often business/economic-related (i.e., not technical). Market uncertainty already exists for agriculture-based alternative fuels such as corn ethanol, for which production has increased by ~33%, and prices have correspondingly lowered since 2015. Passing the Section 45Q tax credit program improves economic feasibility for CCS but may require external investors for a small business to achieve maximum benefits. Public–private partnerships with NDIC and DOE have resulted in foundational technical and regulatory knowledge, growing stakeholder confidence, and a pathway to implementation that enables similar industrial CCS projects in the region to advance. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by NDIC and RTE. The authors would also like to thank CMG, ESRI, IHS, Neuralog, and Schlumberger for allowing the use of their software packages in support of this work.

42 ENGINEERING↗

Hydropower Geotechnical Foundations: Current Practice and Innovation Opportunities for Low-Head Applications

Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to challenges frequently encountered with traditional approaches. With this understanding, it is critically important to understand and document the current state of practice for hydropower geotechnical foundations, identify key challenges, and define opportunities for innovative solutions. To this end, this report documents the current state of practice across the three main phases of geotechnical foundation development: (1) geotechnical site assessment, (2) design, and (3) construction for hydropower systems. It also describes the major challenges with conventional approaches and identifies opportunities for innovation to reduce hydropower foundations costs, timelines, and risks. Key takeaways from this report include the following: Approximately 80% of available low-head sites are expected to have foundations on soil beds rather than rock beds, suggesting that rockfill and earthfill dams may be the most cost-effective conventional dam type for new projects.; Geotechnical and geologic investigation activities are time-consuming and expensive but are essential to define the parameters and criteria needed for foundation design.; Certain riverbed soil and bedrock types present significant technical challenges or require expensive foundation construction, which can prove financially prohibitive for low-head project development.; Modular hydropower design and prefabricated modular foundations represent a promising but unproven paradigm for new hydropower development. Design and construction approaches using optimized and highly repeatable, reliable components would benefit project cost, time, and risk but require additional research and development.; Temporary construction features for foundations, including cofferdams, water diversion, and water control systems, can prove costly and have inherent construction risk.; For economically viable development, hydropower geotechnical foundations should be limited to 4 to 15% of the project’s total initial capital costs. Many proposed projects have experienced cost overruns attributable to foundation difficulties or surprises during construction. These overruns may have been due to inadequate investigations, lack of adequate engineering effort to tailor the structures to site geology and topography, and/or contractual terms, among other considerations.; Challenges for hydropower foundations and opportunities for innovative technology solutions are identified in the following areas (consistent with the three main phases of foundation development): Geotechnical site assessment, Foundation design and materials, Construction methods and technology. Ultimately, this report aims to provide information about geotechnical foundations for low-head hydropower and to motivate transformative technologies to support hydropower growth.

13 HYDRO ENERGY↗