Uranophane dissolution and growth in CaCl2-SiO2(aq) test solutions

The dissolution and growth of uranophane [Ca(UO₂)₂(SiO₃OH)₂·5H₂O] have been examined in Ca- and Si-rich test solutions at low temperatures (20.5 ± 2.0 °C) and near-neutral pH (∼6.0). Uranium-bearing experimental solutions undersaturated and supersaturated with uranophane were prepared in matrices of ∼10⁻² M CaCl₂ and ∼10⁻³ M SiO₂(aq). The experimental solutions were reacted with synthetic uranophane and analyzed periodically over 10 weeks. Interpretation of the aqueous solution data permitted extraction of a solubility constant for the uranophane dissolution reaction and standard state Gibbs free energy of formation for uranophane ( kJ mol⁻¹).     

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.     

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:

     

Redox equilibria of iron and silicate melt structure: Implications for olivine/melt element partitioning

Olivine/melt partitioning of ΣFe, Fe²⁺, Mg²⁺, Ca²⁺, Mn²⁺, Co²⁺, and Ni²⁺ has been determined in the systems CaO-MgO-FeO-Fe₂O₃-SiO₂ (FD) and CaO-MgO-FeO-Fe₂O₃-Al₂O₃-SiO₂ (FDA3) as a function of oxygen fugacity (f_O2) at 0.1 MPa pressure. Total iron oxide content of the starting materials was ∼20 wt%. The f_O2 was to used to control the Fe³⁺/ΣFe (ΣFe: total iron) of the melts. The Fe³⁺/ΣFe and structural roles of Fe²⁺ and Fe³⁺ were determined with ⁵⁷Fe resonant absorption Mössbauer spectroscopy. Changes in melt polymerization, NBO/T, as a function of f_O2 was estimated from the Mössbauer data and existing melt structure information. It varies by ∼100% in melts coexisting with olivine in the FDA3 system and by about 300% in the FD system in the Fe³⁺/ΣFe range of the experiments (0.805-0.092). The partition coefficients ( in olivine/wt% in melt) are systematic functions of f_O2 and, therefore, NBO/T of the melt. There is a -minimum in the FDA3 system at NBO/T-values corresponding to intermediate Fe³⁺/ΣFe (0.34-0.44). In the Al-free system, FD, where the NBO/T values of melts range between ∼1 and ∼2.9, the partition coefficients are positively correlated with NBO/T (decreasing Fe³⁺/ΣFe). These relationships are explained by consideration of solution behavior in the melts governed by Qⁿ-unit distribution and structural changes of the divalent cations in the melts (coordination number, complexing with Fe³⁺, and distortion of the polyhedra).     

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.     

Thermodynamics of multicomponent pyroxenes: I. Formulation of a general model

A model is proposed for the thermodynamic properties of multicomponent pyroxenes in the composition space defined by the end-member component CaMgSi₂O₆ and the exchange components Fe(Mg)_-1, TiAl₂(MgSi₂)_-1, Fe³⁺(Al)_-1, Fe³⁺Al(MgSi)_-1, and Mg(Ca)_-1. It is formulated for the simplifying assumptions that: (1) a molecular mixing type approximation describes changes in the molar configurational entropy associated with the coupled exchange substitutions TiAl₂MgSi₂, Fe³⁺AlMgSi, and Al₂MgSi (and their ferroan equivalents), and (2) Fe²⁺ and Mg²⁺, and Al³⁺ and Fe³⁺ display long-range non-convergent ordering between M2 and octahedral M1 sites, and octahedral M1 and tetrahedral sites, respectively. The molar vibrational Gibbs energy is described by a Taylor expansion of second degree in seven linearly independent composition and ordering variables, which is extended to third degree to account for asymmetry in the mixing of Ca and Mg, and Ca and Fe on the M2 site, and is further modified for the assumption that the standard state properties of Ca end-member components of clinopyroxenes are linearly dependent on the coordination number of Ca²⁺ on the M2 site. The model is shown to be consistent with miscibility gap feaures of pyroxenes in the system CaMgSi₂O₆–CaTiAl₂O₆–CaAl₂SiO₆. In subsequent papers, the model is calibrated for the simplifying assumptions that: (1) all regular-solution-type parameters are constants independent of temperature, (2) Pbca and C2/c end-members have identical heat capacities and coefficients of thermal expansion and compressibility, and (3) the heat capacities and coefficients of thermal expansion and compressibility are zero for all reciprocal reactions relating Pbca and pigeonite or high-calcium pyroxene C2/c endmember components.     

