Stable boron isotope fractionation between dissolved B(OH)3 and B(OH)4

The stable boron isotope ratio (¹¹B/¹⁰B) in marine carbonates is used as a paleo-pH recorder and is one of the most promising paleo-carbonate chemistry proxies. Understanding the thermodynamic basis of the proxy is of fundamental importance, including knowledge on the equilibrium fractionation factor between dissolved boric acid, B(OH)₃, and borate ion, B(OH)₄⁻ (, hereafter α_(B3-B4)). However, this factor has hitherto not been determined experimentally and a theoretically calculated value (Kakihana and Kotaka, 1977, hereafter KK77) has therefore been widely used. I examine the calculations underlying this value. Using the same spectroscopic data and methods as KK77, I calculate the same α_(B3−B4) = 1.0193 at 300 K. Unfortunately, it turns out that in general the result is sensitive to the experimentally determined vibrational frequencies and the theoretical methods used to calculate the molecular forces. Using analytical techniques and ab initio molecular orbital theory, the outcome for α_(B3-B4) varies between ∼1.020 and ∼1.050 at 300 K. However, several arguments suggest that α_(B3-B4) ? 1.030. Measured isotopic shifts in various ¹⁰B-, ²D-, and ¹⁸O-labeled isotopomers do not provide a constraint on stable boron isotope fractionation. I conclude that in order to anchor the fundamentals of the boron pH proxy, experimental work is required. The critics of the boron pH proxy should note, however, that uncertainties in α_(B3-B4) do not bias pH reconstructions provided that organism-specific calibrations are used.     

Melt/biotite B/B isotopic fractionation and the boron local environment in the structure of volcanic glasses

This paper is focused on the role of boron coordination in determining the ¹¹B/¹⁰B isotopic fractionation between melt/glass and biotite at magmatic temperatures. For this purpose, three evolved volcanic rocks from Roccastrada, Mt. Amiata, and Mt. Cimini belonging to the Neogene-Quaternary magmatism of central Italy were studied. In these samples, the measured boron biotite-glass partition coefficient ranges between 0.004 and 0.011, indicating that boron behaves as an incompatible element during biotite crystallization. The ¹¹B magic-angle spinning nuclear magnetic resonance (NMR) spectra reveal the presence of trigonal BO_3/2 units, tetrahedral BO_4/2⁻ sites, and three-coordinated BO_2/2O⁻ species containing one nonbridging oxygen. The relative contributions of these different boron sites were estimated by spectral deconvolution, and it was observed that the fraction of trigonally coordinated boron decreases with increasing K₂O concentration in the glass. The ¹¹B/¹⁰B isotopic fractionation between biotite and melt/glass was observed to be large even at magmatic temperatures and was found to be 1.0066 (Roccastrada sample), 1.00535 (Mt. Amiata sample), and 1.00279 (Mt. Cimini sample). Fractionation is mostly related to the relative amount of trigonal and tetrahedral boron sites in the glass network rather than to other processes, including the speciation of hydrous species in the glass structure. The measured α values are significantly higher than the calculated ones obtained using the reduced partition function ratios (RPFRs) for B(OH)₃ and B(OH)₄⁻ as reported by Kakihana et al. (1977) and the abundance of trigonal and tetrahedral boron obtained by ¹¹B NMR spectra. Furthermore, a nonlinear relationship is observed between the percentage of BO₄ in the glass structure and the measured 1000lnα, suggesting that the approximation of monomeric B(OH)₃ and B(OH)₄⁻ species contributions through ideal mixing in calculating the RPFRs in polyanions (Oi et al., 1989) probably does not apply to silicate glasses.The large B isotopic fractionation measured between glass and biotite and its dependence on the boron coordination in the glass are a limitation to the use of δ¹¹B in the mineral to characterize magmas. Nonetheless, the high incompatible behavior of boron in the most common magmatic minerals rules out that fractional crystallization significantly modified the B isotopic composition of the melt.     

Controlling factors of the δ 11B-pH proxy and its research direction

