Oxygen isotope fractionation in synthetic magnesian calcite

Mg-bearing calcite was precipitated at 25°C in closed system free-drift experiments from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The chemical and isotope composition of the solution and precipitate were investigated during time course experiments of 24-h duration. Monohydrocalcite and calcite precipitated early in the experiments (<8 h), while Mg-calcite was the predominant precipitate (>95%) thereafter. Solid collected at the end of the experiments displayed compositional zoning from pure calcite in crystal cores to up to 23 mol% MgCO₃ in the rims. Smaller excursions in Mg were superimposed on this chemical record, which is characteristic of oscillatory zoning observed in synthetic and natural solid-solution carbonates of differing solubility. Magnesium also altered the predominant morphology of crystals over time from the {104} to {100} and {110} growth forms.The oxygen isotope fractionation factor for the magnesian-calcite-water system (as 10³lnα_Mg-cl-H2O) displayed a strong dependence on the mol% MgCO₃ in the solid phase, but quantification of the relationship was difficult due to the heterogeneous nature of the precipitate. Considering only the Mg-content and δ¹⁸O values for the bulk solid, 10³lnα_Mg-cl-H2O increased at a rate of 0.17 ± 0.02 per mol% MgCO₃; this value is a factor of three higher than the single previous estimate (Tarutani T., Clayton R.N., and Mayeda T. K. (1969) The effect of polymorphims and magnesium substitution on oxygen isotope fractionation between calcium carbonate and water. Geochim. Cosmochim. Acta 33, 987-996). Nevertheless, extrapolation of our relationship to the pure calcite end member yielded a value of 27.9 ± 0.02, which is similar in magnitude to published values for the calcite-water system. Although no kinetic effect was observed on 10³lnα_Mg-cl-H2O for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol · m⁻² · h⁻¹, it was impossible to disentangle the potential effect(s) of precipitation rate and Mg-content on 10³lnα_Mg-cl-H2O due to the heterogeneous nature of the solid.The results of this study suggest that paleotemperatures inferred from the δ¹⁸O values of high magnesian calcite (>10 mol% MgCO₃) may be significantly underestimated. Also, the results underscore the need for additional experiments to accurately characterize the effect of Mg coprecipitation on the isotope systematics of calcite from a chemically homogeneous precipitate or a heterogeneous material that is analyzed at the scale of chemical and isotopic zonation.     

Grain growth kinetics of dolomite,magnesite and calcite: a comparative study

The rates of grain growth of stoichiometric dolomite [CaMg(CO₃)₂] and magnesite (MgCO₃) have been measured at temperatures T of 700–800°C at a confining pressure P _c of 300 MPa, and compared with growth rates of calcite (CaCO₃). Dry, fine-grained aggregates of the three carbonates were synthesized from high purity powders by hot isostatic pressing (HIP); initial mean grain sizes of HIP-synthesized carbonates were 1.4, 1.1, and 17 μm, respectively, for CaMg(CO₃)₂, MgCO₃, and CaCO₃, with porosities of 2, 28, and 0.04% by volume. Grain sizes of all carbonates coarsened during subsequent isostatic annealing, with mean values reaching 3.9, 5.1, and 27 μm for CaMg(CO₃)₂, MgCO₃, and CaCO₃, respectively, in 1 week. Grain growth of dolomite is much slower than the growth rates of magnesite or calcite; assuming normal grain growth and n = 3 for all three carbonates, the rate constant K for dolomite (≃5 × 10⁻⁵ μm³/s) at T = 800°C is less than that for magnesite by a factor of ~30 and less than that for calcite by three orders of magnitude. Variations in carbonate grain growth may be affected by differences in cation composition and densities of pores at grain boundaries that decrease grain boundary mobility. However, rates of coarsening correlate best with the extent of solid solution; K is the largest for calcite with extensive Mg substitution for Ca, while K is the smallest for dolomite with negligible solid solution. Secondary phases may nucleate at advancing dolomite grain boundaries, with implications for deformation processes, rheology, and reaction kinetics of carbonates.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(III) and hematite

