Kinetic and thermodynamic factors controlling the distribution of SO32− and Na+ in calcites and selected aragonites

Significant amounts of SO₄^2?, Na⁺, and OH^? are incorporated in marine biogenic calcites. Biogenic high Mg-calcites average about 1 mole percent SO₄^2?. Aragonites and most biogenic low Mg-calcites contain significant amounts of Na⁺, but very low concentrations of SO₄^2?. The SO₄^2? content of non-biogenic calcites and aragonites investigated was below 100 ppm. The presence of Na⁺ and SO₄^2? increases the unit cell size of calcites. The solid-solutions show a solubility minimum at about 0.5 mole percent SO₄^2? beyond which the solubility rapidly increases. The solubility product of calcites containing 3 mole percent SO₄^2? is the same as that of aragonite. Na⁺ appears to have very little effect on the solubility product of calcites. The amounts of Na⁺ and SO₄^2? incorporated in calcites vary as a function of the rate of crystal growth. The variation of the distribution coefficient ($\text{D}$) of SO₄^2? in calcite at 25.0°C and 0.50 molal NaCl is described by the equation $\text{D = k}{}_0\text{+ k}{}_1\text{R}$ where $\text{k}{}_0$ and $\text{k}{}_1$ are constants equal to $\text{6.16 × 10}{}^{\mathrm{?6}}$ and $\text{3.941 × 10}{}^{\mathrm{?6}}$, respectively, and $\text{R}$ is the rate of crystal growth of calcite in mg·min^?1·g^?1 of seed. The data on Na⁺ are consistent with the hypothesis that a significant amount of Na⁺ occupies interstitial positions in the calcite structure. The distribution of Na⁺ follows a Freundlich isotherm and not the Berthelot-Nernst distribution law. The numerical value of the Na⁺ distribution coefficient in calcite is probably dependent on the number of defects in the calcite structure. The Na⁺ contents of calcites are not very accurate indicators of environmental salinities.     

Complex formation in the ternary system Mg(II)-CO2-H2O

The carbonato and hydrogencarbonato complexes of Mg²⁺ were investigated at 25 and 50° in solutions of the constant ClO₄^? molality (3 M) consisting preponderantly of NaClO₄. The experimental data could be explained assuming the following equilibria: Mg²⁺ + CO_2B + H₂O ag MgHCO⁺₃ + H⁺, $\text{log}{}^1\text{β}{}_1\text{= ?7.64}{}_4\text{± 0.01}{}_7\text{(25°), ?7.46}{}_2\text{± 0.01}{}_1\text{(50°)}$, Mg²⁺ + 2 CO_2g + 2 H₂Oag Mg(HCO₃)⁰₂ ± 2 H⁺, $\text{log}{}^1\text{β}{}_2\text{= ?15.00 ± 0.14 (25°)}$, ?15.37 ± 0.39 (50°), Mg²⁺ + CO_2g + H₂Oag MgCO⁰₃ + 2 H⁺, $\text{log}{}^1\text{k}{}_1\text{= ?15.64 ± 0.06 (25°)}$,?15.23 ± 0.02 (50°), with the assumption γ_MgCO3⁰ = γ_Mg(HCO3)⁰2, ΔG⁰(I = 0) for the reaction MgCO⁰₃ + CO_2g + H₂O = Mg(HCO₃)⁰₂ was estimated to be ?3.9₁ ± 0.8₆ and 0.6 ± 2.4 kJ/mol at 25 and 50°C, respectively. The abundance of carbonate linked Mg(II) species in fresh water systems is discussed.     

The solubilities of calcite,aragonite and vaterite in CO2-H2O solutions between 0 and 90°C,and an evaluation of the aqueous model for the system CaCO3-CO2-H2O

