Gehlenite stability in the system CaO-Al2O3-SiO2-H2O-CO2

P, T, \(X_{{\text{CO}}_{\text{2}} }\) relations of gehlenite, anorthite, grossularite, wollastonite, corundum and calcite have been determined experimentally at P _f =1 and 4 kb. Using synthetic starting minerals the following reactions have been demonstrated reversibly

1. 2 anorthite+3 calcite=gehlenite+grossularite+3 CO₂.
2. anorthite+corundum+3 calcite=2 gehlenite+3 CO₂.
3. 3anorthite+3 calcite=2 grossularite+corundum+3CO₂.
4. grossularite+2 corundum+3 calcite=3 gehlenite+3 CO₂.
5. anorthite+2 calcite=gehlenite+wollastonite+2CO₂.
6. anorthite+wollastonite+calcite=grossularite+CO₂.
7. grossularite+calcite=gehlenite+2 wollastonite+CO₂.

In the T, \(X_{{\text{CO}}_{\text{2}} }\) diagram at P _f =1 kb two isobaric invariant points have been located at 770±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.27 and at 840±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.55. Formation of gehlenite from low temperature assemblages according to (4) and (2) takes place at 1 kb and 715–855° C, \(X_{{\text{CO}}_{\text{2}} }\) =0.1–1.0. In agreement with experimental results the formation of gehlenite in natural metamorphic rocks is restricted to shallow, high temperature contact aureoles.     

Experimental investigation of the mechanism of the reaction: 1 tremolite+11 dolomite ⇌ 8 forsterite+13 calcite+9 CO2+1H2O

Stoichiometric mixtures of tremolite and dolomite were heated to 50° C above equilibrium temperatures to form forsterite and calcite. The pressure of the CO₂-H₂O fluid was 5 Kb and \(X_{{\text{CO}}_{\text{2}} }\) varied from 0.1 to 0.6. The extent of the conversion was determined by the amount of CO₂ produced. The resulting mixtures of unreacted tremolite and dolomite and of newly-formed forsterite and calcite were examined with a scanning electron microscope. All tremolite and dolomite grains showed obvious signs of dissolution. At fluid compositions with \(X_{{\text{CO}}_{\text{2}} }\) less than about 0.4, the forsterite and calcite crystals are randomly distributed throughout the charges, indicating that surfaces of the reactants are not a controlling factor with respect to the sites of nucleation of the products. A change is observed when \(X_{{\text{CO}}_{\text{2}} }\) is greater than about 0.4; the forsterite and calcite crystals now nucleate and grow at the surface of the dolomite grains, thus indicating a change in mechanism at medium CO₂ concentrations. As the reaction progresses, the dolomite grains become more and more surrounded by forsterite and calcite, finally forming armoured relics of dolomite. Under experimental conditions this characteristic texture can only be formed if the CO₂-concentration is greater than about 40 mole %. These findings make it possible to estimate the CO₂-concentration from the texture of the dolomite+tremolite+forsterite+calcite assemblage. The results suggest a dissolution-precipitation mechanism for the reaction investigated. In a simplified form it consists of the following 4 steps:

1. Dissolution of the reactants tremolite and dolomite.
2. Diffusion of the dissolved constituents in the fluid.
3. Heterogeneous nucleation of the product minerals.
4. Growth of forsterite and calcite from the fluid.

Two possible explanations are discussed for the development of the amoured texture at \(X_{{\text{CO}}_{\text{2}} }\) above 0.4. The first is based upon the assumption that dolomite has a lower rate of dissolution than tremolite at high \(X_{{\text{CO}}_{\text{2}} }\) values resulting in preferential calcite and forsterite nucleation and growth on the dolomite surface. An alternative explanation is the formation of a raised CO₂ concentration around the dolomite grains at high \(X_{{\text{CO}}_{\text{2}} }\) values, leading to product precipitation on the dolomite crystals.     

The ionic conductivity in high-pressure polymorphs of calcite

