The Kharamai kimberlite field, Siberia: modification of the lithospheric mantle by the Siberian Trap event

The kimberlites of the Kharamai field intruded through the Siberian Traps shortly after their eruption in Permo-Triassic time. The composition and thermal state of the subcontinental lithospheric mantle (SCLM) beneath the Kharamai field in lower Triassic time have been reconstructed using major- and trace-element analyses of 345 Cr-pyrope garnet xenocrysts from six of the kimberlites, supplemented by a small suite of mantle-derived peridotite xenoliths. The data define a geotherm lying near a 38 mW/m² conductive model to a depth of ca 170 km, where the base of the depleted lithosphere is defined by a marked increase in melt-related metasomatism and by an inflected geotherm. Compared to the SCLM sampled by Devonian (pre-Trap) kimberlites in the same and adjacent terranes, the Kharamai SCLM in Triassic time was warmer and was cooling from a previous thermal high. It was also thinner than the SCLM beneath the Daldyn and Alakit kimberlite fields, and had been strongly metasomatised. The metasomatism lowered the mean Fo content of olivine (from ≥Fo₉₃ to Fo₉₂), greatly reduced the proportion of subcalcic harzburgites, and increased the proportion of fertile lherzolites, especially in the depth range of 80–130 km. The overall pattern of metasomatism is similar to that observed in the SCLM sampled by the Group I kimberlites of the SW Kaapvaal Craton, and inferred to be related to the Karoo thermal event. These observations suggest that events such as the eruption of the Karoo basalts and Siberian Traps change the composition of the SCLM, but do not necessarily destroy it, at distances of several hundred kilometres from the main eruption centres.     

Lithosphere mapping beneath the North American plate

Major- and trace-element analyses of garnets from heavy-mineral concentrates have been used to derive the compositional and thermal structure of the subcontinental lithospheric mantle (SCLM) beneath 16 areas within the core of the ancient Laurentian continent and 11 areas in the craton margin and fringing mobile belts. Results are presented as stratigraphic sections showing variations in the relative proportions of different rock types and metasomatic styles, and the mean Fo content of olivine, with depth. Detailed comparisons with data from mantle xenoliths demonstrate the reliability of the sections.

In the Slave Province, the SCLM in most areas shows a two-layer structure with a boundary at 140–160 km depth. The upper layer shows pronounced lateral variations, whereas the lower layer, after accounting for different degrees of melt-related metasomatism, shows marked uniformity. The lower layer is interpreted as a subcreted plume head, added at ca. 3.2 Ga; this boundary between the layers rises to <100 km depth toward the northern and southern edges of the craton. Strongly layered SCLM suggests that plume subcretion may also have played a role in the construction of the lithosphere beneath Michigan and Saskatchewan.

Outside the Slave Province, most North American Archon SCLM sections are less depleted than similar sections in southern Africa and Siberia; this may reflect extensive metasomatic modification. In E. Canada, the degree of modification increases toward the craton margin, and the SCLM beneath the Kapuskasing Structural Zone is typical of that beneath Proterozoic to Phanerozoic mobile belts.

SCLM sections from several Proterozoic areas around the margin of the Laurentian continental core (W. Greenland, Colorado–Wyoming district, Arkansas) show discontinuities and gaps that are interpreted as the effects of lithosphere stacking during collisional orogeny. Some areas affected by Proterozoic orogenesis (Wyoming Craton, Alberta, W. Greenland) appear to retain buoyant, modified Archean SCLM. Possible juvenile Proterozoic SCLM beneath the Colorado Plateau is significantly less refractory. The SCLM beneath the Kansas kimberlite field is highly melt-metasomatised, reflecting its proximity to the Mid-Continent Rift System.

A traverse across the continent shows that the upper part of the cratonic SCLM is highly magnesian; the decrease in mg# with depth is interpreted as the cumulative effect of metasomatic modification through time. The relatively small variations in seismic velocity within the continental core largely reflect the thickness of this depleted layer. The larger drop in seismic velocity in the surrounding Proton and Tecton belts reflects the closely coupled changes in SCLM composition and geotherm.     

Origin of layered continental mantle (Karelian craton, Finland): Geochemical and Re–Os isotope constraints

Studies of mantle xenolith and xenocryst studies have indicated that the subcontinental lithospheric mantle (SCLM) at the Karelian Craton margin (Fennoscandian Shield) is stratified into at least three distinct layers cited A, B, and C. The origin and age of this layering has, however, remained unconstrained. In order to address this question, we have determined Re–Os isotope composition and a comprehensive set of major and trace elements, from xenoliths representing all these three layers. These are the first Re–Os data from the SCLM of the vast East European Craton.

