Influence of redox conditions on iodide migration through a deep clay formation (Toarcian argillaceous rock,Tournemire, France)

With a half-life of 15.7 Ma, a high mobility and the potential to accumulate in the biosphere, ¹²⁹I is considered, in safety assessment calculations for radioactive waste repositories, to be one of the main contributors to the radiological dose. Several authors have reported that, at low concentration, I⁻ is weakly retained on argillaceous rocks. This process is not yet well-understood and different hypotheses have been put forward as to whether reactive phases or experimental artifacts (e.g. pyrite oxidation) could be the reason for the retention of I⁻ observed at low concentration. The aim of this study was to investigate the effect on I⁻ mobility of (i) the redox conditions and (ii) the amount of pyrite and natural organic matter (NOM) contents of the rock. These questions were addressed by performing batch sorption, through-diffusion and out-diffusion experiments on rock samples of Toarcian argillaceous rock from Tournemire (Aveyron, France). One of the challenges faced during this study was to distinguish actual transport properties from experimental artifacts. A especially elaborate experimental set-up allowed limiting the (i) oxidation of both argillaceous rock and I⁻, and (ii) carbonate precipitation. A comparison of the batch sorption results obtained for two Toarcian clay specimens, that differed in their amount of pyrite and NOM, allowed relating I⁻ sorption to pyrite oxidation. However, no evidence was found to associate the I⁻ behavior to the NOM amounts. While the through-diffusion experiments showed a very slight sorption (distribution ratio (R_d) = 0.016 mL g⁻¹) for the lowest I⁻ concentration under oxic conditions, the out-diffusion tests performed after the through-diffusion experiments on the same cells showed significant sorption under both oxic and anoxic conditions, resulting in R_d ranging from 0.02 mL g⁻¹ to 1.25 mL g⁻¹. The range of R_d values was higher for the upstream reservoir under oxic conditions. The discrepancies observed between the through-diffusion and the out-diffusion experiments suggest a kinetic control of the I⁻ uptake by argillaceous rocks under oxic and anoxic conditions.     

O and H diffusion in uraninite: Implications for fluid-uraninite interactions, nuclear waste disposal, and nuclear forensics

Diffusion coefficients for oxygen and hydrogen were determined from a series of natural uraninite-H₂O experiments between 50 and 700 °C. Under hydrous conditions there are two diffusion mechanisms: (1) an initial extremely fast-path diffusion mechanism that overprinted the oxygen isotopic composition of the entire crystals regardless of temperature and (2) a slower volume-diffusive mechanism dominated by defect clusters that displace or eject nearest neighbor oxygen atoms to form two interstitial sites and two partial vacancies, and by vacancy migration. Using the volume diffusion coefficients in the temperature range of 400-600 °C, diffusion coefficients for oxygen can be represented by D = 1.90e⁻⁵ exp (−123,382 J/RT) cm²/s and for temperatures between 100 and 300 °C the diffusion coefficients can be represented by D = 1.95e⁻¹⁰ exp (−62484 J/RT) cm²/s, where the activation energies for uraninite are 123.4 and 62.5 kJ/mol, respectively. Hydrogen diffusion in uraninite appears to be controlled by similar mechanisms as oxygen. Using the volume diffusion coefficients for temperatures between 50 and 700 °C, diffusion coefficients for hydrogen can be represented by D = 9.28e⁻⁶ exp (−156,528 J/RT) cm²/s for temperatures between 450 and 700 °C and D = 1.39e⁻¹⁴ exp (−34518 J/RT) cm²/s for temperatures between 50 and 400 °C, where the activation energies for uraninite are 156.5 and 34.5 kJ/mol, respectively.Results from these new experiments have implications for isotopic exchange during natural UO₂-water interactions. The exceptionally low δ¹⁸O values of natural uraninites (i.e. −32‰ to −19.5‰) from unconformity-type uranium deposits in Saskatchewan, in conjunction with theoretical and experimental uraninite-water and UO₃-water fractionation factors, suggest that primary uranium mineralization is not in oxygen isotopic equilibrium with coeval clay and silicate minerals. The low δ¹⁸O values have been interpreted as resulting from the low temperature overprinting of primary uranium mineralization in the presence of relatively modern meteoric fluids having δ¹⁸O values of ca. −18‰, despite petrographic and U-Pb isotope data that indicate limited alteration. Our data show that the anomalously low oxygen isotopic composition of the uraninite from the Athabasca Basin can be due to meteoric water overprinting under reducing conditions, and meteoric water or groundwater can significantly affect the oxygen isotopic composition of spent nuclear fuel in a geologic repository, with minimal change to the chemical composition or texture. Moreover, the rather fast oxygen and hydrogen diffusion coefficients for uraninite, especially at low temperatures, suggest that oxygen and hydrogen diffusion may impart characteristic isotopic signals that can be used to track the route of fissile material.     

