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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Enhanced Groundwater Flow on and Below Vera Rubin Ridge, the Murray Formation, Gale Crater: Evidence from Thermochemical Modeling.

NASA’s Mars Science Laboratory Curiosity rover has been exploring Vera Rubin ridge (VRR), part of the Murray formation in Gale crater, Mars, between sol 1809 and 2302. Evidence for Fe-oxides and phyllosilicates in mineralogical and geochemical data for this region was returned by Curiosity [1-5]. We applied thermochemical modeling to con-strain the formation conditions of the phyllosilicate-hematite assemblage identified on and below VRR. Average alteration compositions for the Murray formation on and below VRR were derived using CheMin and APXS data. These compositions were reacted with Gale Portage Water (GPW) between 25–100 °C and for 10% and 50% Fe3+/Fetot of the host rock [6]. Here we summarize models run at 50 °C and 10% Fe3+/Fetot for alteration compositions derived from Murray host rock compositions.

Turner, S. M. R.↗

NASA NTRS1991

Sources and geochemical evolution of cyanide and formaldehyde

The major source of cyanide has, in current paleoatmospheric models, been assumed to be the reaction of photodissociated thermospheric nitrogen with a limiting supply of stratospheric methane. Formaldehyde may be produced with more ease from an atmosphere of carbon dioxide as the dominant carbon species, and from carbonate in solution or sorbed in double layer hydroxide minerals. Potentially more important sources for cyanide and other carbon containing molecules are the partially photoprotected northern and southern auroral ovals where continuous currents reaching several mega-amperes induce ion-molecule reactions, extending into the lower stratosphere. In simulated environments of this kind, the cyanide ion is known to be produced from oxidized carbon species potentially more abundant than methane. Rainout of cyanide and formaldehyde place them in two different geochemical reaction reservoirs. In the anoxic Archean hydrosphere, about 1mM in Fe2(+), the cyanide ion would have been efficiently converted to the stable ferrocyanide complex Fe(CN) sub 6(4-), protecting it from the commonly considered fate of decomposition by hydrolysis, and eventually incorporating it in pyroaurite type minerals, most efficiently in green rust where it converts to insoluble ferriferrocyanide, prussian blue.

Arrhenius, G.↗

DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Timing and Geochemical Evolution of Triassic and Jurassic Magmatism during the Break-up of Pangea along the Eastern North American Margin

While our understanding of continental break-up has seen dramatic advances over the past decades, the impact of syn-rift melt generation on the transition from mechanical stretching to seafloor spreading remains enigmatic. Early Mesozoic continental break-up along the U.S. Eastern North American Margin (ENAM) is characterized by protracted Triassic lithospheric stretching and syn-rift deposition interrupted by voluminous Central Atlantic Magmatic Province (CAMP) magmatism. This magmatic phase is followed by ~50 Myr of protracted rifting prior to steady-state seafloor spreading. Mesozoic magmatism recognized in outcrop and aeromagnetic datasets along the ENAM has been summarily attributed to CAMP, although extensive regional dikes preserved in the ENAM proximal rift domain remain largely undated. This study provides new constraints on the timing, duration, and geochemical evolution of syn-rift ENAM magmatism using in-situ apatite U-Pb geochronology and whole-rock geochemistry and integrates these data into a tectonic model. We present in-situ apatite U-Pb ages and whole-rock geochemical data for 16 diabase dikes and sills from NJ, PA, VA, NC, and SC. A Palisades Sill sample yielded an age of 201 ± 8 Ma and agrees with previous 40 Ar/39 Ar and U-Pb ages. Our new in-situ apatite U-Pb ages fall into three clusters and suggest three main pulses of magmatism along the ENAM, with an initial pulse of magmatism at ~205-195 Ma attributed to CAMP. Whole-rock geochemistry of this pulse shows silica (46-53 wt.%), total alkalis (2-4 wt.%), and enrichment of LREEs relative to HREEs, comparable to previously reported geochemical data for CAMP-related rocks along the ENAM. A second magmatic pulse occurred between ~181-175 Ma and shows slightly higher silica and total alkalis, compared to the older CAMP-age phase. The third pulse between ~166-145 Ma generally shows slightly higher silica and total alkalis. LOI and Sr/Yb values of the two younger pulses are similar to CAMP-age samples suggesting no post-crystallization alteration. These new data demonstrate long-lived syn-rift off-axis magmatism persisted ~50 Myr after the onset of CAMP. This post-CAMP off-axis magmatism is likely related to slow magma-limited spreading during the protracted transition to symmetric seafloor spreading along the ENAM.

Z. Foster-Baril↗

DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗

DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗

DOE OSTI2020

Evaluation of Timber Mountain Recharge and Groundwater Flow in Relation to Pahute Mesa-Oasis Valley Flow System

The purpose of this study was to evaluate the effect of Timber Mountain recharge on groundwater flow in the surrounding Pahute Mesa-Oasis Valley groundwater basin using several different methodologies. These methodologies included examining climatic and soil moisture data, estimating groundwater recharge by chloride mass balance, interpreting water chemical and isotopic data, and modeling groundwater flow along faults.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

Accretion process of the moon

Recent geochemical and geophysical data suggest that the initial temperature of the moon was strongly peaked toward the lunar surface. To explain such an initial temperature distribution, a simple model of accretion process of the moon is presented. The model assumes that the moon was formed from the accumulation of the solid particles or gases in the isolated, closed cloud. Two equations are derived to calculate the accretion rate and surface temperature of the accreting moon. Numerical calculations are made for a wide range of the parameters particle concentration and particle velocity in the cloud. A limited set of the parameters gives the initial temperature profiles as required by geochemical and geophysical data. These models of the proto-moon cloud indicate that the lunar outermost shell, about 400 km thick, was partially or completely molten just after the accretion of the moon and that the moon should have been formed in a period shorter than 1000 years. If the moon formed at a position nearer to the earth than its present one, the moon might have been formed in a period of less than one year.

Mizutani, H.↗

NASA NTRS2004

Gas Hydrate Stability at Low Temperatures and High Pressures with Applications to Mars and Europa

Gas hydrates are implicated in the geochemical evolution of both Mars and Europa [1- 3]. Most models developed for gas hydrate chemistry are based on the statistical thermodynamic model of van der Waals and Platteeuw [4] with subsequent modifications [5-8]. None of these models are, however, state-of-the-art with respect to gas hydrate/electrolyte interactions, which is particularly important for planetary applications where solution chemistry may be very different from terrestrial seawater. The objectives of this work were to add gas (carbon dioxide and methane) hydrate chemistries into an electrolyte model parameterized for low temperatures and high pressures (the FREZCHEM model) and use the model to examine controls on gas hydrate chemistries for Mars and Europa.

Marion, G. M.↗

NASA NTRS2019

Geochemistry of Diagenetic Fluids and Sedimentary Protolith of Gale Crater

The geochemical composition of the sediments of Gale crater can be modeled using two component mixing between a weathered basalt and a sulfate rich endmember. This suggests that the amount of in-situ weathering occurring after deposition was minimal, and most of the variation in weathering indi-ces can be explained simply through addition of Ca-sulfate minerals during diagenesis. Evidence for burial of Gale sediments can be explained through the pres-ence of large ice/dust deposits.

Niles, P. B.↗

NASA NTRS2024

Two mechanisms for the spontaneous emergence, execution, and reprogramming of chemical logic circuits

A key question in origin-of-life research is how primitive life or life-like systems obtained reusable and reprogrammable machinery of processing information to respond to various environmental and internal stimuli. One strategy to approach this question is to view living systems as universal chemical computers. This approach focuses upon uncovering how reusable and reprogrammable chemical NOT, AND, and OR logic gates could appear under prebiotically possible scenarios where simple chemicals and second- or lower-order reactions are more frequent than macromolecules and higher-order reactions. In this article, we theoretically model two generic mechanisms of forming Boolean gates based on geochemically plausible reaction systems that could have preceded complex biopolymers. The mechanisms exploit simple bistable autocatalytic systems to encode binary variables and use heterogeneous dilution rates to differentially activate or deactivate autocatalysis in different reactors. We found that NOT, AND, and OR gates, as well as more complex logic circuits, can spontaneously arise from connected, well-mixed flow reactors based on (i) catalysis by a heterodimer between mutually inhibitory autocatalysts or (ii) selectively permeable tunnels connecting reactors. Our results suggest that simple machinery can support the prebiotic emergence of complex chemical logic circuits and provide new guides for experimentally actualizing chemical computation.

Zhen Peng↗

DOE OSTI2023

Environmental controls on observed spatial variability of soil pore water geochemistry in small headwater catchments underlain with permafrost

