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.     

Chemistry, petrology and bulk oxygen isotope compositions of chondrules from the Mokoia CV3 carbonaceous chondrite

We report bulk chemical compositions and physical properties for a suite of 94 objects, mostly chondrules, separated from the Mokoia CV3_ox carbonaceous chondrite. We also describe mineralogical and petrologic information for a selected subset of the same suite of chondrules. The data are used to examine the range of chondrule bulk compositions, and to investigate the relationships between chondrule mineralogy, texture and bulk compositions, as well as oxygen isotopic properties that we reported previously. Most of the chondrules show minimal metamorphism, corresponding to petrologic subtype <3.2. In general, elemental fractionations observed in chondrule bulk compositions are reflected in the compositions of constituent minerals. For chondrules, mean bulk compositions and compositional ranges are very similar for large (>2 mg) and small (<2 mg) size fractions. Two of the objects studied are described as matrix-rich clasts. These have similar bulk compositions to the chondrule mean, and are potential chondrule precursors. One of these clasts has a similar bulk oxygen isotopic composition to Mokoia chondrules, but the other has an anomalously high value of Δ¹⁷O (+3.60‰).Chondrules are diverse in bulk chemical composition, with factor of 10 variations in most major element abundances that cannot be attributed to secondary processes. The chondrules examined show evidence for extensive secondary oxidation, and possible sulfidization, as expected for an oxidized CV chondrite, but minimal aqueous alteration. Some of the bulk chondrule compositional variation might be the result of chemical (e.g. volatilization or condensation) or physical (e.g. metal loss) processes during chondrule formation. However, we suggest that it is mainly the result of significant variations in the assembly of particles that constituted chondrule precursors. Precursor material likely included a refractory component, possibly inherited from disaggregated CAIs, an FeO-poor ferromagnesian component such as olivine or pyroxene, an oxidized ferromagnesian component, and a metal component. Bulk oxygen isotope ratios of chondrules can be explained if refractory and ferromagnesian precursor materials initially shared similar oxygen isotopic compositions of δ¹⁷O, δ¹⁸O around −50‰, and then significant exchange occurred between the chondrule and surrounding ¹⁶O-poor gas during melting.     

Hydrogen isotope evidence for the origin and evolution of the carbonaceous chondrites

We present new hydrogen isotope data for separated matrix, hydrated chondrules, and other hydrated coarse silicate fragments from nine carbonaceous chondrites. These data were generated using a micro-analytical method involving stepped combustion of tens to hundreds of micrograms of hydrous solids. We also re-evaluate hydrogen isotope data from previous conventional stepped combustion experiments on these and other carbonaceous chondrites.Hydrogen isotope compositions of matrix and whole-rock samples of CM chondrites are correlated with oxygen isotope indices, major and minor-element abundances, and abundance and isotope ratios of other highly volatile elements. These correlations include a monotonic decrease in δD with increasing extent of aqueous alteration and decreasing abundances of highly volatile elements (including C, N and Ar), between extremes of ∼0‰ (least altered, most volatile rich) and −200‰ (most altered, least volatile rich). In plots involving only abundances and/or isotope ratios of highly volatile elements, CI chondrites fall on the high-δD, volatile rich end of the trends defined by CM chondrites; i.e., CI chondrites resemble the least altered CM chondrites in these respects. These trends suggest the protoliths of the CM chondrites (i.e., before aqueous alteration) contained an assemblage of volatiles having many things in common with those in the CI chondrites. If so, then the volatile-element inventory of the CI chondrites was a more widespread component of early solar system objects than suggested by the scarcity of recognized CI meteorites. Differences in volatile-element chemistry between the CI and average CM chondrites can be attributed to aqueous alteration of the latter.Previous models of carbonaceous chondrite aqueous alteration have suggested: (1) the protoliths of the CM chondrites are volatile poor objects like the CO or CV chondrites; and (2) the CI chondrites are more altered products of the same process producing the CM chondrites. Both suggestions appear to be inconsistent with hydrogen isotope data and other aspects of the volatile-element geochemistry of these rocks. We present a model for aqueous alteration of the CM chondrites that reconciles these inconsistencies and suggests revised relationships among the major subtypes of carbonaceous chondrites. Our model requires, among other things, that the water infiltrating CM chondrites had a δD value of ∼−158‰, consistent with initial accretion of CM parent bodies at ∼4 AU.     

