Chemical and oxygen isotopic compositions of accretionary rim and matrix olivine in CV chondrites: Constraints on the evolution of nebular dust

A correlation of petrography, mineral chemistry and in situ oxygen isotopic compositions of fine-grained olivine from the matrix and of fine- and coarse-grained olivine from accretionary rims around Ca-Al-rich inclusions (CAIs) and chondrules in CV chondrites is used here to constrain the processes that occurred in the solar nebula and on the CV parent asteroid. The accretionary rims around Leoville, Vigarano, and Allende CAIs exhibit a layered structure: the inner layer consists of coarse-grained, forsteritic and ¹⁶O-rich olivine (Fa_1-40 and Δ¹⁷O = −24‰ to −5‰; the higher values are always found in the outer part of the layer and only in the most porous meteorites), whereas the middle and the outer layers contain finer-grained olivines that are more fayalitic and ¹⁶O-depleted (Fa_15-50 and Δ¹⁷O = −18‰ to +1‰). The CV matrices and accretionary rims around chondrules have olivine grains of textures, chemical and isotopic compositions similar to those in the outer layers of accretionary rims around CAIs. There is a correlation between local sample porosity and olivine chemical and isotopic compositions: the more compact regions (the inner accretionary rim layer) have the most MgO- and ¹⁶O-rich compositions, whereas the more porous regions (outer rim layers around CAIs, accretionary rims around chondrules, and matrices) have the most MgO- and ¹⁶O-poor compositions. In addition, there is a negative correlation of olivine grain size with fayalite contents and Δ¹⁷O values. However, not all fine-grained olivines are FeO-rich and ¹⁶O-poor; some small (<1 μm in Leoville and 5-10 μm in Vigarano and Allende) ferrous (Fa_>20) olivine grains in the outer layers of the CAI accretionary rims and in the matrix show significant enrichments in ¹⁶O (Δ¹⁷O = −20‰ to −10‰). We infer that the inner layer of the accretionary rims around CAIs and, at least, some olivine grains in the finer portions of accretionary rims and CV matrices formed in an ¹⁶O-rich gaseous reservoir, probably in the CAI-forming region. Grains in the outer layers of the CAI accretionary rims and in the rims around chondrules as well as matrix may have also originated as ¹⁶O-rich olivine. However, these olivines must have exchanged O isotopes to variable extents in the presence of an ¹⁶O-poor reservoir, possibly the nebular gas in the chondrule-forming region(s) and/or fluids in the parent body. The observed trend in isotopic compositions may arise from mixtures of ¹⁶O-rich forsterites with grain overgrowths or newly formed grains of ¹⁶O-poor fayalitic olivines formed during parent body metamorphism. However, the observed correlations of chemical and isotopic compositions of olivine with grain size and local porosity of the host meteorite suggest that olivine accreted as a single population of ¹⁶O-rich forsterite and subsequently exchanged Fe-Mg and O isotopes in situ in the presence of aqueous solutions (i.e., fluid-assisted thermal metamorphism).     

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.     

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.     

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.     

Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa

Fourier transform infrared spectrometry (FTIR) analyses of olivines from peridotite xenoliths found in southern African kimberlites indicate 0 to 80 ppm H₂O concentrations. OH absorbance profiles across olivine grains show homogeneous H contents from core to edge for most samples. In one sample the olivines are H-free, while another has olivines characterized by lower H contents at the grain edges compared to the cores, indicating H loss during transport of the xenolith to the surface. Flat or near-flat H profiles place severe constraints on the duration of H loss from olivine grains, with implications for kimberlite magma ascent rates. Diffusion equations were used to estimate times of H loss of about 4 h for the sample with heterogeneous olivine H contents. Resulting kimberlite ascent rates are calculated to be 5-37 m s⁻¹ minimum, although these estimates are highly dependent on volatile contents and degassing behavior of the host kimberlite magma. Xenolithic olivines from alkali basalts generally have lower H contents and more pronounced H diffusion profiles than do those from kimberlites. This difference is likely caused by higher magma temperatures and lower ascent rates of alkali basalts compared to kimberlites.     

Lithium elemental and isotopic variations in rock-melt interaction

