Characterization of elemental release during microbe-granite interactions at T = 28 °C

This study used batch reactors to characterize the mechanisms and rates of elemental release (Al, Ca, K, Mg, Na, F, Fe, P, Sr, and Si) during interaction of a single bacterial species (Burkholderia fungorum) with granite at T = 28 °C for 35 days. The objective was to evaluate how actively metabolizing heterotrophic bacteria might influence granite weathering on the continents. We supplied glucose as a C source, either NH₄ or NO₃ as N sources, and either dissolved PO₄ or trace apatite in granite as P sources. Cell growth occurred under all experimental conditions. However, solution pH decreased from ∼7 to 4 in NH₄-bearing reactors, whereas pH remained near-neutral in NO₃-bearing reactors. Measurements of dissolved CO₂ and gluconate together with mass-balances for cell growth suggest that pH lowering in NH₄-bearing reactors resulted from gluconic acid release and H⁺ extrusion during NH₄ uptake. In NO₃-bearing reactors, B. fungormum likely produced gluconic acid and consumed H⁺ simultaneously during NO₃ utilization.Over the entire 35-day period, NH₄-bearing biotic reactors yielded the highest release rates for all elements considered. However, chemical analyses of biomass show that bacteria scavenged Na, P, and Sr during growth. Abiotic control reactors followed different reaction paths and experienced much lower elemental release rates compared to biotic reactors. Because release rates inversely correlate with pH, we conclude that proton-promoted dissolution was the dominant reaction mechanism. Solute speciation modeling indicates that formation of Al-F and Fe-F complexes in biotic reactors may have enhanced mineral solubilities and release rates by lowering Al and Fe activities. Mass-balances further reveal that Ca-bearing trace phases (calcite, fluorite, and fluorapatite) provided most of the dissolved Ca, whereas more abundant phases (plagioclase) contributed negligible amounts. Our findings imply that during the incipient stages of granite weathering, heterotrophic bacteria utilizing glucose and NH₄ only moderately elevate silicate weathering reactions that consume atmospheric CO₂. However, by enhancing the dissolution of non-silicate, Ca-bearing trace minerals, they could contribute to high Ca/Na ratios commonly observed in granitic watersheds.     

Characterization of elemental release during microbe-basalt interactions at T = 28 °C

This study used batch reactors to characterize the rates and mechanisms of elemental release during the interaction of a single bacterial species (Burkholderia fungorum) with Columbia River Flood Basalt at T = 28 °C for 36 days. We primarily examined the release of Ca, Mg, P, Si, and Sr under a variety of biotic and abiotic conditions with the aim of evaluating how actively metabolizing bacteria might influence basalt weathering on the continents. Four days after inoculating P-limited reactors (those lacking P in the growth medium), the concentration of viable planktonic cells increased from ∼10⁴ to 10⁸ CFU (Colony Forming Units)/mL, pH decreased from ∼7 to 4, and glucose decreased from ∼1200 to 0 μmol/L. Mass-balance and acid-base equilibria calculations suggest that the lowered pH resulted from either respired CO₂, organic acids released during biomass synthesis, or H⁺ extrusion during uptake. Between days 4 and 36, cell numbers remained constant at ∼10⁸ CFU/mL and pH increased to ∼5. Purely abiotic control reactors as well as control reactors containing inert cells (∼10⁸ CFU/mL) showed constant glucose concentrations, thus confirming the absence of biological activity in these experiments. The pH of all control reactors remained near-neutral, except for one experiment where the pH was initially adjusted to 4 but rapidly rose to 7 within 2 days. Over the entire 36 day period, P-limited reactors containing viable bacteria yielded the highest Ca, Mg, Si, and Sr release rates. Release rates inversely correlate with pH, indicating that proton-promoted dissolution was the dominant reaction mechanism. Both biotic and abiotic P-limited reactors displayed low P concentrations. Chemical analyses of bacteria collected at the end of the experiments, combined with mass-balances between the biological and fluid phases, demonstrate that the absence of dissolved P in the biotic reactors resulted from microbial P uptake. The only P source in the basalt is a small amount of apatite (∼1.2%), which occurs as needles within feldspar grains and glass. We therefore conclude that B. fungorum utilized apatite as a P source for biomass synthesis, which stimulated elemental release from coexisting mineral phases via pH lowering. The results of this study suggest that actively metabolizing bacteria have the potential to influence elemental release from basalt in continental settings.     

