Experimental determination of the role of diffusion on Li isotope fractionation during basaltic glass weathering

In order to use lithium isotopes as tracers of silicate weathering, it is of primary importance to determine the processes responsible for Li isotope fractionation and to constrain the isotope fractionation factors caused by each process as a function of environmental parameters (e.g. temperature, pH). The aim of this study is to assess Li isotope fractionation during the dissolution of basalt and particularly during leaching of Li into solution by diffusion or ion exchange. To this end, we performed dissolution experiments on a Li-enriched synthetic basaltic glass at low ratios of mineral surface area/volume of solution (S/V), over short timescales, at various temperatures (50 and 90 °C) and pH (3, 7, and 10). Analyses of the Li isotope composition of the resulting solutions show that the leachates are enriched in ⁶Li (δ⁷Li = +4.9 to +10.5‰) compared to the fresh basaltic glass (δ⁷Li = +10.3 ± 0.4‰). The δ⁷Li value of the leachate is lower during the early stages of the leaching process, increasing to values close to the fresh basaltic glass as leaching progresses. These low δ⁷Li values can be explained in terms of diffusion-driven isotope fractionation. In order to quantify the fractionation caused by diffusion, we have developed a model that couples Li diffusion with dissolution of the glassy silicate network. This model calculates the ratio of the diffusion coefficients of both isotopes (a = D₇/D₆), as well as its dependence on temperature, pH, and S/V. a is mainly dependent on temperature, which can be explained by a small difference in activation energy (0.10 ± 0.02 kJ/mol) between ⁶Li⁺ and ⁷Li⁺. This temperature dependence reveals that Li isotope fractionation during diffusion is low at low temperatures (T < 20 °C), but can be significant at high temperatures. However, concerning hydrothermal fluids (T > 120 °C), the dissolution rate of basaltic glass is also high and masks the effects of diffusion. These results indicate that the high δ⁷Li values of river waters, in particular in basaltic catchments, and the fractionated values of hydrothermal fluids are mainly controlled by precipitation of secondary phases.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.     

Li isotope fractionation in peridotites and mafic melts

We have measured the Li isotope ratios of a range of co-existing phases from peridotites and mafic magmas to investigate high-temperature fractionations of ⁷Li/⁶Li. The Li isotopic compositions of seven mantle peridotites, reconstructed from analyses of mineral separates, show little variation (δ⁷Li 3.2-4.9‰) despite a wide range in fertility and radiogenic isotopic compositions. The most fertile samples yield a best estimate of δ⁷Li ∼ 3.5‰ for the upper mantle. Bulk analyses of olivine separates from the xenoliths are typically ∼1.5‰ isotopically lighter than co-existing orthopyroxenes, suggestive of a small, high-temperature equilibrium isotope fractionation. On the other hand, bulk analyses of olivine phenocrysts and their host melts are isotopically indistinguishable. Given these observations, equilibrium mantle melting should generate melts with δ⁷Li little different from their sources (<0.5‰ lighter). In contrast to olivine and orthopyroxene, that dominate peridotite Li budgets, bulk clinopyroxene analyses are highly variable (δ⁷Li = 6.6‰ to −8.1‰). Phlogopite separated from a modally metasomatised xenolith yielded an extreme δ⁷Li of −18.9‰. Such large Li isotope variability is indicative of isotopic disequilibrium. This inference is strongly reinforced by in situ, secondary ion mass-spectrometry analyses which show Li isotope zonation in peridotite minerals. The simplest zoning patterns show isotopically light rims. This style of zoning is also observed in the phenocrysts of holocrystalline Hawaiian lavas. More dramatically, a single orthopyroxene crystal from a San Carlos xenolith shows a W-shaped Li isotope profile with a 40‰ range in δ⁷Li, close to the isotope variability seen in all terrestrial whole rock analyses. We attribute Li isotope zonation in mineral phases to diffusive fractionation of Li isotopes, within mineral phases and along melt pathways that pervade xenoliths. Given the high diffusivity of Li, the Li isotope profiles we observe can persist, at most, only a few years at magmatic temperatures. Our results thus highlight the potential of Li isotopes as a high-resolution geospeedometer of the final phases of magmatic activity and cooling.     

Oxygen isotope fractionation between synthetic aragonite and water: Influence of temperature and Mg concentration

