Calcium isotope record of Phanerozoic oceans: Implications for chemical evolution of seawater and its causative mechanisms

A total of 280 brachiopods of Ordovician to Cretaceous age, complemented by published data from belemnites and planktonic foraminifera, are used to reconstruct the evolution of calcium isotope composition of seawater (δ^44/40Ca_SW) over the Phanerozoic. The compiled δ^44/40Ca_SW record shows a general increase from ∼1.3‰ (NIST SRM 915a) at the beginning of the Ordovician to ∼2‰ at present. Superimposed on this trend is a major long-term positive excursion from the Early Carboniferous to Early Permian as well as several short-term, mostly negative, oscillations.A numerical model of the global cycles of calcium, carbon, magnesium and strontium was used to estimate whether the recorded δ^44/40Ca_SW variations can be explained by varying magnitudes of input and output fluxes of calcium to the oceans. The model uses the record of marine ⁸⁷Sr/⁸⁶Sr ratios as proxy for seafloor spreading rates, a record of oceanic Mg/Ca ratios to estimate rates of dolomite formation, and reconstructed atmospheric CO₂, discharge and erosion rates to estimate continental weathering fluxes.The model results indicate that varying magnitudes of the calcium input and output fluxes cannot explain the observed δ^44/40Ca_SW trends, suggesting that the isotope signatures of these fluxes must also have changed. As a possible mechanism we suggest variable isotope fractionation in the sedimentary output flux controlled by the dominant mineralogy in marine carbonate deposits, i.e. the oscillating ‘calcite-aragonite seas’. The ultimate control of the calcium isotope budget of the Phanerozoic oceans appears to have been tectonic processes, specifically variable rates of oceanic crust production that modulated the hydrothermal calcium flux and the oceanic Mg/Ca ratio, which in turn controlled the dominant mineralogy of marine carbonates, hence the δ^44/40Ca_SW. As to the causes of the short-term oscillations recorded in the secular δ^44/40Ca_SW trend, we tentatively propose that these are related to variable rates of dolomite formation and/or to changing chemical composition of the riverine flux, in particular and ratios, induced by variable proportions of silicate vs. carbonate weathering rates on the continents.     

Experimental investigation of the effects of temperature and ionic strength on Mo isotope fractionation during adsorption to manganese oxides

The Mo stable isotope system is being applied to study changes in ocean redox. Such applications implicitly assume that Mo isotope fractionation in aqueous systems is relatively insensitive to frequently changing environmental variables such as temperature (T) and ionic strength (I). A major driver of fractionation is the adsorption of Mo to Mn oxyhydroxide surfaces [Barling J. and Anbar A. D. (2004) Molybdenum isotope fractionation during adsorption by manganese oxides. Earth Planet. Sci. Lett.217(3-4), 315-329]. Here, we report the results of experiments that determine the extent to which Mo isotope fractionation during adsorption of Mo to the Mn oxyhydroxide mineral birnessite is sensitive to T and I. The results are compared to new predictions from quantum chemical computations. We measured fractionation from 1 to 50 °C at I = 0.1 m and found that Δ^97/95Mo_{dissolved-adsorbed} varies from 1.9‰ to 1.6‰ over this temperature range. Experiments were also performed at 25 °C in synthetic seawater (I = 0.7); fractionation at this condition was the same within analytical error as in low ionic strength experiments. These findings confirm that the Mo isotope fractionation during adsorption to Mn oxyhydroxides is relatively insensitive to variations and T and I over environmentally relevant ranges. To relate these findings to potential mechanisms of Mo isotope fractionation, we also report results for density functional theory computations of the fractionation between and various possible structures of molybdic acid as a function of temperature. Because no plausible species fractionates from with a magnitude matching the experiments, we are left with three possibilities to explain the fractionation: (1) solvation effects on the vibrational frequencies of aqueous species considered thus far are significant, such that our calculations in vacuo yield inaccurate fractionations; (2) a trace aqueous species not yet considered fractionates from and then adsorbs to birnessite; or (3) a surface complex not present in solution forms on birnessite in which Mo is not tetrahedrally coordinated. Our findings help validate assumptions underlying paleoceanographic applications of the Mo isotope system and also lead us closer to understanding the mechanism of isotope fractionation during adsorption of Mo to Mn oxyhydroxides.     