Significant boron isotope fractionation occurs in nature (?70 ‰ to +75 ‰) due to the high geochemical reactivity of boron and the large relative mass difference between ¹⁰B and ¹¹B. Since the 1990s, reconstruction of ancient seawater pH using the isotopic composition of boron in bio-carbonates (δ ¹¹B_carb), and then calculation of the past pCO₂ have become important issues for the international isotope geochemistry community, and are called the δ ¹¹B-pH proxy. Although many achievements have been made by this proxy, various aspects of boron systematics require rigorous evaluation. Based on the previous researches, mechanism of boron isotope fractionation, variation of boron isotope (δ ¹¹B) in nature (especially in bio-carbonates) and controlling factors of the δ ¹¹B-pH proxy, such as the dissociation constant of B(OH)₃ in seawater (pK_a), the δ ¹¹B of seawater (δ ¹¹B_SW), the boron isotopic fractionation factor between B(OH) ₄ ^? and B(OH)₃ (α _4–3), and the incorporated species of boron into bio-carbonates, are reviewed in detail and the research directions of this proxy are proposed. Generally, the controversy about pK_a, δ ¹¹B_sw, and α _4–3 is relatively less, but whether boron incorporated into bio-carbonates only in the form of B(OH) ₄ ^? remains doubtful. In the future, it is required that the physicochemical processes that control boron incorporation into carbonates be rigorously characterized and that the related chemical and isotopic fractionation be quantified. It is also necessary and important to establish a “best-fit empirically equation” between δ ¹¹B_carb and pH of seawater based on the precipitation experiments of inorganic or culture experiments of corals or foraminifera. In addition, extended application of the δ ¹¹B-pH proxy to the earlier part of the Phanerozoic relying on the Brachiopods is worthy of studying. Like other geochemical indicators, there are limiting factors of δ ¹¹B; however, it remains a very powerful tool in the reconstruction of past seawater pH at present.     

A critical evaluation of the boron isotope-pH proxy: The accuracy of ancient ocean pH estimates

The boron isotope-pH technique is founded on a theoretical model of carbonate δ¹¹B variation with pH that assumes that the boron isotopic composition of carbonates mirrors the boron isotopic composition of borate in solution (δ¹¹B_carb = δ¹¹B_borate). Knowledge of the fractionation factor for isotope exchange between boric acid and borate in solution (α_4-3), the equilibrium constant for the dissociation of boric acid (pK_B*), as well as the isotopic composition of boron in seawater (δ¹¹B_sw) are required parameters of the model.The available data suggests that both the value of α_4-3 and the history of δ¹¹B_sw are poorly constrained. However, if one assumes that δ¹¹B_carb = δ¹¹B_borate, an empirical value for α_4-3 can be estimated from the results of inorganic carbonate precipitation experiments. This exercise yields an α_4-3 value of ∼0.974 in accordance with recent theoretical estimates, but substantially deviates from the theoretical value of 0.981 often used to estimate paleo-ocean pH. Re-evaluation of ocean pH using an α_4-3 value of 0.974 and published foraminiferal δ¹¹B values for the Cenozoic yield pH estimates that are relatively invariant, but unrealistically high (∼8.4-8.6). Uncertainty increases as foraminiferal ‘vital effects’ are considered and different models for secular changes in seawater δ¹¹B are applied.The inability to capture realistic ocean pH possibly reflects on our understanding of the isotopic relationship between carbonate and borate, as well as the mechanism of boron incorporation in carbonates. Given the current understanding of boron systematics, pH values estimated using this technique have considerable uncertainty, particularly when reconstructions exceed the residence time of boron in the ocean.     

硼同位素分馏的实验理论认识和矿床地球化学研究进展

硼是一种中等挥发性元素,具有¹¹B和¹⁰B两个稳定同位素。两个同位素间高达10%的相对质量差使其在地质过程中引起高达-70‰至+75‰的硼同位素变化。硼在自然界主要与氧键合形成三配位(BO₃)和四配位(BO₄)结构,因而¹¹B和¹⁰B间同位素分馏主要受控于三配体(BO₃)和四面体(BO₄)间配分。本文综述了低温和高温地质过程的硼同位素分馏的理论和实验研究进展。在溶液中B(OH)₃和${B(OH)^{-}_{4}}$间硼同位素分馏受pH和热力学p-T条件控制,实验和理论表征获得常温常压条件下的B(OH)₃和$B(OH)^{-}_{4}$间同位素分馏系数(α_3-4)变化范围为1.019 4至1.033 3。低温条件下矿物(如碳酸盐、黏土矿物(蒙脱石和伊利石)、针铁矿、水锰矿、硼酸盐)与溶液间硼同位素分馏行为除了受p-T-pH影响外,矿物表面吸附引起的分馏效应十分显著。在中高温过程(蒙脱石伊利石化、富硼电气石和白云母矿物与热液流体,以及硅酸盐熔体与流体)中硼同位素分馏行为受到硼配位构型、化学成分以及物理化学条件的控制。随着硼同位素分馏机理研究的深入以及越来越完善的地质储库硼同位素端员特征表征,硼同位素地球化学指标可以灵敏示踪成矿物质来源、探究成矿作用与成因模式和重建成矿过程物理化学条件。目前矿床硼同位素地球化学研究的难点在于实现不同赋存相(如流体、矿物和熔体)中硼配位键合结构和硼同位素组成的精细化表征。     