Kaolinite dissolution and precipitation kinetics at 22 °C and pH 4

Dissolution and precipitation rates of low defect Georgia kaolinite (KGa-1b) as a function of Gibbs free energy of reaction (or reaction affinity) were measured at 22 °C and pH 4 in continuously stirred flowthrough reactors. Steady state dissolution experiments showed slightly incongruent dissolution, with a Si/Al ratio of about 1.12 that is attributed to the re-adsorption of Al on to the kaolinite surface. No inhibition of the kaolinite dissolution rate was apparent when dissolved aluminum was varied from 0 and 60 μM. The relationship between dissolution rates and the reaction affinity can be described well by a Transition State Theory (TST) rate formulation with a Temkin coefficient of 2

     

An experimental study of calcite dissolution rates at acidic conditions and 25 °C in the presence of NaPO3 and MgCl2

Dissolution rates of single calcite crystals were determined from sample weight loss using free-drift rotating disk techniques. Experiments were performed at 25 °C in aqueous HCl solutions over the bulk solution pH range −1 to 3 and in the presence of trace concentrations of aqueous NaPO₃ and MgCl₂. These salts were chosen for this study because aqueous magnesium and phosphate are known to strongly inhibit calcite dissolution at neutral to basic pH. Reactive solutions were undersaturated with respect to possible secondary phases. Neither an inhibition or enhancement of calcite dissolution rates was observed in the presence of aqueous MgCl₂ at pH 1 and 3. The presence of trace quantities of NaPO₃, which dissociates in solution to Na⁺ and H₂PO₄⁻, decreased the overall calcite dissolution rate at pH≤2. This contrasting behavior could be attributed to the different adsorption behavior of these dissolved species. As calcite surfaces are positively charged in acidic solutions, aqueous Mg²⁺ may not adsorb, whereas aqueous phosphate, present as either the anion H₂PO₄⁻ or the neutral species H₃PO₄⁰, readily adsorbs on calcite surfaces leading to significant dissolution inhibition.     

An experimental study of the solubility of molybdenum in H2O and KCl-H2O solutions from 500 °C to 800 °C, and 150 to 300 MPa

The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H₂O and pure H₂O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H₂O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H₂O aqueous solutions. Similarly, in the pure H₂O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS₂).     

The effect of CO2 on the speciation of RbBr in solution at temperatures to 579 °C and pressures to 0.26 GPa

Carbon dioxide- and salt-bearing solutions are common in granulite, ore-forming and magmatic environments. The presence of CO₂ affects mineral solubilities, fluid miscibility, and viscosity and wetting properties, and is expected to affect salt speciation. EXAFS measurements of RbBr-H₂O-CO₂ fluids contained in corundum-osed synthetic fluid inclusions (SFLINCs) have been used to investigate the effect of CO₂ on salt speciation at temperatures to 579 °C and pressures to around 0.26 GPa.Forward modelling indicates that solute dehydration is difficult to distinguish from up to around 40% of Rb-Br ion-pairing, so results refer to the total number of nearest neighbours, which are likely to be mostly O present in waters of hydration, but may also include Br, if ion pairing is present. Additionally, results relate to the number of well-ordered neighbours in the first shell, because nearest neighbours with a high degree of disorder may be present but contribute minimally to the EXAFS signal. Analysis of the EXAFS results at the Rb edge for the CO₂-free solution is consistent with previous work and shows that the number of nearest neighbours for Rb in CO₂-free solutions decreases from 6 ± 0.6 to 1.4 ± 0.1 as temperature increases from 20 to 534 °C. The decrease is accompanied by a decrease in Rb-x bondlengths of 0.05 Å, where x is the first shell scatterer. Results for the CO₂-bearing solution are different to those for the CO₂-free solution. The number of nearest neighbours is 16 and 22% less than for the CO₂-bearing solution at 312 and 445 °C respectively. Changes in the numbers of nearest neighbours correlate well with calculated changes in the bulk solution dielectric constant; CO₂-bearing and CO₂-free solutions lie on the same trend, which suggests that it may be possible to calculate the number of nearest neighbours from dielectric constant. Rb-x bondlengths for the CO₂-bearing solution are statistically indistinguishable to those for the CO₂-free inclusions. Results for Br are worse quality than for Rb so EXAFS analysis could not be completed, however XANES spectra for CO₂-free and CO₂-bearing solutions are consistent with solute dehydration similar to that recorded by the Rb spectra. The conclusions of this study provide support for the notion that CO₂ has a fundamental effect on the mechanics of solubility, and that these effects should be incorporated into conceptual and quantitative thermodynamic models.     