Hydroxamate siderophore-promoted reactions between iron(II) and nitroaromatic groundwater contaminants

Recent studies show that ferrous iron (Fe^II), which is often abundant in anaerobic soil and groundwater, is capable of abiotically reducing many subsurface contaminants. However, studies also demonstrate that Fe^II redox reactivity in geochemical systems is heavily dependent upon metal speciation. This contribution examines the influence of hydroxamate ligands, including the trihydroxamate siderophore desferrioxamine B (DFOB), on Fe^II reactions with nitroaromatic groundwater contaminants (NACs). Experimental results demonstrate that ring-substituted NACs are reduced to the corresponding aniline products in aqueous solutions containing Fe^II complexes with DFOB and two monohydroxamate ligands (acetohydroxamic acid and salicylhydroxamic acid). Reaction rates are heavily dependent upon solution conditions and the identities of both the Fe^II-complexing hydroxamate ligand and the target NAC. Trends in the observed pseudo-first-order rate constants for reduction of 4-chloronitrobenzene (k_obs, s⁻¹) are quantitatively linked to the formation of Fe^II species with standard one-electron reduction potentials, (Fe^III/Fe^II), below −0.3 V. Linear free energy relationships correlate reaction rates with the (Fe^III/Fe^II) values of different electron-donating Fe^II complexes and with the apparent one-electron reduction potentials of different electron-accepting NACs, (ArNO₂). Experiments describing a redox auto-decomposition mechanism for Fe^II-DFOB complexes that occurs at neutral pH and has implications for the stability of hydroxamate siderophores in anaerobic environments are also presented. Results from this study indicate that hydroxamates and other Fe^III-stabilizing organic ligands can form highly redox-active Fe^II complexes that may contribute to the natural attenuation and remediation of subsurface contaminants.     

The activity coefficients of Fe(III) hydroxide complexes in NaCl and NaClO4 solutions

The osmotic coefficients of FeCl₃ at 25 °C from 0.15 to 1.7 m [Rumyantsev et al., Z. Phys. Chem., 218, 1089-1127, 2004] have been used to determine the Pitzer parameters (β⁽⁰⁾, β⁽¹⁾ and C^?) for FeCl₃. Since the differences in the Pitzer coefficients of rare earths in NaCl and NaClO₄ are small, the values of Fe(ClO₄)₃ have been estimated using the differences between La(ClO₄)₃ and LaCl₃. The Pitzer coefficients for FeCl₃ combined with enthalpy and heat capacity data for the rare earths can be used to estimate the activity coefficients of Fe³⁺ in NaCl over a wide range of temperatures (0 to 50 °C) and ionic strength (0 to 6 m).The activity coefficients of Fe³⁺ in NaCl and NaClO₄ solutions have been used to determine the activity coefficients of Fe(OH)²⁺ in these solutions from the measured first hydrolysis constants of Fe³⁺ [Byrne et al., Mar. Chem., 97, 34-48, 2005]. The activity coefficients of , Fe(OH)₃ and from 0 to 50 °C have also been determined from the solubility measurements of Fe(III) in NaCl solutions [Liu and Millero, Geochim. Cosmochim Acta, 63, 3487-3497, 1999]. These activity coefficients have been fitted to the Pitzer equations. These results can be used to estimate the speciation of Fe(III) with OH⁻ in natural waters with high concentrations of NaCl from 0 to 50 °C.     

