SEARCH · Search NASA
Results for “Mineral dissolution”
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.
Atmospheric Dissolved Iron Depostiion to the Global Oceans: Effects of Oxalate-Promoted Fe Dissolution, Photochemical Redox Cycling, and Dust Mineralogy
Mineral dust deposition is suggested to be a significant atmospheric supply pathway of bioavailable iron (Fe) to Fe-depleted surface oceans. In this study, mineral dust and dissolved Fe (Fed) deposition rates are predicted for March 2009 to February 2010 using the 3-D chemical transport model GEOS-Chem implemented with a comprehensive dust-Fe dissolution scheme. The model simulates Fed production during the atmospheric transport of mineral dust taking into account inorganic and organic (oxalate)-promoted Fe dissolution processes, photochemical redox cycling between ferric (Fe(III)) and ferrous (Fe(II)) forms of Fe, dissolution of three different Fe-containing minerals (hematite, goethite, and aluminosilicates), and detailed mineralogy of windblown dust from the major desert regions. Our calculations suggest that during the yearlong simulation is approximately 0.26 Tg (1 Tg = 1012 g) of Fed was deposited to global oceanic regions. Compared to simulations only taking into account proton-promoted Fe dissolution, the addition of oxalate to the dust-Fe mobilization scheme increased total annual model-predicted Fed deposition to global oceanic regions by approximately 75%. The implementation of Fe(II)/Fe(III) photochemical redox cycling in the model allows for the distinction between different oxidation states of deposited Fed. Our calculations suggest that during the daytime, large fractions of Fed deposited to the global oceans is likely to be in Fe(II) form, while nocturnal fluxes of Fed are largely in Fe(III) form. Model simulations also show that atmospheric fluxes of Fed can be strongly influenced by the mineralogy of Fe-containing compounds. This study shows that Fed deposition to the oceans is controlled by total dust-Fe mass concentrations, mineralogy, the surface area of dust particles, atmospheric chemical composition, cloud processing, and meteorological parameters and exhibits complex and spatiotemporally variable patterns. Our study suggests that the explicit model representation of individual processes leading to Fed production within mineral dust are needed to improve the understanding of the atmospheric Fe cycle, and quantify the effect of dust-Fe on ocean biological productivity, carbon cycle, and climate.
Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale
Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.
Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface
Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.
Coatings on Atacama Desert Basalt: A Possible Analog for Coatings on Gusev Plains Basalt
Surface coatings on Gusev Plains basalt have been observed and may contain hematite and nanophase Fe-oxides along with enrichments in P, S, Cl, and K relative to the underlying rock. The Gusev coatings may be derived from the dissolution of adhering soil and/or parent rock along with the addition of S and Cl from outside sources. Transient water for dissolution could be sourced from melting snow during periods of high obliquity, acid fog, and/or ground water (Haskin et al., 2005). Coatings on basalt in the hyper-arid (less than 2mm y(sup -1)) Atacama Desert may assist in understanding the chemistry, mineralogy and formation mechanisms of the Gusev basalt coatings. The Atacama Desert climate is proposed to be analogous to a paleo-Mars climate that was characterized by limited aqueous activity when the Gusev coatings could have formed. The objectives of this work are to (i) determine the chemical nature and extent of surface coatings on Atacama Desert basalt, and (ii) assess coating formation mechanisms in the Atacama Desert. Preliminary backscattered electron imaging of Atacama basalt thin-sections indicated that the coatings are as thick as 20 m. The boundary between the coating and the basalt labradorite, ilmenite, and augite grains was abrupt indicating that the basalt minerals underwent no chemical dissolution. The Atacama coatings have been added to the basalt instead of being derived from basalt chemical weathering. Semi-quantitative energy dispersive spectroscopy shows the coatings to be chemically homogeneous. The coating is depleted in Ca (0.9 wt% CaO) and enriched in K (1.3 wt.% K2O) and Si (69.1 wt.% SiO2) relative to the augite and labradorite grains. A dust source enriched in Si (e.g., poorly crystalline silica) and K and depleted in Ca appears to have been added to the basalt surface. Unlike the Gusev coatings, no P, S, and Cl enrichment was observed. However, Fe (3.2 wt.% FeO) was present in the Atacama coatings suggesting the present of Fe-oxides. While the chemistry of Atacama coating does not mirror the Gusev coating, the coating formation mechanism may be similar. The Atacama coatings of surface basalt are derived completely from exogenous sources. If surface Mars rocks have experienced limited wetting conditions as in the Atacama, then Mars coatings may be derived only from dissolution of material adhering to rock.
Lithological controls on soil geochemistry and clay mineralogy across Spodosols in the coastal temperate rainforest of southeast Alaska
The interaction of similar soil forming factors with different parent materials determines soil chemical weathering, influencing soil processes and properties. In the coastal temperate rainforests of southeast Alaska, Spodosols is the dominant soil order in well-drained parent materials regardless of the lithology. Yet, the role of lithology on chemical weathering, base cation depletion, and mineral transformation during pedogenesis remains elusive. Here, we established a lithosequence comprised of soils derived from tonalite, phyllite, slate, and metavolcanic rocks to test the hypothesis that despite the soils in southeast Alaska presenting similar taxonomy and morphology, the influence of lithology in chemical weathering intensity and mineralogy can be detected. We evaluated physicochemical properties, clay mineralogy, Fe-Al oxides, chemical weathering, and elemental mass balance on eleven Spodosols sampled along the lithosequence. We also propose a new weathering index, the Weathering Index for Spodosols (WISP), suitable to evaluate the weathering degree of soils with Al mobility and base-rich parent materials. This index fills a gap in traditional weathering indices that uses Al as an immobile element and/or does not evaluate the leaching of main base cations. We found the pedons across the lithosequence expressed similar physicochemical properties, with predominantly thick profiles, presenting andic properties, relatively thick spodic horizons, and thin E horizons. Podzolization imposed similar mineral transformation trends, mainly dissolution of chlorite, mica alteration to interstratified mica-vermiculite, and formation of smectite in E horizons. Kaolinite was detected in the soils from all lithologies, but only in trace amounts in phyllite soils. Fe oxides depth distribution was similar but with a higher concentration in the slate (44.5 ± 4.2 g/kg), followed by metavolcanic (31.6 ± 4.0 g/kg), phyllite (27.9 ± 3.8 g/kg), and tonalite soils (19.9 ± 5.7 g/kg) and predominantly as organometallic complexes (extracted by Na-pyrophosphate). The stronger weathering of chlorite and mica in the E horizons was reflected in larger losses of Mg and K, while Ca and Na were more depleted in the spodic horizons, suggesting plagioclase weathering. Slate soils were the most depleted in base cations, followed by phyllite, tonalite, and metavolcanic soils. The more depleted status of slate soils was supported by the higher WISP values (WISP = 62), followed by phyllite (WISP = 55), tonalite (WISP = 53), and metavolcanic soils (WISP = 45). Our results demonstrated that lithology controlled elemental depletion intensities, Fe oxide concentrations, and partially the clay fraction mineralogy in addition to podzolization acting as a dominant pedogenic process across all sites. These findings advance our understanding of the role of lithology on soil mineralogy/geochemistry that impacts critical soil functions, such as the soil carbon cycle, elemental fluxes, organo-mineral interactions, and solid-solution reactions.
Organic Ligand-Mediated Dissolution and Fractionation of Rare-Earth Elements (REEs) from Carbonate and Phosphate Minerals
Not Available
Potential Chemical Impacts of Subsurface CO 2 : An Integrated Experimental and Numerical Assessment for a Case Study of the Ogallala Aquifer
Leakage from geologic CO 2 sequestration (GCS) reservoirs to overlying underground sources of drinking water (USDW) is a tangible risk. This study is an integrated assessment that combines column experiments and reactive transport simulations of sediments sampled from the Ogallala aquifer above an active commercial-scale GCS site (the Farnsworth Unit in northern Texas). Experimental and simulation results suggest that carbonate mineral (calcite and dolomite) dissolution is the most significant reaction following CO 2 intrusion, and is also the dominant source of trace metal release. Cation exchange is another key mechanism controlling trace metal release by cation interference. Most of the trace metals, including Ba, Sr, As, Pb, and Zn, show a short-term release and quickly drop to the baseline values, suggesting low risk to the overlying USDW quality. Other trace metals, such as Mn and U, exhibit a tangible increase of their concentrations in the beginning, and drop to a higher level compared to the baseline, which may become a potential concern of long-term USDW quality change with CO2introduction. This study provides a comprehensive example of combining laboratory experiments and simulations for assessment of CO 2 -sediment interactions with combined release mechanisms in shallow groundwater aquifers. Data presented here provides useful insights for quantitative risk assessment and effective public education regarding CO 2 geological sequestration, and trace metal reactive transport studies in shallow groundwater aquifers.
Overexpression of sulfide:quinone reductase (SQR) in Acidithiobacillus ferrooxidans enhances sulfur, pyrite, and pyrrhotite oxidation
ABSTRACT Hydrogen sulfide is produced during the dissolution of some sulfidic minerals and during the microbial metabolism of reduced sulfur compounds. The sulfide:quinone reductase (SQR) enzyme is able to oxidize H 2 S, and the bioleaching cells Acidithiobacillus ferrooxidans have two SQR genes, only one of which has been characterized. We cloned and overexpressed the two SQR genes in A. ferrooxidans and show that they both have SQR activity. Both AFE_0267 and AFE_1792 are active under anaerobic conditions, but only AFE_1792 is active under aerobic conditions. The effect of the SQR overexpression and the expression of related genes on sulfur metabolism was investigated. The overexpression of SQR improved cell growth and sulfur oxidation, suggesting enhanced SQR activity led to a reduction in H 2 S toxicity as well as providing additional energy through H 2 S oxidation. Additionally, the impact on the oxidation of pyrite and pyrrhotite was investigated. The rate of oxidation of pyrite by the engineered cells was enhanced, and, furthermore, the rate of pyrrhotite oxidation was more than doubled. IMPORTANCE H 2 S is a toxic sulfur intermediate, and the SQR enzyme has evolved to oxidize H 2 S in A. ferrooxidans . In addition to detoxification, H 2 S oxidation provides energy, and overexpression of SQR enhanced aerobic and anaerobic growth on sulfur. The SQR overexpression also enhanced pyrite and pyrrhotite oxidation, which may facilitate the pyrometallurgical processing of a number of critical materials including copper, nickel, and the platinum group metals.
Geochemically-Driven Petrophysical Transformations in Caney Shale and Their Impact on Reservoir Productivity
Geochemically-Driven Petrophysical Transformations in Caney Shale and Their Impact on Reservoir Productivity
Resolving Nanoscale Processes during Carbon Mineralization Using Identical Location Transmission Electron Microscopy
Basalt reservoirs offer the potential for carbon mineralization, aiding in achieving net-zero emissions. However, debates persist about microscopic crystallization mechanisms due to limited characterization techniques under high-temperature and pressure conditions. By using Identical Location Transmission Electron Microscopy (IL-TEM) and cryo-TEM, this study reveals nanoscale interfacial carbonation processes of forsterite and diopside nanoparticles in water-saturated supercritical carbon dioxide under realistic reservoir conditions. Further, both minerals undergo preferential metal cation dissolution into a thin water film, forming porous Si-rich amorphous layers, supporting the leached layer mechanism as the dominant mineral reactivity process. Diopside’s amorphous layer has lower porosity and growth rate relative to forsterite, likely related to the connectivity of silicate tetrahedra. Kinetically favorable nesquehonite and aragonite nanocrystals form on the amorphous layers. These findings support the development of accurate reservoir simulations and help enable commercial-scale carbon storage deployment.
Morphological and crystallographic controls in the replacement of calcite and aragonite by cerussite and otavite
Mineral replacement reactions are essential solid-fluid interactions in natural and industrial processes including metasomatism, diagenesis, and metamorphism, and sequestration of metal toxins from polluted water. In this paper, we explore the morphological evolution of aragonite and calcite (polymorphs of CaCO 3 ) during their replacement by cerussite (PbCO 3 , an isomorph of aragonite) and otavite (CdCO 3 , an isomorph of calcite) in acidic aqueous solutions (with initial pH 3 to 4). Observations using scanning electron and synchrotron-based transmission X-ray microscopies reveal the formation of pseudomorphic shells (~5 to 10 µm thick) of cerussite and otavite when Pb 2+ is reacted with either calcite or aragonite and when Cd 2+ is reacted with aragonite. The formation of pore space at the substrate-precipitate interface and within the precipitate is found to be key in promoting chemical exchange between the dissolving phase and the solution. In contrast, minimal reactivity is observed when Cd 2+ is reacted with calcite. These results demonstrate that the lack of substrate-precipitate epitaxy (either because of inconsistent symmetries or due to large lattice constant mismatch) is a critical factor that enables replacement of primary CaCO 3 minerals by Pb- or Cd-bearing secondary minerals while the presence of epitaxial relationships inhibit those reactions. When calcite and aragonite coexist, the replacement of aragonite by cerussite and otavite occurs preferentially over calcite. We postulate that the relative reactivity of calcite and aragonite to Pb 2+ and Cd 2+ can be determined by interplay between morphological and crystallographic properties of CaCO 3 materials as well as solution chemistry. Our results demonstrate that coupled dissolution-precipitation processes of carbonate minerals play an important role in controlling mobility and bioavailability of lead and cadmium in environments where acidic metal-contaminated water interacts with carbonate-rich rocks and soils.
Replacement of Calcium Carbonate Polymorphs by Cerussite
Calcium carbonate (CaCO 3 ) polymorphs, calcite, aragonite, and vaterite, serve as a major sink to retain various metal ions in natural and engineered systems. Here, we visualize the systematic trends in reactivities of calcite, vaterite, and aragonite to Pb 2+ dissolved in acidic aqueous solutions using in situ optical microscopy combined with ex situ scanning electron and transmission X-ray microscopies. All three polymorphs undergo pseudomorphic replacement by cerussite (PbCO 3 ) but with distinct differences in the evolution of their morphologies. The replacement of calcite and aragonite occurs through the formation of a pseudomorphic cerussite shell (typically 5–10 μm thick) followed by a slower inward propagation of reaction fronts through a thin solution gap (~0.1 μm wide) between the shell and the CaCO 3 substrate. The replacement of vaterite is characterized by the formation of a thinner cerussite shell (≤1 μm thick) and a larger cavity between the shell and the host mineral. These systematic differences in cerussite morphology for different CaCO 3 polymorphs are explained by the relative dissolution and precipitation rates of the reactant and product minerals, coupled with the role of ion transport through the cerussite shells. We also find that the replacement of calcite by cerussite is the slowest when all three polymorphs coexisted. Furthermore, our results provide mechanistic insights into the growth mode of cerussite on dissolving calcium carbonate and demonstrate these CaCO 3 polymorphs as promising substrate materials for removal and recycling of Pb from acidic polluted water and industrial effluents.
Rates of oxidative weathering on the surface of Mars
Implicit in the mnemonic 'MSATT' (Mars surface and atmosphere through time) is that rates of surface processes on Mars through time should be investigated, including studies of the kinetics and mechanism of oxidative weathering reactions occurring in the Martian regolith. Such measurements are described. Two major elements analyzed in the Viking Lander XRF experiment that are most vulnerable to atmospheric oxidation are iron and sulfur. Originally, they occurred as Fe(2+)-bearing silicate and sulfide minerals in basaltic rocks on the surface of Mars. However, chemical weathering reactions through time have produced ferric- and sulfate-bearing assemblages now visible in the Martian regolith. Such observations raise several question about: (1) when the oxidative weathering reactions took place on Mars; (2) whether or not the oxidized regolith is a fossilized remnant of past weathering processes; (3) deducting chemical interactions of the ancient Martian atmosphere with its surface from surviving phases; (4) possible weathering reactions still occurring in the frozen regolith; and (5) the kinetics and mechanism of past and present-day oxidative reactions on Mars. These questions may be addressed experimentally by studying reaction rates of dissolution and oxidation of basaltic minerals, and by identifying reaction products forming on the mineral surfaces. Results for the oxidation of pyrrhotite and dissolved ferrous iron are reported.
Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes
Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.
Impact of aquifer properties on the extent and timeline of CO 2 trapping
Abstract Geologic CO 2 sequestration in porous saline aquifers is a promising approach to reducing atmospheric concentrations of CO 2 . Reactive transport simulations provide the opportunity to analyze which factors influence geochemical reactivity in the reservoir, understand those most important for promoting CO 2 trapping, and assess individual sites. Field‐scale aquifer characterization is time and resource intensive such that here, reactive transport simulations are leveraged to enhance understanding of selected aquifer properties including porosity, permeability, depth of storage, and carbonate mineralogy on the overall CO 2 trapping potential to better select sites promoting geochemical reactivity for CO 2 trapping. There are different mechanisms for sequestrating CO 2 . Once injected, CO 2 will dissolve into the brine to create an acidic environment, resulting in the dissolution of pre‐injection formation minerals. Released ions can reprecipitate as secondary minerals. The dissolved CO 2 and mineralized CO 2 are considered as a more secure form of CO 2 trapping in this study compared to the free supercritical CO 2 . Here, a framework leveraging a controlled set of field scale simulations is developed to facilitate rapid, optimized site selection. Simulations vary aquifer properties to understand the impact of each unique property on CO 2 trapping, tracking, and comparing the amount of supercritical, aqueous, and mineralized CO 2 . The rate at which the CO 2 injected into the aquifer is converted to aqueous or mineralized CO 2 is newly defined here as the sequestration efficiency and used to compare simulation results. The reservoir depth and fraction of carbonate minerals in the formation are shown to be more important factors than reservoir porosity and permeability in affecting CO 2 trapping. However, the impact of aquifer properties on the evolution of injected CO 2 depends on the stage of the sequestration project. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.
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.