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.     

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.     

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

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

Effect of temperature and mass transport on transition metal isotope fractionation during electroplating

Transition metal stable isotope signatures can be useful for tracing both natural and anthropogenic signals in the environment, but only if the mechanisms responsible for fractionation are understood. To investigate isotope fractionations due to electrochemistry (or redox processes), we examine the stable isotope behavior of iron and zinc during the reduction reaction  + 2e⁻ = M_metal as a function of electrochemical driving force, temperature, and time. In all cases light isotopes are preferentially electroplated, following a mass-dependent law. Generally, the extent of fractionation is larger for higher temperatures and lower driving forces, and is roughly insensitive to amount of charge delivered. The maximum fractionations are δ^56/54Fe = −4.0‰ and δ^66/64Zn = −5.5‰, larger than observed fractionations in the natural environment and larger than those predicted due to changes in speciation. All the observed fractionation trends are interpreted in terms of three distinct processes that occur during an electrochemical reaction: mass transport to the electrode, chemical speciation changes adjacent to the electrode, and electron transfer at the electrode. We show that a large isotope effect adjacent the electrode surface arises from the charge-transfer kinetics, but this effect is attenuated in cases where diffusion of ions to the electrode surface becomes the rate-limiting step. Thus while a general increase in fractionation is observed with increasing temperature, this appears to be a result of thermally enhanced mass transport to the reacting interface rather than an isotope effect associated with the charge-transfer kinetics. This study demonstrates that laboratory experiments can successfully distinguish isotopic signatures arising from mass transport, chemical speciation, and electron transfer. Understanding how these processes fractionate metal isotopes under laboratory conditions is the first step towards discovering what role these processes play in fractionating metal isotopes in natural systems.     

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.     

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.     

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).     

Zinc and sulfur isotope variation in sphalerite from carbonate-hosted zinc deposits,Cantabria, Spain

We studied zinc and sulfur isotopes and the chemical composition of sphalerite samples from Picos de Europa (Aliva mine) and sphalerite and hydrozincite samples from La Florida mine, two carbonate-hosted Mississippi Valley-type (MVT) deposits located in northern Spain; despite being close, they are hosted in carbonatic rocks of different ages, Lower Carboniferous and Lower Cretaceous, respectively. The two generations of sphalerite at Picos de Europa show different δ⁶⁶Zn values (stage 1 sphalerite +0.24 per mil and stage 2 sphalerite from ?0.75 to +0.08 per mil). Both generations also differ in the sulfur isotope composition (stage 1 has δ³⁴S?=?+6.6 and stage 2 has δ³⁴S?=??0.9 to +2.9 per mil) and the chemical composition (stage 1 sphalerite, compared to stage 2 sphalerite, is significantly enriched in Pb, As, Mn, Sb, slightly enriched in Ag, Ni, and Cu and depleted in Co, Ga, Tl, Te, Ge, and Sn). We suggest that Zn isotope fractionation was controlled predominantly by pH and T changes. High Zn isotope values reflect rapid precipitation of sphalerite from higher-temperature acidic fluids that carried Zn mostly as chloride species after interaction with carbonate rocks while lower Zn isotope values most likely resulted from a longer precipitation process from fluid at higher pH and decreasing T that carried dominantly Zn sulfide species. At La Florida, sphalerite samples show light ⁶⁶Zn-depleted signatures with δ⁶⁶Zn values from ?0.80 to ?0.01 per mil (mostly between ?0.80 and ?0.24 per mil) and δ³⁴S values from +10.7 to +15.7 per mil without any relationship between the δ⁶⁶Zn and δ³⁴S values. Here, the variation in Zn isotope values is interpreted as related to mixing of fluids from two reservoirs. The Zn was carried by a single deep-seated and higher T (~250–320 °C) fluid, and precipitation took place after mixing with a connate S-rich fluid in a system with mH₂S?>?mZn²⁺ as a result of change in pH, T, and Zn predominant species. The light δ⁶⁶Zn accompanied by heavy δ³⁴S values resulted from fractionation of Zn aqueous sulfides at near-neutral pH and decreasing T. Hydrozincite samples show much heavier δ⁶⁶Zn values (+0.21 to +0.33 per mil), consistent with fractionation during supergene processes.     

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.     

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.     

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

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

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

Speleothem calcite farmed in situ: Modern calibration of δO and δC paleoclimate proxies in a continuously-monitored natural cave system

