Argon and CO2 on the race track in silicate melts: A tool for the development of a CO2 speciation and diffusion model

We have analysed the kinetics of Argon and CO₂ diffusion in simplified iron free rhyolitic to hawaiitic melts using the diffusion couple technique. The concentration distance profiles of Ar and CO₂ were measured with electron microprobe analysis and Fourier Transform Infrared Spectroscopy, respectively. Error functions were fitted to the symmetrical concentration distance profiles to extract the diffusion coefficients.In the temperature range 1373 to 1773 K the activation energies for Ar diffusion range from 169 ± 20 to 257 ± 62 kJ mol⁻¹. Ar diffusivity increases exponentially with the degree of depolymerisation. In contrast, the mobility of total CO₂, that is identical to Ar mobility in rhyolitic melt, keeps constant with changing bulk composition from rhyolite to hawaiite. CO₂ speciation at 1623 K and 500 MPa was modeled for the range of compositions studied using the diffusion data of Ar and total CO₂ in combination with network former diffusion calculated from viscosity data. Within error this model is in excellent agreement with CO₂ speciation data extrapolated from temperatures near the glass transition temperature for dacitic melt composition. This model shows that even in highly depolymerised hawaiitic and tholeiitic melts molecular CO₂ is a stable species and contributes 70 to 80% to the total CO₂ diffusion, respectively.     

Surface simulation studies of the hydration of white rust Fe(OH)2, goethite α-FeO(OH) and hematite α-Fe2O3

Computer modelling techniques were used to elucidate the hydration behaviour of three iron (hydr)oxide minerals at the atomic level: white rust, goethite and hematite. A potential model was first adapted and tested against the bulk structures and properties of eight different iron oxides, oxyhydroxides and hydroxides, followed by surface simulations of Fe(OH)₂, α-FeO(OH) and α-Fe₂O₃. The major interaction between the adsorbing water molecules and the surface is through interaction of their oxygen ions with surface iron ions, followed by hydrogen-bonding to surface oxygen ions. The energies released upon the associative adsorption of water range from 1 to 17 kJ mol⁻¹ for Fe(OH)₂, 26 to 80 kJ mol⁻¹ for goethite and 40 to 85 kJ mol⁻¹ for hematite, reflecting the increasing oxidation of the iron mineral. Dissociative adsorption at goethite and hematite surfaces releases larger hydration energies, ranging from 120 to 208 kJ mol⁻¹ for goethite and 76 to 190 kJ mol⁻¹ for hematite.The thermodynamic morphologies of the minerals, based on the calculated surface energies, agree well with experimental morphologies, where these are available. When the partial pressures required for adsorption of water from the gas phase are plotted against temperature for the goethite and hematite surfaces, taking into account experimental entropies for water, it appears that these minerals may well be instrumental in the retention of water during the cyclic variations in the atmosphere of Mars.     

Transport and thermodynamics constrain belowground carbon turnover in a northern peatland

Rates of anaerobic respiration are of central importance for the long-term burial of carbon (C) in peatlands, which are a relevant sink in the global C cycle. To identify constraints on anaerobic peat decomposition, we determined detailed concentration depth profiles of decomposition end-products, i.e. methane (CH₄) and dissolved inorganic carbon (DIC), along with concentrations of relevant decomposition intermediates at an ombrotrophic Canadian peat bog. The magnitude of in situ net production rates of DIC and CH₄ was estimated by inverse pore-water modeling. Vertical transport in the peat was slow and dominated by diffusion leading to the buildup of DIC and CH₄ with depth (5500 μmol L⁻¹ DIC, 500 μmol L⁻¹ CH₄). Highest DIC and CH₄ production rates occurred close to the water table (decomposition constant k_d ∼ 10⁻³-10⁻⁴ a⁻¹) or in some distinct zones at depth (k_d ∼ 10⁻⁴ a⁻¹). Deeper into the peat, decomposition proceeded very slowly at about k_d = 10⁻⁷ a⁻¹. This pattern could be related to thermodynamic and transport constraints. The accumulation of metabolic end-products diminished in situ energy yields of acetoclastic methanogenesis to the threshold for microbially mediated processes (−20 to −25 kJ mol⁻¹ CH₄). The methanogenic precursor acetate also accumulated (150 μmol L⁻¹). In line with these findings, CH₄ was formed by hydrogenotrophic methanogenesis at Gibbs free energies of −35 to −40 kJ mol⁻¹ CH₄. This was indicated by an isotopic fractionation α_CO2-CH4 of 1.069-1.079. Fermentative degradation of acetate, propionate and butyrate attained Gibbs free energies close to 0 kJ mol⁻¹ substrate. Although methanogenesis was apparently limited by some other factor in some peat layers, transport and thermodynamic constraints likely impeded respiratory processes in the deeper peat. Constraints on the removal of DIC and CH₄ may thus slow decomposition and contribute to the sustained burial of C in northern peatlands.     

