An experimental study of Cobalt (II) complexation in Cl and H2S-bearing hydrothermal solutions

The speciation of cobalt (II) in Cl⁻ and H₂S-bearing solutions was investigated spectrophotometrically at temperatures of 200, 250, and 300 °C and a pressure of 100 bars, and by measuring the solubility of cobaltpentlandite at temperatures of 120-300 °C and variable pressures of H₂S. From the results of these experiments, it is evident that CoHS⁺ and predominate in the solutions except at 150 °C, for which the dominant chloride complex is CoCl₃⁻. The logarithms of the stability constant for CoHS⁺ show moderate variation with temperature, decreasing from 6.24 at 120 °C to 5.84 at 200 °C, and increasing to 6.52 at 300 °C. Formation constants for chloride species increase smoothly with temperature and at 300°C their logarithms reach 8.33 for , 6.44 for CoCl₃⁻, 4.94 to 5.36 for , and 2.42 for CoCl⁺. Calculations based on the composition of a model hydrothermal fluid (Ksp-Mu-Qz, KCl = 0.25 m, NaCl = 0.75 m, ΣS = 0.3 m) suggest that at temperatures ?200 °C, cobalt occurs dominantly as CoHS⁺, whereas at higher temperatures the dominant species is .     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

Arsenic, iron and sulfur co-diagenesis in lake sediments

Profiles of porewater pH and dissolved As, Fe, Mn, sulfate, total sulfide (ΣS^−II), total zero-valent sulfur (ΣS⁰), organic carbon and major ion concentrations, as well as those of solid As, acid-volatile sulfide (AVS), total S, Fe, Mn, Al, organic C, ²¹⁰Pb and ¹³⁷Cs were determined in the sediment of four lakes spanning a range of redox and geochemical conditions. An inverse modeling approach, based on a one-dimensional transport-reaction equation assuming steady-state, was applied to the porewater As profiles and used to constrain the net rates of reactions involving As (). The model defines depth intervals where As is either released to (positive ) or removed from (negative ) the porewaters.At two of the sites, whose bottom water were oxygenated at sampling time, a production zone ( = 12 × 10⁻¹⁸ mol cm⁻³ s⁻¹-71 × 10⁻¹⁸ mol cm⁻³ s⁻¹) is inferred a few cm below the sediment-water interface, coincident with sharp porewater As and Fe peaks that indicate an intense coupled recycling of As and Fe. This process is confirmed by solid As and Fe maxima just below the sediment surface. In these two lakes a zone of As consumption ( = −5 × 10⁻¹⁸ mol cm⁻³ s⁻¹ to −53 × 10⁻¹⁸ mol cm⁻³ s⁻¹), attributed to the slow adsorption of As to authigenic Fe oxyhydroxides, occurs just above the production zone. A second-order rate constant of 0.12 ± 0.03 cm³ mol⁻¹ s⁻¹ is estimated for this adsorption reaction.Such features in the porewater and solid profiles were absent from the two other lakes that develop a seasonally anoxic hypolimnion. Thermodynamic calculations indicate that the porewaters of the four lakes, when sulfidic (i.e., ΣS^−II ? 0.1 μM), were undersaturated with respect to all known solid As sulfides; the calculation also predicts the presence of As^V oxythioanions in the sulfidic waters, as suggested by a recent study. In the sulfidic waters, the removal of As ( = −1 × 10⁻¹⁸ mol cm⁻³ s⁻¹ to −23 × 10⁻¹⁸ mol cm⁻³ s⁻¹) consistently occurred when saturation, with respect to FeS_(s), was reached and when As^V oxythioanions were predicted to be significant components of total dissolved As. This finding has potential implications for As transport in other anoxic waters and should be tested in a wider variety of natural environments.     

Grain boundary partitioning of Ar and He

An experimental procedure has been developed that permits measurement of the partitioning of Ar and He between crystal interiors and the intergranular medium (ITM) that surrounds them in synthetic melt-free polycrystalline diopside aggregates. ³⁷Ar and ⁴He are introduced into the samples via neutron irradiation. As samples are crystallized under sub-solidus conditions from a pure diopside glass in a piston cylinder apparatus, noble gases diffusively equilibrate between the evolving crystal and intergranular reservoirs. After equilibration, ITM Ar and He is distinguished from that incorporated within the crystals by means of step heating analysis. An apparent equilibrium state (i.e., constant partitioning) is reached after about 20 h in the 1450 °C experiments. Data for longer durations show a systematic trend of decreasing ITM Ar (and He) with decreasing grain boundary (GB) interfacial area as would be predicted for partitioning controlled by the network of planar grain boundaries (as opposed to ITM gases distributed in discrete micro-bubbles or melt). These data yield values of GB-area-normalized partitioning, , with units of (Ar/m³ of solid)/(Ar/m² of GB) of 6.8 × 10³-2.4 × 10⁴ m⁻¹. Combined petrographic microscope, SEM, and limited TEM observation showed no evidence that a residual glass phase or grain boundary micro-bubbles dominated the ITM, though they may represent minor components. If a nominal GB thickness (δ) is assumed, and if the density of crystals and the grain boundaries are assumed equal, then a true grain boundary partition coefficient , may be determined. For reasonable values of δ, is at least an order of magnitude lower than the Ar partition coefficient between diopside and melt. Helium partitioning data provide a less robust constraint with between 4 × 10³ and 4 × 10⁴ cm⁻¹, similar to the Ar partitioning data. These data suggest that an ITM consisting of nominally melt free, bubble free, tight grain boundaries can constitute a significant but not infinite reservoir, and therefore bulk transport pathway, for noble gases in fine grained portions of the crust and mantle where aqueous or melt fluids are non-wetting and of very low abundance (i.e., <0.1% fluid). Heterogeneities in grain size within dry equilibrated systems will correspond to significant differences in bulk rock noble gas content.     