The molecular mechanism of Mo isotope fractionation during adsorption to birnessite

Fractionation of Mo isotopes during adsorption to manganese oxides is a primary control on the global ocean Mo isotope budget. Previous attempts to explain what drives the surprisingly large isotope effect δ^97/95Mo_dissolved-δ^97/95Mo_adsorbed=1.8‰ have not successfully resolved the fractionation mechanism. New evidence from extended X-ray absorption fine structure analysis and density functional theory suggests that Mo forms a polymolybdate complex on the surfaces of experimental and natural samples. Mo in this polynuclear structure is in distorted octahedral coordination, while Mo remaining in solution is predominantly in tetrahedral coordination as . Our results indicate that the difference in coordination environment between dissolved Mo and adsorbed Mo is the cause of isotope fractionation. The molecular mechanism of metal isotope fractionation in this system should enable us to explain and possibly predict metal isotope effects in other systems where transition metals adsorb to mineral surfaces.     

Controls on ostracod valve geochemistry: Part 2. Carbon and oxygen isotope compositions

The stable carbon and oxygen isotope compositions of fossil ostracods are powerful tools to estimate past environmental and climatic conditions. The basis for such interpretations is that the calcite of the valves reflects the isotopic composition of water and its temperature of formation. However, calcite of ostracods is known not to form in isotopic equilibrium with water and different species may have different offsets from inorganic precipitates of calcite formed under the same conditions. To estimate the fractionation during ostracod valve calcification, the oxygen and carbon isotope compositions of 15 species living in Lake Geneva were related to their autoecology and the environmental parameters measured during their growth. The results indicate that: (1) Oxygen isotope fractionation is similar for all species of Candoninae with an enrichment in ¹⁸O of more than 3‰ relative to equilibrium values for inorganic calcite. Oxygen isotope fractionation for Cytheroidea is less discriminative relative to the heavy oxygen, with enrichments in ¹⁸O for these species of 1.7 to 2.3‰. Oxygen isotope fractionations for Cyprididae are in-between those of Candoninae and Cytheroidea. The difference in oxygen isotope fractionation between ostracods and inorganic calcite has been interpreted as resulting from a vital effect. (2) Comparison with previous work suggests that oxygen isotope fractionation may depend on the total and relative ion content of water. (3) Carbon isotope compositions of ostracod valves are generally in equilibrium with DIC. The specimens’ δ¹³C values are mainly controlled by seasonal variations in δ¹³C_DIC of bottom water or variation thereof in sediment pore water. (4) Incomplete valve calcification has an effect on carbon and oxygen isotope compositions of ostracod valves. Preferential incorporation of at the beginning of valve calcification may explain this effect. (5) Results presented here as well as results from synthetic carbonate growth indicate that different growth rates or low pH within the calcification site cannot be the cause of oxygen isotope ‘vital effects’ in ostracods. Two mechanisms that might enrich the ¹⁸O of ostracod valves are deprotonation of that may also contribute to valve calcification, and effects comparable to salt effects with high concentrations of Ca and/or Mg within the calcification site that may also cause a higher temperature dependency of oxygen isotope fractionation.     

An integrated sulfur isotope model for Namibian shelf sediments

