Experimental study of nickel(II) interaction with calcite: Adsorption and coprecipitation

Partitioning of Ni in calcite, CaCO₃, was evaluated with the aim of collecting data on partition and distribution coefficients and to enhance understanding about the interaction of Ni with the calcite surface and further incorporation into the bulk. This information will aid in the interpretation of geological processes for safety assessment of waste repositories and contamination of groundwater. Coprecipitation experiments were carried out by the constant addition method at 25 °C and pCO₂ = 1 and 10^−3.5 atm. Ni was moderately partitioned from solution into calcite. For dilute solid solutions (X_Ni < 0.001), Ni partition coefficients were estimated to be ∼1 and found to be weakly dependent on calcite precipitation rate in the range of 3-230 nmol m⁻² s⁻¹. Ni molar fraction in the solid is directly proportional to Ni concentration in the solution. The fit of the data to such a model is good evidence that Ni is taken up as a true solid solution, not simply by physical trapping.     

Solid solutions of trace Eu(III) in calcite: Thermodynamic evaluation of experimental data over a wide range of pH and pCO2

The thermodynamics of dilute Eu-calcite solid solutions formed under widely different pH-pCO₂ conditions at T = 25°C and p = 1 bar were investigated using three sets of Eu(III) uptake experiments, two of which were taken from the literature: (a) recrystallization in synthetic cement pore water at pH ∼ 13 and pCO₂ ∼ 10⁻¹³ bar (this work); (b) coprecipitation in 0.1 M NaClO₄ at pH ∼ 6 and pCO₂ ∼ 1 bar; (c) coprecipitation in synthetic seawater at pH ∼ 8 and pCO₂ ranging from 3 × 10⁻⁴ to 0.3 bar.Solid solution formation was modeled using the Gibbs energy minimization (GEM) method. In a first step (“forward” modeling), we tested ideal binary solid solution models between calcite and the Eu end-members Eu₂(CO₃)_3, EuNa(CO₃)₂, Eu(OH)CO₃ or Eu(OH)₃, for which solids with independently measured solubility products exist. None of these four binary solid solutions was capable of reproducing all three experimental datasets simultaneously. In a second step (“inverse” modeling), ideal binary solid solutions were constructed between calcite and the candidate Eu end-members EuO(OH), EuH(CO₃)₂ and EuO(CO₃)_0.5, for which no independent solubility products are available. For each single data point and each of these end-members, a free energy of formation with inherent activity coefficient term ( = G_α^o + RT lnγ_α) was estimated from “dual thermodynamic” GEM calculations. The statistical mean of was then calculated for each of the three datasets. A specific end-member was considered to be acceptable if a standard deviation of ± 2 kJ mol⁻¹ or less resulted for each single dataset, and if the mean -values calculated for the three datasets coincided. No binary solid solution with any of the seven above mentioned end-members proved to satisfy these criteria.The third step in our analysis involved consideration of ternary solid solutions with CaCO₃ as the major end-member and any two of the seven considered Eu trace end-members. It was found that the three datasets can only be reproduced simultaneously with the ternary ideal solid solution EuH(CO₃)₂ - EuO(OH) - CaCO₃, setting = −1773 kJ mol⁻¹ and = −955 kJ mol⁻¹, whereas all other end-member combinations failed. Our results are consistent with time-resolved laser fluorescence data for Cm(III) and Eu(III) indicating that two distinct species are incorporated in calcite: one partially hydrated, the other completely dehydrated. In conclusion, our study shows that substitution of trivalent for divalent cations in carbonate crystal structures is a more complex process than the classical isomorphic divalent-divalent substitution and may need consideration of multicomponent solid solution models.     

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

Sorption of Eu on Na- and Ca-montmorillonites: experimental investigations and modelling with cation exchange and surface complexation

