Kimberlite Melt Research Articles

Cr-poor and Cr-rich clinopyroxene and garnet megacrysts from southern African Group 1 and Group 2 kimberlites: Clues to megacryst origins and their relationship to kimberlites

Controversies surround the origin of kimberlite megacrysts, including whether and how they are genetically related to their host kimberlites, the relationship between the Cr-poor and Cr-rich suites and the dominant processes responsible for elemental and isotopic variations of megacrysts from a given kimberlite. We present new in-situ major and trace element and Sr isotopic results for clinoyroxene and garnet megacrysts from four southern African kimberlites: Colossus and Orapa (Group 1 kimberlites on the Zimbabwe craton), and Kalkput and Bellsbank (Group 2 kimberlites on the western Kaapvaal craton), that include both Cr-poor and Cr-rich megacryst varieties. Cr-poor megacrysts are present at Colossus, Orapa and Kalkput and the data exhibit tight, well-defined trends on major element diagrams as well as incompatible and rare earth element abundances similar to those previously reported for Cr-poor megacrysts. Cr-rich megacrysts, which are also present at Orapa and are the only variety present at Bellsbank, generally have higher Mg# values, lack well-defined major element trends and show stronger incompatible element enrichments as well as more radiogenic 87Sr/86Sri ratios than Cr-poor megacrysts from the same kimberlite group. Thermobarometry indicates that the Cr-poor megacrysts equilibrated at temperatures of ≈1200 to 1450 °C and pressures of 4.5 to 7.5 GPa. Cr-rich megacrysts, in contrast, extend to temperatures and pressures as low as 700 °C and 3 GPa, respectively. This indicates that, in the studied suites, Cr-poor megacrysts equilibrated at high temperatures in the lower lithosphere (>135 km), whereas Cr-rich megacrysts typically equilibrated at lower temperatures and pressures. Within the Cr-poor megacrysts from Group 1 and Group 2 kimberlites, there is a clear correspondence between kimberlite group, diagnostic incompatible element ratios (e.g., Nb/La) and Sr isotope ratios that parallel the differences noted between whole-rock Group 1 and Group 2 kimberlites. In the case of Cr-poor megacrysts, similar calculated melt compositions in equilibrium with garnet and clinopyroxene from the same kimberlite were obtained using recent high-pressure mineral‑carbonated melt partition coefficients. This suggests formation in conditions close to trace element equilibrium, and is consistent with crystallization from primitive melts with kimberlite-like trace element compositions. In the case of Cr-rich megacrysts, differences in the compositions of melts in equilibrium with clinopyroxene and garnet tend to be larger, and melts in equilibrium with Cr-rich clinopyroxene tend to show significantly greater incompatible element enrichments than those of estimated near-primary kimberlite melts. This could be due to the different behaviour of clinopyroxene and garnet during metasomatic melt-rock interaction, but the apparent disequilibrium between clinopyroxene and garnet could also be due to some of the Cr-rich megacrysts actually being peridotitic xenocrysts. We propose a model for the origin of southern African megacrysts in which carbonated protokimberlite melts formed stockwork-like bodies of variable size in the deep lithosphere (>130 km), which fed networks of melt-filled veins extending into the surrounding and overlying mantle. Crystallization of larger melt bodies resulted in megacryst assemblages dominated by Cr-poor megacrysts, and the incompatible element and isotopic characteristics of these dominantly reflect those of the protokimberlite melt. In contrast, crystallization of smaller melt bodies and their vein networks resulted in megacryst assemblages dominated by Cr-rich megacrysts, which formed as a result of extensive assimilation and metasomatic melt-rock interaction between protokimberlite and peridotite wallrock at low melt/rock ratios, particularly in the middle to shallow lithosphere where pre-existing potassic metasomatic heterogeneities are prevalent. The Cr-rich nature and enrichments in incompatible elements and radiogenic Sr in the Cr-rich megacrysts reflect extensive interaction of their parental magmas with this metasomatized peridotite.

Multistage Mantle Metasomatism during the Generation of Kimberlite Melts: Evidence from Mantle Xenoliths and Megacrysts of the Grib Kimberlite, Arkhangelsk, Russia

Multistage Mantle Metasomatism during the Generation of Kimberlite Melts: Evidence from Mantle Xenoliths and Megacrysts of the Grib Kimberlite, Arkhangelsk, Russia

The effect of oxygen fugacity on diamond resorption in ascending kimberlite melt