In this study the sulfur cycle in the organic-rich mud belt underlying the highly productive upwelling waters of the Namibian shelf is quantified using a 1D reaction-transport model. The model calculates vertical concentration and reaction rate profiles in the top 500 cm of sediment which are compared to a comprehensive dataset which includes carbon, sulfur, nitrogen and iron compounds as well as sulfate reduction (SR) rates and stable sulfur isotopes (³²S, ³⁴S). The sulfur dynamics in the well-mixed surface sediments are strongly influenced by the activity of the large sulfur bacteria Thiomargaritanamibiensis which oxidize sulfide (H₂S) to sulfate () using sea water nitrate () as the terminal electron acceptor. Microbial sulfide oxidation (SOx) is highly efficient, and the model predicts intense cycling between and H₂S driven by coupled SR and SOx at rates exceeding 6.0 mol S m⁻² y⁻¹. More than 96% of the SR is supported by SOx, and only 2-3% of the pool diffuses directly into the sediment from the sea water. A fraction of the produced by Thiomargarita is drawn down deeper into the sediment where it is used to oxidize methane anaerobically, thus preventing high methane concentrations close to the sediment surface. Only a small fraction of total H₂S production is trapped as sedimentary sulfide, mainly pyrite (FeS₂) and organic sulfur (S_org) (∼0.3 wt.%), with a sulfur burial efficiency which is amongst the lowest values reported for marine sediments (<1%). Yet, despite intense SR, FeS₂ and S_org show an isotope composition of ∼5 ‰ at 500 cm depth. These heavy values were simulated by assuming that a fraction of the solid phase sulfur exchanges isotopes with the dissolved sulfide pool. An enrichment in H₂S of ³⁴S towards the sediment-water interface suggests that Thiomargarita preferentially remove H₂³²S from the pore water. A fractionation of 20-30‰ was estimated for SOx (ε_SOx) with the model, along with a maximum fractionation for SR (ε_SR-max) of 100‰. These values are far higher than previous laboratory-based estimates for these processes. Mass balance calculations indicate negligible disproportionation of autochthonous elemental sulfur; an explanation routinely cited in the literature to account for the large fractionations in SR. Instead, the model indicates that repeated multi-stepped sulfide oxidation and intracellular disproportionation by Thiomargarita could, in principle, allow the measured isotope data to be simulated using much lower fractionations for ε_SOx (5‰) and ε_SR (78‰).     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

The origin of Zn isotope fractionation in sulfides

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

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

Inverse kinetic isotope fractionation during bacterial nitrite oxidation

Natural abundance stable isotopes in nitrate (), nitrite (), and nitrous oxide (N₂O) have been used to better understand the cycling of nitrogen in marine and terrestrial environments. However, in order to extract the greatest information from the distributions of these isotopic species, the kinetic isotope effects for each of the relevant microbial reactions are needed. To date, kinetic isotope effects for nitrite oxidation and anaerobic ammonium oxidation (anammox) have not been reported. In this study, the nitrogen isotope effect was measured for microbial nitrite oxidation to nitrate. Nitrite oxidation is the second step in the nitrification process, and it plays a key role in the regeneration of nitrate in the ocean. Surprisingly, nitrite oxidation occurred with an inverse kinetic isotope effect, such that the residual nitrite became progressively depleted in ¹⁵N as the reaction proceeded. Three potential explanations for this apparent inverse kinetic isotope effect were explored: (1) isotope exchange equilibrium between nitrite and nitrous acid prior to reaction, (2) reaction reversibility at the enzyme level, and (3) true inverse kinetic fractionation. Comparison of experimental data to ab initio calculations and theoretical predictions leads to the conclusion that the fractionation is most likely inverse at the enzyme level. Inverse kinetic isotope effects are rare, but the experimental observations reported here agree with kinetic isotope theory for this simple N-O bond-forming reaction. Nitrite oxidation is therefore fundamentally different from all other microbial processes in which N isotope fractionation has been studied. The unique kinetic isotope effect for nitrite oxidation should help to better identify its role in the cycling of nitrite in ocean suboxic zones, and other environments in which nitrite accumulates.     

Carbon isotope fractionation in synthetic magnesian calcite

Mg-calcite was precipitated at 25 °C in closed system, free-drift experiments, from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The carbon stable isotope composition of bulk solid and solution were analyzed from subsamples collected during time course experiments of 24 h duration. Considering only the Mg-content and δ¹³C values for the bulk solid, the carbon isotope fractionation factor for the Mg-calcite-HCO_3(aq)⁻ system (as ) increased with average mol percentage of Mg (X_Mg) in the solid at a rate of (0.024 ± 0.011) per mol% MgCO₃. Extrapolation of this relationship to the pure calcite end member yields a value of 0.82 ± 0.09, which is similar to published values for the calcite-HCO_3(aq)⁻ system. Although did not vary for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol m⁻² h⁻¹, it was not possible to hold Mg-content of the solid constant, so kinetic effect on 10³ ln α could not be evaluated from these experiments.     