Calcium isotope record of Phanerozoic oceans: Implications for chemical evolution of seawater and its causative mechanisms

A total of 280 brachiopods of Ordovician to Cretaceous age, complemented by published data from belemnites and planktonic foraminifera, are used to reconstruct the evolution of calcium isotope composition of seawater (δ^44/40Ca_SW) over the Phanerozoic. The compiled δ^44/40Ca_SW record shows a general increase from ∼1.3‰ (NIST SRM 915a) at the beginning of the Ordovician to ∼2‰ at present. Superimposed on this trend is a major long-term positive excursion from the Early Carboniferous to Early Permian as well as several short-term, mostly negative, oscillations.A numerical model of the global cycles of calcium, carbon, magnesium and strontium was used to estimate whether the recorded δ^44/40Ca_SW variations can be explained by varying magnitudes of input and output fluxes of calcium to the oceans. The model uses the record of marine ⁸⁷Sr/⁸⁶Sr ratios as proxy for seafloor spreading rates, a record of oceanic Mg/Ca ratios to estimate rates of dolomite formation, and reconstructed atmospheric CO₂, discharge and erosion rates to estimate continental weathering fluxes.The model results indicate that varying magnitudes of the calcium input and output fluxes cannot explain the observed δ^44/40Ca_SW trends, suggesting that the isotope signatures of these fluxes must also have changed. As a possible mechanism we suggest variable isotope fractionation in the sedimentary output flux controlled by the dominant mineralogy in marine carbonate deposits, i.e. the oscillating ‘calcite-aragonite seas’. The ultimate control of the calcium isotope budget of the Phanerozoic oceans appears to have been tectonic processes, specifically variable rates of oceanic crust production that modulated the hydrothermal calcium flux and the oceanic Mg/Ca ratio, which in turn controlled the dominant mineralogy of marine carbonates, hence the δ^44/40Ca_SW. As to the causes of the short-term oscillations recorded in the secular δ^44/40Ca_SW trend, we tentatively propose that these are related to variable rates of dolomite formation and/or to changing chemical composition of the riverine flux, in particular and ratios, induced by variable proportions of silicate vs. carbonate weathering rates on the continents.     

Secular variation in the major-ion chemistry of seawater: Evidence from fluid inclusions in Cretaceous halites

The major-ion (Mg²⁺, Ca²⁺, Na⁺, K⁺, , and Cl⁻) chemistry of Cretaceous seawater was determined from analyses of seawater-derived brines preserved as fluid inclusions in marine halites. Fluid inclusions in primary halite from three evaporite deposits were analyzed by the environmental scanning electron microscopy (ESEM) X-ray energy dispersive spectrometry (EDS) technique: the Early Cretaceous (Aptian, 121.0-112.2 Ma) of the Sergipe basin, Brazil and the Congo basin, Republic of the Congo, and the Early to Late Cretaceous (Albian to Cenomanian, 112.2-93.5 Ma) of the Khorat Plateau, Laos, and Thailand. The fluid inclusions in halite indicate that Cretaceous seawater was enriched several fold in Ca²⁺, depleted in , Na⁺, and Mg²⁺, and had lower Na⁺/Cl⁻, Mg²⁺/Ca²⁺, and Mg²⁺/K⁺ ratios compared to modern seawater. Elevated Ca²⁺ concentrations, with Ca²⁺ >  at the point of gypsum saturation, allowed Cretaceous seawater to evolve into Mg²⁺-Ca²⁺-Na⁺-K⁺-Cl⁻ brines lacking measurable .The major-ion composition of Cretaceous seawater was modeled from fluid inclusion chemistries for the Aptian and the Albian-Cenomanian. Aptian seawater was extreme in its Ca²⁺ enrichment, more than three times higher than present day seawater, with a Mg²⁺/Ca²⁺ ratio of 1.1-1.3. Younger, Albian-Cenomanian seawater had lower Ca²⁺ concentrations, and a higher Mg²⁺/Ca²⁺ ratio of 1.2-1.7. Cretaceous (Aptian) seawater has the lowest Mg²⁺/Ca²⁺ ratios so far documented in Phanerozoic seawater from fluid inclusions in halite, and within the range chemically favorable for precipitation of low-Mg calcite ooids and cements. Results from halite fluid inclusions, together with Mg²⁺/Ca²⁺ ratios measured from echinoderm and rudist calcite, all indicate that Early Cretaceous seawater (Hauterivian, Barremian, Aptian, and Albian) had lower Mg²⁺/Ca²⁺ ratios than Late Cretaceous seawater (Coniacian, Santonian, and Campanian). Low Aptian-Albian Mg²⁺/Ca²⁺ seawater ratios coincide with negative excursions of ⁸⁷Sr/⁸⁶Sr ratios and δ³⁴S_SO4, and peak Cretaceous ocean crust production rates, all of which suggests a link between seawater chemistry and midocean ridge hydrothermal brine flux.     

