The effect of pressure on the solubility of minerals in water and seawater

The effect of presure on the solubility of minerals in water and seawater can be estimated from In

$\text{(}\text{K}{}^{\mathrm{P}}{}_{\mathrm{sp}}\text{K}{}^0{}_{\mathrm{sp}}\text{) +}\text{(?ΔVP + 0.5ΔKP}{}^2\text{)}\text{RT}$

where the volume (ΔV) and compressibility (ΔK) changes at atmospheric pressure (P = 0) are given by

$\text{ΔV =}\text{V}\text{?}\text{(M}{}^{\mathrm{+}}\text{, X}{}^{\mathrm{?}}\text{) ?}\text{V}\text{?}\text{[MX(s)]ΔK =}\text{K}\text{?}\text{(M}{}^{\mathrm{+}}\text{, X}{}^{\mathrm{?}}\text{) ?}\text{K}\text{?}\text{[MX(s)]}$

Values of the partial molal volume ($\text{V}\text{?}$) and compressibilty ($\text{K}\text{?}$) in water and seawater have been tabulated for some ions from 0 to 50°C. The compressibility change is quite large (~10 × 10^?3 cm³ bar^?1 mol^?1) for the solubility of most minerals. This large compressibility change accounts for the large differences observed between values of ΔV obtained from linear plots of In K_sp versus P and molal volume data (Macdonald and North, 1974; North, 1974). Calculated values of $\text{K}{}^{\mathrm{P}}{}_{\mathrm{sp}}\text{K}{}^{\mathrm{o}}{}_{\mathrm{sp}}$ for the solubility of CaCO₃, SrSO₄ and CaF₂ in water were found to be in good agreement with direct measurements (Macdonald and North, 1974). Similar calculations for the solubility of minerals in seawater are also in good agreement with direct measurements (Ingle, 1975) providing that the surface of the solid phase is not appreciably altered.     

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 thermodynamics of the carbonate system in seawater

The apparent constants (K'_i) for the ionization of carbonic acid in seawater at various salinities (S,%.) have been fit to equations of the form ln $\text{K'}{}_{\mathrm{i}}\text{= ln K}{}_{\mathrm{i}}\text{+ A}{}_{\mathrm{i}}\text{S}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{+ B}{}_{\mathrm{i}}\text{S}$whereK_i is the thermodynamic ionization constant in water, A_i, and B_i are adjustable parameters. The temperature dependence (TK) of K_i, A_i and B_i were of the form, a₀ + a₁/T + a₃ ln T. Equations of similar forms have been used to analyze the ionization constants for water and boric acid and the solubility product of calcite in seawater. The effect of pressure on the apparent constants ($\text{K}{}^{\mathrm{p}}{}_{\mathrm{i}}\text{K}{}^{\mathrm{o}}{}_{\mathrm{i}}$) have been fit to equations of the form ln $\text{(}\text{K}{}^{\mathrm{p}}{}_{\mathrm{i}}\text{K}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{) = ? (ΔVP + 0.5 ΔK P}{}^2\text{)/RT}$ where the volume (ΔV) and compressibility (ΔK) changes are polynomial functions of temperature. The equations generated for various açids in seawater have been used to examine the carbonate system in seawater. Equations relating the NBS and Tris pH scales have been derived as well as equations of pH as a function of temperature and pressure. The equations from Hansson (1972, Ph.D. Thesis, University of Göteborg, Sweden) and Mehrbachet al. (1973, Limnol. Oceanogr.18, 897–907) have been used to examine the components of the carbonate system. At a fixed total alkalinity and total carbon dioxide, differences of ±0.01 m-equiv kg^?1 in HCO^?₃ and CO^2?₃ were found; however, the [CO₂] and P_co2 are nearly the same. The contribution of borate ion, B(OH)^?₄ determined from the equations of Hansson (1972, Ph.D. Thesis, University of Göteborg, Sweden) and Lyman (1957, Ph.D. Thesis, University of California, Los Angeles) differ by ±0.01 m-equiv kg^?1 for waters with the same salinity and temperature.     

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.     

