The origin of the fumaroles of La Solfatara (Campi Flegrei, South Italy)

The analysis of gaseous compositions from Solfatara (Campi Flegrei, South Italy) fumaroles since the early 1980s, clearly reveals a double thermobarometric signature. A first signature at temperatures of about 360 °C was inferred by methane-based chemical-isotopic geoindicators and by the H₂/Ar geothermometer. These high temperatures, close to the critical point of water, are representative of a deep zone where magmatic gases flash the hydrothermal liquid, forming a gas plume. A second signature was found to be at around 200-240 °C. At these temperatures, the kinetically fast reactive species (H₂ and CO) re-equilibrate in a pure vapor phase during the rise of the plume. A combination of these observations with an original interpretation of the oxygen isotopic composition of the two dominant species, i.e. H₂O and CO₂, shed light on the origin of fumarolic fluids by showing that effluents are mixture between fluids degassed from a magma body and the vapor generated at about 360 °C by the vaporization of hydrothermal liquids. A typical ‘andesitic’ water type (δD ∼ −20‰, δ¹⁸O ∼10‰) and a CO₂-rich composition (X_CO2∼0.4) has been inferred for the magmatic fluids, while for the hydrothermal component a meteoric origin and a CO₂ fugacity fixed by fluid-rock reaction at high temperatures have been estimated. In the time the fraction of magmatic fluids in the fumaroles increased (up to ∼0.5) at each seismic and ground uplift crisis (bradyseism) which occurred at Campi Flegrei, suggesting that bradyseismic crises are triggered by periodic injections of CO₂-rich magmatic fluids at the bottom of the hydrothermal system.     

Insights into C and H storage in the altered oceanic crust: Results from ODP/IODP Hole 1256D

Carbon and hydrogen concentrations and isotopic compositions were measured in 19 samples from altered oceanic crust cored in ODP/IODP Hole 1256D through lavas, dikes down to the gabbroic rocks. Bulk water content varies from 0.32 to 2.14 wt% with δD values from −64‰ to −25‰. All samples are enriched in water relative to fresh basalts. The δD values are interpreted in terms of mixing between magmatic water and another source that can be either secondary hydrous minerals and/or H contained in organic compounds such as hydrocarbons. Total CO₂, extracted by step-heating technique, ranges between 564 and 2823 ppm with δ¹³C values from −14.9‰ to −26.6‰. As for water, these altered samples are enriched in carbon relative to fresh basalts. The carbon isotope compositions are interpreted in terms of a mixing between two components: (1) a carbonate with δ¹³C = −4.5‰ and (2) an organic compound with δ¹³C = −26.6‰. A mixing model calculation indicates that, for most samples (17 of 19), more than 75% of the total C occurs as organic compounds while carbonates represent less than 25%. This result is also supported by independent estimates of carbonate content from CO₂ yield after H₃PO₄ attack. A comparison between the carbon concentration in our samples, seawater DIC (Dissolved Inorganic Carbon) and DOC (Dissolved Organic Carbon), and hydrothermal fluids suggests that CO₂ degassed from magmatic reservoirs is the main source of organic C addition to the crust during the alteration process. A reduction step of dissolved CO₂ is thus required, and can be either biologically mediated or not. Abiotic processes are necessary for the deeper part of the crust (>1000 mbsf) because alteration temperatures are greater than any hyperthermophilic living organism (i.e. T > 110 °C). Even if not required, we cannot rule out the contribution of microbial activity in the low-temperature alteration zones. We propose a two-step model for carbon cycling during crustal alteration: (1) when “fresh” oceanic crust forms at or close to ridge axis, alteration starts with hot hydrothermal fluids enriched in magmatic CO₂, leading to the formation of organic compounds during Fischer-Tropsch-type reactions; (2) when the crust moves away from the ridge axis, these interactions with hot hydrothermal fluids decrease and are replaced by seawater interactions with carbonate precipitation in fractures. Taking into account this organic carbon, we estimate C isotope composition of mean altered oceanic crust at ∼ −4.7‰, similar to the δ¹³C of the C degassed from the mantle at ridge axis, and discuss the global carbon budget. The total flux of C stored in the altered oceanic crust, as carbonate and organic compound, is 2.9 ± 0.4 × 10¹² molC/yr.     

Hydrogen isotopes in volcanic plumes: Tracers for remote temperature sensing of fumaroles