Abstract. Soil pore water (SPW) chemistry can vary substantially across multiple scales in Arctic permafrost landscapes. The magnitude of these variations and their relationship to scale are critical considerations for understanding current controls on geochemical cycling and for predicting future changes. These aspects are especially important for Arctic change modeling where accurate representation of sub-grid variability may be necessary to predict watershed-scale behaviors. Our research goal is to characterize intra- and inter-watershed soil water geochemical variations at two contrasting locations in the Seward Peninsula of Alaska, USA. We then attempt to identify the key factors controlling concentrations of important pore water solutes in these systems. The SPW geochemistry of 18 locations spanning two small Arctic catchments was examined for spatial variability and its dominant environmental controls. The primary environmental controls considered were vegetation, soil moisture and/or redox condition, water–soil interactions and hydrologic transport, and mineral solubility. The sampling locations varied in terms of vegetation type and canopy height, presence or absence of near-surface permafrost, soil moisture, and hillslope position. Vegetation was found to have a significant impact on SPW NO3- concentrations, associated with the localized presence of nitrogen-fixing alders and mineralization and nitrification of leaf litter from tall willow shrubs. The elevated NO3- concentrations were, however, frequently equipoised by increased microbial denitrification in regions with sufficient moisture to support it. Vegetation also had an observable impact on soil-moisture-sensitive constituents, but the effect was less significant. The redox conditions in both catchments were generally limited by Fe reduction, seemingly well-buffered by a cache of amorphous Fe hydroxides, with the most reducing conditions found at sampling locations with the highest soil moisture content. Non-redox-sensitive cations were affected by a wide variety of water–soil interactions that affect mineral solubility and transport. Identification of the dominant controls on current SPW hydrogeochemistry allows for qualitative prediction of future geochemical trends in small Arctic catchments that are likely to experience warming and permafrost thaw. As source areas for geochemical fluxes to the broader Arctic hydrologic system, geochemical processes occurring in these environments are particularly important to understand and predict with regards to such environmental changes.

Conroy, Nathan Alec (ORCID:0000000305973373)↗

NASA NTRS2016

Testing and Resilience of the Impact Origin of the Moon

The leading hypothesis for the origin of the Moon is the giant impact model, which grew out of the post-Apollo science community. The hypothesis was able to explain the high E-M system angular momentum, the small lunar core, and consistent with the idea that the early Moon melted substantially. The standard hypothesis requires that the Moon be made entirely from the impactor, strangely at odds with the nearly identical oxygen isotopic composition of the Earth and Moon, compositions that might be expected to be different if Moon came from a distinct impactor. Subsequent geochemical research has highlighted the similarity of both geochemical and isotopic composition of the Earth and Moon, and measured small but significant amounts of volatiles in lunar glassy materials, both of which are seemingly at odds with the standard giant impact model. Here we focus on key geochemical measurements and spacecraft observations that have prompted a healthy re-evaluation of the giant impact model, provide an overview of physical models that are either newly proposed or slightly revised from previous ideas, to explain the new datasets.

Righter, K.↗

NASA NTRS2014

Preliminary Report on U-Th-Pb Isotope Systematics of the Olivine-Phyric Shergottite Tissint

Geochemical studies of shergottites suggest that their parental magmas reflect mixtures between at least two distinct geochemical source reservoirs, producing correlations between radiogenic isotope compositions, and trace element abundances.. These correlations have been interpreted as indicating the presence of a reduced, incompatible-element- depleted reservoir and an oxidized, incompatible-element-rich reservoir. The former is clearly a depleted mantle source, but there has been a long debate regarding the origin of the enriched reservoir. Two contrasting models have been proposed regarding the location and mixing process of the two geochemical source reservoirs: (1) assimilation of oxidized crust by mantle derived, reduced magmas, or (2) mixing of two distinct mantle reservoirs during melting. The former clearly requires the ancient martian crust to be the enriched source (crustal assimilation), whereas the latter requires a long-lived enriched mantle domain that probably originated from residual melts formed during solidification of a magma ocean (heterogeneous mantle model). This study conducts Pb isotope and U-Th-Pb concentration analyses of the olivine-phyric shergottite Tissint because U-Th-Pb isotope systematics have been intensively used as a powerful radiogenic tracer to characterize old crust/sediment components in mantle- derived, terrestrial oceanic island basalts. The U-Th-Pb analyses are applied to sequential acid leaching fractions obtained from Tissint whole-rock powder in order to search for Pb isotopic source components in Tissint magma. Here we report preliminary results of the U-Th-Pb analyses of acid leachates and a residue, and propose the possibility that Tissint would have experienced minor assimilation of old martian crust.

Moriwaki, R.↗

NASA NTRS1995

Effect of the nonlinearity of the carbonate system on partial pressure of carbon dioxide in the oceans