Understanding the relationships between speleothem stable isotopes (δ¹³C δ¹⁸O) and in situ cave forcing mechanisms is important to interpreting ancient stalagmite paleoclimate records. Cave studies have demonstrated that the δ¹⁸O of inorganically precipitated (low temperature) speleothem calcite is systematically heavier than the δ¹⁸O of laboratory-grown calcite for a given temperature. To understand this apparent offset, rainwater, cave drip water, groundwater, and modern naturally precipitated calcite (farmed in situ) were grown at multiple locations inside Hollow Ridge Cave in Marianna, Florida. High resolution micrometeorological, air chemistry time series and ventilation regimes were also monitored continuously at two locations inside the cave, supplemented with periodic bi-monthly air gas grab sample transects throughout the cave.Cave air chemistry and isotope monitoring reveal density-driven airflow pathways through Hollow Ridge Cave at velocities of up to 1.2 m s⁻¹ in winter and 0.4 m s⁻¹ in summer. Hollow Ridge Cave displays a strong ventilation gradient in the front of the cave near the entrances, resulting in cave air that is a mixture of soil gas and atmospheric CO₂. A clear relationship is found between calcite δ¹³C and cave air ventilation rates estimated by proxies pCO₂ and ²²²Rn. Calcite δ¹³C decreased linearly with distance from the front entrance to the interior of the cave during all seasons, with a maximum entrance-to-interior gradient of Δδ¹³C_CaCO3 = −7‰. A whole-cave “Hendy test” at multiple contemporaneous farming sites reveals that ventilation induces a +1.9 ± 0.96‰ δ¹³C offset between calcite precipitated in a ventilation flow path and calcite precipitated on the edge or out of flow paths. This interpretation of the “Hendy test” has implications for interpreting δ¹³C records in ancient speleothems. Calcite δ¹³C_CaCO3 may be a proxy not only for atmospheric CO₂ or overlying vegetation shifts but also for changes in cave ventilation due to dissolution fissures and ceiling collapse creating and plugging ventilation windows.Farmed calcite δ¹⁸O was found to exhibit a +0.82 ± 0.24‰ offset from values predicted by both theoretical calculations and laboratory-grown inorganic calcite. Unlike δ¹³C_CaCO3, oxygen isotopes showed no ventilation effects, i.e. Δδ¹⁸O_CaCO3 appears to be a function of growth temperature only although we cannot rule out a small effect of (unmeasured) gradients in relative humidity (evaporation) accompanying ventilation. Our results support the findings of other cave investigators that water-calcite fractionation factors observed in speleothem calcite are higher that those measured in laboratory experiments. Cave and laboratory calcite precipitates may differ mainly in the complex effects of kinetic isotope fractionation. Combining our data with other recent speleothem studies, we find a new empirical relationship for cave-specific water-calcite oxygen isotope fractionation across a range of temperatures and cave environments:

1000lnα=16.1(10³T^-1)-24.6     

The isotopic evolution of atmospheric Pb in central Ontario since AD 1800, and its impacts on the soils, waters, and sediments of a forested watershed, Kawagama Lake

A peat core from an ombrotrophic bog documents the isotopic evolution of atmospheric Pb in central Ontario since AD 1804 ± 53 (²¹⁰Pb dating). Despite the introduction of unleaded gasoline in the mid-1970’s, the ratio ²⁰⁶Pb/²⁰⁷Pb in atmospheric deposition has not increased as expected, but rather continues to decline. In fact, snowpack sampling (2005 and 2009) and rainwater samples (2008) show that the isotopic composition of atmospheric Pb today is often far less radiogenic than the gasoline lead that had been used in Canada in the past. The peat, snow, and rainwater data presented here are consistent with the Pb isotope data for aerosols collected in Dorset in 1984 and 1986 which were traced by Sturges and Barrie (1989) to emissions from the Noranda smelter in northern Quèbec, Canada’s largest single source of atmospheric Pb. Understanding atmospheric Pb deposition in central Ontario, therefore, requires not only consideration of natural sources and past contributions from leaded gasoline, but also emissions from metal smelting and refining.Lead in the streams which enter Kawagama Lake today (²⁰⁶Pb/²⁰⁷Pb = 1.16 − 1.19) represents a mixture between the natural values (1.191 − 1.201 estimated using pre-industrial lake sediments) and the values found in the humus layer of the surrounding forest soils (²⁰⁶Pb/²⁰⁷Pb = 1.15 − 1.19). In the lake itself, however, Pb is much less radiogenic (²⁰⁶Pb/²⁰⁷Pb as low as 1.09) than in the streams, with the dissolved fraction less radiogenic than particulate material. The evolution of Pb isotope ratios within the watershed apparently reflects preferential removal by sedimentation of comparatively dense, radiogenic, terrestrial particles (derived from the mineral fraction of soils) from the humus particles with lower ratios of ²⁰⁶Pb/²⁰⁷Pb (because of atmospheric Pb contamination). Despite the contemporary enrichments of Pb in rain and snow, concentrations of dissolved Pb in the lake are extremely low (sometimes below 10 ng/l), with Pb concentrations and Pb/Sc ratios approaching “natural” values because of efficient binding to particles, and their subsequent removal in the watershed.     

