Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

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
NASA NTRS1989

Sensitivity of climate and atmospheric CO2 to deep-ocean and shallow-ocean carbonate burial

A model of the carbonate-silicate geochemical cycle is presented that distinguishes carbonate masses produced by shallow-ocean and deep-ocean carbonate burial and shows that reasonable increases in deep-ocean burial could produce substantial warmings over a few hundred million years. The model includes exchanges between crust and mantle; transients from burial shifts are found to be sensitive to the fraction of nondegassed carbonates subducted into the mantle. Without the habitation of the open ocean by plankton such as foraminifera and coccolithophores, today's climate would be substantially colder.

Volk, Tyler↗
NASA NTRS2002

SNC Meteorites and Martian Reservoirs

Jones first suggested that the inverse covariation of initial epsilon (Nd-143) and Sr-87/Sr-86 of the shergottites could be explained by interaction between mantle-derived magmas with another isotopic reservoir(s) (i.e., assimilation or contamination). In that model, magmas were generated in a source region that was isotopically similar to the Nakhla source and the second reservoir(s) was presumed to be crust. The text also permitted the second reservoir to be another type of mantle, but I can confirm that a second mantle reservoir was never seriously considered by that author. Other features of this model were that (i) it occurred at a particular time, 180 m.y. ago, and (ii) the interacting reservoirs had been separated at approximately 4.5 b.y. In a later paper Jones explored this mixing model more quantitatively and concluded that magmas from a Nakhla-like source region at 180 m.y. would fall on or near an isotopic Nd-Sr-Pb hyperplane defined by the shergottites. This criterion was a necessary prerequisite for the parent magma(s) of the shergottites to have initially been Nakhla-like isotopically. At this juncture, it is perhaps worthwhile to note that this mixing model was not presented to explain geochemical variations but as a justification for a 180 m.y. crystallization age for the shergottites and a 1.3 b.y. crystallization age for the nakhlites. In the mid-1980's crystallization ages estimated for Nakhla ranged from approximately 1.3 b.y to 4.5 b.y. Similarly, preferred crystallization ages for the shergottites ranged from 360 m.y., to 1.3 b.y., to 4.5 b.y. In all these models, the 180 m.y. event seen in the shergottites was deemed to be metamorphic. The fit between the Nakhla-like source region and the shergottite hyperplane was a validation both of the 1.3 b.y. igneous age of Nakhla and the 180 m.y. igneous age of the shergottites.

Jones, J. H.↗
DOE OSTI2020

Iodine effective diffusion coefficients through volcanic rock: Influence of iodine speciation and rock geochemistry

Accurate prediction of the subsurface transport of iodine species is important for the assessment of long-term nuclear waste repository performance, as well as monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty, given that radioiodine decays into radioxenon. However, the transport of iodine through intact geologic media is not well understood, compromising our ability to assess risk associated with radioiodine migration. The current study's goal is to quantify the matrix diffusion of iodine species through saturated volcanic rock, with particular attention paid to the redox environment and potential speciation changes. Here, diffusion experiments were run for iodide through lithophysae-rich lava, lithophysae-poor lava, and welded tuff, whereas iodate diffusion was studied through welded tuff. Iodine transport was compared with a conservative tracer, HDO, and effective diffusion coefficients were calculated. Likely due to a combination of size and anion exclusion effects, iodine species diffused more slowly than the conservative tracer through all rock types tested. Furthermore, oxidation of iodide to iodate was observed in the lithophysae-poor lava, affecting transport. Results provide much needed data for subsurface transport models that predict radioiodine migration from underground sources, and indicate the pressing need for geochemical and redox interactions to be incorporated into these models.

58 GEOSCIENCES↗
DOE OSTI2021

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

The Role of KREEP in the Production of Mg-Suite Magmas and Its Influence on the Extent of Mg-Suite Magmatism in the Lunar Crust

