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.     

Kinetic isotope effect during reduction of iron from a silicate melt

Iron isotopic compositions measured in chondrules from various chondrites vary between δ⁵⁷Fe/⁵⁴Fe = +0.9‰ and −2.0‰, a larger range than for igneous rocks. Whether these compositions were inherited from chondrule precursors, resulted from the chondrule-forming process itself or were produced by later parent body alteration is as yet unclear. Since iron metal is a common phase in some chondrules, it is important to explore a possible link between the metal formation process and the observed iron isotope mass fractionation. In this experimental study we have heated a fayalite-rich composition under reducing conditions for heating times ranging from 2 min to 6 h. We performed chemical and iron isotope analyses of the product phases, iron metal and silicate glass. We demonstrated a lack of evaporation of Fe from the silicate melt in similar isothermal experiments performed under non-reducing conditions. Therefore, the measured isotopic mass fractionation in the glass, ranging between −0.32‰ and +3.0‰, is attributed to the reduction process. It is explained by the faster transport of lighter iron isotopes to the surface where reduction occurs, and is analogous to kinetic isotope fractionation observed in diffusion couples [Richter, F.M., Davis, A.M., Depaolo, D.J., Watson, E.B., 2003. Isotope fractionation by chemical diffusion between molten basalt and rhyolite. Geochim. Cosmochim. Acta67, 3905-3923]. The metal phase contains 90-99.8% of the Fe in the system and lacks significant isotopic mass fractionation, with values remaining similar to that of the starting material throughout. The maximum iron isotope mass fractionation in the glass was achieved within 1 h and was followed by an isotopic exchange and re-equilibration with the metal phase (incomplete at ∼6 h). This study demonstrates that reduction of silicates at high temperatures can trigger iron isotopic fractionation comparable in its bulk range to that observed in chondrules. Furthermore, if metal in Type I chondrules was formed by reduction of Fe silicate, our observed isotopic fractionations constrain chondrule formation times to approximately 60 min, consistent with previous work.     

Liquid composition-dependence of calcium isotope fractionation during diffusion in molten silicates

Liquid phase diffusion experiments were carried out to determine whether diffusive isotopic fractionation of a major chemical element (Ca) varies with chemical composition in high-temperature molten silicates. The objective was to determine how differences in silicate liquid structure, such as the ratio of bridging to non-bridging oxygen atoms, as well as bulk transport properties such as viscosity, relate to isotope discrimination during diffusion. This information, in turn, may relate to the lifetimes and sizes of multi-atom structures in the liquid. Diffusion couples consisting of juxtaposed natural mafic and felsic liquids were held at T = 1450 °C and P = 1.0 GPa for durations of 12-24 h in a standard piston-cylinder assembly. Experiments were done using different mafic endmember compositions (two tholeiitic basalts and a ugandite) and a single rhyolite composition. Major-element diffusion profiles and Ca isotope profiles were measured on the recovered quenched glasses. The starting materials were isotopically indistinguishable, but ⁴⁴Ca/⁴⁰Ca variations of ca. 5‰ arose due to a mass dependence of the Ca diffusion coefficients. Results indicate that the mass dependence of Ca diffusion coefficients varies with the magnitude and direction of aluminum gradients and the viscosity of the liquid. Some Ca fractionations result mainly from Al gradients.A simplified multicomponent diffusion model was used to model the experimental results. The model allows for diffusion of Ca in response to gradients in the concentrations of both CaO as well as Al₂O₃, and the model results are consistent with the inferred existence of at least two distinct species of Ca. The magnitude of isotopic discrimination during diffusion also appears to be stronger on the rhyolite versus the basalt/ugandite side of diffusion couples. The results can largely be accounted for by an adaptation of the model of Dingwell (1990), whereby in high silica liquids, Ca diffuses largely by site hopping through a quasi-stationary aluminosilicate matrix, producing strong isotopic effects because the Ca diffusion is not strongly correlated with the movement of the framework atoms. In low-silica liquids, Ca diffusion is correlated with the movement of the other components and there is less mass discrimination. Combining our Ca results with Ca, Mg, and Li data from previous studies, we show that this model can explain most of the cation- and composition-dependence of diffusive isotopic fractionations observed thus far. A key parameter controlling isotopic discrimination is the ratio of the elemental (Ca, Mg, Li) diffusivity to the Eyring (or Si) diffusivity. However, all experiments done so far also exhibit isotopic features that are not yet fully explained; some of these may relate to small temperature gradients in the capsules, or to more complex coupling effects that are not captured in simplified diffusion models.     