Calcium and magnesium isotope systematics in rivers draining the Himalaya-Tibetan-Plateau region: Lithological or fractionation control?

In rivers draining the Himalaya-Tibetan-Plateau region, the ²⁶Mg/²⁴Mg ratio has a range of 2‰ and the ⁴⁴Ca/⁴²Ca ratio has a range of 0.6‰. The average δ²⁶Mg values of tributaries from each of the main lithotectonic units (Tethyan Sedimentary Series (TSS), High Himalayan Crystalline Series (HHCS) and Lesser Himalayan Series (LHS)) are within 2 standard deviation analytical uncertainty (0.14‰). The consistency of average riverine δ²⁶Mg values is in contrast to the main rock types (limestone, dolostone and silicate) which range in their average δ²⁶Mg values by more than 2‰. Tributaries draining the dolostones of the LHS differ in their values compared to tributaries from the TSS and HHCS. The chemistry of these river waters is strongly influenced by dolostone (solute Mg/Ca close to unity) and both δ²⁶Mg (−1.31‰) and (0.64‰) values are within analytical uncertainty of the LHS dolostone. These are the most elevated values in rivers and rock reported so far demonstrating that both riverine and bedrock values may show greater variability than previously thought.Although rivers draining TSS limestone have the lowest values at −1.41 and 0.42‰, respectively, both are offset to higher values compared to bedrock TSS limestone. The average δ²⁶Mg value of rivers draining mainly silicate rock of the HHCS is −1.25‰, lower by 0.63‰ than the average silicate rock. These differences are consistent with a fractionation of δ²⁶Mg values during silicate weathering. Given that the proportion of Mg exported from the Himalaya as solute Mg is small, the difference in ²⁶Mg/²⁴Mg ratios between silicate rock and solute Mg reflects the ²⁶Mg/²⁴Mg isotopic fractionation factor () between silicate and dissolved Mg during incongruent silicate weathering. The value of of 0.99937 implies that in the TSS, solute Mg is primarily derived from silicate weathering, whereas the source of Ca is overwhelmingly derived from carbonate weathering. The average value in HHCS rivers is within uncertainty of silicate rock at 0.39‰. The widespread hot springs of the High Himalaya have an average δ²⁶Mg value of −0.46‰ and an average value of 0.5‰, distinct from riverine values for δ²⁶Mg but similar to riverine values. Although rivers draining each major rock type have and δ²⁶Mg values in part inherited from bedrock, there is no correlation with proxies for carbonate or silicate lithology such as Na/Ca ratios, suggesting that Ca and Mg are in part recycled. However, in spite of the vast contrast in vegetation density between the arid Tibetan Plateau and the tropical Lesser Himalaya, the isotopic fractionation factor for Ca and Mg between solute and rocks are not systematically different suggesting that vegetation may only recycle a small amount of Ca and Mg in these catchments.The discrepancy between solute and solid Ca and Mg isotope ratios in these rivers from diverse weathering environments highlight our lack of understanding concerning the origin and subsequent path of Ca and Mg, bound as minerals in rock, and released as cations in rivers. The fractionation of Ca and Mg isotope ratios may prove useful for tracing mechanisms of chemical alteration. Ca isotope ratios of solute riverine Ca show a greater variability than previously acknowledged. The variability of Ca isotope ratios in modern rivers will need to be better quantified and accounted for in future models of global Ca cycling, if past variations in oceanic Ca isotope ratios are to be of use in constraining the past carbon cycle.     

