Combined iron and magnesium isotope geochemistry of pyroxenite xenoliths from Hannuoba,North China Craton: implications for mantle metasomatism

We present high-precision iron and magnesium isotopic data for diverse mantle pyroxenite xenoliths collected from Hannuoba, North China Craton and provide the first combined iron and magnesium isotopic study of such rocks. Compositionally, these xenoliths range from Cr-diopside pyroxenites and Al-augite pyroxenites to garnet-bearing pyroxenites and are taken as physical evidence for different episodes of melt injection. Our results show that both Cr-diopside pyroxenites and Al-augite pyroxenites of cumulate origin display narrow ranges in iron and magnesium isotopic compositions (δ⁵⁷Fe = ?0.01 to 0.09 with an average of 0.03 ± 0.08 (2SD, n = 6); δ²⁶Mg = ? 0.28 to ?0.25 with an average of ?0.26 ± 0.03 (2SD, n = 3), respectively). These values are identical to those in the normal upper mantle and show equilibrium inter-mineral iron and magnesium isotope fractionation between coexisting mantle minerals. In contrast, the garnet-bearing pyroxenites, which are products of reactions between peridotites and silicate melts from an ancient subducted oceanic slab, exhibit larger iron isotopic variations, with δ⁵⁷Fe ranging from 0.12 to 0.30. The δ⁵⁷Fe values of minerals in these garnet-bearing pyroxenites also vary widely (?0.25 to 0.08 in olivines, ?0.04 to 0.25 in orthopyroxenes, ?0.07 to 0.31 in clinopyroxenes, 0.07 to 0.48 in spinels and 0.31–0.42 in garnets). In addition, the garnet-bearing pyroxenite shows light δ²⁶Mg (?0.43) relative to the mantle. The δ²⁶Mg of minerals in the garnet-bearing pyroxenite range from ?0.35 for olivine and orthopyroxene, to ?0.34 for clinopyroxene, 0.04 for spinel and ?0.68 for garnet. These measured values stand in marked contrast to calculated equilibrium iron and magnesium isotope fractionation between coexisting mantle minerals at mantle temperatures derived from theory, indicating disequilibrium isotope fractionation. Notably, one phlogopite clinopyroxenite with an apparent later metasomatic overprint has the heaviest δ⁵⁷Fe (as high as 1.00) but the lightest δ²⁶Mg (as low as ?1.50) values of all investigated samples. Overall, there appears to be a negative co-variation between δ⁵⁷Fe and δ²⁶Mg in the Hannuoba garnet-bearing pyroxenite and in the phlogopite clinopyroxenite xenoliths and minerals therein. These features may reflect kinetic isotopic fractionation due to iron and magnesium inter-diffusion during melt–rock interaction. Such processes play an important role in producing inter-mineral iron and magnesium isotopic disequilibrium and local iron and magnesium isotopic heterogeneity in the subcontinental mantle.     

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.     

Experimental determination of magnesium isotope fractionation during higher plant growth

Two higher plant species (rye grass and clover) were cultivated under laboratory conditions on two substrates (solution, phlogopite) in order to constrain the corresponding Mg isotope fractionations during plant growth and Mg uptake. We show that bulk plants are systematically enriched in heavy isotopes relative to their nutrient source. The Δ²⁶Mg_plant-source range from 0.72‰ to 0.26‰ for rye grass and from 1.05‰ to 0.41‰ for clover. Plants grown on phlogopite display Mg isotope signatures (relative to the Mg source) ∼0.3‰ lower than hydroponic plants. For a given substrate, rye grass display lower δ²⁶Mg (by ∼0.3‰) relative to clover. Magnesium desorbed from rye grass roots display a δ²⁶Mg greater than the nutrient solution. Adsorption experiments on dead and living rye grass roots also indicate a significant enrichment in heavy isotopes of the Mg adsorbed on the root surface. Our results indicate that the key processes responsible for heavy isotope enrichment in plants are located at the root level. Both species also exhibit an enrichment in light isotopes from roots to shoots (Δ²⁶Mg_leaf-root = −0.65‰ and −0.34‰ for rye grass and clover grown on phlogopite respectively, and Δ²⁶Mg_leaf-root of −0.06‰ and −0.22‰ for the same species grown hydroponically). This heavy isotope depletion in leaves can be explained by biological processes that affect leaves and roots differently: (1) organo-Mg complex (including chlorophyll) formation, and (2) Mg transport within plant. For both species, a positive correlation between δ²⁶Mg and K/Mg was observed among the various organs. This correlation is consistent with the link between K and Mg internal cycles, as well as with formation of organo-magnesium compounds associated with enrichment in heavy isotopes. Considering our results together with the published range for δ²⁶Mg of natural plants and rivers, we estimate that a significant change in continental vegetation would induce a change of the mean river δ²⁶Mg that is comparable to analytical uncertainties.     

