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.     

Boron isotopic fractionation between minerals and fluids: New insights from in situ high pressure-high temperature vibrational spectroscopic data

Equilibrium boron isotopic fractionations between trigonal B(OH)₃ and tetragonal B(OH)₄⁻ aqueous species have been calculated at high P-T conditions using measured vibrational spectra (Raman and IR) and force-field modeling to compute reduced partition function ratios for B-isotopic exchange following Urey’s theory. The calculated isotopic fractionation factor at 300 K, α_3/4 = 1.0176(2), is slightly lower than the formerly calculated value of α_3/4 = 1.0193 (Kakihana and Kotaka, 1977), due to differences in the determined vibrational frequencies. The effect of pressure on α_3/4 up to 10 GPa and 723 K is shown to be negligible relative to temperature or speciation (pH) effects. Implications for the interpretation of boron fractionation in experimental and natural systems are discussed. We also show that the relationship between seawater-mineral B isotope fractionation and pH can be expressed using two variables, α_3/4 on one hand, and the pK_a of the boric acid-borate equilibrium on the other hand. This latter value is given by the equilibrium of boron species in water for the carbonate-water exchange, but could be governed by mineral surface properties in the case of clays. This may allow defining intrinsic paleo-pHmeters from B isotope fractionation between carbonate and authigenic minerals. Finally, it is shown that fractionation of boron isotopes can be rationalized in terms of the changes in 1) coordination of B from trigonal to tetrahedral in both fluids and minerals; and 2) the ligand nature around B from OH⁻ in the fluid and some hydrous minerals to non-hydrogenated O in many minerals. Relationships are established that allow predicting the isotopic fractionation factor of B between minerals and fluid.     

Boron isotopic composition of atmospheric precipitations and liquid-vapour fractionations

Boron isotope compositions (δ¹¹B) and B concentrations of rains and snows were studied in order to characterize the sources and fractionation processes during the boron atmospheric cycle. The ¹¹B/¹⁰B ratios of instantaneous and cumulative rains and snows from coastal and continental sites show a large range of variations, from −1.5 ± 0.4 to +26.0 ± 0.5‰ and from −10.2 ± 0.5 to +34.4 ± 0.2‰, respectively. Boron concentrations in rains and snows vary between 0.1 and 3.0 ppb. All these precipitation samples are enriched in ¹⁰B compared to the ocean value (δ¹¹B = +39.5‰). An empirical rain-vapour isotopic fractionation of +31‰ is estimated from three largely independent methods. The deduced seawater-vapour fractionation is +25.5‰, with the difference between the rain and seawater fractionations principally reflecting changes in the speciation of boron in the liquid with ∼100% B(OH)₃ present in precipitations. A boron meteoric water line, δD = 2.6δ¹¹B − 133, is proposed which describes the relationship between δD and δ¹¹B in many, but not all, precipitations. Boron isotopic compositions of precipitations can be related to that of the seawater reservoir by the seawater-vapour fractionation and one or more of (1) the rain-vapour isotopic fractionation, (2) evolution of the δ¹¹B value of the atmospheric vapour reservoir via condensation-precipitation processes (Rayleigh distillation process), (3) any contribution of vapour from the evaporation of seawater aerosols, and (4) any contribution from particulate matter, principally sea salt, continental dust and, perhaps more regionally, anthropogenic sources (burning of biomass and fossil fuels). From the δ¹¹B values of continental precipitations, a sea salt contribution cannot be more than a percent or so of the total B in precipitation over these areas.     

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.     

Effects of changing solution chemistry on Fe/Fe isotope fractionation in aqueous Fe-Cl solutions