Application of the Fe isotope system to studies of natural rocks and fluids requires precise knowledge of equilibrium Fe isotope fractionation factors among various aqueous Fe species and minerals. These are difficult to obtain at the low temperatures at which Fe isotope fractionation is expected to be largest and requires careful distinction between kinetic and equilibrium isotope effects. A detailed investigation of Fe isotope fractionation between [Fe^III(H₂O)₆]³⁺ and hematite at 98°C allows the equilibrium ⁵⁶Fe/⁵⁴Fe fractionation to be inferred, which we estimate at 10³lnα_{Fe(III)-hematite} = −0.10 ± 0.20‰. We also infer that the slope of Fe(III)-hematite fractionation is modest relative to 10⁶/T², which would imply that this fractionation remains close to zero at lower temperatures. These results indicate that Fe isotope compositions of hematite may closely approximate those of the fluids from which they precipitated if equilibrium isotopic fractionation is assumed, allowing inference of δ⁵⁶Fe values of ancient fluids from the rock record. The equilibrium Fe(III)-hematite fractionation factor determined in this study is significantly smaller than that obtained from the reduced partition function ratios calculated for [Fe^III(H₂O)₆]³⁺ and hematite based on vibrational frequencies and Mössbauer shifts by [Polyakov 1997] and [Polyakov and Mineev 2000], and Schauble et al. (2001), highlighting the importance of experimental calibration of Fe isotope fractionation factors. In contrast to the long-term (up to 203 d) experiments, short-term experiments indicate that kinetic isotope effects dominate during rapid precipitation of ferric oxides. Precipitation of hematite over ∼12 h produces a kinetic isotope fractionation where 10³lnα_{Fe(III)-hematite} = +1.32 ± 0.12‰. Precipitation under nonequilibrium conditions, however, can be recognized through stepwise dissolution in concentrated acids. As expected, our results demonstrate that dissolution by itself does not measurably fractionate Fe isotopes.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

Raman spectroscopic studies of carbonates part I: High-pressure and high-temperature behaviour of calcite,magnesite, dolomite and aragonite

The room-temperature Raman spectra of aragonite, magnesite and dolomite have been recorded up to 30 GPa and 25 GPa, respectively and no phase changes were observed during compression, unlike calcite. The effect of temperature on the room-pressure Raman spectra of calcite, aragonite, magnesite and dolomite is reported up to 800–1100 K. The measured relative pressure and temperature-shifts of the Raman lines are greater for the lattice modes than for the internal modes of the CO₃ groups. These shifts are used to calculate the mode anharmonic parameters of the observed Raman modes; they are negative and their absolute values are smaller (close to 0) for the internal CO₃ modes than for the lattice modes (4–17 10^?5 K^?1). The temperature shifts of the lattice modes in calcite are considerably larger than those for dolomite and magnesite, and a marked non-linear increase in linewidth is observed above 400° C for calcite. This is consistent with an increasing relaxational component to the libration of the CO₃ groups about their threefold axes, premonitory to the rotational order-disorder transition at higher temperature. This behaviour is not observed for the other calcite structured minerals in this study. We examine systematic variations in the lattice mode frequencies and linewidths with composition, to begin to understand these differences in their anharmonic behaviour. Finally, we have used a simple Debye-Waller model to calculate atomic displacements in calcite, magnesite and dolomite with increasing temperature from the vibrational frequency data, to provide a direct comparison with atomic positional data from high-temperature structure refinements.     

The partitioning of Fe and Mg between olivine and carbonate and the stability of carbonate under mantle conditions

We have investigated the effect of Fe on the stabilities of carbonate (carb) in lherzolite assemblages by determining the partitioning of Fe and Mg between silicate (olivine; ol) and carbonates (magnesite, dolomite, magnesian calcite) at high pressures and temperatures. Fe enters olivine preferentially relative to magnesite and ordered dolomite, but Fe and Mg partition almost equally between disordered calcic carbonate and olivine. Measurement of K _d (X _Fe ^carb X _Mg ^ol /X _Fe ^ol X _Mg ^carb ) as a function of Fe/ Mg ratio indicates that Fe–Mg carbonates deviate only slightly from ideality. Using the regular solution parameter for olivine W _FeMg ^ol of 3.7±0.8 kJ/mol (Wiser and Wood 1991) we obtain for (FeMg)CO₃ a W _FeMg ^carb of 3.05±1.50 kJ/mol. The effect of Ca–Mg–Fe disordering is to raise K _d substantially enabling us to calculate W _CaMg ^carb -W _CaFe ^carb of 5.3±2.2 kJ/mol. The activity-composition relationships and partitioning data have been used to calculate the effect of Fe/Mg ratio on mantle decarbonation and exchange reactions. We find that carbonate (dolomite and magnesian calcite) is stable to slightly lower pressures (by 1 kbar) in mantle lherzolitic assemblages than in the CaO–MgO–SiO₂(CMS)–CO₂ system. The high pressure breakdown of dolomite + orthopyroxene to magnesite + clinopyroxene is displaced to higher pressures (by 2 kbar) in natural compositions relative to CMS. CO₂. We also find a stability field of magnesian calcite in lherzolite at 15–25 kbar and 750–1000°C.     

'Forbidden zone' subduction of sediments to 150 km depth— the reaction of dolomite to magnesite + aragonite in the UHPM metapelites from western Tianshan, China