Confirmation of mass-independent Ni isotopic variability in iron meteorites

We report high-precision analyses of internally-normalised Ni isotope ratios in 12 bulk iron meteorites. Our measurements of ⁶⁰Ni/⁶¹Ni, ⁶²Ni/⁶¹Ni and ⁶⁴Ni/⁶¹Ni normalised to ⁵⁸Ni/⁶¹Ni and expressed in parts per ten thousand (?) relative to NIST SRM 986 as and , vary by 0.146, 0.228 and 0.687, respectively. The precision on a typical analysis is 0.03?, 0.05? and 0.08? for , and , respectively, which is comparable to our sample reproducibility. We show that this ‘mass-independent’ Ni isotope variability cannot be ascribed to interferences, inaccurate correction of instrumental or natural mass-dependent fractionation, fractionation controlled by nuclear field shift effects, nor the influence of cosmic ray spallation. These results thus document the presence of mass-independent Ni isotopic heterogeneity in bulk meteoritic samples, as previously proposed by Regelous et al. (2008) (EPSL 272, 330-338), but our new analyses are more precise and include determination of ⁶⁴Ni. Intriguingly, we find that terrestrial materials do not yield homogenous internally-normalised Ni isotope compositions, which, as pointed out by Young et al. (2002) (GCA 66, 1095-1104), may be the expected result of using the exponential (kinetic) law and atomic masses to normalise all fractionation processes. The certified Ni isotope reference material NIST SRM 986 defines zero in this study, while appropriate ratios for the bulk silicate Earth are given by the peridotites JP-1 and DTS-2 and, relative to NIST SRM 986, yield deviations in , and of −0.006?, 0.036? and 0.119?, respectively. There is a strong positive correlation between and in iron meteorites analyses, with a slope of 3.03 ± 0.71. The variations of Ni isotope anomalies in iron meteorites are consistent with heterogeneous distribution of a nucleosynthetic component from a type Ia supernova into the proto-solar nebula.     

A new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

Modelling carbon isotopes of carbonates in cave drip water

C isotopes in cave drip water are affected by both the C isotope composition of soil air and host rock carbonate. Furthermore, the C isotope composition of cave drip water strongly depends on the calcite dissolution system, i.e., open, closed and intermediate conditions. Here, we present a calcite dissolution model, which calculates the ¹⁴C activity and δ¹³C value of the dissolved inorganic carbon of the drip water. The model is based on the chemical equations describing calcite dissolution (). The most important improvement, relative to previous models, is the combination of the open and closed system conditions in order to simulate the C isotope composition during intermediate states of calcite dissolution and the application to carbon isotope measurements on cave drip waters from Grotta di Ernesto, Italy. The major changes in the C isotope composition of the drip water occur in response to variations in the open-closed system ratio. Additionally, the ¹⁴C activity and the δ¹³C value of the drip water depend on changes in the partial pressure of soil CO₂. Radiocarbon and δ¹³C values of the Grotta di Ernesto drip water are well reproduced by the model.     

New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions

Several studies have shown that there is a strong relationship between the distribution of crenarchaeotal isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) and sea surface temperature (SST). Based on this, a ratio of certain GDGTs, called TEX₈₆ (TetraEther indeX of tetraethers consisting of 86 carbon atoms), was developed as a SST proxy. In this study, we determined the distribution of crenarchaeotal isoprenoid GDGTs in 116 core-top sediments mostly from (sub)polar oceans and combined these data with previously published core-top data. Using this extended global core-top dataset (n = 426), we re-assessed the relationship of crenarchaeal isoprenoid GDGTs with SST. We excluded data from the Red Sea from the global core-top dataset to define new indices and calibration models, as the Red Sea with its elevated salinity appeared to behave differently compared to other parts of the oceans. We tested our new indices and calibration models on three different paleo datasets, representing different temperature ranges. Our results indicate that the crenarchaeol regio-isomer plays a more important role for temperature adaptation in (sub)tropical oceans than in (sub)polar oceans, suggesting that there may be differences in membrane adaptation of the resident crenarchaeotal communities at different temperatures. We, therefore, suggest to apply two different calibration models. For the whole calibration temperature range (−3 to 30 °C), a modified version of TEX₈₆ with a logarithmic function which does not include the crenarchaeol regio-isomer, called , is shown to correlate best with SST: (r² = 0.86, n=396, p <0.0001). Application of on sediments from the subpolar Southern Ocean results in realistic absolute SST estimates and a similar SST trend compared to a diatom SST record from the same core. , which is defined as the logarithmic function of TEX₈₆, yields the best correlation with SST, when the data from the (sub)polar oceans are removed: (r² = 00.87, n = 255, p < 0.0001). Furthermore, gives the best correlation for mescosm data with temperatures ranging between 10 and 46 °C. For Quaternary sediments from the tropical Arabian Sea, both and yield similar trends and SST estimates. However, the extrapolation of calibration on a sediment record from a greenhouse world ocean predicts more reliable absolute SST estimates and relative SST changes in agreement with estimates based on the δ¹⁸O of planktonic foraminifera. Based on the comparison of and derived SSTs using the core top data, we recommend applying above 15 °C and below 15 °C. In cases where paleorecords encompass temperatures both below and above 15 °C, we suggest to use .     

Phosphoric acid fractionation factors for smithsonite and cerussite between 25 and 72°C

The intramolecular kinetic oxygen isotope fractionation between CO₂ and CO₃²⁻ during reaction of phosphoric acid with natural smithsonite (ZnCO₃) and cerussite (PbCO₃) has been determined between 25 and 72°C. While cerussite decomposes in phosphoric acid within a few hours at 25°C, smithsonite reacts very slowly with the acid at 25°C providing yields of CO₂ < 25% after 2 weeks. The low yields result in a low precision for oxygen isotope measurements of the acid-liberated CO₂ (±1.65‰, 1σ, n = 9). The yield and reproducibility of oxygen isotope values of the acid-liberated CO₂ from smithsonite can be improved, the latter to ∼±0.15‰, by increasing the reaction temperature to 50°C for 12 h or to 72°C for 1 h. Our new phosphoric acid fractionation factor for natural cerussite at 25°C deviates significantly from a previously published value on synthetic material. The temperature dependence of the oxygen isotope factionation factor, α between acid-liberated CO₂ and carbonate at 25 to 72°C is given by the following equations

     

Iron isotope fractionation and atom exchange during sorption of ferrous iron to mineral surfaces