Partial pressure of CO2 is a nonlinear function of several seawater properties. Due to the nonlinearity in this relationship, the partial pressure of a uniform ocean would be different from that of a nonuniform ocean with the same bulk seawater properties. Assuming uniformity of seawater properties at some temporal and spatial scales in carbon models leads to systematic errors in partial pressure of CO2. In this paper we evaluate the magnitude of these errors. We partition the Geochemical Ocean Sections Study and Transient Tracers in the Oceans data according to the horizontal structure of several box models from the literature. Our results suggest that assumption of uniformity at large scales leads to understanding of underestimation of global surface ocean partial pressure of CO2 by at least 3 - 12 microatm. Nonlinear effects also introduce systematic errors in the buffer factor estimated from bulk seawater properties. We find the standard deviation of partial pressure of CO2 to be an indicator of the magnitude of the nonlinear effects. We discuss the implications of these errors for some conclusions drawn from carbon models. Biogeochemical processes, such as mixing, gas exchange, or biological activity, influence the distribution of the seawater properties. A shift in spatial or temporal patterns of these processes can modify the nonuniformity of the seawater properties and thus alter the partial pressure of the surface waters, even if the mean intensities of the processes remain constant.

Trela, Piotr↗

NASA NTRS1993

Melt rock components in KREEPy breccia 15205: Petrography and mineral chemistry of KREEP basalts and quartz-normative mare basalts

Many current models for the origin of lunar highland rocks feature as an essential component the assimilation of KREEPy material by primitive magmas parental to the Mg-rich suite and alkali suite plutonic rocks. Similar models have also been proposed for the origin of various mare basalt suites. However, any model which considers assimilation of KREEP an important petrologic process must sooner-or-later deal with the question: what is KREEP? Because pristine KREEP basalts are rare, and most known samples are small (e.g., 15382/15386), the geochemical variability of KREEP basalts is poorly known. Other KREEP compositions which are commonly used in these models include the hypothetical 'high-K KREEP' component of Warren and Wasson, which is derived from Apollo 14 soil data, and the 'superKREEP' quartz-monzodiorite 15405. Lunar breccia 15205 is a polymict regolith breccia that consists of approximately 20% KREEP basalt clasts and 20% quartz-normative basalt clasts in a KREEP-rich matrix. Bulk rock mixing calculations show that this sample comprises about 84% KREEP. The clasts range up to 1 cm in size, but most are considerably smaller. The primary aim is to characterize pristine KREEP basalts petrographically, to establish the range in chemical compositions of KREEP basalts, and to test models that were proposed for their origin. In addition, we may be able to extend the compositional range recognized in the quartz-normative basalt suite and cast some light on its origin as well. Preliminary whole rock geochemical data on the KREEP basalts are presented in a companion paper by M.M. Lindstrom and co-workers. Concentration is on petrography and mineral chemistry of these clasts, and the implications these data have for the origin of the different melt rock suites.

Shervais, John W.↗

DOE OSTI2020

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Post-landing major element quantification using SuperCam laser induced breakdown spectroscopy

The SuperCam instrument on the Perseverance Mars 2020 rover uses a pulsed 1064 nm laser to ablate targets at a distance and conduct laser induced breakdown spectroscopy (LIBS) by analyzing the light from the resulting plasma. SuperCam LIBS spectra are preprocessed to remove ambient light, noise, and the continuum signal present in LIBS observations. Prior to quantification, spectra are masked to remove noisier spectrometer regions and spectra are normalized to minimize signal fluctuations and effects of target distance. In some cases, the spectra are also standardized or binned prior to quantification. To determine quantitative elemental compositions of diverse geologic materials at Jezero crater, Mars, we use a suite of 1198 laboratory spectra of 334 well-characterized reference samples. The samples were selected to span a wide range of compositions and include typical silicate rocks, pure minerals (e.g., silicates, sulfates, carbonates, oxides), more unusual compositions (e.g., Mn ore and sodalite), and replicates of the sintered SuperCam calibration targets (SCCTs) onboard the rover. For each major element (SiO 2 , TiO 2 , Al 2 O 3 , FeO T , MgO, CaO, Na 2 O, K 2 O), the database was subdivided into five “folds” with similar distributions of the element of interest. One fold was held out as an independent test set, and the remaining four folds were used to optimize multivariate regression models relating the spectrum to the composition. We considered a variety of models, and selected several for further investigation for each element, based primarily on the root mean squared error of prediction (RMSEP) on the test set, when analyzed at 3 m. In cases with several models of comparable performance at 3 m, we incorporated the SCCT performance at different distances to choose the preferred model. Shortly after landing on Mars and collecting initial spectra of geologic targets, we selected one model per element. Subsequently, with additional data from geologic targets, some models were revised to ensure results that are more consistent with geochemical constraints. The calibration discussed here is a snapshot of an ongoing effort to deliver the most accurate chemical compositions with SuperCam LIBS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