Aragonite was precipitated in the laboratory at 0, 5, 10, 25, and 40 °C to determine the temperature dependence of the equilibrium oxygen isotope fractionation between aragonite and water. Forced CO₂ degassing, passive CO₂ degassing, and constant addition methods were employed to precipitate aragonite from supersaturated solutions, but the resulting aragonite-water oxygen isotope fractionation was independent of the precipitation method. In addition, under the experimental conditions of this study, the effect of precipitation rate on the oxygen isotope fractionation between aragonite and water was almost within the analytical error of ±∼0.13‰ and thus insignificant. Because the presence of Mg²⁺ ions is required to nucleate and precipitate aragonite from Na-Ca-Cl-HCO₃ solutions under these experimental conditions, the influence of the total Mg²⁺ concentration (up to ∼0.9 molal) on the aragonite-water oxygen isotope fractionation was examined at 25 °C. No significant Mg²⁺ ion effect, or oxygen isotope salt effect, was detected up to 100 mmolal total Mg²⁺ but a noticeable isotope salt effect was observed at ∼0.9 molal total Mg²⁺.On the basis of results of the laboratory synthesis experiments, a new expression for the aragonite-water fractionation is proposed over the temperature range of 0-40 °C:

1000lnα_{aragonite-water}=17.88±0.13(10³/T)-31.14±0.46     

Carbon and hydrogen isotope fractionation associated with the aerobic microbial oxidation of methane, ethane, propane and butane

Carbon isotope fractionation factors associated with the aerobic consumption of methane (C1), ethane (C2), propane (C3), and n-butane (C4) were determined from incubations of marine sediment collected from the Coal Oil Point hydrocarbon seep field, located offshore Santa Barbara, CA. Hydrogen isotope fractionation factors for C1, C2 and C3 were determined concurrently. Fresh sediment samples from two seep areas were each slurried with sea water and treated with C1, C2, C3 or C4, or with mixtures of all four gases. Triplicate samples were incubated aerobically at 15 °C, and the stable isotope composition and headspace levels of C1-C4 were monitored over the course of the experiment. Oxidation was observed for all C1-C4 gases, with an apparent preference for C3 and C4 over C1 and C2 in the mixed-gas treatments. Fractionation factors were calculated using a Rayleigh model by comparing the δ¹³C and δD of the residual C1-C4 gases to their headspace levels. Carbon isotope fractionation factors (reported in ε or (α-1) × 1000 notation) were consistent between seep areas and were −26.5‰ ± 3.9 for C1, −8.0‰ ± 1.7 for C2, −4.8‰ ± 0.9 for C3 and −2.9‰ ± 0.9 for C4. Fractionation factors determined from mixed gas incubations were similar to those determined from individual gas incubations, though greater variability was observed during C1 consumption. In the case of C1 and C3 consumption, carbon isotope fractionation appears to decrease as substrate becomes limiting. Hydrogen isotope fractionation factors determined from the two seep areas differed for C1 oxidation but were similar for C2 and C3. Hydrogen isotope fractionation factors ranged from −319.9‰ to −156.4‰ for C1 incubations, and averaged −61.9‰ ± 8.3 for C2 incubations and −15.1‰ ± 1.9 for C3 incubations. The fractionation factors presented here may be applied to estimate the extent of C1-C4 oxidation in natural gas samples, and should prove useful in further studying the microbial oxidation of these compounds in the natural environment.     

Lithium isotopes in island arc geothermal systems: Guadeloupe, Martinique (French West Indies) and experimental approach