The 2-site protolysis no electrostatics surface complexation and cation exchange (2SPNE/CE) model used in previous work to model the sorption of Ni and Zn on Na- and Ca-montmorillonites was applied to sorption edges and isotherms measured for Eu on these two montmorillonite forms. The aim was to further test the applicability of the sorption model on a trivalent element with a more complex aqueous chemistry. An additional reason for choosing Eu was that it is considered to be a good chemical analogue for other lanthanides and trivalent actinides. With site types, site capacities, and protolysis constants fixed at the values in the Ni/Zn studies, all of the measured sorption edge data could be modelled using cation exchange and the monodentate surface species, ≡S^SOEu²⁺, ≡S^SOEuOH⁺ and ≡S^SOEu(OH)₃⁻, on the strong site type. However, an additional modelling study showed that the same data were almost equally well described by considering bidentate surface complexes, (≡S^SO)₂Eu⁺ and (≡S^SO)₂Eu(OH)₂⁻, and cation exchange. To model the sorption isotherm measurements up to pH = 7.2, only one additional weak site surface complex was required, ≡S^W1OEu²⁺ for the monodentate case and (≡S^W1O)₂Eu⁺ for the bidentate case. Selectivity coefficients are given for Eu³⁺- Ca²⁺ and Eu³⁺- Na⁺ exchange on the planar sites and surface complexation constants for monodentate and bidentate Eu surface species on the edge sites of montmorillonite.     

Arsenic uptake by natural calcite: An XAS study

As-bearing travertine rocks from Tuscany (Italy), where previous studies suggested the existence of a CO₃²⁻ ⇔ AsO₃³⁻ substitution in the calcite lattice, were investigated with X-ray Absorption Spectroscopy (XAS) at the As K-edge (11,867 eV). In two of the studied samples, XANES indicates that As is in the 5+ oxidation state only, and EXAFS analysis reveals a local environment typical of arsenate species. For these samples, the lack of detectable second shell signals suggests a poorly ordered environment, possibly corresponding to an adsorption onto oxide and/or silicate phases. On the other hand, in the third sample XANES reveals a mixed As oxidation state (III and V). This sample also presents evident next nearest neighbor coordination shells, attributed to As-Ca and As-As contributions. The occurrence of next neighbor shells is evidence that part of As is incorporated in an ordered lattice. Furthermore, the local structure revealed by EXAFS is compatible with As incorporation in the calcite phase, as further supported by DFT simulations. The observation of next neighbors shells only in the As(III)-rich sample suggests the substitution of the arsenite group in place of the carbonate one (CO₃²⁻ ⇔ AsO₃³⁻). The conclusion of this work is that uptake of As by calcite is in general less favored than adsorption onto iron oxhydroxides, but could become environmentally important wherever the latter phenomenon is hindered.     

The role of sulfate groups in controlling CaCO3 polymorphism

The nucleation and growth of CaCO₃ phases from aqueous solutions with SO₄²⁻:CO₃²⁻ ratios from 0 to 1.62 and a pH of ∼10.9 were studied experimentally in batch reactors at 25 °C. The mineralogy, morphology and composition of the precipitates were characterized by X-ray diffraction, Fourier transform infrared spectroscopy, scanning electron microscopy and microanalyses. The solids recovered after short reaction times (5 min to 1 h) consisted of a mixture of calcite and vaterite, with a S content that linearly correlates with the SO₄²⁻:CO₃²⁻ ratio in the aqueous solution. The solvent-mediated transformation of vaterite to calcite subsequently occurred. After 24 h of equilibration, calcite was the only phase present in the precipitate for aqueous solutions with SO₄²⁻:CO₃²⁻ ? 1. For SO₄²⁻:CO₃²⁻ > 1, vaterite persisted as a major phase for a longer time (>250 h for SO₄²⁻:CO₃²⁻ = 1.62). To study the role of sulfate in stabilizing vaterite, we performed a molecular simulation of the substitution of sulfate for carbonate groups into the crystal structure of vaterite, aragonite and calcite. The results obtained show that the incorporation of small amounts (<3 mole%) of sulfate is energetically favorable in the vaterite structure, unfavorable in calcite and very unfavorable in aragonite. The computer modeling provided thermodynamic information, which, combined with kinetic arguments, allowed us to put forward a plausible explanation for the observed crystallization behavior.     

Sorption of Eu(III)/Cm(III) on Ca-montmorillonite and Na-illite. Part 1: Batch sorption and time-resolved laser fluorescence spectroscopy experiments