Liquid-vapor fractionation of boron and boron isotopes: Experimental calibration at 400°C/23 MPa to 450°C/42 MPa

We experimentally determined the boron partitioning and boron isotope fractionation between coexisting liquid and vapor in the system H₂O−NaCl−B₂O₃. Experiments were performed along the 400 and 450°C isotherms. Pressure conditions ranged from 23 to 28 MPa at 400°C and from 38 to 42 MPa at 450°C. Boron partitions preferentially into the liquid. Its overall liquid-vapor fractionation is, however, weak: Calculated boron distribution coefficients D_B^liquid-vapor are < 2.5 at all run conditions. With decreasing pressure (i.e. increasing opening of the solvus) D_B^liquid-vapor increases along the individual isotherms. Extrapolation to salt saturated conditions yields maximum boron liquid-vapor fractionations of D_B^liquid-vapor = 1.8 at 450°C and D_B^liquid-vapor = 2.7 at 400°C. ¹¹B preferentially fractionates into the vapor. Calculated Δ¹¹B_vapor-liquid = {[(¹¹B/¹⁰B)_vapor - (¹¹B/¹⁰B)_liquid]/(¹¹B/¹⁰B)_{NBS 951}}*1000 are small and range from 0.2 (± 0.7) to 0.9 (± 0.5) ‰ at 450°C and from 0.1 (± 0.6) to 0.7 (± 0.6) ‰ at 400°C. The data indicate increasing isotopic fractionation with decreasing pressure (i.e. increasing opening of the solvus). Extrapolation to salt saturated conditions yields maximum boron isotope liquid-vapor fractionations of Δ¹¹B_vapor-liquid = 1.5 (± 0.7) ‰ at 450°C and Δ¹¹B_vapor-liquid = 1.3 (± 0.6) ‰ at 400°C. The weak boron isotope fractionation suggests similar trigonal speciation in liquid and vapor. Although the boron and boron isotope fractionation between liquid and vapor is only weak, mass balance calculations indicate that for high degrees of fractionation liquid-vapor phase separation in an open system can significantly alter the boron and boron isotope signature of low-salinity hydrous fluids in hydrothermal systems. Comparing the model calculations with natural oceanic hydrothermal fluids, however, indicate that other processes than fluid phase separation dominate the boron geochemistry in oceanic hydrothermal fluids.     

Iron isotope fractionation between pyrite (FeS2), hematite (Fe2O3) and siderite (FeCO3): A first-principles density functional theory study

In addition to equilibrium isotopic fractionation factors experimentally derived, theoretical predictions are needed for interpreting isotopic compositions measured on natural samples because they allow exploring more easily a broader range of temperature and composition. For iron isotopes, only aqueous species were studied by first-principles methods and the combination of these data with those obtained by different methods for minerals leads to discrepancies between theoretical and experimental isotopic fractionation factors. In this paper, equilibrium iron isotope fractionation factors for the common minerals pyrite, hematite, and siderite were determined as a function of temperature, using first-principles methods based on the density functional theory (DFT). In these minerals belonging to the sulfide, oxide and carbonate class, iron is present under two different oxidation states and is involved in contrasted types of interatomic bonds. Equilibrium fractionation factors calculated between hematite and siderite compare well with the one estimated from experimental data (ln α⁵⁷Fe/⁵⁴Fe = 4.59 ± 0.30‰ and 5.46 ± 0.63‰ at 20 °C for theoretical and experimental data, respectively) while those for Fe(III)_aq-hematite and Fe(II)_aq-siderite are significantly higher that experimental values. This suggests that the absolute values of the reduced partition functions (β-factors) of aqueous species are not accurate enough to be combined with those calculated for minerals. When compared to previous predictions derived from Mössbauer or INRXS data [Polyakov V. B., Clayton R. N., Horita J. and Mineev S. D. (2007) Equilibrium iron isotope fractionation factors of minerals: reevaluation from the data of nuclear inelastic resonant X-ray scattering and Mössbauer spectroscopy. Geochim. Cosmochim. Acta71, 3833-3846], our iron β-factors are in good agreement for siderite and hematite while a discrepancy is observed for pyrite. However, the detailed investigation of the structural, electronic and vibrational properties of pyrite as well as the study of sulfur isotope fractionation between pyrite and two other sulfides (sphalerite and galena) indicate that DFT-derived β-factors of pyrite are as accurate as for hematite and siderite. We thus suggest that experimental vibrational density of states of pyrite should be re-examined.     

