Search NASA⌕ Search

SEARCH · Search NASA

Results for “ultramafics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Isotopic fractionation accompanying CO 2 hydroxylation and carbonate precipitation from high pH waters at The Cedars, California, USA

The Cedars ultramafic block hosts alkaline springs (pH > 11) in which calcium carbonate forms upon uptake of atmospheric CO 2 and at times via mixing with surface water. Here, these processes lead to distinct carbonate morphologies with “floes” forming at the atmosphere-water interface, “snow” of fine particles accumulating at the bottom of pools and terraced constructions of travertine. Floe material is mainly composed of aragonite needles despite CaCO 3 precipitation occurring in waters with low Mg/Ca (<0.01). Precipitation of aragonite is likely promoted by the high pH (11.5–12.0) of pool waters, in agreement with published experiments illustrating the effect of pH on calcium carbonate polymorph selection.

58 GEOSCIENCES↗
DOE OSTI2025

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

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

08 - HYDROGEN↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2022

Exotic Carbonate Mineralization Recovered from a Deep Basalt Carbon Storage Demonstration

Mitigating climate change require transformational advances for carbon dioxide removal, including geologic carbon sequestration in reactive subsurface environments. The Wallula Basalt Carbon Storage Pilot Project demonstrated that CO 2 injected into >800 m deep Columbia River Basalt Group flow top reservoirs mineralizes on month-year timescales. Herein, we present new optical petrography, micro-computed X-ray tomography, and electron microscopy results sidewall cores collected two years after CO 2 injection. As no other anthropogenic carbonates from geologic carbon storage field studies have been recovered, this world-unique sample suite provides unparalleled insight for subsurface carbon mineralization products and paragenesis. Chemically-zoned nodules with Ca/Mn-rich cores and Fe-dominant outer rims are prominent examples of the neoformed carbonate assemblages with ankerite-siderite compositions and exotic divalent cation correlations. Further, paragenetic insights for the timing of aragonite, silica, and fibrous zeolites are clarified based on mineral texture and spatial relationships, along with time-resolved downhole fluid sampling. Collectively, these results clarify the mineralogy, chemistry, and paragenesis of carbon mineralization, providing insight into the ultimate fate and transport of CO 2 in reactive mafic-ultramafic reservoirs

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dithionite Inhibits Iron(III) (Hydr)oxide Formation during Olivine Dissolution Advancing Simultaneous CO 2 Mineralization and Nickel Recovery

Increasing CO 2 concentration poses significant global challenges, impacting both environmental and human health. As we strive for a carbon-neutral energy technology transition, the demand for critical elements (e.g., nickel and cobalt) continues to increase, while high-grade ores are depleting. Combining CO 2 mineralization with the recovery of critical elements from low-grade ores offers an innovative solution. Olivine, a magnesium-rich ultramafic with trace amounts of critical elements, is a promising mineral; however, its impurities, such as iron (Fe(III)), hinder the dissolution of Ni and Mg from olivine, reducing its carbonation and critical element recovery. Here, to address this challenge, this study examined Ni dissolution from San Carlos olivine at high temperatures and high CO 2 pressure. We found that iron(III) (hydr)oxide layers impeded olivine dissolution; however, with sodium dithionite (Na 2 S 2 O 4 ), a reducing agent, olivine dissolution is significantly improved by preventing iron(III) (hydr)oxide formation. With Na 2 S 2 O 4 , within 24 h, Mg and Ni dissolution from olivine increased 2.85-fold and 2.66-fold, respectively, compared to samples without Na 2 S 2 O 4 . After seven cycles, with solutions replaced every 24 h, 98.9% of the total Mg and 84.6% of the total Ni were recovered. This approach enhances olivine’s CO 2 mineralization and improves the sustainability of critical element supply.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Enhanced Recovery of Critical Minerals and Geological Hydrogen in The Mine of the Future

invited perspective article: The rapid demand for critical minerals (CMs) and hydrogen (H2) necessitates innovative solutions beyond conventional mining. This study explores enhanced mineral recovery (EMR) and geological hydrogen production (GeoH2) by leveraging the Earth’s subsurface as a reactive platform. Using ultramafic rocks and engineered fluids, the approach simultaneously mobilizes critical resources and produces H2 through natural processes like serpentinization. This paradigm offers a transformative pathway to secure essential materials and diversify energy resources, setting a foundation for the Terrestrial Mine of the Future

Yan, Keju↗
DOE OSTI2022

Elastic Anisotropy of Lizardite at Subduction Zone Conditions