The range in ⁵⁶Fe/⁵⁴Fe isotopic compositions measured in naturally occurring iron-bearing species is greater than 5‰. Both theoretical modeling and experimental studies of equilibrium isotopic fractionation among iron-bearing species have shown that significant fractionations can be caused by differences in oxidation state (i.e., redox effects in the environment) as well as by bond partner and coordination number (i.e., nonredox effects due to speciation).To test the relative effects of redox vs. nonredox attributes on total Fe equilibrium isotopic fractionation, we measured changes, both experimentally and theoretically, in the isotopic composition of an Fe²⁺-Fe³⁺-Cl-H₂O solution as the chlorinity was varied. We made use of the unique solubility of FeCl₄⁻ in immiscible diethyl ether to create a separate spectator phase against which changes in the aqueous phase could be quantified. Our experiments showed a reduction in the redox isotopic fractionation between Fe²⁺- and Fe³⁺-bearing species from 3.4‰ at [Cl⁻] = 1.5 M to 2.4‰ at [Cl⁻] = 5.0 M, due to changes in speciation in the Fe-Cl solution. This experimental design was also used to demonstrate the attainment of isotopic equilibrium between the two phases, using a ⁵⁴Fe spike.To better understand speciation effects on redox fractionation, we created four new sets of ab initio models of the ferrous chloride complexes used in the experiments. These were combined with corresponding ab initio models for the ferric chloride complexes from previous work. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 6.39‰ to 5.42‰ for Fe(H₂O)₆²⁺, 5.98‰ to 5.34‰ for FeCl(H₂O)₅⁺, and 5.91‰ to 4.86‰ for FeCl₂(H₂O)₄, depending on the model. The theoretical models predict ferric-ferrous fractionation about half as large (depending on model) as the experimental results.Our results show (1) oxidation state is likely to be the dominant factor controlling equilibrium Fe isotope fractionation in solution and (2) nonredox attributes (such as ligands present in the aqueous solution, speciation and relative abundances, and ionic strength of the solution) can also have significant effects. Changes in the isotopic composition of an Fe-bearing solution will influence the resultant Fe isotopic signature of any precipitates.     

Boron isotopes geochemistry of the Changjiang basin rivers

We report analyses of B isotopic compositions in water and suspended particulate matter collected in the Changjiang and its main tributaries. We showed that four sources control the dissolved boron budget; namely atmospheric deposition, evaporite dissolution, anthropogenic inputs and silicate weathering. The contribution of silicate weathering to the dissolved B load ranges from 40% to 50% for the Changjiang main channel and from 45% to 88% for the main tributaries. The isotopic composition of dissolved boron derived from silicate weathering range from −3‰ up to +9‰ suggesting that isotopic fractionation occurs during silicate weathering. The boron isotopic composition of suspended particulate matter range from −11.4‰ to −6‰. Boron derived from silicate weathering is preferentially carried out by the dissolved load which accounts for 30-96% of the total boron. We show that the isotopic compositions of both the dissolved load and suspended particulate matter are controlled by the competition between boron leaching and boron uptake into secondary phases. The first process is characterized by a loss of boron relative to the bedrock without apparent isotopic fractionation whereas the last one is associated to a large isotopic fractionation which enriches the dissolved boron in heavy isotope.     

Isotopic fractionation and reaction kinetics between Cr(III) and Cr(VI) in aqueous media

The redox-sensitive stable isotope geochemistry of chromium bears the potential to monitor the attenuation of chromate pollution and to investigate changes in environmental conditions in the present and the past. The use of stable Cr isotope data as a geo-environmental tracer, however, necessitates an understanding of the reaction kinetics and Cr fractionation behaviour during redox transition and isotope exchange. Here, we report stable chromium isotope fractionation data for Cr(VI) reduction, Cr(III) oxidation and isotopic exchange between soluble Cr(III) and Cr(VI) in aqueous media. The reduction of Cr(VI) to Cr(III) with H₂O₂ under strongly acidic conditions shows a near-equilibrium isotope fractionation of Δ^53/52Cr_{(Cr(III)-Cr(VI))} of −3.54 ± 0.35‰. At pH neutrality, however, the reduction experiments show a kinetic isotope fractionation Δ^53/52Cr_{(Cr(III)-Cr(VI))} of −5‰ for the extent of reduction of up to 85% of the chromium. The oxidation of Cr(III) to Cr(VI) in alkaline media, using H₂O₂ as the oxidant, cannot be explained by a single, unidirectional reaction. Our experiments indicate that the involvement of the unstable intermediates Cr(IV) and Cr(V) and their disproportionation during redox reactions between Cr(III) and Cr(VI) influence the overall fractionation factor, depending on the prevailing pH conditions and the reaction rates. No detectable isotope exchange between soluble Cr(VI) and Cr(III) species at pH values of 5.5 and 7 was revealed over a timescale of days to weeks. This means that, at least within such a time frame, the isotopic composition of Cr(VI) in a natural system will not be influenced by equilibration with any Cr(III) and thus reveal the true extent of reduction, given that the Cr isotope composition of the source Cr(VI) and the fractionation factor for the prevailing conditions are known.     

