Atmospheric pressure changes due to volcanic eruptions and possible water level fluctuations

The data of Reed (1983) are analysed to produce the following empirical equations for the amplitude p ₀ (overall fluctuation) in Pascals of the air pressure wave associated with a volcanic eruption of volume V km³ or a nuclear explosion of strength M Mt: Here s is the distance from the source in km. $$\begin{gathered} \log _{10} p_0 = 4.44 + \log _{10} V - 0.84\log _{10} s \hfill \\ {\text{ }} = 3.44 + \log _{10} M - 0.84\log _{10} s. \hfill \\ \end{gathered} $$ Garrett's (1970) theory is examined on the generation of water level fluctuations by an air pressure wave crossing a water depth discontinuity such as a continental shelf. The total amplitude of the ocean wave is determined to be where c ₂ ¹ = gh ₁, c ₂ ² = gh ₂, g is acceleration of gravity, h ₁ and h ₂ are the water depths on the ocean and shore side of the depth discontinuity, c is the speed of propagation of the air pressure wave, and ? is the water density. $$B = \left[ {\frac{{c_2^2 }}{{c^2 - c_2^2 }} + \frac{{c^2 (c_1 - c_2 )}}{{(c - c_1 )(c^2 - c_2^2 )}}} \right]\frac{{p_0 }}{{g\varrho }}$$ It is calculated that a 10 km³ eruption at Mount St. Augustine would cause a 460 Pa air pressure wave and a discernible water level fluctuation at Vancouver Island of several cm amplitude.     

Carbon Dioxide in igneous petrogenesis: I

A number of experimental CO₂ solubility data for silicate and aluminosilicate melts at a variety of P- T conditions are consistent with solution of CO₂ in the melt by polymer condensation reactions such as SiO _4(m ^4? +CO_2(v)+Si _n O _3n+1(m) ⁽²ⁿ⁺¹⁾ ?Si _n+1O _3n+4(m) ^(2n+4)? +CO _3(m ^)2? . For various metalsilicate systems the relative solubility of CO₂ should depend markedly on the relative Gibbs free change of reaction. Experimental solubility data for the systems Li₂O-SiO₂, Na₂O-SiO₂, K₂O-SiO₂, CaO-SiO₂, MgO-SiO₂ and other aluminosilicate melts are in complete accord with predictions based on Gibbs Free energies of model polycondesation reactions. A rigorous thermodynamic treatment of published P- T-wt.% CO₂ solubility data for a number of mineral and natural melts suggests that for the reaction CO_2(m) ? CO_2(v)

1. CO₂-melt mixing may be considered ideal (i.e., { \(a_{{\text{CO}}_{\text{2}} }^m = X_{{\text{CO}}_{\text{2}} }^m \) );
2. \(\bar V_{{\text{CO}}_{\text{2}} }^m \) , the partial molal volume of CO₂ in the melt, is approximately equal to 30 cm³ mole^?1 and independent of P and T;
3. Δ C _p ⁰ is approximately equal to zero in the T range 1,400° to 1,650 °C and
4. enthalpies and entropies of the dissolution reaction depend on the ratio of network modifiers to network builders in the melt. Analytic expressions which relate the CO₂ content of a melt to P, T, and \(f_{{\text{CO}}_{\text{2}} } \) for andesite, tholeiite and olivine melilite melts of the form

$$\ln X_{{\text{CO}}_{\text{2}} }^m = \ln f_{{\text{CO}}_{\text{2}} } - \frac{A}{T} - B - \frac{C}{T}(P - 1)$$ have been determined. Regression parameters are (A, B, C): andesite (3.419, 11.164, 0.408), tholeiite (14.040, 5.440,0.393), melilite (9.226, 7.860, 0.352). The solubility equations are believed to be accurate in the range 3<P<30 kbar and 1,100°<T<1,650 °C. A series of CO₂ isopleth diagrams for a wide range of T and P are drawn for andesitic, tholeiitic and alkalic melts.     

Experimental study of the physico-chemical conditions of skarn formation

Experimental results show that skarns can—800°C and 500–1,000 bars. Starting materials include intermediate-acidic igneous rocks, volcanic rocks, metamorphic rocks, carbonates of various purities and chemical reagents of analytical purity-grade. Experimental media are: NaCl, NaCl+CaCl₂, NaCl+CaCl₂+MgCl₂, Na₂CO₃ and Na₂SiO₃ solutions. Experimental results show that skarns can be formed under wide physico-chemical conditions:T=400°–800°C,P=500–1,000 bars,pH=4–11, and \(f_{O_2 } = 10^{ - 23} - 10^{ - 11} \) bar. The mineralogy of skarns and the chemical compositions of skarn minerals are generally controlled by the combined factors: the chemical composition of the original rocks, pH values, redox conditions, temperatures and pressures. Isomorphous substitution may have agreat effect on the temperature of formation andfo ₂ of some major skarn minerals. It is found that skarnization occurs preferentially in NaCl and NaCl+CaCl₂ solutions and subortinately in MaCO₃ and Na₂SiO₃ soiutions.     

Mica-feldspar equilibria in supercritical alkali chloride solutions

The equilibrium position of the reaction $$\begin{gathered} 1.5 KAlSi_3 O_8 + HCl = 0.5 KAl_3 Si_3 O_{10} (OH)_2 \hfill \\ + 3SiO_2 + KCl \hfill \\ \end{gathered} $$ has been located at 1 and 2 kb pressure and temperatures between 600° and 670° C using the Ag-AgCl buffer. These data can be combined with information on the dissociation of KC1, HC1 and H₂O to determine species abundances in supercritical aqueous fluids in equilibrium with muscovite — K-feldspar — quartz assemblages. Chloride species become increasingly associated with increasingT, increasing total molality, (m _tot or \(m_{Cl_{tot} } \) ), and decreasing \(P_{H_2 O} \) . Master variable diagrams indicate that the pH of the solutions may vary from near neutral to quite acid. Published data on the paragonite-albite-quartz reaction and exchange reactions involving feldspars and micas were included to calculate speciation in mica-feldspar-NaCl-KCl-HCl-H₂O fluids at 2kb pressure and temperatures between 300° and 600° C. The data are not accurate enough to distinguish different feldspar structural states. Concentration gradients were calculated for individual species between K-feldspar+quartz, muscovite+quartz and andalusite+quartz assemblages at 500° C, 2 kb. Assuming that the proton diffuses most rapidly and that there are no [H⁺] gradients, the molality of the solution must vary 30-fold, with feldspar+quartz at the more concentrated side. The data on mica-feldspar-chloride equilibria are used to interpret the spacial distribution of micas, feldspar and quartz in microfolds. This distribution can be accounted for by pressure solution, due to the fact that non-hydrostatic pressure affects congruently dissolving minerals, auch as quartz, differently from minerals which dissolve incongruently, such as micas and feldspars. We postulate, that during folding at constant \(P_{H_2 O} \) ,T and \(m_{Cl - } \) , gradients in KC1 and SiO₂ are created by stress differences between hinge and limb of a microfold, such that both migrate to the hinge area where quartz precipitates and muscovite is converted to K-felspar, thus accounting for the observed mineral distribution.     

On the Relationship between Refractive Index and Density for SiO2-polymorphs

Five different refraction formulas were applied to SiO₂ polymorphs in order to determine the most suitable refractive index-density relation. 13 SiO₂ polymorphs with topological different tetrahedral frameworks are used in this study including eight new low density SiO₂ polymorphs — so called “guest free porosils”. These SiO₂ polymorphs cover a density range from 1.76 to 2.92 g/cm³. The mean refractive indices (ovn) of the porosils have been determined by the immersion method, the densities (ρ) were calculated from the unit cell parameters. Assuming the polarizability (α) of all SiO₂ polymorphs to be constant the general refractivity formula $$\{ 2\overline {11} 0\} \langle 0001\rangle $$ turned out to be the most suitable for SiO₂ polymorphs. Regression analysis yields an electronic overlap parameter b=1.2(1).     

The reversed alumina contents of orthopyroxene in equilibrium with spinel and forsterite in the system MgO-Al2O3-SiO2

Equilibrium alumina contents of orthopyroxene coexisting with spinel and forsterite in the system MgO-Al₂O₃-SiO₂ have been reversed at 15 different P-T conditions, in the range 1,030–1,600° C and 10–28 kbar. The present data and three reversals of Danckwerth and Newton (1978) have been modeled assuming an ideal pyroxene solid solution with components Mg₂Si₂O₆ (En) and MgAl₂SiO₆ (MgTs), to yield the following equilibrium condition (J, bar, K): $$\begin{gathered} RT{\text{ln(}}X_{{\text{MgTs}}} {\text{/}}X_{{\text{En}}} {\text{) + 29,190}} - {\text{13}}{\text{.42 }}T + 0.18{\text{ }}T + 0.18{\text{ }}T^{1.5} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [0.013 + 3.34 \times 10^{ - 5} (T - 298) - 6.6 \times 10^{ - 7} P]P. \hfill \\ \end{gathered} $$ The data of Perkins et al. (1981) for the equilibrium of orthopyroxene with pyrope have been similarly fitted with the result: $$\begin{gathered} - RT{\text{ln(}}X_{{\text{MgTs}}} \cdot X_{{\text{En}}} {\text{) + 5,510}} - 88.91{\text{ }}T + 19{\text{ }}T^{1.2} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [ - 0.832 - 8.78{\text{ }} \times {\text{ 10}}^{ - {\text{5}}} (T - 298) + 16.6{\text{ }} \times {\text{ 10}}^{ - 7} P]{\text{ }}P. \hfill \\ \end{gathered} $$ The new parameters are in excellent agreement with measured thermochemical data and give the following properties of the Mg-Tschermak endmember: $$H_{f,970}^0 = - 4.77{\text{ kJ/mol, }}S_{298}^0 = 129.44{\text{ J/mol}} \cdot {\text{K,}}$$ and $$V_{298,1}^0 = 58.88{\text{ cm}}^{\text{3}} .$$ The assemblage orthopyroxene+spinel+olivine can be used as a geothermometer for spinel lherzolites, subject to a choice of thermodynamic mixing models for multicomponent orthopyroxene and spinel. An ideal two-site mixing model for pyroxene and Sack's (1982) expressions for spinel activities provide, with the present experimental calibration, a geothermometer which yields temperatures of 800° C to 1,350° C for various alpine peridotites and 850° C to 1,130° C for various volcanic inclusions of upper mantle origin.     

The lattice energy of forsterite. Charge distribution and formation enthalpy of the SiO 4 4− ion

In the lattice energy expression of forsterite, based on a Born-Mayer (electrostatic+repulsive+dispersive) potential, the oxygen charge z _o, the hardness parameter ρ and the repulsive radii r _Mg and r _Si appear as unknown parameters. These were determined by calculating the first and second partial derivatives of the energy with respect to the cell edges, and equalizing them to quantities related to the crystal elastic constants; the overdetermined system of equations was solved numerically, minimizing the root-mean-square deviation. To test the results obtained, the SiO ₄ ^4? ion was assumed to move in the unit-cell, and the least-energy configuration was sought and compared with the experimental one. By combining the two methods, the optimum set of parameters was: z _o=?1.34, ρ=0.27 Å, r _Mg=0.72 Å, r _Si=0.64 Å. The values ?8565.12 and ?8927.28 kJ mol^?1 were obtained, respectively, for the lattice energy E _Land for its ionic component E _L ⁰ ,which accounts for interactions between Mg²⁺ and SiO ₄ ^4? ions only. The charge distribution calculated on the SiO ₄ ^4? ion was discussed and compared with other results. Using appropriate thermochemical cycles, the formation enthalpy and the binding energy of SiO ₄ ^4? were estimated to be: ΔH _f(SiO ₄ ^4? )=2117.6 and E(SiO ₄ ^4? )=708.6 kJ mol^?1, respectively.     

Experimental study of subsolidus phase relations and mixing properties of pyroxene and plagioclase in the system Na2O-CaO-Al2O3-SiO2

In the system Na₂O-CaO-Al₂O₃-SiO₂ (NCAS), the equilibrium compositions of pyroxene coexisting with grossular and corundum were experimentally determined at 40 different P-T conditions (1,100–1,400° C and 20.5–38 kbar). Mixing properties of the Ca-Tschermak — Jadeite pyroxene inferred from the data are (J, K): $$\begin{gathered} G_{Px}^{xs} = X_{{\text{CaTs}}} X_{{\text{Jd}}} [14,810 - 7.15T - 5,070(X_{{\text{CaTs}}} - X_{{\text{Jd}}} ) \hfill \\ {\text{ }} - 3,350(X_{{\text{CaTs}}} - X_{{\text{Jd}}} )^2 ] \hfill \\ \end{gathered} $$ The excess entropy is consistent with a complete disorder of cations in the M2 and the T site. Compositions of coexisting pyroxene and plagioclase were obtained in 11 experiments at 1,190–1,300° C/25 kbar. The data were used to infer an entropy difference between low and high anorthite at 1,200° C, corresponding to the enthalpy difference of 9.6 kJ/mol associated with the C \(\bar 1\) =I \(\bar 1\) transition in anorthite as given by Carpenter and McConnell (1984). The resulting entropy difference of 5.0 J/ mol · K places the transition at 1,647° C. Plagioclase is modeled as ideal solutions, C \(\bar 1\) and I \(\bar 1\) , with a non-first order transition between them approximated by an empirical expression (J, bar, K): $$\Delta G_T = \Delta G_{1,473} \left[ {1 - 3X_{Ab} \tfrac{{T^4 - 1,473^4 }}{{\left( {1,920 - 0.004P} \right)^4 - 1,473^4 }}} \right],$$ where $$\Delta G_{1,473} = 9,600 - 5.0T - 0.02P$$ The derived mixing properties of the pyroxene and plagioclase solutions, combined with the thermodynamic properties of other phases, were used to calculate phase relations in the NCAS system. Equilibria involving pyroxene+plagioclase +grossular+corundum and pyroxene+plagioclase +grossular+kyani te are suitable for thermobarometry. Albite is the most stable plagioclase.     

On the nonstoichiometry and point defects of olivine

This paper presents the point-defect thermodynamics for fayalite and olivine solid solutions (Fe _x Mg_1?x )₂SiO₄. By means of thermogravimetry, the metal-to-oxygen ratio of these silicates has been determined as a function of oxygen potential, compositionx and temperature. Experiments were performed in the range of 1,000° C≦T≦1,280° C and 0.2≦x≦1.0. It is found that V _Me ^″ , Fe _Me ^· and the associate {Fe′ _Si ^′ Fe _Me ^· } are the majority defects. With this knowledge it is possible to calculate the nonstoichiometry at given temperature as a function of \(p_{O_2 } \) and \(a_{SiO_2 } \) . The cation vacancy concentration shows a \(p_{O_2 }^{1/5} \) -dependence (forx≧0.2) and increases at givenT and \(p_{O_2 } \) almost exponentially with compositionx. In the composition range studied here, the silicates show an oxygen excess, and FeO is more soluble in the olivine than SiO₂.     

Use of FePt alloys to eliminate the iron loss problem in 1 atmosphere gas mixing experiments: Theoretical and practical considerations

Theoretical and practical considerations are combined to place limits on the iron content of an FePt alloy that is in equilibrium with silicate melt, olivine and a gas phase of known \(f_{{\text{O}}_{\text{2}} }\) . Equilibrium constants are calculated for the reactions: (1) $$2{\text{Fe}}^{\text{o}} + {\text{SiO}}_{\text{2}} + {\text{O}}_{\text{2}} \rightleftharpoons {\text{Fe}}_{\text{2}} {\text{SiO}}_{\text{4}}$$ (2) $${\text{Fe}}^{\text{o}} + \frac{1}{2}{\text{O}}_{\text{2}} \rightleftharpoons {\text{FeO}}$$ . These equilibria may be used to choose an appropriate iron activity for the FePt alloy of an experiment. The temperature dependence of the equilibrium constants is calculated from experimental data. The Gibbs free energy of reaction (1) obtained using thermochemical data is in close agreement with ΔG_rxn calculated from the experimental data. Reaction (1) has the advantage that it is independent of the Fe²⁺/Fe³⁺ ratio of the melt, but is limited to applications where olivine is a crystallizing phase and requires a formulation for \(a_{{\text{SiO}}_{\text{2}} }^{{\text{liq}}}\) . Reaction (2) uses an empirical approximation for the FeO/Fe₂O₃ ratio of the liquid, and is independent of olivine saturation. However, it requires a formulation for a _FeO ^liq . Either equilibrium constant may be used to calculate the appropriate FePt alloy in equilibrium with a silicate melt. If experiments are conducted at an \(f_{{\text{O}}_{\text{2}} }\) parallel that of a buffer assemblage, a small range of FePt alloys may be used over a large temperature interval. For example, an alloy containing from 6 % to 9 % Fe by weight is in equilibrium with olivine-saturated tholeiites and komatiites at the quartzfayalite-magnetite buffer over the temperature interval 1,400° C to 1,100° C. Lunar basalt liquids in equilibrium with olivine at 1/2 log unit below the iron-wüstite buffer require an FePt alloy that contains 30–50 wt. % iron over a similar temperature interval.     

Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

An experimental study of the partitioning of Fe and Mg between garnet and orthopyroxene

The partitioning of Fe and Mg between garnet and aluminous orthopyroxene has been experimentally investigated in the pressure-temperature range 5–30 kbar and 800–1,200° C in the FeO-MgO-Al₂O₃-SiO₂ (FMAS) and CaO-FeO-MgO-Al₂O₃-SiO₂ (CFMAS) systems. Within the errors of the experimental data, orthopyroxene can be regarded as macroscopically ideal. The effects of Calcium on Fe-Mg partitioning between garnet and orthopyroxene can be attributed to non-ideal Ca-Mg interactions in the garnet, described by the interaction term:W _CaMg ^ga -W _CaFe ^ga =1,400±500 cal/mol site. Reduction of the experimental data, combined with molar volume data for the end-member phases, permits the calibration of a geothermometer which is applicable to garnet peridotites and granulites: $$T(^\circ C) = \left\{ {\frac{{3,740 + 1,400X_{gr}^{ga} + 22.86P(kb)}}{{R\ln K_D + 1.96}}} \right\} - 273$$ with $$K_D = {{\left\{ {\frac{{Fe}}{{Mg}}} \right\}^{ga} } \mathord{\left/ {\vphantom {{\left\{ {\frac{{Fe}}{{Mg}}} \right\}^{ga} } {\left\{ {\frac{{Fe}}{{Mg}}} \right\}}}} \right. \kern-\nulldelimiterspace} {\left\{ {\frac{{Fe}}{{Mg}}} \right\}}}$$ and $$X_{gr}^{ga} = (Ca/Ca + Mg + Fe)^{ga} .$$ The accuracy and precision of this geothermometer are limited by largerelative errors in the experimental and natural-rock data and by the modest absolute variation inK _D with temperature. Nevertheless, the geothermometer is shown to yield reasonable temperature estimates for a variety of natural samples.     

Clinopyroxene on the join CaMgSi2O6-CaFe3+AlSiO6-CaTiAl2O6 at low oxygen fugacity

Subsolidus phase relations on the join CaMgSi₂O₆-CaFe³⁺ AlSiO₆-CaTiAl₂O₆ were studied by the ordinary quenching method at \(f_{O_2 } = 10^{ - 11} \) atm and 1,100°C. Crystalline phases encountered are clinopyroxene_ss (ss:solid solution) (Cpx_ss), melilite (Mel), perovskite (Pv), spinel_ss (Sp_ss), magnetite_ss (Mt_ss) and anorthite (An). There is no Cpx_ss single phase field, and the following assemblages were found; Cpx_ss+Mel, Cpx_ss+Mel+Sp_ss, Cpx_ss+Mel+Pv, Cpx_ss+Mel+Sp_ss+Pv, Cpx_ss+Pv+Sp_ss+An, Sp_ss+Pv+Mel+An+Cpx_ss, Mel+Mt_ss+An+Sp_ss+Cpx_ss+liquid and Mel+Mt_ss+An+Sp_ss+Cpx_ss+Pv. Mössbauer spectral study revealed that Cpx_ss contains both Fe²⁺ and Fe³⁺ in the octahedral site, and it was confirmed that the CaFe³⁺ AlSiO₆ content in the Cpx_ss at low \(f_{O_2 } \) is considerably less than that in the Cpx_ss crystallized in air, whereas the CaFe²⁺Si₂O₆ component increases. The maximum solubility of CaTlAl₂O₆ component in the Cpx_ss at low \(f_{O_2 } \) is higher than that in air. The decrease of CaFe³⁺ AlSiO₆ in the Cpx_ss at low \(f_{O_2 } \) may cause increase of CaTial₂O₆ in the Cpx_ss.     

Thermodynamic properties of andradite and application to skarn with coexisting andradite and hedenbergite

The enthalpy of formation of andradite (Ca₃Fe₂Si₃O₁₂) has been estimated as-5,769.700 (±5) kJ/mol from a consideration of the calorimetric data on entropy (316.4 J/mol K) and of the experimental phaseequilibrium data on the reactions: 1 $$\begin{gathered} 9/2 CaFeSi_2 O_6 + O_2 = 3/2 Ca_3 Fe_2 Si_3 O_{12} + 1/2 Fe_3 O_4 + 9/2 SiO_2 (a) \hfill \\ Hedenbergite andradite magnetite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 4 CaFeSi_2 O_6 + 2 CaSiO_3 + O_2 = 2 Ca_3 Fe_2 Si_3 O_{12} + 4 SiO_2 (b) \hfill \\ Hedenbergite wollastonite andradite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 18 CaSiO_3 + 4 Fe_3 O_4 + O_2 = 6Ca_3 Fe_2 Si_3 O_{12} (c) \hfill \\ Wollastonite magnetite andradite \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} Ca_3 Fe_2 Si_3 O_{12} = 3 CaSiO_3 + Fe_2 O_3 . (d) \hfill \\ Andradite pseudowollastonite hematite \hfill \\ \end{gathered} $$ and $$log f_{O_2 } = E + A + B/T + D(P - 1)/T + C log f_{O_2 } .$$ Oxygen-barometric scales are presented as follows: $$\begin{gathered} E = 12.51; D = 0.078; \hfill \\ A = 3 log X_{Ad} - 4.5 log X_{Hd} ; C = 0; \hfill \\ B = - 27,576 - 1,007(1 - X_{Ad} )^2 - 1,476(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite (Ad)-hedenbergite (Hd)-magnetite-quartz: $$\begin{gathered} E = 13.98; D = 0.0081; \hfill \\ A = 4 log(X_{Ad} / X_{Hd} ); C = 0; \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-quartz: 1 $$\begin{gathered} E = 13.98;{\text{ }}D = 0.0081; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 0;}} \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-calcitequartz: 1 $$\begin{gathered} E = - 1.69;{\text{ }}D = - 0.199; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 2;}} \hfill \\ B = - 20,441 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-calcite: 1 $$\begin{gathered} E = - 17.36;{\text{ }}D = - 0.403; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 4;}} \hfill \\ B = - 11,720 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 \hfill \\ \end{gathered} $$ The oxygen fugacity of formation of those skarns where andradite and hedenbergite assemblage is typical can be calculated by using the above equations. The oxygen fugacity of formation of this kind of skarn ranges between carbon dioxide/graphite and hematite/magnetite buffers. It increases from the inside zones to the outside zones, and appears to decrease with the ore-types in the order Cu, Pb?Zn, Fe, Mo, W(Sn) ore deposits.     