Sorption of Cm(III) and Eu(III) at trace concentrations onto Ca-montmorillonite (SWy-1) and Na-illite (Illite du Puy) has been studied under anaerobic conditions by batch sorption experiments and time-resolved laser fluorescence spectroscopy (TRLFS). Comparison of the results from spectroscopic and batch sorption experiments with Cm and Eu indicates the existence of outer-sphere complexes at pH <4 in the experiments with Na-illite (0.25 g/L solid; 2.5 × 10⁻⁷ mol/L Cm; 0.1 mol/L NaClO₄). In the case of Ca-montmorillonite, (0.25 g/L solid, 2.5 × 10⁻⁷ mol/L Cm or 10⁻⁶ mol/L Eu, 0.066 mol/L Ca(ClO₄)₂), Cm/Eu outer-sphere complexes do not form at significant levels due to the Ca²⁺ competition for the clay mineral cation-exchange sites. TRLFS spectra indicate the formation of inner-sphere surface complexes at pH >5 for both clay minerals. Five H₂O/OH⁻ molecules remain in the first metal ion coordination sphere of the sorbed Eu/Cm. Measured fluorescence lifetimes of sorbed Eu/Cm and peak deconvolution of Cm-spectra are consistent with the formation of surface complexes of the form ≡S-O-Eu/Cm(OH)_x^(2−x)(H₂O)_5−x. At pH ≥ 12 Cm becomes incorporated into a surface precipitate at the Ca-montmorillonite surface presumably composed of Ca(OH)₂ or calcium silicate hydrate. A dramatic shift of the fluorescence emission band by more than 20 nm and a clear increase in the fluorescence lifetime suggests the almost complete displacement of coordinated H₂O and OH⁻. The pH dependent Eu sorption data obtained in batch experiments are consistent with spectroscopic data on Eu and Cm within experimental uncertainties thus demonstrating the validity of Eu as a homologue for trivalent actinides. Parameterization of a two-site protolysis nonelectrostatic surface complexation and cation exchange model using the batch sorption data and spectroscopic results is discussed in Part 2 of this work.     

Goethite, calcite, and organic matter from Permian and Triassic soils: carbon isotopes and CO2 concentrations

Pedogenic goethites in each of two Early Permian paleosols appear to record mixing of two isotopically distinct CO₂ components—atmospheric CO₂ and CO₂ from in situ oxidation of organic matter. The δ¹³C values measured for the Fe(CO₃)OH component in solid solution in these Permian goethites are −13.5‰ for the Lower Leonardian (∼283 Ma BP) paleosol (MCGoeth) and −13.9‰ for the Upper Leonardian (∼270 Ma BP) paleosol (SAP). These goethites contain the most ¹³C-rich Fe(CO₃)OH measured to date for pedogenic goethites crystallized in soils exhibiting mixing of the two aforementioned CO₂ components. δ¹³C measured for 43 organic matter samples in the Lower Leonardian (Waggoner Ranch Fm.) has an average value of −20.3 ± 1.1‰ (1s). The average value yields a calculated Early Permian atmospheric Pco₂ value of about 1 × PAL, but the scatter in the measured δ¹³C values of organic matter permits a calculated maximum Pco₂ of 11 × PAL (PAL = present atmospheric level). Measured values of the mole fraction of Fe(CO₃)OH in MCGoeth and SAP correspond to soil CO₂ concentrations in the Early Permian paleosol profiles of 54,000 and 50,000 ppmV, respectively. Such high soil CO₂ concentrations are similar to modern soils in warm, wet environments.The average δ¹³C values of pedogenic calcite from 9 paleosol profiles stratigraphically associated with MCGoeth (Waggoner Ranch Fm.) range from −6.5‰ to −4.4‰, with a mean δ¹³C value for all profiles of −5.4‰. Thus, the value of Δ¹³C between the pedogenic calcite data set and MCGoeth is 8.1 (±0.9)‰, which is in reasonable accord with the value of 7.7‰ expected if atmospheric Pco₂ and organic matter δ¹³C values were the same for both paleosol types. Furthermore, the atmospheric Pco₂ calculated for the Early Permian from the average measured carbon isotopic compositions of the paleosol calcite and organic matter is also analytically indistinguishable from 1 × PAL, with a maximum calculated atmospheric Pco₂ (permitted by one standard deviation of the organic matter δ¹³C value) of ∼5 × PAL.If, however, measured average δ¹³C values of the plant organic matter are more positive than the original soil organic matter as a result of diagenetic loss of ¹³C-depleted, labile organic compounds, calculated Permian atmospheric Pco₂ using these ¹³C-enriched organic values would underestimate the actual atmospheric Pco₂ using either goethite or calcite. This is the first stratigraphically constrained, intrabasinal study to compare ancient atmospheric CO₂ concentrations calculated from pedogenic goethite and calcite. These results demonstrate that the two different proxies record the same information about atmospheric CO₂.The Fe(CO₃)OH component in pedogenic goethite from a Triassic paleosol in Utah is significantly enriched in ¹³C relative to Fe(CO₃)OH in goethites from soils in which there are mixtures of two isotopic CO₂ components. Field-relationships and the δ¹³C value (−1.9‰) of the Triassic goethite indicate that this ancient paleosol profile experienced mixing of three isotopically distinct CO₂ components at the time of goethite crystallization. The three components were probably atmospheric CO₂, CO₂ from in situ oxidation of organic matter and CO₂ from in situ dissolution of preexisting calcite. Although mixing of three isotopically distinct CO₂ components, as recorded by Fe(CO₃)OH in goethite, has been described in modern soil, this is the first example from a documented paleosol. Its preservation affirms the need for careful, case-by-case assessment of ancient paleosols to establish that goethite in any particular soil is likely to be a valid proxy of atmospheric Pco₂.     

Partition of Sr,Ba, Ca,Y, Eu2+, Eu3+, and other REE between plagioclase feldspar and magmatic liquid: an experimental study

Plagioclase feldspar/magmatic liquid partition coefficients for Sr, Ba, Ca, Y, Eu²⁺, Eu³⁺ and other REE have been determined experimentally at 1 atm total pressure in the temperature range 1150–1400°C. Natural and synthetic melts representative of basaltic and andesitic bulk compositions were used, crystallizing plagioclase feldspar in the composition range An₃₅–An₈₅. Partition coefficients for Sr are greater than unity at all geologically reasonable temperatures, and for Ba are less than unity above approximately 1060°C. Both are strongly dependent upon temperature. Partition coefficients for the trivalent REE are relatively insensitive to temperature. At fixed temperature they decrease monotonically from La to Lu. The partition of Eu is a strong function of oxygen fugacity. Under extreme reducing conditions D_Eu approaches the value of D_Sr.     

History of carbonate ion concentration over the last 100 million years

Instead of having been more or less constant, as once assumed, it is now apparent that the major ion chemistry of the oceans has varied substantially over time. For instance, independent lines of evidence suggest that calcium concentration ([Ca²⁺]) has approximately halved and magnesium concentration ([Mg²⁺]) approximately doubled over the last 100 million years. On the other hand, the calcite compensation depth, and hence the CaCO₃ saturation, has varied little over the last 100 My as documented in deep sea sediments. We combine these pieces of evidence to develop a proxy for seawater carbonate ion concentration ([CO₃²⁻]) over this period of time. From the calcite saturation state (which is proportional to the product of [Ca²⁺] times [CO₃²⁻], but also affected by [Mg²⁺]), we can calculate seawater [CO₃²⁻]. Our results show that [CO₃²⁻] has nearly quadrupled since the Cretaceous. Furthermore, by combining our [CO₃²⁻] proxy with other carbonate system proxies, we provide calculations of the entire seawater carbonate system and atmospheric CO₂. Based on this, reconstructed atmospheric CO₂ is relatively low in the Miocene but high in the Eocene. Finally, we make a strong case that seawater pH has increased over the last 100 My.     

Measurement of the effects of temperature and partial pressure of oxygen on the oxidation states of europium in silicate glasses

The techniques of electron paramagnetic resonance (EPR) were used to measure the concentration ratio of Eu²⁺ to Eu³⁺ in quenched CaMgSi₂O₆, Ca₃Si₃O₁₂, and CaAl₂Si₂O₈ liquids as functions of partial pressure of oxygen and temperature. The redox equilibrium of the Eu ions was described by the reaction 4Eu³⁺ + 20^2? = 4Eu²⁺ + O₂. The reduction of Eu³⁺ to Eu²⁺ was endothermic, and for CaMgSi₂O₆ and Ca₃Al₂Si₃O₁₂ liquids the mean value of ΔH⁰ and the standard deviation from that mean were 25 ± 7 kcal/mole.The magnitude of the Eu anomaly in the distribution coefficients is discussed in terms of the compositions of the solid and liquid phases.     

Interaction between dissolved silica and calcium carbonate: 1. Spontaneous precipitation of calcium carbonate in the presence of dissolved silica