Constantes de formation des complexes hydroxydés de l'aluminium en solution aqueuse de 20 a 70°C

Stability constants of hydroxocomplexes of Al(III):Al(OH)²⁺ and A1(OH)₄^? have been measured in the 20–70°C temperature range by reactions involving only dissolved species. The stability constant ${}^1\text{K}{}_1$ of the first complex ion is studied by measuring pH of solutions of aluminium salts at several concentrations. ${}^1\text{β}{}_4$ of aluminate ion is deduced from equilibrium constants of the reaction between the trioxalato aluminium (III) complex ion and Al³⁺ in acid medium, and between the same complex ion and A1(OH)₄^? in alkaline medium. The K values and the associated ΔH are ${}^1\text{K}{}_1\text{= 10}{}^{\mathrm{?5.00}}$ and ΔH₁ = 11.8 Kcal; ${}^1\text{β}{}_4\text{= 10}{}^{\mathrm{?22.20}}$ and ΔH₄ = 42.45 Kcal. These last results are not in agreement with the values of recent tables for $\text{ΔG}{}^0\text{?}$ and $\text{ΔH}{}^0\text{?}$ of Al³⁺ and Al(OH)₄^?. We suggest a consistent set of data for dissolved and solid Al species and for some aluminosilicates.     

Gibbsite solubility and thermodynamic properties of hydroxy-aluminum ions in aqueous solution at 25°C

Solubility curves were determined for a synthetic gibbsite and a natural gibbsite (Minas Gerais, Brazil) from pH 4 to 9, in 0.2% gibbsite suspensions in 0.01 M NaNO₃ that were buffered by low concentrations of non-complexing buffer agents. Equilibrium solubility was approached from oversaturation (in suspensions spiked with Al(NO₃)₃ solution), and also from undersaturation in some synthetic gibbsite suspensions. Mononuclear Al ion concentrations and pH values were periodically determined. Within 1 month or less, data from over-and undersaturated suspensions of synthetic gibbsite converged to describe an equilibrium solubility curve. A downward shift of the solubility curve, beginning at pH 6.7, indicates that a phase more stable than gibbsite controls Al solubility in alkaline systems. Extrapolation of the initial portion of the high-pH side of the synthetic gibbsite solubility curve provides the first unified equilibrium experimental model of Al ion speciation in waters from pH 4 to 9.The significant mononuclear ion species at equilibrium with gibbsite are Al³⁺, AlOH²⁺, Al(OH)⁺₂ and Al(OH)^?₄, and their ion activity products are $\text{1K}{}_{50}\text{= 1.29 × 10}{}^8$, $\text{1K}{}_{\mathrm{s1}}\text{= 1.33 × 10}{}^3$, $\text{1K}{}_{\mathrm{s2}}\text{= 9.49 × 10}{}^{\mathrm{?3}}$ and $\text{1K}{}_{\mathrm{s4}}\text{= 8.94 × 10}{}^{\mathrm{?15}}$. The calculated standard Gibbs free energies of formation (ΔG°_f) for the synthetic gibbsite and the A1OH²⁺, Al(OH)⁺₂ and Al(OH)^?₄ ions are ?276.0, ?166.9, ?216.5 and ?313.5 kcal mol^?1, respectively. These ΔG°_f values are based on the recently revised ΔG°_f value for Al³⁺ (?117.0 ± 0.3 kcal mol_?1) and carry the same uncertainty. The ΔG°_f of the natural gibbsite is ?275.1 ± 0.4 kcal mol^?, which suggests that a range of ΔG°_f values can exist even for relatively simple natural minerals.     

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}}$

     

8.

Thermodynamic analyses of equilibria involving olivine,orthopyroxene and garnet     

Toshisuke Kawasaki  Yoshito Matsui 《Geochimica et cosmochimica acta》1983,47(10):1661-1679

     

9.

Theoretical prediction of phase relations among aqueous solutions and minerals: Salton Sea geothermal system     

