Oxygen isotopic fractionation during inorganic calcite precipitation ― Effects of temperature, precipitation rate and pH

Stable oxygen isotopic fractionation during inorganic calcite precipitation was experimentally studied by spontaneous precipitation at various pH (8.3 < pH < 10.5), precipitation rates (1.8 < log R < 4.4 μmol m^− 2 h^− 1) and temperatures (5, 25, and 40 °C) using the CO₂ diffusion technique.The results show that the apparent stable oxygen isotopic fractionation factor between calcite and water (α_{calcite–water}) is affected by temperature, the pH of the solution, and the precipitation rate of calcite. Isotopic equilibrium is not maintained during spontaneous precipitation from the solution. Under isotopic non-equilibrium conditions, at a constant temperature and precipitation rate, apparent 1000lnα_{calcite–water} decreases with increasing pH of the solution. If the temperature and pH are held constant, apparent 1000lnα_{calcite–water} values decrease with elevated precipitation rates of calcite. At pH = 8.3, oxygen isotopic fractionation between inorganically precipitated calcite and water as a function of the precipitation rate (R) can be described by the expressions

at 5, 25, and 40 °C, respectively.The impact of precipitation rate on 1000lnα_{calcite–water} value in our experiments clearly indicates a kinetic effect on oxygen isotopic fractionation during calcite precipitation from aqueous solution, even if calcite precipitated slowly from aqueous solution at the given temperature range. Our results support Coplen's work [Coplen T. B. (2007) Calibration of the calcite–water oxygen isotope geothermometer at Devils Hole, Nevada, a natural laboratory. Geochim. Cosmochim. Acta 71, 3948–3957], which indicates that the equilibrium oxygen isotopic fractionation factor might be greater than the commonly accepted value.     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: I. Sr incorporation

Partitioning of strontium during spontaneous calcite formation was experimentally studied using an advanced CO₂-diffusion technique. Results at different precipitation rates and T = 5, 25, and 40 °C show that at constant temperature Sr incorporation into calcite is controlled by the precipitation rate (R in μmol/m²/h) according to the individual expressions

     

The influence of temperature and seawater composition on calcite crystal growth mechanisms and kinetics: Implications for Mg incorporation in calcite lattice

The composition of carbonate minerals formed in past and present oceans is assumed to be significantly controlled by temperature and seawater composition. To determine if and how temperature is kinetically responsible for the amount of Mg incorporated in calcite, we quantified the influence of temperature and specific dissolved components on the complex mechanism of calcite precipitation in seawater. A kinetic study was carried out in artificial seawater and NaCl-CaCl₂ solutions, each having a total ionic strength of 0.7 M. The constant addition technique was used to maintain [Ca²⁺] at 10.5 mmol kg⁻¹ while [] was varied to isolate the role of this variable on the precipitation rate of calcite.Our results show that the overall reaction of calcite precipitation in both seawater and NaCl-CaCl₂ solutions is dominated by the following reaction:

     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

Sorption of yttrium and rare earth elements by amorphous ferric hydroxide: Influence of solution complexation with carbonate

The influence of solution complexation on the sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide was investigated at 25 °C over a range of pH (4.0-7.1) and carbonate concentrations . Distribution coefficients, defined as , where [MS_i]_T is the total concentration of sorbed YREE, M_T is the total YREE concentration in solution, and [S_i] is the concentration of amorphous ferric hydroxide, initially increased in magnitude with increasing carbonate concentration, and then decreased. The initial increase of is due to sorption of YREE carbonate complexes , in addition to sorption of free YREE ions (M³⁺). The subsequent decrease of , which is more extensive for the heavy REEs, is due to the increasing intensity of YREE solution complexation by carbonate ions. The competition for YREEs between solution complexation and surface complexation was modeled via the equation:

     

Experimental study of gold-hydrosulphide complexing in aqueous solutions at 350-500°C, 500 and 1000 bars using mineral buffers