Xenoliths derived from the middle layer B (at  110–180 km depth), which is the main source of harzburgitic garnets and peridotitic diamonds in these kimberlites, are characterised by unradiogenic Os isotopic composition. ¹⁸⁷Os/¹⁸⁸Os shows a good correlation with indices of partial melting implying an age of  3.3. Ga for melt extraction. This age corresponds with the oldest formation ages of the overlying crust, suggesting that layer B represents the unmodified SCLM stabilised during the Paleoarchean. Underlying layer C (at 180–250 km depths) is the main source of Ti-rich pyropes of megacrystic composition but is lacking harzburgitic pyropes. The osmium isotopic composition of layer C xenoliths is more radiogenic compared to layer B, yielding only Proterozoic T_RD ages. Layer C is interpreted to represent a melt metasomatised equivalent to layer B. This metasomatism most likely occurred at ca. 2.0 Ga when the present craton margin formed following continental break-up. Shallow layer A (at  60–110 km depth) has knife-sharp lower contact against layer B indicative of shear zone and episodic construction of SCLM. Layer A peridotites have “ultradepleted” arc mantle-type compositions, and have been metasomatised by radiogenic ¹⁸⁷Os/¹⁸⁸Os, presumably from slab-derived fluids. Since layer A is absent in the core of the craton, its origin can be related to Proterozoic processes at the craton margin. We interpret it to represent the lithosphere of a Proterozoic arc complex (subduction wedge mantle) that became underthrusted beneath the craton margin crust during continental collision  1.9 Ga ago.     

Mesozoic thermal evolution of the southern African mantle lithosphere

The thermal structure of Archean and Proterozoic lithospheric terranes in southern Africa during the Mesozoic was evaluated by thermobarometry of mantle peridotite xenoliths erupted in alkaline magmas between 180 and 60 Ma. For cratonic xenoliths, the presence of a 150–200 °C isobaric temperature range at 5–6 GPa confirms original interpretations of a conductive geotherm, which is perturbed at depth, and therefore does not record steady state lithospheric mantle structure.

Xenoliths from both Archean and Proterozoic terranes record conductive limb temperatures characteristic of a “cratonic” geotherm (40 mW m⁻²), indicating cooling of Proterozoic mantle following the last major tectonothermal event in the region at 1 Ga and the probability of thick off-craton lithosphere capable of hosting diamond. This inference is supported by U–Pb thermochronology of lower crustal xenoliths [Schmitz and Bowring, 2003. Contrib. Mineral. Petrol. 144, 592–618].

The entire region then suffered a protracted regional heating event in the Mesozoic, affecting both mantle and lower crust. In the mantle, the event is recorded at 150 Ma to the southeast of the craton, propagating to the west by 108–74 Ma, the craton interior by 85–90 Ma and the far southwest and northwest by 65–70 Ma. The heating penetrated to shallower levels in the off-craton areas than on the craton, and is more apparent on the southern margin of the craton than in its western interior. The focus and spatial progression mimic inferred patterns of plume activity and supercontinent breakup 30–100 Ma earlier and are probably connected.

Contrasting thermal profiles from Archean and Proterozoic mantle result from penetration to shallower levels of the Proterozoic lithosphere by heat transporting magmas. Extent of penetration is related not to original lithospheric thickness, but to its more fertile character and the presence of structurally weak zones of old tectonism. The present day distribution of surface heat flow in southern Africa is related to this dynamic event and is not a direct reflection of the pre-existing lithospheric architecture.     

Mantle-derived mafic-ultramafic xenoliths and the nature of Indian sub-continental lithosphere

