Meteoric fluid‐rock interaction in Variscan shear zones

Variscan shear zones in the Armorican Massif represent sites of strong fluid‐rock interaction. The hydrogen isotope composition of muscovite (δD_Ms) from syntectonic leucogranite allows to determine the source of fluids that infiltrated the footwall of three detachment zones and the South Armorican Shear Zone. Using temperatures of hydrogen isotope exchange estimated from microstructural data, we calculate the hydrogen isotope ratios of water (δD_water) present within the shear zones during high‐temperature deformation. A ~40‰ difference in δD_water values from deep to shallow crustal level reveals a mixing relationship between deep crustal fluids with higher δD values that range from ?34 to ?33‰, and meteoric fluids with δD values as low as ?74‰ in the upper part of detachment footwalls.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

Genesis of fumarolic emissions as inferred by isotope mass balances: CO2 and water at Vulcano Island, Italy

We have developed a quantitative model of CO₂ and H₂O isotopic mixing between magmatic and hydrothermal gases for the fumarolic emissions of the La Fossa crater (Vulcano Island, Italy). On the basis of isotope balance equations, the model takes into account the isotope equilibrium between H₂O and CO₂ and extends the recent model of chemical and energy two-end-member mixing by Nuccio et al. (1999). As a result, the H₂O and CO₂ content and the δD, δ¹⁸O, and δ¹³C isotope compositions for both magmatic and hydrothermal end-members have been assessed. Low contributions of meteoric steam, added at a shallow depth, have been also recognized and quantified in the fumaroles throughout the period from 1988 to 1998. Nonequilibrium oxygen isotope exchange also seems to be occurring between ascending gases and wall rocks along some fumarolic conduits.The δ¹³C_CO2 of the magmatic gases varies around −3 to 1‰ vs. Peedee belemnite (PDB), following a perfect synchronism with the variations of the CO₂ concentration in the magmatic gases. This suggests a process of isotope fractionation because of vapor exsolution caused by magma depressurization. The hydrogen isotopes in the magmatic gases (−1 to −‰ vs. standard mean ocean water [SMOW]), as well as the above δ¹³C_CO2 value, are coherent with a convergent tectonic setting of magma generation, where the local mantle is widely contaminated by fluids released from the subducted slab. Magma contamination in the crust probably amplifies this effect.The computed isotope composition of carbon and hydrogen in the hydrothermal vapors has been used to calculate the δD and δ¹³C of the entire hydrothermal system, including mixed H₂O-CO₂ vapor, liquid water, and dissolved carbon. We have computed values of about 10‰ vs. SMOW for water and −2 to −6.5‰ vs. PDB for CO₂. On these grounds, we think that Mediterranean marine water (δD_H2O ≈ 10‰) feeds the hydrothermal system. It infiltrates at depth throughout the local rocks, reaching oxygen isotope equilibrium at high temperatures. Interaction processes between magmatic gases and the evolving seawater also seem to occur, causing the dissolution of isotopically fractionated aqueous CO₂ and providing the source for hydrothermal carbon. These results have important implications concerning fluid circulation beneath Vulcano and address the more convenient routine of geochemical surveillance.     

O/O and D/H ratios of pedogenic kaolinite in a North American Cenomanian laterite: Paleoclimatic implications

Kaolinite, gibbsite and quartz are the dominant minerals in samples collected from two outcrops of a Cenomanian (∼95 Ma) laterite in southwestern Minnesota. A combination of measured yields and isotope ratios permitted mass balance calculations of the δD and δ¹⁸O values of the kaolinite in these samples. These calculations yielded kaolinite δD values of about −73‰ and δ¹⁸O values of about +18.7‰. The δD and δ¹⁸O values appear to preserve information on the ancient weathering system.If formed in hydrogen and oxygen isotope equilibrium with water characterized by the global meteoric water line (GMWL), the kaolinite δD and δ¹⁸O values indicate a crystallization temperature of 22 (±5) °C. A nominal paleotemperature of 22 °C implies a δ¹⁸O value for the corresponding water of −6.3‰. The combination of temperature and meteoric water δ¹⁸O values is consistent with relatively intense rainfall at that mid-paleolatitude location (∼40°N) on the eastern shore of the North American Western Interior Seaway. The inferred Cenomanian paleosol temperature of ∼22 °C is in general accord with published mid-Cretaceous continental mean annual temperatures (MAT) estimated from leaf margin analyses of fossil plants.When compared with results from a published GCM-based Cenomanian climate simulation which specifies a latitudinal sea surface temperature (SST) gradient that was either near modern or smaller-than-modern, the kaolinite paleotemperature of 22 °C is closer to the GCM-predicted MAT for a smaller equator-to-pole temperature difference in the mid-Cretaceous. Moreover, the warm, kaolinite-derived, mid-paleolatitude temperature of 22 °C is associated with proxy estimates of high concentrations of atmospheric CO₂ in the Cenomanian. The overall similarity of proxy and model results suggests that the general features of Cenomanian continental climate in that North American locale are probably being revealed.     

