Experimental hydrogen isotope studies—I. Systematics of hydrogen isotope fractionation in the systems epidote-H2O,zoisite-H2O and AlO(OH)-H2O

$\text{H}\text{D}$ Fractionation factors between epidote minerals and water, and between the AlO(OH) dimorphs boehmite and diaspore and water, have been determined between 150 and 650°C. Small water mineral ratios were used to minimise the effect of incongruent dissolution of epidote minerals. Waters were extracted and analysed directly by puncturing capsules under vacuum. Hydrogen diffusion effects were eliminated by using thick-walled capsules.$\text{H}\text{D}$ Exchange rates are very fast between epidote and water (and between boehmite and water), complete exchange taking only minutes above 450°C but several months at 250°C. Exchange between zoisite and water (and between diaspore and water) is very much slower, and an interpolation method was necessary to determine fractionation factors at 450 and below.For the temperature range 300–650°C, the $\text{H}\text{D}$ equilibrium fractionation factor (α^e) between epidote and water is independent of temperature and Fe content of the epidote, and is given by 1000 In α_epidote-H2O^e = ?35.9 ± 2.5, while below 300°C 1000 In , with a ‘cross-over’ estimated to occur at around 185°C. By contrast, zoisite-water fractionations fit the relationship 1000 In .All studied minerals have hydrogen bonding. Fractionations are consistent with the general relationship: the shorter the O-H -- O bridge, the more depleted is the mineral in D.On account of rapid exchange rates, natural epidotes probably acquired their H-isotope compositions at or below 200°C, where fractionations are near or above 0%.; this is in accord with the observation that natural epidotes tend to concentrate D relative to other coexisting hydrous minerals.     

Methane-producing bacteria: natural fractionations of the stable carbon isotopes

The ${}^{13}\text{C}{}^{12}\text{C}$ fractionation factors ($\text{CO}{}_2\text{CH}{}_4$) for the reduction of CO₂ to CH₄ by pure cultures of methane-producing bacteria are, for Methanosarcina barkeri at 40°C, 1.045 ± 0.002; for Methanobacterium strain M.o.H. at 40°C, 1.061 ± 0.002; and, for Methanobacterium thermoautotrophicum at 65°C, 1.025 ± 0.002. These observations suggest that the acetic acid used by acetate dissimilating bacteria, if they play an important role in natural methane production, must have an intramolecular isotopic fractionation ($\text{CO}{}_2\text{H}\text{CH}{}_3$) approximating the observed $\text{CO}{}_2\text{CH}{}_4$ fractionation.     

Dissociation constants for alkali earth and sodium borate ion pairs from 10 to 50°C

Potentiometric measurements in dilute sodium borate solutions with added alkali earth chlordie salts yield the following expressions for the dissociation constants of alkali earth borate ion pairs from 10 to 50°C:

