Two Neoarchean tectonothermal events on the western edge of the North Atlantic Craton, as revealed by SIMS dating of the Saglek Block, Nain Province, Labrador

The Saglek Block forms the northern part of the Nain Province and underwent widespread metamorphism at c. 2.7 Ga, producing the dominant gneissosity and intercalation of supracrustal sequences. Zircon dating of gneiss samples collected along 80 km of the Labrador coast from Ramah Bay in the north to Hebron Fjord in the south confirms the widespread extent of high-grade metamorphism between 2750 and 2700 Ma. In addition, a distinct event between 2550 and 2510 Ma produced felsic melt with peritectic garnet in metavolcanic gneiss and granoblastic recrystallization in mafic granulite. Ductile deformation of granite emplaced at c. 2550 Ma indicates that this later event involved a degree of tectonism during high-T metamorphism. Such tectonism may be related to a hypothesized post-2.7 Ga juxtaposition of the predominantly Eoarchean Saglek Block against the Mesoarchean Hopedale Block, along a north–south boundary that extends from the coast near Nain to offshore of Saglek Bay. Evidence of reworking of c. 2.7 Ga gneisses by c. 2.5 Ga tectonothermal activity has been found elsewhere on the margins of the North Atlantic Craton, of which the Nain Province represents the western margin. In particular, a recent suggestion that c. 2.5 Ga metamorphic ages along the northern margin of the North Atlantic Craton in SW Greenland may record the final assembly of the craton could also apply to the western margin as represented by the rocks of the Nain Province. Supplementary material: Plots and geochemical data are available at https://doi.org/10.6084/m9.figshare.c.4567934