Stable isotope fractionation during bacterial sulfate reduction is controlled by reoxidation of intermediates

Bacterial sulfate reduction is one of the most important respiration processes in anoxic habitats and is often assessed by analyzing the results of stable isotope fractionation. However, stable isotope fractionation is supposed to be influenced by the reduction rate and other parameters, such as temperature. We studied here the mechanistic basics of observed differences in stable isotope fractionation during bacterial sulfate reduction. Batch experiments with four sulfate-reducing strains (Desulfovibrio desulfuricans, Desulfobacca acetoxidans, Desulfonatronovibrio hydrogenovorans, and strain TRM1) were performed. These microorganisms metabolize different carbon sources (lactate, acetate, formate, and toluene) and showed broad variations in their sulfur isotope enrichment factors. We performed a series of experiments on isotope exchange of ¹⁸O between residual sulfate and ambient water. Batch experiments were conducted with ¹⁸O-enriched (δ¹⁸O_water = +700‰) and depleted water (δ¹⁸O_water = −40‰), respectively, and the stable ¹⁸O isotope shift in the residual sulfate was followed. For Desulfovibrio desulfuricans and Desulfonatronovibrio hydrogenovorans, which are both characterized by low sulfur isotope fractionation (ε_S > −13.2‰), δ¹⁸O values in the remaining sulfate increased by only 50‰ during growth when ¹⁸O-enriched water was used for the growth medium. In contrast, with Desulfobacca acetoxidans and strain TRM1 (ε_S < −22.7‰) the residual sulfate showed an increase of the sulfate δ¹⁸O close to the values of the enriched water of +700‰. In the experiments with δ¹⁸O-depleted water, the oxygen isotope values in the residual sulfate stayed fairly constant for strains Desulfovibrio desulfuricans, Desulfobacca acetoxidans and Desulfonatronovibrio hydrogenovorans. However, strain TRM1, which exhibits the lowest sulfur isotope fractionation factor (ε_S < −38.7‰) showed slightly decreasing δ¹⁸O values.Our results give strong evidence that the oxygen atoms of sulfate exchange with water during sulfate reduction. However, this neither takes place in the sulfate itself nor during formation of APS (adenosine-5′-phosphosulfate), but rather in intermediates of the sulfate reduction pathway. These may in turn be partially reoxidized to form sulfate. This reoxidation leads to an incorporation of oxygen from water into the “recycled” sulfate changing the overall ¹⁸O isotopic composition of the remaining sulfate fraction. Our study shows that such incorporation of ¹⁸O is correlated with the stable isotope enrichment factor for sulfur measured during sulfate reduction. The reoxidation of intermediates of the sulfate reduction pathway does also strongly influence the sulfur stable isotope enrichment factor. This aforesaid reoxidation is probably dependent on the metabolic conversion of the substrate and therefore also influences the stable isotope fractionation factor indirectly in a rate dependent manner. However, this effect is only indirect. The sulfur isotope enrichment factors for the kinetic reactions themselves are probably not rate dependent.     

Liquid-vapor fractionation of boron and boron isotopes: Experimental calibration at 400°C/23 MPa to 450°C/42 MPa