European vegetation during Marine Oxygen Isotope Stage-3

European vegetation during representative “warm” and “cold” intervals of stage-3 was inferred from pollen analytical data. The inferred vegetation differs in character and spatial pattern from that of both fully glacial and fully interglacial conditions and exhibits contrasts between warm and cold intervals, consistent with other evidence for stage-3 palaeoenvironmental fluctuations. European vegetation thus appears to have been an integral component of millennial environmental fluctuations during stage-3; vegetation responded to this scale of environmental change and through feedback mechanisms may have had effects upon the environment. The pollen-inferred vegetation was compared with vegetation simulated using the BIOME 3.5 vegetation model for climatic conditions simulated using a regional climate model (RegCM2) nested within a coupled global climate and vegetation model (GENESIS-BIOME). Despite some discrepancies in detail, both approaches capture the principal features of the present vegetation of Europe. The simulated vegetation for stage-3 differs markedly from that inferred from pollen analytical data, implying substantial discrepancy between the simulated climate and that actually prevailing. Sensitivity analyses indicate that the simulated climate is too warm and probably has too short a winter season. These discrepancies may reflect incorrect specification of sea surface temperature or sea-ice conditions and may be exacerbated by vegetation-climate feedback in the coupled global model.     

Carbon isotope effects in carbonate systems

Global carbon cycle models require a complete understanding of the δ¹³C variability of the Earth’s C reservoirs as well as the C isotope effects in the transfer of the element among them. An assessment of δ¹³C changes during CO₂ loss from degassing magmas requires knowledge of the melt-CO₂ carbon isotope fractionation. In order to examine the potential size of this effect for silicate melts of varying composition, ¹³C reduced partition functions were computed in the temperature range 275 to 4000 K for carbonates of varying bond strengths (Mg, Fe, Mn, Sr, Ba, Pb, Zn, Cd, Li, and Na) and the polymorphs of calcite. For a given cation and a given pressure the ¹³C content increases with the density of the carbonate structure. For a given structure the tendency to concentrate ¹³C increases with pressure. The effect of pressure (‰/10 kbar) on the size of the reduced partition function of aragonite varies with temperature; in the pressure range 1 to 10⁵ bars the change is given by:

(1)     

Manganese-chromium isotope systematics of enstatite meteorites