Subduction zones transport water into Earth's deep interior through slab subduction. Serpentine minerals, the primary hydration product of ultramafic peridotite, are abundant in most subduction zones. Characterization of their high-temperature elasticity, particularly their anisotropy, will help us better estimate the extent of mantle serpentinization and the Earth's deep water cycle. Lizardite, the low-temperature polymorph of serpentine, is stable under the P-T conditions of cold subduction slabs (<260°C at 2 GPa), and its high-temperature elasticity remains unknown. Here we report ab initio elasticity and acoustic wave velocities of lizardite at P-T conditions of subduction zones. Our static results agree with previous studies. Its high-temperature velocities are much higher than previous experimental-based lizardite estimates with chrysotile but closer to antigorite velocities. The elastic anisotropy of lizardite is much larger than that of antigorite and could better account for the observed large shear-wave splitting in some cold slabs such as Tonga.

58 GEOSCIENCES↗
DOE OSTI2022

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2024

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO 2 ) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ 34 S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H 2 O-NaCl solutions to calibrate a high-pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO 2 from FMQ-2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO 2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur-rich and high δ 34 S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34 S-rich subducted lithologies.

58 GEOSCIENCES↗
DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2024

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

58 GEOSCIENCES↗
DOE OSTI2025

Complex carbonate phases drive geologic CO2 mineralization

Abstract Geologic carbon sequestration in mafic and ultramafic reservoirs is a scalable strategy for carbon dioxide removal, offering permanent storage via mineralization as stable carbonates. However, there is limited information on the structure and composition of key mineralization endpoints during sequestration. Here, we unravel the atomic structure, composition, and nanoscale morphology of carbonates recovered from a field-scale demonstration of CO 2 mineralization in basalt. Using transmission electron microscopy, we mapped mineralogical variations from the initial to later stages of subsurface carbonate growth and identified a previously unknown cation-ordered ankerite phase that exerts a primary control over carbonation processes. This study has provided a new understanding of subsurface carbonation pathways which will impact the parameterization of predictive geochemical models for future sequestration efforts in basalt formations.

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

A dynamic microbial sulfur cycle in a serpentinizing continental ophiolite

Summary Serpentinization is the hydration and oxidation of ultramafic rock, which occurs as oceanic lithosphere is emplaced onto continental margins (ophiolites), and along the seafloor as faulting exposes this mantle‐derived material to circulating hydrothermal fluids. This process leads to distinctive fluid chemistries as molecular hydrogen (H 2 ) and hydroxyl ions (OH − ) are produced and reduced carbon compounds are mobilized. Serpentinizing ophiolites also serve as a vector to transport sulfur compounds from the seafloor onto the continents. We investigated hyperalkaline, sulfur‐rich, brackish groundwater in a serpentinizing continental ophiolite to elucidate the role of sulfur compounds in fuelling in situ microbial activities. Here we illustrate that key sulfur‐cycling taxa, including Dethiobacter , Desulfitispora and ‘ Desulforudis ’, persist throughout this extreme environment. Biologically catalysed redox reactions involving sulfate, sulfide and intermediate sulfur compounds are thermodynamically favourable in the groundwater, which indicates they may be vital to sustaining life in these characteristically oxidant‐ and energy‐limited systems. Furthermore, metagenomic and metatranscriptomic analyses reveal a complex network involving sulfate reduction, sulfide oxidation and thiosulfate reactions. Our findings highlight the importance of the complete inorganic sulfur cycle in serpentinizing fluids and suggest sulfur biogeochemistry provides a key link between terrestrial serpentinizing ecosystems and their submarine heritage.

Sabuda, Mary C.↗
DOE OSTI2023

Osmotic Pressure Gradient Effects on Water Diffusion in Porous Rock: Can They Pervert Permeability Tests?

Abstract Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.

Mechanics↗
DOE OSTI2021

Surface-wave tomography of the Emeishan large igneous province (China): Magma storage system, hidden hotspot track, and its impact on the Capitanian mass extinction

Large igneous provinces (LIPs) are commonly associated with mass extinctions. However, the precise relations between LIPs and their impacts on biodiversity is enigmatic, given that they can be asynchronous. It has been proposed that the environmental impacts are primarily related to sill emplacement. Therefore, the structure of LIPs' magma storage system is critical because it dictates the occurrence and timing of mass extinction. We use surface-wave tomography to image the lithosphere under the Permian Emeishan large igneous province (ELIP) in southwestern China. We find a northeast-trending zone of high shear-wave velocity (Vs) and negative radial anisotropy (Vsv > Vsh; v and h are vertically and horizontally polarized S waves, respectively) in the crust and lithosphere. We rule out the possibilities of rifting or orogenesis to explain these seismic characteristics and interpret the seismic anomaly as a mafic-ultramafic, dike-dominated magma storage system of the ELIP. We further propose that the anomaly represents a hidden hotspot track that was emplaced before the ELIP eruption. A zone of higher velocity but less-negative radial anisotropy, on the hotspot track but to the northeast of the eruption center in the Panxi region, reflects an elevated proportion of sills emplaced at the incipient stage of the ELIP. Liberation of poisonous gases by the early sill intrusions explains why the mid-Capitanian global biota crisis preceded the peak ELIP eruption by 2–3 m.y.

58 GEOSCIENCES↗
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↗