We report lithium (Li) isotopic measurements in seawater-derived waters that were discharged from geothermal wells, thermal springs, and sub-marine springs located in volcanic island arc areas in Guadeloupe (the Bouillante geothermal field) and Martinique (Lamentin plain and the Diamant areas). While Li isotopic signatures of the geothermal fluids collected from deep reservoirs were found to be homogeneous for a given site, the δ⁷Li signatures for each of these reservoirs were significantly different. The first low temperature (25-250 °C) experiments of Li isotope exchange during seawater/basalt interaction confirmed that Li isotopic exchange is strongly temperature dependent, as previously inferred from natural studies. Li isotopic fractionation ranged from +19.4‰ (Δ_{solution-solid}) at 25 °C to +6.7‰ at 250 °C. These experiments demonstrated the importance of Li isotopic fractionation during the formation of Li-bearing secondary minerals and allowed us to determine the following empirical relationship between isotopic fractionation and temperature: Δ_{solution-solid} = 7847/T − 8.093. Application of experimental results and literature data to the Bouillante area suggested that geothermal water was in equilibrium at 250-260 °C. It likely has a deep and large reservoir located in the upper sheeted dike complex of the oceanic crust, just below the transition zone between andesite volcanic flows and the basaltic dikes. The upper dike section, from which Li is extracted by hydrothermal fluids, was characterized by light Li isotopic values in the rocks, indicating retention of ⁶Li by the altered rocks. For the Lamentin and Diamant areas, the geothermal fluids appeared to be in equilibrium with reservoir volcano-sedimentary rocks at 90-120 °C and 180 °C, respectively. Further evidence for this argument is provided by the fact that only the Na/Li thermometric relationship determined for sedimentary basins yielded temperature values in agreement with those measured or estimated for the reservoir fluids. This suggests the importance of a sedimentary signature in these reservoir rocks. Altogether, this study highlights that the use of Li isotopic systematics is a powerful tool for characterizing the origin of geothermal waters as well as the nature of their reservoir rocks.     

Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths

We present whole rock Li and Mg isotope analyses of 33 ultramafic xenoliths from the terrestrial mantle, which we compare with analyses of 30 (mostly chondritic) meteorites. The accuracy of our new Mg isotope ratio measurement protocol is substantiated by a combination of standard addition experiments, the absence of mass independent effects in terrestrial samples and our obtaining identical values for rock standards using two different separation chemistries and three different mass-spectrometric introduction systems. Carbonaceous, ordinary and enstatite chondrites have irresolvable mean stable Mg isotopic compositions (δ²⁵Mg = −0.14 ± 0.06; δ²⁶Mg = −0.27 ± 0.12‰, 2SD), but our enstatite chondrite samples have lighter δ⁷Li (by up to ∼3‰) than our mean carbonaceous and ordinary chondrites (3.0 ± 1.5‰, 2SD), possibly as a result of spallation in the early solar system. Measurements of equilibrated, fertile peridotites give mean values of δ⁷Li = 3.5 ± 0.5‰, δ²⁵Mg = −0.10 ± 0.03‰ and δ²⁶Mg = −0.21 ± 0.07‰. We believe these values provide a useful estimate of the primitive mantle and they are within error of our average of bulk carbonaceous and ordinary chondrites. A fuller range of fresh, terrestrial, ultramafic samples, covering a variety of geological histories, show a broad positive correlation between bulk δ⁷Li and δ²⁶Mg, which vary from −3.7‰ to +14.5‰, and −0.36‰ to + 0.06‰, respectively. Values of δ⁷Li and δ²⁶Mg lower than our estimate of primitive mantle are strongly linked to kinetic isotope fractionation, occurring during transport of the mantle xenoliths. We suggest Mg and Li diffusion into the xenoliths is coupled to H loss from nominally anhydrous minerals following degassing. Diffusion models suggest that the co-variation of Mg and Li isotopes requires comparable diffusivities of Li and Mg in olivine. The isotopically lightest samples require ∼5-10 years of diffusive ingress, which we interpret as a time since volatile loss in the host magma. Xenoliths erupted in pyroclastic flows appear to have retained their mantle isotope ratios, likely as a result of little prior degassing in these explosive events. High δ⁷Li, coupled with high [Li], in rapidly cooled arc peridotites may indicate that these samples represent fragments of mantle wedge that has been metasomatised by heavy, slab-derived fluids. If such material is typically stirred back into the convecting mantle, it may account for the heavy δ⁷Li seen in some oceanic basalts.     

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.     

Sources and biological fractionation of Silicon isotopes in the Eastern Equatorial Pacific