The lunar magnesian-suite, or Mg-suite, is a series of ancient plutonic rocks from the lunar crust. They have received a considerable amount of attention from lunar scientists since their discovery for three primary reasons: 1) their ages and geochemistry indicate they represent pristine magmatic samples that crystallized very soon after the formation of the Moon; 2) their ages often overlap with ages of the ferroan anorthosite (FAN) crust; and 3) planetary-scale processes are needed in formation models to account for their unique geochemical features. Taken as a whole, the Mg-suite samples, as magmatic cumulate rocks, approximate a fractional crystallization sequence in the low-pressure forsterite-anorthite-silica system, and thus these samples are generally thought to be derived from layered mafic intrusions which crystallized very slowly from magmas that intruded the anorthositic crust. However, no direct linkages have been established between different Mg-suite samples based either on field relationships or geochemistry.The model for the origin of the Mg-suite, which best fits the limited available data, is one where Mg-suite magmas form from melting of a hybrid cumulate package consisting of deep mantle dunite, crustal anorthosite, and KREEP (potassium-rare earth elements-phosphorus) at the base of the crust under the Procellarum KREEP Terrane (PKT). In this model, these three LMO (Lunar Magma Ocean) cumulate components are brought into close proximity by the cumulate overturn process. Deep mantle dunitic cumulates with an Mg number of approximately 90 rise to the base of the anorthositic crust due to their buoyancy relative to colder, more dense Fe- and Ti-rich cumulates. This hybridized source rock melts to form Mg-suite magmas, saturated in Mg-rich olivine and anorthitic plagioclase, that have a substantial KREEP component.

Elardo, S. M.↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2023

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Rock-Sourced Nitrogen in Semi-Arid, Shale-Derived California Soils

Models suggest that rock-derived nitrogen (N) inputs are of global importance to ecosystem N budgets; however, field studies demonstrating the significance of rock N inputs are rare. We examined rock-derived N fluxes in soils derived from sedimentary rocks along a catena formed under a semi-arid climate. Our measurements demonstrate that there are distinct and traceable pools of N in the soil and bedrock and that the fraction of rock-derived N declines downslope along the catena. We used geochemical mass balance weathering flux measurements to estimate a rock-derived N flux of 0.145 to 0.896 kg ha -1 yr -1 at the ridgecrest. We also developed independent N flux estimates using a 15 N-based isotope mixing model. While geochemical mass-balance-based estimates fell within the 95% confidence range derived from the isotope mixing model (-1.1 to 44.3 kg ha -1 yr -1 ), this range was large due to uncertainty in values for atmospheric 15 N deposition. Along the catena, N isotopes suggest a diminishing effect of rock-derived N downslope. Overall, we found that despite relatively large N pools within the saprolite and bedrock, slow chemical weathering and landscape denudation limit the influence of rock-derived N, letting atmospheric N deposition (7.1 kg ha -1 yr -1 ) and N fixation (0.9-3.1 kg ha -1 yr -1 ) dominate N inputs to this grassland ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Utah FORGE 5-2565: Hydrothermal Evolution of Fracture Properties - Workshop Presentation

This is a presentation on the Evolution of Permeability and Strength Recovery of Shear Fractures Under Hydrothermal Conditions project by the U.S. Geological Survey, presented by Dr. David Lockner. The project's objective was to determine how thermal, hydraulic, mechanical, and chemical processes affect the sustainability of fracture networks in geothermal reservoirs and provide strategies for improved EGS techniques that maximize thermal coupling and increase reservoir longevity. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗
DOE OSTI2021

How much carbon can be added to soil by sorption?

Abstract Quantifying the upper limit of stable soil carbon storage is essential for guiding policies to increase soil carbon storage. One pool of carbon considered particularly stable across climate zones and soil types is formed when dissolved organic carbon sorbs to minerals. We quantified, for the first time, the potential of mineral soils to sorb additional dissolved organic carbon (DOC) for six soil orders. We compiled 402 laboratory sorption experiments to estimate the additional DOC sorption potential, that is the potential of excess DOC sorption in addition to the existing background level already sorbed in each soil sample. We estimated this potential using gridded climate and soil geochemical variables within a machine learning model. We find that mid- and low-latitude soils and subsoils have a greater capacity to store DOC by sorption compared to high-latitude soils and topsoils. The global additional DOC sorption potential for six soil orders is estimated to be 107 $$\pm$$ ± 13 Pg C to 1 m depth. If this potential was realized, it would represent a 7% increase in the existing total carbon stock.

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

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

