The timing of kimberlite magmatism in North America: implications for global kimberlite genesis and diamond exploration

Based on a compilation of more than 100 kimberlite age determinations, four broad kimberlite emplacement patterns can be recognized in North America: (1) a northeast Eocambrian/Cambrian Labrador Sea province (Labrador, Québec), (2) an eastern Jurassic province (Ontario, Québec, New York, Pennsylvania), (3) a Cretaceous central corridor (Nunavut, Saskatchewan, central USA), and (4) a western mixed (Cambrian-Eocene) Type 3 kimberlite province (Alberta, Nunavut, Northwest Territories, Colorado/Wyoming). Ten new U–Pb perovskite/mantle zircon and Rb–Sr phlogopite age determinations are reported here for kimberlites from the Slave and Wyoming cratons of western North America. Within the Type 3 Slave craton, at least four kimberlite age domains exist: I-a southwestern Siluro-Ordovician domain (450 Ma), II-a SE Cambrian domain (540 Ma), III-a central Tertiary/Cretaceous domain (48–74 Ma) and IV-a northern mixed domain consisting of Jurassic and Permian kimberlite fields. New U–Pb perovskite results for the 614.5±2.1 Ma Chicken Park and 408.4±2.6 Ma Iron Mountain kimberlites in the State Line field in Colorado and Wyoming confirm the existence of at least two periods of pre-Mesozoic kimberlite magmatism in the Wyoming craton.

A compilation of robust kimberlite emplacement ages from North America, southern Africa and Russia indicates that a high proportion of known kimberlites are Cenozoic/Mesozoic. We conclude that a majority of these kimberlites were generated during enhanced mantle plume activity associated with the rifting and eventual breakup of the supercontinent Gondwanaland. Within this prolific period of kimberlite activity, there is a good correlation between North America and Yakutia for three distinct short-duration (10 my) periods of kimberlite magmatism at 48–60, 95–105 and 150–160 Ma. In contrast, Cenozoic/Mesozoic kimberlite magmatism in southern Africa is dominated by a continuum of activity between 70–95 and 105–120 Ma with additional less-prolific periods of magmatism in the Eocene (50–53 Ma), Jurassic (150–190) and Triassic (235 Ma). Several discrete episodes of pre-Mesozoic kimberlite magmatism variably occur in North America, southern Africa and Yakutia at 590–615, 520–540, 435–450, 400–410 and 345–360 Ma. One of the surprises in the timing of kimberlite magmatism worldwide is the common absence of activity between about 250 and 360 Ma; this period is even longer in southern Africa. This >110 my period of quiescence in kimberlite magmatism is likely linked to relative crustal and mantle stability during the lifetime of the supercontinent Gondwanaland.

Economic diamond deposits in kimberlite occur throughout the Phanerozoic from the Cambrian (Venetia, South Africa; Snap Lake and Kennady Lake, Canada) to the Tertiary (Mwadui, Tanzania; Ekati and Diavik in Lac de Gras, Canada). There are clearly some discrete periods when economic kimberlite-hosted diamond deposits formed globally. In contrast, the Devonian event, which is such an important source of diamonds in Yakutia, is notably absent in the kimberlite record from both southern Africa and North America.     

Temporal, geomagnetic and related attributes of kimberlite magmatism at Ekati, Northwest Territories, Canada

This paper outlines the development of a multi-disciplinary strategy to focus exploration for economic kimberlites on the Ekati property. High-resolution aeromagnetic data provide an over-arching spatial and magnetostratigraphic framework for exploration and kimberlite discovery at Ekati, and hence also for this investigation. The temporal, geomagnetic, spatial and related attributes of kimberlites with variable diamond content have been constrained by judiciously augmenting the information gathered during routine exploration with detailed, laboratory-based or field-based investigations. The natural remanent magnetisation of 36 Ekati kimberlites has been correlated with their age as determined by isotopic dating techniques, and placed in the context of a well-constrained geomagnetic polarity timescale. Kimberlite magmatism occurred over the period 75 to 45 Ma, in at least five temporally discrete intrusive episodes. Based on current evidence, the older kimberlites (75 to 59 Ma) have low diamond contents and are distributed throughout the property. Younger kimberlites (56 to 45 Ma) have moderate to high diamond contents and occur in three distinct intrusive corridors with NNE to NE orientations. Economic kimberlite pipes erupted at 55.4±0.4 Ma along the A154-Lynx intrusive corridor, which is 7 km wide and oriented at 015°, and at 53.2±0.3 Ma along the Panda intrusive corridor, which is 1 km wide and oriented at 038°. The intrusion ages straddle a paleopole reversal at Chron C24n, consistent with the observation that the older economic kimberlites present as aeromagnetic “low” anomalies while the younger economic pipes are characterised as aeromagnetic “highs”. The aeromagnetic responses for these kimberlites are generally muted because they contain volcaniclastic rock types with low magnetic susceptibility. Kimberlites throughout the Ekati property carry a primary natural magnetic remanence (NRM) vector in Ti-bearing groundmass magnetite, and it dominates over vectors related to induced magnetisation. Magnetostratigraphic correlation of Ekati kimberlites may therefore present a powerful adjunct to existing exploration techniques, mainly because the diamond content of Ekati kimberlites apparently is related more to the age of eruption than to any other parameter investigated in this work.     

