The behavior of apatite during crustal anatexis: Equilibrium and kinetic considerations

The solubility and dissolution kinetics of apatite in felsic melts at 850°–1500°C have been examined experimentally by allowing apatite crystals to partially dissolve into apatite-undersaturated melts containing 0–10 wt% water. Analysis of P and Ca gradients in the crystal/melt interfacial region enables determination of both the diffusivities and the saturation levels of these components in the melt. Phosphorus diffusion was identified as the rate-limiting factor in apatite dissolution. Results of four experiments at 8 kbar run in the virtual absence of water yield an activation energy (E) for P diffusion of 143.6 ± 2.8 kcal-mol^?1 and frequency factor (D₀) of 2.23^+2.88_?1.26 × 10⁹cm²-sec^?1. The addition of water causes dramatic and systematic reduction of both E and D₀ such that at 6 wt% H₂O the values are ~25 kcal-mol^?1 and 10^?5 cm²-sec^?1, respectively. At 1300°C, the diffusivity of P increases by a factor of 50 over the first 2% of water added to the melt, but rises by a factor of only two between 2 and 6%, perhaps reflecting the effect of a concentration-dependent mechanism of H₂O solution. Calcium diffusion gradients do not conform well to simple diffusion theory because the release of calcium at the dissolving crystal surface is linked to the transport rate of phosphorus in the melt, which is typically two orders of magnitude slower than Ca. Calcium chemical diffusion rates calculated from the observed gradients are about 50 times slower than calcium tracer diffusion.Apatite solubilities obtained from these experiments, together with previous results, can be described as a function of absolute temperature (T) and melt composition by the expression: In D^apatite/melt_P = [(8400 + ((SiO₂ ? 0.5)2.64 × 10⁴))/T] ? [3.1 + (12.4(SiO₂ ? 0.5))] where SiO₂ is the weight fraction of silica in the melt. This model appears to be valid between 45% and 75% SiO₂, 0 and 10% water, and for the range of pressures expected in the crust.The diffusivity information extracted from the experiments can be directly applied to several problems of geochemical interest, including I) dissolution times for apatite during crustal anatexis, and 2) pileup of P, and consequent local saturation in apatite, at the surfaces of growing major-mineral phases.     

Solubility of sulfur in some magmas at 1 atmosphere