$\text{p}\text{K(}\text{MgH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.266+0.001204 T}$

$\text{p}\text{K(}\text{CaH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.154+0.002170 T}$

$\text{p}\text{K(}\text{SrH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.033+0.001738 T}$

$\text{p}\text{K(}\text{BaH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.942+0.001850 T}$

where T is in °K. Enthalpies for the dissociation reactions at 25°C are less than 1 kcal./mole for all the alkali earth borate ion pairs.Values for pK(NaH₂BO₃°) from 5 to 55°C computed from the experimental data of Owen and King are in good agreement with those determined potentiometrically. The average value from both methods is 0.22 ± 0.1 at 25°C.Application to seawater of computed pK's for MgH₂BO₃⁺, CaH₂BO₃⁺ and NaH₂BO₃⁰ yields an apparent dissociation constant for boric acid of 8.73 vs. 8.70 measured by Lyman, 8.68 by Buch and 8.73 by Byrne and Kester.     

Carbonate equilibria in hydrothermal systems: First ionization of carbonic acid in NaCl media to 300°C

The ionization quotients of aqueous carbon dioxide (carbonic acid) have been precisely determined in NaCl media to 5 m and from 50° to 300°C using potentiometric apparatus previously developed at Oak Ridge National Laboratory. The pressure coefficient was also determined to 250°C in the same media. These results have been combined with selected information in the literature and modeled in two ways to arrive at the best fits and to derive the thermodynamic parameters for the ionization reaction, including the equilibrium constant, activity coefficient quotients, and pressure coefficients. The variation with temperature of the two fundamental quantities $\text{Δ}\text{V}\text{?}{}^{\mathrm{o}}$ and $\text{Δ}\text{C}\text{?}{}^{\mathrm{o}}{}_{\mathrm{p}}$ were examined along the saturation vapor pressure curve and at constant density. The results demonstrated again that for reactions with minimal electrostriction changes the magnitudes and variations of $\text{Δ}\text{C}\text{?}{}^{\mathrm{o}}{}_{\mathrm{p}}$ and $\text{Δ}\text{V}\text{?}{}^{\mathrm{o}}$ with temperature are small and, in addition, $\text{Δ}\text{C}\text{?}{}_{\mathrm{p}}$ and $\text{Δ}\text{V}\text{?}$ are approximately independent of salt concentration.The results have also been applied to an examination of the solubility of calcite as a function of pH (in a given NaCl medium) for the neutral to acidic region both for systems with fixed CO₂ pressure and systems where the calcium ion concentration equals the concentration of carbon. The pH of saturated solutions of calcite with P_CO2 of 12 bars increases from 5.1 to 5.5 between 100° and 300°C.     

Medium pressure crystallization of a monchiquitic magma — evidence from megacrysts of Drever''s block, Ubekendt Ejland, West Greenland

Megacrysts and polymineralic fragments of extraordinary diversity from a Tertiary monchiquitic dyke of Ubekendt Ejland comprise three groups: (1) Cr-diopside-fassaitic diopside + olivine, Fo_90.5?81.5 + CrAl spinels. (II) Fassaitic salite-ferrisalite + KTi-pargasite-ferropargasite + apatite + AlTi-magnetite, (III) Scapolite + hyalophane + potassium feldspar + nepheline + analcime. By comparison with mineralogy and phase relations in the host rock and experimental data from alkaline rocks the megacrysts are related to a sequence of crystallization from primitive monchiquitic to potassic phonolitic magmas rich in H₂O and CO₂ at 5–11 kb. Group I megacrysts formed at temperatures of 1300-1150°C and group II between ? bar at the latter temperature. High P_co2 may have stabilized the scapolite in the more evolved liquid and K-feldspar and nepheline began to crystallize at ca. 800°C possibly together with the ferrisalite.     

Some mineral stability relations in the system CaOMgOSiO2H2OHCl

Mineral-aqueous solution equilibria for the assemblages talc-quartz, tremolite-talc-quartz, diopside-tremolite-quartz, wollastonite-diopside-quartz and wollastonite-quartz have been studied at 2 kb total pressure, 500° to 700°C and chloride concentrations from 0.03 to 6.0 molal. Most work was at 1 m chloride. Both buffered and unbuffered data were obtained and a recalibration of the Ag-AgCl buffer is presented. Log equilibrium quotients at 500°, 600° and 700°C are respectively: Ta-Qz () 2.57, 1.71, 0.73; Tr-Ta-Qz and Di-Tr-Qz () 4.98, 3.99, 2.21 and 7.29, 5.30, 3.56; WoDi-Qz () 3.30, 3.00, 2.79: Wo-Qz () 5.15, 3.95, 2.68. Mineral stability fields plotted in terms of these concentration data more tangibly represent the compositional character of real systems and the mass transfer capabilities of their fluids than do the analogous theoretical activity diagrams.Overall dissociation constants of MgCl₂ and CaCl₂ were calculated from the experimental data using the calculated ionic activity constants for the reactions and the established dissociation constants of HCl. The negative log values are respectively: 3.88. 6.63, 9.20 for CaCl₂ and 4.60, 7.54, 10.37 for MgCl₂ at 500°, 600° and 700°C, 2 kb. The Ca values are about an order of magnitude more positive than the conductance-derived values by Frantz and Marshall (1982).The phase relations developed in this study have application to the genesis of talc, tremolite, and diopside-bearing assemblages in some regional metamorphic rocks, but more specifically to the calcsilicate skarn assemblages of many metasomatic aureoles. The equilibrium fluids are characterized by high concentrations of Ca relative to Mg and increasing $\text{Ca}\text{Mg}$ ratios with decreasing temperatures. The stability fields of talc, tremolite, and quartz expand relative to those of diopside and wollastonite with decreasing temperature, hence their more common appearance as retrograde products in skarn systems.     

Oxygen and sulfur isotope equilibria in the BaSO4HSO4−H2O system from 110 to 350°C and applications

Oxygen isotope exchange between BaSO₄ and H₂O from 110 to 350°C was studied using 1 m H₂SO₄-1 m NaCl and 1 m NaCl solutions to recrystallize the barite. The slow exchange rate (only 7% exchange after 1 yr at 110°C and 91% exchange after 22 days at 350°C in 1 m NaCl solution) prompted the use of the partial equilibrium technique. However, runs at 300 and 350°C were checked by complete exchange experiments. The temperature calibration curve for the isotope exchange is calculated giving most weight to the high temperature runs where the partial equilibrium technique can be tested. Oxygen isotope fractionation factors (α) in 1 m NaCl solution (110–350°C), assuming a value of 1.0407 for α_CO2H2O at 25°C, are:

.These data, when corrected for ion hydration effects in solution (Truesdell, 1974), give the fractionation factors in pure water:

.In the 1 m H₂SO₄-1 m NaCl runs, sulfur isotope fractionation between HSO^?₄ and BaSO₄ is less than the detection limit of 0.4%. A barite-sulfide geothermometer is obtained by combining HSO^?₄H₂S and sulfide-H₂S calibration data.Barite in the Derbyshire ore field, U.K., appears to have precipitated in isotopic equilibrium with water and sulfur in the ore fluid at temperatures less than 150°C. At the Tui Mine, New Zealand, the barite-water geothermometer indicates temperatures of late stage mineralization in the range 100–200°C. A temperature of 350 ± 20°C is obtained from the barite-pyrite geothermometer at the Yauricocha copper deposit, Peru, and oxygen isotope analyses of the barite are consistent with a magmatic origin for the ore fluids.     

Rutile solubility and titanium coordination in silicate melts

The solubility of rutile has been determined in a series of compositions in the K₂O-Al₂O₃-SiO₂ system ($\text{K}{}^1\text{=}\text{K}{}_2\text{O}\text{(K}{}_2\text{O}$ + Al₂O₃) = 0.38–0.90), and the CaO-Al₂O₃-SiO₂ system ($\text{C}{}^1\text{=}\text{CaO}\text{(CaO + Al}{}_2\text{O}{}_3\text{)}\text{= 0.47–0.59}$). Isothermal results in the KAS system at 1325°C, 1400°C, and 1475°C show rutile solubility to be a strong function of the $\text{K}{}^1$ ratio. For example, at 1475°C the amount of TiO₂ required for rutile saturation varies from 9.5 wt% ($\text{K}{}^1\text{= 0.38}$) to 11.5 wt% ($\text{K}{}^1\text{= 0.48}$) to 41.2 wt% ($\text{K}{}^1\text{= 0.90}$). In the CAS system at 1475°C, rutile solubility is not a strong function of $\text{C}{}^1$. The amount of TiO₂ required for saturation varies from 14 wt% ($\text{C}{}^1\text{= 0.48}$) to 16.2 wt% ($\text{C}{}^1\text{= 0.59}$).The solubility changes in KAS melts are interpreted to be due to the formation of strong complexes between Ti and K⁺ in excess of that needed to charge balance Al³⁺. The suggested stoichiometry of this complex is K₂Ti₂O₅ or K₂Ti₃O₇. In CAS melts, the data suggest that Ca²⁺ in excess of A1³⁺ is not as effective at complexing with Ti as is K⁺. The greater solubility of rutile in CAS melts when $\text{C}{}^1$ is less than 0.54 compared to KAS melts of equal $\text{K}{}^1$ ratio results primarily from competition between Ti and Al for complexing cations (Ca vs. K).TiK_β x-ray emission spectra of KAS glasses ($\text{K}{}^1\text{= 0.43–0.60}$) with 7 mole% added TiO₂, rutile, and Ba₂TiO₄, demonstrate that the average Ti-O bond length in these glasses is equal to that of rutile rather than Ba₂TiO₄, implying that Ti in these compositions is 6-fold rather than 4-fold coordinated. Re-examination of published spectroscopic data in light of these results and the solubility data, suggests that the 6-fold coordination polyhedron of Ti is highly distorted, with at least one Ti-O bond grossly undersatisfied in terms of Pauling's rules.     

Use of Pitzer's equations to determine the media effect on the formation of lead chloro complexes

The spectrophotometric measurements of chloro complexes of lead in aqueous HCl, NaCl, MgCl₂ and CaCl₂ solutions at 25°C have been analyzed using Pitzer's specific interaction equations. Parameters for activity coefficients of the complexes PbCl⁺, PbCl₂⁰ and PbCl₃^? have been determined for the various media. Values of $\text{K}{}_1\text{= 30.0 ± 0.6}$, $\text{K}{}_2\text{= 106.7 ± 2.1}$ and $\text{K}{}_3\text{= 73.0 ± 1.5}$ were obtained for the cumulative formation constants. [$\text{Pb}{}^{\mathrm{2+}}\text{+ nCl}{}^{\mathrm{?}}\text{→ PbCl}{}_{\mathrm{n}}{}^{\mathrm{2?n}}\text{)}$]. These values are in reasonable agreement with literature data. The Pitzer parameters for the PbCl ion pairs in various media were used to calculate the speciation of Pb²⁺ in an artificial seawater solution.     

The system pargasite-H2O-CO2: a model for melting of a hydrous mineral with a mixed-volatile fluid—I. Experimental results to 8 kbar

The stability of the amphibole pargasite [NaCa₂Mg₄Al(Al₂Si₆))O₂₂(OH)₂] in the melting range has been determined at total pressures (P) of 1.2 to 8 kbar. The activity of H₂O was controlled independently of P by using mixtures of H₂O + CO₂ in the fluid phase. The mole fraction of H₂O in the fluid ($\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$) ranged from 1.0 to 0.2.At P < 4 kbar the stability temperature (T) of pargasite decreases with decreasing $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at constant P. Above P ? 4 kbar stability T increases as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ is decreased below one, passes through a T maximum and then decreases with a further decrease in $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$. This behavior is due to a decrease in the H₂O content of the silicate liquid as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ decreases. The magnitude of the T maximum increases from about 10°C (relative to the stability T for $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{= 1}$) at P = 5 kbar to about 30°C at P = 8 kbar, and the position of the maximum shifts from $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.6}$ at P = 5 kbar to $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.4}$ at P = 8 kbar.The H₂O content of liquid coexisting with pargasite has been estimated as a function of at 5 and 8 kbar P, and can be used to estimate the H₂O content of magmas. Because pargasite is stable at low values of $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at high P and T, hornblende can be an important phase in igneous processes even at relatively low H₂O fugacities.     

Solubility product of galena at 298°K: A possible explanation for apparent supersaturation in nature

The synthetic chelating agent ethylenediaminetetraacetic acid (EDTA) has been used to evaluate the stoichiometric solubility product of galena (PbS) at 298°K: $\text{K}{}_{\mathrm{s2}}\text{=}{}^{\mathrm{a}}\text{Pb}{}^{\mathrm{2+}}{}^{\mathrm{a}}\text{HS}{}^{\mathrm{?}}{}^{\mathrm{a}}\text{H}{}^{\mathrm{+}}$ This method circumvents the possible uncertainties in the stoichiometry and stability of lead sulfide complexes. At infinite dilution, $\text{Log K}{}_{\mathrm{s2}}\text{= ?12.25 ±0.17}$, and at an ionic strength corresponding to seawater ($\text{I = 0.7 M}$), $\text{Log K}{}_{\mathrm{s2}}\text{= ?11.73 ± 0.05}$. Using the value of $\text{K}{}_{\mathrm{s2}}$ at infinite dilution, and the free energies of formation of HS^? and Pb²⁺ at 298°K (literature values), the free energy of formation of PbS at 298°K is computed to be $\text{?79.1 ± 0.8 KJ/mol (?18.9 Kcal/mol)}$. Galena is shown to be more than two orders of magnitude more soluble than indicated by calculations based on previous thermodynamic data.     

Diffusion of ions in sea water and in deep-sea sediments

The tracer-diffusion coefficient of ions in water, D_j⁰, and in sea water, $\text{D}{}_{\mathrm{j}}{}^1$, differ by no more than zero to 8 per cent. When sea water diffuses into a dilute solution of water, in order to maintain the electro-neutrality, the average diffusion coefficients of major cations become greater but of major anions smaller than their respective $\text{D}{}_{\mathrm{j}}{}^1$ or D_j⁰ values. The tracer diffusion coefficients of ions in deep-sea sediments, D_j,sed., can be related to $\text{D}{}_{\mathrm{j}}{}^1$ by $\text{D}{}_{\mathrm{j,sed.}}\text{= D}{}_{\mathrm{j}}{}^1\text{·}\text{α}\text{θ}{}^2$, where θ is the tortuosity of the bulk sediment and a a constant close to one.     

Transfer and exchange equilibria in a portion of the pyroxene quadrilateral as deduced from natural and experimental data

Differences in the chemical composition of metamorphic and igneous pyroxene minerals may be attributed to a transfer reaction, which determines the Ca content of the minerals, and an exchange reaction, which determines the relative Mg:Fe²⁺ ratios. Natural data for associated Ca pyroxene (Cpx) and orthopyroxene (Opx) or pigeonite are combined with experimental data for Fe-free pyroxenes, to produce the following equations for the Cpx slope of the solvus surface: $\text{> 1080°C: T =}\text{1000}\text{(0.468 + 0.246X}{}^{\mathrm{Cpx}}\text{? 0.123 ln (1–2 [Ca]))}$$\text{< 1080°C: T =}\text{1000}\text{(0.054 + 0.608X}{}^{\mathrm{Cpx}}\text{? 0.304 ln (1–2 [Ca]))}$, and the following equation for the temperature-dependence of the Mg-Fe distribution coefficient: $\text{T =}\text{1130}\text{(}\text{ln}\text{K}{}_{\mathrm{p}}\text{+ 0.505)}$, where T is absolute temperature, X is $\text{Fe}{}^{\mathrm{2+}}\text{(Mg + Fe}{}^{\mathrm{2+}}\text{))}$, [Ca] is $\text{Ca}\text{(Ca + Mg + Fe}{}^{\mathrm{2+}}\text{)}$ in Cpx, and K_D is the distribution coefficient, defined as $\text{X}{}^{\mathrm{Opx}}\text{/(1 ? X}{}^{\mathrm{Opx}}\text{) ÷ X}{}_{\mathrm{Cpx}}\text{/(1 ?}{}^{\mathrm{Cpx}}\text{)}$.The transfer and exchange equations form useful temperature indicators, and when applied to 9 sets of well-studied rocks, yield pairs of temperatures that are in good agreement. For example, temperatures obtained for the Bushveld Complex are 1020°C (solvus equation) and 980°C (exchange equation), based on 7 specimens. The uncertainty in these numbers, due to precision and accuracy errors, is estimated to be ±60°.     