The temporal evolution of North American kimberlites

North American kimberlite magmatism spans a period of time in excess of 1 billion years from Mesoproterozoic kimberlites in the Lake Superior and James Bay Lowlands region of Ontario to Eocene kimberlites in the Lac de Gras field, N.W.T. Based on a compilation of more than 150 robust radiometric age determinations, several distinct kimberlite emplacement patterns are recognized. In general, the temporal pattern of kimberlite emplacement in North America can be broadly subdivided into five domains: (1) a Mesoproterozoic kimberlite province in central Ontario, (2) an Eocambrian/Cambrian Labrador Sea Province in northern Québec and Labrador, (3) an eastern Jurassic Province, (4) a central Cretaceous corridor and (5) a western mixed domain that includes two Type-3 kimberlite provinces (i.e. multiple periods of kimberlite emplacement preserved in the Slave and Wyoming cratons). For some provinces the origin of kimberlite magmatism can be linked to known mantle heat sources such as mantle plume hotspots and upwelling asthenosphere attendant with continental rifting. For example, the timing and location of Mesoproterozoic kimberlites in North America coincides with and slightly precedes the timing of 1.1 Ga intracontinental rifting that culminated in the Midcontinent Rift centered in the Lake Superior region. Many of the kimberlites in the Eocambrian/Cambrian Labrador Sea province were emplaced soon after the opening of the Iapetus Ocean at about 615 Ma and may also be linked to mantle upwelling associated with continental rifting. The eastern Jurassic kimberlites record an age progression where magmatism youngs in a southeast direction from the 200 Ma Rankin Inlet kimberlites to the 155–126 Ma Timiskaming kimberlites. The location of several kimberlite fields and clusters in Ontario and Québec lie along a continental extension of the Great Meteor hotspot track and represents one of the best examples in the world of kimberlite magmatism triggered by mantle plumes. The central Cretaceous (103–94 Ma) corridor extends for more than 4000 km from Somerset Island in northern Canada through the Fort à la Corne field in Saskatchewan to the kimberlites in central USA. This is the first recognized corridor of kimberlite magmatism of this magnitude. The possible westward younging of Cretaceous to Eocene corridors of kimberlite magmatism could reflect major changes in plate geometry during subduction of the Kula–Farallon plate.     

Space–time distribution of interplate moment release including slow earthquakes and the seismo-geodetic coupling in the Sanriku-oki region along the Japan trench

Along the Japan trench where some Mw8 class interplate earthquakes occurred in the past century such as the 1896 Sanriku tsunami earthquake (M6.8, Mt8.6, 12×10²⁰ N m) and the 1968 Tokachi-oki earthquake (Mw8.2, 28×10²⁰ N m), the Pacific plate is subducting under northeast Japan at a rate of around 8 cm/year. The seismic coupling coefficient in this region has been estimated to be 20–40%. In the past decade, three ultra-slow earthquakes have occurred in the Sanriku-oki region (39°N–42°N): the 1989 Sanriku-oki (Mw7.4), the 1992 Sanriku-oki (Mw6.9), and the 1994 Sanriku-oki (Mw7.7) earthquakes. Integrating their interplate moments released both seismically and aseismically, we have the following conclusions. (1) The sum of the seismic moments of the three ultra-slow earthquakes was (4.8–6.6)×10²⁰ N m, which was 20–35% of the accumulated moment (18.6–23.0)×10²⁰ N m, in the region (39°N–40.6°N, 142°E–144°E) for the 21–26 years since the 1968 Mw8.2 Tokachi-oki earthquake. This is consistent with the previous estimates of the seismic coupling coefficient of 20–40%. On the other hand, the sum of the interplate moments including aseismic faulting is (11–16)×10²⁰ N m, leading to a “seismo-geodetic coupling coefficient” of 50–85%, which is an extension of the seismic coupling coefficient to include slow events. (2) The time constants showed a large range from 1 min (10² s) for the 1968 Tokachi-oki earthquake to 10–20 min (10³ s) for the 1896 Sanriku tsunami earthquake, to one day (10⁵ s) for the 1992 Sanriku-oki ultra-slow earthquake, to on the order of one year (10⁷ s) for the 1994 Sanriku-oki ultra-slow earthquakes. (3) Based on the space–time distribution, three “gaps of moment release,” (40.6°N–42°N, 142°E–144°E) 39°N–40°N, 142°E–143°E) and (39°N–40°N, 142°E–144°E), are identified, instead of the gaps of seismicity.     