The application of stable Fe isotopes as a tracer of the biogeochemical Fe cycle necessitates a mechanistic knowledge of natural fractionation processes. We studied the equilibrium Fe isotope fractionation upon sorption of Fe(II) to aluminum oxide (γ-Al₂O₃), goethite (α-FeOOH), quartz (α-SiO₂), and goethite-loaded quartz in batch experiments, and performed continuous-flow column experiments to study the extent of equilibrium and kinetic Fe isotope fractionation during reactive transport of Fe(II) through pure and goethite-loaded quartz sand. In addition, batch and column experiments were used to quantify the coupled electron transfer-atom exchange between dissolved Fe(II) (Fe(II)_aq) and structural Fe(III) of goethite. All experiments were conducted under strictly anoxic conditions at pH 7.2 in 20 mM MOPS (3-(N-morpholino)-propanesulfonic acid) buffer and 23 °C. Iron isotope ratios were measured by high-resolution MC-ICP-MS. Isotope data were analyzed with isotope fractionation models. In batch systems, we observed significant Fe isotope fractionation upon equilibrium sorption of Fe(II) to all sorbents tested, except for aluminum oxide. The equilibrium enrichment factor, , of the Fe(II)_sorb-Fe(II)_aq couple was 0.85 ± 0.10‰ (±2σ) for quartz and 0.85 ± 0.08‰ (±2σ) for goethite-loaded quartz. In the goethite system, the sorption-induced isotope fractionation was superimposed by atom exchange, leading to a δ^56/54Fe shift in solution towards the isotopic composition of the goethite. Without consideration of atom exchange, the equilibrium enrichment factor was 2.01 ± 0.08‰ (±2σ), but decreased to 0.73 ± 0.24‰ (±2σ) when atom exchange was taken into account. The amount of structural Fe in goethite that equilibrated isotopically with Fe(II)_aq via atom exchange was equivalent to one atomic Fe layer of the mineral surface (∼3% of goethite-Fe). Column experiments showed significant Fe isotope fractionation with δ^56/54Fe(II)_aq spanning a range of 1.00‰ and 1.65‰ for pure and goethite-loaded quartz, respectively. Reactive transport of Fe(II) under non-steady state conditions led to complex, non-monotonous Fe isotope trends that could be explained by a combination of kinetic and equilibrium isotope enrichment factors. Our results demonstrate that in abiotic anoxic systems with near-neutral pH, sorption of Fe(II) to mineral surfaces, even to supposedly non-reactive minerals such as quartz, induces significant Fe isotope fractionation. Therefore we expect Fe isotope signatures in natural systems with changing concentration gradients of Fe(II)_aq to be affected by sorption.     

Experimental studies of oxygen isotope fractionation in the carbonic acid system at 15°, 25°, and 40°C

In light of recent studies that show oxygen isotope fractionation in carbonate minerals to be a function of HCO₃⁻ and CO₃²⁻ concentrations, the oxygen isotope fractionation and exchange between water and components of the carbonic acid system (HCO₃⁻, CO₃²⁻, and CO_2(aq)) were investigated at 15°, 25°, and 40°C. To investigate oxygen isotope exchange between HCO₃⁻, CO₃²⁻, and H₂O, NaHCO₃ solutions were prepared and the pH was adjusted over a range of 2 to 12 by the addition of small amounts of HCl or NaOH. After thermal, chemical, and isotopic equilibrium was attained, BaCl₂ was added to the NaHCO₃ solutions. This resulted in immediate BaCO₃ precipitation; thus, recording the isotopic composition of the dissolved inorganic carbon (DIC). Data from experiments at 15°, 25°, and 40°C (1 atm) show that the oxygen isotope fractionation between HCO₃⁻ and H₂O as a function of temperature is governed by the equation:

     

An experimental study of oxygen isotope fractionation between inorganically precipitated aragonite and water at low temperatures

To determine oxygen isotope fractionation between aragonite and water, aragonite was slowly precipitated from Ca(HCO₃)₂ solution at 0 to 50°C in the presence of Mg²⁺ or SO₄²⁻. The phase compositions and morphologies of synthetic minerals were detected by X-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques. The effects of aragonite precipitation rate and excess dissolved CO₂ gas in the initial Ca(HCO₃)₂ solution on oxygen isotope fractionation between aragonite and water were investigated. For the CaCO₃ minerals slowly precipitated by the CaCO₃ or NaHCO₃ dissolution method at 0 to 50°C, the XRD and SEM analyses show that the rate of aragonite precipitation increased with temperature. Correspondingly, oxygen isotope fractionations between aragonite and water deviated progressively farther from equilibrium. Additionally, an excess of dissolved CO₂ gas in the initial Ca(HCO₃)₂ solution results in an increase in apparent oxygen isotope fractionations. As a consequence, the experimentally determined oxygen isotope fractionations at 50°C indicate disequilibrium, whereas the relatively lower fractionation values obtained at 0 and 25°C from the solution with less dissolved CO₂ gas and low precipitation rates indicate a closer approach to equilibrium. Combining the lower values at 0 and 25°C with previous data derived from a two-step overgrowth technique at 50 and 70°C, a fractionation equation for the aragonite-water system at 0 to 70°C is obtained as follows: