Comment on the studies of oxygen isotope fractionation between calcium carbonates and water at low temperatures by Zhou and Zheng (2003; 2005)

Experimental and theoretical aspects of oxygen isotope fractionation in the system calcite-water at low temperatures were critically examined. Contrary to the claim made by Zhou and Zheng [Zhou G.-T., and Zheng Y.-F. (2003) An experimental study of oxygen isotope fractionation between inorganically precipitated aragonite and water at low temperatures. Geochim. Cosmochim. Acta67, 387-399], there is excellent agreement between fractionation factors that were experimentally determined by means of slow, inorganic precipitation of calcite from solutions and those obtained largely from theoretical, statistical-mechanical calculations of the reduced partition function ratios. This agreement strongly suggests that calcite was precipitated from a solution very close to isotopic equilibrium. However, recently Zhou and Zheng [Zhou G.-T., and Zheng Y.-F. (2005) Effect of polymorphic transition on oxygen isotope fractionation between aragonite, calcite and water: a low-temperature experimental study. Am. Miner.90, 1121-1130] presented, without any explanation, conclusions on these major aspects that contradict the previous statements of Zhou and Zheng (2003). The apparent discrepancy in calcite-water oxygen isotope fractionation between experimental and theoretical studies discussed by Zhou and Zheng (2003) originates from the “mineral-water interaction” term in the modified increment method, which was developed by one of the authors (Y.-F. Zheng). We call for evidence for the theoretical nature of the modified increment method, which has never been presented in any of Zheng’s papers. Without such evidence, great caution must be exercised in using fractionation factors derived from the modified increment method.     

Oxygen isotope salt effects at high pressure and high temperature and the calibration of oxygen isotope geothermometers

The influence of NaCl, CaCl₂, and dissolved minerals on the oxygen isotope fractionation in mineral-water systems at high pressure and high temperature was studied experimentally. The salt effects of NaCl (up to 37 molal) and 5-molal CaCl₂ on the oxygen isotope fractionation between quartz and water and between calcite and water were measured at 5 and 15 kbar at temperatures from 300 to 750°C. CaCl₂ has a larger influence than NaCl on the isotopic fractionation between quartz and water. Although NaCl systematically changes the isotopic fractionation between quartz and water, it has no influence on the isotopic fractionation between calcite and water. This difference in the apparent oxygen isotope salt effects of NaCl must relate to the use of different minerals as reference phases. The term oxygen isotope salt effect is expanded here to encompass the effects of dissolved minerals on the fractionations between minerals and aqueous fluids. The oxygen isotope salt effects of dissolved quartz, calcite, and phlogopite at 15 kbar and 750°C were measured in the three-phase systems quartz-calcite-water and phlogopite-calcite-water. Under these conditions, the oxygen isotope salt effects of the three dissolved minerals range from ∼0.7 to 2.1‰. In both three-phase hydrothermal systems, the equilibrium fractionation factors between the pairs of minerals are the same as those obtained by anhydrous direct exchange between each pair of minerals, proving that the use of carbonate as exchange medium provides correct isotopic fractionations for a mineral pair.When the oxygen isotope salt effects of two minerals are different, the use of water as an indirect exchange medium will give erroneous fractionations between the two minerals. The isotope salt effect of a dissolved mineral is also the main reason for the observation that the experimentally calibrated oxygen isotope fractionations between a mineral and water are systematically 1.5 to 2‰ more positive than the results of theoretical calculations. Dissolved minerals greatly affect the isotopic fractionation in mineral-water systems at high pressure and high temperature. If the presence of a solute changes the solubility of a mineral, the real oxygen isotope salt effect of the solute at high pressure and high temperature cannot be correctly derived by using the mineral as reference phase.     

Water and magmas: clarification of a controversy on applications of the Gibbs-Duhem equation

It is shown mathematically that if the activity coefficient of water in ternary water-magma (aluminosilicate) systems is constant or varies only with the mole fraction of water, it is not necessary that the binary magmas form ideal solutions contrary to the claims by Burnham et al. (1978, Geochim. Cosmochim. Acta42, 275–276). A molecular viewpoint is presented to support this argument. The properties of analytical equations capable of representing the activity coefficients of usual and unusual systems are discussed. The correct form of the Gibbs-Dunhem equation for dissociative dissolution processes is presented to disprove the claims by Burnham (1975, Fortschr. Mineral.52, 101–118; 1975, Geochim. Cosmochim. Acta39, 1077–1084), and by Burnhamet al. (1978, Geochim. Cosmochim. Acta42, 275–276).     

Sulfur and oxygen isotope study of sulfate reduction in experiments with natural populations from Fællestrand, Denmark

This study investigates the sulfur and oxygen isotope fractionations of dissimilatory sulfate reduction and works to reconcile the relationships between the oxygen and sulfur isotopic and elemental systems. We report results of experiments with natural populations of sulfate-reducing bacteria using sediment and seawater from a marine lagoon at Fællestrand on the northern shore of the island of Fyn, Denmark. The experiments yielded relatively large magnitude sulfur isotope fractionations for dissimilatory sulfate reduction (up to approximately 45‰ for ³⁴S/³²S) with higher δ¹⁸O accompanying higher δ³⁴S, similar to that observed in previous studies. The seawater used in the experiments was spiked by addition of ¹⁷O-labeled water and the ¹⁷O content of residual sulfate was found to depend on the fraction of sulfate reduced in the experiments. The ¹⁷O data provides evidence for recycling of sulfur from metabolic intermediates and for an ¹⁸O/¹⁶O fractionation of ∼25-30‰ for dissimilatory sulfate reduction. The close correlation between the ¹⁷O data and the sulfur isotope data suggests that isotopic exchange between cell water and external water (reactor water) was rapid under experimental conditions. The molar ratio of oxygen exchange to sulfate reduction was found to be about 2.5. This value is slightly lower than observed in studies of natural ecosystems [e.g., Wortmann U. G., Chernyavsky B., Bernasconi S. M., Brunner B., Böttcher M. E. and Swart P. K. (2007) Oxygen isotope biogeochemistry of pore water sulfate in the deep biosphere: dominance of isotope exchange reactions with ambient water during microbial sulfate reduction (ODP Site 1130). Geochim. Cosmochim. Acta71, 4221-4232]. Using recent models of sulfur isotope fractionations we find that our combined sulfur and oxygen isotopic data places constraints on the proportion of sulfate recycled to the medium (78-96%), the proportion of sulfur intermediate sulfite that was recycled by way of APS to sulfate and released back to the external sulfate pool (∼70%), and also that a fraction of the sulfur intermediates between sulfite and sulfide were recycled to sulfate. These parameters can be constrained because of the independent information provided by δ¹⁸O, δ³⁴S, δ¹⁷O labels, and Δ³³S.     

Element partitioning during precipitation of aragonite from seawater: A framework for understanding paleoproxies

This study presents the results from precipitation experiments carried out to investigate the partitioning of the alkaline earth cations Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺ between abiogenic aragonite and seawater as a function of temperature. Experiments were carried out at 5 to 75 °C, using the protocol of Kinsman and Holland [Kinsman, D.J.J., Holland, H.D., 1969. The coprecipitation of cations with CaCO₃ IV. The coprecipitation of Sr²⁺ with aragonite between 16 and 96 °C. Geochim. Cosmochim. Acta33, 1-17.] The concentrations of Mg Sr and Ba were determined in the fluid from each experiment by inductively coupled plasma-mass spectrometry, and in individual aragonite grains by secondary ion mass spectrometry. The experimentally produced aragonite grains are enriched in trace components (“impurities”) relative to the concentrations expected from crystal-fluid equilibrium, indicating that kinetic processes are controlling element distribution. Our data are not consistent with fractionations produced kinetically in a boundary layer adjacent to the growing crystal because Sr²⁺, a compatible element, is enriched rather than depleted in the aragonite. Element compatibilities, and the systematic change in partitioning with temperature, can be explained by the process of surface entrapment proposed by Watson and Liang [Watson, E.B., Liang, Y., 1995. A simple model for sector zoning in slowly grown crystals: implications for growth rate and lattice diffusion, with emphasis on accessory minerals in crustal rocks. Am. Mineral.80, 1179-1187] and Watson [Watson, E.B., 1996. Surface enrichment and trace-element uptake during crystal growth. Geochim. Cosmochim. Acta60, 5013-5020; Watson, E.B., 2004. A conceptual model for near-surface kinetic controls on the trace-element and stable isotope composition of abiogenic calcite crystals. Geochim. Cosmochim. Acta68, 1473-1488]. This process is thought to operate in regimes where the competition between crystal growth rate and diffusivity in the near-surface region limits the extent to which the solid can achieve partitioning equilibrium with the fluid. A comparison of the skeletal composition of Diploria labyrinthiformis (brain coral) collected on Bermuda with results from precipitation calculations carried out using our experimentally determined partition coefficients indicate that the fluid from which coral skeleton precipitates has a Sr/Ca ratio comparable to that of seawater, but is depleted in Mg and Ba, and that there are seasonal fluctuations in the mass fraction of aragonite precipitated from the calcifying fluid (“precipitation efficiency”). The combined effects of surface entrapment during aragonite growth and seasonal fluctuations in “precipitation efficiency” likely forms the basis for the temperature information recorded in the aragonite skeletons of Scleractinian corals.     

Iron isotope fractionation between pyrite (FeS2), hematite (Fe2O3) and siderite (FeCO3): A first-principles density functional theory study

In addition to equilibrium isotopic fractionation factors experimentally derived, theoretical predictions are needed for interpreting isotopic compositions measured on natural samples because they allow exploring more easily a broader range of temperature and composition. For iron isotopes, only aqueous species were studied by first-principles methods and the combination of these data with those obtained by different methods for minerals leads to discrepancies between theoretical and experimental isotopic fractionation factors. In this paper, equilibrium iron isotope fractionation factors for the common minerals pyrite, hematite, and siderite were determined as a function of temperature, using first-principles methods based on the density functional theory (DFT). In these minerals belonging to the sulfide, oxide and carbonate class, iron is present under two different oxidation states and is involved in contrasted types of interatomic bonds. Equilibrium fractionation factors calculated between hematite and siderite compare well with the one estimated from experimental data (ln α⁵⁷Fe/⁵⁴Fe = 4.59 ± 0.30‰ and 5.46 ± 0.63‰ at 20 °C for theoretical and experimental data, respectively) while those for Fe(III)_aq-hematite and Fe(II)_aq-siderite are significantly higher that experimental values. This suggests that the absolute values of the reduced partition functions (β-factors) of aqueous species are not accurate enough to be combined with those calculated for minerals. When compared to previous predictions derived from Mössbauer or INRXS data [Polyakov V. B., Clayton R. N., Horita J. and Mineev S. D. (2007) Equilibrium iron isotope fractionation factors of minerals: reevaluation from the data of nuclear inelastic resonant X-ray scattering and Mössbauer spectroscopy. Geochim. Cosmochim. Acta71, 3833-3846], our iron β-factors are in good agreement for siderite and hematite while a discrepancy is observed for pyrite. However, the detailed investigation of the structural, electronic and vibrational properties of pyrite as well as the study of sulfur isotope fractionation between pyrite and two other sulfides (sphalerite and galena) indicate that DFT-derived β-factors of pyrite are as accurate as for hematite and siderite. We thus suggest that experimental vibrational density of states of pyrite should be re-examined.     

New high-precision measurements of the isotopic composition of atmospheric argon

Earth’s atmosphere is used as a standard reference gas for mass spectrometric determinations of argon (Ar) isotopes used principally in geochronology. There are three published independent determinations of the Ar isotope composition of modern atmosphere that differ subtly. We have made new high-precision measurements of Ar isotope ratios of five different sources of air using a high-sensitivity multi-collector noble gas mass spectrometer in order to distinguish between them. The isotope ratios, corrected only for backgrounds, reside on a inverse square-root mass law fractionation line that passes through the air value proposed by Lee et al. [Lee J. Y., Marti K., Severinghaus J. P., Kawamura K., Yoo H. S., Lee J. B. and Kim J. S. (2006) A redetermination of the isotopic abundances of atmospheric Ar. Geochim. Cosmochim. Acta70, 4507-4512]. They are distinct from both the other proposed compositions and provide the first independent confirmation of the atmospheric Ar isotope composition. We suggest that the revised values should now be in routine employment.     

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

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

The role of reaction affinity and secondary minerals in regulating chemical weathering rates at the Santa Cruz Soil Chronosequence, California

In order to explore the reasons for the apparent discrepancy between laboratory and field weathering rates and to determine the extent to which weathering rates are controlled by the approach to thermodynamic equilibrium, secondary mineral precipitation, and flow rates, a multicomponent reactive transport model (CrunchFlow) was used to interpret soil profile development and mineral precipitation and dissolution rates at the 226 ka Marine Terrace Chronosequence near Santa Cruz, CA. Aqueous compositions, fluid chemistry, transport, and mineral abundances are well characterized [White A. F., Schulz M. S., Vivit D. V., Blum A., Stonestrom D. A. and Anderson S. P. (2008) Chemical weathering of a Marine Terrace Chronosequence, Santa Cruz, California. I: interpreting the long-term controls on chemical weathering based on spatial and temporal element and mineral distributions. Geochim. Cosmochim. Acta72 (1), 36-68] and were used to constrain the reaction rates for the weathering and precipitating minerals in the reactive transport modeling. When primary mineral weathering rates are calculated with either of two experimentally determined rate constants, the nonlinear, parallel rate law formulation of Hellmann and Tisserand [Hellmann R. and Tisserand D. (2006) Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar. Geochim. Cosmochim. Acta70 (2), 364-383] or the aluminum inhibition model proposed by Oelkers et al. [Oelkers E. H., Schott J. and Devidal J. L. (1994) The effect of aluminum, pH, and chemical affinity on the rates of aluminosilicate dissolution reactions. Geochim. Cosmochim. Acta58 (9), 2011-2024], modeling results are consistent with field-scale observations when independently constrained clay precipitation rates are accounted for. Experimental and field rates, therefore, can be reconciled at the Santa Cruz site.Additionally, observed maximum clay abundances in the argillic horizons occur at the depth and time where the reaction fronts of the primary minerals overlap. The modeling indicates that the argillic horizon at Santa Cruz can be explained almost entirely by weathering of primary minerals and in situ clay precipitation accompanied by undersaturation of kaolinite at the top of the profile. The rate constant for kaolinite precipitation was also determined based on model simulations of mineral abundances and dissolved Al, SiO₂(aq) and pH in pore waters. Changes in the rate of kaolinite precipitation or the flow rate do not affect the gradient of the primary mineral weathering profiles, but instead control the rate of propagation of the primary mineral weathering fronts and thus total mass removed from the weathering profile. Our analysis suggests that secondary clay precipitation is as important as aqueous transport in governing the amount of dissolution that occurs within a profile because clay minerals exert a strong control over the reaction affinity of the dissolving primary minerals. The modeling also indicates that the weathering advance rate and the total mass of mineral dissolved is controlled by the thermodynamic saturation of the primary dissolving phases plagioclase and K-feldspar, as is evident from the difference in propagation rates of the reaction fronts for the two minerals despite their very similar kinetic rate laws.     

Silicon isotopic fractionation of CAI-like vacuum evaporation residues

Calcium-, aluminum-rich inclusions (CAIs) are often enriched in the heavy isotopes of magnesium and silicon relative to bulk solar system materials. It is likely that these isotopic enrichments resulted from evaporative mass loss of magnesium and silicon from early solar system condensates while they were molten during one or more high-temperature reheating events. Quantitative interpretation of these enrichments requires laboratory determinations of the evaporation kinetics and associated isotopic fractionation effects for these elements. The experimental data for the kinetics of evaporation of magnesium and silicon and the evaporative isotopic fractionation of magnesium is reasonably complete for Type B CAI liquids (Richter F. M., Davis A. M., Ebel D. S., and Hashimoto A. (2002) Elemental and isotopic fractionation of Type B CAIs: experiments, theoretical considerations, and constraints on their thermal evolution. Geochim. Cosmochim. Acta66, 521-540; Richter F. M., Janney P. E., Mendybaev R. A., Davis A. M., and Wadhwa M. (2007a) Elemental and isotopic fractionation of Type B CAI-like liquids by evaporation. Geochim. Cosmochim. Acta71, 5544-5564.). However, the isotopic fractionation factor for silicon evaporating from such liquids has not been as extensively studied. Here we report new ion microprobe silicon isotopic measurements of residual glass from partial evaporation of Type B CAI liquids into vacuum. The silicon isotopic fractionation is reported as a kinetic fractionation factor, α_Si, corresponding to the ratio of the silicon isotopic composition of the evaporation flux to that of the residual silicate liquid. For CAI-like melts, we find that α_Si = 0.98985 ± 0.00044 (2σ) for ²⁹Si/²⁸Si with no resolvable variation with temperature over the temperature range of the experiments, 1600-1900 °C. This value is different from what has been reported for evaporation of liquid Mg₂SiO₄ (Davis A. M., Hashimoto A., Clayton R. N., and Mayeda T. K. (1990) Isotope mass fractionation during evaporation of Mg₂SiO₄. Nature347, 655-658.) and of a melt with CI chondritic proportions of the major elements (Wang J., Davis A. M., Clayton R. N., Mayeda T. K., and Hashimoto A. (2001) Chemical and isotopic fractionation during the evaporation of the FeO-MgO-SiO₂-CaO-Al₂O₃-TiO₂-REE melt system. Geochim. Cosmochim. Acta65, 479-494.). There appears to be some compositional control on α_Si, whereas no compositional effects have been reported for α_Mg. We use the values of α_Si and α_Mg, to calculate the chemical compositions of the unevaporated precursors of a number of isotopically fractionated CAIs from CV chondrites whose chemical compositions and magnesium and silicon isotopic compositions have been previously measured.     

Modelling δC and δO in the solution layer on stalagmite surfaces

We derive equations describing the evolution of the carbon and oxygen isotope composition of the bicarbonate in a calcite precipitating solution on the surface of a stalagmite using a classical Rayleigh approach. The combined effects of calcite precipitation, degassing of CO₂ and the buffering effect of the water reservoir are taken into account. Whereas δ¹³C shows a progressive increase to a final constant value, δ¹⁸O shows an initial isotopic enrichment, which exponentially decays due to the buffering effect of the water reservoir. The calculated evolution is significantly different for both carbon and oxygen isotopes than derived in a recent paper [Dreybrodt W. (2008) Evolution of the isotopic composition of carbon and oxygen in a calcite precipitating H₂O-CO₂-CaCO₃ solution and the related isotopic composition of calcite in stalagmites. Geochim. Cosmochim. Acta72, 4712-4724.].Furthermore, we discuss the isotopic evolution of the bicarbonate in the solution for long residence times on the stalagmite surface, i.e., for t→∞. The equilibrium isotope ratio of the bicarbonate is then determined by isotopic exchange between the cave atmosphere and the bicarbonate in the solution and can be calculated by equilibrium isotope fractionation. For strongly ventilated caves exchange with the cave atmosphere will result in higher δ¹³C and δ¹⁸O values than those observed in a pure Rayleigh distillation scenario, for sparsely ventilated caves it will result in lower δ¹³C and δ¹⁸O values.     