Dissolution kinetics of fosteritic olivine at 90-150 °C including effects of the presence of CO2

The steady state dissolution rate of San Carlos olivine [Mg_1.82Fe_0.18 SiO₄] in dilute aqueous solutions was measured at 90, 120, and 150 °C and pH ranging from 2 to 12.5. Dissolution experiments were performed in a stirred flow-through reactor, under either a nitrogen or carbon dioxide atmosphere at pressures between 15 and 180 bar. Low pH values were achieved either by adding HCl to the solution or by pressurising the reactor with CO₂, whereas high pH values were achieved by adding LiOH. Dissolution was stoichiometric for almost all experiments except for a brief start-up period. At all three temperatures, the dissolution rate decreases with increasing pH at acidic to neutral conditions with a slope of close to 0.5; by regressing all data for 2 ? pH ? 8.5 and 90 °C ? T ? 150 °C together, the following correlation for the dissolution rate in CO₂-free solutions is obtained:

     

A simple predictive model of quartz solubility in water-salt-CO2 systems at temperatures up to 1000 °C and pressures up to 1000 MPa

Knowledge of the solubility of quartz over a broad spectrum of aqueous fluid compositions and T-P conditions is essential to our understanding of water-rock interaction in the Earth’s crust. We propose an equation to compute the molality of aqueous silica, m_SiO2(aq), mol·(kg H₂O)⁻¹, in equilibrium with quartz and water-salt-CO₂ fluids, as follows:

     

In situ Atomic Force Microscopy study of dissolution of the barite (0 0 1) surface in water at 30 °C

The dissolution behavior of the barite (0 0 1) surface in pure water at 30 °C was investigated using in situ Atomic Force Microscopy (AFM), to better understand the dissolution mechanism and the microtopographical changes that occur during the dissolution, such as steps and etch pits. The dissolution of the barite (0 0 1) surface started with the slow retreat of steps, after which, about 60 min later, the <hk0> steps of one unit cell layer or multi-layers became two-step fronts (fast “f” and slow “s” steps) with a half-unit cell layer showing different retreat rates. The “f” step had a fast retreat rate (≈(14 ± 1) × 10⁻² nm/s) and tended to have a jagged step edge, whereas the “s” step (≈(1.8 ± 0.1) × 10⁻² nm/s) had a relatively straight front. The formation of the “f” steps led to the formation of a new one-layer step, where the front of the “s” step was overtaken by that of the immediate underlying “f” step. The “f” steps also led to the decrease of the <hk0> steps and the increase in the percentage of stable steps parallel to the [0 1 0] direction during the dissolution.Etch pits, which could be observed after about 90 min, were of three types: triangular etch pits with a depth of a half-unit cell, shallow etch pits, and deep etch pits. The triangular etch pits were bounded by the step edges parallel to [0 1 0], [1 2 0], and [] and had opposite orientations in the upper half and lower half layers. Shallow etch pits that had a depth of two or more half-unit cell layers had any two consecutive pits pointing in the opposite direction of each other. The triangular etch pit appeared to grow by simultaneously removal of a row of ions parallel to each direction from the three step edges. At first, deep etch pits were elongated in the [0 1 0] direction with a curved outline and then gradually developed to an angular form bounded by the {1 0 0}, {3 1 0}, and (0 0 1) faces. The retreat rate of the (0 0 1) face was much slower than those of the {1 0 0} and {3 1 0} and tended to separate into two rates ((0.13 ± 0.01) × 10⁻² nm/s for the deep etch pits derived from a screw dislocation and (0.07 ± 0.01) × 10⁻² nm/s for those from other line defects).The changes in the dissolution rate of a barite (0 0 1) surface during the dissolution were estimated using the retreat rates and densities of the various steps as well as the growth rates, density, and areas of the lateral faces of the deep etch pits that were obtained from this AFM analysis. Our results showed that the dissolution rate of the barite (0 0 1) surface gradually increased and approached the bulk dissolution rate because of the change in the main factor determining the dissolution rate from the density of the steps to the growth and the density of the deep etch pits on the surface.     

