Glacial/interglacial temperature variations in Soreq cave speleothems as recorded by ‘clumped isotope’ thermometry

‘Clumped isotope’ thermometry is based on analyzing mass 47 in CO₂ extracted from carbonates and uses the tracer mass 47 anomaly (Δ₄₇). Δ₄₇ is defined as the deviation of R⁴⁷ from that expected for a random distribution of isotopologues and reflects a temperature dependent preference of ¹³C and ¹⁸O to create a bond with each other in CO₂ or in the carbonate lattice. Being an internal characteristic of the carbonate mineral, it is independent of the isotopic composition of the water in which equilibrium precipitation of the carbonate occurs and can therefore be used to independently determine carbonate growth temperatures. This work provides a first examination of the applicability of ‘clumped isotopes’ thermometry to reconstructing the growth temperatures of speleothems, by examining the glacial/interglacial variations of the Δ₄₇ values of speleothem carbonates from Soreq cave, Israel. The results indicate that the last glacial maximum temperatures were 6-7 °C colder than modern day temperature and a sample at 56 Ky BP was 3 °C colder than the modern. Early Holocene temperatures were slightly above modern day, and late Holocene temperatures were slightly below modern day. These temperature variations are similar to those previously estimated for Eastern Mediterranean sea surface water. Cave water was ¹⁸O depleted in the Holocene compared to modern day (by 0.6-1‰) and 1.1‰ more enriched in the last glacial maximum. Comparison of these cave water δ¹⁸O values with fluid inclusion δD values indicated a late Holocene d-excess value within the range of modern rainfall, implying ∼45% relative humidity. Last glacial maximum and early Holocene d-excess values were significantly lower, suggesting relative humidity of ∼60% and ∼70%, respectively. The temperatures reported in this study were empirically corrected for a non-equilibrium artifact observed in a modern speleothem. The similarity of the temperature variations obtained here to other, independent, records in the region suggests that the Δ₄₇-temperature calibration slope observed in inorganic synthetic calcite and marine organisms may also be applied in speleothems. But the offset observed in modern temperature suggests that the intercept is different so that a separate calibration is needed for accurate absolute temperature reconstruction using speleothem ‘clumped isotopes’. Similar examination of additional caves would be necessary to determine whether such empirical correction can be generally applied or is it a unique characteristic of Soreq cave.     

Stable nitrogen and carbon isotope (δN and δC) variability in shallow tropical Pacific soft coral and black coral taxa and implications for paleoceanographic reconstructions

Soft corals and black corals are useful proxy tools for paleoceanographic reconstructions. However, most work has focused on deep-water taxa and few studies have used these corals as proxy organisms in shallow water (<200 m). To facilitate the use of stable nitrogen and carbon isotope (δ¹⁵N and δ¹³C) records from shallow-water soft coral and black coral taxa for paleoceanographic reconstructions, quantification of the inherent variability in skeletal isotope values between sites, across depth, and among taxa is needed. Here, skeletal δ¹⁵N and δ¹³C values were measured in multiple colonies from eleven genera of soft corals and two genera of black corals from across a depth transect (5-105 m) at two sites in Palau located in the tropical western Pacific Ocean. Overall, no difference in skeletal δ¹⁵N and δ¹³C values between sites was present. Skeletal δ¹⁵N values significantly increased and δ¹³C values decreased with depth. This is consistent with changes in isotope values of suspended particulate organic matter (POM) across the photic zone, suggesting that the primary food source to these corals is suspended POM and that the stable isotopic composition of POM controls the skeletal isotopic composition of these corals. Thus, to compare the isotope records of corals collected across a depth range in the photic zone, first order depth corrections of −0.013‰ m⁻¹ and +0.023‰ m⁻¹ are recommended for δ¹⁵N and δ¹³C, respectively. Average depth-corrected δ¹⁵N values were similar between black corals and soft corals, indicating that corals in these orders feed at a similar trophic level. In contrast, average depth-corrected δ¹³C values of black corals were significantly lower than that of soft corals, potentially resulting from metabolic processes associated with differing skeletal compositions among the orders (i.e., gorgonin vs. chitin based). Thus, a correction of +1.0‰ is recommended for black corals when comparing their δ¹³C-based proxy records to soft corals. After correcting for both the depth and order effects, variability in δ¹⁵N values among corals within each genera was low (standard deviation (SD) of the mean <±0.5‰), with the exception of Acanthorgorgia. The calculated SD of <±0.5‰ provides a first order guideline for the amount of variability that could be expected in a δ¹⁵N record, and suggests that these corals may be useful for δ¹⁵N-based paleoceanographic reconstructions. Variability in δ¹³C values among corals within genera was also low (standard deviation of the mean <±0.5‰) with the exception of Rhipidipathes and Villogorgia. Similar to δ¹⁵N, records from the genera studied here with the exception of Rhipidipathes and Villogorgia may be useful for δ¹³C-based paleoceanographic reconstructions. Overall, using the recommendations developed here, stable isotope records from multiple sites, depths and taxa of these corals can be more rigorously compared.     

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

