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 new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

Ammonium loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

Unsteady diagenetic processes and sulfur biogeochemistry in tropical deltaic muds: Implications for oceanic isotope cycles and the sedimentary record

Sedimentary S cycling is usually conceptualized and interpreted within the context of steadily accreting (1-D) transport-reaction regimes. Unsteady processes, however, are common in many sedimentary systems and can result in dramatically different S reaction balances and diagenetic products than steady conditions. Globally important common examples include tropical deltaic topset and inner shelf muds such as those extending from the Amazon River ∼1600 km along the Guianas coast of South America. These deposits are characterized by episodic reworking of the surface seabed over vertical depths of ∼0.1-3 m. Reworked surface sediments act as unsteady, suboxic batch reactors, unconformably overlying relict anoxic, often methanic deposits, and have diagenetic properties largely decoupled from net accumulation of sediment. Despite well-oxygenated water and an abundant reactive organic matter supply, physical disturbance inhibits macrofauna, and benthic communities are dominated by microbial biomass across immense areas. In the surficial suboxic layer, molecular biological analyses, tracer experiments, sediment C/S/Fe compositions, and δ³⁴S, δ¹⁸O of pore water indicate close coupling of anaerobic C, S, and Fe cycles. δ¹⁸O- can increase by 2-3‰ during anaerobic recycling without net change in δ³⁴S-, demonstrating reduction coupled to complete anaerobic reoxidation to and a δ¹⁸O- reduction + reoxidation fractionation factor?12‰ (summed magnitudes). S reoxidation must be coupled to Fe-oxide reduction, contributing to high dissolved Fe²⁺ (∼1 mM) and Fe mobilization-export. The reworking of Amazon-Guianas shelf muds alone may isotopically alter δ¹⁸O- equivalent in mass to?25% of the annual riverine delivery of to the global ocean. Unsteady conditions result in preservation of unusually heavy δ³⁴S isotopic compositions of residual Cr reducible S, ranging from 0‰ to >30‰ in physically reworked deposits. In contrast, bioturbated facies adjacent to physically reworked regions accumulate isotopically light S (δ³⁴S to −20‰) in otherwise similar decomposition regimes. The isotopic patterns of both physically and biologically reworked regions can be simulated with simple diagenetic models. Heavy S isotopic signatures are largely a consequence of unsteady diffusion and progressive anaerobic burndown into underlying deposits, whereas isotopically depleted bioturbated deposits predominantly reflect biogenic diffusive scaling and isotopic distillation/diffusive pumping associated with reoxidation in burrow walls immediately adjacent to reduced zones. The S isotopic transition from unsteady physically controlled regions of the Amazon delta moving laterally into bioturbated facies mimics the transition of S isotopic patterns temporally in the geologic record during the rise of bioturbation. No special role for S disproportionation is required to explain these differences. The potential role of unsteady, suboxic diagenesis and dynamic reworking of sediments has been largely ignored in models of the evolution of surficial elemental cycling and interpretations of the geologic record.     

New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions

Several studies have shown that there is a strong relationship between the distribution of crenarchaeotal isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) and sea surface temperature (SST). Based on this, a ratio of certain GDGTs, called TEX₈₆ (TetraEther indeX of tetraethers consisting of 86 carbon atoms), was developed as a SST proxy. In this study, we determined the distribution of crenarchaeotal isoprenoid GDGTs in 116 core-top sediments mostly from (sub)polar oceans and combined these data with previously published core-top data. Using this extended global core-top dataset (n = 426), we re-assessed the relationship of crenarchaeal isoprenoid GDGTs with SST. We excluded data from the Red Sea from the global core-top dataset to define new indices and calibration models, as the Red Sea with its elevated salinity appeared to behave differently compared to other parts of the oceans. We tested our new indices and calibration models on three different paleo datasets, representing different temperature ranges. Our results indicate that the crenarchaeol regio-isomer plays a more important role for temperature adaptation in (sub)tropical oceans than in (sub)polar oceans, suggesting that there may be differences in membrane adaptation of the resident crenarchaeotal communities at different temperatures. We, therefore, suggest to apply two different calibration models. For the whole calibration temperature range (−3 to 30 °C), a modified version of TEX₈₆ with a logarithmic function which does not include the crenarchaeol regio-isomer, called , is shown to correlate best with SST: (r² = 0.86, n=396, p <0.0001). Application of on sediments from the subpolar Southern Ocean results in realistic absolute SST estimates and a similar SST trend compared to a diatom SST record from the same core. , which is defined as the logarithmic function of TEX₈₆, yields the best correlation with SST, when the data from the (sub)polar oceans are removed: (r² = 00.87, n = 255, p < 0.0001). Furthermore, gives the best correlation for mescosm data with temperatures ranging between 10 and 46 °C. For Quaternary sediments from the tropical Arabian Sea, both and yield similar trends and SST estimates. However, the extrapolation of calibration on a sediment record from a greenhouse world ocean predicts more reliable absolute SST estimates and relative SST changes in agreement with estimates based on the δ¹⁸O of planktonic foraminifera. Based on the comparison of and derived SSTs using the core top data, we recommend applying above 15 °C and below 15 °C. In cases where paleorecords encompass temperatures both below and above 15 °C, we suggest to use .     

