Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-07T21:22:44.084Z Has data issue: false hasContentIssue false

The Pre-Grenvillian assembly of the southeastern Laurentian margin through the U–Pb–Hf detrital zircon record of Mesoproterozoic supracrustal sequences (Central Grenville Province, Quebec, Canada)

Published online by Cambridge University Press:  15 November 2021

K. Papapavlou*
Affiliation:
GEOTOP, Université du Québec à Montréal, Montréal, QC, H3C 3P8, Canada
A. Moukhsil
Affiliation:
Ministère de l’énergie et des ressources naturelles, Direction du bureau de la connaissance géoscientifique du Québec, Québec, QC, G1H 6R1, Canada
A. Poirier
Affiliation:
GEOTOP, Université du Québec à Montréal, Montréal, QC, H3C 3P8, Canada
J.H.F.L. Davies
Affiliation:
GEOTOP, Université du Québec à Montréal, Montréal, QC, H3C 3P8, Canada
*
Author for correspondence: Konstantinos Papapavlou, Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The detrital zircon perspective on the pre-collisional crustal evolution of the Grenville Province remains poorly explored. In this study, we conducted in situ laser ablation U–Pb–Hf isotopic microanalysis on detrital zircon grains from three pre-orogenic (>1 Ga) supracrustal sequences that crop out in the Central Grenville Province (Lac Saint-Jean region, QC, CA). Detrital zircon grains from vestiges of these sequences record three dominant age peaks at c. 1.46 Ga, 1.62 Ga, 1.85 Ga, and a subordinate peak at 2.7 Ga. The 1.46 Ga and 1.62 Ga age peaks are recorded in detrital zircon grains from a quartzite associated with a metavolcanic sequence (i.e. Montauban Group) with a maximum depositional age of c. 1.44 Ga. In contrast, the c.1.85 Ga age peak is observed from recycled zircon grains in metasediments with maximum depositional ages between 1.2 and 1.3 Ga. The suprachondritic Hf isotope composition in detrital zircon grains of the 1.46 Ga and 1.62 Ga age populations records juvenile crustal growth during peri-Laurentian accretionary orogenesis related to the Pinwarian (1.4–1.5 Ga) and Mazatzalian–Labradorian (1.6–1.7 Ga) events. The detrital zircon grains associated with Penokean–Makkovikian (1.8–1.9 Ga) source rocks record reworking of c. 2.7 Ga continental crust derived from a near-chondritic mantle reservoir. Overall, crust-forming and basement reworking events associated with accretionary orogenesis in southeastern Laurentia are retained in the detrital zircon load of Precambrian basins even after the terminal Grenvillian collision and assembly of Rodinia.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re- use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press

1. Introduction

The Grenville orogen (c. 1090–980 Ma; Rivers, Reference Rivers2008, Reference Rivers2021) is a global-scale collisional orogenic system that terminated the long-lived accretionary growth of Laurentia, leading to the formation of the supercontinent Rodinia (Hoffman, Reference Hoffman1988; Whitmeyer and Karlstrom, Reference Whitmeyer and Karlstrom2007). In the southeastern Canadian shield (Grenville Province; Fig. 1a), c. 1.15 Ga (Shawinigan) and 1.05 Ga (Ottawan) high-grade metamorphism and penetrative deformation concealed the traits of accretionary tectonism, leaving behind gneiss terranes of ambiguous origin. Regional-scale Sm–Nd isotope mapping along the Grenville Province has attempted to untangle the pre-collisional crustal formation history of these large tracts of gneisses (Dickin & McNutt, Reference Dickin and McNutt1989; Dickin, Reference Dickin2000; Dickin et al., Reference Dickin, McNutt, Martin and Guo2009). However, the perspective of detrital zircon on the pre-collisional crustal growth of the Grenville Province remains underexplored despite zircon grains of Laurentian affinity being documented in pre-, syn- and post-collisional metasedimentary sequences (Sager-Kinsman and Parrish, Reference Sager-Kinsman and Parrish1993; Van Breemen and Corriveau, Reference Van Breemen and Corriveau2005; David et al. Reference David, Moukhsil and Dion2010; Kamo et al. Reference Kamo, Heaman and Gower2011; Lasalle et al. Reference Lasalle, Fisher, Indares and Dunning2013; Spencer et al. Reference Spencer, Cawood, Hawkesworth, Prave, Roberts, Horstwood and Whitehouse2015; Groulier et al. Reference Groulier, Indares, Dunning, Moukhsil and Wälle2018; Peck et al. Reference Peck, Quinan and Selleck2019). Together, these studies have documented, through detrital zircon U–Pb geochronology, the presence of 1.4–1.9 Ga sources but, further insights about their crustal production and evolution are missing.

Fig. 1. Simplified geological map of the Lac-Saint-Jean region. The inset (a) depicts the crustal terranes of the Grenville Province sensu Dickin, Reference Dickin2000). The white dashed rectangle depicts the study area. The red circles in (b) denote the four sampling sites.

In this regard, the interrogation of the detrital zircon record, from the Hadean to Phanerozoic, using U–Pb–Hf ± O isotope systematics has placed critical new constraints on the secular evolution and growth of the continental crust (Stevenson & Patchett, Reference Stevenson and Patchett1990; Kemp et al. Reference Kemp, Hawkesworth, Paterson and Kinny2006; Belousova et al. Reference Belousova, Kostitsyn, Griffin, Begg, O’Reilly and Pearson2010; Hawkesworth et al. Reference Hawkesworth, Dhuime, Pietranik, Cawood, Kemp and Storey2010; Parman, Reference Parman2015). Crucially, the initial Hf isotopic composition of single zircon grains monitors crust–mantle differentiation events while remaining resistant to high-grade metamorphism and/or hydrothermal fluid flow that can compromise the whole-rock Sm–Nd isotope system (Moorbath et al. Reference Moorbath, Whitehouse and Kamber1997; Vervoort & Blichert-Toft, Reference Vervoort and Blichert-Toft1999; Hoffmann et al., Reference Hoffmann, Münker, Polat, Rosing and Schulz2011). The advantage of the U–Pb and Lu–Hf isotopic couplet has been demonstrated in many crustal evolution studies (e.g. Condie et al., Reference Condie, Beyer, Belousova, Griffin and O’Reilly2005; Kemp et al. Reference Kemp, Hickman, Kirkland and Vervoort2015; Fisher & Vervoort, Reference Fisher and Vervoort2018) but remains underexploited in the metasedimentary sequences of the Grenville Province. The scarcity of zircon U–Pb–Hf isotopic data from these sequences has inhibited our ability to ground-truth crustal evolution models, which are heavily based on whole-rock Sm–Nd isotopic data from felsic orthogneisses, and in some cases interpreted without U–Pb age constraints from the same samples (e.g. Dickin, Reference Dickin2000).

In this contribution we combine, for the first time, detrital zircon U–Pb (n = 623) and Lu–Hf (n = 366) isotopic data from rare outcrops of quartzites in three Mesoproterozoic supracrustal sequences, from the Lac Saint-Jean region in the central part of the Grenville Province (Fig. 1b). The main aims of our study are to: (a) provide new constraints on the pre-collisional geodynamic evolution and crustal growth of the Grenville Province through examination of the U–Pb–Hf detrital zircon record, (b) compare the zircon U–Pb–Hf isotope signatures with the existing whole-rock Nd isotopic record of crustal growth in the central Grenville Province and (c) better constrain the depositional ages and sources of these poorly documented supracrustal sequences.

2. Pre-collisional accretionary terranes of the Grenville Province

The Grenville Province is divided into two major lithotectonic belts, the Parautochthonous and the Allochthonous Belts (Rivers, Reference Rivers1997). The Parautochthonous consists of basement rocks with Laurentian isotopic and age characteristics, whereas the Allochthonous Belt comprises translated nappes of amphibolite to granulite-facies gneisses (Rivers, Reference Rivers2008). Biotite gneisses in the western part of the Parautochthonous belt exhibit depleted-mantle Nd model ages from 2.5 to 2.9 Ga (Guo & Dickin, Reference Guo and Dickin1996) although Archaean zircon U–Pb crystallization ages have rarely been reported from these areas (Krogh, Reference Krogh1994; Jordan et al. Reference Jordan, Indares and Dunning2006).

The Allochthonous Belt has been characterized, through Nd isotope mapping of granitoid orthogneisses, as an amalgamation of Palaeo- to Mesoproterozoic accretionary terranes (Dickin, Reference Dickin2000). These terranes, sensu Dickin, (Reference Dickin2000), are the Makkovikia–Barilia, Labradoria and Quebecia terranes (Fig. 1a) with crustal formation ages (i.e. depleted-mantle model ages) of c. 1.5 Ga, 1.7–1.8 Ga and 1.9 Ga respectively. In more detail, the Makkovikia has been characterized as a terrane of reworked crust located at the easternmost part of the Grenville Province (i.e. Labrador), and the Barilia as a 1.8–2.0 Ga juvenile Palaeoproterozoic arc in Ontario (Dickin et al., Reference Dickin, McNutt, Martin and Guo2009). These terranes are regarded as products of an arc–continent collisional tectonothermal event known as the Penokean–Makkovikian orogeny, which is well documented in the Great Lakes region and took place between c. 1.8 and 1.9 Ga (Van Schmus & Bickford, Reference Van Schmus and Bickford1981; Barovich et al. Reference Barovich, Patchett, Peterman and Sims1989). U–Pb age-peaks in detrital zircon spectra at c. 1.8–1.9 Ga, from different metasedimentary sequences, have been reported along the Grenville Province and are commonly ascribed to Penokean–Makkovikian sources (Kamo et al. Reference Kamo, Heaman and Gower2011; Moukhsil et al. Reference Moukhsil, Solgadi, Clark, Blouin, Indares and Davis2013; Culshaw et al. Reference Culshaw, Foster, Marsh, Slagstad and Gerbi2016; Peck et al. Reference Peck, Quinan and Selleck2019).

Labradoria was formed during the Mazatzalian–Labradorian orogeny between 1.6 and 1.7 Ga (Wasteneys et al. Reference Wasteneys, Kamo, Moser, Krogh, Gower and Owen1997; Papapavlou et al. Reference Papapavlou, Darling, Storey, Lightfoot, Moser and Lasalle2017) and crops out mostly in Labrador, in the eastern part of the Grenville Province (Connelly & Heaman, Reference Connelly and Heaman1993; Gower, Reference Gower2019). Labradorian-age metaplutonic and mafic/ultramafic suites have also been documented in the central and western parts of the province (Dunning & Indares, Reference Dunning and Indares2010, Culshaw et al. Reference Culshaw, Foster, Marsh, Slagstad and Gerbi2016). In eastern Labrador, c. 1.7 Ga tonalitic and granitic rocks exhibit near-chondritic isotope signatures (ϵNd = −0.13 to +0.42; Scharer, Reference Schärer1991). Later studies suggested that the Labradorian crust-forming event was mostly juvenile but also involved magmatic reworking (i.e. remelting) of older Makkovikian crust in the eastern part of the province (Dickin, Reference Dickin2000). As with detrital zircon grains of Penokean age, grains of Labradorian age are ubiquitous in the pre- and syn-Grenvillian metasedimentary record (Corrigan & Van Breemen, Reference Corrigan and Van Breemen1997; Lasalle et al. Reference Lasalle, Fisher, Indares and Dunning2013; Culshaw et al. Reference Culshaw, Foster, Marsh, Slagstad and Gerbi2016; Peck et al. Reference Peck, Quinan and Selleck2019).

