Enthalpies of formation at 970 K of compounds in the system MgO-Al2O3-SiO2 from high temperature solution calorimetry

The enthalpies of solution of petrologically important phases in the system MgO-Al₂O₃-SiO ₂ were measured in a melt of composition 2PbO · B₂O₃ at 970 ± 2K. The substances investigated included synthetic and natural (meteoritic) enstatite (MgSiO₃), synthetic aluminous enstatite (MgSiO_30.9Al₂O_30.1), synthetic and natural cordierite (Mg₂Al₄Si₅O₁₈), synthetic and natural sapphirine (approx. 7MgO·9Al₂O₃ · 3SiO₂), synthetic spinel (MgAl₂O₄), natural sillimanite (Al₂SiO₅), synthetic forsterite (Mg₂SiO₄), synthetic pyrope (Mg₃Al₂Si₃O₁₂), natural quartz (SiO₂), synthetic periclase (MgO) and corundum (Al₂O₃). Improvement in standardization of the calorimeter solvent made possible greater precision in this study than obtainable in former work in this laboratory on some of the same substances.The enthalpies of formation of enstatite, synthetic cordierite, forsterite and spinel are in reasonable agreement with values previously determined by solution calorimetry. The enthalpy of formation of enstatite is about 0.7 kcal less negative than the value for clinoenstatite resulting from the HF calorimetry of Torgesen and Sahama (J. Amer. Chem. Soc.70. 2156–2160, 1948), and is in accord with predictions based on analysis of published pyroxene equilibrium work. Aluminous enstatite with 10 wt.% Al₂O₃ shows an enthalpy of solution markedly lower than pure MgSiO₃: the measurements lead to an estimate of the enthalpy of formation at 970 K for MgAl₂SiO₆ (Mg-Tschermak) orthopyroxene of + 9.4 ± 1.5 kcal/mole from MgSiO₃ and Al₂O₃.Comparison of the enthalpies of formation of synthetic cordierite and anhydrous natural low-iron cordierite shows that they are energetically quite similar and that the synthetic cordierite is not likely to have large amounts of (Al, Si) tetrahedral disorder. Comparison of the enthalpies of formation of synthetic sapphirine and natural low-iron sapphirine shows, on the other hand, that the former is not a good stability model for the latter. The lower enthalpy of formation of the high-temperature synthetic sample is undoubtedly a consequence of cation disordering.The enthalpy of formation of natural sillimanite is considerably less negative than given by the tables of Robie andWaldbaum (U.S. Geol. Surv. Bull.1259 1968).The measured enthalpy of formation of synthetic pyrope is consistent with that deduced from published equilibrium diagrams in conjunction with the present measured enthalpy of formation of aluminous enstatite. Calculation of the entropy of synthetic pyrope from the present data yields surprisingly high values and suggests that synthetic pyrope is not a good stability model for natural pyrope-rich garnets. Hence, considerable doubt exists about the direct quantitative application of experimental diagrams involving pyropic garnet to discussions of the garnet stability field in the Earth's outer regions.     

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

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

New experimental data for the system CaO-MgO-SiO2-H2O and a synthesis of inferred phase relations

Subsolidus and vapor-saturated liquidus phase relations for a portion of the system CaO-MgO-SiO₂-H₂O, as inferred from experimental data for the composition regions CaMgSi₂O₆-Mg₂SiO₄-SiO₂-H₂O and CaMgSi₂O₆-Mg₂SiO₄-Ca₃MgSi₂O₈ (merwinite)-H₂O, are presented in pressure-temperature projection. Sixteen invariant points and 39 univariant reactions are defined on the basis of the 1 atm and 10 kbar (vapor-saturated) liquidus diagrams. Lack of experimental control over many of the reactions makes the depicted relations schematic in part.An invariant point involving orthoenstatite, protoenstatite, pigeonite, and diopside (all solid solutions) occurs at low pressure (probably between 1 and 2 kbar). At pressures below this invariant point, orthoenstatite breaks down at high temperature to the assemblage diopside + protoenstatite; with increasing temperature, the latter assemblage reacts to form pigeonite. At pressures above the invariant point, pigeonite forms according to the reaction diopside + orthoenstatite = pigeonite, and the assemblage diopside + protoenstatite is not stable. At 1 atm, both pigeonite and protoenstatite occur as primary liquidus phases, but at pressures above 6–7 kbar orthoenstatite is the only Ca-poor pyroxene polymorph which appears on the vapor-saturated liquidus surface.At pressures above approximately 10.8 kbar, only diopside, forsterite, and merwinite occur as primary liquidus phases in the system CaMgSi₂O₆-Mg₂SiO₄-Ca₃MgSi₂O₈-H₂O, in the presence of an aqueous vapor phase. At pressures between 1 atm and 10.2 kbar, both akermanite and monticellite also occur as primary liquidus phases. Comparison of the 1 atm and 10 kbar vapor-saturated liquidus diagrams suggests that melilite basalt bears a low pressure, or shallow depth, relationship to monticellite-bearing ultrabasites.     