Calculations based on approximately 350 new measurements (Ca_T-PCO₂) of the solubilities of calcite, aragonite and vaterite in CO₂-H₂O solutions between 0 and 90°C indicate the following values for the log of the equilibrium constants K_C, K_A, and K_V respectively, for the reaction CaCO₃(s) = Ca²⁺ + CO^2?₃: $\text{Log K}{}_{\mathrm{C}}\text{= ?171.9065 ? 0.077993T +}\text{2839.319}\text{T}\text{+ 71.595 log T}$$\text{Log K}{}_{\mathrm{A}}\text{= ?171.9773 ? 0.077993T +}\text{2903.293}\text{T}\text{+71.595 log T}$$\text{Log K}{}_{\mathrm{V}}\text{= ?172.1295 ? 0.077993T +}\text{3074.688}\text{T}\text{+ 71.595 log T}$ where T is in ^oK. At 25°C the logarithms of the equilibrium constants are ?8.480 ± 0.020, ?8.336 ± 0.020 and ?7.913 ± 0.020 for calcite, aragonite and vaterite, respectively.The equilibrium constants are internally consistent with an aqueous model that includes the CaHCO⁺₃ and CaCO⁰₃ ion pairs, revised analytical expressions for CO₂-H₂O equilibria, and extended Debye-Hückel individual ion activity coefficients. Using this aqueous model, the equilibrium constant of aragonite shows no PCO₂-dependence if the CaHCO⁺₃ association constant is between 0 and 90°C, corresponding to the value logK_Cahco⁺3 = 1.11 ± 0.07 at 25°C. The CaCO⁰₃ association constant was measured potentiometrically to be between 5 and 80°C, yielding logK_CaCO⁰3 = 3.22 ± 0.14 at 25°C.The CO₂-H₂O equilibria have been critically evaluated and new empirical expressions for the temperature dependence of K_H, K₁ and K₂ are $\text{log K}{}_{\mathrm{H}}\text{= 108.3865 + 0.01985076T ?}\text{6919.53}\text{T}\text{? 40.45154 log T +}\text{669365.}\text{T}{}^2$, $\text{log K}{}_1\text{= ?356.3094 ? 0.06091964T +}\text{21834.37}\text{T}\text{+ 126.8339 log T —}\text{1684915.}\text{T}{}^2$ and logK₂ = ?107.8871 ? 0.03252849T + 5151.79/T + 38.92561 logT ? 563713.9/T² which may be used to at least 250°C. These expressions hold for 1 atm. total pressure between 0 and 100°C and follow the vapor pressure curve of water at higher temperatures.Extensive measurements of the pH of Ca-HCO₃ solutions at 25°C and 0.956 atm PCO₂ using different compositions of the reference electrode filling solution show that measured differences in pH are closely approximated by differences in liquid-junction potential as calculated by the Henderson equation. Liquid-junction corrected pH measurements agree with the calculated pH within 0.003-0.011 pH.Earlier arguments suggesting that the CaHCO⁺₃ ion pair should not be included in the CaCO₃-CO₂-H₂O aqueous model were based on less accurate calcite solubility data. The CaHCO⁺₃ ion pair must be included in the aqueous model to account for the observed PCO₂-dependence of aragonite solubility between 317 ppm CO₂ and 100% CO₂.Previous literature on the solubility of CaCO₃ polymorphs have been critically evaluated using the aqueous model and the results are compared.     

The incorporation of Mg2+ and Sr2+ into calcite overgrowths: influences of growth rate and solution composition

The concentrations of Mg²⁺ and Sr²⁺ incorporated within calcite overgrowths precipitated from seawater and related solutions, determined at 25°C, were independent of the precipitation rate over approximately an order of magnitude. The saturation states used to produce this range of precipitation rates varied from 3 to 17 depending on the composition of the solution.The amount of Mg²⁺ incorporated in the overgrowths was not directly proportional to $\text{Mg}{}^{\mathrm{2+}}\text{Ca}{}^{\mathrm{2+}}$ in solution over the entire range (1–20) of ratios studied. Below a ratio of 7.5, the overgrowth was enriched in MgCO₃ relative to what is predicted by the constant distribution coefficient measured above a ratio of 7.5. This increased MgCO₃ correlates with the relative enrichment of adsorbed Mg²⁺. Above a ratio of 7.5 the concentration of MgCO₃ in the calcite overgrowths followed a classical thermodynamic behavior characterized by a constant distribution coefficient of 0.0123 (±0.008 std dev).The concentration of SrCO₃ incorporated in the overgrowths was linearly related to the MgCO₃ content of the overgrowths, and is attributed to increased solubility of SrCO₃ in calcite due to the incorporation of the smaller Mg²⁺ ions.The kinetic data indicate that the growth mechanism involves the adsorption of the cations on the surface of the calcite prior to dehydration and final incorporation. It is suggested that dehydration of cations at the surface is the rate controlling step.     

Mineral-solution equilibria—III. The system Na2OAl2O3SiO2H2OHCl

Chemical equilibrium between sodium-aluminum silicate minerals and chloride bearing fluid has been experimentally determined in the range 500–700°C at 1 kbar, using rapid-quench hydrothermal methods and two modifications of the Ag + AgCl acid buffer technique. The temperature dependence of the thermodynamic equilibrium constant (K) for the reaction $\text{NaAlSi}{}_3\text{O}{}_8\text{+ HCl}{}^{\mathrm{o}}\text{= NaCl}{}^{\mathrm{o}}\text{+}\text{1}\text{2Al}{}_2\text{SiO}{}_5\text{, +}\text{5}\text{2SiO}{}_2\text{+}\text{1}\text{2H}{}_2\text{O}$ Albite Andalusite Qtz. can be described by the following equation: log k = ?4.437 + 5205.6/T(K) The data from this study are consistent with experimental results reported by Montoya and Hemley (1975) for lower temperature equilibria defined by the assemblages albite + paragonite + quartz + fluid and paragonite + andalusite + quartz + fluid. Values of the equilibrium constants for the above reactions were used to estimate the difference in Gibbs free energy of formation between NaCl^o and HCl^o in the range 400–700°C and 1–2 kbar. Similar calculations using data from phase equilibrium studies reported in the literature were made to determine the difference in Gibbs free energy of formation between KCl^o and HCl^o. These data permit modelling of the chemical interaction between muscovite + kspar + paragonite + albite + quartz assemblages and chloride-bearing hydrothermal fluids.     

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.     

Pyroxene geothermometry and geobarometry: experimental and thermodynamic evaluation of some subsolidus phase relations involving pyroxenes in the system CaO-MgO-Al2O3-SiO2

One petrogenetic grid for plagioclase-, spinel- and garnet-lherzolite analogues in the system CaO-MgO-Al₂O₃-SiO₂ is presented from 1 bar to 30 kbar and 400 to 1500°C. Another grid for olivine-gabbro, spinel-gabbro and garnet-pyroxenite analogues in the same system is presented from 1 bar to 25 kbar and 500 to 1500°C. Both grids show the distribution of the mineral assemblages and the variations in the composition of clinopyroxene with temperature and pressure. They were developed by applying simple thermodynamic mixing models of clinopyroxene to experimentally determined clino-pyroxene compositions.Calcium tschermak's pyroxene (CaAl₂SiO₆) in complex CaMgSi₂O₆-CaAl₂SiO₆-Mg₂Si₂O₆ clinopyroxenes is best represented by a local charge balance mixing model where Enthalpy and entropy changes of subsolidus reactions involving variations in the CaAl₂SiO₆ and Mg₂Si₂O₆ content of clinopyroxene are interdependent due to nonideal mixing of these two end-members. CaAl₂SiO₆ can strongly reduce the mutual solubility of clinopyroxene and orthopyroxene at moderate pressures and high temperatures. Failure to take this into account can result in temperature underestimates (up to 150°C) of spinel-lherzolites, garnet-pyroxenites, low pressure garnet-lherzolites, spinel-gabbros, and high pressure plagioclase-lherzolites and olivine-gabbros. However, at temperatures and pressures where the Al₂O₃ content of clinopyroxene is low (e.g. garnet-lherzolite nodules in kimberlite), the mutual solubility is adequantely represented by experimental results in the system CaO-MgO-SiO₂.     

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.     