We experimentally determined the boron partitioning and boron isotope fractionation between coexisting liquid and vapor in the system H₂O−NaCl−B₂O₃. Experiments were performed along the 400 and 450°C isotherms. Pressure conditions ranged from 23 to 28 MPa at 400°C and from 38 to 42 MPa at 450°C. Boron partitions preferentially into the liquid. Its overall liquid-vapor fractionation is, however, weak: Calculated boron distribution coefficients D_B^liquid-vapor are < 2.5 at all run conditions. With decreasing pressure (i.e. increasing opening of the solvus) D_B^liquid-vapor increases along the individual isotherms. Extrapolation to salt saturated conditions yields maximum boron liquid-vapor fractionations of D_B^liquid-vapor = 1.8 at 450°C and D_B^liquid-vapor = 2.7 at 400°C. ¹¹B preferentially fractionates into the vapor. Calculated Δ¹¹B_vapor-liquid = {[(¹¹B/¹⁰B)_vapor - (¹¹B/¹⁰B)_liquid]/(¹¹B/¹⁰B)_{NBS 951}}*1000 are small and range from 0.2 (± 0.7) to 0.9 (± 0.5) ‰ at 450°C and from 0.1 (± 0.6) to 0.7 (± 0.6) ‰ at 400°C. The data indicate increasing isotopic fractionation with decreasing pressure (i.e. increasing opening of the solvus). Extrapolation to salt saturated conditions yields maximum boron isotope liquid-vapor fractionations of Δ¹¹B_vapor-liquid = 1.5 (± 0.7) ‰ at 450°C and Δ¹¹B_vapor-liquid = 1.3 (± 0.6) ‰ at 400°C. The weak boron isotope fractionation suggests similar trigonal speciation in liquid and vapor. Although the boron and boron isotope fractionation between liquid and vapor is only weak, mass balance calculations indicate that for high degrees of fractionation liquid-vapor phase separation in an open system can significantly alter the boron and boron isotope signature of low-salinity hydrous fluids in hydrothermal systems. Comparing the model calculations with natural oceanic hydrothermal fluids, however, indicate that other processes than fluid phase separation dominate the boron geochemistry in oceanic hydrothermal fluids.     

Hydrogen isotope ratios of recent lacustrine sedimentary n-alkanes record modern climate variability

Hydrogen isotope ratios were measured on n-alkanes (n-C₁₂ to n-C₃₁) extracted from recent lake surface sediments along a N-S European transect to test if modern climate variability is recorded in these biomarkers. δD values of the n-alkanes are compared to δD values of meteoric water from the IAEA-GNIP database spanning a range from −119‰ in northern Sweden to −41‰ in southern Italy, to lake water δD values, and to mean annual temperatures, varying between −2.0°C in the north and 13.7°C in the south.δD values of the short-chained n-alkanes n-C₁₂ to n-C₂₀, excluding algal derived n-C₁₇ and n-C₁₉, are higher in the north and lower in the south. The isotopic fractionation ε for hydrogen between meteoric water and the short-chained n-alkanes is increasing from N to S by more than 100‰ and is significantly correlated to mean annual temperature for n-C₁₆ and n-C₁₈. This suggests that these n-alkanes may originate from a different source in the northern lakes, possibly due to petroleum contamination, or are synthesized using a different biochemical pathway.The n-C₁₇ and n-C₁₉ alkanes of algal origin, the n-C₂₁ and n-C₂₃ alkanes originating from water plants, and the long-chain n-alkanes n-C₂₅, n-C₂₇, n-C₂₉, and n-C₃₁ of terrestrial origin, clearly correlate with δD values of meteoric water, lake water, and mean annual temperature, indicating that they excellently record the δD value of meteoric water. The mean hydrogen isotope fractionation ε_C17/w of −157‰ (SD = 13) between n-C₁₇ and meteoric water is fairly constant over the wide range of different climates and lake environments, suggesting only minor influence of environmental factors on this biochemical fractionation. This suggests that δD values of n-C₁₇ are suitable to reconstruct the isotopic composition of source water. The mean fractionation between the long-chain n-alkanes and water is −128‰ (SD = 12). The mean difference of 31‰ between both ε values is likely due to evaporative enrichment of deuterium in the leaf water. If this is the only influence on the enrichment, the difference between the δD values of terrestrial and aquatic compounds might be suitable to reconstruct terrestrial evapotranspiration of the lake environment.     

Water drives the deuterium content of the methane emitted from plants