Ionic conductivity of polycrystalline calcite containing varying amounts of PO ₄ ^3? ions was measured in the pressure range of 1–6 GPa and at room temperature. Electrical conductivity increased with pressure corresponding to the phase transition of calcite I to calcite II. The conductivity in calcite III decreased exponentially with pressure. Calculated activation volumes of the conductivity varied with PO ₄ ^3? content in the range of 0.94–5.34 cm³/mol. This variation corresponded to the lattice parameter change of calcite I due to PO ₄ ^3? incorporation and indicated the contribution of CO ₃ ^2? -vacancies associated with PO ₄ ^3? ions to the conductivity.     

Experimentelle Untersuchung der Metamorphose von kieselig dolomitischen Sedimenten

During an experimental investigation of the metamorphism of siliceous dolomites the equilibrium data of the heterogeneous bivariant reaction 1 $$3{\text{ dolomite + 4 quartz + 1 H}}_{\text{2}} O \rightleftharpoons + 3 calcite + 3 CO_2 $$ were determined for the total fluid pressures of 1,000, 3,000 and 5,000 bars. The equilibrium conditions were found by experiments in which dolomite, quartz and water react to form talc, calcite and CO₂, as well as by experiments with reversible reaction direction. Results are shown on the temperature- \(X_{CO_2 } \) -diagram of Fig. 3. The temperature of formation of talc and calcite depends to a considerable extent on the composition of the CO₂-H₂O-gas phase; this can be read straight off the isobaric (P _f =const.) equilibrium curves in Fig. 3. In addition a strong dependence of the equilibrium temperature on the total pressure P _f was established (see Fig. 5). At a total gas pressure of 1,000 bars dolomite and quartz can react, according to the composition of the CO₂-H₂O-gas phase, to talc and calcite over the whole of the temperature range between about 350° and 490° C. This indicates that at low pressures the formation of talc and calcite takes place in the field of the albite-epidote-hornfels facies. At a pressure of 3,000 bars dolomite and quartz are stable up to about 550° C if the fluid phase is rich in carbon dioxide and correspondingly poor in water. Thus, this paragenesis can occur up to the stability field of staurolite [see annotation (5)] if the partial pressure of CO₂ is large. At the higher total gas pressure of 5,000 bars dolomite and quartz react even at medium CO₂-concentrations only at about 580° C to give talc and calcite. Therefore it is expected that in regional metamorphism at about 5,000 bars pressure or more the paragenesis dolomite plus quartz exists up to and within the stability field of staurolite and reacts only here to form talc and calcite after reaction (1) or tremolite and calcite after the following reaction (2)¹: $$5 dolomite + 8 quartz + 1 H_2 O \rightleftharpoons 1 tremolite + 3 calcite + 7 CO_2 $$ . The exact physico-chemical conditions under which dolomite, quartz and water react on the one hand to form talc, calcite and CO₂, and on the other hand to form tremolite, calcite and carbon dioxide, will be discussed later when our experimental investigations on the formation of tremolite are completed. First results were already published in a short note by Metz, Puhan and Winkler (1968).     

Origin and diagenetic evolution of Ca–Mg–Fe carbonates in some coalfields of Japan

Carbonate concretions, lenses and bands in the Pleistocene, Palaeogene and Upper Triassic coalfields of Japan consist of various carbonate minerals with varied chemical compositions. Authigenic carbonates in freshwater sediments are siderite > calcite > ankerite > dolomite >> ferroan magnesite; in brackish water to marine sediments in the coal measures, calcite > dolomite > ankerite > siderite >> ferroan magnesite; and in the overlying marine deposits, calcite > dolomite >> siderite. Most carbonates were formed progressively during burial within a range of depths between the sediment-water interface and approximately 3 km. The mineral species and the chemical composition of the carbonates are controlled primarily by the initial sedimentary facies of the host sediments and secondarily by the diagenetic evolution of pore water during burial. Based on the regular sequence and burial depth of precipitation of authigenic carbonates in a specific sedimentary facies, three diagenetic stages of carbonates are proposed. Carbonates formed during Stage I (< 500 m) strongly reflect the initial sedimentary facies, e.g. low Ca-Mg siderite in freshwater sediments which are initially rich in iron derived from lateritic soil on the nearby landmass, and Mg calcite and dolomite in brackish-marine sediments whose pore waters abound in Ca²⁺ and Mg²⁺ originating in seawater and calcareous shells. Carbonates formed during Stage II (500–2000 m) include high Ca-Mg siderite, ankerite, Fe dolomite and Fe–Mg calcite in freshwater sediments. The assemblage of Stage II carbonates in brackish-marine sediments in the coal measures is similar to that in freshwater sediments. This suggests similar diagenetic environments owing to an effective migration and mixing of pore water due to the compaction of host sediments. Carbonates formed during Stage III (> 2000 m) are Fe calcite and extremely high Ca-Mg siderite; the latter is exclusively in marine mudstones. The supply of Ca is partly from the alteration of silicates in the sediments at elevated burial temperatures. After uplift, calcite with low Mg content precipitates from percolating groundwater and fills extensional cracks.     

