Episodic diamond genesis at Jwaneng, Botswana, and implications for Kaapvaal craton evolution

Major element and Re–Os isotope analysis of single sulfide inclusions in diamonds from the 240 Ma Jwaneng kimberlite has revealed the presence of at least two generations of eclogitic diamonds at this locality, one Proterozoic (ca. 1.5 Ga) and the other late Archean (ca. 2.9 Ga). The former generation is considered to be the same as that of eclogitic garnet and clinopyroxene inclusion bearing diamonds from Jwaneng with a Sm–Nd isochron age of 1.54 Ga. The latter is coeval with the 2.89 Ga subduction-related generation of eclogitic sulfide inclusion bearing diamonds from Kimberley formed during amalgamation of the western and eastern Kaapvaal craton near the Colesberg magnetic lineament.

The Kimberley, Jwaneng, and Premier kimberlites are key localities for characterizing the relationship between episodic diamond genesis and Kaapvaal craton evolution. Kimberley has 3.2 Ga harzburgitic diamonds associated with creation of the western Kaapvaal cratonic nucleus, and 2.9 Ga eclogitic diamonds resulting from its accretion to the eastern Kaapvaal. Jwaneng has two main eclogitic diamond generations (2.9 and 1.5 Ga) reflecting both stabilization and subsequent modification of the craton. Premier has 1.9 Ga lherzolitic diamonds that postdate Bushveld–Molopo magmatism (but whose precursors have Archean Sm–Nd model ages), as well as 1.2 Ga eclogitic diamonds. Thus, Jwaneng provides the overlap between the dominantly Archean vs. Proterozoic diamond formation evident in the Kimberley and Premier diamond suites, respectively. In addition, the 1.5 Ga Jwaneng eclogitic diamond generation is represented by both sulfide and silicate inclusions, allowing for characterization of secular trends in diamond type and composition. Results for Jwaneng and Kimberley eclogitic sulfides indicate that Ni- and Os-rich end members are more common in Archean diamonds compared to Proterozoic diamonds. Similarly, published data for Kimberley and Premier peridotitic silicates show that Ca-rich (lherzolitic) end members are more likely to be found in Proterozoic diamonds than Archean diamonds. Thus, the available diamond distribution, composition, and age data support a multistage process to create, stabilize, and modify Archean craton keels on a billion-year time scale and global basis.     

Diamonds from Jagersfontein (South Africa): messengers from the sublithospheric mantle

The diamond population from the Jagersfontein kimberlite is characterized by a high abundance of eclogitic, besides peridotitic and a small group of websteritic diamonds. The majority of inclusions indicate that the diamonds are formed in the subcratonic lithospheric mantle. Inclusions of the eclogitic paragenesis, which generally have a wide compositional range, include two groups of eclogitic garnets (high and low Ca) which are also distinct in their rare earth element composition. Within the eclogitic and websteritic suite, diamonds with inclusions of majoritic garnets were found, which provide evidence for their formation within the asthenosphere and transition zone. Unlike the lithospheric garnets all majoritic garnet inclusions show negative Eu-anomalies. A narrow range of isotopically light carbon compositions (δ¹³C −17 to −24 ‰) of the host diamonds suggests that diamond formation in the sublithospheric mantle is principally different to that in the lithosphere. Direct conversion from graphite in a subducting slab appears to be the main mechanism responsible for diamond formation in this part of the Earth’s mantle beneath the Kaapvaal Craton. The peridotitic inclusion suite at Jagersfontein is similar to other diamond deposits on the Kaapvaal Craton and characterized by harzburgitic to low-Ca harzburgitic compositions.     

The evolution of refertilized lithospheric mantle beneath the northeastern Siberian craton: Links between mantle metasomatism,thermal state and diamond potential