In order to understand the distribution of sulfur in igneous rooks, we determined the solubility of sulfur in volcanic rock melts (tholeiite basalt, hawaiite and rhyodacite from Hawaii) at various gas compositions and at 1250° and 1300°C and 1 atm total pressure. The solubility of sulfur in the melt passes through a minimum with change in oxygen partial pressure, if other factors are held constant. For the basaltic liquid at 1200°C, most sulfur in the melt is as dissolved sulfide (S^?2) at oxygen partial pressures below 10^?8 atm and as dissolved sulfate at oxygen partial pressures above 10^?8 atm. Based on the present solubility data, 5 per cent is inferred for volcanic gas at 1 atm total pressure in equilibrium with subaerially extruded Hawaiian tholeiite basalt (Pele's hair with 180 ppm S) at 1200°C and 10^?8 atm P_O2.     

Surface chemistry and dissolution kinetics of glassy rocks at 25°C

The weathering rates and mechanisms of three types of glassy rocks were investigated experimentally at 25 °C, pH 1.0 to 6.2, and reaction times as much as to 3 months. Changes in major element chemistry were monitored concurrently as a function of time in the aqueous solution and within the near surface region of the glass. Leach profiles, obtained by a HF leaching technique, displayed near-surface zones depleted in major cations. These zones increased in depth with increasing time and decreasing pH of reactions. Release rates into the aqueous solution were parabolic for Na and K and linear for Si and Al. A coupled weathering model, involving surface dissolution with concurrent diffusion of Na, K, and Al, produced a mass balance between the aqueous and glass phases. Steady state conditions are reached at pH 1.0 after approximately 3 weeks of reaction. Steady-state is not reached even after 3 months at pH 6.2.An interdiffusion model describes observed changes in Na diffusion profiles for perlite at pH 1.0. The calculated Na self-diffusion coefficient of 5 × 10^?19 cm²·s^?1 at 25°C approximates coefficients extrapolated from previously reported high temperature data for obsidian. The self-diffusion coefficient for H₃O⁺, 1.2 × 10^?20 cm²·s^?1, is similar to measured rates of water diffusion during hydration of obsidian to form perlite.     

Solubility of manganotantalite, zircon and hafnon in highly fluxed peralkaline to peraluminous pegmatitic melts

The behavior of tantalum and zirconium in pegmatitic systems has been investigated through the determination of Ta and Zr solubilities at manganotantalite and zircon saturation from dissolution and crystallization experiments in hydrous, Li-, F-, P- and B-bearing pegmatitic melts. The pegmatitic melts are synthetic and enriched in flux elements: 0.7–1.3 wt% Li₂O, 2–5.5 wt% F, 2.8–4 wt% P₂O₅ and 0–2.8 wt% B₂O₃, and their aluminum saturation index ranges from peralkaline to peraluminous (ASI_Li = Al/[Na + K + Li] = 0.8 to 1.3) with various K/Na ratios. Dissolution and crystallization experiments were conducted at temperatures varying between 700 and 1,150°C, at 200 MPa and nearly water-saturated conditions. For dissolution experiments, pure synthetic, end member manganotantalite and zircon were used in order to avoid problems with slow solid-state kinetics, but additional experiments using natural manganotantalite and zircon of relatively pure composition (i.e., close to end member composition) displayed similar solubility results. Zircon and manganotantalite solubilities considerably increase from peraluminous to peralkaline compositions, and are more sensitive to changes in temperature or ASI of the melt than to flux content. A model relating the enthalpy of dissolution of manganotantalite to the ASI_Li of the melt is proposed: ∆H _diss (kJ/mol) = 304 × ASI_Li − 176 in the peralkaline field, and ∆H _diss (kJ/mol) = −111 × ASI_Li + 245 in the peraluminous field. The solubility data reveal a small but detectable competitivity between Zr and Ta in the melt, i.e., lower amounts of Zr are incorporated in a Ta-bearing melt compared to a Ta-free melt under the same conditions. A similar behavior is observed for Hf and Ta. The competitivity between Zr (or Hf) and Ta increases from peraluminous to peralkaline compositions, and suggests that Ta is preferentially bonded to non-bridging oxygens (NBOs) with Al as first-neighbors, whereas Zr is preferentially bonded to NBOs formed by excess alkalies. As a consequence Zr/Ta ratios, when buffered by zircon and manganotantalite simultaneously, are higher in peralkaline melts than in peraluminous melts.     

Diffusion of 40Ar in biotite: Temperature,pressure and compositional effects

The measured radiogenic ⁴⁰Ar loss from sized biotite (56% annite) samples following isothermalhydrothermal treatment have provided model diffusion coefficients in the temperature interval 600°C to 750°C, calculated on the assumption that Ar transport proceeds parallel to cleavage. These data yield an array on an Arrhenius plot with a slope corresponding to an activation energy 47.0 ± 2 kcal-mol^?1 and a frequency factor of 0.077^+0.21_?0.06 cm²-sec^?1. Together with previous diffusion data for micas in the annitephlogopite series, our results indicate a strong compositional effect, with increasing $\text{Fe}\text{Mg}$ ratio corresponding to an increase in diffusivity. An effective diffusion radius of about 150 μm for biotite is inferred from the experimental data which compares favorably with that estimated from geological studies. A pressure effect on activation energy corresponding to an activation volume of about 14 cm³-mol^?1 is observed. These data yield closure temperature estimates for this biotite composition cooling at rates of 100°C-Ma^?1, 10°C-Ma^?1 and 1°C-Ma^?1 of 345°C, 310°C and 280°C, respectively. ${}^{40}\text{Ar}{}^{39}\text{Ar}$ age-spectrum analysis of a hydrothermally treated biotite yields a complex release pattern casting doubt on the general usefulness of such measurements for geochronological purposes.     

The effect of halogens on Zr diffusion and zircon dissolution in hydrous metaluminous granitic melts

Diffusion of Zr and zircon solubility in hydrous, containing approximately 4.5 wt% H₂O, metaluminous granitic melts with halogens, either 0.35 wt% Cl (LCl) or 1.2 wt% F (MRF), and in a halogen-free melt (LCO) were measured at 1.0 GPa and temperatures between 1,050 and 1,400 °C in a piston-cylinder apparatus using the zircon dissolution technique. Arrhenius equations for Zr diffusion in each hydrous melt composition are, for LCO with 4.4ǂ.4 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaI4aaocqaI4aaocqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaIZaWmcqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabiIda4aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIXaqmcqaI0aancqaIWaamcqGGUaGlcqaIXa % qmcqGHXcqScqaIZaWmcqaIZaWmcqGGUaGlcqaI5aqoaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!571F! D = 2.88 ±0.03x10 ^{- 8} exp( [( - 140.1 ±33.9)/(RT)] )D = 2.88 \pm 0.03x10^{ - 8} \exp \left( {{{ - 140.1 \pm 33.9} \over {RT}}} \right) , for LCl with 4.5ǂ.5 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaIZaWmcqaIZaWmcqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaI1aqncqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabisda0aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIYaGmcqaI1aqncqaI0aancqGGUaGlcqaI4a % aocqGHXcqScqaI2aGncqaI0aancqGGUaGlcqaIXaqmaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!5719! D = 2.33 ±0.05x10 ^{- 4} exp( [( - 254.8 ±64.1)/(RT)] )D = 2.33 \pm 0.05x10^{ - 4} \exp \left( {{{ - 254.8 \pm 64.1} \over {RT}}} \right) and for MRF with 4.9ǂ.3 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaI1aqncqaI0aancqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaIZaWmcqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabiwda1aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIYaGmcqaIYaGmcqaIZaWmcqGGUaGlcqaI4a % aocqGHXcqScqaIXaqmcqaI1aqncqGGUaGlcqaI1aqnaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!5715! D = 2.54 ±0.03x10 ^{- 5} exp( [( - 223.8 ±15.5)/(RT)] )D = 2.54 \pm 0.03x10^{ - 5} \exp \left( {{{ - 223.8 \pm 15.5} \over {RT}}} \right) . Solubilities determined by the dissolution technique were reversed for LCO +4.5ǂ.5 wt% H₂O by crystallization of a Zr-enriched glass of LCO composition at 1,200 and 1,050 °C at 1.0 GPa. The solubility data were used to calculate partition coefficients of Zr between zircon and hydrous melt, which are given by the following expressions: for LCO % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGOnayJa % eG4mamZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 1aqncqGGUaGlcqaI4aaocqaI3aWnaaa!5924! lnD_Zr^zircon/melt = 1.63( [10000/(T)] ) - 5.87\ln D_{Zr}^{zircon/melt} = 1.63\left( {{{10000} \over T}} \right) - 5.87 , for LCl % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGinaqJa % eG4naCZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 0aancqGGUaGlcqaI3aWncqaI1aqnaaa!5920! lnD_Zr^zircon/melt = 1.47( [10000/(T)] ) - 4.75\ln D_{Zr}^{zircon/melt} = 1.47\left( {{{10000} \over T}} \right) - 4.75 and, for MRF by % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGinaqJa % eG4naCZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 0aancqGGUaGlcqaI5aqocqaIXaqmaaa!591C! lnD_Zr^zircon/melt = 1.47( [10000/(T)] ) - 4.91\ln D_{Zr}^{zircon/melt} = 1.47\left( {{{10000} \over T}} \right) - 4.91 . Experiments on the same compositions, but with water contents down to 0.5 wt%, demonstrated reductions in both the diffusion coefficient of Zr and zircon solubility in the melt. The addition of halogens at the concentration levels studied to metaluminous melts has a small effect on either the diffusion of Zr in the melt, or the solubility of zircon at all water concentrations and temperatures investigated. At 800 °C, the calculated diffusion coefficient of Zr is lowest in LCl, 9᎒^-17 m² s^-1, and is highest in LCO, 4᎒^-15 m² s^-1. Extrapolation of the halogen-free solubility data to a magmatic temperature of 800 °C yields solubilities of approximately one-third of those directly measured in similar compositions, predicted by earlier studies of zircon dissolution and based upon analyses of natural rocks. This discrepancy is attributed to the higher oxygen fugacity of the experiments of this study compared with previous studies and nature, and the effect of oxygen fugacity on the structural role of iron in the melt, which, in turn, affects zircon solubility, but does not significantly affect Zr diffusion.     

In situ observation of crystal growth in a basalt melt and the development of crystal size distribution in igneous rocks

The speciation of CO₂ in dacite, phonolite, basaltic andesite, and alkali silicate melt was studied by synchrotron infrared spectroscopy in diamond anvil cells to 1,000 °C and more than 200 kbar. Upon compression to 110 kbar at room temperature, a conversion of molecular CO₂ into a metastable carbonate species was observed for dacite and phonolite glass. Upon heating under high pressure, molecular CO₂ re-appeared. Infrared extinction coefficients of both carbonate and molecular CO₂ decrease with temperature. This effect can be quantitatively modeled as the result of a reduced occupancy of the vibrational ground state. In alkali silicate (NBO/t = 0.98) and basaltic andesite (NBO/t = 0.42) melt, only carbonate was detected up to the highest temperatures studied. For dacite (NBO/t = 0.09) and phonolite melts (NBO/t = 0.14), the equilibrium CO₂ + O^2? = CO₃ ^2? in the melt shifts toward CO₂ with increasing temperature, with ln K = ?4.57 (±1.68) + 5.05 (±1.44) 10³ T ^?1 for dacite melt (ΔH = ?42 kJ mol^?1) and ln K = ?6.13 (±2.41) + 7.82 (±2.41) 10³ T ^?1 for phonolite melt (ΔH = ?65 kJ mol^?1), where K is the molar ratio of carbonate over molecular CO₂ and T is temperature in Kelvin. Together with published data from annealing experiments, these results suggest that ΔS and ΔH are linear functions of NBO/t. Based on this relationship, a general model for CO₂ speciation in silicate melts is developed, with ln K = a + b/T, where T is temperature in Kelvin and a = ?2.69 ? 21.38 (NBO/t), b = 1,480 + 38,810 (NBO/t). The model shows that at temperatures around 1,500 °C, even depolymerized melts such as basalt contain appreciable amounts of molecular CO₂, and therefore, the diffusion coefficient of CO₂ is only slightly dependent on composition at such high temperatures. However, at temperatures close to 1,000 °C, the model predicts a much stronger dependence of CO₂ solubility and speciation on melt composition, in accordance with available solubility data.     

The importance of defining chemical potentials,substitution mechanisms and solubility in trace element diffusion studies: the case of Zr and Hf in olivine

The diffusion, substitution mechanism and solubility limits of Zr and Hf in synthetic forsterite (Mg₂SiO₄) and San Carlos olivine (Mg_0.9Fe_0.1)₂SiO₄ have been investigated between 1,200 and 1,500 °C as a function of the chemical potentials of the components in the system MgO(FeO)–SiO₂–ZrO₂(HfO₂). The effect of oxygen fugacity and crystallographic orientation were also investigated. The solubilities of Zr in forsterite are highest and diffusion fastest when the coexisting three-phase source assemblage includes ZrSiO₄ (zircon) or HfSiO₄ (hafnon), and lower and slower, respectively, when the source assemblage includes MgO (periclase). This indicates that Zr and Hf substitute on the octahedral sites in olivine, charge balanced by magnesium vacancies. Diffusion is anisotropic, with rates along the crystal axes increasing in the order a < b < c. The generalized diffusion relationship as a function of chemical activity (as \(a_{{{\text{SiO}}_{2} }}\)), orientation and temperature is: \(logD_{\text{Zr}} = \frac{1}{4}loga_{{{\text{SiO}}_{2} }} + logD_{0} - \left( {\frac{{368 \pm 17\;{\text{kJ}}\;{\text{mol}}^{ - 1} }}{{2.303\;{\text{RT}}}}} \right)\) where the values of log D ₀ are ?3.8(±0.5), ?3.4(±0.5) and ?3.1(±0.5) along the a, b and c axes, respectively. Most experiments were conducted in air (fO₂ = 10^?0.68 bars), but one at fO₂ = 10^?11.2 bars at 1,400 °C shows no resolvable effect of oxygen fugacity on Zr diffusion. Hf is slightly more soluble in olivine than Zr, but diffuses slightly slower. Diffusivities of Zr in experiments in San Carlos olivine at 1,400 °C, fO₂ = 10^?6.6 bars are similar to those in forsterite at the same conditions, showing that the controls on diffusivities are adequately captured by the simple system (nominally iron-free) experiments. Diffusivities are in good agreement with those measured by Spandler and O’Neill (Contrib Miner Petrol 159:791–818, 2010) in San Carlos olivine using silicate melt as the source at 1,300 °C, and fall within the range of most measurements of Fe–Mg inter-diffusion in olivine at this temperature. Forsterite–melt partitioning experiments in the CaO–MgO–Al₂O₃–SiO₂–ZrO₂/HfO₂ show that the interface concentrations from the diffusion experiments represent true equilibrium solubilities. Another test of internal consistency is that the ratios of the interface concentrations between experiments buffered by Mg₂SiO₄ + Mg₂Si₂O₆ + ZrSiO₄ or Mg₂SiO₄ + ZrSiO₄ + ZrO₂ (high silica activity) to those buffered by Mg₂SiO₄ + MgO + ZrO₂ (low silica activity) agree well with the ratios calculated from thermodynamic data. This study highlights the importance of buffering chemical potentials in diffusion experiments to provide constraints on the interface diffusant concentrations and hence validate the assumption of interface equilibrium.     

Calcite Dissolution in Deionized Water from 50°C to 250°C at 10 MPa: Rate Equation and Reaction Order

Carbonate minerals and water （or geofluids） reactions are important for modeling of geochemical processes and have received considerable attention over the past decades. The calcite dissolution rates from 50℃ to 250℃ at 10 MPa in deionized water with a flow rate varying from 0.2 to 5 mL/min were experimentally measured in a continuous flow column pressure vessel reactor. The dissolution began near the equilibrium with c/ceq 〉 0.3 and finally reached the equilibrium at 100℃-250℃, so the corresponding solubility was also determined as 1.87, 2.02, 2.02 and 1.88×10^-4.mol/L at 100℃, 150℃, 200℃ and 250℃ respectively, which was first increasing and then switching to decreasing with temperature and the maximum value might occur between 150℃ and 200℃. The experimental dissolution rate not only increased with temperature, but also had a rapid increase between 150℃ and 200℃ at a constant flow rate of 4 mL/min. The measured dissolution rates can be described using rate equations of R = k（1-c/ceq）n or R = kc-n. In these equations the reaction order n changed with temperature, which indicates that n was a variable rather than a constant, and the activation energy was 13.4 kJ/mol calculated with R = k（1-c/ceq）n or 18.0 kJ/mol with R = kc^-n, which is a little lower than the surface controlled values. The varied reaction order and lower activation energy indicates that calcite dissolution in this study is a complex interplay of diffusion controlled and surface controlled processes.     

Dissolution sequestration mechanism of CO<Subscript>2</Subscript> at the Shiqianfeng saline aquifer in the Ordos Basin,northwestern China

This study focused on the target injection layers of deep saline aquifers in the Shiqianfeng Fm. in the Carbon Capture and Sequestration (CCS) Demonstration Projects in the Ordos Basin, northwestern China. The study employed a combination method of experiments and numerical simulation to investigate the dissolution mechanism and impact factors of CO₂ in these saline aquifers. The results showed (1) CO₂ solubility in different types of water chemistry were shown in ascending order: MgCl₂-type water < CaCl₂-type water < Na₂SO₄-type water < NaCl-type water < Na₂CO₃-type water < distilled water. These results were consistent with the calculated results undertaken by TOUGHREACT with about 5% margin of error. CO₂ solubility of Shiqianfeng Fm. saline was 1.05 mol/L; (2) compared with distilled water, the more complex the water’s chemical composition, the greater the increase in HCO₃ ^?concentration. While the water’s composition was relatively simple, the tested water’s HCO₃ ^?concentrations were in close accord with the calculated value undertaken by the TOUGHREACT code, and the more complex the water’s composition, the poorer the agreement was, probably due to the complex and unstable HCO₃ ^? complicating matters when in an aqueous solution system including both tested HCO₃ ^?concentration and calculated HCO₃ ^?concentration; (3) the CO₂ solubility in the saline at the temperature conditions of 55 °C and 70 °C were 1.17 and 1.02 mol/L. When compared with the calculated value of 1.20 and 1.05 mol/L, they were almost the same with only 1 and 3% margin of error; concentrations of HCO₃ ^? were 402.73 mg/L (0.007 mol/L) and 385.65 mg/L (0.006 mol/L), while the simulation results were 132.16 mg/L (0.002 mol/L) and 128.52 mg/L (0.002 mol/L). From the contrast between the tested data and the calculated data undertaken by the TOUGHREACT code, it was shown that TOUGHRACT code could better simulate the interaction between saline and CO₂ in the dissolution sequestration capacity. Therefore, TOUGHREACT code could be used for the inter-process prediction of CO₂ long-term geological storage of CO₂; (4) The Ca²⁺ concentration and SO₄ ^2?concentration in saline water had less effect on the solubility of CO₂ and HCO₃ ^?concentration. In addition, TDS and pH values of saline affected not only the solubility of CO₂, but also the conversion of CO₂ to HCO₃ ^? due to that they can affect the activity and acid-base balance. So in fact, we just need to consider that the TDS and pH values are main impact factors in the dissolution sequestration capacity of CO₂ geological sequestration in deep saline aquifers.     

Concentration dependence of water diffusion in obsidian and dacitic melts at high-pressures

The diffusion of water in natural obsidian and model dacitic melts (Ab90Di8Wo2, mol %) has been studied at water vapor pressure up to 170 MPa, temperatures of 1200°C, H₂O contents in melts up to ～6 wt % using a high gas pressure apparatus equipped with a unique internal device. The experiments were carried out by a new low-gradient technique with application of diffusion hydration of a melt from fluid phase. The water solubility in the melts and its concentration along $C_{H_2 O} $ diffusion profiles were determined using IR microspectrometry by application of the modified Bouguer-Beer-Lambert equation. A structural-chemical model was proposed to calculate and predict the concentration dependence of molar absorption coefficients of the hydroxyl groups (OH^?) and water molecules (H₂O) in acid polymerized glasses (quenched melts) in the obsidian-dacite series. The water diffusion coefficients $D_{H_2 O} $ were obtained by the mathematical analysis of concentration profiles and the analytical solution of the second Fick diffusion law using the Boltzman-Matano method. It was shown experimentally that $D_{H_2 O} $ exponentially increases with a growth of water concentration in the obsidian and dacitic melts within the entire range of diffusion profiles. Based on the established quantitative correlation between $D_{H_2 O} $ and viscosity of such melts, a new method was developed to predict and calculate the concentration, temperature, and pressure dependences of $D_{H_2 O} $ in acid magmatic melts (obsidian, rhyolite, albite, granite, dacite) at crustal T, P parameters and water concentrations up to 6 wt %.     

New Zr–Hf Geothermometer for Magmatic Zircons

A geothermometer equation \(T = \frac{{1531}} {{\ln K_d + 0.883}}\), where \(K_{\dot d} = \frac{{X_{Zr}^S X_{Hf}^m }} {{X_{Zr}^m X_{Hf}^s }}\) [X _j ⁱ is the concentration (in ppm) of component i in phase j] is the Zr and Hf distribution coefficient between melt and zircon, and T is temperature in K, was derived by thermodynamic processing of literature experimental data on Zr and Hf distribution between acid melts (m) and zircon (s) and on the solubility of zircon and hafnon in the melts with variable silica content. In calculations with this equations, we assumed the Zr concentration in zircon to be constant: 480000 ppm. It is shown that the commonly observed increase in Hf concentration from the cores to margins of magmatic zircon crystals is caused by the fractional crystallization of zircon. For differentiated acid magmatic series, the initial crystallization temperature of zircon in the least silicic cumulates should be evaluated using the cores of large zircon grains with the highest Zr/Hf ratio. Application of the geothermometer for mafic and intermediate rocks may be hampered due to simultaneous crystallization of zircon with some other ore and mafic minerals relatively enriched in Zr and Hf. The newly derived geothermometer has some advantages over other indicators of the crystallization temperature of magmatic zircon based on the zircon saturation index (Watson and Harrison, 1983; Boehnke et al., 2013) and on Ti concentration in this mineral (Ferry and Watson, 2007) as it does not depend on the major-oxide melt composition and on the accuracy of the estimated SiO₂ and TiO₂ activities in the melts. Calculations of the Zr and Hf fractionation trends in the course of zircon crystallization in granitoid melts allow one to evaluate the temperature at which more evolved melt portions were segregated.     

Thio complexes of gold and the transport of gold in hydrothermal ore solutions

The solubility of gold in aqueous sulphide solutions has been determined from pH_20°C ≈ 4 to pH_20°C ≈ 9.5 in the presence of a pyrite-pyrrhotite redox buffer at temperatures from 160 to 300°C and 1000 bar pressure. Maximum solubilities were obtained in the neutral region of pH as, for example, with m_NaHS = 0.15 m, pH_20°C = 5.96, T = 309°C, P = 1000 bar where a gold solubility of 225 mg/kg was obtained. It was concluded that three thio gold complexes contributed to the solubility. The complex Au₂(HS)₂S^2? predominated in alkaline solution, the Au(HS)₂^? complex occurred in the neutral pH region, and in the acid pH region, it was concluded with less certainty that the Au(HS)° complex was present. Formation constants calculated forAu₂(HS)₂S^2? and Au (HS)₂^? emphasize their high stability. In the temperature range from 175 to 250°C, values of pβ for Au₂(HS)₂S^2? vary from ?53.0 to 47.9 (±1.6) and from ?23.1 to ?19.5 ( ± 1.5) for Au(HS)₂^?. Equilibrium constante for the dissolution reactions, $\text{Au° + H}{}_2\text{S + HS}{}^{\mathrm{?}}\text{? Au(HS)}{}_2{}^{\mathrm{?}}\text{+}\text{1}\text{2}\text{H}{}_2$ and 2Au° + H₂S + 2H8^? ? Au(HS)₂^? + H₂ vary from pK_m = +2.4 to +2.55 (±0.10) for Au₂(HS)₂S^2? and from pK_n = + 1.29 to + 1.19 (±0.10) for Au(HS)₂^? over the temperature range 175 to 250°C. Enthalpies of these dissolution reactions were calculated to be ΔH_m° = ?5.2 ±2.0 kcal/mol and ΔH_n° = +1.7 ±2.0 kcal/mol respectively. It was concluded that gold is probably transported in hydrothermal ore solutions as both thio and chloro complexes and may be deposited in response to changes in temperature, pressure, pH, oxidation potential of the system and total sulphur concentration.     

Celestite (SrSO4(s)) solubility in water,seawater and NaCl solution

Celestite solubility measurements have been conducted in pure water at temperatures from 10 to 90°C. Equilibrium was achieved with respect to a crystalline solid phase from both undersaturated and supersaturated solutions. The measurements show that the solubility undergoes a maximum near 20°C. LogK values for the solubility reaction are adequately described by the following expression over the temperature range 283.15 to 363.15 K: −logK= −35.3106+0.00422837T+318312/T²+14.99586 logT.The following thennodynamic values for the dissolution reaction of SrSO_4(s), at 25°C have been derived: ΔG_R⁰ = 37852 ± 30 Jmol⁻¹ΔH_R⁰ = −1668±920Jmol⁻¹ΔS_R⁰= −132.6±3.2JK^−1mol−1Celestite solubility measurements were also determined in NaCl solutions up to 5 m concentration and from 10 to 40°C. These data are in good agreement with the work of StrÜbel (1966), who reports solubility measurements to temperatures of 100°C.The application of the Pitzer relations and the solubility constants determined in this study to calculate celestite solubility in NaCl solutions yields excellent agreement between predicted values and experimental measurements over the entire range of temperature and NaCl concentration conditions. For the limited number of solubility measurements in seawater-type solutions and mixed-salt brines, the agreement using the Pitzer relations is within three percent of the measured solubility.     

Small-scale Hf isotopic variability in the Peninsula pluton (South Africa): the processes that control inheritance of source 176Hf/177Hf diversity in S-type granites

The ¹⁷⁶Hf/¹⁷⁷Hf composition of inherited and magmatic zircon in the 538 Ma S-type Peninsula pluton (South Africa) has been determined at different scales. In the smallest rock samples investigated (<0.5 dm³), as well as within individual thin sections, magmatic zircon crystals exhibit the same wide range in ε_Hf(538) as the pluton (8ε units). In addition, across a significant range of bulk-rock compositions, both the range and average of the magmatic zircon Hf isotopic composition do not vary significantly with compositional parameters that are expected to scale with the proportion of mantle-derived magma addition (e.g., Mg# and Ca). At all scales, the ε_Hf variability in the magmatic zircon fraction matches well with that portrayed by the time-evolved inherited zircon population [i.e., with the ε_Hf(538) range of the inherited zircon cores]. This evidence suggests that the ε_Hf heterogeneity of magmatic zircon is directly inherited from the source. However, the analysis of zircon core–rim pairs reveals that the ¹⁷⁶Hf/¹⁷⁷Hf composition of the inherited crystals does not directly transfer to their magmatic overgrowths. Small-scale modeling of zircon dissolution and re-precipitation in a static magma generates sub-mm melt domains having variable Zr content and Hf isotope composition. The composition of these domains is controlled by the size and isotope composition of the nearest dissolving zircon crystals and the cooling rate of the magma. These results suggest that in magma systems with a substantial inherited zircon load, zircon crystals within the same rock should record variable ¹⁷⁶Hf/¹⁷⁷Hf in the magmatic zircon fraction.     

Determination of the first ionization constant of silicic acid from quartz solubility in borate buffer solutions to 350°C

The solubility of quartz has been determined in borax buffer solutions having total boron concentrations of 0.10, 0.20, 0.40 and 0.60 mol kg^?1 and over the temperature range 130–350°C at the saturated vapour pressure of the system. The first ionization constant of silicic acid was calculated from the solubility data and varied from ?logK₁ = 8.88 (± 0.15) at 130°C to ?logK₁ = 10.06 (± 0.20) at 350°C. The solubility of quartz in these solutions was due to the presence of the three species, H₄SiO₄, H₃SiO₄^? and NaH₃SiO₄°. The equilibrium constant for the reaction, Na⁺ + H₃SiO₄^? = NaH₃SiO₄° extended from log K_as = 1.18?1.40 (± 0.20) over the temperature interval 135–301°C. The formation of NaH₃SiO₄° ion pairs was concluded to contribute significantly to the solubility of quartz in alkaline hydrothermal solutions when pH > 8 and sodium concentration exceeds 0.10 mol kg^?1.     

DIffusion of Eu and Gd in basalt and obsidian

Tracer diffusion coefficients of ¹⁵³Gd and ¹⁵²Eu in olivine tholeiite have been determined at temperatures between 1150 and 1440°C. The results are identical for both tracers within experimental error. Between 1440 and 1320°C the diffusion coefficients are given by D(Eu, Gd) = 0.058 exp(?40,600/ RT). Between 1320 and 1210°C, the diffusion coefficients are constant at D = (1.4 ± 0.4) × 10^?7 cm²s^?1 and between 1210 and 1150°C, the D values drop irregularly to 4 × 10^?9 cm²s^?1. The liquidus temperature (1270°C) lies within the region of constant D. Such anomalous behavior has not been encountered in previous studies of Ca, Sr, Ba and Co diffusion in basalt. To explain the constant D value near the liquidus, we speculate that the structure of the melt changes as a function of temperature in such a way that the normal temperature dependence of the diffusivity is compensated. For example, the rare earth ions may be displaced from their (high temperature) octahedral coordination sites to other sites where they are more readily dissociated and therefore become progressively more mobile. The behavior below 1210°C may be the result of relatively stable complexes or molecules in the melt or of the formation of a REE bearing crystalline phase that has so far escaped detection. Preliminary results for Eu diffusion in obsidian are D (Eu, 800°C) = 5 × 10^?13 cm² s^?1 and D (Eu, 950°C) = 1.5 × 10^?11 cm² s^?1. These data are consistent with an activation energy of 59 Kcal mole^?1. These low diffusivities indicate that the partitioning of REE in crystallizing intermediate and acidic melts may be controlled by diffusion in the melt rather than equilibrium between the crystal surface and the bulk melt.The diffusion data are applied to partial melting in the mantle, in an attempt to explain how LREE enriched tholeiites may be derived from a LREE depleted mantle source. In this model LREE diffuse from garnet bearing regions that have small melt fractions into garnet free regions that have relatively large melt fractions. REE diffusion is so slow that this process is quantitatively significant only in small partially molten bodies (diameter ～1 km or less) or in larger, but strongly flattened bodies. Internal convective motion during diapiric rise would also increase the efficiency of the process.     

Amorphous silica solubilities—I. Behavior in aqueous sodium nitrate solutions; 25–300°C, 0–6 molal

The solubility of amorphous silica was obtained in aqueous sodium nitrate solutions up to six molal and at temperatures from 25 to 300°C. It was expected that solubilities in aqueous sodium chloride solutions would be similar. At 25°C, the solubility of amorphous silica is lowered from that in water to 0.00086 m in 6.12 m sodium nitrate, or a decrease of 60%. At 300°C, the corresponding decrease is only 27% from a solubility of 0.0269 m in H₂O. From the change in solubility with temperature at a given constant molality of sodium nitrate, the molal heat of solution over the range, 100 to 300°C, increases from + 2.93 kcal mol^?1 in water to + 3.64 kcal mol^?1 in 6m sodium nitrate. The value approaches a constant of +3.8 kcal mol^?1 as sodium nitrate approaches saturation at 10.8 molal.     

The effect of composition on the ferric-ferrous ratio in basaltic liquids at atmospheric pressure

The effects on the ferric-ferrous ratio of varying individual components in a dry basaltic liquid have been determined at atmospheric pressure and constant oxygen fugacity (fO₂). Experiments were conducted by suspending 100 mg samples from pt loops at 1200°C (fO₂ = 10^?8atm) and 1360°C (fO₂ = 10^?6atm) in an atmosphere controlled by mixtures of CO₂ and H₂. A microanalytical wetchemical technique and the electron microprobe were used to determine the composition of the resulting basaltic glasses. In order of decreasing significance, the addition of oxides of K, Na, Si, Al, or Ca produces an increase in the ferric-ferrous ratio of the melt at 1200°C. The change in the ferric-ferrous ratio produced by component addition is less at 1360°C than at 1200°C.     

Pressure dependence ot SrSO4 solubility

The solubilities of SrSO₄ in seawater, 0.65 M NaCl and and distilled water were measured as a function of pressure at 2°C. The thermodynamic solubility product was determined from the distilled water measurements and stoichiometric solubility products were determined from the seawater and Nad measurements. The equilibrium quotient for SrSO₄ dissolution at ionic strength of 0.65 was calculated from the NaCl measurements, using the known NaSO₄^? ionpairing association constant. For each of the solubility products values of $\text{Θ V}$ were determined. These experimental values were all $\text{11.0 ± 0.3}\text{ml}$ mole_? lower than the theoretical values based on anhydrous SrSO₄. This difference may be due to the equilibrating solid phase being a hydrated form of SrSO₄.