Phosphoric acid fractionation factors for smithsonite and cerussite between 25 and 72°C

The intramolecular kinetic oxygen isotope fractionation between CO₂ and CO₃²⁻ during reaction of phosphoric acid with natural smithsonite (ZnCO₃) and cerussite (PbCO₃) has been determined between 25 and 72°C. While cerussite decomposes in phosphoric acid within a few hours at 25°C, smithsonite reacts very slowly with the acid at 25°C providing yields of CO₂ < 25% after 2 weeks. The low yields result in a low precision for oxygen isotope measurements of the acid-liberated CO₂ (±1.65‰, 1σ, n = 9). The yield and reproducibility of oxygen isotope values of the acid-liberated CO₂ from smithsonite can be improved, the latter to ∼±0.15‰, by increasing the reaction temperature to 50°C for 12 h or to 72°C for 1 h. Our new phosphoric acid fractionation factor for natural cerussite at 25°C deviates significantly from a previously published value on synthetic material. The temperature dependence of the oxygen isotope factionation factor, α between acid-liberated CO₂ and carbonate at 25 to 72°C is given by the following equations

     

Boron isotopic fractionation related to boron sorption on humic acid and the structure of surface complexes formed

Boron isotopic fractionation during adsorption onto Ca-flocculated Aldrich humic acid (HA) has been investigated experimentally as a function of solution pH at 25°C and I = 0.15 M. Boron aqueous concentration and isotopic composition were determined by Cs₂BO₂⁺ Positive Thermal Ionization Mass Spectrometry analysis, while the structure of B surface complexes on HA was characterized using ¹¹B Magic Angle Spinning Nuclear Magnetic Resonance (MAS NMR). Significant B sorption on HA was observed at 6 < pH < 12 with a maximum value of Kd, the partition coefficient between adsorbed and aqueous boron, equal to 40 at pH = 9.5-10. Combined ¹¹B MAS NMR analysis and FITEQL modeling of B sorption on HA showed that this element forms tetrahedrally coordinated five- or six-membered ring chelates, most likely 1,2-diol and 1,3-diol complexes at alkaline pH (8 < pH < 11) and dicarboxylic complexes at near neutral conditions (6 < pH < 9). Results of this study demonstrate for the first time that boron sorption on HA induces a strong pH-dependent isotope fractionation—with ¹¹B depleted at the surface of HA—that reaches a maximum at 5 < pH < 9 (α = 0.975, Δ = −25‰) and decreases sharply at pH >9. The measured isotope fractionation cannot be modeled assuming that the isotopic composition of the sorbed borate species is identical to that of B(OH)₄^- species in the parent solution. It is shown that the extent of isotopic fractionation depends not only on B aqueous speciation but also on the distribution and structure of the borate surface complexes formed. In agreement with energetic constrains, calculation of the isotope fractionation factors between aqueous boric acid and boron surface complexes suggests that the formation of the strained six-membered ring 1,3-diol complex yields a much higher fractionation (α_BLP1−III = 0.954-0.960, Δ = −41/-47‰) than that of the very stable five-membered ring 1,2-diol (α_BLP2−III = 0.983, Δ = −18‰). The results of this study open new perspectives to understand and model boron biogeochemical cycle. It is predicted that boron sorption onto organic matter can have important consequences for the boron isotopic composition of surface water reservoirs (seawater, groundwater, soil waters) in which either abundant organic surfaces or significant boron concentrations are available. In addition, the large isotope fractionation between aqueous boric acid and surface boron-organic complexes found in the present work makes boron a promising tracer of biologic activity.     

Chemical and isotopic constraints on water/rock interactions at the Lost City hydrothermal field, 30°N Mid-Atlantic Ridge