Although the diamond potential of cratons is linked mainly to thick and depleted Archean lithospheric keels, there are examples of craton-edge locations and circum-cratonic Proterozoic terranes underlain by diamondiferous mantle. Here, we use the results of comprehensive major and trace-element studies of detrital garnets from diamond-rich Late Triassic (Carnian) sedimentary rocks in the northeastern Siberia to constrain the thermal and chemical state of the pre-Triassic mantle and its ability to sustain the diamond storage. The studied detrital mantle-derived garnets are dominated by low- to medium-Cr lherzolitic (～45%) and low-Cr megacrystic (～39%) chemistries, with a significant proportion of eclogitic garnets (～11%), and only subordinate contribution from harzburgitic garnets (～5%) with variable Cr₂O₃ contents (1.2–8.4 wt.%). Low-Cr megacrysts display uniform, “normal” rare-earth element (REE) patterns with no Eu/Eu* anomalies, systematic Zr and Ti enrichment (mainly within 2.5–5), which are evidence of their crystallization from deep metasomatic melts. Lherzolitic (G9) garnets exhibit normal or humped to MREE-depleted sinusoidal REE patterns and elevated Nd/Y (up to 0.33–0.41) and Zr/Y ratios (up to 7.62). Rare low- to high-Cr harzburgitic (G10) garnets have primarily “depleted”, sinusoidal REE-patterns, low Ti, Y and HREE, but vary significantly in Zr-Hf, Ti and MREE-HREE contents, Nd/Y (within 0.1–2.4) and Zr/Y (1.53–19.9) ratios. The observed trends of chemical enrichment from the most depleted, harzburgitic garnets towards lherzolitic (including high-Ti high-Cr G11-type) garnets and megacrysts result from either voluminous high-temperature metasomatism by plume-derived silicate melts or recurrent mobilization of less voluminous kimberlitic or related carbonated mantle melts, rather than the initially primitive, fertile nature of the Proterozoic SCLM. Calculated Ni-in-garnet temperatures (primarily within ～1150–1250 °C) indicate their derivation from at least ～220 km thick Cr-undersaturated lithosphere at the relevant Devonian to Triassic thermal flow of ～45 mW/m² or cooler. We suggest the existence of rare harzburgitic domains in the primarily lherzolitic diamond-facies SCLM beneath the northeastern Siberian craton at least by Triassic, whereas the abundance of eclogitic garnets, predominance of E-type inclusions in placer diamonds and specific morphologies argue for diamondiferous eclogites occurring within a ～50–65 kbar diamond window of the Olenek province by the same time.     

Rare and unusual mineral inclusions in diamonds from Mwadui, Tanzania

Syngenetic diamond inclusions from the Mwadui kimberlite reveal that an unusually fertile section of lithospheric mantle beneath the Central African Craton was sampled. This is shown by a very high ratio of lherzolitic to harzburgitic garnet inclusions (1:2) and low Mg/Fe-ratios in olivine and orthopyroxene. Geothermometry applied to the peridotitic inclusions indicates disequilibrium between non-touching inclusion pairs to be common. Disequilibrium between garnet-olivine and garnet-orthopyroxene pairs suggests successive iron enrichment during diamond formation, e.g. leading to the presence of harzburgitic garnet and lherzolitic olivine in the same diamond. Apart from the dominant peridotitic inclusion suite (88%), rare eclogitic inclusions occur (2%) and a number of uncertain paragenesis. Two diamonds, one with eclogitic garnets with moderate pyroxene solid solution and the other with a single ferro-periclase inclusion, suggest the contribution of a small sub-lithospheric component. The finding of the association Fe-FeO-Fe₃O₄ in one single diamond indicates diamond formation over a large range of f _O2 conditions, possibly along redox fronts. Steep compositional gradients may also be reflected by the joint occurrence of harzburgitic garnet and a SiO₂-phase in the same diamond. Alternatively the formation of the SiO₂-phase may be due to extreme carbonation of the peridotitic source. Further unusual findings include the exsolution of a silicate phase from magnetite inclusions, (i.e. primary solution of γ-olivine) and an ilmenite inclusion with an eskolaite (Cr₂O₃) component of 14.5 mol%, the latter together with harzburgitic paragenesis silicate inclusions. Received: 23 August 1997 / Accepted: 7 January 1998     

Nitrogen isotope fractionations during progressive metamorphism: A case study from the Paleozoic Cooma metasedimentary complex, southeastern Australia