Copper isotope fractionation during its interaction with soil and aquatic microorganisms and metal oxy(hydr)oxides: Possible structural control

This work is aimed at quantifying the main environmental factors controlling isotope fractionation of Cu during its adsorption from aqueous solutions onto common organic (bacteria, algae) and inorganic (oxy(hydr)oxide) surfaces. Adsorption of Cu on aerobic rhizospheric (Pseudomonas aureofaciens CNMN PsB-03) and phototrophic aquatic (Rhodobacter sp. f-7bl, Gloeocapsa sp. f-6gl) bacteria, uptake of Cu by marine (Skeletonema costatum) and freshwater (Navicula minima, Achnanthidium minutissimum and Melosira varians) diatoms, and Cu adsorption onto goethite (FeOOH) and gibbsite (AlOOH) were studied using a batch reaction as a function of pH, copper concentration in solution and time of exposure. Stable isotopes of copper in selected filtrates were measured using Neptune multicollector ICP-MS. Irreversible incorporation of Cu in cultured diatom cells at pH 7.5-8.0 did not produce any isotopic shift between the cell and solution (Δ^65/63Cu(solid-solution)) within ±0.2‰. Accordingly, no systematic variation was observed during Cu adsorption on anoxygenic phototrophic bacteria (Rhodobacter sp.), cyanobacteria (Gloeocapsa sp.) or soil aerobic exopolysaccharide (EPS)-producing bacteria (P. aureofaciens) in circumneutral pH (4-6.5) and various exposure times (3 min to 48 h): Δ⁶⁵Cu(solid-solution) = 0.0 ± 0.4‰. In contrast, when Cu was adsorbed at pH 1.8-3.5 on the cell surface of soil the bacterium P. aureofacienshaving abundant or poor EPS depending on medium composition, yielded a significant enrichment of the cell surface in the light isotope (Δ⁶⁵Cu (solid-solution) = −1.2 ± 0.5‰). Inorganic reactions of Cu adsorption at pH 4-6 produced the opposite isotopic offset: enrichment of the oxy(hydr)oxide surface in the heavy isotope with Δ⁶⁵Cu(solid-solution) equals 1.0 ± 0.25‰ and 0.78 ± 0.2‰ for gibbsite and goethite, respectively. The last result corroborates the recent works of Mathur et al. [Mathur R., Ruiz J., Titley S., Liermann L., Buss H. and Brantley S. (2005) Cu isotopic fractionation in the supergene environment with and without bacteria. Geochim. Cosmochim. Acta69, 5233-5246] and Balistrieri et al. [Balistrieri L. S., Borrok D. M., Wanty R. B. and Ridley W. I. (2008) Fractionation of Cu and Zn isotopes during adsorption onto amorhous Fe(III) oxyhydroxide: experimental mixing of acid rock drainage and ambient river water. Geochim. Cosmochim. Acta72, 311-328] who reported heavy Cu isotope enrichment onto amorphous ferric oxyhydroxide and on metal hydroxide precipitates on the external membranes of Fe-oxidizing bacteria, respectively.Although measured isotopic fractionation does not correlate with the relative thermodynamic stability of surface complexes, it can be related to their structures as found with available EXAFS data. Indeed, strong, bidentate, inner-sphere complexes presented by tetrahedrally coordinated Cu on metal oxide surfaces are likely to result in enrichment of the heavy isotope on the surface compared to aqueous solution. The outer-sphere, monodentate complex, which is likely to form between Cu²⁺ and surface phosphoryl groups of bacteria in acidic solutions, has a higher number of neighbors and longer bond distances compared to inner-sphere bidentate complexes with carboxyl groups formed on bacterial and diatom surfaces in circumneutral solutions. As a result, in acidic solution, light isotopes become more enriched on bacterial surfaces (as opposed to the surrounding aqueous medium) than they do in neutral solution.Overall, the results of the present study demonstrate important isotopic fractionation of copper in both organic and inorganic systems and provide a firm basis for using Cu isotopes for tracing metal transport in earth-surface aquatic systems. It follows that both adsorption on oxides in a wide range of pH values and adsorption on bacteria in acidic solutions are capable of producing a significant (up to 2.5-3‰ (±0.1-0.15‰)) isotopic offset. At the same time, Cu interaction with common soil and aquatic bacteria, as well as marine and freshwater diatoms, at 4 < pH < 8 yields an isotopic shift of only ±0.2-0.3‰, which is not related to Cu concentration in solution, surface loading, the duration of the experiment, or the type of aquatic microorganisms.     

Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Effects of speciation on equilibrium fractionations and rates of oxygen isotope exchange between (PO4)aq and H2O