We present results of a study of the ⁵³Mn-⁵³Cr isotope systematics in the enstatite chondrites and achondrites (aubrites). The goal of this study was to explore the capabilities of this isotope system to obtain chronological information on these important classes of meteorites and to investigate the original distribution in the inner solar system of the short-lived radionuclide ⁵³Mn. Our earlier work (Lugmair and Shukolyukov, 1998; Shukolyukov and Lugmair, 2000a) has shown that the asteroid belt bodies are characterized by essentially the same initial ⁵³Mn abundance. However, we have found the presence of a gradient in the abundance of the radiogenic ⁵³Cr between the earth-moon system, Mars, and the asteroid Vesta. If this gradient is considered as a function of the heliocentric distance a linear radial dependence is indicated. This can be explained either by an early, volatility controlled Mn/Cr fractionation in the nebula or by an original radially heterogeneous distribution of ⁵³Mn. The enstatite chondrites are suggested to form in the inner zones of the solar nebula, much closer to the Sun than the ordinary chondrites. Therefore, their investigation may be an important test on the hypothesis on a radial heterogeneity in the initial ⁵³Mn.We have studied the bulk samples of the EH4-chondrites Indarch and Abee and the EL6-chondrite Khairpur. Although these meteorites have essentially the same Mn/Cr ratio as the ordinary chondrites, the relative abundance of the radiogenic ⁵³Cr is three times smaller than in the ordinary chondrites. Because these meteorites are primitive (undifferentiated) and no Mn/Cr fractionation had occurred within their parent bodies, this difference is a strong argument in favor of an initially heterogeneous distribution of ⁵³Mn in the early inner solar system. This finding is also consistent with formation of the enstatite chondrites in the inner zones of the solar nebula. Using the characteristic ⁵³Cr excess of the enstatite chondrites and the observed gradient, their place of origin falls at about 1.4 AU or somewhat closer to the Sun (i.e. >1.0-1.4 AU).We also present chronological results for the enstatite chondrites and achondrites. The ‘absolute’ ⁵³Mn-⁵³Cr ages of the EH4-chondrites are old: ∼4565 Ma. The EL6-chondrite Khairpur is ∼4.5 Ma younger, which is in good agreement with the ¹²⁹I-¹²⁹Xe data from the literature. The age of the aubrite Peña Blanca Spring appears to be similar to those of the enstatite chondrites while that of the aubrite Bishopville is at least ∼10 Ma younger, which is also in agreement with the ¹²⁹I-¹²⁹Xe data. The results from bulk samples of aubrites indicate that the last Mn/Cr fractionation in their parent body occurred ∼ 4563 Ma ago and imply an evolution of the Mn-Cr isotope system in an environment with an higher than chondritic Mn/Cr ratio for several millions of years.     

Indicators of parent-body processes: Hydrated chondrules and fine-grained rims in the Mokoia CV3 carbonaceous chondrite

A petrographic and electron microscopic study of the Mokoia CV3 carbonaceous chondrite shows that all of the chondrules and inclusions (>400 μm in diameter) and most of their fine-grained rims studied (referred to as chondrules/rims) contain various amounts of hydrous phyllosilicates (mostly saponite) formed by aqueous alteration of anhydrous silicates. The rims mainly consist of fine-grained olivine and saponite in varying proportions and contain crosscutting veins of Fe-rich olivine. The boundaries between the chondrules and their rims are irregular and show abundant evidence of aqueous alteration interactions between them. In contrast, the host matrix contains very minor amounts of saponite and shows no evidence of such extensive aqueous alteration. The boundaries between the chondrules/rims and the matrix are sharp and show no traces of the matrix having been involved in the alteration of the chondrules/rims. These observations indicate that the aqueous alteration in the chondrules/rims did not occur in the present setting.We suggest that the chondrules/rims are actually clasts transported from a location on the meteorite parent body different from where the Mokoia meteorite was from. The aqueous alteration of the chondrules/rims probably occurred there. The veins in the rims were originally fractures produced in an interchondrule matrix by impacts; these were later filled by Fe-rich olivine during aqueous activity. This location was then involved in impact brecciation, and individual chondrules were ejected as clasts with remnants of the matrix surrounding them. During the continuing brecciation, those chondrule/rim clasts were transported, mixed with anhydrous matrix grains, and finally lithified to the present meteorite. Therefore, the rims are fragmented remnants of a former matrix.Textures characterized by fine-grained rims surrounding chondrules in chondrites have been widely thought to have formed in the solar nebula before they accreted into their parent bodies. However, our results suggest that some textures may not be explained by such an accretionary model; instead, the multi-stage parent-body process modeled for the Mokoia rim formation may be a more plausible explanation.     

The formation of chondrules by open-system melting of nebular condensates