Effect of pH and organic ligands on the kinetics of smectite dissolution at 25 °C

Forward dissolution rates of Na-Montmorillonite (Wyoming) SWy-2 smectite (Ca_0.06Na_0.56)[Al_3.08Fe(III)_0.38Mg_0.54] [Si_7.93 Al_0.07]O₂₀(OH)₄ were measured at 25 °C in a mixed-flow reactor equipped with interior dialysis compartment (6-8 kDa membrane) as a function of pH (1-12), dissolved carbonate (0.5-10 mM), phosphate (10⁻⁵ to 0.03 M), and nine organic ligands (acetate, oxalate, citrate, EDTA, alginate, glucuronic acid, 3,4-dihydroxybenzoic acid, gluconate, and glucosamine) in the concentration range from 10⁻⁵ to 0.03 M. In organic-free solutions, the Si-based rates decrease with increasing pH at 1 ? pH ? 8 with a slope close to −0.2. At 9 ? pH ? 12, the Si-based rates increase with a slope of ∼0.3. In contrast, non-stoichiometric Mg release weakly depends on pH at 1 ? pH ? 12 and decreases with increasing pH. The empirical expression describing Si-release rates [R, mol/cm²/s] obtained in the present study at 25 °C, I = 0.01 M is given by

     

Kinetics and mechanism of natural fluorapatite dissolution at 25 °C and pH from 3 to 12

The dissolution rates of natural fluorapatite (FAP), Ca₁₀(PO₄)₆F₂, were measured at 25 °C in mixed-flow reactors as a function of pH from 3.0 to 11.7, and aqueous calcium, phosphorus, and fluoride concentration. After an initial preferential Ca and/or F release, stoichiometric Ca, P, and F release was observed. Measured FAP dissolution rates decrease with increasing pH at 3 ? pH ? 7, FAP dissolution rates are pH independent at 7 ? pH ? 10, and FAP dissolution rates again decrease with increasing pH at pH ? 10. Measured FAP dissolution rates are independent of aqueous Ca, P, and F concentration at pH ≈ 3 and pH ≈ 10.Apatite dissolution appears to be initiated by the relatively rapid removal from the near surface of F and the Ca located in the M1 sites, via proton for Ca exchange reactions. Dissolution rates are controlled by the destruction of this F and Ca depleted surface layer. The destruction of this layer is facilitated by the adsorption/penetration of protons into the surface at acidic conditions, and by surface hydration at neutral and basic conditions. Taking into account these two parallel mechanisms, measured fluorapatite forward dissolution rates can be accurately described using

     

Experimental study of the effect of pH on the kinetics of montmorillonite dissolution at 25 °C

The effect of pH on the kinetics of smectite (K-montmorillonite) dissolution was investigated at 25 °C in batch and stirred flow-through reactors over the pH range of 1-13.5, in KNO₃ solutions. Dissolution rates were obtained based on the release of Si and Al at steady-state under far from equilibrium conditions. Dissolution was non-stoichiometric between pH 5 and 10, due to adsorption/precipitation of Al. Dissolution rates computed from batch and flow-through experiments were consistent, irrespective of the Si and Al concentrations. Sample pre-treatment and the interlayer cation do not affect the steady-state dissolution rate or stoichiometry of cation release. The rate dependence on pH can be described by:

     

