Basal slip and texture development in calcite: new results from torsion experiments

The deformation behavior of calcite has been of longstanding interest. Through experiments on single crystals, deformation mechanisms were established such as mechanical twinning on in the positive sense and slip on and both in the negative sense. More recently it was observed that at higher temperatures slip in both senses becomes active and, based on slip line analysis, it was suggested that slip may occur. So far there had been no direct evidence for basal slip, which is the dominant system in dolomite. With new torsion experiments on calcite single crystals at 900 K and transmission electron microscopy, this study identifies slip unambiguously by direct imaging of dislocations and diffraction contrast analysis. Including this slip system in polycrystal plasticity simulations, enigmatic texture patterns observed in compression and torsion of calcite rocks at high temperature can now be explained, resolving a long-standing puzzle.     

Structure and chemistry of (111) twin boundaries in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel crystals from Mogok

The atomic scale structure and chemistry of (111) twins in MgAl₂O₄ spinel crystals from the Pinpyit locality near Mogok (Myanmar, formerly Burma) were analysed using complementary methods of transmission electron microscopy (TEM). To obtain a three-dimensional information on the atomic structure, the twin boundaries were investigated in crystallographic projections and Using conventional electron diffraction and high-resolution TEM (HRTEM) analysis we have shown that (111) twins in spinel can be crystallographically described by 180° rotation of the oxygen sublattice normal to the twin composition plane. This operation generates a local hcp stacking in otherwise ccp lattice and maintains a regular sequence of kagome and mixed layers. In addition to rotation, no other translations are present in (111) twins in these spinel crystals. Chemical analysis of the twin boundary was performed by energy-dispersive X-ray spectroscopy (EDS) using a variable beam diameter (VBD) technique, which is perfectly suited for analysing chemical composition of twin boundaries on a sub-nm scale. The VBD/EDS measurements indicated that (111) twin boundary in spinel is Mg-deficient. Quantitative analyses of HRTEM (phase contrast) and HAADF-STEM (Z-contrast) images of (111) twin boundary have confirmed that Mg²⁺ ions are replaced with Be²⁺ ions in boundary tetrahedral sites. The Be-rich twin boundary structure is closely related to BeAl₂O₄ (chrysoberyl) and BeMg₃Al₈O₁₆ (taaffeite) group of intermediate polysomatic minerals. Based on these results, we conclude that the formation of (111) twins in spinel is a preparatory stage of polytype/polysome formation (taaffeite) and is a result of thermodynamically favourable formation of hcp stacking due to Be incorporation on the {111} planes of the spinel structure in the nucleation stage of crystal growth. The twin structure grows as long as the surrounding geochemical conditions allow its formation. The incorporation of Be induces a 2D-anisotropy and exaggerated growth of the crystal along the (111) twin boundary.     

High temperature calorimetry of sulfide systems

The standard enthalpies of formation of FeS (troilite), FeS₂ (pyrite), Co_0.9342S, Co₃S₄ (linnaeite), Co₉S₈ (cobalt pentlandite), CoS₂ (cattierite), CuS (covellite), and Cu₂S (chalcocite) have been determined by high temperature direct reaction calorimetry at temperatures between 700 K and 1021 K. The following results are reported: $$\Delta {\rm H}_{f,FeS}^{tr} = - 102.59 \pm 0.20kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,FeS}^{py} = - 171.64 \pm 0.93kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_{0.934} S} = - 99.42 \pm 1.52kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_9 S_8 }^{ptl} = - 885.66 \pm 16.83kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_3 S_4 }^{In} = - 347.47 \pm 7.27kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,CoS_2 }^{ct} = - 150.94 \pm 4.85kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Cu_2 S}^{cc} = - 80.21 \pm 1.51kJ mol^{ - 1} ,$$ and $$\Delta {\rm H}_{f,CuS}^{cv} = - 53.14 \pm 2.28kJ mol^{ - 1} ,$$ The enthalpy of formation of CuFeS₂ (chalcopyrite) from (CuS+FeS) and from (Cu+FeS₂) was determined by solution calorimetry in a liquid Ni_0.60S_0.40 melt at 1100 K. The results of these measurements were combined with the standard enthalpies of formation of CuS, FeS, and FeS₂, to calculate the standard enthalpy of formation of CuFeS₂. We found \(\Delta {\rm H}_{f,CuFeS_2 }^{ccp} = - 194.93 \pm 4.84kJ mol^{ - 1}\) . Our results are compared with earlier data given in the literature; generally the agreement is good and our values agree with previous estimates within the uncertainties present in both.     