The heat capacity of a natural monticellite and phase equilibria in the system CaO-MgO-SiO2-CO2

The heat capacity of a natural monticellite (Ca_1.00Mg_.09Fe_.91Mn_.01Si_0.99O_3.99) measured between 9.6 and 343 K using intermittent-heating, adiabatic calorimetry yields C_p⁰(298) and S₂₉₈⁰ of 123.64 ± 0.18 and 109.44 ± 0.16 J · mol⁻¹K⁻¹ respectively. Extrapolation of this entropy value to end-member monticellite results in an S⁰₂₉₈ = 108.1 ± 0.2 J · mol⁻¹K⁻¹. High-temperature heat-capacity data were measured between 340–1000 K with a differential scanning calorimeter. The high-temperature data were combined with the 290–350 K adiabatic values, extrapolated to 1700 K, and integrated to yield the following entropy equation for end-member monticellite (298–1700 K): S_T⁰(J · mol⁻¹K⁻¹) = S₂₉₈⁰ + 164.79 In T + 15.337 · 10⁻³T + 22.791 · 10⁵T⁻² − 968.94. Phase equilibria in the CaO-MgO-SiO₂ system were calculated from 973 to 1673 K and 0 to 12 kbar with these new data combined with existing data for akermanite (Ak), diopside (Di), forsterite (Fo), merwinite (Me) and wollastonite (Wo). The location of the calculated reactions involving the phases Mo and Fo is affected by their mutual solid solution. A best fit of the thermodynamically generated curves to all experiments is made when the S⁰₂₉₈ of Me is 250.2 J · mol⁻¹ K⁻¹ less than the measured value of 253.2 J · mol⁻¹ K⁻¹.A best fit to the reversals for the solid-solid and decarbonation reactions in the CaO-MgO-SiO₂-CO₂ system was obtained with the ΔG⁰₂₉₈ (kJ · mole⁻¹) for the phases Ak(−3667), Di(−3025), Fo(−2051), Me(−4317) and Mo(−2133). The two invariant points − Wo and −Fo for the solid-solid reactions are located at 1008 ± 5 K and 6.3 ± 0.1 kbar, and 1361 ± 10 K and 10.2 ± 0.2 kbar respectively. The location of the thermodynamically generated curves is in excellent agreement with most experimental data on decarbonation equilibria involving these phases.     

Holtedahlite, a new magnesium phosphate from Modum, Norway

Holtedahlite occurs in a serpentine-magnesite deposit at Modum in association with althausite and (OH,F) apatite. It is colourless, transparent, has a vitreous lustre, occurs in massive form and shows no cleavage. The Mohs' hardness is $\text{4}\text{1}\text{2}\text{?5}$; specific gravity 2.94(2) (Berman balance), calculated density 2.936 g/cm³. It is uniaxial negative, ω = 1.599(1), ? = 1.597(1).Chemical analysis with the electron microprobe and separate determinations of CO₂ and H₂O give a formula close to Mg₂PO₄OH or more specifically (Mg,Na)₂(PO₄, CO₃OH)(OH,F). Holtedahlite is hexagonal, $\text{α = 11.188(2), c = 4.975(1)}\text{A}\text{?}\text{, V = 539.3(3)}\text{A}\text{?}{}^3\text{Z = 6}$, space group $\text{P321, P3m}\text{l}\text{or}\text{P}\text{3}\text{m}\text{l}$. The strongest lines in the X-ray powder pattern are (in Å, with intensities and indices): 3.722(90) ($\text{11}\text{2}\text{1}$), 3.475(50) ($\text{20}\text{2}\text{1}$), 3.234(30) ($\text{30}\text{3}\text{0}$), 2.796(30), 2.438(100) ($\text{22}\text{4}\text{1}$), 2.177(30) (2.177(30) ($\text{40}\text{4}\text{1}$), 1.859(30) ($\text{22}\text{4}\text{2}$). The infrared spectrum shows the presence of OH^?, CO₃^2?, and PO₄^3? and indicates an O-H… O distance of ca. 3.02 Å.     

