Experimental determination of corundum solubilities in pure water between 400–700°C and 1–3 kbar

Experimentally reversed corundum solubilities in pure water at 400° to 700°C and 0.7 to 3 kbar yield values of dissolved aluminum that range from 1–4 ppm Al. At constant pressure the solubility shows a sigmoidal behavior with a slight maximum at 500°C and minimum at 600°C. Corundum solubility increases with increasing pressure at constant temperature. The dissolved aluminum appears to form an uncharged, but polar species under these conditions probably of the form Al(OH)₃⁰.     

Phase relations in the CH4-H2O-NaCl system at 2 kbar, 300 to 600°C as determined using synthetic fluid inclusions

Synthesis of fluid inclusions in the CH₄-H₂O-NaCl system was accomplished by subjecting fractured quartz or fluorite, along with known quantities of CH₄, H₂O, and NaCl, to a pressure of 2 kbar and temperatures of 300, 400, 500, or 600°C, in sealed Au capsules. Under the elevated P-T conditions, some of the fractures healed, trapping fluids as inclusions. Microthermometric measurements conducted on the fluid inclusions show that at 2 kbar and 400 to 600°C, there are very broad regions of fluid unmixing in the CH₄-H₂O-NaCl system. For those bulk fluid compositions that lie in the two-phase (i.e., immiscible fluids) field, the high-density phase is enriched in NaCl, whereas the low-density phase is enriched in CH₄. For any given bulk composition, the degree of NaCl enrichment in the high-density phase increases, whereas the degree of CH₄ enrichment in the low-density phase decreases, as temperature increases from 400 to 600°C. Our experimental constraints on the size of the two-phase field are generally consistent with results generated using the equation-of-state GEOFLUIDS (available at http://geotherm.ucsd.edu/geofluids/). However, when comparing the compositions of coexisting immiscible fluids, as determined experimentally vs. calculated using GEOFLUIDS, we find that some relatively small but probably significant differences exist between our experiments and this equation of state.     

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.     

Solubilities of corundum, wollastonite and quartz in H2O-NaCl solutions at 800 °C and 10 kbar: Interaction of simple minerals with brines at high pressure and temperature

Solubilities of corundum (Al₂O₃) and wollastonite (CaSiO₃) were measured in H₂O-NaCl solutions at 800 °C and 10 kbar and NaCl concentrations up to halite saturation by weight-loss methods. Additional data on quartz solubility at a single NaCl concentration were obtained as a supplement to previous work. Single crystals of synthetic corundum, natural wollastonite or natural quartz were equilibrated with H₂O and NaCl at pressure (P) and temperature (T) in a piston-cylinder apparatus with NaCl pressure medium and graphite heater sleeves. The three minerals show fundamentally different dissolution behavior. Corundum solubility undergoes large enhancement with NaCl concentration, rising rapidly from Al₂O₃ molality (m_Al2O3) of 0.0013(1) (1σ error) in pure H₂O and then leveling off to a maximum of ∼0.015 at halite saturation (X_NaCl ≈ 0.58, where X is mole fraction). Solubility enhancement relative to that in pure H₂O, , passes through a maximum at X_NaCl ≈ 0.15 and then declines towards halite saturation. Quenched fluids have neutral pH at 25 °C. Wollastonite has low solubility in pure H₂O at this P and T(m_CaSiO3=0.0167(6)). It undergoes great enhancement, with a maximum solubility relative to that in H₂O at X_NaCl ≈ 0.33, and solubility >0.5 molal at halite saturation. Solute silica is 2.5 times higher than at quartz saturation in the system H₂O-NaCl-SiO₂, and quenched fluids are very basic (pH 11). Quartz shows monotonically decreasing solubility from m_SiO2=1.248 in pure H₂O to 0.202 at halite saturation. Quenched fluids are pH neutral. A simple ideal-mixing model for quartz-saturated solutions that requires as input only the solubility and speciation of silica in pure H₂O reproduces the data and indicates that hydrogen bonding of molecular H₂O to dissolved silica species is thermodynamically negligible. The maxima in for corundum and wollastonite indicate that the solute products include hydrates and Na⁺ and/or Cl⁻ species produced by molar ratios of reactant H₂O to NaCl of 6:1 and 2:1, respectively. Our results imply that quite simple mechanisms may exist in the dissolution of common rock-forming minerals in saline fluids at high P and T and allow assessment of the interaction of simple, congruently soluble rock-forming minerals with brines associated with deep-crustal metamorphism.     

Liquid-vapor fractionation of boron and boron isotopes: Experimental calibration at 400°C/23 MPa to 450°C/42 MPa

We experimentally determined the boron partitioning and boron isotope fractionation between coexisting liquid and vapor in the system H₂O−NaCl−B₂O₃. Experiments were performed along the 400 and 450°C isotherms. Pressure conditions ranged from 23 to 28 MPa at 400°C and from 38 to 42 MPa at 450°C. Boron partitions preferentially into the liquid. Its overall liquid-vapor fractionation is, however, weak: Calculated boron distribution coefficients D_B^liquid-vapor are < 2.5 at all run conditions. With decreasing pressure (i.e. increasing opening of the solvus) D_B^liquid-vapor increases along the individual isotherms. Extrapolation to salt saturated conditions yields maximum boron liquid-vapor fractionations of D_B^liquid-vapor = 1.8 at 450°C and D_B^liquid-vapor = 2.7 at 400°C. ¹¹B preferentially fractionates into the vapor. Calculated Δ¹¹B_vapor-liquid = {[(¹¹B/¹⁰B)_vapor - (¹¹B/¹⁰B)_liquid]/(¹¹B/¹⁰B)_{NBS 951}}*1000 are small and range from 0.2 (± 0.7) to 0.9 (± 0.5) ‰ at 450°C and from 0.1 (± 0.6) to 0.7 (± 0.6) ‰ at 400°C. The data indicate increasing isotopic fractionation with decreasing pressure (i.e. increasing opening of the solvus). Extrapolation to salt saturated conditions yields maximum boron isotope liquid-vapor fractionations of Δ¹¹B_vapor-liquid = 1.5 (± 0.7) ‰ at 450°C and Δ¹¹B_vapor-liquid = 1.3 (± 0.6) ‰ at 400°C. The weak boron isotope fractionation suggests similar trigonal speciation in liquid and vapor. Although the boron and boron isotope fractionation between liquid and vapor is only weak, mass balance calculations indicate that for high degrees of fractionation liquid-vapor phase separation in an open system can significantly alter the boron and boron isotope signature of low-salinity hydrous fluids in hydrothermal systems. Comparing the model calculations with natural oceanic hydrothermal fluids, however, indicate that other processes than fluid phase separation dominate the boron geochemistry in oceanic hydrothermal fluids.     

Experimental determination of the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions to 400°C

The solubility of silver sulphide (acanthite/argentite) has been measured in aqueous sulphide solutions between 25 and 400°C at saturated water vapour pressure and 500 bar to determine the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions. The experiments were carried out in a flow-through autoclave, connected to a high-performance liquid chromatographic pump, titanium sampling loop, and a back-pressure regulator on line. Samples for silver determination were collected via the titanium sampling loop at experimental temperatures and pressures. The solubilities, measured as total dissolved silver, were in the range 1.0 × 10⁻⁷ to 1.30 × 10⁻⁴ mol kg⁻¹ (0.01 to 14.0 ppm), in solutions of total reduced sulphur between 0.007 and 0.176 mol kg⁻¹ and pH_T,p of 3.7 to 12.7. A nonlinear least squares treatment of the data demonstrates that the solubility of silver sulphide in aqueous sulphide solutions of acidic to alkaline pH is accurately described by the reactions0.5Ag₂S(s) + 0.5H₂S(aq) = AgHS(aq) K_s,1110.5Ag₂S(s) + 0.5H₂S(aq) + HS⁻ = Ag(HS)2− K_s,122Ag₂S(s) + 2HS⁻ = Ag₂S(HS)22− K_s,232where AgHS(aq) is the dominant species in acidic solutions, Ag(HS)2− under neutral pH conditions and Ag₂S(HS)22− in alkaline solutions. With increasing temperature the stability field of Ag(HS)2− increases and shifts to more alkaline pH in accordance with the change in the first ionisation constant of H₂S(aq). Consequently, Ag₂S(HS)22− is not an important species above 200°C. The solubility constant for the first reaction is independent of temperature to 300°C, with values in the range logK_s,111 = −5.79 (±0.07) to −5.59 (±0.09), and decreases to −5.92 (±0.16) at 400°C. The solubility constant for the second reaction increases almost linearly with inverse temperature from logK_s,122 = −3.97 (±0.04) at 25°C to −1.89 (±0.03) at 400°C. The solubility constant for the third reaction increases with temperature from logK_s,232 = −4.78 (±0.04) at 25°C to −4.57 (±0.18) at 200°C. All solubility constants were found to be independent of pressure within experimental uncertainties. The interaction between Ag⁺ and HS⁻ at 25°C and 1 bar to form AgHS(aq) has appreciable covalent character, as reflected in the exothermic enthalpy and small entropy of formation. With increasing temperature, the stepwise formation reactions become progressively more endothermic and are accompanied by large positive entropies, indicating greater electrostatic interaction. The aqueous speciation of silver is very sensitive to fluid composition and temperature. Below 100°C silver(I) sulphide complexes predominate in reduced sulphide solutions, whereas Ag⁺ and AgClOH⁻ are the dominant species in oxidised waters. In high-temperature hydrothermal solutions of seawater salinity, chloride complexes of silver(I) are most important, whereas in dilute hydrothermal fluids of meteoric origin typically found in active geothermal systems, sulphide complexes predominate. Adiabatic boiling of dilute and saline geothermal waters leads to precipitation of silver sulphide and removal of silver from solution. Conductive cooling has insignificant effects on silver mobility in dilute fluids, whereas it leads to quantitative loss of silver for geothermal fluids of seawater salinity.     

The solubility of platinum and gold in NaCl brines at 1.5 kbar, 600 to 800°C: A laser ablation ICP-MS pilot study of synthetic fluid inclusions

The concentration and distribution of Pt and Au in a fluid-melt system has been investigated by reacting the metals with S-free, single-phase aqueous brines (20, 50, 70 wt% eq. NaCl) ± peraluminous melt at a confining pressure of 1.5 kbar and temperatures of 600 to 800 °C, trapping the fluid in synthetic fluid inclusions (quartz-hosted) and vesicles (silicate melt-hosted), and quantifying the metal content of the trapped fluid and glass by laser ablation ICP-MS. HCl concentration was buffered using the assemblage albite-andalusite-quartz and fO₂ was buffered using the assemblage Ni-NiO. Over the range of experimental conditions, measured concentrations of Pt and Au in the brines (, ) are on on the order of 1-10³ ppm. Concentrations of Pt and Au in the melt (, ) are ∼35-100 ppb and ∼400-1200 ppb, respectively. Nernst partition coefficients (, ) are on the order of 10²-10³ and vary as a function of (non-Henry’s Law behavior). Trapped fluids show a significant range of metal concentrations within populations of inclusions from single experiments (∼ 1 log unit variability for Au; ∼2-3 log unit variability for Pt). Variability in metal concentration within single inclusion groups is attributed to premature brine entrapment (prior to metal-fluid-melt equilibrium being reached); this allows us to make only minimum estimates of metal solubility using metal concentrations from primary inclusions. The data show two trends: (i) maximum and average values of and in inclusions decrease ∼2 orders of magnitude as fluid salinity () increases from ∼4 to 40 molal (20 to 70 wt % eq. NaCl) at a constant temperature; (ii) maximum and average values of increase approximately 1 order of magnitude for every 100°C increase temperature at a fixed . The observed behavior may be described by the general expression:

     

An experimental study of the solubility of molybdenum in H2O and KCl-H2O solutions from 500 °C to 800 °C, and 150 to 300 MPa

The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H₂O and pure H₂O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H₂O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H₂O aqueous solutions. Similarly, in the pure H₂O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS₂).     