Experiments were conducted under canonical nebular conditions to see whether the chemical compositions of the various chondrule types can be derived from a single CI-like starting material by open-system melting and evaporation. Experimental charges, produced at 1580 °C and P_H2 of 1.31×10⁻⁵ atm over 1 to 18 hours, consisted of only two phases, porphyritic olivine crystals in glass. Sulfur, metallic-iron and alkalis were completely evaporated in the first minutes of the experiments and subsequently the main evaporating liquid oxides were FeO and SiO₂. Olivines from short runs (2-4 hours) have compositions of Fo₈₃-Fo₈₉, as in Type IIA chondrules, while longer experimental runs (12-18 hours) produce ∼Fo₉₉ olivine, similar to Type IA chondrules. The concentration of CaO in both olivine (up to 0.6 wt.%) and glass, and their Mg#, increased with increasing heating duration. Natural chondrules also show increasing CaO with decreasing S, alkalis, FeO and SiO₂. The similarities in bulk chemistry, mineralogy and textures between Type IIA and IA chondrules and the experimental charges demonstrate that these chondrules could have formed by the evaporation of CI precursors. The formation of silica-rich chondrules (IIB and IB) by evaporation requires a more pyroxene-rich precursor.Based on the FeO evaporation rates measured here, Type IIA and IA chondrules, were heated for at least ∼0.5 and ∼3.5 h, respectively, if formed at 1580 °C and P_H2 of 1.31×10⁻⁵ atm. Type II chondrules may have experienced higher cooling-rates and less evaporation than Type I.The experimental charges experienced free evaporation and exhibited heavy isotopic enrichments in silicon, as well as zero concentrations of S, Na and K, which are not observed in natural chondrules. However, experiments on potassium-rich melts at the same pressure but in closed capsules showed less evaporation of K, and less K isotopic mass fractionation, than expected as a function of decreasing cooling rate. Thus the environment in which chondrules formed is as important as the kinetic processes they experienced. If chondrule formation occurred under conditions in which evaporated gases remained in the vicinity of the residual melts, the extent of evaporation would be reduced and back reaction between the gas and the melt could contribute to the suppression of isotopic mass fractionation. Hence chondrule formation could have involved evaporative loss without Rayleigh fractionation. Volatile-rich Type II and volatile-poor Type I chondrules may have formed in domains with high and low chondrule concentrations, and high partial pressures of lithophile elements, respectively.     

Iron isotope fractionation by Fe(II)-oxidizing photoautotrophic bacteria

Photoautotrophic bacteria that oxidize ferrous iron (Fe[II]) under anaerobic conditions are thought to be ancient in origin, and the ferric (hydr)oxide mineral products of their metabolism are likely to be preserved in ancient rocks. Here, two enrichment cultures of Fe(II)-oxidizing photoautotrophs and a culture of the genus Thiodictyon were studied with respect to their ability to fractionate Fe isotopes. Fe isotope fractionations produced by both the enrichment cultures and the Thiodictyon culture were relatively constant at early stages of the reaction progress, where the ⁵⁶Fe/⁵⁴Fe ratios of poorly crystalline hydrous ferric oxide (HFO) metabolic products were enriched in the heavier isotope relative to aqueous ferrous iron (Fe[II]_aq) by ∼1.5 ± 0.2‰. This fractionation appears to be independent of the rate of photoautotrophic Fe(II)-oxidation, and is comparable to that observed for Fe isotope fractionation by dissimilatory Fe(III)-reducing bacteria. Although there remain a number of uncertainties regarding how the overall measured isotopic fractionation is produced, the most likely mechanisms include (1) an equilibrium effect produced by biological ligands, or (2) a kinetic effect produced by precipitation of HFO overlaid upon equilibrium exchange between Fe(II) and Fe(III) species. The fractionation we observe is similar in direction to that measured for abiotic oxidation of Fe(II)_aq by molecular oxygen. This suggests that the use of Fe isotopes to identify phototrophic Fe(II)-oxidation in the rock record may only be possible during time periods in Earth’s history when independent evidence exists for low ambient oxygen contents.     

Oxygen isotope stage 3 fluvial and eolian successions in Europe compared with climate model results