Dennis K. Bird  Denis L. Norton 《Geochimica et cosmochimica acta》1981,45(9):1479-1494

     

10.

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

E.J. Reardon 《Chemical Geology》1976,18(4):309-325

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.     

11.

The distribution of divalent trace elements between sulfides,oxides, silicates and hydrothermal solutions: I. Thermodynamic basis     

Dimitrl A. Sverjensky 《Geochimica et cosmochimica acta》1985,49(3):853-864

Experimental data for the standard Gibbs free energies of formation from the elements of a wide variety of metal sulfides and oxides, spinels, olivines and pyroxenes at 25°C and 1 bar define linear correlations, within about ±900 cal·mole^?1, with the corresponding conventional standard partial molal Gibbs free energies of formation of the aqueous M²⁺ cations of the form where a_MaZ and b_MaZ are empirically determined constants characteristic of the structure M_nZ. The only exceptions to correlations of this type are compounds of the heavy alkaline earths Ca, Sr and Ba, which appear to follow correlations with cation radius instead. The linear free energy correlations enable prediction of standard Gibbs free energies of formation of compositional end-members of a particular structure M_nZ provided that a_MaZ and b_MaZ are known accurately. When only the free energy of the Mg end-member is known, the standard Gibbs free energy of formation at 25°C and 1 bar of the Fe endmember, and hence a_MaZ and b_MaZ Can be predicted from the temperature independence of and estimated entropies and heat capacities for the Fe end-member. Using this approach, the free energies of ferrosilite, hedenbergite and annite at 25°C and 1 bar were predicted to within ±1000 cal·mole^?1 of the helgesonet al. (1978) values. Free energies of formation of talc (M₃Si₄O₁₀(OH)₂), clinchlore (M₅Al₂Si₃O₁₀(OH)₈), and tremolite (Ca₂M₅(Si₄O₁₁)₂(OH)₂)-type compounds where M is Mg, Mn, Zn, Fe, Co, or Ni were then predicted at 25°C and 1 bar.Calculation of the equilibrium distribution of Mg, Zn and Sr between galena and hydrothermal solution, and Zn, Mg, Fe and Mn between chlorite and hydrothermal solution demonstrates: (1) that the Sr contents of low temperature galenas (e.g. Mississippi Valley-type) should be negligible (reported analyses of Sr content and Sr isotopic composition of such galenas are probably attributable to fluid inclusions or carbonate inclusions); and (2), that the Zn contents of hydrothermal chlorites in a model of the midoceanic ridge hydrothermal systems are sensitive to temperature, to complexing in the aqueous phase, and to the overall Fe/Mg ratio of the chlorite.     

12.

Rutile solubility and titanium coordination in silicate melts     

James E Dickinson  Paul C Hess 《Geochimica et cosmochimica acta》1985,49(11):2289-2296

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.     

13.

The solubilities of calcite,aragonite and vaterite in CO₂-H₂O solutions between 0 and 90°C,and an evaluation of the aqueous model for the system CaCO₃-CO₂-H₂O     

L.Niel Plummer  Eurybiades Busenberg 《Geochimica et cosmochimica acta》1982,46(6):1011-1040

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.     

14.

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

Frank.J. Millero  Robert H. Byrne 《Geochimica et cosmochimica acta》1984,48(5):1145-1150

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.     

15.

Thermodynamics of brine-salt equilibria — II. The system NaCl-KCl-H₂O from 0 to 200°C     

James R Wood 《Geochimica et cosmochimica acta》1976,40(10):1211-1220