This closeout presentation summarizes a Department of Energy–funded study that evaluated whether large amounts of hydrogen can be safely stored underground in the North Dakota portion of the Williston Basin. The project combined lab testing, computer simulations, and basin‑wide analysis to assess saline formations, depleted oil and gas reservoirs, and salt formations for hydrogen storage capacity, recovery, and risks. Results show that underground hydrogen storage is technically feasible across multiple formation types, with depleted oil and gas reservoirs offering higher recovery and saline formations providing large long‑term storage potential. The study also identifies key challenges—such as wellbore material durability and gas purity management—and recommends pilot projects and further site‑specific studies to support future commercialization.

08 HYDROGEN↗
NASA NTRS2015

Lead Isotopes in Olivine-Phyric Shergottite Tissint: Implications for the Geochemical Evolution of the Shergottite Source Mantle

Geochemically-depleted shergottites are basaltic rocks derived from a martian mantle source reservoir. Geochemical evolution of the martian mantle has been investigated mainly based on the Rb-Sr, Sm-Nd, and Lu-Hf isotope systematics of the shergottites [1]. Although potentially informative, U-Th- Pb isotope systematics have been limited because of difficulties in interpreting the analyses of depleted meteorite samples that are more susceptible to the effects of near-surface processes and terrestrial contamination. This study conducts a 5-step sequential acid leaching experiment of the first witnessed fall of the geochemically-depleted olivinephyric shergottite Tissint to minimize the effect of low temperature distrubence. Trace element analyses of the Tissint acid residue (mostly pyroxene) indicate that Pb isotope compositions of the residue do not contain either a martian surface or terrestrial component, but represent the Tissint magma source [2]. The residue has relatively unradiogenic initial Pb isotopic compositions (e.g., 206Pb/204Pb = 10.8136) that fall within the Pb isotope space of other geochemically-depleted shergottites. An initial μ-value (238U/204Pb = 1.5) of Tissint at the time of crystallization (472 Ma [3]) is similar to a time-integrated mu- value (1.72 at 472 Ma) of the Tissint source mantle calculated based on the two-stage mantle evolution model [1]. On the other hand, the other geochemically-depleted shergottites (e.g., QUE 94201 [4]) have initial μ-values of their parental magmas distinctly lower than those of their modeled source mantle. These results suggest that only Tissint potentially reflects the geochemical signature of the shergottite mantle source that originated from cumulates of the martian magma ocean

Moriwaki, R.↗
DOE OSTI2021

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

This data package contains input files for TOUGHREACT for a modeling study examining the effects of managed aquifer recharge on agricultural lands on nitrate cycling and transport in the Central Valley of CA near Modesto. The files contain all the geochemical species, reactions, and hydrological parameters for the model. The files are text files used for the TOUGH family of code created by LBNL. To use the files a license is required. https://tough.lbl.gov/licensing-download/toughreact-licensing-download/Accompanying Paper Abstract: Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO3-) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying water all-at-once rather than in increments for a fixed volume of recharge transports higher concentrations of NO3- deeper into the profile, which has the potential to exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as riparian corridors, floodplains, wetlands, and other managed landscapes like dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Community Data Mining Approach for Surface Complexation Database Development

This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. Here the framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Understanding Contaminant Migration within a Dynamic River Corridor through Field Experiments and Reactive Transport Modeling