The effects of phosphate speciation on both rates of isotopic exchange and oxygen isotope equilibrium fractionation factors between aqueous phosphate and water were examined over the temperature range 70 to 180°C. Exchange between phosphate and water is much faster at low pH than at high pH, an observation that is similar to what has been observed in the analogous sulfate-water system. Oxygen isotope fractionations between protonated species like H₃PO₄ and H₂PO₄⁻ that are dominant at relatively low pH and species like PO₄³⁻ and ion pairs like KHPO₄⁻ that are dominant at relatively high pH, range between 5 and 8‰ at the temperatures of the experiments. In aqueous phosphate systems at equilibrium, ¹⁸O/¹⁶O ratios increase with increasing degree of protonation of phosphate. This effect can be explained in part by the relative magnitudes of the dissociation constants of the protonated species. Under equilibrium conditions, carbonate in solution or in solid phases concentrates ¹⁸O relative to orthophosphate in solution or in solid phases at all temperatures, supporting the traditional view that biogenic phosphate is precipitated in near oxygen isotope equilibrium with body/ambient aqueous fluids with no attendant vital effects.     

High precision iron isotope measurements of meteoritic material by cold plasma ICP-MS

The first cold plasma ICP-MS (inductively coupled plasma mass spectrometer) Fe isotope study is described. Application of this technique to the analyses of Fe isotopes in a number of meteorites is also reported. The measurement technique relies on reduced temperature operation of the ICP source to eliminate pervasive molecular interferences from Ar complexes associated with conventional ICP-MS. Instrumental mass bias corrections are performed by sample-standard bracketing and using Cu as an external mass bias drift monitor. Repeated measurements of a terrestrial basalt reference sample indicate an external reproducibility of ± 0.06 ‰ for δ⁵⁶Fe and ± 0.25 ‰ for δ⁵⁸Fe (1 σ). The measured iron isotopic compositions of various bulk meteorites, including irons, chondrites and pallasites are identical, within error, to the composition of our terrestrial basalt reference sample suggesting that iron mass fractionation during planet formation and differentiation was non-existent. Iron isotope compositions measured for eight chondrules from the unequilibrated ordinary chondrite Tieschitz range from −0.5 ‰ < δ⁵⁶Fe_chondrules < 0.0 ‰ relative to the terrestrial/meteorite average. Mechanisms for fractionating iron in these chondrules are discussed.     

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

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

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.     

Mg isotope constraints on soil pore-fluid chemistry: Evidence from Santa Cruz, California