Uranophane dissolution and growth in CaCl2-SiO2(aq) test solutions

The dissolution and growth of uranophane [Ca(UO₂)₂(SiO₃OH)₂·5H₂O] have been examined in Ca- and Si-rich test solutions at low temperatures (20.5 ± 2.0 °C) and near-neutral pH (∼6.0). Uranium-bearing experimental solutions undersaturated and supersaturated with uranophane were prepared in matrices of ∼10⁻² M CaCl₂ and ∼10⁻³ M SiO₂(aq). The experimental solutions were reacted with synthetic uranophane and analyzed periodically over 10 weeks. Interpretation of the aqueous solution data permitted extraction of a solubility constant for the uranophane dissolution reaction and standard state Gibbs free energy of formation for uranophane ( kJ mol⁻¹).     

Experimental determination of the solubility of natural wollastonite in pure water up to pressures of 5 GPa and at temperatures of 400-800 °C

The solubility of natural, near-end-member wollastonite-I (>99.5% CaSiO₃) has been determined at temperatures from 400 to 800 °C and pressures between 0.8 and 5 GPa in piston-cylinder apparatus with the weight-loss method. Chemical analysis of quench products and optical monitoring in a hydrothermal diamond anvil cell demonstrates that no additional phases form during dissolution. Wollastonite-I, therefore, dissolves congruently in the pressure-temperature range investigated. The solubility of CaSiO₃ varies between 0.175 and 13.485 wt% and increases systematically with both temperature and pressure up to 3.0 GPa. Above 3.0 GPa wollastonite-I reacts rapidly to the high-pressure modification wollastonite-II. No obvious trends are evident in the solubility of wollastonite-II, with values between 1.93 and 10.61 wt%. The systematics of wollastonite-I solubility can be described well by a composite polynomial expression that leads to isothermal linear correlation with the density of water. The molality of dissolved wollastonite-I in pure water is then

log(m_woll)=2.2288-3418.23×T^-1+671386.84×T^-2+logρ_H2O×(5.4578+2359.11×T^-1).     

Experimental determination of synthetic NdPO4 monazite end-member solubility in water from 21°C to 300°C: implications for rare earth element mobility in crustal fluids

The solubility of synthetic NdPO₄ monazite end-member has been determined experimentally from 21 to 300°C in aqueous solutions at pH = 2, and at 21°C and pH = 2 for GdPO₄. Measurements were performed in batch reactors, with regular solution sampling for pH measurement, rare earths and phosphorous analysis by inductively coupled plasma mass spectrometry (ICP-MS) coupled with a desolvation system. Scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS) were employed to check that no reprecipitation of secondary phases occurred and that the mineral surfaces remained those of a monazite. Coupled with speciation calculations, measured solution compositions permitted the determination of NdPO₄ and GdPO₄ solubility products which are in general agreement with previous experimental determination on rhabdophane at 25°C, but showing that monazite is more than two orders of magnitude less soluble than inferred on the basis of previous thermodynamic estimates. The temperature evolution from 21 to 300°C of the equilibrium constant (K) of the NdPO₄ monazite end-member dissolution reaction given by:

     

CO2-H2O mixtures in the geological sequestration of CO2. II. Partitioning in chloride brines at 12-100°C and up to 600 bar

Correlations presented by Spycher et al. (2003) to compute the mutual solubilities of CO₂ and H₂O are extended to include the effect of chloride salts in the aqueous phase. This is accomplished by including, in the original formulation, activity coefficients for aqueous CO₂ derived from several literature sources, primarily for NaCl solutions. Best results are obtained when combining the solubility correlations of Spycher et al. (2003) with the activity coefficient formulation of Rumpf et al. (1994) and Duan and Sun (2003), which can be extended to chloride solutions other than NaCl. This approach allows computing mutual solubilities in a noniterative manner with an accuracy typically within experimental uncertainty for solutions up to 6 molal NaCl and 4 molal CaCl₂.     

Solubility of KFe(CrO4)2·2H2O at 4–75°C

The solubility of KFe(CrO₄)₂·2H₂O, a precipitate recently identified in a Cr(VI)-contaminated soil, was studied in dissolution and precipitation experiments. Ten dissolution experiments were conducted at 4–75°C and initial pH values between 0.8 and 1.2 using synthetic KFe(CrO₄)₂·2H₂O. Four precipitation experiments were conducted at 25°C with final pH values between 0.16 and 1.39. The log K_SP for the reaction

相似文献   

Experimental determination of the diffusion coefficient for calcium in olivine between 900°C and 1500°C

Experiments have been carried out to determine the temperature, oxygen fugacity (fO₂) and compositional dependence of the tracer diffusion coefficient (D) of calcium in olivine. These data constrain the diffusion coefficient over the temperature range 900 to 1500°C for the three principal crystallographic axes. Well constrained linear relationships between the reciprocal of the absolute temperature and log(D) exist at any given oxygen fugacity. There is a strong dependence of the diffusion coefficient on oxygen fugacity with D ∝ fO₂^(1/3). This makes a knowledge of the T-fO₂ path followed by geological samples a prerequisite for modelling Ca diffusion in olivine. The best fitting preexponential factor (Do) and activation energy (E) to the Arrhenius equation log (D) = log [Do exp(−E/RT)] + 0.31Δ log fO₂ for Ca diffusion in olivine at a given oxygen fugacity (fO₂*) are given by:diffusion along [100]: log [Do (m²/s)] = −10.78 ± 0.43; E = 193 ± 11 kJ/moldiffusion along [010]: log [Do (m²/s)] = −10.46 ± 0.37; E = 201 ± 10 kJ/moldiffusion along [001]: log [Do (m²/s)] = −10.02 ± 0.29; E = 207 ± 8 kJ/molwhere Δ log fO₂ = log[fO₂*] − log[10⁻¹²] with fO₂* in units of bars. There is no measurable compositional dependence of the diffusion coefficient between Fo₈₃ and Fo₉₂. Diffusion in Fo₁₀₀ has a much higher activation energy than in Fe-bearing olivine and has a weaker fO₂ dependence.     

Solubility of grossular, Ca3Al2Si3O12, in H2O-NaCl solutions at 800 °C and 10 kbar, and the stability of garnet in the system CaSiO3-Al2O3-H2O-NaCl

