Liquid immiscibility in deep-seated magmas and the generation of carbonatite melts

This paper reviews the results of investigations of melt inclusions in minerals of carbonatites and spatially associated silicate rocks genetically related to various deep-seated undersaturated silicate magmas of alkaline ultrabasic, alkaline basic, lamproitic, and kimberlitic compositions. The analysis of this direct genetic information showed that all the deep magmas are inherently enriched in volatile components, the most abundant among which are carbon dioxide, alkalis, halides, sulfur, and phosphorus. The volatiles probably initially served as agents of mantle metasomatism and promoted melting in deep magma sources. The derived magmas became enriched in carbon dioxide, alkalis, and other volatile components owing to the crystallization and fractionation of early high-magnesium minerals and gradually acquired the characteristics of carbonated silicate liquids. When critical compositional parameters were reached, the accumulated volatiles catalyzed immiscibility, the magmas became heterogeneous, and two-phase carbonate-silicate liquid immiscibility occurred at temperatures of ≥1280–1250°C. The immiscibility was accompanied by the partitioning of elements: the major portion of fluid components partitioned together with Ca into the carbonate-salt fraction (parental carbonatite melt), and the silicate melt was correspondingly depleted in these components and became more silicic. After spatial separation, the silicate and carbonate-silicate melts evolved independently during slow cooling. Differentiation and fractionation were characteristic of silicate melts. The carbonatite melts became again heterogeneous within the temperature range from 1200 to 800–600°C and separated into immiscible carbonate-salt fractions of various compositions: alkali-sulfate, alkali-phosphate, alkali-fluoride, alkali-chloride, and Fe-Mg-Ca carbonate. In large scale systems, polyphase silicate-carbonate-salt liquid immiscibility is usually manifested during the slow cooling and prolonged evolution of deeply derived melts in the Earth’s crust. It may lead to the formation of various types of intrusive carbonatites: widespread calcite-dolomite and rare alkali-sulfate, alkali-phosphate, and alkali-halide rocks. The initial alkaline carbonatite melts can retain their compositions enriched in P, S, Cl, and F only at rapid eruption followed by instantaneous quenching.     

Crystallization of authigenic carbonates in mud volcanoes at Lake Baikal

This paper presents data on authigenic siderite first found in surface sediments from mud volcanoes in the Central (K-2) and Southern (Malen’kii) basins of Lake Baikal. Ca is the predominant cation, which substitutes Fe in the crystalline lattice of siderite. The enrichment of the carbonates in the ¹³C isotope (from +3.3 to +6.8‰ for the Malen’kii volcano and from +17.7 to +21.9‰ for K-2) results from the crystallization of the carbonates during methane generation via the bacterial destruction of organic matter (acetate). The overall depletion of the carbonates in ¹⁸O is mainly inherited from the isotopic composition of Baikal water.     

Experimental study of fluorine and chlorine contents in mica (biotite) and their partitioning between mica,phonolite melt,and fluid

This paper reports the results of a study of the composition of mica (biotite) crystallizing in the system of phonolite melt-Cl- and F-bearing aqueous fluid at T ～ 850°C, P = 200 MPa, and \(f_{O_2 } \) = Ni-NiO, as well as data on F and Cl partitioning between coexisting phases. It was established that Cl content in mica is significantly lower than in phonolite melt and, especially, in fluid. Fluorine shows a different behavior in this system: its content in mica is always higher than in phonolite melt but lower than in fluid. The mica-melt partition coefficients of Cl and F also behave differently. The Cl partition coefficient gradually increases from 0.17 to 0.33 with increasing Cl content in the system, whereas the partition coefficient of F sharply decreases from 3.0 to 1.0 with increasing total F content. The apparent partition coefficients of F between biotite and groundmass (melt) in various magmatic rocks are usually significantly higher than the experimental values. It was supposed that the higher ^Bt/glassD_F values in natural samples could be related to the influence of later oxidation reactions, reequilibration of biotite at continuously decreasing \(f_{H_2 O} \)/f _HF ratio, and an increase in this coefficients with decreasing total F content in the system.     

Inclusions of silicate and sulfate melts in chrome diposide from the Inagli deposit,Yakutia, Russia