Fluvial and eolian successions of oxygen isotope stage 3 are compared with global (GCM) and regional climate (RCM) modeling experiments of the stage 3 and last glacial maximum climate in Europe. Differences in precipitation between stage-3 stades and interstades were minor, which is confirmed by the fluvial successions. The fluvial response to climate variation is non-uniform, and in southern Europe more pronounced than in northern Europe. The model simulations indicate a strong western winter circulation over Europe during stage 3, which is supported by the eolian deposits data. Wind speeds in the last glacial maximum simulation appear modest compared with those of stage 3, which contrasts with the abundance of eolian deposits. This suggests that during glacial climates the stabilizing effect of vegetation determines eolian sedimentation rates, rather than wind speed. Stage 3 can be divided into an older part (>45,000 cal yr B.P.) with a relatively stable landscape and moist climate and a younger part with more frequent climate change and decreasing landscape stability.     

Ca,Al-rich inclusions, amoeboid olivine aggregates, and Al-rich chondrules from the unique carbonaceous chondrite Acfer 094: I. mineralogy and petrology

Based on their mineralogy and petrography, ∼200 refractory inclusions studied in the unique carbonaceous chondrite, Acfer 094, can be divided into corundum-rich (0.5%), hibonite-rich (1.1%), grossite-rich (8.5%), compact and fluffy Type A (spinel-melilite-rich, 50.3%), pyroxene-anorthite-rich (7.4%), and Type C (pyroxene-anorthite-rich with igneous textures, 1.6%) Ca,Al-rich inclusions (CAIs), pyroxene-hibonite spherules (0.5%), and amoeboid olivine aggregates (AOAs, 30.2%). Melilite in some CAIs is replaced by spinel and Al-diopside and/or by anorthite, whereas spinel-pyroxene assemblages in CAIs and AOAs appear to be replaced by anorthite. Forsterite grains in several AOAs are replaced by low-Ca pyroxene. None of the CAIs or AOAs show evidence for Fe-alkali metasomatic or aqueous alteration. The mineralogy, textures, and bulk chemistry of most Acfer 094 refractory inclusions are consistent with their origin by gas-solid condensation and may reflect continuous interaction with SiO and Mg of the cooling nebula gas. It appears that only a few CAIs experienced subsequent melting. The Al-rich chondrules (ARCs; >10 wt% bulk Al₂O₃) consist of forsteritic olivine and low-Ca pyroxene phenocrysts, pigeonite, augite, anorthitic plagioclase, ± spinel, FeNi-metal, and crystalline mesostasis composed of plagioclase, augite and a silica phase. Most ARCs are spherical and mineralogically uniform, but some are irregular in shape and heterogeneous in mineralogy, with distinct ferromagnesian and aluminous domains. The ferromagnesian domains tend to form chondrule mantles, and are dominated by low-Ca pyroxene and forsteritic olivine, anorthitic mesostasis, and Fe,Ni-metal nodules. The aluminous domains are dominated by anorthite, high-Ca pyroxene and spinel, occasionally with inclusions of perovskite; have no or little FeNi-metal; and tend to form cores of the heterogeneous chondrules. The cores are enriched in bulk Ca and Al, and apparently formed from melting of CAI-like precursor material that did not mix completely with adjacent ferromagnesian melt. The inferred presence of CAI-like material among precursors for Al-rich chondrules is in apparent conflict with lack of evidence for melting of CAIs that occur outside chondrules, suggesting that these CAIs were largely absent from chondrule-forming region(s) at the time of chondrule formation. This may imply that there are several populations of CAIs in Acfer 094 and that mixing of “normal” CAIs that occur outside chondrules and chondrules that accreted into the Acfer 094 parent asteroid took place after chondrule formation. Alternatively, there may have been an overlap in the CAI- and chondrule-forming regions, where the least refractory CAIs were mixed with Fe-Mg chondrule precursors. This hypothesis is difficult to reconcile with the lack of evidence of melting of AOAs which represent aggregates of the least refractory CAIs and forsterite grains.     

Pyrite oxidation in moist air