Boron isotopic fractionation related to boron sorption on humic acid and the structure of surface complexes formed

Boron isotopic fractionation during adsorption onto Ca-flocculated Aldrich humic acid (HA) has been investigated experimentally as a function of solution pH at 25°C and I = 0.15 M. Boron aqueous concentration and isotopic composition were determined by Cs₂BO₂⁺ Positive Thermal Ionization Mass Spectrometry analysis, while the structure of B surface complexes on HA was characterized using ¹¹B Magic Angle Spinning Nuclear Magnetic Resonance (MAS NMR). Significant B sorption on HA was observed at 6 < pH < 12 with a maximum value of Kd, the partition coefficient between adsorbed and aqueous boron, equal to 40 at pH = 9.5-10. Combined ¹¹B MAS NMR analysis and FITEQL modeling of B sorption on HA showed that this element forms tetrahedrally coordinated five- or six-membered ring chelates, most likely 1,2-diol and 1,3-diol complexes at alkaline pH (8 < pH < 11) and dicarboxylic complexes at near neutral conditions (6 < pH < 9). Results of this study demonstrate for the first time that boron sorption on HA induces a strong pH-dependent isotope fractionation—with ¹¹B depleted at the surface of HA—that reaches a maximum at 5 < pH < 9 (α = 0.975, Δ = −25‰) and decreases sharply at pH >9. The measured isotope fractionation cannot be modeled assuming that the isotopic composition of the sorbed borate species is identical to that of B(OH)₄^- species in the parent solution. It is shown that the extent of isotopic fractionation depends not only on B aqueous speciation but also on the distribution and structure of the borate surface complexes formed. In agreement with energetic constrains, calculation of the isotope fractionation factors between aqueous boric acid and boron surface complexes suggests that the formation of the strained six-membered ring 1,3-diol complex yields a much higher fractionation (α_BLP1−III = 0.954-0.960, Δ = −41/-47‰) than that of the very stable five-membered ring 1,2-diol (α_BLP2−III = 0.983, Δ = −18‰). The results of this study open new perspectives to understand and model boron biogeochemical cycle. It is predicted that boron sorption onto organic matter can have important consequences for the boron isotopic composition of surface water reservoirs (seawater, groundwater, soil waters) in which either abundant organic surfaces or significant boron concentrations are available. In addition, the large isotope fractionation between aqueous boric acid and surface boron-organic complexes found in the present work makes boron a promising tracer of biologic activity.     

Ab initio molecular orbital calculations for boron isotope fractionations on boric acids and borates

Ab initio quantum chemistry calculations have been performed on the isotopic exchange reaction between B(OH)₃ and B(OH)₄⁻. Several calculation methods have been carefully compared and evaluated. The “water-droplet” method is chosen to investigate this isotope exchange reaction using cluster models with up to 34 water molecules surrounding the solute. HF/6-31G* level calculations coupled with a 0.920 scaling factor are used for the frequency calculations. A larger K value (1.027) is obtained from this study compared to the commonly used 1.0194 (Kakihana et al., 1977).The fractionations for several boric acid polymers and boron minerals are also studied. Our results suggest that assuming the BO₄ bonding in B(OH)₄⁻ is identical to that in borosilicates is wrong. Tetrahedral boron in silicates has a significantly smaller reduced isotopic partition function ratio (RPFR) and hence will be much isotopically lighter than in B(OH)₄⁻.The new theoretical curve of pH vs. δ¹¹B composition of B(OH)₄⁻ using our calculated 1.027 can be used to predict pH values for equilibrium cases such as incorporation into inorganic calcite. We also find that the shape of this curve is very sensitive to both K and pK_a value, giving the possibility of also predicting salinity from the different shapes of the curve.     

Inter-mineral Fe isotope variations in mantle-derived rocks and implications for the Fe geochemical cycle

