Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geochemical”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Decadal-Scale Hydrothermal Alteration of Basalt: Implications for Long-Term Reactivity in Subsurface Technologies

Global demands to secure critical mineral supply chains and managing water resources require innovative resource management strategies. Subsurface basalt reservoirs offer a dual-purpose application, where CO2 enhanced critical metal recovery can be coupled with water remediation strategies for the extraction of economically important and environmentally relevant elements. This study utilizes static batch experiments conducted on Columbia River Basalt Group samples to evaluate long-term fluid–rock interactions relevant to coupled subsurface applications. Samples were reacted under elevated temperature and pressure conditions, for 296 and 4,061 days, providing a unique dataset to assess geochemical evolution, element partitioning, and aqueous chemistry over decadal timescales. The alteration was dominated by silicate, aluminosilicate, and oxide phases. Progressive hydrothermal alteration increased the cation exchange capacity through the formation of secondary minerals (e.g. clays and zeolites), influencing solute retention, ion exchange behavior, and aqueous chemistry. Grain-scale and bulk analyses revealed selective cation release and retention patterns, establishing a mechanistic baseline for predicting element mobility. These geochemical constraints inform the prospective use of CO2-enhanced metal recovery for critical resources (e.g. Co, Ni, Mn, Mg, Ti, Fe) and associated water-quality impacts, providing a foundation for the design and management of engineered fluid–rock systems in basalt reservoirs.

Gibson, Marie N.↗
DOE OSTI2019

A Closer Look at Fe(II) Passivation of Goethite

Our understanding of how Fe(II) reacts with Fe(III) oxides has evolved significantly based on evidence for electron transfer at the oxide-water interface and Fe(II)-catalyzed recrystallization. Recent experiments, as well as computational calculations, further suggest that defects at oxide surfaces can facilitate electron transfer. There is, however, some evidence that these processes cease after continued contact with Fe(II). Here, we explore the sustainability of these processes over longer time frames by measuring whether exposure to Fe(II) inhibits Fe(II)-goethite electron transfer, and whether this inhibition is reversible. To quantify the extent of electron transfer, we used selective isotope labeling with 57Fe Mössbauer spectroscopy. The results revealed that pre-reacting goethite with Fe(II) decreases the extent of electron transfer during subsequent exposures, and that less electron transfer occurs with increasing reaction time. Similarly, X-ray absorption spectroscopy and magnetic circular dichroism suggest that reaction with Fe(II) consumes surface defects and forms a passivation layer, making goethite less susceptible to electron transfer. When Fe(II) is removed by either wet chemical extraction or oxidation, electron transfer is mostly restored, suggesting that inhibition of electron transfer is reversible. Our results show that in environmental systems, which are commonly subjected to transient geochemical fluctuations, electron transfer between Fe(II) and Fe oxides, and processes linked to it, are likely to be relevant beyond just short time scales.

Andrade de Notini, Luiza↗
DOE OSTI2019

Collapse of Reacted Fracture Surface Decreases Permeability and Frictional Strength

Abstract Geochemical and geomechanical perturbations of the subsurface caused by the injection of fluids present risks of leakage and seismicity. This study investigated how acidic fluid flow affects hydraulic and frictional properties of fractures using experiments with 3.8‐cm‐long specimens of Eagle Ford shale, a laminated shale with carbonate‐rich strata. In low‐pressure flow cells, one set of samples was exposed to acidic brine and another set was exposed to neutral brine. X‐ray computed tomography and energy‐dispersive X‐ray spectroscopy revealed that samples exposed to acidic brine were calcite‐depleted and had developed a porous altered layer, while the other set showed no evidence of alteration. After reaction, samples were compressed and sheared in a triaxial cell that supplied normal stress and differential pore pressure at prescribed sliding velocities, independently measuring friction and permeability. During the initial compression, the porous altered layer collapsed into fine particles that filled the fracture. This effectively impeded flow and sealed the fracture, resulting in fracture permeability to decrease 1 to 2 orders of magnitude relative to the unaltered fractures. This is a favorable outcome in subsurface applications where the goal is to reduce leakage risks. However, during shear the reacted fracture had lower frictional strength because the fine‐grained particles in the collapsed layer prevented the formation of interlocking microasperities. Therefore, coupled geochemical and geomechanical processes that could favorably seal fractures could also increase the likelihood of induced seismicity. These findings have important implications for geological carbon sequestration, pressurized fluid energy storage, geothermal energy, and other subsurface technologies.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic Visibility of Eclogite in the Earth’s Upper Mantle—Implications From High Pressure‐Temperature Single‐Crystal Elastic Properties of Omphacite

Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite: K S0 ’ = 4.5(1), G 0 ’ = 1.53(5), ∂K S0 /∂ T = −0.029(5) GPa/K, ∂G 0 /∂ T = −0.013(5) GPa/K, with K S0 = 123(3) GPa, G 0 = 74(2) GPa, and ρ 0 = 3.34(1) g/cm 3 . We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle.