Comment on “Hydrolysis of neptunium(V) at variable temperatures (10-85 °C)” by L. Rao, T.G. Srinivasan, A.Yu. Garnov, P. Zanonato, P. Di Bernardo, and A. Bismondo

In a recent study [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830.] the hydrolysis of Np(V) was investigated at 10-85 °C by absorption spectroscopy, potentiometry, and microcalorimetry along the titration of Np(V) solutions with tetramethylammonium hydroxide up to pH 10. However, there is strong evidence that the precautions to avoid competing reactions with carbonate were not sufficient and that the measured effects are not caused by the formation of Np(V) hydroxide complexes but primarily by the formation of Np(V) carbonate complexes. The reported equilibrium constants, enthalpies, entropies, and heat capacities for the complexes NpO₂OH(aq) and are severely in error and must not be used for the geochemical modeling of neptunium. If the hydrolysis constants reported by Rao et al. [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830] are used to calculate neptunium solubilities in alkaline solutions relevant for nuclear waste repositories, the Np(V) concentrations are overestimated by orders of magnitude.     

Response to the comment by Z.F. Zhao and T.S. Gao (2007) on “Paleozoic ages and excess Ar in garnets from the Bixiling eclogite in dabieshan, China: New insights from Ar/Ar dating by stepwise crushing”

In our paper we supposed that the Paleozoic ages of the Bixiling eclogitic garnets by ⁴⁰Ar/³⁹Ar crushing in vacuo were related to the UHP metamorphism. Zhao and Gao [Zhao Z.-F. and Gao T.-S. (2007) Comment on “Paleozoic ages and excess ⁴⁰Ar in garnets from the Bixiling eclogite in Dabieshan, China:Newinsights from ⁴⁰Ar/³⁹Ar dating by stepwise crushing” by Qiu and Wijbrans (2006). Geochim. Cosmochim. Acta71(24), 6046-6050] dispute our conclusions and re-interpret them as “the age of garnet growth by low-T/HP blueschist/eclogite-facies metamorphism of the UHP eclogite precursor during arc-continent collision in the early Paleozoic” without presenting any new petrological evidence. Based on our present understanding of ⁴⁰Ar/³⁹Ar crushing technique and petrological observations, we do not agree with their re-interpretation and still prefer our original interpretation.     

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.     

Cooling rate based on schreibersite growth for the Emery mesosiderite

The Emery mesosiderite contains large Ni-rich grains of the phosphide schreibersite, which have exsolved from kamacite. Computer simulation of diffusion-controlled growth of this schreibersite indicates that exsolution occurred during cooling at the rate of 0.1°C/Myr. This determination supports the cooling rate estimated for mesosiderites by Powell (Geochim. Cosmochim. Acta33, 789–810), using taenite-kamacite data.     

Meteoritic component and impact melt composition at the lac à l'Eau Claire (clearwater) impact structures,Quebec

The major and selected trace element composition of the melt at the eastern Lac à l'Eau Claire impact structure can be modelled by 52.3% granodiorite, 15.1% granulite, 25.7% quartz-monzonite plus 6.9% Orgueil (C1-chondrite). A satisfactory fit to melt at the western structure requires 57.4% granulite and 42.6% quartz-monzonite. These models agree with the suggestions of Palmeet al. (1978, Geochim. Cosmochim. Acta42, 313–323) on the type and extent of meteoritic component in the melts.     

Kinetic isotope effect during reduction of iron from a silicate melt