The well-studied Paleozoic Cooma metamorphic complex in southeastern Australia is characterized by a uniform siliciclastic protolith, of uniform age, with a continuous range of metamorphic grade from subgreenschist- to upper amphibolite-facies, and migmatite-grade in an annular pattern around the Cooma granodiorite. Those conditions are optimal for investigating variations of N concentrations and δ¹⁵N values during progressive metamorphism. Nitrogen concentrations decrease and δ¹⁵N increases with increasing metamorphic grade (sub-chlorite zone: 120 ppm N, δ¹⁵N = 2.3‰; chlorite zone: 110 ppm N, δ¹⁵N = 3.0‰; biotite and andalusite zone: 85 ppm N, δ¹⁵N = 3.8 ‰; sillimanite and migmatite zones: 40 ppm N, δ¹⁵N = 10.7‰). Covariation of K and N contents is consistent with N substituting for K as NH₄⁺ in micas. Observed trends of increasing δ¹⁵N values with decreasing nitrogen concentrations can be explained by a continuous release of nitrogen depleted in ¹⁵N with progressive metamorphism, which causes an enrichment of ¹⁵N in the residual nitrogen of the rock. Equilibrium models for Rayleigh distillation and batch volatilisation for data of the greenschist and amphibolite facies metasedimentary rocks can be explained by N₂-NH₄⁺ exchange at temperatures of 300-600 °C, whereas observed large fractionations for the upper amphibolite-facies and melt products in the migmatite-grade samples may be interpreted as NH₃-NH₄⁺ exchanges at temperature of 650-730 °C. Lower values in the highest grade zones may also stem in part from input of ¹⁵N-depleted fluids from the granodiorite.The magnitude of isotope fractionation of nitrogen is about 1-2‰ during progressive metamorphism of metasedimentary rocks from sub-chlorite zone to biotite-andalusite zone, which is consistent with previous studies. Consequently, the large spread of δ¹⁵N values in Archean greenschist-facies metasedimentary rocks of −6‰ to 30‰ can be accounted for by variable mixtures of mantle plume-dominated volatiles with a δ¹⁵N of −5‰, and a ¹⁵N-enriched marine sedimentary kerogen component inherited from a CI chondrite veneer having δ¹⁵N of 30‰ to 42‰.     

Regional patterns in the paragenesis and age of inclusions in diamond, diamond composition, and the lithospheric seismic structure of Southern Africa

The Archean lithospheric mantle beneath the Kaapvaal–Zimbabwe craton of Southern Africa shows ±1% variations in seismic P-wave velocity at depths within the diamond stability field (150–250 km) that correlate regionally with differences in the composition of diamonds and their syngenetic inclusions. Seismically slower mantle trends from the mantle below Swaziland to that below southeastern Botswana, roughly following the surface outcrop pattern of the Bushveld-Molopo Farms Complex. Seismically slower mantle also is evident under the southwestern side of the Zimbabwe craton below crust metamorphosed around 2 Ga. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa, and Jwaneng have Re–Os isotopic ages that range from circa 2.9 Ga to the Proterozoic and show little correspondence with these lithospheric variations. However, silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane do show some regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity correlates with a greater proportion of eclogitic versus peridotitic silicate inclusions in diamond, a greater incidence of younger Sm–Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds whereas the converse is true for diamonds from higher velocity mantle. The oldest formation ages of diamonds indicate that the mantle keels which became continental nuclei were created by middle Archean (3.2–3.3 Ga) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of sulfide inclusions that are eclogitic in the 2.9 Ga age population links late Archean (2.9 Ga) subduction-accretion events involving an oceanic lithosphere component to craton stabilization. These events resulted in a widely distributed younger Archean generation of eclogitic diamonds in the lithospheric mantle. Subsequent Proterozoic tectonic and magmatic events altered the composition of the continental lithosphere and added new lherzolitic and eclogitic diamonds to the already extensive Archean diamond suite.     

Cl,N, C isotopic analysis of common agro-chemicals for identifying non-point source agricultural contaminants

The isotopic compositions of commercially available herbicides were analyzed to determine their respective ¹⁵N, ¹³C and ³⁷Cl signatures for the purposes of developing a discrete tool for tracing and identifying non-point source contaminants in agricultural watersheds. Findings demonstrate that of the agrochemicals evaluated, chlorine stable isotopes signatures range between δ³⁷Cl = −4.55‰ and +3.40‰, whereas most naturally occurring chlorine stable isotopes signatures, including those of road salt, sewage sludge and fertilizers, vary in a narrow range about the Standard Mean Ocean Chloride (SMOC) between −2.00‰ and +1.00‰. Nitrogen stable isotope values varied widely from δ¹⁵N = −10.86‰ to +1.44‰ and carbon stable isotope analysis gave an observed range between δ¹³C = −37.13‰ and −21.35‰ for the entire suite of agro-chemicals analyzed. When nitrogen, carbon and chlorine stable isotope analyses were compared in a cross-correlation analysis, statistically independent isotopic signatures exist suggesting a new potential tracer tool for identifying herbicides in the environment.     

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

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

The oxygen-isotope composition of chondrules and isolated forsterite and olivine grains from the Tagish Lake carbonaceous chondrite