Hao, Ming↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗
DOE OSTI2024

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Linking the Spin Transition of Ferric Iron in δ‐(Al,Fe)OOH to Water Storage in the Lower Mantle

As the most massive geochemical reservoir, the lower mantle affects the Earth's budget of volatile elements, including hydrogen or H 2 O. The properties of minerals in the lower mantle are further affected by changes in the electronic configurations of iron cations, that is, by spin transitions. The feedback between spin transitions and potential storage of H 2 O in solid hydrous phases in the lower mantle, however, remains unexplored. By combining high‐pressure nuclear resonant inelastic X‐ray scattering and high‐pressure high‐temperature X‐ray diffraction experiments, we constrained the thermal equation of state of δ‐(Al,Fe)OOH, a member of the phase H solid solution. Based on the derived thermal equation of state of δ‐(Al,Fe)OOH and the underlying thermodynamic model, we calculate the excess Gibbs free energy that arises from the spin transition of ferric iron in this compound and evaluate the effect on phase equilibria. The results of our analysis show that the spin transition of ferric iron in phase H may significantly reduce the thermodynamic activity and hence the concentration of H 2 O in a coexisting hydrous melt. As a consequence, nominally anhydrous minerals of the lower mantle may become dehydrated in the presence of phase H. Our analysis further suggests that, under certain conditions, the spin transition may expand the thermal stability of Fe 3+ ‐bearing phase H and create a geochemical link between the storage of H 2 O in phase H and ferric iron in the lower mantle.