The rate of pyrite oxidation in moist air was determined by measuring, over time, the pressure difference between a sealed chamber containing pyrite plus oxygen and a control. The experiments carried out at 25°C, 96.7% fixed relative humidity, and oxygen partial pressures of 0.21, 0.61, and 1.00 atm showed that the rate of oxygen consumption is a function of oxygen partial pressure and time. The rates of oxygen consumption (r, mol/m²sec) fit the expression

(A)     

Oxygen-exchange pathways in aluminum polyoxocations

Using molecular dynamics simulations and electronic structure methods, we postulate a mechanism to explain the complicated reactivity trends that are observed for oxygen isotope exchange reactions between sites in aluminum polyoxocations of the ε-Keggin type and bulk solution. Experimentally, the molecules have four nonequivalent oxygens that differ considerably in reactivity both within a molecule, and between molecules in the series: Al₁₃, GaAl₁₂, and GeAl₁₂ [MO₄Al₁₂(OH)₂₄(H₂O)₁₂ⁿ⁺(aq); with M = Al(III) for Al₁₃, n = 7; M = Ga(III) for GaAl₁₂, n = 7; M = Ge(IV) for GeAl₁₂, n = 8]. We find that a partly dissociated, metastable intermediate molecule of expanded volume is necessary for exchange of both sets of μ₂-OH and that the steady-state concentration of this intermediate reflects the bond strengths between the central metal and the μ₄-O. Thus the central metal exerts extraordinary control over reactions at hydroxyl bridges, although these are three bonds away.This mechanism not only explains the reactivity trends for oxygen isotope exchange in μ₂-OH and η-OH₂ sites in the ε-Keggin aluminum molecules, but also explains the observation that the reactivities of minerals tend to reflect the presence of highly coordinated oxygens, such as the μ₄-O in boehmite, α-, and γ-Al₂O₃ and their Fe(III) analogs. The partial dissociation of these highly coordinated oxygens, coupled with simultaneous activation and displacement of neighboring metal centers, may be a fundamental process by which metals atoms undergo ligand exchanges at mineral surfaces.     

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

Carbon and hydrogen isotope fractionation by moderately thermophilic methanogens

A series of laboratory studies were conducted to increase understanding of stable carbon (¹³C/¹²C) and hydrogen (D/H) isotope fractionation arising from methanogenesis by moderately thermophilic acetate- and hydrogen-consuming methanogens. Studies of the aceticlastic reaction were conducted with two closely related strains of Methanosaeta thermophila. Results demonstrate a carbon isotope fractionation of only 7‰ (α = 1.007) between the methyl position of acetate and the resulting methane. Methane formed by this process is enriched in ¹³C when compared with other natural sources of methane; the magnitude of this isotope effect raises the possibility that methane produced at elevated temperature by the aceticlastic reaction could be mistaken for thermogenic methane based on carbon isotopic content. Studies of H₂/CO₂ methanogenesis were conducted with Methanothermobacter marburgensis. The fractionation of carbon isotopes between CO₂ and CH₄ was found to range from 22 to 58‰ (1.023 ≤ α ≤ 1.064). Greater fractionation was associated with low levels of molecular hydrogen and steady-state metabolism. The fractionation of hydrogen isotopes between source H₂O and CH₄ was found to range from 127 to 275‰ (1.16 ≤ α ≤ 1.43). Fractionation was dependent on growth phase with greater fractionation associated with later growth stages. The maximum observed fractionation factor was 1.43, independent of the δD-H₂ supplied to the culture. Fractionation was positively correlated with temperature and/or metabolic rate. Results demonstrate significant variability in both hydrogen and carbon isotope fractionation during methanogenesis from H₂/CO₂. The relatively small fractionation associated with deuterium during H₂/CO₂ methanogenesis provides an explanation for the relatively enriched deuterium content of biogenic natural gas originating from a variety of thermal environments. Results from these experiments are used to develop a hypothesis that differential reversibility in the enzymatic steps of the H₂/CO₂ pathway gives rise to variability in the observed carbon isotope fractionation. Results are further used to constrain the overall efficiency of electron consumption by way of the hydrogenase system in M. marburgensis, which is calculated to be less than 55%.     