An experimental study of alkali metal distributions in feldspars and micas

K and Rb distributions between aqueous alkali chloride vapour phase (0.7 molar) and coexisting phlogopites and sanidines have been investigated in the range 500 to 800°C at 2000 kg/cm² total pressure.Complete solid solution of RbMg₃AlSi₃O₁₀(OH)₂ in KMg₃AlSi₃O₁₀(OH)₂ exists at and above 700°C. At 500°C a possible miscibility gap between approximately 0.2 and 0.6 mole fraction of the Rb end-member is indicated.Only limited solid solution of Rb AlSi₃O₈ in KAlSi₃O₈ has been found at all temperatures investigated.Distribution coefficients, expressed as $\text{K}{}_{\mathrm{d}}\text{= (}\text{Rb/K}\text{)}$ in solid/(Rb/K) in vapour, are appreciably temperature-dependent but at each temperature are independent of composition for low Rb end-member mole fractions in the solids. The determined $\text{K}{}_{\mathrm{D}}$ values and their approximate Rb end-member mole fraction ($\text{X}{}_{\mathrm{RM}}$) ranges of constancy are summarized as follows: (°C)$\text{T}$$\text{K}{}_{\mathrm{D}}{}^{\mathrm{Phlog/Vap.}}$$\text{X}{}_{\mathrm{RM}}$$\text{K}{}_{\mathrm{D}}{}^{\mathrm{Sandi/Vap.}}$$\text{X}{}_{\mathrm{rm}}$

     