Jinshajiangite—a new Ba-Mn-Fe-Ti-bearing silicate mineral

Jinshajiangite, (Na, K)₅(Ba, Ca)₄(Fe²⁺, Mn)₁₅(Ti, Fe³⁺, Nb)₈(SiO₄)₁₅ (F, O, OH)₁₀, is a new mineral which occurs as thin tabular crystals in an arfedsonite-rich albite dyke near the Jinshajiang River winding through Sichuan Province, China. Jinshajiangite is monoclinic. Possible space groupC2/m, Cm orC2. Unit cell:a=10.732,b=13.847,c=20.817 Å,β=95°3′,z=2. The strongest lines in the X-ray diffraction pattern [d in (hkl)] are: 10.2 (7) (002), 3.44 (10) ( \(\bar 311\) , 310, \(\bar 205\) , 006), 3.15 (8) (205), 2.630 (8) (136, 243), 2.570 (8) ( \(\bar 430\) ), 2.202 (4) (405), 1.755 (4) ( \(\bar 536\) ), and 1.715 (5b) (3. 1. 10). The new mineral is black red brownish red or golden red in color. Lustervitreous. Streak light yellow. Cleavges {010} and {100} perfect. H. (Vicker) 430 kg/mm². Specific gravity 3.61 (meas.). Density 3.56 (cale.). It is optically biaxial positive, 2V=72,r<v; refractive indices:N _x =1.792,N _y =1.801,N _z =1.852. Oblique extinction anglec?X=13°, sign positive. The new mineral is strongly pleochroic:X=libht golden yellow,Y=brownish yellow,Z=brownish red. The DTA shows two endothermic peaks at about 795°C and 995°C. Infrared spectrum absorp- tion bands are observed at 308, 380, 492, 520, 620, 965, and 1,000 cm^?1.     