The oxygen-isotope compositions (obtained by laser fluorination) of hand-picked separates of isolated forsterite, isolated olivine and chondrules from the Tagish Lake carbonaceous chondrite describe a line (δ¹⁷O = 0.95 * δ¹⁸O − 3.24; R² = 0.99) similar to the trend known for chondrules from other carbonaceous chondrites. The isolated forsterite grains (Fo_99.6-99.8; δ¹⁸O = −7.2‰ to −5.5‰; δ¹⁷O = −9.6‰ to −8.2‰) are more ¹⁶O-rich than the isolated olivine grains (Fo_39.6-86.8; δ¹⁸O = 3.1‰ to 5.1‰; δ¹⁷O = −0.3‰ to 2.2‰), and have chemical and isotopic characteristics typical of refractory forsterite. Chondrules contain olivine (Fo_97.2-99.8) with oxygen-isotope compositions (δ¹⁸O = −5.2‰ to 5.9‰; δ¹⁷O = −8.1‰ to 1.2‰) that overlap those of isolated forsterite and isolated olivine. An inverse relationship exists between the Δ¹⁷O values and Fo contents of Tagish Lake isolated forsterite and chondrules; the chondrules likely underwent greater exchange with ¹⁶O-poor nebular gases than the forsterite. The oxygen-isotope compositions of the isolated olivine grains describe a trend with a steeper slope (1.1 ± 0.1, R² = 0.94) than the carbonaceous chondrite anhydrous mineral line (CCAM_slope = 0.95). The isolated olivine may have crystallized from an evolving melt that exchanged with ¹⁶O-poor gases of somewhat different composition than those which affected the chondrules and isolated forsterite. The primordial components of the Tagish Lake meteorite formed under conditions similar to other carbonaceous chondrite meteorite groups, especially CMs. Its alteration history has its closest affinities to CI carbonaceous chondrites.     

Diamondiferous lithospheric roots along the western margin of the Kalahari Craton—the peridotitic inclusion suite in diamonds from Orapa and Jwaneng

The Orapa and Jwaneng kimberlites are located along the western margin of the Kalahari Craton and the prevalence of eclogitic over peridotitic diamonds in both mines has recently been linked to lower P-wave velocities in the deep mantle lithosphere (relative to the bulk of the craton) to suggest a diamond formation event prompted by mid-Proterozoic growth and modification of preexisting Archean lithosphere (Shirey et al. 2002). Here we study peridotitic diamonds from both mines, with an emphasis on the style of metasomatic source enrichment, to evaluate their relationship with this major eclogitic diamond formation event. In their major element chemistry, the peridotitic inclusions compare well with a world-wide database but reveal differences to diamond sources located in the interior of the Western Terrane of the Kaapvaal block, where the classical mines in the Kimberley region are located. The most striking difference is the relative paucity of low-Ca (<2 wt% CaO in garnet) harzburgites and a low ratio of harzburgitic to lherzolitic garnets (2:1). This suggests that lithospheric mantle accreted to the rim of the Zimbabwe and Kaapvaal blocks was overall chemically less depleted. Alternatively, this more fertile signature may be assigned to stronger metasomatic re-enrichment but the trace element signature of garnet inclusions is not in favor of strong enrichment in major elements. For both mines the majority of lherzolitic and harzburgitic garnet inclusions are characterized by moderately sinusoidal REE_N patterns and low Ti, Zr and Y contents, indicative of a metasomatic agent with very high LREE/HREE and low HFSE. This is consistent with metasomatism by a CHO-fluid or, as modeled by Burgess and Harte (2003), a highly fractionated, low-volume silicate melt from the MORB-source. In both cases, changes in the major element chemistry of the affected rocks will be limited. In a few garnets from Orapa preferential MREE enrichment is observed, suggesting that the percolating fluid/melt fractionated a LREE-phyllic phase (such as crichtonite). The overall moderate degree of metasomatism reflected by the inclusion chemistry is in stark contrast to lithospheric sections for Orapa and Jwaneng based on mantle xenocrysts and xenoliths, revealing extensive mantle metasomatism (Griffin et al. 2003). This suggests that the formation of peridotitic diamonds predates the intensive modification of the subcratonic lithosphere during Proterozoic rifting and compression, implying that diamonds may survive major tectonothermal events.Editorial responsibility: J. Hoefs     

Sulfur and carbon isotope records from 1700 to 800 Ma carbonates of the Jixian section, northern China: Implications for secular isotope variations in Proterozoic seawater and relationships to global supercontinental events