Iron isotope and major- and minor-element compositions of coexisting olivine, clinopyroxene, and orthopyroxene from eight spinel peridotite mantle xenoliths; olivine, magnetite, amphibole, and biotite from four andesitic volcanic rocks; and garnet and clinopyroxene from seven garnet peridotite and eclogites have been measured to evaluate if inter-mineral Fe isotope fractionation occurs in high-temperature igneous and metamorphic minerals and if isotopic fractionation is related to equilibrium Fe isotope partitioning or a result of open-system behavior. There is no measurable fractionation between silicate minerals and magnetite in andesitic volcanic rocks, nor between olivine and orthopyroxene in spinel peridotite mantle xenoliths. There are some inter-mineral differences (up to 0.2‰ in ⁵⁶Fe/⁵⁴Fe) in the Fe isotope composition of coexisting olivine and clinopyroxene in spinel peridotites. The Fe isotope fractionation observed between clinopyroxene and olivine appears to be a result of open-system behavior based on a positive correlation between the Δ⁵⁶Fe_{clinopyroxene-olivine} fractionation and the δ⁵⁶Fe value of clinopyroxene and olivine. There is also a significant difference in the isotopic compositions of garnet and clinopyroxene in garnet peridotites and eclogites, where the average Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation is +0.32 ± 0.07‰ for six of the seven samples. The one sample that has a lower Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation of 0.08‰ has a low Ca content in garnet, which may reflect some crystal chemical control on Fe isotope fractionation. The Fe isotope variability in mantle-derived minerals is interpreted to reflect subduction of isotopically variable oceanic crust, followed by transport through metasomatic fluids. Isotopic variability in the mantle might also occur during crystal fractionation of basaltic magmas within the mantle if garnet is a liquidus phase. The isotopic variations in the mantle are apparently homogenized during melting processes, producing homogenous Fe isotope compositions during crust formation.     

Re-evaluating boron speciation in biogenic calcite and aragonite using B MAS NMR

Understanding the partitioning of aqueous boron species into marine carbonates is critical for constraining the boron isotope system for use as a marine pH proxy. Previous studies have assumed that boron was incorporated into carbonate through the preferential uptake of tetrahedral borate B(OH)₄⁻. In this study we revisit this assumption through a detailed solid state ¹¹B magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopic study of boron speciation in biogenic and hydrothermal carbonates. Our new results contrast with those of the only previous NMR study of carbonates insofar as we observe both trigonal and tetrahedral coordinated boron in almost equal abundances in our biogenic calcite and aragonite samples. In addition, we observe no strict dependency of boron coordination on carbonate crystal structure. These NMR observations coupled with our earlier re-evaluation of the magnitude of boron isotope fractionation between aqueous species suggest that controls on boron isotope composition in marine carbonates, and hence the pH proxy, are more complex that previously suggested.     

Intra-shell boron isotope ratios in the symbiont-bearing benthic foraminiferan Amphistegina lobifera: Implications for δB vital effects and paleo-pH reconstructions

The boron isotope composition of marine carbonates is considered to be a seawater pH proxy. Nevertheless, the use of δ¹¹B has some limitations such as the knowledge of the fractionation factor (α_4-3) between boric acid and the borate ion and the amplitude of “vital effects” on this proxy that are not well constrained. Using secondary ion mass spectrometry (SIMS) we have examined the internal variability of the boron isotope ratio in the shallow water, symbionts bearing foraminiferan Amphistegina lobifera. Specimens were cultured at constant temperature (24 ± 0.1 °C) in seawater with pH ranging between 7.90 and 8.45. Intra-shell boron isotopes showed large variability with an upper limit value of ≈30‰. Our results suggest that the fractionation factor α_4-3 of 0.97352 (Klochko et al., 2006) is in better agreement with our experiments and with direct pH measurements in seawater vacuoles associated with the biomineralization process in these foraminifera. Despite the large variability of the skeletal pH values in each cultured specimen, it is possible to link the lowest calculated pH values to the experimental culture pH values while the upper pH limit is slightly below 9. This variability can be interpreted as follows: foraminifera variably increase the pH at the biomineralization site to about 9. This increase above ambient seawater pH leads to a range in δ¹¹B (Δ¹¹B) for each seawater pH. This Δ¹¹B is linearly correlated with the culture seawater pH with a slope of −13.1 per pH unit, and is independent of the fractionation factor α_4-3, or the δ¹¹B_sw through time. It may also be independent of the pK_B (the dissociation constant of boric acid) value. Therefore, Δ¹¹B in foraminifera can potentially reconstruct paleo-pH of seawater.     

Hydrogen isotope fractionation between OH-bearing minerals and water

