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.     

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.     

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.     

Distribution coefficients of Eu and Sr for plagioclase-liquid and clinopyroxene-liquid equilibria in oceanic ridge basalt: an experimental study

The distribution coefficients of Eu and Sr for plagioclase-liquid and clinopyroxene-liquid pairs as a function of temperature and oxygen fugacity were experimentally investigated using an oceanic ridge basalt enriched with Eu and Sr as the starting material. Experiments were conducted between 1190° and 1140°C over a range of oxygen fugacities between 10^?8 and 10^?14 atm.The molar distribution coefficients are given by the equations: log^PL_Sr = 7320/T ? 4.62logK^CPX_Sr = 18020/T ? 13.10. Similarly, the weight fraction distribution coefficients are given by the equations: logD^PL_Sr = 6570/T ? 4.30logD^CPX_Sr = 18434/T ? 13.62.Although the mole fraction distribution coefficients have a smaller dependence on bulk composition than do the weight fraction distribution coefficients, they are not independent of bulk composition, thereby restricting the application of these experimental results to rocks similar to oceanic ridge basalts in bulk composition.Because the Sr distribution coefficients are independent of oxygen fugacity, they may be used as geothermometers. If the temperature can be determined independently — for example, with the Sr distribution coefficients, the Eu distribution coefficients may be used as oxygen geobarometers. Throughout the range of oxygen fugacities ascribed to terrestrial and lunar basalts, plagioclase concentrates Eu but clinopyroxene rejects Eu.     

The role of CO2 in the chemical modification of deep continental crust

Compositional differences between granulite facies rocks and equivalent amphibolite facies rocks and the observation of CO₂-rich fluid inclusions in granulites, have led to the suggestion that CO₂ must play a role in modifying the composition of deep continental crust. How CO₂ effects this change has remained unclear. Using the thermodynamic properties of aqueous ions in a fluid of evolving $\text{CO}{}_2\text{H}{}_2\text{O}$ ratio, it is possible to model the incongruent dissolution of feldspars under conditions appropriate for granulite facies metamorphism. The results demonstrate that dissolution will be strongly enhanced at high $\text{CO}{}_2\text{H}{}_2\text{O}$ ratios, with ion solubilities being Na⁺ >K⁺ ? Ca⁺⁺. This enhancement is compatible with the reported compositional contrasts between granulite and amphibolite facies rock, but requires large fluid volumes.To test the dissolution model, a detailed field and petrologic study was conducted in a well exposed granulite facies terrane in West Greenland. Strong correlation between fluid composition and bulk rock chemistry can be documented; CO₂-rich regions contain rocks which consistently have low ratios, while H₂O-rich regions consistently have high ratios. Magnetite rims on sulfide grains are ubiquitous in high regions and are absent in high regions, and they provide evidence that CO₂ was introduced into the region. These correlations and observations are predictable from the properties of the dissolution process. These considerations, along with observations regarding graphite petrogenesis, provide strong arguments that the total fluid volume interacting with the rock during metamorphism was very large, in some cases equaling or exceeding total rock volume. Such large fluid volumes can lead to significant compositional modification of the crust, and will mask the original protolith chemistry. Such processes should lead to Ca- and Al-enriched, Na-, K-, S- and Si-depleted residues in the deep crust.     

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.     

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.     

Calcium diffusion in a simple silicate melt to 30 kbar

