Phase relations of carbonate-bearing eclogite assemblages from 2.5 to 5.5 GPa: implications for petrogenesis of carbonatites

We have experimentally investigated the phase and melting relations of garnet + clinopyroxene + carbonate assemblages at 2.5–5.5 GPa, to assess the feasibility of carbonated eclogite as a source for some crustally emplaced carbonatites. The solidus of our composition was at 1,125 °C at 2.5 GPa, 1,225 °C at 3.5 GPa and 1,310 °C at 5.0 GPa. Melts were sodic calcio-dolomitic carbonatites, and were markedly more calcic than the dolomitic melts produced by partial melting of carbonated peridotite. Na contents of the experimental carbonatites decreased with increasing pressure when compared at similar degrees of melting, and SiO₂ contents increased with degree of melting. Experiments on a second composition with enhanced Na₂O demonstrated its strong effect in lowering melting temperatures in carbonate eclogite. Natural carbonated eclogite bodies in the peridotitic upper mantle will have a range of solidus temperatures. In many cases, carbonate will be molten in the upper 250 km. Carbonate melt would segregate from its source eclogite at very low melt fractions and infiltrate surrounding peridotitic wall rock. This would result in metasomatic enrichment of the peridotitic wall rock, but its exact nature will depend on the relative P–T positions of the eclogite + CO₂ and peridotite + CO₂ solidii. As a result of these inevitable metasomatic interactions, it is considered unlikely that carbonatite melts derived from carbonated eclogite in the upper mantle could be emplaced into the crust unmodified. However, they may have a role in metasomatically enriching and carbonating parts of the upper mantle, producing sources suitable for subsequent production of silica undersaturated silicate liquids and carbonatites ultimately emplaced in the crust.Editorial responsibility: J. Hoefs     

Compressional behavior of omphacite to 47 GPa

Omphacite is an important mineral component of eclogite. Single-crystal synchrotron X-ray diffraction data on natural (Ca, Na) (Mg, Fe, Al)Si₂O₆ omphacite have been collected at the Advanced Photon Source beamlines 13-BM-C and 13-ID-D up to 47 GPa at ambient temperature. Unit cell parameter and crystal structure refinements were carried out to constrain the isothermal equation of state and compression mechanism. The third-order Birch–Murnaghan equation of state (BM3) fit of all data gives V ₀ = 423.9(3) Å³, K _T0 = 116(2) GPa and K _T0′ = 4.3(2). These elastic parameters are consistent with the general trend of the diopside–jadeite join. The eight-coordinated polyhedra (M2 and M21) are the most compressible and contribute to majority of the unit cell compression, while the SiO₄ tetrahedra (Si1 and Si2) behave as rigid structural units and are the most incompressible. Axial compressibilities are determined by fitting linearized BM3 equation of state to pressure dependences of unit cell parameters. Throughout the investigated pressure range, the b-axis is more compressible than the c-axis. The axial compressibility of the a-axis is the largest among the three axes at 0 GPa, yet it quickly drops to the smallest at pressures above 5 GPa, which is explained by the rotation of the stiffest major compression axis toward the a-axis with the increase in pressure.     

The strength of ruby from X-ray diffraction under non-hydrostatic compression to 68 GPa

Polycrystalline ruby (α-Al₂O₃:Cr³⁺), a widely used pressure calibrant in high-pressure experiments, was compressed to 68.1 GPa at room temperature under non-hydrostatic conditions in a diamond anvil cell. Angle-dispersive X-ray diffraction experiments in a radial geometry were conducted at beamline X17C of the National Synchrotron Light Source. The stress state of ruby at high pressure and room temperature was analyzed based on the measured lattice strain. The differential stress of ruby increases with pressure from ~3.4 % of the shear modulus at 18.5 GPa to ~6.5 % at 68.1 GPa. The polycrystalline ruby sample can support a maximum differential stress of ~16 GPa at 68.1 GPa under non-hydrostatic compression. The results of this study provide a better understanding of the mechanical properties of this important material for high-pressure science. From a synthesis of existing data for strong ceramic materials, we find that the high-pressure yield strength correlates well with the ambient pressure Vickers hardness.     

Textures in deforming forsterite aggregates up to 8 GPa and 1673 K