Buchen, Johannes [Univ. of Bayreuth (Germany); Cal↗
DOE OSTI2026

Soil Organic Matter Stabilization by Polyvalent Cations in a Buried Alkaline Soil

Buried paleosols can store large quantities of organic carbon (C), much of which persists for millennia due to isolation from surface processes that promote decomposition. Subsoil organic matter (SOM) persistence is often enhanced by mineral associations and ionic conditions—particularly high clay content and polyvalent cations—that limit microbial degradation and leaching. However, the vulnerability of these deep C stocks under erosion or environmental change remains poorly understood. This study investigates controls on SOM stabilization in the Brady paleosol and overlying modern soils across contrasting geomorphic settings in the Great Plains of Nebraska where Late Quaternary loess deposition and erosion created a sequence of buried and exposed paleosols. We sampled soils along burial and erosional toposequences and analyzed their physicochemical properties and radiocarbon‐based persistence of occluded particulate organic matter (oPOM) and mineral fractions (MF). Brady paleosol showed greater persistence (lower Fm) of oPOM and MF than modern soils, particularly under burial. This was linked to higher silt and clay content, elevated electrical conductivity, and increased exchangeable calcium and magnesium content, supporting roles for organo‐mineral interactions, flocculation, and carbonate cementation. In modern soils, SOM persistence and C content were more strongly tied to pH and cation exchange capacity. Erosional exposure reduced SOM stability and promoted geochemical convergence toward modern surface soils. These findings show that burial enhances SOM persistence via multiple stabilization mechanisms, while erosion increases subsoil C vulnerability. Our results underscore the importance of geomorphic and geochemical context in predicting soil C stability under environmental change.

Environmental sciences↗
DOE OSTI2026

Carbon monoxide oxidation expands the known metabolic capacity in anaerobic methanotrophic consortia

Consortia of anaerobic methane-oxidizing archaea (ANME-2) and sulphate-reducing bacteria (SRB) represent globally relevant syntrophic associations capable of growing with minimal amounts of free energy and can persist when methane becomes limiting. Carbon monoxide (CO) has been reported in seep environments and represents a thermodynamically favourable alternative electron donor due to its low reduction potential. Here, we show that environmental ANME-SRB consortia can oxidize CO in the absence of methane, in anoxic microcosm experiments using a combination of stable isotope geochemical tracers, metatranscriptomics, and single cell activity measurements (FISH–nanoSIMS). The oxidation of CO was coupled with sulphate-reduction by syntrophic consortia, and, in the absence of sulphate, through CO 2 reduction to methane by ANME-2. Under these conditions, the production of methane was one ninth the rate of methanotrophy coupled to sulphate-reduction. Paired single cell FISH-nanoSIMS analysis of anabolic activity indicates that CO respiration appears to support cell maintenance rather than active growth, consistent with the observed down-regulation of energy generating complexes in ANME (e.g., mtr, rnf, etc.). The versatile capability of CO oxidation by anaerobic methanotrophic consortia broadens our understanding of carbon cycling in methane seeps and highlights potential mechanisms of resilience by methanotrophic archaea under changing geochemical regimes.

03 NATURAL GAS↗
DOE OSTI2022

Assessment of chemo-mechanical impacts of CO 2 sequestration on the caprock formation in Farnsworth oil field, Texas

This study evaluates the chemo-mechanical influence of injected CO 2 on the Morrow B sandstone reservoir and the upper Morrow shale caprock utilizing data from the inverted 5-spot pattern centered on Well 13-10A within the Farnsworth unit (FWU). This study also seeks to evaluate the integrity of the caprock and the long-term CO 2 storage capability of the FWU. The inverted 5-spot pattern was extracted from the field-scale model and tuned with the available field observed data before the modeling work. Two coupled numerical simulation models were utilized to continue the study. First, a coupled hydro-geochemical model was constructed to simulate the dissolution and precipitation of formation minerals by modeling three intra-aqueous and six mineral reactions. In addition, a coupled hydro-geomechanical model was constructed and employed to study the effects of stress changes on the caprock’s porosity, permeability, and ground displacement. The Mohr–Coulomb circle and failure envelope were used to determine caprock failure. In this work, the CO 2 -WAG injection is followed by the historical field-observed strategy. During the forecasting period, a Water Alternating Gas (WAG) injection ratio of 1:3 was utilized with a baseline bottom-hole pressure constraint of 5500 psi for 20 years. A post-injection period of 1000 years was simulated to monitor the CO 2 plume and its effects on the CO 2 storage reservoir and caprock integrity. The simulation results indicated that the impacts of the geochemical reactions on the porosity of the caprock were insignificant as it experienced a decrease of about 0.0003% at the end of the 1000-year post-injection monitoring. On the other hand, the maximum stress-induced porosity change was about a 1.4% increase, resulting in about 4% in permeability change. It was estimated that about 3.3% of the sequestered CO 2 in the formation interacted with the caprock. Despite these petrophysical property alterations and CO 2 interactions in the caprock, the caprock still maintained its elastic properties and was determined to be far from its failure.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

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

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
DOE OSTI2021

Soil metabolome response to whole-ecosystem warming at the Spruce and Peatland Responses under Changing Environments experiment

In this study, a suite of complementary environmental geochemical analyses, including NMR and gas chromatography-mass spectrometry (GC-MS) analyses of central metabolites, Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) of secondary metabolites, and lipidomics, was used to investigate the influence of organic matter (OM) quality on the heterotrophic microbial mechanisms controlling peatland CO2, CH4, and CO2:CH4 porewater production ratios in response to climate warming. Our investigations leverage the Spruce and Peatland Responses under Changing Environments (SPRUCE) experiment, where air and peat warming were combined in a whole-ecosystem warming treatment. We hypothesized that warming would enhance the production of plant-derived metabolites, resulting in increased labile OM inputs to the surface peat, thereby enhancing microbial activity and greenhouse gas production. Because shallow peat is most susceptible to enhanced warming, increases in labile OM inputs to the surface, in particular, are likely to result in significant changes to CO2 and CH4 dynamics and methanogenic pathways. In support of this hypothesis, significant correlations were observed between metabolites and temperature consistent with increased availability of labile substrates, which may stimulate more rapid turnover of microbial proteins. An increase in the abundance of methanogenic genes in response to the increase in the abundance of labile substrates was accompanied by a shift toward acetoclastic and methylotrophic methanogenesis. Our results suggest that as peatland vegetation trends toward increasing vascular plant cover with warming, we can expect a concomitant shift toward increasingly methanogenic conditions and amplified climate–peatland feedbacks.