Experimental delineation of the “e” ⇌ I\bar 1 and “e” ⇌ C\bar 1 transformations in intermediate plagioclase feldspars

Approximately 125 hydrothermal annealing experiments have been carried out in an attempt to bracket the stability fields of different ordered structures within the plagioclase feldspar solid solution. Natural crystals were used for the experiments and were subjected to temperatures of ～650°C to ～1,000°C for times of up to 370 days at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =600 bars, or \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =1,200 bars. The structural states of both parent and product materials were characterised by electron diffraction, with special attention being paid to the nature of type e and type b reflections (at h+k=(2n+1), l=(2n+1) positions). Structural changes of the type C \(\bar 1\) → I \(\bar 1\) , C \(\bar 1\) → “e” structure, I \(\bar 1\) → “e” and “e” structure → I \(\bar 1\) have been followed. There are marked differences between the ordering behaviour of crystals with compositions on either side of the C \(\bar 1\) ? I \(\bar 1\) transition line. In the composition range ～ An₅₀ to ～ An₇₀ the e structure appears to have a true field of stability relative to I \(\bar 1\) ordering, and a transformation of the type I \(\bar 1\) ? e has been reversed. It is suggested that the e structure is the more stable ordered state at temperatures of ～ 800°C and below. For compositions more albite-rich than ～ An₅₀ the upper temperature limit for long range e ordering is lower than ～ 750°C, and there is no evidence for any I \(\bar 1\) ordering. The evidence for a true stability field for “e” plagioclase, which is also consistent with calorimetric data, necessitates reanalysis both of the ordering behaviour of plagioclase crystals in nature and of the equilibrium phase diagram for the albite-anorthite system. Igneous crystals with compositions of ～ An₆₅, for example, probably follow a sequence of structural states C \(\bar 1\) → I \(\bar 1\) → e during cooling. The peristerite, Bøggild and Huttenlocher miscibility gaps are clearly associated with breaks in the albite, e and I \(\bar 1\) ordering behaviour but their exact topologies will depend on the thermodynamic character of the order/disorder transformations.     

Deformation twinning in calcite,dolomite, and other rhombohedral carbonates

Transmission electron microscope (tem) observations of single and multiple twins in calcite and dolomite are presented, and the results are analysed by means of selected area diffraction and trace analysis. Simple twinning in rhodochrosite and kutnahorite is also analysed. It is shown that the ordered carbonates, such as dolomite, have a common twinning plane {01 \(\bar 1\) 2} and this appears to be their only mode of deformation twinning. The carbonates with higher symmetry, such as calcite, have {01 \(\bar 1\) 8} as the primary twinning plane but calcite itself has other twinning mechanisms, of which the most important is illustrated. Crossing and stopping twins are also discussed. It is shown that twinning in calcite, which occurs predominantly at low temperatures, is characterized by the generation of large numbers of glide dislocations.     

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.     

Die Reaktion 2 Zoisit+1 CO2 ⇌ 3 Anorthit+1 Calcit+1 H2O