16.

Trace element partitioning and melt structure: An experimental study at 1 atm pressure     

Bjørn O. Mysen  David Virgo 《Geochimica et cosmochimica acta》1980,44(12):1917-1930

Diopside-melt and forsterite-melt rare earth (REE) and Ni partition coefficients have been determined as a function of bulk compositions of the melt. Available Raman spectroscopic data have been used to determine the structures of the melts coexisting with diopside and forsterite. The compositional dependence of the partition coefficients is then related to the structural changes of the melt.The melts in all experiments have a ratio of nonbridging oxygens to tetrahedral cations ($\text{NBO}\text{T}$) between 1 and 0. The quenched melts consist of structural units that have, on the average, 2 (chain), 1 (sheet) and 0 (three-dimensional network) nonbridging oxygens per tetrahedral cation. The proportions of these structural units in the melts, as well as the overall $\text{NBO}\text{T}$, change as a function of the bulk composition of the melt.It has been found that Ce, Sm, Tm and Ni crystal-liquid partition coefficients ($\text{K}{}^{\mathrm{crystal?liq}}{}_{\mathrm{i}}\text{=}\text{C}{}^{\mathrm{crystal}}{}_{\mathrm{i}}\text{C}{}^{\mathrm{liq}}{}_{\mathrm{i}}$) decrease linearly with increasing $\text{NBO}\text{T}$. The values of the individual REE crystal-liquid trace element partition coefficients have different functional relations to $\text{NBO}\text{T}$, so that the degree of light REE enrichment of the melts would depend on their $\text{NBO}\text{T}$.The solution mechanisms of minor oxides such as CO₂, H₂O, TiO₂, P₂O₅ and Fe₂O₃ in silicate melts are known. These data have been recast as changes of $\text{NBO}\text{T}$ of the melts with regard to the type of oxide and its concentration in the melt. From such data the dependence of crystal-liquid partition coefficients on concentration and type of minor oxide in melt solution has been calculated.     