An integrated sulfur isotope model for Namibian shelf sediments

In this study the sulfur cycle in the organic-rich mud belt underlying the highly productive upwelling waters of the Namibian shelf is quantified using a 1D reaction-transport model. The model calculates vertical concentration and reaction rate profiles in the top 500 cm of sediment which are compared to a comprehensive dataset which includes carbon, sulfur, nitrogen and iron compounds as well as sulfate reduction (SR) rates and stable sulfur isotopes (³²S, ³⁴S). The sulfur dynamics in the well-mixed surface sediments are strongly influenced by the activity of the large sulfur bacteria Thiomargaritanamibiensis which oxidize sulfide (H₂S) to sulfate () using sea water nitrate () as the terminal electron acceptor. Microbial sulfide oxidation (SOx) is highly efficient, and the model predicts intense cycling between and H₂S driven by coupled SR and SOx at rates exceeding 6.0 mol S m⁻² y⁻¹. More than 96% of the SR is supported by SOx, and only 2-3% of the pool diffuses directly into the sediment from the sea water. A fraction of the produced by Thiomargarita is drawn down deeper into the sediment where it is used to oxidize methane anaerobically, thus preventing high methane concentrations close to the sediment surface. Only a small fraction of total H₂S production is trapped as sedimentary sulfide, mainly pyrite (FeS₂) and organic sulfur (S_org) (∼0.3 wt.%), with a sulfur burial efficiency which is amongst the lowest values reported for marine sediments (<1%). Yet, despite intense SR, FeS₂ and S_org show an isotope composition of ∼5 ‰ at 500 cm depth. These heavy values were simulated by assuming that a fraction of the solid phase sulfur exchanges isotopes with the dissolved sulfide pool. An enrichment in H₂S of ³⁴S towards the sediment-water interface suggests that Thiomargarita preferentially remove H₂³²S from the pore water. A fractionation of 20-30‰ was estimated for SOx (ε_SOx) with the model, along with a maximum fractionation for SR (ε_SR-max) of 100‰. These values are far higher than previous laboratory-based estimates for these processes. Mass balance calculations indicate negligible disproportionation of autochthonous elemental sulfur; an explanation routinely cited in the literature to account for the large fractionations in SR. Instead, the model indicates that repeated multi-stepped sulfide oxidation and intracellular disproportionation by Thiomargarita could, in principle, allow the measured isotope data to be simulated using much lower fractionations for ε_SOx (5‰) and ε_SR (78‰).     

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.     

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.     

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.     

Time-scales for magmatic differentiation at the Snaefellsjökull central volcano, western Iceland: Constraints from U-Th-Pa-Ra disequilibria in post-glacial lavas

We present major and trace element and Sr-Nd-Pb and U-Th-Pa-Ra isotope data for a small sample suite of primarily post-glacial, mildly alkalic volcanic rocks from the Snaefellsjökull central volcano situated off the main rift systems in western Iceland. The volcanic rocks are crystal-poor and range from olivine alkali basalt to trachyte and show tight correlations of major and trace elements that are explained by fractional crystallization involving removal of olivine, plagioclase, clinopyroxene, Fe-Ti oxide and apatite. Sr-Nd-Pb isotopes are practically invariant, consistent with derivation from the same source region. During fractionation from primitive basalt to evolved trachyte, (²³⁰Th/²³²Th), (²³⁰Th/²³⁸U) and (²³¹Pa/²³⁵U) decrease progressively at broadly constant (²³⁸U/²³²Th). A continuous closed-system fractionation model that assumes constant initial (²³⁰Th/²³²Th) in the basaltic precursor melt indicates that hawaiite was derived from olivine basalt by ∼50% fractional crystallization within and trachyte by ∼80% fractionation within . An overrepresentation of evolved basalts and hawaiites with young inferred magma ages in the dataset is consistent with the parental precursor to these magmas intruded into the sub-volcanic magma plumbing system as a consequence of lithospheric rebound caused by deglaciation. Lavas affected by apatite removal have higher (²³¹Pa/²³⁵U) than predicted for simple radioactive decay, suggesting apatite significantly fractionates U from Pa. The proposed fractionation model consistently explains our U-series data assuming and and . If applicable, these D values would indicate that the effect of apatite fractionation must be adequately considered when assessing differentiation time scales using (²³¹Pa/²³⁵U) disequilibria data.     