The youngest accretionary event in the pre-collisional history of the Grenville Province is the c. 1.46–1.52 Ga Pinwarian orogeny (Tucker & Gower, Reference Tucker and Gower1994; Gower & Krogh, Reference Gower and Krogh2002) that is expressed through titanite U–Pb dating in reactivated mylonitic shear zones in adjacent orogens (e.g. Southern Province; Papapavlou et al. Reference Papapavlou, Darling, Lightfoot, Lasalle, Gibson, Storey and Moser2018). In the eastern part of the Grenville Province, Pinwarian magmatism reworked Labradoria whereas further west it led to the construction of Quebecia (Dickin, Reference Dickin2000). Recent studies focusing on the Pinwarian magmatism suggest that it took place above two subduction zones with opposing polarities: a north-dipping subduction zone below Laurentia and a south-dipping intra-oceanic subduction zone, which led to the construction of a nascent magmatic arc (Martin & Dickin, Reference Martin and Dickin2005; Vautour & Dickin, Reference Vautour and Dickin2019).

The Mid- to Late-Proterozoic Elzevirian (c. 1.2–1.3 Ga) and Shawinigan (c. 1.19–1.14 Ga) orogenic pulses constitute the final accretionary events (Rivers, Reference Rivers2008) before the collisional Grenville orogeny that lasted until c. 960 Ma (Jannin et al. Reference Jannin, Gervais, Moukhsil, Augland and Crowley2018, Turlin et al. Reference Turlin, Vanderhaeghe, Gervais, André-Mayer, Moukhsil and Zeh2019) and was ended by late collisional alkaline magmatism (Groulier et al. Reference Groulier, Turlin, André-Mayer, Ohnenstetter, Crépon, Boulvais, Poujol, Rollion-Bard, Zeh, Moukhsil, Solgadi and El basbas2020b).

3. Regional geology of the Saguenay – Lac Saint-Jean region

The Lac Saint-Jean region is located in the central part of the Grenville Province within the Allochthonous Belt (Fig. 1). The oldest lithostratigraphic unit in the area is the Pinwarian Cap de la Mer amphibolite, with a U–Pb zircon age of 1506 ± 13 Ma, which is included within a group of basement units called the Saguenay Gneiss Complex (Hébert & Van Breemen, Reference Hébert and Van Breemen2004). This area belongs to the Portneuf – St Maurice lithotectonic domain that consists of two main units, the La Bostonnais complex, which is a 1.4–1.37 Ga suite of metamorphosed calc-alkaline plutonic rocks, and the Montauban group, a c. 1.45 Ga metamorphosed volcano-sedimentary sequence commonly observed as enclaves within the La Bostonnais complex (Fig. 1b; L Nadeau et al., unpub. data, 1992; Corrigan & Van Breemen, Reference Corrigan and Van Breemen1997). The area is dominated by anorthosite–mangerite–charnockite–granite (AMCG) magmatic suites and hosts the >20 000 km2 Lac Saint-Jean Anorthositic Suite (Woussen et al. Reference Woussen, Dimroth, Corriveau and Archer1981). U–Pb zircon dating (isotope dilution thermal ionization mass spectrometry (ID-TIMS)) has shown that the Lac Saint-Jean anorthosite suite was emplaced between c. 1160 and 1140 Ma and was followed by two less voluminous AMCG magmatic pulses at c. 1082–1050 Ma and 1020–1010 Ma (Higgins & Van Breemen, Reference Higgins and Van Breemen1996). Regional Grenvillian metamorphism in the area attained granulite facies with PT conditions of 6.2 kbar and 780 °C followed by a retrograde amphibolite facies event reaching 4.7 kbar and 690 °C , respectively (Hébert & Van Breemen, Reference Hébert and Van Breemen2004).

3.a. Supracrustal lithological complexes