Macrocrystal phlogopite Rb–Sr dates for the Ekati property kimberlites, Slave Province, Canada: evidence for multiple intrusive episodes in the Paleocene and Eocene

New Rb–Sr age determinations using macrocrystal phlogopite are presented for 27 kimberlites from the Ekati property of the Lac de Gras region, Slave Province, Canada. These new data show that kimberlite magmatism at Ekati ranges in age from at least Late Paleocene (61 Ma) to Middle Eocene time (45 Ma). Older, perovskite-bearing kimberlites from Ekati extend this age range to Late Cretaceous time (74 Ma). Within this age range, emplacement episodes at 48, 51–53, 55–56 and 59–61 Ma can be recognized. Middle Eocene kimberlite magmatism of the previously dated Mark kimberlite (47.5 Ma) is shown to include four other pipes from the east-central Ekati property. A single kimberlite (Aaron) may be younger than the 47.5 Ma Mark kimberlite. The economically important Panda kimberlite is precisely dated in this study to be 53.3±0.6 Ma using the phlogopite isochron method, and up to six additional kimberlites from the central Ekati property have Early Eocene ages indistinguishable from that of Panda, including the Koala and Koala North occurrences. Late Paleocene 55–56 Ma kimberlite magmatism, represented by the Diavik kimberlite pipes adjacent to the southeastern Ekati property, is shown to extend onto the southeastern Ekati property and includes three, and possibly four, kimberlites. A precise eight-point phlogopite isochron for the Cobra South kimberlite yields an emplacement age of 59.7±0.4 Ma; eight other kimberlites from across the Ekati property have similar Late Paleocene Rb–Sr model ages. The addition of 27 new emplacement ages for kimberlites from the Ekati property confirms that kimberlite magmatism from the central Slave Province is geologically young, despite ages ranging back to Cambrian time from elsewhere in the Slave Province. With the available geochronologic database, Lac de Gras kimberlites with the highest diamond potential are currently restricted to the 51–53 and 55–56 Ma periods of kimberlite magmatism.     

Age and emplacement of late-Variscan granites of the western Bohemian Massif with main focus on the Hauzenberg granitoids (European Variscides, Germany)

The Variscan Hauzenberg pluton consists of granite and granodiorite that intruded late- to postkinematically into HT-metamorphic rocks of the Moldanubian unit at the southwestern margin of the Bohemian Massif (Passauer Wald). U–Pb dating of zircon single-grains and monazite fractions, separated from medium- to coarse-grained biotite-muscovite granite (Hauzenberg granite II), yielded concordant ages of 320 ± 3 and 329 ± 7 Ma, interpreted as emplacement age. Zircons extracted from the younger Hauzenberg granodiorite yielded a ²⁰⁷Pb–²⁰⁶Pb mean age of 318.6 ± 4.1 Ma. The Hauzenberg granite I has not been dated. The pressure during solidification of the Hauzenberg granite II was estimated at 4.6 ± 0.6 kbar using phengite barometry on magmatic muscovite, corresponding to an emplacement depth of 16-18 km. The new data are compatible with pre-existing cooling ages of biotite and muscovite which indicate the Hauzenberg pluton to have cooled below T = 250–400 °C in Upper Carboniferous times. A compilation of age data from magmatic and metamorphic rocks of the western margin of the Bohemian Massif suggests a west- to northwestward shift of magmatism and HT/LP metamorphism with time. Both processes started at > 325 Ma within the South Bohemian Pluton and magmatism ceased at ca. 310 Ma in the Bavarian Oberpfalz. The slight different timing of HT metamorphism in northern Austria and the Bavarian Forest is interpreted as being the result of partial delamination of mantle lithosphere or removal of the thermal boundary layer.     

Temperature–time paths from phosphate accessory phase paragenesis in the Honey Brook Upland and associated cover sequence, SE Pennsylvania, USA

Analysis of monazite-bearing lithologies from the Precambrian Honey Brook Upland (HBU) and overlying metasedimentary Paleozoic Chester Valley Sequence (CVS) (SE PA, USA) reveals overprinting of primary major and accessory phase parageneses by texturally and compositionally disparate secondary accessory phase parageneses. Two-pyroxene temperatures of 915–945 °C for reconstituted pyroxene reflect emplacement temperatures of felsic plutonic rocks (opdalite, charnockite) prior to Mesoproterozoic metamorphism. Monazite in metavolcanic felsic gneiss yields three age domains at 1009 ± 4 Ma (2 s.e.), 965 ± 6, and 876 ± 10 Ma. The first two domains record metamorphism of the HBU after anorthosite intrusion; peak monazite–xenotime temperatures for the monazite core domain are 700 °C, and high Th/U values in the second (overgrowth) age domain likely reflect a second high-T monazite growth episode. Formation of cummingtonite coronas on orthopyroxene in opdalite constrains maximum 1010 Ma metamorphic temperatures in the “granulite-facies” terrane to 730–740 °C. Evidence of increased Cl fluid activity in the 965 Ma metamorphism includes higher Cl content of matrix apatite relative to garnet-included apatite (metavolcanics), and Cl-bearing K-hornblende succeeding cummingtonite in coronal overgrowths (opdalite). Extreme monazite Th/U values (75–250) in the rim domain suggest growth during low-T hydrothermal alteration. In the opdalite, secondary singe-grain monazite and monazite + xenotime metasomites in apatite yield ages of 714 ± 24 and 586 ± 88 Ma, temperatures of 325–425 °C, and are interpreted to reflect thermal disturbances associated with late Proterozoic plutonic and volcanic activity in the Upland. This thermal disturbance may be recorded by Rb–Sr age of 567 Ma for biotite from a HBU gneiss. Monazite age domains in metaquartzite (378 ± 28, 272 ± 44 Ma) suggest that low-grade metamorphism (260–320 °C, Mnz–Xno thermometry) of the CVS is not a result of Taconian orogenesis.     

Conditions and timing of high-pressure Variscan metamorphism in the South Carpathians, Romania

High-pressure (HP) metamorphic rocks, including garnet peridotite, eclogite, HP granulite, and HP amphibolite, are important constituents of several tectonostratigraphic units in the pre-Alpine nappe stack of the Getic–Supragetic (GS) basement in the South Carpathians. A Variscan age for HP metamorphism is firmly established by Sm–Nd mineral–whole-rock isochrons for garnet amphibolite, 358±10 Ma, two samples of eclogite, 341±8 and 344±7 Ma, and garnet peridotite, 316±4 Ma.

A prograde history for many HP metamorphic rocks is documented by the presence of lower pressure mineral inclusions and compositional zoning in garnet. Application of commonly accepted thermobarometers to eclogite (grt+cpx±ky±phn±pg±zo) yields a range in “peak” pressures and temperatures of 10.8–22.3 kbar and 545–745 °C, depending on tectonostratigraphic unit and locality. Zoisite equilibria indicate that activity of H₂O in some samples was substantially reduced, ca. 0.1–0.4. HP granulite (grt+cpx+hb+pl) and HP amphibolite (grt+hbl+pl) may have formed by retrogression of eclogites during high-temperature decompression. Two types of garnet peridotite have been recognized, one forming from spinel peridotite at ca. 1150–1300 °C, 25.8–29.0 kbar, and another from plagioclase peridotite at 560 °C, 16.1 kbar.

The Variscan evolution of the pre-Mesozoic basement in the South Carpathians is similar to that in other segments of the European Variscides, including widespread HP metamorphism, in which P–T–t characteristics are specific to individual tectonostratigraphic units, the presence of diverse types of garnet peridotite, diachronous subduction and accretion, nappe assembly in pre-Westphalian time due to collision of Laurussia, Gondwana, and amalgamated terranes, and finally, rapid exhumation, cooling, and deposition of eroded debris in Westphalian to Permian sedimentary basins.     

Lithosphere development in the Slave craton: a linked crustal and mantle perspective

The late tectonic evolution of the Slave craton involves extensive magmatism, deformation, and high temperature-low pressure (HT-LP) metamorphism. We argue that the nature of these tectonic events is difficult to reconcile with early, pre-2.7 Ga development and preservation of a thick tectosphere, and suggest that crust–mantle coupling and stabilization occurred only late in the orogenic development of the craton. The extent and repetitiveness of the tectonic reworking documented within the Mesoarchean basement complex of the western Slave, together with the development of large-volume, extensional mafic magmatism at 2.7 Ga within the basement complex argue against preservation of a widespread, thick, cool Mesoarchean tectosphere beneath the western Slave craton prior to Neoarchean tectonism. Broad-scale geological and geophysical features of the Slave craton, including orientation of an early F1 fold belt, distribution of ca. 2.63–2.62 Ga plutonic rocks, and the distribution of geochemical, petrological and geophysical domains within the mantle lithosphere collectively highlight the importance of an NE–SW structural grain to the craton. These trends are oblique to the earlier, ca. 2.7 Ga north–south trending boundary between Mesoarchean and Neoarchean crustal domains, and are interpreted to represent a younger structural feature imposed during northwest or southeast-vergent tectonism at ca. 2.64–2.61 Ga. Extensive plutonism, in part mantle-derived, crustal melting and associated HT-LP metamorphism argue for widespread mantle heat input to the crust, a feature most consistent with thin (<100 km) lithosphere at that time. We propose that the mantle lithosphere developed by tectonic imbrication of one or more slabs subducted beneath the craton at the time of development of the D1 structural grain, producing the early 2.63–2.62 Ga arc-like plutonic rocks. Subsequent collision (external to the present craton boundaries) possibly accompanied by partial delamination of some of the underthrust lithosphere, produced widespread deformation (D2) and granite plutonism throughout the province at 2.6–2.58 Ga. An implication of this model is that diamond formation in the Slave should be Neoarchean in age.     