Isotopic fractionation of the major elements of molten basalt by chemical and thermal diffusion

Samples produced in piston cylinder experiments were used to document the thermal isotopic fractionation of all the major elements of basalt except for aluminum and the fractionation of iron isotopes by chemical diffusion between a natural basalt and rhyolite. The thermal isotopic fractionations are summarized in terms of a parameter Ω_i defined as the fractionation in per mil per 100 °C per atomic mass units difference between the isotopes. For molten basalt we report Ω_Ca = 1.6, Ω_Fe = 1.1, Ω_Si = 0.6, Ω_O = 1.5. In an earlier paper we reported Ω_Mg = 3.6. These fractionations represent a steady state balance between thermal diffusion and chemical diffusion with the mass dependence of the thermal diffusion coefficient being significantly larger than the mass dependence of the chemical diffusion coefficients for isotopes of the same element. The iron isotopic measurements of the basalt-rhyolite diffusion couple showed significant fractionation that are parameterized in terms of a parameter β_Fe = 0.03 when the ratio of the diffusion coefficients D₅₄ and D₅₆ of ⁵⁴Fe and ⁵⁶Fe is expressed in terms of the atomic mass as D₅₄/D₅₆ = (56/54)^β_Fe. This value of β_Fe is smaller than what we had measured earlier for lithium, magnesium and calcium (i.e., β_Li = 0.215, β_Ca = 0.05, β_Mg = 0.05) but still significant when one takes into account the high precision with which iron isotopic compositions can be measured (i.e., ±0.03‰) and that iron isotope fractionations at magmatic temperatures from other causes are extremely small. In a closing section we discuss technological and geological applications of isotopic fractionations driven by either or both chemical and thermal gradients.     

Lithium isotopes in the system Qz-Ms-fluid:An experimental study

The fractionation of lithium isotopes among quartz, muscovite, and a chloride-bearing aqueous fluid has been investigated experimentally at 400°-500°C and 50-100 MPa. Experiments were performed for 15-60 days in cold seal vessels with natural mineral specimens. Lithium was introduced primarily through the fluid, which also contained KCl and HCl. In most runs, the fluid was prepared with the L-SVEC standard (δ⁷Li = 0) and was 1 M in total chloride with K/Li/H = 100/10/1. In two experiments, a ⁶Li spike was employed. The experiments demonstrate that quartz and muscovite are susceptible to pronounced, rapid shifts in Li isotopic composition by diffusion through interaction with a Li-bearing fluid, particularly at 500°C. At 500°C, fractionation factors were determined to be Δ_{quartz-aqueous} ≅ +8 to + 12‰ and Δ_{muscovite-aqueous} ≅ +18 to + 20‰. An intermineral fractionation factor is given by Δ_{muscovite-quartz} ≈ +9‰. At 400°C, the results suggest Δ_{quartz-aqueous} ≈ +4 to + 6‰. The study provides evidence of systematic fractionation in lithium isotopes at the temperatures of some magmatic processes, such as those associated with porphyry-type ore systems and pegmatites.     

Lithium isotopic disequilibrium of minerals in the spinel lherzolite xenoliths from Boeun,Korea