Mineral-solution equilibria—IV. Solubilities and the thermodynamic properties of FeCl20 in the system Fe2O3-H2-H2O-HCl

The solubility of hematite in chloride-bearing hydrothermal fluids was determined in the temperature range 400–600°C and at 1000 and 2000 bars using double-capsule, rapid-quench hydrothermal techniques and a modification of the Ag + AgCl buffer method (Frantz and Popp, 1979). The changes in the molalities of associated hydrogen chloride () as a function of the molality of total iron in the fluid at constant temperature and pressure were used to identify the predominant species of iron in the hydrothermal fluid. The molality of associated HCl varied from 0.01 to 0.15. Associated FeCl₂⁰ was found to be the most abundant species in equilibrium with hematite. Determination of Cl/Fe in the fluid in equilibrium with hematite yields values approximately equal to 2.0 suggesting that ferrous iron is the dominant oxidation state.The equilibrium constant for the reaction Fe₂O₃ + 4HCl⁰ + H₂ = 2FeCl₂⁰ + 3H₂O was calculated and used to estimate the difference in Gibbs free energy between FeCl₂⁰ and HCl⁰ in the temperature range 400–600°C at 1000 and 2000 bars pressure.     

40Ar39Ar dating,Ar diffusion properties,and cooling rate determinations of severely shocked chondrites

Determinations of ${}^{40}\text{Ar}{}^{39}\text{Ar}$ ages are reported for seven severely shock-heated chondrites. Shaw gives a plateau age of 4.29 Gyr. Louisville, Farmington, and Wickenburg give well-defined intercept ages of 0.5–0.6 Gyr. Orvinio, Arapahoe, and Lubbock show complex ${}^{40}\text{Ar}{}^{39}\text{Ar}$ release curves, with age minima of 0.7–1.0 Gyr. Degassing times of 0.5–1.0 Gyr are suggested for these meteorites. Most severely shocked chondrites were apparently not totally degassed of ⁴⁰Ar by the event, but retained from ~ 2 to ~45% of their ⁴⁰Ar. When calculated values of the diffusion parameter, $\text{D}\text{a}{}^2$, for Ar are examined in Arrhenius plots, they show two distinct linear relationships, which apparently correspond to the degassing of different mineral phases with distinct $\text{K}\text{Ca}$ ratios and different average temperatures for Ar release. The experimentally determined values of $\text{D}\text{a}{}^2$ for the high temperature phase of several severely shocked chondrites are ~10^?7 to 10^?5sec^?1 for their determined shock-heating temperatures of ~950°C to ~ 1200°C. The inferred reheating temperatures, $\text{D}\text{a}{}^2$ values, and fraction of ⁴⁰Ar loss during the reheating event for these seven chondrites suggest post-shock cooling rates and burial depth of ~ 10^?2 10^?4°C/sec and ~0.5–2m, respectively. For three chondrites these cooling rates agree with those determined from Ni diffusion in metal grains: for five chondrites the cooling rates derived from ⁴⁰Ar and Ni disagree by a factor of ~10⁵. It is suggested that five of these severely shocked chondrites were part of large ejecta blankets containing hot material and cold clasts with a distribution of sizes and that the cooling rate of this ejecta appreciably decreased as a function of time.     

Partitioning of Ni2+, Co2+, Fe2+, Mn2+ and Mg2+ between olivine and silicate melts: compositional dependence of partition coefficient