Experimental investigation of the metamorphism of siliceous dolomites

For the reaction: 1 diopside+3 dolomite ?2 forsterite+4 calcite+2 CO₂ (14) the following P _total?T-brackets have been determined experimentally in the presence of a gasphase consisting of 90 mole%CO₂ and 10 mole%H₂O∶1 kb, 544°±20° C; 3kb, 638°±15° C; 5kb, 708°±10° C; 10kb, 861°±10° C. The determination was carried out with well defined synthetic minerals in the starting mixture. The MgCO₃-contents of the magnesian calcites formed by the reaction in equilibrium with dolomite agree very well with the calcite-dolomite miscibility gap, which can be recalculated from the activities and the activity coefficients of MgCO₃ as given by Gordon and Greenwood (1970). The equilibrium constant K _14b was calculated with respect to the reference pressure P ₀=1 bar using the experimentally determined \(P_{total} TX_{CO_2 }\) brackets, the activities of MgCO₃ and CaCO₃ (Gordon and Greenwood 1970; Skippen 1974) and the fugacities of CO₂ Holloway (1977) considering the correction of Flowers (1979). Results are plotted as function of the absolute reciprocal temperature in Fig. 1. For the temperature range of 530° to 750° C the following linear expression can be given for the natural logarithm of K_14b: (g) $$[ln K_{14b} ]_T^P = - \frac{{18064.43}}{{T\left( {^\circ K} \right)}} + 38.58 + \frac{{0.308(P - 1 bar)}}{{T\left( {^\circ K} \right)}}$$ where P is the total pressure in bars and T the temperature in degrees Kelvin. Combining Equation (g) with the activities of MgCO₃ and CaCO₃ gives the equilibrium fugacity \(f_{CO_2 }\) : (i) $$[ln f_{CO_2 } ]_T^P = - \frac{{11635.44}}{{T\left( {^\circ K} \right)}} + 21.09 + \frac{{0.154(P - 1 bar)}}{{T\left( {^\circ K} \right)}}$$ Equation (i) and the fugacities of CO₂ permit to calculate the equilibrium data in terms of \(P_{CO_2 }\) and T (see Fig. 3) or P _total, T and \(X_{CO_2 }\) (see Fig. 5). Combining the \(P_{total} TX_{CO_2 }\) equilibrium data of the above reaction with those of the previously investigated reaction (Metz 1976): 1 tremolite+11 dolomite ?8 forsterite+13 calcite+9 CO₂+1 H₂O yields the stability conditions of the four-mineral assemblage: diopside+calcian dolomite+forsterite +magnesian calcite and the stability conditions of the five-mineral assemblage: tremolite+calcian dolomite+forsterite +magnesian calcite+diopside both shown in Fig. 6. Since these assemblages are by no means rare in metamorphic siliceous dolomites (Trommsdorff 1972; Suzuki 1977; Puhan 1979) the data of Fig. 6 can be used to determine the pressure of metamorphism and to estimate the composition of the CO₂-H₂O fluid provided the temperature of the metamorphic event was determined using the calcite-dolomite geothermometer.     