Thermochemistry of the high structural state plagioclases

The enthalpies of solution of a suite of 19 high-structural state synthetic plagioclases were measured in a Pb₂B₂O₅ melt at 970 K. The samples were crystallized from analyzed glasses at 1200°C and 20 kbar pressure in a piston-cylinder apparatus. A number of runs were also made on Amelia albite and Amelia albite synthetically disordered at 1050–1080°C and one bar for one month and at 1200°C and 20 kbar for 10 hr. The component oxides of anorthite, CaO, Al₂O₃ and SiO₂, were remeasured.The ΔH of disorder of albite inferred in the present study from albite crystallized from glass is 3.23 kcal, which agrees with the 3.4 found by Holm and Kleppa (1968). It is not certain whether this value includes the ΔH of a reversible displacive transition to monoclinic symmetry, as suggested by Helgesonet al. (1978) for the Holm-Kleppa results. The enthalpy of solution value for albite accepted for the solid solution series is based on the heat-treated Amelia albite and is 2.86 kcal less than for untreated Amelia albite.The enthalpy of formation from the oxides at 970 K of synthetic anorthite is ?24.06 ± 0.31 kcal, significantly higher than the ?23.16 kcal found by Charluet al. (1978), and in good agreement with the value of ?23.89 ± 0.82 given by Robieet al. (1979), based on acid calorimetry.The excess enthalpy of mixing in high plagioclase can be represented by the expression, valid at 970 K: ΔH^ex(±0.16 kcal) = 6.7461 X_abX²_An + 2.0247 X_AnX²_Ab where X_Ab and X_An are, respectively, the mole fractions of NaAlSi₃O₈ and CaAl₂Si₂O₈. This ΔH^ex, together with the mixing entropy of Kerrick and Darken's (1975) Al-avoidance model, reproduces almost perfectly the free energy of mixing found by Orville (1972) in aqueous cation-exchange experiments at 700°C. It is likely that Al-avoidance is the significant stabilizing factor in the high plagioclase series, at least for X_An≥ 0.3. At high temperatures the plagioclases have nearly the free energies of ideal one-site solid solutions. The Al-avoidance model leads to the following Gibbs energy of mixing for the high plagioclase series: $\text{ΔG}{}^{\mathrm{mix}}\text{= ΔH}{}^{\mathrm{ex}}\text{+ RT X}{}_{\mathrm{Ab}}\text{ln[X}{}^2{}_{\mathrm{Ab}}\text{(2 ? X}{}_{\mathrm{Ab}}\text{)]+ X}{}_{\mathrm{An}}\text{ln}\text{[X}{}_{\mathrm{An}}\text{(1+X}{}_{\mathrm{An}}\text{)}{}^2\text{]}\text{4}$. The entropy and enthalpy of mixing should be very nearly independent of temperature because of the unlikelihood of excess heat capacity in the albite-anorthite join.     

Sorption of noble gases by solids,with reference to meteorites. I. Magnetite and carbon