Thermodynamic functions describing salt solubilities and activities in the system NaCl-KCl-H₂O have been obtained over the temperature range 0–200°C using existing isopiestic and solubility data. The two-parameter equation used previously at 25°C was sufficient to characterize the activity coefficient of the aqueous binary systems from infinite dilution to saturation, and the mixing behavior could be described by Harned's Rule over this temperature range. Activity coefficients have been calculated for NaCl-KCl-H₂O mixtures as a function of temperature using the expression: $\text{log γ}{}_{\mathrm{\pm A}}\text{= ? Aγ√I/(1+}\text{a}\text{?}\text{Bγ√I) + B}{}^{\mathrm{.}}\text{I + α}{}_{\mathrm{AB}}\text{I}$. Calculated solubility curves are in good agreement with the experimental data up to 200°C and estimates of ΔH⁰_(sat) and ΔC⁰_p(sat) calculated from the first and second temperature derivatives of log K are consistent with measured standard state enthalpies and heat capacities for both NaCl and KCl. The absolute value of the mixing parameters, α_AB, decrease with increasing temperature and approach zero for both salts in the vicinity of 200°C. The excess free energy of mixing thus approaches zero and may be due to the transition from a dissociated electrolyte solution to one in which associated ions or molecular species are beginning to dominate over ‘free’ ions. The mixing parameter decreases more rapidly for KCl than for NaCl and suggests that KCl-rich systems become associated at lower temperatures than NaCl-rich systems.     

16.

The use of the specific interaction model to estimate the partial molal volumes of electrolytes in seawater     

Frank J Millero 《Geochimica et cosmochimica acta》1977,41(2):215-223

The specific interaction model has been used to determine the partial molal volume of electrolytes in 0.725 m NaCl and 35‰ salinity seawater solutions at 25°C. The partial molal volumes of electrolytes (MX) were estimated at a given ionic strength (I) from

$\text{V}\text{(}\text{MX}\text{) =}\text{V}{}^0\text{(}\text{MX}\text{) +}\text{S}{}_{\mathrm{v}}\text{I}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{(1 + I}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}\text{+ v}{}_{\mathrm{M}}\text{∑}\text{X}\text{B}{}_{\mathrm{<mtext>MX</mtext>}}\text{[X] + v}{}_{\mathrm{X}}\text{∑}\text{M}\text{B}{}_{\mathrm{<mtext>MX</mtext>}}\text{[X]}$

, where S_V is the Debye-Hückel limiting law slope, v_i is the number of ions i formed when MX dissociated, [i] is the total molality of ion i and B_MX is a specific interaction parameter that varies slowly with ionic strength. The values of $\text{V}\text{(}\text{MX}\text{)}$ estimated by using this equation were found to agree very well with experimental values in NaCl and seawater providing there are not strong interactions between M and X. For electrolytes that form ion pairs (i.e. MX°) corrections must be made. Methods are discussed for making these corrections.     

17.

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

Ralph Kretz 《Geochimica et cosmochimica acta》1982,46(3):411-421

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

18.

Complex formation in the ternary system Mg(II)-CO₂-H₂O     

Walter Riesen  Heinz Gamsjäger  Paul W. Schindler 《Geochimica et cosmochimica acta》1977,41(9):1193-1200

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.     

19.

The effect of pressure on the ionization of boric acid in sodium chloride and seawater from molal volume data at 0 and 25°C     

Gary K Ward  Frank J Millero 《Geochimica et cosmochimica acta》1975,39(12):1595-1604