Diffusion and sorption of Cs, I and HTO in samples of the argillaceous Wakkanai Formation from the Horonobe URL, Japan: Clay-based modeling approach

Diffusion and sorption behaviors of cationic Cs⁺, anionic I⁻ and neutral HTO in samples of the Wakkanai Formation from the Horonobe underground research laboratory (URL), Japan, were investigated as a function of ionic strength (I) of synthetic groundwater by through-diffusion and batch sorption experiments and mechanistic modeling. The effective diffusivities (D_e) measured by through-diffusion experiments showed cation excess and anion exclusion effects, which were strongly dependent on I; D_e for Cs⁺ decreased as I increased, D_e for I⁻ showed the opposite dependency and D_e for HTO showed no dependence. The sorption of Cs⁺ measured by through-diffusion and batch sorption experiments were described by Freundlich isotherms with consistent parameters and decreased with I as a result of competitive ion exchange.Diffusion and sorption behaviors were interpreted by assuming the clay components of illite and smectite control diffusion and sorption mechanisms. The component additive (CA) sorption model, which includes illite and smectite contents and their ion exchange constants, provided a reasonable account of the Cs⁺ sorption trends measured as functions of I and Cs concentration. The diffusion model was developed by coupling the electrical double layer (EDL) model, describing the change of ionic concentrations (cation excess and anion deficit) and viscoelectric effects caused by electrostatic interaction at negatively charged clay surfaces, and a simplified pore model assuming one type of pore shape and includes their size distribution. When averaging the electrostatic effects by using the pore surface area distribution, the model could predict the cation excess and anion exclusion effects, and its dependence on I reasonably well. This result implies the nanoscale pores dominating the pore surface area can strongly impact on ionic diffusion in argillaceous rocks. The clay-based modeling approach described here provides a useful tool to predict ionic diffusion and sorption in argillaceous rocks.     

Tracing interactions between natural argillites and hyper-alkaline fluids from engineered cement paste and concrete: Chemical and isotopic monitoring of a 15-years old deep-disposal analogue

Samples of Toarcian argillite were collected both next to and far from a CEM II cement paste and a CEM II concrete, within the specific context of a 15-a old borehole located in the Tournemire Experimental Platform (Aveyron, France). The objectives were evaluation of the mineralogical and geochemical changes of the claystone at the contact with the cementitious materials and determination of the spatial extent of the interactions. The approach includes the examination of the mineralogical (XRD, SEM, TEM), chemical (major, trace, rare earth elements) and isotopic (Sr, C, O) compositions of argillite whole-rocks and of various soluble phases, at two scales: in the rock matrix (P1 scale) and along micro-cracks (P2 scale). The two study scales outline nearly similar mineralogical modifications, shown by the presence of Ca silicate hydrates (C–S–H) and newly-formed CaCO₃ within 10–15 mm of the cement paste and concrete. Chemical data from whole-rock argillites indicate few changes in a slightly thicker zone (18–20 mm), mainly consisting of an increase in the CaO wt.%, and a decrease in Sr contents. The other elementary contents remained quite constant except for MgO, which suggests redistribution with precipitation of a Mg-rich mineral phase at 20 mm from cement paste/concrete interface. Acetic acid leachates had more pronounced variations, including a decrease of the total elementary content in the same ‘geochemical disturbed zone’ (GDZ), together with a significant increase of the Sr isotopic ratios. A combination of Sr and C/O isotopic patterns was used to distinguish the behavior of secondary cementitious phases in the clay-rich rock: (i) calcite dissolution and re-precipitation is supported by C/O isotopic data and (ii) C–S–H neoformation is evidenced by the ⁸⁷Sr/⁸⁶Sr ratios; this tool also contributes to determine the origin of the fluids. The proportion of newly-formed C–S–H in the matrix and in the micro-cracks of the argillite is modeled.     

Permian volcanic rocks from the Apuseni Mountains (Romania): Geochemistry and tectonic constraints