The solid‐state reaction magnesite (MgCO₃) + calcite (aragonite) (CaCO₃) = dolomite (CaMg(CO₃)₂) has been identified in metapelites from western Tianshan, China. Petrological studies show that two metamorphic stages are recorded in the metapelites: (1) the peak mineral assemblage of magnesite and calcite pseudomorphs after aragonite which is only preserved as inclusions within dolomite; and (2) the retrograde glaucophane‐chloritoid facies mineral assemblage of glaucophane, chloritoid, dolomite, garnet, paragonite, chlorite and quartz. The peak metamorphic temperatures and pressures are calculated to be 560–600 °C, 4.95–5.07 GPa based on the calcite–dolomite geothermometer and the equilibrium calculation of the reaction dolomite = magnesite + aragonite, respectively. These give direct evidence in UHP metamorphic rocks from Tianshan, China, that carbonate sediments were subducted to greater than 150 km depth. This UHP metamorphism represents a geotherm lower than any previously estimated for subduction metamorphism (< 3.7 °C km^?1) and is within what was previously considered a ‘forbidden’ condition within Earth. In terms of the carbon cycle, this demonstrates that carbonate sediments can be subducted to at least 150 km depth without releasing significant CO₂ to the overlying mantle wedge.     

Thermochemistry of carbonate-pyroxene equilibria

The enthalpies of drop solution of calcite, magnesite, dolomite, wollastonite and diopside have been measured in a lead borate solvent at 977 K in a Calvettype microcalorimeter. The carbonate calorimetry was done under flowing gas atmosphere. Both natural and synthetic samples were used. From these calorimetric data, the enthalpies of several reactions of carbonate with quartz were calculated. The enthalpies of these reactions (kJ/mol) at 298 K are: calcite+quartzwollastonite+CO₂, 92.3±1.0; magnesite+quartzenstatite+CO₂, 82.9±2.8; dolomite+quartzdiopside+CO₂, 163.0±1.9. These values generally are in agreement with those calculated from Robie et al., Helgeson et al., Berman and Holland and Powell. The enthalpy of dolomite-quartz reaction overlaps marginally with those from Berman and Holland and Powell. The enthalpy of formation of dolomite from magnesite and calcite (-11.1±2.5 kJ/mol) was also derived from the measured enthalpies, and this value is consistent with that from acid solution calorimetric measurements as shown by Navrotsky and Capobianco, but different from values in the earlier literature. These results support the premise that drop-solution of carbonates into molten lead borate results in a well-defined final state consisting of dissolved oxide and evolved CO₂. This was also confirmed by weight change experiments. Thus, oxide melt calorimetry is applicable to carbonates.     

Origin of ultramafic-hosted magnesite on Margarita Island,Venezuela

Ultramafic-hosted deposits of magnesite (MgCO₃) have been studied on Margarita Island, Venezuela, to elucidate the source of carbon and conditions of formation for this type of ore. Petrographic, mineralogic, and δ¹⁸O data indicate that magnesite precipitated on Margarita in near-surface environments at low P and T. δ¹³C ranges from −9 to −16‰ PDB within the magnesite and −8 to −10‰ PDB within some calcite and dolomite elsewhere on the island. The isotopically light dolomite fills karst and the calcite occurs as stock-work veins which resemble the magnesite deposits. These carbon isotopic ratios are consistent with a deep-seated source rather than an overlying source from a zone of surficial weathering. However, there is not much enrichment of precious metals and no enrichment of heavy rare-earth elements, as would be expected if the carbon had migrated upward as aqueous carbonate ions. The carbon probably has risen as a gaseous mixture of CO₂ and CH₄ which partially dissolved in near-surface water before leaching cations and precipitating as magnesite and other carbonates. The process probably is ongoing, given regional exhalation of carbonaceous gases.     

Calcium isotope fractionation in calcite and aragonite