The equilibrium conditions of the following reaction 2 zoisite +1 CO₂?3 anorthite+1 calcite+1 H₂O 2 Ca₂Al₃[O/OH/SiO₄/Si₂O₇]+1 CO₂?3 CaAl₂Si₂O₈+1 CaCO₃+1 H₂O have been determined experimentally at total pressures of P _j= 2000 bars, P _f =5000 bars, and P _f =7000 bars. Owing to the vertical position of the equilibrium curves in isobaric T- \(X_{{\text{CO}}_{\text{2}} }\) diagrams, the composition of the binary H₂O-CO₂ fluid phase coexisting with zoisite is independent of temperature in the temperature interval investigated. According to our experiments, orthorhombic zoisite is only stable in equilibrium with a fluid phase at a concentration of CO₂ which is less than, respectively, ca. 2 Mol% CO₂ at P _f =2000 bars, ea. 6 Mol% at P _f =5000 bars, and ca. 10 Mol% at P _f =7000 bars. Thus, the fluid phase coexisting with zoisite is rich in H₂O. While this is independent of temperature the experimental data demonstrate that the influence of pressure cannot be neglected: With increasing pressure the concentration of CO₂ of the fluid phase coexisting with zoisite can rise a little. The position of the reaction studied, which is independent of temperature and exhibits small values of \(X_{{\text{CO}}_{\text{2}} }\) ,leads to two important petrogenetic conclusions:

1. The occurrence of zoisite is an indicator for a CO₂-poor and H₂O-rich fluid composition during metamorphism of marly calcsilicates.
2. If the concentration of CO₂ of the fluid phase coexisting with zoisite exceeds the equilibrium value of \(X_{{\text{CO}}_{\text{2}} }\) calcite+anorthite+H₂O is formed from zoisite+CO₂. Thus, a considerable increase in the anorthite-content of plagioelase is possible.

     

Transmission electron microscopy on meteoritic troilite

Phase transitions and associated domains of meteoritic troilite (FeS) have been studied by means of transmission electron microscopy (TEM). Three polymorphs have been found, two of which can be described by superstructures of the NiAs-type structure (A, C subcell). The P \(\overline 6\) 2c (√3A, 2C) polymorph, stable at room temperature, displays antiphase domains with the displacement vector 1/3< \(\overline {\text{1}}\) 10>. In situ heating experiments showed that the P \(\overline 6\) 2c polymorph changes at temperatures of 115°–150° C into an orthorhombic pseudohexagonal transitional phase with the probable space group Pmcn (A,√3A, C). It contains antiphase domains with the displacement vector 1/2 [110] and twins with a threefold twin-axis parallel c. When heated above 210° C the transitional phase transforms into the high-temperature modification with NiAs structure (P6 ₃/mmc). All observed phase transitions are reversible. The occurrence of antiphase and twin domains, respectively, agrees with the symmetry reductions involved in the subsolidus phase transitions. This is demonstrated by group-subgroup relationships among the space groups P6 ₃/mmc, Pmcn, and P \(\overline 6\) 2c.     

Suitability of Existing Numerical Model Codes and Thermodynamic Databases for the Prognosis of Calcite Dissolution Processes in Near-Surface Sediments Due to a CO2 Leakage Investigated by Column Experiments