The Homestead kimberlite, central Montana, USA: mineralogy, xenocrysts, and upper-mantle xenoliths

The Homestead kimberlite was emplaced in lower Cretaceous marine shale and siltstone in the Grassrange area of central Montana. The Grassrange area includes aillikite, alnoite, carbonatite, kimberlite, and monchiquite and is situated within the Archean Wyoming craton. The kimberlite contains 25–30 modal% olivine as xenocrysts and phenocrysts in a matrix of phlogopite, monticellite, diopside, serpentine, chlorite, hydrous Ca–Al–Na silicates, perovskite, and spinel. The rock is kimberlite based on mineralogy, the presence of atoll-textured groundmass spinels, and kimberlitic core-rim zoning of groundmass spinels and groundmass phlogopites.

Garnet xenocrysts are mainly Cr-pyropes, of which 2–12% are G10 compositions, crustal almandines are rare and eclogitic garnets are absent. Spinel xenocrysts have MgO and Cr₂O₃ contents ranging into the diamond inclusion field. Mg-ilmenite xenocrysts contain 7–11 wt.% MgO and 0.8–1.9 wt.% Cr₂O₃, with (Fe⁺³/Fe^tot) from 0.17–0.31. Olivine is the only obvious megacryst mineral present. One microdiamond was recovered from caustic fusion of a 45-kg sample.

Upper-mantle xenoliths up to 70 cm size are abundant and are some of the largest known garnet peridotite xenoliths in North America. The xenolith suite is dominated by dunites, and harzburgites containing garnet and/or spinel. Granulites are rare and eclogites are absent. Among 153 xenoliths, 7% are lherzolites, 61% are harzburgites, 31% are dunites, and 1% are orthopyroxenites. Three of 30 peridotite xenoliths that were analysed are low-Ca garnet–spinel harzburgites containing G10 garnets. Xenolith textures are mainly coarse granular, and only 5% are porphyroclastic.

Xenolith modal mineralogy and mineral compositions indicate ancient major-element depletion as observed in other Wyoming craton xenolith assemblages, followed by younger enrichment events evidenced by tectonized or undeformed veins of orthopyroxenite, clinopyroxenite, websterite, and the presence of phlogopite-bearing veins and disseminated phlogopite. Phlogopite-bearing veins may represent kimberlite-related addition and/or earlier K-metasomatism.

Xenolith thermobarometry using published two-pyroxene and Al-in-opx methods suggest that garnet–spinel peridotites are derived from 1180 to 1390 °C and 3.6 to 4.7 GPa, close to the diamond–graphite boundary and above a 38 mW/m² shield geotherm. Low-Ca garnet–spinel harzburgites with G10 garnets fall in about the same T and P range. Most spinel peridotites with assumed 2.0 GPa pressure are in the same T range, possibly indicating heating of the shallow mantle. Four of 79 Cr diopside xenocrysts have P–T estimates in the diamond stability field using published single-pyroxene P–T calculation methods.     

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.     

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.     

The Late Cretaceous diamondiferous pyroclastic kimberlites from the Fort à la Corne (FALC) field,Saskatchewan craton,Canada: Petrology,geochemistry and genesis