Iron isotopic compositions measured in chondrules from various chondrites vary between δ⁵⁷Fe/⁵⁴Fe = +0.9‰ and −2.0‰, a larger range than for igneous rocks. Whether these compositions were inherited from chondrule precursors, resulted from the chondrule-forming process itself or were produced by later parent body alteration is as yet unclear. Since iron metal is a common phase in some chondrules, it is important to explore a possible link between the metal formation process and the observed iron isotope mass fractionation. In this experimental study we have heated a fayalite-rich composition under reducing conditions for heating times ranging from 2 min to 6 h. We performed chemical and iron isotope analyses of the product phases, iron metal and silicate glass. We demonstrated a lack of evaporation of Fe from the silicate melt in similar isothermal experiments performed under non-reducing conditions. Therefore, the measured isotopic mass fractionation in the glass, ranging between −0.32‰ and +3.0‰, is attributed to the reduction process. It is explained by the faster transport of lighter iron isotopes to the surface where reduction occurs, and is analogous to kinetic isotope fractionation observed in diffusion couples [Richter, F.M., Davis, A.M., Depaolo, D.J., Watson, E.B., 2003. Isotope fractionation by chemical diffusion between molten basalt and rhyolite. Geochim. Cosmochim. Acta67, 3905-3923]. The metal phase contains 90-99.8% of the Fe in the system and lacks significant isotopic mass fractionation, with values remaining similar to that of the starting material throughout. The maximum iron isotope mass fractionation in the glass was achieved within 1 h and was followed by an isotopic exchange and re-equilibration with the metal phase (incomplete at ∼6 h). This study demonstrates that reduction of silicates at high temperatures can trigger iron isotopic fractionation comparable in its bulk range to that observed in chondrules. Furthermore, if metal in Type I chondrules was formed by reduction of Fe silicate, our observed isotopic fractionations constrain chondrule formation times to approximately 60 min, consistent with previous work.     

New data on the kinetics of the dissolution of alkali feldspars at 200°C in CO2 charged water

The kinetics of the dissolution of a (Na. K) sanidine in CO₂-charged water, were studied experimentally at 200°C as a function of the surface area s and of the time t. The molalities of Na⁺ and K⁺ have been plotted against the product st. As for albite and adularia (Lagache, Bull. Soc. Fr. Minéral. Cristallogr. 88, 223–253, 1965), the rates of dissolution of Na² and K⁺ are smooth continuous functions of st, which implies that they are controlled by the composition of the solution.The comparison between the dissolution of pure sodic or potassic feldspars and that of an intermediate feldspar shows that the sanidine dissolves as if it were composed of albite and K-feldspar grains in the proportions corresponding to its composition.Theoretical considerations presented by Helgeson (Geochim. Cosmochim. Acta35, 421–169, 1971; The Feldspars, pp. 184–217, 1972) and Pa?es (Geochim. Cosmochim. Acta37, 2641–2663, 1973) are reviewed: both had suggested, arguing from my first experimental results, that the dissolution could be described by a process of diffusional mass transfer through a surface layer of reaction products.The present experiments do not agree with such an interpretation of the mechanism of dissolution.     

GEOCARBSULF: A combined model for Phanerozoic atmospheric O2 and CO2