Among the risks of CO₂ storage is the potential of CO₂ leakage into overlaying formations and near-surface potable aquifers. Through a leakage, the CO₂ can intrude into protected groundwater resources, which can lead to groundwater acidification followed by potential mobilisation of heavy metals and other trace metals through mineral dissolution or ion exchange processes. The prediction of pH buffer reactions in the formations overlaying a CO₂ storage site is essential for assessing the impact of CO₂ leakages in terms of trace metal mobilisation. For buffering the pH-value, calcite dissolution is one of the most important mechanisms. Although calcite dissolution has been studied for decades, experiments conducted under elevated CO₂ partial pressures are rare. Here, the first study for column experiments is presented applying CO₂ partial pressures from 6 to 43 bars and realising a near-natural flow regime. Geochemical calculations of calcite dissolution kinetics were conducted using PHREEQC together with different thermodynamic databases. Applying calcite surface areas, which were previously acquired by N₂-BET or calculated based on grain diameters, respectively, to the rate laws according to Plummer et al. (Am J Sci 278:179–216, doi:10.2475/ajs.278.2.179, 1978) or Palandri and Kharaka (US Geol Surv Open file Rep 2004–1068:71, 2004) in the numerical simulations led to an overestimation of the calcite dissolution rate by up to three orders of magnitude compared to the results of the column experiments. Only reduction of the calcite surface area in the simulations as a fitting procedure allowed reproducing the experimental results. A reason may be that the diffusion boundary layer (DBL), which depends on the groundwater flow velocity and develops at the calcite grain surface separating it from the bulk of the solution, has to be regarded: The DBL leads to a decrease in the calcite dissolution rate under natural laminar flow conditions compared to turbulent mixing in traditional batch experiments. However, varying the rate constants by three orders of magnitudes in a field scale PHREEQC model simulating a CO₂ leakage produced minor variations in the pH buffering through calcite dissolution. This justifies the use of equilibrium models when calculating the calcite dissolution in CO₂ leakage scenarios for porous aquifers and slow or moderate groundwater flow velocities. However, the selection of the thermodynamic database has an impact on the dissolved calcium concentration, leading to an uncertainty in the simulation results. The resulting uncertainty, which applies also to the calculated propagation of an aquifer zone depleted in calcite through dissolution, seems negligible for shallow aquifers of approximately 60 m depth, but amounts to 35 % of the calcium concentration for aquifers at a depth of approximately 400 m.     

Carbon and oxygen diffusion in calcite: Effects of Mn content and P H2O

The diffusion rates of carbon and oxygen in two calcite crystals of different Mn contents have been studied between 500° and 800° C in a CO₂-H₂O atmosphere (P _{CO 2}=1?5 bars, P _H2_O=0.02?24 bars) labeled with ¹³C and ¹⁸O. Isotope concentration gradients within annealed specimens were measured using a secondary ion microprobe by depth profiling parallel and perpendicular to the c axis. Despite the anisotropic structure of calcite, the diffusion of carbon and oxygen are both very nearly isotropic. Least-squares fitting of the carbon data to an Arrhenius relation gives an activation energy of 87±2 kcal/mole, with D ₀ terms dependent only slightly upon direction: 1 $$D_{\text{0}} {\text{(}}\parallel c{\text{) = }}\left( {9\frac{{ + 12}}{{ - 5}}} \right){\text{x10}}^{\text{2}} cm^2 /s$$ , 2 $$D_{\text{0}} {\text{(}} \bot c{\text{) = }}\left( {5\frac{{ + 6}}{{ - 3}}} \right){\text{x10}}^{\text{2}} cm^2 /s$$ . These results are in close agreement with previous determinations. Results for oxygen diffusion, however, give D values much larger than those previously reported for dry conditions; at 650° to 800° C the D values are two orders of magnitude larger. The diffusion of oxygen, unlike carbon, is strongly dependent on water pressure, as well as Mn content, and does not fit an Arrhenius relation over the entire temperature range. On the basis of these observations and considerations of the defect chemistry of calcite, it is proposed that carbon migrates as a Frenkel pair. The diffusion of oxygen, however, appears to be more complicated and may depend upon several simultaneous mechanisms.     

The high-temperature heat capacity of natural calcite (CaCO3)

The heat capacity (C _P) of a natural sample of calcite (CaCO₃) has been measured from 350 to 775 K by differential scanning calorimetry (DSC). Heat capacities determined for a powdered sample and a single-crystal disc are in close agreement and have a total uncertainty of ±1 percent. The following equation for the heat capacity of calcite from 298 to 775 K was fit by least squares to the experimental data and constrained to join smoothly with the low-temperature heat capacity data of Staveley and Linford (1969) (C _P in J mol^?1 K^?1, T in K): $$\begin{gathered} C_p = - 184.79 + 0.32322T - 3,688,200T^{ - 2} \hfill \\ {\text{ }} - (1.2974{\text{ }} \times {\text{ 10}}^{ - {\text{4}}} )T^2 + 3,883.5T^{ - 1/2} \hfill \\ \end{gathered} $$ Combining this equation with the S ₂₉₈ ⁰ value from Staveley and Linford (1969), entropies for calcite are calculated and presented to 775 K. A simple method of extrapolating the heat capacity function of calcite above 775 K is presented. This method provides accurate entropies of calcite for high-temperature thermodynamic calculations, as evidenced by calculation of the equilibrium: CaCO_{3 (s)}=CaO_(s)+CO_{2 (s)}.     