Microbial ammonia oxidation and enhanced nitrogen cycling in the Endeavour hydrothermal plume

Ammonium was injected from the subseafloor hydrothermal system at the Endeavour Segment, Juan de Fuca Ridge, into the deep-sea water column resulting in an -rich (?177 nM) neutrally buoyant hydrothermal plume. This was quickly removed by both autotrophic ammonia oxidation and assimilation. The former accounted for at least 93% of total net removal, with its maximum rate in the neutrally buoyant plume (?53 nM d⁻¹) up to 10-fold that in background deep water. Ammonia oxidation in this plume potentially added 26-130 mg into the deep-sea water column. This oxidation process was heavily influenced by the presence of organic-rich particles, with which ammonia-oxidizing bacteria (AOB) were often associated (40-68%). AOB contributed up to 10.8% of the total microbial communities within the plume, and might constitute a novel lineage of β-proteobacterial AOB based on 16S rRNA and amoA phylogenetic analyses. Meanwhile, assimilation rates were also substantially enhanced within the neutrally buoyant plume (?26.4 nM d⁻¹) and accounted for at least 47% of total net removal rates. The combined oxidation and assimilation rates always exceeded total net removal rates, suggesting active in situregeneration rates of at least an order of magnitude greater than the particulate nitrogen flux from the euphotic zone. Ammonia oxidation is responsible for turnover of 0.7-13 days and is probably the predominant in situ organic carbon production process (0.6-13 mg C m⁻² d⁻¹) at early stages of Endeavour neutrally buoyant plumes.     

The problem of inherited Ar* in dating impact glass by the Ar/Ar method: Evidence from the Tswaing impact crater (South Africa)

The Tswaing meteorite impact crater is a 1.13 km diameter structure located in the 2.05 Ga Nebo granite of the Bushveld Complex. The impact age had previously been determined by fission track dating to 220 ± 104 ka. ⁴⁰Ar/³⁹Ar step-heating and total fusion experiments performed on single- and multi-grain impact glass aliquots gave apparent ages ranging from 1.0 ± 0.3 Ma to 204 ± 6 Ma. These “ages” indicate that the radiogenic Ar derived from the target rocks has not been completely degassed as a result of the impact process, despite fusion of the target material. Results of step-heating experiments imply that the trapped within the glass is located in two distinct reservoirs thought to be the glass matrix and fluid/vapor inclusions (or un-melted residual clasts). Calculations assuming an age of 0.2 ± 0.1 Ma for Tswaing (fission track data) reveal that the amount of inherited ⁴⁰Ar*() relative to the pre-impact concentration varies from 0.015% to 4.15%. The spread defined by likely reflects the various quench rates experienced by the glass, most certainly due to the pre-impact position of the sample relative to the center of the crater. We compare the influence of on the apparent ⁴⁰Ar/³⁹Ar age determination of five impact structures. Our calculations show that the main characteristic controlling the age offset (for a given proportion of ) is the age difference between the impact and the target rocks (i.e., the ⁴⁰Ar* concentration in the target rock). The buffer effect for a given crater structure can be predicted knowing the age of the basement and having a rough estimation of the age of the crater structure itself. The occurrence of is likely influenced by (1) the degree of polymerization (i.e., silicate structure complexity) of the target rock and presumably related to the diffusivity of Ar in the melt and glass, (2) the Ar partial pressure at the grain boundary, (3) the quantity of energy involved in the impact, and (4) the porosity of the target rocks. For glass that inevitably suffers inherited and/or excess ⁴⁰Ar*, the use of the inverse isochron technique can be appropriate but should be applied with careful statistical treatment.     

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.     

The effect of sulfur on the partitioning of Ni and other first-row transition elements between olivine and silicate melt

