The effect of early diagenesis on the Fe isotope compositions of porewaters and authigenic minerals in continental margin sediments

Iron isotope compositions in marine pore fluids and sedimentary solid phases were measured at two sites along the California continental margin, where isotope compositions range from δ⁵⁶Fe = −3.0‰ to +0.4‰. At one site near Monterey Canyon off central California, organic matter oxidation likely proceeds through a number of diagenetic pathways that include significant dissimilatory iron reduction (DIR) and bacterial sulfate reduction, whereas at our other site in the Santa Barbara basin DIR appears to be comparatively small, and production of sulfides (FeS and pyrite) was extensive. The largest range in Fe isotope compositions is observed for Fe(II)_aq in porewaters, which generally have the lowest δ⁵⁶Fe values (minimum: −3.0‰) near the sediment surface, and increase with burial depth. δ⁵⁶Fe values for FeS inferred from HCl extractions vary between ∼−0.4‰ and +0.4‰, but pyrite is similar at both stations, where an average δ⁵⁶Fe value of −0.8 ± 0.2‰ was measured. We interpret variations in dissolved Fe isotope compositions to be best explained by open-system behavior that involves extensive recycling of Feflux. This study is the first to examine Fe isotope variations in modern marine sediments, and the results show that Fe isotopes in the various reactive Fe pools undergo isotopic fractionation during early diagenesis. Importantly, processes dominated by sulfide formation produce high-δ⁵⁶Fe values for porewaters, whereas the opposite occurs when Fe(III)-oxides are present and DIR is a major pathway of organic carbon respiration. Because shelf pore fluids may carry a negative δ⁵⁶Fe signature it is possible that the Fe isotope composition of ocean water reflects a significant contribution of shelf-derived iron to the open ocean. Such a signature would be an important means for tracing iron sources to the ocean and water mass circulation.     

Isotopic evidence (Sr/Sr, δLi) for alteration of the oceanic crust at deep-rooted mud volcanoes in the Gulf of Cadiz, NE Atlantic Ocean

The chemical and isotopic composition of pore fluids is presented for five deep-rooted mud volcanoes aligned on a transect across the Gulf of Cadiz continental margin at water depths between 350 and 3860 m. Generally decreasing interstitial Li concentrations and ⁸⁷Sr/⁸⁶Sr ratios with increasing distance from shore are attributed to systematically changing fluid sources across the continental margin. Although highest Li concentrations at the near-shore mud volcanoes coincide with high salinities derived from dissolution of halite and late-stage evaporites, clayey, terrigenous sediments are identified as the ultimate Li source to all pore fluids investigated. Light δ⁷Li values, partly close to those of hydrothermal vent fluids (δ⁷Li: +11.9‰), indicate that Li has been mobilized during high-temperature fluid/sediment or fluid/rock interactions in the deep sub-surface. Intense leaching of terrigenous clay has led to radiogenic ⁸⁷Sr/⁸⁶Sr ratios (∼0.7106) in pore fluids of the near-shore mud volcanoes. In contrast, non-radiogenic ⁸⁷Sr/⁸⁶Sr ratios (∼0.7075) at the distal locations are attributed to admixing of a basement-derived fluid component, carrying an isotopic signature from interaction with the basaltic crust. This inference is substantiated by temperature constraints from Li isotope equilibrium calculations suggesting exchange processes at particularly high temperatures (>200 °C) for the least radiogenic pore fluids of the most distal location.Advective pore fluids in the off-shore reaches of the Gulf of Cadiz are influenced by successive exchange processes with both oceanic crust and terrigenous, fine-grained sediments, resulting in a chemical and isotopic signature similar to that of fluids in near-shore ridge flank hydrothermal systems. This suggests that deep-rooted mud volcanoes in the Gulf of Cadiz represent a fluid pathway intermediate between mid-ocean ridge hydrothermal vent and shallow, marginal cold seep. Due to the thicker sediment coverage and slower fluid advection rates, the overall geochemical signature is shifted towards the sediment-diagenetic signal compared to ridge flank hydrothermal environments.     

Li isotope fractionation in peridotites and mafic melts

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

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.     

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.     

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.     