In high-temperature volcanic fumaroles (>400 °C), the isotopic composition of molecular hydrogen (H₂) reaches equilibrium with that of the fumarolic H₂O. In this study, we used this hydrogen isotope exchange equilibrium of fumarolic H₂ as a tracer for the remote temperature at volcanic fumaroles. In this remote sensing, we deduced the hydrogen isotopic composition (δD value) of fumarolic H₂ from those in the volcanic plume. To ascertain that we can estimate the δD value of fumarolic H₂ from those in a volcanic plume, we estimated the values in three fumaroles with outlet temperatures of 630 °C (Tarumae), 203 °C (Kuju), and 107 °C (E-san). For this we measured the concentration and δD value of H₂ in each volcanic plume, along with those determined directly at each fumarole. The average and maximum mixing ratios of fumarolic H₂ within a plume’s total H₂ were 97% and 99% (at Tarumae), 89% and 96% (at Kuju), and 97% and 99% (at E-san). We found a linear relationship between the depletion in the δD values of H₂, with the reciprocal of H₂ concentration. Furthermore, the estimated end-member δD value for each H₂-enriched component (−260 ± 30‰ vs. VSMOW in Tarumae, −509 ± 23‰ in Kuju, and −437 ± 14‰ in E-san) coincided well with those observed at each fumarole (−247.0 ± 0.6‰ in Tarumae, −527.7 ± 10.1‰ in Kuju, and −432.1 ± 2.5‰ in E-san). Moreover, the calculated isotopic temperatures at the fumaroles agreed to within 20 °C with the observed outlet temperature at Tarumae and Kuju. We deduced that the δD value of the fumarolic H₂ was quenched within the volcanic plume. This enabled us to remotely estimate these in the fumarole, and thus the outlet temperature of fumaroles, at least for those having the outlet temperatures more than 400 °C. By applying this methodology to the volcanic plume emitted from the Crater 1 of Mt. Naka-dake (the volcano Aso) where direct measurement on fumaroles was impractical, we estimated that the δD value of the fumarolic H₂ to be −172 ± 16‰ and the outlet temperature to be 868 ± 97 °C. The remote temperature sensing using hydrogen isotopes developed in this study is widely applicable to many volcanic systems.     

Geochemical characterization of secondary microbial gas occurrence in the Songliao Basin, NE China

A suite of natural gases from the northern Songliao Basin in NE China were characterized for their molecular and carbon isotopic composition. Gases from shallow reservoirs display clear geochemical evidence of alteration by biodegradation, with very high dryness (C₁/C₂₊ > 100), high C₂/C₃ and i-C₄/n-C₄ ratios, high nitrogen content and variable carbon dioxide content. Isotopic values show wide range variations (δ¹³C_CH4 from −79.5‰ to −45.0‰, δ¹³C_C2H6 from −53.7‰ to −32.2‰, δ¹³C_C3H8 from −36.5‰ to −20.1‰, δ¹³C_nC4H10 from −32.7‰ to −24.5‰, and δ¹³C_CO2 from −21.6‰ to +10.5‰). A variety of genetic types can be recognized on the basis of chemical and isotopic composition together with their geological occurrence. Secondary microbial gas generation was masked by primary microbial gas and the mixing of newly generated methane with thermogenic methane already in place in the reservoir can cause very complicated isotopic signatures. System openness also was considered for shallow biodegraded gas accumulations. Gases from the Daqing Anticline are relatively wet with ¹³C enriched methane and ¹³C depleted CO₂, representing typically thermogenic origin. Gases within the Longhupao-Da’an Terrace have variable dryness, ¹³C enriched methane and variable δ¹³C of CO₂, suggesting dominant thermogenic origin and minor secondary microbial methane augment. The Puqian-Ao’nan Uplift contains relatively dry gas with ¹³C depleted methane and ¹³C enriched CO₂, typical for secondary microbial gas with a minor part of thermogenic methane. Gas accumulations in the Western Slope are very dry with low carbon dioxide concentrations. Some gases contain ¹³C depleted methane, ethane and propane, indicating low maturity/primary microbial origin. Recognition of varying genetic gas types in the Songliao Basin helps explain the observed dominance of gas in the shallow reservoir and could serve as an analogue for other similar shallow gas systems.     

Subsurface processes at the lucky strike hydrothermal field, Mid-Atlantic ridge: evidence from sulfur, selenium, and iron isotopes