Despite the occurrence of highly variable lithium (Li) elemental distribution and isotopic fractionation in mantle mineral, the mechanism of Li heterogeneity and fractionation remains a controversial issue. We measured Li contents and isotopic compositions of olivine and clinopyroxene xenocrysts and phenocrysts from kamafugite host lavas, as well as minerals in melt pockets occurring as metasomatic products in peridotite xenoliths from the Western Qinling, central China. The olivine xenocrysts in the kamafugites show compositional zonation. The cores have high Mg# (100 × Mg/(Mg+Fe); 91.0–92.2) and Li abundances (5.63–21.7 ppm), low CaO contents (≤0.12 wt%) and low δ⁷Li values (−39.6 to −6.76‰), which overlap with the compositional ranges of the olivines in the melt pockets as well as those in peridotite xenoliths. The rims of the olivine xenocrysts display relatively low Mg# (85.9–88.2), high CaO contents (0.19–0.38 wt%) and high δ⁷Li values (18.3–26.9‰), which are comparable to the olivine phenocrysts (Mg#: 86.4–87.1; CaO: 0.20–0.28 wt%; Li: 12.4–36.8 ppm; δ⁷Li: 18.1–26.0‰) and the silicate-melt metasomatized olivines. The clinopyroxene phenocrysts and clinopyroxenes in the melt pockets have no distinct characteristics with respect to the Li abundances and δ⁷Li values, but show higher and lower CaO contents, respectively, than the clinopyroxenes from silicate and carbonatite metasomatized samples. These features indicate that Li concentration and isotopic signatures of the cores of the xenocrysts recorded carbonatite melt-peridotite reaction (carbonatite metasomatism) at mantle depth, and the variations in the rims probably resulted from xenocryst–host magma interaction during ascent. Our results reveal that the interaction with carbonatite and silicate melts gave rise to an increase in Li abundance in minerals of peridotite xenoliths at mantle depth or during transportation. In terms of δ⁷Li, the carbonatite and silicate melts produced remarkably contrasting δ⁷Li variations in olivine. Based on the systematic variations of Li abundances and Li isotopes in olivines, we suggest that the δ⁷Li value of olivine is a more important indicator than that of clinopyroxene in discriminating carbonatite and silicate melt interaction agents with peridotites.     

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 enrichment and isotopic variation in minerals from peridotite xenoliths from northwestern Ethiopian plateau

We report Lithium (Li) concentrations and isotopic compositions for co-existing olivine, orthopyroxene (opx), and clinopyroxene (cpx) mineral separates from depleted and metasomatised peridotite xenoliths hosted by basaltic lavas from northwestern Ethiopian plateau (Gundeweyn area). The peridotites contain five lherzolites and one harzburgite and are variably depleted and enriched in LREE relative to HREE. In both depleted and enriched lherzolites, Li is preferentially incorporated into olivine (2.4-3.3 ppm) compared to opx (1.4-2.1 ppm) and cpx (1.4-2.0 ppm) whereas the Li contents of olivines (5.4 ppm) from an enriched harzburgiteare higher than those of lherzolites. Olivines from the samples show higher Li abundances than normal mantle olivines (1.6-1.9 ppm) indicating the occurrence of Li enrichments through melt-preroditite interaction. The average δ⁷ Li values range from +2.2 to +6.0‰ in olivine, from -0.1 to +2.0‰ in opx and from -4.4 to -0.9‰ in cpx from the lherzolites. The Li isotopic composition (3.5‰) of olivines from harzburgite fall within the range of olivine from lherzolites but the opxs show low in δ⁷Li (-2.0‰). Overall Li isotopic compositions of olivines from the peridotites fall within the range of normal mantle olivine, δ⁷Li values of ~+4±2‰ within uncertainty, reflecting metasomatism (enrichment) of the peridotites by isotopically heavy Li-rich asthenospheric melt. Li isotope zonation is also observed in most peridotite minerals. Majority of olivine grains display isotopically heavy cores and light rims and the reverse case is observed for some olivine grains. Orthopyroxene and clinopyroxene grains show irregular distribution in δ⁷Li. These features of Li isotopic compositions within and between grains in the samples reflect the effect of diffusion-driven isotopic fractionation during meltperidotite interaction and cooling processes.     

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.     