Use of Pitzer Equations to Examine the Formation of Mercury(II) Hydroxide and Chloride Complexes in NaClO<Subscript>4</Subscript> Media

The thermodynamic stability constants for the hydrolysis and formation of mercury (Hg²⁺) chloride complexes

Twinning induced by the rhombohedral to orthorhombic phase transition in lanthanum gallate (LaGaO3)

Phase-transformation-induced twins in pressureless-sintered lanthanum gallate (LaGaO₃) ceramics have been analysed using the transmission electron microscopy (TEM). Twins are induced by solid state phase transformation upon cooling from the rhombohedral to orthorhombic (o, Pnma) symmetry at ∼145°C. Three types of transformation twins {101} _o , {121} _o , and {123} _o were found in grains containing multiple domains that represent orientation variants. Three orthorhombic orientation variants were distinguished from the transformation domains converged into a triple junction. These twins are the reflection type as confirmed by tilting experiment in the microscope. Although not related by group–subgroup relation, the transformation twins generated by phase transition from rhombohedral to orthorhombic are consistent with those derived from taking cubic aristotype of the lowest common supergroup symmetry as an intermediate metastable structure. The r→ o phase transition of first order in nature may have occurred by a diffusionless, martensitic-type or discontinuous nucleation and growth mechanism.     

Foundations of cable car towers upon alpine glaciers

This paper presents a design approach for strip footings upon glacier ice. Safety against ultimate limit state is proved by the geotechnical slip-line field solution by Prandtl. Glacier ice at 0°C can be modelled as purely cohesive material. Statistical evaluation of uniaxial compression tests with high strain rate revealed a mean value of the cohesion of 600 kPa and a characteristic value c _k = 355 kPa (5% fractile). With a coefficient of variation V _c = 0.3, the partial safety factor turns out to be γ _c = 1.9. An approximate solution for estimating the creep settlement rate is presented to check the serviceability limit state: with the width b of the strip foundation, p the foundation pressure and for ice at 0°C. Experiences on Stubai glacier with grate shaped footings showed that creep settlements occurring per year due to maximum foundation pressures 250 kPa did not influence the operation and the maintenance of the cable cars.     

Thermal expansion and spontaneous strain associated with the normal-incommensurate phase transition in melilites

The linear thermal expansions of åkermanite (Ca₂MgSi₂O₇) and hardystonite (Ca₂ZnSi₂O₇) have been measured across the normal-incommensurate phase transition for both materials. Least-squares fitting of the high temperature (normal phase) data yields expressions linear in T for the coefficients of instantaneous linear thermal expansion, $$\alpha _1 = \frac{1}{l}\frac{{dl}}{{dT}}$$ for åkermanite: $$\begin{gathered} \alpha _{[100]} = 6.901(2) \times 10^{ - 6} + 1.834(2) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 2.856(1) \times 10^{ - 6} + 11.280(1) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ for hardystonite: $$\begin{gathered} \alpha _{[100]} = 15.562(5) \times 10^{ - 6} - 1.478(3) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 11.115(5) \times 10^{ - 6} + 11.326(3) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ Although there is considerable strain for temperatures within 10° C of the phase transition, suggestive of a high-order phase transition, there appears to be a finite ΔV of transition, and the phase transition is classed as “weakly first order”.     

H<Subscript>2</Subscript>O diffusion in peralkaline to peraluminous rhyolitic melts