Mantle derived xenoliths in India are known to occur in the Proterozoic ultrapotassic rocks like kimberlites from Dharwar and Bastar craton and Mesozoic alkali igneous rocks like lamrophyres, nephelinites and basanites. The xenoliths in kimberlites are represented by garnet harzburgites, lherzolites, wehrlite, olivine clinopyroxenites and kyaniteeclogite varieties. The PT conditions estimated for xenoliths from the Dharwar craton suggest that the lithosphere was at least 185 km thick during the Mid-Proterozoic period. The ultrabasic and eclogite xenoliths have been derived from depths of 100–180 km and 75–150 km respectively. The Kalyandurg and Brahmanpalle clusters have sampled the typical Archaean subcontinental lithospheric mantle (SCLM) with a low geotherm (35 mW/m²) and harzburgitic to lherzolitic rocks with median X_mg ^olivine > 0.93. The base of the depleted lithosphere at 185–195 km depth is marked by a 10–15 km layer of strongly metasomatised peridotites (X_mg ^olivine > ∼0.88). The Anampalle and Wajrakarur clusters 60 km to the NW show a distinctly different SCLM; it has a higher geotherm (37.5 to 40 mW/m²) and contains few subcalcic harzburgites, and has a median X_mg ^olivine = 0.925. In contrast, the kimberlites of the Uravakonda and WK-7 clusters sampled quite fertile (median X_mg ^olivine ∼0.915) SCLM with an elevated geotherm (> 40 mW/m²). The lamrophyres, basanites and melanephelinites associated with the Deccan Volcanic Province entrain both ultramafic and mafic xenoliths. The ultramafic group is represented by (i) spinel lherzolites, harzburgites, and (ii) pyroxenites. Single pyroxene granulite and two pyroxene granulites constitutes the mafic group. Temperature estimates for the West Coast xenoliths indicate equilibration temperatures of 500–900°C while the pressure estimates vary between 6–11 kbar corresponding to depths of 20–35 km. This elevated geotherm implies that the region is characterized by abnormally high heat flow, which is also supported by the presence of linear array of hot springs along the West Coast. Spinel peridotite xenoliths entrained in the basanites and melanephelinites from the Kutch show low equilibrium temperatures (884–972°C). The estimated pressures obtained on the basis of the absence of both plagioclase and garnet in the xenoliths and by referring the temperatures to the West Coast geotherm is ∼ 15 kbar (40–45 km depth). The minimum heat flow of 60 to 70 mW/m² has been computed for the Kutch xenolith (Bhujia hill), which is closely comparable to the oceanic geotherm. Xenolith studies from the West Coast and Kutch indicate that the SCLM beneath is strongly metasomatised although the style of metasomatism is different from that below the Dharwar Craton.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

The Siberian lithosphere traverse: mantle terranes and the assembly of the Siberian Craton

The kimberlite fields scattered across the NE part of the Siberian Craton have been used to map the subcontinental lithospheric mantle (SCLM), as it existed during Devonian to Late Jurassic time, along a 1000-km traverse NE–SW across the Archean Magan and Anabar provinces and into the Proterozoic Olenek Province. 4100 garnets and 260 chromites from 65 kimberlites have been analysed by electron probe (major elements) and proton microprobe (trace elements). These data, and radiometric ages on the kimberlites, have been used to estimate the position of the local (paleo)geotherm and the thickness of the lithosphere, and to map the detailed distribution of specific rock types and mantle processes in space and time. A low geotherm, corresponding approximately to the 35 mW/m² conductive model of Pollack and Chapman [Tectonophysics 38, 279–296, 1977], characterised the Devonian lithosphere beneath the Magan and Anabar crustal provinces. The Devonian geotherm beneath the northern part of the area was higher, rising to near a 40 mW/m² conductive model. Areas intruded by Mesozoic kimberlites are generally characterised by this higher, but still ‘cratonic' geotherm. Lithosphere thickness at the time of kimberlite intrusion varied from ca. 190 to ca. 240 km beneath the Archean Magan and Anabar provinces, but was less (150–180 km) beneath the Proterozoic Olenek Province already in Devonian time. Thinner Devonian lithosphere (140 km) in parts of this area may be related to Riphean rifting. Near the northern end of the traverse, differences in geotherm, lithosphere thickness and composition between the Devonian Toluopka area and the nearby Mesozoic kimberlite fields suggest thinning of the lithosphere by ca. 50–60 km, related to Devonian rifting and Triassic magmatism. A major conclusion of this study is that the crustal terrane boundaries defined by geological mapping and geophysical data (extended from outcrops in the Anabar Shield) represent major lithospheric sutures, which continue through the upper mantle and juxtapose lithospheric domains that differ significantly in composition and rock-type distribution between 100 and 250 km depth. The presence of significant proportions of harzburgitic and depleted lherzolitic garnets beneath the Magan and Anabar provinces is concordant with their Archean surface geology. The lack of harzburgitic garnets, and the chemistry of the lherzolitic garnets, beneath most of the other fields are consistent with the Proterozoic surface rocks. Mantle sections for different terranes within the Archean portion of the craton show pronounced differences in bulk composition, rock-type distribution, metasomatic overprint and lithospheric thickness. These observations suggest that individual crustal terranes, of both Archean and Proterozoic age, had developed their own lithospheric roots, and that these differences were preserved during the Proterozoic assembly of the craton. Data from kimberlite fields near the main Archean–Proterozoic suture (the Billyakh Shear Zone) suggest that reworking and mixing of Archean and Proterozoic mantle was limited to a zone less than 100 km wide.     

Kimberlites of the Man craton, West Africa

The Man craton in West Africa is an Archaean craton formerly joined to the Guyana craton (South America) that was rifted apart in the Mesozoic. Kimberlites of the Man craton include three Jurassic-aged clusters in Guinea, two Jurassic-aged clusters in Sierra Leone, and in Liberia two clusters of unknown age and one Neoproterozoic cluster recently dated at 800 Ma.

All of the kimberlites irrespective of age occur as small pipes and prolific dykes. Some of the Banankoro cluster pipes in Guinea, the Koidu pipes in Sierra Leone and small pipes in the Weasua cluster in Liberia contain hypabyssal-facies kimberlite and remnants of the so-called transitional-facies and diatreme-facies kimberlite. Most of the Man craton kimberlites are mineralogically classified as phlogopite kimberlites, although potassium contents are relatively low. They are chemically similar to mica-poor Group 1A Southern African examples.

The Jurassic kimberlites are considered to represent one province of kimberlites that track from older bodies in Guinea (Droujba 153 Ma) to progressively younger kimberlites in Sierra Leone (Koidu, 146 Ma and Tongo, 140 Ma). The scarcity of diatreme-facies kimberlites relative to hypabyssal-facies kimberlites and the presence of the so-called transitional-facies indicate that the pipes have been eroded down to the interface between the root and diatreme zones. From this observation, it is concluded that extensive erosion (1–2 km) has occurred since the Jurassic. In addition to erosion, the presence of abundant early crystallizing phlogopite is considered to have had an effect on the relatively small sizes of the Man craton kimberlites.     

Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low ε _Nd and high ⁸⁷Sr/⁸⁶Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low ⁸⁷Sr/⁸⁶Sr and intermediate ε _Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.     

Samples of Proterozoic iron-enriched mantle from the Premier kimberlite

Two populations of mantle xenoliths from the Proterozoic Premier kimberlite show an absence of potassic metasomatism common in Phanerozoic kimberlites. The Premier samples are relatively enriched in Fe and Ti, and contain Fe mica and aluminous amphibole instead of Mg-phlogopite and K-richterite. These features are consistent with a recently identified ρ wave anomaly beneath this part of the Kaapvaal craton ascribed to refertilization of the mantle. Upwelling of sublithospheric mantle to produce the Bushveld Igneous Complex is considered to be the source of silicate melt available for metasomatism. The resultant refertilized Fe-, Ti-, and Al-enriched mantle composition resembles that which is required to form Proterozoic troctolitic magmas.     

The crust-mantle boundary beneath cratons and craton margins: a transect across the south-west margin of the Kaapvaal craton

“Lower-crustal suite” xenoliths occur in “on-craton” and “off-craton” kimberlites located across the south-western margin of the Kaapvaal craton, southern Africa.

Rock types include mafic granulite (plagioclase-bearing assemblages), eclogite (plagioclase-absent assemblages with omphacitic clinopyroxene) and garnet pyroxenite (“orthopyroxene-bearing eclogite”). The mafic granulites are subdivided into three groups: garnet granulites (cpx + grt + plag + qtz); two pyroxene garnet granulites (cpx + opx + grt + plag); kyanite granulites (cpx + grt + ky + plag + qtz). Reaction microstructures preserved in many of the granulite xenoliths involve the breakdown of plagioclase by a combination of reactions: (1) cpx + plag → grt + qtz; (2) plag → grt + ky + qtz; (3) plag → cpx (jd-rich) + qtz. Compositional zoning in minerals associated with these reactions records the continuous transition from granulite facies mineral assemblages and pressure (P) — temperature (T) conditions to those of eclogite facies.

Two distinct P-T arrays are produced: (1) “off-craton” granulites away from the craton margin define a trend from 680 °C, 7.5 kbar to 850 °C, 12 kbar; (2) granulite xenoliths from kimberlites near the craton margin and “on-craton” granulites produce a trend with similar geothermal gradient but displaced to lower T by ˜ 100 °C. Both P-T fields define higher geothermal gradients than the model steady state conductive continental geotherm (40 mWm²) and are not consistent with the paleogeotherm constructed from mantle-derived garnet peridotite xenoliths.

A model involving intrusion of basic magmas around the crust/mantle boundary followed by isobaric cooling is proposed to explain the thermal history of the lower crust beneath the craton margin. The model is consistent with the thermal evolution of the exposed Namaqua-Natal mobile belt low-pressure granulites and the addition of material from the mantle during the Namaqua thermal event (c. 1150 Ma). The xenolith P-T arrays are not interpreted as representing paleogeotherms at the time of entrainment in the host kimberlite. They most likely record P-T conditions “frozen-in” during various stages of the tectonic juxtaposition of the Namaqua Mobile Belt with the Kaapvaal craton.     