Energy gap and density in SiO2 polymorphs

A relationship between the energy gap (E _G) and density (ρ) for pure SiO₂ polymorphs is derived from atomic weights and first ionization potentials of free silicon and oxygen atoms. Theoretical considerations are based on the Lorentz electron theory of solids. The eigenfrequency v₀ of elementary electron oscillators, in energy units h v₀, is identified with the energy gap of a solid. The numerical relation is expressed as \(E_G = \sqrt {139.24 - 13.8327\rho } \) is in eV. For low-quartz with a density of 2.65 g/cm³ and also for stishovite with a density of 4.28 g/cm³, the energy gap E _G=10.1 eV and 8.9 eV, respectively. From laboratory measurements for low-quartz E _G=10.2 eV. The energy gap-density relation suggests a critical density value of ρ_x ≈ 10.1 g/cm³ for an SiO₂ phase when the energy gap vanishes (E _G=0), which is consistent with estimated densities for a high pressure silica polymorph with the fluorite structure.     

Ca-Fe-Mg olivines: phase relations and a solution model

Reversed phase equilibrium experiments in the system (Ca, Mg, Fe)₂SiO₄ provide four tielines at P?1 bar and 1 kbar and 800° C–1,100° C. These tielines have been used to model the solution properties of the olivine quadrilateral following the methods described by Davidson et al. (1981) for quadrilateral clinopyroxenes. The discrepancy between the calculated phase relations and the experimentally determined tielines is within the uncertainty of the experiments. The solution properties of quadrilateral olivines can be described by a non-convergent site-disorder model that allows for complete partitioning of Ca on the M2 site, highly disordered Fe-Mg cation distributions and limited miscibility between high-Ca and low-Ca olivines. The ternary data presented in this paper together with binary solution models for the joins Fo-Mo and Fa-Kst have been used to evaluate two solution parameters: $$\begin{gathered} F^0 \equiv 2(\mu _{{\rm M}o}^0 - \mu _{{\rm K}st}^0 ) + \mu _{Fa}^0 - \mu _{Fo}^0 = 12.660 (1.6) kJ, \hfill \\ \Delta G_*^0 \equiv \mu _{{\rm M}gFe}^0 + \mu _{FeMg}^0 - \mu _{Fo}^0 - \mu _{Fa}^0 = 7.030 (3.9) kJ. \hfill \\ \end{gathered} $$ . Ternary phase quilibrium data for olivines tightly constrain the value of F⁰, but not that for ΔG _* ⁰ which describes nonideality in Fe-Mg mixing. From this analysis, we infer a function for the apparent standard state energy of Kst: $$\begin{gathered} \mu _{{\rm K}st}^0 = - 102.79 \pm 0.8 - (T - 298)(0.137026) \hfill \\ + (T - 298 - T1n(T/298))(0.155519) \hfill \\ + (T - 298)^2 (2.8242E - 05)/2 \hfill \\ + (T - 298)^2 (2.9665E + 03)/(2T(298)^2 ) kJ \hfill \\ \end{gathered} $$ where T is in Kelvins and the 298 K value is relative to oxides.     