This study is a comprehensive, stable isotope survey of the marine carbonate-dominated, upper Paleo- to lower Neoproterozoic stratigraphy of Jixian County, China. Carbonate-associated sulfate (CAS) was extracted and measured for δ³⁴S_CAS using the same samples analyzed for δ¹³C_carbonate. This integrated proxy approach is a step towards a more comprehensive picture of secular variation in the composition of Proterozoic seawater. We specifically sampled marine carbonate intervals from the lower section of the Chuanlinggou Formation, Changcheng Group (ca. 1700 Ma) to the top of the Jingeryu Formation, Qingbaikou Group (ca. 800 Ma). δ¹³C_carbonate values are mostly negative in the upper Paleoproterozoic Changcheng Group, with an ascending trend from −3‰ to 0‰. We observed variation of approximately 0 ± 1‰ in the Mesoproterozoic Jixian Group, and positive values of +2 ± 2‰ characterize the lower Neoproterozoic Qingbaikou Group. Stratigraphic variations in δ³⁴S_CAS are more remarkable in their ranges and magnitudes, including conspicuously high values exceeding +30‰ in the three intervals at ca. 1700 Ma, 1300-1100 Ma, and 1000-900 Ma. In the Changcheng Group, δ³⁴S_CAS values are typically higher than +25‰, with only a few values of less than +15‰. In contrast, most of the data spanning from the Mesoproterozoic Tieling Formation of the Jixian Group to the lower Neoproterozoic Jingeryu Formation of the Qingbaikou Group are highly variable between +10‰ and +25‰, with some values exceeding +25‰.In the late Paleoproterozoic (1700-1600 Ma), a >10‰ decrease in δ³⁴S_CAS and ∼3‰ increase in δ¹³C_carbonate are coincident with, and likely related to, the breakup of Columbia, a supercontinent that predated Rodinia. Carbon and sulfur isotope data from the Mesoproterozoic, when global tectonic activity was comparatively weaker, fall mostly in the ranges of +15 ± 10‰ and 0 ± 1‰, respectively, but fluctuations of >20‰ for δ³⁴S_CAS and >3‰ for the δ¹³C_carbonate at ca. 1450-1400 Ma may reflect subduction and large-scale magmatic activity in island arcs marking the end of Columbia breakup. From the late Mesoproterozoic (ca. 1300-1100 Ma) to the early Neoproterozoic (ca. 800 Ma), the δ¹³C and δ³⁴S of seawater increased gradually with increasing variability. Most impressive areδ³⁴S_CAS values that exceed +30‰ in two intervals at ca. 1300-1100 Ma and ca. 1000-900 Ma, which may reflect the assembly and early breakup of Rodinia. Although gaps in the record remain, and studies of even higher resolution are warranted, our results suggest that changes in paleoceanographic conditions linked to global tectonics strongly influenced the biogeochemical cycles of C and S. Furthermore, periods of the Proterozoic previously noted for their isotopic invariability show clear isotopic expressions of this tectonic activity.     

Quadruple sulfur isotope analysis of ca. 3.5 Ga Dresser Formation: New evidence for microbial sulfate reduction in the early Archean

Multiple sulfur isotope system is a powerful new tracer for atmospheric, volcanic, and biological influences on sulfur cycles in the anoxic early Earth. Here, we report high-precision quadruple sulfur isotope analyses (³²S/³³S/³⁴S/³⁶S) of barite, pyrite in barite, and sulfides in related hydrothermal and igneous rocks occurring in the ca. 3.5 Ga Dresser Formation, Western Australia. Our results indicate that observed isotopic variations are mainly controlled by mixing of mass-dependently (MD) and non-mass-dependently fractionated (non-MD) sulfur reservoirs. Based on the quadruple sulfur isotope systematics (δ³⁴S-Δ³³S-Δ³⁶S) for these minerals, four end-member sulfur reservoirs have been recognized: (1) non-MD sulfate (δ³⁴S = −5 ± 2‰; Δ³³S = −3 ± 1‰); (2) MD sulfate (δ³⁴S = +10 ± 3‰); (3) non-MD sulfur (δ³⁴S > +6‰; Δ³³S > +4‰); and (4) igneous MD sulfur (δ³⁴S = Δ³³S = 0‰). The first and third components show a clear non-MD signatures, thus probably represent sulfate and sulfur aerosol inputs. The MD sulfate component (2) is enriched in ³⁴S (+10 ± 3‰) and may have originated from microbial and/or abiotic disproportionation of volcanic S or SO₂. Our results reconfirm that the Dresser barites contain small amounts of pyrite depleted in ³⁴S by 15-22‰ relative to the host barite. These barite-pyrite pairs exhibit a mass-dependent relationship of δ³³S/δ³⁴S with slope less than 0.512, which is consistent with that expected for microbial sulfate reduction and is significantly different from that of equilibrium fractionation (0.515). The barite-pyrite pairs also show up to 1‰ difference in Δ³⁶S values and steep Δ³⁶S/Δ³³S slopes, which deviate from the main Archean array (Δ³⁶S/Δ³³S = −0.9) and are comparable to isotope effects exhibited by sulfate reducing microbes (Δ³⁶S/Δ³³S = −5 to −11). These new lines of evidence support the existence of sulfate reducers at ca. 3.5 Ga, whereas microbial sulfur disproportionation may have been more limited than recently suggested.     