Mg isotope ratios (²⁶Mg/²⁴Mg) are reported in soil pore-fluids, rain and seawater, grass and smectite from a 90 kyr old soil, developed on an uplifted marine terrace from Santa Cruz, California. Rain water has an invariant ²⁶Mg/²⁴Mg ratio (expressed as δ²⁶Mg) at −0.79 ± 0.05‰, identical to seawater δ²⁶Mg. Detrital smectite (from the base of the soil profile, and therefore unweathered) has a δ²⁶Mg value of 0.11‰, potentially enriched in ²⁶Mg by up to 0.3‰ compared to the bulk silicate Earth Mg isotope composition (although within the range of all terrestrial silicates). The soil pore-waters show a continuous profile with depth for δ²⁶Mg, ranging from −0.99‰ near the surface to −0.43‰ at the base of the profile. Shallow pore-waters (<1 m) have δ²⁶Mg values that are similar to, or slightly lower than the rain waters. This implies that the degree of biological cycling of Mg in the pore-waters is relatively small and is quantified as <32%, calculated using the average Mg isotope enrichment factor between grass and rain (δ²⁶Mg_grass-δ²⁶Mg_rain) of 0.21‰. The deep pore-waters (1-15 m deep) have δ²⁶Mg values that are intermediate between the smectite and rain, ranging from −0.76‰ to −0.43‰, and show a similar trend with depth compared to Sr isotope ratios. The similarity between Sr and Mg isotope ratios confirms that the Mg in the pore-waters can be explained by a mixture between rain and smectite derived Mg, despite the fact that Mg and Sr concentrations may be buffered by the exchangeable reservoir. However, whilst Sr isotope ratios in the pore-waters span almost the complete range between mineral and rain inputs, Mg isotopes compositions are much closer to the rain inputs. If Mg and Sr isotope ratios are controlled uniquely by a mixture, the data can be used to estimate the mineral weathering inputs to the pore-waters, by correcting for the rain inputs. This isotopic correction is compared to the commonly used chloride correction for precipitation inputs. A consistent interpretation is only possible if Mg isotope ratios are fractionated either by the precipitation of a secondary Mg bearing phase, not detected by conventional methods, or selective leaching of ²⁴Mg from smectite. There is therefore dual control on the Mg isotopic composition of the pore-waters, mixing of two inputs with distinct isotopic compositions, modified by fractionation. The data provide (1) further evidence for Mg isotope fractionation at the surface of the Earth and (2) the first field evidence of Mg isotope fractionation during uptake by natural plants. The coherent behaviour of Mg isotope ratios in soil environments is encouraging for the development of Mg isotope ratios as a quantitative tracer of both weathering inputs of Mg to waters, and the physicochemical processes that cycle Mg, a major cation linked to the carbon cycle, during continental weathering.     

Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments

The chemical and isotopic composition of speleothem calcite and particularly that of stalagmites and flowstones is increasingly exploited as an archive of past environmental change in continental settings. Despite intensive research, including modelling and novel approaches, speleothem data remain difficult to interpret. A possible way foreword is to apply a multi-proxy approach including non-conventional isotope systems. For the first time, we here present a complete analytical dataset of magnesium isotopes (δ²⁶Mg) from a monitored cave in NW Germany (Bunker Cave). The data set includes δ²⁶Mg values of loess-derived soil above the cave (−1.0 ± 0.5‰), soil water (−1.2 ± 0.5‰), the carbonate hostrock (−3.8 ± 0.5‰), dripwater in the cave (−1.8 ± 0.2‰), speleothem low-Mg calcite (stalactites, stalagmites; −4.3 ± 0.6‰), cave loam (−0.6 ± 0.1‰) and runoff water (−1.8 ± 0.1‰) in the cave, respectively. Magnesium-isotope fractionation processes during weathering and interaction between soil cover, hostrock and solute-bearing soil water are non-trivial and depend on a number of variables including solution residence times, dissolution rates, adsorption effects and potential neo-formation of solids in the regolith and the carbonate aquifer. Apparent Mg-isotope fractionation between dripwater and speleothem low-Mg calcite is about 1000lnα_Mg-cc-Mg(aq) = −2.4‰. A similar Mg-isotope fractionation (1000lnα_Mg-cc-Mg(aq) ≈ −2.1‰) is obtained by abiogenic precipitation experiments carried out at aqueous Mg/Ca ratios and temperatures close to cave conditions. Accordingly, ²⁶Mg discrimination during low-Mg calcite formation in caves is highly related to inorganic fractionation effects, which may comprise dehydration of Mg²⁺ prior to incorporation into calcite, surface entrapment of light isotopes and reaction kinetics. Relevance of kinetics is supported by a significant negative correlation of Mg-isotope fractionation with the precipitation rate for inorganic precipitation experiments.     

Chemical and isotopic constraints on water/rock interactions at the Lost City hydrothermal field, 30°N Mid-Atlantic Ridge