The Late Cretaceous (ca. 100 Ma) diamondiferous Fort à la Corne (FALC) kimberlite field in the Saskatchewan (Sask) craton, Canada, is one of the largest known kimberlite fields on Earth comprising essentially pyroclastic kimberlites. Despite its discovery more than two decades ago, petrological, geochemical and petrogenetic aspects of the kimberlites in this field are largely unknown. We present here the first detailed petrological and geochemical data combined with reconnaissance Nd isotope data on drill-hole samples of five major kimberlite bodies. Petrography of the studied samples reveals that they are loosely packed, clast-supported and variably sorted, and characterised by the presence of juvenile lapilli, crystals of olivine, xenocrystal garnet (peridotitic as well as eclogitic paragenesis) and Mg-ilmenite. Interclast material is made of serpentine, phlogopite, spinel, carbonate, perovskite and rutile. The mineral compositions, whole-rock geochemistry and Nd isotopic composition (Nd: + 0.62 to − 0.37) are indistinguishable from those known from archetypal hypabyssal kimberlites. Appreciably lower bulk-rock CaO (mostly < 5 wt%) and higher La/Sm ratios (12–15; resembling those of orangeites) are a characteristic feature of these rocks. Their geochemical composition excludes any effects of significant crustal and mantle contamination/assimilation. The fractionation trends displayed suggest a primary kimberlite melt composition indistinguishable from global estimates of primary kimberlite melt, and highlight the dominance of a kimberlite magma component in the pyroclastic variants. The lack of Nb-Ta-Ti anomalies precludes any significant role of subduction-related melts/fluids in the metasomatism of the FALC kimberlite mantle source region. Their incompatible trace elements (e.g., Nb/U) have OIB-type affinities whereas the Nd isotope composition indicates a near-chondritic to slightly depleted Nd isotope composition. The Neoproterozoic (~ 0.6–0.7 Ga) depleted mantle (T_DM) Nd model ages coincide with the emplacement age (ca. 673 Ma) of the Amon kimberlite sills (Baffin Island, Rae craton, Canada) and have been related to upwelling protokimberlite melts during the break-up of the Rodinia supercontinent and its separation from Laurentia (North American cratonic shield). REE inversion modelling for the FALC kimberlites as well as for the Jericho (ca. 173 Ma) and Snap Lake (ca. 537 Ma) kimberlites from the neighbouring Slave craton, Canada, indicate all of their source regions to have been extensively depleted (~ 24%) before being subjected to metasomatic enrichment (1.3–2.2%) and subsequent small-degree partial melting. These findings are similar to those previously obtained on Mesozoic kimberlites (Kaapvaal craton, southern Africa) and Mesoproterozoic kimberlites (Dharwar craton, southern India). The striking similarity in the genesis of kimberlites emplaced over broad geological time and across different supercontinents of Laurentia, Gondwanaland and Rodinia, highlights the dominant petrogenetic role of the sub-continental lithosphere. The emplacement of the FALC kimberlites can be explained both by the extensive subduction system in western North America that was established at ca. 150 Ma as well as by far-field effects of the opening of the North Atlantic ocean during the Late Cretaceous.     

U–Pb zircon and monazite geochronology of post-collisional Hercynian granitoids from the Central Iberian Zone (Northern Portugal)

In the Central Iberian Zone (CIZ) of the Iberian Massif large volumes of granitoids were emplaced during the post-collisional stage of the Hercynian orogeny (syn- to post-D3, the last ductile deformation phase). Twelve granitic units and a quartz monzodiorite were selected for a U–Pb zircon and monazite geochronological study. They represent successive stages of the D3 event. The Ucanha-Vilar, Lamego, Sameiro and Refoios do Lima plutons are coeval (313±2 Ma, 319±4 Ma, 316±2 Ma and 314±2 Ma, respectively) and belong to the earliest stage. Later on the Braga massif was emplaced, its different units yielding the same age: 309±3 Ma for the Braga granite, 309±1 Ma for the Gonça granite and 311±5 Ma for a related quartz monzodiorite. The Braga massif is subcontemporaneous with the Agrela and Celeirós plutons (307±3.5 Ma and 306±2 Ma, respectively), in agreement with field data. The Briteiros granite is younger (300±1 Ma), followed by the emplacement of the Peneda–Gerês massif (Gerês, Paufito, Illa and Carris granites). The Gerês granite, emplaced at 296±2 Ma, seems to represent a first magmatic pulse immediately followed by the intrusion of the Paufito granite at 290±2.5 Ma. For the Carris granite a minimum emplacement age of 280±5 Ma was obtained. Based on these results the following chronology is proposed: (1) syn-D3 biotite granitoids, 313–319 Ma; (2) late-D3 biotite-dominant granitoids, 306–311 Ma; (3) late- to post-D3 granitoids, ca. 300 Ma; (4) post-D3 granitoids, 290–296 Ma. These chronological data indicate that successive granitic intrusions were emplaced in the CIZ during a short time span of about 30 Ma that corresponds to the latest stages of the Hercynian orogeny. A rapid and drastic change occurred at about 300 Ma, between a compressive ductile tectonic regime (D3, ca. 300–320 Ma) associated to calc-alkaline, monzonitic and aluminopotassic plutonism and a fragile phase of deformation (D4) which controlled the emplacement of the subalkaline ferro-potassic plutonism at 290–296 Ma.     

A subsurface criterion for predictive exploration of kimberlites from Bouguer gravity in the eastern Dharwar craton,India