Spinel lherzolite xenoliths found in Boeun, Korea, have protogranular to porphyroclastic textures and are enclosed in a Miocene alkali basalt. The lithium concentration and isotopic compositions of olivine, clinopyroxene, and orthopyroxene separates from the spinel lherzolite, and whole rocks of the spinel lherzolites and alkali basalt were determined by inductively coupled plasma mass spectrometry (ICP-MS) and multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). The lithium concentrations of the olivines and orthopyroxenes range from 2.2 to 5.0 ppm and from 2.1 to 6.4 ppm, respectively. In contrast, the clinopyroxenes have larger lithium concentrations, from 2.0 to 8.4 ppm, which reflect their preferential lithium enrichment. The lithium isotopic compositions (δ⁷Li) of olivines (-5.4 to + 3.5‰), orthopyroxenes (-11.4 to -0.1‰), and clinopyroxenes (-14.4 to -4.7‰) range far beyond the average mantle composition of + 4 ± 2‰. The lithium isotopic composition of the host rock, alkali basalt (3.4‰), is within the range of the intraplate and oceanic island basalts. The spinel lherzolites from Boeun exhibits strong elemental and isotopic disequilibria due to the different lithium and lithium isotope diffusion velocities in the olivine, orthopyroxene, and clinopyroxene minerals after eruption and magma cooling.     

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.     

Lithium isotopic composition and concentration of the upper continental crust

The Li isotopic composition of the upper continental crust is estimated from the analyses of well-characterized shales, loess, granites and upper crustal composites (51 samples in total) from North America, China, Europe, Australia and New Zealand. Correlations between Li, δ⁷Li, and chemical weathering (as measured by the Chemical Index of Alteration (CIA)), and δ⁷Li and the clay content of shales (as measured by Al₂O₃/SiO₂), reflect uptake of heavy Li from the hydrosphere by clays. S-type granites from the Lachlan fold belt (-1.1 to -1.4‰) have δ⁷Li indistinguishable from their associated sedimentary rocks (-0.7 to 1.2‰), and show no variation in δ⁷Li throughout the differentiation sequence, suggesting that isotopic fractionation during crustal anatexis and subsequent differentiation is less than analytical uncertainty (±1‰, 2σ). The isotopically light compositions for both I- and S-type granites from the Lachlan fold belt (-2.5 to + 2.7 ‰) and loess from around the world (-3.1 to + 4.5‰) reflect the influence of weathering in their source regions. Collectively, these lithologies possess a limited range of Li isotopic compositions (δ⁷Li of −5‰ to + 5‰), with an average (δ⁷Li of 0 ± 2‰ at 1σ) that is representative of the average upper continental crust. Thus, the Li isotopic composition of the upper continental crust is lighter than the average upper mantle (δ⁷Li of + 4 ± 2‰), reflecting the influence of weathering on the upper crustal composition. The concentration of Li in the upper continental crust is estimated to be 35 ± 11 ppm (2σ), based on the average loess composition and correlations between insoluble elements (Ti, Nb, Ta, Ga and Al₂O₃, Th and HREE) and Li in shales. This value is somewhat higher than previous estimates (∼20 ppm), but is probably indistinguishable when uncertainties in the latter are accounted for.     

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.     

Temperature-dependent isotopic fractionation of lithium between clinopyroxene and high-pressure hydrous fluids

The fractionation of lithium isotopes between synthetic spodumene as representative of Li-bearing clinopyroxene and Cl- and OH-bearing aqueous fluids was experimentally determined between 500 and 900°C at 2.0 GPa. In all the experiments, ⁷Li was preferentially partitioned into the fluid. The fractionation is temperature dependent and approximated by the equation Δ⁷Li_{(clinopyroxene–fluid)}=−4.61×(1,000/T [K]) + 2.48; R ²=0.86. Significant Li isotopic fractionation of about 1.0‰ exists even at high temperatures of 900°C. Using neutral and weakly basic fluids revealed that the amount of fractionation is not different. The Li isotopic fractionation between altered basalt and hot spring water (350°C) in natural samples is in good agreement with our experimentally determined fractionation curve. The data confirm earlier speculations drawn from the Li isotopic record of dehydrated metamorphic rocks that fluids expelled from a dehydrating slab carry heavier Li into the mantle wedge, and that a light Li component is introduced into the deeper mantle. Li and Li isotopes are redistributed among wedge minerals as fluids travel across the wedge into hotter regions of arc magma production. This modifies the Li isotopic characteristics of slab-derived fluids erasing their source memory, and explains the absence of cross-arc variations of Li isotopes in arc basalts.     

Diffusion induced Li isotopic fractionation during the cooling of magmatic rocks: The case of pyroxene phenocrysts from nakhlite meteorites