We report results from axisymmetric deformation experiments carried out on forsterite aggregates in the deformation-DIA apparatus, at upper mantle pressures and temperatures (3.1–8.1 GPa, 1373–1673 K). We quantified the resulting lattice preferred orientations (LPO) and compare experimental observations with results from micromechanical modeling (viscoplastic second-order self-consistent model—SO). Up to 6 GPa (~185-km depth in the Earth), we observe a marked LPO consistent with a dominant slip in the (010) plane with one observation of a dominant [100] direction, suggesting that [100](010) slip system was strongly activated. At higher pressures (deeper depth), the LPO becomes less marked and more complex with no evidence of a dominant slip system, which we attribute to the activation of several concurrent slip systems. These results are consistent with the pressure-induced transition in the dominant slip system previously reported for olivine and forsterite. They are also consistent with the decrease in the seismic anisotropy amplitude observed in the Earth’s mantle at depth greater than ~200 km.     

Phase relation of CaSO<Subscript>4</Subscript> at high pressure and temperature up to 90 GPa and 2300 K

Calcium sulfate (CaSO₄), one of the major sulfate minerals in the Earth’s crust, is expected to play a major role in sulfur recycling into the deep mantle. Here, we investigated the crystal structure and phase relation of CaSO₄ up to ~90 GPa and 2300 K through a series of high-pressure experiments combined with in situ X-ray diffraction. CaSO₄ forms three thermodynamically stable polymorphs: anhydrite (stable below 3 GPa), monazite-type phase (stable between 3 and ~13 GPa) and barite-type phase (stable up to at least 93 GPa). Anhydrite to monazite-type phase transition is induced by pressure even at room temperature, while monazite- to barite-type transition requires heating at least to 1500 K at ~20 GPa. The barite-type phase cannot always be quenched from high temperature and is distorted to metastable AgMnO₄-type structure or another modified barite structure depending on pressure. We obtained the pressure–volume data and density of anhydrite, monazite- and barite-type phases and found that their densities are lower than those calculated from the PREM model in the studied P–T conditions. This suggests that CaSO₄ is gravitationally unstable in the mantle and fluid/melt phase into which sulfur dissolves and/or sulfate–sulfide speciation may play a major role in the sulfur recycling into the deep Earth.     

Stability of phlogopite in ultrapotassic kimberlite-like systems at 5.5–7.5 GPa

Hydrous K-rich kimberlite-like systems are studied experimentally at 5.5–7.5 GPa and 1200–1450?°C in terms of phase relations and conditions for formation and stability of phlogopite. The starting samples are phlogopite–carbonatite–phlogopite sandwiches and harzburgite–carbonatite mixtures consisting of Ol?+?Grt?+?Cpx?+?L (±Opx), according to the previous experimental results obtained at the same P–T parameters but in water-free systems. Carbonatite is represented by a K- and Ca-rich composition that may form at the top of a slab. In the presence of carbonatitic melt, phlogopite can partly melt in a peritectic reaction at 5.5 GPa and 1200–1350?°C, as well as at 6.3–7.0 GPa and 1200?°C: 2Phl?+?CaCO₃ (L)?Cpx?+?Ol?+?Grt?+?K₂CO₃ (L)?+?2H₂O (L). Synthesis of phlogopite at 5.5 GPa and 1200–1350?°C, with an initial mixture of H₂O-bearing harzburgite and carbonatite, demonstrates experimentally that equilibrium in this reaction can be shifted from right to left. Therefore, phlogopite can equilibrate with ultrapotassic carbonate–silicate melts in a?≥?150?°C region between 1200 and 1350?°C at 5.5 GPa. On the other hand, it can exist but cannot nucleate spontaneously and crystallize in the presence of such melts in quite a large pressure range in experiments at 6.3–7.0 GPa and 1200?°C. Thus, phlogopite can result from metasomatism of peridotite at the base of continental lithospheric mantle (CLM) by ultrapotassic carbonatite agents at depths shallower than 180–195 km, which creates a mechanism of water retaining in CLM. Kimberlite formation can begin at 5.5 GPa and 1350?°C in a phlogopite-bearing peridotite source generating a hydrous carbonate–silicate melt with 10–15 wt% SiO₂, Ca# from 45 to 60, and high K enrichment. Upon further heating to 1450?°C due to the effect of a mantle plume at the CLM base, phlogopite disappears and a kimberlite-like melt forms with SiO₂ to 20 wt% and Ca#?=?35–40.     