Silicon isotopes in dissolved silicic acid were measured in the upper four kilometers between 4°N and 3°S latitude at 110°W longitude in the eastern Equatorial Pacific. Silicon isotopes became progressively heavier with silicic acid depletion of surface water as expected from biological fractionation. The value of ε estimated by applying a steady-state isotope fractionation model to data from all stations between 4°N and 3°S was −0.77 ± 0.12‰ (std. err.). When the analysis was restricted to those stations whose temperature and salinity profiles indicated that they were directly influenced by upwelling of the Equatorial Undercurrent (EUC), the resulting value of ε was −1.08 ± 0.27‰ (std. err.) similar to the value established in culture studies (−1.1‰). When the non steady state Rayleigh model was applied to the same restricted data set the resulting value of ε was significantly more positive, −0.61 ± 0.16‰ (std. err.). To the extent that the equatorial system approximates a steady state these results support a value of −1.1‰ for the fractionation factor for isotopes of Si in the sea. Without the assumption of steady state the value of ε can only be constrained to be between −0.6 and −1.1‰. Silicic acid in Equatorial Pacific Deep Water below 2000 m had a near constant δ³⁰Si of +1.32 ± 0.05‰. That value is significantly more positive than obtained for North Pacific Deep Water at similar depths at stations to the northwest of our study area (0.9-1.0‰) and it is slightly less positive than new measures of the δ³⁰Si of silicic acid from the silicic acid plume centered over the Cascadia basin in the Northeast Pacific (Si(OH)₄ > 180  μM, δ³⁰Si = +1.46 ± 0.12‰ (SD, n = 4). We show that the data from the equator and Cascadia basin fit a general trend of increasing δ³⁰Si(OH)₄ with increasing silicic acid concentration in the deep sea, but that the isotope values from the Northeast Pacific are anomalously light. The observed level of variation in the silicon isotope composition of deep waters from this single ocean basin is considerably larger than that predicted by current models based on fractionation during opal formation with no isotope effect during dissolution. Confirmation of such high variability in deep water δ³⁰Si(OH)₄ within individual ocean basins will require reassessment of the mechanisms controlling the distribution of isotopes of silicon in the sea.     

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.     

Determination of carbon fractionation factor between aqueous carbonate and CO2(g) in two-direction isotope equilibration

The experiments were conducted in the open CO₂ system to find out the equilibrium fractionation between the carbonate ion and CO₂(g). The existence of isotopic equilibrium was checked using the two-direction approach by passing the CO₂−N₂ gases with different δ¹³C compositions (− 1.5‰ and − 23‰) through the carbonate solution with δ¹³C = − 4.2‰. The Δ_{CO3T²⁻−CO2(g)} equilibrium fractionation is given as 6.03 ± 0.17‰ at 25 °C. Discussion is provided about the significance of carbonate complexing in determination of Δ_{CO3T²⁻−CO2(g)} and Δ_{HCO3T⁻−CO2(g)} fractionations. Finally, an isotope numerical model of flow and kinetics of hydration and dehydroxylation is built to predict the isotopic behaviour of the system with time.     

Calcium isotope fractionation in calcite and aragonite

Calcium isotope fractionation was measured on skeletal aragonite and calcite from different marine biota and on inorganic calcite. Precipitation temperatures ranged from 0 to 28°C. Calcium isotope fractionation shows a temperature dependence in accordance with previous observations: 1000 · ln(α_cc) = −1.4 + 0.021 · T (°C) for calcite and 1000 · ln(α_ar) = −1.9 + 0.017 · T (°C) for aragonite. Within uncertainty the temperature slopes are identical for the two polymorphs. However, at all temperatures calcium isotopes are more fractionated in aragonite than in calcite. The offset in δ^44/40Ca is about 0.6‰. The underlying mechanism for this offset may be related to the different coordination numbers and bond strengths of the calcium ions in calcite and aragonite crystals, or to different Ca reaction behavior at the solid-liquid interface. Recently, the observed temperature dependence of the Ca isotope fractionation was explained quantitatively by the temperature control on precipitation rates of calcium carbonates in an experimental setting (Lemarchand et al., 2004). We show that this mechanism can in principle also be applied to CaCO₃ precipitation in natural environments in normal marine settings. Following this model, Ca isotope fractionation in marine Ca carbonates is primarily controlled by precipitation rates. On the other hand the larger Ca isotope fractionation of aragonite compared to calcite can not be explained by different precipitation rates. The rate control model of Ca isotope fractionation predicts a strong dependence of the Ca isotopic composition of carbonates on ambient CO₃²⁻ concentration. While this model is in general accordance with our observations in marine carbonates, cultured specimens of the planktic foraminifer Orbulina universa show no dependence of Ca-isotope fractionation on the ambient CO₃²⁻ concentration. The latter observation implies that the carbonate chemistry in the calcifying vesicles of the foraminifer is independent from the ambient carbonate ion concentration of the surrounding water.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

Magnesium isotope fractionation in silicate melts by chemical and thermal diffusion

Two types of laboratory experiments were used to quantify magnesium isotopic fractionations associated with chemical and thermal (Soret) diffusion in silicate liquids. Chemical diffusion couples juxtaposing a molten natural basalt (SUNY MORB) and a molten natural rhyolite (Lake County Obsidian) were run in a piston cylinder apparatus and used to determine the isotopic fractionation of magnesium as it diffused from molten basalt to molten rhyolite. The thermal diffusion experiments were also run in a piston cylinder apparatus but with a sample made entirely of molten SUNY MORB displaced from the hotspot of the assembly furnace so that the sample would have a temperature difference of about 100-200 °C from one end to the other. The chemical diffusion experiments showed fractionations of ²⁶Mg/²⁴Mg by as much as 7‰, which resulted in an estimate for the mass dependence of the self-diffusion coefficients of the magnesium isotopes corresponding to D_26Mg/D_24Mg=(24/26)^β with β = 0.05. The thermal diffusion experiments showed that a temperature difference of about 100 °C resulted in the MgO, CaO, and FeO components of the basalt becoming slightly enriched by about 1 wt% in the colder end while SiO₂ was enriched by several wt% in the hotter end. The temperature gradient also fractionated the magnesium isotopes. A temperature difference of about 150 °C produced an 8‰ enrichment of ²⁶Mg/²⁴Mg at the colder end relative to the hotter end. The magnesium isotopic fractionation as a function of temperature in molten basalt corresponds to 3.6 × 10⁻²‰/°C/amu.     

Oxygen isotope fractionation of dissolved oxygen during reduction by ferrous iron

The oxygen isotope fractionation factor of dissolved oxygen gas has been measured during inorganic reduction by aqueous FeSO₄ at 10−54 °C under neutral (pH 7) and acidic (pH 2) conditions, with Fe(II) concentrations ranging up to 0.67 mol L⁻¹, in order to better understand the geochemical behavior of oxygen in ferrous iron-rich groundwater and acidic mine pit lakes. The rate of oxygen reduction increased with increasing temperature and increasing Fe(II) concentration, with the pseudo-first-order rate constant k ranging from 2.3 to 82.9 × 10⁻⁶ s⁻¹ under neutral conditions and 2.1 to 37.4 × 10⁻⁷ s⁻¹ under acidic conditions. The activation energy of oxygen reduction was 30.9 ± 6.6 kJ mol⁻¹ and 49.7 ± 13.0 kJ mol⁻¹ under neutral and acidic conditions, respectively. Oxygen isotope enrichment factors (ε) become smaller with increasing temperature, increasing ferrous iron concentration, and increasing reaction rate under acidic conditions, with ε values ranging from −4.5‰ to −11.6‰. Under neutral conditions, ε does not show any systematic trends vs. temperature or ferrous iron concentration, with ε values ranging from −7.3 to −10.3‰. Characterization of the oxygen isotope fractionation factor associated with O₂ reduction by Fe(II) will have application to elucidating the process or processes responsible for oxygen consumption in environments such as groundwater and acidic mine pit lakes, where a number of possible processes (e.g. biological respiration, reduction by reduced species) may have taken place.     

Stable isotope fractionation during bacterial sulfate reduction is controlled by reoxidation of intermediates

Bacterial sulfate reduction is one of the most important respiration processes in anoxic habitats and is often assessed by analyzing the results of stable isotope fractionation. However, stable isotope fractionation is supposed to be influenced by the reduction rate and other parameters, such as temperature. We studied here the mechanistic basics of observed differences in stable isotope fractionation during bacterial sulfate reduction. Batch experiments with four sulfate-reducing strains (Desulfovibrio desulfuricans, Desulfobacca acetoxidans, Desulfonatronovibrio hydrogenovorans, and strain TRM1) were performed. These microorganisms metabolize different carbon sources (lactate, acetate, formate, and toluene) and showed broad variations in their sulfur isotope enrichment factors. We performed a series of experiments on isotope exchange of ¹⁸O between residual sulfate and ambient water. Batch experiments were conducted with ¹⁸O-enriched (δ¹⁸O_water = +700‰) and depleted water (δ¹⁸O_water = −40‰), respectively, and the stable ¹⁸O isotope shift in the residual sulfate was followed. For Desulfovibrio desulfuricans and Desulfonatronovibrio hydrogenovorans, which are both characterized by low sulfur isotope fractionation (ε_S > −13.2‰), δ¹⁸O values in the remaining sulfate increased by only 50‰ during growth when ¹⁸O-enriched water was used for the growth medium. In contrast, with Desulfobacca acetoxidans and strain TRM1 (ε_S < −22.7‰) the residual sulfate showed an increase of the sulfate δ¹⁸O close to the values of the enriched water of +700‰. In the experiments with δ¹⁸O-depleted water, the oxygen isotope values in the residual sulfate stayed fairly constant for strains Desulfovibrio desulfuricans, Desulfobacca acetoxidans and Desulfonatronovibrio hydrogenovorans. However, strain TRM1, which exhibits the lowest sulfur isotope fractionation factor (ε_S < −38.7‰) showed slightly decreasing δ¹⁸O values.Our results give strong evidence that the oxygen atoms of sulfate exchange with water during sulfate reduction. However, this neither takes place in the sulfate itself nor during formation of APS (adenosine-5′-phosphosulfate), but rather in intermediates of the sulfate reduction pathway. These may in turn be partially reoxidized to form sulfate. This reoxidation leads to an incorporation of oxygen from water into the “recycled” sulfate changing the overall ¹⁸O isotopic composition of the remaining sulfate fraction. Our study shows that such incorporation of ¹⁸O is correlated with the stable isotope enrichment factor for sulfur measured during sulfate reduction. The reoxidation of intermediates of the sulfate reduction pathway does also strongly influence the sulfur stable isotope enrichment factor. This aforesaid reoxidation is probably dependent on the metabolic conversion of the substrate and therefore also influences the stable isotope fractionation factor indirectly in a rate dependent manner. However, this effect is only indirect. The sulfur isotope enrichment factors for the kinetic reactions themselves are probably not rate dependent.     

Lithium isotopic systematics of peridotite xenoliths from Hannuoba, North China Craton: Implications for melt-rock interaction in the considerably thinned lithospheric mantle

Li concentrations and isotopic compositions of coexisting minerals (ol, opx, and cpx) from peridotite xenoliths entrained in the Hannuoba Tertiary basalts, North China Craton, provide insight into Li isotopic fractionation between mantle minerals during melt-rock interaction in the considerably thinned lithospheric mantle. Bulk analyses of mineral separates show significant enrichment of Li in cpx (2.4-3.6 ppm) relative to olivine (1.2-1.8 ppm), indicating that these peridotites have been affected by mantle metasomatism with mafic silicate melts. Bulk olivine separates (δ⁷Li ∼ +3.3‰ to +6.4‰) are isotopically heavier than coexisting pyroxenes (δ⁷Li ∼ −3.3‰ to −8.2‰ in cpx, and −4.0‰ to −6.7‰ in opx). Such large variation suggests Li elemental and isotopic disequilibrium. This conclusion is supported by results from in situ SIMS analyses of mineral grains where significant Li elemental and isotopic zonations exist. The olivine and opx have lower Li concentrations and heavier Li isotopes in the rims than in the cores. This reverse correlation of δ⁷Li with Li concentrations indicates diffusive fractionation of Li isotopes. However, the zoning patterns in coexisting cpx show isotopically heavier rims with higher Li abundances. This positive correlation between δ⁷Li and Li concentrations suggests a melt mixing trend. We attribute Li concentration and isotope zonation in minerals to the effects of two-stage diffusive fractionation coupled with melt-rock interaction. The earliest melts may have been derived from the subducted oceanic slab with low δ⁷Li values produced by isotopic fractionation during the dehydration of the seawater-altered slab. Melts at later stages were derived from the asthenosphere and interacted with the peridotites, producing the Li elemental and isotopic zoning in mineral grains. These data thus provide evidence for multiple-stage peridotite-melt interaction in the lithospheric mantle beneath the northern North China Craton.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

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.     

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.