An Early Permian volcanic assemblage is well exposed in the central-western part of the Apuseni Mountains (Romania). The rocks are represented by rhyolites, basalts and subordinate andesites suggesting a bimodal volcanic activity that is intimately associated with a post-orogenic (Variscan) syn-sedimentary intra-basinal continental molasse sequences. The mafic and mafic-intermediate rocks belong to sub-alkaline tholeiitic series were separated in three groups (I–III) showing a high Th and Pb abundances, depletion in Nb, Ta and Sr, and slightly enriched in LREE patterns (La_N/Yb_N = 1.4–4.4). Isotopically, the rocks of Group I have the initial ratios ⁸⁷Sr/⁸⁶Sr_(i) = 0.709351–0.707112, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512490–0.512588 and high positive ?_Nd270 values from 3.9 to 5.80; the rocks of Group II present for the initial ratios values ⁸⁷Sr/⁸⁶Sr_(i) = 0.709434–0.710092, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512231–0.512210 and for ?_Nd270 the negative values from −1.17 to −1.56; the rocks of Group III display for the initial ratios the values ⁸⁷Sr/⁸⁶Sr_(i) = 0.710751–0.709448, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512347–0.512411 and for ?_Nd270 the positive values from 1.64 to 2.35. The rocks resembling continental tholeiites, suggest a mantle origin and were further affected by fractionation and crustal contamination. In addition, the REE geochemistry (1 > Sm_N/Yb_N < 2.5; 0.9 > La_N/Sm_N < 2.5) suggests that these rocks were generated by high percentage partial melting of a metasomatized mantle in the garnet peridotite facies. The felsic rocks are enriched in Cs, Rb Th and U and depleted in Nb, Ta, Sr, Eu, and Ti. The REE fractionation patterns show a strong negative Eu anomaly (Eu/Eu* = 0.23–0.40). The felsic rocks show the initial ratios the values: ⁸⁷Sr/⁸⁶Sr_(i) = 0.704096–0.707805, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512012–0.512021 and for ?_Nd270 the negative values from −5.27 to −5.44. They suggest to be generated within the lower crust during the emplacement of mantle-derived magmas that provided necessary heat to crustal partial melting.     

Effect of confining pressure on the diffusion of HTO, 36Cl− and 125I− in a layered argillaceous rock (Opalinus Clay): diffusion perpendicular to the fabric

The through- and out-diffusion of HTO, ³⁶Cl⁻ and ¹²⁵I⁻ in Opalinus Clay, an argillaceous rock from the northern part of Switzerland, was studied under different confining pressures between 4 and 15 MPa. The direction of diffusion and the confining pressure were perpendicular to the bedding. Confining pressure had only a small effect on diffusion. An increase in pressure from 4 to 15 MPa resulted in a decrease of the effective diffusion coefficient of ∼20%. Diffusion accessible porosities were not measurably affected. The values of the effective diffusion coefficients, D_e, ranged between (5.6±0.4)×10⁻¹² and (6.7±0.4)×10⁻¹² m² s⁻¹ for HTO, (7.1±0.5)×10⁻¹³ and (9.1±0.6)×10⁻¹³ m² s⁻¹ for ³⁶Cl⁻ and (4.5±0.3)×10⁻¹³ and (6.6±0.4)×10⁻¹³ m² s⁻¹ for ¹²⁵I⁻. The rock capacity factors, α, measured were circa 0.14 for HTO, 0.040 for ³⁶Cl⁻ and 0.080 for ¹²⁵I⁻. Because of anion exclusion effects, anions diffuse slower and exhibit smaller diffusion accessible porosities than the uncharged HTO. Unlike ³⁶Cl⁻, ¹²⁵I⁻ sorbs weakly on Opalinus Clay resulting in a larger rock capacity factor. The sorption coefficient, K_d, for ¹²⁵I⁻ is of the order of 1–2×10⁻⁵ m³ kg⁻¹. The effective diffusion coefficient for HTO is in good agreement with values measured in other sedimentary rocks and can be related to the porosity using Archie's Law with exponent m=2.5.     

The kinetics of iodide oxidation by the manganese oxide mineral birnessite