Melt inclusions were studied in chrome diopside from the Inagli deposit of gemstones in the Inagli massif of alkaline ultrabasic rocks of potassic affinity in the northwestern Aldan shield, Yakutia, Russia. The chrome diopside is highly transparent and has an intense green color. Its Cr₂O₃ content varies from 0.13 to 0.75 wt %. Primary and primary-secondary polyphase inclusions in chrome diopside are dominated by crystal phases (80–90 vol %) and contain aqueous solution and a gas phase. Using electron microprobe analysis and Raman spectroscopy, the following crystalline phases were identified. Silicate minerals are represented by potassium feldspar, pectolite [NaCa₂Si₃O₈(OH)], and phlogopite. The most abundant minerals in the majority of inclusions are sulfates: glaserite (aphthitalite) [K₃Na(SO₄)₂], glauberite [Na₂Ca(SO₄)₂], aluminum sulfate, anhydrite (CaSO₄), gypsum (CaSO₄ × 2H₂O), barite (BaSO₄), bloedite [Na₂Mg(SO₄)₂ × 4H₂O], thenardite (NaSO₄), polyhalite [K₂Ca₂Mg(SO₄)₄ × 2H₂O], arcanite (K₂SO₄), and celestite (SrSO₄). In addition, apatite was detected in some inclusions. Chlorides are probably present among small crystalline phases, because some analyses of aggregates of silicate and sulfate minerals showed up to 0.19–10.3 wt % Cl. Hydrogen was identified in the gas phase of polyphase inclusions by Raman spectroscopy. The composition of melt from which the chrome diopside crystallized was calculated on the basis of the investigation of silicate melt inclusions. This melt contains 53.5 wt % SiO₂, considerable amounts of CaO (16.3 wt %), K₂O (7.9 wt %), Na₂O (3.5 wt %), and SO₃ (1.4 wt %) and moderate amounts of Al₂O₃ (7.5 wt %), MgO (5.8 wt %), FeO (1.1 wt %), and H₂O (0.75 wt %). The content of Cr₂O₃ in the melt was 0.13 wt %. Many inclusions were homogenized at 770–850°C, when all of the crystals and the gas phase were dissolved. The material of inclusions heated up to the homogenization temperature became heterogeneous even during very fast quenching (two seconds) producing numerous small crystals. This fact implies that most of the inclusions contained a salt (rather than silicate) melt of sulfate-dominated composition. Such inclusions were formed from salt globules (with a density of about 2.5 g/cm³) occurring as an emulsion in the denser (2.6 g/cm³) silicate melt from which the chrome diopside crystallized.     

Late Palaeozoic hydrocarbon migration through the Clair field, West of Shetland, UK Atlantic margin

Geochemical analysis of bitumen- and hydrocarbon-bearing fluid inclusions from the Devonian-Carboniferous Clair field indicates that the reservoirs contain a mixture of oils from different marine and lacustrine sources. Reconstruction of the Clair field oil-charge history using fluid inclusion petrography show that oil-charging occurred at times of K-feldspar, quartz and calcite cementation. Temperature–composition–time data yielded from the integration of fluid inclusion microthermometry with high-resolution Ar–Ar dating, date hydrocarbon-bearing K-feldspar overgrowths at 247 ± 3.3 Ma. These data show that in order for oil to be trapped within primary fluid inclusions in K-feldspar overgrowths, hydrocarbon migration throughout the UK Atlantic margin must have been taking place during the Late Palaeozoic and as such, current industry oil-play models based solely on oil charging from Jurassic-Cretaceous marine sources are clearly incomplete and need revision. Apatite fission track analysis and vitrinite reflectance data were used to reconstruct thermal burial histories and assess potential oil generation from Middle Devonian lacustrine source rocks. Thermal history data from wells along The Rona Ridge adjacent to the Clair field show that the Palaeozoic section was heated to greater than 100 °C at some time between 270 and 230 Ma, confirming that Devonian source rocks were mature and expelling oil during the Late Palaeozoic at the time that authigenic K-feldspar overgrowths were growing in the Clair field.     

Water release patterns of heated speleothem calcite and hydrogen isotope composition of fluid inclusions

Speleothem fluid inclusions are a potential paleo-precipitation proxy to reconstruct past rainwater isotopic composition (δ¹⁸O, δD). To get a better insight in the extraction of inclusion water from heated speleothem calcite, we monitored the water released from crushed and uncrushed speleothem calcite, heated to 900 °C at a rate of 300 °C/h, with a quadrupole mass spectrometer. Crushed calcite released water in three not well individualised peaks between 25 and 360 °C, 360 and 650 °C and between 650 and 800 °C while uncrushed calcite released water in two distinct temperature intervals: between 25 and 550 °C and between 550 and 900 °C.Water from two speleothems from the Han-sur-Lesse cave was recovered using three different techniques: i) the crushing and heating to 360 °C technique, ii) the decrepitation by heating to 550 °C and iii) the decomposition by heating to 900 °C technique. Measurements of the δD of water recovered by the decomposition of Han-sur-Lesse calcite heated to 900 °C did not show a 20 to 30‰ offset as found by previous authors. However a difference of 7‰ was observed between water released before and after decomposition of the calcite. Water recovery from the Han-sur-Lesse samples suggests that a simple heating technique (up to 550 °C) without crushing could both (a) recover water with δD representative of that of the drip water and (b) double the water yield as compared to the crushing and heating method.Our study warns for possible contamination of the recovered inclusion water with hydration water of lime, responsible for the recovery of water with very negative δD values.     

Li-Sr-Nd isotope signatures of the plume and cratonic lithospheric mantle beneath the margin of the rifted Tanzanian craton (Labait)