Calcium isotope fractionation was measured on skeletal aragonite and calcite from different marine biota and on inorganic calcite. Precipitation temperatures ranged from 0 to 28°C. Calcium isotope fractionation shows a temperature dependence in accordance with previous observations: 1000 · ln(α_cc) = −1.4 + 0.021 · T (°C) for calcite and 1000 · ln(α_ar) = −1.9 + 0.017 · T (°C) for aragonite. Within uncertainty the temperature slopes are identical for the two polymorphs. However, at all temperatures calcium isotopes are more fractionated in aragonite than in calcite. The offset in δ^44/40Ca is about 0.6‰. The underlying mechanism for this offset may be related to the different coordination numbers and bond strengths of the calcium ions in calcite and aragonite crystals, or to different Ca reaction behavior at the solid-liquid interface. Recently, the observed temperature dependence of the Ca isotope fractionation was explained quantitatively by the temperature control on precipitation rates of calcium carbonates in an experimental setting (Lemarchand et al., 2004). We show that this mechanism can in principle also be applied to CaCO₃ precipitation in natural environments in normal marine settings. Following this model, Ca isotope fractionation in marine Ca carbonates is primarily controlled by precipitation rates. On the other hand the larger Ca isotope fractionation of aragonite compared to calcite can not be explained by different precipitation rates. The rate control model of Ca isotope fractionation predicts a strong dependence of the Ca isotopic composition of carbonates on ambient CO₃²⁻ concentration. While this model is in general accordance with our observations in marine carbonates, cultured specimens of the planktic foraminifer Orbulina universa show no dependence of Ca-isotope fractionation on the ambient CO₃²⁻ concentration. The latter observation implies that the carbonate chemistry in the calcifying vesicles of the foraminifer is independent from the ambient carbonate ion concentration of the surrounding water.     

First-principles investigation of the concentration effect on equilibrium fractionation of Ca isotopes in forsterite

Previous theoretical studies have found that the concentration variations within a certain range have a prominent effect on inter-mineral equilibrium isotope fractionation (10³lnα). Based on the density functional theory, we investigated how the average Ca–O bond length and the reduced partition function ratios (10³lnβ) and 10³lnα of ⁴⁴Ca/⁴⁰Ca in forsterite (Fo) are affected by its Ca concentration. Our results show that Ca–O bond length in forsterite ranges from 2.327 to 2.267 Å with the Ca/(Ca + Mg) varying between a narrow range limited by an upper limit of 1/8 and a lower limit of 1/64. However, outside this narrow range, i.e., Ca/(Ca + Mg) is lower than 1/64 or higher than 1/8, Ca–O bond length becomes insensitive to Ca concentration and maintains to be a constant. Because the 10³lnβ is negatively correlated with Ca–O bond length, the 10³lnβ significantly increases with decreasing Ca/(Ca + Mg) when 1/64 < Ca/(Ca + Mg) < 2/16. As a consequence, the 10³lnα between forsterite and other minerals also strongly depend on the Ca content in forsterite. Combining previous studies with our results, the heavier Ca isotopes enrichment sequence in minerals is: forsterite > orthopyroxene > clinopyroxene > calcite ≈ diopside > dolomite > aragonite. Olivine and pyroxenes are enriched in heavier Ca isotope compared to carbonates. The 10³lnα between forsterite with a Ca/(Ca + Mg) of 1/64 and clinopyroxene (Ca/Mg = 1/1, i.e., diopside) is up to ~ 0.64‰ at 1200 K. The large 10³lnα_Fo-diopside relative to the current analytical precision for Ca isotope measurements suggests that the dependence of 10³lnα_Fo-diopside on temperature can be used as a thermometer, similar to the one based on the 10³lnα of ⁴⁴Ca/⁴⁰Ca between orthopyroxene and diopside. These two Ca isotope thermometers both have a precision approximate to that of elemental thermometers and provide independent constraints on temperature.

     

Ca isotopes in carbonate sediment and pore fluid from ODP Site 807A: The Ca(aq)-calcite equilibrium fractionation factor and calcite recrystallization rates in Pleistocene sediments

The calcium isotopic compositions (δ⁴⁴Ca) of 30 high-purity nannofossil ooze and chalk and 7 pore fluid samples from ODP Site 807A (Ontong Java Plateau) are used in conjunction with numerical models to determine the equilibrium calcium isotope fractionation factor (α_s−f) between calcite and dissolved Ca²⁺ and the rates of post-depositional recrystallization in deep sea carbonate ooze. The value of α_s−f at equilibrium in the marine sedimentary section is 1.0000 ± 0.0001, which is significantly different from the value (0.9987 ± 0.0002) found in laboratory experiments of calcite precipitation and in the formation of biogenic calcite in the surface ocean. We hypothesize that this fractionation factor is relevant to calcite precipitation in any system at equilibrium and that this equilibrium fractionation factor has implications for the mechanisms responsible for Ca isotope fractionation during calcite precipitation. We describe a steady state model that offers a unified framework for explaining Ca isotope fractionation across the observed precipitation rate range of ∼14 orders of magnitude. The model attributes Ca isotope fractionation to the relative balance between the attachment and detachment fluxes at the calcite crystal surface. This model represents our hypothesis for the mechanism responsible for isotope fractionation during calcite precipitation. The Ca isotope data provide evidence that the bulk rate of calcite recrystallization in freshly-deposited carbonate ooze is 30-40%/Myr, and decreases with age to about 2%/Myr in 2-3 million year old sediment. The recrystallization rates determined from Ca isotopes for Pleistocene sediments are higher than those previously inferred from pore fluid Sr concentration and are consistent with rates derived for Late Pleistocene siliciclastic sediments using uranium isotopes. Combining our results for the equilibrium fractionation factor and recrystallization rates, we evaluate the effect of diagenesis on the Ca isotopic composition of marine carbonates at Site 807A. Since calcite precipitation rates in the sedimentary column are many orders of magnitude slower than laboratory experiments and the pore fluids are only slightly oversaturated with respect to calcite, the isotopic composition of diagenetic calcite is likely to reflect equilibrium precipitation. Accordingly, diagenesis produces a maximum shift in δ⁴⁴Ca of +0.15‰ for Site 807A sediments but will have a larger impact where sedimentation rates are low, seawater circulates through the sediment pile, or there are prolonged depositional hiatuses.     