O and H diffusion in uraninite: Implications for fluid-uraninite interactions, nuclear waste disposal, and nuclear forensics

Diffusion coefficients for oxygen and hydrogen were determined from a series of natural uraninite-H₂O experiments between 50 and 700 °C. Under hydrous conditions there are two diffusion mechanisms: (1) an initial extremely fast-path diffusion mechanism that overprinted the oxygen isotopic composition of the entire crystals regardless of temperature and (2) a slower volume-diffusive mechanism dominated by defect clusters that displace or eject nearest neighbor oxygen atoms to form two interstitial sites and two partial vacancies, and by vacancy migration. Using the volume diffusion coefficients in the temperature range of 400-600 °C, diffusion coefficients for oxygen can be represented by D = 1.90e⁻⁵ exp (−123,382 J/RT) cm²/s and for temperatures between 100 and 300 °C the diffusion coefficients can be represented by D = 1.95e⁻¹⁰ exp (−62484 J/RT) cm²/s, where the activation energies for uraninite are 123.4 and 62.5 kJ/mol, respectively. Hydrogen diffusion in uraninite appears to be controlled by similar mechanisms as oxygen. Using the volume diffusion coefficients for temperatures between 50 and 700 °C, diffusion coefficients for hydrogen can be represented by D = 9.28e⁻⁶ exp (−156,528 J/RT) cm²/s for temperatures between 450 and 700 °C and D = 1.39e⁻¹⁴ exp (−34518 J/RT) cm²/s for temperatures between 50 and 400 °C, where the activation energies for uraninite are 156.5 and 34.5 kJ/mol, respectively.Results from these new experiments have implications for isotopic exchange during natural UO₂-water interactions. The exceptionally low δ¹⁸O values of natural uraninites (i.e. −32‰ to −19.5‰) from unconformity-type uranium deposits in Saskatchewan, in conjunction with theoretical and experimental uraninite-water and UO₃-water fractionation factors, suggest that primary uranium mineralization is not in oxygen isotopic equilibrium with coeval clay and silicate minerals. The low δ¹⁸O values have been interpreted as resulting from the low temperature overprinting of primary uranium mineralization in the presence of relatively modern meteoric fluids having δ¹⁸O values of ca. −18‰, despite petrographic and U-Pb isotope data that indicate limited alteration. Our data show that the anomalously low oxygen isotopic composition of the uraninite from the Athabasca Basin can be due to meteoric water overprinting under reducing conditions, and meteoric water or groundwater can significantly affect the oxygen isotopic composition of spent nuclear fuel in a geologic repository, with minimal change to the chemical composition or texture. Moreover, the rather fast oxygen and hydrogen diffusion coefficients for uraninite, especially at low temperatures, suggest that oxygen and hydrogen diffusion may impart characteristic isotopic signals that can be used to track the route of fissile material.     

FTIR spectroscopy of OH in olivine: A new tool in kimberlite exploration

Our study of olivines from Canadian kimberlites shows that the application of FTIR spectroscopy significantly improves the reliability of olivine as a kimberlite indicator mineral (KIM). We have developed an algorithm that yields the water concentration and the normalized intensity of the OH IR absorption band at 3572 cm⁻¹ from unpolished olivine grains of unknown thickness. For 80% of kimberlitic olivines these two parameters are significantly higher than those for olivines from non-kimberlitic magmas and consequently, olivines with water concentrations >60 ppm and a strong absorption band at 3572 cm⁻¹ can be reliably classified as being kimberlitic.We have identified two major spectral features in the OH absorption bands of kimberlitic olivines that allow for a more detailed classification: (a) the presence of three types of high-requency OH absorption bands (Group 1A, 1B and 1C) and (b) the proportion of low-frequency OH absorption bands (Group 2) relative to high-frequency bands (Group 1). Comparison of our results with experimental studies suggests that differences within Group 1 OH absorption bands are due to different pressures of crystallization or hydrogenation. The three identified types of Group 1 OH absorption bands approximately correspond to high (P > 2 GPa, Group 1A), moderate (2-1 GPa, Group 1B), and low (<1 GPa, Group 1C) pressures of hydrogenation. Group 2 OH IR absorption bands in olivines with NiO > 3500 ppm are interpreted to reflect olivine-orthopyroxene equilibria and hence are indicative of xenocrystic olivine derived from lherzolitic or harzburgitic mantle sources. Interaction of xenocrystic olivine with hydrous kimberlitic melts with low silica activity likely will cause a gradual increase in Group 1 absorption bands. Therefore, FTIR spectra of olivine can be used to obtain qualitative estimates of the duration of interaction between mantle material and a kimberlitic melt.In addition to applications in kimberlite and diamond exploration, FTIR spectra of olivine phenocrysts, combined with mineral chemical data, may also provide insights into kimberlite evolution. Our data suggest that in some instances the ascent of kimberlitic magmas could have been interrupted at or near the Moho, followed by olivine crystallization and exsolution of aqueous fluids.     