Water speciation and diffusion in haploandesitic melts at 743-873 K and 100 MPa

Water is an important volatile component in andesitic eruptions and deep-seated andesitic magma chambers. We report an investigation of H₂O speciation and diffusion by dehydrating haploandesitic melts containing ?2.5 wt.% water at 743-873 K and 100 MPa in cold-seal pressure vessels. FTIR microspectroscopy was utilized to measure species [molecular H₂O (H₂O_m) and hydroxyl group (OH)] and total H₂O (H₂O_t) concentration profiles on the quenched glasses from the dehydration experiments. The equilibrium constant of the H₂O speciation reaction H₂O_m+O?2OH, K = (X_OH)²/(X_H2OmX_O) where X means mole fraction on a single oxygen basis, in this Fe-free andesite varies with temperature as ln K = 1.547-2453/T where T is in K. Comparison with previous speciation data on rhyolitic and dacitic melts indicates that, for a given water concentration, Fe-free andesitic melt contains more hydroxyl groups. Water diffusivity at the experimental conditions increases rapidly with H₂O concentration, contrary to previous H₂O diffusion data in an andesitic melt at 1608-1848 K. The diffusion profiles are consistent with the model that molecular H₂O is the diffusion species. Based on the above speciation model, H₂O_m and H₂O_t diffusivity (in m²/s) in haploandesite at 743-873 K, 100 MPa, and H₂O_t ? 2.5 wt.% can be formulated as

     

CO2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.     

Water solubility in haplogranitic melts coexisting with H2O-H2 fluids

The water solubility in haplogranitic melts (normative composition Ab₃₉Or₃₂Qz₂₉) coexisting with H₂O-H₂ fluids at 800 and 950 °C and 1, 2 and 3 kbar vapour pressure has been determined using IR spectroscopy. The experiments were performed in internally heated pressure vessels and the hydrogen fugacity (f _H2) was controlled using the double capsule technique and oxygen buffer assemblages (WM and IW). Due to the limited lifetimes of these oxygen buffers the water solubility was determined from diffusion profiles (concentration-distance profiles) measured with IR spectroscopy in the quenched glasses. The reliability of the experimental strategy was demonstrated by comparing the results of short- and long-duration experiments performed with pure H₂O fluids. The water solubility in Ab₃₉Or₃₂Qz₂₉ melts equilibrated with H₂O-H₂ fluids decreases progressively with decreasing f _H2O, as f _H2 (or X _H2) increases in the fluid phase. The effect of H₂ on the evolution of the water solubility is similar to that of CO₂ or another volatile with a low solubility in the melt and can be calculated in a first approximation with the Burnham water solubility model. Recalculation of high temperature water speciation for AOQ melts coexisting with H₂O-H₂ fluids at 800 °C, 2 kbar suggests that the concentrations of molecular H₂O are proportional to f _H2O (calculated using available mixing models), indicating Henrian behaviour for the solubility of molecular H₂O in haplogranitic melts. Received: 29 June 1998 / Accepted: 10 March 1999     

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.     

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.     

Carbon dioxide in silicate melts: A molecular dynamics simulation study