In this study we explore the detrital zircon record from three of the four supracrustal complexes that crop out in the wider area of the Lac Saint-Jean region. These are the Montauban, Wabash, Barrois and Saint-Onge complexes (Fig. 1b). The Montauban complex consists of quartzofeldspathic gneisses interlayered with amphibolite, pillowed tholeiitic metabasalts with boudinage structure, volcanic and volcaniclastic rocks (i.e. lappili tuff), metalliferous calc-silicate rocks, garnetite, and minor quartzite which was sampled in this study (sample 17GC1086; Stamatelopoulou-Seymour & Maclean, Reference Stamatelopoulou-Seymour and MacLean1977; Rondot, Reference Rondot1978). Supracrustal lithologies of the Montauban group are of economic importance since they host metamorphosed exhalative massive sulphide (Zn–Pb–Cu–Au–Ag) and stratiform gold occurrences (Bernier et al. Reference Bernier, Pouliot and MacLean1987). The Wabash group is a metamorphosed volcano-sedimentary sequence that consists of four tectonostratigraphic units that mainly comprise migmatitic paragneisses, calcitic and dolomitic marbles, minor impure quartzite layers (sample 18YD2075), amphibolites and quartzofeldspathic rocks. The Barrois group comprises chiefly Opx-bearing orthogneisses, Bt-bearing paragneisses, charnockitic to enderbitic gneisses, amphibolite and garnetite. The Barrois group is divided into the Barrois I unit that contains metasedimentary lithologies (i.e. paragneisses) and the Barrois II unit that consists of orthogneisses which likely constitute the basement onto which the Barrois I unit was deposited. The biotite-bearing paragneisses consists of interbedded, decimetre- to metre-scale, quartzite horizons, which were sampled here (samples 19AM02 and 19FS7085). Aeromagnetic maps suggest the juxtaposition of the Barrois and Wabash supracrustal packages most likely along a stratigraphic contact. The Saint-Onge group is composed of quartzofeldspathic paragneisses, calc-silicate rocks, marbles and minor quartzite. Lithological members of this group show extensive alteration in outcrop scale and were not sampled in this study. Four quartzite samples from the Lac Saint-Jean region were collected: these samples are the 17GC1086 (Montauban Group), 18YD2075 (Wabash Complex), 19FS7085 and 19AM02 (Barrois Complex). Sample coordinates and petrographic descriptions are provided in supplementary material file 1 (in the Supplementary Material available online at https://doi.org/10.1017/S0016756821001023).

4. Analytical techniques and methodology

Sample processing and microanalytical work were conducted at the facilities of GEOTOP-UQAM (Montreal, QC, CA) using conventional zircon separation procedures involving crushing with a jaw crusher and disc mill followed by Wilfley table, heavy liquid, and magnetic separation using a Franz isodynamic separator. The zircon grains were annealed at 1000 °C for 48 hours and were then mounted in epoxy for imaging and microanalysis. Cathodoluminescence (CL) imaging of the zircon grains was conducted using a Hitachi S-3400N scanning electron microscope (SEM). The U–Pb and Lu–Hf isotopic analyses were performed during separate sessions on overlapping or adjacent analytical spots within the same CL growth domain assuming no differences in the isotopic composition or age of the different volumes. The U–Pb and Lu–Hf analytical sessions were performed using a Photon Machines Analyte G2 excimer laser (193 nm) coupled with either a single-collector high-resolution inductively coupled plasma mass spectrometer (Nu Attom, HR-ICP-MS, for U–Pb microanalysis) or a multicollector ICP-MS (Nu Plasma II, for Hf isotope microanalysis). The ablation took place in a two-volume HELEX II sample cell with He as the carrier gas which was mixed with the Ar make-up gas downstream before entering the ICP torch. Instrument parameters of the analytical sessions and secondary reference material data are provided in supplementary material file 1 (in the supplementary material available online at https://doi.org/10.1017/S0016756821001023). The U–Pb and Lu–Hf time-resolved isotopic data were reduced and processed using iolite 3.63 and IsoplotR (Paton et al. Reference Paton, Hellstrom, Paul, Woodhead and Hergt2011; Vermeesch, Reference Vermeesch2018).

4.a. U–Pb isotopic analysis

The U–Pb isotopic data were collected using a spot size of 25 μm at a repetition rate of 5 Hz and fluence of 3 J cm−2. The 91500 zircon (206Pb/238U age of 1062.4 ± 0.4 Ma; Wiedenbeck et al., Reference Wiedenbeck, Allé, Corfu, Griffin, Meier, Oberli, Von Quadt, Roddick and Spiegel1995) was used as the primary reference material for the correction of laser-induced elemental fractionation and instrument drift. Secondary zircon reference materials were interspersed every 10–15 unknowns and yielded the following weighted mean 207Pb/206Pb ages: BB9 (567 ± 16 Ma, 2s, MSWD = 2, n = 8), Mudtank (713 ± 17 Ma, 2s, MSWD = 1.5, n = 35), GJ-1 (604 ± 11 Ma, 2s, MSWD = 2.7, n = 38), UQZ8 (1119 ± 9 Ma, 2s, MSWD = 5, n = 15) and Pentecote (1367 ± 39 Ma, 2s, MSWD = 3, n = 10). These dates corroborate within uncertainty the published reference values (Martignole et al. Reference Martignole, Machado and Nantel1993; Machado & Gauthier, Reference Machado and Gauthier1996; Jackson et al. Reference Jackson, Pearson, Griffin and Belousova2004; Santos et al. Reference Santos, Lana, Scholz, Buick, Schmitz, Kamo, Gerdes, Corfu, Tapster, Lancaster, Storey, Basei, Tohver, Alkmin, Nalini, Krambrock, Fantini and Wiedenbeck2017; Gain et al. Reference Gain, Gréau, Henry, Belousova, Dainis, Griffin and O’Reilly2019).

The U–Pb isotopic data (see also U–Pb supplementary material available online at https://doi.org/10.1017/S0016756821001023) that are presented from the cores of zircon grains have ≤10% discordance and Th/U > 0.1 following the rationale of previous studies (Papapavlou et al. Reference Papapavlou, Strachan, Storey and Bullen2021). The filtered data are presented as 207Pb/206Pb ages in Kernel Density Estimation (KDE) diagrams with an adaptive bandwidth using the Java application DensityPlotter (Vermeesch, Reference Vermeesch2012). Maximum depositional ages (MDAs) were calculated using the maximum likelihood age (MLA) estimation algorithm of Vermeesch (Reference Vermeesch2021). Grains that do not show igneous zoning in CL were not considered during the calculation of MDAs but are presented in the U–Pb supplementary material (available online at https://doi.org/10.1017/S0016756821001023).

4.b. Lu–Hf isotopic analysis

The Lu–Hf isotopic data were collected using a 50 μm spot with a laser frequency of 15 Hz and fluence of c. 9 Jcm−2. A typical Lu–Hf analysis consisted of 25 s of background measurement and 20 s of data acquisition. The 171Yb, 173Yb, 174Yb, 175Lu, 176(Yb + Lu + Hf), 177Hf, 178Hf, 179Hf, 180Hf, 181Ta and 182W beams were collected using 11 Faraday cups in fixed positions attached to 1011 Ω feedback resistors, with the mass separation being controlled by an optical zoom lens between the magnetic sector and the detectors. The 91500 zircon was used as the primary reference material to correct instrument drift during the analytical session. Mass bias and isobaric interference corrections of 176Yb and 176Lu on 176Hf were made following the rationale and guidelines of previous studies (Woodhead et al., Reference Woodhead, Hergt, Shelley, Eggins and Kemp2004; Fisher et al. Reference Fisher, Vervoort and Hanchar2014). Specifically, the β Yb (Yb mass bias) was determined by measuring the 173Yb/171Yb ratio relative to 173Yb/171Yb = 1.132685 (Chu et al., Reference Chu, Taylor, Chavagnac, Nesbitt, Boella, Milton, German, Bayon and Burton2002). Following the calculation of β Yb, the contribution of 176Yb on 176Hf was subtracted using the intensity of the interference-free 173Yb and a 176Yb/173Yb = 0.79618 (Chu et al., Reference Chu, Taylor, Chavagnac, Nesbitt, Boella, Milton, German, Bayon and Burton2002). The contribution of the 176Lu on the 176Hf signal was estimated assuming that β Lu = β Yb using the intensity of the interference-free 175Lu with a 176Lu/175Lu = 0.02656 (Blichert-Toft et al. Reference Blichert-Toft, Chauvel and Albarède1997). The mass bias correction on the 176Hf/177Hf ratios was estimated using 179Hf/177Hf = 0.7325 (Patchett & Tatsumoto, Reference Patchett and Tatsumoto1981). Mass bias corrections for Yb, Hf and Lu were made using an exponential law. For the Yb mass bias corrections of grains with low Yb/Hf ratios (e.g. Mudtank and Munzirc02a, which have 173Yb beam intensities below 5mV), the β Hf was used instead of the poorly characterized β Yb. For quality control on the mass bias corrections, the 178Hf/177Hf ratio on the analysed reference materials yielded a value of 1.46731 ± 52 (1s), in agreement with the published value of 1.46735 (Thirlwall & Anczkwiewicz, Reference Thirlwall and Anczkiewicz2004).

Secondary zircon reference materials with variable 176Yb/177Hf ratios were analysed during the Lu–Hf session to ensure that our 176Yb correction yielded accurate results (see also Lu–Hf supplementary material file, available online at https://doi.org/10.1017/S0016756821001023). The following weighted mean 176Hf/177Hf values were obtained: BB9 zircon (0.281666 ± 26, MSWD = 1.5, n = 9, 2s, ϵHf value of −27.2 ± 1.2), Maniitsoq (0.280865 ± 15; MSWD = 0.9, n = 10, 2s, ϵHf value of −0.28 ± 0.9), Munzirc12c (0.282129 ± 12; MSWD = 1.3, n = 17, 2s), Munzirc02a (0.282160 ± 8, MSWD = 1.5, n = 18, 2s), Munzirc32c (0.282150 ± 19; MSWD = 2, n = 14, 2s), Munzirc42b (0.282164 ± 14; MSWD = 0.3, n = 11, 2s), Plešovice (0.282475 ± 6; MSWD = 0.6, n = 30, 2s, ϵHf value of −3.5 ± 0.6), SA01 (0.282296 ± 24; MSWD = 1.5, n = 11, 2s, ϵHf value of −5.7 ± 1.5), GJ-1 (0.282006 ± 16; MSWD = 0.6, n = 17, 2s, ϵHf value of −14.3 ± 1.1) and Mudtank (0.282528 ± 11; MSWD = 1.5, n = 16, 2s, ϵHf value of +7.1 ± 1) (Morel et al. Reference Morel, Nebel, Nebel-Jacobsen, Miller and Vroon2008; Sláma et al., Reference Sláma, Košler, Condon, Crowley, Gerdes, Hanchar, Horstwood, Morris, Nasdala, Norberg, Schaltegger, Schoene, Tubrett and Whitehouse2008; Fisher et al. Reference Fisher, Hanchar, Samson, Dhuime, Blichert-Toft, Vervoort and Lam2011; Santos et al. Reference Santos, Lana, Scholz, Buick, Schmitz, Kamo, Gerdes, Corfu, Tapster, Lancaster, Storey, Basei, Tohver, Alkmin, Nalini, Krambrock, Fantini and Wiedenbeck2017; Gain et al. Reference Gain, Gréau, Henry, Belousova, Dainis, Griffin and O’Reilly2019; Marsh et al. Reference Marsh, Jorgensen, Petrus, Hamilton and Mole2019; Huang et al., Reference Huang, Wang, Yang, Ramezani, Yang, Zhang, Yang, Xia, Feng, Lin, Wang, Ma, He, Xie and Wu2020). The reproducibility (2s) of the secondary reference materials varies from 0.004 to 0.029 %. The analytical errors are quadratic additions of the in-run error (2SE) and the reproducibility (2s) of the Plešovice standard and are typically ∼± 2 ϵHf units. The initial 176Hf/177Hf isotope ratios and ϵHf values were calculated using the 176Lu decay constant of 1.867 × 10−11 (Söderlund et al. Reference Söderlund, Patchett, Vervoort and Isachsen2004) and the Lu–Hf Chondritic Uniform Reservoir (CHUR) parameters of Bouvier et al. (Reference Bouvier, Vervoort and Patchett2008). The initial 176Hf/177Hf isotope ratios and ϵHf values for each detrital zircon were calculated using the 207Pb/206Pb date of each analytical spot.

5. Results

5.a. Sample 17GC1086 (Montauban Group)

The examined zircon grains are rounded to sub-rounded with aspect ratios 1:1 to 4:1, exhibiting oscillatory, sector-zoned or homogeneous zones in CL with brighter cores commonly overgrown by thin (<15 μm) darker rims (Fig. 2; spot 1086R7.3). Planar banding and domains with faint variations in CL response are also observed. Four dominant age peaks at 1.46 Ga, 1.62 Ga, 1.87 Ga and 2.7 Ga are recorded in this sample (Fig. 3; sample 17GC1086). The dominant age components at c. 1.46 Ga and 1.62 Ga indicate the prevalence of sources with Pinwarian (34 % of the analyses) and Labradorian (29 % of the analyses) affinities, with subordinate amounts of grains with Penokean ages (22 % of the analyses). The MLA algorithm yields a date of 1442 ± 8 Ma (Fig. 4a) that is interpreted as the MDA of the sample. This date agrees with the youngest cluster of nine U–Pb analyses from grains with oscillatory zoning in CL and Th/U values from 0.3 up to 1.1 that yield a weighted mean 207Pb/206Pb date of 1436 ± 11 Ma (n = 9, MSWD = 0.9, 2s). A cluster of younger 207Pb/206Pb dates from 1399 Ma to 1334 Ma is also observed, but the zoning patterns in CL of the grains are not diagnostic of igneous growth and were not considered in the estimation of the MDA (see also Supplementary Material, available online at https://doi.org/10.1017/S0016756821001023).

Fig. 2. Representative cathodoluminescence images from detrital zircon grains of the four quartzite samples. Red and blue circles show the location of the 25 μm and 50 μm beam location for U–Pb and Lu–Hf microanalysis, respectively.

Fig. 3. Kernel density estimation (KDE) 207Pb/206Pb age diagrams and weighted mean 207Pb/206Pb age diagrams from the youngest zircon grains in the four quartzite samples. The numbers adjacent to the peaks denote the age maxima in Ga. The weighted mean plots next to the KDE diagrams represent MDA estimates based on the youngest acceptable dates for each supracrustal complex.

Fig. 4. Radial plots and MDA estimates of the (a) Montauban Group quartzite, (b) Wabash Complex quartzite and (c, d) Barrois Complex quartzites using isoplotR (Vermeesch, Reference Vermeesch2018) and the Maximum Likelihood Age (MLA) algorithm (Vermeesch, Reference Vermeesch2021). The radial plots display heteroscedastic data (i.e. data points with unequal uncertainties) where each data point represents a 207Pb/206Pb date. These plots are scatter plots of standardized estimates (y axis) against precisions (x axis).The angular position of each data point denotes the date, whereas the horizontal distance from the point of origin varies with the precision of the analysis (i.e. imprecise measurements on the left side and more precise on the right side of the plot). For further details on age calculation see Vermeesch (Reference Vermeesch2021) and Galbraith (Reference Galbraith1988).

The ϵHf values of the analysed zircon grains vary from −18.4 to +10.7 (Fig. 5). The majority of zircon grains with 207Pb/206Pb ages between c.1.4 and 1.56 Ga (92 % of the analyses) exhibit suprachrondritic ϵHf values (ϵHf > 0) consistent with the juvenile signature of Pinwarian sources. Similarly, 84 % of the analyses from the 1.6–1.7 Ga age population also yield suprachondritic ϵHf values (ϵHf = +0.2 to +5.8) with 16 % recording ϵHf values between −1.1 and −8. The third major age component, with ages between c. 1.8 and 1.9 Ga, shows an increasing contribution from sources that define a vertical mixing array between sub- and suprachondritic components (ϵHf between −1.2 and +7.0). The minor 2.7 Ga age component consists chiefly of grains with ϵHf values that vary from −10.3 to +3.8.

Fig. 5. Epsilon Hf versus age with Lu–Hf isotopic data from the detrital zircon grains of the four quartzites reported in this study. Black and white coloured are Lu–Hf data points from the dataset of Spencer et al. (Reference Spencer, Cawood, Hawkesworth, Prave, Roberts, Horstwood and Whitehouse2015). The black-filled data points represent data from the pre-orogenic Siamarnekh Formation (SF). The white-filled data points depict, for comparative purposes, Lu–Hf isotopic data from pre-, syn- and post-orogenic formations from Labrador and Scotland. The lines with different slopes (DM = depleted mantle, MC = mafic crust, CC = continental crust) correspond to different Lu–Hf isotope evolution curves for the depleted mantle (176Lu/177Hf = 0.03976; Vervoort et al. Reference Vervoort, Kemp and Fisher2018), continental crust (176Lu/177Hf = 0.0113; Rudnick & Gao, Reference Rudnick and Gao2003) and mafic crust (176Lu/177Hf = 0.022; Nebel et al. Reference Nebel, Nebel-Jacobsen, Mezger and Berndt2007) reservoirs. The CHUR reference line (ϵHf = 0) involves an uncertainty of ±0.4 ϵHf units (Bouvier et al. Reference Bouvier, Vervoort and Patchett2008). The grey band is drawn parallel to the 176Lu/177Hf = 0.022 line and covers c. 70 % of the 1.8–1.9 Ga zircon grains.

5.b. Sample 18YD2075 (Wabash complex)

Two grain-size populations consisting of euhedral zircon grains with aspect ratios varying from 1:1 up to 3:1 are observed in this sample. The examined grains exhibit bright cores in CL, with oscillatory or sector zoning, overgrown by darker euhedral rims of possible metamorphic origin (Fig. 2). Zircon grains with homogeneous grey to dark grey response in CL are also common. The region between the brighter rounded to sub-rounded cores and the darker rims often has a brighter response in CL, indicating a separate growth event (Fig. 2; 2075R1.3). The four age components recorded in the detrital spectrum are c. 1.61 Ga, 1.80 Ga, 1.90 Ga and 2.68 Ga (Fig. 3; 18YD2075). The MLA algorithm yields a date of 1309 ± 38 Ma (Fig. 4b) that is interpreted as the MDA of the sample, corroborating previous MDA constraints of this complex (Moukhsil et al. Reference Moukhsil, Solgadi, Belkacim, Augland and David2015). In comparison, the two youngest dates with Th/U values of 0.14 and 0.19 agree with the MDA yielding a weighted average 207Pb/206Pb date of 1316 ± 150 (2s).

The ϵHf values from the analysed zircon grains vary from −26.9 to +10.5 (Fig. 5). Specifically, the Lu–Hf analyses from the minor 1.4–1.5 Ga and c.1.6–1.7 Ga age populations have ϵHf values varying from +1.2 to +10.5 and −12.4 to +8.1, respectively. The majority of the analysed grains belong to the 1.7–2.0 Ga age range (56% of the analyses) and exhibit a spread in ϵHf values from −26.9 to +7.2. Zircon grains with ages >2 Ga (30% of the analyses) have ϵHf values that vary from −10.6 to +7.5.

5.c. Sample 19AM02B (Barrois complex)

The zircon grains of the 19AM02 quartzite are rounded to sub-rounded with aspect ratios 1:1 to 2:1. The examined grains exhibit bright cores in CL, commonly with sector zoning, overgrown by rims (≤20 μm thick) that are dark grey in CL (Fig. 2; spot 02R10.15). Two dominant age peaks are recorded in the detrital spectrum at 1.85 Ga and 2.65 Ga (Fig. 3; 19AM02B). The MLA algorithm yields a date of 1202 ± 40 Ma (Fig. 4c) that is interpreted as the MDA of the sample. In comparison, the three youngest dates from zircon grains with Th/U values between 0.4 and 1, that overlap at 2s, yield a weighted average 207Pb/206Pb date of 1303 ± 11 Ma.

The ϵHf values on the examined grains vary from −12.3 to +10. In more detail, two analyses of the youngest age population at c. 1.2 Ga and 1.3 Ga yield ϵHf values of −11.3 and −12, respectively. The Lu–Hf analyses on grains of the c. 1.8–2.0 Ga and 2.3–3.0 Ga age populations exhibit ϵHf values that vary from −12.3 to +10.0 and −8.6 to +7.5 (Fig. 5).

5.d. Sample 19FS7085 (Barrois complex)

The examined zircon grains are rounded to sub-rounded with aspect ratios 1:1 to 3:1. The imaged grains exhibit bright cores in CL overgrown by oscillatory zoned mantle domains (Fig. 2; 7085R3.11); domains with dark grey response in CL that exhibit oscillatory or no zoning are also observed. A dominant age peak is recorded in the sample at 1.83 Ga, with two subordinate peaks at 1.44 Ga and 2.63 Ga (Fig. 3; 19FS7085). The MLA algorithm yields a date of 1281 ± 48 Ma (Fig. 4d) that is interpreted as the MDA of the sample. The three youngest U–Pb dates with Th/U values from 0.47 to 0.84 yield a weighted average 207Pb/206Pb date of 1294 ± 64 (2s) which is consistent with the calculated MDA.

The ϵHf values on the examined grains vary from −17.8 to +10.5 (Fig. 5). Specifically, the Lu–Hf analyses in the 1.7–1.9 Ga age population show the widest spread in ϵHf values from −17.8 to +10.4 with a prevalence of negative ϵHf values (84% of the analyses). Four U–Pb analyses that fall within the 1.6–1.7 Ga age population show a smaller spread in ϵHf values from −3.66 to +5.56. Seven Lu–Hf analyses from grains of the c. 1.4–1.5 Ga population vary from −17.5 to +8.8, with a prevalence of positive ϵHf values (72 % of the analyses), whereas those with ages >2.3 Ga yielded ϵHf values varying from −15.3 to +3.7.

6. Discussion

6.a. Timeline of crustal evolution

6.a.1 Juvenile crustal growth at c. 1.4–1.5 Ga and 1.6–1.7 Ga

The suprachondritic Hf isotopic values in the c. 1.4–1.5 Ga age population support juvenile crustal additions during the accretionary growth of southeastern Laurentia or reworking of a juvenile source with short crustal residence time. In the case of the 1.6–1.7 Ga age population the presence of subchondritic ϵHf values is consistent with an older crustal component assimilated during mantle-derived magmatism. Juvenile whole-rock Nd isotopic compositions of Pinwarian and Labradorian orthogneisses exhibit ϵNd values that vary from +2.2 to +5.8 and +2.3 to +3.6, respectively (Dickin, Reference Dickin2000; Zelek & Dickin, Reference Zelek and Dickin2013). Geochemical and zircon Lu–Hf isotopic data from other Pinwarian metaplutonic suites in the Grenville Province also show juvenile signatures with ϵHf values varying from 6.3 ± 0.5 to 10.7 ± 0.7 (Augland et al., Reference Augland, Moukhsil, Solgadi and Indares2015), supporting the Lu–Hf detrital zircon dataset of this study.

It has been proposed that juvenile plutonic rocks have low preservation potential during supercontinent assembly (Hawkesworth et al. Reference Hawkesworth, Dhuime, Pietranik, Cawood, Kemp and Storey2010; Spencer et al. Reference Spencer, Cawood, Hawkesworth, Prave, Roberts, Horstwood and Whitehouse2015). However, in our study the age peaks at c. 1.45 Ga and 1.65 Ga in the detrital spectra record accretionary orogenesis that survived the effects of the protracted collisional Grenvillian orogenesis. These peri-Laurentian accretionary events are members of the Great Proterozoic Accretionary Orogen of Laurentia that contains the largest fraction of juvenile continental crust produced globally during the Palaeo- and Mesoproterozoic (Condie, Reference Condie2013). Thus, here we show that, in Central Grenville Province, the record of juvenile crust production prior to collision is preserved in vestiges of pre-orogenic basins, supporting the few U–Pb and Lu–Hf zircon studies from pristine meta-igneous lithologies (Augland et al., Reference Augland, Moukhsil, Solgadi and Indares2015).

6.a.2 Crustal reworking at 1.8–1.9 Ga

The vertical mixing array in ϵHf-time space is consistent with the assimilation of mantle-derived magmas with pre-existing crust during the Penokean–Makkovikian orogeny (Fig. 5). Alternatively, the suprachondritic ϵHf values at 1.8–1.9 Ga age could represent a distinct juvenile isotopic source crystallized between 1.8 and 1.9 Ga, with the subchondritic ϵHf data reflecting pure crustal reworking. Neither scenario can be excluded since juvenile intra-oceanic and continental arcs that reworked pre-existing crust were constructed between 1.8 and 1.9 Ga (e.g. Labrador; Culshaw et al. Reference Culshaw, Ketchum and Barr2000) and could have shed detritus into the examined pre-collisional basins.

In our dataset, a 176Lu/177Hf evolution trajectory with slope of 0.011 between a near-chondritic 2.7 Ga reservoir and the end-member ϵHf value of the sub-chondritic source at 1.8–1.9 Ga (arrow line; Fig. 5) is consistent with the reworking of Archaean (i.e. Superior Province) continental crust. The ϵHf-time pattern could alternatively be interpreted as the reworking of variably contaminated mafic sources extracted from a depleted-mantle reservoir between c. 2.8 Ga and 3.6 Ga (176Lu/177Hf = 0.022; grey band in Fig. 5). We lean towards the former scenario, considering the development of an Andean-type convergent margin during the Penokean orogeny at the southern margin of Laurentia (Schulz & Cannon., Reference Schulz and Cannon2007) and the presence of 2.7 Ga chondritic ϵHf detrital zircon values. The only previous Lu–Hf isotopic data from detrital zircon grains of a pre-orogenic formation come from the Siamarnekh formation with MDA at c.1.25 Ga (Spencer et al. Reference Spencer, Cawood, Hawkesworth, Prave, Roberts, Horstwood and Whitehouse2015) that crops out in northern Labrador. The zircon grains of the Siamarnekh formation exhibit two dominant age peaks at 1.9 Ga and 2.15 Ga and show a vertical spread in ϵHf-time space. Interestingly, the 2.15 Ga grains fall along the crustal reworking trajectory (176Lu/177Hf = 0.011; Fig. 5), but further constraints about their origin and provenance have not been reported. Also, the most proximal Penokean-affinity sources, from the Western Grenville Province, that could have shed detritus towards the sequences examined here are gneisses of the Barillia terrane with Nd model ages in the 1.8–2.0 Ga range. These gneisses contain no Nd isotopic evidence of reworked Archaean crustal components and have been characterized as juvenile (Dickin, Reference Dickin2000). However, lead isotopic data from leached feldspar crystals from various lithologies in the western Grenville Province suggest mixing of Archaean crust with 1.8 or 1.5 Ga juvenile magmatic rocks, with the latter preferred in order to be compatible with the 1.8 Ga Nd model ages in the area (DeWolf & Mezger, Reference DeWolf and Mezger1994). More distal terranes, in the Great Lakes region and Eastern Labrador, exhibit stronger evidence for Penokean–Makkovikian reworking of Archaean basement crust. In the Great Lakes region, synorogenic Penokean plutons that intrude the Archaean Marshfield terrane show ϵNd values that vary from −1 to −7 resulting from the contamination of magmas with Archaean crust, supporting previous isotopic studies (Van Schmus & Anderson, Reference Van Schmus and Anderson1977; Van Wyck & Johnson, Reference Van Wyck and Johnson1997). Further eastwards, in the Makkovik Province of eastern Labrador, the Island Harbour Bay Plutonic Suite constitutes a calc-alkaline magmatic arc that reworked Archaean crust between c. 1.89 and 1.87 Ga (Ketchum et al. Reference Ketchum, Culshaw and Barr2002). Similarly, Makkovikian granitoids in this area show negative ϵNd values that range from −7.2 to −2.5 as a result of variable assimilation of granitic magmas with Archaean crust of the Nain Province (Barr et al. Reference Barr, White, Culshaw and Ketchum2001). Taking into consideration that the 1.8–1.9 Ga grains are not first cycle and there is no U–Pb isotopic evidence for 1.8–1.9 Ga magmatism proximally to the study area, we attribute the parental sources of the 1.8–1.9 Ga grains to reworked Archaean basement terranes of the Great Lakes area or Makkovik Province.

6.b. Provenance and tectonic constraints

6.b.1. Montauban quartzite

The suprachondritic Hf isotopic compositions in the zircon grains of the c.1.46 Ga age population indicate they were derived from juvenile sources. Potential sources are the c. 1.45 Ga Montauban arc (L Nadeau et al., unpub. data, 1992), an intraoceanic arc that accreted to Laurentia at c. 1.39 Ga (Sappin et al. Reference Sappin, Constantin, Clark and van Breemen2009), and juvenile Pinwarian granitoids (Augland et al., Reference Augland, Moukhsil, Solgadi and Indares2015). The detrital spectrum in Montauban quartzite and the c. 5–10 Ma age gap between the dominant age peak and the MDA would be consistent with deposition in a convergent margin (possibly a back-arc) setting (e.g. Cawood et al. Reference Cawood, Hawkesworth and Dhuime2012). The mingling of clastic, pyroclastic and pillow basalts in the Montauban Group likely records the early stages of back-arc basin formation at c. 1.44 Ga (Fig. 6a). The Escoumins supracrustal belt is a laterally equivalent supracrustal sequence to the Montauban Group, with a similar MDA of 1476 ± 10 Ma (Groulier et al. Reference Groulier, Indares, Dunning, Moukhsil and Wälle2018). The two sequences have age spectra with dominance from Pinwarian sources, but in the Montauban quartzite there is also a significant contribution from Labradorian sources (Fig. 6b). The minor age peak in both sequences between c. 1.8 and 1.9 Ga indicates the limited contribution from Penokean–Makkovikian sources. Comparatively, the detrital signal of Montauban and Escoumins sequences is markedly different from quartzites of the Pinwarian, Wakeham intracontinental rift basin, in eastern Grenville Province (Van Breemen & Corriveau, Reference Van Breemen and Corriveau2005). The unimodal detrital spectra with an age peak at 1.8–1.9 Ga in Wakeham basin quartzites (Fig. 6c) indicate provenance most probably from the Palaeoproterozoic Torngat orogen or the Makkovikian terrane. The disparate detrital age distributions demonstrate the different tectonic setting and catchment areas of the Pinwarian basins along the orogen.

Fig. 6. (a) Tectonic setting (not in scale) of the examined quartzite from the supracrustal sequence of the Montauban Group; (b) KDE plots of U–Pb dates from the Montauban and Escoumins supracrustal group (Groulier et al. Reference Groulier, Indares, Dunning, Moukhsil and Wälle2018); (c) KDE plots of U–Pb dates from the Wakeham basin (U–Pb data from Wodicka et al. Reference Wodicka, David, Parent, Gobeil and Verpaelst2003 and Van Breemen & Corriveau., Reference Van Breemen and Corriveau2005); (d) tectonic setting (not in scale) of the Wabash–Barrois supracrustal groups; (e) KDE plots of U–Pb dates from the Wabash–Barrois and Saint-Boniface supracrustal groups (U–Pb data from this study and Peck et al. Reference Peck, Quinan and Selleck2019); (f) KDE plots of U–Pb dates from supracrustal sequences of the Adirondack lowlands (U–Pb data from Chiarenzelli et al. Reference Chiarenzelli, Kratzmann, Selleck and de Lorraine2015; and pers. comm. with J. Chiarenzelli).

6.b.2. Wabash and Barrois quartzites

The age gap between the MDA and the major age peaks at 1.8–1.9 Ga in Barrois and Wabash complexes is more than 500 Ma. Specifically, c. 97 % of the analyses record an age gap of more than 150 Ma between the dominant age peak and the MDA (Fig. 6d), suggesting that both sequences were deposited in intracratonic extensional sedimentary basins (e.g. Cawood et al. Reference Cawood, Hawkesworth and Dhuime2012). In terms of proto-sources, Penokean–Makkovikian (c. 1.8–1.9 Ga) lithologies to our knowledge are not exposed in the study area. Therefore, the 1.8–1.9 Ga detrital zircon grains, hosted in the Wabash and Barrois quartzites, either constitute distal detritus deposited in intermediate sedimentary reservoirs and recycled during Elzevirian rifting (1.2–1.3 Ga) or were transported via orogen-parallel drainage (present coordinates), possibly in large fluvial systems. U–Pb detrital zircon spectra with dominant age peaks between c. 1.8 and 1.9 Ga and MDA at 1.2–1.3 Ga, similarly with the Wabash and Barrois complexes, are observed within different crustal domains of the Grenville Province. Specifically, southwestwards of the study area, in the Morin terrane, the Saint-Boniface metasedimentary sequence (Corrigan & Van Breemen, Reference Corrigan and Van Breemen1997) with MDA at c. 1.25 Ga (Peck et al. Reference Peck, Quinan and Selleck2019) shows dominant age-peaks at c. 1.25 and 1.82 Ga (Fig. 6e). These age peaks also indicate the feeding of these basins dominantly from Penokean and Elzevirian sources. Further west, in the Adirondack region, three quartzites from metasedimentary sequences with MDAs at c. 1.25–1.3 Ga contain a higher proportion of Elzevirian and a lower proportion of Penokean detrital zircon grains than those found in Quebec (Fig. 6f; Chiarenzelli et al. Reference Chiarenzelli, Kratzmann, Selleck, de Lorraine and Lupulescu2017). These 1.2–1.3 Ga basins support inferences from bimodal volcanic and AMCG plutonic suites which indicate that extensional tectonism was occurring during that period in the Grenville Province and on a global scale at the core of Nuna supercontinent (Evans & Mitchell, Reference Evans and Mitchell2011; Indares & Moukhsil, Reference Indares and Moukhsil2013).

7. Conclusions

From the examination of three poorly exposed Mesoproterozoic supracrustal complexes from Precambrian basins of the Central Grenville Province we show that:

  • Detrital zircon U–Pb data from a quartzite hosted within a metavolcanic sequence (Montauban complex) indicates the development of a Pinwarian back-arc basin with maximum depositional age of 1442 ± 8 Ma (2s).

  • The detrital zircon grains of the Montauban quartzite record two dominant age peaks at c. 1.46 Ga and 1.62 Ga that exhibit suprachondritic ϵHf values recording juvenile crust formation during the Pinwarian (c. 1.4–1.5 Ga) and Labradorian (c. 1.6–1.7 Ga) orogenies.

  • The U–Pb detrital zircon spectra from two supracrustal complexes in the Lac Saint-Jean region (Wabash and Barrois complexes) with MDA of c. 1.2–1.3 Ga suggest deposition of 1.8–1.9 Ga recycled or far-travelled detrital zircon grains in an intracontinental extensional basin.

  • The sub- and suprachondritic ϵHf values found in the 1.8–1.9 Ga detrital zircon grains manifest the reworking of a near-chondritic 2.7 Ga continental basement during the Andean-type, Penokean–Makkovikian orogeny.

  • The accretionary stages of supercontinent assembly are retained in the detrital zircon archive of pre-collisional basins in Central Grenville Province even after the terminal Grenvillian collision at c. 1 Ga.

Acknowledgements

KP thanks R Stevenson, and A Indares for comments and suggestions on early versions of the manuscript, and J Gogot for technical assistance. The comments from anonymous reviewers and François Turlin enhanced the manuscript and are appreciated. This is Ministère de l’Energie et des Ressources naturelles publication number 8449-2016-01.

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756821001023

References

Augland, LE, Moukhsil, A, Solgadi, F and Indares, A (2015) Pinwarian to Grenvillian magmatic evolution in the central Grenville Province: new constraints from ID–TIMS U–Pb ages and coupled Lu–Hf S–MC–ICP–MS data. Canadian Journal of Earth Sciences 52(9), 701–21. doi: 10.1139/cjes-2014-0232.CrossRefGoogle Scholar
Barovich, KM, Patchett, PJ, Peterman, ZE and Sims, PK (1989) Nd isotopes and the origin of 1.9-1.7 Ga Penokean continental crust of the Lake Superior region. Geological Society of America Bulletin 101(3), 333–8. doi: 10.1130/0016-7606(1989)101<0333:NIATOO>2.3.CO;2.2.3.CO;2>CrossRefGoogle Scholar
Barr, SM, White, CE, Culshaw, NG and Ketchum, JWF (2001) Geology and tectonic setting of Paleoproterozoic granitoid suites in the Island Harbour Bay area, Makkovik Province, Labrador. Canadian Journal of Earth Sciences 38(3), 441–63. doi: 10.1139/cjes-38-3-441.CrossRefGoogle Scholar
Belousova, EA, Kostitsyn, YA, Griffin, WL, Begg, GC, O’Reilly, SY and Pearson, NJ (2010) The growth of the continental crust: constraints from zircon Hf-isotope data. Lithos 119, 457–66. doi: 10.1016/j.lithos.2010.07.024.CrossRefGoogle Scholar
Bernier, L, Pouliot, G and MacLean, WH (1987) Geology and metamorphism of the Montauban north gold zone: a metamorphosed polymetallic exhalative deposit, Grenville province, Quebec. Economic Geology 82, 2076–90. doi: 10.2113/gsecongeo.82.8.2076.CrossRefGoogle Scholar
Blichert-Toft, J, Chauvel, C and Albarède, F (1997) Separation of Hf and Lu for high-precision isotope analysis of rock samples by magnetic sector-multiple collector ICP-MS. Contributions to Mineralogy and Petrology 127, 248–60. doi: 10.1007/s004100050278.CrossRefGoogle Scholar
Bouvier, A, Vervoort, JD and Patchett, PJ (2008) The Lu-Hf and Sm-Nd isotopic composition of CHUR: constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth and Planetary Science Letters 273, 48–57. doi: 10.1016/j.epsl.2008.06.010.CrossRefGoogle Scholar
Cawood, PA, Hawkesworth, CJ and Dhuime, B (2012) Detrital zircon record and tectonic setting. Geology 40(10), 875–78. doi: 10.1130/G32945.1.CrossRefGoogle Scholar
Chiarenzelli, J, Kratzmann, D, Selleck, B and de Lorraine, W (2015) Age and provenance of Grenville supergroup rocks, Trans-Adirondack Basin, constrained by detrital zircons. Geology 43(2), 183–6. doi: 10.1130/G36372.1.CrossRefGoogle Scholar
Chiarenzelli, J, Kratzmann, D, Selleck, B, de Lorraine, W and Lupulescu, M (2017) Utility of detrital zircon grains from upper amphibolite facies rocks of the Grenville Supergroup, Adirondack Lowlands, Northeastern United States. In Sediment Provenance: Influences on Compositional Change from Source to Sink (ed. R Mazumder), 371–401. doi: 10.1016/B978-0-12-803386-9.00014-9.Google Scholar
Chu, NC, Taylor, RN, Chavagnac, V, Nesbitt, RW, Boella, RM, Milton, JA, German, CR, Bayon, G and Burton, K (2002). Hf isotope ratio analysis using multi-collector inductively coupled plasma mass spectrometry: an evaluation of isobaric interference corrections. Journal of Analytical Atomic Spectrometry 17, 1567–74. doi: 10.1039/b206707b.CrossRefGoogle Scholar
Condie, KC (2013) Preservation and recycling of crust during accretionary and collisional phases of proterozoic orogens: a bumpy road from Nuna to Rodinia. Geosciences (Switzerland) 3(2), 240–61. doi: 10.3390/geosciences3020240.Google Scholar
Condie, KC, Beyer, E, Belousova, E, Griffin, WL and O’Reilly, SY (2005) U-Pb isotopic ages and Hf isotopic composition of single zircons: the search for juvenile Precambrian continental crust. Precambrian Research 139, 42–100. doi: 10.1016/j.precamres.2005.04.006.CrossRefGoogle Scholar
Connelly, JN and Heaman, LM (1993) U-Pb geochronological constraints on the tectonic evolution of the Grenville Province, western Labrador. Precambrian Research 63, 123–42. doi: 10.1016/0301-9268(93)90008-P.CrossRefGoogle Scholar
Corrigan, D and Van Breemen, O (1997) U-Pb age constraints for the lithotectonic evolution of the Grenville Province along the Mauricie transect, Quebec. Canadian Journal of Earth Sciences 34(3), 229–316. doi: 10.1139/e17-027.CrossRefGoogle Scholar
Culshaw, N, Foster, J, Marsh, J, Slagstad, T and Gerbi, C (2016) Kiosk domain, Central Gneiss Belt, Grenville Province, Ontario: a Labradorian palimpsest preserved in the ductile deep crust. Precambrian Research 280, 249–78. doi: 10.1016/j.precamres.2016.05.002.CrossRefGoogle Scholar
Culshaw, N, Ketchum, J and Barr, S (2000) Structural evolution of the Makkovik Provinces, Labrador, Canada: tectonic processes during 200 Myr at a Paleoproterozoic active margin. Tectonics 19(5), 961–77. doi: 10.1029/1999TC001156.CrossRefGoogle Scholar
David, J, Moukhsil, A and Dion, C (2010) Datations U-Pb effectuées dans la Province de Grenville en 2008–2009. RP 2010-10. Québec: Ministère des Ressources naturelles, 17 pp.Google Scholar
DeWolf, CP and Mezger, K (1994) Lead isotope analyses of leached feldspars: constraints on the early crustal history of the Grenville Orogen. Geochimica et Cosmochimica Acta 58(24), 5537–50. doi: 10.1016/0016-7037(94)90248-8.CrossRefGoogle Scholar
Dickin, AP and McNutt, RH (1989). Nd model age mapping of the southeast margin of the Archean foreland in the Grenville province of Ontario. Geology 17(4), 299302. doi: 10.1130/0091-7613(1989)017<0299:NMAMOT>2.3.CO;2 2.3.CO;2>CrossRefGoogle Scholar
Dickin, AP (2000) Crustal formation in the Grenville Province: Nd-isotope evidence. Canadian Journal of Earth Sciences 37(2–3), 165–81. doi: 10.1139/e99-039.CrossRefGoogle Scholar
Dickin, AP, McNutt, RH, Martin, C and Guo, A (2009) The extent of juvenile crust in the Grenville Province: Nd isotope evidence. Geological Society of America Bulletin 122, 870–83. doi: 10.1130/B26381.1.CrossRefGoogle Scholar
Dunning, G and Indares, A (2010) New insights on the 1.7-1.0 Ga crustal evolution of the central Grenville Province from the Manicouagan – Baie Comeau transect. Precambrian Research 180(3–4), 204–26. doi: 10.1016/j.precamres.2010.04.005.CrossRefGoogle Scholar
Easton, RM (1992) The Grenville Province and the Proterozoic history of central and southern Ontario. In Geology of Ontario Special, vol. 4, part 2 (eds Thurston, PC, Williams, HR, Sutcliffe, RH and Stott, GM), pp. 715904. Ontario: Ontario Ministry of Northern Development and Mines.Google Scholar
Evans, DAD and Mitchell, RN (2011) Assembly and breakup of the core of Paleoproterozoic-Mesoproterozoic supercontinent Nuna. Geology 39(5), 443–6. doi: 10.1130/G31654.1.CrossRefGoogle Scholar
Fisher, CM, Hanchar, JM, Samson, SD, Dhuime, B, Blichert-Toft, J, Vervoort, JD and Lam, R (2011) Synthetic zircon doped with hafnium and rare earth elements: a reference material for in situ hafnium isotope analysis. Chemical Geology 286(1–2), 32–47. doi: 10.1016/j.chemgeo.2011.04.013.CrossRefGoogle Scholar
Fisher, CM and Vervoort, JD (2018) Using the magmatic record to constrain the growth of continental crust: —the Eoarchean zircon Hf record of Greenland. Earth and Planetary Science Letters 488, 79–91. doi: 10.1016/j.epsl.2018.01.031.CrossRefGoogle Scholar
Fisher, CM, Vervoort, JD and Hanchar, JM (2014) Guidelines for reporting zircon Hf isotopic data by LA-MC-ICPMS and potential pitfalls in the interpretation of these data. Chemical Geology 363, 125–33. doi: 10.1016/j.chemgeo.2013.10.019.CrossRefGoogle Scholar
Gain, SEM, Gréau, Y, Henry, H, Belousova, E, Dainis, I, Griffin, WL and O’Reilly, SY (2019) Mud Tank Zircon: long-term evaluation of a reference material for U-Pb dating, Hf-isotope analysis and trace element analysis. Geostandards and Geoanalytical Research 43, 339–54. doi: 10.1111/ggr.12265.CrossRefGoogle Scholar
Galbraith, RF (1988) Graphical display of estimates having differing standard errors. Technometrics 30, 271. doi: 10.2307/1270081.CrossRefGoogle Scholar
Gower, CF (2019) Regional Geology of Eastern Labrador (Eastern Makkovik and Grenville Provinces). St. John’s: Government of Newfoundland and Labrador, Department of Natural Resources, Geological Survey, Memoir 4, 654 pp. (Geofile LAB/1752.).Google Scholar
Gower, CF and Krogh, TE (2002) A U-Pb geochronological review of the Proterozoic history of the eastern Grenville Province. Canadian Journal of Earth Sciences 39(5), 795–829. doi: 10.1139/e01-090.CrossRefGoogle Scholar
Groulier, PA, Indares, A, Dunning, G and Moukhsil, A (2020a) Andean style 1.50–1.35 Ga arc dynamics in the Southeastern Laurentian margin: the rifting and reassembly of Quebecia. Terra Nova 32(6), 450–7. doi: 10.1111/ter.12482.CrossRefGoogle Scholar
Groulier, PA, Indares, A, Dunning, G, Moukhsil, A and Wälle, M (2018) Peri-Laurentian, Pinwarian-age oceanic arc crust preserved in the Grenville Province: insights from the Escoumins supracrustal belt. Precambrian Research 311, 37–64. doi: 10.1016/j.precamres.2018.04.001.CrossRefGoogle Scholar
Groulier, PA, Turlin, F, André-Mayer, AS, Ohnenstetter, D, Crépon, A, Boulvais, P, Poujol, M, Rollion-Bard, C, Zeh, A, Moukhsil, A, Solgadi, F and El basbas, A (2020b) Silicate-carbonate liquid immiscibility: insights from the Crevier alkaline intrusion (Quebec). Journal of Petrology 61(3). doi: 10.1093/petrology/egaa033.CrossRefGoogle Scholar
Guo, A and Dickin, AP (1996) The southern limit of Archean crust and significance of rocks with Paleoproterozoic model ages: Nd model age mapping in the Grenville Province of western Quebec. Precambrian Research 77(3–4), 231–41. doi: 10.1016/0301-9268(95)00095-x.CrossRefGoogle Scholar
Hawkesworth, CJ, Dhuime, B, Pietranik, AB, Cawood, PA, Kemp, AIS and Storey, CD (2010) The generation and evolution of the continental crust. Journal of the Geological Society 167, 229–48. doi: 10.1144/0016-76492009-072.CrossRefGoogle Scholar
Hébert, C and Van Breemen, O (2004) Mesoproterozoic basement of the Lac St. Jean Anorthosite Suite and younger Grenvillian intrusions in the Saguenay region, Québec: structural relationships and U-Pb geochronology. Memoir of the Geological Society of America 197, 65–79. doi: 10.1130/0-8137-1197-5.65.Google Scholar
Higgins, MD and Van Breemen, O (1996) Three generations of anorthosite-mangerite-charnockite-granite (AMCG) magmatism, contact metamorphism and tectonism in the Saguenay-Lac-Saint-Jean region of the Grenville Province, Canada. Precambrian Research 79(3–4), 327–46. doi: 10.1016/0301-9268(95)00102-6.CrossRefGoogle Scholar
Hoffman, P (1988) United plates of America, the birth of a craton: early Proterozoic assembly and growth of Laurentia. Annual Review of Earth and Planetary Sciences 16, 543–603. doi: 10.1146/annurev.earth.16.1.543.CrossRefGoogle Scholar
Hoffmann, JE, Münker, C, Polat, A, Rosing, MT and Schulz, T (2011) The origin of decoupled Hf-Nd isotope compositions in Eoarchean rocks from southern West Greenland. Geochimica et Cosmochimica Acta 75(21), 6610–28. doi: 10.1016/j.gca.2011.08.018.CrossRefGoogle Scholar
Huang, C, Wang, H, Yang, JH, Ramezani, J, Yang, C, Zhang, SB, Yang, YH, Xia, XP, Feng, LJ, Lin, J, Wang, TT, Ma, Q, He, HY, Xie, LW and Wu, ST (2020) SA01 – a proposed Zircon reference material for Microbeam U-Pb age and Hf-O isotopic determination. Geostandards and Geoanalytical Research 44(1), 103–23. doi: 10.1111/ggr.12307.CrossRefGoogle Scholar
Indares, A, Dunning, G, Cox, R, Gale, D and Connelly, J (1998) High-pressure, high-temperature rocks from the base of thick continental crust: geology and age constraints from the Manicouagan Imbricate Zone, eastern Grenville Province. Tectonics 17(3), 426–40. doi: 10.1029/98TC00373.CrossRefGoogle Scholar
Indares, A and Moukhsil, A (2013) Geon 12 crustal extension in the central Grenville Province, implications for the orogenic architecture, and potential influence on the emplacement of anorthosites. Canadian Journal of Earth Sciences 50(9), 955–66. doi: 10.1139/cjes-2012-0161.CrossRefGoogle Scholar
Jackson, SE, Pearson, NJ, Griffin, WL and Belousova, EA (2004) The application of laser ablation-inductively coupled plasma-mass spectrometry to in situ U-Pb zircon geochronology. Chemical Geology 211, 4769. doi: 10.1016/j.chemgeo.2004.06.017.CrossRefGoogle Scholar
Jannin, S, Gervais, F, Moukhsil, A, Augland, LE and Crowley, JL (2018) Déformations tardi-grenvilliennes dans la ceinture parautochtone (province de grenville centrale): contraintes géochronologiques par couplage de méthodes U-Pb de haute résolution spatiale et de haute précision. Canadian Journal of Earth Sciences 55, 406–35. doi: 10.1139/cjes-2017-0129.CrossRefGoogle Scholar
Jordan, SL, Indares, A and Dunning, G (2006) Partial melting of metapelites in the Gagnon terrane below the high-pressure belt in the Manicouagan area (Grenville Province): pressure-temperature (P-T) and U-Pb age constraints and implications. Canadian Journal of Earth Sciences 43(9), 1309–29. doi: 10.1139/E06-038.CrossRefGoogle Scholar
Kamo, SL, Heaman, LM and Gower, CF (2011) Evidence for post-1200 Ma – pre-Grenvillian supracrustal rocks in the Pinware terrane, eastern Grenville Province at Battle Harbour, Labrador. Canadian Journal of Earth Sciences 48, 371–87. doi: 10.1139/E10-052.CrossRefGoogle Scholar
Kemp, AIS, Hawkesworth, CJ, Paterson, BA and Kinny, PD (2006) Episodic growth of the Gondwana supercontinent from hafnium and oxygen isotopes in zircon. Nature 439, 580–83. doi: 10.1038/nature04505.CrossRefGoogle Scholar
Kemp, AIS, Hickman, AH, Kirkland, CL and Vervoort, JD (2015) Hf isotopes in detrital and inherited zircons of the Pilbara Craton provide no evidence for Hadean continents. Precambrian Research 261, 112–26. doi: 10.1016/j.precamres.2015.02.011.CrossRefGoogle Scholar
Ketchum, JWF, Culshaw, NG and Barr, SM (2002) Anatomy and orogenic history of a Paleoproterozoic accretionary belt: the Makkovik Province, Labrador, Canada. Canadian Journal of Earth Sciences 39(5), 711–30. doi: 10.1139/e01-099.CrossRefGoogle Scholar
Krogh, TE (1994) Precise U-Pb ages for Grenvillian and pre-Grenvillian thrusting of Proterozoic and Archean metamorphic assemblages in the Grenville Front tectonic zone, Canada. Tectonics 13(4), 963–82. doi: 10.1029/94TC00801.CrossRefGoogle Scholar
Lasalle, S, Fisher, CM, Indares, A and Dunning, G (2013) Contrasting types of Grenvillian granulite facies aluminous gneisses: insights on protoliths and metamorphic events from zircon morphologies and ages. Precambrian Research 228, 117–30. doi: 10.1016/j.precamres.2013.01.014.CrossRefGoogle Scholar
Machado, N and Gauthier, G (1996) Determination of 207Pb/206Pb ages on zircon and monazite by laser-ablation ICPMS and application to a study of sedimentary provenance and metamorphism in southeastern Brazil. Geochimica et Cosmochimica Acta 60, 5063–73. doi: 10.1016/S0016-7037(96)00287-6.CrossRefGoogle Scholar
Marsh, JH, Jorgensen, TRC, Petrus, JA, Hamilton, MA and Mole, DR (2019) U-Pb, trace element, and Hafnium isotopic composition of the Maniitsoq zircon: a potential new zircon Archean reference material. Goldschmidt 2019, 2161, Abstracts. Barcelona.Google Scholar
Martignole, J, Machado, N and Nantel, S (1993) Timing of intrusion and deformation of the Riviere-Pentecote anorthosite (Grenville Province). Journal of Geology 101(5). doi: 10.1086/648256.CrossRefGoogle Scholar
Martin, C and Dickin, AP (2005) Styles of Proterozoic crustal growth on the southeast margin of Laurentia: evidence from the central Grenville Province northwest of Lac St.-Jean, Quebec. Canadian Journal of Earth Sciences 42, 1643–52. doi: 10.1139/e05-052.CrossRefGoogle Scholar
Moorbath, S, Whitehouse, MJ and Kamber, BS (1997) Extreme Nd-isotope heterogeneity in the early Archaean – fact or fiction? Case histories from northern Canada and West Greenland. Chemical Geology 135(3–4), 213–31. doi: 10.1016/S0009-2541(96)00117-9.CrossRefGoogle Scholar
Morel, MLA, Nebel, O, Nebel-Jacobsen, YJ, Miller, JS and Vroon, PZ (2008) Hafnium isotope characterization of the GJ-1 zircon reference material by solution and laser-ablation MC-ICPMS. Chemical Geology 225(1–2), 231–5. doi: 10.1016/j.chemgeo.2008.06.040.Google Scholar
Moukhsil, A, Solgadi, F, Belkacim, S, Augland, LE and David, J (2015) Géologie de la région de Parent, Haut-Saint-Maurice (partie ouest du Grenville). Québec: Ministère de l’Énergie et des Ressources naturelles. RG 2015-04, 59 pp. + 1 map (scale 1/125 000).Google Scholar
Moukhsil, A, Solgadi, F, Clark, T, Blouin, S, Indares, A and Davis, DW (2013) Géologie du nord-ouest de la région du barrage Daniel-Johnson (Manic 5), Côte-Nord (22K14, 22K15, 22K16, 22N02 et 22N03). Québec: Ministère des Ressources naturelles. RG 2013-01, 44 pp. + 1 plan.Google Scholar
Nebel, O, Nebel-Jacobsen, Y, Mezger, K and Berndt, J (2007) Initial Hf isotope compositions in magmatic zircon from early Proterozoic rocks from the Gawler Craton, Australia: a test for zircon model ages. Chemical Geology 241(1–2), 23–37. doi: 10.1016/j.chemgeo.2007.02.008.CrossRefGoogle Scholar
Papapavlou, K, Darling, JR, Lightfoot, PC, Lasalle, S, Gibson, L, Storey, CD and Moser, D (2018) Polyorogenic reworking of ore-controlling shear zones at the South Range of the Sudbury impact structure: a telltale story from in situ U–Pb titanite geochronology. Terra Nova 30(3), 254–61. doi: 10.1111/ter.12332.CrossRefGoogle Scholar
Papapavlou, K, Darling, JR, Storey, CD, Lightfoot, PC, Moser, DE and Lasalle, S (2017) Dating shear zones with plastically deformed titanite: new insights into the orogenic evolution of the Sudbury impact structure (Ontario, Canada). Precambrian Research 291, 220–35. doi: 10.1016/j.precamres.2017.01.007.CrossRefGoogle Scholar
Papapavlou, K, Strachan, RA, Storey, CD and Bullen, D (2021) Tectonic significance of a supra-ophiolitic sedimentary cover succession, Unst, Shetland, Scottish Caledonides: insights from the U–Pb–Hf detrital zircon record. Journal of the Geological Society 178(5). doi: 10.1144/jgs2020-169.CrossRefGoogle Scholar
Parman, SW (2015) Time-lapse zirconography: imaging punctuated continental evolution. Geochemical Perspectives Letters 1, 43–52. doi: 10.7185/geochemlet.1505.Google Scholar
Patchett, PJ and Tatsumoto, M (1981). A routine high-precision method for Lu-Hf isotope geochemistry and chronology. Contributions to Mineralogy and Petrology 75(3), 263–67. doi: 10.1007/BF01166766.CrossRefGoogle Scholar
Paton, C, Hellstrom, J, Paul, B, Woodhead, J and Hergt, J (2011) Iolite: freeware for the visualisation and processing of mass spectrometric data. Journal of Analytical Atomic Spectrometry 26, 2508–18. doi: 10.1039/c1ja10172b.CrossRefGoogle Scholar
Peck, WH, Quinan, MP and Selleck, BW (2019) Detrital zircon constraints on Grenville sedimentation at the margin of Laurentia. Precambrian Research 331, 105342. doi: 10.1016/j.precamres.2019.105342.CrossRefGoogle Scholar
Rivers, T (1997) Lithotectonic elements of the Grenville Province: review and tectonic implications. Precambrian Research 86(3–4), 117–54. doi: 10.1016/s0301-9268(97)00038-7.CrossRefGoogle Scholar
Rivers, T (2008) Assembly and preservation of lower, mid, and upper orogenic crust in the Grenville Province: implications for the evolution of large hot long-duration orogens. Precambrian Research 167, 237–59. doi: 10.1016/j.precamres.2008.08.005.CrossRefGoogle Scholar
Rivers, T (2021) The Grenvillian Orogeny and Rodinia. In Encyclopedia of Geology (eds-in-chief D Alderton and SA Elias), pp. 187201. Cambridge: Academic Press. doi: 10.1016/b978-0-12-409548-9.12506-0.CrossRefGoogle Scholar
Rondot, J (1978) Région du Saint-Maurice. Québec: Ministère de l’Énergie et des Ressources naturelles. DPV 594, 91 pp. + 2 maps.Google Scholar
Rudnick, RL and Gao, S (2003) Composition of the continental crust. In Treatise on Geochemistry (eds-in-chief HD Holland and KK Turekian) 3, 1–64. doi: 10.1016/B0-08-043751-6/03016-4.CrossRefGoogle Scholar
Sager-Kinsman, EA and Parrish, RR (1993) Geochronology of detrital zircons from the Elzevir and Frontenac terranes, Central Metasedimentary Belt, Grenville Province, Ontario. Canadian Journal of Earth Sciences 30, 465–73. doi: 10.1139/e93-034.CrossRefGoogle Scholar
Santos, MM, Lana, C, Scholz, R, Buick, I, Schmitz, MD, Kamo, SL, Gerdes, A, Corfu, F, Tapster, S, Lancaster, P, Storey, CD, Basei, MAS, Tohver, E, Alkmin, A, Nalini, H, Krambrock, K, Fantini, C and Wiedenbeck, M (2017) A new appraisal of Sri Lankan BB zircon as a reference material for LA-ICP-MS U-Pb geochronology and Lu-Hf isotope tracing. Geostandards and Geoanalytical Research 41, 335–58. doi: 10.1111/ggr.12167.CrossRefGoogle Scholar
Sappin, AA, Constantin, M, Clark, T and van Breemen, O (2009) Geochemistry, geochronology, and geodynamic setting of Ni-Cu ± PGE mineral prospects hosted by mafic and ultramafic intrusions in the Portneuf-Mauricie Domain, Grenville Province, Quebec. Canadian Journal of Earth Sciences 46(5), 331–53. doi: 10.1139/E09-022.Google Scholar
Schärer, U (1991) Rapid continental crust formation at 1.7 Ga from a reservoir with chondritic isotope signatures, eastern Labrador. Earth and Planetary Science Letters 102(2), 110–33. doi: 10.1016/0012-821X(91)90002-Y.CrossRefGoogle Scholar
Schulz, KJ and Cannon, WF (2007) The Penokean orogeny in the Lake Superior region. Precambrian Research 157, 425. doi: 10.1016/j.precamres.2007.02.022.CrossRefGoogle Scholar
Sláma, J, Košler, J, Condon, DJ, Crowley, JL, Gerdes, A, Hanchar, JM, Horstwood, MSA, Morris, GA, Nasdala, L, Norberg, N, Schaltegger, U, Schoene, B, Tubrett, MN and Whitehouse, MJ (2008) Plešovice zircon: a new natural reference material for U-Pb and Hf isotopic microanalysis. Chemical Geology 249(1–2), 1–35. doi: 10.1016/j.chemgeo.2007.11.005.CrossRefGoogle Scholar
Söderlund, U, Patchett, PJ, Vervoort, JD and Isachsen, CE (2004) The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions. Earth and Planetary Science Letters 219(3–4), 311–24. doi: 10.1016/S0012-821X(04)00012-3.CrossRefGoogle Scholar
Spencer, CJ, Cawood, PA, Hawkesworth, CJ, Prave, AR, Roberts, NMW, Horstwood, MSA and Whitehouse, MJ (2015) Generation and preservation of continental crust in the Grenville Orogeny. Geoscience Frontiers 6(3), 357–72. doi: 10.1016/j.gsf.2014.12.001.CrossRefGoogle Scholar
Stamatelopoulou-Seymour, K and MacLean, WH (1977) Geochemistry of possible metavolcanic rocks and their relationship to mineralization at Montauban-les-mines, Quebec. Canadian Journal of Earth Sciences 14(11), 2440–52. doi: 10.1139/e77-212.CrossRefGoogle Scholar
Stevenson, RK and Patchett, PJ (1990) Implications for the evolution of continental crust from Hf isotope systematics of Archean detrital zircons. Geochimica et Cosmochimica Acta 54, 1683–97. doi: 10.1016/0016-7037(90)90400-F.CrossRefGoogle Scholar
Thirlwall, MF and Anczkiewicz, R (2004) Multidynamic isotope ratio analysis using MC-ICP-MS and the causes of secular drift in Hf, Nd and Pb isotope ratios. International Journal of Mass Spectrometry 235(1), 59–81. doi: 10.1016/j.ijms.2004.04.002.CrossRefGoogle Scholar
Tucker, RD and Gower, CF (1994) A U-Pb geochronological framework for the Pinware Terrane, Grenville Province, southeast Labrador. Journal of Geology 102(1), 67–78. doi: 10.1086/629648.CrossRefGoogle Scholar
Turlin, F, Vanderhaeghe, O, Gervais, F, André-Mayer, AS, Moukhsil, A, Zeh, A and Ion Probe Team Nancy (I.P.T.N) (2019) Petrogenesis of LREE-rich pegmatitic granite dykes in the central Grenville Province by partial melting of Paleoproterozoic-Archean metasedimentary rocks: evidence from zircon U-Pb-Hf-O isotope and trace element analyses. Precambrian Research 327, 327–60. doi: 10.1016/j.precamres.2019.02.009.CrossRefGoogle Scholar
Van Breemen, O and Corriveau, L (2005) U-Pb age constraints on arenaceous and volcanic rocks of the Wakeham Group, eastern Grenville Province. Canadian Journal of Earth Sciences 42, 1677–97. doi: 10.1139/e05-079.CrossRefGoogle Scholar
Van Schmus, WR and Anderson, JL (1977) Gneiss and migmatite of Archean age in the Precambrian basement of central Wisconsin. Geology 5, 45–48. doi: 10.1130/0091-7613(1977)5<45:GAMOAA>2.0.CO;2.2.0.CO;2>CrossRefGoogle Scholar
Van Schmus, WR and Bickford, ME (1981) Proterozoic chronology and evolution of the Midcontinent region, North America. Developments in Precambrian Geology 4, 261–96. doi: 10.1016/S0166-2635(08)70016-1.CrossRefGoogle Scholar
Van Wyck, N and Johnson, CM (1997) Common lead, Sm-Nd, and U-Pb constraints on petrogenesis, crustal architecture, and tectonic setting of the Penokean orogeny (Paleoproterozoic) in Wisconsin. Bulletin of the Geological Society of America 109(7), 799–808. doi: 10.1130/0016-7606(1997)109<0799:CLSNAU>2.3.CO;2.2.3.CO;2>CrossRefGoogle Scholar
Vautour, S and Dickin, A (2019) Nd isotope mapping of a Pinwarian-age composite arc belt in the Quebecia terrane of the central Grenville Province, Canada. Precambrian Research 332(15). doi: 10.1016/j.precamres.2019.105409.CrossRefGoogle Scholar
Vermeesch, P (2012) On the visualisation of detrital age distributions. Chemical Geology 312–313(18), 190–4. doi: 10.1016/j.chemgeo.2012.04.021.Google Scholar
Vermeesch, P (2018) IsoplotR: a free and open toolbox for geochronology. Geoscience Frontiers 9(5), 1479–93. doi: 10.1016/j.gsf.2018.04.001.CrossRefGoogle Scholar
Vermeesch, P (2021) Maximum depositional age estimation revisited. Geoscience Frontiers 12, 843–50. doi: 10.1016/j.gsf.2020.08.008.CrossRefGoogle Scholar
Vervoort, JD and Blichert-Toft, J (1999) Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochimica et Cosmochimica Acta. doi: 10.1016/S0016-7037(98)00274-9.CrossRefGoogle Scholar
Vervoort, JD, Kemp, AI and Fisher, CM (2018) Hf Isotope Constraints on Evolution of the Depleted Mantle and Growth of Continental Crust. American Geophysical Union, Fall meeting 2018, abstract #V23A-07.Google Scholar
Wasteneys, HA, Kamo, SL, Moser, D, Krogh, TE, Gower, CF and Owen, JV (1997) U-Pb geochronological constraints on the geological evolution of the Pinware terrane and adjacent areas, Grenville Province, southeast Labrador, Canada. Precambrian Research. doi: 10.1016/S0301-9268(96)00030-7.CrossRefGoogle Scholar
Whitmeyer, SJ and Karlstrom, KE (2007) Tectonic model for the Proterozoic growth of North America. Geosphere 3, 220–59. doi: 10.1130/GES00055.1.CrossRefGoogle Scholar
Wiedenbeck, M, Allé, P, Corfu, F, Griffin, WL, Meier, M, Oberli, F, Von Quadt, A, Roddick, JC and Spiegel, W (1995) Three natural zircon standards for U-Th-Pb, Lu-Hf, trace element and REE analyses. Geostandards Newsletter 19(1), 1–23. doi: 10.1111/j.1751-908X.1995.tb00147.x.CrossRefGoogle Scholar
Wodicka, N, David, J, Parent, M, Gobeil, A and Verpaelst, P (2003) Géochronologie U-Pb et Pb-Pb de la région de Sept-iles – Natashquan, Province de Grenville, moyenne Cote-Nord. In Geologie et ressources minerales de la partie est de la Province de Grenville (eds D Brisebois and T Clark). Montréal: Ministère des Ressources naturelles, 59–117.Google Scholar
Woodhead, J, Hergt, J, Shelley, M, Eggins, S and Kemp, R (2004) Zircon Hf-isotope analysis with an excimer laser, depth profiling, ablation of complex geometries, and concomitant age estimation. Chemical Geology 209(1–2), 121–35. doi: 10.1016/j.chemgeo.2004.04.026.CrossRefGoogle Scholar
Woussen, G, Dimroth, E, Corriveau, L and Archer, P (1981) Crystallization and emplacement of the Lac St-Jean anorthosite massif (Quebec, Canada). Contributions to Mineralogy and Petrology 76, 343–50. doi: 10.1007/BF00375461.CrossRefGoogle Scholar
Zelek, M and Dickin, A (2013) Nd isotope mapping of crustal terranes in the parent-clova area, Quebec: implications for the evolution of the Laurentian margin in the Central Grenville Province. Geosciences (Switzerland) 3(3), 448–65. doi: 10.3390/geosciences3030448.Google Scholar
Figure 0

Fig. 1. Simplified geological map of the Lac-Saint-Jean region. The inset (a) depicts the crustal terranes of the Grenville Province sensu Dickin, 2000). The white dashed rectangle depicts the study area. The red circles in (b) denote the four sampling sites.

Figure 1

Fig. 2. Representative cathodoluminescence images from detrital zircon grains of the four quartzite samples. Red and blue circles show the location of the 25 μm and 50 μm beam location for U–Pb and Lu–Hf microanalysis, respectively.

Figure 2

Fig. 3. Kernel density estimation (KDE) 207Pb/206Pb age diagrams and weighted mean 207Pb/206Pb age diagrams from the youngest zircon grains in the four quartzite samples. The numbers adjacent to the peaks denote the age maxima in Ga. The weighted mean plots next to the KDE diagrams represent MDA estimates based on the youngest acceptable dates for each supracrustal complex.

Figure 3

Fig. 4. Radial plots and MDA estimates of the (a) Montauban Group quartzite, (b) Wabash Complex quartzite and (c, d) Barrois Complex quartzites using isoplotR (Vermeesch, 2018) and the Maximum Likelihood Age (MLA) algorithm (Vermeesch, 2021). The radial plots display heteroscedastic data (i.e. data points with unequal uncertainties) where each data point represents a 207Pb/206Pb date. These plots are scatter plots of standardized estimates (y axis) against precisions (x axis).The angular position of each data point denotes the date, whereas the horizontal distance from the point of origin varies with the precision of the analysis (i.e. imprecise measurements on the left side and more precise on the right side of the plot). For further details on age calculation see Vermeesch (2021) and Galbraith (1988).

Figure 4

Fig. 5. Epsilon Hf versus age with Lu–Hf isotopic data from the detrital zircon grains of the four quartzites reported in this study. Black and white coloured are Lu–Hf data points from the dataset of Spencer et al. (2015). The black-filled data points represent data from the pre-orogenic Siamarnekh Formation (SF). The white-filled data points depict, for comparative purposes, Lu–Hf isotopic data from pre-, syn- and post-orogenic formations from Labrador and Scotland. The lines with different slopes (DM = depleted mantle, MC = mafic crust, CC = continental crust) correspond to different Lu–Hf isotope evolution curves for the depleted mantle (176Lu/177Hf = 0.03976; Vervoort et al.2018), continental crust (176Lu/177Hf = 0.0113; Rudnick & Gao, 2003) and mafic crust (176Lu/177Hf = 0.022; Nebel et al.2007) reservoirs. The CHUR reference line (ϵHf = 0) involves an uncertainty of ±0.4 ϵHf units (Bouvier et al.2008). The grey band is drawn parallel to the 176Lu/177Hf = 0.022 line and covers c. 70 % of the 1.8–1.9 Ga zircon grains.

Figure 5

Fig. 6. (a) Tectonic setting (not in scale) of the examined quartzite from the supracrustal sequence of the Montauban Group; (b) KDE plots of U–Pb dates from the Montauban and Escoumins supracrustal group (Groulier et al.2018); (c) KDE plots of U–Pb dates from the Wakeham basin (U–Pb data from Wodicka et al.2003 and Van Breemen & Corriveau., 2005); (d) tectonic setting (not in scale) of the Wabash–Barrois supracrustal groups; (e) KDE plots of U–Pb dates from the Wabash–Barrois and Saint-Boniface supracrustal groups (U–Pb data from this study and Peck et al.2019); (f) KDE plots of U–Pb dates from supracrustal sequences of the Adirondack lowlands (U–Pb data from Chiarenzelli et al.2015; and pers. comm. with J. Chiarenzelli).

Supplementary material: File

Papapavlou et al. supplementary material

Papapavlou et al. supplementary material

Download Papapavlou  et al. supplementary material(File)
File 821.6 KB