Low temperature vent fluids (<91 °C) issuing from the ultramafic-hosted hydrothermal system at Lost City, 30°N Mid-Atlantic Ridge, are enriched in dissolved volatiles (H₂,CH₄) while attaining elevated pH values, indicative of the serpentization processes that govern water/rock interactions deep in the oceanic crust. Here, we present a series of theoretical models to evaluate the extent of hydrothermal alteration and assess the effect of cooling on the systematics of pH-controlled B aqueous species. Peridotite-seawater equilibria calculations indicate that the mineral assemblage composed of diopside, brucite and chrysotile likely dictates fluid pH at moderate temperature serpentinization processes (<300 °C), by imposing constraints on the aCa⁺⁺/a²H⁺ ratios and the activity of dissolved SiO₂. Based on Sr abundances and the ⁸⁷Sr/⁸⁶Sr isotope ratios of vent fluids reported from Lost City, estimated water/rock mass ratios (w/r = 2-4) are consistent with published models involving dissolved CO₂ and alkane concentrations. Combining the reported δ¹⁸O values of vent fluids (0.7‰) with such w/r mass ratios, allows us to bracket subseafloor reaction temperatures in the vicinity of 250 °C. These estimates are in agreement with previous theoretical studies supporting extensive conductive heat loss within the upflow zones. Experimental studies on peridotite-seawater alteration suggest that fluid pH increases during cooling which then rapidly enhances boron removal from solution and incorporation into secondary phases, providing an explanation for the highly depleted dissolved boron concentrations measured in the low temperature but alkaline Lost City vent fluids. Finally, to account for the depleted ¹¹B composition (δ¹¹B ∼25-30‰) of vent fluids relative to seawater, isotopic fractionation between tetrahedrally coordinated aqueous boron species with BO₃-bearing mineral sites (e.g. in calcite, brucite) is proposed.     

Iron isotope fractionation and atom exchange during sorption of ferrous iron to mineral surfaces

The application of stable Fe isotopes as a tracer of the biogeochemical Fe cycle necessitates a mechanistic knowledge of natural fractionation processes. We studied the equilibrium Fe isotope fractionation upon sorption of Fe(II) to aluminum oxide (γ-Al₂O₃), goethite (α-FeOOH), quartz (α-SiO₂), and goethite-loaded quartz in batch experiments, and performed continuous-flow column experiments to study the extent of equilibrium and kinetic Fe isotope fractionation during reactive transport of Fe(II) through pure and goethite-loaded quartz sand. In addition, batch and column experiments were used to quantify the coupled electron transfer-atom exchange between dissolved Fe(II) (Fe(II)_aq) and structural Fe(III) of goethite. All experiments were conducted under strictly anoxic conditions at pH 7.2 in 20 mM MOPS (3-(N-morpholino)-propanesulfonic acid) buffer and 23 °C. Iron isotope ratios were measured by high-resolution MC-ICP-MS. Isotope data were analyzed with isotope fractionation models. In batch systems, we observed significant Fe isotope fractionation upon equilibrium sorption of Fe(II) to all sorbents tested, except for aluminum oxide. The equilibrium enrichment factor, , of the Fe(II)_sorb-Fe(II)_aq couple was 0.85 ± 0.10‰ (±2σ) for quartz and 0.85 ± 0.08‰ (±2σ) for goethite-loaded quartz. In the goethite system, the sorption-induced isotope fractionation was superimposed by atom exchange, leading to a δ^56/54Fe shift in solution towards the isotopic composition of the goethite. Without consideration of atom exchange, the equilibrium enrichment factor was 2.01 ± 0.08‰ (±2σ), but decreased to 0.73 ± 0.24‰ (±2σ) when atom exchange was taken into account. The amount of structural Fe in goethite that equilibrated isotopically with Fe(II)_aq via atom exchange was equivalent to one atomic Fe layer of the mineral surface (∼3% of goethite-Fe). Column experiments showed significant Fe isotope fractionation with δ^56/54Fe(II)_aq spanning a range of 1.00‰ and 1.65‰ for pure and goethite-loaded quartz, respectively. Reactive transport of Fe(II) under non-steady state conditions led to complex, non-monotonous Fe isotope trends that could be explained by a combination of kinetic and equilibrium isotope enrichment factors. Our results demonstrate that in abiotic anoxic systems with near-neutral pH, sorption of Fe(II) to mineral surfaces, even to supposedly non-reactive minerals such as quartz, induces significant Fe isotope fractionation. Therefore we expect Fe isotope signatures in natural systems with changing concentration gradients of Fe(II)_aq to be affected by sorption.     

