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.     

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.     

Lithium isotopes in Guatemalan and Franciscan HP-LT rocks: Insights into the role of sediment-derived fluids during subduction

High-pressure, low-temperature (HP-LT) rocks from a Cretaceous age subduction complex occur as tectonic blocks in serpentinite mélange along the Motagua Fault (MF) in central Guatemala. Eclogite and jadeitite among these are characterized by trace element patterns with enrichments in fluid mobile elements, similar to arc lavas. Eclogite is recrystallized from MORB-like altered oceanic crust, presumably at the boundary between the down-going plate and overlying mantle wedge. Eclogite geochemistry, mineralogy and petrography suggest a two step petrogenesis of (1) dehydration during prograde metamorphism at low temperatures (<500 °C) followed by (2) partial rehydration/fertilization at even lower T during exhumation. In contrast, Guatemalan jadeitites are crystallized directly from low-T aqueous fluid as veins in serpentinizing mantle during both subduction and exhumation. The overall chemistry and mineralogy of Guatemalan eclogites are similar to those from the Franciscan Complex, California, implying similar P-T-x paths.Li concentrations (?90 ppm) in mineral separates and whole rocks (WR) from Guatemalan and Franciscan HP-LT rocks are significantly higher than MORB (4-6 ppm), but similar to HP-LT rocks globally. Li isotopic compositions range from −5‰ to +5‰ for Guatemalan HP-LT rocks, and −4‰ to +1‰ for Franciscan eclogites, overlapping previous findings for other HP-LT suites. The combination of Li concentrations greater than MORB, and Li isotopic values lighter than MORB are inconsistent with a simple dehydration model. We prefer a model in which Li systematics in Guatemalan and Franciscan eclogites reflect reequilibration with subduction fluids during exhumation. Roughly 5-10% of the Li in these fluids is derived from sediments.Model results predict that the dehydrated bulk ocean crust is isotopically lighter (δ⁷Li ? +1 ± 3‰) than the depleted mantle (∼+3.5 ± 0.5‰), while the mantle wedge beneath the arc is the isotopic complement of the bulk crust. A subduction fluid with an AOC-GLOSS composition over the full range of model temperatures (50-600 °C) gives an average fluid δ⁷Li (∼+7 ± 5‰ 1σ) that is isotopically heavier than the depleted mantle. If the lowest temperature steps are excluded (50-260 °C) as too cold to participate in circulation of the mantle wedge, then the average subduction fluid (δ⁷Li = +4 ± 2.3‰ 1σ, is indistinguishable from depleted mantle. Because of the relatively compatible nature of Li in metamorphic minerals, the most altered part of the crust (uppermost extrusives), may retain a Li isotopic signature (∼+5 ± 3‰) heavier than the bulk crust. The range of Li isotopic values for OIB, IAB and MORB overlap, making it is difficult to resolve which of these components may contribute to the recycled component in the mantle using δ⁷Li alone.     

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.     

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.     

Geothermal waters from the Taupo Volcanic Zone,New Zealand: Li,B and Sr isotopes characterization

Chemical and isotopic data for 23 geothermal water samples collected in New Zealand within the Taupo Volcanic Zone (TVZ) are reported. Major and trace elements including Li, B and Sr and their isotopic compositions (δ⁷Li, δ¹¹B, ⁸⁷Sr/⁸⁶Sr) were determined in high temperature geothermal waters collected from deep boreholes in different geothermal fields (Ohaaki, Wairakei, Mokai, Kawerau and Rotokawa geothermal systems). Lithium concentrations are high (from 4.5 to 19.9 mg/L) and Li isotopic compositions (δ⁷Li) are homogeneous, ranging between −0.5‰ and +1.4‰. In particular, it is noteworthy that, except for the samples from the Kawerau geothermal field having slightly higher δ⁷Li values (+1.4%), the other geothermal waters have a near constant δ⁷Li signature around a mean value of 0‰ ± 0.6 (2σ, n = 21). Boron concentrations are also high and relatively homogeneous for the geothermal samples, falling between 17.5 and 82.1 mg/L. Boron isotopic compositions (δ¹¹B) are all negative, and display a range between −6.7‰ and −1.9‰. These B isotope compositions are in agreement with those of the Ngawha geothermal field in New Zealand. Lithium and B isotope signatures are in a good agreement with a fluid signature mainly derived from water/rock interaction involving magmatic rocks with no evidence of seawater input. On the other hand, Sr concentrations are lower and more heterogeneous and fall between 2 and 165 μg/L. The ⁸⁷Sr/⁸⁶Sr ratios range from 0.70549 to 0.70961. These Sr isotope compositions overlap those of the Rotorua geothermal field in New Zealand, confirming that some geothermal waters (with more radiogenic Sr) have interacted with bedrocks from the metasedimentary basement. Each of these isotope systems on their own reveals important information about particular aspects of either water source or water/rock interaction processes, but, considered together, provide a more integrated understanding of the geothermal systems from the TVZ in New Zealand.     

Carbon and hydrogen isotope variation of plant biomarkers in a plant-soil system

Stable carbon (δ¹³C) and hydrogen (δD) isotopic compositions of n-alkanes, anteiso-alkanes, n-alkanoic acids, n-alkanols, phytol and sterols in raw leaves of Acer argutum and Acer carpinifolium, their fallen leaves, mold and soils from a natural Acer forest were measured in order to: (1) understand isotopic variation of the plant biomarkers in a plant-soil system and (2) evaluate which biomarker is the most effective recorder of soil vegetation. Long-chain (> C₂₄) n-alkanes, n-alkanoic acids and n-alkanols are gradually enriched in ¹³C up to 12.9‰ (average of 4.3‰) and depleted in D up to 94‰ (average of 55‰) from raw leaves to soils. However, anteiso-alkanes, phytol and sterols show little variation in both δ¹³C (< ± 1‰) and δD (< ± 2‰) from raw leaves to soils. These isotope signatures in a plant-soil system indicate that isoprenoid plant biomarkers such as sterols in soils faithfully preserve the isotopic compositions of dominant higher plants growing on the soils without a diagenetic effect upon the isotopic compositions. In contrast, long-chain n-alkyl molecules in soils undergo specific isotopic modification during biodegradation associated with early diagenesis and/or a significant contribution from heterotrophic reworking.     

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.     

Use of B isotopes as a tracer of anthropogenic emissions in the atmosphere of Paris,France

In order to test the potential of B isotopes as a tracer of contamination of the atmosphere, the B isotopic composition of rainwater samples monitored over a year in the centre of Paris, France were determined. Boron concentrations range from 19 nmol/L to 500 nmol/L and δ¹¹B range from 0‰ to +38‰. Mean annual values are 148 nmol/L and +25‰, respectively. The results suggest that variability in B isotopic compositions is mainly caused by mixing of two main sources, although isotopic fractionation during the evaporation–condensation processes may also be important. One source is a marine component, which exhibits a heavy B isotopic composition. The decrease of δ¹¹B in rainwater with increasing NO₃/B and SO₄/B molar ratios suggests that a second source may be anthropogenic emissions. To constrain this end-member, B was determined in urban particulates, which were enriched in the light isotope and the lowest values were consistent with a B-rich fossil fuel composition. These results confirm the great sensitivity of B to anthropogenic sources and the ability of B isotopic ratios to decipher the origin of B in the atmosphere.     

Stable isotope chemistry of fossil bone as a new paleoclimate indicator

During fossilization, bone is thought to recrystallize and alter chemically on timescales of kyr to a few tens of kyr, i.e., similar to the timescale for formation of soils. Therefore, C- and O-isotope compositions of bone apatite should correlate with trends in soil water composition and aridity, and serve as paleoclimate indicators. This hypothesis was tested by analyzing C- and O-isotope compositions of the CO₃ component of fossil bone apatite from mid-Oligocene through late Pleistocene units in Oregon and western Idaho, including the John Day (19.4-30.0 Ma), Mascall (15.2-15.8 Ma), and Rattlesnake (7.2-7.8 Ma) Formations, whose paleosol sequences have been studied in detail, and the Juntura (10-11 Ma), Hagerman (3.2 Ma), and Fossil Lake (<23-650 ka) fossil localities. Tooth enamel δ¹⁸O values provide a baseline of meteoric water compositions. Stable isotope compositions of bone CO₃ do change in response to broad climatic trends, but show poor correlation with compositions of corresponding paleosol CO₃ at specific horizons. Instead, compositional deviations between bone and paleosol CO₃ correlate with compositional deviations with the next higher paleosol; this suggests that the timescale for fossilization exceeds one paleosol cycle. Based on stratigraphic evidence and simple alteration models, fossilization timescales are estimated at 20-50 kyr, indicating that bone CO₃ will prove most useful for sequences spanning >100 kyr. C-isotopes show negative and strong positive deviations during wet and dry climates respectively, and short-term trends correspond well with changes in aridity within the Mascall and Rattlesnake Formations, as inferred from paleosols. A proposed correction to δ¹⁸O values based on δ¹³C anomalies implies a small, ∼1.5‰ increase in meteoric water δ¹⁸O during the late Oligocene global warming event, consistent with a minimum temperature increase of ∼4 °C. A strong inferred decrease in δ¹⁸O of 4-5‰ after 7 Ma closely parallels compositional changes in tooth enamel, and reflects a doubling in the height of the Cascade Range.     

Lithium isotope geochemistry of marine pore waters - Insights from cold seep fluids

Lithium concentration and isotope data (δ⁷Li) are reported for pore fluids from 18 cold seep locations together with reference fluids from shallow marine environments, a sediment-hosted hydrothermal system and two Mediterranean brine basins. The new reference data and literature data of hydrothermal fluids and pore fluids from the Ocean Drilling Program follow an empirical relationship between Li concentration and δ⁷Li (δ⁷Li = −6.0(±0.3) · ln[Li] + 51(±1.2)) reflecting Li release from sediment or rocks and/or uptake of Li during mineral authigenesis. Cold seep fluids display δ⁷Li values between +7.5‰ and +45.7‰, mostly in agreement with this general relationship. Ubiquitous diagenetic signals of clay dehydration in all cold seep fluids indicate that authigenic smectite-illite is the major sink for light pore water Li in deeply buried continental margin sediments. Deviations from the general relationship are attributed to the varying provenance and composition of sediments or to transport-related fractionation trends. Pore fluids on passive margins receive disproportionally high amounts of Li from intensely weathered and transported terrigenous matter. By contrast, on convergent margins and in other settings with strong volcanogenic input, Li concentrations in pore water are lower because of intense Li uptake by alteration minerals and, most notably, adsorption of Li onto smectite. The latter process is not accompanied by isotope fractionation, as revealed from a separate study on shallow sediments. A numerical transport-reaction model was applied to simulate Li isotope fractionation during upwelling of pore fluids. It is demonstrated that slow pore water advection (order of mm a⁻¹) suffices to convey much of the deep-seated diagenetic Li signal into shallow sediments. If carefully applied, Li isotope systematics may, thus, provide a valuable record of fluid/mineral interaction that has been inherited several hundreds or thousands of meters below the actual seafloor fluid escape structure.     

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.     

Note on the isotopic geochemistry of fossil-lacustrine tufas in carbonate plateau—A study from Dungul region (SW Egypt)

Lithological, chemical, and stable isotope data are used to characterize lacustrine tufas dating back to pre-late Miocene and later unknown times, capping different surfaces of a Tertiary carbonate (Sinn el-Kedab) plateau in Dungul region in the currently hyperarid south-western Egypt. These deposits are composed mostly of calcium carbonate, some magnesium carbonate and clastic particles plus minor amounts of organic matter. They have a wide range of (Mg/Ca)_molar ratios, from 0.03 to 0.3. The bulk-tufa carbonate has characteristic isotope compositions: (δ¹³C_mean = −2.49 ± 0.99‰; δ¹⁸O_mean = −9.43 ± 1.40‰). The δ¹³C values are consistent with a small input from C4 vegetation or thinner soils in the recharge area of the tufa-depositing systems. The δ¹⁸O values are typical of fresh water carbonates. Covariation between δ¹³C and δ¹⁸O values probably is a reflection of climatic conditions such as aridity. The tufas studied are isotopically similar to the underlying diagenetic marine chalks, marls and limestones (δ¹³C_mean = −2.06 ± 0.84‰; δ¹⁸O_mean = −10.06 ± 1.39‰). The similarity has been attributed to common meteoric water signatures. This raises large uncertainties in using tufas (Mg/Ca)_molar, δ¹³C and δ¹⁸O records as proxies of paleoclimatic change and suggests that intrinsic compositional differences in material sources within the plateau may mask climatic changes in the records.     

The triple isotopic composition of oxygen in leaf water

The isotopic composition of atmospheric O₂ depends on the rates of oxygen cycling in photosynthesis, respiration, photochemical reactions in the stratosphere and on δ¹⁷O and δ¹⁸O of ocean and leaf water. While most of the factors affecting δ¹⁷O and δ¹⁸O of air O₂ have been studied extensively in recent years, δ¹⁷O of leaf water—the substrate for all terrestrial photosynthesis—remained unknown. In order to understand the isotopic composition of atmospheric O₂ at present and in fossil air in ice cores, we studied leaf water in field experiments in Israel and in a European survey. We measured the difference in δ¹⁷O and δ¹⁸O between stem and leaf water, which is the result of isotope enrichment during transpiration. We calculated the slopes of the lines linking the isotopic compositions of stem and leaf water. The obtained slopes in ln(δ¹⁷O + 1) vs. ln(δ¹⁸O + 1) plots are characterized by very high precision (∼0.001) despite of relatively large differences between duplicates in both δ¹⁷O and δ¹⁸O (0.02-0.05‰). This is so because the errors in δ¹⁸O and δ¹⁷O are mass-dependent. The slope of the leaf transpiration process varied between 0.5111 ± 0.0013 and 0.5204 ± 0.0005, which is considerably smaller than the slope linking liquid water and vapor at equilibrium (0.529). We further found that the slope of the transpiration process decreases with atmospheric relative humidity (h) as 0.522-0.008 × h, for h in the range 0.3-1. This slope is neither influenced by the plant species, nor by the environmental conditions where plants grow nor does it show strong variations along long leaves.     

Nature of volatile depletion and genetic relationships in enstatite chondrites and aubrites inferred from Zn isotopes

Enstatite meteorites include the undifferentiated enstatite chondrites and the differentiated enstatite achondrites (aubrites). They are the most reduced group of all meteorites. The oxygen isotope compositions of both enstatite chondrites and aubrites plot along the terrestrial mass fractionation line, which suggests some genetic links between these meteorites and the Earth as well.For this study, we measured the Zn isotopic composition of 25 samples from the following groups: aubrites (main group and Shallowater), EL chondrites, EH chondrites and Happy Canyon (impact-melt breccia). We also analyzed the Zn isotopic composition and elemental abundance in separated phases (metal, silicates, and sulfides) of the EH4, EL3, and EL6 chondrites. The different groups of meteorites are isotopically distinct and give the following values (‰): aubrite main group (−7.08 < δ⁶⁶Zn < −0.37); EH3 chondrites (0.15 < δ⁶⁶Zn < 0.31); EH4 chondrites (0.15 < δ⁶⁶Zn < 0.27); EH5 chondrites (δ⁶⁶Zn = 0.27 ± 0.09; n = 1); EL3 chondrites (0.01 < δ⁶⁶Zn < 0.63); the Shallowater aubrite (1.48 < δ⁶⁶Zn < 2.36); EL6 chondrites (2.26 < δ⁶⁶Zn < 7.35); and the impact-melt enstatite chondrite Happy Canyon (δ⁶⁶Zn = 0.37).The aubrite Peña Blanca Spring (δ⁶⁶Zn = −7.04‰) and the EL6 North West Forrest (δ⁶⁶Zn = 7.35‰) are the isotopically lightest and heaviest samples, respectively, known so far in the Solar System. In comparison, the range of Zn isotopic composition of chondrites and terrestrial samples (−1.5 < δ⁶⁶Zn < 1‰) is much smaller ( [Luck et al., 2005] and [Herzog et al., 2009]).EH and EL3 chondrites have the same Zn isotopic composition as the Earth, which is another example of the isotopic similarity between Earth and enstatite chondrites. The Zn isotopic composition and abundance strongly support that the origin of the volatile element depletion between EL3 and EL6 chondrites is due to volatilization, probably during thermal metamorphism. Aubrites show strong elemental depletion in Zn compared to both EH and EL chondrites and they are enriched in light isotopes (δ⁶⁶Zn down to −7.04‰). This is the opposite of what would be expected if Zn elemental depletion was due to evaporation, assuming the aubrites started with an enstatite chondrite-like Zn isotopic composition. Evaporation is therefore not responsible for volatile loss from aubrites. On Earth, Zn isotopes fractionate very little during igneous processes, while differentiated meteorites show only minimal Zn isotopic variability. It is therefore very unlikely that igneous processes can account for the large isotopic fractionation of Zn in aubrites. Condensation of an isotopically light vapor best explains Zn depletion and isotopically light Zn in these puzzling rocks. Mass balance suggests that this isotopically light vapor carries Zn lost by the EL6 parent body during thermal metamorphism and that aubrites evolved from an EL6-like parent body. Finally, Zn isotopes suggest that Shallowater and aubrites originate from distinct parent bodies.     

Heterogeneous magnesium isotopic composition of the upper continental crust

High-precision Mg isotopic data are reported for ∼100 well-characterized samples (granites, loess, shales and upper crustal composites) that were previously used to estimate the upper continental crust composition. Magnesium isotopic compositions display limited variation in eight I-type granites from southeastern Australia (δ²⁶Mg = −0.25 to −0.15) and in 15 granitoid composites from eastern China (δ²⁶Mg = −0.35 to −0.16) and do not correlate with SiO₂ contents, indicating the absence of significant Mg isotope fractionation during differentiation of granitic magma. Similarly, the two S-type granites, which represent the two end-members of the S-type granite spectrum from southeastern Australia, have Mg isotopic composition (δ²⁶Mg = −0.23 and −0.14) within the range of their potential source rocks (δ²⁶Mg = −0.20 and +0.15) and I-type granites, suggesting that Mg isotope fractionation during crustal anatexis is also insignificant. By contrast, δ²⁶Mg varies significantly in 19 A-type granites from northeastern China (−0.28 to +0.34) and may reflect source heterogeneity.Compared to I-type and S-type granites, sedimentary rocks have highly heterogeneous and, in most cases, heavier Mg isotopic compositions, with δ²⁶Mg ranging from −0.32 to +0.05 in nine loess from New Zealand and the USA, from −0.27 to +0.49 in 20 post-Archean Australian shales (PAAS), and from −0.52 to +0.92 in 20 sedimentary composites from eastern China. With increasing chemical weathering, as measured by the chemical index of alternation (CIA), δ²⁶Mg values show a larger dispersion in shales than loess. Furthermore, δ²⁶Mg correlates negatively with δ⁷Li in loess. These characteristics suggest that chemical weathering significantly fractionates Mg isotopes and plays an important role in producing the highly variable Mg isotopic composition of sedimentary rocks.Based on the estimated proportions of major rock units within the upper continental crust and their average MgO contents, a weighted average δ²⁶Mg value of −0.22 is derived for the average upper continental crust. Our studies indicate that Mg isotopic composition of the upper crust is, on average, mantle-like but highly heterogeneous, with δ²⁶Mg ranging from −0.52 to +0.92. Such large isotopic variation mainly results from chemical weathering, during which light Mg isotopes are lost to the hydrosphere, leaving weathered products (e.g., sedimentary rocks) with heavy Mg isotopes.     

Transfer and early diagenesis of biogenic silica oxygen isotope signals during settling and sedimentation of diatoms in a temperate freshwater lake (Lake Holzmaar, Germany)

We have investigated the transfer of oxygen isotope signals of diatomaceous silica (δ¹⁸O_diatom) from the epilimnion (0-7 m) through the hypolimnion to the lake bottom (∼20 m) in freshwater Lake Holzmaar, Germany. Sediment-traps were deployed in 2001 at depths of 7 and 16 m to harvest fresh diatoms every 28 days. The 7 m trap collected diatoms from the epilimnion being the main zone of primary production, while the 16 m trap collected material already settled through the hypolimnion. Also a bottom sediment sample was taken containing diatom frustules from approximately the last 25 years. The δ¹⁸O_diatom values of the 7 m trap varied from 29.4‰ in spring/autumn to 26.2‰ in summer according to the temperature dependence of oxygen isotope fractionation and represent the initial isotope signal in this study. Remarkably, despite the short settling distance δ¹⁸O_diatom values of the 7 and the 16 m trap were identical only during spring and autumn seasons while from April to September δ¹⁸O_diatom values of the 16 m trap were roughly ∼1.5‰ enriched in ¹⁸O compared to those of the 7 m trap. Isotopic exchange with the isotopically lighter water of the hypolimnion would shift the δ¹⁸O_diatom value to lower values during settling from 7 to 16 m excluding this process as a cause for the deviation. Dissolution of opal during settling with intact organic coatings of the diatom cells and near neutral pH of the water should only cause a minor enrichment of the 16 m values. Nevertheless, opal from the bottom sediment was found to be 2.5‰ enriched in ¹⁸O compared to the weighted average of the opal from the 7 m trap. Thus, resuspension of bottom material must have contributed to the intermediate δ¹⁸O_diatom signal of the 16 m trap during summer. Dissolution experiments allowed further investigation of the cause for the remarkably enriched δ¹⁸O_diatom value of the bottom sediment. Experiments with different fresh diatomaceous materials show an increase of opaline ¹⁸O at high pH values which is remarkably reduced when organic coatings of the cells still exist or at near neutral pH. In contrast, high pH conditions do not affect the δ¹⁸O_diatom values of sub-fossil and even fossil opal. IR analyses show that the ¹⁸O enrichment of the sedimentary silica is associated with a decrease in Si-OH groups and the formation of Si-O-Si linkages. This indicates a silica dehydroxylation process as cause for the isotopic enrichment of the bottom sediment. Silica dissolution and dehydroxylation clearly induce a maturation process of the diatom oxygen isotope signal presumably following an exponential behaviour with a rapid initial phase of signal alteration. The dynamics of this process is of particular importance for the quantitative interpretation of sedimentary δ¹⁸O_diatom values in terms of palaeothermometry.     

Insights into Li and Li isotope cycling and sub-arc metasomatism from veined mantle xenoliths,Kamchatka

Harzburgitic xenoliths cut by pyroxenite veins from Avachinsky volcano, Kamchatka, are derived from the sub-arc mantle and record element transfer from the slab to the arc. Olivine and orthopyroxene in the harzburgites have Li isotopic compositions (δ⁷Li = +2.8 to +5.6) comparable to estimates of the upper mantle (δ⁷Li ~ +4 ± 2). The pyroxenite veins, which represent modal metasomatism and may therefore provide information about the metasomatic agent, have mantle-normalized trace element characteristics that suggest overprinting of their mantle source by an aqueous, slab-derived fluid. These include relative enrichments of Pb over Ce, U over Th and Sr over Nd. Li is enriched relative to the HREE, and ortho- and clinopyroxene from the veins are in Li elemental and isotopic equilibrium with each other and the surrounding harzburgite. Vein samples (δ⁷Li = +3.0 to +5.0) do not record a significant slab-derived δ⁷Li signature. These observations can be reconciled if slab Li diffusively re-equilibrates in the mantle wedge. Modeling demonstrates that Li equilibration of small (1–2 cm width) veins or melt conduits is achieved at mantle wedge temperatures within 10¹–10⁵ years. We conclude that strongly fractionated Li isotopic signatures cannot be sustained for long periods in the sub-arc mantle, at least at shallow (<70 km) depths. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Climatic dependence of stable carbon and oxygen isotope signals recorded in speleothems: From soil water to speleothem calcite

Understanding the relationship between stable isotope signals recorded in speleothems (δ¹³C and δ¹⁸O) and the isotopic composition of the carbonate species in the soil water is of great importance for their interpretation in terms of past climate variability. Here the evolution of the carbon isotope composition of soil water on its way down to the cave during dissolution of limestone is studied for both closed and open-closed conditions with respect to CO₂.The water entering the cave flows as a thin film towards the drip site. CO₂ degasses from this film within approx. 10 s by molecular diffusion. Subsequently, chemical and isotopic equilibrium is established on a time scale of several 10-100 s. The δ¹³C value of the drip water is mainly determined by the isotopic composition of soil CO₂. The evolution of the δ¹⁸O value of the carbonate species is determined by the long exchange time T_ex, between oxygen in carbonate and water of several 10,000 s. Even if the oxygen of the CO₂ in soil water is in isotopic equilibrium with that of the water, dissolution of limestone delivers oxygen with a different isotopic composition changing the δ¹⁸O value of the carbonate species. Consequently, the δ¹⁸O value of the rainwater will only be reflected in the drip water if it has stayed in the rock for a sufficiently long time.After the water has entered the cave, the carbon and oxygen isotope composition of the drip water may be altered by CO₂-exchange with the cave air. Exchange times, , of about 3000 s are derived. Thus, only drip water, which drips in less than 3000 s onto the stalagmite surface, is suitable to imprint climatic signals into speleothem calcite deposited from it.Precipitation of calcite proceeds with time constants, τ_p, of several 100 s. Different rate constants and equilibrium concentrations for the heavy and light isotopes, respectively, result in isotope fractionation during calcite precipitation. Since T_ex ? τ_p, exchange with the oxygen in the water can be neglected, and the isotopic evolution of carbon and oxygen proceed analogously. For drip intervals T_d < 0.1τ_p the isotopic compositions of both carbon and oxygen in the solution evolve linearly in time. The calcite precipitated at the apex of the stalagmite reflects the isotopic signal of the drip water.For long drip intervals, when calcite is deposited from a stagnant water film, long drip intervals may have a significant effect on the isotopic composition of the DIC. In this case, the isotopic composition of the calcite deposited at the apex must be determined by averaging over the drip interval. Such processes must be considered when speleothems are used as proxies of past climate variability.     

Semi-automatic determination of the carbon and oxygen stable isotope compositions of calcite and dolomite in natural mixtures

A semi-automatic, on-line method was developed to determine the δ¹³C and δ¹⁸O values of coexisting calcite and dolomite. An isotopic mass balance is used to calculate the compositions of dolomite after having measured that of calcite and of the “bulk” sample. The limit of validity of this method is established by performing isotopic measurements of artificial mixtures made of precisely weighted and isotopically-characterised dolomite and calcite. The accuracy and repeatability of the calculation of dolomite δ¹³C and δ¹⁸O are statistically determined with a Monte-Carlo procedure of error propagation. Stable isotope ratios are determined by using an automated MultiPrep™ system on-line with an isotope-ratio mass-spectrometer (IRMS). The reaction time and the temperature of reaction were optimised by comparing the results with the isotopic composition of known mixtures. The best results were obtained by phosphoric acid digestion after 20 min at 40 °C for calcite and 45 min at 90 °C for dolomite. This procedure allows an accurate determination of the isotopic ratios from small samples (300 μg). Application of this protocol to natural mixtures of calcite and dolomite requires the accurate determination of the relative abundance of calcite and dolomite by combining Mélières manocalcimetry (MMC) and X-ray diffractometry (XRD).