When transported by magmas to the Earth's surface, diamond crystals underwent resorption, the intensity of which significantly differed in various kimberlite pipes. We experimentally simulated diamond resorption at different oxygen fugacities (fO2) in ascending kimberlite magma enriched in CO2 and H2O. The experiments were carried out using specially prepared unaltered Group I kimberlite from the Udachnaya East pipe (Yakutia) and model carbonatite at 3.0 GPa, 1200–1400 °C, and fO2 in a range of NNO–2 to NNO + 3.2 log units (where NNO is Ni-NiO buffer). Over the investigated range of conditions, resorption of octahedral diamond crystals is found to occur according to a single scenario. Negative trigons and shield-shaped laminae develop on the {111} faces and crystal edges are truncated by the surfaces of tetrahexahedroids. The rate of diamond resorption increases in all studied systems as fO2 and temperature are raised. In this case, water-enriched melts are the most aggressive media in the investigated T–fO2 interval. Among the most oxidized high-temperature melts, it is carbonatite melts depleted in SiO2 that provide the maximum rate of diamond resorption. Furthermore, the rates of diamond resorption we obtained are an order of magnitude higher than those previously measured in silicate melts containing CO2 and H2O, at fO2 values from the NNO buffer to NNO-2. Therefore, high oxygen fugacity, a temperature of ~1400 °C, and essentially carbonate composition of water-containing magma could provide a high intensity of diamond resorption at the mantle stage of magma ascent to the surface. Apparently, this process primarily influenced the formation of the appearance and preservation of natural diamond crystals in kimberlite pipes.

Magnetotelluric 3-D full tensor inversion in the Dharwar craton, India: Mapping of subduction polarity and kimberlitic melt

Complex geological structures and processes that took place in the Dharwar craton formation make it difficult to understand the evolution history. 3-D magnetotelluric inversion is a challenging task for the imaging of sub-surface structures. Data at 40 stations in a gridded fashion are used in this study for inversion. A controversy exists regarding the subduction polarity between the eastern and western Dharwar craton. Based on the conductivity anomalies mapped in the sub-surface, the lithosphere can be divided into the shallower and deeper lithosphere. The study delineated several crustal and lithospheric upper mantle conductors. In the crustal region, several conductive features (~10 Ω-m) are imaged in the western part, central, and eastern part of the profile. A new finding of this 3-D study is a conductor in the eastern Dharwar craton in the depth range of 65–140 km. The base of this conductor shows the graphite diamond stability field and is correlated with the kimberlites/lamproites present in the region. An uppermost mantle conductor is present at the depth range of 80–200 km in the central part of the study area. Sulphides and carbon-rich fluids could be one cause of the conductors mapped in the crust. The low electrical resistivity imaged in the deeper lithosphere could be due to the refertilization of the mantle scar in the Cretaceous age by the passage of several hotspots. The lithospheric thickness estimated beneath the Dharwar craton in this study is more than 200 km. This study reveals geophysical evidence for the eastward subduction polarity in the Dharwar craton.

The spatial and temporal evolution of primitive melt compositions within the Lac de Gras kimberlite field, Canada: Source evolution vs lithospheric mantle assimilation

Kimberlites are the deepest-derived magmas and provide unique insights into the composition and evolution of the deep Earth. However, this is hindered by contamination of kimberlite melts by mantle and crustal material during ascent. Previous work has shown that the Nd and Hf isotope compositions of the 47 – 75 Ma Lac de Gras (LDG) kimberlites (western Canada) evolve with time from geochemically enriched compositions towards more radiogenic values, which could result from either radiogenic ingrowth or a decreasing contribution of deeply subducted crustal material to the kimberlite source. While previously published TiO2-rich and TiO2-poor whole-rock compositions indicate that at least two distinct primitive melt compositions occur in the LDG field, it is uncertain whether temporal evolution of the kimberlite source influenced primitive melt compositional variability. In this study, the compositions of early formed olivine (rims overgrowths on xenocrysts) and chromite, together with xenocrystic olivine (cores) in 28 LDG kimberlites are combined with new and existing Nd and Hf isotopic compositions and geochronological constraints to examine the temporal and spatial variability of primitive kimberlite melt compositions and their relation to source evolution.Olivine rim Mg# [=100 × Mg/(Mg + Fe2+); 89.6 ± 0.2 to 91.8 ± 0.2] and chromite Ti# [=100 × Ti/(Ti + Al + Cr); 2.5 ± 0.3 to 10.1 ± 1.5], used as proxies for primitive kimberlite melt compositions, indicate that a range of melt compositions were parental to the LDG kimberlites. A positive correlation between olivine rim Mg# and commonly resorbed xenocrystic olivine cores (Mg# = 89.4 ± 1.9 to 92.1 ± 1.3) indicate that assimilation of laterally heterogeneous lithospheric mantle contributed to primitive melt variability. Spatially proximal kimberlites in the Panda, Falcon and Arnie clusters in the centre of the field feature olivine rim Mg# and chromite Ti# with moderate to strong correlations with emplacement ages (47 to 75 Ma), and initial 176Hf/177Hf and 143Nd/144Nd ratios. These correlations indicate that once the effects of lithospheric mantle assimilation are screened out by considering proximal kimberlites, which traversed a similar lithosphere, the compositions of primitive melts at LDG are controlled by the temporal evolution of the sub-lithospheric source. The temporal variability of primitive melt compositions at LDG (i.e., increasingly radiogenic Nd and Hf isotopes and increasing melt Mg# with decreasing age) suggests that the source of these kimberlites evolved via progressive consumption of a deeply subducted, enriched component that was previously mixed with a geochemically more depleted mantle. This work highlights that the compositions of kimberlites are profoundly influenced by both their asthenospheric sources, including deeply subducted material, and assimilation of laterally heterogeneous lithospheric mantle during ascent.