Wilson, Rachel M.↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2023

Tectonic and geological setting influence hot spring microbiology

Hydrothermal systems form at divergent and convergent boundaries of lithospheric plates and within plates due to weakened crust and mantle plumes, playing host to diverse microbial ecosystems. Little is known of how differences in tectonic setting influence the geochemical and microbial compositions of these hydrothermal ecosystems. Here, for this work, coordinated geochemical and microbial community analyses were conducted on 87 high-temperature (>65°C) water and sediment samples from hot springs in Yellowstone National Park, Wyoming, USA ( n = 41; mantle plume setting), Iceland ( n = 41, divergent boundary), and Japan ( n = 5; convergent boundary). Region-specific variation in geochemistry and sediment-associated 16S rRNA gene amplicon sequence variant (ASV) composition was observed, with 16S rRNA gene assemblages being nearly completely distinguished by region and pH being the most explanatory parameter within regions. Several low abundance ASVs exhibited cosmopolitan distributions across regions, while most high-abundance ASVs were only identified in specific regions. The presence of some cosmopolitan ASVs across regions argues against dispersal limitation primarily shaping the distribution of taxa among regions. Rather, the results point to local tectonic and geologic characteristics shaping the geochemistry of continental hydrothermal systems that then select for distinct microbial assemblages. These results provide new insights into the co-evolution of hydrothermal systems and their microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sulfur cycling at natural hydrocarbon and sulfur seeps in Santa Paula Creek, CA

Biogeochemical cycling of sulfur is relatively understudied in terrestrial environments compared to marine environments. However, the comparative ease of access, observation, and sampling of terrestrial settings can expand our understanding of organisms and processes important in the modern sulfur cycle. Furthermore, these sites may allow for the discovery of useful process analogs for ancient sulfur-metabolizing microbial communities at times in Earth's past when atmospheric O 2 concentrations were lower and sulfide was more prevalent in Earth surface environments. Here, we identified a new site at Santa Paula Creek (SPC) in Ventura County, CA—a remarkable freshwater, gravel-bedded mountain stream charged with a range of oxidized and reduced sulfur species and heavy hydrocarbons from the emergence of subsurface fluids within the underlying sulfur- and organic-rich Miocene-age Monterey Formation. SPC hosts a suite of morphologically distinct microbial biofacies that form in association with the naturally occurring hydrocarbon seeps and sulfur springs. We characterized the geology, stream geochemistry, and microbial facies and diversity of the Santa Paula Creek ecosystem. Using geochemical analyses and 16S rRNA gene sequencing, we found that SPC supports a dynamic sulfur cycle that is largely driven by sulfide-oxidizing microbial taxa, with contributions from smaller populations of sulfate-reducing and sulfur-disproportionating taxa. This preliminary characterization of SPC revealed an intriguing site in which to study geological and geochemical controls on microbial community composition and to expand our understanding of sulfur cycling in terrestrial environments.

58 GEOSCIENCES↗