Partitioning of strontium between calcite,dolomite and liquids: An experimental study under higher temperature diagenetic conditions,and a model for the prediction of mineral pairs for geothermometry

The partitioning of Sr between calcite, dolomite and liquids is essentially independent of temperature between 150° and 350° C. The partition coefficients corrected for number of cation sites are b _calc=0.096 and b _dol= 0.048 for 1 mol cations/6 mol H₂O liquid. Upon dilution the partition coefficients increase, but their ratio stays constant at about 2∶1. This ratio is due to the fact that calcite has twice as many Ca-sites for Sr-substitution as dolomite. The 2∶1 relationship is also observed in natural calcite and dolomite which have undergone diagenesis. The temperature independence of partitioning is caused by the relatively small thermal expansion of calcite and dolomite. Thermal expansion between 25° and 400° C was found to follow the equations V _calc=7.0·10⁻⁴ T(°C)+36.95 and V _dol=6.9·10⁻⁴ T(°C)+32.24, V: cm³/mol. Therefore calcite and dolomite cannot serve as a temperature indicator. To have an ideal geothermometer a mineral pair with high and low thermal expansion is required. Literature date demonstrate that wurtzite, sphalerite, and galena are such minerals.     

Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments

The chemical and isotopic composition of speleothem calcite and particularly that of stalagmites and flowstones is increasingly exploited as an archive of past environmental change in continental settings. Despite intensive research, including modelling and novel approaches, speleothem data remain difficult to interpret. A possible way foreword is to apply a multi-proxy approach including non-conventional isotope systems. For the first time, we here present a complete analytical dataset of magnesium isotopes (δ²⁶Mg) from a monitored cave in NW Germany (Bunker Cave). The data set includes δ²⁶Mg values of loess-derived soil above the cave (−1.0 ± 0.5‰), soil water (−1.2 ± 0.5‰), the carbonate hostrock (−3.8 ± 0.5‰), dripwater in the cave (−1.8 ± 0.2‰), speleothem low-Mg calcite (stalactites, stalagmites; −4.3 ± 0.6‰), cave loam (−0.6 ± 0.1‰) and runoff water (−1.8 ± 0.1‰) in the cave, respectively. Magnesium-isotope fractionation processes during weathering and interaction between soil cover, hostrock and solute-bearing soil water are non-trivial and depend on a number of variables including solution residence times, dissolution rates, adsorption effects and potential neo-formation of solids in the regolith and the carbonate aquifer. Apparent Mg-isotope fractionation between dripwater and speleothem low-Mg calcite is about 1000lnα_Mg-cc-Mg(aq) = −2.4‰. A similar Mg-isotope fractionation (1000lnα_Mg-cc-Mg(aq) ≈ −2.1‰) is obtained by abiogenic precipitation experiments carried out at aqueous Mg/Ca ratios and temperatures close to cave conditions. Accordingly, ²⁶Mg discrimination during low-Mg calcite formation in caves is highly related to inorganic fractionation effects, which may comprise dehydration of Mg²⁺ prior to incorporation into calcite, surface entrapment of light isotopes and reaction kinetics. Relevance of kinetics is supported by a significant negative correlation of Mg-isotope fractionation with the precipitation rate for inorganic precipitation experiments.     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: II. Ca isotopes