Ion-microprobe was used to measure Li abundances and isotopic compositions in pyroxenes from three Martian meteorites belonging to the nakhlite family. The profiles performed across augite crystals from Northwest Africa 817 show a large isotopic zoning from crystal cores (δ⁷Li ∼ 0‰) to rims (δ⁷Li ∼ +20‰) while Li abundances are almost constant (∼9.2 μg/g). Unlike NWA 817, the pyroxene studied in the Miller Range 03346 nakhlite shows a zoning in Li abundance, with concentrations increasing from ∼2.5 μg/g in the core to ∼9 μg/g in the rim. The augite rim (δ⁷Li = +7‰) is slightly enriched in ⁷Li with regard to the core (δ⁷Li = +4‰), but most of the isotopic variations observed occur at an intermediate position along the profile, where δ⁷Li falls down to ∼−11‰. In the case of Nakhla, Li concentrations in augite increase from cores (∼3.5 μg/g) to rims (∼6.5 μg/g), while the δ⁷Li variation is restricted (i.e., between δ⁷Li = +6.0 and +12.6‰). For the three meteorites the Li abundances were also measured in the groundmass, which was found to be enriched in lithium (∼10 μg/g). Conventional magmatic and post-magmatic processes such as alteration and fractional crystallization, fail to explain the dataset obtained on nakhlites. Degassing processes, which were previously proposed to explain the Li distribution in shergottite crystals, cannot result in the strong decoupling between Li abundances and isotopic composition observed in nakhlites. We suggest that the original magmatic Li distributions (concentrations and isotopic compositions) in nakhlites have been modified by diffusion of Li from the Li-rich groundmass towards the pyroxene crystals during sub-solidus cooling. Diffusion appears to have been efficient for NWA 817 and MIL 03346 but, apparently, did not produce a significant migration of Li in Nakhla, possibly because of the lower abundance of groundmass in the latter. Diffusion induced Li redistributions may also affect terrestrial porphyric rocks but very specific cooling rates are required to quench the diffusion profiles as observed in two of the present nakhlites.     

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.     

Carbon isotopic fractionation during a mesocosm bloom experiment dominated by Emiliania huxleyi: Effects of CO2 concentration and primary production

We investigated the effect of CO₂ and primary production on the carbon isotopic fractionation of alkenones and particulate organic matter (POC) during a natural phytoplankton bloom dominated by the coccolithophore Emiliania huxleyi. In nine semi-closed mesocosms (∼11 m³ each), three different CO₂ partial pressures (pCO₂) in triplicate represented glacial (∼180 ppmv CO₂), present (∼380 ppmv CO₂), and year 2100 (∼710 ppmv CO₂) CO₂ conditions. The largest shift in alkenone isotopic composition (4-5‰) occurred during the exponential growth phase, regardless of the CO₂ concentration in the respective treatment. Despite the difference of ∼500 ppmv, the influence of pCO₂ on isotopic fractionation was marginal (1-2‰). During the stationary phase, E. huxleyi continued to produce alkenones, accumulating cellular concentrations almost four times higher than those of exponentially dividing cells. Our isotope data indicate that, while alkenone production was maintained, the interaction of carbon source and cellular uptake dynamics by E. huxleyi reached a steady state. During stationary phase, we further observed a remarkable increase in the difference between δ¹³C of bulk organic matter and of alkenones spanning 7-12‰. We suggest that this phenomenon is caused mainly by a combination of extracellular release of ¹³C-enriched polysaccharides and subsequent particle aggregation induced by the production of transparent exopolymer particles (TEP).     

Li isotopic variations in single pyroxenes from the Northwest Africa 480 shergottite (NWA 480): a record of degassing of Martian magmas?