The spatial distribution of the deuterium content of precipitation has a well-established latitudinal variation that is reflected in organic molecules in plants growing at different locations. Some laboratory and field studies have already shown that the deuterium content of methane emitted from methanogens can be partially related to δD variations of the water in the surrounding environment. Here we present a similar relation for the methane emitted from plant biomass under UV radiation. To show this relation, we determined the hydrogen isotopic composition of methane released from leaves of a range of plants grown with water of different deuterium content (δD = −130‰ to +115‰). The plant leaves were irradiated with UV light and the CH₄ isotopic composition was measured by continuous flow isotope ratio mass spectrometry (CF-IRMS). Furthermore, the deuterium content of bulk biomass and of the methoxyl (OCH₃) groups of the biomass was measured. The D/H ratio successively decreases from bulk biomass (δD = −106‰ to −50‰) via methoxyl groups (δD = −310‰ to −115‰) to the CH₄ emitted (δD = −581‰ to −196‰). The range of isotope ratios in bulk biomass and OCH₃ groups is smaller than in the water used to grow the plants. Methoxyl groups, which contain only non-exchangeable hydrogen, can be used to assess the fraction of external water that was incorporated before OCH₃ groups were formed. Surprisingly, the CH₄ formed under UV irradiation has a wider isotopic range than the OCH₃ groups. Although the precise production pathway cannot be fully determined, the presented experiments indicate that methoxyl groups are not the only source substrate for CH₄, but that other sources, including very depleted ones, must contribute. The main limitation to the interpretation of the data is the possible influence of exchangeable water, which could not be quantified. Future studies should include measurements of leaf water and avoid interaction between different plants via the gas phase. Despite these deficiencies, the results suggest that the deuterium content of the methane generated from plants under UV irradiation is closely linked to δD in precipitation. This dependency, which should also exist for other biogenic methane sources could be evaluated with global isotope models.     

Pressure effects on the reduced partition function ratio for hydrogen isotopes in water

We have developed a simple, yet accurate theoretical method for calculating the reduced isotope partition function ratio (RIPFR) for hydrogen of water at elevated pressures. This approach requires only accurate equations of state (EOS) for pure isotopic end-members (H₂O and D₂O), which are available in the literature. The effect of pressure or density on the RIPFR of water was calculated relative to that of ideal-gas water at infinitely low pressure for the temperature range from 0 to 527 °C. For gaseous and low-pressure (ca. ?15 MPa) supercritical phases of water, the RIPFR increases slightly (1-1.3‰) with pressure or density in a fashion similar to those of many other geologic materials. However, in liquid and high-pressure (>20 MPa) supercritical phases, the RIPFR of water decreases (0.5-6‰) with increasing pressure (or density) to 100 MPa. This rather unique phenomenon is ascribed to the inverse molar volume isotope effects (MVIE) of liquid and high-density supercritical waters, V (D₂O) > V (H₂O), while other substances including minerals show the normal MVIE. These theoretical predictions were experimentally confirmed by Horita et al. [Horita, J., Cole, D.R., Polyakov, V.B., Driesner, T., 2002. Experimental and theoretical study of pressure effects on hydrogen isotope fractionation in the system brucite-water at elevated temperatures. Geochim. Cosmochim. Acta66, 3769 - 3788.] for the system brucite-water. Although the P-T ranges for the EOS of normal and heavy waters are rather limited, our modeling indicates that the RIPFR of water continues to decrease with pressure above 100 MPa. The method developed here can be applied to any other geologic fluids, if accurate EOS for their isotopic end-members is available. These results have important implications for the interpretation of high-pressure isotopic partitioning in the Earth, the outer planets, and their moons.     

Isotopic constraints on ice age fluids in active geothermal systems: Reykjanes, Iceland

The Reykjanes geothermal system is located on the landward extension of the Mid-Atlantic Ridge in southwest Iceland, and provides an on-land proxy to high-temperature hydrothermal systems of oceanic spreading centers. Previous studies of elemental composition and salinity have shown that Reykjanes geothermal fluids are likely hydrothermally modified seawater. However, δD values of these fluids are as low as −23‰, which is indicative of a meteoric water component. Here we constrain the origin of Reykjanes hydrothermal solutions by analysis of hydrogen and oxygen isotope compositions of hydrothermal epidote from geothermal drillholes at depths between 1 and 3 km. δD_EPIDOTE values from wells RN-8, -9, -10 and -17 collectively range from −60 to −78‰, and δ¹⁸O_EPIDOTE in these wells are between −3.0 and 2.3‰. The δD values of epidote generally increase along a NE trend through the geothermal field, whereas δ¹⁸O values generally decrease, suggesting a southwest to northeast migration of the geothermal upflow zone with time that is consistent with present-day temperatures and observed hydrothermal mineral zones. For comparative analysis, the meteoric-water dominated Nesjavellir and Krafla geothermal systems, which have a δD_FLUID of ∼ −79‰ and −89‰, respectively, show δD_EPIDOTE values of −115‰ and −125‰. In contrast, δD_EPIDOTE from the mixed meteoric-seawater Svartsengi geothermal system is −68‰; comparable to δD_EPIDOTE from well RN-10 at Reykjanes.Stable isotope compositions of geothermal fluids in isotopic equilibrium with the epidotes at Reykjanes are computed using published temperature dependent hydrogen and oxygen isotope fractionation curves for epidote-water, measured isotope composition of the epidotes and temperatures approximated from the boiling point curve with depth. Calculated δD and δ¹⁸O of geothermal fluids are less than 0‰, suggesting that fluids of meteoric or glacial origin are a significant component of the geothermal solutions. Additionally, δD_FLUID values in equilibrium with geothermal epidote are lower than those of modern-day fluids, whereas calculated δ¹⁸O_FLUID values are within range of the observed fluid isotope composition. We propose that modern δD_EPIDOTE and δD_FLUID values are the result of diffusional exchange between hydrous alteration minerals that precipitated from glacially-derived fluids early in the evolution of the Reykjanes system and modern seawater-derived geothermal fluids. A simplified model of isotope exchange in the Reykjanes geothermal system, in which the average starting δD_ROCK value is −125‰ and the water to rock mass ratio is 0.25, predicts a δD_FLUID composition within 1‰ of average measured values. This model resolves the discrepancy between fluid salinity and isotope composition of Reykjanes geothermal fluids, explains the observed disequilibrium between modern fluids and hydrothermal epidote, and suggests that rock-fluid interaction is the dominant control over the evolution of fluid isotope composition in the hydrothermal system.     