Ca isotope fractionation during inorganic calcite formation was experimentally studied by spontaneous precipitation at various precipitation rates (1.8 < log R < 4.4 μmol/m²/h) and temperatures (5, 25, and 40 °C) with traces of Sr using the CO₂ diffusion technique.Results show that in analogy to Sr/Ca [see Tang J., Köhler S. J. and Dietzel M. (2008) Sr²⁺/Ca²⁺ and ⁴⁴Ca/⁴⁰Ca fractionation during inorganic calcite formation: I. Sr incorporation. Geochim. Cosmochim. Acta] the ⁴⁴Ca/⁴⁰Ca fractionation during calcite formation can be followed by the Surface Entrapment Model (SEMO). According to the SEMO calculations at isotopic equilibrium no fractionation occurs (i.e., the fractionation coefficient α_calcite-aq = (⁴⁴Ca/⁴⁰Ca)_s/(⁴⁴Ca/⁴⁰Ca)_aq = 1 and Δ^44/40Ca_calcite-aq = 0‰), whereas at disequilibrium ⁴⁴Ca is fractionated in a primary surface layer (i.e., the surface entrapment factor of ⁴⁴Ca, F_44Ca < 1). As a crystal grows at disequilibrium, the surface-depleted ⁴⁴Ca is entrapped into the newly formed crystal lattice. ⁴⁴Ca depletion in calcite can be counteracted by ion diffusion within the surface region. Our experimental results show elevated ⁴⁴Ca fractionation in calcite grown at high precipitation rates due to limited time for Ca isotope re-equilibration by ion diffusion. Elevated temperature results in an increase of ⁴⁴Ca ion diffusion and less ⁴⁴Ca fractionation in the surface region. Thus, it is predicted from the SEMO that an increase in temperature results in less ⁴⁴Ca fractionation and the impact of precipitation rate on ⁴⁴Ca fractionation is reduced.A highly significant positive linear relationship between absolute ⁴⁴Ca/⁴⁰Ca fractionation and the apparent Sr distribution coefficient during calcite formation according to the equation

Δ^44/40Ca_calcite-aq=(−1.90±0.26)·logD_Sr−2.83±0.28     

Dissolution Kinetics of Calcite in 0.1 M NaCl Solution at Room Temperature: An Atomic Force Microscopic (AFM) Study

Atomic force microscopy (AFM) was used to study the rates of migration of the (10¯1 4) plane of a single-crystal of calcite dissolving in 0.1 M NaCl aqueous solutions at room temperature. The solution pH and P_{CO 2} controlled in the ranges 4.4 < pH < 12.2 and 0 < P_{CO 2} < 10^-3.5 atm (ambient), respectively. Measured step velocities were compared with the mineral dissolution rates determined from the calcium fluxes. The step velocity is defined as the average of the velocities of the obtuse and acute steps. Rates of step motion increased gradually from 1.4(±0.2) at pH 5.3 to 2.4(±0.3) nm s^-1 at pH 8.2, whereas the rates inverted and decreased to the minimum value of 0.69(±0.18) nm s^-1 at pH 10.8. For pH > 10.8, only the velocity of the obtuse steps increased as pH increased, whereas that of acute steps gradually decreased.The dissolution rate of the mineral can be calculated from the measured step velocities and average slope, which is proportional to the concentration of exposed monomolecular steps on the surface. The average slope of the dissolving mineral, measured at pH 5.6 and 9.7, was 0.026 (±0.015). Using this slope, we calculate bulk dissolution rates for 5.3 < pH < 12.2 of 4.9(±3.0) × 10^-11 to 1.8(±1.0) × 10^-10 mol cm^-2 s^-1. The obtained dissolution rate can be expressed by the following empirical equation:R_dss = 10^{-4.66(±0.13)}[H⁺] + 10^{-3.87(±0.06)}[HCO₃ ^-] + 10^{-7.99(plusmn; 0.08)}[OH^-]We propose that calcite dissolution in these solutions is controlled by elementary reactions that are similar to those that control the dissolution of other amphoteric solids, such as oxides. The mechanisms include the proton-enhanced hydration and detachment of calcium-carbonate ion pairs. The detachments are enhanced by the presence of adsorbed nucleophiles, such as hydroxyl and bicarbonate ions, and by protons adsorbed to key oxygens. A molecular model is proposed that illustrates these processes.     

Experimental studies of oxygen isotope fractionation between rhodochrosite (MnCO3) and water at low temperatures

Rhodochrosite crystals were precipitated from Na-Mn-Cl-HCO₃ parent solutions following passive, forced and combined passive-to-forced CO₂ degassing methods. Forced and combined passive-to-forced CO₂ degassing produced rhodochrosite crystals with a small non-equilibrium oxygen isotope effect whereas passive CO₂ degassing protocols yielded rhodochrosite in apparent isotopic equilibrium with water. On the basis of the apparent equilibrium isotopic data, a new temperature-dependent relation is proposed for the oxygen isotope fractionation between rhodochrosite and water between 10 and 40 °C:

1000lnα_{rhodochrosite-water}=17.84±0.18(10³/T)-30.24±0.62     

An experimental and theoretical determination of oxygen isotope fractionation in the system magnetite-H2O from 300 to 800°C