Hydrogen isotope fractionation factors between hydroxyl-bearing minerals and water were determined at temperatures ranging between 400 and 850°C. The hydrogen isotope exchange rates for the mineral-water pairs examined were very slow. In most cases it was necessary to use an interpolation method for the determination of the hydrogen isotope equilibrium fractionation factor, α_e.For the temperature range of 450–850°C the hydrogen isotope fractionation factors for the mica-water and amphibole-water systems are simply expressed as a function of temperature and the molar fractions of the six-fold coordinated cations in the crystal, regardless of mineral species, as follows: 10³ In α_{e(mineral-water)} = ? 22.4 (10⁶T^?2) + 28.2 + (2X_Al ? 4X_Mg ? 68X_Fe), where X is the molar fraction of the cations. As the equation indicates, for any specific composition of the OH-bearing minerals, the change of α_e with temperature, over the temperature range investigated, is the same for all minerals studied. Thus for any specified values of X_Al, X_Mg, and X_Fe for these minerals, the relationship between α_e and T is 10³ In α_e = αT^?2 + k. Consequently, hydrogen isotope fractionation among coexisting minerals is temperature independent and cannot be used as a hydrogen isotope geothermometer.Some exceptions to the above general observations exist for minerals such as boehmite and kaolinite. In these minerals hydrogen bonding modifies the equilibrium hydrogen isotopic fractionation between mineral and water.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(III) and hematite

Application of the Fe isotope system to studies of natural rocks and fluids requires precise knowledge of equilibrium Fe isotope fractionation factors among various aqueous Fe species and minerals. These are difficult to obtain at the low temperatures at which Fe isotope fractionation is expected to be largest and requires careful distinction between kinetic and equilibrium isotope effects. A detailed investigation of Fe isotope fractionation between [Fe^III(H₂O)₆]³⁺ and hematite at 98°C allows the equilibrium ⁵⁶Fe/⁵⁴Fe fractionation to be inferred, which we estimate at 10³lnα_{Fe(III)-hematite} = −0.10 ± 0.20‰. We also infer that the slope of Fe(III)-hematite fractionation is modest relative to 10⁶/T², which would imply that this fractionation remains close to zero at lower temperatures. These results indicate that Fe isotope compositions of hematite may closely approximate those of the fluids from which they precipitated if equilibrium isotopic fractionation is assumed, allowing inference of δ⁵⁶Fe values of ancient fluids from the rock record. The equilibrium Fe(III)-hematite fractionation factor determined in this study is significantly smaller than that obtained from the reduced partition function ratios calculated for [Fe^III(H₂O)₆]³⁺ and hematite based on vibrational frequencies and Mössbauer shifts by [Polyakov 1997] and [Polyakov and Mineev 2000], and Schauble et al. (2001), highlighting the importance of experimental calibration of Fe isotope fractionation factors. In contrast to the long-term (up to 203 d) experiments, short-term experiments indicate that kinetic isotope effects dominate during rapid precipitation of ferric oxides. Precipitation of hematite over ∼12 h produces a kinetic isotope fractionation where 10³lnα_{Fe(III)-hematite} = +1.32 ± 0.12‰. Precipitation under nonequilibrium conditions, however, can be recognized through stepwise dissolution in concentrated acids. As expected, our results demonstrate that dissolution by itself does not measurably fractionate Fe isotopes.     

Oxygen isotope salt effects at high pressure and high temperature and the calibration of oxygen isotope geothermometers

The influence of NaCl, CaCl₂, and dissolved minerals on the oxygen isotope fractionation in mineral-water systems at high pressure and high temperature was studied experimentally. The salt effects of NaCl (up to 37 molal) and 5-molal CaCl₂ on the oxygen isotope fractionation between quartz and water and between calcite and water were measured at 5 and 15 kbar at temperatures from 300 to 750°C. CaCl₂ has a larger influence than NaCl on the isotopic fractionation between quartz and water. Although NaCl systematically changes the isotopic fractionation between quartz and water, it has no influence on the isotopic fractionation between calcite and water. This difference in the apparent oxygen isotope salt effects of NaCl must relate to the use of different minerals as reference phases. The term oxygen isotope salt effect is expanded here to encompass the effects of dissolved minerals on the fractionations between minerals and aqueous fluids. The oxygen isotope salt effects of dissolved quartz, calcite, and phlogopite at 15 kbar and 750°C were measured in the three-phase systems quartz-calcite-water and phlogopite-calcite-water. Under these conditions, the oxygen isotope salt effects of the three dissolved minerals range from ∼0.7 to 2.1‰. In both three-phase hydrothermal systems, the equilibrium fractionation factors between the pairs of minerals are the same as those obtained by anhydrous direct exchange between each pair of minerals, proving that the use of carbonate as exchange medium provides correct isotopic fractionations for a mineral pair.When the oxygen isotope salt effects of two minerals are different, the use of water as an indirect exchange medium will give erroneous fractionations between the two minerals. The isotope salt effect of a dissolved mineral is also the main reason for the observation that the experimentally calibrated oxygen isotope fractionations between a mineral and water are systematically 1.5 to 2‰ more positive than the results of theoretical calculations. Dissolved minerals greatly affect the isotopic fractionation in mineral-water systems at high pressure and high temperature. If the presence of a solute changes the solubility of a mineral, the real oxygen isotope salt effect of the solute at high pressure and high temperature cannot be correctly derived by using the mineral as reference phase.     