Calcium-45 was used as a radiotracer to measure self-diffusion coefficients for Ca in a sodium-calcium-aluminosilicate melt (29% Na₂O, 5% CaO, 10% Al₂O₃, 56% SiO₂) at temperatures in the range 1100–1400°C and pressures to 30 kbar. Calcium diffusivity (D_Ca) was found to depend upon both temperature and pressure in a complex but systematic manner: $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?P}\text{)}{}_{\mathrm{T}}$ is always negative and has a larger absolute value at lower temperatures; $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?T}\text{)}{}_{\mathrm{P}}$ is positive and increases with increasing pressure. The overall dependence of D_Ca upon T and P is given approximately by $\text{D}{}_{\mathrm{caT.P}}\text{= [0.0025 exp(}\text{-23,107}\text{RT}\text{)] exp [P}\text{(0.7297T ? 1261.32)}\text{RT}\text{]}$. When expressed in terms of volume (V_a) and energy (E) of activation, the results are as follows: V_a ranges from 2.2 cm³/mole at 1400°C to 11.9 cm³/mole at 1100°C. and E ranges from 25.4 kcal/mole (1 kban to 49.8 kcal/mole (20 kbar).From the systematic dependence of D_Ca upon T and P, it is concluded that diffusion of Ca²⁺ in silicate melts does not take place by means of a vacant site mechanism, but is controlled instead by the amount and distribution of free volume in the melt structure.If it is assumed that the viscosity of the melt used in this study decreases with increasing pressure (Kushiro, 1976, J. Geophys. Res.81, 6351–6356) as D_Ca does, then the Stokes-Einstein inverse relation between viscosity and diffusivity is clearly violated, and its validity for silicate melts must be questioned. Thus, it appears that in silicate melts, unlike many liquids, viscous flow and diffusion are fundamentally different transport processes, involving different structural units.The effect of pressure on calcium diffusion is too small to invalidate kinetic models of upper mantle processes that have been based upon diffusivity values measured at 1 atm. Pressure may, however, induce significant reductions in the diffusion rates of large ions such as Rb⁺ or SiO^4?₄ in silicate melts.     

The significance of metamorphic fluorite in the Adirondacks

Thermodynamic calculations for selected silicate-oxide-fluorite assemblages indicate that several commonly occurring fluorite-bearing assemblages are restricted to relatively narrow fields at constant P?T. The presence of fayalite-ferrohedenbergite-fluorite-quartz ± magnetite and ferrosalite-fluorite-quartz-magnetite assemblages in orthogneisses from Au Sable Forks, Wanakena and Lake Pleasant, New York, buffered fluorine and oxygen fugacities during the granulite facies metamorphism in the Adirondack Highlands. These buffering assemblages restrict$\text{?}{}_{\mathrm{F}}{}_2$ to 10^{?29 ± 1} bar and $\text{?}{}_0{}_2$ to 10^{?16 ± 1} bar at the estimated metamorphic temperature of 1000K and pressure of 7 kbar. The assemblage biotite-magnetite-ilmenite-K-feldspar, found in the same Au Sable Forks outcrop as the fayalite-fluorite-ferrohedenbergite-quartz-magnetitie assemblage, restricts H₂O fugacities to less than 10^3·3 bar. These fugacities limit H₂ and HF fugacities to less than 10¹ bar for the Au Sable outcrop. The data indicate that relative to H₂O, O₂, H₂, F₂ and HF are not major species in the fluid equilibrated with Adirondack orthogneisses. The calculated F₂ fugacilies are similar to the upper limits possible for plagioclase-bearing rocks and probably represent the upper $\text{?}{}_{\mathrm{F}}{}_2$ limit for metamorphism in the Adirondacks and in other granulite facies terranes.     

The ionization of acids in estuarine waters

The stoichiometric, $\text{K}{}_{\mathrm{HA}}{}^1$, and apparent, K'_HA, constants for the ionization of a number of weak acids (NH₄⁺, HSO₄^?, HF, H₂O, B(OH)₃, H₂CO₃, HCO₃^?, H₃PO₄, H₂PO₄^?, HPO₄², H₃AsO₄ H₂AsO₄^? and HAsO₄^2?) in seawater at 25°C diluted with water have been fitted to equations of the form (Millero, 1979). In $\text{K}{}_{\mathrm{HA}}{}^1\text{= In K}{}_{\mathrm{HA}}\text{+ AS}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{+ BS}$ where In K_HA is the thermodynamic constant in water, S is the salinity, A and B are adjustable parameters. The validity of this equation in estuarine waters has been examined by using an ion pairing model (Millero and Schreiber, 1981). The calculated values of $\text{K}{}_{\mathrm{HA}}{}^1$ and K'_HA at S = 35%. are in good agreement with the measured values for all the systems examined. The equation used to extrapolate the measured values to pure water K_HA predicted values that agreed with those determined by using the ion pairing model. The exception was the ionization of HPO₄^2? due to the strong interactions of Ca²⁺ and Mg²⁺ with PO₄^3?. The differences in the predicted values of $\text{K}{}_{\mathrm{HA}}{}^1$ in seawater diluted with pure water and average river water were very small for all the acids except HPO₄^2? (the maximum ΔpK = 0.96 in average river water). The larger difference in the $\text{K}{}_{\mathrm{HA}}{}^1$ for HPO₄^2? in river waters is due to the strong interactions of Ca²⁺ and PO₄^3?.     

Estimation of the dielectric constant of H2O from experimental solubilities of quartz,and calculation of the thermodynamic properties of aqueous species to 900°C at 2 kb

Experimental quartz solubilities in H₂O (Anderson and Burnham, 1965, 1967) were used together with equations of state for quartz and aqueous species (Helgesonet al., 1978; Walther and Helgeson, 1977) to calculate the dielectric constant of H₂O () at pressures and temperatures greater than those for which experimental measurements (Heger, 1969; Lukashovet al., 1975) are available ($\text{0.001 ? P ? 5}\text{kb}$ and $\text{0 ? T ? 600°C}$). Estimates of computed in this way for 2 kb (which are the most reliable) range from 9.6 at 600°C to 5.6 at 800°C. These values are 0.5 and 0.8 units greater, respectively, than corresponding values estimated by Quist and Marshall (1965), but they differ by <0.3 units from extrapolated values computed from Pitzer's (1983) adaptation of the Kirkwood (1939) equation. The estimates of generated from quartz solubilities at 2 kb were fit with a power function of temperature, which was then used together with equations and data given by Helgeson and Kirkham (1974a,b, 1976) Helgesonet al. (1981), and Helgeson (1982b, 1984) to calculate Born functions, Debye Hückel parameters, and the thermodynamic properties of Na⁺, K⁺, Mg⁺⁺, Ca⁺⁺, and other aqueous species of geologic interest at temperatures to 900°C.     

Evaluation of deep temperatures of hydrothermal systems by a new gas geothermometer

The chemical composition of gas mixtures emerging in thermal areas can be used to evaluate the deep thermal temperatures. Chemical analyses of the gas compositions for 34 thermal systems were considered and an empirical relationship developed between the relative concentrations of H₂S, H₂, CH₄ and CO₂ and the reservoir temperature. The evaluated temperatures can be expressed by: $\text{t°C =}\text{24775}\text{α + β + 36.05}\text{?273}$ where $\text{α = 2 log}\text{CH}{}_4\text{CO}{}_2\text{?log}\text{H}{}_2\text{CO}{}_2\text{?3 log}\text{H}{}_2\text{S}\text{CO}{}_2$ (concentrations in % by volume) and β = 7 logP_co2     

The kinetics of silica-water reactions

A differential rate equation for silica-water reactions from 0–300°C has been derived based on stoichiometry and activities of the reactants in the reaction SiO₂(s) + 2H₂O(l) = H₄SiO₄(aq)

where ($\text{A}\text{M}$) = (the relative interfacial area between the solid and aqueous phases/the relative mass of water in the system), and k₊ and k_? are the rate constants for, respectively, dissolution and precipitation. The rate constant for precipitation of all silica phases is $\text{log k}{}_{\mathrm{?}}\text{= ? 0.707 ?}\text{2598}\text{T}\text{(T, K)}$ and E_act for this reaction is 49.8 kJ mol^?1. Corresponding equilibrium constants for this reaction with quartz, cristobalite, or amorphous silica were expressed as $\text{log K = a + bT +}\text{c}\text{T}$. Using $\text{K =}\text{k}{}_{\mathrm{+}}\text{k}{}_{\mathrm{?}}$, k was expressed as $\text{log k}{}_{\mathrm{+}}\text{= a + bT +}\text{c}\text{T}$ and a corresponding activation energy calculated:

     

15.

Methane solubilities in multisalt solutions     

Patricia A. Byrne  Ronald K. Stoessell 《Geochimica et cosmochimica acta》1982,46(11):2395-2397

Solubilities of methane in multisalt solutions at 550 psia and 25°C can be predicted from single-salt salting coefficients. The ionic strength contribution of the ith salt, I_i, is multiplied by its molal salting coefficient, k_mi, in the following summation over all salts:

$\text{log}\text{M}{}_{\mathrm{o}}\text{M}{}_{\mathrm{s}}\text{= ∑}{}_{\mathrm{i}}\text{k}{}_{\mathrm{m}}{}_{\mathrm{i}}\text{I}{}_{\mathrm{i}}$

where m_o and m_s are molal methane solubilities in distilled water and the salt solution, respectively, at the T, P and methane fugacity of interest.This equation predicts methane solubility in multisalt brines containing Na⁺, K⁺, Mg⁺², Ca⁺², Cl^?, SO₄^?2 and CO₃^?2 ions. k_mi values reported by Stoessell and Byrne (1982b) can be used in solubility predictions in brines at earth surface conditions. Prediction in reservoir brines would require determination of k_mi, for the different salts at reservoir temperatures and pressures.     

16.

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

Yves Couturier  Gil Michard  Gérard Sarazin 《Geochimica et cosmochimica acta》1984,48(4):649-659

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.     

17.

Precise age determinations and petrogenetic studies using the KCa method     

B.D. Marshall  D.J. DePaolo 《Geochimica et cosmochimica acta》1982,46(12):2537-2545

High precision mass spectrometric determination of calcium isotope ratios allows the ⁴⁰K → ⁴⁰Ca radioactive decay to be used for dating a much broader range of geologic materials than is suggested by previous work. ${}^{40}\text{Ca}{}^{42}\text{Ca}$ is used to monitor enrichments in ⁴⁰Ca and can be measured to ±0.01% (2σ) using an exponential mass discrimination correction (Russell et al., 1978) and large ion currents. The earth's mantle has such a low $\text{K}\text{Ca}$ (～0.01) that it has retained “primordial” ${}^{40}\text{Ca}{}^{42}\text{Ca}\text{= 151.016}$ ± 0.011 (normalized to ${}^{42}\text{Ca}{}^{44}\text{Ca}\text{= 0.31221}$), as determined by measurements on two meteorites, pyroxene from an ultramafic nodule, metabasalt, and carbonatite. ${}^{40}\text{Ca}{}^{42}\text{Ca}$ ratios can be conveniently expressed relative to this value as ?_Ca in units of 10^?4. To test the method for age dating, a mineral isochron has been obtained on a sample of Pikes Peak granite, which has been shown to have concordant KAr, RbSr, and UPb ages. Plagioclase, K-feldspar, biotite, and whole rock yield an age of 1041 ± 32 m.y. (2σ) in agreement with previous age determinations (λ_K = 0.5543 b.y.^?1, $\text{λ}{}_{\mathrm{\beta ?}}\text{λ}{}_{\mathrm{<mtext>K</mtext>}}\text{= 0.8952}$, ⁴⁰K = 0.01167%). The initial ${}^{40}\text{Ca}{}^{42}\text{Ca}$ of 151.024 ± 0.016 (?_Ca = +0.5 ± 1.0), indicates that assimilation of high K/Ca crust was insufficient to affect calcium isotopes. Measurements on two-mica granite from eastern Nevada indicate that the magma sources had K/Ca ≈ 1, similar to intermediate-composition crustal rocks. These results show that the KCa system can be used as a precise geochromometer for common felsic igneous and metamorphic rocks, and may prove applicable to sedimentary rocks containing authigenic K minerals. The relatively short half-life of ⁴⁰K, the non-volatile daughter, and the fact that potassium and calcium are stoichiometric constituents of many minerals, make the KCa system complementary to other dating methods, and potentially applicable to a variety of geologic problems.     

18.

Density calculations for silicate liquids. I. Revised method for aluminosilicate compositions     

Y Bottinga  D Weill  P Richet 《Geochimica et cosmochimica acta》1982,46(6):909-919

Equations are developed for calculating the density of aluminosilicate liquids as a function of composition and temperature. The mean molar volume at reference temperature T_r, is given by $\text{V}{}_{\mathrm{r}}\text{= ∑X}{}_{\mathrm{i}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{+ X}{}_{\mathrm{A}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$, where the summation is taken over all oxide components except A1₂O₃, X stands for mole fraction, terms are constants derived independently from an analysis of volume-composition relations in alumina-free silicate liquids, and $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is the composition-dependent apparent partial molar volume of Al₂O₃. The thermal expansion coefficient of aluminosilicate liquids is given by $\text{α = ∑X}{}_{\mathrm{i}}\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}\text{+ X}{}_{\mathrm{A}}\text{}\text{?}\text{ga}{}_{\mathrm{A}}{}^{\mathrm{o}}$, where $\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}$ terms are constants independent of temperature and composition, and $\text{}\text{?}\text{ga}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is a composition-dependent term representing the effect of Al₂O₃ on the thermal expansion. Parameters necessary to calculate the volume of silicate liquids at any temperature T according to V(T) = V_rexp[α(T-T_r)], where T_r = 1400°C have been evaluated by least-square analysis of selected density measurements in aluminosilicate melts. Mean molar volumes of aluminosilicate liquids calculated according to the model equation conform to experimentally measured volumes with a root mean square difference of 0.28 $\text{cc}\text{mole}$ and an average absolute difference of 0.90% for 248 experimental observations. The compositional dependence of $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is discussed in terms of several possible interpretations of the structural role of Al³⁺ in aluminosilicate melts.     

19.

Intrinsic oxygen fugacity measurements: techniques and results for spinels from upper mantle peridotites and megacryst assemblages     

Richard J. Arculus  John W. Delano 《Geochimica et cosmochimica acta》1981,45(6):899-913

A new technique for the determination of intrinsic oxygen fugacities () of single and polyphase geological samples with solid ZrO₂, oxygen-specific electrolytes is described. Essentially the procedure involves isolating the emf signal from the sample from that unavoidably imposed by the residual atmosphere inside the sample-bearing sensor. By varying the of the residual atmosphere, it is possible to determine a ‘plateau’ value of constant recorded from the sensor which represents a reversed intrinsic measurement for the sample alone, and where the extent of the plateau reflects the innate buffering capability of the sample. A measure of the precision and accuracy of the data obtained is the fact that identical values are obtained whether on a heating or cooling cycle of the sample + compatible atmosphere system.These techniques have been applied to measurements of the intrinsic of spinels from peridotites and megacryst assemblages from Australia, West Germany and the U.S.A. Oxidation states range from ～- 0.25 log₁₀ units more oxidized to 1 log₁₀ unit more reduced than the iron-wüstite (IW) buffer. The overall reduced nature of the spinels and the range of 's obtained are striking features of the data. One implication of the results is that the majority of mantle-derived magmas are initially highly reduced, and the relatively oxidized values observed at surface (～- 4–5 log₁₀ orders more oxidized than IW) reflect late-stage alteration, perhaps by H₂ loss (Sato, 1978).     

20.

The thermodynamics of the carbonate system in seawater     

Frank J Millero 《Geochimica et cosmochimica acta》1979,43(10):1651-1661

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.     

	a	b	c	Eact(kJ mol -1)
Quarts	1.174	-2.028 x 103	-4158	67.4–76.6
α-Cristobalite	-0.739	0	-3586	68.7
β-Cristobalite	-0.936	0	-3392	65.0
Amorphous silica	-0.369	-7.890 x 10-4	3438	60.9–64.9

Methane solubilities in multisalt solutions

Solubilities of methane in multisalt solutions at 550 psia and 25°C can be predicted from single-salt salting coefficients. The ionic strength contribution of the ith salt, I_i, is multiplied by its molal salting coefficient, k_mi, in the following summation over all salts:

$\text{log}\text{M}{}_{\mathrm{o}}\text{M}{}_{\mathrm{s}}\text{= ∑}{}_{\mathrm{i}}\text{k}{}_{\mathrm{m}}{}_{\mathrm{i}}\text{I}{}_{\mathrm{i}}$

where m_o and m_s are molal methane solubilities in distilled water and the salt solution, respectively, at the T, P and methane fugacity of interest.This equation predicts methane solubility in multisalt brines containing Na⁺, K⁺, Mg⁺², Ca⁺², Cl^?, SO₄^?2 and CO₃^?2 ions. k_mi values reported by Stoessell and Byrne (1982b) can be used in solubility predictions in brines at earth surface conditions. Prediction in reservoir brines would require determination of k_mi, for the different salts at reservoir temperatures and pressures.     

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.     

Precise age determinations and petrogenetic studies using the KCa method

High precision mass spectrometric determination of calcium isotope ratios allows the ⁴⁰K → ⁴⁰Ca radioactive decay to be used for dating a much broader range of geologic materials than is suggested by previous work. ${}^{40}\text{Ca}{}^{42}\text{Ca}$ is used to monitor enrichments in ⁴⁰Ca and can be measured to ±0.01% (2σ) using an exponential mass discrimination correction (Russell et al., 1978) and large ion currents. The earth's mantle has such a low $\text{K}\text{Ca}$ (～0.01) that it has retained “primordial” ${}^{40}\text{Ca}{}^{42}\text{Ca}\text{= 151.016}$ ± 0.011 (normalized to ${}^{42}\text{Ca}{}^{44}\text{Ca}\text{= 0.31221}$), as determined by measurements on two meteorites, pyroxene from an ultramafic nodule, metabasalt, and carbonatite. ${}^{40}\text{Ca}{}^{42}\text{Ca}$ ratios can be conveniently expressed relative to this value as ?_Ca in units of 10^?4. To test the method for age dating, a mineral isochron has been obtained on a sample of Pikes Peak granite, which has been shown to have concordant KAr, RbSr, and UPb ages. Plagioclase, K-feldspar, biotite, and whole rock yield an age of 1041 ± 32 m.y. (2σ) in agreement with previous age determinations (λ_K = 0.5543 b.y.^?1, $\text{λ}{}_{\mathrm{\beta ?}}\text{λ}{}_{\mathrm{<mtext>K</mtext>}}\text{= 0.8952}$, ⁴⁰K = 0.01167%). The initial ${}^{40}\text{Ca}{}^{42}\text{Ca}$ of 151.024 ± 0.016 (?_Ca = +0.5 ± 1.0), indicates that assimilation of high K/Ca crust was insufficient to affect calcium isotopes. Measurements on two-mica granite from eastern Nevada indicate that the magma sources had K/Ca ≈ 1, similar to intermediate-composition crustal rocks. These results show that the KCa system can be used as a precise geochromometer for common felsic igneous and metamorphic rocks, and may prove applicable to sedimentary rocks containing authigenic K minerals. The relatively short half-life of ⁴⁰K, the non-volatile daughter, and the fact that potassium and calcium are stoichiometric constituents of many minerals, make the KCa system complementary to other dating methods, and potentially applicable to a variety of geologic problems.     

Density calculations for silicate liquids. I. Revised method for aluminosilicate compositions

Equations are developed for calculating the density of aluminosilicate liquids as a function of composition and temperature. The mean molar volume at reference temperature T_r, is given by $\text{V}{}_{\mathrm{r}}\text{= ∑X}{}_{\mathrm{i}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{+ X}{}_{\mathrm{A}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$, where the summation is taken over all oxide components except A1₂O₃, X stands for mole fraction, terms are constants derived independently from an analysis of volume-composition relations in alumina-free silicate liquids, and $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is the composition-dependent apparent partial molar volume of Al₂O₃. The thermal expansion coefficient of aluminosilicate liquids is given by $\text{α = ∑X}{}_{\mathrm{i}}\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}\text{+ X}{}_{\mathrm{A}}\text{}\text{?}\text{ga}{}_{\mathrm{A}}{}^{\mathrm{o}}$, where $\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}$ terms are constants independent of temperature and composition, and $\text{}\text{?}\text{ga}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is a composition-dependent term representing the effect of Al₂O₃ on the thermal expansion. Parameters necessary to calculate the volume of silicate liquids at any temperature T according to V(T) = V_rexp[α(T-T_r)], where T_r = 1400°C have been evaluated by least-square analysis of selected density measurements in aluminosilicate melts. Mean molar volumes of aluminosilicate liquids calculated according to the model equation conform to experimentally measured volumes with a root mean square difference of 0.28 $\text{cc}\text{mole}$ and an average absolute difference of 0.90% for 248 experimental observations. The compositional dependence of $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is discussed in terms of several possible interpretations of the structural role of Al³⁺ in aluminosilicate melts.     

Intrinsic oxygen fugacity measurements: techniques and results for spinels from upper mantle peridotites and megacryst assemblages

A new technique for the determination of intrinsic oxygen fugacities () of single and polyphase geological samples with solid ZrO₂, oxygen-specific electrolytes is described. Essentially the procedure involves isolating the emf signal from the sample from that unavoidably imposed by the residual atmosphere inside the sample-bearing sensor. By varying the of the residual atmosphere, it is possible to determine a ‘plateau’ value of constant recorded from the sensor which represents a reversed intrinsic measurement for the sample alone, and where the extent of the plateau reflects the innate buffering capability of the sample. A measure of the precision and accuracy of the data obtained is the fact that identical values are obtained whether on a heating or cooling cycle of the sample + compatible atmosphere system.These techniques have been applied to measurements of the intrinsic of spinels from peridotites and megacryst assemblages from Australia, West Germany and the U.S.A. Oxidation states range from ～- 0.25 log₁₀ units more oxidized to 1 log₁₀ unit more reduced than the iron-wüstite (IW) buffer. The overall reduced nature of the spinels and the range of 's obtained are striking features of the data. One implication of the results is that the majority of mantle-derived magmas are initially highly reduced, and the relatively oxidized values observed at surface (～- 4–5 log₁₀ orders more oxidized than IW) reflect late-stage alteration, perhaps by H₂ loss (Sato, 1978).     

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.