Oxygen isotope fractionations have been determined between magnetite and water from 300 to 800°C and pressures between 10 and 215MPa. We selected three reaction pathways to investigate fractionation: (a) reaction of fine-grained magnetite with dilute aqueous NaCl solutions; (b) reduction of fine-grained hematite through reaction with dilute acetic acid; and (c) oxidation of fine iron power in either pure water or dilute NaCl solutions. Effective use of acetic acid was limited to temperatures up to about 400°C, whereas oxide-solution isotope exchange experiments were conducted at all temperatures. Equilibrium ¹⁸O/¹⁶O fractionation factors were calculated from the oxide-water experiments by means of the partial isotope exchange method, where generally four isotopically different waters were used at any given temperature. Each run product was characterized by X-ray diffraction (XRD), scanning electron microscopy (SEM), and on a limited basis, high-resolution transmission electron microscopy (HRTEM) and Mössbauer spectroscopy. Results from the microscopic examinations indicate the formation of well-crystallized octahedra and dodecahedra of magnetite where the extent of crystallization, grain size, and grain habit depend on the initial starting material, P, T, solution composition, and duration of the run.The greatest amount of oxygen isotope exchange (∼90% or greater) was observed in experiments where magnetite either recrystallized in the presence of 0.5 m NaCl from 500 to 800°C or formed from hematite reacted with 0.5 m acetic acid at 300, 350 and 400°C. Fractionation factors (10³ ln α_mt-H2O) determined from these partial exchange experiments exhibit a steep decrease (to more negative values) with decreasing temperature down to about 500°C, followed by shallower slope. A least-squares regression model of these partial exchange data, which accounts for analytical errors and errors generated by mass balance calculations, gives the following expression for fractionation that exhibits no minimum: 1000lnα_lmt-lw=−8.984(±0.3803)x+3.302(±0.377)x²—0.426(±0.092)x³ with an R² = 0.99 for 300 ≤ T≤ 800°C (x = 10⁶/T²). The Fe oxidation results also exhibit this type of temperature dependence but shifted to slightly more negative 10³ ln α values; there is the suggestion that a kinetic isotope effect may contribute to these fractionations. A theoretical assessment of oxygen isotope fractionation using β-factors derived from heat capacity and Mössbauer temperature (second-order Doppler) shift measurements combined with known β-factors for pure water yield fractionations that are somewhat more negative compared to those determined experimentally. This deviation may be due to the combined solute effects of dissolved magnetite plus NaCl (aq), as well as an underestimation of β_mt at low temperatures. The new magnetite-water experimental fractionations agree reasonably well with results reported from other experimental studies for temperatures ≥ 500°C, but differ significantly with estimates based on quasi-theoretical and empirical approaches. Calcite-magnetite and quartz-magnetite fractionation factors estimated from the combination of magnetite β’s calculated in this study with those for calcite and quartz reported by Clayton and Kieffer (1991) agree very closely with experimentally determined mineral-pair fractionations.     

Origin of the Rubian carbonate-hosted magnesite deposit, Galicia, NW Spain: mineralogical, REE, fluid inclusion and isotope evidence