Low temperature vent fluids (<91 °C) issuing from the ultramafic-hosted hydrothermal system at Lost City, 30°N Mid-Atlantic Ridge, are enriched in dissolved volatiles (H₂,CH₄) while attaining elevated pH values, indicative of the serpentization processes that govern water/rock interactions deep in the oceanic crust. Here, we present a series of theoretical models to evaluate the extent of hydrothermal alteration and assess the effect of cooling on the systematics of pH-controlled B aqueous species. Peridotite-seawater equilibria calculations indicate that the mineral assemblage composed of diopside, brucite and chrysotile likely dictates fluid pH at moderate temperature serpentinization processes (<300 °C), by imposing constraints on the aCa⁺⁺/a²H⁺ ratios and the activity of dissolved SiO₂. Based on Sr abundances and the ⁸⁷Sr/⁸⁶Sr isotope ratios of vent fluids reported from Lost City, estimated water/rock mass ratios (w/r = 2-4) are consistent with published models involving dissolved CO₂ and alkane concentrations. Combining the reported δ¹⁸O values of vent fluids (0.7‰) with such w/r mass ratios, allows us to bracket subseafloor reaction temperatures in the vicinity of 250 °C. These estimates are in agreement with previous theoretical studies supporting extensive conductive heat loss within the upflow zones. Experimental studies on peridotite-seawater alteration suggest that fluid pH increases during cooling which then rapidly enhances boron removal from solution and incorporation into secondary phases, providing an explanation for the highly depleted dissolved boron concentrations measured in the low temperature but alkaline Lost City vent fluids. Finally, to account for the depleted ¹¹B composition (δ¹¹B ∼25-30‰) of vent fluids relative to seawater, isotopic fractionation between tetrahedrally coordinated aqueous boron species with BO₃-bearing mineral sites (e.g. in calcite, brucite) is proposed.     

Intra-shell boron isotope ratios in the symbiont-bearing benthic foraminiferan Amphistegina lobifera: Implications for δB vital effects and paleo-pH reconstructions

The boron isotope composition of marine carbonates is considered to be a seawater pH proxy. Nevertheless, the use of δ¹¹B has some limitations such as the knowledge of the fractionation factor (α_4-3) between boric acid and the borate ion and the amplitude of “vital effects” on this proxy that are not well constrained. Using secondary ion mass spectrometry (SIMS) we have examined the internal variability of the boron isotope ratio in the shallow water, symbionts bearing foraminiferan Amphistegina lobifera. Specimens were cultured at constant temperature (24 ± 0.1 °C) in seawater with pH ranging between 7.90 and 8.45. Intra-shell boron isotopes showed large variability with an upper limit value of ≈30‰. Our results suggest that the fractionation factor α_4-3 of 0.97352 (Klochko et al., 2006) is in better agreement with our experiments and with direct pH measurements in seawater vacuoles associated with the biomineralization process in these foraminifera. Despite the large variability of the skeletal pH values in each cultured specimen, it is possible to link the lowest calculated pH values to the experimental culture pH values while the upper pH limit is slightly below 9. This variability can be interpreted as follows: foraminifera variably increase the pH at the biomineralization site to about 9. This increase above ambient seawater pH leads to a range in δ¹¹B (Δ¹¹B) for each seawater pH. This Δ¹¹B is linearly correlated with the culture seawater pH with a slope of −13.1 per pH unit, and is independent of the fractionation factor α_4-3, or the δ¹¹B_sw through time. It may also be independent of the pK_B (the dissociation constant of boric acid) value. Therefore, Δ¹¹B in foraminifera can potentially reconstruct paleo-pH of seawater.     

First experimental determination of iron isotope fractionation between hematite and aqueous solution at hydrothermal conditions