Experimental investigation of the solubility of albite and jadeite in H2O, with paragonite + quartz at 500 and 600 °C, and 1-2.25 GPa

The solubilities of the assemblages albite + paragonite + quartz and jadeite + paragonite + quartz in H₂O were determined at 500 and 600 °C, 1.0-2.25 GPa, using hydrothermal piston-cylinder methods. The three minerals are isobarically and isothermally invariant in the presence of H₂O, so fluid composition is uniquely determined at each pressure and temperature. A phase-bracketing approach was used to achieve accurate solubility determinations. Albite + quartz and jadeite + quartz dissolve incongruently in H₂O, yielding residual paragonite which could not be retrieved and weighed. Solution composition fixed by the three-mineral assemblage at a given pressure and temperature was therefore bracketed by adding NaSi₃O_6.5 glass in successive experiments, until no paragonite was observed in run products. Solubilities derived from experiments bounding the appearance of paragonite thus constrain the equilibrium fluid composition. Results indicate that, at a given pressure, Na, Al, and Si concentrations are higher at 600 °C than at 500 °C. At both 500 and 600 °C, solubilities of all three elements increase with pressure in the albite stability field, to a maximum at the jadeite-albite-quartz equilibrium. In the jadeite stability field, element concentrations decline with continued pressure increase. At the solubility maximum, Na, Al, and Si concentrations are, respectively, 0.16, 0.05, and 0.48 molal at 500 °C, and 0.45, 0.27, and 1.56 molal at 600 °C. Bulk solubilities are 3.3 and 10.3 wt% oxides, respectively. Observed element concentrations are everywhere greater than those predicted from extrapolated thermodynamic data for simple ions, monomers, ion pairs, and the silica dimer. The measurements therefore require the presence of additional, polymerized Na-Al-Si-bearing species in the solutions. The excess solubility is >50% at all conditions, indicating that polymeric structures are the predominant solutes in the P-T region studied. The solubility patterns likely arise from combination of the large solid volume change associated with the albite-jadeite-quartz equilibrium and the rise in Na-Al-Si polymerization with approach to the hydrothermal melting curves of albite + quartz and jadeite + quartz. Our results indicate that polymerization of Na-Al-Si solutes is a fundamental aspect of fluid-rock interaction at high pressure. In addition, the data suggest that high-pressure metamorphic isograds can impose unexpected controls on metasomatic mass transfer, that significant metasomatic mass transfer prior to melting should be considered in migmatitic terranes, and that polymeric complexes may be an important transport agent in subduction zones.     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.     

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

Solubility of grossular, Ca3Al2Si3O12, in H2O-NaCl solutions at 800 °C and 10 kbar, and the stability of garnet in the system CaSiO3-Al2O3-H2O-NaCl