To simulate trapping of meteoritic noble gases by solids, 18 samples of Fe₃O₄ were synthesized in a noble gas atmosphere at 350–720 K by the reactions: 3Fe + 4H₂O → Fe₃O₄ + 4H₂ (Ne, Ar, Kr, Xe) 3Fe + 4CO → Fe₃O₃ + 4C + carbides (Xe only) Phases were separated by selective solvents (HgCl₂, HCl). Noble gas contents were analyzed by mass spectrometry, or, in runs where 36 d Xe¹²⁷ tracer was used, by γ-counting. Surface areas, as measured by the BET method, ranged from 1 to 400 m²/g. Isotopic fractionations were below the detection limit of 0.5%/m.u.Sorption of Xe on Fe₃O₄ and C obeys Henry's Law between $\text{1 × 10}{}^{\mathrm{?8}}$ and $\text{4 × 10}{}^{\mathrm{?5}}$ atm, but shows only a slight temperature dependence between 650 and 720 K ($\text{ΔH}{}_{\mathrm{sol}}\text{= ?4 ± 2}\text{kcal/mole}$). The mean distribution coefficient K_Xe is $\text{0.28 ± 0.09}$ cc STP/g atm for Fe₃O₄ and only a factor of $\text{1.2 ± 0.4}$ greater for C; such similarity for two cogenetic phases was predicted by Lewis et al. (1977). Stepped heating and etching experiments show that 20–50% of the total Xe is physically adsorbed and about 20% is trapped in the solid. The rest is chemisorbed with $\text{ΔH}{}_{\mathrm{s}}\text{? ?13}\text{kcal/mole}$. The desorption or exchange half-time for the last two components is >10² yr at room temperature.Etching experiments showed a possible analogy to “Phase $\text{Q}$” in meteorites. A typical carbon + carbide sample, when etched with HNO₃, lost 47% of its Xe but only 0.9% of its mass, corresponding to a ~0.6 Å layer. Though this etchable, surficial gas component was more thermolabile than $\text{Q}$ (release $\text{T}$ below 1000°C, compared to 1200–1600°C), another experiment shows that the proportion of chemisorbed Xe increases upon moderate heating (1 hr at 450°C). Apparently adsorbed gases can become “fixed” to the crystal, by processes not involving volume diffusion (recrystallization, chemical reaction, migration to traps, etc.). Such mechanisms may have acted in the solar nebula, to strengthen the binding of adsorbed gases.Adsorbed atmospheric noble gases are present in all samples, and dominate whenever the noble gas partial pressure in the atmosphere is greater than that in the synthesis. Many of the results of Lancet and Anders (1973) seem to have been dominated by such an atmospheric component; others are suspect for other reasons, whereas still others seem reliable. When the doubtful samples of Lancet and Anders are eliminated or corrected, the fractionation pattern—as in our samples—no longer peaks at Ar, but rises monotonically from Ne to Xe. No clear evidence remains for the strong temperature dependence claimed by these authors.     

Thermochemistry of some pyroxenes and related compounds

The enthalpies of formation of a number of crystalline silicates from the oxides at 986 K were determined by oxide melt solution calorimetry. The values of ΔH°_{f, 986}, in kcal/mol, are as follows: MgCaSi₂O₆, ? 34.3 ± 0.4; CoCaSi₂O₆, ? 26.7 ± 0.5; NiCaSi₂O₆, ? 27.1 ± 0.5; MnSiO₃, ? 6.3 ± 0.3; Mn₂SiO₄, ? 12.2 ± 0.3. In addition, for MnSiO₃ (rhodonite)→ MnSiO₃ (pyroxmangite), ΔH°₉₈₆ = + 0.06 ± 0.3₃kcal/mol and for MgCaSi₂O₆ (diopside) = MgCaSi₂O₆ (glass), ΔH°₉₈₆ = + 21.0 ± 0.3 kcal/ mol. For hedenbergite, FeCaSi₂O₆, ΔG°₁₃₅₀ = ?25.6 ± 1.5 kcal/mol. In terms of pyroxene phase equilibria and crystal chemistry, our thermochemical data support the generally accepted crystallographic arguments that (a) the C2/c clinopyroxene structure increases in stability with decreasing size of the ion occupying the Ml site in the MCaSi₂O₆ series, and (b) the energy (and enthalpy) differences between orthopyroxene, clinopyroxene, and pyroxenoid structures are generally quite small and often less than 500 cal/mol in magnitude.     

An evaluation of rate equations for calcite precipitation kinetics at pCO2 less than 0.01 atm and pH greater than 8