The kinetics of iodide (I⁻) and molecular iodine (I₂) oxidation by the manganese oxide mineral birnessite (δ-MnO₂) was investigated over the pH range 4.5-6.25. I⁻ oxidation to iodate proceeded as a two-step reaction through an I₂ intermediate. The rate of the reaction varied with both pH and birnessite concentration, with faster oxidation occurring at lower pH and higher birnessite concentration. The disappearance of I⁻ from solution was first order with respect to I⁻ concentration, pH, and birnessite concentration, such that −d[I⁻]/dt = k[I⁻][H⁺][MnO₂], where k, the third order rate constant, is equal to 1.08 ± 0.06 × 10⁷ M⁻² h⁻¹. The data are consistent with the formation of an inner sphere I⁻ surface complex as the first step of the reaction, and the adsorption of I⁻ exhibited significant pH dependence. Both I₂, and to a lesser extent, sorbed to birnessite. The results indicate that iodine transport in mildly acidic groundwater systems may not be conservative. Because of the higher adsorption of the oxidized I species I₂ and , as well as the biophilic nature of I₂, redox transformations of iodine must be taken into account when predicting I transport in aquifers and watersheds.     

Tissue N content and N natural abundance in epilithic mosses for indicating atmospheric N deposition in the Guiyang area,SW China

Tissue N contents and δ¹⁵N signatures in 175 epilithic mosses were investigated from urban to rural sites in Guiyang (SW China) to determine atmospheric N deposition. Moss N contents (0.85–2.97%) showed a significant decrease from the urban area (mean = 2.24 ± 0.32%, 0–5 km) to the rural area (mean = 1.27 ± 0.13%, 20–25 km), indicating that the level of N deposition decreased away from the urban environment, while slightly higher N contents re-occurred at sites beyond 30 km, suggesting higher N deposition in more remote rural areas. Moss δ¹⁵N ranged from −12.50‰ to −1.39‰ and showed a clear bimodal distribution (−12‰ to −6‰ and −5‰ to −2‰), suggesting that there are two main sources for N deposition in the Guiyang area. More negative δ¹⁵N (mean = −8.87 ± 1.65‰) of urban mosses mainly indicated NH₃ released from excretory wastes and sewage, while the less negative δ¹⁵N (from −3.83 ± 0.82‰ to −2.48 ± 0.95‰) of rural mosses were mainly influenced by agricultural NH₃. With more negative values in the urban area than in the rural area, the pattern of moss δ¹⁵N variation in Guiyang was found to be opposite to cities where N deposition is dominated by NO_x–N. Therefore, NH_x–N is the dominant N form deposited in the Guiyang area, which is supported by higher NH_x–N than NO_x–N in local atmospheric deposition. From the data showing that moss is responding to NH_x–N/NO_x–N in deposition it can be further demonstrated that the variation of moss δ¹⁵N from the Guiyang urban to rural area was more likely controlled by the ratio of urban-NH_x/agriculture-NH_x than the ratio of NH_x–N/NO_x–N. The results of this study have extended knowledge of atmospheric N sources in city areas, showing that urban sewage discharge could be important in cities co-generic to Guiyang.     

The distribution of short-lived radioisotopes in the early solar system and the chronology of asteroid accretion, differentiation, and secondary mineralization