Lithium concentrations and isotopic compositions of olivine and ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd of coexisting clinopyroxene from peridotite xenoliths from the Quaternary Labait volcano, Tanzania, document the influence of rift-related metasomatism on the ancient cratonic mantle. Olivines show negative correlations between Fo content and both δ⁷Li and Li concentrations. Olivines in iron-rich peridotites (Fo_85–87) have high Li concentrations (3.2–4.8 ppm) and heavy δ⁷Li (+5.2 to +6.6). In contrast, olivines in ancient, refractory peridotites have lower Li concentrations (∼2 ppm) and relatively light δ⁷Li (+2.6 to +3.5). This reflects mixing between ancient, refractory cratonic lithosphere and asthenosphere-derived rift magmas. A uniquely fertile, deformed, high-temperature garnet lherzolite, interpreted to be from the base of the lithosphere, has a ⁸⁷Sr/⁸⁶Sr of 0.7029 and ¹⁴³Nd/¹⁴⁴Nd of 0.51286, similar to HIMU oceanic basalts. It provides the best estimate of the Sr–Nd isotope composition of the upwelling mantle (i.e., plume, sensu lato) underlying this portion of the East African Rift, and is slightly less radiogenic compared to previous estimates of the plume that were based on rift basalts. Although elevated δ⁷Li are not exclusive to HIMU source regions, the data collectively indicate that the plume beneath Labait has HIMU characteristics in Sr, Nd and Li isotope composition. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Chemical and boron-isotope variations in tourmalines from an S-type granite and its source rocks: the Erongo granite and tourmalinites in the Damara Belt,Namibia

Tourmaline is widespread in metapelites and pegmatites from the Neoproterozoic Damara Belt, which form the basement and potential source rocks of the Cretaceous Erongo granite. This study traces the B-isotope variations in tourmalines from the basement, from the Erongo granite and from its hydrothermal stage. Tourmalines from the basement are alkali-deficient schorl-dravites, with B-isotope ratios typical for continental crust (δ¹¹B average −8.4‰ ± 1.4, n = 11; one sample at −13‰, n = 2). Virtually all tourmaline in the Erongo granite occurs in distinctive tourmaline-quartz orbicules. This “main-stage” tourmaline is alkali-deficient schorl (20–30% X-site vacancy, Fe/(Fe + Mg) 0.8–1), with uniform B-isotope compositions (δ¹¹B −8.7‰ ± 1.5, n = 49) that are indistinguishable from the basement average, suggesting that boron was derived from anatexis of the local basement rocks with no significant shift in isotopic composition. Secondary, hydrothermal tourmaline in the granite has a bimodal B-isotope distribution with one peak at about −9‰, like the main-stage tourmaline, and a second at −2‰. We propose that the tourmaline-rich orbicules formed late in the crystallization history from an immiscible Na–B–Fe-rich hydrous melt. The massive precipitation of orbicular tourmaline nearly exhausted the melt in boron and the shift of δ¹¹B to −2‰ in secondary tourmaline can be explained by Rayleigh fractionation after about 90% B-depletion in the residual fluid. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

The effects of metamorphism on O and Fe isotope compositions in the Biwabik Iron Formation, northern Minnesota

The Biwabik Iron Formation of Minnesota (1.9 Ga) underwent contact metamorphism by intrusion of the Duluth Complex (1.1 Ga). Apparent quartz–magnetite oxygen isotope temperatures decrease from ∼700°C at the contact to ∼375°C at 2.6 km distance (normal to the contact in 3D). Metamorphic pigeonite at the contact, however, indicates that peak temperatures were greater than 825°C. The apparent O isotope temperatures, therefore, reflect cooling, and not peak metamorphic conditions. Magnetite was reset in δ¹⁸O as a function of grain size, indicating that isotopic exchange was controlled by diffusion of oxygen in magnetite for samples from above the grunerite isograd. Apparent quartz–magnetite O isotope temperatures are similar to calculated closure temperatures for oxygen diffusion in magnetite at a cooling rate of ∼5.6°C/kyr, which suggests that the Biwabik Iron Formation cooled from ∼825 to 400°C in ∼75 kyr at the contact with the Duluth Complex. Isotopic exchange during metamorphism also occurred for Fe, where magnetite–Fe silicate fractionations decrease with increasing metamorphic grade. Correlations between quartz–magnetite O isotope fractionations and magnetite–iron silicate Fe isotope fractionations suggest that both reflect cooling, where the closure temperature for Fe was higher than for O. The net effect of metamorphism on δ¹⁸O–δ⁵⁶Fe variations in magnetite is a strong increase in δ¹⁸O_Mt and a mild decrease in δ⁵⁶Fe with increasing metamorphic grade, relative to the isotopic compositions that are expected at the low temperatures of initial magnetite formation. If metamorphism of Iron Formations occurs in a closed system, bulk O and Fe isotope compositions may be preserved, although re-equilibration among the minerals may occur for both O and Fe isotopes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.