Thermodynamic properties of almandine-grossular garnet solid solutions

The mixing properties of Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnet solid solutions have been studied in the temperature range 850–1100° C. The experimental method involves measuring the composition of garnet in equilibrium with an assemblage in which the activity of the Ca₃Al₂Si₃O₁₂ component is fixed. Experiments on the assemblage garnet solid solution, anorthite, Al₂SiO₅ polymorph and quartz at known pressure and temperature fix the activity of the Ca₃Al₂Si₃O₁₂ component through the equilibrium: 1 $$\begin{gathered} {\text{3CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} \rightleftarrows {\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} \hfill \\ {\text{Anorthite garnet}} \hfill \\ {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + SiO}}_{\text{2}} \hfill \\ {\text{ sillimanite/kyanite quartz}}{\text{.}} \hfill \\ \end{gathered}$$ This equilibrium, with either sillimanite or kyanite as the aluminosilicate mineral, was used to control \({\text{a}}_{{\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} }^{{\text{gt}}} \) . The compositions of the garnet solutions produced were determined by measurement of their unit cell edges. At 1 bar Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnets exhibit negative deviations from ideality at the Fe-rich end of the series and positive deviations at the calcium end. With increasing pressure the activity coefficients for the Ca₃Al₂Si₃O₁₂ component increase because the partial molar volume of this component is greater than the molar volume of pure grossular. Previous studies indicate that the activity coefficients for the Ca₃Al₂Si₃O₁₂ component also increase with increasing (Mg/Mg+Fe) ratio of the garnet. The region of negative deviation from ideality implies a tendency towards formation of a stable Fe-Ca garnet component. Evidence in support of this conclusion has been found in a natural Fe-rich garnet which was found to contain two different garnet phases of distinctly different compositions.     