The diffusion of water in a peralkaline and a peraluminous rhyolitic melt was investigated at temperatures of 714–1,493 K and pressures of 100 and 500 MPa. At temperatures below 923 K dehydration experiments were performed on glasses containing about 2 wt% H₂O _t in cold seal pressure vessels. At high temperatures diffusion couples of water-poor (<0.5 wt% H₂O _t ) and water-rich (~2 wt% H₂O _t ) melts were run in an internally heated gas pressure vessel. Argon was the pressure medium in both cases. Concentration profiles of hydrous species (OH groups and H₂O molecules) were measured along the diffusion direction using near-infrared (NIR) microspectroscopy. The bulk water diffusivity () was derived from profiles of total water () using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between and Both methods consistently indicate that is proportional to in this range of water contents for both bulk compositions, in agreement with previous work on metaluminous rhyolite. The water diffusivity in the peraluminous melts agrees very well with data for metaluminous rhyolites implying that an excess of Al₂O₃ with respect to alkalis does not affect water diffusion. On the other hand, water diffusion is faster by roughly a factor of two in the peralkaline melt compared to the metaluminous melt. The following expression for the water diffusivity in the peralkaline rhyolite as a function of temperature and pressure was obtained by least-squares fitting:

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

Dissolution Kinetics of Dolomite in Water at Elevated Temperatures

Kinetic experiments of dolomite dissolution in water over a temperature range from 25 to 250°C were performed using a flow through packed bed reactor. Authors chose three different size fractions of dolomite samples: 18–35 mesh, 35–60 mesh, and 60–80 mesh. The dissolution rates of the three particle size samples of dolomite were measured. The dissolution rate values are changed with the variation of grain size of the sample. For the sample through 20–40 mesh, both the release rate of Ca and the release rate of Mg increase with increasing temperature until 200°C, then decrease with continued increasing temperature. Its maximum dissolution rate occurs at 200°C. The maximum dissolution rates for the sample through 40–60 mesh and 60–80 mesh happen at 100°C. Experimental results indicate that the dissolution of dolomite is incongruent in most cases. Dissolution of fresh dolomite was non-stoichiometric, the Ca/Mg ratio released to solution was greater than in the bulk solid, and the ratio increases with rising temperatures from 25 to 250°C. Observations on dolomite dissolution in water are presented as three parallel reactions, and each reaction occurs in consecutive steps as

A Coupled Riverine-Marine Fractionation Model for Dissolved Rare Earths and Yttrium

Fractionation of yttrium (Y) and the rare earth elements (REEs) begins in riverine systems and continues in estuaries and the ocean. Models of yttrium and rare earth (YREE) distributions in seawater must therefore consider the fractionation of these elements in both marine and riverine systems. In this work we develop a coupled riverine/marine fractionation model for dissolved rare earths and yttrium, and apply this model to calculations of marine YREE fractionation for a simple two-box (riverine/marine) geochemical system. Shale-normalized YREE concentrations in seawater can be expressed in terms of fractionation factors ( _ij ) appropriate to riverine environments ( ) and seawater ( ):

The role of iron in tetrahedrite and tennantite determined by Rietveld refinement of neutron powder diffraction data