At Lucky Strike near the Azores Triple Junction, the seafloor setting of the hydrothermal field in a caldera system with abundant low-permeability layers of cemented breccia, provides a unique opportunity to study the influence of subsurface geological conditions on the hydrothermal fluid evolution. Coupled analyses of S isotopes performed in conjunction with Se and Fe isotopes have been applied for the first time to the study of seafloor hydrothermal systems. These data provide a tool for resolving the different abiotic and potential biotic near-surface hydrothermal reactions. The δ³⁴S (between 1.5‰ and 4.6‰) and Se values (between 213 and 1640 ppm) of chalcopyrite suggest a high temperature end-member hydrothermal fluid with a dual source of sulfur: sulfur that was leached from basaltic rocks, and sulfur derived from the reduction of seawater sulfate. In contrast, pyrite and marcasite generally have lower δ³⁴S within the range of magmatic values (0 ± 1‰) and are characterized by low concentrations of Se (<50 ppm). For ⁸²Se/⁷⁶Se ratios, the δ⁸²Se values range from basaltic values of near −1.5‰ to −7‰. The large range and highly negative values of hydrothermal deposits observed cannot be explained by simple mixing between Se leached from igneous rock and Se derived from seawater. We interpret the Se isotope signature to be a result of leaching and mixing of a fractionated Se source located beneath hydrothermal chimneys in the hydrothermal fluid. At Lucky Strike we consider two sources for S and Se: (1) the “end-member” hydrothermal fluid with basaltic Se isotopic values (−1.5‰) and typical S isotope hydrothermal values of 1.5‰; (2) a fractionated source hosted in subsurface environment with negative δ³⁴S values, probably from bacterial reduction of seawater sulfate and negative δ⁸²Se values possibly derived from inorganic reduction of Se oxyanions. Fluid trapped in the subsurface environment is conductively cooled and has restricted mixing and provide favorable conditions for subsurface microbial activity which is potentially recorded by S isotopes. Fe isotope systematic reveals that Se-rich high temperature samples have δ⁵⁷Fe values close to basaltic values (∼0‰) whereas Se-depleted samples precipitated at medium to low temperature are systematically lighter (δ⁵⁷Fe values between −1 to −3‰). An important implication of our finding is that light Fe isotope composition down to −3.2‰ may be explained entirely by abiotic fractionation, in which a reservoir effect during sulfide precipitation was able to produce highly fractionated compositions.     

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

S-33 constraints on the seawater sulfate contribution in modern seafloor hydrothermal vent sulfides

Sulfide sulfur in mid-oceanic ridge hydrothermal vents is derived from leaching of basaltic-sulfide and seawater-derived sulfate that is reduced during high temperature water rock interaction. Conventional sulfur isotope studies, however, are inconclusive about the mass-balance between the two sources because ³⁴S/³²S ratios of vent fluid H₂S and chimney sulfide minerals may reflect not only the mixing ratio but also isotope exchange between sulfate and sulfide. Here, we show that high-precision analysis of S-33 can provide a unique constraint because isotope mixing and isotope exchange result in different Δ³³S (≡δ³³S-0.515 δ³⁴S) values of up to 0.04‰ even if δ³⁴S values are identical. Detection of such small Δ³³S differences is technically feasible by using the SF₆ dual-inlet mass-spectrometry protocol that has been improved to achieve a precision as good as 0.006‰ (2σ).Sulfide minerals (marcasite, pyrite, chalcopyrite, and sphalerite) and vent H₂S collected from four active seafloor hydrothermal vent sites, East Pacific Rise (EPR) 9-10°N, 13°N, and 21°S and Mid-Atlantic Ridge (MAR) 37°N yield Δ³³S values ranging from −0.002 to 0.033 and δ³⁴S from −0.5‰ to 5.3‰. The combined δ³⁴S and Δ³³S systematics reveal that 73 to 89% of vent sulfides are derived from leaching from basaltic sulfide and only 11 to 27% from seawater-derived sulfate. Pyrite from EPR 13°N and marcasite from MAR 37°N are in isotope disequilibrium not only in δ³⁴S but also in Δ³³S with respect to associated sphalerite and chalcopyrite, suggesting non-equilibrium sulfur isotope exchange between seawater sulfate and sulfide during pyrite precipitation. Seafloor hydrothermal vent sulfides are characterized by low Δ³³S values compared with biogenic sulfides, suggesting little or no contribution of sulfide from microbial sulfate reduction into hydrothermal sulfides at sediment-free mid-oceanic ridge systems. We conclude that ³³S is an effective new tracer for interplay among seawater, oceanic crust and microbes in subseafloor hydrothermal sulfur cycles.     

Reactions governing the chemistry of crater fumaroles from Vulcano Island,Italy, and implications for volcanic surveillance

More than 200 chemical and isotope analyses of fumarolic fluids collected at the Fossa Grande crater, Vulcano Island, during the 1980s show that the main process controlling these fluids is mixing between the gas released by a magma body and the vapour produced through evaporation of brines of marine origin. Large variations in the relative contribution of these two sources have been observed during the last 10 a. The main species (H₂O and CO₂), the inert gases (He and N₂), and the D content of steam are fixed by the mixing processes; they are therefore the best tracers the fraction of the deep magmatic component in the fumarolic fluids discharged at the surface. In contrast, the “fast” species (H₂ and CO) equilibrate at T,P values close to the outlet temperature and atmospheric pressure, and under redox conditions governed by the SO₂H₂S buffer, as indicated by thermodynamic calculations.Acid gases (HCl, HF, H₂S and SO₂) are partly contributed by the magmatic component and partly produced by the reactions between hot rocks, steam and salts which take place in the “dry” zones surrounding the central magmatic gas column, as suggested by the good agreement between their analytical and theoretical contents.     

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.     

