SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗

DOE OSTI2021

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Characterization of ultrafast-laser ablation of micro-structures in Li-ion battery anode and cathode materials: Morphology, rate, and efficiency

Characterization of the rate and quality of ultrafast-laser ablation of Li-ion battery (LIB) electrode materials is presented for a collection of common and next-generation electrodes. Laser ablated micro-structures on the surface of LIB electrodes have been shown to provide dramatic enhancement of high-rate capability and electrode wetting. However, industrial adoption is hampered by a lack of data enabling informed choice of laser parameters and predicting process throughput. Here this work bridges this gap by providing characterization of the ablation process at more laser parameters (laser fluence and number of pulses used) than are currently available in the literature. Further, we expand on previous graphite and lithium iron phosphate (LFP) ablation work by extending ablation characterization to new LIB materials, providing high data resolution and adopting new characterization metrics which are relevant for industrial application of this technology. Ablated pores are characterized by their ablated depth, volume, and how the depth and volume ablation rate changes as a function of pore depth. Finally, we provide a detailed characterization of the morphology of laser ablated micro-structures which informs how material and laser parameters affect the quality of laser-processed electrodes.

25 ENERGY STORAGE↗

DOE OSTI2023

The influence of submarine canyons-related processes on recent benthic foraminiferal distribution, Espírito Santo Basin, Southeastern Brazil

The complex topography of submarine canyons may result in different composition of benthic foraminifera assemblages. To understand how trophic, hydrological and sedimentological conditions in submarine canyons can influence the distribution of benthic foraminifera, and to use this information to corroborate paleoenvironmental interpretations for the Holocene, we investigated recent benthic foraminiferal assemblages (total fauna >63 μm) and sedimentological data in two canyons (CANWN and CAND) in the Espírito Santo Basin (ESB) between 18°20' and 21°20' S. Surface sediment samples (0–2 cm) were collected inside the canyons (150 to 1300 m water depth) and in shelf-slope adjacent transects (50 to 1300 m water depth). The density, taxonomic diversity, and assemblage composition of benthic foraminifera change with depth and location. The distinct ecological preferences of the most abundant taxa allowed us to recognize five benthic foraminiferal groups. Three groups (V, III, and I) are present in different bathymetric sectors; Group V: outer shelf (50 m), Group III: upper, and middle – lower slope (150 to 400 m), and Group I: middle – lower slope (1000 to 1300 m). Groups II and IV show no characteristic bathymetric distribution and are present only in CAND and in CANWN, respectively. Group V consists of Hanzawaia boueana, Peneroplis planatus, and Quinqueloculina lamarckiana; Group III is dominated by Globocassidulina rossensis and Trifarina spp.; Group I consists of Globocassidulina crassa, Bolivina lowmani, Gavelinopsis versiformis, Alabaminella weddellensis, and Epistominella exigua. The main species in Group II (CAND, 150, 1000 to 1300 m) are Trifarina angulosa, Globocassidulina subglobosa, and Discorbis vilardeboanus. Group IV (middle – lower CANWN, 1000 to 1300 m), consists mainly of agglutinated species Glomospira charoides, Rhabdammina abyssorum, and Psammosphaera fusca. Further, our data suggest that the quantity (and quality) of food supply, hydrodynamic conditions and sediment properties are the main drivers controlling the bathymetric distribution of benthic foraminiferal assemblages in both canyons. The middle – lower CANNW revealed unstable trophic conditions, related to terrigenous sediment input due to turbidity currents. In CAND, the foraminiferal assemblages ecology indicated sufficient organic matter supply that favors species establishment and diversity, indicating a more productive and less unstable environment than in CANWN.

58 GEOSCIENCES↗

DOE OSTI2023

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A methodology for customizing implantation profiles of light ions using a single thin foil energy degrader

We report a method was developed to quantify the spatial distribution and implantation depth of energy-degraded light ions with a thin foil rotating energy degrader for use during multiple ion beam irradiation. The methodology covers three physical phenomena: ions passing through the thin foil, ions travelling through the vacuum to the target, and ion implantation into the target, and accounts for the distribution of ions both in depth and in plane. The processes of energy straggling and scattering were calculated using SRIM. The effects of raster-scanning, and the geometry of the system were implemented in scripts handling the SRIM output files. Elastic backscattering (EBS) using 2.38 MeV H+ protons was used to measure the helium depth profiles after implantation with and without thin foil energy degradation. Defect analysis with transmission electron microscopy confirmed the implantation profiles measured with EBS and calculated with SRIM.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Characterization of Gulf Stream energy resources: A 30-year high-resolution simulation study

The strength and proximity of the Gulf Stream (GS) along the U.S. coast offers opportunities for marine energy extraction. Leveraging high-resolution (up to 400 m) long-term (30-year) hindcasts based on an unstructured grid model, this study characterizes the energy potential of the GS in two viable locations: the Florida Straits (FS) and Cape Hatteras (CH). Results reveal distinct characteristics in the kinetic energy flux (KEF) and kinetic energy density (KED) for the two regions. The 30-year mean KEF at CH is 27.50 GW, surpassing the 19.74 GW estimated for the FS, which is attributed to the wider cross-section of the GS in CH. Additionally, CH exhibits higher KEF variability due to GS’s larger path shifts and meanders. Estimated KEDs at ADCP site locations reach up to 2,908 W/m 2 in the FS and 1,512 W/m 2 in the CH at a depth of 20 m. High-energy zones are closer to shore in the FS, whereas energy-rich areas at CH are located further offshore. A depth-based KED analysis shows significant decreases across typical energy extraction depths (20 to 100 m), with reductions up to 44%. Finally, this study provides high-resolution resource characterization, forming a foundation for optimized site selection and marine hydrokinetic system design.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

Optimizing aluminum oxide passivation layers—Laser-induced plasmas for layer formation, real-time diagnostics, and layer evaluation

Molten salts have beneficial thermophysical and electrolytic properties, but the aggressive nature of molten salts requires corrosion mitigation strategies, such as structural material surface treatment. In this study, a laser-induced breakdown spectroscopy (LIBS) system was used to simultaneously ablate and monitor de-excitation spectra to form and identify dense alumina passivation layers. Additionally, high-frequency (kHz) LIBS imaging allowed rapid elemental mapping and depth profiling for surface O/Al ratios via minimally destructive analysis. The paired t-test rejected the null hypothesis that air and Ar cover gases were equivalent (p = 0.00012), and the analysis of variance (ANOVA) showed that overlap between shots (p = 0.03112) and cover gas flow rate (p = 0.02729) were statistically significant. At 75% overlap, the O/Al ratio increased to 10% ± 2%, and 90% overlap rose further to 31% ± 2%. Regardless of overlap, tuning the gas flow from 0.5 to 2.5 L min −1 enhanced the O/Al ratio by 24% ± 4%. LIBS depth profiling showed a 3× increase in layer thickness from 75% to 90% overlap. Molecular band peaks of AlO, a precursor to Al 2 O 3 found during treatment, suggest that band intensity could be used for real-time Al 2 O 3 optimization. Ultimately, by adjusting laser spot overlap and cover gas flow rate, a 266 nm laser was shown to form Al 2 O 3 layers on an Al6061 sample, with cover gas flow rate primarily affecting the O/Al ratio and layer thickness correlating with laser spot overlap. This combination of molecular spectra collection, imaging, and depth profiling demonstrated the robust layer analysis capabilities of LIBS.

Corrosion↗

DOE OSTI2022

Sink or source? Insights into the behavior of copper and zinc in the sediment porewater of a constructed wetland by peepers

The H-02 free water surface constructed wetland has been applied to remove heavy metals, mainly copper (Cu) and zinc (Zn), from wastewater on the Savannah River Site (Aiken, SC, USA). More and more studies focus on the metal behavior between the sediment and the overlying water, which directly reflects the stability of metals after sedimentation in constructed wetlands. This study focused on the biogeochemical pathways in metal bioavailability and remobilization in the sediment after metals were removed from the overlying water. The dialysis sampling devices (peepers) were used to collect porewater samples from eight depths in the sediment for the measurement of Cu, Zn, dissolved organic carbon (DOC), and major anions (sulfate and chloride). Surface water samples were also collected for the measurement of Cu, Zn, DOC, and anions. Different temporal trends were observed for dissolved Cu between the surface and bottom waters, but not for dissolved Zn. There were no obvious changes in porewater metal concentrations with increasing depths in the sediment. Sediment served as a sink for Cu as only 3% of porewater samples showed higher labile Cu concentrations than the surface water during the entire year. Furthermore, sediment served as a source for Zn in summer and winter as 32% of porewater samples showed higher labile Zn concentrations than the surface water, and these samples were collected at all sediment depths. We think the seasonal changes in the behaviors of Cu and Zn were primarily controlled by the sulfur dynamics and the metal removal processes in the constructed wetland, as well as the different complexing chemistry between Cu and Zn. Also, our study supports that peeper is a powerful tool for studying the biogeochemistry of metals in the sediment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Quasi-In-Situ Analysis of Electrode Top Atomic Layers via High-Sensitivity Low-Energy Ion Scattering and Potential-Controlled Sample Transfer