Silicon isotope and trace element constraints on the origin of ∼3.5 Ga cherts: Implications for Early Archaean marine environments

Silicon (Si) isotope variability in Precambrian chert deposits is significant, but proposed explanations for the observed heterogeneity are incomplete in terms of silica provenance and fractionation mechanisms involved. To address these issues we investigated Si isotope systematics, in conjunction with geochemical and mineralogical data, in three well-characterised and approximately contemporaneous, ∼3.5 Ga chert units from the Pilbara greenstone terrane (Western Australia).We show that Si isotope variation in these cherts is large (−2.4‰ to +1.3‰) and was induced by near-surface processes that were controlled by ambient conditions. Cherts that formed by chemical precipitation of silica show the largest spread in δ³⁰Si (−2.4‰ to +0.6‰) and are characterised by positive Eu, La and Y anomalies and overall depletions in lithophile trace elements. Silicon isotope systematics in these orthochemical deposits are explained by (1) mixing between hydrothermal fluids and seawater, and/or (2) fractionation of hydrothermal fluids by subsurface losses of silica due to conductive cooling. Rayleigh-type fractionation of hydrothermal fluids was largely controlled by temperature differences between these fluids and seawater. Lamina-scale Si isotope heterogeneity within individual chemical chert samples up to 2.2‰ is considered to reflect the dynamic nature of hydrothermal activity. Silicified volcanogenic sediments lack diagnostic REE+Y anomalies, are enriched in lithophile elements, and exhibit a much more restricted range of positive δ³⁰Si (+0.1‰ to +1.1‰), which points to seawater as the dominant source of silica.The proposed model for Si isotope variability in the Early Archaean implies that chemical cherts with the most negative δ³⁰Si formed from pristine hydrothermal fluids, whereas silicified or chemical sediments with positive δ³⁰Si are closest to pure seawater deposits. Taking the most positive value found in this study (+1.3‰), and assuming that the Si isotope composition of seawater is governed by input of fractionated hydrothermal fluids, we infer that the temperature of ∼3.5 Ga seawater was below ∼55 °C.     

Ca isotopes in carbonate sediment and pore fluid from ODP Site 807A: The Ca(aq)-calcite equilibrium fractionation factor and calcite recrystallization rates in Pleistocene sediments

The calcium isotopic compositions (δ⁴⁴Ca) of 30 high-purity nannofossil ooze and chalk and 7 pore fluid samples from ODP Site 807A (Ontong Java Plateau) are used in conjunction with numerical models to determine the equilibrium calcium isotope fractionation factor (α_s−f) between calcite and dissolved Ca²⁺ and the rates of post-depositional recrystallization in deep sea carbonate ooze. The value of α_s−f at equilibrium in the marine sedimentary section is 1.0000 ± 0.0001, which is significantly different from the value (0.9987 ± 0.0002) found in laboratory experiments of calcite precipitation and in the formation of biogenic calcite in the surface ocean. We hypothesize that this fractionation factor is relevant to calcite precipitation in any system at equilibrium and that this equilibrium fractionation factor has implications for the mechanisms responsible for Ca isotope fractionation during calcite precipitation. We describe a steady state model that offers a unified framework for explaining Ca isotope fractionation across the observed precipitation rate range of ∼14 orders of magnitude. The model attributes Ca isotope fractionation to the relative balance between the attachment and detachment fluxes at the calcite crystal surface. This model represents our hypothesis for the mechanism responsible for isotope fractionation during calcite precipitation. The Ca isotope data provide evidence that the bulk rate of calcite recrystallization in freshly-deposited carbonate ooze is 30-40%/Myr, and decreases with age to about 2%/Myr in 2-3 million year old sediment. The recrystallization rates determined from Ca isotopes for Pleistocene sediments are higher than those previously inferred from pore fluid Sr concentration and are consistent with rates derived for Late Pleistocene siliciclastic sediments using uranium isotopes. Combining our results for the equilibrium fractionation factor and recrystallization rates, we evaluate the effect of diagenesis on the Ca isotopic composition of marine carbonates at Site 807A. Since calcite precipitation rates in the sedimentary column are many orders of magnitude slower than laboratory experiments and the pore fluids are only slightly oversaturated with respect to calcite, the isotopic composition of diagenetic calcite is likely to reflect equilibrium precipitation. Accordingly, diagenesis produces a maximum shift in δ⁴⁴Ca of +0.15‰ for Site 807A sediments but will have a larger impact where sedimentation rates are low, seawater circulates through the sediment pile, or there are prolonged depositional hiatuses.     

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

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

