Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “soil depth”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Optimal load attachment of a deeply embedded ring anchor in clay

A Deeply Embedded Ring Anchor (DERA) system has been developed as a cost-effective solution for mooring arrays of floating offshore wind turbines (FOWTs) to the seabed. The DERA boasts several key features, including its versatility in various soil types, compact size, compatibility with diverse mooring systems, multi-line potential, and robust performance even under unintentional loading conditions. While prior preliminary studies have provided valuable insights into how the DERA can enhance cost-effectiveness by offering a high load capacity, these studies have predominantly focused on optimizing anchor performance under translational horizontal and vertical loading. However, to design the DERA optimally, we must also consider its ability to handle inclined loading conditions in addition to lateral and axial loadings. Due to its shorter length compared to a conventional caisson, the DERA has less resistance to moments, making it more sensitive to horizontal load capacity and the optimal load attachment depth concerning load angle. For this reason, our study introduces an analytical approach to evaluate the effects of inclined loading on anchor performance, utilizing the previously validated upper bound plastic limit analysis (PLA) method. In investigating the optimal load attachment of the DERA, this paper conducts a parametric study to analyze how factors such as load attachment depth, anchor aspect ratio, and load inclination affect the DERA’s load capacity. Our findings indicate that PLA can serve as a valuable analytical tool for assessing the ultimate load capacity of the DERA, particularly under inclined loading conditions.

42 ENGINEERING↗
DOE OSTI2022

An approach to sulfate geoengineering with surface emissions of carbonyl sulfide

Sulfate geoengineering (SG) methods based on lower stratospheric tropical injection of sulfur dioxide (SO 2 ) have been widely discussed in recent years, focusing on the direct and indirect effects they would have on the climate system. Here a potential alternative method is discussed, where sulfur emissions are located at the surface or in the troposphere in the form of carbonyl sulfide (COS) gas. There are two time-dependent chemistry–climate model experiments designed from the years 2021 to 2055, assuming a 40 Tg-Syr -1 artificial global flux of COS, which is geographically distributed following the present-day anthropogenic COS surface emissions (SG-COS-SRF) or a 6 Tg-Syr -1 injection of COS in the tropical upper troposphere (SG-COS-TTL). The budget of COS and sulfur species is discussed, as are the effects of both SG-COS strategies on the stratospheric sulfate aerosol optical depth (~Δτ=0.080 in the years 2046–2055), aerosol effective radius (0.46 µm), surface SO x deposition (+8.9 % for SG-COS-SRF; +3.3 % for SG-COS-TTL), and tropopause radiative forcing (RF; ~-1.5 W m –2 in all-sky conditions in both SG-COS experiments). Indirect effects on ozone, methane and stratospheric water vapour are also considered, along with the COS direct contribution. According to our model results, the resulting net RF is –1.3 W m –2 , for SG-COS-SRF, and –1.5 W m –2 , for SG-COS-TTL, and it is comparable to the corresponding RF of –1.7 W m –2 obtained with a sustained injection of 4 Tg-Syr -1 in the tropical lower stratosphere in the form of SO 2 (SG-SO2, which is able to produce a comparable increase of the sulfate aerosol optical depth). Significant changes in the stratospheric ozone response are found in both SG-COS experiments with respect to SG-SO2 (~5 DU versus +1.4 DU globally). According to the model results, the resulting ultraviolet B (UVB) perturbation at the surface accounts for –4.3 % as a global and annual average (versus –2.4 % in the SG-SO2 case), with a springtime Antarctic decrease of –2.7 % (versus a +5.8 % increase in the SG-SO2 experiment). Overall, we find that an increase in COS emissions may be feasible and produce a more latitudinally uniform forcing without the need for the deployment of stratospheric aircraft. However, our assumption that the rate of COS uptake by soils and plants does not vary with increasing COS concentrations will need to be investigated in future work, and more studies are needed on the prolonged exposure effects to higher COS values in humans and ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Peatland hydrology across scale: a probabilistic framework for confronting variability, heterogeneity, and uncertainty (Final Report)