Stable isotope and petrologic evidence for open-system degassing during the climactic and pre-climactic eruptions of Mt. Mazama, Crater Lake, Oregon

Evaluation of the extent of volatile element recycling in convergent margin volcanism requires delineating likely source(s) of magmatic volatiles through stable isotopic characterization of sulfur, hydrogen and oxygen in erupted tephra with appropriate assessment of modification by degassing. The climactic eruption of Mt. Mazama ejected approximately 50 km³ of rhyodacitic magma into the atmosphere and resulted in formation of a 10-km diameter caldera now occupied by Crater Lake, Oregon (lat. 43°N, long. 122°W). Isotopic compositions of whole-rocks, matrix glasses and minerals from Mt. Mazama climactic, pre-climactic and postcaldera tephra were determined to identify the likely source(s) of H₂O and S. Integration of stable isotopic data with petrologic data from melt inclusions has allowed for estimation of pre-eruptive dissolved volatile concentrations and placed constraints on the extent, conditions and style of degassing.Sulfur isotope analyses of climactic rhyodacitic whole rocks yield δ³⁴S values of 2.8-14.8‰ with corresponding matrix glass values of 2.4-13.2‰. δ³⁴S tends to increase with stratigraphic height through climactic eruptive units, consistent with open-system degassing. Dissolved sulfur concentrations in melt inclusions (MIs) from pre-climactic and climactic rhyodacitic pumices varies from 80 to 330 ppm, with highest concentrations in inclusions with 4.8-5.2 wt% H₂O (by FTIR). Up to 50% of the initial S may have been lost through pre-eruptive degassing at depths of 4-5 km. Ion microprobe analyses of pyrrhotite in climactic rhyodacitic tephra and andesitic scoria indicate a range in δ³⁴S from −0.4‰ to 5.8‰ and from −0.1‰ to 3.5‰, respectively. Initial δ³⁴S values of rhyodacitic and andesitic magmas were likely near the mantle value of 0‰. Hydrogen isotope (δD) and total H₂O analyses of rhyodacitic obsidian (and vitrophyre) from the climactic fall deposit yielded values οf −103 to −53‰ and 0.23-1.74 wt%, respectively. Values of δD and wt% H₂O of obsidian decrease towards the top of the fall deposit. Samples with depleted δD, and mantle δ¹⁸O values, have elevated δ³⁴S values consistent with open-system degassing. These results imply that more mantle-derived sulfur is degassed to the Earth’s atmosphere/hydrosphere through convergent margin volcanism than previously attributed. Magmatic degassing can modify initial isotopic compositions of sulfur by >14‰ (to δ³⁴S values of 14‰ or more here) and hydrogen isotopic compositions by 90‰ (to δD values of −127‰ in this case).     

Silicon isotopes in meteorites and planetary core formation

The silicon (Si) isotope compositions of 42 meteorite and terrestrial samples have been determined using MC-ICPMS with the aim of resolving the current debate over their compositions and the implications for core formation. No systematic δ³⁰Si differences are resolved between chondrites (δ³⁰Si = −0.49 ± 0.15‰, 2σ_SD) and achondrites (δ³⁰Si = −0.47 ± 0.11‰, 2σ_SD), although enstatite chondrites are consistently lighter (δ³⁰Si = −0.63 ± 0.07‰, 2σ_SD) in comparison to other meteorite groups. The data reported here for meteorites and terrestrial samples display an average difference Δ³⁰Si_{BSE−meteorite∗} = 0.15 ± 0.10‰, which is consistent within uncertainty with previous studies. No effect from sample heterogeneity, preparation, chemistry or mass spectrometry can be identified as responsible for the reported differences between current datasets. The heavier composition of the bulk silicate Earth is consistent with previous conclusions that Si partitioned into the metal phase during metal-silicate equilibration at the time of core formation. Fixing the temperature of core formation to the peridotite liquidus and using an appropriate metal silicate fractionation factor (ε ∼0.89), the Δ³⁰Si_{BSE−meteorite∗} value from this study indicates that the Earth core contains at least 2.5 and possibly up to 16.8 wt% Si.     

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.     

Mineral inclusions in diamonds from the River Ranch kimberlite, Zimbabwe

More than 99% of mineral inclusions in diamonds from the River Ranch pipe in the Late Archean Limpopo Mobile Belt (Zimbabwe), are phases of harzburgitic paragenesis, namely olivine (Fo_92–93), orthopyroxene (Mg# = 93), G10 garnets and chromites. The diamond inclusion (DI) chemistry demonstrates a limited overlap with River Ranch kimberlite macrocrysts: the DI garnets are more Ca-undersaturated, and DI spinel and garnet are more Mg-rich. Most River Ranch diamond inclusions were equilibrated at T = 1080–1320 °C, P = 47–61 kbar, and f _O2 between IW and WM buffers. The P/T profile beneath the Limpopo Mobile Belt (LMB) is consistent with a paleo-heat flow of 41–42 mW/m², similar to calculations for Roberts Victor, but hotter than for the Finsch, Kimberley, Koffiefontein and Premier Mines. This is ascribed to the younger tectonothermal age of the LMB and its proximity to Late Archean oceans. Like diamond inclusions from all other kimberlites studied, the River Ranch DI have a lithospheric affinity and therefore indicate that an ancient, chemically depleted, thick (at least 200 km) mantle root existed beneath the Limpopo Mobile Belt 530–540 Ma ago. The mantle root might have developed beneath the continental Central Zone of the LMB as early as the Archean, and could be alien to the overthrust allochthonous sheet of the Limpopo Belt. Oxygen fugacity estimates for diamond inclusions at River Ranch are similar to other diamondiferous harzburgites beneath the Kaapvaal craton, indicating that the Kaapvaal mantle as a whole was well buffered and homogeneous with respect to f _O2 at the time of peridotitic diamond crystallization. Received: 11 January 1995 / Accepted: 10 June 1997     

Melting of enriched Archean subcontinental lithospheric mantle: Evidence from the ca. 1760 Ma volcanic rocks of the Xiong’er Group,southern margin of the North China Craton

The Xiong’er Group is an important geologic unit in the southern margin of the North China Craton. It is dominated by the volcanic rocks, dated at 1763 ± 15 Ma, that have SiO₂ contents ranging from 52.10 wt% to 73.51 wt%. These volcanic rocks are sub-alkaline and can be classified into three subgroups: basaltic andesites, andesites and rhyolites. They unexceptionally show enrichment of light rare earth elements (LREE) and share similar trace element patterns. Depletions in Nb, Ta, Sr, P and Ti relative to the adjacent elements are evident for all the samples. The volcanic rocks are evolved with low MgO contents (0.29–5.88 wt%) and accordingly low Mg^# values of 11–53. The Nd isotopes are enriched and show a weak variation with ?_Nd(t) = −7.12 to −9.63. Zircon Hf isotopes are also enriched with ?_Hf(t) = −12.02 ± 0.45. The volcanic rocks of the Xiong’er Group are interpreted to represent fractional crystallization of a common mantle source. The volcanic rocks might have been generated by high-degree partial melting of a lithospheric mantle that was originally modified by the oceanic subduction in the Late Archean. This brings a correlation with the subduction-modified lithospheric mantle in an extensional setting during breakup of the Columbia supercontinent in the late Paleoproterozoic, rather than in an arc setting. The elevated SiO₂ contents and evolved radiogenic isotope features indicate the possible incorporation into their source of lower crustal materials that have similar Nd isotopic characteristics to the subcontinental lithospheric mantle. The existence of extensive Xiong’er volcanic rocks (60,000 km²) indicates an early large-scale subduction-related metasomatism in the area and probably suggest a flat subduction model for the plate-margin magmatism in the Late Archean.     

Copper isotopic zonation in the Northparkes porphyry Cu-Au deposit, SE Australia

Significant, systematic Cu isotopic variations have been found in the Northparkes porphyry Cu-Au deposit, NSW, Australia, which is an orthomagmatic porphyry Cu deposit. Copper isotope ratios have been measured in sulfide minerals (chalcopyrite and bornite) by both solution and laser ablation multi-collector inductively coupled plasma-mass spectrometry (MC-ICP-MS). The results from both methods show a variation in δ⁶⁵Cu of hypogene sulfide minerals of greater than 1‰ (relative to NIST976). Significantly, the results from four drill holes through two separate ore bodies show strikingly similar patterns of Cu isotope variation. The patterns are characterized by a sharp down-hole decrease from up to 0.8‰ (0.29 ± 0.56‰, 1σ, n = 20) in the low-grade peripheral alteration zones (phyllic-propylitic alteration zone) to a low of ∼−0.4‰ (−0.25 ± 0.36‰, 1σ, n = 30) at the margins of the most mineralized zones (Cu grade >1 wt%). In the high-grade cores of the systems, the compositions are more consistent at around 0.2‰ (0.19 ± 0.14‰, 1σ, n = 40). The Cu isotopic zonation may be explained by isotope fractionation of Cu between vapor, solution and sulfides at high temperature, during boiling and sulfide precipitation processes. Sulfur isotopes also show an isotopically light shell at the margins of the high-grade ore zones, but these are displaced from the low δ⁶⁵Cu shells, such that there is no correlation between the Cu and S isotope signatures. Fe isotope data do not show any discernable variation along the drill core. This work demonstrates that Cu isotopes show a large response to high-temperature porphyry mineralizing processes, and that they may act as a vector to buried mineralization.     

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ⁶⁵Cu values (based on ⁶⁵Cu/⁶³Cu) of enargite (δ⁶⁵Cu = −0.01 ± 0.10‰; 2σ) and chalcopyrite (δ⁶⁵Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ⁶⁵Cu < 1‰). These mineral samples show lower δ⁶⁵Cu values than stream waters (1.38‰ ? δ⁶⁵Cu ? 1.69‰). The average isotopic fractionation (Δ_aq-min = δ⁶⁵Cu_aq − δ⁶⁵Cu_min, where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in ⁶⁵Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite (Δ_aq-min^o=-0.57±0.14‰, where min^o refers to the starting mineral) and no apparent fractionation for enargite (Δ_aq-min^o=0.14±0.14‰). Abiotic fractionation is attributed to preferential oxidation of ⁶⁵Cu⁺ at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of ⁶⁵Cu_aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ⁶⁵Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.     

Tissue N content and N natural abundance in epilithic mosses for indicating atmospheric N deposition in the Guiyang area,SW China

Tissue N contents and δ¹⁵N signatures in 175 epilithic mosses were investigated from urban to rural sites in Guiyang (SW China) to determine atmospheric N deposition. Moss N contents (0.85–2.97%) showed a significant decrease from the urban area (mean = 2.24 ± 0.32%, 0–5 km) to the rural area (mean = 1.27 ± 0.13%, 20–25 km), indicating that the level of N deposition decreased away from the urban environment, while slightly higher N contents re-occurred at sites beyond 30 km, suggesting higher N deposition in more remote rural areas. Moss δ¹⁵N ranged from −12.50‰ to −1.39‰ and showed a clear bimodal distribution (−12‰ to −6‰ and −5‰ to −2‰), suggesting that there are two main sources for N deposition in the Guiyang area. More negative δ¹⁵N (mean = −8.87 ± 1.65‰) of urban mosses mainly indicated NH₃ released from excretory wastes and sewage, while the less negative δ¹⁵N (from −3.83 ± 0.82‰ to −2.48 ± 0.95‰) of rural mosses were mainly influenced by agricultural NH₃. With more negative values in the urban area than in the rural area, the pattern of moss δ¹⁵N variation in Guiyang was found to be opposite to cities where N deposition is dominated by NO_x–N. Therefore, NH_x–N is the dominant N form deposited in the Guiyang area, which is supported by higher NH_x–N than NO_x–N in local atmospheric deposition. From the data showing that moss is responding to NH_x–N/NO_x–N in deposition it can be further demonstrated that the variation of moss δ¹⁵N from the Guiyang urban to rural area was more likely controlled by the ratio of urban-NH_x/agriculture-NH_x than the ratio of NH_x–N/NO_x–N. The results of this study have extended knowledge of atmospheric N sources in city areas, showing that urban sewage discharge could be important in cities co-generic to Guiyang.