Although iron isotopes provide a new powerful tool for tracing a variety of geochemical processes, the unambiguous interpretation of iron isotope ratios in natural systems and the development of predictive theoretical models require accurate data on equilibrium isotope fractionation between fluids and minerals. We investigated Fe isotope fractionation between hematite (Fe₂O₃) and aqueous acidic NaCl fluids via hematite dissolution and precipitation experiments at temperatures from 200 to 450 °C and pressures from saturated vapor pressure (P_sat) to 600 bar. Precipitation experiments at 200 °C and P_sat from aqueous solution, in which Fe aqueous speciation is dominated by ferric iron (Fe^III) chloride complexes, show no detectable Fe isotope fractionation between hematite and fluid, Δ⁵⁷Fe_{fluid-hematite} = δ⁵⁷Fe_fluid − δ⁵⁷Fe_hematite = 0.01 ± 0.08‰ (2 × standard error, 2SE). In contrast, experiments at 300 °C and P_sat, where ferrous iron chloride species (FeCl₂ and FeCl⁺) dominate in the fluid, yield significant fluid enrichment in the light isotope, with identical values of Δ⁵⁷Fe_{fluid-hematite} = −0.54 ± 0.15‰ (2SE) both for dissolution and precipitation runs. Hematite dissolution experiments at 450 °C and 600 bar, in which Fe speciation is also dominated by ferrous chloride species, yield Δ⁵⁷Fe_{fluid-hematite} values close to zero within errors, 0.15 ± 0.17‰ (2SE). In most experiments, chemical, redox, and isotopic equilibrium was attained, as shown by constancy over time of total dissolved Fe concentrations, aqueous Fe^II and Fe^III fractions, and Fe isotope ratios in solution, and identical Δ⁵⁷Fe values from dissolution and precipitation runs. Our measured equilibrium Δ⁵⁷Fe_{fluid-hematite} values at different temperatures, fluid compositions and iron redox state are within the range of fractionations in the system fluid-hematite estimated using reported theoretical β-factors for hematite and aqueous Fe species and the distribution of Fe aqueous complexes in solution. These theoretical predictions are however affected by large discrepancies among different studies, typically ±1‰ for the Δ⁵⁷Fe _{Fe(aq)-hematite} value at 200 °C. Our data may thus help to refine theoretical models for β-factors of aqueous iron species. This study provides the first experimental calibration of Fe isotope fractionation in the system hematite-saline aqueous fluid at elevated temperatures; it demonstrates the importance of redox control on Fe isotope fractionation at hydrothermal conditions.     

The effect of aqueous diffusion on the fractionation of chlorine and bromine stable isotopes

Diffusive isotopic fractionation factors are important in order to understand natural processes and have practical application in radioactive waste storage and carbon dioxide sequestration. We determined the isotope fractionation factors and the effective diffusion coefficients of chloride and bromide ions during aqueous diffusion in polyacrylamide gel. Diffusion was determined as functions of temperature, time and concentration. The effect of temperature is relatively large on the diffusion coefficient (D) but only small on isotope fractionation. For chlorine, the ratio, D_³⁵Cl/D_³⁷Cl varied from 1.00128 ± 0.00017 (1σ) at 2 °C to 1.00192 ± 0.00015 at 80 °C. For bromine, D_⁷⁹Br/D_⁸¹Br varied from 1.00098 ± 0.00009 at 2 °C to 1.0064 ± 0.00013 at 21 °C and 1.00078 ± 0.00018 (1σ) at 80 °C. There were no significant effects on the isotope fractionation due to concentration. The lack of sensitivity of the diffusive isotope fractionation to anything at the most common temperatures (0 to 30 °C) makes it particularly valuable for application to understanding processes in geological environments and an important natural tracer in order to understand fluid transport processes.     

Mercury isotope fractionation during liquid-vapor evaporation experiments