Lithium abundances and isotopic compositions were measured by ion microprobe in individual grains of pyroxene, and in a few maskelynites and Ca-phosphates grains, from the Martian meteorite Northwest Africa 480 (NWA 480).In pyroxenes Li abundances are nearly constant from core to rim with concentrations ranging between 3 and 4 μg/g. In contrast, a significant isotopic zoning is observed with δ⁷Li values increasing within single crystals from ∼ −17‰ in the core to ∼ +10‰ in the rim, most of the variability being observed in the core. Plagioclase (now maskelynite) and phosphate crystals, which co-crystallized with the pyroxene rims, display similar δ⁷Li values. Because of the incompatible behavior of Li, the present constancy of Li concentrations within zoned pyroxenes rules out any simple crystallization model in a closed system for Li. The large Li isotopic variations observed within pyroxenes support this conclusion. There is no evidence in support of secondary alteration of NWA 480 to explain the Li isotopic variations, which thus most likely reflect magmatic processes on Mars. Degassing might explain the Li systematics observed in NWA 480 pyroxenes. Because Li has a strong affinity with water-rich fluids, a significant loss of Li from NWA 480 parental melt can happen upon melt emplacement and cooling. Such a Li loss could compensate the effect of crystal fractionation and thus help to maintain constant the Li content of the melt. Li isotopic fractionation is anticipated to accompany this process, ⁷Li being depleted relative to ⁶Li in the volatile phase. The magnitude of the isotopic change of the fractionating melts is difficult to predict because it depends on the value of the Li isotopic fractionation and on the amount of Li loss, but at first glance it seems consistent with the increase of δ⁷Li values observed in NWA 480 pyroxenes with increasing fractionation. The present data suggest that degassing prevailed not only during the crystallization of shergottites like Zagami and Shergotty, but also during the crystallization of the other types of basaltic shergottites.     

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.     

Quantifying Li isotope fractionation during smectite formation and implications for the Li cycle

Tri-octahedral Li-Mg smectites (hectorites) were synthesized at temperatures ranging from 25 to 250 °C, in the presence of solutions highly enriched in lithium. After removing all the exchangeable lithium from the synthesized clays, Li isotope fractionation (Δ⁷Li_{clay-solution}) was determined. This fractionation was linked to Li incorporation into the structural octahedral site, substituting for Mg²⁺. As predicted, experimental Δ⁷Li_{clay-solution} inversely correlates with temperature, and ranges from −1.6‰ ± 1.3‰ at 250 °C to −10.0‰ ± 1.3‰ at 90 °C, and then stays relatively constant down to 25 °C. The relatively constant isotope fractionation factor below 90 °C may be due to high concentrations of edge octahedra in low crystallinity smectites. The isotopic fractionation factor (α), for a given temperature, does not depend on the solution matrix, nor on the amount of structural Li incorporated into the clay. Empirical linear laws for α as a function of 1/T (K) were inferred. Smectite Li contents and smectite-solution distribution coefficients (D_Li/Mg) increase with temperature, as expected for a substitution process. The fractions of dissolved Li incorporated into the smectite octahedral sites are small and do not depend on the duration of the experiment. In a seawater-like matrix solution, less Li is incorporated into the smectites, probably as a result of competition with dissolved Mg²⁺ ions for incorporation into the octahedral sites. The high Li contents observed in marine smectites are therefore best explained either by a significant contribution from basalts, by adsorption processes, or by the influence of seawater chemical composition on distribution coefficients. We also calculate, using present-day estimates of hydrothermal water and river fluxes, that a steady-state ocean would require a relatively large global clay-water Li isotope fractionation (−12‰ to −21‰). This study demonstrates the ability of laboratory experiments to quantify the impact of secondary phases on the Li geochemical cycle and associated isotope fractionations.     

Ion-exchange fractionation of copper and zinc isotopes