The compressibility of Fe- and Al-bearing phase D to 30 GPa

High-pressure in situ X-ray diffraction experiment of Fe- and Al-bearing phase D (Mg_0.89Fe_0.14Al_0.25Si_1.56H_2.93O₆) has been carried out to 30.5 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields values of V ₀ = 86.10 ± 0.05 ų; K ₀ = 136.5 ± 3.3 GPa and K′ = 6.32 ± 0.30. If K′ is fixed at 4.0 K ₀ = 157.0 ± 0.7 GPa, which is 6% smaller than Fe–Al free phase D reported previously. Analysis of axial compressibilities reveals that the c-axis is almost twice as compressible (K _c  = 93.6 ± 1.1 GPa) as the a-axis (K _a  = 173.8 ± 2.2 GPa). Above 25 GPa the c/a ratio becomes pressure independent. No compressibility anomalies related to the structural transitions of H-atoms were observed in the pressure range to 30 GPa. The density reduction of hydrated subducting slab would be significant if the modal amount of phase D exceeds 10%.     

Structural evolution of hemimorphite at high pressure up to 4.2 GPa

The high-pressure structural evolution of hemimorphite, Zn₄Si₂O₇(OH)₂·H₂O, a = 8.3881(13), b = 10.7179(11), c = 5.1311(9) Å, V = 461.30(12) Å³, space group Imm2, Z = 2, was studied by single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions up to 4.2 GPa. In the pressure range of 0.0001–2.44 GPa, the unit-cell parameters change almost linearly. The phase transition (probably of the second order) with symmetry reduction from Imm2 (hemimorphite-I) to Pnn2 (hemimorphite-II) was found near 2.5 GPa. The structure compressibility increases somewhat above the phase transition. Namely, the initial unit-cell volume decreases by 3.6% at 2.44 GPa and by 7.2% at 4.20 GPa. The hemimorphite framework can be described as built up of secondary building units (SBU) Zn₄Si₂O₇(OH)₂. These blocks are combined to form the rods arranged along the c-axis; these rods are multiplied by basic and I-translations of orthorhombic unit cell. The symmetry reduction is caused by the rotation of the rods along their axis. In hemimorphite-I, the compression affects mainly the SBU dimensions, whereas a rectangular section of the channels having mm2 symmetry remains practically unchanged. An appreciable decrease in this section in hemimorphite-II is determined by its oblique distortion with the loss of m planes. It results from opposite rotation of adjacent SBU, which also leads into the loss of I-translation. In hemimorphite-I, the coordination of H₂O molecules is fourfold planar; the hydrogen-bonded hydroxyls and H₂O molecules form infinite ribbons along the c-axis. In hemimorphite-II, an additional short H₂O–O contact appears as a result of asymmetric deformation of the channels. The appearance of this new contact provides the possibility for re-orientation of hydrogen bonds. The planar coordination of H₂O molecules changes to tetrahedral and the ribbons are transformed to islands (OH)₂–H₂O.     

Stability of Kerogen Classification with Regard to Image Segmentation

This paper investigates the stability of an automatic system for classifying kerogen material from images of sieved rock samples. The system comprises four stages: image acquisition, background removal, segmentation, and classification of the segmented kerogen pieces as either inertinite or vitrinite. Depending upon a segmentation parameter d, called “overlap”, touching pieces of kerogen may be split differently. The aim of this study is to establish how robust the classification result is to variations of the segmentation parameter. There are two issues that pose difficulties in carrying out an experiment. First, even a trained professional may be uncertain when distinguishing between isolated pieces of inertinite and vitrinite, extracted from transmitted-light microscope images. Second, because manual labelling of large amount of data for training the system is an arduous task, we acquired the true labels (ground truth) only for the pieces obtained at overlap d=0.5. To construct ground truth for various values of d we propose here label-inheritance trees. With thus estimated ground truth, an experiment was carried out to evaluate the robustness of the system to changes in the segmentation through varying the overlap value d. The average system accuracy across values of d spanning the range from 0 to 1 was 86.5%, which is only slightly lower than the accuracy of the system at the design value of d=0.5 (89.07%).     