Liquid-vapor mercury isotope fractionation was investigated under equilibrium and dynamic conditions. Equilibrium evaporation experiments were performed in a closed glass system under atmospheric pressure between 0 and 22 °C, where vapor above the liquid was sampled at chemical equilibrium. Dynamic evaporation experiments were conducted in a closed glass system under 10⁻⁵ bar vacuum conditions varying (1) the fraction of liquid Hg evaporated at 22 °C and (2) the temperature of evaporation (22-100 °C). Both, residual liquid and condensed vapor fractions were analyzed using stannous chloride CV-MC-ICP-MS.Equilibrium evaporation showed a constant liquid-vapor fractionation factor (α_202/198) of 1.00086 ± 0.00022 (2SD, n = 6) within the 0-22 °C range. The 22 °C dynamic evaporations experiments displayed Rayleigh distillation fractionation behavior with liquid-vapor α_202/198 = 1.0067 ± 0.0011 (2SD), calculated from both residual and condensed vapor fractions. Our results confirm historical data (1920s) from Brönsted, Mulliken and coworkers on mercury isotopes separation using evaporation experiments, for which recalculated δ²⁰²Hg′ showed a liquid-vapor α_202/198 of 1.0076 ± 0.0017 (2SD). This liquid-vapor α_202/198 is significantly different from the expected kinetic α_202/198 value ((202/198)^0.5 = 1.0101). A conceptual evaporation model of back condensation fluxes within a thin layer at the liquid-vapor interface was used to explain this discrepancy. The δ²⁰²Hg′ of condensed vapor fractions in the 22-100 °C temperature range experiments showed a negative linear relationship with 10⁶/T², explained by increasing rates of exchange within the layer with the increase in temperature.Evaporation experiments also resulted in non-mass-dependent fractionation (NMF) of odd ¹⁹⁹Hg and ²⁰¹Hg isotopes, expressed as Δ¹⁹⁹Hg′ and Δ²⁰¹Hg′, the deviation in ‰ from the mass fractionation relationship with even isotopes. Liquid-vapor equilibrium yielded Δ¹⁹⁹Hg′/Δ²⁰¹Hg′ relationship of 2.0 ± 0.6 (2SE), which is statistically not different from the one predicted for the nuclear field shift effect (Δ¹⁹⁹Hg/Δ²⁰¹Hg ≈ 2.47). On the other hand, evaporation under dynamic conditions at 22 °C led to negative anomalies in the residual liquid fractions that are balanced by positive anomalies in condensed vapors with lower Δ¹⁹⁹Hg′/Δ²⁰¹Hg′ ratios of 1.2 ± 0.4 (2SD). This suggests that either magnetic isotope effects may have occurred without radical chemistry or an unknown NMF process on odd isotopes operated during liquid mercury evaporation.     

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

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

Boron speciation in aqueous fluids at 22 to 600°C and 0.1 MPa to 2 GPa

The speciation of boron in H₂O+H₃BO₃±NaCl and H₂O+Na₂B₄O₇ fluids was studied in situ at temperatures between 22 and 600°C and pressures from 0.1 MPa to ∼2 GPa using Raman spectroscopy and a hydrothermal diamond anvil cell. Additionally, we determined the frequency shifts of the 877 cm⁻¹ Raman line of [B(OH)₃]⁰ in aqueous fluids with temperature (∂ν₈₇₇/∂T)_{p = 0.1 MPa} = −0.02532 cm⁻¹K⁻¹ and pressure (∂ν₈₇₇/∂P)_{T = 22°C} = 4.06 cm⁻¹GPa⁻¹. The observed species in acidic fluids were [B(OH)₃]⁰ and smaller amounts of a four-coordinated boron species which may be attributed to dissolved metaboric acid HBO₂(aq). The ratio of this B^[4]-O species to [B(OH)₃]⁰ increases with temperature and decreases slightly with addition of NaCl. In alkaline solutions, polyboric ions depolymerize rapidly with temperature. Thus, [B(OH)₃]⁰ and [B(OH)₄]⁻ were the only remaining detectable species at 500 and 600°C. The Raman spectra showed an increase of [B(OH)₃]⁰ relative to [B(OH)₄]⁻ with temperature and an increase of [B(OH)₄]⁻ relative to [B(OH)₃]⁰ with pressure.The general trend in the boron speciation is a higher stability of simpler complexes with temperature. The experimental observations strongly indicate that planar three-coordinated [B(OH)₃]⁰ is the predominant boron species in the aqueous phase over a wide range of P-T-pH conditions. This supports the validity of previous assumptions on boron coordination in crustal and mantle wedge fluids.     

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.     

