Fractionation of Cu and Zn isotopes during adsorption onto amorphous Fe(III) oxyhydroxide: Experimental mixing of acid rock drainage and ambient river water

Fractionation of Cu and Zn isotopes during adsorption onto amorphous ferric oxyhydroxide is examined in experimental mixtures of metal-rich acid rock drainage and relatively pure river water and during batch adsorption experiments using synthetic ferrihydrite. A diverse set of Cu- and Zn-bearing solutions was examined, including natural waters, complex synthetic acid rock drainage, and simple NaNO₃ electrolyte. Metal adsorption data are combined with isotopic measurements of dissolved Cu (⁶⁵Cu/⁶³Cu) and Zn (⁶⁶Zn/⁶⁴Zn) in each of the experiments. Fractionation of Cu and Zn isotopes occurs during adsorption of the metal onto amorphous ferric oxyhydroxide. The adsorption data are modeled successfully using the diffuse double layer model in PHREEQC. The isotopic data are best described by a closed system, equilibrium exchange model. The fractionation factors (α_soln-solid) are 0.99927 ± 0.00008 for Cu and 0.99948 ± 0.00004 for Zn or, alternately, the separation factors (Δ_soln-solid) are −0.73 ± 0.08‰ for Cu and −0.52 ± 0.04‰ for Zn. These factors indicate that the heavier isotope preferentially adsorbs onto the oxyhydroxide surface, which is consistent with shorter metal-oxygen bonds and lower coordination number for the metal at the surface relative to the aqueous ion. Fractionation of Cu isotopes also is greater than that for Zn isotopes. Limited isotopic data for adsorption of Cu, Fe(II), and Zn onto amorphous ferric oxyhydroxide suggest that isotopic fractionation is related to the intrinsic equilibrium constants that define aqueous metal interactions with oxyhydroxide surface sites. Greater isotopic fractionation occurs with stronger metal binding by the oxyhydroxide with Cu > Zn > Fe(II).     

Isotopic variations of dissolved copper and zinc in stream waters affected by historical mining

Zinc and Cu play important roles in the biogeochemistry of natural systems, and it is likely that these interactions result in mass-dependent fractionations of their stable isotopes. In this study, we examine the relative abundances of dissolved Zn and Cu isotopes in a variety of stream waters draining six historical mining districts located in the United States and Europe. Our goals were to (1) determine whether streams from different geologic settings have unique or similar Zn and Cu isotopic signatures and (2) to determine whether Zn and Cu isotopic signatures change in response to changes in dissolved metal concentrations over well-defined diel (24-h) cycles.Average δ⁶⁶Zn and δ⁶⁵Cu values for streams varied from +0.02‰ to +0.46‰ and −0.7‰ to +1.4‰, respectively, demonstrating that Zn and Cu isotopes are heterogeneous among the measured streams. Zinc or Cu isotopic changes were not detected within the resolution of our measurements over diel cycles for most streams. However, diel changes in Zn isotopes were recorded in one stream where the fluctuations of dissolved Zn were the largest. We calculate an apparent separation factor of ∼0.3‰ (^66/64Zn) between the dissolved and solid Zn reservoirs in this stream with the solid taking up the lighter Zn isotope. The preference of the lighter isotope in the solid reservoir may reflect metabolic uptake of Zn by microorganisms. Additional field investigations must evaluate the contributions of soils, rocks, minerals, and anthropogenic components to Cu and Zn isotopic fluxes in natural waters. Moreover, rigorous experimental work is necessary to quantify fractionation factors for the biogeochemical reactions that are likely to impact Cu and Zn isotopes in hydrologic systems. This initial investigation of Cu and Zn isotopes in stream waters suggests that these isotopes may be powerful tools for probing biogeochemical processes in surface waters on a variety of temporal and spatial scales.     

Zn and Cu isotopic variations in chondrites and iron meteorites: Early solar nebula reservoirs and parent-body processes