The solubility of gold was measured in KCl solutions (0.001-0.1 m) at near-neutral to weakly acidic pH in the presence of the K-feldspar-muscovite-quartz, andalusite-muscovite-quartz, and pyrite-pyrrhotite-magnetite buffers at temperatures 350 to 500°C and pressures 0.5 and 1 kbar. These mineral buffers were used to simultaneously constrain pH, f(S₂), and f(H₂). The experiments were performed using a CORETEST flexible Ti-cell rocking hydrothermal reactor enabling solution sampling at experimental conditions. Measured log m(Au) (mol/kg H₂O) ranges from −7.5 at weakly acid pH to −5.9 in near-neutral solutions, and increases slightly with temperature. Gold solubility in weakly basic and near-neutral solutions decreases with decreasing pH at all temperatures, which implies that Au(HS)₂ is the dominant Au species in solution. In more acidic solutions, solubility is independent of pH. Comparison of the experimentally measured solubilities with literature values for Au hydrolysis constants demonstrates that at 350°C dominates Au aqueous speciation at the weakly acidic pH and f(S₂)/f(H₂) conditions imposed by the pyrite-pyrrhotite-magnetite buffer. In contrast, at temperatures >400°C becomes less important and predominates in weakly acid solutions. Solubility data collected in this study were used to calculate the following equilibrium reaction constants:

     

The ionic strength dependence of lead (II) carbonate complexation in perchlorate media

Lead speciation in many aqueous geochemical systems is dominated by carbonate complexation. However, direct observations of Pb²⁺ complexation by carbonate ions are few in number. This work represents the first investigation of the equilibrium over a range of ionic strength. Through spectrophotometric observations of formation at 25 °C in NaHCO₃-NaClO₄ solutions, formation constants of the form were determined between 0.001 and 5.0 molal ionic strength. Formation constant results were well represented by the equation:

     

A potentiometric study of the stability of aqueous yttrium-acetate complexes from 25 to 175 °C and 1-1000 bar

The stability of yttrium-acetate (Y-Ac) complexes in aqueous solution was determined potentiometrically at temperatures 25-175 °C (at P_s) and pressures 1-1000 bar (at 25 and 75 °C). Measurements were performed using glass H⁺-selective electrodes in potentiometric cells with a liquid junction. The species YAc²⁺ and were found to dominate yttrium aqueous speciation in experimental solutions at 25-100 °C (log [Ac⁻] < −1.5, pH < 5.2), whereas at 125, 150 and 175 °C introduction of into the Y-Ac speciation model was necessary. The overall stability constants β_n were determined for the reaction

     

Comparative carbonate and hydroxide complexation of cations in seawater

Comparative concentrations of carbonate and hydroxide complexes in natural solutions can be expressed in terms of reactions with bicarbonate that have no explicit pH dependence (). Stability constants for this reaction with n = 1 were determined using conventional formation constant data expressed in terms of hydroxide and carbonate. Available data indicate that stability constants appropriate to seawater at 25 °C expressed in the form are on the order of 10^4.2 for a wide range of cations (M^z+) with z = +1, +2 and +3. Φ₁ is sufficiently large that species appear to substantially dominate MOH^z−1 species in seawater. Evaluations of comparative stepwise carbonate and hydroxide stability constant behavior leading to the formation of n = 2 and n = 3 complexes suggest that carbonate complexes generally dominate hydroxide complexes in seawater, even for cations whose inorganic speciation schemes in seawater are currently presumed to be strongly dominated by hydrolyzed forms (). Calculated stability constants, and , indicate that the importance of carbonate complexation is sufficiently large that carbonate and hydroxide complexes would be generally comparable even if calculated Φ₂ and Φ₃ values are overestimated by two or more orders of magnitude. Inclusion of mixed ligand species in carbonate-hydroxide speciation models allows cation complexation intensities (M_T/[M^z+]) to be expressed in the following form:

     

Free energy of formation for green rust sodium sulphate (NaFe6Fe3(OH)18(SO4)2(s))

In a recent study, sulphate-bearing green rust (GR_SO4) was shown to incorporate Na⁺ in its structure (NaFe^II₆Fe^III₃(OH)₁₈(SO₄)_2(s); GR_Na,SO4). The compound was synthesised by aerial oxidation of Fe(OH)_2(s) in the presence of NaOH. This paper reports on its free energy of formation .Freshly synthesised GR_Na,SO4 was titrated with 0.5 M H₂SO₄ in an inert atmosphere at 25 °C, producing dissolved Fe²⁺ and magnetite or goethite. Solution concentrations, PHREEQC and the MINTEQ database were used to calculate reaction constants for the reactions:

     

Molybdic acid ionisation under hydrothermal conditions to 300 °C

This UV spectrophotometric study was aimed at providing precise, experimentally derived thermodynamic data for the ionisation of molybdic acid (H₂MoO₄) from 30 to 300 °C and at equilibrium saturated vapour pressures. The determination of the equilibrium constants and associated thermodynamic parameters were facilitated by spectrophotometric measurements using a specially designed high temperature optical Ti-Pd flow-through cell with silica glass windows.The following van’t Hoff isochore equations describe the temperature dependence of the first and second ionisation constants of molybdic acid up to 300 °C:

     