Petrology and U–Pb geochronology of lower crustal xenoliths and the development of a craton, Slave Province, Canada

Lower crustal xenoliths recovered from Eocene to Cambrian kimberlites in the central and southern Slave craton are dominated by mafic granulites (garnet, clinopyroxene, plagioclase±orthopyroxene), with subordinate metatonalite and peraluminous felsic granulites. Geothermobarometry indicates metamorphic conditions of 650–800 °C at pressures of 0.9–1.1 GPa. The metamorphic conditions are consistent with temperatures expected for the lower crust of high-temperature low-pressure (HT-LP) metamorphic belts characteristic of Neoarchean metamorphism in the Slave craton. U–Pb geochronology of zircon, rutile and titanite demonstrate a complex history in the lower crust. Mesoarchean protoliths occur beneath the central Slave supporting models of an east-dipping boundary between Mesoarchean crust in the western and Neoarchean crust in the eastern Slave. At least, two episodes of igneous and metamorphic zircon growth occurred in the interval 2.64–2.58 Ga that correlate with the age of plutonism and metamorphism in the upper crust, indicating magmatic addition to the lower crust and metamorphic reworking during this period. In addition, discrete periods of younger zircon growth at ca. 2.56–2.55 and 2.51 Ga occurred 20–70 my after the cessation of ca. 2.60–2.58 Ga regional HT-LP metamorphism and granitic magmatism in the upper crust. This pattern of younger metamorphic events in the deep crust is characteristic of the Slave as well as other Archean cratons (e.g., Superior). The high temperature of the lower crust immediately following amalgamation of the craton, coupled with evidence for continued metamorphic zircon growth for >70 my after ‘stabilization’ of the upper crust, is difficult to reconcile with a thick (200 km), cool lithospheric mantle root beneath the craton prior to this event. We suggest that thick tectosphere developed synchronously or after these events, most likely by imbrication of mantle beneath the craton at or after ca. 2.6 Ga. The minimum age for establishing a cratonic like geotherm is given by lower crustal rutile ages of ca. 1.8 Ga in the southern Slave. Transient heating and possible magmatic additions to the lower crust continued through the Proterozoic, with possible additional growth of the tectosphere.     

Regional patterns in the paragenesis and age of inclusions in diamond, diamond composition, and the lithospheric seismic structure of Southern Africa

The Archean lithospheric mantle beneath the Kaapvaal–Zimbabwe craton of Southern Africa shows ±1% variations in seismic P-wave velocity at depths within the diamond stability field (150–250 km) that correlate regionally with differences in the composition of diamonds and their syngenetic inclusions. Seismically slower mantle trends from the mantle below Swaziland to that below southeastern Botswana, roughly following the surface outcrop pattern of the Bushveld-Molopo Farms Complex. Seismically slower mantle also is evident under the southwestern side of the Zimbabwe craton below crust metamorphosed around 2 Ga. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa, and Jwaneng have Re–Os isotopic ages that range from circa 2.9 Ga to the Proterozoic and show little correspondence with these lithospheric variations. However, silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane do show some regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity correlates with a greater proportion of eclogitic versus peridotitic silicate inclusions in diamond, a greater incidence of younger Sm–Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds whereas the converse is true for diamonds from higher velocity mantle. The oldest formation ages of diamonds indicate that the mantle keels which became continental nuclei were created by middle Archean (3.2–3.3 Ga) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of sulfide inclusions that are eclogitic in the 2.9 Ga age population links late Archean (2.9 Ga) subduction-accretion events involving an oceanic lithosphere component to craton stabilization. These events resulted in a widely distributed younger Archean generation of eclogitic diamonds in the lithospheric mantle. Subsequent Proterozoic tectonic and magmatic events altered the composition of the continental lithosphere and added new lherzolitic and eclogitic diamonds to the already extensive Archean diamond suite.     

Petrology and geochemistry of diamondiferous Mesoproterozoic kimberlites from Wajrakarur kimberlite field,Eastern Dharwar craton,southern India: genesis and constraints on mantle source regions