The kinetics of spontaneous precipitation of CaCO₃ from aqueous solution in the presence of dissolved silica was investigated by recording pH as a function of time. The presence of dissolved silica, at concentrations below saturation with respect to the amorphous phase, decreases induction time for CaCO₃ nucleation, but does not affect CaCO₃ polymorphism. For a “pure” system without silica, the surface free energy, σ, determined from classical nucleation theory is 42 mJ m⁻². This agrees well with values reported in the literature for vaterite and indicates some degree of heterogeneous nucleation, which can occur because of the relatively low degree of supersaturation used for the experiments. In the presence of 1 and 2 mM silica, σ is 37 and 34 mJ m⁻², indicating an increasing degree of heterogeneous nucleation as the amount of polymeric silica increases. The ratio of Ca²⁺ to CO₃²⁻ activity was a governing parameter for determining which CaCO₃ polymorph precipitated. At high Ca²⁺ to CO₃²⁻ activity ratios, almost all initial solid was vaterite, whereas at low ratios, a mixture of vaterite and calcite was observed. In solutions with low Ca²⁺ to CO₃²⁻ activity ratios, the presence of silica at concentrations above saturation with respect to amorphous silica led to formation of only calcite and strongly influenced the crystalline structure and morphology of the precipitates. At high Ca²⁺ to CO₃²⁻ ratios, system behaviour did not differ from that without silica.     

Modeling of subsurface calcite dissolution, including the respiration and reoxidation processes of marine sediments in the region of equatorial upwelling off Gabon

Mineralization of organic matter and the subsequent dissolution of calcite were simulated for surface sediments of the upper continental slope off Gabon by using microsensors to measure O₂, pH, pCO₂ and Ca²⁺ (in situ), pore-water concentration profiles of NO₃⁻, NH₄⁺, Fe²⁺, and Mn²⁺ and SO₄²⁻ (ex situ), as well as sulfate reduction rates derived from incubation experiments. The transport and reaction model CoTReM was used to simulate the degradation of organic matter by O₂, NO₃⁻, Fe(OH)₃ and SO₄²⁻, reoxidation reactions involving Fe²⁺ and Mn²⁺, and precipitation of FeS. Model application revealed an overall rate of organic matter mineralization amounting to 50 μmol C cm⁻² yr⁻¹, of which 77% were due to O₂, 17% to NO₃⁻ and 3% to Fe(OH)₃ and 3% to SO₄²⁻. The best fit for the pH profile was achieved by adapting three different dissolution rate constants of calcite ranging between 0.01 and 0.5% d⁻¹ and accounting for different calcite phases in the sediment. A reaction order of 4.5 was assumed in the kinetic rate law. A CaCO₃ flux to the sediment was estimated to occur at a rate of 42 g m⁻² yr⁻¹ in the area of equatorial upwelling. The model predicts a redissolution flux of calcite amounting to 36 g m⁻² yr⁻¹, thus indicating that ∼90% of the calcite flux to the sediment is redissolved.     

Laser Ablation ICPMS study of trace element partitioning between plagioclase and basaltic melts: an experimental approach

Plagioclase-melt partition coefficients (D) for 34 trace elements at natural concentration levels were determined experimentally in a natural MORB composition at atmospheric pressure using thin Pt-wire loops. Experiments were carried out at three temperatures (1,220, 1,200, and 1,180°C), and at three different oxygen fugacities (fO₂ = IW, QFM, air) in order to assess the effect of fO₂ on the partitioning of elements with multiple valence (Fe, Eu, Cr). Run products were analyzed by laser-ablation ICP-MS. Most trace element Ds increase slightly as temperature decreases, except for D _Zr, D _Fe, D _Eu and D _Cr that vary systematically with fO₂. Applying the Lattice Strain Model to our data suggests the presence of Fe²⁺ entirely in the octahedral site at highly to moderate reducing conditions, while Fe³⁺ was assigned wholly to the tetrahedral site of the plagioclase structure. Furthermore, we provide a new quantitative framework for understanding the partitioning behaviour of Eu, which occurs as both 2+ and 3+ cations, depending on fO₂and confirm the greater compatibility of Eu²⁺, which has an ionic radius similar to Sr, relative to Eu³⁺ in plagioclase and the higher Eu²⁺/ Eu³⁺ under reducing conditions. For petrogenetic basaltic processes, a combined fractionation of Eu²⁺–Sr and Fe–Mg by plagioclase has considerable potential as an oxybarometer for natural magmatic rocks.     

Diagenesis and geochemistry of porites corals from Papua New Guinea: Implications for paleoclimate reconstruction