17.

An electrochemical study of Ni²⁺, Co²⁺, and Zn²⁺ ions in melts of composition CaMgSi₂O₆     

Krystyna W. Semkow  Ronald A. Rizzo  Larry A. Haskin  David J. Lindstrom 《Geochimica et cosmochimica acta》1982,46(10):1879-1889

Cyclic voltammetry has been done for Ni²⁺, Co²⁺, and Zn²⁺ in melts of diopside composition in the temperature range 1425 to 1575°C. Voltammetric curves for all three ions excellently match theoretical curves for uncomplicated, reversible charge transfer at the Pt electrode. This implies that the neutral metal atoms remain dissolved in the melt. The reference electrode is a form of oxygen electrode. Relative to that reference assigned a reduction potential of 0.00 volt, the values of standard reduction potential for the ions are $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.32 ± .01}\text{V}$, $\text{E}{}^1\text{(}\text{Co}{}^{\mathrm{2+}}\text{Co}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.45 ± .02}\text{V}$, and $\text{E}{}^1\text{(}\text{Zn}{}^{\mathrm{2+}}\text{Zn}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.53 ± .01}\text{V}$. The electrode reactions are rapid, with first order rate constants of the order of 10^?2 cm/sec. Diffusion coefficients were found to be 2.6 × 10^?6 cm²/sec for Ni²⁺, 3.4 × 10^?6 cm²/sec for Co²⁺, and 3.8 × 10^?6 cm²/sec for Zn²⁺ at 1500°C. The value of $\text{E1}$ ($\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0$, diopside) is a linear function of temperature over the range studied, with values of ?0.35 V at 1425°C and ?0.29 V at 1575°C. At constant temperature the value of $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{, 1525°C}$) was not observed to vary with composition over the range CaO · MgO · 2SiO₂ to CaO·MgO·3SiO₂ or from 1.67 CaO·0.33MgO·2SiO₂ to 0.5 CaO·1.5MgO·2SiO₂. The value for the diffusion coefficient for Ni²⁺ decreased by an order of magnitude at 1525°C over the compositional range CaO · MgO · 1.25SiO₂ to CaO · MgO · 3SiO₂. This is consistent with a mechanism by which Ni²⁺ ions diffuse by moving from one octahedral coordination site to another in the melt, with the same Ni²⁺ species discharging at the cathode regardless of the SiO₂ concentration in the melt.     

18.

Diffusivity of oxygen in jadeite and diopside melts at high pressures     

N. Shimizu  I. Kushiro 《Geochimica et cosmochimica acta》1984,48(6):1295-1303

The diffusivity of oxygen was determined in melts of Jadeite (NaAlSi₂O₆) and diopside (CaMgSi₂O₆) compositions using diffusion couples with ¹⁸O as a tracer. In the Jadeite melt, the diffusivity of oxygen increases from $\text{6.87}{}_{\mathrm{?0.25}}{}^{\mathrm{+0.28}}\text{× 10}{}^{\mathrm{?10}}\text{cm}{}^2\text{/}\text{sec}$ at 5 Kb to $\text{1.32 ± 0.08 × 10}{}^{\mathrm{?9}}\text{cm}{}^2\text{/}\text{sec}$ at 20 Kb at constant temperature (1400°C), whereas in the diopside melt at 1650°C, the diffusivity decreases from $\text{7.30}{}_{\mathrm{?0.18}}{}^{0.29}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 10 Kb to $\text{5.28}{}_{\mathrm{?0.55}}{}^{\mathrm{+0.60}}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 17 Kb. These results demonstrate that the diffusivity is inversely correlated with the viscosity of the melt. For the jadeite melt, in particular, the inverse correlation is very well approximated by the Eyring equation using the diameter of oxygen ions as a unit distance of translation, suggesting that the viscous flow is rate-limited by the diffusion of individual oxygen ions. In the diopside melt, the activation volume is slightly greater than the molar volume of oxygen ion, indicating that the individual oxygen ion is the diffusion unit. The negative activation volume obtained for the jadeite melt is interpreted as the volume decrease associated with a diffusive jump of an oxygen ion due to local collapse of the network structure.     