Clumped isotope thermometry of carbonatites as an indicator of diagenetic alteration

We measure the clumped isotopic signature of carbonatites to assess the integrity of the clumped isotope paleothermometer over long timescales (10⁷-10⁹ years) and the susceptibility of the proxy to closed system re-equilibration of isotopes during burial diagenesis. We find pristine carbonatites that have primary oxygen isotope signatures, along with a Carrara marble standard, do not record clumped isotope signatures lighter than 0.31‰ suggesting atoms of carbon and oxygen freely exchange within the carbonate lattice at temperatures above 250-300 °C. There is no systematic trend in the clumped isotope signature of pristine carbonatites with age, although partial re-equilibration to lower temperatures can occur if a carbonatite has been exposed to burial temperatures for long periods of time. We conclude that the solid-state re-ordering of carbon and oxygen atoms is sufficiently slow to enable the use of clumped isotope paleothermometry on timescales of 10⁸ years, but that diagenetic resetting can still occur, even without bulk recrystallization. In addition to the carbonatite data, an inorganic calibration of the clumped isotope paleothermometer for low temperature carbonates (7.5-77 °C) is presented and highlights the need for further inter-lab standardization.     

Deep-sea coral geochemistry: Implication for the vital effect

Deep-sea corals hold a great potential as a key to important aspects of paleoceanography for at least two reasons, 1) they offer temporal high resolution records of deep-sea environment, because they have growth banding structures, 2) and they are well suited for studying vital effects, because the deep-sea environment does not change over short time scales. However, the relationship between the chemical composition of deep-sea coral skeletons and environmental factors is not well understood. In this study, the chemical composition of deep-sea corals was measured in bulk individuals and along skeletal micro-structures. Among the bulk individuals, δ¹⁸O value and Sr / Ca ratio show a negative but weak correlation with ambient temperature. On the other hand, the Mg / Ca ratio has a positive, weak correlation with the temperature. Large variations were found among samples collected from similar temperature. The variation is up to 3.8‰ for δ¹⁸O, 0.9 mmol/mol for Sr / Ca ratios, and 0.78 mmol/mol for Mg / Ca ratios among samples with ambient average temperature within 1 °C. This variation may be due to a large vital effect. The centers of calcification (COCs), which were formed at high calcification rate, have lower Sr / Ca, U / Ca and higher Mg / Ca ratios than surrounding fasciculi. This chemical distribution supports the model that elemental incorporation depends on calcification rate. This suggests that calcification rate is a very important factor for the chemical composition in deep-sea corals and is one of the most significant mechanisms of the vital effect. Because of this large vital effect, further investigations are essential to use the deep-sea coral as a temperature proxy.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Stable isotopes in deep-sea corals and a new mechanism for “vital effects”