We have performed a series of molecular dynamics simulations aimed at the evaluation of the solubility of CO₂ in silicate melts of natural composition (from felsic to ultramafic). In making in contact within the simulation cell a supercritical CO₂ phase with a silicate melt of a given composition, we have been able to evaluate (i) the solubility of CO₂ in the P-T range 1473-2273 K and 20-150 kbar, (ii) the density change experienced by the CO₂-bearing melt, (iii) the respective concentrations of CO₂ and species in the melt, (iv) the lifetime and the diffusivity of these species and (v) the structure of the melt around the carbonate groups. The main results are the following:(1) The solubility of CO₂ increases markedly with the pressure in the three investigated melts (a rhyolite, a mid-ocean ridge basalt and a kimberlite) from about ∼2 wt% CO₂ at 20 kbar to ∼25 wt% at 100 kbar and 2273 K. The solubility is found to be weakly dependent on the melt composition (as far as the present compositions are concerned) and it is only at very high pressure (above ∼100 kbar) that a clear hierarchy between solubilities occurs (rhyolite < MORB < kimberlite). Furthermore at a given pressure the calculated solubility is negatively correlated with the temperature.(2) In CO₂-saturated melts, the proportion of carbonate ions is positively correlated with the pressure at isothermal condition and is negatively correlated with the temperature at isobaric condition (and vice versa for molecular CO₂). Furthermore, at fixed (P, T) conditions the proportion of carbonate ions is higher in CO₂-undersaturated melts than in the CO₂-saturated melt. Although the proportion of molecular CO₂ decreases when the degree of depolymerization of the melt increases, it is still significant in CO₂-saturated basic and ultrabasic compositions at high temperatures. This finding is at variance with experimental data on CO₂-bearing glasses which show no evidence of molecular CO₂ as soon as the degree of depolymerization of the melt is high (e.g. basalt). These conflicting results can be reconciled with each other by noticing that a simple low temperature extrapolation of the simulation data predicts that the proportion of molecular CO₂ in basaltic melts might be negligible in the glass at room temperature.(3) The carbonate ions are found to be transient species in the liquid phase, with a lifetime increasing exponentially with the inverse of the temperature. Contrarily to a usual assumption, the diffusivity of carbonate ions into the liquid silicate is not vanishingly small with respect to that of CO₂ molecules: in MORB they differ from each other by a factor of ∼6 at 1473 K and only a factor of ∼2 at 2273 K. Although the bulk diffusivity of CO₂ is governed primarily by the diffusivity of CO₂ molecules, the carbonate ions contribute significantly to the diffusivity of CO₂ in depolymerized melts.(4) Concerning the structure of the CO₂-bearing silicate melt, the carbonate ions are found to be preferentially associated with NBO’s of the melt, with an affinity for NBOs which exceeds that for BOs by almost one order of magnitude. This result explains why the concentration in carbonate ions is positively correlated with the degree of depolymerization of the melt and diminishes drastically in fully polymerized melts where the number of NBO’s is close to zero. Furthermore, the network modifier cations are not randomly distributed in the close vicinity of carbonate groups but exhibit a preferential ordering which depends at once on the nature of the cation and on the melt composition. However at the high temperatures investigated here, there is no evidence of long lived complexes between carbonate groups and metal cations.     

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 .     

Dissolution of Columbia River Basalt under mildly acidic conditions as a function of temperature: Experimental results relevant to the geological sequestration of carbon dioxide

Increasing attention is being focused on the rapid rise of CO₂ levels in the atmosphere, which many believe to be the major contributing factor to global climate change. Sequestering CO₂ in deep geological formations has been proposed as a long-term solution to help stabilize CO₂ levels. However, before such technology can be developed and implemented, a basic understanding of H₂O–CO₂ systems and the chemical interactions of these fluids with the host formation must be obtained. Important issues concerning mineral stability, reaction rates, and carbonate formation are all controlled or at least significantly impacted by the kinetics of rock–water reactions in mildly acidic, CO₂-saturated solutions. Basalt has recently been identified as a potentially important host formation for geological sequestration. Dissolution kinetics of the Columbia River Basalt (CRB) were measured for a range of temperatures (25–90 °C) under mildly acidic to neutral pH conditions using the single-pass flow-through test method. Under anaerobic conditions, the normalized dissolution rates for CRB decrease with increasing pH (3 ? pH ? 7) with a slope, η, of −0.15 ± 0.01. Activation energy, E_a, has been estimated at 32.0 ± 2.4 kJ mol⁻¹. Dissolution kinetics measurements like these are essential for modeling the rate at which CO₂-saturated fluids react with basalt and ultimately drive conversion rates to carbonate minerals in situ.     