We evaluate initial (²⁶Al/²⁷Al)_I, (⁵³Mn/⁵⁵Mn)_I, and (¹⁸²Hf/¹⁸⁰Hf)_I ratios, together with ²⁰⁷Pb/²⁰⁶Pb ages for igneous differentiated meteorites and chondrules from ordinary chondrites for consistency with radioactive decay of the parent nuclides within a common, closed isotopic system, i.e., the early solar nebula. The relative initial isotopic abundances of ²⁶Al, ⁵³Mn, and ¹⁸²Hf in differentiated meteorites and chondrules are consistent with decay from common solar system initial values, here denoted by I(Al)_SS, I(Mn)_SS, and I(Hf)_SS, respectively. I(Mn)_SS and I(Hf)_SS = 9.1 ± 1.7 × 10⁻⁶ and 1.07 ± 0.08 × 10⁻⁴, respectively, correspond to “canonical” I(Al)_SS = 5.1 × 10⁻⁵. I(Hf)_SS so determined is consistent with I(Hf)_SS = 9.72 ± 0.44 × 10⁻⁵ directly determined from an internal Hf-W isochron for CAI minerals. I(Mn)_SS is within error of the lowest value directly measured for CAIs. We suggest that erratically higher values measured for CAIs in carbonaceous chondrites may reflect proton irradiation of unaccreted CAIs by the early Sun after other asteroids destined for melting by ²⁶Al decay had already accreted. The ⁵³Mn incorporated within such asteroids would have been shielded from further “local” spallogenic contributions from within the solar system. The relative initial isotopic abundances of the short-lived nuclides are less consistent with the ²⁰⁷Pb/²⁰⁶Pb ages of the corresponding materials than with one another. The best consistency of short- and long-lived chronometers is obtained for (¹⁸²Hf/¹⁸⁰Hf)_I and the ²⁰⁷Pb/²⁰⁶Pb ages of angrites. (¹⁸²Hf/¹⁸⁰Hf)_I decreases with decreasing ²⁰⁷Pb/²⁰⁶Pb ages at the rate expected from the 8.90 ± 0.09 Ma half-life of ¹⁸²Hf. The model solar system age thus determined is T_SS,Hf-W = 4568.3 ± 0.7 Ma. (²⁶Al/²⁷Al)_I and (⁵³Mn/⁵⁵Mn)_I are less consistent with ²⁰⁷Pb/²⁰⁶Pb ages of the corresponding meteorites, but yield T_SS,Mn-Cr = 4568.2 ± 0.5 Ma relative to I(Al)_SS = 5.1 × 10⁻⁵ and a ²⁰⁷Pb/²⁰⁶Pb age of 4558.55 ± 0.15 Ma for the LEW86010 angrite. The Mn-Cr method with I(Mn)_SS = 9.1 ± 1.7 × 10⁻⁶ is useful for dating accretion (if identified with chondrule formation), primary igneous events, and secondary mineralization on asteroid parent bodies. All of these events appear to have occurred approximately contemporaneously on different asteroid parent bodies. For I(Mn)_SS = 9.1 ± 1.7 × 10⁻⁶, parent body differentiation is found to extend at least to ∼5 Ma post-T_SS, i.e., until differentiation of the angrite parent body ∼4563.5 Ma ago, or ∼4564.5 Ma ago using the directly measured ²⁰⁷Pb/²⁰⁶Pb ages of the D’Orbigny-clan angrites. The ∼1 Ma difference is characteristic of a remaining inconsistency for the D’Orbigny-clan between the Al-Mg and Mn-Cr chronometers on one hand, and the ²⁰⁷Pb/²⁰⁶Pb chronometer on the other. Differentiation of the IIIAB iron meteorite and ureilite parent bodies probably occurred slightly later than for the angrite parent body, and at nearly the same time as one another as shown by the Mn-Cr ages of IIIAB irons and ureilites, respectively. The latest recorded episodes of secondary mineralization are for carbonates on the CI carbonaceous chondrite parent body and fayalites on the CV carbonaceous chondrite parent body, both extending to ∼10 Ma post-T_SS.     

Elemental and isotopic fractionation of Type B CAI-like liquids by evaporation

Vacuum evaporation experiments with Type B CAI-like starting compositions were carried out at temperatures of 1600, 1700, 1800, and 1900 °C to determine the evaporation kinetics and evaporation coefficients of silicon and magnesium as a function of temperature as well as the kinetic isotope fractionation factor for magnesium. The vacuum evaporation kinetics of silicon and magnesium are well characterized by a relation of the form J = J_oe^−E/RT with J_o = 4.17 × 10⁷ mol cm⁻² s⁻¹, E = 576 ± 36 kJ mol⁻¹ for magnesium, J_o = 3.81 × 10⁶ mol cm⁻² s⁻¹, E = 551 ± 63 kJ mol⁻¹ for silicon. These rates only apply to evaporation into vacuum whereas the actual Type B CAIs were almost certainly surrounded by a finite pressure of a hydrogen-dominated gas. A more general formulation for the evaporation kinetics of silicon and magnesium from a Type B CAI-like liquid that applies equally to vacuum and conditions of finite hydrogen pressure involves combining our determinations of the evaporation coefficients for these elements as a function of temperature (γ = γ₀e^−E/RT with γ₀ = 25.3, E = 92 ± 37 kJ mol⁻¹ for γ_Si; γ₀ = 143, E = 121 ± 53 kJ mol⁻¹ for γ_Mg) with a thermodynamic model for the saturation vapor pressures of Mg and SiO over the condensed phase. High-precision determinations of the magnesium isotopic composition of the evaporation residues from samples of different size and different evaporation temperature were made using a multicollector inductively coupled plasma mass spectrometer. The kinetic isotopic fractionation factors derived from this data set show that there is a distinct temperature effect, such that the isotopic fractionation for a given amount of magnesium evaporated is smaller at lower temperature. We did not find any significant change in the isotope fractionation factor related to sample size, which we interpret to mean that recondensation and finite chemical diffusion in the melt did not affect the isotopic fractionations. Extrapolating the magnesium kinetic isotope fractionations factors from the temperature range of our experiments to temperatures corresponding to partially molten Type B CAI compositions (1250-1400 °C) results in a value of α_Mg ≈ 0.991, which is significantly different from the commonly used value of .     