The major ion, δCa, δCa, and δMg geochemistry of granite weathering at pH = 1 and T = 25 °C: power-law processes and the relative reactivity of minerals

We dissolved Boulder Creek Granodiorite in a plug flow reactor for 5794 h at pH = 1 and T = 25 °C. The primary purpose of the experiment was to identify controls on dissolved δ^44/40Ca, δ^44/42Ca, and δ^26/24Mg values during granite weathering. Herein, we also examine the origin of Ca and Mg isotopic variability among minerals composing the Boulder Creek Granodiorite, and we constrain fundamental characteristics of granite weathering important for quantifying the elemental and isotopic geochemistry of the reactor output. Nine Ca-bearing minerals display an 8.80‰ range of δ^44/40Ca values and a 0.51‰ range of δ^44/42Ca values. Three Mg-bearing minerals display a 1.53‰ range of δ^26/24Mg values. These ranges expressed at the mineralogical scale are higher than the ranges thus far reported for bulk igneous rocks. Most of the δ^44/40Ca variability reflects ⁴⁰Ca enrichment in K-feldspar, and to a lesser extent, biotite, due to the radioactive decay of ⁴⁰K over the 1.7 Ga age of the rock, whereas the entire range of δ^44/42Ca values reflects mass-dependent isotope fractionation during igneous differentiation and crystallization. The range of δ^26/24Mg values may represent either fractionation during the chloritization of biotite or interaction of the Boulder Creek Granodiorite with Mg-rich metamorphic fluids having low δ^26/24Mg values.The elemental and isotopic composition of the reactor output varied substantially during the experiment. We synthesize the mineralogical and fluid data using coupled mass-conservation equations solved at non-steady-state. Model calculations reveal an intricate balance between increasing specific surface area and decreasing mineral concentrations. While surface area normalized dissolution rate constants were time-invariant, specific surface area increased as a power-law function of time through positive feedbacks between mechanical disaggregation, chemical dissolution, and mineral depletion. Variations in dissolved δ^44/40Ca, δ^44/42Ca, and δ^26/24Mg values reflect conservative mixing rather than fractionation. Apatite and calcite initially control δ^44/40Ca and δ^44/42Ca values, followed by biotite, titanite, epidote, hornblende, and plagioclase. The release of radiogenic ⁴⁰Ca clearly defines the period where biotite dissolution dominates. The brucite layer of chlorite initially controls δ^26/24Mg values, followed by biotite, the TOT layer of chlorite, and hornblende. Through direct isotopic tracking, these results demonstrate that trace minerals, such as apatite and calcite in the case of Ca and brucite in the case of Mg, dominate elemental release during the incipient stages of granite weathering. The results further show that biotite dissolution dominates the middle stages of granite weathering and that plagioclase dissolution only becomes important during relatively late stages. The Ca and Mg isotope variations associated with these stages are distinct and potentially resolvable in soil mineral weathering studies.     

Natural variations of δSi ratios during progressive basalt weathering, Hawaiian Islands