Cu and Fe chalcopyrite leach activation energies and the effect of added Fe

The leaching kinetics of chalcopyrite (CuFeS₂) concentrate in sulfuric acid leach media with and without the initial addition of Fe³⁺ under carefully controlled solution conditions (E_h 750 mV SHE, pH 1) at various temperatures from 55 to 85 °C were measured. Kinetic analyses by (i) apparent rate (not surface area normalised), and rate dependence using (ii) a shrinking core model and (iii) a shrinking core model in conjunction with Fe³⁺ activity, were performed to estimate the activation energies (E_a) for Cu and Fe dissolution.The E_a values determined for Cu and Fe leaching in the absence of added Fe³⁺ are within experimental error, 80 ± 10 kJ mol⁻¹ and 84 ± 10 kJ mol⁻¹, respectively (type iii analyses E_a are quoted unless stated otherwise), and are indicative of a chemical reaction controlled process. On addition of Fe³⁺ the initial Cu leach rate (up to 10 h) was increased and Cu was released to solution preferentially over Fe, with the E_a value of 21 ± 5 kJ mol⁻¹ (type ii analysis) suggestive of a transport controlled rate determining process. However, the rate of leaching rapidly decreased until it was consistently slower than for the equivalent leaches where Fe³⁺ was not added. The resulting E_a value for this leach regime of 83 ± 10 kJ mol⁻¹ is within experimental error of that determined in the absence of added Fe³⁺. In contrast to Cu release, Fe release to solution was consistent with a chemical reaction controlled leach rate throughout. The Fe release E_a of 76 ± 10 kJ mol⁻¹ is also within experimental error of that determined in the absence of added Fe³⁺. Where type (ii) and (iii) analyses were both successfully carried out (in all cases except for Cu leaching with added Fe³⁺, <10 h) the E_a derived are within experimental error. However, the type (iii) analyses of the leaches in the presence of added Fe³⁺ (>10 h), as compared to in the absence of added Fe³⁺, returned a considerably smaller pre-exponential factors for both Cu and Fe leach analyses commensurate with the considerably slower leach rate, suggestive of a more applicable kinetic analysis.XPS examination of leached chalcopyrite showed that the surface concentration of polysulfide and sulfate was significantly increased when Fe³⁺ was added to the leach liquor. Complementary SEM analysis revealed the surface features of chalcopyrite, most likely due to the nature of the polysulfide formed, are subtly different with greater surface roughness upon leaching in the absence of added Fe³⁺ as compared to a continuous smooth surface layer formed in the presence of added Fe³⁺. These observations suggest that the effect of Fe³⁺ addition on the rate of leaching is not due to the change in the chemical reaction controlled mechanism but due to a change in the available surface area for reaction.     

Activation energy for diffusion of helium in water-saturated, compacted Na-montmorillonite

It is important to clarify the migration behavior of hydrogen gas dissolved in water-saturated, compacted bentonite, which is a promising material for geologic disposal of high-level waste and TRU waste disposal. The diffusion coefficients of helium, which can be detected under extremely low background conditions, in water-saturated, compacted Na-montmorillonite were determined as a function of temperature by a transient diffusion method. The activation energies for diffusion of helium were then obtained. The activation energies were from 6.9 ± 4.8 to 19 ± 2.8 kJ mol^− 1 and were regarded to be independent of dry density. The activation energies of helium in water-saturated Na-montmorillonite were roughly equal to those in bulk water, 14.9 kJ mol^− 1, and in ice, from 11 to 13 kJ mol^− 1. It is possible that helium diffuses not only in pore water but also in interlayer water.     

Oxygen isotope fractionation of dissolved oxygen during reduction by ferrous iron