Offsets from isotopic equilibrium in biogenic carbonates have complicated paleoclimate reconstructions for decades. A new archive of climate, deep-sea corals, is used to evaluate the calcification processes, independent of photosynthesis, that contribute to these offsets. Carbon and oxygen stable isotope data from six modern deep-sea corals show strong linear trends between δ¹³C and δ¹⁸O. Slopes of these trends between samples are similar and range between 1.9 to 2.6 for Δδ¹³C/Δδ¹⁸O. Linear trends intersect isotopic equilibrium for δ¹⁸O and are slightly depleted for δ¹³C. Variations in the isotopic ratios are strongly correlated with the density banding structure. Isotopically depleted aragonite is associated with light, quickly precipitating bands, whereas isotopically enriched points correspond to slowly accumulating, less dense aragonite. The densest white band at the trabecular center is furthest from isotopic equilibrium for both carbon and oxygen. Data from this region fall off the linear trend between δ¹⁸O and δ¹³C. This deviation, where δ¹³C remains constant while the δ¹⁸O continues to decrease, does not support “vital effect” mechanisms that call upon kinetic fractionation to explain offsets from isotopic equilibrium. We propose a new mechanism for vital effects in these deep-sea corals that is based on a thermodynamic response to a biologically induced pH gradient in the calcifying region.     

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.     

Magnesium stable isotope fractionation in marine biogenic calcite and aragonite

This survey of magnesium stable isotope compositions in marine biogenic aragonite and calcite includes samples from corals, sclerosponges, benthic porcelaneous and planktonic perforate foraminifera, coccolith oozes, red algae, and an echinoid and brachiopod test. The analyses were carried out using MC-ICP-MS with an external repeatability of ±0.22‰ (2SD for δ²⁶Mg; n = 37), obtained from a coral reference sample (JCp-1).Magnesium isotope fractionation in calcitic corals and sclerosponges agrees with published data for calcitic speleothems with an average Δ²⁶Mg_{calcite-seawater} = −2.6 ± 0.3‰ that appears to be weakly related to temperature. With one exception (Vaceletia spp.), aragonitic corals and sclerosponges also display uniform Mg isotope fractionations relative to seawater with Δ²⁶Mg_{biogenic aragonite-seawater} = −0.9 ± 0.2.Magnesium isotopes in high-Mg calcites from red algae, echinoids and perhaps some porcelaneous foraminifera as well as in all low-Mg calcites (perforate foraminifera, coccoliths and brachiopods) display significant biological influences. For planktonic foraminifera, the Mg isotope data is consistent with the fixation of Mg by organic material under equilibrium conditions, but appears to be inconsistent with Mg removal from vacuoles. Our preferred model, however, suggests that planktonic foraminifera synthesize biomolecules that increase the energetic barrier for Mg incorporation. In this model, the need to remove large quantities of Mg from vacuole solutions is avoided. For the high-Mg calcites from echinoids, the precipitation of amorphous calcium carbonate may be responsible for their weaker Mg isotope fractionation.Disregarding superimposed biological effects, it appears that cation light isotope enrichments in CaCO₃ principally result from a chemical kinetic isotope effect, related to the incorporation of cations at kink sites. In this model, the systematics of cation isotope fractionations in CaCO₃ relate to the activation energy required for cation incorporation, which probably reflects the dehydration of the cation and the crystal surface and bond formation at the incorporation site. This kinetic incorporation model predicts (i) no intrinsic dependence on growth rate, unless significant back reaction upon slow growth reduces the isotope fractionation towards that characteristic for equilibrium isotope partitioning (this may be observed for Ca isotopes in calcites), (ii) a small decrease of isotope fractionation with increasing temperature that may be amplified if higher temperatures promote back reaction and (iii) a sensitivity to changes in the activation barrier caused by additives such as anions or biomolecules or by the initial formation of amorphous CaCO₃.     

Intensified mid-Holocene Asian monsoon recorded in corals from Kikai Island, subtropical northwestern Pacific