Modulated structures in calcian dolomite: A study by electron microscopy

Calcian dolomite from the Devonian Lost Burro formation has been investigated with electron microscopy techniques. Electron diffraction shows evidence for “c” and “d” type reflections which may occur independently and are indicative of ordered superstructures. High resolution electron microscopy combined with selected area optical diffraction is the basis for models to explain the superstructures in calcian dolomite. It is proposed that “c” reflections are due to ordered substitution of Mg by Ca in basal cation layers. “d” reflections result when the rhombohedral stacking of basal layers is interrupted by intercalation of additional Ca layers. During electron irradiation at 1 MeV the Mg-Ca distribution becomes disordered and the crystal structure attains calcite symmetry. The arrangement of CO₃ groups remains ordered.     

The thermodynamics of the I \overline{1} –P \overline{1} phase transition in Ca-rich plagioclase from an assessment of the spontaneous strain

In-situ powder diffraction measurements between 90 and 935?K on four anorthite-rich plagioclase samples (An₁₀₀, An₉₆Ab₄, An₈₉Ab₁₁ and An₇₈Ab₂₂) were used to determine the detailed evolution of these samples through the $I \overline{1} $ – $P \overline{1} $ phase transition. The c-type reflections indicative of $P \overline{1} $ symmetry were detected only in An₁₀₀, An₉₆Ab₄, whereas deviations in the evolution of the unit-cell parameters with temperature were observed in all samples, most prominently in the β unit-cell angle. The c-type reflections disappear at ~510 and ~425?K in An₁₀₀ and An₉₆Ab₄ respectively, and their intensity decreases according to a tricritical trend $ I^{2} \propto \left( {T - T_{\text{c}} } \right) $ . The cell parameter changes were used to determine the spontaneous strains arising from the transition which were modelled with Landau theory, allowing for low-temperature quantum saturation, in order to determine the thermodynamic behaviour. In An₁₀₀ tricritical behaviour was observed [T _c?=?512.7(4)?K; θ_s?=?394(4)] in good agreement with previous studies, and the c-type superlattice reflections indicative of $P \overline{1} $ symmetry persist up to the T _c determined from the spontaneous strain, and then disappear. The evolution of the spontaneous strain in An₉₆Ab₄ is tricritical at low temperatures [T _c?=?459(1) K, θ_s?=?396(5)] up to the temperature of disappearance of c-type reflections, but becomes second order beyond ~440?K. In An₈₉Ab₁₁ the strain displays second-order behaviour throughout [T _c?=?500(1) and θ_s?=?212(5)], and the c-type reflections are not detected in the powder diffraction patterns at any temperature. The apparent discrepancy between the absence of c-type reflections in temperature ranges where the cell parameters display significant spontaneous strain is resolved through consideration of the sizes of the anti-phase domains within the crystals. It is deduced that the tricritical phase transition occurs in well-ordered crystals with large domains in which the behavior of individual domains is dominant (i.e. in pure anorthite) or where the $P \overline{1} $ distortions within the domains are large enough to dominate the structural coherency strains between the domains. When both the magnitude of the $P \overline{1} $ pattern of displacements of the tetrahedral framework become smaller and the influence of the structural coherency between anti-phase domains becomes significant, the thermodynamic behavior becomes 2nd-order in character, the c-type reflections disappear, and the orientation of the spontaneous strain changes.     

A note on the significance of the assemblage calcite-quartz-plagioclase-paragonite-graphite