High-precision Zn isotopic variations are reported for carbonaceous chondrites (CC), equilibrated (EOC) and unequilibrated (UOC) ordinary chondrites, iron meteorites from the IAB-IIICD (nonmagmatic) and IIIA (magmatic) groups, and metal from the Brenham pallasite. For irons, δ⁶⁵Cu values are also reported. Data have also been obtained on a coarse-grained type-B calcium-, aluminum-rich refractory inclusion (CAI) from Allende and on acid leaches of Allende (CV3), Krymka (LL3), and Charsonville (H6). Variations expressed as δ⁶⁶Zn (deviation in parts per thousand of ⁶⁶Zn/⁶⁴Zn in samples relative to a standard) spread over a range of 0.3‰ for carbonaceous chondrites, 2‰ for ordinary chondrites, and 4‰ for irons.The measured ⁶⁶Zn/⁶⁴Zn, ⁶⁷Zn/⁶⁴Zn, and ⁶⁸Zn/⁶⁴Zn ratios vary linearly with mass difference and define a common isotope fractionation line with terrestrial samples, which demonstrates that Zn was derived from an initially single homogeneous reservoir. The δ⁶⁶Zn values are correlated with meteorite compositions and slightly decrease in the order CI, CM, CV-CO, and to UOC. The isotopically light Zn of Allende CAI and the acid-resistant residues of Allende and Krymka show that the light component is associated with refractory material, presumably minerals from the spinel-group. This, together with the reverse correlation between relative abundances of light Zn isotopes and volatile element abundances, suggests that Zn depletion in planetary bodies with respect to CI cannot be ascribed to devolatilization of CI-like material. These observations rather suggest that refractory material reacted with a gas phase enriched in the lighter Zn isotopes. Alternatively, chondrules with their associated rims should carry a light Zn isotopic signature. The δ⁶⁶Zn values of unequilibrated chondrites are rather uniform, whereas equilibrated chondrites show distinctly more isotopic variability.The values of δ⁶⁵Cu-δ⁶⁶Zn in irons define two trends. The moderate and positively correlated Cu and Zn isotope variations in IIIA and pallasite samples probably reflect crystallization of silicate, sulfide, and solid metal from the liquid metal. The range of δ⁶⁶Zn values of the IAB-IIICD group is large (>3‰) and contrasts with the moderate fractionation of Cu isotopes. We interpret this feature and the negative δ⁶⁶Zn-δ⁶⁵Cu correlation as reflecting mixing, possibly achieved by percolation, between metals from a regolith devolatilized at low temperature (enriched in heavy zinc) and metallic liquids formed within the parent body.     

Atmospheric deposition and isotope biogeochemistry of zinc in ombrotrophic peat

Zinc isotope ratios were measured in the top sections of dated ombrotrophic peat cores in Finland to investigate their potential as proxies for atmospheric sources and to constrain post depositional processes affecting the geochemical record. The peat deposits were located in Hietajärvi, a background site well away from any point pollution source and representing ‘background’ conditions, in Outokumpu, next to a mining site, and in Harjavalta, next to a smelter. Measured total concentrations, calculated excess concentrations and mass balance considerations suggest that zinc is subjected to important biogeochemical cycling within the peat. Significant isotopic variability was found in all three peat bogs, with heavier zinc in the deeper and lighter zinc in the upper sections. Isotope ratios and concentrations correlated in the two peats located next to dominant point sources, i.e. the smelting and mining site, suggesting that zinc isotopes trace pollution sources. Concentration and isotope peaks were offset from the period of mining and smelting activity, supporting migration of zinc down the profile. The δ⁶⁶Zn_JMC (where δ⁶⁶Zn = [(⁶⁶Zn/⁶⁴Zn)_sample/(⁶⁶Zn/⁶⁴Zn)_JMC-standard − 1] × 10³) of the top section sample at the remote Hietajärvi site was 0.9‰ and we suggest this represents the regional background isotope signature of atmospheric zinc. The deeper sections of the peat cores show isotopically heavier zinc than any potential atmospheric source, indicating that post depositional processes affected the isotopic records. The large variations encountered (up to 1.05‰ for δ⁶⁶Zn) and Rayleigh modelling imply that multiple fractionation of zinc during diagenetic alterations occurs and nutrient recycling alone cannot explain the fractionation pattern.We propose that zinc isotopes are amenable to identify different atmospheric zinc sources, including zinc derived from anthropogenic activities such as mining and smelting, but multiple biogeochemical processes seriously affect the record and they need to be evaluated and assessed carefully if zinc isotopes are used in terrestrial paleorecords.     

Ion-exchange fractionation of copper and zinc isotopes

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

Mo isotope fractionation during adsorption to Fe (oxyhydr)oxides