Rietveld refinement of neutron powder diffraction data on four samples of synthetic, iron-bearing tetrahedrite (Cu_12?xFe_xSb₄S₁₃) with x = 0.28, 0.69, 0.91, 2.19 and four samples of synthetic tennantite (Cu_12?xFe_xAs₄S₁₃) with x = 0.33, 0.38, 0.86, 1.5 indicate unambiguously that iron is incorporated into tetrahedral M1 (12d) sites and not into triangular M2 (12e) sites in the cubic crystal structure (space group I $ \ifmmode\expandafter\bar\else\expandafter\=\fi{4} Rietveld refinement of neutron powder diffraction data on four samples of synthetic, iron-bearing tetrahedrite (Cu_12−xFe_xSb₄S₁₃) with x = 0.28, 0.69, 0.91, 2.19 and four samples of synthetic tennantite (Cu_12−xFe_xAs₄S₁₃) with x = 0.33, 0.38, 0.86, 1.5 indicate unambiguously that iron is incorporated into tetrahedral M1 (12d) sites and not into triangular M2 (12e) sites in the cubic crystal structure (space group I 3 m). The refinement results also confirm that M2 is a split (24g), flat-pyramidal site situated statistically on both sides of the S1−S1–S2 triangle. In tetrahedrite, this split is about 0.6 ?, in tennantite about 0.7 ?. Trends in bond lengths and magnitude of the M2 split were evaluated by means of linear regression with Fe concentration as the independent variable.     

Opening and resetting temperatures in heating geochronological systems

We present a theoretical model for diffusive daughter isotope loss in radiochronological systems with increasing temperature. It complements previous thermochronological models, which focused on cooling, and allows for testing opening and resetting of radiochronometers during heating. The opening and resetting temperatures are, respectively,

Influence of hydrogen on Fe–Mg interdiffusion in (Mg,Fe)O and implications for Earth’s lower mantle

Interdiffusion of Fe and Mg in (Mg,Fe)O has been investigated experimentally under hydrous conditions. Single crystals of MgO in contact with (Mg_0.73Fe_0.27)O were annealed hydrothermally at 300 MPa between 1,000 and 1,250°C and using a Ni–NiO buffer. After electron microprobe analyses, the dependence of the interdiffusivity on Fe concentration was determined using a Boltzmann–Matano analysis. For a water fugacity of ∼300 MPa, the Fe–Mg interdiffusion coefficient in Fe _x Mg_1−x O with 0.01 ≤ x ≤ 0.25 can be described by with and C = −80 ± 10 kJ mol⁻¹. For x = 0.1 and at 1,000°C, Fe–Mg interdiffusion is a factor of ∼4 faster under hydrous than under anhydrous conditions. This enhanced rate of interdiffusion is attributed to an increased concentration of metal vacancies resulting from the incorporation of hydrogen. Such water-induced enhancement of kinetics may have important implications for the rheological properties of the lower mantle.

Interaction of Freshly Precipitated Silica Gel with Aqueous Silicic Acid Solutions under Ambient and Near Neutral pH-conditions: A Detailed Analysis of Linear Rate Law

Interaction of freshly precipitated silica gel with aqueous solutions was studied at laboratory batch experiments under ambient and near neutral pH-conditions. The overall process showed excellent reversibility: gel growth could be considered as an opposite process to dissolution and a linear rate law could be applied to experimental data. Depending on the used rate law form, the resulting rate constants were sensitive to errors in parameters/variables such as gel surface area, equilibrium constants, Si-fluxes, and reaction quotients. The application of an Integrated Exponential Model appeared to be the best approach for dissolution data evaluation. It yielded the rate constants k _dissol ∼ (4.50 ± 0.68) × 10⁻¹² and k _growth ∼ (2.58 ± 0.39) × 10⁻⁹ mol m⁻² s⁻¹ for zero ionic strength. In contrast, a Differential Model gave best results for growth data modeling. It yielded the rate constants k _dissol ∼ (1.14 ± 0.44) × 10⁻¹¹ and k _growth ∼ (6.08 ± 2.37) × 10⁻⁹ mol m⁻² s⁻¹ for higher ionic strength (I ∼ 0.04 to 0.11 mol L⁻¹). The found silica gel solubility at zero ionic strength was somewhat lower than the generally accepted value. Based on the and standard Gibbs free energy of silica gel formation was calculated as and −850,318 ± 20 J mol⁻¹, respectively. Activation energies for silica gel dissolution and growth were determined as and respectively. An universal value for growth of any silica polymorph, is not consistent with the value for silica gel growth, which questions the hypothesis about one unique activated complex controlling the silica polymorph growth.     

New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers

The models recognize that ZrSiO₄, ZrTiO₄, and TiSiO₄, but not ZrO₂ or TiO₂, are independently variable phase components in zircon. Accordingly, the equilibrium controlling the Zr content of rutile coexisting with zircon is ZrSiO₄ = ZrO₂ (in rutile) + SiO₂. The equilibrium controlling the Ti content of zircon is either ZrSiO₄ + TiO₂ = ZrTiO₄ + SiO₂ or TiO₂ + SiO₂ = TiSiO₄, depending whether Ti substitutes for Si or Zr. The Zr content of rutile thus depends on the activity of SiO₂ as well as T, and the Ti content of zircon depends on and as well as T. New and published experimental data confirm the predicted increase in the Zr content of rutile with decreasing and unequivocally demonstrate that the Ti content of zircon increases with decreasing . The substitution of Ti in zircon therefore is primarily for Si. Assuming a constant effect of P, unit and that and are proportional to ppm Zr in rutile and ppm Ti in zircon, [log(ppm Zr-in-rutile) + log] = A₁ + B₁/T(K) and [log(ppm Ti-in-zircon) + log − log] = A₂ + B₂/T, where the A and B are constants. The constants were derived from published and new data from experiments with buffered by either quartz or zircon + zirconia, from experiments with defined by the Zr content of rutile, and from well-characterized natural samples. Results are A₁ = 7.420 ± 0.105; B₁ = −4,530 ± 111; A₂ = 5.711 ± 0.072; B₂ = −4,800 ± 86 with activity referenced to α-quartz and rutile at P and T of interest. The zircon thermometer may now be applied to rocks without quartz and/or rutile, and the rutile thermometer applied to rocks without quartz, provided that and are estimated. Maximum uncertainties introduced to zircon and rutile thermometry by unconstrained and can be quantitatively assessed and are ≈60 to 70°C at 750°C. A preliminary assessment of the dependence of the two thermometers on P predicts that an uncertainty of ±1 GPa introduces an additional uncertainty at 750°C of ≈50°C for the Ti-in-zircon thermometer and of ≈70 to 80°C for the Zr-in-rutile thermometer.     

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

In situ study of the $$ R\overline{3} c \to R\overline{3} m $$ orientational disorder in calcite

The temperature dependences of the crystal structure and intensities of the (113) and (211) reflections in calcite, CaCO₃, were studied using Rietveld structure refinements based on synchrotron powder X-ray diffraction data. Calcite transforms from to at about T _c = 1,240 K. A CO₃ group occupies, statistically, two positions with equal frequency in the disordered phase, but with unequal frequency in the partially ordered phase. One position for the CO₃ group is rotated by 180° with respect to the other. The unequal occupancy of the two orientations in the partially ordered phase is obtained directly from the occupancy factor, x, for the O1 site and gives rise to the order parameter, S = 2x − 1. The a cell parameter shows a negative thermal expansion at low T, followed by a plateau region at higher T, then a steeper contraction towards T _c, where the CO₃ groups disorder in a rapid process. Using a modified Bragg–Williams model, fits were obtained for the order parameter S, and for the intensities of the (113) and (211) reflections. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Fluid–rock interaction at a carbonatite-gneiss contact, Alnö, Sweden

We evaluate balanced metasomatic reactions and model coupled reactive and isotopic transport at a carbonatite-gneiss contact at Alnö, Sweden. We interpret structurally channelled fluid flow along the carbonatite-gneiss contact at ～640°C. This caused (1) metasomatism of the gneiss, by the reaction: ${\hbox{biotite} + \hbox{quartz} + \hbox{oligoclase} + \hbox{K}_{2} \hbox{O} +\,\hbox{Na}_{2}\hbox{O} \pm \hbox{CaO} \pm \hbox{MgO} \pm \hbox{FeO} = \hbox{albite} + \hbox{K-feldspar} + \hbox{arfvedsonite} + \hbox{aegirene-}\hbox{augite} + \hbox{H}_{2} \hbox{O} + \hbox{SiO}_{2}}We evaluate balanced metasomatic reactions and model coupled reactive and isotopic transport at a carbonatite-gneiss contact at Aln?, Sweden. We interpret structurally channelled fluid flow along the carbonatite-gneiss contact at ∼640°C. This caused (1) metasomatism of the gneiss, by the reaction: , (2) metasomatism of carbonatite by the reaction: calcite + SiO₂ = wollastonite + CO₂, and (3) isotopic homogenization of the metasomatised region. We suggest that reactive weakening caused the metasomatised region to widen and that the metasomatic reactions are chemically (and possibly mechanically) coupled. Spatial separation of reaction and isotope fronts in the carbonatite conforms to a chromatographic model which assumes local calcite–fluid equilibrium, yields a timescale of 10²–10⁴ years for fluid–rock interaction and confirms that chemical transport towards the carbonatite interior was mainly by diffusion. We conclude that most silicate phases present in the studied carbonatite were acquired by corrosion and assimilation of ijolite, as a reactive by-product of this process and by metasomatism. The carbonatite was thus a relatively pure calcite–H₂O−CO₂–salt melt or fluid.