Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geoscience”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Onset of magma ocean solidification on Mars inferred from Mn-Cr chronometry

The mantle of Mars probably differentiated through the crystallization of a magma ocean during the first tens of million years (Ma) of Solar System evolution. However, the exact timescale of large-scale silicate differentiation of the martian mantle is debated, and in particular, it remains unclear when differentiation commenced. In this paper, we applied the short-lived 53 Mn- 53 Cr system to martian meteorites in order to date the onset of large-scale mantle differentiation on Mars. The new Cr isotope data demonstrate that martian meteorites exhibit no resolvable radiogenic 53 Cr variations, and instead have a uniform +20.2±1.2 (95% conf.) parts-per-million excess in 53 Cr/ 52 Cr relative to the terrestrial mantle. The investigated groups of martian meteorites are lithologically varied and derive from diverse mantle sources that probably had variable Mn/Cr. Hence, the lack of 53 Cr variability among martian meteorites demonstrates that silicate differentiation on Mars occurred after the extinction of 53 Mn. Provided that the sources of the martian meteorites have Mn/Cr variations that are typical of the terrestrial planets, this result implies that the onset of large-scale silicate differentiation must have occurred later than 20±5 Ma after Solar System formation. The onset of silicate differentiation on Mars inferred here is significantly later than time estimates for segregation of the martian core which conservatively occurred within <10 Ma after Solar System formation. Thus, the new Mn-Cr data imply that there was a small, but resolvable, time gap of at least 5 Ma between core formation and magma ocean solidification on Mars. If the age of core segregation is taken at face value, our results imply that the martian magma ocean remained mostly molten over several Ma. This inferred longevity of the magma ocean is inconsistent with thermal models predicting rapid (<1 Ma) solidification of the martian magma ocean. Although there is currently no unique solution to this conundrum, our results can potentially be explained by a protracted history of impact bombardment that delayed differentiation in a shallow magma ocean on Mars, or perhaps more readily, by the presence of an early and dense atmosphere that acted as an insulator and prevented the magma ocean from cooling quickly.

58 GEOSCIENCES↗
DOE OSTI2022

Tungsten and molybdenum isotopic evidence for an impact origin of pallasites

The origin of pallasites—stony-iron meteorites mainly composed of olivine and Fe-Ni metal—is debated and proposed formation scenarios broadly range from models that explain pallasite formation by internal processes in the mantle of a differentiated planetesimal to those that involve impact–induced mixing of core and mantle materials. Here, the origin of pallasites is examined by studying the nebular source regions of their precursor material using Mo isotopes and their history of metal-silicate segregation using Hf-W chronometry. We report new Mo and W isotopic data for a large suite of pallasite metal samples, alongside Pt isotope data to quantify superimposed cosmic ray exposure effects. Most main-group pallasites exhibit uniform pre-exposure 182 W and Mo isotopic compositions that bear an excellent similarity to those of IIIAB iron meteorites. Four main-group pallasites and the IIIAB iron Thunda have more radiogenic pre-exposure 182 W compositions, but display the same Mo isotopic composition as other main-group pallasites and IIIAB irons. This strong chronological and genetic link strongly suggests that main-group pallasite metal originated in the IIIAB parent body core. This, combined with prior Pd-Ag chronometric evidence for an early collisional disruption of the IIIAB parent body, implies that main-group pallasites formed by impact–induced mixing of metal and silicates rather than by an internal process on the IIIAB parent body. This mixing led to elevated 182 W compositions in some pallasites, which are best accounted for by partial re-equilibration of IIIAB metal with radiogenic 182 W from the colliding body. Altogether, our results support models that explain main-group pallasite formation by injection of pallasite metal into the mantle of another differentiated body, implying that pallasite silicates did not primarily derive from the IIIAB mantle, but instead from that of the colliding body.

58 GEOSCIENCES↗
DOE OSTI2020

Citation for the 2019 Victor M. Goldschmidt award to Donald J. DePaolo

Perhaps the comments most commonly heard when a colleague receives professional recognition are “well-deserved” and “long-overdue”. As such, it is my pleasure to introduce Donald DePaolo as the recipient of this year’s Goldschmidt Award as it is both long-overdue and so well-deserved. While Don has received many honors over the course of his career, the Goldschmidt best recognizes his broad impact on geochemistry. There are few who have contributed as much to geochemistry as Don DePaolo. Finally, his scientific contributions are voluminous and broad, ranging from pioneering efforts to establish the Sm/Nd isotopic system, and its applications to mantle evolution and basalt petrogenesis, to studies of magma mixing, assimilation and crystal fractionation, to evolution of the lower crust and granite petrogenesis, to chemical oceanography, diagenesis and seafloor alteration, to reactive transport accompanying rock-water interaction, and to kinetic isotopic fractionation associated with crystal growth, diffusion and evaporation.

58 GEOSCIENCES↗
DOE OSTI2022

Cadmium isotope systematics in sedimentary carbonate: Extending the utility of the cadmium isotope palaeo-productivity proxy

The cadmium (Cd) isotope compositions of seven Holocene (<11.7 ka) carbonate sediments, including the bulk carbonate, as well as the > 140 μm (zooplankton), <63 μm and < 20 μm (phytoplankton) fractions, and individual species of foraminifera isolated from the bulk carbonate were obtained for sub-samples processed using the Cd-cleaning and Mg/Ca cleaning methods. The respective Cd isotope systematics of each carbonate archive were then compared with the Cd isotope signature within the ocean using modern seawater data. We show that a constant Cd isotope fractionation factor of -0.25 ± 0.03 ‰ (2SE, n = 4) describes the partitioning of Cd from seawater into the carbonate phase ($Δ^{114}Cd_{CdC O_3 - seawater}$), allowing for accurate reconstructions of the past seawater nutrient regime. The Cd isotope signature of the < 63 µm and < 20 µm (phytoplankton) fractions can display comparable values to those recorded in the modern surface ocean, and can be used to model the magnitude of Cd uptake and utilisation in the past ocean. Individual species of planktic and benthic foraminifera, which inhabit different depths of the water column, can provide a full surface to deep water column reconstruction of the Cd nutrient regime, while infaunal benthic foraminifera can be used to interrogate the Cd isotope systematics of sediment porewaters. Bulk sediment and the > 140 µm (zooplankton) fraction also retain important information relating to the overlying water column Cd isotope systematics and can serve as essential archives for palaeo-productivity reconstructions where the separation of individual components from the bulk sediment is not possible. This new Cd isotope palaeo-productivity toolkit can considerably aid reconstruction efforts and vastly improves the ability to quantify the role of phytoplankton in modulating carbon sequestration through transitions in Earth’s climate evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Chronology of alkali anorthosite 14304 clast “b” records basin forming impact at ~3.95 Ga

Alkali-suite rocks constitute one of three major suites of lunar crustal rocks. As such, constraining their formation timescales and petrogenesis is important for understanding the earliest magmatic history of the Moon. However, the magmatic history of alkali-suite rocks is partly obscured by superimposed effects of major basin-forming impact events on the lunar nearside. Consequently, unambiguous crystallization ages of samples from this suite of rocks have not been determined. Here, the aim of this study is to better understand the petrogenetic history of the alkali-suite and the potential superimposed effects of impact metamorphism by determining Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages for an alkali anorthosite clast from Apollo 14 lithic breccia 14304 termed clast “b”. The new chronologic measurements of clast “b” yield concordant Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages of 3947 ± 13 Ma, 3975 ± 34 Ma, and 3937 ± 37 Ma respectively, resulting in a weighted mean age of 3949 ± 11 Ma. This age is not interpreted to date an igneous event related to production of the lunar highlands crust and instead the chronology and petrography of clast “b” are most readily explained by an impact event at ~3.95 Ga that caused near-complete isotopic re-equilibration of the Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar chronometers. The weighted mean age of 3949 ± 11 Ma of clast “b” is several hundred Ma younger than 4.3–4.4 Ga ages typically determined for samples of other crustal rock suites but in very good agreement with independent estimates for the formation of Imbrium basin ejecta and other marginally older impact events which are thought to have been sampled at the Apollo 14 landing site. Thus, although petrologic and geochemical examination suggest that clast “b” is a pristine igneous clast, its age likely records an impact event at the Apollo 14 landing site. Whereas the various determined ages of clast “b” do not reflect the timescales of alkali-suite magmatism, the relatively low initial Sr and Nd isotopic compositions of clast “b” indicate that its protolith evolved with very low 147 Sm/ 144 Nd and 87 Rb/ 86 Sr that are distinct from estimates for urKREEP but similar to that of lunar plagioclase. This implies that the igneous protolith of clast “b” derived from a plagioclase-dominated KREEP-rich source that must have formed after the formation of the urKREEP source at ~4.35 Ga but well before the impact event recorded by clast “b” at ~3.95 Ga.

58 GEOSCIENCES↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2026

The solubility of NdPO 4 and DyPO 4 and stability of Nd and Dy chloride and hydroxyl complexes as a function of pH and salinity to 450 °C