Partitioning of Ni²⁺, Co²⁺, Fe²⁺, Mn²⁺ and Mg²⁺ between olivine and silicate melts has been determined near the join (Mg_0.5·-Fe_0.5)₂SiO₄-K₂O 4SiO₂ and for seven different basaltic compositions. The experiments were made at 1 atm total pressure, 1500-1100°C, and under moderate to reducing oxygen fugacities. The concentration factor, defined as K^MO = (MO)^ol/(MO)^liq (molar ratio), increases markedly for all the cations studied as the olivine component of the liquid decreases. Much of the increase in K^MO is considered as due to the compositional effect of the coexisting liquid: the temperature effect on K^MO is probably opposite to the compositional effect (K^MO decreases as temperature decreases).The partition coefficient K^MO-MgO = (MO/MgO)^ol/(MO/MgO)^liq for the reaction, M_ol²⁺ + Mg_liq²⁺ = M_liq²⁺ + Mg_ol²⁺. is relatively constant over a wide range of SiO₂ content of the liquid, except in the case of Ni²⁺. The partition coefficients have similar ranges both in synthetic and natural rock systems: K^NiO-MgO = 1.8–3.0, K^CoO-MgO = 0.6–0.8, K^FeO-MgO = 0.27–0.38, and K^MnO-MgO = 0.23–0.32. There is a systematic variation in the partition coefficient K^MO-MgO with the composition of liquid; K^MO-MgO increases with increasing SiO₂ content of melt. The partition coefficients for Co²⁺, Fe²⁺ and Mn²⁺ are useful to test the equilibration of olivine with magma of a wide compositional range.     

Propriétés thermodynamiques des solutions (Na,K)ClH2O entre 400 et 800°C, 1–2 Kb: une revue des équilibres d'échange avec les silicates sodipotassiques

Several independant determinations of the difference in Gibbs free energy of formation (from the elements at 25°, 1 bar) between NaCl⁰ and KCl⁰ in aqueous solutions (molality > 0.5) are derived from equilibrium data between alkali feldspars, feldspathoids (nepheline-kalsilite), micas (muscoviteparagonite) and hydrothermal (Na, K)Cl-H₂O solutions. These results along with other data from the literature are discussed. The relation: $\text{Δ}{}^0\text{?,}\text{KCl}{}^0\text{? Δ}{}^0\text{/tf,}\text{NaCl}{}^0\text{(}\text{J}\text{) = ?16500(± 2500) ? 18(± 4) T(}\text{K}\text{)}$. is proposed from 400 to 800°C and 1 to 2 Kbar.     

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.     

Mineral-solution equilibria—I. An experimental study of complexing and thermodynamic properties of aqueous MgCl2 in the system MgO-SiO2-H2O-HCl

Speciation of aqueous magnesium in the system MgO-SiO₂-H₂O-HCl in supercritical aqueous fluids has been investigated using standard rapid-quench hydrothermal techniques and a modification of the Ag + AgCl buffer method (Frantz and Eugster, 1973. Am. J. Sci.267, 268–286). A concentric double-capsule charge was utilized. The outer gold capsule contained the assemblage talc + quartz + Ag + AgCl + H₂O-MgCl₂ fluid; the inner platinum capsule, Ag + AgCl + H₂O-HCl fluid. During the experiments, and thus $\text{?}{}_{\mathrm{HCl}}$ equilibrated between the two capsules. After quenching, measurement of the chloride concentration in the fluid in the inner capsule and total magnesium in the fluid in the outer capsule defines the concentrations of HCl and Mg that coexist with talc + quartz in the outer capsule. Changes in the measured molality of HCl as a function of the total magnesium concentration at constant P and T were used to identify the predominant species of magnesium in the hydrothermal fluid. Experimental results showed that at 2000 bar, MgCl°₂ is the predominant species above 550°C and Mg²⁺, below 400°C. Data at intermediate temperatures when combined with the dissociation constant for HCl were used to obtain the dissociation constant for MgCl°₂. The results of these experiments were combined with results from experiments using Ag + AgCl in conjunction with the oxygen buffer, hematite-magnetite, to obtain the equilibrium constant for the reaction $\text{1}\text{3}\text{Talc + 2HC1° H}{}_2\text{O MgCl°}{}_2\text{+}\text{4}\text{3}\text{Quartz +}\text{4}\text{3}\text{H}{}_2\text{O}$ from which the difference in Gibbs free energy of MgCl°₂ and HC1° was obtained as a function of temperature at 1000, 1500 and 2000 bar pressure, Solubility constants for brucite. forsterite, chrysotile, and talc were calculated.     