The apparent molal volume, φ_V of boric acid, B(OH)₃ and sodium borate, NaB(OH)₄, have been determined in 35%. salinity seawater and 0·725 molal NaCl solutions at 0 and 25°C from precise density measurements. Similar to the behavior of nonelectrolytes and electrolytes in pure water, the φ_V of B(OH)₃ is a linear function of added molality and the φ_V of NaB(OH)₄ is a linear function of the square root of added molarity in seawater and NaCl solutions. The partial molal volumes, $\text{V}\text{?}\text{1}$, of B(OH)₃ and NaB(OH)₄ in seawater and NaCl were determined from the φ_V's by extrapolating to infinite dilution in the medium. The $\text{V}\text{?}\text{1}$ of B(OH)₃ is larger in NaCl and seawater than pure water apparently due to the ability of electrolytes to dehydrate the nonelectrolyte B(OH)₃. The $\text{V}\text{?}\text{1}$ of NaB(OH)₄ in itself, NaCl and seawater is larger than the expected value at 0·725 molal ionic strength due to ion pair formation [Na⁺ + B(OH)₄^? → NaB(OH)₄⁰]. The volume change for the formation of NaB(OH)₄⁰ in itself and NaCl was found to be equal to 29·4 ml mol^?1 at 25°C and 0·725 molal ionic strength. These large $\text{Δ}\text{V}\text{?}\text{1}$'s indicate that at least one water molecule is released when the ion pair is formed [Na⁺ + B(OH)₄^? → H₂O + NaOB(OH)₂⁰]. The observed $\text{V}\text{?}\text{1}$ in seawater and the $\text{Δ}\text{V}\text{?}\text{1}$ (NaB⁰) in water and NaCl were used to estimate $\text{Δ}\text{V}\text{?}\text{1 (MgB}{}^{\mathrm{+}}\text{) = Δ}\text{V}\text{?}\text{1 (CaB}{}^{\mathrm{+}}\text{) = 38·4 ml mol}{}^{\mathrm{?1}}$ for the formation of MgB⁺ and CaB⁺. The volume change for the ionization of B(OH)₃ in NaCl and seawater was determined from the molal volume data. Values of $\text{Δ}\text{V}\text{?}\text{1 = ?29·2 and ?25·9 ml mol}{}^{\mathrm{?1}}$ were found in seawater and $\text{Δ}\text{V}\text{?}\text{1 = ?21·6 and ?26·4}$ in NaCl, respectively, at 0 and 25°C. The effect of pressure on the ionization of B(OH)₃ in NaCl and seawater at 0 and 25°C determined from the volume change is in excellent agreement with direct measurements in artificial seawater (culberson and Pytkowicz, 1968; Disteche and Disteche, 1967) and natural seawater (Culberson and Pytkowicz, 1968).     

20.

Sodium and potassium coprecipitation in aragonite     

Art F White 《Geochimica et cosmochimica acta》1977,41(5):613-625

The coprecipitation of Na and K was experimentally investigated in aragonite. The distribution functions were determined at pH 6.8 and 8.8 over aqueous Na and K concentrations of between 5 × 10^?4and 2.0 M and temperatures of between 25 and 75°C.The mole fractions of Na and K in aragonite are related to the aqueous ratios of Na and Ca by a function of the form

$\text{log}\text{X}{}_{\mathrm{Na}}{}_2\text{CO}{}_3\text{,K}{}_2\text{CO}{}_3\text{= C}{}_0\text{+ C}{}_1\text{log}\text{a}{}^2{}_{\mathrm{Na}}\text{? ,K}{}^{\mathrm{?}}\text{a}{}_{\mathrm{Ca}}{}^{\mathrm{2+}}$

where C₀ and C₁ are constants at a given temperature. This equation was derived by a statistical model assuming a heterogeneous energy distribution for the sites of incorporation. The independence of the coprecipitation process from aqueous anion activities suggests that carbonate is the only anionic species in the solid solution.     

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

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.     

The distribution of divalent trace elements between sulfides,oxides, silicates and hydrothermal solutions: I. Thermodynamic basis