Halogen diffusion in a basaltic melt

The diffusion of the halogens fluorine, chlorine and bromine was measured in a hawaiitic melt from Mt. Etna at 500 MPa and 1.0 GPa, 1250 to 1450 °C at anhydrous conditions; the diffusion of F and Cl in the melt was also studied with about 3 wt% of dissolved water. Experiments were performed using the diffusion-couple technique in a piston cylinder. Most experiments were performed with only one halogen diffusing between the halogen-enriched and halogen-poor halves of the diffusion couple, but a few experiments with a mixture of halogens (F, Cl and Br) were also performed in order to investigate the possibility of interactions between the halogens during diffusion. Fluorine and chlorine diffusivity show a very similar behavior, slightly diverging at low temperature. Bromine diffusion is a factor of about 2-5 lower than the other halogens in this study. Diffusion coefficients for fluorine range between 2.3 × 10⁻¹¹ and 1.4 × 10⁻¹⁰ m² s⁻¹, for chlorine between 1.1 × 10⁻¹¹ and 1.3 × 10⁻¹⁰ and for bromine between 9.4 × 10⁻¹² and 6.8 × 10⁻¹¹ m² s⁻¹. No pressure effect was detected at the conditions investigated. In experiments involving mixed halogens, the diffusivities appear to decrease slightly (by a factor of ∼3), and are more uniform among the three elements. However, activation energies for diffusion do not appear to differ between experiments with individual halogens or when they are all mixed together. The effect of water increases the diffusion coefficients of F and Cl by no more than a factor of 3 compared to the anhydrous melt (D_F = 4.0 × 10⁻¹¹ to 1.6 × 10⁻¹⁰ m² s⁻¹; D_Cl = 3.0 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ m² s⁻¹). Comparing our results to the diffusion coefficients of other volatiles in nominally dry basaltic melts, halogen diffusivities are about one order of magnitude lower than H₂O, similar to CO₂, and a factor of ∼5 higher than S. The contrasting volatile diffusivities may affect the variable extent of volatile degassing upon melt depressurization and vesiculation, and can help our understanding of the compositions of rapidly grown magmatic bubbles.     

Isotopic fractionation of the major elements of molten basalt by chemical and thermal diffusion

Samples produced in piston cylinder experiments were used to document the thermal isotopic fractionation of all the major elements of basalt except for aluminum and the fractionation of iron isotopes by chemical diffusion between a natural basalt and rhyolite. The thermal isotopic fractionations are summarized in terms of a parameter Ω_i defined as the fractionation in per mil per 100 °C per atomic mass units difference between the isotopes. For molten basalt we report Ω_Ca = 1.6, Ω_Fe = 1.1, Ω_Si = 0.6, Ω_O = 1.5. In an earlier paper we reported Ω_Mg = 3.6. These fractionations represent a steady state balance between thermal diffusion and chemical diffusion with the mass dependence of the thermal diffusion coefficient being significantly larger than the mass dependence of the chemical diffusion coefficients for isotopes of the same element. The iron isotopic measurements of the basalt-rhyolite diffusion couple showed significant fractionation that are parameterized in terms of a parameter β_Fe = 0.03 when the ratio of the diffusion coefficients D₅₄ and D₅₆ of ⁵⁴Fe and ⁵⁶Fe is expressed in terms of the atomic mass as D₅₄/D₅₆ = (56/54)^β_Fe. This value of β_Fe is smaller than what we had measured earlier for lithium, magnesium and calcium (i.e., β_Li = 0.215, β_Ca = 0.05, β_Mg = 0.05) but still significant when one takes into account the high precision with which iron isotopic compositions can be measured (i.e., ±0.03‰) and that iron isotope fractionations at magmatic temperatures from other causes are extremely small. In a closing section we discuss technological and geological applications of isotopic fractionations driven by either or both chemical and thermal gradients.     

Trace/minor element:calcium ratios in cultured benthic foraminifera. Part II: Ontogenetic variation

Ontogenetic (developmental stage) measurements of Mg/Ca and Sr/Ca were made on the benthic foraminifer Bulimina aculeata, which were cultured under controlled physicochemical conditions of temperature, pH, alkalinity, salinity, and trace- and minor-element concentrations. We utilized two methods of ontogenetic sampling—whole specimens progressively increasing in length and laser microdissection of a single specimen with subsequent analysis of dissected portions. A novel high-resolution laser-microdissection (HRLM) method allowed for precise (10 μm) cuts of the foraminiferal tests (shells) along the geometrically complex sutures distinguishing individual chambers. This new microdissection method limited sample loss and cross-contamination between foraminiferal chambers. Little or no variation in D_Sr was observed at different foraminiferal developmental stages. Conversely, D_Mg was enriched during a mid-developmental stage of whole-specimen samples (150-225 μm D_Mg = 1.6 × 10⁻³) compared to earlier and later stages (<150 μm, >225 μm D_Mg = 8.3 × 10⁻⁴). Further analysis of HRLM ontogenetic samples showed a larger, age-dependent D_Mg signature variation. This increase in shell Mg/Ca may contribute substantially to the measured inter-individual variability in Mg/Ca temperature prediction for cultured B. aculeata. Due to relatively large Mg/Ca inter- and intra-individual variability, measuring similar-size foraminiferal samples may improve the precision of paleotemperature prediction. Additionally, partial dissolution of the highest ontogenetically Mg-enriched calcite (D_Mg = 1.3 × 10⁻²-1.6 × 10⁻²) may occur in undersaturated bottom-water environments or during reductive cleaning procedures. Thus, the calcite phases remaining after partial dissolution by either natural or laboratory cleaning processes may not accurately represent the calcification environment.     

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     

Hydration of rhyolitic glass during weathering as characterized by IR microspectroscopy

The mechanism and rate of hydration of rhyolitic glass during weathering were studied. Doubly polished thin sections of two rhyolites with different duration of weathering (Ohsawa lava: 26,000 yr, Awanomikoto lava: 52,000 yr) were prepared. Optical microscope observation showed that altered layers had developed along the glass surfaces. IR spectral line profile analysis was conducted on the glass sections from the surface to the interior for a length of 250 μm and the contents of molecular H₂O (H₂O_m), OH species (OH) and total water (H₂O_t) were determined. The diffusion profile of H₂O_m in Ohsawa lava extends beyond the layer observed by optical microscope. The content of H₂O_m in the hydrated region is much higher than that of OH species. Thus, the reaction from H₂O_m to OH appears to be little and H₂O_m is the dominant water species moving into the glass during weathering. Based on the concentration profiles, the diffusion coefficients of H₂O_m(D_H2Om) and H₂O_t(D_H2Ot) were determined to be 2.8 × 10⁻¹⁰ and 3.4 × 10⁻¹⁰ μm² s⁻¹ for Ohsawa lava, and 5.2 × 10⁻¹¹ and 4.1 × 10⁻¹¹ μm² s⁻¹ for Awanomikoto lava, respectively. The obtained D_H2Om during weathering are more than 2-3 orders of magnitude larger than the diffusion coefficient at ∼20 °C that is extrapolated from the diffusivity data for >400 °C. This might suggest that the mechanism of water transport is different at weathering conditions and >400 °C.     

Experimental constraints on the mobility of Rhenium in silicate liquids

The volatization of Rhenium (Re) from melts of natural basalt, dacite and a synthetic composition in the CaO-MgO-Al₂O₃-SiO₂ system has been investigated at 0.1 MPa and 1250-1350 °C over a range of fO₂ conditions from log fO₂ = −10 to −0.68. Experiments were conducted using open top Pt crucibles doped with Re and Yb. Analysis of quenched glasses by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) normal to the melt/gas interface showed concentration profiles for Re, to which a semi-infinite one-dimensional diffusion model could be applied to extract diffusion coefficients (D). The results show Re diffusivity in basalt at 1300 °C in air is log D_Re = −7.2 ± 0.3 cm²/s and increases to log D_Re = −6.6 ± 0.3 cm²/s when trace amounts of Cl were added to the starting material. At fO₂ conditions below the nickel-nickel oxide (NNO) buffer Re diffusivity decreases to and to in dacitic melt. In the CMAS composition, . The diffusivity of Re is comparable to Ar and CO₂ in basalt at 500 MPa favoring its release as a volatile. Our results support the contention that subaerial degassing is the cause of lower Re concentrations in arc-type and ocean island basalts compared to mid-ocean ridge basalts.     