19.

Prediction of Gibbs free energies of calcite-type carbonates and the equilibrium distribution of trace elements between carbonates and aqueous solutions     

Dimitri A Sverjensky 《Geochimica et cosmochimica acta》1984,48(5):1127-1134

A linear correlation exists between the standard Gibbs free energies of formation of calcite-type carbonates (MCO3) and the corresponding conventional standard Gibbs free energies of formation of the aqueous divalent cations (M2+) at 25 °C and 1 bar ΔGMCO30 = m(ΔGf,M2+0) ? 141,200 cal · mole?1 where m is equal to 0.9715. This relationship enables prediction of the standard free energies of formation of numerous hypothetical carbonates with the calcite structure. Associated uncertainties typically range from about ± 250 to 600 cal · mole?1. An important consequence of the above correlation is that the thermodynamic equilibrium constant for the distribution of two trace elements M and N between carbonate mineral and aqueous solution at 25 °C and 1 bar is proportional to the free energy difference between the corresponding two aqueous ions: In Combination of predicted standard free energies, entropies and volumes of carbonate minerals at 25°C and 1 bar with standard free energies of aqueous ions and the equation of state in Helgesonet al. (1981) enables prediction of the thermodynamic equilibrium constant for trace element distribution between carbonates and aqueous solutions at elevated temperatures and pressures. Interpretation of the thermodynamic equilibrium constant in terms of concentration ratios in the aqueous phase is considerably simplified if pairs of divalent trace elements are considered that have very similar ionic radii (e.g., $\text{Sr}{}^{\mathrm{2+}}\text{Pb}{}^{\mathrm{2+}}$, $\text{Mg}{}^{\mathrm{2+}}\text{Zn}{}^{\mathrm{2+}}$). In combination with data for the stabilities of complex ions in aqueous solutions, the above calculations enable useful limits to be placed on the concentrations of trace elements in hydrothermal solutions.     

20.

The solubility of noble gases in water and in NaCl brine     

S.P. Smith  B.M. Kennedy 《Geochimica et cosmochimica acta》1983,47(3):503-515

A direct-sampling, mass-spectrometric technique has been used to measure simultaneously the solubilities of He, Ne, Ar, Kr, and Xe in fresh water and NaCl brine (0 to 5.2 molar) from 0° to 65 °C, and at 1 atm total pressure of moist air. The argon solubility in the most concentrated brines is 4 to 7 times less than in fresh water at 65 °C and 0°C, respectively. The salt effect is parameterized using the Setschenow equation.