Boron isotopes as pH proxy: A new look at boron speciation in deep-sea corals using B MAS NMR and EELS

Dissolved boron in modern seawater occurs in the form of two species, trigonal boric acid B(OH)₃ and tetrahedral borate ion . One of the key assumption in the use of boron isotopic compositions of carbonates as pH proxy is that only borate ions, , are incorporated into the carbonate. Here, we investigate the speciation of boron in deep-sea coral microstructures (Lophelia pertusa specimen) by using high field magic angle spinning nuclear magnetic resonance (¹¹B MAS NMR) and electron energy-loss spectroscopy (EELS). We observe both boron coordination species, but in different proportions depending on the coral microstructure, i.e. centres of calcification versus fibres. These results suggest that careful sampling is necessary before performing boron isotopic measurements in deep-sea corals. By combining the proportions of B(OH)₃ and determined by NMR and our previous ion microprobe boron isotope measurements, we propose a new equation for the relation between seawater pH and boron isotopic composition in deep-sea corals.     

Modelling the geochemical cycle of boron: Implications for the long-term δ11B evolution of seawater and oceanic crust

The boron geochemical cycle has been simulated using a time-dependent geochemical box model that was coupled to a one-dimension model of seawater–oceanic crust interactions. Boron elemental and isotopic compositions of oceanic rocks as a function of depth were calculated by mass balance, using the temperature and porosity profiles of the crust as well as the available experimental and empirical distribution coefficients and fractionation factors between mineral and water. Ranges of boron elemental and isotopic variations of seawater were calculated for crust–seawater interactions that take place from the ridge-axis to the off-axis closure of the hydrothermal system. The present-day δ¹¹B of seawater (40‰) could represent a steady-state value. However, depending on crustal permeability, lifetime of water–rock interactions, and expansion rate of the oceanic ridge, the δ¹¹B of seawater may vary from 30‰ to 50‰ at the 10 million year scale.Some boron isotope compositions of Cretaceous biogenic carbonates and ophiolitic serpentinites from Oman are comparable to modern rock samples, suggesting that the δ¹¹B of Cretaceous seawater was close to the present-day value. Low δ¹¹B values of some biogenic carbonates cannot be attributed to low pH values of past seawater, but more probably to δ¹¹B variations of seawater or diagenetic alteration by crustal aqueous fluids. Boron isotope composition of hydrothermally altered serpentines could be considered as a promising proxy of the seawater composition.     

Boron isotope geochemistry of Paleozoic brachiopod calcite: Implications for a secular change in the boron isotope geochemistry of seawater over the Phanerozoic

Boron isotope composition of marine carbonates has been proposed as a paleo-pH proxy and potential tool to reconstruct atmospheric pCO₂. The precise knowledge of the boron isotopic composition of ancient seawater represents the fundamental prerequisite for any paleo-pH reconstruction. This contribution presents boron isotope values for Silurian to Permian brachiopod calcite that might be used to reconstruct pH or boron isotope composition of past oceans. All brachiopod shells were screened for diagenetic recrystallization by means of cathodoluminescence microscopy, trace element geochemistry (B, Fe, Mn, Sr) as well as SEM. Only nonluminescent shells revealing well-preserved microstructures, high strontium and boron concentrations as well as low iron and manganese contents were accepted for boron isotope analysis. The boron isotope ratios of Silurian, Devonian, Pennsylvanian and Permian brachiopod calcite range from 6.8 to 11.0‰, 7.3 to 14.9‰, 12.4 to 15.8‰ and 10.1 to 11.7‰, respectively. These δ¹¹B values are significantly lower in comparison to δ¹¹B values of modern biogenic carbonates and indicate that the Paleozoic oceans were depleted in ¹¹B by up to 10‰. Box modeling of the boron geochemical cycle suggests that the significant depletion of ¹¹B in the oceanic reservoir may have been initiated by an enhanced continental boron discharge. Our data support the earlier made conclusion that boron isotopes may not be used in the geological past as reliable paleo-pH proxy unless the boron isotopic composition of ancient oceans can be constrained by further studies.     

Theoretical investigation of iron isotope fractionation between Fe(H2O)6 and Fe(H2O)6: Implications for iron stable isotope geochemistry