Melting of portlandite up to 6 GPa

 Using the high-pressure differential thermal analysis (HP-DTA) system in a cubic multianvil high-pressure apparatus, we measured the melting points of portlandite, Ca(OH)₂, up to 6 GPa and 1000 °C. We detected endothermic behavior at the temperature and pressure conditions of 800 °C and 2.5 GPa, 769 °C and 3.5 GPa, 752 °C and 4.0 GPa, 686 °C and 5.0 GPa, and 596 °C and 6.0 GPa, respectively, due to melting of portlandite. By in situ X-ray studies under pressure, the melting of portlandite was observed at 730 °C and 4.32 GPa and at 640 °C and 5.81 GPa, respectively. Results of both HP-DTA and X-ray studies were consistent within experimental error. The melting is congruent and has a negative Clapeyron slope, indicating that liquid Ca(OH)₂ has higher densities than crystalline portlandite in this pressure range. Received: 19 June 1999 / Revised, accepted: 11 September 1999     

Equation of state of adamite up to 11 GPa: a synchrotron X-ray diffraction study

The compression behavior of natural adamite [Zn₂AsO₄OH] has been investigated up to 11.07 GPa at room temperature utilizing in situ angle-dispersive X-ray diffraction and a diamond anvil cell. No phase transition has been observed within the pressure range investigated. A third-order Birch–Murnaghan equation of state fitted to all of the data points yielded V ₀ = 430.1(4) Å³, K ₀ = 80(3) GPa, K′ ₀ = 1.9(5). The K ₀ was obtained as 69(1) GPa when K′ ₀ was fixed at 4. Analysis of axial compressible moduli shows the intense compression anisotropy of adamite: K _a0 = 37(3) GPa, K _b0 = 153(6) GPa, K _c0 = 168(8) GPa; hence, a axis is the most compressible and the compressibility of b and c axis is comparable. Furthermore, the comparisons among the compressional properties of adamite, libethenite (Cu₂PO₄OH, also belongs to olivenite group), and andalusite (Al₂SiO₄O has the similar structure with adamite) at high pressure were made.     

Electrical conductivity of NaCl-bearing aqueous fluids to 600 °C and 1 GPa

The electrical conductivity of aqueous fluids containing 0.01, 0.1, and 1 M NaCl was measured in an externally heated diamond cell to 600 °C and 1 GPa. These measurements therefore more than double the pressure range of previous data and extend it to higher NaCl concentrations relevant for crustal and mantle fluids. Electrical conductivity was generally found to increase with pressure and fluid salinity. The conductivity increase observed upon variation of NaCl concentration from 0.1 to 1 M was smaller than from 0.01 to 0.1 M, which reflects the reduced degree of dissociation at high NaCl concentration. Measured conductivities can be reproduced (R ² = 0.96) by a numerical model with log \(\sigma\) = ?1.7060– 93.78/T + 0.8075 log c + 3.0781 log \(\rho\) + log \(\varLambda\) ₀(T, \(\rho\)), where \(\sigma\) is the conductivity in S m^?1, T is temperature in K, c is NaCl concentration in wt%, \(\rho\) is the density of pure water (in g/cm³) at given pressure and temperature, and \(\varLambda\) ₀ (T, \(\rho\)) is the molar conductivity of NaCl in water at infinite dilution (in S cm² mol^?1), \(\varLambda\) ₀ = 1573–1212 \(\rho\) + 537 062/T–208 122 721/T ². This model allows accurate predictions of the conductivity of saline fluids throughout most of the crust and upper mantle; it should not be used at temperatures below 100 °C. In general, the data show that already a very small fraction of NaCl-bearing aqueous fluid in the deep crust is sufficient to enhance bulk conductivities to values that would be expected for a high degree of partial melting. Accordingly, aqueous fluids may be distinguished from hydrous melts by comparing magnetotelluric and seismic data. H₂O–NaCl fluids may enhance electrical conductivities in the deep crust with little disturbance of v _p or v _p/v _s ratios. However, at the high temperatures in the mantle wedge above subduction zones, the conductivity of hydrous basaltic melts and saline aqueous fluids is rather similar, so that distinguishing these two phases from conductivity data alone is difficult. Observed conductivities in forearc regions, where temperatures are too low to allow melting, may be accounted for by not more than 1 wt% of an aqueous fluid with 5 wt% NaCl, if this fluid forms a continuous film or fills interconnected tubes.     