Studying the Stability of the K/Ar Isotopic System of Phlogopites in Conditions of High T, P: 40Ar/39Ar Dating, Laboratory Experiment, Numerical Simulation

Typically, 40Ar/39Ar dating of phlogopites from deep-seated xenoliths of kimberlite pipes produces estimates that suggest much older ages than those when these pipes were intruded. High-pressure (3 GPa) laboratory experiments enabled the authors to explore the behaviour of argon in the phlogopite structure under the conditions that correspond to the mantle, at the temperatures (from 700 to 1000 °С), far exceeding closure temperature of the K/Ar isotopic system. “Volume diffusion” remains foremost for describing the mobility of argon in phlogopite at high pressures. The mantle material age can be estimated through the dating of the phlogopites from deep-seated xenoliths of kimberlites, employing the 40Ar/39Ar method, subject to correction for a partial loss of radiogenic 40Ar when xenolith moves upwards to the Earth’s surface. The obtained data served as the basis for proposing the behaviour model of the K/Ar isotopic system of minerals in conditions of great depths (lower crust, mantle), and when transporting xenoliths in the kimberlite melt.

Silicate Inclusions in Metamorphic Diamonds from the Ultra-High Pressure Kokchetav Complex (Kazakhstan)

Mineral inclusions in cubic diamonds from garnet–clinopyroxene rock of the Kokchetav massif were studied. The coexistence of fluid and silicate inclusions in the central part of the diamond of the G0 sample was revealed by means of transmission electron microscopy. Silicate inclusions are represented by intergrowths of garnet and mica, which are spatially related with the carbonate and fluid inclusions. The first finding of silicate inclusions in the cubic diamonds from the UHP complex discovered over 50 years of their study is apparently due to a selective capture of the silicate minerals in the process of the diamond crystallization from the carbonate-bearing C–O–H fluid. The processes of diamond crystallization in the metamorphic deeply subducted rocks and upper mantle rocks, which are carried to the surface as xenoliths by kimberlite melts, have much in common.

MAGLAB: A computing platform connecting geophysical signatures to melting processes in Earth's mantle

Decompression melting of the upper mantle produces magmas and volcanism at the Earth's surface. Experimental petrology demonstrates that the presence of CO2 and H2O enhances peridotite melting anywhere within the upper mantle down to approximately 200–300 km depth. The presence of mantle melts with compositions ranging from carbonate-rich to silicate-rich unavoidably affects the geophysical signals retrieved from Earth's mantle. Geochemical investigations of erupted intraplate magmas along with geophysical surveys allow for constraining the nature and volume of primary melts, and a sound formalism is required to integrate these diverse datasets into a realistic model for the upper mantle including melting processes. Here, we introduce MAGLAB, a model developed to calculate the composition and volume fraction of melts in the upper mantle, together with the corresponding electrical conductivity of partially molten mantle peridotites at realistic pressure-temperature conditions and volatile contents. We use MAGLAB to show how the compositions of intraplate magmas relate to variations in lithosphere thickness. Progressive partial melting of a homogeneous peridotitic mantle source can in theory create the diversity of compositions observed among the spectrum of intraplate magma types, with kimberlite melts beneath thick continental shields, alkaline magmas such as melilitite, nephelinite and basanite beneath thinner continents and relatively old plus thick oceanic lithospheres, and ‘regular’ basalts beneath the youngest and thinnest oceanic lithospheres as well as beneath significantly thinned continental lithospheres. MAGLAB calculations support recent experimental findings about the role of H2O in the upper mantle on producing primary kimberlitic melts in addition to CO2. We demonstrate the robustness of MAGLAB calculations by reproducing the compositions of erupted melts as well as associated mantle electrical conductivities beneath the Society hotspot in the Pacific Ocean. A comparison of our simulations with magnetotelluric surveys at various oceanic settings shows that the heterogeneities in electrical conductivity of Earth's upper mantle are related to variations in volatile content via the presence of small (generally <<1 wt%) and heterogeneously distributed fractions of CO2-H2O-bearing melts.

Magma Mingling in Kimberlites: Evidence from the Groundmass Cocrystallization of Two Spinel-Group Minerals

We present the results of a detailed petrographic study of fresh coherent samples of the Menominee kimberlite sampled at site 73, located in Menominee County, MI, USA. Our objective is to account for its unusual and complex paragenetic sequence. Several generations of olivine, ilmenite, and spinel-group minerals are described. Early olivine and ilmenite are xenocrystic and were replaced or overgrown by primary minerals. Zoned microcrysts of olivine have a xenocrystic core mantled by a first rim in which rutile, geikielite, and spinel s.s. (spinel sensu stricto) cocrystallized. The in situ U–Pb dating of a microcryst of primary rutile yielded 168.9 ± 4.4 Ma, which was interpreted as the age of emplacement. The groundmass consists of olivine, spinel s.s., a magnesian ulvöspinel–ulvöspinel–magnetite (MUM) spinel, calcite, and dolomite. An extremely low activity of Si is suggested by the crystallization of spinel s.s. instead of phlogopite in the groundmass. The presence of djerfisherite microcrysts indicates high activities of Cl and S during the late stages of melt crystallization. The occurrence of two distinct spinel-group minerals (spinel s.s. and qandilite-rich MUM) in the groundmass is interpreted as clear evidence of the mingling of a magnesiocarbonatitic melt with a dominant kimberlitic melt.

Sr–Nd–O isotopic evidence of variable sources of mantle metasomatism in the subcratonic lithospheric mantle beneath the Grib kimberlite, northwestern Russia

To provide new insights into the type and extent of mantle metasomatism in the subcratonic lithospheric mantle, we examined the Sr–Nd–O isotopic compositions of orthopyroxene, clinopyroxene, garnet, ilmenite and phlogopite from sheared garnet lherzolite, granular garnet harzburgites and lherzolites and clinopyroxene-phlogopite rocks from the Grib kimberlite in the Arkhangelsk diamond province in northwestern Russia.Clinopyroxene and orthopyroxene from sheared garnet lherzolite initially have a close value of 87Sr/86Sr(t) (~0.7034) and close weak positive εNd. Orthopyroxene and clinopyroxene are in oxygen isotope equilibrium with coexisting olivine. Clinopyroxene from a garnet harzburgite has a low 87Sr/86Sr(t) isotope ratio of 0.70266. Clinopyroxene from granular garnet lherzolites has a relatively narrow variation in 87Sr/86Sr(t) (0.70456–0.70582) and considerably larger variations in εNd (−4.3 − +1.0) isotope ratios. Garnet displays elevated initial 87Sr/86Sr(t) values (0.70540–0.70633). Ilmenite shows a narrow range in 87Sr/86Sr(t) (0.70497–0.70522) coupled with εNd values of +0.4 and +3.5. These isotopic data suggest granular garnet lherzolite of mantle metasomatism took place during the interaction of kimberlite melts with SCLM that contained mica-amphibole-rutile-ilmenite-diopside (MARID)-type metasomes. Clinopyroxenes from clinopyroxene-phlogopite (phlogopite wehrlite) xenoliths display a broader range in 87Sr/86Sr(t) (0.70486–0.70813) that is significantly higher than the kimberlite values and a circa-chondritic εNd (−0.1 − +1.3) with a restricted δ18O range (5.11‰–5.33‰). More radiogenic Sr isotopic composition decoupled from Nd isotopes could have been induced by metasomatic melt/fluid related to a subducted material.The isotopic compositions of mantle minerals preserve Sr–Nd isotopic evidence of pre-kimberlite metasomatic events that were probably due to incomplete reequilibration with ultramafic carbonated melt. Based on mineral pairs Rb–Sr isochrons and a clinopyroxene-based Sm–Nd errochron, these mantle metasomatic events correspond to ~550–600 Ma and ~1200 Ma episodes of magmatic–thermal activity.

Thermobarometry of CO2-rich, silica-undersaturated melts constrains cratonic lithosphere thinning through time in areas of kimberlitic magmatism

Cratonic lithosphere is believed to have been chemically buoyant and mechanically resistant to destruction over billions of years. Yet the absence of cratonic roots at some Archean terrains casts doubt on the craton stability and longevity on a global scale. As unique mantle-derived melts at ancient continents, silica-poor, kimberlitic melts are ideal tools to constrain the temporal variation of lithosphere thickness and the processes affecting the lithosphere root. However, no reliable thermobarometer exists to date for strongly silica-undersaturated, mantle-derived melts. Here we develop a new thermobarometer for silica-poor, CO2-rich melts using high-temperature, high-pressure experimental data. Our barometer is calibrated based on a new observation of pressure-dependent variation of Al2O3 in partial melts saturated with garnet and olivine, while our thermometer is calibrated based on the well-known olivine-melt Mg-exchange. For applications to natural magmas, we also establish a correction scheme to estimate their primary melt compositions.Applying this liquid-based thermobarometer to the estimated primary melt compositions for a global kimberlite dataset, we show that the equilibration depths between primary kimberlite melts and mantle peridotites indicate a decrease of up to ∼150 km in cratonic lithosphere thickness globally during the past ∼2 Gyr. Together with the temporal coupling between global kimberlite frequency and cold subduction flux since ∼2 Gyr ago, our results imply a causal link between lithosphere thinning and supply of CO2-rich melts enhanced by deep subduction of carbonated oceanic crusts. While hibernating at the lithosphere root, these melts chemically metasomatize and rheologically weaken the rigid lithosphere and consequently facilitate destruction through convective removal in the ambient mantle or thermo-magmatic erosion during mantle plume activities.

Evolution of kimberlite magmas in the crust: A case study of groundmass and mineral-hosted inclusions in the Mark kimberlite (Lac de Gras, Canada)