A model for the combined long-term cycles of carbon and sulfur has been constructed which combines all the factors modifying weathering and degassing of the GEOCARB III model [Berner R.A., Kothavala Z., 2001. GEOCARB III: a revised model of atmospheric CO₂ over Phanerozoic time. Am. J. Sci. 301, 182-204] for CO₂ with rapid recycling and oxygen dependent carbon and sulfur isotope fractionation of an isotope mass balance model for O₂ [Berner R.A., 2001. Modeling atmospheric O₂ over Phanerozoic time. Geochim. Cosmochim. Acta65, 685-694]. New isotopic data for both carbon and sulfur are used and new feedbacks are created by combining the models. Sensitivity analysis is done by determining (1) the effect on weathering rates of using rapid recycling (rapid recycling treats carbon and sulfur weathering in terms of young rapidly weathering rocks and older more slowly weathering rocks); (2) the effect on O₂ of using different initial starting conditions; (3) the effect on O₂ of using different data for carbon isotope fractionation during photosynthesis and alternative values of oceanic δ¹³C for the past 200 million years; (4) the effect on sulfur isotope fractionation and on O₂ of varying the size of O₂ feedback during sedimentary pyrite formation; (5) the effect on O₂ of varying the dependence of organic matter and pyrite weathering on tectonic uplift plus erosion, and the degree of exposure of coastal lands by sea level change; (6) the effect on CO₂ of adding the variability of volcanic rock weathering over time [Berner, R.A., 2006. Inclusion of the weathering of volcanic rocks in the GEOCARBSULF model. Am. J. Sci.306 (in press)]. Results show a similar trend of atmospheric CO₂ over the Phanerozoic to the results of GEOCARB III, but with some differences during the early Paleozoic and, for variable volcanic rock weathering, lower CO₂ values during the Mesozoic. Atmospheric oxygen shows a major broad late Paleozoic peak with a maximum value of about 30% O₂ in the Permian, a secondary less-broad peak centered near the Silurian/Devonian boundary, variation between 15% and 20% O₂ during the Cambrian and Ordovician, a very sharp drop from 30% to 15% O₂ at the Permo-Triassic boundary, and a more-or less continuous rise in O₂ from the late Triassic to the present.     

Europium retention onto clay minerals from 25 to 150 °C: Experimental measurements, spectroscopic features and sorption modelling

The sorption of Eu(III) onto kaolinite and montmorillonite was investigated up to 150 °C. The clays were purified samples, saturated with Na in the case of montmorillonite. Batch experiments were conducted at 25, 40, 80 and 150 °C in 0.5 M NaClO₄ solutions to measure the distribution coefficients (Kd) of Eu as a trace element (<10⁻⁶ mol/L) between the solution and kaolinite. For the Na-montmorillonite, we used Kd results from a previous study [Tertre, E., Berger, G., Castet, S., Loubet, M., Giffaut, E., 2005. Experimental study of adsorption of Ni²⁺, Cs⁺ and Ln³⁺ onto Na-montmorillonite up to 150 °C. Geochim. Cosmochim. Acta69, 4937-4948] obtained under exactly the same conditions. The number and nature of the Eu species sorbed onto both clay minerals were investigated by time resolved laser fluorescence spectroscopy (TRLFS) in specific experiments in the same temperature range. We identified a unique inner-sphere complex linked to the aluminol sites in both clays, assumed to be AlOEu²⁺ at the edge of the particles, and a second exchangeable outer-sphere complex for montmorillonite, probably in an interlayer position. The Kd values were used to adjust the parameters of a surface complexation model (DLM: diffuse layer model) from 25 to 150 °C. The number of Eu complexes and the stoichiometry of reactions were constrained by TRLFS. The acidity constants of the amphoteric aluminol sites were taken from another study [Tertre, E., Castet, S., Berger, G., Loubet, M., Giffaut, E. Acid/base surface chemistry of kaolinite and Na-montmorillonite at 25 and 60 °C: experimental study and modelling. Geochim. Cosmochim. Acta, in press], which integrates the influence of the negative structural charge of clays on the acid/base properties of edge sites as a function of temperature and ionic strength. The results of the modelling show that the observed shift of the sorption edge towards low pH with increasing temperature results solely from the contribution of the AlOEu²⁺ edge complexes. Finally, we successfully tested the performance of our model by confronting the predictions with experimental Kd data. We used our own data obtained at lower ionic strength (previous study) or higher suspension density and higher starting concentration (TRLFS runs, this study), as well as published data from other experimental studies [Bradbury, M.H., Baeyens, B., 2002. Sorption of Eu on Na and Ca-montmorillonite: experimental investigations and modeling with cation exchange and surface complexation. Geochim. Cosmochim. Acta66, 2325-2334; Kowal-Fouchard, A., 2002. Etude des mécanismes de rétention des ions U(IV) et Eu(III) sur les argiles: influence des silicates. Ph.D. Thesis, Université Paris Sud, France, 330p].