The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H₂O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca₃Al₂Si₃O₁₂ with increasing NaCl mole fraction (X_NaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H₂O value over the investigated range, indicating strong solute interaction with NaCl. The Ca₃Al₂Si₃O₁₂ mole fraction versus NaCl mole fraction curve has a broad plateau between X_NaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H₂O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO₃-Al₂O₃-H₂O-NaCl. At high Al₂O₃/CaSiO₃ bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H₂O, but disappears once grossular is nucleated. At X_NaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO₃: Al₂O₃ molality increases from 0.0013 in initially pure H₂O to near 0.15 at X_NaCl = 0.4 in CaSiO₃-saturated solutions, a >100-fold enhancement. In contrast, addition of Al₂O₃ to wollastonite-saturated NaCl solutions increases CaSiO₃ molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al₂O₃ solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.     

A simple predictive model of quartz solubility in water-salt-CO2 systems at temperatures up to 1000 °C and pressures up to 1000 MPa

Knowledge of the solubility of quartz over a broad spectrum of aqueous fluid compositions and T-P conditions is essential to our understanding of water-rock interaction in the Earth’s crust. We propose an equation to compute the molality of aqueous silica, m_SiO2(aq), mol·(kg H₂O)⁻¹, in equilibrium with quartz and water-salt-CO₂ fluids, as follows:

     

Kinetics of the reaction MgO + Al2O3 → MgAl2O4 and Al-Mg interdiffusion in spinel at 1200 to 2000°C and 1.0 to 4.0 GPa

The rate of spinel (MgAl₂O₄) growth at the interface between MgO and Al₂O₃ was investigated systematically at temperatures of 1200° to ∼2000°C and pressures between 1.0 and 4.0 GPa with a solid-media, piston-cylinder apparatus. As reported in previous 1-atm studies, the thickness (ΔX) of the spinel layer increases linearly with the square root of time for experiments differing only in duration, irrespective of pressure-temperature (P-T) conditions. The reaction rate constant (k = ΔX²/2t) is log-linear in 1/T and also in pressure. The apparent activation energy of 410 kJ/mol is independent of pressure; the apparent activation volume increases systematically with increasing temperature. Electron microprobe traverses across the spinel layer reveal a significant Al excess and charge-compensating Mg deficit near the spinel/corundum interface. This nonstoichiometry is promoted by high temperatures (>1500°C), suppressed by high pressures and varies linearly across the spinel to a near-stoichiometric composition at the interface with periclase. The Al and Mg composition gradients can be used to extract interdiffusion coefficients for Al ↔ Mg exchange through the spinel, which are described by D?=2.5×10⁻⁶ exp(−28200/T) m²sThese diffusivities differ substantially from the reaction rate constant k, reflecting the fact that k is a combination of the diffusivity and the reaction potential as indicated by the difference in spinel composition across the spinel layer (i.e., coexisting with corundum vs. coexisting with periclase). A simple model can be used to separate the two effects and show that the reaction potential (i.e., the MgO-Al₂O₃ phase diagram) is sensitive to changes in both temperature and pressure, whereas the governing diffusivity depends only on temperature.     

Controls on gold solubility in arc magmas: An experimental study at 1000 °C and 4 kbar

In order to (1) explain the worldwide association between epithermal gold-copper-molybdenum deposits and arc magmas and (2) test the hypothesis that adakitic magmas would be Au-specialized, we have determined the solubility of Au at 4 kbar and 1000 °C for three intermediate magmas (two adakites and one calc-alkaline composition) from the Philippines. The experiments were performed over a fO₂ range corresponding to reducing (∼NNO−1), moderately oxidizing (∼NNO+1.5) and strongly oxidizing (∼NNO+3) conditions as measured by solid Ni-Pd-O sensors. They were carried out in gold containers, the latter serving also as the source of gold, in presence of variable amounts of H₂O and, in a few additional experiments, of S. Concentrations of Au in glasses were determined by LA-ICPMS. Gold solubility in melt is very low (30-240 ppb) but increases with fO₂ in a way consistent with the dissolution of gold as both Au¹⁺ and Au³⁺ species. In the S-bearing experiments performed at ∼NNO−1, gold solubility reaches much higher values, from ∼1200 to 4300 ppb, and seems to correlate with melt S content. No systematic difference in gold solubility is observed between the adakitic and the non-adakitic compositions investigated. Oxygen fugacity and the sulfur concentration in melt are the main parameters controlling the incorporation and concentration of gold in magmas. Certain adakitic and non-adakitic magmas have high fO₂ and magmatic S concentrations favorable to the incorporation and transport of gold. Therefore, the cause of a particular association between some arc magmas and Au-Cu-Mo deposits needs to be searched in the origin of those specialized magmas by involvement of Au- and S-rich protoliths. The subducted slab, which contains metal-rich massive sulfides, may constitute a potentially favorable protolith for the genesis of magmas specialized with respect to gold.     

CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100°C and up to 600 bar

Evaluating the feasibility of CO₂ geologic sequestration requires the use of pressure-temperature-composition (P-T-X) data for mixtures of CO₂ and H₂O at moderate pressures and temperatures (typically below 500 bar and below 100°C). For this purpose, published experimental P-T-X data in this temperature and pressure range are reviewed. These data cover the two-phase region where a CO₂-rich phase (generally gas) and an H₂O-rich liquid coexist and are reported as the mutual solubilities of H₂O and CO₂ in the two coexisting phases. For the most part, mutual solubilities reported from various sources are in good agreement. In this paper, a noniterative procedure is presented to calculate the composition of the compressed CO₂ and liquid H₂O phases at equilibrium, based on equating chemical potentials and using the Redlich-Kwong equation of state to express departure from ideal behavior. The procedure is an extension of that used by King et al. (1992), covering a broader range of temperatures and experimental data than those authors, and is readily expandable to a nonideal liquid phase. The calculation method and formulation are kept as simple as possible to avoid degrading the performance of numerical models of water-CO₂ flows for which they are intended. The method is implemented in a computer routine, and inverse modeling is used to determine, simultaneously, (1) new Redlich-Kwong parameters for the CO₂-H₂O mixture, and (2) aqueous solubility constants for gaseous and liquid CO₂ as a function of temperature. In doing so, mutual solubilities of H₂O from 15 to 100°C and CO₂ from 12 to 110°C and up to 600 bar are generally reproduced within a few percent of experimental values. Fugacity coefficients of pure CO₂ are reproduced mostly within one percent of published reference data.