Molybdenum (Mo) isotopes have great potential as a paleoredox indicator, but this potential is currently restricted by an incomplete understanding of isotope fractionations occurring during key (bio)geochemical processes. To address one such uncertainty we have investigated the isotopic fractionation of Mo during adsorption to a range of Fe (oxyhydr)oxides, under variable Mo/Fe-mineral ratios and pH. Our data confirm that Fe (oxyhydr)oxides can readily adsorb Mo, highlighting the potential importance of this removal pathway for the global Mo cycle. Furthermore, adsorption of Mo to Fe (oxyhydr)oxides is associated with preferential uptake of the lighter Mo isotopes. Fractionations between the solid and dissolved phase (Δ⁹⁸Mo) increase at higher pH, and also vary with mineralogy, increasing in the order magnetite (Δ⁹⁸Mo = 0.83 ± 0.60‰) < ferrihydrite (Δ⁹⁸Mo = 1.11 ± 0.15‰) < goethite (Δ⁹⁸Mo = 1.40 ± 0.48‰) < hematite (Δ⁹⁸Mo = 2.19 ± 0.54‰). Small differences in isotopic fractionation are also seen at varying Mo/Fe-mineral ratios for individual minerals. The observed isotopic behaviour is consistent with both fractionation during adsorption to the mineral surface (a function of vibrational energy) and adsorption of different Mo species/structures from solution. The different fractionation factors determined for different Fe (oxyhydr)oxides suggests that these minerals likely exert a major control on observed natural Mo isotope compositions during sediment deposition beneath suboxic through to anoxic (but non-sulfidic) bottom waters. Our results confirm that Mo isotopes can provide important information on the spatial extent of different paleoredox conditions, providing they are used in combination with other techniques for evaluating the local redox environment and the mineralogy of the depositing sediments.     

Contrasting isotopic signatures between anthropogenic and geogenic Zn and evidence for post-depositional fractionation processes in smelter-impacted soils from Northern France

Zinc isotopes have been studied along two smelter-impacted soil profiles sampled near one of the largest Pb and Zn processing plants in Europe located in northern France, about 50 km south of Lille. Mean δ⁶⁶Zn values along these two soil profiles range from +0.22 ± 0.17‰ (2σ) to +0.34 ± 0.17‰ (2σ) at the lowest horizons and from +0.38 ± 0.45‰ (2σ) to +0.76 ± 0.14‰ (2σ) near the surface. The δ⁶⁶Zn values in the lowest horizons of the soils are interpreted as being representative of the local geochemical background (mean value +0.31 ± 0.38‰), whereas heavier δ⁶⁶Zn values near the surface of the two soils are related to anthropogenic Zn. This anthropogenic Zn occurs in the form of franklinite (ZnFe₂O₄)-bearing slag grains originating from processing wastes at the smelter site and exhibiting δ⁶⁶Zn values of +0.81 ± 0.20‰ (2σ). The presence of franklinite is indicated by EXAFS analysis of the topsoil samples from both soil profiles as well as by micro-XANES analysis of the surface horizon of a third smelter-impacted soil from a distant site. These results indicate that naturally occurring Zn and smelter-derived Zn exhibit significantly different δ⁶⁶Zn values, which suggests that zinc isotopes can be used to distinguish between geogenic and anthropogenic sources of Zn in smelter-impacted soils. In addition to a possible influence of additional past sources of light Zn (likely Zn-sulfides and Zn-sulfates directly emitted by the smelter), the light δ⁶⁶Zn values in the surface horizons compared to smelter-derived slag materials are interpreted as resulting mainly from fractionation processes associated with biotic and/or abiotic pedological processes (Zn-bearing mineral precipitation, Zn complexation by organic matter, and plant uptake of Zn). This conclusion emphasizes the need for additional Zn isotopic studies before being able to use Zn isotopes to trace sources and pathways of this element in surface environments.     

Iron isotope fractionation and atom exchange during sorption of ferrous iron to mineral surfaces

The application of stable Fe isotopes as a tracer of the biogeochemical Fe cycle necessitates a mechanistic knowledge of natural fractionation processes. We studied the equilibrium Fe isotope fractionation upon sorption of Fe(II) to aluminum oxide (γ-Al₂O₃), goethite (α-FeOOH), quartz (α-SiO₂), and goethite-loaded quartz in batch experiments, and performed continuous-flow column experiments to study the extent of equilibrium and kinetic Fe isotope fractionation during reactive transport of Fe(II) through pure and goethite-loaded quartz sand. In addition, batch and column experiments were used to quantify the coupled electron transfer-atom exchange between dissolved Fe(II) (Fe(II)_aq) and structural Fe(III) of goethite. All experiments were conducted under strictly anoxic conditions at pH 7.2 in 20 mM MOPS (3-(N-morpholino)-propanesulfonic acid) buffer and 23 °C. Iron isotope ratios were measured by high-resolution MC-ICP-MS. Isotope data were analyzed with isotope fractionation models. In batch systems, we observed significant Fe isotope fractionation upon equilibrium sorption of Fe(II) to all sorbents tested, except for aluminum oxide. The equilibrium enrichment factor, , of the Fe(II)_sorb-Fe(II)_aq couple was 0.85 ± 0.10‰ (±2σ) for quartz and 0.85 ± 0.08‰ (±2σ) for goethite-loaded quartz. In the goethite system, the sorption-induced isotope fractionation was superimposed by atom exchange, leading to a δ^56/54Fe shift in solution towards the isotopic composition of the goethite. Without consideration of atom exchange, the equilibrium enrichment factor was 2.01 ± 0.08‰ (±2σ), but decreased to 0.73 ± 0.24‰ (±2σ) when atom exchange was taken into account. The amount of structural Fe in goethite that equilibrated isotopically with Fe(II)_aq via atom exchange was equivalent to one atomic Fe layer of the mineral surface (∼3% of goethite-Fe). Column experiments showed significant Fe isotope fractionation with δ^56/54Fe(II)_aq spanning a range of 1.00‰ and 1.65‰ for pure and goethite-loaded quartz, respectively. Reactive transport of Fe(II) under non-steady state conditions led to complex, non-monotonous Fe isotope trends that could be explained by a combination of kinetic and equilibrium isotope enrichment factors. Our results demonstrate that in abiotic anoxic systems with near-neutral pH, sorption of Fe(II) to mineral surfaces, even to supposedly non-reactive minerals such as quartz, induces significant Fe isotope fractionation. Therefore we expect Fe isotope signatures in natural systems with changing concentration gradients of Fe(II)_aq to be affected by sorption.     

Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed

Here the hydrogeochemical constraints of a tracer dilution study are combined with Fe and Zn isotopic measurements to pinpoint metal loading sources and attenuation mechanisms in an alpine watershed impacted by acid mine drainage. In the tested mountain catchment, δ⁵⁶Fe and δ⁶⁶Zn isotopic signatures of filtered stream water samples varied by ∼3.5‰ and 0.4‰, respectively. The inherent differences in the aqueous geochemistry of Fe and Zn provided complimentary isotopic information. For example, variations in δ⁵⁶Fe were linked to redox and precipitation reactions occurring in the stream, while changes in δ⁶⁶Zn were indicative of conservative mixing of different Zn sources. Fen environments contributed distinctively light dissolved Fe (<−2.0‰) and isotopically heavy suspended Fe precipitates to the watershed, while Zn from the fen was isotopically heavy (>+0.4‰). Acidic drainage from mine wastes contributed heavier dissolved Fe (∼+0.5‰) and lighter Zn (∼+0.2‰) isotopes relative to the fen. Upwelling of Fe-rich groundwater near the mouth of the catchment was the major source of Fe (δ⁵⁶Fe ∼ 0‰) leaving the watershed in surface flow, while runoff from mining wastes was the major source of Zn. The results suggest that given a strong framework for interpretation, Fe and Zn isotopes are useful tools for identifying and tracking metal sources and attenuation mechanisms in mountain watersheds.     

Interaction between zinc and freshwater and marine diatom species: Surface complexation and Zn isotope fractionation

This work is devoted to characterization of zinc interaction in aqueous solution with two marine planktonic (Thalassiosira weissflogii = TW, Skeletonema costatum = SC) and two freshwater periphytic species (Achnanthidium minutissimum = AMIN, Navicula minima = NMIN) by combining adsorption and electrophoretic measurements with surface complexation modeling and by assessing Zn isotopes fractionation during both long term uptake and short term adsorption on diatom cells and their frustules. Reversible adsorption experiments were performed at 25 and 5 °C as a function of exposure time (5 min to 140 h), pH (2 to 10), zinc concentration in solution (10 nM to 1 mM), ionic strength (I = 0.001 to 1.0 M) and the presence of light. While the shape of pH-dependent adsorption edge is almost the same for all four species, the constant-pH adsorption isotherm and maximal Zn binding capacities differ by an order of magnitude. The extent of adsorption increases with temperature from 5 to 25 °C and does not depend on light intensity. Zinc adsorption decreases with increase of ionic strength suggesting competition with sodium for surface sites. Cell number-normalized concentrations of sorbed zinc on whole cells and their silica frustules demonstrated only weak contribution of the latter (10-20%) to overall zinc binding by diatom cell wall. Measurements of electrophoretic mobilities (μ) revealed negative diatoms surface potential in the full range of zinc concentrations investigated (0.15-760 μmol/L), however, the absolute value of μ decreases at [Zn] > 15 μmol/L suggesting a change in surface speciation. These observations allowed us to construct a surface complexation model for Zn binding by diatom surfaces that postulates the constant capacitance of the electric double layer and considers Zn complexation with carboxylate and silanol groups. Thermodynamic and structural parameters of this model are based on previous acid-base titration and spectroscopic results and allow quantitative reproduction of all adsorption experiments. Although Zn adsorption constants on carboxylate groups are almost the same, Zn surface adsorption capacities are very different among diatom species which is related to the systematic differences in their cell wall composition and thickness. Measurements of Zn isotopic composition (⁶⁶Zn/(⁶⁴Zn)) performed using a multicollector ICP MS demonstrated that irreversible incorporation of Zn in cultured diatom cells produces enrichment in heavy isotope compared to growth media (Δ⁶⁶Zn(solid-solution) = 0.27 ± 0.05, 0.08 ± 0.05, 0.21 ± 0.05, and 0.19 ± 0.05‰ for TW, SC, NMIN, and AMIN species, respectively). Accordingly, an enrichment of cells in heavy isotopes (Δ⁶⁶Zn(solid-solution) = 0.43 ± 0.1 and 0.27 ± 0.1‰ for NMIN and AMIN, respectively) is observed following short-term Zn sorption on freshwater cells in nutrient media at pH ∼ 7.8. Finally, diatoms frustules are enriched in heavy isotopes compared to solution during Zn adsorption on silica shells at pH ∼ 5.5 (Δ⁶⁶Zn(solid-solution) = 0.35 ± 0.10‰). Measured isotopes fractionation can be related to the structure and stability of Zn complexes formed and they provide a firm basis for using Zn isotopes for biogeochemical tracing.     

Isotopic fractionation of Cu in tektites

Tektites are terrestrial natural glasses of up to a few centimeters in size that were produced during hypervelocity impacts on the Earth’s surface. It is well established that the chemical and isotopic composition of tektites is generally identical to that of the upper terrestrial continental crust. Tektites typically have very low water content, which has generally been explained by volatilization at high temperature; however, the exact mechanism is still debated. Because volatilization can fractionate isotopes, comparing the isotopic composition of volatile elements in tektites with those of their source rocks may help to understand the physical conditions during tektite formation.Interestingly, volatile chalcophile elements (e.g., Cd and Zn) seem to be the only elements for which isotopic fractionation is known so far in tektites. Here, we extend this study to Cu, another volatile chalcophile element. We have measured the Cu isotopic composition for 20 tektite samples from the four known different strewn fields. All of the tektites (except the Muong Nong-types) are enriched in the heavy isotopes of Cu (1.98 < δ⁶⁵Cu < 6.99) in comparison to the terrestrial crust (δ⁶⁵Cu ≈ 0) with no clear distinction between the different groups. The Muong Nong-type tektites and a Libyan Desert Glass sample are not fractionated (δ⁶⁵Cu ≈ 0) in comparison to the terrestrial crust. To refine the Cu isotopic composition of the terrestrial crust, we also present data for three geological reference materials (δ⁶⁵Cu ≈ 0).An increase of δ⁶⁵Cu with decreasing Cu abundance probably reflects that the isotopic fractionation occurred by evaporation during heating. A simple Rayleigh distillation cannot explain the Cu isotopic data and we suggest that the isotopic fractionation is governed by a diffusion-limited regime. Copper is isotopically more fractionated than the more volatile element Zn (δ^66/64Zn up to 2.49‰). This difference of behavior between Cu and Zn is predicted in a diffusion-limited regime, where the magnitude of the isotopic fractionation is regulated by the competition between the evaporative flux and the diffusive flux at the diffusion boundary layer. Due to the difference of ionic charge in silicates (Zn²⁺ vs. Cu⁺), Cu has a diffusion coefficient that is larger than that of Zn by at least two orders of magnitude. Therefore, the larger isotopic fractionation in Cu than in Zn in tektites is due to the significant difference in their respective chemical diffusivity.     

Kinetic fractionation of Fe isotopes during transport through a porous quartz-sand column

Sorption and desorption processes are an important part of biological and geochemical metallic isotope cycles. Here, we address the dynamic aspects of metallic isotopic fractionation in a theoretical and experimental study of Fe sorption and desorption during the transport of aqueous Fe(III) through a quartz-sand matrix. Transport equations describing the behavior of sorbing isotopic species in a water saturated homogeneous porous medium are presented; isotopic fractionation of the system (Δ_{sorbedmetal-soln}) being defined in terms of two parameters: (i) an equilibrium fractionation factor, α_e; and (ii) a kinetic sorption factor, α₁. These equations are applied in a numerical model that simulates the sorption-desorption of Fe isotopes during injection of a Fe(III) solution pulse into a quartz matrix at pH 0-2 and explores the effects of the kinetic and equilibrium parameters on the Fe-isotope evolution of porewater. The kinetic transport theory is applied to a series of experiments in which pulses of Na and Fe(III) chloride solutions were injected into a porous sand grain column. Fractionation factors of α_e = 1.0003 ± 0.0001 and α₁ = 0.9997 ± 0.0004 yielded the best fit between the transport model and the Fe concentration and δ⁵⁶Fe data. The equilibrium fractionation (Δ⁵⁶Fe_{sorbedFe-soln}) of 0.3‰ is comparable with values deduced for adsorption of metallic cations on iron and manganese oxide surfaces and suggests that sandstone aquifers will fractionate metallic isotopes during sorption-desorption reactions. The ability of the equilibrium fractionation factor to describe a natural system, however, depends on the proximity to equilibrium, which is determined by the relative time scales of mass transfer and chemical reaction; low fluid transport rates should produce a system that is less dependent on kinetic effects. The results of this study are applicable to Fe-isotope fractionation in clastic sediments formed in highly acidic conditions; such conditions may have existed on Mars where acidic oxidizing ground and surface waters may have been responsible for clastic sedimentation and metallic element transport.     

Calculation of equilibrium stable isotope partition function ratios for aqueous zinc complexes and metallic zinc

The goal of this study is to determine reduced partition function ratios for a variety of species of zinc, both as a metal and in aqueous solutions in order to calculate equilibrium stable isotope partitioning. We present calculations of the magnitude of Zn stable-isotope fractionation (^66,67,68Zn/⁶⁴Zn) between aqueous species and metallic zinc using measured vibrational spectra (fit from neutron scattering studies of metallic zinc) and a variety of electronic structure models. The results show that the reduced metal, Zn(0), will be light in equilibrium with oxidized Zn(II) aqueous species, with the best estimates for the Zn(II)-Zn(0) fractionation between hexaquo species and metallic zinc being Δ^66/64Zn_aq-metal ∼ 1.6‰ at 25 °C, and Δ^66/64Zn_aq-metal ∼ 0.8‰ between the tetrachloro zinc complex and metallic zinc at 25 °C using B3LYP/aug-cc-pVDZ level of theory and basis set. To examine the behavior of zinc in various aqueous solution chemistries, models for Zn(II) complex speciation were used to determine which species are thermodynamically favorable and abundant under a variety of different conditions relevant to natural waters, experimental and industrial solutions. The optimal molecular geometries for [Zn(H₂O)₆]²⁺, [Zn(H₂O)₆]·SO₄, [ZnCl₄]²⁻ and [Zn(H₂O)₃(C₃H₅O(COO)₃)]⁻ complexes in various states of solvation, protonation and coordination were calculated at various levels of electronic structure theory and basis set size. Isotopic reduced partition function ratios were calculated from frequency analyses of these optimized structures. Increasing the basis set size typically led to a decrease in the calculated reduced partition function ratios of ∼0.5‰ with values approaching a plateau using the aug-cc-pVDZ basis set or larger. The widest range of species were studied at the B3LYP/LAN2DZ/6-31G^∗ level of theory and basis-set size for comparison. Aqueous zinc complexes where oxygen is bound to the metal center tended to have the largest reduced partition function ratios, with estimated fractionations ranging from 2.2 to 2.9‰ (⁶⁶Zn/⁶⁴Zn) at 25 °C relative to metallic zinc. The tetrahedrally coordinated tetrachloro zinc complex, where zinc is bound exclusively to chloride, had the lowest reduced partition function ratio for a Zn(II) species (Δ^66/64Zn_aq-metal ∼ 1-1.3‰ at 25 °C). Increasing the number of waters in the second shell of solvation of the above complexes led to variable results, most commonly leading to a decrease of ∼0.2 to 0.3‰ in calculated Δ^66/64Zn_aq-metal at 25 °C.These estimates are useful in the interpretation of observed fractionations during the electrochemical deposition of zinc, where aqueous-metal fractionations of up to 5.5‰ are observed. The models show these are not caused by an equilibrium fractionation process. These results suggest that the redox cycle of zinc during industrial processing may be responsible for isotopically distinct reservoirs of zinc observed in polluted environments. The leaching of metallic zinc or zinc tailings from industrial sites could lead to the observed heavy signature in river systems, the magnitude of which will be reliant on the source material and the aqueous species that form.     

Copper isotope fractionation during its interaction with soil and aquatic microorganisms and metal oxy(hydr)oxides: Possible structural control

This work is aimed at quantifying the main environmental factors controlling isotope fractionation of Cu during its adsorption from aqueous solutions onto common organic (bacteria, algae) and inorganic (oxy(hydr)oxide) surfaces. Adsorption of Cu on aerobic rhizospheric (Pseudomonas aureofaciens CNMN PsB-03) and phototrophic aquatic (Rhodobacter sp. f-7bl, Gloeocapsa sp. f-6gl) bacteria, uptake of Cu by marine (Skeletonema costatum) and freshwater (Navicula minima, Achnanthidium minutissimum and Melosira varians) diatoms, and Cu adsorption onto goethite (FeOOH) and gibbsite (AlOOH) were studied using a batch reaction as a function of pH, copper concentration in solution and time of exposure. Stable isotopes of copper in selected filtrates were measured using Neptune multicollector ICP-MS. Irreversible incorporation of Cu in cultured diatom cells at pH 7.5-8.0 did not produce any isotopic shift between the cell and solution (Δ^65/63Cu(solid-solution)) within ±0.2‰. Accordingly, no systematic variation was observed during Cu adsorption on anoxygenic phototrophic bacteria (Rhodobacter sp.), cyanobacteria (Gloeocapsa sp.) or soil aerobic exopolysaccharide (EPS)-producing bacteria (P. aureofaciens) in circumneutral pH (4-6.5) and various exposure times (3 min to 48 h): Δ⁶⁵Cu(solid-solution) = 0.0 ± 0.4‰. In contrast, when Cu was adsorbed at pH 1.8-3.5 on the cell surface of soil the bacterium P. aureofacienshaving abundant or poor EPS depending on medium composition, yielded a significant enrichment of the cell surface in the light isotope (Δ⁶⁵Cu (solid-solution) = −1.2 ± 0.5‰). Inorganic reactions of Cu adsorption at pH 4-6 produced the opposite isotopic offset: enrichment of the oxy(hydr)oxide surface in the heavy isotope with Δ⁶⁵Cu(solid-solution) equals 1.0 ± 0.25‰ and 0.78 ± 0.2‰ for gibbsite and goethite, respectively. The last result corroborates the recent works of Mathur et al. [Mathur R., Ruiz J., Titley S., Liermann L., Buss H. and Brantley S. (2005) Cu isotopic fractionation in the supergene environment with and without bacteria. Geochim. Cosmochim. Acta69, 5233-5246] and Balistrieri et al. [Balistrieri L. S., Borrok D. M., Wanty R. B. and Ridley W. I. (2008) Fractionation of Cu and Zn isotopes during adsorption onto amorhous Fe(III) oxyhydroxide: experimental mixing of acid rock drainage and ambient river water. Geochim. Cosmochim. Acta72, 311-328] who reported heavy Cu isotope enrichment onto amorphous ferric oxyhydroxide and on metal hydroxide precipitates on the external membranes of Fe-oxidizing bacteria, respectively.Although measured isotopic fractionation does not correlate with the relative thermodynamic stability of surface complexes, it can be related to their structures as found with available EXAFS data. Indeed, strong, bidentate, inner-sphere complexes presented by tetrahedrally coordinated Cu on metal oxide surfaces are likely to result in enrichment of the heavy isotope on the surface compared to aqueous solution. The outer-sphere, monodentate complex, which is likely to form between Cu²⁺ and surface phosphoryl groups of bacteria in acidic solutions, has a higher number of neighbors and longer bond distances compared to inner-sphere bidentate complexes with carboxyl groups formed on bacterial and diatom surfaces in circumneutral solutions. As a result, in acidic solution, light isotopes become more enriched on bacterial surfaces (as opposed to the surrounding aqueous medium) than they do in neutral solution.Overall, the results of the present study demonstrate important isotopic fractionation of copper in both organic and inorganic systems and provide a firm basis for using Cu isotopes for tracing metal transport in earth-surface aquatic systems. It follows that both adsorption on oxides in a wide range of pH values and adsorption on bacteria in acidic solutions are capable of producing a significant (up to 2.5-3‰ (±0.1-0.15‰)) isotopic offset. At the same time, Cu interaction with common soil and aquatic bacteria, as well as marine and freshwater diatoms, at 4 < pH < 8 yields an isotopic shift of only ±0.2-0.3‰, which is not related to Cu concentration in solution, surface loading, the duration of the experiment, or the type of aquatic microorganisms.     

Stable Cu and Zn isotope ratios as tracers of sources and transport of Cu and Zn in contaminated soil

Copper and Zn metals are produced in large quantities for different applications. During Cu production, large amounts of Cu and Zn can be released to the environment. Therefore, the surroundings of Cu smelters are frequently metal-polluted. We determined Cu and Zn concentrations and Cu and Zn stable isotope ratios (δ⁶⁵Cu, δ⁶⁶Zn) in three soils at distances of 1.1, 3.8, and 5.3 km from a Slovak Cu smelter and in smelter wastes (slag, sludge, ash) to trace sources and transport of Cu and Zn in soils. Stable isotope ratios were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICP-MS) in total digests. Soils were heavily contaminated with concentrations up to 8087 μg g⁻¹ Cu and 2084 μg g⁻¹ Zn in the organic horizons. The δ⁶⁵Cu values varied little (−0.12‰ to 0.36‰) in soils and most wastes and therefore no source identification was possible. In soils, Cu became isotopically lighter with increasing depth down to 0.4 m, likely because of equilibrium reactions between dissolved and adsorbed Cu species during transport of smelter-derived Cu through the soil. The δ⁶⁶Zn_IRMM values were isotopically lighter in ash (−0.41‰) and organic horizons (−0.85‰ to −0.47‰) than in bedrock (−0.28‰) and slag (0.18‰) likely mainly because of kinetic fractionation during evaporation and thus allowed for separation of smelter-Zn from native Zn in soil. In particular in the organic horizons large variations in δ⁶⁶Zn values occur, probably caused by biogeochemical fractionation in the soil-plant system. In the mineral horizons, Zn isotopes showed only minor shifts to heavier δ⁶⁶Zn values with depth mainly because of the mixing of smelter-derived Zn and native Zn in the soils. In contrast to Cu, Zn isotope fractionation between dissolved and adsorbed species was probably only a minor driver in producing the observed variations in δ⁶⁶Zn values. Our results demonstrate that metal stable isotope ratios may serve as tracer of sources, vertical dislocation, and biogeochemical behavior in contaminated soil.     