The behavior of a persistent uranium plume within an extended river corridor at the DOE Hanford site is dominantly controlled by river stage fluctuations in the adjacent Columbia River. The plume behavior is further complicated by substantial heterogeneity in physical and geochemical properties of the host aquifer sediments. Multi-scale field and laboratory experiments and reactive transport modeling were integrated to understand the complex plume behavior influenced by highly variable hydrologic and geochemical conditions in time and space. In this paper, we (1) describe multiple data sets from field-scale uranium adsorption and desorption experiments performed at our experimental well-field, (2) develop a reactive transport model that incorporates hydrologic and geochemical heterogeneities characterized from multi-scale and multi-type datasets and a surface complexation reaction network based on laboratory studies, and (3) compare the modeling and observation results to provide insights on how to refine the conceptual model and reduce prediction uncertainties. The experimental results revealed significant spatial variability in uranium adsorption/desorption behavior, while modeling demonstrated that ambient hydrologic and geochemical conditions and heterogeneities in sediment physical and chemical properties both contributed to complex plume behavior and its persistence. This research underscores the great challenges in adequately characterizing this type of site to model the reactive transport processes over scales of 10 m or more. Our analysis provides important insights into the characterization, understanding, modeling, and remediation of groundwater contaminant plumes influenced by dynamic surface water and groundwater interactions.

Chen, Xingyuan↗
NASA NTRS1988

Martian regolith geochemistry and sampling techniques

Laboratory study of samples of the intermediate and fine-grained regolith, including duricrust peds, is a fundamental prerequisite for understanding the types of physical and chemical weathering processes on Mars. The extraordinary importance of such samples is their relevance to understanding past changes in climate, availability (and possible physical state) of water, eolian forces, the thermal and chemical influences of volcanic and impact processes, and the inventory and fates of Martian volatiles. Fortunately, this regolith material appears to be ubiquitous over the Martian surface, and should be available at many different landing sites. Viking data has been interpreted to indicate a smectite-rich regolith material, implying extensive weathering involving aqueous activity and geochemical alteration. An all-igneous source of the Martian fines has also been proposed. The X-ray fluorescence measurement data set can now be fully explained in terms of a simple two-component model. The first component is silicate, having strong geochemical similarities with Shergottites, but not other SNC meteorites. The second component is salt. Variations in these components could produce silicate and salt-rich beds, the latter being of high potential importance for microenvironments in which liquid water (brines) could exist. It therefore would be desirable to scan the surface of the regolith for such prospects.

Clark, B. C.↗
DOE OSTI2022

Early diagenetic processes in an iron-dominated marine depositional system

The early diagenetic interplay between reactive iron, sulfur, and organic matter in the bathymetrically isolated Santa Monica Basin (SMB) sediments are investigated in this study. Here, we explore solid-phase and porewater profiles from the basin, supplemented with a transect from 71 to 907 m water depth that includes oxygenated (>60 μM O 2 ) bottom waters near the coast and oxygen-deficient waters (~4 µM O 2 ) in the basin. The geochemical data of the basin sediments are further scrutinized by means of reactive transport modeling. The results show that the basin sediments do not follow the traditional geochemical signatures of oxygen-deficient settings. A lack of dissolved sulfide accumulation and sulfurized iron persists despite the sediments being deposited under reducing conditions (without bioturbation/bioirrigation), strong organic carbon input (TOC up to 5.0 wt%), and active dissimilatory sulfate reduction. Not only did we find an exceptional enrichment in highly reactive Fe in the surface sediments (~45 % of total Fe), but the enrichment of reactive Fe, including ferrihydrite, persists downcore and coexists with high levels of dissolved Fe. The enhanced preservation of Fe oxides and lack of iron-sulfide precipitation is in part explained by detection via Mössbauer spectra of iron oxides bounded to organic matter (Fe[III]-OM coprecipitates). The modeled Fe budget shows that most of the Fe oxides in the surface sediments are internally recycled by upward diffusion and subsequent oxidation of Fe 2+ . Sulfide oxidation coupled to Fe reduction effectively precludes sulfide accumulation while enhancing build-up of dissolved Fe, fueling the Fe cycle within the first 5 cm depth. Continuous reoxidation of Fe 2+ enhances the formation of Fe(III)-OM coprecipitates, limiting the amount of reactive organic matter. In the unavailability of labile organic matter, other than within the uppermost layers, the organic-rich sediment profiles are dominated by Fe cycling that limits the production and preservation of sulfides and enhances the preservation of Fe oxides and organic carbon. Finally, this study highlights key local controls on Fe availability in marginal basins and describes an intricate biogeochemical C-Fe-S cycling in modern and possibly ancient marine systems with important implications for Fe availability in the marine realm.

58 GEOSCIENCES↗