First experimental determination of iron isotope fractionation between hematite and aqueous solution at hydrothermal conditions

Although iron isotopes provide a new powerful tool for tracing a variety of geochemical processes, the unambiguous interpretation of iron isotope ratios in natural systems and the development of predictive theoretical models require accurate data on equilibrium isotope fractionation between fluids and minerals. We investigated Fe isotope fractionation between hematite (Fe₂O₃) and aqueous acidic NaCl fluids via hematite dissolution and precipitation experiments at temperatures from 200 to 450 °C and pressures from saturated vapor pressure (P_sat) to 600 bar. Precipitation experiments at 200 °C and P_sat from aqueous solution, in which Fe aqueous speciation is dominated by ferric iron (Fe^III) chloride complexes, show no detectable Fe isotope fractionation between hematite and fluid, Δ⁵⁷Fe_{fluid-hematite} = δ⁵⁷Fe_fluid − δ⁵⁷Fe_hematite = 0.01 ± 0.08‰ (2 × standard error, 2SE). In contrast, experiments at 300 °C and P_sat, where ferrous iron chloride species (FeCl₂ and FeCl⁺) dominate in the fluid, yield significant fluid enrichment in the light isotope, with identical values of Δ⁵⁷Fe_{fluid-hematite} = −0.54 ± 0.15‰ (2SE) both for dissolution and precipitation runs. Hematite dissolution experiments at 450 °C and 600 bar, in which Fe speciation is also dominated by ferrous chloride species, yield Δ⁵⁷Fe_{fluid-hematite} values close to zero within errors, 0.15 ± 0.17‰ (2SE). In most experiments, chemical, redox, and isotopic equilibrium was attained, as shown by constancy over time of total dissolved Fe concentrations, aqueous Fe^II and Fe^III fractions, and Fe isotope ratios in solution, and identical Δ⁵⁷Fe values from dissolution and precipitation runs. Our measured equilibrium Δ⁵⁷Fe_{fluid-hematite} values at different temperatures, fluid compositions and iron redox state are within the range of fractionations in the system fluid-hematite estimated using reported theoretical β-factors for hematite and aqueous Fe species and the distribution of Fe aqueous complexes in solution. These theoretical predictions are however affected by large discrepancies among different studies, typically ±1‰ for the Δ⁵⁷Fe _{Fe(aq)-hematite} value at 200 °C. Our data may thus help to refine theoretical models for β-factors of aqueous iron species. This study provides the first experimental calibration of Fe isotope fractionation in the system hematite-saline aqueous fluid at elevated temperatures; it demonstrates the importance of redox control on Fe isotope fractionation at hydrothermal conditions.     