We used a reproducible seeded growth technique with a pH-stat to study the kinetics of calcite precipitation at 25°C. We performed different experiments at initial Ca²⁺ and HCO₃^? concentrations ranging from 0.7–2 and 4–7 mmol L^?1, pH values ranging from 8.25 to 8.70, $\text{pCO}{}_2$ values ranging from 0.0006 to 0.01 atm, and ionic strengths ranging from 0.015 to 0.10 mol L^?1. With this experimental data set, we used initial rate measurements and integral methods to test several precipitation rate equations. Rate equations that possess a disequilibrium functional dependence, such as the BURTON et al. (1951) dislocation model, forms of the Davies and Jones (1955) model, and the model used by Reddy and Nancollas (1973), did not adequately describe the kinetics of calcite precipitation at pH greater than 8 and $\text{pCO}{}_2$ less than 0.01 atm. Rate equations that describe independent dissolution and precipitation mechanisms with elementary reactions, such as the equation presented by Plummeret al. (1978), and nancollas and Reddy (1971) were more successful. However, Plummer's model did not adequately describe the rate of all experiments due to the presence of an OH^? surface term in the precipitation rate equation. The elementary reaction of the Nancollas and Reddy model is written in terms of bulk Ca²⁺ and CO₃^? concentrations, and appears to be the most successful model which describes calcite precipitation at pH > 8 and $\text{pCO}{}_2\text{< 0.01}$ atm. The Nancollas and Reddy model, altered to account for varying ionic strengths, adequately described the rate of all experiments and yielded a precipitation rate constant of $\text{118.2 ± 13.9}$ dm⁶ mol^?1 m^?2 s^?1, with an apparent Arrhenius activation energy of 48.1 kJ mol^?1.     

Lead orthophosphates—II. Stability of cholopyromophite at 25°c

The conversion of secondary lead orthophosphate [PbHPO₄] into chloropyromorphite [Pb₅(PO₄)₃Cl] in ca. 10^?1 M NaCl solutions has been investigated at 25°C. From the composition of the supernatant solutions, the solubility product constant for Pb₅(PO₄)₃Cl has been calculated to be 10^?84.4±0.1, corresponding to $\text{ΔG}{}_{\mathrm{?}}\text{°}$ of ?906.2 kcal mol^?1. The solution equilibria and phase relationships in the system PbCl₂-PbO-P₂O₈-H₂O are discussed along with the geological implications.     

Thermochemical study of glasses in the system NaAlSi3O8-KAlSi3O8-Si4O8 and the join Na1.6Al1.6Si2.4O8-K1.6Al1.6Si2.4O8

Enthalpies of solution in 2PbO · B₂O₃ at 981 K have been measured for glasses in the system albite-orthoclase-silica and along the join Na_1.6Al_1.6Si_2.4O₈-K_1.6Al_1.6Si_2.4O₈. The join KAlSi₃O₈-Si₄O₈ shows zero heat of mixing similar to that found previously for NaAlSi₃O₈-Si₄O₈ glasses. Albite-orthoclase glasses show negative heats of mixing symmetric about Ab₅₀Or₅₀ (W_n = ? 2.4 ± 0.8 kcal). Negative heats of (Na, K) mixing are also found at $\text{Si}\text{(Si + Al)}\text{= 0.6.}$ Ternary excess enthalpies of mixing in the glassy system Ab-Or-4Q are positive but rarely exceed 1 kcal mol^?1.Using earlier studies of the thermodynamic properties of the crystals, the present calorimetric data and the “two-lattice” entropy model, the albite-orthoclase phase diagram is calculated in good agreement with experimental data. Attempts to calculate albite-silica and orthoclase-silica phase diagrams reveal complexities probably related to significant (but unknown) mutual solid solubility between cristobalite and alkali feldspar and to the very small heat and entropy of fusion of SiO₂.     

High pressure phase transformations in the joins Mg2SiO4-Ca2SiO4 and MgO-CaSiO3

High pressure phase transformations for all the mineral phases along the joins Mg₂SiO₄-Ca₂-SiO₄ and MgO-CaSiO₃ in the system MgO-CaO-SiO₂ were investigated in the pressure range between 100 and 300 kbar at about 1,000 °C, by means of the technique involving a diamond-anvil press coupled with laser heating. In addition to the four end-members, there are three stable intermediate mineral components in these two joins. Phase behaviour of all the end-member components at high pressure have been reported earlier and are reviewed here. Results of this study reveal that the three intermediate components are all unstable relative to the end-members at pressures greater than 200 kbar. Ultimately, monticellite (CaMgSiO₄) decomposes into CaSiO₃ (perovskite-type)+MgO; merwinite (Ca₃MgSi₂O₈) decomposes into Ca₂SiO₄(K₂NiF₄-type)+CaSiO₃ (perovskite-type)+MgO; and akermanite (Ca₂MgSi₂O₇) decomposes into CaSiO₃ (perovskite-type)+MgO. Note that the decomposition reactions of all phases studied here result in the formation of MgO. Intermediate Ca-Mg silicates transform to pure Ca-silicates plus MgO, while pure Mg₂SiO₄ transforms to MgSiO₃+MgO.     

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.     