Crystallographic preferred orientation in wüstite (FeO) through the cubic-to-rhombohedral phase transition

Magnesiowüstite, (Mg_0.08Fe_0.88)O, and wüstite, Fe_0.94O, were compressed to ~36?GPa at ambient temperature in the diamond anvil cell (DAC) at the Advanced Light Source. X-ray diffraction patterns were taken in situ in radial geometry in order to study the evolution of crystallographic preferred orientation through the cubic-to-rhombohedral phase transition. Under uniaxial stress in the DAC, {100}_c planes aligned perpendicular to the compression direction. The {100}_c in cubic became { $\left\{ {10\bar 14} \right\}$ }_r in rhombohedral and remained aligned perpendicular to the compression direction. However, the {101}_c and {111}_c planes in the cubic phase split into { ${10{\bar{1}}4}$ }_r and { ${11{\bar{2}}0}$ }_r, and (0001)_r and { ${10{\bar{1}}1}$ }_r, respectively, in the rhombohedral phase. The { ${11{\bar{2}}0}$ }_r planes preferentially aligned perpendicular to the compression direction while { ${10{\bar{1}}4}$ }_r oriented at a low angle to the compression direction. Similarly, { ${10{\bar{1}}1}$ }_r showed a slight preference to align more closely perpendicular to the compression direction than (0001)_r. This variant selection may occur because the 〈 ${10{\bar{1}}4}$ 〉_r and [0001]_r directions are the softer of the two sets of directions. The rhombohedral texture distortion may also be due to subsequent deformation. Indeed, polycrystal plasticity simulations indicate that for preferred { ${10{\bar{1}}4}$ }〈 ${1{\bar{2}}10}$ 〉_r and { ${11{\bar{2}}0}$ }〈 ${{\bar{1}}101}$ 〉_r slip and slightly less active { ${10{\bar{1}}1}$ }〈 ${{\bar{1}}2{\bar{1}}0}$ 〉_r slip, the observed texture pattern can be obtained.