Silicon stable isotopes can be used to trace the biogeochemical pathways of Si as it moves from its continental sources to its sink in ocean sediments. Along the way, Si is incorporated into clay minerals, taken up by plants where it forms plant opal, and leached into rivers, the major land-to-ocean conduit. Compared to igneous rocks, the waters that drain continents are enriched in heavy Si isotopes, but the mechanisms that control fractionation have not been elucidated. We studied Si isotope fractionation along a 4 million yr basaltic soil chronosequence on the Hawaiian Islands. Using the natural context of these samples in combination with laboratory experiments, we demonstrate that the isotopic composition of dissolved Si in weathering systems is determined by the combined effects of rock disintegration, clay mineral neosynthesis, and Si biocycling. Weathering preferentially releases ²⁸Si into solution, whereas secondary mineral formation preferentially removes ²⁸Si from solution. In humid environments, leached soils have lost large amounts of this soluble Si, thus creating a net loss of ³⁰Si from the entire soil system. As soils develop and greater fractions of Si reside in neoformed clay minerals, δ³⁰Si_{bulk soil} values change progressively toward more negative values; basalt δ³⁰Si values are about −0.5‰, but older soils have δ³⁰Si values up to −2.5‰. The difference between the solid and solution δ³⁰Si values remains more or less constant with progressive weathering, and therefore, soil water from older soils has a more negative δ³⁰Si composition. In the upper horizons of the Hawaiian soils, this weathering-driven δ³⁰Si shift is modified by the addition of unweathered primary minerals via dust, carrying δ³⁰Si values of about −0.5‰, and by biocycling of Si via plants, producing negative δ³⁰Si values in phytoliths and positive δ³⁰Si values in soil solutions derived from upper horizons. Due to the high concentrations of dissolved Si in these near-surface layers, rivers have more positive δ³⁰Si values than predicted based on the weathering status of the lower horizons. When combined with published δ³⁰Si values from large rivers worldwide, we find that the results from Hawaii point to weathering control of Si isotopes delivered to the oceans, and thus, to an important continent-ocean linkage that warrants further investigation.     

Mg isotope heterogeneity in the Allende meteorite measured by UV laser ablation-MC-ICPMS and comparisons with O isotopes

First results from a new UV laser ablation MC-ICPMS method for measuring Mg isotope ratios in situ in meteoritical materials show that there are mass-dependent variations in δ²⁵Mg and δ²⁶Mg up to 1.5 ‰ per amu in chondrules and 0.3‰ per amu in a CAI from the Allende meteorite. In both cases the mass-dependent fractionation is associated with alteration. Comparisons with laser ablation O isotope data indicate that incorporation of pre-existing grains of forsterite with distinct Mg and O isotopic compositions and post-formation alteration both contributed to the variability in Mg isotope ratios in the chondrules, resulting in a correlation between high δ²⁵Mg and low Δ¹⁷O. The laser ablation analyses of the CAI show that high-precision determinations of both δ²⁵Mg and δ²⁶Mg can be used to discriminate features of the ²⁶Al-²⁶Mg isotope system that are relevant to chronology from those that result from element mobility.     

Magnesium and iron isotopes in 2.7 Ga Alexo komatiites: Mantle signatures, no evidence for Soret diffusion, and identification of diffusive transport in zoned olivine

Komatiites from Alexo, Canada, are well preserved and represent high-degree partial mantle melts (∼50%). They are thus well suited for investigating the Mg and Fe isotopic compositions of the Archean mantle and the conditions of magmatic differentiation in komatiitic lavas. High precision Mg and Fe isotopic analyses of 22 samples taken along a 15-m depth profile in a komatiite flow are reported. The δ²⁵Mg and δ²⁶Mg values of the bulk flow are −0.138 ± 0.021‰ and −0.275 ± 0.042‰, respectively. These values are indistinguishable from those measured in mantle peridotites and chondrites, and represent the best estimate of the composition of the silicate Earth from analysis of volcanic rocks. Excluding the samples affected by secondary Fe mobilization, the δ⁵⁶Fe and δ⁵⁷Fe values of the bulk flow are +0.044 ± 0.030‰, and +0.059 ± 0.044‰, respectively. These values are consistent with a near-chondritic Fe isotopic composition of the silicate Earth and minor fractionation during komatiite magma genesis. In order to explain the early crystallization of pigeonite relative to augite in slowly cooled spinifex lavas, it was suggested that magmas trapped in the crystal mush during spinifex growth differentiated by Soret effect, which should be associated with large and coupled variations in the isotopic compositions of Mg and Fe. The lack of variations in Mg and Fe isotopic ratios either rules out the Soret effect in the komatiite flow or the effect is effaced as the solidification front migrates downward through the flow crust. Olivine separated from a cumulate sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg values relative to the bulk flow, which modeling shows can be explained by kinetic isotope fractionation associated with Fe-Mg inter-diffusion in olivine. Such variations can be used to identify diffusive processes involved in the formation of zoned minerals.     

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.     

Magnesium isotope systematics of the lithologically varied Moselle river basin, France