gneisses from outcrops around Saglek Bay differ in age across the Handy Fault (Fig. 1), with c. 3.6 Ga protoliths to the west and c. 3.3 Ga to the east. However, recent dating (Komiya et al. 2015;Kusiak et al. 2018;Sałacinśka et al. 2018Sałacinśka et al. , 2019 revealed complicated and tectonized relationships between Eoarchean and younger TTG gneisses on both sides of the fault. Metamorphosed supracrustal assemblages of sedimentary and volcanic rocks, with associated ultramafic gneisses, occur as discontinuous layers within gneissosity typically less than 100 m thick, and have been divided into sparsely distributed pre-or syn-Uivak supracrustal rocks (the Nulliak assemblage) and post-Uivak Meso-to Neoarchean Upernavik supracrustal rocks (Bridgwater & Schiøtte 1991). Isotopic U-Pb and Hf data from detrital zircon from Upernavik metasediments indicate deposition after c. 3.0 Ga, and it has been suggested that this unit includes unrelated supracrustal packages with various ages of deposition . Similarly, there is some uncertainty about the age of deposition of volcanic rocks and sediments that in part formed the Nulliak assemblage. Detrital zircons with an age of c. 3850 Ma (Nutman & Collerson 1991) supported deposition after that time, but Komiya et al. (2015) and Shimojo et al. (2016) favoured a greater antiquity, namely, >3.9 Ga, for these rocks. Whitehouse et al. (2019) have recently questioned this, suggesting that there is confusion over the assignment of metasedimentary and mafic gneisses to the Nulliak or Upernavik 'assemblages'. Graphite-bearing metasediments, claimed by Komiya et al. (2015) to belong to the Nulliak assemblage, contain detrital zircon that demonstrate a much younger provenance (see Whitehouse et al. 2019, after Schiøtte et al. 1992. Also, mafic tectonic enclaves near Nulliak Is. that were assigned to the Nulliak assemblage on the map by Ryan & Martineau (2012) have Sm-Nd isotopic signatures that show that some of these are much younger than the enclosing Uivak gneiss (Morino et al. 2017). In the absence of more extensive direct dating of tectonic enclaves and belts, the true age of many of the supracrustal rocks in the Saglek Block remains in question.
The assembly of the Saglek Block, comprising the Eoarchean Uivak gneiss, Mesoarchean tonalitic to gabbroic gneisses, and supracrustal packages, may be attributed to c. 2.7 Ga tectonism during the widespread high-T metamorphism Krogh & Kamo 2006;Kusiak et al. 2018;Sałacinśka et al. 2018). The above-cited studies were mostly focused around Saglek Bay, but c. 2.7 Ga metamorphic zircon was also identified as far south as Drachart Island , and as far north as the Avayalik Islands near the tip of the eastern Labrador pensinsula, where Scott (1995) proposed an extension of the Nain Province as reworked crust within the Torngat Orogen. A corresponding c. 2.7 Ga craton-forming event is recognized in similar gneisses on the conjugate section of southwestern Greenland (e.g. Nutman et al. 2004Nutman et al. , 2013Kirkland et al. 2018).
Post-dating the formation of the gneisses in the Saglek Block, 2.5 Ga mineral ages (U-Pb zircon, monazite and titanite, and K-Ar hornblende; see Discussion for references) have been attributed to the thermal and hydrothermal effects of post-tectonic granitic magmatism (Baadsgaard et al. 1979;Schiøtte et al. 1992). Alternatively, c. 2.5 Ga monazite and titanite ages from offshore drilling samples collected by Wasteneys et al. (1996), along with c. 2.7 Ga ages of detrital zircon, led Connelly & Ryan (1996) to infer a north-south-trending tectonic boundary between the Saglek and Hopedale blocks. The Hopedale Block comprises late Paleoarchean to Neoarchean protoliths metamorphosed at c. 3.0-2.8 Ga. The inferred boundary with the Saglek Block extends offshore north of Nain; but has not been directly observed, being obscured by the extensive Proterozoic Nain Plutonic Suite. Connelly & Ryan (1996) suggested a link between the Saglek-Hopedale boundary and the Okak shear zone (van Kranendonk & Helmstaedt 1990), to the south. The Okak shear deforms a granitic pluton on Okak Island (Fig. 1), which was inferred by Schiøtte et al. (1992) to have an age of c. 2.5 Ga (based on unpublished data of Roddick and van Kranendonk and the age of metamorphic monazite in adjacent metasediments).
More recent monazite dating by Kusiak et al. (2018) has increased the known extent of high-temperature metamorphism at both c. 2.7 and c. 2.5 Ga in the Saglek Block between Ramah Bay and Hebron Fjord. However, it is not clear whether these ages represent a prolonged period of high-T metamorphism, a gneissforming event with subsequent passive thermal activity and granite emplacement, or two discrete tectonothermal events. The purpose of this paper is to re-evaluate the timing and significance of deformation and metamorphism in the Saglek Block, utilizing new U-Pb isotopic dating of zircon and monazite.

Field relationships
This study investigates the sequence of deformation events in the Saglek Block, based on coastal field work between Ramah Bay and Hebron Fjord (Fig. 1), conducted by our team in 2014 and 2017, and the scheme by van Kranendonk & Helmstaedt (1990) for the North River-Nutak area, 100 km south of Saglek Bay. In almost all localities, pre-deformation relationships between the diverse rock types have been transposed into a high-strain gneissosity and a multistage deformational history has long been recognized (Morgan 1975;Schiøtte et al. 1990). Because gneissosity affects both Eoarchean and younger Archean TTG protoliths, as well as Mesoarchean Upernavik supracrustal rocks, and because published age data indicate widespread zircon and monazite growth during metamorphism at c. 2.7 Ga (Bridgwater & Schiøtte 1991), this may have been the tectonothermal event during which the Saglek Block was assembled. Commonly observed intrafolial folds of gneissic laminations are indications of high-strain ductile deformation (D 1 ) prior to that which produced the dominant gneissosity (D 2 ). The dominant gneissosity is in most places vertical to steeply dipping, with lesser domains of low-angle layering, such as that observed on the cliff face at Cape Uivak (Fig. 2a). The difference between D 1 and D 2 structures can be observed where S 2 flattens leucosome in gneisses with S 1 gneissosity and/or mineral foliation ( Fig. 2b-d). Elsewhere, S 1 has been transposed into a highstrain, moderately dipping to subvertical gneissic S 2 , which trends predominantly north-south. It is unknown whether D 2 significantly postdates D 1 , but because it transposes D 1 fabrics in late Mesoarchean to early Neoarchean Upernavik supracrustal rocks, as well as older TTG gneisses, and no intrusive rocks separate D 1 and D 2 , it is likely that they represent stages of a single tectonothermal event. Dating of gneisses with Eoarchean protoliths has established an additional, much older episode of high-grade metamorphism at c. 3.6 Ga (Sałacinśka et al. , 2019; however, no large-scale structures have been distinguished in the field that relate to this earlier event. In outcrops c. 100 km to the south of Saglek Bay, van Kranendonk & Helmstaedt (1990) described ductile thrusting (F 0 ) in the Upernavik supracrustal rocks and recumbent folding (F n+1 ) in both supracrustal rocks and TTG gneisses during high-T, high-P metamorphism. These are low-angle features, in contrast to the predominantly steep north-southtrending nature of S 2 gneissosity in most of the Saglek area. Although it is possible that these relate to the low-angle macrofold at Cape Uivak, the latter, along with recumbent folds on nearby Big Is., have been attributed either to nappes produced during late Archean, regional, asymmetric folding that generated the present map pattern (Bridgwater et al. 1975) or else to a separate recumbent folding event superimposed on the gneiss pattern produced by the aforementioned asymmetric folding (B. Ryan, pers. comm. 2019).
Post-D 2 structures tend to be localized. Upright minor folds are found with axial planes parallel to dominant gneissosity, and are interpreted as recording the waning stages of D 2 tectonism. There are abundant granitoid stocks, sills and dykes that cut the dominant gneissosity. Such granitoids have been classified by previous researchers (Bridgwater & Schiøtte 1991;Schiøtte et al. 1992) as 'post-tectonic', with metamorphism attributed to late syntectonic magmatism in the waning stages of Neoarchean tectonism. However, in many localities between Saglek Bay and Hebron Fjord, granitoid stocks and dykes are strongly deformed, especially on the islands and in eastern coastal regions. On Dog Is., coarse metagranite that intrudes Uivak gneiss has a steep S 3 foliation that is axial planar to open F 3 folds where granitic melt has pooled in fold noses (Fig. 2e). Intense L 3 defined by stretching of recrystallized fabrics in both pre-D 2 gneisses and post-D 2 granitoids indicates high rates of simple shear during the D 3 event. Dynamic recrystallization of granitoid produced augen gneiss 1 km to the east of St John's Harbour (Fig. 2f ), and the augen show alignment with coarse-grained biotite in quartz and orthoclase. This alignment is parallel to that in the matrix, where quartz, feldspar and biotite have recrystallized into an anastomosing S 3 foliation (Fig. 2g). This is consistent with progressive crystallization of the granitoid during a high-strain ductile event. Mylonite micro-shear zones cut across S 3 foliation. There is an increase in D 3 strain eastwards and southeastwards from Saglek Bay to the coast, with increasing development of F 3 meso-to-macro folding with variably plunging fold axes and intense L 3 stretching and recrystallization parallel to fold axes. Such features are possibly related to D n+3 structures described by van Kranendonk & Helmstaedt (1990) further south, which they attributed to a large amphibolite-facies shear zone that runs north-south along the coastal fringe of the Saglek Block south of Saglek Bay to Okak Island, where it deforms syntectonic granitic plutons assumed to have intruded at c. 2.5 Ga ). However, this shear zone involves retrogression of granulites to amphibolite-facies gneisses, whereas no such shear zone-related retrogression in association with D 3 structures is observed around Saglek Bay or in granulites around Hebron Fjord.

Sample selection and description
Several samples were collected for age determination between Ramah Bay and Hebron Fjord (Fig. 1b). Metamorphic grade at these localities varies from amphibolite to granulite facies (Ryan & Martineau 2012), albeit with varying degrees of later lower-grade overprinting, especially around Ramah Bay. The samples include felsic orthogneisses (L1414, L1488, L1489, L1491 and L1493), intermediate orthogneisses having the chemical characteristics of altered volcanic rocks (L1458 and L1487), mafic granulites (L1453 and L1490) and metapelitic gneiss (L1492). A sample of syn-D 3 granitoid (L1412) was also collected. Classification of orthogneisses is based on whole-rock geochemistry, using the total alkali v. silica diagram (Middlemost 1994) for igneous protoliths with <65 wt% SiO 2 , and the ternary normative feldspar classification of Barker (1979;after O'Connor 1965). Plots are provided with geochemical data in the supplementary material. Here, samples are briefly described according to structural relationship and locality (Fig. 1b). Mineral modes and major element geochemistry are presented in Tables 1 and 2, respectively.
Fine-to medium-grained grey felsic orthogneisses matching the description of Uivak I gneiss were collected from Reichel Head (L1491, L1493), Little Ramah Bay (L1488, L1489) and Dog Island (L1414). All have granoblastic fabrics with S 2 defined by millimetre-to centimetre-scale quartzofeldspathic layers (leucosome) and aligned biotite with or without hornblende. Leucosome consisting of intergrown quartz and feldspar (Fig. 3a) or quartz and mesoperthite (Fig. 2d), which is coarser than the granoblastic fabric in the host gneiss, has resulted from crystallization of partial melt, with minor recrystallization on grain margins providing evidence of limited subsequent deformation. Samples L1488 and L1489 are representatives of trondhjemitic orthogneisses from Little Ramah Bay that have, respectively, an abundance and a scarcity of nebulitic leucosome. These two samples were combined for the purpose of dating. The leucosome of sample L1414 felsic orthogneiss from Dog Island (Fig. 2b) is stromatic and the host gneiss varies from a patchily heterogeneous texture, interpreted as the recrystallization of a coarse-grained granitoid (L1414A), to a finer-grained, homogeneous pale grey gneiss with few laminations of leucosome (L1414B). Parts A and B were therefore dated separately.
Mafic samples were collected from S 2 layers hosted by TTG orthogneisses at Little Ramah Bay (L1490) and the south shore of Hebron Fjord (L1453). The latter was tentatively assigned by Ryan & Martineau (2012) to the Nulliak supracrustal assemblage; however, the outcrop also contains aluminous metasediments more typical of the Upernavik supracrustal rocks. Both samples have granoblastic texture with two-pyroxene-and hornblendebearing assemblages typical of mafic granulite generated from basaltic protoliths (Fig. 3b), but L1453 has a stronger foliation, with stromatic leucosome and associated garnet-biotite-rich selvages. A subvertical NW-trending S 2 gneissosity at the Hebron Fjord locality is crenulated by open to tight F 3 folds with SW-dipping axial planar S 3 defined by aligned biotite and NW-plunging axes. Some patches of garnet-leucosome truncate S 2 but are deformed by D 3 structures, indicating partial melting of mafic orthogneiss during both events.
The S 2 gneissosity and stromatic leucosome found in the orthogneisses is also present in sample L1492 of metapelite from Reichel Head (Fig. 2c). The sample is rich in graphite, similar to metasedimentary rocks that have been claimed to be early Eoarchean in age by Tashiro et al. (2017). Garnet poikiloblasts enclose S 2 -aligned flakes of biotite, graphite and sillimanite (Fig. 3c). Leucosome rich in K-feldspar and quartz is also flattened into S 2 , and the mineral assemblage is characteristic of granulitefacies metamorphism.
Pyroxene-quartz-bearing samples were taken from Little Ramah Bay (L1487) and Upernavik Island (L1458). The former has been described by Kusiak et al. (2018), whereas the latter is a typical orthopyroxene-garnet gneiss (Fig. 3d) found interlayered with aluminous and mafic gneisses that form the Upernavik 'assemblage' (Ryan & Martineau 2012). It is more aluminous than typical andesite, but unlike metapelite from the same locality, it contains abundant orthopyroxene, and is chemically characteristic of altered metavolcanic or volcanoclastic rocks, similar to Mesoarchean deposits at Qussuk and Storø in southwestern Greenland (Szilas et al. 2016(Szilas et al. , 2017. Composition grades across S 2 gneissic layers from orthopyroxene-plagioclase (L1458A) to garnet-biotite (L1458B) gneiss; however, the difference is in modal proportion only, and all phases are present in both rock types. Abundant leucosome is present in L1458, as S 2 -cutting layers with biotite-rich selvages and euhedral garnet poikiloblasts, as are commonly formed through incongruent melting of pelitic rocks (L1458C; Fig. 3d). Garnet is anhedral and slightly poikilitic in both types, with large inclusions of rutile and quartz (Fig. 3e). The leucosome is quartz-rich and moderately deformed, with warped and recrystallized quartz grains wrapping garnet poikiloblasts that contain S 3aligned grains of biotite and monazite (Fig. 2h). The presence of coarse biotite flakes aligned with the foliation in biotite-rich selvages, and with the sub-mylonitic recrystallization of quartz in the leucosome, supports the interpretation that garnet formed through incongruent melting of the host gneiss, and that the melt crystallized under stress.

Methods
Detailed analytical protocols and data reduction procedures are provided in the Appendix. Determination of bulk-rock geochemistry for major elements was undertaken by Acme Labs in Vancouver, Canada, through Bureau Veritas, Poland. For Zr-inrutile thermometry, electron microprobe (EMP) analysis was carried out on a Cameca SX-100 instrument at the Electron Microprobe Laboratory, State Geological Institute of Dionýza Štúra, Bratislava, Slovakia. For ion microprobe analysis, plugs drilled from polished thin sections, and monazite and zircon mineral grains separated from crushed samples, were mounted in epoxy, polished and imaged by scanning electron microscope (SEM) with backscattered electron (BSE) and cathodoluminescence (CL) detectors at the John de Laeter Centre, Curtin University, Western Australia. Isotopic U-Pb dating of zircon and monazite grains was by sensitive highresolution ion microprobe (SHRIMP II) at the John de Laeter Centre, Curtin University in Perth, Western Australia, and by CAMECA IMS 1280 ion microprobe at the NordSIMS facility, Swedish Museum of Natural History, Stockholm. All ion microprobe data are quoted with 1σ analytical errors, whereas weighted mean and discordia intercept ages are quoted at 95% confidence levels, and include the decay-constant error of the concordia curve.

Rutile thermometry
The abundance of rutile in sample L1458 allowed for in situ EMP analysis of zirconium contents to estimate temperatures of mineral growth. The formulation of Watson et al. (2006) was used; data are    provided in Table 3. Three grains of rutile included in garnet from the garnet-biotite-rich part of the sample (L1458B, Fig. 3f

Zircon dating
Zircon grains from all samples are subhedral to anhedral, with cathodoluminescence imaging (Fig. 4) revealing rims with anhedral, graduated or sector zoning typical of growth under     Values corrected for common Pb. Disc. % = (1-( high-grade metamorphic conditions. Cores having euhedral, graduated and/or oscillatory growth zoning, typical of crystallization from an evolving magma, were found in all samples except L1453 (Hebron Fjord) and L1490 (Little Ramah Bay), the latter two having rounded or irregular cores without distinct growth zoning. Sub-grain domains of zircon with features typical of growth during metamorphism were targeted for spot analysis in all samples, and zircon grains with magmatic growth features were targeted in samples L1412 (near St Johns Harbour), L1414A/B (Dog Island), L1487 (Little Ramah Bay) and L1491 (Reichel Head). Isotopic U-Pb data (Table 4) are presented in Tera-Wasserburg concordia plots (Fig. 5) along with 207 Pb/ 206 Pb mean ages for concordant populations and Model 1 discordia intercept ages for linear arrays. Older outliers from rounded or irregular cores, which are interpreted as xenocrystic or inherited zircon, were not included in the calculation of ages and statistics from the identified populations. For cores and grains with growth zoning characteristic of igneous zircon, Model 1 discordia chords were calculated with forced lower intercepts of 2720 ± 50 Ma, approximating the time period within which granulite-facies gneisses were estimated to have formed from older magmatic protoliths in the Saglek Block , and references therein). Statistical test values (mean square of weighted deviates; MSWD) and other details are provided with the concordia plots in Figure 5. For andesitic orthogneiss L1487 and trondhjemitic orthogneiss L1491, discordia chords yield upper intercept ages of 3664 ± 35 Ma and 3715 ± 26 Ma, respectively. The latter includes five concordant data with a mean 207 Pb/ 206 Pb age of 3714 ± 11 Ma. Data from zircon in trondhjemitic orthogneiss L1414B spread between c. 3650 and 3590 Ma. Mean 207 Pb/ 206 Pb ages were also derived from igneous zircon in meta-trondhjemite layer L1414A (2749 ± 3 Ma) and meta-monzonite L1412 (2547 ± 3 Ma). In all cases, the estimates are interpreted as the time of crystallization of igneous protoliths, with the exception of L1414B; in that sample, the cluster of analyses at c. 3590 Ma can be interpreted as a maximum age only for the protolith, assuming that they were not disturbed by later metamorphism. Igneous zircon from metapelite L1492 yielded scattered ages between c. 3280 and 2950 Ma, which are interpreted as dating detrital sources for the metasediment, although here again, the possibility of disturbance at c. 2.7 Ga cannot be discounted.
Zircon with metamorphic morphologies, either as rims with discordant boundaries to cores or as distinctly equant rounded and sector-zoned grains, yielded statistically valid (MSWD ≤ 1.3) mean 207 Pb/ 206 Pb ages for samples of andesitic orthogneiss L1487 (2742 ± 8 Ma), trondhjemitic orthogneiss L1488 (2750 ± 7 Ma) and mafic granulite L1490 (2739 ± 9 Ma). Slightly more scattered data were derived from metapelite L1492 (c. 2750-2720 Ma) and metagranite L1493 (c. 2740L1493 (c. -2710. Two data from light-CL rims in zircon from sample L1414A yielded c. 2710 Ma ages. Together, these data from six samples are interpreted as dating zircon growth during high-T metamorphism between c. 2750 and 2710 Ma. The dataset from mafic granulite sample L1453 is more complex, with 51 analyses from unzoned, concentric and sectorzoned grains and cores ranging between c. 2740 and 2680 Ma (group 1 ages, Fig. 5), and 25 data from unzoned or gradationally zoned rims ranging between c. 2570 and 2510 Ma (group 2 ages). Analyses in group 1 record variable U contents (Fig. 5), whereas those from group 2 have uniformly low U contents. The groups represent periods of zircon growth during two separate metamorphic events. To better define the gap in time between the events, subsets of statistically equivalent data were extracted from the youngest ages in group 1 and the oldest ages in group 2. The 42 youngest out of 51 data in group 1 yield a mean 207 Pb/ 206 Pb age of 2702 ± 2 Ma, and the 20 oldest data out of 25 in group 2 yield a mean 207 Pb/ 206 Pb age of 2551 ± 6 Ma. These mean ages provide statistically robust estimates for the minimum age of zircon growth in the first metamorphic event, and the maximum age of growth in the second, respectively.

In situ monazite dating
To constrain the timing of mineral growth in high-grade metamorphic assemblages, plugs containing monazite and surrounding minerals were drilled from polished thin sections, mounted and analysed by secondary ion mass spectrometry (SIMS; Table 5, Fig. 5). Monazite in andesitic orthogneiss L1487 occurs as xenoblastic grains in a granoblastic assemblage that is strongly retrogressed, as described by Kusiak et al. (2018). Owing to the marginal alteration of monazite grains in thin section (Fig. 6a), data were also taken from unaltered fragments of monazite separated from the orthogneiss and mounted in a polished epoxy plug. Excluding three slightly younger, discordant data points, 12 analyses from a combination of grains in drilled thin sections and separates yielded a mean 207 Pb/ 206 Pb age of 2709 ± 14 Ma. Monazite in metapelite sample L1492 is unaltered and has polygonal grain boundaries with other metamorphic phases (Fig. 6b), and yields a mean 207 Pb/ 206 Pb age of 2727 ± 6 Ma. Ages were also collected from each part of metavolcanic rock L1458 (A, B and C). Two 10 μm wide monazite inclusions in garnet porphyroblasts from the garnet-biotite-rich part (L1458B) yield spot ages of c. 2680 and c. 2670 Ma. Monazite occurs more abundantly in association with garnet-leucosome L1458C, in which millimetre-scale preferentially aligned inclusions in garnet poikiloblasts are parallel to S 3 defined by biotite inclusions (Fig. 6c). Four inclusions of monazite yield ages that range from 2550 to 2510 Ma. Excluding the two oldest analyses, nine data yield a mean 207 Pb/ 206 Pb age of 2522 ± 7 Ma, which can be considered as a robust statistical minimum age for the period of metamorphism. Age data were also obtained from monazite separated from the orthopyroxene-rich part (L1458A) and yielded a mean 207 Pb/ 206 Pb age of 2551 ± 5 Ma. Mean ages from samples L1487 and L1492 are attributed to monazite growth during the first period of high-T metamorphism. Monazite from sample L1458, dated as inclusions in garnet from parts B and C, indicates multiple stages of mineral growth. Those grains present in cross-cutting garnet-leucosome (L1458C) fall within the second period of mineral growth at 2.5 Ga identified in zircon from other samples in this study, as does monazite in the matrix of part A. The two monazite inclusions in the garnet-biotite-rich part (B) fall between the two stages of zircon growth in other samples, but agree with some monazite age estimates obtained by Kusiak et al. (2018). This may be an indication of monazite growth and/or disturbance continuing after 2700 Ma, but as a separate generation from the second stage of growth from 2550 to 2510 Ma.

Significance and correlation between Labrador and Greenland
The new results from monazite and zircon associated with metamorphic assemblages and deformation fabrics, especially where supported by dating structurally constrained meta-granitoids, provide evidence of two distinct high-temperature tectonothermal events: at 2750-2700 Ma and 2550-2510 Ma (Fig. 7). The separation of the two episodes of high-T mineral growth is clearer than that observed in EMP monazite dating by Kusiak et al. (2018). There is evidence of partial melting and crystallization of anatectic melt in both the earlier and later stages of each of the events. Therefore, the growth of zircon and monazite after 2550 Ma is probably not due to the 'thermal effects' of granitic emplacement, as suggested by Schiøtte et al. (1992); rather, it is more likely that Fig. 5. Tera-Wasserburg concordia plots of SIMS age data from zircon and monazite. Spot ages represented by error ellipses are colour-coded according to type: green, c. 2.7 Ga metamorphic; blue, c. 2.5 Ga metamorphic; red or pink, igneous; grey, xenocrystic or inherited; tan, detrital. Values corrected for common Pb. Disc. % = (1-( tectonothermal activity is the progenitor of granitic melts that were emplaced both during and after high-strain deformation. These include the c. 2530 Ma major granitic stockworks described by Baadsgaard et al. (1979) on the coast and islands outside Saglek Bay. A re-examination of localities in which c. 2.5 Ga magmatism and mineral growth occurs shows that such ages are scattered along the Saglek Block from Saglek Bay to Nain, and on both sides of the Handy Fault (Fig. 8). Further north there is a lack of data; however, zircon growth during metamorphic events at c. 2.7 and 2.5 Ga have been recognized by Scott (1995) in meta-tonalites at Home and Avayalik Islands, which may be part of the Nain Province. Nevertheless, it is likely that the effects of the 2.5 Ga event increase towards the south and east, as such ages were also obtained from zircon and monazite in drill-core samples taken c. 40 km outside Saglek Bay see Fig. 8 Wasteneys et al. 1996). However, the presence of c. 2680 Ma granitoid sheets cutting folded and metamorphosed Lister gneiss constrains gneiss formation to c. 2.7 Ga (Schiøtte et al. 1989). This, along with the evidence for c. 2.7 and c. 2.5 Ga metamorphism in samples from the Saglek area in our study and in that of Kusiak et al. (2018), does not contradict the terrane boundary proposed by Wasteneys et al. (1996), but does suggest that the assembly of the Saglek and Hopedale blocks was earlier than c.
2.5 Ga. The late Archean metamorphic events and the assembly of two different crustal blocks in northern Labrador may be analogous to the juxtaposition of terranes having differing structural and metamorphic histories in the Archean of southwestern Greenland    (Nutman & Friend 2007;Friend & Nutman 2019). Major late Archean terrane boundaries along the coast of southwestern Greenland tend to run NE into the glacial cap, rather than following the general north-south trend of gneisses in the Saglek area. Circa 2.8-2.7 Ga high-grade metamorphism that strongly affects the Nain Province also does so in the vicinity of Nuuk, with grade decreasing towards the east (e.g. Nutman & Friend 2007;Dziggel et al. 2017). This part of the North Atlantic Craton contains a complex mixture of Eoarchean and Paleoarchean terranes, and Mesoarchean arc assemblages, and has many similarities in timing and composition to the gneisses of the Saglek Block. The extensive c. 2560 Ma Qôrqut Granite Complex that intrudes gneisses inland from Nuuk Naeraa et al. 2014) is a potential correlative of syn-to late-D 3 magmatism in the Saglek block, as marginal tectonic reworking of the Qôrqut has been observed ). However, no significant granite metamorphic event at c. 2.5 Ga is recognized in this part of Greenland. Some 100 km to the north, dating of c. 2.5 Ga metamorphic monazite in gneisses near Maniitsoq and inland led Dyck et al. (2015) to define the Majorqaq Belt, a NE-trending mobile belt between the main 2.7 Ga assembled part of the North Atlantic Craton and the Mesoarchean Maniitsoq block further north. Dyck et al. (2015) suggested that the belt resulted from the collision of the Maniitsoq block subsequent to southward subduction of an ocean basin beneath the North Atlantic Craton, and that the Qôrqut Granite Complex is the product of slab dewatering. The Majorqaq Belt might well correlate with c. 2.5 Ga tectonothermal activity in the Saglek Block. In this case, the Qôrqut Granite Complex would correlate well with large plutons of the same age found to the south of Saglek in the Okak area . Indeed, the Maniitsoq block itself was subjected to marginal reworking to the north by Paleoproterozoic tectonism, similar to the northern and western margins of the Nain Province (St-Onge et al. 2009). However, there is a lack of data from the Labrador coast north of Ramah Bay that would allow any clear correlation to be made with the hypothesized 'Majorqaq' Belt. In addition, this would dissociate the Eoarchean Uivak gneiss from terranes of similar age in the Itsaq Gneiss Complex (Fig. 8). The latter has complicated ductile structural relationships with Paleo-to Mesoarchean terranes (Hoffmann et al. 2014). If such relationships are similarly complicated in the Saglek Block, the intercalation of older and younger crust there may also be a product of amalgamation at c. showing the location of magmatic and metamorphic mineral ages (U-Pb zircon, monazite and titanite ages, and K-Ar hornblende ages) that fall within the c. 2.7 and 2.5 Ga events, and the approximate known extent of c. 2.7 Ga metamorphism in the Nain Province and southwestern Greenland. A hypothetical boundary between the Hopedale and Saglek blocks, slightly modified from that proposed by Connelly & Ryan (1996) and Wasteneys et al. (1996), is shown by dashed line, parallel to the Handy Fault. Terranes in southwestern Greenland (bound by black lines) and the extent of Proterozoic orogenic fronts (thick black lines) are from Henriksen et al. (2009). Age data are from Wanless et al. (1970Wanless et al. ( , 1974, Barton in zircon and/or micro-inclusions. Very low amounts of common Pb were detected during the spot analyses with (<0.1% of total 206 Pb), in many cases below detection limit for 204 Pb based on the electron multiplier background. Where common Pb corrections were deemed necessary on the basis of measurable 204 Pb (>3× standard deviation on the average background), these were small and therefore insensitive to the precise composition of common Pb. Data reduction was performed using the NordSIMS-developed suite of software of M. J. Whitehouse. All ion microprobe data are quoted with 1σ analytical errors, whereas weighted mean and discordia intercept ages are quoted at 95% confidence levels, and include the decay-constant error of the concordia curve.
For in situ (i.e. within polished thin section) and grain mount monazite analysis, the SHRIMP II was operated with a primary beam of O 2 ions focused through a 50 μm Köhler aperture to produce an oval 10 μm wide spot with a surface current of 0.2-0.4 nA. Secondary ionization was measured without energy filtering on a single electron multiplier on 13 mass stations from 202 (LaPO 2 ) to 270 (UO 2 ), with a mass resolution of >5200 for the latter. Secondary ion retardation was used to eliminate ion scatter. Mass stations 202 (LaPO 2 ), 203 (CePO 2 ), 205.9 (NdPO 2 ), 232 (Th), 244.8 (YCeO) and 264 (ThO 2 ) were analysed for matrix corrections and interference on 204 Pb, following the protocols outlined by Fletcher et al. (2010). Mass stations were measured through six cycles, with typical count times of 10 s per cycle for 204 Pb, background (at 204.04 a.m.u.) and 206 Pb, 30 s for 207 Pb and 5 s for 208 Pb. Reduction of raw data for standards and samples was performed using the SQUID 2.5 and Isoplot 3.70 add-ins for Microsoft Excel 2003 (Ludwig 2001(Ludwig , 2009). Age ( 206 Pb/ 238 U) and abundance of U were calibrated against reference monazite French (514 Ma; 1000 ppm U). High La and high Y-Nd-U standards Z2234 and Z2908, respectively, were used for matrix and interference corrections, following the method described by Fletcher et al. (2010). Corrections for common Pb on isotopic U/ Pb values and ages were carried out with common Pb estimated from 204 Pb counts and the composition of Broken Hill lead.
For Zr-in-rutile thermometry, electron microprobe analysis was undertaken at the Electron Microprobe Laboratory, State Geological Institute of Dionýza Štúra, Bratislava, Slovakia, utilizing a Cameca SX-100 electron microprobe equipped with four wavelengthdispersive spectrometers. Large high-sensitivity, LPET and LLIF crystals and a conventional TAP crystal were used for analysis. Analytical conditions were chosen to balance the best analytical conditions against reasonable acquisition times. An accelerating voltage of 15 kV was used, with a probe current of 200 nA. Zirconium contents were calibrated against an in-house standard.

Scientific editing by Randall Parrish
Correction notice: The copyright has been updated to Open Access.