Coupled records of Sr/Ca and oxygen isotope ratios (δ¹⁸O) of coral skeletons have been used to produce quantitative estimates of paleo-sea surface temperature (SST) and δ¹⁸O of surface seawater that can in some cases be converted to sea surface salinity (SSS). Two fossil corals from Kikai Island in the subtropical northwestern Pacific, a location affected by East Asian summer and winter monsoons, were analyzed to investigate differences between mid-Holocene and present-day SST and SSS. At 6180 cal yr BP, SSTs were roughly the same as today, both in summer and winter; δ¹⁸O_seawater and SSS values were higher both in summer (+ 0.5‰, +1.1 psu) and in winter (+ 0.2‰, + 0.6 psu) than modern values. At 7010 cal yr BP, SSTs were slightly cooler both in summer and winter (−0.8 and −0.6 °C), whereas δ¹⁸O_seawater and SSS had higher values in summer (+ 0.3‰, + 0.6 psu) and in winter (+ 0.8‰, + 1.9 psu) than present-day values. These results are consistent with other marine records for the mid-Holocene of the low and midlatitudes in the northwestern Pacific. Such regional conditions indicate that the East Asian summer and winter monsoons were more intense in the mid-Holocene, which was likely a function of the mid-Holocene insolation regime.     

Environmental memory and a possible seasonal bias in the stable isotope composition of (U-Th)/He-dated goethite from the Canadian Arctic

Goethite (Ax-2) from Axel Heiberg Island (∼80°N) on the margin of the Arctic Ocean is the dominant mineral in a sample of “petrified” Eocene wood, but U, Th, and He measurements suggest that the goethite (α-FeOOH) crystallized in the latest Miocene/Pliocene (ca. 5.5 to 2.8 Ma). Measured δD and δ¹⁸O values of Ax-2 are −221 (±6)‰ and −9.6 (±0.5)‰, respectively. The inferred δD and δ¹⁸O values of the ancient water were about −139‰ and −18.6‰, respectively, with a calculated temperature of crystallization of 3 (±5)°C, which compares with the modern summer (J-J-A) temperature of 3 °C and contrasts with a modern MAT of −19 °C. Published results from various biological proxies on nearby Ellesmere Island indicate a Pliocene (∼4 Ma) MAT of either −6 or −0.4 °C and corresponding seasonal amplitudes of about 18 or 13 °C. A conductive heat flow model suggests that a temperature of 3 °C could represent goethite crystallization at depths of ∼100-200 cm (for MAT = −6 °C) or ∼250-450 cm (for MAT = −0.4 °C) over seasonally restricted intervals of time.The δ¹⁸O value of the Ax-2 water (−18.6‰) is more positive than the modern J-J-A precipitation (−22‰). In combination, the paleotemperatures and δ¹⁸O values of ancient waters (from Ax-2 and published results from three Eocene or Pliocene proxy sites on Axel Heiberg and Ellesmere Islands) are consistent with a warm season bias in those isotopic proxies. The results are also consistent with higher proportions of J-J-A precipitation in the annual total. If so, this emphasizes the importance of seasonality at high latitudes even in times of warmer global climates, and suggests that the Arctic hydrologic cycle, as expressed in the seasonal distribution and isotopic composition of precipitation (perhaps modified by a warmer Arctic Ocean), differed from modern.The δ¹³C value of the Fe(CO₃)OH component in the Ax-2 goethite is +6.6‰, which is much more positive than expected if crystallizing goethite incorporated CO₂ derived primarily from oxidation of relict Eocene wood with δ¹³C values of about −24‰. This apparent paradox may be resolved if the goethite is a product of oxidation of ¹³C-rich siderite, which had previously replaced wood in an Eocene methanogenic burial environment. Thus, the goethite retains a carbon isotope “memory” of a diagenetic Eocene event, but a δD and δ¹⁸O record of the latest Miocene/Pliocene Arctic climate.     

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.     

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

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.     

First experimental determination of iron isotope fractionation between hematite and aqueous solution at hydrothermal conditions