Experimental studies of equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Here we compare new experimental studies with theoretical predictions of equilibrium iron isotopic fractionation among aqueous ferric chloride complexes (Fe(H₂O)₆³⁺, FeCl(H₂O)₅²⁺, FeCl₂(H₂O)₄⁺, FeCl₃ (H₂O)₃, and FeCl₄^-), using the Fe-Cl-H₂O system as a simple, easily-modeled example of the larger variety of iron-ligand compounds, such as chlorides, sulfides, simple organic acids, and siderophores. Isotopic fractionation (⁵⁶Fe/⁵⁴Fe) among naturally occuring iron-bearing species at Earth surface temperatures (up to ∼3‰) is usually attributed to redox effects in the environment. However, theoretical modeling of reduced isotopic partition functions among iron-bearing species in solution also predicts fractionations of similar magnitude due to non-redox changes in speciation (i.e., ligand bond strength and coordination number). In the present study, fractionations are measured in a series of low pH ([H⁺] = 5 M) solutions of ferric chloride (total Fe = 0.0749 mol/L) at chlorinities ranging from 0.5 to 5.0 mol/L. Advantage is taken of the unique solubility of FeCl₄^- in immiscible diethyl ether to create a separate spectator phase, used to monitor changing fractionation in the aqueous solution. Δ⁵⁶Fe_aq-eth = δ⁵⁶Fe (total Fe remaining in aqueous phase)−δ⁵⁶Fe (FeCl₄^- in ether phase) is determined for each solution via MC-ICPMS analysis.Both experiments and theoretical calculations of Δ⁵⁶Fe_aq-eth show a downward trend with increasing chlorinity: Δ⁵⁶Fe_aq-eth is greatest at low chlorinity, where FeCl₂(H₂O)₄⁺ is the dominant species, and smallest at high chlorinity where FeCl₃(H₂O)₃ is dominant. The experimental Δ⁵⁶Fe_aq-eth ranges from 0.8‰ at [Cl^-] = 0.5 M to 0.0‰ at [Cl^-] = 5.0 M, a decrease in aqueous-ether fractionation of 0.8‰. This is very close to the theoretically predicted decreases in Δ⁵⁶Fe_aq-eth, which range from 1.0 to 0.7‰, depending on the ab initio model.The rate of isotopic exchange and attainment of equilibrium are shown using spiked reversal experiments in conjunction with the two-phase aqueous-ether system. Equilibrium under the experimental conditions is established within 30 min.The general agreement between theoretical predictions and experimental results points to substantial equilibrium isotopic fractionation among aqueous ferric chloride complexes and a decrease in ⁵⁶Fe/⁵⁴Fe as the Cl^-/Fe³⁺ ion ratio increases. The effects on isotopic fractionation shown by the modeling of this simple iron-ligand system imply that ligands present in an aqueous environment are potentially important drivers of fractionation, are indicative of possible fractionation effects due to other speciation effects (such as iron-sulfide systems or iron bonding with organic ligands), and must be considered when interpreting iron isotope fractionation in the geological record.     

Lithium isotope systematics in a forested granitic catchment (Strengbach, Vosges Mountains, France)

Over the last decade it has become apparent that Li isotopes may be a good proxy to trace silicate weathering. However, the exact mechanisms which drive the behaviour of Li isotopes in surface environments are not totally understood and there is a need to better calibrate and characterize this proxy. In this study, we analysed the Li concentrations and isotopic compositions in the various surface reservoirs (soils, rocks, waters and plants) of a small forested granitic catchment located in the Vosges Mountains (Strengbach catchment, France, OHGE http://ohge.u-strasbg.fr). Li fluxes were calculated in both soil profiles and at the basin scale and it was found that even in this forested basin, atmospheric inputs and litter fall represented a minor flux compared to input derived from the weathering of rocks and soil minerals (which together represent a minimum of 70% of dissolved Li). Li isotope ratios in soil pore waters show large depth dependent variations. Average dissolved δ⁷Li decreases from −1.1‰ to −14.4‰ between 0 and −30 cm, but is +30.7‰ at −60 cm. This range of Li isotopic compositions is very large and it encompasses almost the entire range of terrestrial Li isotope compositions that have been previously reported. We interpret these variations to result from both the dissolution and precipitation of secondary phases. Large isotopic variations were also measured in the springs and stream waters, with δ⁷Li varying from +5.3‰ to +19.6‰. δ⁷Li increases from the top to the bottom of the basin and also covaries with discharge at the outlet. These variations are interpreted to reflect isotopic fractionations occurring during secondary phase precipitation along the water pathway through the rocks. We suggest that the dissolved δ⁷Li increases with increasing residence time of waters through the rocks, and so with increasing time of interaction between waters and solids. A dissolution precipitation model was used to fit the dissolved Li isotopic compositions. It was found that the isotopic compositions of springs and stream waters are explicable by an isotopic fractionation of −5‰ to −14‰ (best fit −10.8‰), in agreement with Li incorporation into clay. In soil solutions, it was found that isotopic fractionation during secondary precipitation is larger (at least −23‰), suggesting a major role for different secondary phases, such as iron oxides that maybe incorporate Li with a higher isotope fractionation.