Coral proxy records of sea surface temperature (SST) and hydrological balance have become important tools in the field of tropical paleoclimatology. However, coral aragonite is subject to post-depositional diagenetic alteration in both the marine and vadose environments. To understand the impact of diagenesis on coral climate proxies, two mid-Holocene Porites corals from raised reefs on Muschu Island, Papua New Guinea, were analysed for Sr/Ca, δ¹⁸O, and δ¹³C along transects from 100% aragonite to 100% calcite. Thin-section analysis showed a characteristic vadose zone diagenetic sequence, beginning with leaching of primary aragonite and fine calcite overgrowths, transitional to calcite void filling and neomorphic, fabric selective replacement of the coral skeleton. Average calcite Sr/Ca and δ¹⁸O values were lower than those for coral aragonite, decreasing from 0.0088 to 0.0021 and −5.2 to −8.1‰, respectively. The relatively low Sr/Ca of the secondary calcite reflects the Sr/Ca of dissolving phases and the large difference between aragonite and calcite Sr/Ca partition coefficients. The decrease in δ¹⁸O of calcite relative to coral aragonite is a function of the δ¹⁸O of precipitation. Carbon-isotope ratios in secondary calcite are variable, though generally lower relative to aragonite, ranging from −2.5 to −10.4%. The variability of δ¹³C in secondary calcite reflects the amount of soil CO₂ contributing ¹³C-depleted carbon to the precipitating fluids. Diagenesis has a greater impact on Sr/Ca than on δ¹⁸O; the calcite compositions reported here convert to SST anomalies of 115°C and 14°C, respectively. Based on calcite Sr/Ca compositions in this study and in the literature, the sensitivity of coral Sr/Ca-SST to vadose-zone calcite diagenesis is 1.1 to 1.5°C per percent calcite. In contrast, the rate of change in coral δ¹⁸O-SST is relatively small (−0.2 to 0.2°C per percent calcite). We show that large shifts in δ¹⁸O, reported for mid-Holocene and Last Interglacial corals with warmer than present Sr/Ca-SSTs, cannot be caused by calcite diagenesis. Low-level calcite diagenesis can be detected through X-ray diffraction techniques, thin section analysis, and high spatial resolution sampling of the coral skeleton and thus should not impede the production of accurate coral paleoclimate reconstructions.     

Interdiffusion of H2O and CO2 in metamorphic fluids at ∼490 to 690°C and 1 GPa

Interdiffusion coefficients have been determined for H₂O-CO₂ mixtures by quantifying the flux of CO₂ between two fluid-filled chambers in a specially designed piston-cylinder cell. The two chambers, which are maintained at 1.0 GPa and at temperatures differing by ∼100°C, each contain the X_CO2-buffering assemblage calcite + quartz + wollastonite, in H₂O. The positive dependence of X_CO2 on temperature results in a down-temperature, steady-state flux of CO₂ through a capillary tube that connects the two chambers. This flux drives the wollastonite = calcite + quartz equilibrium to the right in the cooler chamber, producing a measurable amount of calcite that is directly related to CO₂-H₂O interdiffusion rates. Diffusivities calculated from seven experiments range from 1.0 × 10⁻⁸ to 6.1 × 10⁻⁸ m²/s for mean capillary temperatures between ∼490 and 690°C. The data set can be approximated by an Arrhenius-type relation:

     

Rate-controlled calcium isotope fractionation in synthetic calcite