Selenium isotopes trace the source and redox processes in the black shale-hosted Se-rich deposits in China

We analyzed the Se isotopic composition of black shales and related kerogen and sulfide fractions from the Zunyi Ni-Mo-Se deposit, the La’erma Se-Au deposit and the Yutangba Se deposit in southern China to constrain metal sources and accumulation processes, both subjects of disagreement in the scientific community. Se at the Zunyi Ni-Mo-Se polymetallic deposit displayed a restricted range of δ⁸²Se values (−1.6‰ to 2.4‰ with a mean of 0.6‰) suggesting a major hydrothermal origin where aqueous Se was probably transported as H₂Se, along with H₂S, and precipitated directly as selenides or in sulfides. Se at the La’erma Se-Au deposit covers a larger range in δ⁸²Se values (−3.8‰ to 5.4‰ with a mean of 0.3‰), suggesting Se redistribution following redox transformations, leading to kinetic isotopic fractionation. The largest Se isotopic variation so far in natural terrestrial samples was found in the Yutangba Se deposit, with δ⁸²Se values varying from −12.77‰ to 4.93‰. On the basis of variations in Se isotopes in the deposit, along with other geological and geochemical evidence, the “redox model” (supergene alteration) explains the occurrence of native Se in the deposit. Overall, hydrothermal systems may be a potentially important Se source to form economic deposits in comparison to seawater sources. Significantly, our study indicates that either secondary hydrothermal or supergene alteration is a key factor in Se enrichment in black shales. Redistribution of Se, and probably other redox-sensitive metals like Mo, Cr and V, leads to isotopic fractionation which may be used to fingerprint such alteration/precipitation processes.     

Boron isotopic fractionation related to boron sorption on humic acid and the structure of surface complexes formed

Boron isotopic fractionation during adsorption onto Ca-flocculated Aldrich humic acid (HA) has been investigated experimentally as a function of solution pH at 25°C and I = 0.15 M. Boron aqueous concentration and isotopic composition were determined by Cs₂BO₂⁺ Positive Thermal Ionization Mass Spectrometry analysis, while the structure of B surface complexes on HA was characterized using ¹¹B Magic Angle Spinning Nuclear Magnetic Resonance (MAS NMR). Significant B sorption on HA was observed at 6 < pH < 12 with a maximum value of Kd, the partition coefficient between adsorbed and aqueous boron, equal to 40 at pH = 9.5-10. Combined ¹¹B MAS NMR analysis and FITEQL modeling of B sorption on HA showed that this element forms tetrahedrally coordinated five- or six-membered ring chelates, most likely 1,2-diol and 1,3-diol complexes at alkaline pH (8 < pH < 11) and dicarboxylic complexes at near neutral conditions (6 < pH < 9). Results of this study demonstrate for the first time that boron sorption on HA induces a strong pH-dependent isotope fractionation—with ¹¹B depleted at the surface of HA—that reaches a maximum at 5 < pH < 9 (α = 0.975, Δ = −25‰) and decreases sharply at pH >9. The measured isotope fractionation cannot be modeled assuming that the isotopic composition of the sorbed borate species is identical to that of B(OH)₄^- species in the parent solution. It is shown that the extent of isotopic fractionation depends not only on B aqueous speciation but also on the distribution and structure of the borate surface complexes formed. In agreement with energetic constrains, calculation of the isotope fractionation factors between aqueous boric acid and boron surface complexes suggests that the formation of the strained six-membered ring 1,3-diol complex yields a much higher fractionation (α_BLP1−III = 0.954-0.960, Δ = −41/-47‰) than that of the very stable five-membered ring 1,2-diol (α_BLP2−III = 0.983, Δ = −18‰). The results of this study open new perspectives to understand and model boron biogeochemical cycle. It is predicted that boron sorption onto organic matter can have important consequences for the boron isotopic composition of surface water reservoirs (seawater, groundwater, soil waters) in which either abundant organic surfaces or significant boron concentrations are available. In addition, the large isotope fractionation between aqueous boric acid and surface boron-organic complexes found in the present work makes boron a promising tracer of biologic activity.     

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

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

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.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.