$\text{ln}\text{[}\text{β}{}^{\mathrm{i}}{}_{\mathrm{o}}\text{(T)}\text{β}{}^{\mathrm{i}}\text{(T)}\text{= MK}{}^{\mathrm{i}}{}_{\mathrm{M}}\text{(T)}$

where M is NaCl moiarity, βⁱ_o(T) and βⁱ(T) the Bunsen solubility coefficients for gas i in fresh water and brine, and Kⁱ_M(T) the empirical salting coefficient. Values of Kⁱ_M(T) are calculated using volumetric concentration units for noble gas and NaCl content and are independent of NaCl molarity. Below about 40°C, temperature coefficients of all Kⁱ_M are negative. The value of K^He_M is a minimum at 40°C. K^Ar_M decreases from about 0.40 at 0°C to 0.28 at 65 °C. The absolute magnitudes of the differences in salting coefficients (relative to K^Ar_M) decrease from 0° to 65°C. Over the range of conditions studied, all noble gases are salted out, and K^He_M ? K^Ne_M < K^Ar_M < K^Kr_M < K^Xe_M.From the solubility data, we calculated $\text{Δ}\text{G}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{H}{}^0{}_{\mathrm{tr}}$ and $\text{Δ}\text{C}{}^{\mathrm{O}}{}_{\mathrm{p,tr}}$ for the transfer of noble gases from fresh water to 1 molar NaCl solutions. At low temperatures $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, is positive, but decreases and becomes negative at temperatures ranging from about 25°C for He to 45°C for Xe. At low temperatures, the dissolved electrolyte apparently interferes with the formation of a cage of solvent molecules about the noble gas atom. At higher temperatures, the local environment of the gas atom in the brine appears to be slightly more ordered than in pure water, possibly reflecting the longer effective range of the ionic fields at higher temperature.The measured solubilities can be used to model noble gas partitioning in two-phase geothermal systems at low temperatures. The data can also be used to estimate the temperature and concentration dependence of the salt effect for other alkali halides. Extrapolation of the measured data is not possible due to the incompletely-characterized minima in the temperature dependence of the salting coefficients. The regularities in the data observed at low temperatures suggest relatively few high-temperature data will be required to model the behavior of noble gases in high-temperature geothermal brines.     

(°C)T	$\text{K}{}_{\mathrm{D}}{}^{\mathrm{Phlog/Vap.}}$	$\text{X}{}_{\mathrm{RM}}$	$\text{K}{}_{\mathrm{D}}{}^{\mathrm{Sanid/Vap.}}$	$\text{X}{}_{\mathrm{RM}}$
500	$\text{0.64 ± 0.11}$	0–0.2	$\text{0.17 ± 0.04}$	0–0.07
700	$\text{1.11 ± 0.11}$	0–0.2	$\text{0.33 ± 0.04}$	0–0.1
800	$\text{1.28 ± 0.03}$	0–0.2	$\text{0.45 ± 0.06}$	0–0.1

Trace element partitioning and melt structure: An experimental study at 1 atm pressure

Diopside-melt and forsterite-melt rare earth (REE) and Ni partition coefficients have been determined as a function of bulk compositions of the melt. Available Raman spectroscopic data have been used to determine the structures of the melts coexisting with diopside and forsterite. The compositional dependence of the partition coefficients is then related to the structural changes of the melt.The melts in all experiments have a ratio of nonbridging oxygens to tetrahedral cations ($\text{NBO}\text{T}$) between 1 and 0. The quenched melts consist of structural units that have, on the average, 2 (chain), 1 (sheet) and 0 (three-dimensional network) nonbridging oxygens per tetrahedral cation. The proportions of these structural units in the melts, as well as the overall $\text{NBO}\text{T}$, change as a function of the bulk composition of the melt.It has been found that Ce, Sm, Tm and Ni crystal-liquid partition coefficients ($\text{K}{}^{\mathrm{crystal?liq}}{}_{\mathrm{i}}\text{=}\text{C}{}^{\mathrm{crystal}}{}_{\mathrm{i}}\text{C}{}^{\mathrm{liq}}{}_{\mathrm{i}}$) decrease linearly with increasing $\text{NBO}\text{T}$. The values of the individual REE crystal-liquid trace element partition coefficients have different functional relations to $\text{NBO}\text{T}$, so that the degree of light REE enrichment of the melts would depend on their $\text{NBO}\text{T}$.The solution mechanisms of minor oxides such as CO₂, H₂O, TiO₂, P₂O₅ and Fe₂O₃ in silicate melts are known. These data have been recast as changes of $\text{NBO}\text{T}$ of the melts with regard to the type of oxide and its concentration in the melt. From such data the dependence of crystal-liquid partition coefficients on concentration and type of minor oxide in melt solution has been calculated.     

An electrochemical study of Ni2+, Co2+, and Zn2+ ions in melts of composition CaMgSi2O6

Cyclic voltammetry has been done for Ni²⁺, Co²⁺, and Zn²⁺ in melts of diopside composition in the temperature range 1425 to 1575°C. Voltammetric curves for all three ions excellently match theoretical curves for uncomplicated, reversible charge transfer at the Pt electrode. This implies that the neutral metal atoms remain dissolved in the melt. The reference electrode is a form of oxygen electrode. Relative to that reference assigned a reduction potential of 0.00 volt, the values of standard reduction potential for the ions are $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.32 ± .01}\text{V}$, $\text{E}{}^1\text{(}\text{Co}{}^{\mathrm{2+}}\text{Co}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.45 ± .02}\text{V}$, and $\text{E}{}^1\text{(}\text{Zn}{}^{\mathrm{2+}}\text{Zn}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.53 ± .01}\text{V}$. The electrode reactions are rapid, with first order rate constants of the order of 10^?2 cm/sec. Diffusion coefficients were found to be 2.6 × 10^?6 cm²/sec for Ni²⁺, 3.4 × 10^?6 cm²/sec for Co²⁺, and 3.8 × 10^?6 cm²/sec for Zn²⁺ at 1500°C. The value of $\text{E1}$ ($\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0$, diopside) is a linear function of temperature over the range studied, with values of ?0.35 V at 1425°C and ?0.29 V at 1575°C. At constant temperature the value of $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{, 1525°C}$) was not observed to vary with composition over the range CaO · MgO · 2SiO₂ to CaO·MgO·3SiO₂ or from 1.67 CaO·0.33MgO·2SiO₂ to 0.5 CaO·1.5MgO·2SiO₂. The value for the diffusion coefficient for Ni²⁺ decreased by an order of magnitude at 1525°C over the compositional range CaO · MgO · 1.25SiO₂ to CaO · MgO · 3SiO₂. This is consistent with a mechanism by which Ni²⁺ ions diffuse by moving from one octahedral coordination site to another in the melt, with the same Ni²⁺ species discharging at the cathode regardless of the SiO₂ concentration in the melt.     

Diffusivity of oxygen in jadeite and diopside melts at high pressures