Coupled lithium- and iron isotope fractionation during magmatic differentiation

In this study, we investigated Fe and Li isotope fractionation between mineral separates of olivine pheno- and xenocrysts (including one clinopyroxyene phenocryst) and their basaltic hosts. Samples were collected from the Canary Islands (Teneriffa, La Palma) and some German volcanic regions (Vogelsberg, Westerwald and Hegau). All investigated bulk samples fall in a tight range of Li and Fe isotope compositions (δ⁵⁶Fe_wr = 0.06–0.17‰ and δ⁷Li_ma = 2.5–5.2‰, assuming δ⁷Li of the olivine-free matrix is virtually identical to that of the bulk sample for mass balance reasons). In contrast, olivine phenocrysts display highly variable, but generally light Fe and mostly light Li isotope compositions compared to their respective olivine-free basaltic matrix, which was considered to represent the melt (with δ⁵⁶Fe_ol = ? 0.24 to 0.14‰ and δ⁷Li_ol = ? 10.5 to + 6.5‰, respectively). Single olivine crystals from one sample display even a larger range of δ⁵⁶Fe_ol between ? 0.7 and + 0.1‰. One single clinopyroxene phenocryst displays the lightest Li isotope composition (δ⁷Li_cpx = ? 17.7‰), but no Fe isotope fractionation relative to melt. The olivine phenocrysts show variable Mg# and Ni (correlated in most cases) that range between 0.89 and 0.74 and between 300 and 3000 μg/g, respectively. These olivines likely grew by fractional crystallization in an evolving magma. One sample from the Vogelsberg volcano contained olivine xenocrysts (Mg# > 0.89 and Ni > 3000 μg/g), in addition to olivine phenocrysts. This sample displays the highest Li- and the second highest Fe-isotope fractionation between olivine and melt (Δ⁷Li_ol-melt = ? 13; Δ⁵⁶Fe_ol-melt = ? 0.29).Our data, i.e. the variable olivine- at constant whole rock and matrix isotope compositions, strongly indicate disequilibrium, i.e. kinetic Fe and Li isotope fractionation between olivine and melt (for Li also between cpx and melt) during fractional crystallization. Δ⁷Li_ol-melt is correlated with the Li partitioning between olivine and melt (i.e. with Li_ol/Li_melt), indicating Li isotope fractionation due to preferential (faster) diffusion of ⁶Li into olivine during fractional crystallization. Olivine with low Δ⁷Li_ol-melt, also have low Δ⁵⁶Fe_ol-melt, indicating that Fe isotope fractionation is also driven by diffusion of isotopically light Fe into olivine, potentially, as Fe–Mg inter-diffusion. The lowest Δ⁵⁶Fe_ol-melt (? 0.40) was observed in a sample from Westerwald (Germany) with abundant magnetite, indicating relatively oxidizing conditions during magma differentiation. This may have enhanced equilibrium Fe isotope fractionation between olivine and melt or fine dispersed magnetite in the basalt matrix may have shifted its Fe isotope composition towards higher δ⁵⁶Fe. The decoupling of Li- and Fe isotope fractionation in cpx is likely due to faster diffusion of Li relative to Fe in cpx, implying that the large investigated cpx phenocryst resided in the magma for only a short period of time which was sufficient for Li- but not for Fe diffusion. The absence of any equilibrium Fe isotope fractionation between the investigated cpx phenocryst and its basaltic host may be related to the similar Fe^3 +/Fe^2 + of cpx and melt. In contrast to cpx, the generally light Fe isotope composition of all investigated olivine separates implies the existence of equilibrium- (in addition to diffusion-driven) isotope fractionation between olivine and melt, on the order of 0.1‰.     

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.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

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.     