Effect of pressure on the diffusivity of network-forming cations in melts of jadeitic compositions

The diffusivities of network-forming cations (Si⁴⁺, Al³⁺, Ge⁴⁺ and Ga³⁺) in melts of the jadeitic composition NaAl(Si, Ge)₂O₆ and Na(Al, Ga)Si₂O₆ have been measured at pressures between 6 and 20 kbar at 1400°C. The rates of interdiffusion of Si⁴⁺-Ge⁴⁺ and Al³⁺-Ge³⁺ increase with increasing pressure at constant temperature. The results are consistent with the ion-dynamics computer simulations of Jadeite melt by Angellet al. (1982, 1983). The coefficient measured for the Si⁴⁺-Ge⁴⁺ interdiffusion is between 8 × 10^?10 and 2.5 × 10^?8$\text{cm}{}^2\text{sec}$ at 6 kbar, depending on the composition of the melt, whereas at 20 kbar it is between 7 × 10^?9 and 2 × 10^?7$\text{cm}{}^2\text{sec}$. The effect of pressure is greater for more Si-rich compositions (i.e., closer to NaAlSi₂O₆ composition). The coefficient measured for the Al³⁺-Ga³⁺ inter- diffusion is between 9 × 10^?10 and 3 × 10^?9 cm²/sec at 6 kbar and between 3 × 10^?9 and 1 × 10^?8$\text{cm}{}^2\text{sec}$ at 20 kbar. The rate of increase in diffusivity with pressure of Al³⁺-Ga³⁺ (a factor of 3–4) is smaller than that of Si⁴⁺-Ge⁴⁺ (a factor of 7–17).The Si⁴⁺-Ge⁴⁺ interdiffusion in melts of Na₂O · 4(Si, Ge)O₂ composition has also been measured at 8 and 15 kbar for comparison. The effect of pressure on the diffusivity in this melt is significantly smaller than that for the jadeitic melts. The increase in diffusivity of the network-forming cations in jadeitic melts with increasing pressure may be related to the decrease in viscosity of the same melt. The present results, as well as the ion-dynamics simulations, suggest that the homogenization of partial melts and mixing of magmas would be more efficient at greater depths.     

Geochronologic investigations of the Sudbury Nickel Irruptive and the Superior Province granites north of Sudbury

RbSr ($\text{λ}\text{Rb}\text{= 1.39 × 10}{}^{\mathrm{?11}}\text{yr}{}^{\mathrm{?1}}\text{)}$ and U-Pb ($\text{λ 238 = 1.54 × 10}{}^{\mathrm{?10}}\text{yr}{}^{\mathrm{?1}}\text{, λ235 = 9.72 × 10}{}^{\mathrm{?10}}\text{yr}{}^{\mathrm{?1}}$) measurements were undertaken in the Sudbury area, Sudbury, Ontario to determine the ages of the Sudbury Nickel Irruptive, Superior Province granites north of Sudbury, Sudbury Breccia and subsequent metamorphism. The Sudbury Nickel Irruptive norite whole rock Rb-Sr data yield an age of $\text{1883 ± 136}\text{Myr}$ ($\text{I.R}\text{. = 0.7071 ± 0.0005}$; all results quoted at 2π level) while the Nickel Irruptive micropegmatite Rb-Sr system has been disturbed and does not yield an isochron. A plagioclase-whole rock pair from the norite near the norite-micropegmatite transition yields an age of 1866 Myr, which when taken in conjunction with field (Stevenson and Colgrove, 1968) and geochemical (Naldrettet al., 1970, 1972) data does not support the conclusion of gibbins and McNurr (1972) that the micropegmatite is a later intrusion rather than a differentiate of the magma which produced the norite. Rb-Sr studies of the Superior Province granites north of Sudbury yield an age of $\text{2698 ± 162}\text{Myr}$ ($\text{I.R.}\text{= 0.7019 ± 0.0012}$). U-Pb zircon studies of these granites and granitic clasts within the Sudbury Breccia yield an age of $\text{2.71 ± 0.05}$ Byr and suggest the breccia granitic clasts were derived from the Superior Province granites. The granitic rocks ~150 km north of Sudbury have been undisturbed for ~ 2.6 Byr based on Rb-Sr mineral studies, whereas the granites and Sudbury Breccia within ~ 15 km of the Nickel Irruptive, as well as the Sudbury norite at the perimeter of the Irruptive have been disturbed by the Penokean Orogeny 1.7–1.75 Byr ago. The Penokean event appears to have overprinted isotopic evidence of the Sudbury impact event at least in the area studied.     