The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H₂O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca₃Al₂Si₃O₁₂ with increasing NaCl mole fraction (X_NaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H₂O value over the investigated range, indicating strong solute interaction with NaCl. The Ca₃Al₂Si₃O₁₂ mole fraction versus NaCl mole fraction curve has a broad plateau between X_NaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H₂O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO₃-Al₂O₃-H₂O-NaCl. At high Al₂O₃/CaSiO₃ bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H₂O, but disappears once grossular is nucleated. At X_NaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO₃: Al₂O₃ molality increases from 0.0013 in initially pure H₂O to near 0.15 at X_NaCl = 0.4 in CaSiO₃-saturated solutions, a >100-fold enhancement. In contrast, addition of Al₂O₃ to wollastonite-saturated NaCl solutions increases CaSiO₃ molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al₂O₃ solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.     

Dissolution of illite in saline-acidic solutions at 25 °C

The dissolution rate of illite, a common clay mineral in Australian soils, was studied in saline-acidic solutions under far from equilibrium conditions. The clay fraction of Na-saturated Silver Hill illite (K_1.38Na_0.05)(Al_2.87Mg_0.46Fe³⁺_0.39Fe²⁺_0.28Ti_0.07)[Si_7.02Al_0.98]O₂₀(OH)₄ was used for this study. The dissolution rates were measured using flow-through reactors at 25 ± 1 °C, solution pH range of 1.0-4.25 (H₂SO₄) and at two ionic strengths (0.01 and 0.25 M) maintained using NaCl solution. Illite dissolution rates were calculated from the steady state release rates of Al and Si. The dissolution stoichiometry was determined from Al/Si, K/Si, Mg/Si and Fe/Si ratios. The release rates of cations were highly incongruent during the initial stage of experiments, with a preferential release of Al and K over Si in majority of the experiments. An Al/Si ratio >1 was observed at pH 2 and 3 while a ratio close to the stoichiometric composition was observed at pH 1 and 4 at the higher ionic strength. A relatively higher K⁺ release rate was observed at I = 0.25 in 2-4 pH range than at I = 0.01, possibly due to ion exchange reaction between Na⁺ from the solution and K⁺ from interlayer sites of illite. The steady state release rates of K, Fe and Mg were higher than Si over the entire pH range investigated in the study. From the point of view of the dominant structural cations (Si and Al), stoichiometric dissolution of illite occurred at pH 1-4 in the higher ionic strength experiments and at pH ?3 for the lower ionic strength experiments. The experiment at pH 4.25 and at the lower ionic strength exhibited lower R_Al (dissolution rate calculated from steady state Al release) than R_Si (dissolution rate calculated from steady state Si release), possibly due to the adsorption of dissolved Al as the output solutions were undersaturated with respect to gibbsite. The dissolution of illite appears to proceed with the removal of interlayer K followed by the dissolution of octahedral cations (Fe, Mg and Al), the dissolution of Si is the limiting step in the illite dissolution process. A dissolution rate law showing the dependence of illite dissolution rate on proton concentration in the acid-sulfate solutions was derived from the steady state dissolution rates and can be used in predicting the impact of illite dissolution in saline acid-sulfate environments. The fractional reaction orders of 0.32 (I = 0.25) and 0.36 (I = 0.01) obtained in the study for illite dissolution are similar to the values reported for smectite. The dissolution rate of illite is mainly controlled by solution pH and no effect of ionic strength was observed on the dissolution rates.     

Dissolution kinetics of calcite at 50-70 °C: An atomic force microscopic study under near-equilibrium conditions

Direct measurements of calcite faces were performed using in situ atomic force microscopy (AFM) to reveal the dissolution processes as a function of solution saturation state and temperature. Time-sequential AFM images demonstrated that step velocities at constant temperature increased with increasing undersaturation. The anisotropy of obtuse and acute step velocities appeared to become more significant as solutions approached equilibrium and temperature increased. At saturation state Ω > 0.02, a curvilinear boundary was formed at the intersection of two acute steps and the initially rhombohedral etch pit exhibited a nearly triangular shape. This suggests that the and steps may not belong to the calcite-aqueous solution equilibrium system. Further increase in the saturation state (Ω ? 0.3) led to a lack of etch pit formation and dissolution primarily occurred at existing steps, in accordance with Teng (2004). Analysis of step kinetics at different temperatures yielded activation energies of 25 ± 6 kJ/mol and 14 ± 13 kJ/mol for obtuse and acute steps, respectively. The inconsistencies in etch pit morphology, step anisotropy, and step activation energies from the present study with those of studies far-from-equilibrium can be explained by increased influence of the backward reaction, or growth, near-equilibrium. We propose that the backward reaction occurs preferentially at the acute-acute kink sites. The kinetics and effective activation energies of near-equilibrium calcite dissolution presented in this work provide accurate experimental data under likely CO₂ sequestration conditions, and thus are crucial to the development of robust geochemical models that predict the long-term performance of mineral-trapped CO₂.