The origin and mechanisms of salinization of the lower Jordan river

The chemical and isotopic (⁸⁷Sr/⁸⁶Sr, δ¹¹B, δ³⁴S_sulfate, δ¹⁸O_water, δ¹⁵N_nitrate) compositions of water from the Lower Jordan River and its major tributaries between the Sea of Galilee and the Dead Sea were determined in order to reveal the origin of the salinity of the Jordan River. We identified three separate hydrological zones along the flow of the river:

(1)

A northern section (20 km downstream of its source) where the base flow composed of diverted saline and wastewaters is modified due to discharge of shallow sulfate-rich groundwater, characterized by low ⁸⁷Sr/⁸⁶Sr (0.7072), δ³⁴S_sulfate (−2‰), high δ¹¹B (∼36‰), δ¹⁵N_nitrate (∼15‰) and high δ¹⁸O_water (−2 to-3‰) values. The shallow groundwater is derived from agricultural drainage water mixed with natural saline groundwater and discharges to both the Jordan and Yarmouk rivers. The contribution of the groundwater component in the Jordan River flow, deduced from mixing relationships of solutes and strontium isotopes, varies from 20 to 50% of the total flow.

(2)

A central zone (20-50 km downstream from its source) where salt variations are minimal and the rise of ⁸⁷Sr/⁸⁶Sr and SO₄/Cl ratios reflects predominance of eastern surface water flows.

(3)

A southern section (50-100 km downstream of its source) where the total dissolved solids of the Jordan River increase, particularly during the spring (70-80 km) and summer (80-100 km) to values as high as 11.1 g/L. Variations in the chemical and isotopic compositions of river water along the southern section suggest that the Zarqa River (⁸⁷Sr/⁸⁶Sr∼0.70865; δ¹¹B∼25‰) has a negligible affect on the Jordan River. Instead, the river quality is influenced primarily by groundwater discharge composed of sulfate-rich saline groundwater (Cl^-=31-180 mM; SO₄/Cl∼0.2-0.5; Br/Cl∼2-3×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.70805; δ¹¹B∼30‰; δ¹⁵N_nitrate ∼17‰, δ³⁴S_sulfate=4-10‰), and Ca-chloride Rift valley brines (Cl^-=846-1500 mM; Br/Cl∼6-8×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.7080; δ¹¹B>40‰; δ³⁴S_sulfate=4-10‰). Mixing calculations indicate that the groundwater discharged to the river is composed of varying proportions of brines and sulfate-rich saline groundwater. Solute mass balance calculations point to a ∼10% contribution of saline groundwater (Cl⁻=282 to 564 mM) to the river. A high nitrate level (up to 2.5 mM) in the groundwater suggests that drainage of wastewater derived irrigation water is an important source for the groundwater. This irrigation water appears to leach Pleistocene sediments of the Jordan Valley resulting in elevated sulfate contents and altered strontium and boron isotopic compositions of the groundwater that in turn impacts the water quality of the lower Jordan River.

     

Paleogene paleoclimate reconstruction using oxygen isotopes from land and freshwater organisms: the use of multiple paleoproxies