The Archaean-Proterozoic Dharwar craton has many recorded occurrences of diamondiferous kimberlites. Reports of kimberlite emplacement in parts of the tectonically complex eastern Dharwar craton and a significant density contrast between kimberlites and the host peninsular gneisses motivated us to conduct gravity studies in the Narayanpet-Irladinne area of the eastern Dharwar craton. This region is contiguous with the Maddur-Narayanpet kimberlite that lies to its north, while the river Krishna lies to its south. From observed association of reported kimberlites in the Maddur-Narayanpet field with subsurface topography of the assumed three-layer earth section obtained by Bouguer gravity modelling, we developed a subsurface criterion for occurrence of kimberlites in the present study area. Using this criterion, five potential zones for kimberlite localization were identified in the Narayanpet-Irladinne region, eastern Dharwar craton.     

Punctuated,long-lived emplacement history of the Renard 2 kimberlite,Canada, revealed by new high precision U-Pb groundmass perovskite dating

Kimberlites are rare volatile-rich ultramafic magmas thought to erupt in short periods of time (<1 Myr) but there is a growing body of evidence that the emplacement history of a kimberlite can be significantly more protracted. In this study we report a detailed geochronology investigation of a single kimberlite pipe from the Renard cluster in north-central Québec. Ten new high precision ID-TIMS (isotope dilution – thermal ionization mass spectrometry) U-Pb groundmass perovskite dates from the main pipe-infilling kimberlites and several small hypabyssal kimberlites from the Renard 2 pipe indicate kimberlite magmatism lasted at least ~20 Myr. Two samples of the main pipe-infilling kimberlites yield identical weighted mean ²⁰⁶Pb/²³⁸U perovskite dates with a composite date of 643.8 ± 1.0 Myr, interpreted to be the best estimate for main pipe emplacement. In contrast, six hypabyssal kimberlite samples yielded a range of weighted mean ²⁰⁶Pb/²³⁸U perovskite dates between ~652-632 Myr. Multiple dates determined from these early-, syn- and late-stage small hypabyssal kimberlites in the Renard 2 pipe demonstrate this rock type (commonly used to date kimberlites) help to constrain the duration of kimberlite intrusion history within a pipe but do not necessarily reliably record the emplacement age of the main diatreme in the Renard cluster. Our results provide the first robust geochronological data on a single kimberlite that confirms the field relationships initially observed by Wagner (1914) and Clement (1982); the presence of antecedent (diatreme precursor) intrusions, contemporaneous (syn-diatreme) intrusions, and consequent (post-diatreme) cross-cutting intrusions. The results of this detailed U-Pb geochronology study indicate a single kimberlite pipe can record millions of years of magmatism, much longer than previously thought from the classical viewpoint of a rapid and short-duration emplacement history.

     

1214 ± 5 Ma dyke from the Darling Range,southwestern Yilgarn Craton,Western Australia

High thorium euhedral, twinned and elongate zircons from the felsic part of a mafic dyke located in the Archaean Yilgarn Craton approximately 30 km northeast of Perth and approximately 2 km east of the Darling Fault, have consistent 207 Pb/ 206 Pb ages of 1214 ± 5 Ma. This age is interpreted as the age of dyke emplacement and is identical, within the uncertainties, with other U–Pb dyke ages reported for the southwest Yilgarn Craton. The present result extends the known occurrence of ca 1210 Ma dykes to the western margin of the Yilgarn Craton and confirms earlier conclusions that a major mafic dyke emplacement occurred throughout the southern Yilgarn Craton during a short‐lived magmatic pulse at ca 1210 Ma.     

Diamonds and their mineral inclusions from the A154 South pipe, Diavik Diamond Mine, Northwest territories, Canada

Mineral inclusions recovered from 100 diamonds from the A154 South kimberlite (Diavik Diamond Mines, Central Slave Craton, Canada) indicate largely peridotitic diamond sources (83%), with a minor (12%) eclogitic component. Inclusions of ferropericlase (4%) and diamond in diamond (1%) represent “undetermined” parageneses.

Compared to inclusions in diamonds from the Kaapvaal Craton, overall higher CaO contents (2.6 to 6.0 wt.%) of harzburgitic garnets and lower Mg-numbers (90.6 to 93.6) of olivines indicate diamond formation in a chemically less depleted environment. Peridotitic diamonds at A154 South formed in an exceptionally Zn-rich environment, with olivine inclusions containing more than twice the value (of  52 ppm) established for normal mantle olivine. Harzburgitic garnet inclusions generally have sinusoidal rare earth element (REE_N) patterns, enriched in LREE and depleted in HREE. A single analyzed lherzolitic garnet is re-enriched in middle to heavy REE resulting in a “normal” REE_N pattern. Two of the harzburgitic garnets have “transitional” REE_N patterns, broadly similar to that of the lherzolitic garnet. Eclogitic garnet inclusions have normal REE_N patterns similar to eclogitic garnets worldwide but at lower REE concentrations.