Magnesium and strontium isotope signatures were determined during different seasons for the main rivers of the Moselle basin, northeastern France. This small basin is remarkable for its well-constrained and varied lithology on a small distance scale, and this is reflected in river water Sr isotope compositions. Upstream, where the Moselle River drains silicate rocks of the Vosges mountains, waters are characterized by relatively high ⁸⁷Sr/⁸⁶Sr ratios (0.7128-0.7174). In contrast, downstream of the city of Epinal where the Moselle River flows through carbonates and evaporites of the Lorraine plateau, ⁸⁷Sr/⁸⁶Sr ratios are lower, down to 0.70824.Magnesium in river waters draining silicates is systematically depleted in heavy isotopes (δ²⁶Mg values range from −1.2 to −0.7‰) relative to the value presently estimated for the continental crust and a local diorite (−0.5‰). In comparison, δ²⁶Mg values measured in soil samples are higher (∼0.0‰). This suggests that Mg isotope fractionation occurs during mineral leaching and/or formation of secondary clay minerals. On the Lorraine plateau, tributaries draining marls, carbonates and evaporites are characterized by low Ca/Mg (1.5-3.2) and low Ca/Sr (80-400) when compared to local carbonate rocks (Ca/Mg = 29-59; Ca/Sr = 370-2200), similar to other rivers draining carbonates. The most likely cause of the Mg and Sr excesses in these rivers is early thermodynamic saturation of groundwater with calcite relative to magnesite and strontianite as groundwater chemistry progressively evolves in the aquifer. δ²⁶Mg of the dissolved phases of tributaries draining mainly carbonates and evaporites are relatively low and constant throughout the year (from −1.4‰ to −1.6‰ and from −1.2‰ to −1.4‰, respectively), within the range defined for the underlying rocks. Downstream of Epinal, the compositions of the Moselle River samples in a δ²⁶Mg vs. ⁸⁷Sr/⁸⁶Sr diagram can be explained by mixing curves between silicate, carbonate and evaporite waters, with a significant contribution from the Vosgian silicate lithologies (>70%). Temporal co-variation between δ²⁶Mg and ⁸⁷Sr/⁸⁶Sr for the Moselle River throughout year is also observed, and is consistent with a higher contribution from the Vosges mountains in winter, in terms of runoff and dissolved element flux. Overall, this study shows that Mg isotopes measured in waters, rocks and soils, coupled with other tracers such as Sr isotopes, could be used to better constrain riverine Mg sources, particularly if analytical uncertainties in Mg isotope measurements can be improved in order to perform more precise quantifications.     

Geochemical fingerprints of dolomitization in Bahamian carbonates: Evidence from sulphur,calcium, magnesium and clumped isotopes

In an effort to constrain the mechanism of dolomitization in Neogene dolomites in the Bahamas and improve understanding of the use of chemostratigraphic tracers in shallow‐water carbonate sediments the δ³⁴S, Δ₄₇, δ¹³C, δ¹⁸O, δ^44/40Ca and δ²⁶Mg values and Sr concentrations have been measured in dolomitized intervals from the Clino core, drilled on the margin of Great Bahama Bank and two other cores (Unda and San Salvador) in the Bahamas. The Unda and San Salvador cores have massively dolomitized intervals that have carbonate associated sulphate δ³⁴S values similar to those found in contemporaneous seawater and δ^44/40Ca, δ²⁶Mg values, Sr contents and Δ₄₇ temperatures (25 to 30°C) indicating relatively shallow dolomitization in a fluid‐buffered system. In contrast, dolomitized intervals in the Clino core have elevated values of carbonate associated sulphate δ³⁴S values indicating dolomitization in a more sediment‐buffered diagenetic system where bacterial sulphate reduction enriches the residual in ³⁴S, consistent with high sediment Sr concentrations and low δ^44/40Ca and high δ²⁶Mg values. Only dolomites associated with hardgrounds in the Clino core have carbonate associated δ³⁴S values similar to seawater, indicating continuous flushing of the upper layers of the sediment by seawater during sedimentary hiatuses. This interpretation is supported by changes to more positive δ^44/40Ca values at hardground surfaces. All dolomites, whether they formed in an open fluid‐buffered or closed sediment‐buffered diagenetic system have similar δ²⁶Mg values suggesting that the HMC transformed to dolomite. The clumped isotope derived temperatures in the dolomitized intervals in Clino yield temperatures that are higher than normal, possibly indicating a kinetic isotope effect on dolomite Δ₄₇ values associated with carbonate formation through bacterial sulphate reduction. The findings of this study highlight the utility of applying multiple geochemical proxies to disentangle the diagenetic history of shallow‐water carbonate sediments and caution against simple interpretations of stratigraphic variability in these geochemical proxies as indicating changes in the global geochemical cycling of these elements in seawater.     

Magnesium isotope fractionation in silicate melts by chemical and thermal diffusion

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

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.     

Iron isotope fractionation during skarn-type metallogeny: A case study of Xinqiao Cu–S–Fe–Au deposit in the Middle–Lower Yangtze valley

Fe isotope compositions of mineral separates and bulk samples from Xinqiao Cu–S–Fe–Au skarn type deposit were investigated. An overall variation in δ⁵⁷Fe values from − 1.22‰ to + 0.73‰ has been observed, which shows some regularity. The δ⁵⁷Fe values of endoskarn and the earliest formed Fe-mineral phase magnetite are ca.1.2‰ and ca. 0.3‰ lower, respectively, relative to the quartz–monzodiorite stock, indicating that fluid exsolved from the stock is enriched in light Fe isotopes. Moreover, spatial and temporal variations in δ⁵⁷Fe values are observed, which suggest iron isotope fractionation during fluid evolution. Precipitation of Fe-bearing minerals results in the Fe isotope composition of residual fluids evolving with time. Precipitation of Fe (III) minerals incorporating heavy iron isotopes preferentially leaves the remaining fluid enriched in light isotopes, while precipitation of Fe (II) minerals preferentially taking-up light iron isotopes, and makes the Fe isotopic composition of the fluid progressively heavier. The regularity of Fe isotope variations occurred during fluid exsolution and evolution indicates that the dominant Fe source of Xinqiao deposit is magmatic. Overall, this study demonstrates that Fe isotope composition has great potential in unraveling ore-forming processes, as well as constraining the metal sources of ore deposits.     

Variations of δSi and Ge/Si with weathering and biogenic input in tropical basaltic ash soils under monoculture

In soils, silicon released by mineral weathering can be retrieved from soil solution through clay formation, Si adsorption onto secondary oxides and plant uptake, thereby impacting the Si-isotopic signature and Ge/Si ratio of dissolved Si (DSi) exported to rivers. Here we use these proxies to study the contribution of biogenic Si (BSi) in a soil-plant system involving basaltic ash soils differing in weathering degree under intensive banana cropping. δ³⁰Si and Ge/Si ratios were determined in bulk soils (<2 mm), sand (50-2000 μm), silt (2-50 μm), amorphous Si (ASi, 2-50 μm) and clay (<2 μm) fractions: δ³⁰Si by MC-ICP-MS Nu Plasma in medium resolution, operating in dry plasma with Mg doping (δ³⁰Si vs. NBS28 ± 0.12‰ ± 2σ_SD), Ge/Si computed after determination of Ge and Si concentrations by HR-ICP-MS and ICP-AES, respectively. Components of the ASi fraction were quantified by microscopic counting (phytoliths, diatoms, ashes). Compared to fresh ash (δ³⁰Si = −0.38‰; Ge/Si = 2.21 μmol mol⁻¹), soil clay fractions (<2 μm) were enriched in light Si isotopes and Ge: with increasing weathering degree, δ³⁰Si decreased from −1.19 to −2.37‰ and Ge/Si increased from 4.10 to 5.25 μmol mol⁻¹. Sand and silt fractions displayed δ³⁰Si values close to fresh ash (−0.33‰) or higher due to saharian dust quartz deposition, whose contribution was evaluated by isotopic mass balance calculation. Si-isotopic signatures of bulk soils (<2 mm) were strongly governed by the relative proportions of primary and secondary minerals: the bulk soil Si-isotopic budget could be closed indicating that all the phases involved were identified. Microscopic counting highlighted a surface accumulation of banana phytoliths and a stable phytolith pool from previous forested vegetation. δ³⁰Si and Ge/Si values of clay fractions in poorly developed volcanic soils, isotopically heavier and Ge-depleted in surface horizons, support the occurrence of a DSi source from banana phytolith dissolution, available for Si sequestration in clay-sized secondary minerals (clay minerals formation and Si adsorption onto Fe-oxide). In the soil-plant system, δ³⁰Si and Ge/Si are thus highly relevant to trace weathering and input of DSi from phytoliths in secondary minerals, although not quantifying the net input of BSi to DSi.     

The behaviour of Li and Mg isotopes during primary phase dissolution and secondary mineral formation in basalt

This study presents lithium (Li) and magnesium (Mg) isotope data from experiments designed to assess the effects of dissolution of primary phases and the formation of secondary minerals during the weathering of basalt. Basalt glass and olivine dissolution experiments were performed in mixed through-flow reactors under controlled equilibrium conditions, at low pH (2-4) in order to keep solutions undersaturated (i.e. far-from equilibrium) and inhibit the formation of secondary minerals. Combined dissolution-precipitation experiments were performed at high pH (10 and 11) increasing the saturation state of the solutions (moving the system closer to equilibrium) and thereby promoting the formation of secondary minerals.At conditions far from equilibrium saturation state modelling and solution stoichiometry suggest that little secondary mineral formation has occurred. This is supported by the similarity of the dissolution rates of basalt glass and olivine obtained here compared to those of previous experiments. The δ⁷Li isotope composition of the experimental solution is indistinguishable from that of the initial basalt glass or olivine indicating that little fractionation has occurred. In contrast, the same experimental solutions have light Mg isotope compositions relative to the primary phases, and the solution becomes progressively lighter with time. In the absence of any evidence for secondary mineral formation the most likely explanation for these light Mg isotope compositions is that there has been preferential loss of light Mg during primary phase dissolution.For the experiments undertaken at close to equilibrium conditions the results of saturation state modelling and changes in solution chemistry suggest that secondary mineral formation has occurred. X-ray diffraction (XRD) measurements of the reacted mineral products from these experiments confirm that the principal secondary phase that has formed is chrysotile. Lithium isotope ratios of the experimental fluid become increasingly heavy with time, consistent with previous experimental work and natural data indicating that ⁶Li is preferentially incorporated into secondary minerals, leaving the solution enriched in ⁷Li. The behaviour of Mg isotopes is different from that anticipated or observed in natural systems. Similar to the far from equilibrium experiments initially light Mg is lost during olivine dissolution, but with time the δ²⁶Mg value of the solution becomes increasingly heavy. This suggests either preferential loss of light, and then heavy Mg from olivine, or that the secondary phase preferentially incorporates light Mg from solution. Assuming that the secondary phase is chrysotile, a Mg-silicate, the sense of Mg fractionation is opposite to that previously associated with silicate soils and implies that the fractionation of Mg isotopes during silicate precipitation may be mineral specific. If secondary silicates do preferentially remove light Mg from solution then this could be a possible mechanism for the relatively heavy δ²⁶Mg value of seawater. This study highlights the utility of experimental studies to quantify the effects of natural weathering reactions on the Li and Mg geochemical cycles.     

Tracing mechanisms controlling the release of dissolved silicon in forest soil solutions using Si isotopes and Ge/Si ratios

The terrestrial biogenic Si (BSi) pool in the soil-plant system is ubiquitous and substantial, likely impacting the land-ocean transfer of dissolved Si (DSi). Here, we consider the mechanisms controlling DSi in forest soil in a temperate granitic ecosystem that would differ from previous works mostly focused on tropical environments. This study aims at tracing the source of DSi in forest floor leachates and in soil solutions under various tree species at homogeneous soil and climate conditions, using stable Si isotopes and Ge/Si ratios. Relative to granitic bedrock, clays minerals were enriched in ²⁸Si and had high Ge/Si ratios, while BSi from phytoliths was also enriched in ²⁸Si, but had a low Ge/Si ratio. Such a contrast is useful to infer the relative contribution of silicate weathering and BSi dissolution in the shallow soil on the release of DSi in forest floor leachate solutions. The δ³⁰Si values in forest floor leachates (−1.38‰ to −2.05‰) are the lightest ever found in natural waters, and Ge/Si ratios are higher in forest floor leachates relative to soil solution. These results suggest dissolution of ²⁸Si and Ge-enriched secondary clay minerals incorporated by bioturbation in organic-rich horizons in combination with an isotopic fractionation releasing preferentially light Si isotopes during this dissolution process. Ge/Si ratios in soil solutions are governed by incongruent weathering of primary minerals and neoformation of secondary clays minerals. Tree species influence Si-isotopic compositions and Ge/Si ratios in forest floor leachates through differing incorporation of minerals in organic horizons by bioturbation and, to a lesser extent, through differing Si recycling.     

Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments

The chemical and isotopic composition of speleothem calcite and particularly that of stalagmites and flowstones is increasingly exploited as an archive of past environmental change in continental settings. Despite intensive research, including modelling and novel approaches, speleothem data remain difficult to interpret. A possible way foreword is to apply a multi-proxy approach including non-conventional isotope systems. For the first time, we here present a complete analytical dataset of magnesium isotopes (δ²⁶Mg) from a monitored cave in NW Germany (Bunker Cave). The data set includes δ²⁶Mg values of loess-derived soil above the cave (−1.0 ± 0.5‰), soil water (−1.2 ± 0.5‰), the carbonate hostrock (−3.8 ± 0.5‰), dripwater in the cave (−1.8 ± 0.2‰), speleothem low-Mg calcite (stalactites, stalagmites; −4.3 ± 0.6‰), cave loam (−0.6 ± 0.1‰) and runoff water (−1.8 ± 0.1‰) in the cave, respectively. Magnesium-isotope fractionation processes during weathering and interaction between soil cover, hostrock and solute-bearing soil water are non-trivial and depend on a number of variables including solution residence times, dissolution rates, adsorption effects and potential neo-formation of solids in the regolith and the carbonate aquifer. Apparent Mg-isotope fractionation between dripwater and speleothem low-Mg calcite is about 1000lnα_Mg-cc-Mg(aq) = −2.4‰. A similar Mg-isotope fractionation (1000lnα_Mg-cc-Mg(aq) ≈ −2.1‰) is obtained by abiogenic precipitation experiments carried out at aqueous Mg/Ca ratios and temperatures close to cave conditions. Accordingly, ²⁶Mg discrimination during low-Mg calcite formation in caves is highly related to inorganic fractionation effects, which may comprise dehydration of Mg²⁺ prior to incorporation into calcite, surface entrapment of light isotopes and reaction kinetics. Relevance of kinetics is supported by a significant negative correlation of Mg-isotope fractionation with the precipitation rate for inorganic precipitation experiments.     

CaO-MgO-Al2o3-SiO2 liquids: chemical and isotopic effects of Mg and Si evaporation in a closed system of solar composition

A method is shown for calculating vapor pressures over a CMAS droplet in a gas of any composition. It is applied to the problem of the evolution of the chemical and Mg and Si isotopic composition of a completely molten droplet having the composition of a likely refractory inclusion precursor during its evaporation into the complementary, i.e. modified solar, gas from which it originally condensed, a more realistic model than previous calculations in which the ambient gas is pure H_2(g). Because the loss rate of Mg is greater than that of Si, the vapor pressure of Mg_(g) falls and its ambient pressure rises faster than those of SiO_(g) during isothermal evaporation, causing the flux of Mg_(g) to approach zero faster and MgO to approach its equilibrium concentration sooner than SiO₂. As time passes, δ²⁵Mg and δ²⁹Si increase in the droplet and decrease in the ambient gas. The net flux of each isotope crossing the droplet/gas interface is the difference between its outgoing and incoming flux. δ²⁵Mg and δ²⁹Si of this instantaneous gas become higher, first overtaking their values in the ambient gas, causing them to increase with time, and later overtaking their values in the droplet itself, causing them to decrease with time, ultimately reaching their equilibrium values. If the system is cooling during evaporation and if mass transfer ceases at the solidus temperature, 1500 K, final MgO and SiO₂ contents of the droplet are slightly higher in modified solar gas than in pure H_2(g), and the difference increases with decreasing cooling rate and increasing ambient pressure. During cooling under some conditions, net fluxes of evaporating species become negative, causing reversal of the evaporation process into a condensation process, an increase in the MgO and/or SiO₂ content of the droplet with time, and an increase in their final concentrations with increasing ambient pressure and/or dust/gas ratio. At cooling rates <∼3 K/h, closed-system evaporation at P^tot ∼ 10⁻³ bar in a modified solar gas, or at lower pressure in systems with enhanced dust/gas ratio, can yield the same δ²⁵Mg in a residual CMAS droplet for vastly different evaporated fractions of Mg. The δ²⁵Mg of a refractory residue may thus be insufficient to determine the extent of Mg loss from its precursor. Evaporation of Mg into an Mg-bearing ambient gas causes δ²⁶Mg and δ²⁵Mg of the residual droplet to fall below values expected from Rayleigh fractionation for the amount of ²⁴Mg evaporated, with the degree of departure increasing with increasing fraction evaporated and ambient pressure of Mg. δ²⁶Mg and δ²⁵Mg do not depart proportionately from Rayleigh fractionation curves, with δ²⁵Mg being less than expected on the basis of δ²⁶Mg by up to ∼1.2‰. Such departures from Rayleigh fractionation could be used in principle to distinguish heavily from lightly evaporated residues with the same δ²⁵Mg.     