Carbon isotope fractionation in synthetic magnesian calcite

Mg-calcite was precipitated at 25 °C in closed system, free-drift experiments, from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The carbon stable isotope composition of bulk solid and solution were analyzed from subsamples collected during time course experiments of 24 h duration. Considering only the Mg-content and δ¹³C values for the bulk solid, the carbon isotope fractionation factor for the Mg-calcite-HCO_3(aq)⁻ system (as ) increased with average mol percentage of Mg (X_Mg) in the solid at a rate of (0.024 ± 0.011) per mol% MgCO₃. Extrapolation of this relationship to the pure calcite end member yields a value of 0.82 ± 0.09, which is similar to published values for the calcite-HCO_3(aq)⁻ system. Although did not vary for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol m⁻² h⁻¹, it was not possible to hold Mg-content of the solid constant, so kinetic effect on 10³ ln α could not be evaluated from these experiments.     

Determination of Sr and Ba partition coefficients between apatite and water from 5°C to 60°C: a potential new thermometer for aquatic paleoenvironments

Apparent partition coefficients of Sr and Ba between calcium phosphate and water were measured experimentally for temperature ranging from 5°C to 60°C. Calcium phosphates were precipitated from an aqueous mixture of Na₂HPO₄ · 2H₂O (10⁻² M) and CaCl₂ · 2H₂O (10⁻² M). Spiked solutions of Sr or Ba were introduced into the CaCl₂ · 2H₂O solution at Sr/Ca and Ba/Ca ratios of 0.1. The experiment consisted in sampling the liquid and solid phases after 1, 6, 48, and 96 h of interaction. The amorphous calcium phosphate (ACP) precipitated early in the experiment was progressively replaced by hydroxylapatite (HAP), except at 5°C where brushite (di-calcium phosphate di-hydrate or DCPD) was formed. We observed that the crystallinity of the solid phase increased with time for a given temperature and increased with temperature for a given time of reaction. With the exception of the experiment at 5°C, yield (R%) and apparent partition coefficients (K_a-w^Sr/Ca and K_a-w^Ba/Ca) both decreased with increasing reaction time. After 96 h, R%, K_a-w^Sr/Ca and K_a-w^Ba/Ca were observed to be constant, suggesting that the solid phases were at steady-state with respect to the aqueous solutions. The thermodependence of Sr and Ba partitioning between apatite and water at low temperature could therefore be calculated:

     

Synthesis, characterization and thermochemistry of a Pb-jarosite

The enthalpy of formation from the elements of a well-characterized synthetic Pb-jarosite sample corresponding to the chemical formula (H₃O)_0.74Pb_0.13Fe_2.92(SO₄)₂(OH)_5.76(H₂O)_0.24 was measured by high temperature oxide melt solution calorimetry. This value ( = −3695.9 ± 9.7 kJ/mol) is the first direct measurement of the heat of formation for a lead-containing jarosite. Comparison to the thermochemical properties of hydronium jarosite and plumbojarosite end-members strongly suggests the existence of a negative enthalpy of mixing possibly related to the nonrandom distribution of Pb²⁺ ions within the jarosite structure. Based on these considerations, the following thermodynamic data are proposed as the recommended values for the enthalpy of formation from the elements of the ideal stoichiometric plumbojarosite Pb_0.5Fe₃(SO₄)₂(OH)₆:  = −3118.1 ± 4.6 kJ/mol,  = −3603.6 ± 4.6 kJ/mol and S° = 376.6 ± 4.5 J/(mol K). These data should prove helpful for the calculation of phase diagrams of the Pb-Fe-SO₄-H₂O system and for estimating the solubility product of pure plumbojarosite. For illustration, the evolution of the estimated solubility product of ideal plumbojarosite as a function of temperature in the range 5-45 °C was computed (Log(K_sp) ranging from −24.3 to −26.2). An Eh-pH diagram is also presented.     

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.     