The petrology and geochemistry of some new occurrences of Mesoproterozoic diamondiferous hypabyssal-facies kimberlites from the Chigicherla, Wajrakarur-Lattavaram and Kalyandurg clusters of the Wajrakarur kimberlite field (WKF), Eastern Dharwar craton (EDC), southern India, are reported. The kimberlites contain two generations of olivine, and multiple groundmass phases including phlogopite, spinel, calcite, dolomite, apatite, perovskite, apatite and rare titanite, and xenocrysts of eclogitic garnet and picro-ilmenite. Since many of the silicate minerals in these kimberlites have been subjected to carbonisation and alteration, the compositions of the groundmass oxide minerals play a crucial role in their characterisation and in understanding melt compositions. While there is no evidence for significant crustal contamination in these kimberlites, some limited effects of ilmenite entrainment are evident in samples from the Kalyandurg cluster. Geochemical studies reveal that the WKF kimberlites are less differentiated and more primitive than those from the Narayanpet kimberlite field (NKF), Eastern Dharwar craton. Highly fractionated (La/Yb = 108–145) chondrite-normalised distribution patterns with La abundances of 500–1,000 × chondrite and low heavy rare earth elements (HREE) abundances of 5–10 × chondrite are characteristic of these rocks. Metasomatism by percolating melts from the convecting mantle, rather than by subduction-related processes, is inferred to have occurred in their source regions based on incompatible element signatures. While the majority of the Eastern Dharwar craton kimberlites are similar to the Group I kimberlites of southern Africa in terms of petrology, geochemistry and Sr–Nd isotope systematics, others show the geochemical traits of Group II kimberlites or an overlap between Group I and II kimberlites. Rare earth element (REE)-based semi-quantitative forward modelling of batch melting of southern African Group I and II kimberlite source compositions involving a metasomatised garnet lherzolite and very low degrees of partial melting demonstrate that (1) WKF and NKF kimberlites display a relatively far greater range in the degree of melting than those from the on-craton occurrences from southern Africa and are similar to that of world-wide melilitites, (2) different degrees of partial melting of a common source cannot account for the genesis of all the EDC kimberlites, (3) multiple and highly heterogeneous kimberlite sources involve in the sub-continental lithospheric mantle (SCLM) in the Eastern Dharwar craton and (4) WKF and NKF kimberlites generation is a resultant of complex interplay between the heterogeneous sources and their different degrees of partial melting. These observations are consistent with the recent results obtained from inversion modelling of REE concentrations from EDC kimberlites in that both the forward as wells as inverse melting models necessitate a dominantly lithospheric, and not asthenospheric, mantle source regions. The invading metasomatic (enriching) melts percolating from the convecting (asthenosphere) mantle impart an OIB-like isotopic signature to the final melt products.     

Integrated models of diamond formation and craton evolution

Two decades of diamond research in southern Africa allow the age, average N content and carbon composition of diamonds, and the dominant paragenesis of their syngenetic silicate and sulfide inclusions to be integrated on a cratonwide scale with a model of craton formation. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa and Jwaneng have Re–Os isotopic ages that range from circa 2.9 Ga to the mid-Proterozoic and display little correspondence with the prominent variations in the P-wave velocity (±1%) that the mantle lithosphere shows at depths within the diamond stability field (150–225 km). Silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane show a regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity relative to the craton average correlates with a greater proportion of eclogitic vs. peridotitic silicate inclusions in diamond, a greater incidence of younger Sm–Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds. The oldest formation ages of diamonds support a model whereby mantle that became part of the continental keel of cratonic nuclei first was created by middle Archean (3.2–3.3 Ga or older) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of eclogitic sulfide inclusions in the 2.9 Ga age population links late Archean (2.9 Ga) subduction–accretion events to craton stabilization. These events resulted in a widely distributed, late Archean generation of eclogitic diamonds in an amalgamated craton. Subsequent Proterozoic tectonic and magmatic events altered the composition of the continental lithosphere and added new lherzolitic and eclogitic diamonds to the already extensive Archean diamond suite. Similar age/paragenesis systematics are seen for the more limited data sets from the Slave and Siberian cratons.     

Isotopic characteristics and petrogenesis of the lamproites and kimberlites of central west Greenland

Kimberlites which intruded the Sisimiut (formerly Holsteinsborg) region of central west Greenland during the Early Palaeozoic have initial ⁸⁷Sr/⁸⁶Sr between 0.7028 and 0.7033 and ε_Nd between + 1.3 and + 3.9. Mid-Proterozoic potassic lamproites from the same region have initial ⁸⁷Sr/⁸⁶Sr between 0.7045 and 0.7060, ε_Nd between −13 and −10 and unradiogenic initial Pb isotopic compositions. The isotopic data favour an asthenospheric mantle source for the kimberlite magmas, in common with “basaltic” kimberlites from other localities, whereas the lamproite magma sources evolved in isolation from the convecting mantle for > 1000 Ma, probably within the subcontinental lithospheric mantle of the Greenland craton, prior to emplacement of the lamproites.     