The mobility of rare earth elements (REE) in geological systems is often controlled by the stability of monazite and xenotime, which are important hosts for the light and heavy REE. These REE phosphates constitute important resources and provide information on ore formation conditions and timing due to their uses as geothermometers and geochronometers. While the thermodynamic properties for these minerals are well-established up to high temperature and pressure, the properties for the aqueous REE complexes are not well constrained with limited information up to 300 °C. In this study, the solubility of synthetic NdPO 4 and DyPO 4 endmembers were measured in hydrothermal sub- to supercritical NaCl-bearing aqueous solutions from 350 at saturated water vapor pressure to 450 °C at 700 bar. The speciation of Nd and Dy was determined in acidic to alkaline solutions (pH 25 °C values from 2 to 10), and indicates that the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 predominate at low salinity (0.01 mol/kg NaCl) with some contribution of REE chloride species REECl 2+ and REECl 2 + in acidic fluids. The REE phosphate solubility measured in these experiments is up to two orders of magnitude higher than predicted using existing thermodynamic properties from literature for aqueous species extrapolated from lower temperature data to supercritical conditions. To address this discrepancy, the thermodynamic properties of the REE aqueous species were optimized using GEMSFITS to derive the formation constants for the REE chloride (β Cl ) and hydroxyl (β OH ) complexes in the studied temperature range. This study highlights the importance of the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 over a wide range of pH and temperature. As a result, the revised thermodynamic properties for Dy and Nd chloride and hydroxyl species more accurately predict the solubility of NdPO 4 and DyPO 4 and provide insight into the speciation and mobility of the light and heavy REE in hydrothermal supercritical fluids in the crust.

58 GEOSCIENCES↗
DOE OSTI2025

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗
DOE OSTI2026

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2025

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Contrasting Effects of Catecholate and Hydroxamate Siderophores on Molybdenite Dissolution

Molybdenum (Mo) is essential for many enzymes but is often sequestered within minerals, rendering it not readily bioavailable. Metallophores, metabolites secreted by microorganisms and plants, promote mineral dissolution to increase metal bioavailability. However, interactions between metallophores and Mo-bearing minerals remain unclear. Here, in this study, catecholate protochelin and hydroxamate desferrioxamine B (DFOB) were utilized to examine their effects on dissolution of the common Mo-bearing mineral, molybdenite (MoS 2 ), under both oxic and anoxic conditions. Protochelin promoted molybdenite dissolution under oxic conditions with the formation of MoO 3 on the surface and Mo-siderophore complexes in solution. This was attributed to air-oxidation of both molybdenite and protochelin, as evidenced by lack of dissolution under anoxic conditions but enhanced dissolution by either pre-oxidized protochelin or pre-oxidized molybdenite. Liquid chromatography-mass spectroscopy (LC-MS), X-ray photoelectron spectroscopy (XPS), and Time-of-flight secondary ion mass spectrometry (ToF-SIMS) analyses revealed degradation of protochelin and adsorptions of its byproducts on molybdenite surface to promote dissolution. Conversely, DFOB inhibited molybdenite dissolution under both oxic and anoxic conditions, likely attributed to surface adsorption of DFOB and its weak complexation with Mo(VI) at the circumneutral pH. This work highlights the need to consider the balance between promoting and inhibitory effects of different metallophores on Mo-mineral dissolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Using SF 6 and Xe to Monitor Gas Migration Through Explosion-Generated Fracture Networks

Here we describe a field study where tracer gas was injected into a subsurface cavity created by a small chemical explosion beneath the water table. The main objective of the study is to compare the migration of sulfur hexafluoride (SF 6 ) and xenon (Xe) through an explosion-generated fracture network and to study the influence of ground water on gas transport. A mixture of tracer gases (50% of SF 6 and 50% of Xe) was injected on 31 October 2018 and gas sampling continued until 8 November 2018. We observe similar trends in SF 6 and Xe concentrations at four ground surface sampling sites. The changes in the SF 6 /Xe ratios with time show that more SF 6 than Xe is observed during the barometric pressure lows when the absolute measured concentrations are highest. Conversely, the ratio SF 6 /Xe is less than 1 during the high-pressure intervals when absolute measured concentrations are low. The results of the experiment suggest that during barometric pressure lows the tracer is migrating to the surface primarily by advective gas phase transport, whereas during barometric pressure highs, advection is suppressed and near-surface evaporation of interstitial pore fluid with tracer dissolved in it becomes more important. Thus, the results of the experiment show that the gas concentrations at the surface are controlled by the combined effects of the gas dissolution into pore water and the barometric pressure fluctuations.

58 GEOSCIENCES↗
DOE OSTI2020

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗
DOE OSTI2020

Semiquantitative Estimates of Rainfall Variability During the 8.2 kyr Event in California Using Speleothem Calcium Isotope Ratios

A multiproxy record from a fast-growing stalagmite reveals variable hydroclimate on the California coast across the 8.2 kyr event and a precursor event likely caused by initial drainage of proglacial Lake Agassiz. Using speleothem δ 44 Ca, we develop the first semiquantitative estimates of paleorainfall variability for California through calibration with measurements of the modern climate and cave environment. We find that the magnitude of rainfall variability during the 8.2 kyr event approached the multiyear variability observable in the recent past (1950–2019) and the magnitude of variability during the precursor event likely exceeded this range. Additionally, we observe other instances of multidecadal variability comparable in magnitude to the precursor event during the record. Finally, our work suggests that speleothem calcium isotope ratios are a powerful semiquantitative means to reconstruct paleorainfall, although numerous factors must be assessed in each cave system before applying this approach.

58 GEOSCIENCES↗