Natural fractionation of U/U

The isotopic composition of U in nature is generally assumed to be invariant. Here, we report variations of the ²³⁸U/²³⁵U isotope ratio in natural samples (basalts, granites, seawater, corals, black shales, suboxic sediments, ferromanganese crusts/nodules and BIFs) of ∼1.3‰, exceeding by far the analytical precision of our method (≈0.06‰, 2SD). U isotopes were analyzed with MC-ICP-MS using a mixed ²³⁶U-²³³U isotopic tracer (double spike) to correct for isotope fractionation during sample purification and instrumental mass bias. The largest isotope variations found in our survey are between oxidized and reduced depositional environments, with seawater and suboxic sediments falling in between. Light U isotope compositions (relative to SRM-950a) were observed for manganese crusts from the Atlantic and Pacific oceans, which display δ²³⁸U of −0.54‰ to −0.62‰ and for three of four analyzed Banded Iron Formations, which have δ²³⁸U of −0.89‰, −0.72‰ and −0.70‰, respectively. High δ²³⁸U values are observed for black shales from the Black Sea (unit-I and unit-II) and three Kupferschiefer samples (Germany), which display δ²³⁸U of −0.06‰ to +0.43‰. Also, suboxic sediments have slightly elevated δ²³⁸U (−0.41‰ to −0.16‰) compared to seawater, which has δ²³⁸U of −0.41 ± 0.03‰. Granites define a range of δ²³⁸U between −0.20‰ and −0.46‰, but all analyzed basalts are identical within uncertainties and slightly lighter than seawater (δ²³⁸U = −0.29‰).Our findings imply that U isotope fractionation occurs in both oxic (manganese crusts) and suboxic to euxinic environments with opposite directions. In the first case, we hypothesize that this fractionation results from adsorption of U to ferromanganese oxides, as is the case for Mo and possibly Tl isotopes. In the second case, reduction of soluble U^VI to insoluble U^IV probably results in fractionation toward heavy U isotope compositions relative to seawater. These findings imply that variable ocean redox conditions through geological time should result in variations of the seawater U isotope compositions, which may be recorded in sediments or fossils. Thus, U isotopes might be a promising novel geochemical tracer for paleo-redox conditions and the redox evolution on Earth. The discovery that ²³⁸U/²³⁵U varies in nature also has implications for the precision and accuracy of U-Pb dating. The total observed range in U isotope compositions would produce variations in ²⁰⁷Pb/²⁰⁶Pb ages of young U-bearing minerals of up to 3 Ma, and up to 2 Ma for minerals that are 3 billion years old.     

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.     

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

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.     

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

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

Boron isotopic composition of atmospheric precipitations and liquid-vapour fractionations

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

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.     

The behavior of molybdenum and its isotopes across the chemocline and in the sediments of sulfidic Lake Cadagno, Switzerland