Experimental data for the standard Gibbs free energies of formation from the elements of a wide variety of metal sulfides and oxides, spinels, olivines and pyroxenes at 25°C and 1 bar define linear correlations, within about ±900 cal·mole^?1, with the corresponding conventional standard partial molal Gibbs free energies of formation of the aqueous M²⁺ cations of the form where a_MaZ and b_MaZ are empirically determined constants characteristic of the structure M_nZ. The only exceptions to correlations of this type are compounds of the heavy alkaline earths Ca, Sr and Ba, which appear to follow correlations with cation radius instead. The linear free energy correlations enable prediction of standard Gibbs free energies of formation of compositional end-members of a particular structure M_nZ provided that a_MaZ and b_MaZ are known accurately. When only the free energy of the Mg end-member is known, the standard Gibbs free energy of formation at 25°C and 1 bar of the Fe endmember, and hence a_MaZ and b_MaZ Can be predicted from the temperature independence of and estimated entropies and heat capacities for the Fe end-member. Using this approach, the free energies of ferrosilite, hedenbergite and annite at 25°C and 1 bar were predicted to within ±1000 cal·mole^?1 of the helgesonet al. (1978) values. Free energies of formation of talc (M₃Si₄O₁₀(OH)₂), clinchlore (M₅Al₂Si₃O₁₀(OH)₈), and tremolite (Ca₂M₅(Si₄O₁₁)₂(OH)₂)-type compounds where M is Mg, Mn, Zn, Fe, Co, or Ni were then predicted at 25°C and 1 bar.Calculation of the equilibrium distribution of Mg, Zn and Sr between galena and hydrothermal solution, and Zn, Mg, Fe and Mn between chlorite and hydrothermal solution demonstrates: (1) that the Sr contents of low temperature galenas (e.g. Mississippi Valley-type) should be negligible (reported analyses of Sr content and Sr isotopic composition of such galenas are probably attributable to fluid inclusions or carbonate inclusions); and (2), that the Zn contents of hydrothermal chlorites in a model of the midoceanic ridge hydrothermal systems are sensitive to temperature, to complexing in the aqueous phase, and to the overall Fe/Mg ratio of the chlorite.     

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.     

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.     

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.     

Thermodynamics of brine-salt equilibria — II. The system NaCl-KCl-H2O from 0 to 200°C

Thermodynamic functions describing salt solubilities and activities in the system NaCl-KCl-H₂O have been obtained over the temperature range 0–200°C using existing isopiestic and solubility data. The two-parameter equation used previously at 25°C was sufficient to characterize the activity coefficient of the aqueous binary systems from infinite dilution to saturation, and the mixing behavior could be described by Harned's Rule over this temperature range. Activity coefficients have been calculated for NaCl-KCl-H₂O mixtures as a function of temperature using the expression: $\text{log γ}{}_{\mathrm{\pm A}}\text{= ? Aγ√I/(1+}\text{a}\text{?}\text{Bγ√I) + B}{}^{\mathrm{.}}\text{I + α}{}_{\mathrm{AB}}\text{I}$. Calculated solubility curves are in good agreement with the experimental data up to 200°C and estimates of ΔH⁰_(sat) and ΔC⁰_p(sat) calculated from the first and second temperature derivatives of log K are consistent with measured standard state enthalpies and heat capacities for both NaCl and KCl. The absolute value of the mixing parameters, α_AB, decrease with increasing temperature and approach zero for both salts in the vicinity of 200°C. The excess free energy of mixing thus approaches zero and may be due to the transition from a dissociated electrolyte solution to one in which associated ions or molecular species are beginning to dominate over ‘free’ ions. The mixing parameter decreases more rapidly for KCl than for NaCl and suggests that KCl-rich systems become associated at lower temperatures than NaCl-rich systems.     

The use of the specific interaction model to estimate the partial molal volumes of electrolytes in seawater

The specific interaction model has been used to determine the partial molal volume of electrolytes in 0.725 m NaCl and 35‰ salinity seawater solutions at 25°C. The partial molal volumes of electrolytes (MX) were estimated at a given ionic strength (I) from

$\text{V}\text{(}\text{MX}\text{) =}\text{V}{}^0\text{(}\text{MX}\text{) +}\text{S}{}_{\mathrm{v}}\text{I}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{(1 + I}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}\text{+ v}{}_{\mathrm{M}}\text{∑}\text{X}\text{B}{}_{\mathrm{<mtext>MX</mtext>}}\text{[X] + v}{}_{\mathrm{X}}\text{∑}\text{M}\text{B}{}_{\mathrm{<mtext>MX</mtext>}}\text{[X]}$

, where S_V is the Debye-Hückel limiting law slope, v_i is the number of ions i formed when MX dissociated, [i] is the total molality of ion i and B_MX is a specific interaction parameter that varies slowly with ionic strength. The values of $\text{V}\text{(}\text{MX}\text{)}$ estimated by using this equation were found to agree very well with experimental values in NaCl and seawater providing there are not strong interactions between M and X. For electrolytes that form ion pairs (i.e. MX°) corrections must be made. Methods are discussed for making these corrections.     

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

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 effect of pressure on the ionization of boric acid in sodium chloride and seawater from molal volume data at 0 and 25°C

The apparent molal volume, φ_V of boric acid, B(OH)₃ and sodium borate, NaB(OH)₄, have been determined in 35%. salinity seawater and 0·725 molal NaCl solutions at 0 and 25°C from precise density measurements. Similar to the behavior of nonelectrolytes and electrolytes in pure water, the φ_V of B(OH)₃ is a linear function of added molality and the φ_V of NaB(OH)₄ is a linear function of the square root of added molarity in seawater and NaCl solutions. The partial molal volumes, $\text{V}\text{?}\text{1}$, of B(OH)₃ and NaB(OH)₄ in seawater and NaCl were determined from the φ_V's by extrapolating to infinite dilution in the medium. The $\text{V}\text{?}\text{1}$ of B(OH)₃ is larger in NaCl and seawater than pure water apparently due to the ability of electrolytes to dehydrate the nonelectrolyte B(OH)₃. The $\text{V}\text{?}\text{1}$ of NaB(OH)₄ in itself, NaCl and seawater is larger than the expected value at 0·725 molal ionic strength due to ion pair formation [Na⁺ + B(OH)₄^? → NaB(OH)₄⁰]. The volume change for the formation of NaB(OH)₄⁰ in itself and NaCl was found to be equal to 29·4 ml mol^?1 at 25°C and 0·725 molal ionic strength. These large $\text{Δ}\text{V}\text{?}\text{1}$'s indicate that at least one water molecule is released when the ion pair is formed [Na⁺ + B(OH)₄^? → H₂O + NaOB(OH)₂⁰]. The observed $\text{V}\text{?}\text{1}$ in seawater and the $\text{Δ}\text{V}\text{?}\text{1}$ (NaB⁰) in water and NaCl were used to estimate $\text{Δ}\text{V}\text{?}\text{1 (MgB}{}^{\mathrm{+}}\text{) = Δ}\text{V}\text{?}\text{1 (CaB}{}^{\mathrm{+}}\text{) = 38·4 ml mol}{}^{\mathrm{?1}}$ for the formation of MgB⁺ and CaB⁺. The volume change for the ionization of B(OH)₃ in NaCl and seawater was determined from the molal volume data. Values of $\text{Δ}\text{V}\text{?}\text{1 = ?29·2 and ?25·9 ml mol}{}^{\mathrm{?1}}$ were found in seawater and $\text{Δ}\text{V}\text{?}\text{1 = ?21·6 and ?26·4}$ in NaCl, respectively, at 0 and 25°C. The effect of pressure on the ionization of B(OH)₃ in NaCl and seawater at 0 and 25°C determined from the volume change is in excellent agreement with direct measurements in artificial seawater (culberson and Pytkowicz, 1968; Disteche and Disteche, 1967) and natural seawater (Culberson and Pytkowicz, 1968).     

Sodium and potassium coprecipitation in aragonite

The coprecipitation of Na and K was experimentally investigated in aragonite. The distribution functions were determined at pH 6.8 and 8.8 over aqueous Na and K concentrations of between 5 × 10^?4and 2.0 M and temperatures of between 25 and 75°C.The mole fractions of Na and K in aragonite are related to the aqueous ratios of Na and Ca by a function of the form

$\text{log}\text{X}{}_{\mathrm{Na}}{}_2\text{CO}{}_3\text{,K}{}_2\text{CO}{}_3\text{= C}{}_0\text{+ C}{}_1\text{log}\text{a}{}^2{}_{\mathrm{Na}}\text{? ,K}{}^{\mathrm{?}}\text{a}{}_{\mathrm{Ca}}{}^{\mathrm{2+}}$

where C₀ and C₁ are constants at a given temperature. This equation was derived by a statistical model assuming a heterogeneous energy distribution for the sites of incorporation. The independence of the coprecipitation process from aqueous anion activities suggests that carbonate is the only anionic species in the solid solution.