Understanding past climate change is critical to the interpretation of earth history. Even though relative temperature change has been readily assessed in the marine record, it has been more difficult in the terrestrial record due to restricted taxonomic distribution and isotopic fractionation. This problem could be overcome by the use of multiple paleoproxies. Therefore, the δ¹⁸O isotopic composition of five paleoproxies (rodent tooth enamel, δ¹⁸O_Phosphate = +17.7 ± 2.0‰ n = 74 (VSMOW); fish scale ganoine δ¹⁸O_Phosphate = +19.7 ± 0.7‰ n = 20 (VSMOW); gastropod shell δ¹⁸O_Calcite = −1.7 ± 1.3‰ n = 50 (VPDB); charophyte gyrogonite δ¹⁸O_Calcite = −2.4 ± 0.5‰ n = 20 (VPDB); fish otolith δ¹⁸O_Aragonite = δ¹⁸O = −3.6 ± 0.6‰ n = 20 (VPDB)) from the Late Eocene (Priabonian) Osborne Member (Headon Hill Formation, Solent Group, Hampshire Basin, UK) were determined. Because diagenetic alteration was shown to be minimal the phosphate oxygen component of rodent tooth enamel (as opposed to enamel carbonate oxygen) was used to calculate an initial δ¹⁸O_{Local water} value of 0.0 ± 3.4‰. However, a skewed distribution, most likely as a result of the ingestion of evaporating water, necessitated the calculation of a corrected δ¹⁸O_{Local water} value of −1.3 ± 1.7‰ (n = 62). This δ¹⁸O_{Local water} value corresponds to an approximate mean annual temperature of 18 ± 1°C. Four other mean paleotemperatures can also be calculated by combining the δ¹⁸O_{Local water} value with four independent freshwater paleoproxies. The calculated paleotemperature using the fish scale thermometry equations most likely represents the mean temperature (21 ± 2°C) of the entire length of the growing season. This should be concordant with the paleotemperature calculated using the Lymnaea shell thermometry equation (23 ± 2°C). The lack of concordance is interpreted to be the result of diagenetic alteration of the originally aragonitic Lymnaea shell to calcite. The mean paleotemperature calculated using the charophyte gyrogonite thermometry equation (21 ± 2°C), on the other hand, most likely represents the mean temperature of a single month toward the end of the growing season. The fish otolith mean paleotemperature (28 ± 2°C) most likely represents the mean temperature of the warmest months of the growing season. An approximate mean annual temperature of 18 ± 1°C, in addition to a mean growing season paleotemperature of 21 ± 2°C (using fish scale only) with a warmest month temperature of 28 ± 2°C, and high associated standard deviations suggest that a subtropical to warm temperate seasonal climate existed during the deposition of the Late Eocene Osborne Member.     

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

Comment on “A critical evaluation of the boron isotope-pH proxy: The accuracy of ancient ocean pH estimates” by M. Pagani, D. Lemarchand, A. Spivack and J. Gaillardet

Pagani et al. [Pagani M., Lemarchand D., Spivack A., and Gaillardet J. (2005). A critical evaluation of the boron isotope-pH proxy: the accuracy of ancient ocean pH estimates. Geochim. Cosmochim. Acta69(4), 953-961] use data from previous boron isotope studies to suggest that the fractionation between boric acid and borate in seawater as well as the history of δ¹¹B in seawater are poorly understood, thus limiting our ability to capture realistic ocean pH with this proxy. Although we agree with the authors that the long recognized uncertainty in the secular variation of δ¹¹B_seawater imposes a temporal limit on paleo-pH reconstructions, their evaluation of the δ¹¹B/pH relationship in carbonates is flawed. Potential complications from vital, temperature and dissolution effects reported in that paper are based on studies that are experimentally and/or analytically poorly constrained. Using published validation studies we will demonstrate that many of the problems outlined by Pagani et al. have already been addressed, or are based on misinterpretations of previous work. Most importantly, statistical evaluation suggests empirical data are best described by a fractionation of ∼20‰. Recent paleoreconstructions confirm that the boron isotope proxy can be used with confidence, if sample selection and analyses are done carefully.     

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.     

Natural fractionation of U/U