Molybdenum (Mo) isotope studies in black shales can provide information about the redox evolution of the Earth’s oceans, provided the isotopic consequences of Mo burial into its major sinks are well understood. Previous applications of the Mo isotope paleo-ocean redox proxy assumed quantitative scavenging of Mo when buried into sulfidic sediments. This paper contains the first complete suite of Mo isotope fractionation observations in a sulfidic water column and sediment system, the meromictic Lake Cadagno, Switzerland, a small alpine lake with a pronounced oxygen-sulfide transition reaching up to H₂S ∼ 200 μM in the bottom waters (or about 300 μM total sulfide: ΣS²⁻ = H₂S + HS⁻ + S²⁻). We find that Mo behaves conservatively in the oxic zone and non-conservatively in the sulfidic zone, where dissolved Mo concentrations decrease from 14 nM to 2-8 nM across this transition. Dissolved Mo in the upper oxic waters has a δ⁹⁸Mo_oxic = 0.9 ± 0.1‰, which matches that of the riverine input, δ⁹⁸Mo_river = 0.9 ± 0.1‰. In the deeper sulfidic waters, a subaquatic source delivers Mo at 1.55 ± 0.1‰, but the dissolved Mo is even heavier at δ⁹⁸Mo_sulfidic = 1.8‰. Sediment traps in the sulfidic zone of the lake collect particles increasingly enriched in Mo with depth, with δ⁹⁸Mo values significantly fractionated at −0.8‰ to −1.2‰ both near the chemocline and in the deepest trap. Suspended particulates in the sulfidic waters carry lighter Mo than the ambient dissolved Mo pool by ∼0.3-1.5‰. Sedimentary Mo concentrations correlate with total organic carbon and yield Mo levels which are two orders of magnitude higher than typical crustal values found in rocks from the catchment area. Solid-phase Mo in the sediment shows a slightly positive δ⁹⁸Mo trend with depth, from δ⁹⁸Mo = 1.2‰ to 1.4‰ while the pore waters show dramatic enrichments of Mo (>2000 nM) with a relatively light isotope signature of δ⁹⁸Mo = 0.9-1.0‰.These data are explained if Mo is converted to particle-reactive oxythiomolybdates in the sulfidic waters and is fractionated during removal from solution onto particles. Isotope fractionation is expressed in the water column, despite the high sulfide concentrations, because the rate of Mo removal is fast compared to the slow reaction kinetics of thiomolybdate formation. However, elemental and isotopic mass balances show that Mo is indeed quantitatively removed to the lake sediments and thus the isotopic composition of the sediments reflects sources to the sulfidic water. This efficient Mo drawdown is expected to occur in settings where H₂S is very much in excess over Mo or in a restricted setting where the water renewal rate is slow compared to the Mo burial rate. We present a model for the Mo isotope fractionation in sulfidic systems associated with the slow reaction kinetics and conclude that quantitative removal will occur in highly sulfidic and restricted marine systems.     

Oxygen and iron isotope constraints on near-surface fractionation effects and the composition of lunar mare basalt source regions