Magmatic fluids of Tatun volcanic group,Taiwan

Two distinctive magmatic fluids were recognized in the Tatun volcanic group (TVG), Taiwan. One is a relatively reduced fluid represented by the fumarolic gases at Hsiao-you-ken (HYK) geothermal field. Another is an oxidized fluid containing high concentrations of HCl represented by the fumarolic gases at Da-you-ken (DYK). An intermediate gas was recognized at Gung-tze-ping (GTP) and She-hung-ping (SHP). The fumarolic gases at HYK and GTP possess the features of so-called primary steam generated on mixing of magmatic gas and meteoric groundwater. The fumarolic gases at DYK are a simple mixture between magmatic gas and water vapor of meteoric origin. The CO₂/H₂O molar ratio of the magmatic component in the fumarolic gases at DYK was estimated to be 0.018, meanwhile it was estimated to be 0.027 for the fumarolic gases at HYK and GTP, suggesting the magma beneath DYK is depleted in volatiles relative to the magma beneath HYK and GTP. The estimated CO₂/H₂O ratio for the magmatic component is comparable to that of some active volcanoes in Japan, suggesting the enrichment of volatiles in the magmas beneath TVG.     

Geochemical features and origin of natural gas in heavy oil area of the Western Slope,Songliao Basin,China

The Western Slope of the Songliao Basin is rich in heavy oil resources (>70 × 10⁸ bbl), around which there are shallow gas reservoirs (∼1.0 × 10¹² m³). The gas is dominated by methane with a dryness over 0.99, and the non-hydrocarbon component being overwelmingly nitrogen. Carbon isotope composition of methane and its homologs is depleted in ¹³C, with δ¹³C₁ values being in the range of −55‰ to −75‰, δ¹³C₂ being in the range of −40‰ to −53‰ and δ¹³C₃ being in the range of −30‰ to −42‰, respectively. These values differ significantly from those solution gases source in the Daqing oilfield. This study concludes that heavy oils along the Western Slope were derived from mature source rocks in the Qijia-Gulong Depression, that were biodegraded. The low reservoir temperature (30–50 °C) and low salinity of formation water with neutral to alkaline pH (NaHCO₃) appeared ideal for microbial activity and thus biodegradation. Natural gas along the Western Slope appears mainly to have originated from biodegradation and the formation of heavy oil. This origin is suggested by the heavy δ¹³C of CO₂ (−18.78‰ to 0.95‰) which suggests that the methane was produced via fermentation as the terminal decomposition stage of the oil.     

Sulfur isotope signatures in gypsiferous sediments of the Estancia and Tularosa Basins as indicators of sulfate sources, hydrological processes, and microbial activity

In order to reconstruct paleo-environmental conditions for the saline playa lakes of the Rio Grande Rift, we investigated sediment sulfate sources using sulfur isotope compositions of dissolved ions in modern surface water, groundwater, and precipitated in the form of gypsum sediments deposited during the Pleistocene and Holocene in the Tularosa and Estancia Basins. The major sulfate sources are Lower and Middle Permian marine evaporites (δ³⁴S of 10.9-14.4‰), but the diverse physiography of the Tularosa Basin led to a complex drainage system which contributed sulfates from various sources depending on the climate at the time of sedimentation. As inferred from sulfur isotope mass balance constraints, weathering of sulfides of magmatic/hydrothermal and sedimentary origin associated with climate oscillations during Last Glacial Maximum contributed about 35-50% of the sulfates and led to deposition of gypsum with δ³⁴S values of −1.2‰ to 2.2‰ which are substantially lower than Permian evaporates. In the Estancia Basin, microbial sulfate reduction appears to overprint sulfur isotopic signatures that might elucidate past groundwater flows. A Rayleigh distillation model indicates that about 3-18% of sulfates from an inorganic groundwater pool (δ³⁴S of 12.6-13.8‰) have been metabolized by bacteria and preserved as partially to fully reduced sulfur-bearing minerals species (elemental sulfur, monosulfides, disulfides) with distinctly negative δ³⁴S values (−42.3‰ to −20.3‰) compared to co-existing gypsum (−3.8‰ to 22.4‰). For the Tularosa Basin microbial sulfate reduction had negligible effect on δ³⁴S value of the gypsiferous sediments most likely because of higher annual temperatures (15-33 °C) and lower organic carbon content (median 0.09%) in those sediments leading to more efficient oxidation of H₂S and/or smaller rates of sulfate reduction compared to the saline playas of the Estancia Basin (5-28 °C; median 0.46% of organic carbon).The White Sands region of the Tularosa Basin is frequently posited as a hydrothermal analogue for Mars. High temperatures of groundwater (33.3 °C) and high δ¹⁸O(H₂O) values (1.1‰) in White Sands, however, are controlled predominantly by seasonal evaporation rather than the modern influx of hydrothermal fluids. Nevertheless, it is possible that some of the geochemical processes in White Sands, such as sulfide weathering during climate oscillations and upwelling of highly mineralized waters, might be considered as valid terrestrial analogues for the sulfate cycle in places such as Meridiani Planum on Mars.     