Cu isotopic fractionation in the supergene environment with and without bacteria

The isotopic composition of dissolved Cu and solid Cu-rich minerals [δ⁶⁵Cu (‰) = (⁶⁵Cu/⁶³Cu_sample/⁶⁵Cu/⁶³Cu_std) - 1)*1000] were monitored in batch oxidative dissolution experiments with and without Thiobacillus ferrooxidans. Aqueous copper in leach fluids released during abiotic oxidation of both chalcocite and chalcopyrite was isotopically heavier (δ⁶⁵Cu = 5.34‰ and δ⁶⁵Cu = 1.90‰, respectively, [±0.16 at 2σ]) than the initial starting material (δ⁶⁵Cu = 2.60 ± 0.16‰ and δ⁶⁵Cu = 0.58 ± 0.16‰, respectively). Isotopic mass balance between the starting material, aqueous copper, and secondary minerals precipitated in these experiments explains the heavier isotopic values of aqueous copper. In contrast, aqueous copper from leached chalcocite and chalcopyrite inoculated with Thiobacillus ferrooxidans was isotopically similar to the starting material. The lack of fractionation of the aqueous copper in the biotic experiments can best be explained by assuming a sink for isotopically heavy copper present in the bacteria cells with δ⁶⁵Cu = 5.59 ± 0.16‰. Consistent with this inference, amorphous Cu-Fe oxide minerals are observed surrounding cell membranes of Thiobacillus grown in the presence of dissolved Cu and Fe.Extrapolating these experiments to natural supergene environments implies that release of isotopically heavy aqueous Cu from oxidative leach caps, especially under abiotic conditions, should result in precipitates in underlying enrichment blankets that are isotopically heavy. Where iron-oxidizing cells are involved, isotopically heavy oxidized Cu entrained in cellular material may become associated with leach caps, causing the released aqueous Cu to be less isotopically enriched in the heavy isotope than predicted for the abiotic system. Rayleigh fractionation trends with fractionation factors calculated from our experiments for both biotic and abiotic conditions are consistent with large numbers of individual abiotic or biotic leaching events, explaining the supergene chalcocites in the Morenci and Silver Bell porphyry copper deposits.     

Early diagenetic processes during Mn-carbonate formation: evidence from the isotopic composition of authigenic Ca-rhodochrosites of the Baltic Sea

The formation of authigenic Ca-rich rhodochrosite (ACR) in sapropelic sediments of the Gotland Basin, Baltic Sea, is governed by deepwater renewal processes whereby saline water from the North Atlantic flushes the brackish anoxic Baltic Deeps. The carbon and oxygen isotopic compositions of these Mn-carbonates suggest that ACR formation takes place just below the sediment surface and that dissolved compounds from the deepwater column, such as water and bicarbonate molecules, were incorporated in ACR during authigenesis. Porewaters near the sediment surface display δ¹⁸O values of −5.4‰ (VSMOW) and are generally depleted in ¹⁸O, compared to the oxygen isotopic composition of water in equilibrium with Mn-carbonate solid solutions (ACR δ¹⁸O values are −4.6‰). This suggests that early burial diagenetic processes significantly modify the initial isotopic composition of water during Mn-carbonate formation. The reduction of sulfate having δ¹⁸O values of +8.4‰ accounts for a permanent enrichment of porewater ¹⁸O and observed δ¹⁸O values at depth equal to −4.6‰. However, this process does not explain the observed disequilibrium in the oxygen isotopic composition between water and ACR close to the sediment surface where Mn-carbonate formation takes place. Based on isotopic mass balance calculations, we suggest that MnO₂ with δ¹⁸O values of +8.9‰ released oxygen enriched in ¹⁸O into the anoxic porewaters close below the sediment surface. This process should occur after oxygenation events during deepwater renewal when MnO₂ accumulates at the surface of anoxic sediments. Manganese carbonates formed in these waters display δ¹⁸O values of ∼1.0‰ heavier than values expected solely from the initial deepwater composition. This quantitatively explains the discrepancy between paleosalinities calculated from ACR δ¹⁸O based on Mn-carbonate/water isotopic equilibrium fractionation and direct observations for the same period. Our results emphasize the important role of microbial MnO₂ reduction during rhodochrosite authigenesis and suggest that Mn(II) activity, rather than alkalinity, is the limiting component for sedimentary Mn-carbonate formation.