Quantifying Li isotope fractionation during smectite formation and implications for the Li cycle

Tri-octahedral Li-Mg smectites (hectorites) were synthesized at temperatures ranging from 25 to 250 °C, in the presence of solutions highly enriched in lithium. After removing all the exchangeable lithium from the synthesized clays, Li isotope fractionation (Δ⁷Li_{clay-solution}) was determined. This fractionation was linked to Li incorporation into the structural octahedral site, substituting for Mg²⁺. As predicted, experimental Δ⁷Li_{clay-solution} inversely correlates with temperature, and ranges from −1.6‰ ± 1.3‰ at 250 °C to −10.0‰ ± 1.3‰ at 90 °C, and then stays relatively constant down to 25 °C. The relatively constant isotope fractionation factor below 90 °C may be due to high concentrations of edge octahedra in low crystallinity smectites. The isotopic fractionation factor (α), for a given temperature, does not depend on the solution matrix, nor on the amount of structural Li incorporated into the clay. Empirical linear laws for α as a function of 1/T (K) were inferred. Smectite Li contents and smectite-solution distribution coefficients (D_Li/Mg) increase with temperature, as expected for a substitution process. The fractions of dissolved Li incorporated into the smectite octahedral sites are small and do not depend on the duration of the experiment. In a seawater-like matrix solution, less Li is incorporated into the smectites, probably as a result of competition with dissolved Mg²⁺ ions for incorporation into the octahedral sites. The high Li contents observed in marine smectites are therefore best explained either by a significant contribution from basalts, by adsorption processes, or by the influence of seawater chemical composition on distribution coefficients. We also calculate, using present-day estimates of hydrothermal water and river fluxes, that a steady-state ocean would require a relatively large global clay-water Li isotope fractionation (−12‰ to −21‰). This study demonstrates the ability of laboratory experiments to quantify the impact of secondary phases on the Li geochemical cycle and associated isotope fractionations.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.     

Boron isotope geochemistry of Paleozoic brachiopod calcite: Implications for a secular change in the boron isotope geochemistry of seawater over the Phanerozoic

Boron isotope composition of marine carbonates has been proposed as a paleo-pH proxy and potential tool to reconstruct atmospheric pCO₂. The precise knowledge of the boron isotopic composition of ancient seawater represents the fundamental prerequisite for any paleo-pH reconstruction. This contribution presents boron isotope values for Silurian to Permian brachiopod calcite that might be used to reconstruct pH or boron isotope composition of past oceans. All brachiopod shells were screened for diagenetic recrystallization by means of cathodoluminescence microscopy, trace element geochemistry (B, Fe, Mn, Sr) as well as SEM. Only nonluminescent shells revealing well-preserved microstructures, high strontium and boron concentrations as well as low iron and manganese contents were accepted for boron isotope analysis. The boron isotope ratios of Silurian, Devonian, Pennsylvanian and Permian brachiopod calcite range from 6.8 to 11.0‰, 7.3 to 14.9‰, 12.4 to 15.8‰ and 10.1 to 11.7‰, respectively. These δ¹¹B values are significantly lower in comparison to δ¹¹B values of modern biogenic carbonates and indicate that the Paleozoic oceans were depleted in ¹¹B by up to 10‰. Box modeling of the boron geochemical cycle suggests that the significant depletion of ¹¹B in the oceanic reservoir may have been initiated by an enhanced continental boron discharge. Our data support the earlier made conclusion that boron isotopes may not be used in the geological past as reliable paleo-pH proxy unless the boron isotopic composition of ancient oceans can be constrained by further studies.     

Hydrogen isotope fractionation factor for mixed-layer illite/smectite at 60° to 150°C: new data from the northeast Texas Gulf Coast

