skip to main content


Title: Effects of Charged Solute-Solvent Interaction on Reservoir Temperature during Subsurface CO2 Injection
A short-term side-effect of CO2 injection is a developing low-pH front that forms ahead of the bulk water injectant, due to differences in solute diffusivity. Observations of downhole well temperature show a reduction in aqueous-phase temperature with the arrival of a low-pH front, followed by a gradual rise in temperature upon the arrival of a high concentration of bicarbonate ion. In this work, we model aqueous-phase transient heat advection and diffusion, with the volumetric energy generation rate computed from solute-solvent interaction using the Helgeson–Kirkham–Flowers (HKF) model, which is based on the Born Solvation model, for computing specific molar heat capacity and the enthalpy of charged electrolytes. A computed injectant water temperature profile is shown to agree with the actual bottom hole sampled temperature acquired from sensors. The modeling of aqueous-phase temperature during subsurface injection simulation is important for the accurate modeling of mineral dissolution and precipitation because forward dissolution rates are governed by a temperature-dependent Arrhenius model.  more » « less
Award ID(s):
2019194
NSF-PAR ID:
10354628
Author(s) / Creator(s):
Date Published:
Journal Name:
Minerals
Volume:
12
Issue:
6
ISSN:
2075-163X
Page Range / eLocation ID:
752
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Secondary minerals in martian nakhlites provide a powerful tool for investigating the nature, composition, and duration of aqueous activity in the martian crust. Northwest Africa (NWA) 998 crystallized early from the nakhlite magmatic source and has evidence of minimal signatures of the late hydrothermal alteration event that altered the nakhlites. Using FIB‐TEM techniques to study a cumulus apatite grain in NWA 998, we report the first evidence of a submicron‐scale vein consisting of fluorapatite and an SiO2‐rich phase. Fluorapatite grew epitaxially on the walls of an opened cleavage plane of host F‐bearing chlorapatite and the SiO2‐rich phase filled the center of the vein. The presence of nanoporosity and nanometer‐scale amorphous material and the sharp interface between the vein and the host apatite indicate the vein represents a coupled dissolution–reprecipitation process that generated apatite of a different composition that was more stable with the fluid. Using experimental data and diffusion coefficients of Cl in apatite from the literature, we conclude that the vein was caused by a low temperature (~300°C), slightly acidic, F‐, Si‐rich, aqueous fluid that acted as a closed system. Based on the characteristics of the vein (formation by rapid injection of fluid) and the fluid (composition, temperature, pH), and the lack of terrestrial weathering products in our SEM and TEM images, we infer that the vein is pre‐terrestrial in origin. Our observations support the hypothesis that the heat source triggering a hydrothermal system was a low‐shock velocity impact and rule out a magmatic origin. Finally, the vein could have formed from a late‐stage fluid different from that reported in other nakhlites, but formation during the same magmatic event by, for example, a less evolved fluid might also be plausible.

     
    more » « less
  2. null (Ed.)
    Abstract The McMurdo Dry Valley region is the largest ice-free area of Antarctica. Ephemeral streams flow here during the austral summer, transporting glacial meltwater to perennially ice-covered, closed basin lakes. The chemistry of 24 Taylor Valley streams was examined over the two-decade period of monitoring from 1993 to 2014, and the geochemical behavior of two streams of contrasting physical and biological character was monitored across the seven weeks of the 2010–2011 flow season. Four species dominate stream solute budgets: HCO3–, Ca2+, Na+, and Cl–, with SO42–, Mg2+, and K+ present in significantly lesser proportions. All streams contain dissolved silica at low concentrations. Across Taylor Valley, streams are characterized by their consistent anionic geochemical fingerprint of HCO3 > Cl > SO4, but there is a split in cation composition between 14 streams with Ca > Na > Mg > K and 10 streams with Na > Ca > Mg > K. Andersen Creek is a first-order proglacial stream representative of the 13 short streams that flow <1.5 km from source to gage. Von Guerard is representative of 11 long streams 2–7 km in length characterized by extensive hyporheic zones. Both streams exhibit a strong daily cycle for solute load, temperature, dissolved oxygen, and pH, which vary in proportion to discharge. A well-expressed diurnal co-variation of pH with dissolved oxygen is observed for both streams that reflects different types of biological control. The relative consistency of Von Guerard composition over the summer flow season reflects chemostatic regulation, where water in transient storage introduced during times of high streamflow has an extended opportunity for water-sediment interaction, silicate mineral dissolution, and pore-water exchange. 
    more » « less
  3. Abstract

    Projections of future conditions within the critical zone—earthcasts—can be used to understand the potential effects of changes in climate on processes affecting landscapes. We are developing an approach to earthcast how weathering will change in the future using scenarios of climate change. As a first step here, we use the earthcasting approach to model aspect‐related effects on soil water chemistry and weathering on hillsides in a well‐studied east‐west trending watershed (Shale Hills, Pennsylvania, USA). We completed model simulations of solute chemistry in soil water with and without the effect of aspect for comparison to catchment observations. With aspect included, aqueous weathering fluxes were higher on the sunny side of the catchment. But the effect of aspect on temperature (0.8 °C warmer soil on sunny side) and recharge (100 mm/year larger on shaded side) alone did not explain the magnitude of the observed higher weathering fluxes on the sunny side. Modeled aspect‐related differences in weathering fluxes only approach field observations when we incorporated the measured differences in clay content observed in augered soils on the two hillslopes. We also had to include a biolifting module to accurately describe cation concentrations in soil water versus depth. Biolifting lowered some mineral dissolution rates while accelerating kaolinite precipitation. These short‐duration simulations also highlighted that the inherited differences in particle size on the two sides of the catchment might in themselves be explained by weathering under different microclimates caused by aspect—over longer durations than simulated with our models.

     
    more » « less
  4. Lead (Pb) solubility is commonly limited by dissolution–precipitation reactions of secondary mineral phases in contaminated soils and water. In the research described here, Pb solubility and free Pb2+ ion activities were measured following the precipitation of Pb minerals from aqueous solutions containing sulfate or carbonate in a 1:5 mole ratio in the absence and presence of phosphate over the pH range 4.0–9.0. Using X-ray diffraction and Fourier-transform infrared spectroscopic analysis, we identified anglesite formed in sulfate-containing solutions at low pH. At higher pH, Pb carbonate and carbonate-sulfate minerals, hydrocerussite and leadhillite, were formed in preference to anglesite. Precipitates formed in the Pb-carbonate systems over the pH range of 6 to 9 were composed of cerussite and hydrocerussite, with the latter favored only at the highest pH investigated. The addition of phosphate into the Pb-sulfate and Pb-carbonate systems resulted in the precipitation of Pb3(PO4)2 and structurally related pyromorphite minerals and prevented Pb sulfate and carbonate mineral formation. Phosphate increased the efficiency of Pb removal from solution and decreased free Pb2+ ion activity, causing over 99.9% of Pb to be precipitated. Free Pb2+ ion activities measured using the ion-selective electrode revealed lower values than predicted from thermodynamic constants, indicating that the precipitated minerals may have lower KSP values than generally reported in thermodynamic databases. Conversely, dissolved Pb was frequently greater than predicted based on a speciation model using accepted thermodynamic constants for Pb ion-pair formation in solution. The tendency of the thermodynamic models to underestimate Pb solubility while overestimating free Pb2+ activity in these systems, at least in the higher pH range, indicates that soluble Pb ion-pair formation constants and KSP values need correction in the models. 
    more » « less
  5. Hydrogen is a rapidly diffusing monovalent cation in nominally anhydrous minerals (NAMs, such as olivine, orthopyroxene, and clinopyroxene), which is potentially re-equilibrated during silicate melt-rock and aqueous fluid-rock interactions in massif and abyssal peridotites. We apply a 3D numerical diffusion modeling technique to provide first-order timescales of complete hydrogen re-equilibration in olivine, clinopyroxene, and orthopyroxene over the temperature range 600-1200°C. Model crystals are 1-3 mm along the c-axis and utilize H+ diffusion coefficients appropriate for Fe-bearing systems. Two sets of models were run with different boundary compositions: 1) “low-H models” are constrained by mineral-melt equilibrium partitioning with a basaltic melt that has 0.75 wt% H2O and 2) “high-H models,” which utilize the upper end of the estimated range of mantle water solubility for each phase. Both sets of models yield re-equilibration timescales that are identical and are fast for all phases at a given temperature. These timescales have strong log-linear trends as a function of temperature (R2 from 0.97 to 0.99) that can be used to calculate expected re-equilibration time at a given temperature and grain size. At the high end of the model temperatures (1000-1200°C), H+ completely re-equilibrates in olivine, orthopyroxene, and clinopyroxene within minutes to hours, consistent with previous studies. These short timescales indicate that xenolith NAM mantle water contents are likely to be overprinted prior to eruption. The models also resolve the decoupled water-trace element relationship in Southwest Indian Ridge peridotites, in which peridotite REE abundances are reproduced by partial melting models whereas the relatively high NAM H2O contents require later re-equilibration with melt. At temperatures of 600-800°C, which correspond to conditions of hydrothermal alteration of pyroxene to amphibole and talc, H+ re-equilibration typically occurs over a range of timescales spanning days to years. These durations are well within existing estimates for the duration of fluid flow in oceanic hydrothermal systems, suggesting that peridotite NAM water contents are susceptible to diffusive overprinting during higher temperature hydrothermal alteration. Thus, diffusion during aqueous fluid-rock interactions may also explain NAM H2O contents that are too high to reflect residues of melting. These relatively short timescales at low temperatures suggest that the origin of water contents measured in peridotite NAMs requires additional constraints on sample petrogenesis, including petrographic and trace element analyses. Our 3D model results also hint that H+ may diffuse appreciably during peridotite serpentinization, but diffusion coefficients at low temperature are unconstrained and additional experimental investigations are needed. 
    more » « less