Component mobility at 900 °C and 18 kbar from experimentally grown coronas in a natural gabbro

Several approximately 100-μm-wide reaction zones were grown under experimental conditions of 900 °C and 18 kbar along former olivine-plagioclase contacts in a natural gabbro. The reaction zone comprises two distinct domains: (i) an irregularly bounded zone with idiomorphic grains of zoisite and minor corundum and kyanite immersed in a melt developed at the plagioclase side and (ii) a well-defined reaction band comprising a succession of mineral layers forming a corona structure around olivine. Between the olivine and the plagioclase reactant phases we observe the following layer sequence: olivine|pyroxene|garnet|partially molten domain|plagioclase. Within the pyroxene layer two micro-structurally distinct layers comprising enstatite and clinopyroxene can be discerned. Chemical potential gradients persisted for the CaO, Al₂O₃, SiO₂, MgO and FeO components, which drove diffusion of Ca, Al and Si bearing species from the garnet-matrix interface to the pyroxene-olivine interface and diffusion of Mg- and Fe-bearing species in the opposite direction. The systematic mineralogical organization and chemical zoning across the corona suggest that the olivine corona was formed by a “diffusion-controlled” reaction. We estimate a set of diffusion coefficients and conclude that L_AlAl < L_CaCa < (L_SiSi, L_FeFe) < L_MgMg during reaction rim growth.     

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.     

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.     

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.     

Petrogenesis of olivine-phyric shergottite Yamato 980459, revisited

Primitive magmas provide critical information on mantle sources, but most Martian meteorites crystallized from fractionated melts. An olivine-phyric shergottite, Yamato 980459 (Y-980459), has been interpreted to represent a primary melt, because its olivine megacrysts have magnesian cores (Fo84-86) that appear to be in equilibrium with the Y-980459 whole-rock composition based on Fe-Mg partitioning. However, crystal size distribution (CSD) plots for Y-980459 olivines show a size gap, suggesting a cumulus origin for some megacrysts. Because melting experiments using the Y-980459 whole-rock composition have been used to infer the thermal structure and volatile contents of the Martian mantle, the interpretation that this rock is primitive should be scrutinized.We report major, minor and trace element compositions of Y-980459 olivines and compare them with results from melting experiments (both hydrous and anhydrous) and thermodynamic calculations. Cores of the olivine megacrysts have major and minor element contents identical to those of the most magnesian olivines from the experiments, but they differ slightly from those of thermodynamic calculations. This is probably because the Y-980459 whole-rock composition lies near the limit of the range of liquids used to calibrate these models. The megacryst cores (Fo80-85) exhibit minor and trace element (Mn-Ni-Co-Cr-V) characteristics distinct from other olivines (megacryst rims and groundmass olivines, Fo < 80), implying that the megacryst cores crystallized under more reduced conditions (∼IW + 1).Y-980459 contains pyroxenes with orthopyroxene cores mantled by pigeonite and augite. We also found some reversely zoned pyroxenes that have augite cores (low-Mg#) mantled by orthopyroxenes (high-Mg#), although they are uncommon. These reversely zoned pyroxenes are interpreted to have grown initially as atoll-like crystals with later crystallization filling in the hollow centers, implying disequilibrium crystallization at a moderate cooling rate (3-7 °C/h). The calculated REE pattern of a melt in equilibrium with normally zoned pyroxene is parallel to those of glass and the Y-980459 whole-rock as well as other depleted olivine-phyric shergottites, suggesting that Y-980459 was derived from a depleted mantle reservoir.Considering the CSD patterns of Y-980459 olivines, we propose that the olivine megacrysts are cumulus crystals which probably formed in a feeder conduit by continuous melt replenishment, and the parent melt composition was indistinguishable from the Y-980459 whole-rock with 0-2 wt% of H₂O and 0-5 wt% of CO₂. The final magma pulse entrained these cumulus olivines and then crystallized groundmass olivines and pyroxenes. Although Y-980459 contains small amounts of cumulus olivine (<∼6 vol%), we conclude that the Y-980459 whole-rock composition closely approximates a Martian primary melt composition.