Whether transition element isotopes can be fractionated at equilibrium in nature is still uncertain. Standard solutions of Cu and Zn were eluted on an anion-exchange resin, and the isotopic compositions of Cu (with respect to Zn) of the eluted fractions were measured by multiple-collector inductively coupled plasma mass spectrometry. It was found that for pure Cu solutions, the elution curves are consistent with a ⁶³Cu/⁶⁵Cu mass fractionation coefficient of 0.46‰ in 7 mol/L HCl and 0.67‰ in 3 mol/L HCl between the resin and the solution. Batch fractionation experiments confirm that equilibrium fractionation of Cu between resin and 7 mol/L HCl is ∼0.4‰ and therefore indicates that there is no need to invoke kinetic fractionation during the elution. Zn isotope fractionation is an order of magnitude smaller, with a ⁶⁶Zn/⁶⁸Zn fractionation factor of 0.02‰ in 12 mol/L HCl. Cu isotope fractionation results determined from a chalcopyrite solution in 7 mol/L HCl give a fractionation factor of 0.58‰, which indicates that Fe may interfere with Cu fractionation.Comparison of Cu and Zn results suggests that the extent of Cu isotopic fractionation may signal the presence of so far unidentified polynuclear complexes in solution. In contrast, we see no compelling reason to ascribe isotope fractionation to the coexistence of different oxidation states. We further suggest that published evidence for iron isotopic fractionation in nature and in laboratory experiments may indicate the distortion of low-spin Fe tetrahedral complexes.The isotope geochemistry of transition elements may shed new light on their coordination chemistry. Their isotopic fractionation in the natural environment may be interpreted using models of thermodynamic fractionation.     

Kinetic isotopic fractionation during diffusion of ionic species in water

Experiments specifically designed to measure the ratio of the diffusivities of ions dissolved in water were used to determine . The measured ratio of the diffusion coefficients for Li and K in water (D_Li/D_K = 0.6) is in good agreement with published data, providing evidence that the experimental design being used resolves the relative mobility of ions with adequate precision to also be used for determining the fractionation of isotopes by diffusion in water. In the case of Li, we found measurable isotopic fractionation associated with the diffusion of dissolved LiCl (D_⁷Li/D_⁶Li=0.99772±0.00026). This difference in the diffusion coefficient of ⁷Li compared to ⁶Li is significantly less than that reported in an earlier study, a difference we attribute to the fact that in the earlier study Li diffused through a membrane separating the water reservoirs. Our experiments involving Mg diffusing in water found no measurable isotopic fractionation (D_²⁵Mg/D_²⁴Mg=1.00003±0.00006). Cl isotopes were fractionated during diffusion in water (D_³⁷Cl/D_³⁵Cl=0.99857±0.00080) whether or not the co-diffuser (Li or Mg) was isotopically fractionated. The isotopic fractionation associated with the diffusion of ions in water is much smaller than values we found previously for the isotopic fractionation of Li and Ca isotopes by diffusion in molten silicate liquids. A major distinction between water and silicate liquids is that water surrounds dissolved ions with hydration shells, which very likely play an important but still poorly understood role in limiting the isotopic fractionation associated with diffusion.     

Carbon and hydrogen isotopic compositions of products of open-system catalytic hydrogenation of CO2: Implications for abiogenic hydrocarbons in Earth’s crust

This paper reports the isotope effects in an open-system Fischer-Tropsch type (FTT) synthesis, with implications for the origin of natural abiogenic hydrocarbons. The starting form of carbon was CO₂, with carbon and hydrogen isotopic compositions measured for products of catalytic hydrogenation of CO₂ on iron and cobalt catalysts (FTCO₂-Fe and FTCO₂-Co) at 350 and 245 °C, respectively, and 10 MPa. The carbon isotopic composition of the resulting saturated hydrocarbons (alkanes) as a function of carbon number shows a positive trend for both FTCO₂-Fe and FTCO₂-Co, with a fractionation of 2-4‰ and 3-6‰ between CH₄ and C₂H₆ over the Fe and Co catalysts, respectively. The unsaturated hydrocarbons (alkenes) do not show any trend. A strong kinetic isotope fractionation (>40‰) occurred between CO₂ and CH₄ in both experiments. The hydrogen isotope fractionation between alkanes appeared to be similar to that found in natural (thermogenic and biogenic) gases, with enrichment in deuterium of longer hydrocarbon chains; the dominant H/D fractionation occurred between CH₄ and C₂H₆. Alkenes in the products of the FTCO₂-Fe reaction are enriched in deuterium (∼50‰) and do not show any trend versus carbon number. We suggest that other than FTT reactions or a simple mixing are responsible for the occurrence of the inverse isotopic trends in both δ¹³C and δD found in light hydrocarbons in some terrestrial environments and meteorites.     

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.