Kimberlites are the surface manifestation of deeply-derived (>150 km) and rapidly ascended magmas. Fresh kimberlite rocks are exceptionally rare, as most of them are invariably modified by pervasive deuteric and/or post-magmatic fluids that overprint the original mineralogy. In this study, we examined fresh archetypal kimberlite from the Mark pipe (Lac de Gras, Canada), which is characterised by well-preserved olivine and groundmass minerals. The sequence of crystallisation of the parental melt and its major compositional features, including oxygen fugacity, were reconstructed using textural relationships between magmatic minerals, their zoning patterns and crystal/melt/fluid inclusions. Crystal and multiphase primary, pseudosecondary and secondary melt/fluid inclusions in olivine, Cr-diopside, spinel, perovskite, phlogopite/kinoshitalite, apatite and calcite preserve a record of different stages of kimberlite melt evolution. Melt/fluid inclusions are generally more depleted in silica and more enriched in alkalis (K, Na), alkali-earth (Ba, Sr) and halogens (Cl, F) relative to the whole-rock composition of the Mark kimberlite. These melt/fluid inclusion compositions, in combination with presence of elevated CaO (up to 1.73 wt%), in Mg-rich olivine rinds, crystallisation of groundmass kinoshitalite, carbonates (calcite, Sr-Ba-bearing) and alkali-enriched rims around apatite suggest that there was progressive enrichment in CO2, alkalis and halogens in the evolving parental melt.The Mark kimberlite groundmass is characterised by the following stages of in-situ crystallisation: (1) olivine rims around xenocrystic cores + Cr-spinel/TIMAC. (2) Mg-rich olivine rinds around olivine rims/cores + MUM-spinel (followed by pleonaste and Mg-magnetite) + monticellite (+ partial resorption of olivine, along with the formation of ferropericlase and CO2 as a result of decarbonation reactions) + perovskite + apatite. (3) Olivine outmost rinds, which are coeval with phlogopite/kinoshitalite + apatite + sulphides + carbonate (calcite, Ba-Sr-Na-bearing varieties). In addition, oxygen fugacity of the Mark kimberlite was constrained by olivine-chromite, perovskite and monticellite oxygen barometry and showed that the parental melt became progressively more oxidised in response to fractional crystallisation. (4) Deuteric (i.e. late-stage magmatic) and/or post-magmatic (i.e. external fluids) alteration of magmatic minerals (e.g., olivine, monticellite, ferropericlase) and crystallisation of mesostasis serpentine, K-bearing chlorite and brucite (i.e. replacement of ferropericlase). The absence of any alkali (Na, K) and halogen (F, Cl) rich groundmass minerals in the Mark kimberlite may be attributed to these elements becoming concentrated in the late-stage melt where they potentially formed unstable, water-soluble carbonates (such as those observed in melt inclusions). Consequently, these minerals were most likely removed from the groundmass by deuteric and/or post-magmatic alteration.

Geochemical Evidence for Participation of the Subducted Crust in the Process of Transformation of the Subcontinental Mantle in the Yakutian Diamondiferous Province

The data available indicate the complex evolution of deformed peridotites of mantle xenoliths, the P–T parameters of which indicate that they are fragments of the metasomatized lower part of the cratonic lithosphere. The zoning established in garnets from xenoliths in kimberlite pipes is interpreted as a result of metasomatism that occurred shortly before xenoliths reached the surface. Metasomatic alterations in xenoliths of deformed harzburgites were manifested not only in the development of zoning of minerals. The study results show that there is a discrepancy between the data calculated based on the contents of incompatible elements in minerals of xenoliths and those obtained due to direct measurements of the bulk composition of xenoliths. To determine the balance of incompatible elements, we have carried out experiments on leaching xenoliths of deformed lherzolites from the Udachnaya kimberlite pipe. It was established that a significant part of LREEs in the studied xenoliths occurs in the intergranular space. The distribution pattern of incompatible elements and, in particular, the presence of a positive Eu anomaly indicate that the appearance of the intergranular component is not associated with contamination of xenoliths by the kimberlite melt. Quite a few xenoliths demonstrate a positive Eu anomaly, which indicates the influence of the subducted crustal component at one of the modification stages of xenoliths.

High S and high CO2 contents in haplokimberlite: An experimental and Raman spectroscopic study

Sulfur is an important element present in natural kimberlites and along with CO2, S can play a role in the kimberlite degassing. We have investigated experimentally the change in S content and CO2 solubility in synthetic kimberlitic melts in response to a range of pressure (0.5 to 2.0 GPa) and temperature (1500 to 1525 °C). Several initial S concentrations were investigated ranging from 0 to 24000 ppm. The dissolved CO2 and S were determined by Raman spectroscopy and Electron Probe Micro-Analyses. Under the investigated oxidizing conditions (ΔFMQ + 1), S is dissolved in the glass only as S6+ forming sulfate molecular groups (SO42−). The measured S concentration in the glasses increases from 2900 to 22000 ppm. These results suggest that the experimental conditions were below saturation with respect to S and that the S solubility is higher than 22000 ppm for kimberlitic melts; regardless of the experimental conditions considered here. CO2 is dissolved as CO32− molecular groups. The CO2 solubility ranges from 3.0 to 11.3 wt% between 0.5 and 2.0 GPa. CO2 solubility is not affected by the presence of S; which suggests that SO42− and CO32− clusters have two distinct molecular environments not interacting together. This result implies that both CO2 and S are efficiently transported by kimberlitic melt from the upper mantle towards the atmosphere.

Kimberlite Metasomatism of the Lithosphere and the Evolution of Olivine in Carbonate-rich Melts — Evidence from the Kimberley Kimberlites (South Africa)

AbstractOlivine is the most abundant phase in kimberlites and is stable throughout most of the crystallization sequence, thus providing an extensive record of kimberlite petrogenesis. To better constrain the composition, evolution, and source of kimberlites we present a detailed petrographic and geochemical study of olivine from multiple dyke, sill, and root zone kimberlites in the Kimberley cluster (South Africa). Olivine grains in these kimberlites are zoned, with a central core, a rim overgrowth, and occasionally an external rind. Additional ‘internal’ and ‘transitional’ zones may occur between the core and rim, and some samples of root zone kimberlites contain a late generation of high-Mg olivine in cross-cutting veins. Olivine records widespread pre-ascent (proto-)kimberlite metasomatism in the mantle including the following features: (1) relatively Fe-rich (Mg# &amp;lt;89) olivine cores interpreted to derive from the disaggregation of kimberlite-related megacrysts (20 % of cores); (2) Mg–Ca-rich olivine cores (Mg# &amp;gt;89; &amp;gt;0·05 wt% CaO) suggested to be sourced from neoblasts in sheared peridotites (25 % of cores); (3) transitional zones between cores and rims probably formed by partial re-equilibration of xenocrysts (now cores) with a previous pulse of kimberlite melt (i.e. compositionally heterogeneous xenocrysts); (4) olivine from the Wesselton water tunnel sills, internal zones (I), and low-Mg# rims, which crystallized from a kimberlite melt that underwent olivine fractionation and stalled within the shallow lithospheric mantle. Magmatic crystallization begins with internal olivine zones (II), which are common but not ubiquitous in the Kimberley olivine. These zones are euhedral, contain rare inclusions of chromite, and have a higher Mg# (90·0 ± 0·5), NiO, and Cr2O3 contents, but are depleted in CaO compared with the rims. Internal olivine zones (II) are interpreted to crystallize from a primitive kimberlite melt during its ascent and transport of olivine toward the surface. Their compositions suggest assimilation of peridotitic material (particularly orthopyroxene) and potentially sulfides prior to or during crystallization. Comparison of internal zones (II) with liquidus olivine from other mantle-derived carbonate-bearing magmas (i.e. orangeites, ultramafic lamprophyres, melilitites) shows that low (100×) Mn/Fe (∼1·2), very low Ca/Fe (∼0·6), and moderate Ni/Mg ratios (∼1·1) appear to be the hallmarks of olivine in melts derived from carbonate-bearing garnet-peridotite sources. Olivine rims display features indicative of magmatic crystallization, which are typical of olivine rims in kimberlites worldwide; that is, primary inclusions of chromite, Mg-ilmenite and rutile, homogeneous Mg# (88·8 ± 0·3), decreasing Ni and Cr, and increasing Ca and Mn. Rinds and high-Mg olivine are characterized by extreme Mg–Ca–Mn enrichment and Ni depletion, and textural relationships indicate that these zones represent replacement of pre-existing olivine, with some new crystallization of rinds. These zones probably precipitated from evolved, oxidized, and relatively low-temperature kimberlite fluids after crustal emplacement. In summary, this study demonstrates the utility of combined petrography and olivine geochemistry to trace the evolution of kimberlite magmatic systems from early metasomatism of the lithospheric mantle by (proto-)kimberlite melts, to crystallization at different depths en route to surface, and finally late-stage deuteric or hydrothermal fluid alteration after crustal emplacement.

Isotopic analyses of clinopyroxenes demonstrate the effects of kimberlite melt metasomatism upon the lithospheric mantle

The trace element and radiogenic isotope systematics of clinopyroxene have frequently been used to characterise mantle metasomatic processes, because it is the main host of most lithophile elements in the lithospheric mantle. To further our understanding of mantle metasomatism, both solution-mode Sr-Nd-Hf-Pb and in situ trace element and Sr isotopic data have been acquired for clinopyroxene grains from a suite of peridotite (lherzolites and wehrlites), MARID (Mica-Amphibole-Rutile-Ilmenite-Diopside), and PIC (Phlogopite-Ilmenite-Clinopyroxene) rocks from the Kimberley kimberlites (South Africa). The studied mantle samples can be divided into two groups on the basis of their clinopyroxene trace element compositions, and this subdivision is reinforced by their isotopic ratios. Type 1 clinopyroxene, which comprises PIC, wehrlite, and some sheared lherzolite samples, is characterised by low Sr (~100–200 ppm) and LREE concentrations, moderate HFSE contents (e.g., ~40–75 ppm Zr; La/Zr < 0.04), and restricted isotopic compositions (e.g., 87Sr/86Sri = 0.70369–0.70383; εNdi = +3.1 to +3.6) resembling those of their host kimberlite magmas. Available trace element partition coefficients can be used to show that Type 1 clinopyroxenes are close to being in equilibrium with kimberlite melt compositions, supporting a genetic link between kimberlites and these metasomatised lithologies. Thermobarometric estimates for Type 1 samples in this study indicate equilibration depths of 135–160 km within the lithosphere, thus showing that kimberlite melt metasomatism is prevalent in the deeper part of the lithosphere beneath Kimberley. In contrast, Type 2 clinopyroxenes occur in MARID rocks and coarse granular lherzolites in this study, which derive from shallower depths (<135 km), and have higher Sr (~350–1000 ppm) and LREE contents, corresponding to higher La/Zr of > ~ 0.05. The isotopic compositions of Type 2 clinopyroxenes are more variable and extend from compositions resembling the “enriched mantle” towards those of Type 1 rocks (e.g., εNdi = −12.7 to −4.4). To constrain the source of these variations, in situ Sr isotope analyses of clinopyroxene were undertaken, including zoned grains in Type 2 samples. MARID and lherzolite clinopyroxene cores display generally radiogenic but variable 87Sr/86Sri values (0.70526–0.71177), which are correlated with Sr contents and La/Zr ratios, and which might be explained by the interaction between peridotite and melts from different enriched sources within the lithospheric mantle. Most notably, the rims of these Type 2 clinopyroxenes trend towards compositions similar to those of the host kimberlite and Type 1 clinopyroxene from PIC and wehrlites. These results are interpreted to represent clinopyroxene overgrowth during late-stage (shortly before/during entrainment) metasomatism by kimberlite magmas. Our study shows that a pervasive, alkaline metasomatic event caused MARID to be generated and harzburgites to be converted to lherzolite in the lithospheric mantle beneath the Kimberley area, which was followed by kimberlite metasomatism during Cretaceous magmatism. This latter event is the time at which discrete PIC, wehrlite, and sheared lherzolite lithologies were formed, and MARID and granular lherzolites were partly modified.

Reconciling seismic structures and Late Cretaceous kimberlite magmatism in northern Alberta, Canada

AbstractThe Late Cretaceous kimberlites in northern Alberta, Canada, intruded into the Paleoproterozoic crust and represent a nonconventional setting for the discovery of diamonds. Here, we examined the origin of kimberlite magmatism using a multidisciplinary approach. A new teleseismic survey reveals a low-velocity (−1%) corridor that connects two deep-rooted (&amp;gt;200 km) quasi-cylindrical anomalies underneath the Birch Mountains and Mountain Lake kimberlite fields. The radiometric data, including a new U-Pb perovskite age of 90.3 ± 2.6 Ma for the Mountain Lake intrusion, indicate a northeast-trending age progression in kimberlite magmatism, consistent with the (local) plate motion rate of North America constrained by global plate reconstructions. Taken together, these observations favor a deep stationary (relative to the lower mantle) source region for kimberlitic melt generation. Two competing models, mantle plume and slab subduction, can satisfy kinematic constraints and explain the exhumation of ultradeep diamonds. The plume hypothesis is less favorable due to the apparent age discrepancy between the oldest kimberlites (ca. 90 Ma) and the plume event (ca. 110 Ma). Alternatively, magma generation may have been facilitated by decompression of hydrous phases (e.g., wadsleyite and ringwoodite) within the mantle transition zone in response to thermal perturbations by a cold slab. The three-dimensional lithospheric structures largely controlled melt migration and intrusion processes during the Late Cretaceous kimberlite magmatism in northern Alberta.

Apatite compositions and groundmass mineralogy record divergent melt/fluid evolution trajectories in coherent kimberlites caused by differing emplacement mechanisms

Kimberlites are pipe-like igneous bodies, consisting of a pyroclastic crater and diatreme, commonly underlain by coherent root-zone rocks, and with associated dyke/sill complexes. The processes that control the different modes of coherent kimberlite emplacement remain uncertain. In addition, late evolution of kimberlite melts during emplacement into the upper crust remains poorly constrained. Therefore, it is unclear whether there is a link between melt composition/evolution and the emplacement mechanism of coherent kimberlites (i.e. planar dykes/sills vs. irregular bodies in the root zone). An absence of comparative studies on late-stage magmatic phases across the different emplacement modes of coherent kimberlite from the same locality hamper resolution of these issues. Therefore, we report petrographic and mineral chemical data for groundmass apatite in samples of dyke, sill, and root-zone kimberlites from the Kimberley cluster (South Africa). Early crystallised phases (olivine, spinel, Mg-ilmenite) in dyke/sill and root-zone kimberlites have indistinguishable compositions, and hence crystallised from similar primitive melts. Conversely, apatite compositions are generally distinct in dyke/sill (low Sr, high and variable Si) and root-zone kimberlites (high and variable Sr, low Si). The Si enrichment of apatite in dykes/sills is attributed to the coupled incorporation of CO32− and SiO44− for PO43−, reflecting higher CO2 contents in their parental melts, and potentially higher Si contents due to the preferential crystallisation of carbonates over mica/monticellite. The low Sr contents of apatite in dyke/sill kimberlites reflect equilibrium with a (kimberlite) melt (i.e. DSr is close to unity for carbonate and silicate melts), whereas the higher Sr contents of apatite in root-zone kimberlites require crystallisation from, or overprinting by a H2O ± CO2 fluid (significantly higher DSr). The relative enrichment of CO2 in kimberlite dykes/sills is evident from the abundance of carbonates, the presence of mesostasis dolomite and calcite phenocrysts in some samples, and concomitant reduced proportions of other groundmass phases (e.g. serpentine, mica, monticellite). During late alteration of kimberlite dykes/sills, monticellite is typically replaced by carbonates, whereas olivine and pleonaste are relatively stable, indicating the melts which form dykes/sills evolve to higher CO2/H2O ratios. It is unlikely that these two distinct evolutionary paths were caused by crustal contamination before or during near surface magma emplacement, because crustal assimilation is not recorded in the O and Sr isotopic composition of late crystallising olivine rinds or carbonates, respectively. We suggest that higher concentrations of CO2 are retained in kimberlite dykes/sills due to higher confining pressures (i.e. lack of breakthrough to the surface). In contrast, exsolution of CO2 from root-zone kimberlites increased melt H2O/CO2 ratios and promoted the crystallisation of mica and monticellite at the expense of dolomite and calcite. Apatite compositions have the potential to aid in the discrimination of kimberlites from lamproites (higher LREE, Sr, F, and S, lower Si contents) and carbonatites (higher LREE, F, Cl and S, lower Fe contents). However, the compositions of kimberlitic apatite overlap those from aillikites, probably due to similar late-stage melt compositions.

Metasomatic Evolution of Coesite-Bearing Diamondiferous Eclogite from the Udachnaya Kimberlite

A coesite-bearing diamondiferous eclogite from the Udachnaya kimberlite (Daldyn field, Siberian craton) has been studied to trace its complex evolution recorded in rock-forming and minor mineral constituents. The eclogite sample is composed of rock-forming omphacite (60 vol%), garnet (35 vol%) and quartz/coesite (5 vol%) and contains intergranular euhedral zoned olivine crystals, up to 200 µm long, coexisting with phlogopite, orthopyroxene, clinopyroxene (secondary), K-feldspar, plagioclase, spinel, sodalite and djerfisherite. Garnet grains are zoned, with a relatively homogeneous core and a more magnesian overgrowth rim. The rim zones further differ from the core in having higher Zr/Y (6 times that in the cores), ascribed to interaction with, or precipitation from, a kimberlite-related melt. Judging by pressure-temperature estimates (~1200 °C; 6.2 GPa), the xenolith originated at depths of ~180–200 km at the base of the continental lithosphere. The spatial coexistence of olivine, orthopyroxene and coesite/quartz with K-Na-Cl minerals in the xenolith indicates that eclogite reacted with a deep-seated kimberlite melt. However, Fe-rich olivine, orthopyroxene and low-pressure minerals (sodalite and djerfisherite) likely result from metasomatic reaction at shallower depths during transport of the eclogite by the erupting kimberlite melt. Our results demonstrate that a mixed eclogitic-peridotitic paragenesis, reported previously from inclusions in diamond, can form by interaction of eclogite and a kimberlite-related melt.

Multi-stage evolution of kimberlite melt as inferred from inclusions in garnet megacrysts in the Grib kimberlite (Arkhangelsk region, Russia)

To provide new insights into the origin of garnet megacrysts and evolution of kimberlite melts, we studied in detail the polymineralic and monomineralic inclusions and their host garnets from the Grib kimberlite (Arkhangelsk diamond province, Russia). Low-Cr and high-Cr garnet megacrysts and eclogitic garnets contain abundant polymineralic and rare monomineralic inclusions. Monomineralic inclusions presented by clinopyroxene, ilmenite, olivine replaced by serpentine, were found exclusively within the low-Cr megacrysts. The composition of clinopyroxene exhibits geochemical equilibrium with the host garnet, indicating its primary origin during the formation of the megacryst assemblage. The low-Cr garnet–clinopyroxene mineral assemblage formed as a result of high-temperature, melt-associated mantle metasomatism by failed kimberlite within the lithospheric mantle (T = 1150 °C and P = 5.5 GPa). Polymineralic inclusions are characterised by a silicate or silicate-sulphate matrix. The central part of the silicate inclusions is filled by serpentine and contains ilmenite, spinel, perovskite, calcite and apatite. At the contact with host garnets, phlogopite, spinel and amphibole occur as reaction minerals. Composition of spinel and other minerals within inclusions with silicate matrix suggests that kimberlite melt was trapped at mantle pressures. Inclusions with silicate matrix were found in all garnets. The matrix of silicate-sulphate inclusions consists of silicate cryptocrystalline phases and sulphate minerals (celestine–barite) and contains calcite grains. The inclusions are distributed in some low-Cr garnet megacrysts and eclogitic garnet. The silicate-sulphate inclusions were crystallised from the late-stage kimberlite melt. Diverse reaction textures are evidences of disequilibrium between the host crystals and polymineralic inclusions and indicate that garnet and the hosted inclusions reacted with the ascending kimberlite melt. The silicate-sulphate inclusions with a thin rim of epidote within eclogitic garnets indicate that a kimberlite melt invaded the garnet and induced partial melting. The studied inclusions allow us to propose three stages of the Grib kimberlite evolution: 1) generation of garnet megacrysts and primary inclusions due to melt metasomatism, 2) reaction of the high-Ti kimberlite melt with garnet megacrysts (including their dissolution) and 3) alteration of the inclusions in garnet after kimberlite ascent.