Sulfur isotope systematics of the 1982 El Chichón trachyandesite: an ion microprobe study

Sulfur isotopic compositions were determined by ion microprobe for 36 spots on anhydrite crystals in trachyandesitic pumices erupted from El Chichón Volcano in 1982. Individual anhydrite crystals are homogeneous in δ³⁴S, within the ±1‰ (2σ) uncertainty of the method, but crystal-to-crystal variations are large (+2.5 to +10.9‰). The mean δ³⁴S for anhydrite (+6.4 ± 2.1‰, 1σ) is significantly lower than earlier results for bulk anhydrite separates (+9.0 to +9.2‰). The difference between the mean δ³⁴S values in these two populations may reflect a grain-size effect, with heavier sulfur concentrated in smaller anhydrite crystals, few of which were analyzed by ion microprobe. Variations in δ³⁴S show no correlation with complex textures in anhydrite revealed by cathodoluminescence color. Ion-microprobe analyses of δ³⁴S were also obtained on six ovoid-shaped inclusions of pyrrhotite, chalcopyrite, and/or intermediate sulfide solid solution hosted by silicate or oxide crystals, interpreted to be magmatic (δ³⁴S = −0.1 to +2.7‰; mean +0.7‰), and on four irregularly shaped multiphase sulfide fragments in the matrix, interpreted as xenocrystic, which range widely in δ³⁴S (−3.7 to +5.5‰). We evaluate four different mixing scenarios involving (1) magmatic anhydrite and sedimentary sulfate, (2) magmatic anhydrite and hydrothermal anhydrite, and anhydrite and coexisting sulfide crystals precipitated in different domains of a common magma reservoir that were affected by (3) different degrees of degassing or (4) different degrees of crustal sulfur contamination. The model involving physical contamination of sedimentary sulfate is considered untenable. The other three models are considered to be viable, but none of them can explain all observations. The results of this study and other recent investigations prompt a re-evaluation of the sulfur budget for the 1982 El Chichón eruption. We estimate that 2.2 × 10¹³ g of S was emitted, and that 58 wt.% of the sulfur was present as anhydrite prior to eruption, with the remainder in a vapor phase, with H₂S/SO₂ ≈ 9. The bulk magmatic δ³⁴S value for the 1982 El Chichón trachyandesite is estimated as +4.1 to +5.8‰, typical of the relatively heavy sulfur isotopic compositions that characterize subduction-related magmas.     

Sr, C, and O isotope geochemistry of Ordovician brachiopods: a major isotopic event around the Middle-Late Ordovician transition

Here we present Sr, C, and O isotope curves for Ordovician marine calcite based on analyses of 206 calcitic brachiopods from 10 localities worldwide. These are the first Ordovician-wide isotope curves that can be placed within the newly emerging global biostratigraphic framework. A total of 182 brachiopods were selected for C and O isotope analysis, and 122 were selected for Sr isotope analysis. Seawater ⁸⁷Sr/⁸⁶Sr decreased from 0.7090 to 0.7078 during the Ordovician, with a major, quite rapid fall around the Middle-Late Ordovician transition, most probably caused by a combination of low continental erosion rates and increased submarine hydrothermal exchange rates. Mean δ¹⁸O values increase from −10‰ to −3‰ through the Ordovician with an additional short-lived increase of 2 to 3‰ during the latest Ordovician due to glaciation. Although diagenetic alteration may have lowered δ¹⁸O in some samples, particularly those from the Lower Ordovician, maximum δ¹⁸O values, which are less likely to be altered, increase by more than 3‰ through the Ordovician in both our data and literature data. We consider that this long-term rise in calcite δ¹⁸O records the effect of decreasing tropical seawater temperatures across the Middle-Late Ordovician transition superimposed on seawater δ¹⁸O that was steadily increasing from ≤−3‰ standard mean ocean water (SMOW). By contrast, δ¹³C variation seems to have been relatively modest during most of the Ordovician with the exception of the globally documented, but short-lived, latest Ordovician δ¹³C excursion up to +7‰. Nevertheless, an underlying trend in mean δ¹³C can be discerned, changing from moderately negative values in the Early Ordovician to moderately positive values by the latest Ordovician. These new isotopic data confirm a major reorganization of ocean chemistry and the surface environment around 465 to 455 Ma. The juxtaposition of the greatest recorded swings in Phanerozoic seawater ⁸⁷Sr/⁸⁶Sr and δ¹⁸O at the same time as one of the largest marine transgressions in Phanerozoic Earth history suggests a causal link between tectonic and climatic change, and emphasizes an endogenic control on the O isotope budget during the Early Paleozoic. Better isotopic and biostratigraphic constraints are still required if we are to understand the true significance of these changes. We recommend that future work on Ordovician isotope stratigraphy focus on this outstanding Middle-Late Ordovician event.     