Determination of the 87Rb decay constant

The decay constant ⁸⁷Rb has been redetermined by measuring the amount of radiogenic ⁸⁷Sr produced over a period of 19 years, in 20 g samples of purified RbClO₄, using isotope dilution techniques. The rubidium sample was spiked with ⁸⁴Sr and the nanogram quantities of strontium separated by coprecipitation with Ba(NO₃)₂. Analyses were carried out on a 25cm, 90° sector mass spectrometer equipped with a Spiraltron electron multiplier. Measurement of three independent ratios permitted continuous monitoring of the ion beam fractionation. The average of nine determinations gives a value for the decay constant of $\text{1.419(±0.012) × 10}{}^{\mathrm{?11}}\text{yr}{}^{\mathrm{?1}}\text{(2σ)}$. [$\text{τ}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{= 4.89(±0.04) × 10}{}^{10}\text{yr}$.]     

Alkali-free tourmaline in the system MgO-Al2O3-B2O3-SiO2-H2O

An end member of the tourmaline series with a structural formula □(Mg₂Al)Al₆(BO₃)₃[Si₆O₁₈](OH)₄ has been synthesized in the system MgO-Al₂O₃-B₂O₃-SiO₂-H₂O where it represents the only phase with a tourmaline structure. Our experiments provide no evidence for the substitutions Al → Mg + H, Mg → 2H, B + H → Si, and AlAl → MgSi and we were not able to synthesize a phase “Mg-aluminobuergerite” characterized by Mg in the (3a)-site and a strong (OH)-deficiency reported by Rosenberg and Foit (1975). The alkali-free tourmaline has a vacant (3a)-site and is related to dravite by the □ + Al for Na + Mg substitution. It is stable from at least 300°C to about 800°C at low fluid pressures and 100% excess B₂O₃, and can be synthesized up to a pressure of 20 kbars. At higher temperatures the tourmaline decomposes into grandidierite or a boron-bearing phase possibly related to mullite (“B-mullite”), quartz, and unidentified solid phases, or the tourmaline melts incongruently into corundum + liquid, depending on pressure. In the absence of excess B₂O₃ tourmaline stability is lowered by about 60°C. Tourmaline may coexist with the other MgO-Al₂O₃-B₂O₃-SiO₂-H₂O phases forsterite, enstatite, chlorite, talc, quartz, grandidierite, corundum, spinel, “B-mullite,” cordierite, and sinhalite depending on the prevailing PTX-conditions.The (3a)-vacant tourmaline has the space group R3m with $\text{a =15.90}\text{A}\text{?}$, $\text{c = 7.115}\text{A}\text{?}$, and $\text{V = 1557.0}\text{A}\text{?}{}^3$. However, these values vary at room temperature with the pressure-temperature conditions of synthesis by $\text{±0.015}\text{A}\text{?}\text{in a}$, $\text{±0.010}\text{A}\text{?}\text{in c}$, and $\text{±4.0}\text{A}\text{?}{}^3\text{in V}$, probably as a result of MgAl order/disorder relations in the octahedral positions. Despite these variations intensity calculations support the assumed structural formula. Refractive indices are n_o = 1.631(2), n_E = 1.610(2), Δn = 0.021. The infrared spectrum is intermediate between those of dravite and elbaite. The common alkali and calcium deficiencies of natural tourmalines may at least partly be explained by miscibilities towards (3a)-vacant end members. The apparent absence of (3a)-vacant tourmaline in nature is probably due to the lack of fluids that carry boron but no Na or Ca.     