Electrocatalytic reactions involve interfacial interactions between the surfaces of electrodes and reactive species at an electrolyte interface. There are presently no universal or unambiguous methods to directly assay the active top atomic layer composition that influences the reactivity of these electrodes under relevant operating conditions. Low-energy ion scattering (LEIS) spectroscopy is a surface characterization technique that yields compositional analysis of the outermost atomic layer of a material, but it must be performed in ultrahigh vacuum (UHV). Application of LEIS measurements to electrochemical materials that are removed from ambient liquid-phase environments thus leaves an open question as to whether the surface that is transferred to UHV is truly the surface that manifested during the electrochemical reaction. Toward the goal of preserving the active surface state, we developed a sample transfer workflow for LEIS enabling air-free removal and drying of an electrode from an electrochemical cell while maintaining control of the potential using an auxiliary electrode. The potential-controlled emersion method was demonstrated to give distinct potential-dependent surface compositions for a Cu−Pd alloy relative to removal after uncontrolled return to open-circuit potential. A Cu-enriched surface was found at anodic potential and a Pd-enriched surface at cathodic potential, suggesting that the approach can be used to retain representative atomic configurations during transfer. Since adsorbates will often persist from the reaction environment, conventional sample pretreatment methods for removal, including atomic O and atomic H exposure, were also contrasted. Both methods were found to differ with results from incidental low-dose depth profiling by the LEIS primary ion source, which removes adventitious species and surface atoms during the course of repeated measurements. These depth profiles were found to be sensitive to sample history and thus qualitatively informative, despite the possible changes induced by ion damage. The results exhibit (i) the need for complete control over the polarization state of the sample at all times (no excursions to open circuit during transfer) and (ii) the utility of low-dose depth profiling to capture changes in the near-surface composition.

Alloys↗

DOE OSTI2023

Molecular and Dual-Isotopic Profiling of the Microbial Controls on Nitrogen Leaching in Agricultural Soils under Managed Aquifer Recharge

Nitrate (NO 3 – ) leaching is a serious health and ecological concern in global agroecosystems, particularly those under the application of agricultural-managed aquifer recharge (Ag-MAR); however, there is an absence of information on microbial controls affecting NO 3 – leaching outcomes. We combine natural dual isotopes of NO 3 – ( 15 N/ 14 N and 18 O/ 16 O) with metagenomics, quantitative polymerase chain reaction (PCR), and a threshold indicator taxa analysis (TITAN) to investigate the activities, taxon profiles, and environmental controls of soil microbiome associated with NO 3 – leaching at different depths from Californian vineyards under Ag-MAR application. The isotopic signatures demonstrated a significant priming effect (P < 0.01) of Ag-MAR on denitrification activities in the topsoil (0–10 cm), with a 12–25-fold increase of 15 N–NO 3 – and 18 O–NO 3 – after the first 24 h of flooding, followed by a sharp decrease in the enrichment of both isotopes with ~80% decline in denitrification activities thereafter. In contrast, deeper soils (60–100 cm) showed minimal or no denitrification activities over the course of Ag-MAR application, thus resulting in 10–20-fold of residual NO 3 – being leached. Metagenomic profiling and laboratory microcosm demonstrated that both nitrifying and denitrifying groups, responsible for controlling NO 3 – leaching, decreased in abundance and potential activity rates with soil depth. TITAN suggested that Nitrosocosmicus and Bradyrhizobium, as the major nitrifier and denitrifier, had the highest and lowest tipping points with regard to the NO 3 – changes (P < 0.05), respectively. Overall, our study provides new insight into specific depth limitations of microbial controls on soil NO 3 – leaching in agroecosystems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Reorganization of Water at Aqueous Aluminum Chloride (AlCl 3 ) Interfaces: Vibrational Sum Frequency Generation and Molecular Dynamics Simulations

AlCl 3 hydration states and complexation are not well understood both in solutions and at the air–aqueous interface despite their potential significance in natural waters and their industrial and energy-related applications. Here, we investigated Al 3+ and Cl – ion behaviors in an AlCl 3 aqueous bulk solution and at the air–aqueous interface using interface-selective vibrational sum frequency generation (SFG), Raman and infrared spectroscopies, molecular dynamics (MD) simulation, as well as molecular-informed reduced modeling. Our reduced modeling reveals relatively long-range effects for Al 3+ as compared to monovalent ions such as Na + indicating that the interfacial depth of trivalent ions can be significantly larger than that of monovalent ions at the air–water interface. MD simulations reveal interfacial stratification and multiple layering of the ions. Compression of the Al 3+ and Cl – distributions with increasing concentrations from 0.5 to 2.5 m is also observed in the subsurface regions. Significant SSP- and PPP-polarized SFG OH spectral intensity increases are observed from 0.5 to 1.5 m and 0.5 to 2.5 m, respectively, indicative of interfacial depth increases and a change in average orientation above 1.5 m. Extensive evaluation of SFG spectra, Fresnel-corrected using several approaches, shows the same trends. The nonmonotonic trend points to a changing structure in surface and subsurface water orientation and hydrogen bonding environment generally consistent with the MD simulation of stratification and water orientation changes. Furthermore, solvent-shared ion pairing is implicated with MD simulation radial distribution analysis and consistent with infrared spectral identification of the hexaaqua aluminum ion in the solution phase. Spectral evidence of a strong Al 3+ hydration shell and the acidic behavior of the Al 3+ ions is obvious in the Raman and infrared spectra of the bulk solution. In conclusion, we show that the MD dipole potential is directly related to the MD second-order susceptibility of the interface, χ SFG–MD (2) , both of which correlate up to ∼35 Å with the spectral observations of increasing and then saturating intensities, suggesting that both ion stratification and interfacial depth determine the water orientations at an air–water interface of 1-3 electrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Quantitative Measurement of Li-Ion Concentration and Diffusivity in Solid-State Electrolyte