The isotopic composition of U in nature is generally assumed to be invariant. Here, we report variations of the ²³⁸U/²³⁵U isotope ratio in natural samples (basalts, granites, seawater, corals, black shales, suboxic sediments, ferromanganese crusts/nodules and BIFs) of ∼1.3‰, exceeding by far the analytical precision of our method (≈0.06‰, 2SD). U isotopes were analyzed with MC-ICP-MS using a mixed ²³⁶U-²³³U isotopic tracer (double spike) to correct for isotope fractionation during sample purification and instrumental mass bias. The largest isotope variations found in our survey are between oxidized and reduced depositional environments, with seawater and suboxic sediments falling in between. Light U isotope compositions (relative to SRM-950a) were observed for manganese crusts from the Atlantic and Pacific oceans, which display δ²³⁸U of −0.54‰ to −0.62‰ and for three of four analyzed Banded Iron Formations, which have δ²³⁸U of −0.89‰, −0.72‰ and −0.70‰, respectively. High δ²³⁸U values are observed for black shales from the Black Sea (unit-I and unit-II) and three Kupferschiefer samples (Germany), which display δ²³⁸U of −0.06‰ to +0.43‰. Also, suboxic sediments have slightly elevated δ²³⁸U (−0.41‰ to −0.16‰) compared to seawater, which has δ²³⁸U of −0.41 ± 0.03‰. Granites define a range of δ²³⁸U between −0.20‰ and −0.46‰, but all analyzed basalts are identical within uncertainties and slightly lighter than seawater (δ²³⁸U = −0.29‰).Our findings imply that U isotope fractionation occurs in both oxic (manganese crusts) and suboxic to euxinic environments with opposite directions. In the first case, we hypothesize that this fractionation results from adsorption of U to ferromanganese oxides, as is the case for Mo and possibly Tl isotopes. In the second case, reduction of soluble U^VI to insoluble U^IV probably results in fractionation toward heavy U isotope compositions relative to seawater. These findings imply that variable ocean redox conditions through geological time should result in variations of the seawater U isotope compositions, which may be recorded in sediments or fossils. Thus, U isotopes might be a promising novel geochemical tracer for paleo-redox conditions and the redox evolution on Earth. The discovery that ²³⁸U/²³⁵U varies in nature also has implications for the precision and accuracy of U-Pb dating. The total observed range in U isotope compositions would produce variations in ²⁰⁷Pb/²⁰⁶Pb ages of young U-bearing minerals of up to 3 Ma, and up to 2 Ma for minerals that are 3 billion years old.     

Silicon isotopic composition of dissolved silicon and suspended particulate matter in the Yellow River, China, with implications for the global silicon cycle

The silicon isotopic composition of dissolved silicon and suspended particulate matter (SPM) were systematically investigated in water samples from the mainstem of the Yellow River and 4 major tributaries. The SPM content of the Yellow River varied from 1.4 to 38,560 mg/L, averaging 3568 mg/L, and the δ³⁰Si of suspended particulate matter (δ³⁰Si_SPM) varied from 0.3‰ to −0.4‰, averaging −0.02‰. The major factors affecting the SPM content and the δ³⁰Si_SPM values in the Yellow River were inferred to be the mineralogical, chemical and isotopic characteristics of the sediments from the Loess Plateau and a combination of the climate and the flow discharge of the river.The major ions in the Yellow River water were Na⁺, Ca²⁺, Mg²⁺, HCO₃⁻, SO₄²⁻ and Cl⁻. High salt concentration was observed in samples from the middle and lower reaches, likely reflecting the effects of evaporation and irrigation because the Na⁺, Mg²⁺, SO₄²⁻ and K⁺ concentrations were correlated with the Cl⁻ concentration. The dissolved Si concentration (D_Si) increased downstream, varying from 0.016 to 0.323 mM. The δ³⁰Si of dissolved Si (δ³⁰Si_Diss) varied from 0.4‰ to 2.5‰, averaging 1.28‰. The major processes controlling the D_Si and δ³⁰Si_Diss of the Yellow River are (a) the weathering of silicate rocks, (b) the formation of phytoliths in plants, (c) the evaporation of water from and the addition of meteoric water to the river system, which only affects concentrations, (d) the adsorption and desorption of aqueous monosilicic acid on iron oxide, and (e) the dissolution of phytoliths in soils.The D_Si and δ³⁰Si_Diss values of global rivers vary spatially and temporally in response to changes in climate, chemical weathering intensity and biological activity. The moderately positive δ³⁰Si_Diss values observed in the Yellow River may be attributed to the higher rates of chemical weathering and biological activities that have been observed in this catchment in comparison with those of other previously studied catchments, excluding the Yangtze River. Human activities may also potentially influence chemical weathering and biological activities and affect the D_Si and δ³⁰Si_Diss values of the major rivers of the world. Further river studies should be performed to gain a better understanding of the global Si isotope budget.     

Changing moisture sources over the last 330,000 years in Northern Oman from fluid-inclusion evidence in speleothems