Kinetics and mechanism of polythionate oxidation to sulfate at low pH by O2 and Fe

Polythionates (S_xO₆²⁻) are important in redox transformations involving many sulfur compounds. Here we investigate the oxidation kinetics and mechanisms of trithionate and tetrathionate oxidation between pH 0.4 and pH 2 in the presence of Fe³⁺ and/or oxygen. In these solutions, Fe³⁺ plus oxygen oxidizes tetrathionate and trithionate at least an order of magnitude faster than oxygen alone. Kinetic measurements, coupled with density functional calculations, suggest that the rate-limiting step for tetrathionate oxidation involves Fe³⁺ attachment, followed by electron density shifts that result in formation of a sulfite radical and S₃O₃⁰ derivatives. The overall reaction orders for trithionate and tetrathionate are fractional due to rearrangement reactions and side reactions between reactants and intermediate products. The pseudo-first order rate coefficients for tetrathionate range from 10⁻¹¹ s⁻¹ at 25°C to 10⁻⁸ s⁻¹ at 70°C, compared to 2 × 10⁻⁷ s⁻¹ at 35 °C for trithionate. The apparent activation energy (E_A) for tetrathionate oxidation at pH 1.5 is 104.5 ± 4.13 kJ/mol. A rate law at pH 1.5 and 70°C between 0.5 and 5 millimolar [Fe³⁺] is of the form:

     

The effect of calcium on aqueous uranium(VI) speciation and adsorption to ferrihydrite and quartz

Recent studies of uranium(VI) geochemistry have focused on the potentially important role of the aqueous species, CaUO₂(CO₃)₃²⁻ and Ca₂UO₂(CO₃)₃⁰(aq), on inhibition of microbial reduction and uranium(VI) aqueous speciation in contaminated groundwater. However, to our knowledge, there have been no direct studies of the effects of these species on U(VI) adsorption by mineral phases. The sorption of U(VI) on quartz and ferrihydrite was investigated in NaNO₃ solutions equilibrated with either ambient air (430 ppm CO₂) or 2% CO₂ in the presence of 0, 1.8, or 8.9 mM Ca²⁺. Under conditions where the Ca₂UO₂(CO₃)₃⁰(aq) species predominates U(VI) aqueous speciation, the presence of Ca in solution lowered U(VI) adsorption on quartz from 77% in the absence of Ca to 42% and 10% at Ca concentrations of 1.8 and 8.9 mM, respectively. U(VI) adsorption to ferrihydrite decreased from 83% in the absence of Ca to 57% in the presence of 1.8 mM Ca. Surface complexation model predictions that included the formation constant for aqueous Ca₂UO₂(CO₃)₃⁰(aq) accurately simulated the effect of Ca²⁺ on U(VI) sorption onto quartz and ferrihydrite within the thermodynamic uncertainty of the stability constant value. This study confirms that Ca²⁺ can have a significant impact on the aqueous speciation of U(VI), and consequently, on the sorption and mobility of U(VI) in aquifers.     

Seawater nutrient and carbonate ion concentrations recorded as P/Ca, Ba/Ca, and U/Ca in the deep-sea coral Desmophyllum dianthus