Our project aimed to guide the representations of peatland hydrology and water-carbon feedbacks within Earth System Models, by improving process understanding and developing parsimonious models that elucidate the effects of hydroclimatic variability and spatial heterogeneity. Our research questions centered on the three key drivers of peatland hydrology: (1) seasonal and interannual hydroclimatic fluctuations, (2) spatial heterogeneity of peatland microtopography, and (3) hydrological connectivity across landscape units within a peatland watershed. Our analyses and model testing used existing long-term datasets within the Marcell Experimental Forest (MEF), with new data collected whenever necessary. To date, we have demonstrated the strong influence of water table elevations on the temperature sensitivity of CH 4 emissions, using a newly developed eddy covariance dataset spanning eleven years at MEF. Specifically, higher water tables dampen the increase in CH 4 emissions in the spring as well as their decrease in the fall, resulting in hysteresis. These results imply that any hydroclimatological changes in peatlands that shift seasonal water availability from winter to summer will increase annual CH 4 emissions, even if soil temperature remains unchanged. To further investigate water table elevation change within the peatland microtopography, we installed automated water table gauges across four bog-forest transects (the “lagg”) in two peatland watersheds at MEF, which are hotspots of intense biogeochemical activity. These measurements will give us new information about the extent of lagg expansion and contraction during snowmelt and high intensity rainfall events, as well as the directionality of water and nutrient flow into the bog from the surrounding forests. Finally, we investigated the hydrological connectivity across the peatland watershed complex during early spring – a critical period where water table is elevated from snowmelt – using extended hydrological records from MEF (e.g., snow and frost depths, water table elevations, streamflow). The importance of thermal and wintertime processes within the Energy Exascale Earth System Model (E3SM) was determined through sensitivity analyses. Results from data analyses show that (i) streamflow has decreased over decades likely due to increased evapotranspiration rates (despite no detectable trends in precipitation), (ii) spring streamflow generation is controlled by a sequence of fill-and-spill pathways from the snowpack to the stream through peatland water stores, and (iii) frost depth is a key explanatory variable for the timing and magnitude of streamflow. These results suggest that frost plays an important role in connecting surface water storage in peatlands to stream outlets, and that the hydrological connectivity across the peatland watershed complex can mediate the sensitivity of hydrological responses to climate variations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Origin of tungsten and geochemical controls on its occurrence and mobilization in shallow sediments from Fallon, Nevada, USA

In this study, tungsten (W) occurrence and speciation was investigated in sediments collected from Fallon, Nevada where previous studies have linked elevated W levels in human body fluids to an unusual cluster of childhood leukemia cases. The speciation of sedimentary W was determined by μ-XRF mapping and μ-XANES. The W content of the analyzed surface sediments ranged between 81 and 25,908 mg/kg, which is significantly higher than the W content in deeper sediments which ranged from 37 to 373 mg/kg at 30 cm depth. The μ-XANES findings reveal that approximately 20–50% of the total W in the shallow sediment occurs in the metallic form (W 0 ); the rest occurs in the oxide form (W VI O 3 ). Because W 0 does not occur naturally, its elevated concentrations in surface sediments point toward a possible local anthropogenic origin. The oxidation of metallic W 0 with meteoric waters likely leads to the formation of W VI O 3 . The chief water-soluble W species was identified as WO 4 2– by chromatographic separation and speciation modeling. These results led us to postulate that W 0 particles from a currently unknown but local source(s) is (are) deposited onto the soils and/or surface sediments. The W 0 in interaction with meteoric water is oxidized to W VI O 3 , and as these sediment-water interactions progress, WO 4 2– is formed in the water at pH ~7. Under pH < 7, and sufficient W concentrations, tungstate tends to polymerize, and polymerized species are less likely to adsorb onto sediments. Polymerized species have lower affinity than monomers, which leads to enhanced mobility of W.

58 GEOSCIENCES↗
DOE OSTI2025

Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data

ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.