Stable isotope compositions of tourmalines from granites and related hydrothermal rocks of the Karagwe-Ankolean belt,northwest Tanzania

Tourmaline is a ubiquitous mineral in the Mid-Proterozoic, peraluminous, syn- to post-tectonic granites and aplites and the related hydrothermal rocks of the Karagwe-Ankolean belt in northwest Tanzania. Electron microprobe analysis indicates that tourmalines from all of the intrusive and hydrothermal lithologies: (1) belong to the schorl-dravite solid-solution series; and (2) plot within the field occupied by tourmaline from Li-poor granitoids on the Fe-Al-Mg classification diagram. Oxygen isotope compositions range from +12.2 to +11.6‰ (SMOW) for magmatic tourmalines and from +10.8 to +9.8‰ for those of hydrothermal origin. Hydrogen isotope compositions vary from −79 to − 65‰ (SMOW) for magmatic tourmalines and from −99 to −84‰ for hydrothermal tourmalines. Water contents measured by manometry are constant at 3.0–3.2 wt.%. Within the broad grouping there arc systematic variations in both chemical [particularly Fe_tot/(Fe_tot + Mg ratio)] and isotopic composition that relate to evolving magmatic and hydrothermal conditions. Igneous differentiation [increasing Fe_tot/(Fe_tot + Mg) in magmatic tourmaline] has produced trends with higher δ¹⁸O in quartz, lower δ¹⁸O in tourmaline, and larger Δ_QTZ.−TOUR.-values, that reflect a combination of a reduction of crystallization temperature and an increase of Fe_tot/ (Fe_tot + Mg) in the residual melt. Subsequent cooling and interaction of an exsolved, B-rich magmatic fluid with the pelitic country rocks, resulted in the deposition of hydrothermal tourmaline with increasing Fe_tot/(Fe_tot + Mg) ratios, and progressively lower δ¹⁸O and δD -values.     

Extreme oxygen isotope signature of meteoric water in magmatic zircon from metagranite in the Sulu orogen, China: Implications for Neoproterozoic rift magmatism

Unusual ¹⁸O depletion, with δ¹⁸O values as negative as −10‰ to −4‰ relative to VSMOW, was reported in zircons from ultrahigh-pressure eclogite-facies metamorphic rocks in the Dabie-Sulu orogenic belt, China. But it is critical for the negative δ¹⁸O zircons to be distinguished between magmatic and metamorphic origins, because the ¹⁸O depletion can be acquired by high-T eclogite-facies metamorphism of meteoric-hydrothermally altered low δ¹⁸O rocks. While zircon O diffusion kinetics has placed a reasonable constraint on this, zircon trace element compositions can provide a straightforward distinction between the magmatic and metamorphic origins. This paper reports our finding of unusual ¹⁸O depletion in zircon from granitic gneiss in the northeastern end of the Sulu orogen. Zircon δ¹⁸O values vary from −7.8‰ to −3.1‰ along a profile of 50 m length at Zaobuzhen. They are close to extremely low δ¹⁸O values of −9.0‰ to −5.9‰ for metagranite at Qinglongshan and adjacent areas in the southwestern end of the Sulu orogen. CL imaging suggests that the low δ¹⁸O zircons at Zaobuzhen are primarily of magmatic origin, but underwent different degrees of metamorphic modification. Zircon U-Pb dating yields middle Neoproterozoic ages of 751 ± 27 to 779 ± 25 Ma for protolith crystallization and Triassic ages of 214 ± 10 to 241 ± 33 Ma for metamorphic resetting. However, no metamorphic modification occurs in zircon REE patterns that only indicate magmatic recrystallization and hydrothermal alteration, respectively. Thus, the negative δ¹⁸O zircons are interpreted as crystallizing from negative δ¹⁸O magmas due to melting of meteoric-hydrothermally altered negative δ¹⁸O rocks in an active rift setting at about 780 Ma. The variation in zircon δ¹⁸O values indicates considerable O isotope heterogeneity in its granitic protolith. Zircon Lu-Hf isotope analyses give positive ε_Hf(t) values of 1.6-4.1 and Hf model ages of 1.18-1.30 Ga. This suggests that the granitic protolith was derived from the mid-Neoproterozoic reworking of late Mesoproterozoic juvenile crust. The metagranites at Zaobuzhen and Qinglongshan, about 450 km apart, are two known occurrences of the unusually low δ¹⁸O zircons below −6‰ so far reported in the Sulu orogen. They are similar to each other in both protolith and metamorphic ages, so that they share the same nature of both Neoproterozoic protolith and Triassic metamorphism. Therefore, the locally negative δ¹⁸O zircons may register centers of low δ¹⁸O magmatism during the supercontinental rifting.     