Solubility measurements of synthetic neodymium monazite as a function of temperature at 2 kbars, and aqueous neodymium speciation in equilibrium with monazite

The solubility of synthetic NdPO₄ monazite end-member was experimentally determined from 300 up to 800 °C, at 2000 bars in pure water, and in aqueous chloride or phosphate solutions. Both the classical weight-loss method and a new method based on isotope dilution coupled with thermal ionization mass spectrometer were used. In the range of temperature studied monazite showed a prograde solubility from 10^−5.4 m at 300 °C up to 10^−2.57 m at 800 °C. Experiments in H₂O-H₃PO₄-NaCl-HCl solutions suggested Nd(OH)₃⁰ was the major species that was formed at high temperature and pressure. The equilibrium constants (log K) for the reaction:

     

Chemical controls on the solubility, speciation and mobility of lanthanum at near surface conditions: A geochemical modeling study

Recent experimental determinations of the solubility products of common rare earth minerals such as monazite and xenotime and stability constants for chloride, sulfate, carbonate and hydroxide complexes provide a basis to model quantitatively the solubility, and therefore the mobility, of rare earth elements (REE) at near surface conditions. Data on the mobility of REE and stabilities of REE complexes at near-neutral conditions are of importance to safe nuclear waste disposal, and environmental monitoring. The aim of this study is to understand REE speciation and solubility of a given REE in natural environments. In this study, a series of formation constants for La aqueous complexes are recommended by using the specific interaction theory (SIT) for extrapolation to infinite dilution. Then, a thermodynamic model has been employed for calculation of the solubility and speciation of La in soil solutions reacted with the La end-member of mineral monazite (LaPO₄), and other La-bearing solid phases including amorphous lanthanum hydroxide (La(OH)₃, am) and different La carbonates, as a function of various inorganic and organic ligand concentrations. Calculations were carried out at near-neutral pH (pH 5.5–8.5) and 25 °C at atmospheric CO₂ partial pressure. The model takes account of the species: La³⁺, LaCl²⁺, , , , , , , , , La(OH)²⁺, LaOx⁺, , LaAc²⁺ and (where Ox²⁻ = oxalate and Ac⁻ = acetate).The calculations indicate that the La species that dominate at pH 5.5–8.5 in the baseline model soil solution (BMSS) include La³⁺, LaOx⁺, , and in order of increasing importance as pH rises. The solubility of monazite in the BMSS remains less than 3 × 10⁻⁹ M, exhibiting a minimum of 2 × 10⁻¹² M at pH 7.5. The calculations quantitatively demonstrate that the concentrations of La controlled by the solubility of other La-bearing solid phases are many orders of magnitude higher than those controlled by monazite in the pH range from 5.5 to 8.5, suggesting that monazite is likely to be the solubility-controlling phase at this pH range. The calculations also suggest that significant mobility of La (and other REE) is unlikely because high water–rock ratios on the order of at least 10⁴ (mass ratio) are required to move 50% of the La from a soil. An increase in concentration of oxalate by one order of magnitude from that of the baseline model solution results in the dominance of LaOx⁺ at pH 5.5–7.5. Similarly, the increase in concentration of by one order of magnitude makes the dominant species at pH 5.5–7.5. Above pH 7.5, carbonate complexes are important. The increase in oxalate or concentrations by one order of magnitude can enhance the solubility of monazite by a factor of up to about 6 below neutral pH, in comparison with that in the baseline model soil solution. From pH 7.0 to 8.5, the solubility of monazite in the soil solutions with higher concentrations of oxalate or is similar, or almost identical, to that in the BMSS.     

Pb and rare earth element diffusion in xenotime

Diffusion of Pb and the rare earth elements Sm, Dy and Yb have been characterized in synthetic xenotime under dry conditions. The synthetic xenotime was grown via a Na₂CO₃–MoO₃ flux method. The sources of diffusant for the rare earth diffusion experiments were REE phosphate powders, with experiments run using sources containing a single REE. For Pb, the source consisted a mixture of YPO₄ and PbTiO₃. Experiments were performed by placing source and xenotime in Pt capsules, and annealing capsules in 1 atm furnaces for times ranging from 30 min to several weeks, at temperatures from 1000 to 1500 °C. The REE and Pb distributions in the xenotime were profiled by Rutherford Backscattering Spectrometry (RBS).The following Arrhenius relations are obtained for diffusion in xenotime, normal to (101):