Imaging global chemical and thermal heterogeneity in the subcontinental lithospheric mantle with garnets and xenoliths: Geophysical implications

Suites of mantle-derived xenoliths in volcanic rocks provide estimates of the geothermal gradient and composition of the subcontinental lithospheric mantle (SCLM) at the time of the volcanic eruption. The development of single-grain thermometry and barometry, applied to xenocryst minerals in volcanic rocks, has greatly expanded the number of localities for which such data can be obtained and made it feasible to map the geology of the SCLM on a broader scale, both vertically and laterally. From garnet xenocrysts, it is possible to derive profiles showing mean values of olivine composition, bulk-rock composition, density and seismic velocities, as well as geotherm parameters and constraints on the thickness of the SCLM. Geochemical profiles, coupled with Re–Os dating of peridotites and their enclosed sulfide minerals, show that Archean or Proterozoic SCLM is preserved at shallow levels beneath many areas of younger tectonothermal age; this implies rapid vertical variations in Vs and Vp with depth, which may affect seismic interpretations. Data from several hundred localities worldwide define a secular evolution in the composition of the SCLM, related to the tectonothermal age of the overlying crust. Archean SCLM is typically strongly depleted in basaltic components, highly magnesian and thick (160–250 km), and has low geotherms; Phanerozoic SCLM is typically fertile (rich in basaltic components), Fe-rich, thin (50–100 km) and has a range of high geotherms; Proterozoic SCLM (much of which may be reworked Archean mantle) tends to be intermediate in all respects. The correlated variations in SCLM fertility, lithospheric thickness and geotherm reinforce the effects of each on seismic velocity, and produce more rapid lateral variations in seismic response than would result from thermal effects alone. These correlations are the key to using seismic tomography images to map the lateral extent of different types of SCLM.     

Peridotitic diamonds from the Slave and the Kaapvaal cratons—similarities and differences based on a preliminary data set

A comparison of the diamond productions from Panda (Ekati Mine) and Snap Lake with those from southern Africa shows significant differences: diamonds from the Slave typically are un-resorbed octahedrals or macles, often with opaque coats, and yellow colours are very rare. Diamonds from the Kaapvaal are dominated by resorbed, dodecahedral shapes, coats are absent and yellow colours are common. The first two features suggest exposure to oxidizing fluids/melts during mantle storage and/or transport to the Earth's surface, for the Kaapvaal diamond population.

Comparing peridotitic inclusions in diamonds from the central and southern Slave (Panda, DO27 and Snap Lake kimberlites) and the Kaapvaal indicates that the diamondiferous mantle lithosphere beneath the Slave is chemically less depleted. Most notable are the almost complete absence of garnet inclusions derived from low-Ca harzburgites and a generally lower Mg-number of Slave inclusions.

Geothermobarometric calculations suggest that Slave diamonds originally formed at very similar thermal conditions as observed beneath the Kaapvaal (geothermal gradients corresponding to 40–42 mW/m² surface heat flow), but the diamond source regions subsequently cooled by about 100–150 °C to fall on a 37–38 mW/m² (surface heat flow) conductive geotherm, as is evidenced from touching (re-equilibrated) inclusions in diamonds, and from xenocrysts and xenoliths. In the Kaapvaal, a similar thermal evolution has previously been recognized for diamonds from the De Beers Pool kimberlites. In part very low aggregation levels of nitrogen impurities in Slave diamonds imply that cooling occurred soon after diamond formation. This may relate elevated temperatures during diamond formation to short-lived magmatic perturbations.

Generally high Cr-contents of pyrope garnets (inside and outside of diamonds) indicate that the mantle lithosphere beneath the Slave originally formed as a residue of melt extraction at relatively low pressures (within the stability field of spinelperidotites), possibly during the extraction of oceanic crust. After emplacement of this depleted, oceanic mantle lithosphere into the Slave lithosphere during a subduction event, secondary metasomatic enrichment occurred leading to strong re-enrichment of the deeper (>140 km) lithosphere. Because of the extent of this event and the occurrence of lower mantle diamonds, this may be related to an upwelling plume, but it may equally just reflect a long term evolution with lower mantle diamonds being transported upwards in the course of “normal” mantle convection.     

Timing of Precambrian melt depletion and Phanerozoic refertilization events in the lithospheric mantle of the Wyoming Craton and adjacent Central Plains Orogen