Optimal methods for estimating kinetic isotope effects from different forms of the Rayleigh distillation equation

The interpretation of stable isotope values hinges on precise, accurate estimates of kinetic isotope effects (α), which are equal to ¹k/²k, where ¹k and ²k are the reaction rate constants for the two isotopes. Kinetic isotope effects are commonly determined by monitoring the reactant concentration (C) and isotope ratio (R) relative to their initial values (C₁ and R₁, respectively). Values of α are estimated from the C and R values by using the Rayleigh distillation equation (RDE).

(A)     

Determination of the mass-dependence of cadmium isotope fractionation during evaporation

Mass fractionation laws relate the fractionation factor α_A for one isotope ratio to the fractionation factor α_B for a second isotope ratio of the same element, with a fractionation exponent β such that α_A = α_B^β. The exponent β defines the mass-dependence of the mass fractionation law and thus determines the slope of a mass fractionation line in linearized three isotope space. The generalized power law (GPL) defines β as a function of a variable exponent n. The laws that aim to describe equilibrium and kinetic isotope fractionations are special cases of the GPL with n = −1 and n → 0, respectively.Large isotope fractionations (up to 10% for ¹⁰⁶Cd/¹¹⁴Cd) were found to accompany the evaporation of molten Cd into vacuum at about 180°C. The slopes of the fractionation lines (β-values) were obtained by analyzing the Cd isotope compositions of the evaporation residues relative to the starting material with two different multiple collector-ICPMS instruments. For the most fractionated sample, the difference between the theoretical β-values, that describe kinetic and equilibrium isotope fractionation, is 10 to 20 times larger than the measurement uncertainty. A mass-dependence with n = −0.35 was determined for this sample. This result differs significantly from the value that would be expected for simple kinetic evaporation (n → 0), which is governed by the diffusion of monatomic Cd from the melt into vacuum. The observed “non-kinetic” mass-dependence probably results from partial recondensation (back reaction) of Cd vapor into the melt phase. This interpretation requires that equilibrium evaporation of Cd at about 180°C is associated with significant isotope fractionation.The present study demonstrates that the mechanism of isotope fractionation can be investigated by studying the associated mass-dependence, which can be determined by measuring the isotope ratios of a fractionated product relative to the starting material. The quantification of mass fractionation line slopes with the GPL should aid the interpretation of mass-dependent and small mass-independent isotope effects.     

Isotopic exchange of carbon-bound hydrogen over geologic timescales

The increasing popularity of compound-specific hydrogen isotope (D/H) analyses for investigating sedimentary organic matter raises numerous questions about the exchange of carbon-bound hydrogen over geologic timescales. Important questions include the rates of isotopic exchange, methods for diagnosing exchange in ancient samples, and the isotopic consequences of that exchange. This article provides a review of relevant literature data along with new data from several pilot studies to investigate such issues. Published experimental estimates of exchange rates between organic hydrogen and water indicate that at warm temperatures (50-100°C) exchange likely occurs on timescales of 10⁴ to 10⁸ yr. Incubation experiments using organic compounds and D-enriched water, combined with compound-specific D/H analyses, provide a new and highly sensitive method for measuring exchange at low temperatures. Comparison of δD values for isoprenoid and n-alkyl carbon skeletons in sedimentary organic matter provides no evidence for exchange in young (<1 Ma), cool sediments, but strong evidence for exchange in ancient (>350 Ma) rocks. Specific rates of exchange are probably influenced by the nature and abundance of organic matter, pore-water chemistry, the presence of catalytic mineral surfaces, and perhaps even enzymatic activity.Estimates of equilibrium fractionation factors between organic H and water indicate that typical lipids will be depleted in D relative to water by ∼75 to 140‰ at equilibrium (30°C). Thus large differences in δD between organic molecules and water cannot be unambiguously interpreted as evidence against hydrogen exchange. A better approach may be to use changes in stereochemistry as a proxy for hydrogen exchange. For example, estimated rates of H exchange in pristane are similar to predicted rates for stereochemical inversion in steranes and hopanes. The isotopic consequences of this exchange remain in question. Incubations of cholestene with D₂O indicate that the number of D atoms incorporated during structural rearrangements can be far less than the number of C-H bonds that are broken. Sample calculations indicate that, for steranes in immature sediments, the D/H ratio imparted by biosynthesis may be largely preserved in spite of significant structural changes.     