The diffusivity of oxygen was determined in melts of Jadeite (NaAlSi₂O₆) and diopside (CaMgSi₂O₆) compositions using diffusion couples with ¹⁸O as a tracer. In the Jadeite melt, the diffusivity of oxygen increases from $\text{6.87}{}_{\mathrm{?0.25}}{}^{\mathrm{+0.28}}\text{× 10}{}^{\mathrm{?10}}\text{cm}{}^2\text{/}\text{sec}$ at 5 Kb to $\text{1.32 ± 0.08 × 10}{}^{\mathrm{?9}}\text{cm}{}^2\text{/}\text{sec}$ at 20 Kb at constant temperature (1400°C), whereas in the diopside melt at 1650°C, the diffusivity decreases from $\text{7.30}{}_{\mathrm{?0.18}}{}^{0.29}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 10 Kb to $\text{5.28}{}_{\mathrm{?0.55}}{}^{\mathrm{+0.60}}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 17 Kb. These results demonstrate that the diffusivity is inversely correlated with the viscosity of the melt. For the jadeite melt, in particular, the inverse correlation is very well approximated by the Eyring equation using the diameter of oxygen ions as a unit distance of translation, suggesting that the viscous flow is rate-limited by the diffusion of individual oxygen ions. In the diopside melt, the activation volume is slightly greater than the molar volume of oxygen ion, indicating that the individual oxygen ion is the diffusion unit. The negative activation volume obtained for the jadeite melt is interpreted as the volume decrease associated with a diffusive jump of an oxygen ion due to local collapse of the network structure.     

Prediction of Gibbs free energies of calcite-type carbonates and the equilibrium distribution of trace elements between carbonates and aqueous solutions

A linear correlation exists between the standard Gibbs free energies of formation of calcite-type carbonates (MCO3) and the corresponding conventional standard Gibbs free energies of formation of the aqueous divalent cations (M2+) at 25 °C and 1 bar ΔGMCO30 = m(ΔGf,M2+0) ? 141,200 cal · mole?1 where m is equal to 0.9715. This relationship enables prediction of the standard free energies of formation of numerous hypothetical carbonates with the calcite structure. Associated uncertainties typically range from about ± 250 to 600 cal · mole?1. An important consequence of the above correlation is that the thermodynamic equilibrium constant for the distribution of two trace elements M and N between carbonate mineral and aqueous solution at 25 °C and 1 bar is proportional to the free energy difference between the corresponding two aqueous ions: In Combination of predicted standard free energies, entropies and volumes of carbonate minerals at 25°C and 1 bar with standard free energies of aqueous ions and the equation of state in Helgesonet al. (1981) enables prediction of the thermodynamic equilibrium constant for trace element distribution between carbonates and aqueous solutions at elevated temperatures and pressures. Interpretation of the thermodynamic equilibrium constant in terms of concentration ratios in the aqueous phase is considerably simplified if pairs of divalent trace elements are considered that have very similar ionic radii (e.g., $\text{Sr}{}^{\mathrm{2+}}\text{Pb}{}^{\mathrm{2+}}$, $\text{Mg}{}^{\mathrm{2+}}\text{Zn}{}^{\mathrm{2+}}$). In combination with data for the stabilities of complex ions in aqueous solutions, the above calculations enable useful limits to be placed on the concentrations of trace elements in hydrothermal solutions.     

The solubility of noble gases in water and in NaCl brine

A direct-sampling, mass-spectrometric technique has been used to measure simultaneously the solubilities of He, Ne, Ar, Kr, and Xe in fresh water and NaCl brine (0 to 5.2 molar) from 0° to 65 °C, and at 1 atm total pressure of moist air. The argon solubility in the most concentrated brines is 4 to 7 times less than in fresh water at 65 °C and 0°C, respectively. The salt effect is parameterized using the Setschenow equation.

$\text{ln}\text{[}\text{β}{}^{\mathrm{i}}{}_{\mathrm{o}}\text{(T)}\text{β}{}^{\mathrm{i}}\text{(T)}\text{= MK}{}^{\mathrm{i}}{}_{\mathrm{M}}\text{(T)}$

where M is NaCl moiarity, βⁱ_o(T) and βⁱ(T) the Bunsen solubility coefficients for gas i in fresh water and brine, and Kⁱ_M(T) the empirical salting coefficient. Values of Kⁱ_M(T) are calculated using volumetric concentration units for noble gas and NaCl content and are independent of NaCl molarity. Below about 40°C, temperature coefficients of all Kⁱ_M are negative. The value of K^He_M is a minimum at 40°C. K^Ar_M decreases from about 0.40 at 0°C to 0.28 at 65 °C. The absolute magnitudes of the differences in salting coefficients (relative to K^Ar_M) decrease from 0° to 65°C. Over the range of conditions studied, all noble gases are salted out, and K^He_M ? K^Ne_M < K^Ar_M < K^Kr_M < K^Xe_M.From the solubility data, we calculated $\text{Δ}\text{G}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{H}{}^0{}_{\mathrm{tr}}$ and $\text{Δ}\text{C}{}^{\mathrm{O}}{}_{\mathrm{p,tr}}$ for the transfer of noble gases from fresh water to 1 molar NaCl solutions. At low temperatures $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, is positive, but decreases and becomes negative at temperatures ranging from about 25°C for He to 45°C for Xe. At low temperatures, the dissolved electrolyte apparently interferes with the formation of a cage of solvent molecules about the noble gas atom. At higher temperatures, the local environment of the gas atom in the brine appears to be slightly more ordered than in pure water, possibly reflecting the longer effective range of the ionic fields at higher temperature.The measured solubilities can be used to model noble gas partitioning in two-phase geothermal systems at low temperatures. The data can also be used to estimate the temperature and concentration dependence of the salt effect for other alkali halides. Extrapolation of the measured data is not possible due to the incompletely-characterized minima in the temperature dependence of the salting coefficients. The regularities in the data observed at low temperatures suggest relatively few high-temperature data will be required to model the behavior of noble gases in high-temperature geothermal brines.