Diffusivities among the REE do not differ greatly in xenotime over the investigated temperature range, in contrast to findings for the REE in zircon [Cherniak, D.J., Hanchar, J.M., Watson, E.B., 1997. Rare earth diffusion in zircon. Chem. Geol. 134, 289–301.], where the LREE diffuse more slowly, and with higher activation energies for diffusion, than the heavier rare earths. In zircon, these differences among diffusion of the rare earths are attributed to the relatively large size of the REE with respect to Zr, for which they likely substitute in the zircon lattice. With the systematic increase in ionic radius from the heavy to lighter REE, this size mismatch becomes more pronounced and diffusivities of the LREE are as consequence slower. Although xenotime is isostructural with zircon, the REE are more closely matched in size to Y, so in xenotime this effect appears much smaller and the REE diffuse at similar rates. In addition, the process of diffusion in xenotime likely involves simple REE^+ 3 → Y^+ 3 exchange, without charge compensation as needed for REE^+ 3 → Zr^+ 4 exchange in zircon. This latter factor may also contribute to the large activation energies for diffusion of the REE in zircon (i.e., 691–841 kJ mol^− 1, [Cherniak, D.J., Hanchar, J.M., Watson, E.B., 1997. Rare earth diffusion in zircon. Chem. Geol. 134, 289–301.]), in comparison with those for xenotime.For Pb, the following Arrhenius relation is obtained (also normal to (101)):

These measurements suggest that Pb diffusion in xenotime is quite slow, even slower than Pb diffusion in monazite and zircon, and considerably slower than diffusion of the REE in xenotime. Xenotime may therefore be even more retentive of Pb isotope signatures than either monazite or zircon in cases where Pb isotopes are altered solely by volume diffusion. However, because the activation energy for Pb diffusion in xenotime is lower than those for monazite and zircon, Pb diffusion may be somewhat faster at many temperatures of geologic interest in xenotime than in monazite or zircon.     

Stability constants for the formation of lead chloride complexes as a function of temperature and ionic strength

The stability constants for the formation of lead (Pb²⁺) with chloride

     

A spectrophotometric study of samarium (III) speciation in chloride solutions at elevated temperatures

The speciation of samarium (III) in chloride-bearing solutions was investigated spectrophotometrically at temperatures of 100-250 °C and a pressure of 100 bars. The simple hydrated ion, Sm³⁺, is predominant at ambient temperature, but chloride complexes are the dominant species at elevated temperatures. Cumulative formation constants for samarium chloride species were calculated for the following reactions:

     

Solubility and dissolution rate of silica in acid fluoride solutions

We performed 57 batch reactor experiments in acidic fluoride solutions to measure the dissolution rate of quartz. These rate data along with rate data from published studies were fit using multiple linear regression to produce the following non-unique rate law for quartz

where 10^−5.13 < a_HF < 10^1.60, −0.28 < pH < 7.18, and 298 < T < 373 K. Similarly, 97 amorphous silica dissolution rate data from published studies were fit by multiple linear regression to produce the following non-unique rate law for amorphous silica

where 10^−2.37 < a_HF < 10^1.61, −0.32 < pH < 4.76 and 296 < T < 343 K. Regression of the rates versus other combinations of solution species, e.g.  + H⁺, F⁻ + H⁺, HF + , HF + F⁻, or  + F⁻, produced equally good fits. Any of these rate laws can be interpreted to mean that the rate-determining step for silica dissolution in fluoride solutions involves a coordinated attack of a Lewis acid, on the bridging O atom and a Lewis base on the Si atom. This allows a redistribution of electrons from the Si–O bond to form a O–H group and a Si–FH group.     

H2O diffusion models in rhyolitic melt with new high pressure data

Water diffusion in silicate melts is important for understanding bubble growth in magma, magma degassing and eruption dynamics of volcanos. Previous studies have made significant progress on water diffusion in silicate melts, especially rhyolitic melt. However, the pressure dependence of H₂O diffusion is not constrained satisfactorily. We investigated H₂O diffusion in rhyolitic melt at 0.95–1.9 GPa and 407–1629 °C, and 0.2–5.2 wt.% total water (H₂O_t) content with the diffusion-couple method in a piston-cylinder apparatus. Compared to previous data at 0.1–500 MPa, H₂O diffusivity is smaller at higher pressures, indicating a negative pressure effect. This pressure effect is more pronounced at low temperatures. Assuming H₂O diffusion in rhyolitic melt is controlled by the mobility of molecular H₂O (H₂O_m), the diffusivity of H₂O_m (D_H2Om) at H₂O_t ≤ 7.7 wt.%, 403–1629 °C, and ≤ 1.9 GPa is given by

D_H2Om=D₀exp(aX),