The magnitude of equilibrium iron isotope fractionation between Fe(H₂O)₆³⁺ and Fe(H₂O)₆²⁺ is calculated using density functional theory (DFT) and compared to prior theoretical and experimental results. DFT is a quantum chemical approach that permits a priori estimation of all vibrational modes and frequencies of these complexes and the effects of isotopic substitution. This information is used to calculate reduced partition function ratios of the complexes (10³ · ln(β)), and hence, the equilibrium isotope fractionation factor (10³ · ln(α)). Solvent effects are considered using the polarization continuum model (PCM). DFT calculations predict fractionations of several per mil in ⁵⁶Fe/⁵⁴Fe favoring partitioning of heavy isotopes in the ferric complex. Quantitatively, 10³ · ln(α) predicted at 22°C, ∼ 3 ‰, agrees with experimental determinations but is roughly half the size predicted by prior theoretical results using the Modified Urey-Bradley Force Field (MUBFF) model. Similar comparisons are seen at other temperatures. MUBFF makes a number of simplifying assumptions about molecular geometry and requires as input IR spectroscopic data. The difference between DFT and MUBFF results is primarily due to the difference between the DFT-predicted frequency for the ν₄ mode (O-Fe-O deformation) of Fe(H₂O)₆³⁺ and spectroscopic determinations of this frequency used as input for MUBFF models (185-190 cm⁻¹ vs. 304 cm⁻¹, respectively). Hence, DFT-PCM estimates of 10³ · ln(β) for this complex are ∼ 20% smaller than MUBFF estimates. The DFT derived values can be used to refine predictions of equilibrium fractionation between ferric minerals and dissolved ferric iron, important for the interpretation of Fe isotope variations in ancient sediments. Our findings increase confidence in experimental determinations of the Fe(H₂O)₆³⁺ − Fe(H₂O)₆²⁺ fractionation factor and demonstrate the utility of DFT for applications in “heavy” stable isotope geochemistry.     

Combined iron and magnesium isotope geochemistry of pyroxenite xenoliths from Hannuoba,North China Craton: implications for mantle metasomatism

We present high-precision iron and magnesium isotopic data for diverse mantle pyroxenite xenoliths collected from Hannuoba, North China Craton and provide the first combined iron and magnesium isotopic study of such rocks. Compositionally, these xenoliths range from Cr-diopside pyroxenites and Al-augite pyroxenites to garnet-bearing pyroxenites and are taken as physical evidence for different episodes of melt injection. Our results show that both Cr-diopside pyroxenites and Al-augite pyroxenites of cumulate origin display narrow ranges in iron and magnesium isotopic compositions (δ⁵⁷Fe = ?0.01 to 0.09 with an average of 0.03 ± 0.08 (2SD, n = 6); δ²⁶Mg = ? 0.28 to ?0.25 with an average of ?0.26 ± 0.03 (2SD, n = 3), respectively). These values are identical to those in the normal upper mantle and show equilibrium inter-mineral iron and magnesium isotope fractionation between coexisting mantle minerals. In contrast, the garnet-bearing pyroxenites, which are products of reactions between peridotites and silicate melts from an ancient subducted oceanic slab, exhibit larger iron isotopic variations, with δ⁵⁷Fe ranging from 0.12 to 0.30. The δ⁵⁷Fe values of minerals in these garnet-bearing pyroxenites also vary widely (?0.25 to 0.08 in olivines, ?0.04 to 0.25 in orthopyroxenes, ?0.07 to 0.31 in clinopyroxenes, 0.07 to 0.48 in spinels and 0.31–0.42 in garnets). In addition, the garnet-bearing pyroxenite shows light δ²⁶Mg (?0.43) relative to the mantle. The δ²⁶Mg of minerals in the garnet-bearing pyroxenite range from ?0.35 for olivine and orthopyroxene, to ?0.34 for clinopyroxene, 0.04 for spinel and ?0.68 for garnet. These measured values stand in marked contrast to calculated equilibrium iron and magnesium isotope fractionation between coexisting mantle minerals at mantle temperatures derived from theory, indicating disequilibrium isotope fractionation. Notably, one phlogopite clinopyroxenite with an apparent later metasomatic overprint has the heaviest δ⁵⁷Fe (as high as 1.00) but the lightest δ²⁶Mg (as low as ?1.50) values of all investigated samples. Overall, there appears to be a negative co-variation between δ⁵⁷Fe and δ²⁶Mg in the Hannuoba garnet-bearing pyroxenite and in the phlogopite clinopyroxenite xenoliths and minerals therein. These features may reflect kinetic isotopic fractionation due to iron and magnesium inter-diffusion during melt–rock interaction. Such processes play an important role in producing inter-mineral iron and magnesium isotopic disequilibrium and local iron and magnesium isotopic heterogeneity in the subcontinental mantle.