The δD values of water and clay collected from 0.6 to 5.5 km deep wells in 9 oil and gas fields within the Houston salt basin located in the northeast Texas Gulf Coast were used to reevaluate the temperature dependence of hydrogen isotope fractionation between mixed layer illite/smectite (I/S) and water, and the equations of Yeh (1980), Geochim. Cosmochim. Acta, 42:140-143) and Capuano (1992), Geochim. Cosmochim. Acta, 56:2547-2554) for the calculation of αI/S-water, both of which are commonly cited in the literature. δD values of water and clay in the normally pressured and geopressured sections are different. In the normally pressured section (<2.6 km), δD_water is constant (−15 ± 1 (1σ)‰) while δD_clay increases linearly from −59 to −43‰ with increasing depth. In contrast, in the geopressured section δD_water decreases linearly from −7 to −26‰ with increasing depth while δD_clay is nearly constant (−36 ± 3 (1σ)‰). The opposite trends are a product of D/H exchange between clay and water in a water-dominated system within the normally pressured section and rock-dominated system in the geopressured section. Assuming hydrogen isotope equilibrium between sediments and enclosed pore water, the new δD_water and δD_clay data were used to calculate the hydrogen isotope fractionation factor between I/S and water (αI/S-water), which showed a good fit to the equation derived by Capuano (1992), but not to the equation derived by Yeh (1980), both for the normally and geopressured data despite that both data sets are strikingly different. This suggests that hydrogen isotope equilibrium has been achieved in these two different regimes and that the equation by Capuano (1992) is valid and has wider applicability.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

Experimental study of arsenic speciation in vapor phase to 500°C: implications for As transport and fractionation in low-density crustal fluids and volcanic gases

The stoichiometry and stability of arsenic gaseous complexes were determined in the system As-H₂O ± NaCl ± HCl ± H₂S at temperatures up to 500°C and pressures up to 600 bar, from both measurements of As^(III) and As^(V) vapor-liquid and vapor-solid partitioning, and X-ray absorption fine structure (XAFS) spectroscopic study of As^(III)-bearing aqueous fluids. Vapor-aqueous solution partitioning for As^(III) was measured from 250 to 450°C at the saturated vapor pressure of the system (P_sat) with a special titanium reactor that allows in situ sampling of the vapor phase. The values of partition coefficients for arsenious acid (H₃AsO₃) between an aqueous solution (pure H₂O) and its saturated vapor (K = mAs_vapor /mAs_liquid) were found to be independent of As^(III) solution concentrations (up to ∼1 to 2 mol As/kg) and equal to 0.012 ± 0.003, 0.063 ± 0.023, and 0.145 ± 0.020 at 250, 300, and 350°C, respectively. These results are interpreted by the formation, in the vapor phase, of As(OH)₃(gas), similar to the aqueous As hydroxide complex dominant in the liquid phase. Arsenic chloride or sulfide gaseous complexes were found to be negligible in the presence of HCl or H₂S (up to ∼0.5 mol/kg of vapor). XAFS spectroscopic measurements carried out on As^(III)-H₂O (±NaCl) solutions up to 500°C demonstrate that the As(OH)₃ complex dominates As speciation both in dense H₂O-NaCl fluids and low-density supercritical vapor. Vapor-liquid partition coefficients for As^(III) measured in the H₂O-NaCl system up to 450°C are consistent with the As speciation derived from these spectroscopic measurements and can be described by a simple relationship as a function of the vapor-to-liquid density ratio and temperature. Arsenic^(III) partitioning between vapor and As-concentrated solutions (>2 mol As/kg) or As₂O₃ solid is consistent with the formation, in the vapor phase, of both As₄O₆ and As(OH)₃. Arsenic^(V) (arsenic acid, H₃AsO₄) vapor-liquid partitioning at 350°C for dilute aqueous solution was interpreted by the formation of AsO(OH)₃ in the vapor phase.The results obtained were combined with the corresponding properties for the aqueous As(III) hydroxide species to generate As(OH)₃(gas) thermodynamic parameters. Equilibrium calculations carried out by using these data indicate that As(OH)₃(gas) is by far the most dominant As complex in both volcanic gases and boiling hydrothermal systems. This species is likely to be responsible for the preferential partition of arsenic into the vapor phase as observed in fluid inclusions from high-temperature (400 to 700°C) Au-Cu (-Sn, -W) magmatic-hydrothermal ore deposits. The results of this study imply that hydrolysis and hydration could be also important for other metals and metalloids in the H₂O-vapor phase. These processes should be taken into account to accurately model element fractionation and chemical equilibria during magma degassing and fluid boiling.     