The isotopic composition of Ca (Δ⁴⁴Ca/⁴⁰Ca) in calcite crystals has been determined relative to that in the parent solutions by TIMS using a double spike. Solutions were exposed to an atmosphere of NH₃ and CO₂, provided by the decomposition of (NH₄)₂CO₃, following the procedure developed by previous workers. Alkalinity, pH and concentrations of CO₃²⁻, HCO₃⁻, and CO₂ in solution were determined. The procedures permitted us to determine Δ(⁴⁴Ca/⁴⁰Ca) over a range of pH conditions, with the associated ranges of alkalinity. Two solutions with greatly different Ca concentrations were used, but, in all cases, the condition [Ca²⁺]>>[CO₃²⁻] was met. A wide range in Δ(⁴⁴Ca/⁴⁰Ca) was found for the calcite crystals, extending from 0.04 ± 0.13‰ to −1.34 ± 0.15‰, generally anti-correlating with the amount of Ca removed from the solution. The results show that Δ(⁴⁴Ca/⁴⁰Ca) is a linear function of the saturation state of the solution with respect to calcite (Ω). The two parameters are very well correlated over a wide range in Ω for each solution with a given [Ca]. The linear correlation extended from Δ(⁴⁴Ca/⁴⁰Ca) = −1.34 ± 0.15‰ to 0.04 ± 0.13‰, with the slopes directly dependent on [Ca]. Solutions, which were vigorously stirred, showed a much smaller range in Δ(⁴⁴Ca/⁴⁰Ca) and gave values of −0.42 ± 0.14‰, with the largest effect at low Ω. It is concluded that the diffusive flow of CO₃²⁻ into the immediate neighborhood of the crystal-solution interface is the rate-controlling mechanism and that diffusive transport of Ca²⁺ is not a significant factor. The data are simply explained by the assumptions that: a) the immediate interface of the crystal and the solution is at equilibrium with Δ(⁴⁴Ca/⁴⁰Ca) ∼ −1.5 ± 0.25‰; and b) diffusive inflow of CO₃²⁻ causes supersaturation, thus precipitating Ca from the regions exterior to the narrow zone of equilibrium. The result is that Δ(⁴⁴Ca/⁴⁰Ca) is a monotonically increasing (from negative values to zero) function of Ω. We consider this model to be a plausible explanation of most of the available data reported in the literature. The well-resolved but small and regular isotope fractionation shifts in Ca are thus not related to the diffusion of very large hydrated Ca complexes, but rather due to the ready availability of Ca in the general neighborhood of the crystal-solution interface. The largest isotopic shift which occurs as a small equilibrium effect is then subdued by supersaturation precipitation for solutions where [Ca²⁺]>>[CO₃²⁻] + [HCO₃⁻]. It is shown that there is a clear temperature dependence of the net isotopic shifts that is simply due to changes in Ω due to the equilibrium “constants” dependence on temperature, which changes the degree of saturation and hence the amount of isotopically unequilibrated Ca precipitated. The effects that are found in natural samples, therefore, will be dependent on the degree of diffusive inflow of carbonate species at or around the crystal-liquid interface in the particular precipitating system, thus limiting the equilibrium effect.     

Immobilization of strontium during iron biomineralization coupled to dissimilatory hydrous ferric oxide reduction

The potential for incorporation of strontium (Sr) into biogenic Fe(II)-bearing minerals formed during microbial reduction of synthetic hydrous ferric oxide (HFO) was investigated in circumneutral bicarbonate-buffered medium containing SrCl₂ at concentrations of 10 μM, 100 μM, or 1.0 mM. CaCl₂ (10 mM) was added to some experiments to simulate a Ca-rich groundwater. In Ca-free systems, 89 to 100% of total Sr was captured in solid-phase compounds formed during reduction of 30 to 40 mmol Fe(III) L⁻¹ over a 1-month period. A smaller fraction of total Sr (25 to 34%) was incorporated into the solid phase in cultures amended with 10 mM CaCl₂. X-ray diffraction identified siderite and ferroan ankerite as major end products of HFO reduction in Ca-free and Ca-amended cultures, respectively. Scanning electron microscopy-energy dispersive x-ray spectroscopy revealed the presence of Sr associated with carbonate phases. Selective extraction of HFO reduction end products indicated that 46 to 100% of the solid-phase Sr was associated with carbonates. The sequestration of Sr into carbonate phases in the Ca-free systems occurred systematically according to a heterogeneous (Doerner-Hoskins) partition coefficient (D_D-H) of 1.81 ± 0.15. This D_D-H value was 2 to 10 times higher than values determined for incorporation of Sr (10 μM) into FeCO₃(s) precipitated abiotically at rates comparable to or greater than rates observed during HFO reduction, and fivefold higher than theoretical partition coefficients for equilibrium Fe(Sr)CO₃ solid solution formation. Surface complexation and entrapment of Sr by rapidly growing siderite crystals (and possibly other biogenic Fe(II) solids) provides an explanation for the intensive scavenging of Sr in the Ca-free systems. The results of abiotic siderite precipitation experiments in the presence and absence of excess Ca indicate that substitution of Ca for Sr at foreign element incorporation sites (mass action effect) on growing FeCO₃(s) surfaces can account for the inhibition of Sr incorporation into the siderite component of ankerite formed in the Ca-amended HFO reduction experiments. Likewise, substitution of Fe(II) for Sr may explain the absence of major Sr partitioning into the calcite component of ankerite. The findings indicate that under appropriate conditions, sequestration of metals in siderite produced during bacterial Fe(III) oxide reduction may provide a mechanism for retarding the migration of Sr and other divalent metal contaminants in anaerobic, carbonate-rich sedimentary environments.     