Graphite-bearing CO2-fluid inclusions in granulites: Insights on graphite precipitation and carbon isotope evolution

Graphite in deep crustal enderbitic (orthopyroxene + garnet + plagioclase + quartz) granulites (740°C, 8.9 kb) of Nilgiri hills, southern India were investigated for their spectroscopic and isotopic characteristics. Four types of graphite crystals were identified. The first type (Gr_I), which is interstitial to other mineral grains, can be grouped into two subtypes, Gr_IA and Gr_IB. Gr_IA is either irregular in shape or deformed, and rough textured with average δ¹³C values of −12.7 ± 0.4‰ (n = 3). A later generation of interstitial graphite (Gr_IB) shows polygonal crystal shapes and highly reflecting smooth surface features. These graphite grains are more common and have δ¹³C values of −11.9 ± 0.3‰ (n = 14). Both subtypes show well-defined Raman shifts suggesting a highly crystalline nature. Cores of interstitial graphite grains have, on average, lower δ¹³C values by ∼0.5‰ compared to that of the rim. The second type of graphite (Gr_II) occurs as solid inclusions in silicate minerals, commonly forming regular hexagonal crystals with a slightly disordered structure. The third type of graphite (Gr_III) is associated with solid inclusions (up to 100 μm) that have decrepitation halos of numerous small (<15 μm) satellite fluid inclusions of pure CO₂ with varying density (1.105 to 0.75 g/cm³). The fourth type of graphite (Gr_IV) is found as daughter crystals within primary type CO₂-fluid inclusions in garnet and quartz. These fluid inclusions have a range of densities (1.05 to 0.90 g/cm³), but in general are significantly less dense than graphite-free primary, pure CO₂ fluid inclusions (1.12 g/cm³). Raman spectral characteristics of graphite inside fluid inclusions suggest graphite crystallization at low temperature (∼ 500°C). The precipitation of graphite probably occurred during the isobaric cooling of CO₂-rich peak metamorphic fluid as a result of oxyexsolution of oxide phases. The oxyexsolution process is evidenced by the magnetite-ilmenite granular exsolution textures and the systematic presence of numerous micron-sized rutile and other oxide inclusions in association with fluid inclusions within garnet, plagioclase, and quartz.The carbon isotope compositions of coexisting CO₂ (in fluid inclusions) and graphite show a fractionation (α₂CO−gr) of ∼6‰ in garnet, consistent with the existing theoretical estimates of α₂CO−gr at 800°C. A subsequent generation of CO₂ inclusions trapped in matrix quartz and quartz segregation have higher δ¹³C values, −4‰ and −2.9‰ respectively. Graphite in quartz segregations also has higher δ¹³C values (−9.8‰) than those in enderbite (−12.7‰). Micro-graphite crystals included in garnet, quartz (enderbite), and quartz (segregation) have average δ¹³C values of −11.1, −10.4, and −8.7‰ respectively, indicating progressive enrichment in ¹³C with a decrease in temperature of recrystallization of respective minerals. This progressive enrichment is also observed in carbon isotope compositions of fluid inclusion CO₂, suggesting isotopic equilibrium during graphite precipitation from CO₂ fluids. Thus, the carbon isotope record preserved in these rocks by the interstitial graphite, CO₂ fluid in enderbite, graphite microcrystals, graphite in quartz segregation, and CO₂ fluid in quartz segregation, suggests a temperature-controlled isotopic evolution. This evolution is in accordance with a closed system Rayleigh-type graphite precipitation process which progressively enriched residual CO₂ in ¹³C.     

Evidence for a nitrogen flux directly derived from the European subcontinental mantle in the Western Eger Rift, Central Europe

In the Czech-German border region of the Vogtland and NW Bohemia (western Eger rift, Central Europe), chemical and isotopic compositions (C, N, He, Ar) of free gas from a thermal water escape (fluorite mine, Schönbrunn), two mineral springs (“Eisenquelle,” Bad Brambach; “Sprudel III,” Bad Elster) and a mofette (Bublak) located along an ∼40-km long traverse are reported. The gases of Bublak and Bad Brambach are CO₂-rich (>99 vol.%) and have δ¹³C values of −1.95 and −4.29‰, respectively. With distance from the center of CO₂ degassing (Bublak) the δ¹³C values decrease, most likely due to physico-chemical fractionation of CO₂ between gaseous and aqueous phases rather than to admixture of organic/biogenic CO₂. The δ¹⁵N values range between −3.2 and −0.6‰, compared to an upper mantle value of −4.0 ± 1.0‰. The four locations are characterized by ³He/⁴He ratios decreasing from 5.9 R_a in the center (Bublak) to 0.8 R_a in the periphery (Schönbrunn) and give evidence for mixing of He from a deep-seated magmatic source with a crustal source. The location with the highest ³He/⁴He ratio (5.9 R_a) is accompanied by the highest ⁴⁰Ar/³⁶Ar (550). We argue that the nitrogen of the Bublak mofette gas is a mixture of predominantly atmospheric and mantle-derived components, whereas at the other three locations crustal nitrogen may also be present. The Bublak δ¹⁵N value of ≈−4.5 ± 1.0‰ represents the first free gas δ¹⁵N reference from the European subcontinental mantle (ESCM) and indicates that, in contrast to the ³He/⁴He ratios, the δ¹⁵N values are equal for ESCM and MORB, respectively.     

Incremental vacuum dehydration-decarbonation experiments on a natural gibbsite (α-Al(OH3)): CO2 abundance and δC values

Incremental vacuum dehydration-decarbonation experiments were performed at 190°C on chemically “cleaned” aliquots of a gibbsite-dominated, Eocene-age bauxite sample with evolution of CO₂ and H₂O. “Plateau” F (CO₂/H₂O ratios) and δ¹³C values of the CO₂ derived from gibbsite were attained over the dehydration interval, X_v(H₂) = 0.16 to 0.67 (i.e., 16 to 67% breakdown of gibbsite). The plateau value of F for gibbsite was 0.0043 ± 0.0003, while the corresponding δ¹³C value of evolved CO₂ was −16.0‰±0.4‰. Additional experiments on chemically cleaned aliquots included (1) treatment with a solution of 0.3M Na-Citrate + 0.1M Na-Dithionite and (2) an exchange experiment with 0.1 bar of ¹³C-depleted CO₂ (−46‰) at 105°C for 64.5 h. Neither of these additional treatments resulted in a measurable perturbation of plateau values of F or δ¹³C for CO₂ evolved from gibbsite during dehydroxylation. These results support published work on Holocene samples which suggested that CO₂ occluded in gibbsite may preserve information on δ¹³C values of CO₂ in ancient terrestrial systems. The plateau values of F observed in the Eocene gibbsite indicate that it may be possible to experimentally calibrate a relationship between the concentration of CO₂ occluded in gibbsite and CO₂ in the environment at the time of crystallization. Such a calibration would significantly enhance the value of gibbsite as a source of information on ancient oxidized carbon systems.     

Chemical and isotopic equilibrium between CO2 and CH4 in fumarolic gas discharges: Generation of CH4 in arc magmatic-hydrothermal systems

The chemical and isotopic composition of fumarolic gases emitted from Nisyros Volcano, Greece, and of a single gas sample from Vesuvio, Italy, was investigated in order to determine the origin of methane (CH₄) within two subduction-related magmatic-hydrothermal environments.Apparent temperatures derived from carbon isotope partitioning between CH₄ and CO₂ of around 340°C for Nisyros and 470°C for Vesuvio correlate well with aquifer temperatures as measured directly and/or inferred from compositional data using the H₂O-H₂-CO₂-CO-CH₄ geothermometer. Thermodynamic modeling reveals chemical equilibrium between CH₄, CO₂ and H₂O implying that carbon isotope partitioning between CO₂ and CH₄ in both systems is controlled by aquifer temperature.N₂/³He and CH₄/³He ratios of Nisyros fumarolic gases are unusually low for subduction zone gases and correspond to those of midoceanic ridge environments. Accordingly, CH₄ may have been primarily generated through the reduction of CO₂ by H₂ in the absence of any organic matter following a Fischer-Tropsch-type reaction. However, primary occurrence of minor amounts of thermogenic CH₄ and subsequent re-equilibration with co-existing CO₂ cannot be ruled out entirely. CO₂/³He ratios and δ¹³C_CO2 values imply that the evolved CO₂ either derives from a metasomatized mantle or is a mixture between two components, one outgassing from an unaltered mantle and the other released by thermal breakdown of marine carbonates. The latter may contain traces of organic matter possibly decomposing to CH₄ during thermometamorphism.