The oxygen isotope fractionation factor of dissolved oxygen gas has been measured during inorganic reduction by aqueous FeSO₄ at 10−54 °C under neutral (pH 7) and acidic (pH 2) conditions, with Fe(II) concentrations ranging up to 0.67 mol L⁻¹, in order to better understand the geochemical behavior of oxygen in ferrous iron-rich groundwater and acidic mine pit lakes. The rate of oxygen reduction increased with increasing temperature and increasing Fe(II) concentration, with the pseudo-first-order rate constant k ranging from 2.3 to 82.9 × 10⁻⁶ s⁻¹ under neutral conditions and 2.1 to 37.4 × 10⁻⁷ s⁻¹ under acidic conditions. The activation energy of oxygen reduction was 30.9 ± 6.6 kJ mol⁻¹ and 49.7 ± 13.0 kJ mol⁻¹ under neutral and acidic conditions, respectively. Oxygen isotope enrichment factors (ε) become smaller with increasing temperature, increasing ferrous iron concentration, and increasing reaction rate under acidic conditions, with ε values ranging from −4.5‰ to −11.6‰. Under neutral conditions, ε does not show any systematic trends vs. temperature or ferrous iron concentration, with ε values ranging from −7.3 to −10.3‰. Characterization of the oxygen isotope fractionation factor associated with O₂ reduction by Fe(II) will have application to elucidating the process or processes responsible for oxygen consumption in environments such as groundwater and acidic mine pit lakes, where a number of possible processes (e.g. biological respiration, reduction by reduced species) may have taken place.     

Water diffusion in dacitic melt

H₂O diffusion in dacitic melt was investigated at 0.48-0.95 GPa and 786-893 K in a piston-cylinder apparatus. The diffusion couple design was used, in which a nominally dry dacitic glass makes one half and is juxtaposed with a hydrous dacitic glass containing up to ∼8 wt.% total water (H₂O_t). H₂O concentration profiles were measured on quenched glasses with infrared microspectroscopy. The H₂O diffusivity in dacite increases rapidly with water content under experimental conditions, similar to previous measurements at the same temperature but at pressure <0.15 GPa. However, compared with the low-pressure data, H₂O diffusion at high pressure is systematically slower. H₂O diffusion profiles in dacite can be modeled by assuming molecular H₂O (H₂O_m) is the diffusing species. Total H₂O diffusivity D_H2Ot within 786-1798 K, 0-1 GPa, and 0-8 wt.% H₂O_t can be expressed as: where D_H2Ot is in m²/s, T is temperature in K, P is pressure in GPa, K = exp(1.49 − 2634/T) is the equilibrium constant of speciation reaction (H₂O_m+O?2OH) in the melt, X = C/18.015/[C/18.015 + (100 − C)/33.82], C is wt.% of H₂O_t, and 18.015 and 33.82 g/mol correspond to the molar masses of H₂O and anhydrous dacite on a single oxygen basis. Compared to H₂O diffusion in rhyolite, diffusivity in dacite is lower at intermediate temperatures but higher at superliquidus temperatures. This general H₂O diffusivity expression can be applied to a broad range of geological conditions, including both magma chamber processes and volcanic eruption dynamics from conduit to the surface.     

The effect of fluorine on liquidus phase relationships in the system Qz-Ab-Or with excess water at 1 kb

Liquidus phase relationships have been determined experimentally for the system Qz-Ab-Or with excess water and 1, 2, and 4 wt.% added fluorine at 1 kb pressure. With increasing fluorine content the position of the quartz-alkali feldspar field boundary moves away from the quartz apex. The position of the minimum melting composition and the minimum liquidus temperature change progressively from Qz₃₇Ab₃₄Or₂₉ and 730° C for the fluorine free system (Tuttle and Bowen 1958) to Qz₁₅Ab₅₈Or₂₇ and 630° C for the system with 4 wt.% added fluorine. Exploratory experiments have been carried out below the liquidus, and have indicated that for certain bulk compositions an assemblage consisting of two alkali feldspars, quartz, melt and vapour can exist at temperatures as low as 550° C at 1 kb.The experimental results suggest that there may be an interaction between fluorine and aluminosilicate complexes present within the melt, to produce aluminofluoride (AlF ₆ ^3– ) complex anions (Manning et al. 1980). The observed changes in liquidus phase relationships with increasing fluorine content indicate that the compositions of certain fluorine-rich granitic rocks are consistent with an origin by crystallisation of residual melts enriched in fluorine by magmatic differentiation. Such residual melts may exist at relatively low temperatures, and may form part of a continuum between granite magmatism and associated hydrothermal activity. Because of the observed preference of fluorine for aluminosilicate phases at the magmatic stage, the presence of fluorine alone is not considered to play a direct part in the generation of residual mineralising hydrothermal fluids.     