Although iron isotopes provide a new powerful tool for tracing a variety of geochemical processes, the unambiguous interpretation of iron isotope ratios in natural systems and the development of predictive theoretical models require accurate data on equilibrium isotope fractionation between fluids and minerals. We investigated Fe isotope fractionation between hematite (Fe₂O₃) and aqueous acidic NaCl fluids via hematite dissolution and precipitation experiments at temperatures from 200 to 450 °C and pressures from saturated vapor pressure (P_sat) to 600 bar. Precipitation experiments at 200 °C and P_sat from aqueous solution, in which Fe aqueous speciation is dominated by ferric iron (Fe^III) chloride complexes, show no detectable Fe isotope fractionation between hematite and fluid, Δ⁵⁷Fe_{fluid-hematite} = δ⁵⁷Fe_fluid − δ⁵⁷Fe_hematite = 0.01 ± 0.08‰ (2 × standard error, 2SE). In contrast, experiments at 300 °C and P_sat, where ferrous iron chloride species (FeCl₂ and FeCl⁺) dominate in the fluid, yield significant fluid enrichment in the light isotope, with identical values of Δ⁵⁷Fe_{fluid-hematite} = −0.54 ± 0.15‰ (2SE) both for dissolution and precipitation runs. Hematite dissolution experiments at 450 °C and 600 bar, in which Fe speciation is also dominated by ferrous chloride species, yield Δ⁵⁷Fe_{fluid-hematite} values close to zero within errors, 0.15 ± 0.17‰ (2SE). In most experiments, chemical, redox, and isotopic equilibrium was attained, as shown by constancy over time of total dissolved Fe concentrations, aqueous Fe^II and Fe^III fractions, and Fe isotope ratios in solution, and identical Δ⁵⁷Fe values from dissolution and precipitation runs. Our measured equilibrium Δ⁵⁷Fe_{fluid-hematite} values at different temperatures, fluid compositions and iron redox state are within the range of fractionations in the system fluid-hematite estimated using reported theoretical β-factors for hematite and aqueous Fe species and the distribution of Fe aqueous complexes in solution. These theoretical predictions are however affected by large discrepancies among different studies, typically ±1‰ for the Δ⁵⁷Fe _{Fe(aq)-hematite} value at 200 °C. Our data may thus help to refine theoretical models for β-factors of aqueous iron species. This study provides the first experimental calibration of Fe isotope fractionation in the system hematite-saline aqueous fluid at elevated temperatures; it demonstrates the importance of redox control on Fe isotope fractionation at hydrothermal conditions.     

Controls, variation, and a record of climate change in detailed stable isotope record in a single bryozoan skeleton

The long-lived (about 20 yr) bryozoan Adeonellopsis sp. from Doubtful Sound, New Zealand, precipitates aragonite in isotopic equilibrium with seawater, exerting no metabolic or kinetic effects. Oxygen isotope ratios (δ¹⁸O) in 61 subsamples (along three branches of a single unaltered colony) range from −0.09 to +0.68‰ PDB (mean = +0.36‰ PDB). Carbon isotope ratios (δ¹³C) range from +0.84 to +2.18‰ PDB (mean = +1.69‰ PDB). Typical of cool-water carbonates, δ¹⁸O-derived water temperatures range from 14.2 to 17.5 °C. Adeonellopsis has a minimum temperature growth threshold of 14 °C, recording only a partial record of environmental variation. By correlating seawater temperatures derived from δ¹⁸O with the Southern Oscillation Index, however, we were able to detect major events such as the 1983 El Niño. Interannual climatic variation can be recorded in skeletal carbonate isotopes. The range of within-colony isotopic variability found in this study (0.77‰ in δ¹⁸O and 1.34 in δ¹³C) means that among-colony variation must be treated cautiously. Temperate bryozoan isotopes have been tested in less than 2% of described extant species — this highly variable phylum is not yet fully understood.     

Kinetic isotopic fractionation during carbonate dissolution in laboratory experiments: Implications for detection of microbial CO2 signatures using δC-DIC