Oxygen and iron isotope analyses of low-Ti and high-Ti mare basalts are presented to constrain their petrogenesis and to assess stable isotope variations within lunar mantle sources. An internally-consistent dataset of oxygen isotope compositions of mare basalts encompasses five types of low-Ti basalts from the Apollo 12 and 15 missions and eight types of high-Ti basalts from the Apollo 11 and 17 missions. High-precision whole-rock δ¹⁸O values (referenced to VSMOW) of low-Ti and high-Ti basalts correlate with major-element compositions (Mg#, TiO₂, Al₂O₃). The observed oxygen isotope variations within low-Ti and high-Ti basalts are consistent with crystal fractionation and match the results of mass-balance models assuming equilibrium crystallization. Whole-rock δ⁵⁶Fe values (referenced to IRMM-014) of high-Ti and low-Ti basalts range from 0.134‰ to 0.217‰ and 0.038‰ to 0.104‰, respectively. Iron isotope compositions of both low-Ti and high-Ti basalts do not correlate with indices of crystal fractionation, possibly owing to small mineral-melt iron fractionation factors anticipated under lunar reducing conditions.The δ¹⁸O and δ⁵⁶Fe values of low-Ti and the least differentiated high-Ti mare basalts are negatively correlated, which reflects their different mantle source characteristics (e.g., the presence or absence of ilmenite). The average δ⁵⁶Fe values of low-Ti basalts (0.073 ± 0.018‰, n = 8) and high-Ti basalts (0.191 ± 0.020‰, n = 7) may directly record that of their parent mantle sources. Oxygen isotope compositions of mantle sources of low-Ti and high-Ti basalts are calculated using existing models of lunar magma ocean crystallization and mixing, the estimated equilibrium mantle olivine δ¹⁸O value, and equilibrium oxygen-fractionation between olivine and other mineral phases. The differences between the calculated whole-rock δ¹⁸O values for source regions, 5.57‰ for low-Ti and 5.30‰ for high-Ti mare basalt mantle source regions, are solely a function of the assumed source mineralogy. The oxygen and iron isotope compositions of lunar upper mantle can be approximated using these mantle source values. The δ¹⁸O and δ⁵⁶Fe values of the lunar upper mantle are estimated to be 5.5 ± 0.2‰ (2σ) and 0.085 ± 0.040‰ (2σ), respectively. The oxygen isotope composition of lunar upper mantle is identical to the current estimate of Earth’s upper mantle (5.5 ± 0.2‰), and the iron isotope composition of the lunar upper mantle overlaps within uncertainty of estimates for the terrestrial upper mantle (0.044 ± 0.030‰).     

Microscale oxygen isotope variations in 1.9 Ga Gunflint cherts: Assessments of diagenesis effects and implications for oceanic paleotemperature reconstructions

Variations in the oxygen isotope composition (δ¹⁸O) of five cherts from the 1.9 Ga Gunflint iron formation (Canada) were studied at the micrometer scale by ion microprobe to try to better understand the processes that control δ¹⁸O values in cherts and to improve seawater paleotemperature reconstructions. Gunflint cherts show clearly different δ¹⁸O values for different types of silica with for instance a difference of ≈15‰ between detrital quartz and microquartz. Microquartz in the five samples is characterized by large intra sample variations in δ¹⁸O values, (δ¹⁸O of quartz varies from 4.6‰ to 6.6‰ at the 20 μm scale and from ≈12‰ to 14‰ at 2 μm scale). Isotopic profiles in microquartz adjacent to hydrothermal quartz veins demonstrate that microquartz more than ≈200 μm away from the veins has preserved its original δ¹⁸O value.At the micrometer spatial resolution of the ion probe, data reveal that microquartz has preserved a considerable δ¹⁸O heterogeneity that must be regarded as a signature inherited from its diagenetic history. Modelling of the δ¹⁸O variations produced during the diagenetic transformation of sedimentary amorphous silica precursors into microquartz allows us to calculate seawater temperature (T_sw at which the amorphous silica precipitated) and diagenesis temperature (T_diagenesis at which microquartz formed) that reproduce the δ¹⁸O distributions (mean, range and shape) measured at micrometer scale in microquartz. The two critical parameters in this modelling are the δ¹⁸O value and the mass fraction of the diagenetic fluid. Under these assumptions, the most likely ranges for T_sw and T_diagenesis are from 37 to 52 °C and from 130 to 170 °C, respectively.     

Sources of fluids and gases expelled at cold seeps offshore Georgia, eastern Black Sea

Four seep sites located within an ∼20 km² area offshore Georgia (Batumi seep area, Pechori Mound, Iberia Mound, and Colkheti Seep) show characteristic differences with respect to element concentrations, and oxygen, hydrogen, strontium, and chlorine isotope signatures in pore waters, as well as impregnation of sediments with petroleum and hydrocarbon potential. All seep sites have active gas seepage, near surface authigenic carbonates and gas hydrates. Cokheti Seep, Iberia Mound, and Pechori Mound are characterized by oil-stained sediments and gas seepage decoupled from deep fluid advection and bottom water intrusion induced by gas bubble release. Pechori Mound is further characterized by deep fluid advection of lower salinity pore fluids. The Pechori Mound pore fluids are altered by mineral/water reactions at elevated temperatures (between 60 and 110 °C) indicated by heavier oxygen and lighter chlorine isotope values, distinct Li and B enrichment, and K depletion. Strontium isotope ratios indicate that fluids originate from late Oligocene strata. This finding is supported by the occurrence of hydrocarbon impregnations within the sediments. Furthermore, light hydrocarbons and high molecular weight impregnates indicate a predominant thermogenic origin for the gas and oil at Pechori Mound, Iberia Mound, and Colkheti Seep. C₁₅₊ hydrocarbons at the oil seeps are allochtonous, whereas those at the Batumi seep area are autochthonous. The presence of oleanane, an angiosperm biomarker, suggests that the hydrocarbon source rocks belong to the Maikopian Formation. In summary, all investigated seep sites show a high hydrocarbon potential and hydrocarbons of Iberia Mound, Colkheti Seep, and Pechori Mound are predominantly of thermogenic origin. However, only at the latter seep site advection of deep pore fluids is indicated.