Speleothems from Hoti Cave in northern Oman provide a record of continental pluvial periods over the last 330,000 yr. Periods of rapid speleothem deposition occurred from 6000 to 10,500, 78,000 to 82,000, 120,000 to 135,000, 180,000 to 200,000, and 300,000 to 330,000 yr ago, with little or no growth during the intervening periods. During each of these five pluvial periods, δD values of water extracted from speleothem fluid inclusions (δD_FI) are between −60 and −20‰ (VSMOW) and δ¹⁸O values of speleothem calcite (δ¹⁸O_C) are between −12 and −4‰ to (VPDB). These values are much more negative than modern rainfall (for δD) or modern stalagmites (for δ¹⁸O). Previous work on the isotopic composition of rainfall in Oman has shown that northern and southern moisture sources are isotopically distinct. Combined measurements of the δD values of fluid-inclusion water with calculated δ¹⁸O values from peak interglacial speleothems indicate that groundwater was predominantly recharged by the southern (Indian Ocean) moisture source, when the monsoon rainfall belt moved northward and reached Northern Oman during each of these periods.     

The effect of early diagenesis on the Fe isotope compositions of porewaters and authigenic minerals in continental margin sediments

Iron isotope compositions in marine pore fluids and sedimentary solid phases were measured at two sites along the California continental margin, where isotope compositions range from δ⁵⁶Fe = −3.0‰ to +0.4‰. At one site near Monterey Canyon off central California, organic matter oxidation likely proceeds through a number of diagenetic pathways that include significant dissimilatory iron reduction (DIR) and bacterial sulfate reduction, whereas at our other site in the Santa Barbara basin DIR appears to be comparatively small, and production of sulfides (FeS and pyrite) was extensive. The largest range in Fe isotope compositions is observed for Fe(II)_aq in porewaters, which generally have the lowest δ⁵⁶Fe values (minimum: −3.0‰) near the sediment surface, and increase with burial depth. δ⁵⁶Fe values for FeS inferred from HCl extractions vary between ∼−0.4‰ and +0.4‰, but pyrite is similar at both stations, where an average δ⁵⁶Fe value of −0.8 ± 0.2‰ was measured. We interpret variations in dissolved Fe isotope compositions to be best explained by open-system behavior that involves extensive recycling of Feflux. This study is the first to examine Fe isotope variations in modern marine sediments, and the results show that Fe isotopes in the various reactive Fe pools undergo isotopic fractionation during early diagenesis. Importantly, processes dominated by sulfide formation produce high-δ⁵⁶Fe values for porewaters, whereas the opposite occurs when Fe(III)-oxides are present and DIR is a major pathway of organic carbon respiration. Because shelf pore fluids may carry a negative δ⁵⁶Fe signature it is possible that the Fe isotope composition of ocean water reflects a significant contribution of shelf-derived iron to the open ocean. Such a signature would be an important means for tracing iron sources to the ocean and water mass circulation.     

A comparison of isotope fractionation of carbon and hydrogen from paddy field rice roots and soil bacterial enrichments during CO2/H2 methanogenesis

To better understand the isotope biogeochemistry of paddy field CH₄, we investigated carbon and hydrogen isotope fractionation during CO₂ reduction by a methanogenic community enriched from California paddy field soil and rice plants. Results from analyses of terminal restriction fragment length polymorphism (T-RFLP) and sequences of the archaeal small-subunit (SSU) rRNA-encoding genes (rDNA) showed a difference in methanogenic community structure between the soil (dominated by Methanobacteriaceae) and roots (dominated by Methanospirillaceae) which was essentially the same for sampling dates 15 and 99 days after flooding (DAF). CO₂/H₂ methanogenesis by these microbial communities produced CH₄ with different isotope ratios and fractionation factors (α factors). The carbon isotope α factors in an open system with a continuous supply of 0.5% H₂ were 1.050 ± 0.002 and 1.057 ± 0.001 for soil and root enrichment cultures at 15 DAF, and 1.052 ± 0.0.002 and 1.059 ± 0.002 for soil and root enrichment cultures at 99 DAF, respectively. These α factors are similar to, but distinct from values previously obtained from cultures of mesophilic methanogens and are larger than calculated values (1.045) for paddy soil. Fractionation of hydrogen isotopes was also studied in a closed system under 80% H₂. The difference in α factors between soil and root enrichment cultures remained clear. The hydrogen isotope fractionations between culture water and the product CH₄ were −327 ± 14‰ and −319 ± 18‰ for soil enrichments, and −389 ± 17‰ and −382 ± 21‰ for root enrichments at 15 DAF and 99 DAF, respectively.