Phase relations to 3 kbar in the systems edenite + H2O and edenite + excess quartz + H2O

Amphiboles approached edenite (NaCa₂Mg₅Si₇AlO₂₂(OH)₂), richterite (Na₂CaMg₅Si₈O₂₂(OH)₂), tremolite (□Ca₂Mg₅Si₈O₂₂(OH)₂) solid solutions were studied by conventional hydrothermal techniques employing the bulk compositions edenite, and edenite + additional quartz, all with excess H₂O. For the stoichiometric edenite bulk composition + excess H₂O, the equilibrium phase assemblage is diopside + Na-phlogopite + forsterite + fluid at, and just above the amphibole high-temperature limit at 850 ± 5°C, 500 bar, and 880 ± 5°C, 1000 bar. The breakdown temperature of sodic phlogopite is 855 ± 3°C at 500 bar, and 890 ± 5°C at 700 bar, producing nepheline + plagioclase (or melt), additional forsterite and fluid. Diopside and Na-phlogopite solid solution coexist over a broad P_fluid-T region, even within the amphibole field, where they are associated with an edenite-richterite (-tremolite) solid solution of approximate composition Ed₃₅Rc₅₀Tr₁₅.In the system edenite + 4 quartz + excess H₂O, nearly pure tremolite and albite coexist stably between 670° and 830°C at 1000 bar and give way to the possibly metastable assemblage diopside + talc + albite below 670°C. In the presence of albite, tremolite reacts to produce diopside + quartz + enstatite + fluid above 830°C at 1000 bar. For the investigated silica-rich bulk composition, amphibole P_fluid-T stability is divided by the albite melting curve into a tremolite + albite field, and a tremolite + aqueous melt field. Substantial equilibrium solid solution of tremolite towards edenite or richterite was not observed for silica-excess bulk compositions. Metastable edenite-rich amphiboles initially synthesized change to tremolite with increasing run length in the presence of free SiO₂.Edenitic amphibole is stable only over a very limited temperature range in silica-undersaturated environments, thus accounting for its rarity in nature. Na-phlogopite solid solutions are also disfavored by high a_SiO2; even for nepheline-normative lithologies, a hypothesized rapid low-temperature conversion to vermiculite or smectite could partly explain the scarcity of sodic phlogopite in rocks.     

Volatilization of sodium from silicate melt spheres and its application to the formation of chondrules

The rates of volatilization of Na from liquid spheres of chondrule compositions have been determined as functions of time, temperature, partial pressure of oxygen, and sizes of the spheres. The Na₂O content in the sphere is uniform in each run. but it decreases with time of the run, indicating that the rate of diffusion of Na in the liquid is greater than that of volatilization, and that the latter is the rate-controlling process. The rate of sodium volatilization becomes greater with increasing temperature and with decreasing P_O2 and size of the spheres. The relation of the Na₂O content in the liquid sphere with time and its size indicate that the amount of Na₂O volatilized from the liquid spheres within unit time is proportional to the surface area of the spheres and the concentration of Na₂O in the liquid. From these relations, the rate of volatilization of sodium can be obtained at constant temperature and P_o2. The rate of volatilization of sodium satisfies the Arrhenius relation within the temperature range from about 1450–1600 C at 10^?9,2 atm p_O2; the activation energy for the sodium volatilization is approximately 100 kcal-mole^?1. The rate is also approximately proportional to within the range of p_O2 from 10^?10.2 to 10^?5.0 atm at about 1500° C. Based on the present results and the Na₂O contents in chondrules. it is suggested that they experienced an instant heating with maximum temperature of 1400–2200° C followed by an immediate cooling.