Variability of the He3 and Ne21 production rates in ordinary chondrites

Al²⁶ and noble gas contents of 6 ordinary chondrites with $\text{He}{}^3\text{Ne}{}^{21}$ ratios above 6.0 or below 4.0 are used to infer the variability of the production rates of He³ and Ne²¹ (P_He³ and P_Ne²¹). The ratio of the observed Al²⁶ content to a calculated, normal value is taken as a measure of the change of P_Ne²¹ from its normal value. The corresponding change in P_He³ is then computed from the observed $\text{He}{}^3\text{Ne}{}^{21}$ ratio and an average value of P_He³.According to these calculations which exclude orbital effects, P_He³ will be near the average value in meteorites with high $\text{He}{}^3\text{Ne}{}^{21}$ ratios, while P_Ne²¹ will be about 30 per cent below normal. In meteorites with low $\text{He}{}^3\text{Ne}{}^{21}$ ratios, P_He³ may be depressed by as much as 25 per cent from normal while P_Ne²¹ may be 15–20 per cent above the average.     

Thermodynamics of brine-salt equilibria—I. The systems NaCl-KCl-MgCl2-CaCl2-H2O and NaCl-MgSO4-H2O at 25°C

A thermodynamic model for concentrated brines has been developed which is capable of predicting the solubilities of many of the common evaporite minerals in chloro-sulfate brines at 25°C and 1 atm. The model assumes that the behaviour of the mean stoichiometric ionic activity coefficient in mixtures of aqueous electrolytes can be described by the Scatchard deviation function and Harned's Rule. In solutions consisting of one salt and H₂O, the activity coefficient is described by the expression where $\text{a}\text{?}$ and B? salt specific parameters obtained from data regression. In a mixture of n electrolytes and H₂O, B? for the ith component is given by $\text{B}{}^{\mathrm{<mtext>i</mtext><mtext>?</mtext>}}{}_{\mathrm{i}}\text{=B}\text{i}\text{?}{}_{\mathrm{i}}\text{+σ α}{}_{\mathrm{ij}}\text{y}{}_{\mathrm{j}}$ where α_ij is a (constant) mixing parameter characterizing the interaction of the i and j components and y_j is the ionic strength fraction of the jth component. The activity of H₂O is obtained from a Gibbs-Duhem integration and does not require any additional parameters or assumptions. In this study, parameters have been obtained for the systems NaCl-KCl-MgCl₂-CaCl₂-H₂O and NaCl-MgSO₄-H₂O at 25°C and 1 atm. Computed solubility curves and solution compositions predicted for invariant points in these systems agree well with the experimental data. The model is flexible and easily extended to other systems and to higher temperatures.     

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°.     

Solubilities of some hydrous REE phosphates with implications for diagenesis and sea water concentrations

Solubility product determinations suggest that the hydrous phosphates of the rare earths, REPO₄ · $\text{xH}{}_2\text{O}$, are important in controlling the sea water REE concentrations. Two of these solids, rhabdophane, (P6₂22) and “hydrous xenotime”, (14₁/amd), have been synthesized at 100°C via the acid hydrolysis of the respective REE pyrophosphate. The solubility products at infinite dilution were determined to be $\text{pK}{}^0\text{= 24.5}$, (La at 25°C); 26.0, (Pr at 100°C); 25.7, (Nd at 100°C): and 25.5, (Er at 100°C). On the basis of calculations involving the reaction of RE³⁺ with apatite to form the hydrous phosphate, the lanthanum concentration in sea water is predicted to be about 140 pmol/L. Laboratory experiments support the hypothesis that apatite is a substrate for reactions with dissolved REE.