Here, we present a quantitative method to measure the concentration and diffusivity of Li ion in a solid-state electrolyte with nanoscale depth resolution. We designed a standard sample with different depths of trenches and used time-of-flight secondary ion mass spectroscopy, inductively coupled plasma optical emission spectroscopy, three-dimensional (3D) optical microscopy, and electrochemical strain microscopy to obtain effective Li-ion concentration and the corresponding diffusivity as a function of depth. As a result, we verified the dependence of electromechanical strain induced by ionic oscillation on mobile-ion concentration. We also showed the diffusivity change as a function of effective concentration and discussed the difference in diffusivity near the surface and inside the bulk. Our method can be widely used for electromechanical strain in various ionic conductors, rendering our technique universal for quantitative analysis of ionic dynamics at a multiscale.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions

Cloud droplet spectral relative dispersion is critical to parameterizations of cloud radiative properties, warm-rain initiation, and aerosol-cloud interactions in models; however, there is no consistent relationship between relative dispersion and volume-mean radius in literature, which hinders improving relative dispersion parameterization and calls for physical explanation. Here we show, by analyzing aircraft observations of cumulus clouds during Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations, that the correlation between relative dispersion and volume-mean radius changes from positive to negative as volume-mean radius increases. With the new observation, we postulate that the sign of the correlation is determined by whether or not condensation (evaporation) occurs simultaneously with significant new activation (deactivation). The hypothesis is validated by simulations of both an adiabatic cloud parcel model and a parcel model accounting for entrainment-mixing. A new quantity, first bin strength, is introduced to quantify this new observation. Theoretical analysis of truncated gamma and modified gamma size distributions further supports the hypothesis and reconciles the contrasting relationships between relative dispersion and volume-mean radius, including the results in polluted fog observations. The results could shed new light on the so-called “twilight zone” between cloudy and cloud-free air, which in turn affects evaluation of aerosol-cloud interactions and retrieval of aerosol optical depth.

Lu, Chunsong↗

DOE OSTI2020

Carbon Thaw Rate Doubles When Accounting for Subsidence in a Permafrost Warming Experiment

Permafrost soils, which are perennially frozen soils found throughout cold regions, contain vast quantities of carbon and ice. When permafrost thaws, carbon can be lost to the atmosphere, contributing to climate change. This means it is important to track permafrost thaw, which is often done using active layer thickness, or the depth of the seasonally thawed surface layer of soil. However, ice volume can be lost from thawing permafrost, causing the soil surface to drop. Conventional measurements do not account for this surface drop, and the rate of thaw could therefore be underestimated. We found that experimentally warmed soils dropped at a rate of 6 cm year-1, mostly due to loss of ice volume and also due to the loss of soil mass. When accounting for the change in soil surface height over time, the full depth of permafrost thaw was 49% greater. The increased depth of thaw resulted in more than twice as much carbon being thawed as was estimated with standard methods that did not account for subsidence. These findings suggest that permafrost is thawing more quickly than long-term records indicate and that this could result in additional carbon release contributing to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Quantifying Dynamic Water Storage in Unsaturated Bedrock with Borehole Nuclear Magnetic Resonance

Quantifying the volume of water that is stored in the subsurface is critical to studies of water availability to ecosystems, slope stability, and water-rock interactions. In a variety of settings, water is stored in fractured and weathered bedrock as rock moisture. However, few techniques are available to measure rock moisture in unsaturated rock, making direct estimates of water storage dynamics difficult to obtain. Here, we use borehole nuclear magnetic resonance (NMR) at two sites in seasonally dry California to quantify dynamic rock moisture storage. Additionally, we show strong agreement between NMR estimates of dynamic storage and estimates derived from neutron logging and mass balance techniques. The depths of dynamic storage are up to 9 m and likely reflect the depth extent of root water uptake. To our knowledge, these data are the first to quantify the volume and depths of dynamic water storage in the bedrock vadose zone via NMR.

58 GEOSCIENCES↗

DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