Stable carbon and oxygen isotope fractionation in bivalve (Placopecten magellanicus) larval aragonite

The relationship between stable isotope composition (δ¹³C and δ¹⁸O) in seawater and in larval shell aragonite of the sea scallop, Placopecten magellanicus, was investigated in a controlled experiment to determine whether isotopes in larval shell aragonite can be used as a reliable proxy for environmental conditions. The linear relationship between δ¹³C_DIC and δ¹³C_aragonite (r² = 0.97, p < 0.0001, RMSE = 0.18) was:

δ¹³C_DIC=1.15(±0.05)∗δ¹³C_aragonite-0.85(±0.04)     

Use of isotope ratios to assess sources of Pb and Zn dispersed in the environment during mining and ore processing within the Orlovka–Spokoinoe mining site (Russia)

Element concentrations, element ratios and Pb and Zn isotope data are reported for different geologic samples (barren and ore-bearing granites and host rocks), technogenic products (ore concentrates and tailings) and biologic samples (lichens and birch leaves) from the Orlovka–Spokoinoe mining district, Eastern Transbaikalia, Russia, with the aim to trace the sources of Pb and Zn at a local level within the mining site. Lichens and birch leaves were used as receptors of contamination within the mining site. Pb/Zr and Zn/Zr values indicated Pb and Zn enrichment relative to host rocks. Zn isotope data of 15 geologic and 11 lichen samples showed different Zn isotopic signatures with the total range for the geologic suite of −0.4‰ to +1.2‰ and for lichens of +0.4‰ to +1.4‰ in δ⁶⁶Zn relative to Lyon JMC Zn standard. The source of isotopically heavy Zn within the Orlovka–Spokoinoe mining site could be potentially associated with long-range atmospheric aerosols that also contributed Pb to the studied mining site. Our results demonstrated that Zn isotopes might be used as new tools for Zn source assessment.     

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.     

Zinc Isotopic Compositions of NIST SRM 683 and Whole‐Rock Reference Materials

In this study the homogeneity of the zinc isotopic composition in the NIST SRM 683 reference material was examined by measuring the Zn isotopic signature in microdrilled sample powders from two metal nuggets. Zinc was purified using AG MP‐1M resin and then measured by MC‐ICP‐MS. Instrumental mass bias was corrected using the “sample‐standard bracketing” method and empirical external normalisation with Cu doping. After evaluating the potential effects of varying acid mass fractions and different matrices, high‐precision Zn isotope data were obtained with an intermediate measurement precision better than ± 0.05‰ (δ⁶⁶Zn, 2s) over a period of 5 months. The δ⁶⁶Zn_JMC‐Lyon mean values of eighty‐four and fourteen drilled powders from two nuggets were 0.11 ± 0.02‰ and 0.12 ± 0.02‰, respectively, indicating that NIST SRM 683 is a good isotopic reference material with homogeneous Zn isotopes. The Zn isotopic compositions of seventeen rock reference materials were also determined, and their δ⁶⁶Zn values were in agreement with most previously published data within 2s. The δ⁶⁶Zn values of most of the rock reference materials analysed were in the range 0.22–0.36‰, except for GSP‐2 (1.07 ± 0.06‰, n = 12), NOD‐A‐1 (0.96 ± 0.03‰, n = 6) and NOD‐P‐1 (0.78 ± 0.03‰, n = 6). These comprehensive data should serve as reference values for quality assurance and interlaboratory calibration exercises.     