REE constraints on fractionation processes of massive-type anorthosites on the Lofoten Islands, Norway

Summary Rare Earth Element (REE) data of 34 samples of magmatic rocks from the Lofoten Islands in Norway lend support to the derivation of anorthosites, ferrodiorites and jotunites by fractionation and cumulus processes from typical basaltic magma. Both REE concentration and Eu anomalies (expressed as Eu/Eu^*) form continuous linear trends from anorthosite towards gabbro, ferrodiorite and jotunite in discrimination diagrams against molar CaO/Al₂O₃ ratios indicating the predominant accumulation of plagioclase. Eu/Eu^* decreases from about 4 in the cumulates (anorthosites) to around 1 in the fine-grained gabbroic dikes and to below 1 in some ferrodiorites and the jotunite. The various types of ferrodiorites and the jotunite are regarded as residual liquids, in some cases with variable amounts of cumulus plagioclase. The whole fractionation series from gabbro towards anorthosites and ferrodiorites can be observed in a single intrusion. With increasing fractionation, the REE patterns generally change from flat, slightly LREE-enriched or LREE-depleted to steep and strongly LREE-enriched. These changes and the REE abundances are mainly controlled by the abundance of apatite. Temporally and spatially related mangerites and charnockites form a trend from low-SiO₂ mangerites with Eu/Eu^* > 1 to intermediate-SiO₂ acidic mangerites with Eu/Eu^* ≈ 1 and charnockites with Eu/Eu^* < 1. Accordingly, the low-SiO₂ mangerites are interpreted as alkali feldspar-rich cumulates and the charnockites as residual liquids derived from the acidic mangerites. The mangerites with Eu/Eu^* around 1 have patterns similar to those of some highly evolved ferrodiorites possibly indicating a genetic link. Received December 12, 1999; revised version accepted November 15, 2000     

Stable isotope variations in modern tropical speleothems: Evaluating equilibrium vs. kinetic isotope effects

Applications of speleothem calcite geochemistry in climate change studies require the evaluation of the accuracy and sensitivity of speleothem proxies to correctly infer paleoclimatic information. The present study of Harrison’s Cave, Barbados, uses the analysis of the modern climatology and groundwater system to evaluate controls on the C and O isotopic composition of modern speleothems. This new approach directly compares the δ¹⁸O and δ¹³C values of modern speleothems with the values for their corresponding drip waters in order to assess the degree to which isotopic equilibrium is achieved during calcite precipitation. If modern speleothems can be demonstrated to precipitate in isotopic equilibrium, then ancient speleothems, suitable for paleoclimatic studies, from the same cave environment may also have been precipitated in isotopic equilibrium. If modern speleothems are precipitated out of isotopic equilibrium, then the magnitude and direction of the C and O isotopic offsets may allow specific kinetic and/or equilibrium isotopic fractionation mechanisms to be identified.Carbon isotope values for the majority of modern speleothem samples from Harrison’s Cave fall within the range of equilibrium values predicted from the combined use of (1) calcite-water fractionation factors from the literature, (2) measured temperatures, and (3) measured δ¹³C values of the dissolved inorganic carbon of drip waters. Calcite samples range from ∼0.8‰ higher to ∼1.1‰ lower than predicted values. The ¹³C depletions are likely caused by kinetically driven departures in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions, caused by rapid calcite growth. ¹³C enrichments can be accounted for by Rayleigh distillation of the HCO₃⁻ (aq) reservoir during degassing of ¹³C-depleted CO₂.Modern speleothems from Harrison’s Cave are not in O isotopic equilibrium with their corresponding drip waters and are 0.2‰ to 2.3‰ enriched in ¹⁸O relative to equilibrium values. δ¹⁸O variations in modern calcite are likely controlled by kinetically driven changes in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions to nonequilibrium conditions, consistent with rapid calcite growth. In contrast to δ¹³C, δ¹⁸O values of modern calcite may not be affected by Rayleigh distillation during degassing because CO₂ hydration and hydroxylation reactions will buffer the O isotopic composition of the HCO₃⁻ (aq) reservoir. If the effects of Rayleigh distillation manifest themselves in the O isotopic system, they will result in ¹⁸O enrichment in the HCO₃⁻ (aq) reservoir and ultimately in the precipitated CaCO₃.