The Rubian magnesite deposit (West Asturian—Leonese Zone, Iberian Variscan belt) is hosted by a 100-m-thick folded and metamorphosed Lower Cambrian carbonate/siliciclastic metasedimentary sequence—the Cándana Limestone Formation. It comprises upper (20-m thickness) and lower (17-m thickness) lens-shaped ore bodies separated by 55 m of slates and micaceous schists. The main (lower) magnesite ore body comprises a package of magnesite beds with dolomite-rich intercalations, sandwiched between slates and micaceous schists. In the upper ore body, the magnesite beds are thinner (centimetre scale mainly) and occur between slate beds. Mafic dolerite dykes intrude the mineralisation. The mineralisation passes eastwards into sequence of bedded dolostone (Buxan) and laminated to banded calcitic marble (Mao). These show significant Variscan extensional shearing or fold-related deformation, whereas neither Rubian dolomite nor magnesite show evidence of tectonic disturbance. This suggests that the dolomitisation and magnesite formation postdate the main Variscan deformation. In addition, the morphology of magnesite crystals and primary fluid inclusions indicate that magnesite is a neoformed hydrothermal mineral. Magnesite contains irregularly distributed dolomite inclusions (<50 μm) and these are interpreted as relics of a metasomatically replaced dolostone precursor. The total rare earth element (REE) contents of magnesite are very similar to those of Buxan dolostone but are depleted in light rare earth elements (LREE); heavy rare earth element concentrations are comparable. However, magnesite REE chondrite normalised profiles lack any characteristic anomaly indicative of marine environment. Compared with Mao calcite, magnesite is distinct in terms of both REE concentrations and patterns. Fluid inclusion studies show that the mineralising fluids were MgCl₂–NaCl–CaCl₂–H₂O aqueous brines exhibiting highly variable salinities (3.3 to 29.5 wt.% salts). This may be the result of a combination of fluid mixing, migration of pulses of variable-salinity brines and/or local dissolution and replacement processes of the host dolostone. Fluid inclusion data and comparison with other N Iberian dolostone-hosted metasomatic deposits suggest that Rubian magnesite probably formed at temperatures between 160 and 200°C. This corresponds, at hydrostatic pressure (500 bar), to a depth of formation of ~~5 km. Mineralisation-related Rubian dolomite yields δ ¹⁸O values (δ ¹⁸O: 12.0–15.4‰, mean: 14.4±1.1‰) depleted by around 5‰ compared with barren Buxan dolomite (δ ¹⁸O: 17.1–20.2‰, mean: 19.4±1.0‰). This was interpreted to reflect an influx of ¹⁸O-depleted waters accompanied by a temperature increase in a fluid-dominated system. Overlapping calculated δ ¹⁸O_fluid values (~+5‰ at 200°C) for fluids in equilibrium with Rubian dolomite and magnesite show that they were formed by the same hydrothermal system at different temperatures. In terms of δ ¹³C values, Rubian dolomite (δ ¹³C: −1.4 to 1.9‰, mean: 0.4±1.3‰) and magnesite (δ ¹³C: −2.3 to 2.4‰, mean: 0.60±1.0‰) generally exhibit more negative δ ¹³C values compared with Buxan dolomite (δ ¹³C: −0.2 to 1.9‰, mean: 0.8±0.6‰) and Mao calcite (δ ¹³C: −0.3 to 1.5‰, mean: 0.6±0.6‰), indicating progressive modification to lower δ ¹³C values through interaction with hydrothermal fluids. ⁸⁷Sr/⁸⁶Sr ratios, calculated at 290 Ma, vary from 0.70849 to 0.70976 for the Mao calcite and from 0.70538 to 0.70880 for the Buxan dolostone. The ⁸⁷Sr/⁸⁶Sr ratios in Rubian magnesite are more radiogenic and range from 0.71123 to 0.71494. The combined δ ¹⁸O–δ ¹³C and ⁸⁷Sr/⁸⁶Sr data indicate that the magnesite-related fluids were modified basinal brines that have reacted and equilibrated with intercalated siliciclastic rocks. Magnesite formation is genetically linked to regional hydrothermal dolomitisation associated with lithospheric delamination, late-Variscan high heat flow and extensional tectonics in the NW Iberian Belt. A comparison with genetic models for the Puebla de Lillo talc deposits suggests that the formation of hydrothermal replacive magnesite at Rubian resulted from a metasomatic column with magnesite forming at higher fluid/rock ratios than dolomite. In this study, magnesite generation took place via the local reaction of hydrothermal dolostone with the same hydrothermal fluids in very high permeability zones at high fluid/rock ratios (e.g. faults). It was also possibly aided by additional heat from intrusive dykes or sub-cropping igneous bodies. This would locally raise isotherms enabling a transition from the dolomite stability field to that of magnesite.Editorial handling: F. Tornos     

Antigorite-ophicarbonates: Contact metamorphism in Valmalenco,Italy

Outside the Bergell tonalite contact aureole, ophicarbonate rocks consist of blocks of antigorite schist embedded in veins of calcite ± tremolite. An antigorite schistosity predates some of these calcite veins. Mono- and bimineralic assemblages occur in reaction zones associated with the veins. Within the aureole, the ophicarbonate veining becomes less distinct and polymineralic assemblages become more frequent. A regular sequence of isobaric univariant assemblages is found, separated by isograds corresponding to isobaric invariant assemblages. In order of increasing grade the invariant assemblages are: antigorite+diopside+olivine+tremolite+calcite antigorite+dolomite+olivine+tremolite+calcite antigorite+olivine+talc+magnesite antigorite+dolomite+olivine+tremolite+talc These assemblages match a previously derived topology in P-T-X_CO2 space for the system CaO-MgO-SiO₂-H₂O-CO₂; the field sequence can be used to adjust the relative locations of calculated invariant points with respect to temperature. Isobaric univariant and invariant assemblages are plotted along a profile map to permit direct comparison with the phase diagram.It is inferred that, during the formation of the ophicarbonate veins, calcite precipitated from fluid introduced into the serpentinite. During contact metamorphism, however, the compositions of pore fluids evolved by reaction in the ophicarbonate rocks were largely buffered by the solid phases. This control occurred on a small scale, because there are local variations in the buffering solid assemblages within a centimeter range.