The biotite zone assemblage: calcite-quartz-plagioclase (An₂₅)-phengite-paragonite-chlorite-graphite, is developed at the contact between a carbonate and a pelite from British Columbia. Thermochemical data for the equilibrium paragonite+calcite+2 quartz=albite+ anorthite+CO₂+H₂O yields: $$\log f{\text{H}}_{\text{2}} {\text{O}} + \log f{\text{CO}}_{\text{2}} = 5.76 + 0.117 \times 10^{ - 3} (P - 1)$$ for a temperature of 700°K and a plagioclase composition of An₂₅. By combining this equation with equations describing equilibria between graphite and gas species in the system C-H-O, the following partial pressures: \(P{\text{H}}_2 {\text{O}} = 2572{\text{b, }}P{\text{CO}}_2 = 3162{\text{b, }}P{\text{H}}_2 = 2.5{\text{b, }}P{\text{CH}}_4 = 52.5{\text{b, }}P{\text{CO}} = 11.0{\text{b}}\) are obtained for \(f{\text{O}}_2 = 10^{ - 26}\) . If total pressure equals fluid pressure, then the total pressure during metamorphism was approximately 6 kb. The total fluid pressure calculated is extremely sensitive to the value of \(f{\text{O}}_2\) chosen.     

Contribution of {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} and 13CTDIC Evaluation to the Identification of CO2 Sources in Volcanic Groundwater Systems: Influence of Hydrometeorological Conditions and Lava Flow Morphologies—Application to the Argnat Basin (Chaîne des Puys, Massif Central, France)

Mineralization of groundwater in volcanic aquifers is partly acquired through silicates weathering. This alteration depends on the dissolution of atmospheric, biogenic, or mantellic gaseous CO₂ whose contributions may depend on substratum geology, surface features, and lava flow hydrological functionings. Investigations of $ {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} $ and δ¹³C_TDIC (total dissolved inorganic carbon) on various spatiotemporal scales in the unsaturated and saturated zones of volcanic flows of the Argnat basin (French Massif Central) have been carried out to identify the carbon sources in the system. Mantellic sources are related to faults promoting CO₂ uplift from the mantle to the saturated zone. The contribution of this source is counterbalanced by infiltration of water through the unsaturated zone, accompanied by dissolution of soil CO₂ or even atmospheric CO₂ during cold periods. Monitoring and modeling of δ¹³C_TDIC in the unsaturated zone shows that both $ {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} $ and δ¹³C_TDIC are controlled by air temperature which influences soil respiration and soil-atmosphere CO₂ exchanges. The internal geometry of volcanic lava flows controls water patterns from the unsaturated zone to saturated zone and thus may explain δ¹³C heterogeneity in the saturated zone at the basin scale.     

Experimentelle Untersuchung der Reaktion Dolomit + Kalifeldspat + H2O ⇌ Phlogopit + Calcit + CO2

The equilibrium curve for the reaction 3 dolomite + 1 K-feldspar + 1 H₂O=1 phlogopite + 3 calcite + 3 CO₂ was determined experimentally at a total gas pressure of 2000 bars using two different methods.

1. In the first case water alone was added to the reactants. The CO₂ component of the gas phase was producted solely by the reaction under favourable P-T conditions. This manner of carrying out the reaction is called the “water method”. With this method sufficient time must be allowed for the gas phase to attain a constant composition (see Fig. 1). Reverse reactions were carried out using reaction products of the forward reaction.
2. In the second case silver oxalate + water were added to the reactants. Breakdown of the silver oxalate leads to formation of a CO₂-H₂O gasphase of definite composition. At constant temperature and gas pressure the \(X_{{\text{CO}}_{\text{2}} } \) determines whether the reaction products will be phlogopite + calcite or dolomite + K-feldspar. In this case it is not necessary to wait for equilibrium to be attained. This method is abbreviated the “oxalate method”. Reactants for reverse reactions are not identical with the products of the forward reaction.

At high temperatures the results of the two different methods agree well (see Tables 1 and 2). Equilibrium was attained in one case at 490° C and \(X_{{\text{CO}}_{\text{2}} } \) of approximately 0.77, and in the other case at 520° C and \(X_{{\text{CO}}_{\text{2}} } \) of 0.90. At lower temperatures there are considerable differences in the results. With the water method an \(X_{{\text{CO}}_{\text{2}} } \) of about 0.25 was reached at 450° C. With the oxalate method dolomite K-feldspar and water still react with each other at even higher \(X_{{\text{CO}}_{\text{2}} } \) values. Phlogopite, calcite and CO₂ are formed together with metastable talc. There are no criteria to indicate which of the methods is the correct one at lower temperatures and in Fig. 2, therefore, both equilibrium curves are plotted.     