The mechanisms of water diffusion in polymerized silicate melts

Diffusion of water was experimentally investigated for melts of albitic (Ab) and quartz-orthoclasic (Qz₂₉Or₇₁, in wt %) compositions with water contents in the range of 0 to 8.5 wt % at temperatures of 1100 to 1200 °C and at pressures of 1.0 and 5.0 kbar. Apparent chemical diffusion coefficients of water (D _water) were determined from concentration-distance profiles measured by FTIR microspectroscopy. Under the same P-T condition and water content the diffusivity of water in albitic, quartz-orthoclasic and haplogranitic (Qz₂₈Ab₃₈ Or₃₄, Nowak and Behrens, this issue) melts is identical within experimental error. Comparison to data published in literature indicates that anhydrous composition only has little influence on the mobility of water in polymerized melts but that the degree of polymerization has a large effect. For instance, D_water is almost identical for haplogranitic and rhyolitic melts with 0.5–3.5 wt % water at 850 °C but it is two orders of magnitude higher in basaltic than in haplogranitic melts with 0.2–0.5 wt % water at 1300 °C. Based on the new water diffusivity data, recently published in situ near-infrared spectroscopic data (Nowak 1995; Nowak and Behrens 1995), and viscosity data (Schulze et al. 1996) for hydrous haplogranitic melts current models for water diffusion in silicate melts are critically reviewed. The NIR spectroscopy has indicated isolated OH groups, pairs of OH groups and H₂O molecules as hydrous species in polymerized silicate melts. A significant contribution of isolated OH groups to the transport of water is excluded for water contents above 10 ppm by comparison of viscosity and water diffusion data and by inspection of concentration profiles from trace water diffusion. Spectroscopic measurements have indicated that the interconversion of H₂O molecules and OH pairs is relatively fast in silicate glasses and melts even at low temperature and it is inferred that this reaction is an active step for migration of water. However, direct jumps of H₂O molecules from one cavity within the silicate network to another one can not be excluded. Thus, we favour a model in which water migrates by the interconversion reaction and, possibly, small sequences of direct jumps of H₂O molecules. In this model, immobilization of water results from dissociation of the OH pairs. Assuming that the frequency of the interconversion reaction is faster than that of diffusive jumps, OH pairs and water molecules can be treated as a single diffusing species having an effective diffusion coefficient . The shape of curves of D_water versus water content implies that increases with water content. The change from linear to exponential dependence of D_water between 2 and 3 wt % water is attributed to the influence of the dissociation reaction at low water content and to the modification of the melt structure by incorporation of OH groups. Received: 26 March 1996 / Accepted: 23 August 1996     

Thermodynamic properties of H4SiO4 in the ideal gas state as evaluated from experimental data

Solid phases of silicon dioxide react with water vapor with the formation of hydroxides and oxyhydroxides of silica. Recent transpiration and mass-spectrometric studies convincingly demonstrate that H₄SiO₄ is the predominant form of silica in vapor phase at water pressure in excess of 10⁻² MPa. Available literature transpiration and solubility data for the reactions of solid SiO₂ phases and low-density water, extending from 424 to 1661 K, are employed for the determination of Δ_fG⁰, Δ_fH⁰ and S⁰ of H₄SiO₄ in the ideal gas state at 298.15 K, 0.1 MPa. In total, there are 102 data points from seven literature sources. The resulting values of the thermodynamic functions of H₄SiO₄(g) are: Δ_fG⁰ = −1238.51 ± 3.0 kJ mol⁻¹, Δ_fH⁰ = −1340.68 ± 3.5 kJ mol⁻¹ and S⁰ = 347.78 ± 6.2 J K⁻¹ mol⁻¹. These values agree quantitatively with one set of ab initio calculations. The relatively large uncertainties are mainly due to conflicting data for H₄SiO₄(g) from various sources, and new determinations of would be helpful. The thermodynamic properties of this species, H₄SiO₄(g), are necessary for realistic modeling of silica transport in a low-density water phase. Applications of this analysis may include the processes of silicates condensation in the primordial solar nebula, the precipitation of silica in steam-rich geothermal systems and the corrosion of SiO₂-containing alloys and ceramics in moist environments.     

Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar

Here we report on an experimental investigation of the relation between the dissolution rate of albite feldspar and the Gibbs free energy of reaction, ΔG_r. The experiments were carried out in a continuously stirred flow-through reactor at 150 °C and pH_(150 °C) 9.2. The dissolution rates R are based on steady-state Si and Al concentrations and sample mass loss. The overall relation between ΔG_r and R was determined over a free energy range of −150 < ΔG_r < −15.6 kJ mol⁻¹. The data define a continuous and highly non-linear, sigmoidal relation between R and ΔG_r that is characterized by three distinct free energy regions. The region furthest from equilibrium, delimited by −150 < ΔG_r < −70 kJ mol⁻¹, represents an extensive dissolution rate plateau with an average rate . In this free energy range the rates of dissolution are constant and independent of ΔG_r, as well as [Si] and [Al]. The free energy range delimited by −70 ? ΔG_r ? −25 kJ mol⁻¹, referred to as the ‘transition equilibrium’ region, is characterized by a sharp decrease in dissolution rates with increasing ΔG_r, indicating a very strong inverse dependence of the rates on free energy. Dissolution nearest equilibrium, defined by ΔG_r > −25 kJ mol⁻¹, represents the ‘near equilibrium’ region where the rates decrease as chemical equilibrium is approached, but with a much weaker dependence on ΔG_r. The lowest rate measured in this study, R = 6.2 × 10⁻¹¹ mol m⁻² s⁻¹ at ΔG_r = −16.3 kJ mol⁻¹, is more than two orders of magnitude slower than the plateau rate. The data have been fitted to a rate equation (adapted from Burch et al. [Burch, T. E., Nagy, K. L., Lasaga, A. C., 1993. Free energy dependence of albite dissolution kinetics at 80 °C and pH 8.8. Chem. Geol.105, 137-162]) that represents the sum of two parallel reactions

R=k₁[1-exp(-ng^m₁)]+k₂[1-exp(-g)]^m₂,     

Solubility and near-equilibrium dissolution rates of quartz in dilute NaCl solutions at 398-473 K under alkaline conditions

The dissolution-precipitation of quartz controls porosity and permeability in many lithologies and may be the best studied mineral-water reaction. However, the rate of quartz-water reaction is relatively well characterized far from equilibrium but relatively unexplored near equilibrium. We present kinetic data for quartz as equilibrium is approached from undersaturation and more limited data on the approach from supersaturated conditions in 0.1 molal NaCl + NaOH + NaSiO(OH)₃ solutions with pH 8.2-9.7 at 398, 423, 448, and 473 K. We employed a potentiometric technique that allows precise determination of solution speciation within 2 kJ mol⁻¹ of equilibrium without the need for to perturb the system through physical sampling and chemical analysis. Slightly higher equilibrium solubilities between 423 and 473 K were found than reported in recent compilations. Apparent activation energies of 29 and 37 kJ mol⁻¹ are inferred for rates of dissolution at two surface sites with different values of connectedness: dissolution at Q¹ or Q² silicon sites, respectively. The dissolution mechanism varies with ΔG such that reactions at both sites control dissolution up until a critical free energy value above which only reactions at Q¹ sites are important. When our near-equilibrium dissolution rates are extrapolated far from equilibrium, they agree within propagated uncertainty at 398 K with a recently published model by Bickmore et al. (2008). However, our extrapolated rates become progressively slower than model predictions with increasing temperature. Furthermore, we see no dependence of the postulated Q¹ reaction rate on pH, and a poorly-constrained pH dependence of the postulated Q² rate. Our slow extrapolated rates are presumably related to the increasing contribution of dissolution at Q³ sites far from equilibrium. The use of the potentiometric technique for rate measurement will yield both rate data and insights into the mechanisms of dissolution over a range of chemical affinity. Such measurements are needed to model the evolution of many natural systems quantitatively.