The effect of sulfur dissolved as sulfide (S²⁻) in silicate melts on the activity coefficients of NiO and some other oxides of divalent cations (Ca, Cr, Mn, Fe and Co) has been determined from olivine/melt partitioning experiments at 1400 °C in six melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), and in derivatives of these compositions at 1370 °C, obtained from the six CMAS compositions by substituting Fe for Mg (FeCMAS). Amounts of S²⁻ were varied from zero to sulfide saturation, reaching 4100 μg g⁻¹ S in the most sulfur-rich silicate melt. The sulfide solubilities compare reasonably well with those predicted from the parameterization of the sulfide capacity of silicate melts at 1400 °C of O’Neill and Mavrogenes (2002), although in detail systematic deviations indicate that a more sophisticated model may improve the prediction of sulfide capacities.The results show a barely discernible effect of S²⁻ in the silicate melt on Fe, Co and Ni partition coefficients, and also surprisingly, a tiny but resolvable effect on Ca partitioning, but no detectable effect on Cr, Mn or some other lithophile incompatible elements (Sc, Ti, V, Y, Zr and Hf). Decreasing Mg# of olivine (reflecting increasing FeO in the system) has a significant influence on the partitioning of several of the divalent cations, particularly Ca and Ni. We find a remarkably systematic correlation between and the ionic radius of M²⁺, where M = Ca, Cr, Mn, Fe, Co or Ni, which is attributable to a simple relationship between size mismatch and excess free energies of mixing in Mg-rich olivine solid solutions.Neither the effect of S²⁻ nor of Mg#^ol is large enough by an order of magnitude to account for the reported variations of obtained from electron microprobe analyses of olivine/glass pairs from mid-ocean ridge basalts (MORBs). Comparing these MORB glass analyses with the Ni-MgO systematics of MORB from other studies in the literature, which were obtained using a variety of analytical techniques, shows that these electron microprobe analyses are anomalous. We suggest that the reported variation of with S content in MORB is an analytical artifact.Mass balance of melt and olivine compositions with the starting compositions shows that dissolved S²⁻ depresses the olivine liquidus of haplobasaltic silicate melts by 5.8 × 10⁻³ (±1.3 × 10⁻³) K per μg g⁻¹ of S²⁻, which is negligible in most contexts. We also present data for the partitioning of some incompatible trace elements (Sc, Ti, Y, Zr and Hf) between olivine and melt. The data for Sc and Y confirm previous results showing that and decrease with increasing SiO₂ content of the melt. Values of average 0.01 with most falling in the range 0.005-0.015. Zr and Hf are considerably more incompatible than Ti in olivine, with and about 10⁻³. The ratio / is well constrained at 0.611 ± 0.016.     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

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:

     

Coupled mobilization of dissolved organic matter and metals (Cu and Zn) in soil columns

Dissolved organic carbon (DOC) is a key component involved in metal displacement in soils. In this study, we investigated the concentration profiles of soil-borne DOC, Cu and Zn at various irrigation rates with synthetic rain water under quasi steady-state conditions, using repacked soil columns with a metal-polluted topsoil and two unpolluted subsoils. Soil solution was collected using suction cups installed at centimeter intervals over depth. In the topsoil the concentrations of DOC, dissolved metals (Zn and Cu), major cations (Ca²⁺ and Mg²⁺) and anions ( and ) increased with depth. In the subsoil, the Cu and Zn concentrations dropped to background levels within 2 cm. All compounds were much faster mobilized in the first 4 cm than in the rest of the topsoil. DOC and Cu concentrations were higher at higher flow rates for a given depth, whereas the concentrations of the other ions decreased with increasing flow rate. The decomposition of soil organic matter resulted in the formation of DOC, , and and was the main driver of the system. Regression analysis indicated that Cu mobilization was governed by DOC, whereas Zn mobilization was primarily determined by Ca and to a lesser extent by DOC. Labile Zn and Cu²⁺ concentrations were well predicted by the NICA-Donnan model. The results highlight the value of high-resolution in-situ measurements of DOC and metal mobilization in soil profiles.     

Dissociation constants and speciation in aqueous Li2SO4 and K2SO4 from measurements of electrical conductance to 673 K and 29 MPa

The electrical conductivities of aqueous solutions of Li₂SO₄ and K₂SO₄ have been measured at 523-673 K at 20-29 MPa in dilute solutions for molalities up to 2 × 10⁻² mol kg⁻¹. These conductivities have been fitted to the conductance equation of Turq, Blum, Bernard, and Kunz with a consensus mixing rule and mean spherical approximation activity coefficients. In the temperature interval 523-653 K, where the dielectric constant, ε, is greater than 14, the electrical conductance data can be fitted by a solution model which includes ion association to form , , and , where M is Li or K. The adjustable parameters of this model are the first and second dissociation constants of the M₂SO₄. For the 673 K and 300 kg m⁻³ state point where the Coulomb interactions are the strongest (dielectric constant, ε = 5), models with more extensive association give good fits to the data. In the case of the Li₂SO₄ model, including the multi-ion associate, , gave an extremely good fit to the conductance data.