Garnet peridotite xenoliths from the Sloan kimberlite (Colorado) are variably depleted in their major magmaphile (Ca, Al) element compositions with whole rock Re-depletion model ages generally consistent with this depletion occurring in the mid-Proterozoic. Unlike many lithospheric peridotites, the Sloan samples are also depleted in incompatible trace elements, as shown by the composition of separated garnet and clinopyroxene. Most of the Sloan peridotites have intermineral Sm–Nd and Lu–Hf isotope systematics consistent with this depletion occurring in the mid-Proterozoic, though the precise age of this event is poorly defined. Thus, when sampled by the Devonian Sloan kimberlite, the compositional characteristics of the lithospheric mantle in this area primarily reflected the initial melt extraction event that presumably is associated with crust formation in the Proterozoic—a relatively simple history that may also explain the cold geotherm measured for the Sloan xenoliths.

The Williams and Homestead kimberlites erupted through the Wyoming Craton in the Eocene, near the end of the Laramide Orogeny, the major tectonomagmatic event responsible for the formation of the Rocky Mountains in the late Cretaceous–early Tertiary. Rhenium-depletion model ages for the Homestead peridotites are mostly Archean, consistent with their origin in the Archean lithospheric mantle of the Wyoming Craton. Both the Williams and Homestead peridotites, however, clearly show the consequences of metasomatism by incompatible-element-rich melts. Intermineral isotope systematics in both the Homestead and Williams peridotites are highly disturbed with the Sr and Nd isotopic compositions of the minerals being dominated by the metasomatic component. Some Homestead samples preserve an incompatible element depleted signature in their radiogenic Hf isotopic compositions. Sm–Nd tie lines for garnet and clinopyroxene separates from most Homestead samples provide Mesozoic or younger “ages” suggesting that the metasomatism occurred during the Laramide. Highly variable Rb–Sr and Lu–Hf mineral “ages” for these same samples suggest that the Homestead peridotites did not achieve intermineral equilibrium during this metasomatism. This indicates that the metasomatic overprint likely was introduced shortly before kimberlite eruption through interaction of the peridotites with the host kimberlite, or petrogenetically similar magmas, in the Wyoming Craton lithosphere.     

A Re–Os isotope and PGE study of kimberlite-derived peridotite xenoliths from Somerset Island and a comparison to the Slave and Kaapvaal cratons

The concentrations of platinum-group elements (PGE; Os, Ir, Ru, Pd and Pt) and Re, and the Os isotopic compositions were determined for 33 lithospheric mantle peridotite xenoliths from the Somerset Island kimberlite field. The Os isotopic compositions are exclusively less radiogenic than estimates of bulk-earth (¹⁸⁷Os/¹⁸⁸Os as low as 0.1084) and require a long-term evolution in a low Re–Os environment. Re depletion model ages (T_RD) indicate that the cratonic lithosphere of Somerset Island stabilised by at least 2.8 Ga, i.e. in the Neoarchean and survived into the Mesozoic to be sampled by Cretaceous kimberlite magmatism. An Archean origin also is supported by thermobarometry (Archean lithospheric keels are characterised by >150 km thick lithosphere), modal mineralogy and mineral chemistry observations. The oldest ages recorded in the lithospheric mantle beneath Somerset Island are younger than the Mesoarchean (>3 Ga) ages recorded in the Slave craton lithospheric mantle to the southwest [Irvine, G.J., et al., 1999. Age of the lithospheric mantle beneath and around the Slave craton: a Rhenium–Osmium isotopic study of peridotite xenoliths from the Jericho and Somerset Island kimberlites. Ninth Annual V.M. Goldschmidt Conf., LPI Cont., 971: 134–135; Irvine, G.J., et al., 2001. The age of two cratons: a PGE and Os-Isotopic study of peridotite xenoliths from the Jericho kimberlite (Slave craton) and the Somerset Island kimberlite field (Churchill Province). The Slave–Kaapvaal Workshop, Merrickville, Ontario, Canada]. Younger, Paleoproterozoic, T_RD model ages for Somerset Island samples are generally interpreted as the result of open system behaviour during metasomatic and/or magmatic processes, with possibly the addition of new lithospheric material during tectono-thermal events related to the Taltson–Thelon orogen. PGE patterns highly depleted in Pt and Pd generally correspond to older Archean T_RD model ages indicating closed system behaviour since the time of initial melt extraction. Younger Proterozoic T_RD model ages generally correspond to more complex PGE patterns, indicating open system behaviour with possible sulfide or melt addition. There is no correlation between the age of the lithosphere and depth, at Somerset Island.