Seawater temperature proxies based on DSr, DMg, and DU from culture experiments using the branching coral Porites cylindrica

In order to investigate the incorporation of Sr, Mg, and U into coral skeletons and its temperature dependency, we performed a culture experiment in which specimens of the branching coral (Porites cylindrica) were grown for 1 month at three seawater temperatures (22, 26, and 30 °C). The results of this study showed that the linear extension rate of P. cylindrica has little effect on the skeletal Sr/Ca, Mg/Ca, and U/Ca ratios. The following temperature equations were derived: Sr/Ca (mmol/mol) = 10.214(±0.229) − 0.0642(±0.00897) × T (°C) (r² = 0.59, p < 0.05); Mg/Ca (mmol/mol) = 1.973(±0.302) + 0.1002(±0.0118) × T (°C) (r² = 0.67, p < 0.05); and U/Ca (μmol/mol) = 1.488(±0.0484) − 0.0212(±0.00189) × T (°C) (r² = 0.78, p < 0.05). We calculated the distribution coefficient (D) of Sr, Mg, and U relative to seawater temperature and compared the results with previous data from massive Porites corals. The seawater temperature proxies based on D calibrations of P. cylindrica established in this study are generally similar to those for massive Porites corals, despite a difference in the slope of D_U calibration. The calibration sensitivity of D_Sr, D_Mg, and D_U to seawater temperature change during the experiment was 0.64%/°C, 1.93%/°C, and 1.97%/°C, respectively. These results suggest that the skeletal Sr/Ca ratio (and possibly the Mg/Ca and/or U/Ca ratio) of the branching coral P. cylindrica can be used as a potential paleothermometer.     

Internally consistent thermodynamic data for magnesium sulfate hydrates

Experimental studies on the stability of several Mg-sulfate hydrates including epsomite (MgSO₄·7H₂O), hexahydrite (MgSO₄·6H₂O), starkeyite (MgSO₄·4H₂O), and kieserite (MgSO₄·H₂O) as a function of temperature and relative humidity are in poor agreement with calculations based on thermodynamic properties of these substances taken from the literature. Therefore, we synthesized four different MgSO₄ hydrates and measured their enthalpies of formation by solution calorimetry at T = 298.15 K. The resulting enthalpies of formation from the elements are:

•

Δ_fH⁰₂₉₈ (epsomite) = −3387.7 ± 1.3 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (hexahydrite) = −3088.1 ± 1.1 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (sanderite, MgSO₄·2H₂O) = −1894.9 ± 1.3 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (kieserite) = −1612.4 ± 1.3 kJmol⁻¹

Using mathematical programming (MAP) techniques, standard thermodynamic values consistent both with our calorimetric data and previously published humidity brackets could be derived:

•

Epsomite: Δ_fH⁰₂₉₈ = −3388.7 kJmol⁻¹, S⁰₂₉₈ = 371.3 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2871.0 kJmol⁻¹

•

Hexahydrite: Δ_fH⁰₂₉₈ = −3087.3 kJmol⁻¹, S⁰₂₉₈ = 348.5 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2632.3 kJmol⁻¹

•

Starkeyite: Δ_fH⁰₂₉₈ = −2496.1 kJmol⁻¹, S⁰₂₉₈ = 259.9 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2153.8 kJmol⁻¹

•

Kieserite: Δ_fH⁰₂₉₈ = −1611.5 kJmol⁻¹, S⁰₂₉₈ = 126.0 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −1437.9 kJmol⁻¹

Additionally, heat capacity measurements and standard entropy determinations of several magnesium sulfate hydrate minerals from the literature are analyzed and judged against estimates obtained from a linear combination of the heat capacities of MgSO₄ and hexagonal ice. The results of the MAP analysis are compared to these estimates to conclude that heat capacity and entropy correlate well with the number of waters of hydration. However, even the good correlation is not good enough to capture the fine variations in these properties. Consequently, their experimental measurement is inevitable if reliable thermodynamic data are sought. Our MAP thermodynamic data show that epsomite, hexahydrite, and kieserite have stability fields in the T-%RH space. Starkeyite is metastable. Although no MAP data could have been derived for pentahydrite (MgSO₄·5H₂O) and sanderite, their transient existence suggest that both of them are metastable as well.