Laboratory experiments on reagent-grade calcium carbonate and carbonate rich glacial sediments demonstrate previously unreported kinetic fractionation of carbon isotopes during the initial hydrolysis and early stages of carbonate dissolution driven by atmospheric CO₂. There is preferential dissolution of Ca¹²CO₃ during hydrolysis, resulting in δ¹³C-DIC values that are significantly lighter isotopically than the bulk carbonate. The fractionation factor for this kinetic isotopic effect is defined as ε_carb. ε_carb is greater on average for glacial sediments (−17.4‰) than for calcium carbonate (−7.8‰) for the < 63 μm size fraction, a sediment concentration of 5 g L⁻¹ and closed system conditions at 5°C. This difference is most likely due to the preferential dissolution of highly reactive ultra-fine particles with damaged surfaces that are common in subglacial sediments. The kinetic isotopic fractionation has a greater impact on δ¹³C-DIC at higher CaCO₃:water ratios and is significant during at least the first 6 h of carbonate dissolution driven by atmospheric CO₂ at sediment concentrations of 5 g L⁻¹. Atmospheric CO₂ dissolving into solution following carbonate hydrolysis does not exhibit any significant equilibrium isotopic fractionation for at least ∼ 6 h after the start of the experiment at 5°C. This is considerably longer than previously reported in the literature. Thus, kinetic fractionation processes will likely dominate the δ¹³C-DIC signal in natural environments where rock:water contact times are short <6-24 h (e.g., glacial systems, headwaters in fluvial catchments) and there is an excess of carbonate in the sediments. It will be difficult apply conventional isotope mass balance techniques in these types of environment to identify microbial CO₂ signatures in DIC from δ¹³C-DIC data.     

Evidence for ocean acidification in the Great Barrier Reef of Australia

Geochemical records preserved in the long-lived carbonate skeleton of corals provide one of the few means to reconstruct changes in seawater pH since the commencement of the industrial era. This information is important in not only determining the response of the surface oceans to ocean acidification from enhanced uptake of CO₂, but also to better understand the effects of ocean acidification on carbonate secreting organisms such as corals, whose ability to calcify is highly pH dependent. Here we report an ∼200 year δ¹¹B isotopic record, extracted from a long-lived Porites coral from the central Great Barrier Reef of Australia. This record covering the period from 1800 to 2004 was sampled at yearly increments from 1940 to the present and 5-year increments prior to 1940. The δ¹¹B isotopic compositions reflect variations in seawater pH, and the δ¹³C changes in the carbon composition of surface water due to fossil fuel burning over this period. In addition complementary Ba/Ca, δ¹⁸O and Mg/Ca data was obtained providing proxies for terrestrial runoff, salinity and temperature changes over the past 200 years in this region. Positive thermal ionization mass spectrometry (PTIMS) method was utilized in order to enable the highest precision and most accurate measurements of δ¹¹B values. The internal precision and reproducibility for δ¹¹B of our measurements are better than ±0.2‰ (2σ), which translates to a precision of better than ±0.02 pH units. Our results indicate that the long-term pre-industrial variation of seawater pH in this region is partially related to the decadal-interdecadal variability of atmospheric and oceanic anomalies in the Pacific. In the periods around 1940 and 1998 there are also rapid oscillations in δ¹¹B compositions equivalent changes in pH of almost 0.5 U. The 1998 oscillation is co-incident with a major coral bleaching event indicating the sensitivity of skeletal δ¹¹B compositions to loss of zooxanthellate symbionts. Importantly, from the 1940s to the present-day, there is a general overall trend of ocean acidification with pH decreasing by about 0.2-0.3 U, the range being dependent on the value assumed for the fractionation factor α_(B3-B4) of the boric acid and borate species in seawater. Correlations of δ¹¹B with δ¹³C during this interval indicate that the increasing trend towards ocean acidification over the past 60 years in this region is the result of enhanced dissolution of CO₂ in surface waters from the rapidly increasing levels of atmospheric CO₂, mainly from fossil fuel burning. This suggests that the increased levels of anthropogenic CO₂ in atmosphere has already caused a significant trend towards acidification in the oceans during the past decades. Observations of surprisingly large decreases in pH across important carbonate producing regions, such as the Great Barrier Reef of Australia, raise serious concerns about the impact of Greenhouse gas emissions on coral calcification.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.     

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.