Rare earth element and Pb isotope variations in a 52 kyr peat core from Lynch’s Crater (NE Queensland, Australia): Proxy development and application to paleoclimate in the Southern Hemisphere

Accurate prediction of future climate scenarios is contingent on our understanding of past and present climate mechanisms. This is done in part through the reconstruction of historical climate changes using proxy records from terrestrial and marine archives. Terrestrial archives covering the Holocene and late Pleistocene are limited, most acutely in the Southern Hemisphere. Here, Rare earth elements (REE) and Pb isotopes are developed as inorganic geochemical proxies of mineral dust source changes and, by extension, climate change. Using a peat core from Lynch’s Crater in NE Queensland, Australia, we present the first long-term (c. 52 kyr) terrestrial record of atmospheric REE and Pb deposition (with the exception of four wet events which represent periods of erosion from the crater itself) in the Southern Hemisphere covering both glacial and interglacial times. Based on a combination of correlation analyses, Al and Ti normalised profiles and elemental patterns, we establish REE are immobile within the peat deposit and not subject to significant post depositional diagenetic changes (important particularly for Ce). This is vital as REE can be mobile under acid and organic rich conditions like those that can occur during the development of a peat deposit. The volcanic provinces of eastern Australia have characteristic Eu anomaly signatures, which allowed their use in a novel way to detect changes in dust source to Lynch’s Crater. Between 41,095 and 52,505 BP the deposit was under the influence of dust carried by long distance transport (>1500 km) from SE Australia. From 8525 to 40,815 BP regional sources (100-1500 km) dominated the deposited signals while between 1740 and 8390 BP the dust signal was controlled by local sources (<100 km). These findings were also confirmed by Pb isotope data. Changepoint modelling refined the timing of these changes in dust source, recognizing concurrent shifts in our tracing tools ((Eu/Eu^∗)_PAAS and ²⁰⁶Pb/²⁰⁷Pb). These changepoints were then compared to other palaeoenvironmental records (pollen, lake levels and dune building) from eastern Australia and found to be similar. Our results demonstrate that REE and Pb isotopes are effective tools for tracing past changes in atmospheric dust sources and to the study of climate change using minerotrophic peat deposits.     

Sources and biological fractionation of Silicon isotopes in the Eastern Equatorial Pacific

Silicon isotopes in dissolved silicic acid were measured in the upper four kilometers between 4°N and 3°S latitude at 110°W longitude in the eastern Equatorial Pacific. Silicon isotopes became progressively heavier with silicic acid depletion of surface water as expected from biological fractionation. The value of ε estimated by applying a steady-state isotope fractionation model to data from all stations between 4°N and 3°S was −0.77 ± 0.12‰ (std. err.). When the analysis was restricted to those stations whose temperature and salinity profiles indicated that they were directly influenced by upwelling of the Equatorial Undercurrent (EUC), the resulting value of ε was −1.08 ± 0.27‰ (std. err.) similar to the value established in culture studies (−1.1‰). When the non steady state Rayleigh model was applied to the same restricted data set the resulting value of ε was significantly more positive, −0.61 ± 0.16‰ (std. err.). To the extent that the equatorial system approximates a steady state these results support a value of −1.1‰ for the fractionation factor for isotopes of Si in the sea. Without the assumption of steady state the value of ε can only be constrained to be between −0.6 and −1.1‰. Silicic acid in Equatorial Pacific Deep Water below 2000 m had a near constant δ³⁰Si of +1.32 ± 0.05‰. That value is significantly more positive than obtained for North Pacific Deep Water at similar depths at stations to the northwest of our study area (0.9-1.0‰) and it is slightly less positive than new measures of the δ³⁰Si of silicic acid from the silicic acid plume centered over the Cascadia basin in the Northeast Pacific (Si(OH)₄ > 180  μM, δ³⁰Si = +1.46 ± 0.12‰ (SD, n = 4). We show that the data from the equator and Cascadia basin fit a general trend of increasing δ³⁰Si(OH)₄ with increasing silicic acid concentration in the deep sea, but that the isotope values from the Northeast Pacific are anomalously light. The observed level of variation in the silicon isotope composition of deep waters from this single ocean basin is considerably larger than that predicted by current models based on fractionation during opal formation with no isotope effect during dissolution. Confirmation of such high variability in deep water δ³⁰Si(OH)₄ within individual ocean basins will require reassessment of the mechanisms controlling the distribution of isotopes of silicon in the sea.