Carbon isotopic values (δ¹³C) range from − 10.5‰ to + 0.7‰, with 94% of diamonds falling between − 6.3‰ and − 4.0‰. Nitrogen concentrations range from below detection (< 10 ppm) to 3800 ppm and aggregation states cover the entire spectrum from poorly aggregated (Type IaA) to fully aggregated (Type IaB). Diamonds without evidence of previous plastic deformation (which may have accelerated nitrogen aggregation) typically have < 25% of their nitrogen in the fully aggregated B-centres. Assuming diamond formation beneath the Central Slave to have occurred in the Archean [Westerlund, K.J., Shirey, S.B., Richardson, S.H., Gurney, J.J., Harris, J.W., 2003b. Re–Os systematics of diamond inclusion sulfides from the Panda kimberlite, Slave craton. VIIIth International Kimberlite Conference, Victoria, Canada, Extended Abstracts, 5p.], such low aggregation states indicate mantle residence at fairly low temperatures (< 1100 °C). Geothermometry based on non-touching inclusion pairs, however, indicates diamond formation at temperatures around 1200 °C. To reconcile inclusion and nitrogen based temperature estimates, cooling by about 100–200 °C shortly after diamond formation is required.     

Formation of continental flood volcanism — The perspective of setting of melting

Models of continental flood basalt (CFB) formation are evaluated by examining their implications for the setting, mainly temperature and depth, of melting which is assumed to result from adiabatic decompression. Most attractive is the model of melting in upwelling bodies (probably plume heads) rising to the base of the continental lithosphere. This constrains the melting to 120–150 km or deeper (continental lithospheric thickness) and thus the plume potential temperatures to ≥ 300 °C higher than ambient mantle. The primary melts should be hot, MgO-rich, modified during ascent, and assimilate components of the lithosphere, which can provide the continental-like geochemical signature of many CFB. Circulation within the upwellings and presence of eclogite patches also influence magma generation and composition. Dehydration melting when plumes heat the lowermost lithosphere can generate CFB only if the source region contains ca. 15% hydrous minerals beneath the entire area covered by flood volcanics, which is difficult to justify. On the other hand, assimilation of “continental” chemical components from large parts of the lithosphere does not require such extreme metasomatism. Decompression melting under strongly thinned rifted lithosphere requires lower potential temperatures of the rising material and lesser modification of the primary magmas than the plume head model of CFB formation. Available observations do not support the contemporaneity of flood volcanism with rifts having the required sizes and histories, but more information is needed to further test this model. On the other hand, magma production can assist rift initiation and lithospheric rupture, so subsequent thinning can explain the common formation of volcanic rifted margins immediately following CFB emplacement. Ancient LIP should record the same processes as seen in young CFB.     

Petrography,Sr-isotope geochemistry and geochronology of the Nxau Nxau kimberlites,north-west Botswana

The Nxau Nxau kimberlites in northwest Botswana belong to the Xaudum kimberlite province that also includes the Sikereti, Kaudom and Gura kimberlite clusters in north-east Namibia. The Nxau Nxau kimberlites lie on the southernmost extension of the Congo Craton, which incorporates part of the Damara Orogenic Belt on its margin. The Xaudum kimberlite province is geographically isolated from other known clusters but occurs within the limits of the NW-SE oriented, Karoo-aged Okavango Dyke Swarm and near NE-SW faults interpreted as the early stages of the East African Rift System. Petrographic, geochronological and isotopic studies were undertaken to characterise the nature of these kimberlites and the timing of their emplacement. The Nxau Nxau kimberlites exhibit groundmass textures, mineral phases and Sr-isotope compositions (⁸⁷Sr/⁸⁶Sr_i of 0.7036 ± 0.0002; 2σ) that are characteristic of archetypal (Group I) kimberlites. U-Pb perovskite, ⁴⁰Ar/³⁹Ar phlogopite and Rb-Sr phlogopite ages indicate that the kimberlites were emplaced in the Cretaceous, with perovskite from four samples yielding a preferred weighted average U-Pb age of 84 ± 4 Ma (2σ). This age is typical of many kimberlites in southern Africa, indicating that the Xaudum occurrences form part of this widespread Late Cretaceous kimberlite magmatic province. This time marks a significant period of tectonic stress reorganisation that could have provided the trigger for kimberlite magmatism. In this regard, the Nxau Nxau kimberlites may form part of a NE-SW oriented trend such as the Lucapa corridor, with implications for further undiscovered kimberlites along this corridor.