Ogle, Kiona [School of Informatics, Computing, and↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2022

Impacts of climate change on future hurricane induced rainfall and flooding in a coastal watershed: A case study on Hurricane Harvey

The warming climate is likely to increase hurricane-associated extreme rainfall and lead to sea-level rise (SLR). Thus, how the floods induced by intense hurricanes respond to these potential changes is of great concern. This study investigates the future warmer climate impacts on hurricane-induced extreme rainfall, and—more importantly—the subsequent compound flooding at the watershed scale (from an event-based analysis perspective). To this goal, a modeling framework is designed based on the Distributed Hydrology Soil Vegetation Model (DHSVM), the Two-Dimensional fvand the Regional Community Earth System Model (R-CESM). The framework was applied to Hurricane Harvey (2017) at the Clear Creek watershed (a coastal watershed in the southern Houston) as a case study. The results show that the projected maximum rainfall totals over the watershed would be exacerbated by 17.7 % and 49.7 % in the 2050s and 2090s (respectively) under Representative Concentration Pathway 8.5 (RCP 8.5). This means a 16.1 % increase in Harvey rainfall over the watershed per degree Celsius increase in Mean Surface Temperature over the Gulf of Mexico region (18°~31° N, 77°~98° W). Meanwhile, the increases in maximum inundation extent would be 11.0 % (2050s) and 19.5 % (2090s). Furthermore, considerable increases in maximum inundation depth and duration in regions along the middle and downstream of Clear Creek (and also those around Clear Lake) are expected. The projected SLR will have little effect on the maximum inundation depth and extent if storm surge changes are not taken into account; meanwhile, it will influence the inundation duration at downstream locations. In conclusion, this modeling framework can be also applied at other coastal watersheds to evaluate the projected climate change impacts on the compound flooding induced by extreme climate events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil carbon change in intensive agriculture after 25 years of conservation management

Changes in soil organic carbon (SOC) and nitrogen (SON) are strongly affected by land management but few long-term comparative studies have surveyed changes throughout the whole soil profile. We quantified 25-year SOC and SON changes to 1 m in 10 replicate ecosystems at an Upper Midwest, USA site. We compared four annual cropping systems in maize (Zea mays)-soybean (Glycine max)-winter wheat (Triticum aestivum) rotations, each managed differently (Conventional, No-till, Reduced input, and Biologically based); in three managed perennial systems (hybrid Poplar (Populus × euramericana), Alfalfa (Medicago sativa), and Conifer (Pinus spp.); and in three successional systems (Early, Mid- and Late succession undergoing a gradual change in species composition and structure over time). Both Reduced input and Biologically based systems included winter cover crops. Neither SOC nor SON changed significantly in the Conventional or Late successional systems over 25 years. All other systems gained SOC and SON to different degrees. SOC accrual was fastest in the Early successional system (0.8 ± 0.1 Mg C ha –1 y –1 ) followed by Alfalfa and Conifer (avg. 0.7 ± 0.1 Mg C ha –1 y –1 ), Poplar, Reduced input, and Biologically based systems (avg. 0.4 ± 0.1 Mg C ha –1 y –1 ), and Mid-successional and No-till systems (0.3 and 0.2 Mg C ha –1 y –1 , respectively). Over the most recent 12 years, rates of SOC accrual slowed in all systems except Reduced input and Mid-successional. There was no evidence of SOC loss at depth in any system, including No-till. Rates of SON accrual ranged from 64.7 to 0.8 kg N ha –1 y –1 in the order Alfalfa ≥ Early successional > Reduced input and Biologically based ≥ Poplar > No-till and Conifer > Mid-successional systems. Pyrogenic C levels in the Conventional, Early, and Late successional systems were similar despite 17 years of annual burning in the Early successional system (~ 15 % of SOC to 50 cm, on average, and ~40 % of SOC from 50 to 100 cm). Results underscore the importance of cover crops, perennial crops, and no-till options for sequestering whole profile C in intensively managed croplands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Code associated with "Assessing the Effect of a Deep-Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights from a Multi-Site US Study"

R code associated with "Assessing the Effect of a Deep-Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights from a Multi-Site US Study" These scripts can be used to analyze the accompanying dataset: https://doi.org/10.5281/zenodo.16620529. (1) resampling_functions.R : Creates functions for resampling data for nonparametric bootstrapping (2) bootstrap_depth_profiles.R : Applies resampling functions to the dataset (3) gapfill_masscalcs.R : Prepares data for equivalent mass calculations (4) get_14C_EM.R : Applies soilR to obtain the 14C end member for mixing model calculations (5) mass_aggregate.R : Performs equivalent mass and mass weighted averaging calculations.

Carbon↗
DOE OSTI2020

Multi-model Hydroclimate Projections for the Alabama-Coosa-Tallapoosa River Basin in the Southeastern United States

Abstract This study uses a high-resolution, process-based modeling framework to assess the impacts of changing climate on water resources for the Alabama-Coosa-Tallapoosa River Basin in the southeastern United States. A 33-member ensemble of hydrologic projections was generated using 3 distributed hydrologic models (Precipitation-Runoff Modeling System, Variable Infiltration Capacity, and Distributed Hydrology Soil Vegetation Model) of different complexity. These hydrologic models were driven by dynamically downscaled and bias-corrected future climate simulations from 11 Coupled Model Intercomparison Project Phase 5 global climate models under Representative Concentration Pathway 8.5 emission scenario, with 40 years each in baseline (1966–2005) and future (2011–2050) periods. The hydroclimate response, in general, projects an increase in mean seasonal precipitation, runoff, and streamflow. The high and low flows are projected to increase and decrease, respectively, in general, suggesting increased likelihood of extreme rainfall events and intensification of the hydrologic cycle. The uncertainty associated with the ensemble hydroclimate response, analyzed through an analysis of variance technique, suggests that the choice of climate model is more critical than the choice of hydrologic model for the studied region. This study provides in-depth insights of hydroclimate response and associated uncertainties to support informed decisions by water resource managers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Biogeochemistry of Pond B (Savannah River Site, South Carolina, USA): Water column and Sediments

Pond B at Savannah River Site (SRS, South Carolina) is a monomictic reservoir that received SRS R reactor cooling water from 1961–1964. Previous studies conducted between the 1980s–1990s on the water column and sediments of Pond B measured trace amounts of Pu (33 MBq 238Pu and 430 MBq 239,240Pu), 241Am, and 137Cs. Since then, the pond has been relatively isolated and the radionuclide concentrations have not been monitored over time. Herein, about 30 years after the last publication on Pond B, we are re-evaluating the geochemistry and radionuclide distribution within Pond B at five locations along a horizontal transect from the inlet to outlet. In addition, we are conducting the first analysis of the microbial community composition. This data package consists of water column measurements (1) using a multi-probe sonde, (2) plutonium isotopes and total concentrations, (3) total organic/inorganic carbon concentrations, (4) trace metals and major ions concentrations, and (6) microbial community composition. Sediment measurements include (1) plutonium isotopes and total concentrations and (2) microbial community composition.For sediment core and porewater, please see other Pond B data package (ess-dive-60d2352a4495472-20200615T203925174).The data in this package will be submitted for peer review in 2023.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

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

Molecular Imaging of Microbially Induced Corrosion of Synthetic Archeological Glasses by a Rhizosphere Bacterium

Microbially induced corrosion (MIC) focuses on the degradation of solid materials, such as glass or metal. Soil microbes are often associated with the corrosion of foreign objects in the rhizosphere. Paenibacillus polymyxa SCE2, a facultative anaerobic bacterium in soil, is of the same genus as bacteria found near nuclear waste disposal sites. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) was used for imaging surface changes induced by P. polymyxa SCE2 cultured on two synthetic glass coupons to represent natural analogs of materials that were studied in relation to the vitrification of nuclear waste. Multimodal imaging was used to verify bacterial coverage across the glass surface after long-term growth. ToF-SIMS spectral analysis showed detection of glass component ions, such as silicon oxide (m/z – 59.96 SiO 2 – ) and aluminum oxide (m/z – 101.95 Al 2 O 3 – ), and biofilm’s extracellular polymeric substance (EPS) components, such as pentadecanoic acid (m/z – 241.22 C 15 H 29 O 2 – ) and sterol lipids (m/z – 311.16 C 20 H 23 O 3 – ). ToF-SIMS spectral, imaging, and depth profiling analyses showed that the glass rich in silica and other light elements (“granite glass”) had more “corrosion related” peaks than the glass that was less silica-rich and contained more iron (“dike glass”). Furthermore, these surface and interface compositional and spatial differences observed in the mass spectra and imaging were attributed to bacterial metabolism and an electron transfer mechanism influenced by morphological and compositional differences between the two types of glasses. ToF-SIMS is effective in studying microbial effects, bringing new molecular insights into MIC in a broader context of materials degradation.

Amorphous materials↗
DOE OSTI2025

Old Woman Creek Wetland Sediment and Electrochemical Sensor Microbial Community, 2023

We are developing a technique to monitor microbiological activities referred to as zero resistance ammetry, which entails the deployment of graphite electrodes in sediments. Measurement of current between electrodes of contrasting redox regimes and/or predominant terminal electron accepting processes can be used as an indicator of the extents of microbiological activity. We deployed an electrode array at depths of 2 mm, 4 mm, 76 mm, 78 mm, 152 mm, 154 mm, 227 mm, and 229 mm below the wetland sediment water interface in the Old Woman Creek National Estuarine Research Center, Huron, OH, USA (Lat. = 41.380833, Long. = -82.508889). A core was collected from adjacent sediment and subsamples were collected from depth intervals of 0 – 25 mm, 25 – 127 mm, 127 – 128 mm, and below 178 mm. To determine if the microbial communities attached to the electrodes were reflective of the adjacent sediment-associated microbial community, we conducted a 16S rRNA gene-based (V4 region) survey of these respective materials. This data package contains the results of these surveys, including metadata on the depths from which samples were collected (samples.csv), DNA extraction and sequencing information (OWC_DEPTH_AMPLICON_SEQUENCING_METADATA), sequence processing information (OWC_DEPTH_BIOINFORMATIC_METADATA.csv), an operational taxonomic unit (OTU) table (OWC_DEPTH_97OTUS_TABLE.csv), and nucleotide sequences of OTUs (OWC_DEPTH_97OTUS_SEQS.fasta). All files can be opened using a text-editing application. The fasta file is compatible with bioinformatics applications.

54 ENVIRONMENTAL SCIENCES↗