A critical evaluation of the boron isotope-pH proxy: The accuracy of ancient ocean pH estimates

The boron isotope-pH technique is founded on a theoretical model of carbonate δ¹¹B variation with pH that assumes that the boron isotopic composition of carbonates mirrors the boron isotopic composition of borate in solution (δ¹¹B_carb = δ¹¹B_borate). Knowledge of the fractionation factor for isotope exchange between boric acid and borate in solution (α_4-3), the equilibrium constant for the dissociation of boric acid (pK_B*), as well as the isotopic composition of boron in seawater (δ¹¹B_sw) are required parameters of the model.The available data suggests that both the value of α_4-3 and the history of δ¹¹B_sw are poorly constrained. However, if one assumes that δ¹¹B_carb = δ¹¹B_borate, an empirical value for α_4-3 can be estimated from the results of inorganic carbonate precipitation experiments. This exercise yields an α_4-3 value of ∼0.974 in accordance with recent theoretical estimates, but substantially deviates from the theoretical value of 0.981 often used to estimate paleo-ocean pH. Re-evaluation of ocean pH using an α_4-3 value of 0.974 and published foraminiferal δ¹¹B values for the Cenozoic yield pH estimates that are relatively invariant, but unrealistically high (∼8.4-8.6). Uncertainty increases as foraminiferal ‘vital effects’ are considered and different models for secular changes in seawater δ¹¹B are applied.The inability to capture realistic ocean pH possibly reflects on our understanding of the isotopic relationship between carbonate and borate, as well as the mechanism of boron incorporation in carbonates. Given the current understanding of boron systematics, pH values estimated using this technique have considerable uncertainty, particularly when reconstructions exceed the residence time of boron in the ocean.     

Isotopic and element exchange during serpentinization and metasomatism at the Atlantis Massif (MAR 30°N): Insights from B and Sr isotope data

The Lost City hydrothermal system at the southern Atlantis Massif (Mid-Atlantic Ridge, 30°N) provides a natural laboratory for studying serpentinization processes, the temporal evolution of ultramafic-hosted hydrothermal systems, and alteration conditions during formation and emplacement of an oceanic core complex. Here we present B, O, and Sr isotope data to investigate fluid/rock interaction and mass transfer during detachment faulting and exhumation of lithospheric sequences within the Atlantis Massif. Our data indicate that extensive serpentinization was a seawater-dominated process that occurred predominately at temperatures of 150-250 °C and at high integrated W/R ratios that led to a marked boron enrichment (34-91 ppm). Boron removal from seawater during serpentinization is positively correlated with changes in δ¹¹B (11-16‰) but shows no correlation with O-isotope composition. Modeling indicates that B concentrations and isotope values of the serpentinites are controlled by transient temperature-pH conditions. In contrast to prior studies, we conclude that low-temperature marine weathering processes are insignificant for boron geochemistry of the Atlantis Massif serpentinites. Talc- and amphibole-rich fault rocks formed within a zone of detachment faulting at temperatures of approximately 270-350 °C and at low W/R ratios. Talc formation in ultramafic domains in the massif was subsequent to an early stage of serpentinization and was controlled by the access of Si-rich fluids derived through seawater-gabbro interactions. Replacement of serpentine by talc resulted in boron loss and significant lowering of δ¹¹B values (9-10‰), which we model as the product of progressive extraction of boron. Our study provides new constraints on the boron geochemical cycle at oceanic spreading ridges and suggests that serpentinization associated with ultramafic-hosted hydrothermal systems may have important implications for the behavior of boron in subduction zone settings.