As paleoceanographic archives, deep sea coral skeletons offer the potential for high temporal resolution and precise absolute dating, but have not been fully investigated for geochemical reconstructions of past ocean conditions. Here we assess the utility of skeletal P/Ca, Ba/Ca and U/Ca in the deep sea coral D. dianthus as proxies of dissolved phosphate (remineralized at shallow depths), dissolved barium (trace element with silicate-type distribution) and carbonate ion concentrations, respectively. Measurements of these proxies in globally distributed D. dianthus specimens show clear dependence on corresponding seawater properties. Linear regression fits of mean coral Element/Ca ratios against seawater properties yield the equations: P/Ca_coral (μmol/mol) = (0.6 ± 0.1) P/Ca_sw(μmol/mol) - (23 ± 18), R² = 0.6, n = 16 and Ba/Ca_coral(μmol/mol) = (1.4 ± 0.3) Ba/Ca_sw(μmol/mol) + (0 ± 2), R² = 0.6, n = 17; no significant relationship is observed between the residuals of each regression and seawater temperature, salinity, pressure, pH or carbonate ion concentrations, suggesting that these variables were not significant secondary dependencies of these proxies. Four D. dianthus specimens growing at locations with Ω_arag ? 0.6 displayed markedly depleted P/Ca compared to the regression based on the remaining samples, a behavior attributed to an undersaturation effect. These corals were excluded from the calibration. Coral U/Ca correlates with seawater carbonate ion: U/Ca_coral(μmol/mol) = (−0.016 ± 0.003) (μmol/kg) + (3.2 ± 0.3), R² = 0.6, n = 17. The residuals of the U/Ca calibration are not significantly related to temperature, salinity, or pressure. Scatter about the linear calibration lines is attributed to imperfect spatial-temporal matches between the selected globally distributed specimens and available water column chemical data, and potentially to unresolved additional effects. The uncertainties of these initial proxy calibration regressions predict that dissolved phosphate could be reconstructed to ±0.4 μmol/kg (for 1.3-1.9 μmol/kg phosphate), and dissolved Ba to ±19 nmol/kg (for 41-82 nmol/kg Ba_sw). Carbonate ion concentration derived from U/Ca has an uncertainty of ±31μmol/kg (for ). The effect of microskeletal variability on P/Ca, Ba/Ca, and U/Ca was also assessed, with emphasis on centers of calcification, Fe-Mn phases, and external contaminants. Overall, the results show strong potential for reconstructing aspects of water mass mixing and biogeochemical processes in intermediate and deep waters using fossil deep-sea corals.     

Molecular and structural investigation of yukonite (synthetic & natural) and its relation to arseniosiderite

This paper describes the first synthesis method of yukonite, its thorough molecular and structural analysis along with natural specimens originating from Tagish Lake (Canada) and Grotta della Monaca (Italy) for comparison, and its structural relation to arseniosiderite. The synthetic and natural yukonites were found to have a range of composition according to the general formula Ca₂Fe_3-5(AsO₄)₃(OH)_4-10·xH₂O where x = 2-11. The synthetic yukonite was found to be equivalent at the atomic, molecular and structural level to the Tagish Lake yukonite. At the molecular level, arseniosiderite, via vibrational spectroscopy, was found to have a H-bonding system as in scorodite and exhibit an extra arsenate mode indicative of groups. Heating experiments along with ATR-IR analysis indicated the presence of structural water and hydroxyl units in arseniosiderite. In yukonite in contrast, a wide diffuse H-bonding environment was observed with only arsenate groups. The presence of both structural water and hydroxyl groups was further verified via ATR-IR spectroscopy. The As K, Fe 2p and Ca 2p XANES spectra of yukonite and arseniosiderite were found to be identical, confirming that the local nature of the As, Fe and Ca atoms in these structures is the same. Diffraction analysis (X-ray and electron) showed that yukonite consist of nano-crystalline domains while in the case of arseniosiderite micro-size single crystal domains exist.