Experimental study of the incorporation of Li, Sc, Al and other trace elements into olivine

We performed a series of synthesis experiments at 1 atm pressure to investigate the substitution mechanisms of 1+ and 3+ ions into olivine. Forsterite crystals were grown from bulk compositions that contained the element of interest (e.g. Li) and different amounts of additional single trace elements. By working at constant (major element) liquid composition and temperature we eliminated all compositional effects other than those due to the trace elements. Mineral-melt pairs were then analysed to determine the compositional-dependence of the partition coefficient (D), which corresponds to , and where [element] refers to weight concentration of the element in the respective phase.We find that Li forms a stable coupled substitution with Sc and, at above ∼500 ppm Sc in the crystal, Li⁺ and Sc³⁺ ions form an ordered neutral complex ([LiSc]). This complex dissociates at lower trace element concentrations and a second, concentration-independent, mechanism begins to dominate. This second solution mechanism is most likely 2Li⁺ ⇔ Mg²⁺ where one of the Li atoms is in an interstitial position in the crystal lattice. Natural olivines show Li contents slightly greater than Sc (on an atomic basis), indicating that both substitution mechanisms are significant. Unlike Sc, Al does not appear to form a stable complex with Li in the olivine structure.Sodium is highly incompatible in olivine with of ∼0.00015-0.03. Olivine-liquid partitioning of Na⁺ is independent of Sc³⁺ or Al³⁺ concentration. This indicates that the coupled substitution of Na⁺ with any 3+ ions is unlikely. Instead, the relevant substitution mechanism appears to be 2Na⁺ ⇔ Mg²⁺. Although independent of 3+ ion concentration, is inversely correlated with the Li concentration of both melts and crystals, implying that Na competes (unsuccessfully) with Li to replace Mg in the olivine structure.Aluminium is highly incompatible in forsterite . Values of are similar for all phase pairs synthesised from starting materials containing between 10 and 100,000 ppm Al. This suggests that Al is principally incorporated in forsterite by replacing one Mg and one Si atom , where the Al atoms on octahedral (Mg) and tetrahedral (Si) sites are dissociated from one another.The incorporation of gallium into forsterite is influenced by the presence of Li. Where Li concentration in the crystal is much greater than that of Ga (on an atomic basis) we find an excellent correlation between and melt Li content. This relationship indicates that Ga³⁺ and Li⁺ replace 2Mg²⁺ on octahedral sites and that the Ga and Li atoms are, like Sc and Li, strongly associated in the crystal structure.The mechanism by which scandium is incorporated into forsterite is strongly governed by the presence Li. As discussed above, ordered complexes form readily in forsterite in Li-rich experiments. Under Li-absent but Sc-rich conditions (Sc in the crystal >∼500 ppm), is proportional to the concentration of Sc in the melt. This indicates that Sc incorporation is charge-balanced by the formation of magnesium vacancies , and that both species are associated . At lower Sc concentrations (<500 ppm in the crystal), the concentration-dependence of partitioning indicates that the complexes dissociate.Our results demonstrate that partitioning of 1+ and 3+ ions into olivine is complex and involves a range of point defects which yield strongly composition-dependent crystal-melt partition coefficients. Since physical and chemical properties of natural olivine, such as diffusion of ⁶Li and ⁷Li and H₂O solubility, depend on the concentrations of the defects identified in this study, our results provide an important insight into how determining substitution mechanisms can improve our understanding of large-scale mantle processes and properties.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.     

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:

     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Experimental studies of metal-silicate partitioning of Sb: Implications for the terrestrial and lunar mantles

The terrestrial mantle has a well defined Sb depletion of ∼7 ± 1 (Jochum and Hofmann, 1997), and the lunar mantle is depleted relative to the Earth by a factor of ∼50 ± 5 (Wolf and Anders, 1980). Despite these well defined depletions, there are few data upon which to evaluate their origin—whether due to volatility or core formation. We have carried out a series of experiments to isolate several variables such as oxygen fugacity, temperature, pressure, and silicate and metallic melt compositions, on the magnitude of . The activity of Sb in FeNi metal is strongly composition dependent such that solubility of Sb as a function of fO₂ must be corrected for the metal composition. When the correction is applied, Sb solubility is consistent with 3+ valence. Temperature series (at 1.5 GPa) shows that decreases by a factor of 100 over 400 °C, and a pressure series exhibits an additional decrease between ambient pressure (100 MPa) and 13 GPa. A strong dependence upon silicate melt composition is evident from a factor of 100 decrease in between nbo/t values of 0.3 and 1.7. Consideration of all these variables indicates that the small Sb depletion for the Earth’s mantle can be explained by high PT equilibrium partitioning between metal and silicate melt . The relatively large lunar Sb depletion can also be explained by segregation of a small metallic core, at lower pressure conditions where is much higher (2500).     

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.     

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 origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.