Crystal structures of Zn2SiO4 III and IV synthesized at 6.5–8 GPa and 1,273 K

We report the crystal structures determined under ambient condition for two Zn₂SiO₄ polymorphs synthesized at 6.5 GPa and 1,273 K (phase III) and 8 GPa and 1,273 K (phase IV) and also compare their ²⁹Si MAS NMR spectroscopic characteristics with those of other Zn₂SiO₄ polymorphs (phases I, II and V). Electron microprobe analysis revealed that both of phases III and IV are stoichiometric like the lower-pressure polymorphs (phases I and II), contrary to previous report. The crystal structures were solved using an ab initio structure determination technique from synchrotron powder X-ray diffraction data utilizing local structural information from ²⁹Si MAS NMR as constraints and were further refined with the Rietveld technique. Phase III is orthorhombic (Pnma) with a = 10.2897(5), b = 6.6711(3), c = 5.0691(2) Å. It is isostructural with the high-temperature (Zn_1.1Li_0.6Si_0.3)SiO₄ phase and may be regarded as a ‘tetrahedral olivine’ type that resembles the ‘octahedral olivine’ structure in the (approximately hexagonally close packed) oxygen arrangement and tetrahedral Si positions, but has Zn in tetrahedral, rather than octahedral coordination. Phase IV is orthorhombic (Pbca) with a = 10.9179(4), b = 9.6728(4), c = 6.1184(2) Å. It also consists of tetrahedrally coordinated Zn and Si and features unique edge-shared Zn₂O₆ dimers. The volumes per formula under ambient condition for phases III and IV are both somewhat larger than that of the lower-pressure polymorph, phase II, suggesting that the two phases may have undergone structural changes during temperature quench and/or pressure release.     

Patterns of Phase Formation in Kimberlite under Reducing Conditions at 4 GPa and 1500°C

Doklady Earth Sciences - Experimental data on the interaction of an iron melt with natural kimberlite at a temperature of 1500 ± 25°C and a pressure of 4.0 ± 0.2 GPa corresponding to...     

Thermo-elastic behaviour of Be2BO3OH (hambergite) up to 7 GPa and 1,100 K

The thermo-elastic behaviour of Be₂BO₃(OH)_0.96F_0.04 (i.e. natural hambergite, Z = 8, a = 9.7564(1), b = 12.1980(2), c = 4.4300(1) Å, V = 527.21(1) ų, space group Pbca) has been investigated up to 7 GPa (at 298 K) and up to 1,100 K (at 0.0001 GPa) by means of in situ single-crystal X-ray diffraction and synchrotron powder diffraction, respectively. No phase transition or anomalous elastic behaviour has been observed within the pressure range investigated. P?V data fitted to a third-order Birch–Murnaghan equation of state give: V ₀ = 528.89(4) ų, K _T0 = 67.0(4) GPa and K′ = 5.4(1). The evolution of the lattice parameters with pressure is significantly anisotropic, being: K _T0(a):K _T0(b):K _T0(c) = 1:1.13:3.67. The high-temperature experiment shows evidence of structure breakdown at T > 973 K, with a significant increase in the full-width-at-half-maximum of all the Bragg peaks and an anomalous increase in the background of the diffraction pattern. The diffraction pattern was indexable up to 1,098 K. No new crystalline phase was observed up to 1,270 K. The diffraction data collected at room-T after the high-temperature experiment showed that the crystallinity was irreversibly compromised. The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α ₀ + α ₁ T ^?1/2. The refined parameters for Be₂BO₃(OH)_0.96F_0.04 are: α ₀ = 7.1(1) × 10^?5 K^?1 and α ₁ = ?8.9(2) × 10^?4 K ^?1/2 for the unit-cell volume, α ₀(a) = 1.52(9) × 10^?5 K^?1 and α ₁(a) = ?1.4(2) × 10^?4 K ^?1/2 for the a-axis, α ₀(b) = 4.4(1) × 10^?5 K^?1 and α ₁(b) = ?5.9(3) × 10^?4 K ^?1/2 for the b-axis, α ₀(c) = 1.07(8) × 10^?5 K^?1 and α ₁(c) = ?1.5(2) × 10^?4 K ^?1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α ₀(a):α ₀(b):α ₀(c) = 1.42:4.11:1. The main deformation mechanisms in response to the applied temperature, based on Rietveld structure refinement, are discussed.