An exotic kind of cosmic material: graphite-containing xenoliths from the Krymka (LL3.1) chondrite

Seven graphite-containing xenoliths were found in the Krymka (LL3.1) chondrite. The xenoliths have the following chemical and mineralogical characteristics which distinguish them from the Krymka host: (1) low totals in bulk chemical analyses obtained by electron microprobe; (2) high bulk Fe abundances; (3) a uniform recrystallized, chondrule-free texture; (4) the presence of euhedral graphite and carbon-rich material; (5) higher quantities of troilite and metal; (6) a relatively homogeneous composition of silicates; (7) a distinctive composition of metal, chromite and phosphate; (8) isotopically heavy C in graphite compared to both bulk Krymka and graphite in other ordinary chondrites. The xenoliths are mineralogically similar, but not identical, to the Krymka carbonaceous clast K1, which bears graphite microcrystals, organic compounds and mysterite. They resemble carbonaceous chondrites, both chemically and isotopically. The mineralogical, chemical and isotopic data for the graphite-containing fragments suggest that this material represents metamorphosed varieties of a previously unknown type of unequilibrated carbonaceous matter. Most likely, the graphite has a metamorphic origin and was crystallized from C-containing precursor materials through the following transformation sequence: organic compounds → C-rich material → graphite.     

Formation and stability of schwertmannite in acidic mining lakes

Schwertmannite (ideal formula: Fe₈O₈(OH)₆SO₄) is typically found as a secondary iron mineral in pyrite oxidizing environments. In this study, geochemical constraints upon its formation are established and its role in the geochemical cycling of iron between reducing and oxidizing conditions are discussed. The composition of surface waters was analyzed and sediments characterized by X-ray diffraction, FTIR spectroscopy and determination of the Fe:S ratio in the oxalate extractable fraction from 18 acidic mining lakes. The lakes are exposed to a permanent supply of pyritegenous ferrous iron from adjacent ground water. In 3 of the lakes the suspended matter was fractionated using ultra filtration and analyzed with respect to their mineral composition. In addition, stability experiments with synthetic schwertmannite were performed. The examined lake surface waters were O₂-saturated and have sulfate concentrations (10.3 ± 5.5 mM) and pH values (3.0 ± 0.6) that are characteristic for the stability window of schwertmannite. Geochemical modeling implied that i) the waters were saturated with respect to schwertmannite, which controlled the activity of Fe³⁺ and sulfate, and ii) a redox equilibrium exists between Fe²⁺ and schwertmannite. In the uppermost sediment layers (1 to 5 cm depth), schwertmannite was detectable in 16 lakes—in 5 of them by all three methods. FTIR spectroscopy also proved its occurrence in the colloidal fraction (1-10 kDa) in all of the 3 investigated lake surface waters. The stability of synthetic schwertmannite was examined as a function of pH (2-7) by a 1-yr experiment. The transformation rate into goethite increased with increasing pH. Our study suggests that schwertmannite is the first mineral formed after oxidation and hydrolysis of a slightly acidic (pH 5-6), Fe(II)-SO₄ solution, a process that directly affects the pH of the receiving water. Its occurrence is transient and restricted to environments, such as acidic mining lakes, where the coordination chemistry of Fe³⁺ is controlled by the competition between sulfate and hydroxy ions (i.e. mildly acidic).