Bimetasomatic zoning in the CaO-MgO-SiO2-H2O-CO2 system: Experiments with the use of natural rock samples

Experiments reproducing the development of bimetasomatic zoning in the CaO-MgO-SiO₂-H₂O-CO₂ system were conducted at elevated P-T parameters with the use of samples of naturally occurring quartzdolomite and calcite-serpentinite rocks. In order to maintain mass transfer exclusively via the diffusion-controlled mechanism, we used the method of the ensured compaction of the cylindrical sample surface with a thin-walled gold tube. In the course of the experiments, a single diopside zone ～2.5 × 10^?5 m thick was obtained at the quartz-dolomite interface at T = 600°C, $P_{H_2 O + CO_2 } $ = 200 MPa, and $X_{CO_2 } $ = 0.5 for 25–40 days and a succession of metasomatic zones at T = 750°C, $P_{H_2 O + CO_2 } $ = 300 MPa, and $X_{CO_2 } $ = 0.4 for 48 days. The metasomatic zones were as follows (listed in order from quartz to dolomite): wollastonite ‖ diopside ‖ tremolite ‖ calcite + forsterite; with the average width of the diopside zone equal to ～1.3 × 10^?5 m and the analogous part of the wollastonite zone equal to ～2.6 × 10^?5 m. Two zones (listed in order from calcite to serpentine) diopside and diopside-forsterite (the average widths of these zones were ～6 × 10^?4 and ～8 × 10^?4 m, respectively) were determined to develop at contact between serpentine and calcite during experiments that lasted 124 days at T = 500°C, $P_{H_2 O + CO_2 } $ = 200 MPa, and $X_{CO_2 } $ = 0.2–0.4. In the former and latter situations, the growth rate of the zoning ranged between 3.1 × 10^?12 and 1.2 × 10^?11 m/s and between 5.6 × 10^?11 and 7.5 × 10^?11 m/s, respectively. The higher growth rate in the latter case can be explained by the higher water mole fraction in the fluid, with this water released during serpentinite decomposition in the experiments. The development of the only diopside zone in the experiments modeling the interaction of quartz and dolomite at T = 600–650°C and $P_{H_2 O + CO_2 } $ = 200 MPa is in conflict with theoretical considerations underlain by the Korzhinskii-Fisher-Joesten model. The interaction of quartz and dolomite in the CaO-MgO-SiO₂-CO₂-H₂O system at the P-T- $X_{CO_2 } $ parameters specified above should be attended by the origin of a number of reaction zones consisting of various proportions of talc, forsterite, tremolite, diopside, and calcite. The saturation of the fluid with respect to these minerals was likely not reached, and this resulted in the degeneration of the respective stability fields in the succession of zones. Conceivably, this was related to the insufficient rates of quartz and dolomite dissolution and the relatively low diffusion rates of the dissolved species in the low-permeable medium. In the experiments with interacting calcite and serpentine, the zoning calcite ‖ diopside ‖ diopside + forsterite ‖ serpentine developed in its complete form, in agreement with the theory. Equilibrium was likely achieved in these experiments due to the higher diffusion coefficients.     

Ore-forming solution and genesis of a rock crystal deposit in south Hunan,indicated by fluid inclusion studies

Si, Al, Ca, Mg, Fe, Na, K, CO ₃ ^?2 , F, etc. are detected from the fluid inclusion leachates. Among these constituents, Si, Na, and CO ₃ ^?2 are predominant, amounting to more than 80 percent. This indicates that the ore-forming solution must be alkaline with Si, Na, and CO ₃ ^?2 as its dominant components. Homogenization temperatures for the solution range from 80 to 360°C. Although rock quartz can crystallize at the above temperature interval, perfect crystals of economic importance are largely formed below 260°C. The temperature of formation increases toward the granite intrusives at a rate of about one degree per meter. It is estimated from the lithostatic load that the salinity of rock quartz is 17–23 (NaCl wt%), while that of vein quartz is relatively high as compared with the former. There is a tendency for the salinity of the ore-forming solution to increase with depth.